Next Article in Journal
A Theoretical and Experimental Investigation on the Fracture Mechanism of Center-Symmetric Closed Crack in Compacted Clay under Compression–Shear Loading
Next Article in Special Issue
Application of Symmetric Explicit Symplectic Integrators in Non-Rotating Konoplya and Zhidenko Black Hole Spacetime
Previous Article in Journal
Grid Harmonics Suppression for Three Phase Dual-Frequency Grid-Connected Inverter Based on Feedforward Compensation
Previous Article in Special Issue
Contrasting the Implicit Method in Incoherent Lagrangian and the Correction Map Method in Hamiltonian
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spin-Gravity Coupling in a Rotating Universe

1
Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA
2
School of Astronomy, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran
3
Department of Physics, Sharif University of Technology, Tehran 11365-9161, Iran
4
Theoretical Physics Laboratory, Nuclear Safety Institute, Russian Academy of Sciences, B. Tulskaya 52, 115191 Moscow, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(8), 1518; https://doi.org/10.3390/sym15081518
Submission received: 15 June 2023 / Revised: 17 July 2023 / Accepted: 28 July 2023 / Published: 1 August 2023
(This article belongs to the Special Issue Symmetry in Gravity Research)

Abstract

:
The coupling of intrinsic spin with the nonlinear gravitomagnetic fields of Gödel-type spacetimes is studied. We work with Gödel-type universes in order to show that the main features of spin-gravity coupling are independent of causality problems of the Gödel universe. The connection between the spin–gravitomagnetic field coupling and Mathisson’s spin-curvature force is demonstrated in the Gödel-type universe. That is, the gravitomagnetic Stern–Gerlach force due to the coupling of spin with the gravitomagnetic field reduces in the appropriate correspondence limit to the classical Mathisson spin-curvature force.

1. Introduction

Inertia is the intrinsic tendency of matter to remain in a given condition. The state of matter in spacetime is determined by its mass and spin; indeed, mass and spin characterize the irreducible unitary representations of the Poincaré group [1]. Therefore, mass and spin determine the inertial properties of a particle. In classical physics, the inertial forces that act on a particle are proportional to its inertial mass; moreover, the moment of inertia is the rotational analogue of mass. The inertial effects of intrinsic spin are independent of the inertial mass of the particle and depend purely on intrinsic spin. The inertia of intrinsic spin is of quantum origin, and its properties, therefore, complement the inertial characteristics of the mass and orbital angular momentum of the particle.
It turns out that the intrinsic spin  S  of a particle couples to the rotation of a noninertial observer, thus resulting in a Hamiltonian of the form  H s r = S · Ω , where  Ω  is the angular velocity of the observer’s local spatial frame with respect to a nonrotating (i.e., Fermi–Walker) transported frame. For an intuitive explanation of this type of coupling, let us consider a noninertial observer that is at rest in Minkowski spacetime but refers its observations to axes that rotate uniformly with angular speed  Ω  in the positive sense about the direction of propagation of a plane electromagnetic wave of frequency  ω > Ω . The Fourier analysis of the electromagnetic field detected by the noninertial observer reveals that the measured frequency of the wave is given by  ω Ω , where the upper (lower) sign refers to positive (negative) helicity radiation. One can understand this result as a kind of “rotational Doppler effect”: In a positive (negative) helicity electromagnetic wave, the electric and magnetic fields rotate in the positive (negative) sense with the wave frequency  ω  about the direction of propagation. The noninertial observer thus realizes that the positive (negative) helicity radiation has electric and magnetic fields that rotate in the positive (negative) sense with frequency  ω Ω  ( ω + Ω ) about the direction of wave propagation. Multiplication of the measured frequency by results in the measured energy by the noninertial observer, namely,  ω Ω , which illustrates the coupling of photon helicity with rotation. A general consequence of spin-rotation coupling should be noted here: There is a certain shift in energy when polarized radiation passes through a rotating spin flipper. To demonstrate this effect within the context of the present discussion, imagine that the noninertial observer is replaced by a uniformly rotating half-wave plate. That is, the electromagnetic radiation of frequency  ω in  is normally incident on the plate rotating with  Ω < ω in . The frequency of the radiation within the stationary medium of the half-wave plate remains constant and approximately equal to  ω in Ω , since we have neglected time dilation for simplicity. The outgoing radiation has an opposite helicity to the incident radiation and frequency  ω out , where  ω in Ω ω out ± Ω  due to helicity–rotation coupling. Therefore,  ω out ω in 2 Ω  and the photon energy in passing through the rotating half-wave plate is shifted by  2 Ω .
A general account of spin-rotation coupling is contained in [2], and more recent discussions of its observational basis can be found in [3,4,5,6]. A similar phenomenon occurs in a gravitational field [7,8,9]. The spin-rotation effect can be theoretically extended to the spin–gravity coupling via the gravitational Larmor theorem [10,11], which is the rotational side of Einstein’s principle of equivalence. Imagine a free test gyroscope with its center of mass held at rest in a gravitational field; then, the locally measured components of the gyroscope’s spin vector undergo a precessional motion with an angular velocity that is given by the locally measured gravitomagnetic field. The Gravity Probe B (GP-B) space experiment has measured the gravitomagnetic field of the Earth [12,13].
According to the gravitational Larmor theorem, the gravitomagnetic field of a rotating system is locally equivalent to a rotation resulting in a Hamiltonian for intrinsic spin–gravity coupling of the form  H s g = S · B , where  B  is the relevant gravitomagnetic field [14]. The spin–gravity coupling is of basic physical significance due to the fundamental nature of the intrinsic spin of the particles and the universality of the gravitational interaction. For prospects regarding the measurement of intrinsic spin–gravity coupling, see [15,16,17,18,19]. In general,  B  depends on the position, and the intrinsic spin–gravity coupling leads to a measured gravitomagnetic Stern–Gerlach force of the form  ( S · B ) . This gravitational force, which acts on a test particle, is completely independent of its inertial mass and depends solely on its intrinsic spin. It has been shown in [20], within the framework of linearized general relativity, that the gravitomagnetic Stern–Gerlach force associated with spin–gravity coupling reduces in the correspondence limit to Mathisson’s classical spin-curvature force [21,22]. It would be interesting to extend this result to the nonlinear regime. The purpose of the present work is to study further the inertial effects of intrinsic spin by investigating the intrinsic spin–gravity coupling for spinning test particles in Gödel-type spacetimes. For background material, ref. [20] and the references cited therein should be consulted for further important information regarding the topic of spin-rotation–gravity coupling and its experimental basis.

2. Gravitomagnetism in the Gödel-Type Universe

With respect to the spacetime coordinates  x μ = ( c t , x , y , z ) , the metric of the Gödel solution [23] of Einstein’s gravitational field equations arises as a special case in the class of the so-called Gödel-type models [24,25,26], which are described by the following line element:
d s 2 = g μ ν d x μ d x ν = d t 2 2 σ e μ x d t d y + d x 2 + κ e 2 μ x d y 2 + d z 2 ,
with arbitrary constant parameters  μ , σ , and  κ . In our conventions, the speed of light  c = 1  and Planck’s constant  = 1 , unless specified otherwise; moreover, the metric signature is +2, and the Greek indices run from 0 to 3, while the Latin indices run from 1 to 3. The system of coordinates in metric (1) is admissible provided
σ + κ > 0 .
Moreover, we assume throughout that  σ > 0 . In general, the Gödel-type universe contains closed timelike curves, which could lead to problems with causality. However, one can demonstrate [24] that closed timelike curves are absent in model (1), provided
κ 0 .
Specifically, for the Gödel universe,  κ = 1  in metric (1); therefore, closed timelike curves do exist in the Gödel universe. To ensure that our considerations regarding spin–gravity coupling are independent of the causality difficulties of the Gödel universe, we use metric (1) for our main calculations in this paper.
The Gödel-type universe is a regular stationary and spatially homogeneous spacetime that contains rotating matter. Consider the class of observers that are all spatially at rest in this spacetime. Each such observer has a velocity 4-vector  u μ = δ 0 μ  that is free of acceleration, expansion, and shear; however, it is rotating in the negative sense about the z axis and its vorticity 4-vector
ω μ = 1 2 η μ ν ρ σ u ν u ρ ; σ ,
is purely spatial  ω μ = ( 0 , ω ) , with the 3-vector
ω = Ω z , Ω = μ 2 σ σ + κ .
For the sake of definiteness, we henceforth assume that  Ω > 0 ; then, Equation (5) implies that  μ > 0  as well. Here,  η α β γ δ = ( g ) 1 / 2 ϵ α β γ δ  is the Levi–Civita tensor, and  ϵ α β γ δ  is the alternating symbol with  ϵ 0123 = 1 . It is interesting to note that, in nonrelativistic fluid mechanics, the vorticity vector  ω N  is defined as  ω N = × v , where  v  is the flow velocity. If the fluid rotates with a spatially uniform angular velocity  Ω  such that  v = Ω × x , then  ω N = 2 Ω . In this paper, we follow the relativistic definition of vorticity.
The geometry of the Gödel-type model has been studied by a number of authors [27,28,29]. The Weyl curvature of Gödel-type spacetime is of type D in the Petrov classification. The Gödel-type universe admits five Killing vector fields, namely,  t y z x μ y y  and [24,30]:
K = 2 σ e μ x σ + κ t 2 μ y x + μ 2 y 2 e 2 μ x σ + κ y .
We are interested in the measurements of an observer that is free and spatially at rest in spacetime with a 4-velocity vector  u μ = d x μ / d τ  and proper time  τ , where  τ = t + constant . The observer carries along its geodesic world line a natural tetrad frame  e μ α ^  that is orthonormal, namely,
g μ ν e μ α ^ e ν β ^ = η α ^ β ^ ,
where  η μ ν = diag ( 1 , 1 , 1 , 1 )  is the Minkowski metric tensor. Indeed,
e 0 ^ = t , e 1 ^ = x , e 2 ^ = σ σ + κ t + e μ x σ + κ y , e 3 ^ = z ,
where the spatial axes of the observer’s frame are primarily along the background coordinate axes. By introducing the dual coframe  ϑ α ^ ,
ϑ 0 ^ = d t + σ e μ x d y , ϑ 1 ^ = d x , ϑ 2 ^ = σ + κ e μ x d y , ϑ 3 ^ = d z ,
such that  e α ^ ϑ β ^ = δ α β , the line element (1) is recast into
d s 2 = d t + σ e μ x d y 2 + d x 2 + ( σ + κ ) e 2 μ x d y 2 + d z 2 .
Let  λ μ α ^  be the orthonormal tetrad frame that is parallel transported along the observer’s geodesic world line such that  D λ μ α ^ / d τ = 0 . We find that
λ μ 1 ^ = e μ 1 ^ cos Ω τ + e μ 2 ^ sin Ω τ , λ μ 2 ^ = e μ 1 ^ sin Ω τ + e μ 2 ^ cos Ω τ ,
while  λ μ 3 ^ = e μ 3 ^ , and naturally  λ μ 0 ^ = e μ 0 ^ = u μ . It is simple to check these results using the Christoffel symbols:
Γ 10 0 = σ σ + κ Ω , Γ 20 1 = σ + κ e μ x Ω , Γ 10 2 = e μ x σ + κ Ω ,
which are the only nonzero components of  Γ ν 0 μ . Therefore, the observer’s natural frame rotates with respect to the parallel-transported frame about their common z axis with frequency  Ω , which is consistent with vorticity (5).
Let us now consider the special case of metric (1) with parameters
μ = 2 Ω , σ = 2 , κ = 1 .
With these parameters, metric (1) reduces to the Gödel line element
d s 2 = d t 2 2 2 e 2 Ω x d t d y + d x 2 e 2 2 Ω x d y 2 + d z 2 .
For the Gödel universe, Einstein’s field equations
R μ ν 1 2 g μ ν R + Λ g μ ν = 8 π G T μ ν
have a perfect fluid source of
T μ ν = ( ρ + p ) u μ u ν + p g μ ν ,
where  ρ  is the energy density, p is the pressure, and  u μ = δ 0 μ  is the 4-velocity vector of the perfect fluid. In this special case,  R μ ν = 2 Ω 2 u μ u ν , and
2 Ω 2 = 8 π G ( ρ + p ) , Λ + Ω 2 = 8 π G p .
In the absence of the cosmological constant  Λ , we have as the source of the Gödel universe a perfect fluid with a stiff equation of state  ρ = p = Ω 2 / ( 8 π G ) . Another possibility is dust ( p = 0 ) with  4 π G ρ = Λ = Ω 2 . It follows from Equation (17) that  Λ = 4 π G ( ρ p ) ; therefore, in any realistic situation, the cosmological constant of the Gödel universe must be negative or zero ( Λ 0 ).
The spinning test particle in the Gödel universe is immersed in the perfect fluid source, and its intrinsic spin couples to the vorticity of the fluid. The nature of the spin–gravity coupling and its connection with Mathisson’s classical spin-curvature force provided the original motivation for the present work.
After this brief digression regarding the Gödel universe, we return to the Gödel-type metric with explicit components:
( g μ ν ) = 1 0 σ W 0 0 1 0 0 σ W 0 κ W 2 0 0 0 0 1 , ( g μ ν ) = κ σ + κ 0 σ σ + κ W 1 0 0 1 0 0 σ σ + κ W 1 0 1 σ + κ W 2 0 0 0 0 1 ,
where  W ( x ) = e μ x , and  g = σ + κ W ( x ) .

