1. Introduction
Spin effects in gravity straddle the boundary between quantum and classical physics. The difference between quantum and classical behavior becomes particularly transparent with spin, which is relevant to low energy approaches to quantum gravity.
Quantum beats are an unequivocal indication that the system considered obeys the laws of quantum mechanics. Quantum systems with a two-dimensional Hilbert space are also called qubits. This definition is borrowed from quantum computing where two-level quantum systems play a predominant role. Qubits provide examples of systems that are genuinely quantum mechanical and, at the same time, simple, can be studied within the confines of first quantization, and are ideal in the study of relativistic gravity close to, or at the quantum level. There are gravitational qubits in the universe. Some of them are disussed below.
In what follows, use is made of the external field approximation [
1,
2] that treats gravity as a classical theory when it interacts with quantum particles. This approximation can be applied successfully to all those problems involving gravitational sources of weak to intermediate strength for which the full-fledged use of general relativity is not required [
3,
4,
5,
6,
7,
8,
9,
10,
11], it is encountered in the solution of relativistic wave equations, and takes different forms according to the statistics obeyed by the particles [
2,
6,
12,
13,
14]. The approximation can also be applied to theories in which acceleration has an upper limit [
15,
16,
17,
18,
19,
20,
21,
22,
23,
24] and that allow for the resolution of astrophysical and cosmological singularities in quantum gravity [
25,
26]. It is of interest to those theories of asymptotically safe gravity that can be expressed as Einstein gravity coupled to a scalar field [
27], and can produce results complementary to those of the method of space-time deformation [
28].
At the same time, theoretical developments by Mashhoon [
19,
20,
29,
30,
31], and by other authors [
12,
32,
33,
34,
35,
36] in the field of spin-gravity coupling require a scheme that involves all components of the metric tensor. Finally, recent experimental observations of important rotation-related classical effects [
37,
38,
39], of spin-rotation coupling for photons [
40] and neutrons [
41], the development of a spin rotator for neutron interferometry [
42], and the generation of spin currents via spin-rotation coupling [
43] indicate the degree of maturity and breadth of scope reached by the field.
In the formalism introduced in references [
1,
2], the effect of gravity on wave functions is contained in a phase factor. If phase differences develop in processes involving the qubits studied below, measurements become possible. This is the common thread that links the various sections of this work. If gravitation produces qubits, then these may be observable and yield useful experimental results. Quantum physics and gravity may meet well before the onset of the quantum gravity regime usually associated with Planck’s length and there still are interesting problems to investigate at lower scales. For instance, in addition to the important classical effect observed and discussed in references [
37,
38,
39], there also is a quantum Lense–Thirring effect, that represents the action of the Lense–Thirring metric on a particle wave function. By applying the procedure of [
1,
2], one finds [
36] that the phase difference produced by a gravitational source of mass
M, radius
R and angular velocity
is
, where
is the effective Lense–Thirring frequency of a gyroscope and
replaces the period of a satellite in the classical calculation. Its observation with neutron interferometers of typical dimension
cm, still seems difficult but would complete nicely what we know at present about rotation in relativity.
For the sake of completeness, some essential points are being repeated. The key player in what follows is the covariant Dirac equation
that determines the behavior of spin-1/2 particles in the presence of a gravitational field
. In Equation (
1),
,
is the covariant derivative,
the spin connection and the matrices
satisfy the relations
. Both
and
can be obtained from the usual constant Dirac matrices
by using the vierbein fields
and the relations
where
. We use units
and the notations are as in [
3].
Equation (
1) can be solved exactly [
12,
44] to first order in the metric deviation
, where the Minkowski metric
has signature -2. This is achieved by first transforming Equation (
1) into the equation
where
By multiplying Equation (
3) on the left by
, we obtain the equation
whose solution
is exact to first order. The operator
is defined as
and
satisfies the usual free Dirac equation
In Equations (
4) and (
7), the path integrals are taken along the classical world line of the particle starting from an arbitrary reference point
P. Only the path to
needs to be known in the integrations indicated because Equation (
4) already is a first order solution. The positive energy solutions of Equation (
8) are given by
where
,
,
,
and
represents the Pauli matrices. In addition,
can take the forms
and
where
, and
.
and
are the angular and linear momentum operators of the particle. It follows from Equations (
6) and (
4) that the solution of Equation (
1) can be written in the form [
6]
and also as
where
is of first order in
. The factor
on the r.h.s. of Equation (
11) appears because both sides of the equation must agree when the gravitational field vanishes.
On multiplying Equation (
1) on the left by
and using the relations
and
we obtain the equation
In Equations (
13) and (
14),
is the linearized Riemann tensor,
R the corresponding Ricci scalar, and
.
By using Equation (
4), we also find
Equation (
14) implies that the gyro-gravitational ratio of a massive Dirac particle is one, as found in [
45,
46,
47].
The transformations of coordinates
, with
small of first order, lead to the “gauge” transformations
. It is therefore necessary to show that
in Equation (
11) is gauge invariant. In fact, on applying Stokes theorem to a closed spacetime path
C and using Equation (
12), we find that
changes by
where
is a surface bound by
C and
is the total angular momentum of the particle. Equation (
16) shows that Equations (
10) and (
11) are gauge invariant and confirms that, to first order in the gravitational field, the gyro-gravitational ratio of a Dirac particle is one. Use of Equation (
10) or Equation (
11) assures the correct treatment of both spin and angular momentum.
The plan of this work is as follows. In
Section 2 we discuss qubits represented by particles in accelerators. In
Section 3 we derive the gravitational deflection of particles propagating in a gravitational background represented by the Lense–Thirring metric and obtain the contribution due to the rotation of the source. The neutrino helicity transitions are derived in
Section 4 and some astrophysical consequences are discussed in
Section 5. Spin currents and spin motion are presented in
Section 6 and are followed by a summary.
2. Spin-Rotation Coupling in Accelerators
The spin-rotation effect described by Mashhoon is conceptually important, since it extends our knowledge of rotational inertia to the quantum level and violates the principle of equivalence [
29,
31] that is well-tested at the classical level.
It has, of course, been argued that the principle of equivalence does not hold true in the quantum world. This is the case for phase shifts in particle interferometers [
44,
48] and wave functions depend on the masses of the particles involved [
49]. In addition, the equivalence principle does not apply in the context of the causal interpretation of quantum mechanics as shown by Holland [
50]. Several models predicting quantum violations of the equivalence principle have also been discussed in the literature [
51], also in connection with neutrino oscillations [
52,
53,
54,
55]. The Mashhoon term, in particular, yields different potentials for different particles and for different spin states and cannot, therefore, be regarded as universal. It plays, nonetheless, an essential role in precise measurements of the
factor of the muon.
The experiment [
56,
57] involves muons in a storage ring. Muons on equilibrium orbits within a small fraction of the maximum momentum are almost completely polarized with spin vectors pointing in the direction of motion. As the muons decay, those electrons projected forward in the muon rest frame are detected around the ring. Their modulated angular distribution reflects the precession of the muon spin along the cyclotron orbits.
Our calculations use the covariant Dirac equation and are performed in the rotating frame of the muon and do not therefore require a relativistic treatment of inertial spin effects [
58]. Then the vierbein formalism yields
and
where
and
are the three-acceleration and three-rotation of the observer and
in the chiral representation of the usual Dirac matrices. The second term in Equation (
17) represents the Mashhoon effect. The first term drops out. The remaining contributions to the Dirac Hamiltonian, to first order in
and
, are [
32,
44]
All quantities in
H are time-independent and are referred to a left-handed set of three axes rotating about the
-axis in the clockwise direction of motion of the muons. The muon momentum is directed along the
-axis which is tangent to the muon orbits. The magnetic field is
. Only the Mashhoon term then couples the helicity states of the muon. The remaining terms contribute to the overall energy
E of the states, and we indicate by
the corresponding part of the Hamiltonian.
Before decay the muon states can be represented as
where
and
are the right and left helicity states of the Hamiltonian
and satisfy the equation
The total effective Hamiltonian is
, where
represents the total magnetic moment of the muon and
is the Bohr magneton. We will neglect the presence of electric fields that also affect the muon spin. Their effects can be controlled in suitable ways [
57].
The coefficients
and
in Equation (
19) evolve in time according to
where
M is the matrix
and
represents the width of the muon and is not particularly relevant to what follows. Equations (
21) and (
22) describe a two-dimensional qubit. The non-diagonal form of
M (when
) implies that rotation does not couple universally to matter.
M has eigenvalues
and eigenstates
The muon states that satisfy Equation (
21), and the condition
at
, are
where
The spin-flip probability is therefore
The
-term in Equation (
24) accounts for the observed exponential decrease in electron counts due to the loss of muons by radioactive decay [
57]. The term in square brackets represents the well known phenomenon of quantum beats which one should expect because muons are quantum systems. It also represents the characteristic behavior of a two-dimensional qubit.
The spin-rotation contribution to
is represented by
which is the cyclotron angular velocity
[
57]. The spin-flip angular frequency is then
which is precisely the observed modulation frequency of the electron counts [
57,
59] and yields the value
of the energy level splitting. This result is independent of the value of the anomalous magnetic moment of the particle. It is therefore the spin-rotation coupling that gives evidence to the
term in
by exactly cancelling, in
, the much larger contribution
that one would get if the fermion had no anomalous magnetic moment. The cancellation is made possible by the non-diagonal form of
M and is therefore a direct consequence of the violation of the equivalence principle.
It is perhaps surprising that spin-rotation coupling as such has almost gone unnoticed for such a long time. It is, however, significant that its effect is observed in an experiment that has already provided crucial tests of quantum electrodynamics and a test of Einstein’s time-dilation formula to better than a 0.1 percent accuracy.
Applications of these ideas to compound spin systems like heavy ions in accelerators can be found in [
60,
61,
62].
3. Geometrical Optics of Spin-1/2 Particles
In this Section we study the propagation of a spin-1/2 particle in the Lense–Thirring metric [
63] represented, in its post-Newtonian form, by
where
and
M,
R,
and
are mass, radius, angular velocity, and angular momentum of the source. The vierbein field to
is
The gravitational contribution in Equation (
28) can be further isolated by writing
. The components of the spin connection can be calculated using Equations (
2) and (
28) and are
and have the explicit form
In what follows, use is made of the Dirac representation of the
, of the first derivative of
with respect to
and of the second derivative
where
are the Christoffell symbols of the second type.
For the Lense–Thirring metric and to
, these are
In the geometrical optics approximation , where k is the momentum of the particle, the geometrical phase is sufficient to reproduce the classical angle of deflection, as it should, but also some effects due to the angular velocity of rotation of the source.
The deflection angle
is defined by
where
is the unperturbed momentum and
, for
.
It follows from Equations (
7) and (
10) that, once
is chosen to be a plane wave solution of the flat spacetime Dirac equation, the geometrical phase of a particle of four-momentum
is given by
where
and
The components of
can be determined from the equation
We consider the two cases of propagation along the z-axis, which is parallel to the angular momentum of the source, and along the x-axis, orthogonal to it. In both instances, the fermions are assumed to be ultrarelativistic, i.e., .
When motion is along the z-direction and (i.e., are eigenstates of ).
We consider fermions starting from
with impact parameter
and propagating along
,
. We find
and
We finally obtain
which is the deflection predicted by general relativity for photons, with corrections due to the fermion mass and to
. In the limit
Equation (
40) reduces to
When the fermions propagate along
x, the deflection angle is
The first term is just that predicted by general relativity.
Contrary to the case of propagation along
z, the contribution of
does not vanish in the limit
. In fact, in this limit we get
4. Neutrino Helicity Transitions
In what follows, it is convenient to write the left and right neutrino wave functions in the form
where
represents the Pauli matrices.
are eigenvectors of
corresponding to negative and positive helicity and
,
. This notation already takes into account the fact that if
are the helicity states, then we have
for relativistic neutrinos. The propagation is
in vacuo.
In general, the spin precesses during the motion of the neutrino. This can be expected because of the presence of in .
We now study the helicity flip of one flavor neutrinos as they propagate in the gravitational field produced by a rotating mass. The neutrino state vector can be written as
where
and
is an affine parameter along the world-line. In order to determine
and
, we can write Equation (
10) as
where
and
is a plane wave solution of Equation (
9). The latter can be written as
should also be normalized. However, this is unnecessary, because it is shown below that
is already of
and can only produce higher order terms. From Equations (
45), (
46), and (
48) we obtain
An equation for
can be derived in an entirely similar way.
If we consider neutrinos which are created in the left-handed state, then
, and we obtain
where
. As remarked in [
64],
need not be a null vector if we assume that the neutrino moves along an “average” trajectory. We also find, to lowest order,
where
are the usual Christoffel symbols, and
In order to solve the evolution equations for
and
and complete the equations describing this two-dimensional qubit, one also needs the terms
and
of the usual qubit matrix
M. In what follows, we compute the probability amplitude Equation (
52) for neutrinos propagating along the
z and the
x directions explicitly.
For propagation along the
z-axis, we have
and
and we choose
. We find
Summing up, and neglecting terms of
, Equation (
52) becomes
As a consequence
and the probability amplitude for the
transition is of
, as expected.
Integrating Equation (
55) from
to
z, yields
It also follows that
In this qubit, the first term in Equation (
57) comes from the mass of the gravitational source. The second from the source’s angular momentum and vanishes for
because the contribution from
to 0 exactly cancels that from 0 to
. In fact, if we consider neutrinos propagating from 0 to
, we obtain
According to semiclassical spin precession equations [
65], there should be no spin motion because spin and
are parallel. The probabilities Equations (
57) and (
58) mark therefore a departure from expected results. They yield however results that are small of second order. Both expressions vanish for
, as it should because helicity is conserved [
66]. It is interesting to observe that spin precession also occurs when
vanishes [
67,
68]. In the case of Equation (
57) the mass contribution is larger when
, which, close to the source, with
, becomes
and is always satisfied. In the case of Equation (
58), the rotational contribution is larger if
which restricts the region of dominance to a strip about the
z-axis in the equatorial plane, if the source is compact and
is relatively large.
In proximity of the source where the gravitational field is stronger and
the evolution equations for
and
are
and
,
is almost constant, and
and
oscillate with frequency
The contribution of the source rotation is therefore
. The helicity oscillations discussed are, in principle, relevant in astrophysics because right-handed neutrinos are considered sterile. This point is discussed in the next section.
When propagation is along
x, we put
,
. The calculation can be simplified by assuming that the motion is in the equatorial plane with
,
. We then have
Summing up, and neglecting terms of
, Equation (
52) becomes
The contributions to
again vanish and we get
Integrating Equation (
62) from
to
x, we obtain
and
Obviously, the contribution of
M is the same as for
z-axis propagation. However, the two cases differ substantially in the behavior of the term containing
. In this case, in fact, the term does not vanish for
. If we consider neutrinos generated at
and propagating to
, we find
The
M term is larger when
. At the poles
and the
M term dominates because the condition
is always satisfied. The angular momentum contribution prevails in proximity of the equatorial plane. The transition probability vanishes at
.
An altogether different type of qubit is represented by neutrino flavor oscillations. They have been discussed in the context of the Lense–Thirring metric in [
6]. The qubit frequency is in this case proportional to
, where
and
are the masses of the neutrino mass eigenstates.
5. Neutrino Conversion in Supernovae
The results of the previous section may be applied to the propagation of a beam of neutrinos in vacuo. The presence of a medium is realized by means of a potential
V. The neutrinos are massive and may therefore have a magnetic moment
. In the presence of an external magnetic field
and of the Mashhoon term proportional to the angular velocity of the source
, the evolution equations become [
33,
69]
where
and
are the
R and
L projection operators. Equation (
66) leads to neutrino oscillations.
The frequency of oscillation is then
where
is the component of the angular velocity normal to the neutrino trajectory. In particular, if a beam of neutrinos consists of
particles at
, the relative numbers of
and
at
z will be
These oscillations are interesting because the
’s, if they exist, do not interact. They would therefore provide an energy dissipation mechanism with possible astrophysical implications. The conversion rate is not large for galaxies and white dwarfs. In fact one can obtain from Equation (
67)
for galaxies of size
L for which
km/s. Similarly, for white dwarfs for which
, one finds
. On the other hand, the
’s diffuse out of a canonical neutron star in a time 1 to 10 s, during which they travel a maximum distance
cm between collisions. This and the fact that for a millisecond pulsar the conversion rate
is
at distances
cm suggest that the dynamics of the star could be affected by such a cooling mechanism. Indeed the star may even cool too rapidly at higher rotational speeds for a pulsar to form.
The magnetic moment of the neutrino does not appear in the calculations because magnetic spin-flip rates of magnitude comparable to Equation (
67) would require magnetic moments in excess of the value
predicted by the standard model.
The behavior of neutrinos in a medium is modified by a potential
V that vanishes for
’s. In the core of a supernova
V can be written as [
70]
where
, the
Y’s represent the lepton fractions present, and
g/cm
3. For supernovae
V can be large, of the order of several electron volts, and rotation may be neglected. Only the acceleration term in Equation (
67) need be considered and the effective Hamiltonian then has the form [
69]
where
and
is the component of the acceleration transverse to the neutrino trajectory. If the initial state is pure
and the number of particles in this state is
at
, then the corresponding numbers of
and
at
z are
where
If
, spin precession is strongly suppressed and the flux of particles at
z consists mainly of
’s. The conversion takes place at resonance if
.
Summarizing, the components of acceleration transverse to the particle path couple to its spin. This and the Mashhoon term applied to massive neutrinos produce oscillations, which may have macroscopic effects if the ’s are sterile, as frequently assumed. In fact, conversion by rotation-spin coupling may help to explain why pulsars of period shorter than a millisecond are relatively rare.
6. Spin Currents
The realization that the flow of spin angular momentum can be separated from that of charge has recently stimulated intense interest in fundamental spin physics [
71], particularly in view of its applications [
43,
72,
73].
In this section we study the generation and control of spin currents by rotation and acceleration [
3]. In this context the fundamental tool still is the covariant Dirac equation.
We use the first order solutions of Equation (
1) that have the form
where
is a solution of Equation (
8) and the operator
is given by Equation (
10) or Equation (
11).
When acceleration and rotation are present,
is given by [
32,
44]
where
and
represent acceleration and rotation respectively. To first order the tetrad is given by
from which the spinorial connection can be calculated in the usual way. The result is
and
.
For electrons,
corresponds to the choice
and
to
. Substituting into Equation (
9), one finds the spinors
and
. These are not eigenspinors of the matrix
and do not, therefore, represent the spin components in the
z-direction. They become however eigenspinors of
when
, or when
(electron rest frame).
The appropriate way to determine whether there is transfer of angular momentum between the external non-inertial field and the electron spin is to use the third rank spin current tensor [
74]
that in Minkowski space satisfies the conservation law
when all
vanish and yields in addition the expected result
in the rest frame of the particle. Writing
, using the relation
and substituting Equations (
72) and (
11) into Equation (
75), one obtains, to
,
It is therefore possible to separate
in inertial and non-inertial parts. The first term on the r.h.s. of Equation (
76) gives the result expected when
and the external field vanishes. From Equation (
76) one finds
where terms containing
and
have been eliminated. It therefore follows that the external field invalidates the simple conservation law
and that there is in this case continual interchange between spin and orbital angular momentum. The result is entirely similar to that found for external electromagnetic fields [
74]. Thus, in principle, one can use non-inertial fields to generate spin currents. In the rest frame of the particle and when
, one finds
and
. In Equation (
78)
corresponds to
. The direct coupling of the non-inertial field to the particle’s spin current violates the law
. Conservation is restored if the parameters
, or
, or both vanish.
Qubits appear in the actual spin motion. In general, transfer of angular momentum between external and non-inertial fields occurs when the operator
has some non-diagonal matrix elements. If in fact at time
a beam of electrons is entirely of the
variety, at time
t the fraction of
is
. The last expression becomes, along the electron world line,
where, as usual,
and
is the affine parameter along the world line. From [
6]
one can see that
A useful way to visualise the spin motion under the action of rotation and acceleration follows from the Mashhoon term
, where
, and from
. The two interaction terms lead to the first order equation of motion [
75]
Note that
, introduced by writing
, where
, does not contribute to Equation (
80) because
. Note also that the terms
drop out, in the particle rest frame, on account of
. No mixed effects of first order in rotation or acceleration and first order in the electromagnetic field are therefore present in this calculation. This applies to all terms containing the magnetic field
, like the Zeeman term, and electric fields, like the spin-orbit interaction, that are present in the lowest order Dirac Hamiltonian that can be derived from Equation (
1) [
44]. To
, contributions to Equation (
80) from the electromagnetic field are present in the actual determination of the electron’s path, as stated above.
From Equation (
80) one obtains
where
. The parameter
corresponds to the electron four-momentum when
and
.
Some general conclusions can be drawn from Equation (
83).
If , the particles move in the -plane. If, in addition, , then only if .
If, however, is also due to rotation, the conditions imply and therefore . This is the relevant case of motion in the -plane with rotation along an axis perpendicular to it. One cannot, therefore, have a rotation-induced spin current in this instance.
Even for
one can have
if one of
and
does not vanish. This is a direct consequence of the spin rotation interaction or Mashhoon term contained in Equation (
1) [
19,
20,
44].
A few examples are discussed below.
Consider an electron wave packet moving along the
x-axis of a frame rotating about the same axis. Then
and the remaining parameters are
and
. While the electron propagates along
x,
and
propagate in opposite directions along
x because of Equations (
80) and (
82). For a beam the spin current generated by rotation is therefore
, with obvious meaning of the symbols. One finds
which holds for
on account of the requirement
.
Consider next a wave packet moving in the plane
itself moving along
z with velocity
, while rotating about the
z-axis with
. The other parameters are
and
. From Equation (
83) one finds
and
, where
R is the radius of the circle described by the wave packet in the plane
. The motion of the center of mass of the wave packet is helical and so is the motion of the spin components which propagate, however, in opposite directions giving rise to a spin current.
In the case of motion occurring in the plane
(
) and
, one also gets
and
where
,
, and
is the cyclotron frequency of the electrons along the circular path determined by the constant magnetic field
B. Substituting Equations (
84) in (
79) one finds
which holds for all
t for which
. While both spin up and down electrons move on a circle of radius
R about the
z-axis, they propagate in opposite directions because of the spin-rotation coupling, thus generating a spin current.
Summarizing, rotation and acceleration can be used to generate and control spin currents. This follows from the covariant Dirac equation and its exact solutions to . To this order, external electromagnetic fields can be taken into account through the particle motion. The transition amplitude for the conversion spin-up to spin-down is proportional to and is expressed as a function of , , and of the electron four-momentum before the onset of rotation and acceleration. The same expression suggests criteria for the generation of spin currents and the transfer of momentum and angular momentum to them. No energy can, of course, be transferred from the non-inertial fields to the spin currents as long as remains stationary.
The particular forms of Equation (
83) discussed above provide additional examples of gravitational qubits.
7. Summary
Gravitational qubits are the simplest quantum systems that can be employed to study gravity. They exist in the laboratory and astrophysical conditions. Attention has been focussed on spin- fermions because of their simple eigenstate structure. Spin-flip transitions occur in nature frequently. We have considered here some of those that are characterized by sustained oscillations.
To the laboratory belong particles rotating in accelerators. Their quantum beats refer to spin oscillations mediated by the Mashhoon spin-rotation interaction. They have been observed and play an important role in important measurement of the anomalous magnetic moment of the muon.
Section 3 is entirely devoted to the calculation of the deflection of fermions in a gravitational background described by the Lense–Thirring metric. The interest is limited here to the action of rotation on the particle spin. The procedure can be also applied to bosons. Basically, the part of the deflection that does not depend on rotation is that predicted by general relativity. Rotation of the source yields in general additional corrections which are due to the particle spin and are therefore quantum mechanical. These, for instance, are present in the deflection Equations (
40) and (
42), but with a noticeable difference: in the case of Equation (
42) the contribution of the source angular momentum does not vanish in the limit
.
The neutrino helicity oscillations in vacuo described in
Section 4 are small, but intriguing because
evolve into
which are sterile. This energy dissipation mode could be relevant to compact astrophysical objects. The introduction in
Section 5 of a medium in which neutrinos propagate increases the
transition probability if the medium potential is attuned to the difference of the mass squared of the neutrino mass eigenstates. The fact that pulsars of period shorter than a millisecond are rare, lends support to the dissipation mechanism implied by the
conversion by rotation-spin coupling.
Spin currents and spin motion are introduced in
Section 6. It is shown that gravitational fields violate the usual formulation of the spin conservation law. Some particular qubits are discussed together with the associated spin currents. We also consider the case of fermions accelerated to the maximal acceleration limits contemplated by Caianiello [
15,
16,
17]. The calculations [
76] confirm that continual interchange between spin and angular momentum can occur in this instance, but only if the acceleration is time-dependent. This requires a transfer of energy from a very compact star, or a black hole and the particle. Even in the case, uniform acceleration produces no observable effects on the particle spin, in agreement with [
77].
Near-neighbor momenta oscillations of a gas of particles in a gravitational field are discussed in [
78].