Non--geodesic circular motion of massive spinning test bodies

Recent interest on studying possible violations of the Equivalence Principle has led to the development of space satellite missions testing it for bodies moving on circular orbits around Earth. This experiment establishes that the validity of the Equivalence Principle is independent of the composition of bodies. However, the internal degrees of freedom of the bodies (such as spin) were not taken into account. In this work, it is shown exactly that the circular orbit motion of test bodies does present a departure from geodesic motion when spin effects are not negligible. Using a Lagrangian theory for spinning massive bodies, an exact solution for their circular motion is found showing that the non--geodesic behavior manifests through different tangential velocities of the test bodies, depending on the orientation of its spin with respect to the total angular momentum of the satellite. Besides, for circular orbits, spinning test bodies present no tangential acceleration. We estimate the difference of the two possible tangential velocities for the case of circular motion of spinning test bodies orbiting Earth.

them as limiting cases of rotating fluids moving in gravitational fields. The Mathisson-Papapetrou equations (MPE) have been used to obtain several exact solutions (see for example Refs. [14][15][16][17][18]). However, the MPE present several problems in the description of spinning massive particles. For example, they are dynamical equations of third order that do not preserve the square root of the Casimir operator of the Poincaré group P µ P µ (formed by the momentum P µ ), among others. These difficulties are analyzed below [19].
Action approaches can be used to describe the dynamics of spinning massive particles from first principles [20]. In this work, we use the Lagrangian theory developed in Res. [21][22][23][24][25][26][27][28][29][30][31][32][33], which also allow us to avoid the several issues of MPE. This theory can be derived rigously from first principles allowing the proper treatment of the crucial lack of parallelism between velocity and momentum, which otherwise cannot be obtained as the canonical momentum cannot be appropriately defined (as in MPE). This Lagrangian theory has been used to study spinning massive particles (tops) in different contexts and gravitational fields [25][26][27][28][29][30][31][32][33], always finding new effects on the dynamics associated to the non-geodesics motion of tops due to spin-gravity coupling. Moreover, in Ref. [21] was shown that this Lagrangian model for tops matches the experimental conclusions of Ref. [2]. The Lagrangian theory exactly predicts the results of 87 Rb atoms experiment [2], rigourously showing that in tops in a vertical free-falling trajectory with spins aligned (with the trajectory), the forces induced by the spin-gravity coupling vanishes, and thereby the top does follow a geodesic [21]. Besides, in Ref. [21], a different and more concrete experimental setting was proposed for tops moving "parabolically" in a non-geodesic orbit, where a measurement could be possibly performed.
This previous success in the agreement of the results which stem from Lagrangian theory for tops with experiment, leads to wonder what experimental settings can be appropriated to measure deviations of from geodesic orbits. It is the purpose of this work to study the circular orbits of tops in a Schwarzschild background using this Lagrangian theory, showing how the spin induces non-geodesic motion of the test bodies. We apply these results to estimate these deviations for a possible circular orbit around Earth. There exist other schemes that allow us to study the motion of tops in gravitational fields, such as post-Newtonian approximation for spinning massive particle [34][35][36][37][38][39]. In particular, the chaotic motions in post-Newtonian systems of spinning compact binaries were investigated in Refs. [40][41][42][43][44][45][46][47][48][49]. Also, circular orbits of tops on a Schwarzschild background have been studied using the MPE approach [50], and the chaotic motions of tops in this spacetime background were explored in Refs. [51,52]. However, in here we restrict ourselves to the Lagrangian formalism of Sec. II, as it allows us to obtain an exact solution for the motion of tops in circular trajectories, without the difficulties introduced by the MPE approach.

II. LAGRANGIAN THEORY FOR TOPS
The Lagrangian model for spinning particles consider tops with mass m, spin J, energy E and total angular momentum j. The full theory is developed in Refs. [21][22][23][24][25][26][27], and we limit ourselves here to highlight its most relevant results.