3. Mathisson’s Spin-Curvature Force

To connect Mathisson’s classical spin-curvature force in the correspondence limit with intrinsic spin that is purely of quantum origin, it proves useful to introduce a classical model of intrinsic spin. To simplify matters, we permanently attach a free spin vector  S  to a Newtonian point particle resulting in a “pole-dipole” particle. The particle thus carries the spin vector along its world line, and the corresponding equations of motion in a gravitational field are the Mathisson–Papapetrou pole–dipole equations [21,31]:
D P μ d ς = 1 2 R μ ν α β U ν S α β ,
D S μ ν d ς = P μ U ν P ν U μ ,
where  U μ = d x μ / d ς  is the 4-velocity of the pole–dipole particle,  U μ U μ = 1 , and  ς  is its proper time. The particle’s 4-momentum is  P μ  and its spin tensor is  S μ ν , which satisfies the Frenkel–Pirani supplementary condition [32,33]
S μ ν U ν = 0 .
In this system, the inertial mass of the particle m m : = P μ U μ , and the magnitude of its spin s s 2 : = 1 2 S μ ν S μ ν , are constants of the motion. Moreover, Pirani has shown that the spin vector  S μ S μ U μ = 0 ,
S μ = 1 2 η μ ν ρ σ U ν S ρ σ , S α β = η α β γ δ U γ S δ .
is Fermi–Walker-transported along the particle’s world line [33]. That is, the Mathisson–Papapetrou equations for a spinning test particle, together with the Frenkel–Pirani supplementary condition, imply that the spin vector of a test pole–dipole particle is nonrotating in this classical model, which is consistent with the inertia of the intrinsic spin. Furthermore, the Mathisson–Papapetrou equations, together with the Frenkel–Pirani supplementary condition, imply that, in the massless limit, the spinning massless test particle follows a null geodesic with the spin vector parallel or is antiparallel to its direction of motion [34]. Hence, our classical model is consistent with physical expectations.
What is the influence of the inertia of the intrinsic spin on the motion of the spinning particle? From Equation (20), we find
P μ = m U μ + S μ ν D U ν d ς ;
thus, in the absence of spin,  P μ = m U μ , and the particle simply follows a timelike geodesic of the background gravitational field. In the presence of spin, on the other hand, the Mathisson spin-curvature force  F μ F μ U μ = 0 ,
F μ = 1 2 R μ ν α β U ν S α β = * R μ ν ρ σ U ν S ρ U σ , * R μ ν ρ σ = 1 2 η μ ν α β R α β ρ σ ,
must be taken into account [22]. It follows from Equation (23) that  P μ m U μ  is of the second order in spin; hence, the Mathisson–Papapetrou equations of motion to first order in spin become [35]
D S μ ν d ς 0 ,
and
m D U μ d ς F μ = 1 2 R μ ν α β U ν S α β .

4. Spin-Vorticity–Gravity Coupling

We now turn to the behavior of spinning test particles in the Gödel-type spacetime. Within the framework of linearized general relativity, it can be shown in general that, in source-free Ricci-flat regions of the gravitational field, the Mathisson force corresponds to the Stern–Gerlach force associated with the spin–gravitomagnetic field coupling [20]. In the Gödel-type universe, on the other hand, the spinning particle is immersed in the source of the gravitational field. Is  F μ = μ ( H s g )  still valid for the Gödel-type spacetime?
Let us consider a spinning test particle held at rest in space at fixed  ( x , y , z )  coordinates in the Gödel-type spacetime. According to the free reference observer with adapted tetrad frames of  e μ α ^  and  λ μ α ^  at the same location, the spin vector to linear order stays fixed with respect to the parallel-propagated frame as a consequence of Equation (25); that is,  S i ^ i = 1 , 2 , 3 , are constants of the motion, where  S α ^ = S μ λ μ α ^ ; hence,
S 0 ^ = 0 , S i ^ = S μ λ μ i ^ .
The motion of the comoving observer has vorticity in accordance with Equation (4), and we therefore expect that the spin should couple to the vorticity resulting in the spin-vorticity Hamiltonian given by
H s v = S · ω = Ω S 3 ^ .
Furthermore, the spin vector precesses with frequency  Ω z  with respect to the observer’s natural frame  e μ i ^  based on the spatial coordinate axes. The Hamiltonian associated with this motion is the spin–gravity Hamiltonian given by
H s g = S · B ,
where  B = Ω = Ω z  is the gravitomagnetic field in this case. The result is
H s g = Ω S 3 ^ .
The spin–gravity coupling is indeed the same as the spin–vorticity coupling in this case, since the spinning particle, while engulfed by the source of the gravitational field, is fixed in space and comoving with the observer. It is clear that in this case  μ ( H s g ) = 0 , so that the Stern–Gerlach force vanishes. To calculate the Mathisson force in this case, we need to find the Riemann curvature tensor for the Gödel-type universe, since the Mathisson force is directly proportional to the spacetime curvature.
In metric (1), the nonzero components of the Riemann tensor can be obtained from
R 0101 = Ω 2 , R 0202 = ( κ + σ ) e 2 μ x Ω 2 , R 0112 = σ e μ x Ω 2 , R 1212 = κ 4 κ σ + 5 e 2 μ x Ω 2 .
We are interested in the components of the curvature tensor projected onto the orthonormal tetrad frame  λ μ α ^  adapted to our fiducial observer, namely,
R α ^ β ^ γ ^ δ ^ = R μ ν ρ σ λ μ α ^ λ ν β ^ λ ρ γ ^ λ σ δ ^ .
The measured components of the Riemann tensor can be expressed via its symmetries as a  6 × 6  matrix in the standard manner with indices that range over the set  { 01 , 02 , 03 , 23 , 31 , 12 } . The end result is of the general form:
E H H T S ,
where  E H , and  S  represent the measured gravitoelectric, gravitomagnetic, and spatial components of the Riemann curvature tensor, respectively, and  E  and  S  are symmetric matrices, while  H  is traceless. In the case of Gödel-type spacetime (1), we find that  H = 0 , and
( E i ^ j ^ ) = Ω 2 0 0 0 Ω 2 0 0 0 0 , ( S i ^ j ^ ) = 0 0 0 0 0 0 0 0 1 + 4 κ σ Ω 2 .
These results are equally valid if the curvature tensor is projected onto the natural frame  e μ α ^  of the reference observer.
We find that the Mathisson force, given by Equation (24), can be expressed as
F μ = λ μ α ^ F α ^ , F 0 ^ = 0 , F i ^ = H j ^ i ^ S j ^ .
However,  H = 0 ; therefore, the measured components of the Mathisson force vanish as well. That is,
F μ = μ ( H s g ) = 0 .
It is important to verify this result in a quasi-inertial Fermi normal coordinate system established about the world line of an arbitrary reference observer that is spatially at rest.

5. Fermi Coordinates in Gödel-Type Spacetimes

To explore spin–gravity coupling in Fermi coordinates, it is convenient to set up a quasi-inertial system of coordinates based on the nonrotating spatial frame adapted to a fiducial geodesic observer that is at rest in space with fixed  ( x , y , z )  coordinates and a 4-velocity vector  u μ = δ 0 μ  in Gödel-type spacetime (1). The reference observer establishes in the neighborhood of its world line a Fermi normal coordinate system based on the parallel-propagated spatial frame  λ μ i ^ i = 1 , 2 , 3 , given by Equation (11). That is, at each event  x ¯ μ ( τ )  on its world line, there is a local hypersurface formed by all spacelike geodesic curves that are orthogonal to the observer’s world line at  x ¯ μ ( τ ) . Consider an event with coordinates  x μ  on this hypersurface that can be connected to  x ¯ μ ( τ )  by a unique spacelike geodesic of proper length . Then, the reference observer can assign Fermi coordinates  X μ = ( T , X i )  to  x μ  such that
T : = τ , X i : = ξ μ λ μ i ^ ( τ ) .
Here,  ξ μ ξ μ u μ = 0 , is a unit spacelike vector tangent to the unique spacelike geodesic at  x ¯ μ ( τ ) .
For the case of Gödel’s universe, one can find the exact Fermi metric coefficients [30]. The previous results are generalized here for Gödel-type spacetime (1). For the spacelike geodesics  x μ ( ) , we use Killing vector fields  t , y , and  z  to derive the equations of motion:
t + σ e μ x y = E , σ e μ x t κ e 2 μ x y = k .
z = h , t 2 2 σ e μ x t y + x 2 + κ e 2 μ x y 2 + z 2 = 1 .
Here,  E , k , and h are integration constants; moreover, a prime denotes the derivative of a spacetime coordinate with respect to proper length , e.g.,  t = d t / d . The condition  ξ μ λ μ 0 ^ = 0 , where  ξ μ = d x μ / d , implies that  E = 0 . Then, with  z = h  and  E = 0 , we find
t = σ e μ x k σ + κ ,
y = e 2 μ x k σ + κ ,
x 2 + e 2 μ x k 2 σ + κ = 1 h 2 .
The ordinary differential Equation (42) has the general solution for  x ( )  given by
e μ x = α 0 cosh a + b ,
where the constant parameters are fixed as
α 0 a = | k | μ σ + κ , a = μ 1 h 2 ,
and the condition
α 0 cosh b = 1
is imposed to satisfy  x ( 0 ) = 0 .
Substituting Equation (43) into Equations (40) and (41), we find the solutions for  t ( )  and  y ( ) :
t τ = 2 μ σ σ + κ k | k | arctan e a + b arctan e b .
y = 1 μ α 0 1 σ + κ k | k | tanh ( a + b ) tanh b .
Then, making use of Equations (8)–(11) and (37), we derive for the Fermi coordinates
T = τ , Z = h .
X cos ( Ω T ) Y sin ( Ω T ) = | k | σ + κ sinh b .
X sin ( Ω T ) + Y cos ( Ω T ) = k σ + κ .
As in [30], we introduce the cylindrical coordinates
X = ρ cos θ , Y = ρ sin θ ,
and recast Equations (49) and (50) into
cos ( θ + Ω T ) = tanh b ,
sin ( θ + Ω T ) = | k | k cosh b ,
μ ρ = a .
As a result, we rewrite the solutions (43), (47), and (46) as
e μ x = cosh ( μ ρ ) + sinh ( μ ρ ) cos ( θ + Ω T ) ,
σ + κ μ y = tanh ( μ ρ ) sin ( θ + Ω T ) 1 + tanh ( μ ρ ) cos ( θ + Ω T ) ,
tan μ 2 σ + κ σ T t = e μ ρ 1 sin ( θ + Ω T ) 1 cos ( θ + Ω T ) + [ 1 + cos ( θ + Ω T ) ] e μ ρ .
Finally, the transformation from  ( t , x , y , z )  to Fermi coordinates  ( T , X , Y , Z )  can be conveniently written in terms of the new variables:
R = μ ρ , F = θ + Ω T ,
as follows:
e μ x = cosh R + sinh R cos F ,
σ + κ μ y = sinh R sin F cosh R + sinh R cos F ,
tan μ 2 σ + κ σ T t = e R 1 sin F 1 cos F + ( 1 + cos F ) e R .
By differentiation, we obtain
d t + σ e μ x d y = d T + 1 μ σ σ + κ ( cosh R 1 ) d F ,
d x 2 + ( κ + σ ) e 2 μ x d y 2 = 1 μ 2 d R 2 + sinh 2 R d F 2 .
It remains to substitute these results into Equation (10) to derive the line element of the Gödel-type universe in terms of the Fermi coordinates. We find
d s 2 = ( 1 + L ) d T 2 2 Ω K d T ( X d Y Y d X ) + d X 2 + d Y 2 + d Z 2 + F X 2 + Y 2 ( X d Y Y d X ) 2 ,
where
L = σ 4 ( σ + κ ) sinh 2 R σ + 2 κ σ + κ ( cosh R 1 ) 2 ,
K = κ σ + κ ( cosh R 1 ) 2 R 2 ,
F = sinh 2 R R 2 1 σ σ + κ ( cosh R 1 ) 2 R 2
are functions of the variable
R = 2 Ω σ + κ σ ( X 2 + Y 2 ) 1 / 2 .

6. Spin–Gravity Coupling in Fermi Coordinates

In general, the spacetime metric in the Fermi system is given by
d s 2 = g ^ μ ν d X μ d X ν ,
where
g ^ 00 = 1 R 0 ^ i ^ 0 ^ j ^ ( T ) X i X j + , g ^ 0 i = 2 3 R 0 ^ j ^ i ^ k ^ ( T ) X j X k + ,
and
g ^ i j = δ i j 1 3 R i ^ k ^ j ^ l ^ ( T ) X k X l + .
In these expansions in powers of spatial Fermi coordinates, the coefficients are, in general, functions of T and consist of components of the Riemann curvature tensor and its covariant derivatives as measured by the reference observer that permanently occupies the spatial origin of the Fermi coordinate system. That is, the metric of the Fermi normal coordinate system established on the basis of a parallel-propagated spatial frame along the world line of a geodesic observer is the Minkowski metric plus perturbations caused by the curvature of spacetime. Fermi coordinates are admissible within a cylindrical spacetime region around the world line of the fiducial observer, and the radius of this cylinder is given by an appropriate radius of the curvature of spacetime [30].
As defined in Equation (32),  R α ^ β ^ γ ^ δ ^  are evaluated at the origin of spatial Fermi coordinates via the projection of the Riemann tensor on the tetrad frame  λ μ α ^  of the fiducial observer; indeed, for the stationary Gödel-type spacetime, the nonzero components of  R α ^ β ^ γ ^ δ ^  are constants and can be obtained from
R 0 ^ 1 ^ 0 ^ 1 ^ = R 0 ^ 2 ^ 0 ^ 2 ^ = Ω 2 , R 1 ^ 2 ^ 1 ^ 2 ^ = 1 + 4 κ σ Ω 2
via the symmetries of the Riemann curvature tensor.
We define the curvature-based gravitoelectric potential  Φ ^  and gravitomagnetic vector potential  A ^  via  g ^ 00 = 1 + 2 Φ ^  and  g ^ 0 i = 2 A ^ i  [36,37]. Indeed,
Φ ^ = 1 2 R 0 ^ i ^ 0 ^ j ^ X i X j + , A ^ i = 1 3 R 0 ^ j ^ i ^ k ^ X j X k + .
The corresponding fields are given by
E ^ = Φ ^ , B ^ = × A ^ ,
as expected; more explicitly,
E ^ i = R 0 ^ i ^ 0 ^ j ^ X j + , B ^ i = 1 2 ϵ i j k R j ^ k ^ 0 ^ l ^ X l + .
To the lowest order, the gravitomagnetic field vanishes in the Gödel-type spacetime; therefore, we need to compute higher-order terms.
In Section 5, we derived the exact Fermi metric coefficients for the Gödel-type universe. They are given explicitly by
g ^ 00 = 1 σ 4 ( σ + κ ) sinh 2 R σ + 2 κ σ + κ ( cosh R 1 ) 2 ,
g ^ 0 i = κ σ + κ Ω ( cosh R 1 ) 2 R 2 ( Y , X , 0 ) ,
and
( g ^ i j ) = 1 + A C 0 C 1 + B 0 0 0 1 ,
where
A = F Y 2 X 2 + Y 2 , B = F X 2 X 2 + Y 2 , C = F X Y X 2 + Y 2 .
The exact Fermi coordinate system has been established around the fiducial observer fixed at  X = Y = Z = 0 .
For  κ 0 , there are no closed timelike curves. In the special case of the Gödel universe with parameters (13), there are no closed timelike curves within a cylindrical region about the Z axis with
R = 2 Ω ( X 2 + Y 2 ) 1 / 2 R max , R max = 2 ln ( 1 + 2 ) .
Indeed, a circle in the  ( X , Y )  plane inside this domain is spacelike; however, it becomes null for  R = R max  and timelike for  R > R max .
The stationary and divergence-free gravitomagnetic vector field of the Gödel-type universe is given by  B ^ 1 = B ^ 2 = 0  and
B ^ 3 = κ σ + κ Ω ( cosh R 1 ) sinh R R .
It is interesting to note that  B ^ 3  and its first derivative with respect to  R  vanish at  R = 0 ; then,  B ^ 3  monotonically increases with increasing  R  and diverges as  R . More explicitly,
B ^ 3 = 2 κ σ Ω 3 ( X 2 + Y 2 ) 1 + Ω 2 σ + κ σ ( X 2 + Y 2 ) + 2 Ω 4 5 σ + κ σ 2 ( X 2 + Y 2 ) 2 + ,
so that the fiducial observer measures a null gravitomagnetic field at its location  ( X = Y = Z = 0 ) . Furthermore, the gravitomagnetic field away from the Z axis points along Z and is cylindrically symmetric; indeed, it vanishes all along Z but increases monotonically away from the Z axis and eventually diverges as the radius of the cylinder about the Z axis approaches infinity.
Within the Fermi coordinate system, it is useful to define the class of fundamental observers that remain at rest in space, each with fixed  ( X , Y , Z )  coordinates. For our present purposes, we concentrate on the set of fundamental observers that occupy a cylindrical region in the neighborhood of the Z axis. Specifically, in this region, we can express the metric tensor in Fermi coordinates as
g ^ μ ν = η μ ν + h ^ μ ν ,
where the nonzero components of the gravitational potentials are given by
h ^ 00 = Ω 2 ( X 2 + Y 2 ) , h ^ 01 = κ σ Ω 3 ( X 2 + Y 2 ) Y , h ^ 02 = κ σ Ω 3 ( X 2 + Y 2 ) X ,
and
h ^ 11 = 1 3 1 + 4 κ σ Ω 2 Y 2 , h ^ 12 = 1 3 1 + 4 κ σ Ω 2 X Y , h ^ 22 = 1 3 1 + 4 κ σ Ω 2 X 2 .
That is, for the sake of simplicity, we confine our considerations to a cylindrical region about the Z axis such that  Ω | X | = Ω | Y | ε , where  0 < ε 1  and all terms of order  ε 4  and higher are neglected in our analysis.
In the cylindrical neighborhood of the fiducial observer under consideration, fundamental observers have access to adapted orthonormal tetrad frames  φ μ α ^ , given in the Fermi coordinate system  ( T , X , Y , Z )  by
φ μ 0 ^ = ( 1 + 1 2 h ^ 00 , 0 , 0 , 0 ) , φ μ 1 ^ = ( h ^ 01 , 1 1 2 h ^ 11 , 0 , 0 ) ,
φ μ 2 ^ = ( h ^ 02 , h ^ 12 , 1 1 2 h ^ 22 , 0 ) , φ μ 3 ^ = ( 0 , 0 , 0 , 1 ) .
These tetrad axes are primarily along the Fermi coordinate directions; indeed, for  X = Y = 0 φ μ α ^ λ μ α ^ . According to these fundamental observers, a spinning particle held at rest in space has a 4-velocity vector in the Fermi system given by  U ^ μ = φ μ 0 ^ ; moreover, its spin vector has the following measured components:
S ^ 0 ^ = 0 , S ^ i ^ = S ^ μ φ μ i ^ ,
since  S ^ μ U ^ μ = 0 . Furthermore, the gravitomagnetic field at the location of the spin is given by
B ^ 1 = 0 , B ^ 2 = 0 , B ^ 3 = 1 2 ( X h ^ 02 Y h ^ 01 ) = 2 κ σ Ω 3 ( X 2 + Y 2 ) ,
in agreement with Equation (82) within our approximation scheme. The Hamiltonian for spin-gravity coupling in the Fermi frame is thus given by
H ^ s g = S ^ · B ^ = 2 κ σ Ω 3 ( X 2 + Y 2 ) S ^ 3 ^ ,
which reduces in our approximation to  2 κ σ Ω 3 ( X 2 + Y 2 ) S 3 ^ , where  S 3 ^  is a constant. The corresponding Stern–Gerlach force is then
μ H ^ s g = 4 κ σ Ω 3 S 3 ^ ( 0 , X , Y , 0 ) .
Next, we need to compute the Mathisson force in the Fermi frame, namely,
F ^ μ = 1 2 R ^ μ ν α β U ^ ν S ^ α β .
For metric (83), the curvature tensor to the first order in the perturbation is given by
R ^ μ ν α β = 1 2 ( h ^ μ β , ν α + h ^ ν α , μ β h ^ ν β , μ α h ^ μ α , ν β ) .
We are interested in the gravitomagnetic components of this curvature tensor as measured by the fundamental observers. Projection of this tensor on the tetrad frame  φ μ α ^  does not affect its components in our approximation scheme. We find in this case
( H ^ i ^ j ^ ) = 0 0 κ 1 0 0 κ 2 0 0 0 ,
where
κ 1 = 1 2 X ( X h ^ 02 Y h ^ 01 ) , κ 2 = 1 2 Y ( X h ^ 02 Y h ^ 01 ) .
Hence,  F ^ 0 ^ = 0 , and  F ^ i ^ = H ^ i ^ j ^ S ^ j ^ = ( κ 1 , κ 2 , 0 ) S 3 ^  at the level of approximation under consideration here. Moreover, Equation (89) implies
κ 1 = 4 κ σ Ω 3 X , κ 2 = 4 κ σ Ω 3 Y .
Therefore,  F ^ μ = μ H ^ s g  as measured by the fundamental observers within the cylindrical domain in the Fermi frame.
We have thus far relied on the classical pole–dipole model for the evaluation of spin–gravity coupling. It is important to demonstrate that our considerations are consistent with the solutions of the Dirac equation in the Gödel-type universe.