A. Equations of motion
It is well-known that the velocity u µ of a spinning particle is not parallel, in general, to the canonical momentum vector P µ . The velocity vector may, under some circumstances, become spacelike [22][23][24]. However, the momentum vector remains always timelike and gives rise to the dynamical conservation law of mass m 2 ≡ P µ P µ > 0 [23,27]. The spin of tops is defined through an antisymmetric tensor S µν (see below). The action S = L dλ associated to the Lagrangian theory for tops is λ-reparametrization invariant, where the Lagrangian L(a 1 , a 2 , a 3 , a 4 ) = (a 1 ) 1/2 L a 2 /a 1 , a 3 /(a 1 ) 2 , a 4 /(a 1 ) 2 is an arbitrary function of four invariants a 1 , a 2 , a 3 , a 4 , and L is an arbitrary function of three variables where a 1 ≡ u µ u µ , a 2 ≡ σ µν σ µν = −tr(σ 2 ), a 3 ≡ u α σ αβ σ βγ u γ , a 4 ≡ det(σ) [21,23,27], where u µ and σ µν are the top's velocity and angular velocity respectively defined in terms of derivatives with respect to the arbitrary parameter λ (see Refs. [21][22][23][24][25][26][27]). The momentum vector P µ and the antisymmetric spin tensor S µν are canonically conjugated to the position and orientation of the top, P µ ≡ ∂L/∂u µ and S µν ≡ ∂L/∂σ µν = −S νµ . Explicit examples of such Lagrangians can be found in Refs. [22,27]. In this way, it is found that the dynamics of a top describes a non-geodesic behavior, seen through the momentum equation [23,26,27] and the equation for the spin tensor The overdot represents the derivative with respect to an arbitrary parameter (λ), in such a way that velocity u µ =ẋ µ is the derivative of coordinates. In addition, Γ ν ρτ are the Christoffel symbols for the metric field g µν (the speed of light is set equal to 1). The six independent components of the antisymmetric spin tensor generate Lorentz transformations, and in order to restrict them to generate three dimensional rotations we impose the Tulczyjew constraint S µν P ν = 0 [19,22,23,53]. This constraint has been shown to be important in the consistency of a theory for spinning massive particles [19], as it can be deduced as a constraint which emerges from the Lagrangian of the theory, and not an external imposition on the top dynamics [22] (there are other Hamiltonian formulations which do not require these constraints [54]). Lastly, in this theory, the (square) top spin J 2 ≡ 1 2 S µν S µν can be shown to be a conserved quantity [23][24][25][26][27].
The non-geodesic behavior of a top moving on a background gravitational field is determined by Eqs. (1) and (2), plus the constraint. As a result, the top can be interpreted as an extended object that is subject to tidal forces due to gravity. Spin gives internal structure to the classical massive particles, and they cannot be longer described as pointlike objects. Due to the fact that any extended object is crossed by infinitely many geodesics (only a pointlike object is traversed by just one geodesic) the averaged motion does not align with any of the constituent geodesics, and the motion is, in general, non-geodesic. Similar effects have been studied for fields [55,56] (which are naturally extended objects) and electromagnetic waves [55,57,58]. Thus, one should expect that the inclusion of spin in the dynamics of massive particles should lead to non-geodesic orbits.

B. Comparison with Mathisson-Papapetrou theory
The above theory does not coincide with the MPE for massive spinning particles [12,13]. As it is discussed in Ref. [19], the MPE present several problems in their description for tops. The MPE formalism is composed by the third-order dynamical equations where P µ is the momentum defined in that theory, and D/Ds is the s-parametrized covariant derivative for the proper time s of the spinning particle. It is not difficult to obtain [19] that the system preserves the timelike behavior of the velocity u µ u µ = 1. However, the MPE establishes that there is no dynamical conservation law for the mass [19]. In fact, it can be proved that the MPE, under the Pirani constraint S µν u ν = 0, the square of momentum P µ P µ is not a constant of motion [19]. On the other hand, under the Tulczyjew constraint, it can be shown that the MPE implies that the mass is not a constant of motion [19]. Both results of the MPE represent a very undesired behavior for a relativistic theory of massive particles. We refer the reader to Ref. [19], where a deep discussion on the difficulties of the MPE formalism is presented. Without a Lagrangian formulation the canonical momentum and the spin tensor cannot be appropriately defined. That is the reason the MPE equations have those difficulties, whereas the above Lagrangian theory does not.