7. Dirac Equation in the Gödel-Type Universe

Let us start with the Dirac equation in the form [38,39]
( i γ α α m ) Ψ = 0 , μ = μ + Γ μ ,
where the fermion wave function  Ψ  is a 4-component spacetime scalar variable composed of the pair of 2-spinors  φ  and  χ :
Ψ = φ χ , φ = φ 1 φ 2 , χ = χ 1 χ 2 .
As before, we assume that the observer in the gravitational field has a natural adapted orthonormal tetrad field and
γ α = e α β ^ γ β ^ , { γ μ , γ ν } = 2 g μ ν ( x ) I 4 ,
where  I n  is the n-dimensional identity matrix, and
γ 0 ^ = I 2 0 0 I 2 , γ i ^ = 0 σ i σ i 0 .
Here,  σ i  are Pauli matrices, namely,
σ 1 = 0 1 1 0 , σ 2 = 0 i i 0 , σ 3 = 1 0 0 1 .
The spin connection  Γ μ  (also known as Fock–Ivanenko coefficients) is given by
Γ μ = i 4 e ν α ^ e ν β ^ ; μ σ α ^ β ^ , σ α ^ β ^ : = i 2 [ γ α ^ , γ β ^ ] .
Making use of tetrad frame (8), we find, after some algebra, the explicit form of the Dirac Equation (97) in Gödel-type spacetime (1):
[ γ 0 ^ σ σ + κ γ 2 ^ i t γ 1 ^ p x e μ x κ + σ γ 2 ^ p y γ 3 ^ p z + i μ 2 γ 1 ^ + μ 4 σ σ + κ γ 0 ^ Σ 3 ^ m ] Ψ = 0 .
Here, as usual, the momentum operator is  p = i  and the spin operator  Σ  is given by the matrix
Σ i ^ = σ i 0 0 σ i .
Next, due to the symmetries of Gödel-type spacetime, we assume a solution of the form
Ψ = ψ ( x ) exp ( i ω t + i k 2 y + i k 3 z ) ,
where the four components of  ψ ( x )  satisfy ordinary differential equations, namely,
d ψ d x = M ψ ,
where  M  is the  4 × 4  matrix
M = A + i k 3 0 i B + i k 3 A i B 0 0 i B + A + i k 3 i B 0 i k 3 A + i m 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 .
Here,  A ±  and  B ±  are given by
A ± = ω σ σ + κ ± Ω σ + κ σ + k 2 e μ x , B ± = ω ± Ω 2 .
The spin–vorticity–gravity coupling is evident in the way the frequency of the radiation is changed by  ± Ω / 2 , in agreement with previous results [37,40,41]. If  k 2 = 0 , the waves can only travel parallel or antiparallel to the rotation axis. In this case, matrix  M  has constant elements, and the general solution of Equation (106) can be expressed in terms of the eigenvalues and eigenfunctions of  M . It turns out that no propagation can occur in this case due to the requirement that the wave amplitude be finite at all times [41]. These general results for the Dirac equation are consistent with the propagation of the scalar and electromagnetic waves in the Gödel-type universe; for brief accounts of these latter topics, see the appendices at the end of this paper.
To deal with the general case, we henceforth assume  k 2 0  and change to  ξ = e μ x  instead of x as the independent variable. Let us recall here that  μ > 0 , since we have explicitly assumed  Ω > 0 . For  > x > , we find that  ξ  goes from zero to  + ; hence,  ξ  is a radial coordinate. In terms of  ξ , Equation (106) takes the form
ξ d ψ d ξ = M ψ ,
where matrix  M  is simply related to  M , namely,
M = A ¯ + i γ 0 i B ¯ + i γ A ¯ i B ¯ 0 0 i B ¯ + A ¯ + i γ i B ¯ 0 i γ A ¯ i m 2 Ω σ σ + κ 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 .
Here,  A ¯ ±  and  B ¯ ±  are given by
A ¯ ± = ω 2 Ω σ σ + κ ± 1 2 + β ξ , B ¯ ± = 1 2 σ σ + κ ω Ω ± 1 2 ,
and we have introduced dimensionless parameters
β = k 2 2 Ω σ σ + κ , γ = k 3 2 Ω σ σ + κ .
To clarify the structure of the resulting system (109)–(110), we note that the 4-spinor (98) can be decomposed into the sum of the left and right spinors,
ψ = ψ L + ψ R , ψ L = 1 2 ( 1 γ 5 ) ψ , ψ R = 1 2 ( 1 + γ 5 ) ψ ,
where  γ 5 : = i γ 0 ^ γ 1 ^ γ 2 ^ γ 3 ^ . By definition, the left and right spinors are eigenstates of the  γ 5  matrix:  γ 5 ψ L = ψ L  and  γ 5 ψ R = ψ R . Furthermore, we decompose the left and right spinors into the eigenstates of the  Σ 3 ^  spin matrix (i.e., the “spin-up” and “spin-down” states):
ψ L = ψ + L + ψ L , ψ R = ψ + R + ψ R , Σ 3 ^ ψ ± L = ± ψ ± L , Σ 3 ^ ψ ± R = ± ψ ± R .
After these steps, we thus have
ψ + L = L + 1 0 1 0 , ψ L = L 0 1 0 1 , ψ + R = R + 1 0 1 0 , ψ R = R 0 1 0 1 ,
where explicitly
L + = 1 2 ( φ 1 + χ 1 ) , L = 1 2 ( φ 2 + χ 2 ) ,
R + = 1 2 ( φ 1 χ 1 ) , R = 1 2 ( φ 2 χ 2 ) .
Taking these definitions into account, we can straightforwardly recast system (109)–(110) into an equivalent but more transparent form:
ξ d d ξ + A ¯ + L + = i ( B ¯ + + γ ) L + i m 2 Ω σ σ + κ R ,
ξ d d ξ A ¯ L = i ( B ¯ γ ) L + + i m 2 Ω σ σ + κ R + ,
ξ d d ξ + A ¯ + R + = i ( B ¯ + γ ) R i m 2 Ω σ σ + κ L ,
ξ d d ξ A ¯ R = i ( B ¯ + γ ) R + i m 2 Ω σ σ + κ L + .
The nontrivial mass mixes the left and right modes. However, for the massless ( m = 0 ) case or in the high-energy approximation ( m c 2 Ω 1 ), we can neglect the last terms on the right-hand sides. As a result, the left modes  L ±  decouple from the right modes  R ±  and the system reduces to
ξ d d ξ + A ¯ + L + = i ( B ¯ + + γ ) L ,
ξ d d ξ A ¯ L = i ( B ¯ γ ) L + ,
ξ d d ξ + A ¯ + R + = i ( B ¯ + γ ) R ,
ξ d d ξ A ¯ R = i ( B ¯ + γ ) R + .
It is interesting to mention that in this approximation scheme Equation (109) can also be solved by a different approach that is briefly described in Appendix A.

7.1. Explicit Solutions