C. Integrability and non-integrability in Schwarzschild background
Several authors have discussed the integrability and non-integrability of different models for spinning objects in a Schwarzschild field background. For example, the non-integrability of system formed by a Schwarzschild black hole orbited by a spinning companion in the Mathisson-Papapetrou formalism under the extreme-mass-ratio limit was proved in Refs. [51,52], showing the chaoticity of the model. Furthermore, in Ref. [42,59] was found the onset of chaos in the post-Newtonian Lagrangian formulation of the two-black hole system with one body spinning, whereas the post-Newtonian Hamiltonian formulation of the two-black hole system with one body spinning was shown to be integrable and regular [60,61]. By using the construction of canonical conjugate spin variables and the relation between Lagrangian and Hamiltonian approaches at the same post-Newtonian order [62,63], it has been proved that these two models are integrable and non-chaotic, regardless of the post-Newtonian order [64]. On the other hand, a system formed by a Schwarzschild black hole orbited by a spinning companion in the extreme-mass-ratio limit, described by a Lagrangian theory, can be chaotic [42].
Here we present only the key steps to obtain the solution for a circular motion of the top with spin perpendicular to the plane of motion. We refer to readers to those references for a full and detailed procedure to get the solutions for the equations of motion derived from the Lagrangian theory.
Assuming a Schwarzschild field background (describing approximatelly the Earth gravitational field), the equatorial motion of a top can be solved exactly, as any equatorial plane can be defined for circular motion to take place. We write the metric in spherical coordinates g tt = 1 − 2r 0 /r, g rr = − (1 − 2r 0 /r) −1 , g θθ = −r 2 , g φφ = −r 2 sin 2 θ, where r 0 = GM with the gravitational constant G and the Earth mass M . The circular motion is defined as such byṙ = 0. Besides, without any loss of generality, we can study the the motion in the plane defined by θ = π/2. If the top is initially in that plane andθ = 0, then it remains in that equatorial plane, where P θ = 0 [21,23,27]. In this solution, spin can be chosen to be orthogonal to the equatorial plane S rθ = S θφ = S 0θ = 0 [21,23,27], being parallel or antiparallel to the total angular momentum of the top along the whole trajectory. This total momentum angular j for the top in this trayectory is a conserved quantity. In Refs. [21,23,27] is shown that the general solutions for the momenta equations (1) are P φ = (−j ± EJ/m)/(1 − η), and P t = [E ∓ jJr 0 /(mr 3 )]/(1 − η), with the dimensionless parameter η = J 2 r 0 /(m 2 r 3 ). Here, the ± stands for two trajectories that depend on the spin orientation, parallel or antiparallel to the total angular momentum of the top, both of them remaining perpendicular to the plane of motion. These two momenta are conserved (Ṗ t = 0 andṖ φ = 0) for circular motion [21,27], and thus, the Eötvös ratio is meaningless for this particular orbit. The non-geodesic motion manifests itself in changes of the velocity, not acceleration.
From the constant of motion P µ P µ = m 2 , we get that P r = 0 = P 2 , in consistency with the circular motion solution, and the relation between the radial momentum and the radial velocitẏ r = (1 − 2r 0 /r) (P r /P t ) = 0, given by solutions of the Lagrangian theory for tops [21,23,27]. This constraint determines the energy of each trajectory of a top moving on a circular orbit of radius r. The solutions of Eq. (4) for energy can be readily obtained as where Circular orbits are obtained by the study of the behavior of the effective potential V ± (r) = E ± (r) 2 [65]. By calculating the zeros of the derivative of the effective potential a condition for the radius of the orbits can be obtained, while their stability can be studied through the positive behavior of its second derivative. If spin is neglected (Λ 1 = 1 and Λ 2 = 0), the effective potential becomes simply V = (1 − 2r 0 /r) j 2 /r 2 + m 2 (1 − η) 2 , and from the zeros of the its derivative we obtain a relation between the allowed radius and the angular momentum of spinless particles [65] From where we obtain that circular orbits are allowed only for r > 3r 0 . By checking the second derivative of the effective potential, the stability of circular orbits are restricted only to r > 6r 0 . When spin is included, the condition for circular orbits can be obtained by the zeros of the first derivative of the complete solution for energies (5). The calculation is not straightforward, however for a small spin contribution J j, and for tops far from the black hole r 0 r, the allowed stable circular orbits are still for r 6r 0 . This is the case for the estimations of the trayectory deviations for tops around Earth studied in next section.
The interplay of its spin with gravity, introduces different corrections in this tangential velocity, which depends on the spin orientation, such that the top with antiparallel spin is faster than the one with the parallel spin to the total angular momentum rφ − > rφ + . A possible measurement of a maximal manifestation of non-geodesic motion can be achieved if two test bodies (with equal composition) are set to rotate in order to have opposite (internal) angular momenta directions, parallel and antiparallel to the total angular momentum of the circular motion of the satellite. In such cases, any deviation from geodesic orbits must be reflected in different measurements of the angular velocities of the test bodies. Thus, the dimensionless ratio is non-zero in the case of circular orbit for the non-geodesic behavior of tops. Notice that for these circular motions described by solution (8), tops do not present tangential accelerationsφ = 0, as the radius remains constant. Therefore, there is no relative acceleration between the two test bodies. Furthermore, when spin is neglected J = 0 (η = 0), a massive particle can only have a unique angular velocityφ = j/(r 2 E) and an unique energy, yielding the usual result δ = 0 for geodesic motion in the Schwarzschild field [66]. The approximately vanishing Eötvös ratio and δ = 0 are the results measured in the MS mission [5], which is in agreement with the Lagrangian theory for tops.