Multiplying Equation (122) by  i ( B ¯ γ )  and Equation (123) by  i ( B ¯ + + γ ) , we derive the second-order equations for the left modes:
ξ d d ξ + A ¯ + ξ d d ξ A ¯ L = B ¯ + B ¯ + γ ( B ¯ + B ¯ ) + γ 2 L ,
ξ d d ξ A ¯ + ξ d d ξ + A ¯ L + = B ¯ + B ¯ + γ ( B ¯ + B ¯ ) + γ 2 L + .
In Equation (111), it is useful to introduce a dimensionless parameter  α ,
α : = ω 2 Ω σ σ + κ , A ¯ ± = α + β ξ ± 1 2 ;
then,
A ¯ + A ¯ = ( α + β ξ ) 2 1 4 , A ¯ + A ¯ = 1 ,
B ¯ + B ¯ = σ σ + κ ω 2 ( 2 Ω ) 2 1 16 , B ¯ + B ¯ = 1 2 σ σ + κ .
Employing the ansatz
L ± = ξ 1 u 1 2 ,
we can recast Equations (126) and (127) into the form
ξ 2 d 2 d ξ 2 u s + 1 4 μ ˜ f 2 β 2 ξ 2 2 β ξ ( α + s ) u s = 0 ,
where  s = ± 1 2  and
μ ˜ f 2 = α 2 + γ 2 B ¯ + B ¯ + γ ( B ¯ + B ¯ ) = 1 μ 2 ω 2 κ σ + κ + k 3 Ω / 2 2 .
With a new independent variable  ξ ˜ = 2 | β | ξ , Equation (132) can be reduced to Whittaker’s equation [42]:
d 2 u s d ξ ˜ 2 + 1 4 + κ ˜ f ξ ˜ + 1 4 μ ˜ f 2 ξ ˜ 2 u s = 0 ,
where
κ ˜ f = β | β | α + s .
The Dirac field is a linear perturbation on the Gödel-type spacetime; therefore,  ψ ( x )  should be bounded. Demanding that  ψ ( x )  be finite everywhere, the acceptable solution of Whittaker’s equation is given via the confluent hypergeometric functions by
u s = u s 0 exp ( 1 2 ξ ˜ ) ξ ˜ 1 2 + μ ˜ f 1 F 1 ( 1 2 + μ ˜ f κ ˜ f , 1 + 2 μ ˜ f ; ξ ˜ ) ,
where
1 2 + μ ˜ f κ ˜ f = n , n = 0 , 1 , 2 , .
In this case, the confluent hypergeometric function can be expressed in terms of the associated Laguerre polynomial.
For  k 2 < 0 β  is negative and  κ ˜ f = α + s = ω 2 Ω σ σ + κ + s , with  s = ± 1 2 . Then, combining Equations (137) and (133), we derive the dispersion relation
ω = ( 2 n + 1 + 2 s ) Ω ± ( k 3 Ω / 2 ) 2 κ σ ( 2 n + 1 + 2 s ) 2 Ω 2 1 / 2 .
Note that solutions with both signs of energy are admissible.
Similarly, multiplying Equation (124) by  i ( B ¯ + γ )  and Equation (125) by  i ( B ¯ + γ ) , we derive the second-order equations for the right modes:
ξ d d ξ + A ¯ + ξ d d ξ A ¯ R = B ¯ + B ¯ γ ( B ¯ + B ¯ ) + γ 2 R ,
ξ d d ξ A ¯ + ξ d d ξ + A ¯ R + = B ¯ + B ¯ γ ( B ¯ + B ¯ ) + γ 2 R + .
Using the ansatz
R ± = ξ 1 v 1 2 ,
we recast Equations (139) and (140) into
ξ 2 d 2 d ξ 2 v s + 1 4 μ ¯ f 2 β 2 ξ 2 2 β ξ ( α + s ) v s = 0 ,
where  s = ± 1 2 , as before, but now we have
μ ¯ f 2 = α 2 + γ 2 B ¯ + B ¯ γ ( B ¯ + B ¯ ) = 1 μ 2 ω 2 κ σ + κ + k 3 + Ω / 2 2 .
With the independent variable  ξ ˜ = 2 | β | ξ , Equation (142) can again be reduced to Whittaker’s equation:
d 2 v s d ξ ˜ 2 + 1 4 + κ ˜ f ξ ˜ + 1 4 μ ¯ f 2 ξ ˜ 2 v s = 0 ,
where  κ ˜ f  is given by Equation (135). The regular solution of Equation (144) is given by
v s = v s 0 exp ( 1 2 ξ ˜ ) ξ ˜ 1 2 + μ ¯ f 1 F 1 ( 1 2 + μ ¯ f κ ˜ f , 1 + 2 μ ¯ f ; ξ ˜ ) ,
where
1 2 + μ ¯ f κ ˜ f = n , n = 0 , 1 , 2 , .
As before, we can combine Equations (146) and (143) to derive the dispersion relation
ω = ( 2 n + 1 + 2 s ) Ω ± ( k 3 + Ω / 2 ) 2 κ σ ( 2 n + 1 + 2 s ) 2 Ω 2 1 / 2 .
The motion of Dirac waves in the Gödel-type universe is in general agreement with the corresponding results for the scalar and electromagnetic wave propagation described in Appendix B and Appendix C.

7.2. Dealing with Subtle Points of Dirac Theory on Curved Spacetimes

In order to have a correct quantum-mechanical interpretation, Dirac Equation (97) should be recast into the form of the Schrödinger equation:
i Ψ t = H Ψ .
In flat spacetime with the Minkowski metric  η μ ν = diag ( 1 , 1 , 1 , 1 ) , the trivial frame  e μ α ^ = δ α μ , and the spin connection  Γ μ = 0 , this is straightforward. Multiplying Equation (97) by  γ 0 ^ , we derive Schrödinger Equation (148) with the Hermitian Hamiltonian
H = β D m + α D · p .
Here we denote, as usual, the matrices
β D : = γ 0 ^ , α D i ^ : = γ 0 ^ γ i ^ = 0 σ i σ i 0 , i = 1 , 2 , 3 .
In addition, one also needs a quantum-probabilistic picture, which is related to the normalization of the wave function. As is well known, a direct consequence of the Dirac Equation (97) is the conservation of the vector current, which, in flat spacetime, can be expressed as
μ J μ = 0 , J μ = Ψ ¯ γ μ Ψ .
Integration over 3-space yields a global conservation law
d 3 x J 0 = d 3 x Ψ Ψ = constant = 1 .
The physical interpretation of the Dirac fermion dynamics is based on Equations (148) and (152), especially when the fermionic particle interacts with external fields.
Dirac theory on curved manifolds, however, involves a number of subtleties. In particular, the differential conservation law (151) is replaced by its curved version:
μ J μ = 1 g μ g J μ = 0 , J μ = e μ α ^ Ψ ¯ γ α ^ Ψ ,
which yields the global conservation law
d 3 x g J 0 = d 3 x g e 0 α ^ Ψ γ 0 ^ γ α ^ Ψ = constant = 1 .
For the natural Gödel-type tetrad frame (8), we have  e 0 α ^ γ 0 ^ γ α ^ = 1 σ σ + κ α D 2 ^ ; therefore, the physical interpretation of the solutions is unclear. In addition, Dirac Equation (103) obviously cannot be directly recast into the form of the Schrödinger wave Equation (148).
Both issues are related to the choice of the tetrad frame, which is defined up to an arbitrary local Lorentz transformation. The choice (8) corresponds to the so-called Landau-Lifshitz gauge with  e 0 i ^ = 0  and  e i 0 ^ = 0 . The situation is essentially improved when one chooses the Schwinger gauge for the frame, where  e i 0 ^ = 0  and  e 0 i ^ = 0 . Then, Equation (154) reduces to an “almost flat” form:
d 3 x g e 0 0 ^ Ψ Ψ = 1 ,
and the Dirac equation is straightforwardly recast into the Schrödinger form [38,39].
This suggests replacing the original tetrad frame (8) with a new one
e ˜ 0 ^ = κ σ + κ t + σ κ e μ x y , e ˜ 1 ^ = x , e ˜ 2 ^ = e μ x κ y , e ˜ 3 ^ = z ,
where we assume  κ > 0 . Obviously, this choice corresponds to the Schwinger gauge  e ˜ i 0 ^ = 0  and  e ˜ 0 i ^ = 0  for  i = 1 , 2 , 3 .
For Gödel-type spacetimes, the two frames (8) and (156) are related by the Lorentz transformation
e ˜ α ^ = Λ β ^ α ^ e β ^ ,
where explicitly
Λ α ^ β ^ = σ + κ κ 0 σ κ 0 0 1 0 0 σ κ 0 σ + κ κ 0 0 0 0 1 .
Interestingly, the transformation with constant matrix elements is global, whereas, in general, only local Lorentz transformations are possible.
The change of a frame on the spacetime affects the fermionic wave function
Ψ Ψ ˜ = L 1 Ψ
via the spinor matrix L that satisfies
L 1 γ α ^ L = Λ α ^ β ^ γ β ^ .
Using a convenient parametrization with  cosh ζ = σ + κ κ  and  sinh ζ = σ κ , we easily derive
L = cosh ( ζ / 2 ) I 4 + sinh ( ζ / 2 ) α D 2 ^ = cosh ( ζ / 2 ) I 2 sinh ( ζ / 2 ) σ 2 sinh ( ζ / 2 ) σ 2 cosh ( ζ / 2 ) I 2 .
The spinor transformation (159) mixes the spin-up and spin-down states ( L ± ) for the left modes (and similarly for the right modes), and an appropriate normalization of the solutions should be fixed for the squares  Ψ ˜ Ψ ˜  of the transformed wave functions.

8. Dirac Equation in Fermi Frame