IV. ESTIMATIONS OF THE TRAJECTORY DEVIATIONS
The inclusion of spin into the the test bodies is essential for experiments carried out to demonstrate the validity of non-geodesic motion of massive spinning bodies. To study the first order corrections to the circular orbit of test bodies due to its spin, let us consider the following approximations. Let us study tops motion with J j. In this case, the total momentum angular j is approximately by the orbital momentum angular of the top.
At first order in spin, the energy solutions of Eq. (4) is With this solution, we can obtain the angular velocity (8) of tops at first order in spin to bė Therefore, at first order in spin, the δ-parameter (9) becomes Now, consider a possible experiment near the Earth surface (where r 0 ∼ 4.4 × 10 −3 [m] r ∼ 7 × 10 6 [m]) in order to measure the non-geodesic behavior of a top in circular trayectory. This kind of experiment have been perfomed by the MS mission [5]. Let us consider a model of two test bodies with the same total angular momentum (which is a conserved quantity for each one) but different spin orientations. Thus, considering a small spin contribution and r 0 r, the top total angular momentum can be approximated by Eq. (7). Notice that this angular momentum is consistent with Kepler's law of motion, as for a general angular momentum j = mrv (with test top's velocity v much smaller than the speed of light) we can write v = rΩ for a circular orbit, where Ω is the angular frequency of the satellite's trajectory. By Kepler's law, this frequency is related to the trajectory radius by Ω 2 = r 0 /r 3 . Using this, the angular momentum becomes j ≈ m √ rr 0 , which is an approximation to Eq. (7).
Thus, the energy of top (10) at first order becomes while the approximated angular velocity (11) at first order iṡ Notice that the gravitational correction factor 3r 0 /2r to the spin-coupling also appears in post-Newtonian theories [67]. Lastly, the parameter (12) at first order becomes simply The behavior of this parameter is shown in Fig. 1, in terms of dimensionless quantities. We have plotted δ/J as function of the normalized radius z = r/r 0 , where J = J/(mr 0 ). The plot is presented for z 1, according to the used approximations. It is straightforward to show from (15) that δ/J = (1 − 3/z)/ √ z, implying that the parameter has the same behavior for any given spin, depending only on the orbit radius r. From the plot, we can see that δ decays as 1/ √ z. We can estimate in a closer form the parameter (15) for a test spinning body. Assuming that the intrinsic spin of each top can be estimated as J ∼ md 2 ω, where d is a characteristic length of the experimental test top body and ω is its internal angular frequency of rotation. Thus, the correction to the circular trajectory of two massive test objects with different spin orientations can be estimated through the ratio (15) to be δ ≈ (d 2 ω/ √ r 0 r)(1 − 3r 0 /r).
This approximated δ-ratio depends on the the test top parameters r, ω and d. Assuming that the test tops have dimensions of the order of d ∼ 10 −2 [m], then the ratio gives of the order where the internal angular frequency of tops is measured in [Hz]. For a higher intrinsic angular velocity, a larger deviation from a geodesic path can be achieved.

V. CONCLUSIONS
The main purpose of this work is to bring attention to the notion of non-geodesic motion, showing that it is the spin (internal dynamics) and not the composition, size or shape of the test bodies, the key for detecting a possible non-geodesic trajectory.
In the simplest case of a circular motion for a spinning massive particle, the complete dynamics can be solved exactly, and it describes a non-geodesics motion whose deviations depend on the magnitude and direction of the particle's spin. Notice that the solution presented above has constant momenta, and thus, the non-geodesic motion cannot be detected by measuring the Eötvös ratio, as the non-geodesic motion is manifested only in the change of the velocities of test bodies. Of course, this is not a general rule, and more complicated motions in different spacetimes may indeed present non-vanishing Eötvös ratios [21,27,29]. As the δ-ratio show in Eqs. (9), when spin is negligible then geodesic motion is expected. Therefore, it comes as no surprise that the EP is repeatedly confirmed for non-spinning test bodies.
Any experimental setup designed to measured non-geodesic motion must be constructed in order to capture the spin-gravity coupling and its effects. In particular, for experiments around Earth the δ-ratio has enough freedom to adjust the parameters of the orbit radius of the satellite and the characteristic length and inner angular velocity of test body, depending on the accuracy of the setup. Thus, any experimental setting should also consider the angular momenta of test bodies in order to prove the validity of non-geodesic motion due to spin.