Let us next consider the Dirac equation in the quasi-inertial Fermi frame of Section 6. We are interested in the propagation of Dirac particles as described by fundamental observers that are all spatially at rest in the Fermi frame and occupy the limited cylindrical region about the Z axis such that  Ω | X | = Ω | Y | ε . As before, we ignore all terms of order  ε 4  and higher. The preferred observers have adapted orthonormal tetrad frames  φ μ α ^  given in Equations (86)–(87). Let us note that  φ μ α ^  can be written in the  ( T , X , Y , Z )  coordinate system as
φ μ 0 ^ = ( 1 + 1 2 h ^ 00 , h ^ 01 , h ^ 02 , 0 ) , φ μ 1 ^ = ( 0 , 1 + 1 2 h ^ 11 , h ^ 12 , 0 ) ,
φ μ 2 ^ = ( 0 , 0 , 1 + 1 2 h ^ 22 , 0 ) , φ μ 3 ^ = ( 0 , 0 , 0 , 1 ) .
We employ perturbations beyond Minkowski spacetime in our Fermi frame; hence, in the absence of  h ^ μ ν , we have  φ μ α ^ δ α μ . To simplify matters even further, we assume henceforth that the deviation from Minkowski spacetime is only due to the gravitomagnetic potentials  h ^ 01 = κ σ Ω 3 Y ( X 2 + Y 2 )  and  h ^ 02 = κ σ Ω 3 X ( X 2 + Y 2 )  that give rise to the gravitomagnetic field  B ^ = ( 0 , 0 , B ^ 3 ) , where  B ^ 3 = 2 κ σ Ω 3 ( X 2 + Y 2 ) .
With these assumptions, the spin connection in (102) can be computed using the tetrad system  φ μ α ^  that is adapted to our reference observers and we find
γ μ Γ ^ μ = i 2 B ^ 3 σ 3 0 0 σ 3 .
That is, the spin connection is proportional to the gravitomagnetic field of the Gödel-type universe in the Fermi frame under consideration here.
For the sake of simplicity, we assume a solution of the Dirac equation that propagates along the Z axis and is of the form
Ψ ^ = ψ ^ ( X , Y ) exp ( i ω T + i k 3 Z ) .
Moreover, it is convenient to define
X ^ = ψ ^ 1 ψ ^ 3 , Y ^ = ψ ^ 2 ψ ^ 4 .
In this case, Dirac’s equation reduces to
X + i Y + κ σ ω Ω 3 ( X 2 + Y 2 ) ( X + i Y ) X ^ = i κ σ Ω 3 ( X 2 + Y 2 ) σ 1 Y ^ + i k 3 ω + m ω m k 3 Y ^ ,
and
X i Y κ σ ω Ω 3 ( X 2 + Y 2 ) ( X i Y ) Y ^ = i κ σ Ω 3 ( X 2 + Y 2 ) σ 1 X ^ + i k 3 ω + m ω m k 3 X ^ .
Here,  X : = / X , etc.; furthermore, we note that
( X ± i Y ) ( X 2 + Y 2 ) 2 = 4 ( X 2 + Y 2 ) ( X ± i Y ) ,
( X ± i Y ) [ ( X 2 + Y 2 ) ( X i Y ) ] = 4 ( X 2 + Y 2 ) .
In the absence of the gravitational perturbation, the positive-frequency plane wave solutions of the free Dirac equation propagating in the Z direction are given by
w ^ ± e i ω T + i k 3 Z ,
where the spin of the Dirac particle is either parallel ( w ^ + ) or antiparallel ( w ^ ) to the Z direction; that is,
w ^ + = N 1 0 ϱ 0 , w ^ = N 0 1 0 ϱ .
Here,  N  and  N  are positive normalization constants,  ω = ( m 2 + k 3 2 ) 1 / 2 , and
ϱ : = k 3 ω + m = ω m k 3 .
With these background states, we solve Equations (167) and (168) to the linear order in the gravitomagnetic perturbation and obtain, after some algebra,
Ψ ^ + = N exp [ 3 κ 8 σ ω Ω 3 ( X 2 + Y 2 ) 2 ] i κ 4 σ Ω 3 ϱ ( X 2 + Y 2 ) ( X + i Y ) ϱ exp [ 3 κ 8 σ ω Ω 3 ( X 2 + Y 2 ) 2 ] i κ 4 σ Ω 3 ( X 2 + Y 2 ) ( X + i Y ) e i ω T + i k 3 Z
Ψ ^ = N i κ 4 σ Ω 3 ϱ ( X 2 + Y 2 ) ( X i Y ) exp [ 3 κ 8 σ ω Ω 3 ( X 2 + Y 2 ) 2 ] i κ 4 σ Ω 3 ( X 2 + Y 2 ) ( X i Y ) ϱ exp [ 3 κ 8 σ ω Ω 3 ( X 2 + Y 2 ) 2 ] e i ω T + i k 3 Z .
These solutions of Dirac’s equation exhibit the coupling of spin with the gravitomagnetic field of a Gödel-type universe and may be compared and contrasted with the results of Appendix C for the propagation of circularly polarized electromagnetic waves along the Z axis in the Fermi frame.
We should note that fermions in Gödel-type universes have been the subject of a number of previous studies; see, for instance, refs. [43,44,45,46,47,48] and the references cited therein.

9. Discussion

Spin–gravity coupling represents a physically important subject matter in view of the basic nature of the intrinsic spin of particles and the universality of the gravitational interaction. We have investigated in detail the coupling of intrinsic spin with the gravitomagnetic fields of a three-parameter class of Gödel-type spacetimes. These stationary and homogeneous rotating universes are characterized by the set of constant parameters  ( κ , σ , μ ) ; for  κ < 0 ; there are closed timelike curves (CTCs) in spacetime, while, for  κ 0 , the CTCs are absent. For  ( κ , σ , μ ) ( 1 , 2 , 2 Ω ) , we recover Gödel’s rotating universe model, where  Ω > 0  is the frequency of rotation. Regarding the background Gödel-type spacetimes, we have studied Dirac’s equation and worked out its solutions; furthermore, we have extended our results to exact Fermi normal coordinate systems in these universes. We have shown that the Stern–Gerlach force due to the coupling of intrinsic spin with the gravitomagnetic field of a Gödel-type spacetime is in agreement with the correspondence limit with the classical Mathisson spin-curvature force. This is a nonlinear generalization of previous work that focused on linearized general relativity [20]. Our main results turn out to be independent of the possible causality difficulties of the Gödel-type spacetimes.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this theoretical study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Alternative Solution of Equation (109)

The purpose of this appendix is to present a different approach to the solution of Equation (109).
We can write Equation (109) in the form
ξ d ( U ψ ) d ξ = U M U 1 ( U ψ ) ,
where  U  is a constant unitary matrix given by
U = 1 2 I 2 I 2 I 2 I 2 .
Under this similarity transformation, we have
U γ 0 ^ U = 0 I 2 I 2 0 = γ 5 , U γ i ^ U = γ i ^ .
That is, the standard representation of the Dirac matrices is thus transformed to the chiral (Weyl) representation. Employing this representation, we find
U M U 1 = A ¯ + i B ¯ + i γ 0 0 i B ¯ + i γ A ¯ 0 0 0 0 A ¯ + i B ¯ + i γ 0 0 i B ¯ + i γ A ¯ i m 2 Ω σ σ + κ 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 ,
where  A ¯ ±  and  B ¯ ±  are given by Equation (111). By expressing  U ψ  in the form
U ψ = 2 R L , R = R + R , L = L + L ,
where  R  and  L  are now right-handed and left-handed two-component Weyl spinors, we recover the system of Equations (118)–(121). The rest of the analysis would follow the treatment presented in Section 7.

Appendix B. Scalar Waves in the Gödel-Type Universe

Consider first a scalar field  ϕ  of inertial mass m propagating on the background Gödel-type spacetime (1). The wave equation is
g μ ν ϕ ; μ ν m 2 c 2 2 ϕ = 0 ,
where  / ( m c )  is the Compton wavelength of the particle. The back reaction is of the second order in the perturbation and can be neglected. The scalar wave equation can be written as
1 g x μ g g μ ν ϕ x ν m 2 c 2 2 ϕ = 0 ,
where, for metric (1),  g = e μ x σ + κ . Moreover,  t y , and  z  are Killing vector fields; therefore, we assume that
ϕ ( x ) = e i ω t + i k 2 y + i k 3 z ϕ ¯ ( ξ ) , ξ : = e μ x ,
where  ξ  increases from 0 to as the x coordinate decreases from  +  to  . In terms of the new radial variable  ξ , the equation for  ϕ ¯  reduces to
d 2 ϕ ¯ d ξ 2 α s 2 + β s ξ + ζ s ( ζ s + 1 ) ξ 2 ϕ ¯ = 0 ,
where
α s = k 2 μ σ + κ , β s = 2 ω k 2 c μ 2 σ σ + κ , ζ s ( ζ s + 1 ) = 1 μ 2 ω 2 c 2 κ σ + κ + k 3 2 + m 2 c 2 2 .
Let us assume  ζ s > 0  and note that for  k 2 = 0 , Equation (A9) for  ϕ ¯  has solutions of the form  ξ ζ s  and  ξ ζ s + 1  that diverge at  ξ = 0  and  ξ = , respectively. However, the scalar perturbation must be finite everywhere; therefore, waves cannot freely propagate parallel or antiparallel to the axis of rotation of the Gödel-type spacetime. Next, for  k 2 0 , we introduce a new variable,  ξ ¯ : = | k 2 | σ Ω ( σ + κ ) ξ , in terms of which Equation (A9) takes the form of Whittaker’s equation [42]:
d 2 ϕ ¯ d ξ ¯ 2 + 1 4 + κ ¯ s ξ ¯ + 1 4 μ ¯ s 2 ξ ¯ 2 ϕ ¯ = 0 ,
where
κ ¯ s = ω 2 Ω k 2 | k 2 | σ σ + κ , μ ¯ s = ± ( ζ s + 1 2 ) .
In terms of the confluent hypergeometric functions, bounded solutions of this equation can be expressed up to proportionality constants by the following:
exp ( 1 2 ξ ¯ ) ξ ¯ ζ s + 1 1 F 1 ( n , 2 ζ s + 2 ; ξ ¯ ) , n = 0 , 1 , 2 , .
Here,  ζ s > 0 μ ¯ s = ζ s + 1 / 2 , and
ζ s + 1 + ω 2 Ω k 2 | k 2 | σ σ + κ = n , ω = ± 2 Ω ( n + ζ s + 1 ) σ + κ σ ,
for  k 2 < 0  (upper plus sign) or  k 2 > 0  (lower minus sign), respectively. Negative frequency in the case of  k 2 > 0  indicates that waves traveling forward in time move backward along the y direction. Finally, we note that only certain frequencies are allowed for the scalar waves; for instance, for   k 2 < 0 , we have  ω n = 2 Ω ( n + ζ s + 1 ) ( σ + κ ) / σ . That is,
ω n ± = ( 2 n + 1 ) Ω ± 4 n ( n + 1 ) κ σ + 1 Ω 2 + k 3 2 + m 2 c 2 2 1 / 2 ,
where  ω n + > 0  for all n by definition, while  ω n > 0  for  n = 1 , 2 , 3 ,  only if
ω n + ω n = 4 n ( n + 1 ) σ + κ σ Ω 2 k 3 2 m 2 c 2 2 > 0 .
For further work on the scalar perturbations of the Gödel-type universe and its extensions, see [49,50,51].

Appendix C. Electromagnetic Waves in the Gödel-Type Universe

The propagation of electromagnetic radiation in the Gödel universe was originally investigated in the search for the coupling of photon helicity with the rotation of matter [52]. In Gödel-type spacetimes, Maxwell’s equations can be reduced to an equation of the form of Equation (A9), where, instead of the quantities in Equation (A10), we find [53]
α s α e m = K 2 ± μ σ + κ , β s β e m = 2 ω K 2 ± c μ 2 σ σ + κ ,
and  ζ s ζ e m , where
ζ e m ( ζ e m + 1 ) = 1 μ 2 ω 2 c 2 κ σ + κ + ( K 3 ± ) 2 2 Ω K 3 ± ,
since photon is massless ( m = 0 ). The helicity coupling evident in Equation (A18) is consistent with the spin–vorticity–gravity coupling described in Section 4. That is, based on the results of Section 4, we would expect the corresponding Hamiltonian for a photon to be proportional to  ± K 3 ± Ω / ω  so that, in terms of frequency, we would have  ± K 3 ± Ω / ω . The effect should disappear in the case of a null geodesic consistent with the eikonal limit  ω . For further extensions and generalizations to Gödel-type universes, see [53,54,55,56,57,58,59].

AEM Waves in the Fermi Frame

We consider the propagation of electromagnetic radiation on the background quasi-inertial Fermi normal coordinate system. In terms of the Faraday tensor  F μ ν , the source-free Maxwell equations can be expressed as
F [ μ ν , ρ ] = 0 , ( g F μ ν ) , ν = 0 .
Using the same approach as in [52], we replace the gravitational field by a hypothetical optical medium that occupies Euclidean space with Cartesian Fermi coordinates  ( X , Y , Z ) . The electromagnetic field Equation (A19) reduces to the traditional form of Maxwell’s equations in a medium with the following decompositions:
F μ ν ( E ˜ , B ˜ ) , g F μ ν ( D ˜ , H ˜ ) .
That is,  F 0 i = E ˜ i  and  F i j = ϵ i j k B ˜ k ; similarly,  g F 0 i = D ˜ i  and  g F i j = ϵ i j k H ˜ k . Here,  ϵ i j k  is the totally antisymmetric symbol with  ϵ 123 = 1 . The corresponding optical medium turns out to be gyrotropic with constitutive relations [60,61,62,63,64]:
D ˜ i = ϵ ^ i j E ˜ j ( G ^ × H ˜ ) i , B ˜ i = μ ^ i j H ˜ j + ( G ^ × E ˜ ) i ,
where the characteristics of the medium are conformally invariant and are given by
ϵ ^ i j = μ ^ i j = g ^ g ^ i j g ^ 00 , G ^ i = g ^ 0 i g ^ 00 .
Expressing electromagnetic fields in the standard complex form and introducing the Riemann–Silberstein vectors,
F ˜ ± = E ˜ ± i H ˜ , S ˜ ± = D ˜ ± i B ˜ ,
the wave propagation equation can be expressed as the Dirac equation for photons in the gravitational field. That is,
× F ˜ ± = ± i S ˜ ± t , · S ˜ ± = 0 ,
where
S ˜ p ± = ϵ ^ p q F ˜ q ± ± i ( G ^ × F ˜ ± ) p .
The Dirac-type equation implies  t ( · S ˜ ± ) = 0 ; therefore, if  · S ˜ ± = 0  initially, then it is valid for all time.
To interpret the physical meaning of these results, it proves useful to consider plane electromagnetic waves of frequency  ω  propagating along the z axis in a global inertial frame with coordinates  x μ = ( t , x )  in Minkowski spacetime. Maxwell’s equations are linear; therefore, we can use complex electric and magnetic fields and use the convention that only the real parts correspond to measurable quantities. The waves can have two independent orthogonal linear polarization states along the  x ^  and  y ^  directions, where  x ^  is a unit vector along the x axis, etc. The circular polarization states are constructed from the linear polarization states via superposition; in this case, the electric ( e ) and magnetic ( b ) fields can be expressed as
e ± = 1 2 a ± ( x ^ ± i y ^ ) e i ω ( t z ) , b ± = i 2 a ± ( x ^ ± i y ^ ) e i ω ( t z ) ,
where  a +  and  a  are constant complex amplitudes. Here, the upper (lower) sign represents waves in which the orthogonal electric and magnetic fields rotate in the positive (negative) sense about the direction of the wave motion. In the case of a photon with positive (negative) circular polarization, the photon has positive (negative) helicity, namely, its spin is  +  ( ) along its direction of propagation. The Riemann–Silberstein vectors have interesting behaviors for helicity states of the photon; in fact, for positive-helicity radiation,
e + + i b + = a + ( x ^ + i y ^ ) e i ω ( t z ) , e + i b + = 0 ,
while for radiation with negative helicity,
e + i b = 0 , e i b = a ( x ^ i y ^ ) e i ω ( t z ) .
Hence,  e + i b  ( e i b ) represents in essence an electromagnetic wave with positive (negative) helicity. It is important to note that Equations (A24) and (A25) that represent the propagation of the electromagnetic test fields in a gravitational field completely decouple for different helicity states.
Imagine the propagation of electromagnetic waves with definite helicity along the Z axis in the Fermi normal coordinate system in the Gödel-type spacetime. The universe rotates in the negative sense about the Z axis. We confine our considerations to the cylindrical region near the rotation axis where the perturbation analysis contained in Equations (83)–(85) is valid. To simplify matters, we take into account only the gravitomagnetic potentials  h ^ 01  and  h ^ 02 , and we ignore the other potentials; therefore, in Equation (A22), we have
ϵ ^ i j = μ ^ i j 1 , G ^ κ σ Ω 3 ( X 2 + Y 2 ) ( Y , X , 0 ) .
It is straightforward to show that, in this case, the field Equations (A24) and (A25) have the solution
F ˜ 1 ± = a ^ ± exp [ i ω ( T Z ) κ 4 σ ω Ω 3 ( X 2 + Y 2 ) 2 ] ,
F ˜ 2 ± = ± i a ^ ± exp [ i ω ( T Z ) κ 4 σ ω Ω 3 ( X 2 + Y 2 ) 2 ] ,
and  F ˜ 3 ± = 0 . Here,  a ^ +  and  a ^  are constant amplitudes for the positive and negative helicity waves in the Fermi frame, respectively. If the wave propagates along the axis of rotation (i.e., the  Z  direction), then, in Equations (A30) and (A31), we have  Z Z  and  ±  in the exponents of these equations, as well as in the coefficient of the latter equation. For  Ω = 0 , the Fermi frame reduces to a global inertial frame in Minkowski spacetime and we recover waves of the form given in Equations (A26)–(A28).
The helicity–gravitomagnetic field coupling is evident in these results and corresponds to Equations (89) and (90) of Section 6; indeed, the form of this coupling is reminiscent of the helicity–twist coupling studied in [65].

References

  1. Wigner, E.P. On unitary representations of the inhomogeneous Lorentz group. Ann. Math. 1939, 40, 149–204. [Google Scholar] [CrossRef]
  2. Mashhoon, B. Nonlocal Gravity; Oxford University Press: Oxford, UK, 2017. [Google Scholar]
  3. Demirel, B.; Sponar, S.; Hasegawa, Y. Measurement of the spin-rotation coupling in neutron polarimetry. New J. Phys. 2015, 17, 023065. [Google Scholar] [CrossRef] [Green Version]
  4. Danner, A.; Demirel, B.; Sponar, S.; Hasegawa, Y. Development and performance of a miniaturised spin rotator suitable for neutron interferometer experiments. J. Phys. Commun. 2019, 3, 035001. [Google Scholar] [CrossRef]
  5. Danner, A.; Demirel, B.; Kersten, W.; Wagner, R.; Lemmel, H.; Sponar, S.; Hasegawa, Y. Spin-rotation coupling observed in neutron interferometry. NPJ Quantum Inf. 2020, 6, 23. [Google Scholar] [CrossRef] [Green Version]
  6. Yu, T.; Luo, Z.; Bauer, G.E.W. Chirality as generalized spin-orbit interaction in spintronics. Phys. Rept. 2023, 1009, 1–115. [Google Scholar] [CrossRef]
  7. de Oliveira, C.G.; Tiomno, J. Representations of Dirac equation in general relativity. Nuovo Cimento 1962, 24, 672–687. [Google Scholar] [CrossRef]
  8. Hehl, F.W.; Ni, W.-T. Inertial effects of a Dirac particle. Phys. Rev. D 1990, 42, 2045–2048. [Google Scholar] [CrossRef]
  9. Damião Soares, I.I.; Tiomno, J. The physics of the Sagnac-Mashhoon effects. Phys. Rev. D 1996, 54, 2808–2813. [Google Scholar] [CrossRef]
  10. Larmor, J. On the theory of the magnetic influence on spectra; and on the radiation from moving ions. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1897, 44, 503–512. [Google Scholar] [CrossRef] [Green Version]
  11. Mashhoon, B. On the gravitational analogue of Larmor’s theorem. Phys. Lett. A 1993, 173, 347–354. [Google Scholar] [CrossRef]
  12. Everitt, C.W.F.; DeBra, D.B.; Parkinson, B.W.; Turneaure, J.P.; Conklin, J.W.; Heifetz, M.I.; Keiser, G.M.; Silbergleit, A.S.; Holmes, T.; Kolodziejczak, J.; et al. Gravity Probe B: Final results of a space experiment to test general relativity. Phys. Rev. Lett. 2011, 106, 221101. [Google Scholar] [CrossRef] [Green Version]
  13. Everitt, C.W.F.; Muhlfelder, B.; DeBra, D.B.; Parkinson, B.W.; Turneaure, J.P.; Silbergleit, A.S.; Acworth, E.B.; Adams, M.; Adler, R.; Bencze, W.J.; et al. The Gravity Probe B test of general relativity. Class. Quantum Gravity 2015, 32, 224001. [Google Scholar] [CrossRef] [Green Version]
  14. Papini, G. Spin-gravity coupling and gravity-induced quantum phases. Gen. Relativ. Gravit. 2008, 40, 1117–1144. [Google Scholar] [CrossRef] [Green Version]
  15. Mashhoon, B. On the coupling of intrinsic spin with the rotation of the Earth. Phys. Lett. A 1995, 198, 9–13. [Google Scholar] [CrossRef]
  16. Mashhoon, B. Gravitational couplings of intrinsic spin. Class. Quantum Gravity 2008, 17, 2399–2409. [Google Scholar] [CrossRef] [Green Version]
  17. Tarallo, M.G.; Mazzoni, T.; Poli, N.; Sutyrin, D.V.; Zhang, X.; Tino, G.M. Test of Einstein Equivalence Principle for 0-spin and half-integer-spin atoms: Search for spin-gravity coupling effects. Phys. Rev. Lett. 2014, 113, 023005. [Google Scholar] [CrossRef] [Green Version]
  18. Fadeev, P.; Wang, T.; Band, Y.B.; Budker, D.; Graham, P.W.; Sushkov, A.O.; Kimball, D.F.J. Gravity Probe Spin: Prospects for measuring general-relativistic precession of intrinsic spin using a ferromagnetic gyroscope. Phys. Rev. D 2021, 103, 044056. [Google Scholar] [CrossRef]
  19. Vergeles, S.N.; Nikolaev, N.N.; Obukhov, Y.N.; Silenko, A.J.; Teryaev, O.V. General relativity effects in precision spin experimental tests of fundamental symmetries. Physics-Uspekhi 2023, 66, 109–147. [Google Scholar] [CrossRef]
  20. Mashhoon, B. Gravitomagnetic Stern-Gerlach force. Entropy 2021, 23, 445. [Google Scholar] [CrossRef]
  21. Mathisson, M. Neue Mechanik materieller Systeme. Acta Phys. Pol. 1937, 6, 163–200, Reprinted: Gen. Relativ. Gravit. 2010, 42, 1011–1048. [Google Scholar] [CrossRef]
  22. Mashhoon, B. Spin-gravity coupling. Acta Phys. Pol. 2008, Suppl. 1, 113–122. [Google Scholar] [CrossRef]
  23. Gödel, K. An example of a new type of cosmological solutions of Einstein’s field equations of gravitation. Rev. Mod. Phys. 1949, 21, 447–450. [Google Scholar] [CrossRef] [Green Version]
  24. Obukhov, Y.N. On physical foundations and observational effects of cosmic rotation. In Colloquium on Cosmic Rotation, Berlin, 1998; Scherfner, M., Chrobok, T., Shefaat, M., Eds.; Wissenschaft und Technik Verlag: Berlin, Germany, 2000; pp. 23–96. [Google Scholar] [CrossRef]
  25. Rebouças, M.J.; Tiomno, J. Homogeneity of Riemannian spacetimes of Gödel type. Phys. Rev. D 1983, 28, 1251–1264. [Google Scholar] [CrossRef]
  26. Rebouças, M.J.; Tiomno, J. A class of inhomogeneous Gödel-type models. Nuovo Cim. B 1985, 90, 204–210. [Google Scholar] [CrossRef]
  27. Hawking, S.W.; Ellis, G.F.R. The Large Scale Structure of Space-Time; Cambridge University Press: Cambridge, UK, 1973. [Google Scholar] [CrossRef]
  28. Stephani, H.; Kramer, D.; MacCallum, M.A.H.; Hoenselaers, C.; Herlt, E. Exact Solutions of Einstein’s Field Equations, 2nd ed.; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar] [CrossRef]
  29. Griffiths, J.B.; Podolsky, J. Exact Space-Times in Einstein’s General Relativity; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar] [CrossRef]
  30. Chicone, C.; Mashhoon, B. Explicit Fermi coordinates and tidal dynamics in de Sitter and Gödel spacetimes. Phys. Rev. D 2006, 74, 064019. [Google Scholar] [CrossRef] [Green Version]
  31. Papapetrou, A. Spinning test-particles in general relativity. I. Proc. R. Soc. A 1951, 209, 248–258. [Google Scholar] [CrossRef]
  32. Frenkel, J. Die Electrodynamik des rotierenden Elektrons. Z. Phys. 1926, 37, 243–262. [Google Scholar] [CrossRef]
  33. Pirani, F.A.E. On the physical significance of the Riemann tensor. Acta Phys. Pol. 1956, 15, 389–405, Reprinted: Gen. Relativ. Gravit. 2009, 41, 1215–1232. [Google Scholar] [CrossRef]
  34. Mashhoon, B. Massless spinning test particles in a gravitational field. Ann. Phys. 1975, 89, 254–257. [Google Scholar] [CrossRef]
  35. Chicone, C.; Mashhoon, B.; Punsly, B. Relativistic motion of spinning particles in a gravitational field. Phys. Lett. A 2005, 343, 1–7. [Google Scholar] [CrossRef] [Green Version]
  36. Mashhoon, B. Gravitoelectromagnetism: A brief review. In The Measurement of Gravitomagnetism: A Challenging Enterprise; Iorio, L., Ed.; Nova Science: New York, NY, USA, 2007; pp. 29–39. [Google Scholar] [CrossRef]
  37. Bini, D.; Mashhoon, B.; Obukhov, Y.N. Gravitomagnetic helicity. Phys. Rev. D 2022, 105, 064028. [Google Scholar] [CrossRef]
  38. Obukhov, Y.N.; Silenko, A.J.; Teryaev, O.V. Spin in an arbitrary gravitational field. Phys. Rev. D 2013, 88, 084014. [Google Scholar] [CrossRef] [Green Version]
  39. Obukhov, Y.N.; Silenko, A.J.; Teryaev, O.V. General treatment of quantum and classical spinning particles in external fields. Phys. Rev. D 2017, 96, 105005. [Google Scholar] [CrossRef] [Green Version]
  40. Mashhoon, B.; Obukhov, Y.N. Spin precession in inertial and gravitational fields. Phys. Rev. D 2013, 88, 064037. [Google Scholar] [CrossRef] [Green Version]
  41. Bini, D.; Chicone, C.; Mashhoon, B. Spacetime splitting, admissible coordinates and causality. Phys. Rev. D 2012, 85, 104020. [Google Scholar] [CrossRef] [Green Version]
  42. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions; National Bureau of Standards: Washington, DC, USA, 1964. [Google Scholar]
  43. Damião Soares, I. Gravitational coupling of neutrinos to matter vorticity: Microscopic asymmetries. Phys. Rev. D 1981, 23, 272–286. [Google Scholar] [CrossRef]
  44. Leahy, D. Scalar and neutrino fields in the Gödel universe. Int. J. Theor. Phys. 1982, 21, 703–753. [Google Scholar] [CrossRef]
  45. Damião Soares, I.; Rodrigues, L.M.C.S. Gravitational coupling of neutrinos to matter vorticity. II. Microscopic asymmetries in angular momentum modes. Phys. Rev. D 1985, 31, 422–424. [Google Scholar] [CrossRef] [PubMed]
  46. Pimentel, L.O.; Macías, A. Klein-Gordon and Weyl equations in the Gödel universe. Phys. Lett. A 1986, 117, 325–327. [Google Scholar] [CrossRef]
  47. Villalba, V.M. Dirac spinor in a nonstationary Gödel-type cosmological universe. Mod. Phys. Lett. A 1993, 8, 3011–3018. [Google Scholar] [CrossRef] [Green Version]
  48. Pimentel, L.O.; Camacho, A.; Macías, A. Weyl equation in Gödel type universes. Mod. Phys. Lett. A 1994, 9, 3703–3706. [Google Scholar] [CrossRef] [Green Version]
  49. Hiscock, W.A. Scalar perturbations in the Gödel universe. Phys. Rev. D 1978, 17, 1497–1500. [Google Scholar] [CrossRef]
  50. Thakurta, S.N.G. Scalar perturbations in some cosmological metrics. Phys. Rev. D 1980, 21, 864–866. [Google Scholar] [CrossRef]
  51. Saibatalov, R.X. Scalar field in causal Gödel-type space-times. Gravit. Cosmol. 2001, 7, 293–296. [Google Scholar]
  52. Mashhoon, B. Influence of gravitation on the propagation of electromagnetic radiation. Phys. Rev. D 1975, 11, 2679–2684. [Google Scholar] [CrossRef]
  53. Korotky, V.A.; Obukhov, Y.N. Electromagnetic waves in rotating universe. Mosc. Univ. Phys. Bull. 1991, 46, 4–6. [Google Scholar]
  54. Cohen, J.M.; Vishveshwara, C.V.; Dhurandhar, S.V. Electromagnetic fields in the Gödel universe. J. Phys. A Math. Gen. 1980, 13, 933–938. [Google Scholar] [CrossRef]
  55. Obukhov, Y.N.; Korotky, V.A. The Weyssenhoff fluid in Einstein-Cartan theory. Class. Quantum Gravity 1987, 4, 1633–1657. [Google Scholar] [CrossRef]
  56. Korotky, V.A.; Obukhov, Y.N. Kinematic analysis of cosmological models with rotation. Sov. Phys. JETP 1991, 72, 11–15. [Google Scholar]
  57. Abd-Eltwab, Y.A. Electromagnetic fields in the homogeneous Gödel-type universes. Nuovo Cim. B 1993, 108, 465–469. [Google Scholar] [CrossRef]
  58. Saibatalov, R.X. Electromagnetic field in causal and acausal Gödel-type space-times. Gen. Relativ. Gravit. 1995, 27, 697–711. [Google Scholar] [CrossRef]
  59. Havare, A.; Yetkin, T. Exact solution of the photon equation in stationary Gödel-type and Gödel spacetimes. Class. Quantum Gravity 2002, 19, 2783–2791. [Google Scholar] [CrossRef] [Green Version]
  60. Skrotskii, G.V. The influence of gravitation on the propagation of light. Sov. Phys. Dokl. 1957, 2, 226–229. [Google Scholar]
  61. Plebanski, J. Electromagnetic waves in gravitational fields. Phys. Rev. 1960, 118, 1396–1408. [Google Scholar] [CrossRef]
  62. de Felice, F. On the gravitational field acting as an optical medium. Gen. Relativ. Gravit. 1971, 2, 347–357. [Google Scholar] [CrossRef]
  63. Volkov, A.M.; Izmest’ev, A.A.; Skrotskii, G.V. The propagation of electromagnetic waves in a Riemannian space. Sov. Phys. JETP 1971, 32, 686–689. [Google Scholar]
  64. Hehl, F.W.; Obukhov, Y.N. Foundations of Classical Electrodynamics: Charge, Flux, and Metric; Birkhäuser: Boston, MA, USA, 2003. [Google Scholar] [CrossRef] [Green Version]
  65. Bini, D.; Chicone, C.; Mashhoon, B.; Rosquist, K. Spinning particles in twisted gravitational wave spacetimes. Phys. Rev. D 2018, 98, 024043. [Google Scholar] [CrossRef] [Green Version]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mashhoon, B.; Molaei, M.; Obukhov, Y.N. Spin-Gravity Coupling in a Rotating Universe. Symmetry 2023, 15, 1518. https://doi.org/10.3390/sym15081518

AMA Style

Mashhoon B, Molaei M, Obukhov YN. Spin-Gravity Coupling in a Rotating Universe. Symmetry. 2023; 15(8):1518. https://doi.org/10.3390/sym15081518

Chicago/Turabian Style

Mashhoon, Bahram, Masoud Molaei, and Yuri N. Obukhov. 2023. "Spin-Gravity Coupling in a Rotating Universe" Symmetry 15, no. 8: 1518. https://doi.org/10.3390/sym15081518

APA Style

Mashhoon, B., Molaei, M., & Obukhov, Y. N. (2023). Spin-Gravity Coupling in a Rotating Universe. Symmetry, 15(8), 1518. https://doi.org/10.3390/sym15081518

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop