Revisiting the gravitomagnetic clock e ﬀ ect

To the ﬁrst post-Newtonian order, if two test particles revolve in opposite directions about a massive, spinning body along two circular and equatorial orbits with the same radius, they take di ﬀ erent times to return to the reference direction relative to which their motion is measured: it is the so-called gravitomagnetic clock e ﬀ ect. The satellite moving in the same sense of the rotation of the primary is slower, and experiences a retardation with respect to the case when the latter does not spin, while the one circling in the opposite sense of the rotation of the source is faster, and its orbital period is shorter than it would be in the static case. The resulting time di ﬀ erence due to the stationary gravitomagnetic ﬁeld of the central spinning body is proportional to the angular momentum per unit mass of the latter through a numerical factor which so far has been found to be 4 π . A numerical integration of the equations of motion of a ﬁctitious test particle moving along a circular path lying in the equatorial plane of a hypothetical rotating object by including the gravitomagnetic acceleration to the ﬁrst post-Newtonian order shows that, actually, the gravitomagnetic corrections to the orbital periods are larger by a factor of 4 in both the prograde and retrograde cases. Such an outcome, which makes the proportionality coe ﬃ cient of the gravitomagnetic di ﬀ erence in the orbital periods of the two counter-revolving orbiters equal to 16 π , conﬁrms an analytical calculation recently published in the literature by the present author.


Introduction
Within the weak-field and slow-motion approximation of the General Theory of Relativity (GTR), by gravitomagnetic clock effect it is usually meant the difference ∆T gvm between the orbital periods of two counter-revolving test particles, i.e. nonspinning, moving along circular orbits of identical radius r 0 in the equatorial plane of a rotating body of mass M and angular momentum J (Mitskevich and Pulido Garcia 1970;Vladimirov et al. 1987;Cohen and Mashhoon 1993;Mashhoon 1997;You 1998;Mashhoon et al. 1999;Tartaglia 2000c;Mashhoon and Santos 2000;Mashhoon et al. 2001;Hackmann and Lämmerzahl 2014;Scheumann et al. 2023).It turns out that ∆T gvm is proportional to J/ M c 2 , where c is the speed of light in vacuum, through a numerical factor that has been calculated in the literature to be equal to 4π.Such an intriguing relativistic feature of motion was the subject of several papers investigating its possible detection as well; see Gronwald et al. (1997); Mashhoon et al. (1999); Lichtenegger et al. (2000); Tartaglia (2000a,b); Iorio (2001b,a); Iorio and Lichtenegger (2005); Lichtenegger et al. (2006).It has also relevant consequences in astrophysical contexts such as Kerr black hole spacetimes (de Felice 1995; Bonnor and Steadman 1999;Bini and Jantzen 2003;Faruque 2004;Bini et al. 2005).For other versions of the gravitomagnetic clock effect involving spinning orbiters in the Kerr spacetime, see, e.g., Bini et al. (2004); Mashhoon and Singh (2006).
The standard approach in deriving the aforementioned form of the gravitomagnetic clock effect is to calculate the time interval T gvm required to a test particle to come back to some fixed reference direction in the orbital plane from which it began its motion, assumed circular throughout the overall variation of the azimuthal angle ϕ reckoned from such a line and spanning an interval of 2π, when the general relativistic gravitomagnetic acceleration is added to the Newtonian inverse-square one.The unit vector Ĵ of the primary's angular momentum is assumed to be known, so that one can align the reference z axis with it, and the reference {x, y} plane coincides with the equatorial one of the source.It is assumed the point of view of a distant observer fixed with respect lorenzo.iorio@libero.itarXiv:1407.5022v5[gr-qc] 14 Dec 2023 to distant stars who uses the coordinate time t as own proper time; the difference of the proper times τ of the counter-orbiting test particles is identical to T gvm up to corrections of order O c −4 .Indeed, it is dτ = √ g 00 dt √ 1 + h 00 dt (1 + h 00 /2) dt, where g 00 is the "00" component of the spacetime metric tensor, and h 00 is a small correction of post-Newtonian (pN) order.
In this paper, we will numerically check such an established result by integrating the equations of motion of a fictitious test particle orbiting a putative spinning body by removing the foregoing restrictions about Ĵ. Stated differently, a circular and equatorial orbit will be considered, but its orientation in the reference frame adopted will be arbitrary.It turns out that, in such a scenario, the line of the nodes, among other things, remains fixed; then, it will be naturally assumed as reference polar axis from which the azimuthal angle ϕ is counted in the orbital plane.Furthermore, the argument of latitude u, which reckons the instantaneous position of the test particle along its orbit just from the unit vector l of the line of the nodes, will play the role of ϕ.Thus, the time interval1 between two consecutive crossings of the fixed line of the nodes by the test particle will be inspected by looking at the time needed to the cosine of the angle between its position vector r and l to assume again its initial value.The numerical integration will be performed for both the senses of motion of the test particle.As a result, a discrepancy larger than the established form of the gravitomagnetic clock effect by a multiplicative factor of2 4 will be found, in agreement with an analytical calculation recently appeared in the literature (Iorio 2016).
The paper is organized as follows.In Section 2, the gravitomagnetic acceleration experienced by a satellite orbiting a rotating primary and some of its main orbital effects are reviewed.The gravitomagnetic features of circular and equatorial orbits arbitrarily oriented in space are outlined in Section 3. Section 4 is devoted to the numerical integrations of the equations of motion in the scenario of Section 3. Section 5 deals with a tentative explanation of the resulting discrepancy with respect to the standard case.The impact of the gravitomagnetic field on the interval between two consecutive passages of the test particle at the periapsis, known as anomalistic period, is numerically investigated in Section 6.The possibility of using the existing Earth's artificial satellites LAGEOS and LARES 2 to detect a certain form of gravitomagnetic clock effect in view of their peculiar orbital configuration is investigated in Section 7. Section 8 summarizes the obtained results and offers concluding remarks.Appendix A includes a list of the symbols and definitions used in the paper.

Gravitomagnetic effects for a circular and equatorial orbit
Let us consider an orbital configuration such that the satellite's angular momentum h is (anti-)parallel to the body's spin axis, irrespectively of the orientation of the latter one in the reference frame adopted, i.e. such that see Figure 1 and Figure 2.
In this case, according to Equations ( 5)-( 6), the satellite's orbital plane stays unchanged in space, being aligned with the equator of the primary.The same occurs also if the latter is oblate since the resulting Newtonian precessions of the node and the inclination are proportional to Ĵ • m and Ĵ • l, respectively (Iorio 2017).In particular, the line of nodes remains fixed, and l is constant; thus, it can be naturally assumed as a reference direction in the orbital plane.Instead, in general, the line of the apsides changes because of Equation ( 7).The argument of latitude u, counted just from l taken as reference polar axis, can be used as polar angle reckoning the instantaneous position of the test particle along its orbit.According to Equations (2)-(4), if Equations ( 8)-(9) hold, the normal component of the LT acceleration vanishes, and the gravitomagnetically perturbed motion is entirely in-plane.
If the orbit is also circular, then Equation ( 1) is only radial, reducing to The "+" and "-" signs in Equation ( 10) correspond to the prograde Ĵ • ĥ = +1 and retrograde Ĵ • ĥ = −1 motion with respect to the sense of rotation of the primary, respectively; Figure 1 shows the first case, while the second one is depicted in Figure 2.
It can be noted that when the test particle orbits the central body along the same sense of the rotation of the latter, i.e. for the "+" sign, Equation ( 10) is radially directly outward, weakening the overall gravitational pull to which the test particle is subjected; thus, an increase in the time required to complete a full orbital revolution, which, in this case, can be naturally assumed as the time interval between two consecutive crossings of the fixed line of the nodes, is expected since the motion is slower.On the contrary, when the satellite moves in the opposite direction with respect to the sense of rotation of its primary, the overall gravitational tug experienced by the former is enhanced since Equation ( 10) is radially directed inward (" − " sign), the motion is faster, and the completion of an orbital revolution is expected to take less time.Mitskevich and Pulido Garcia (1970) obtained, incorrectly, the opposite result, while Vladimirov et al. (1987) agreed with the picture just outlined.
In the next Section, we will quantitatively assess such distinctive features of motion, whose non-Machian nature was discussed by Mashhoon et al. (2001), by numerically integrating the equations of motion of a fictitious test particle orbiting a hypothetical massive, spinning body.Let us consider a fictitious test particle orbiting a massive, spinning body with the same mass and spin axis orientation of, say, Jupiter, whose relevant physical parameters are listed in Table 1.It is assumed that the satellite moves along a circular and equatorial prograde orbit starting from the ascending node along the line of the nodes, i.e., Ĵ • ĥ = +1 and r0 • l = +1.Its equations of motion, in Cartesian coordinates, are numerically integrated with and without the LT acceleration of Equation ( 1) starting from the same initial conditions.The angular momentum J of the fictitious primary is assumed large enough to easily visualize its gravitomagnetic effects remaining, at the same time, to a level compatible with the 1pN approximation.In particular, a value for J is adopted such that the ratio of Equation ( 10) to the Newtonian monopole term is as little as 0.005.Then, a time series for the cosine of the angle between the line of the nodes and the test particle's position vector, i.e. r • l, is numerically produced in both the runs and plotted versus time in order to reckon the instant of time when l is crossed again after the initial instant, i.e., when r • l = +1 again.Figure 3 displays both the Keplerian and the gravitomagnetic time series over one orbital revolution.It turns out that, when Equation ( 1) is included in the equations of motion, the passage of the test particle at the line of the nodes occurs after an amount of time equal to the Keplerian orbital period increased by the gravitomagnetic correction

L. Iorio
L. Iorio which, in the fictional example considered, amounts to the 1% of P Kep ; such a retardation is represented by the shaded area in Figure 3.
Figure 4 depicts the same scenario, but for a retrograde orbit, i.e., with Ĵ • h = −1.In this case, the crossing of the line of the nodes occurs in advance with respect to the purely Keplerian case by the same quantity of Equation ( 11); such an advance is represented by the shaded area in Figure 4. Further numerical integrations showed that by varying the orbital radius does not affect the magnitude of the gravitomagnetic retardation/advance.
Figure 5 deals with the instant of time t when the two counter-orbiting test particles meet each other for the first time after their common start.In the Keplerian case, such an event would occur just after half a revolution at ϕ Kep = π.Instead, when the LT acceleration of Equation ( 1) is present, it is expected that the retrograde particle, which moves faster, reaches the prograde one, which is slower, before than in the static case.In Figure 5, the numerically produced time series of r± • l, obtained by including Equation (1) in the equations of motion of both particles, are plotted together, and the time when they cross is marked with a vertical dashed green line.It turns out that with gvm :=  1) for a circular and equatorial orbit arbitrarily oriented in space starting from the ascending node with Ĵ • ĥ = +1 (prograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the line of the nodes, which is fixed both in the Keplerian and LT cases, differ by +8π J/ M c 2 which, in units of P Kep , amounts to 0.01.The shaded area represents such a retardation with respect to the Keplerian period.
Theory predicts just what we came up with in Section 4, as shown by Eq. ( 74) and Eq.(77) of Iorio (2016).Such an outcome has an impact also on the gravitomagnetic clock effect which, now, is For the Earth, since its angular momentum per unit mass amounts to (Petit and Luzum 2010) J ⊕ /M ⊕ 9 × 10 8 m 2 s −1 , Equation ( 14) returns ∆T ⊕ gvm 5 × 10 −7 s.The angular momentum of Jupiter is (Soffel et al. 2003) J Jup 6.9 × 10 38 kg m 2 s −1 ; then, one has ∆T  1) for a circular and equatorial orbit arbitrarily oriented in space starting from the ascending node with Ĵ • ĥ = −1 (retrograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the line of the nodes, which is fixed both in the Keplerian and LT cases, differ by −8π J/ M c 2 which, in units of P Kep , amounts to −0.01.The shaded area represents such an advance with respect to the Keplerian period.
plane of the supermassive black hole in Sgr A * at the Galactic Center, Equation ( 14), calculated for (Daly et al. 2023) χ • 0.90 and (Gillessen et al. 2009) M • 4.1 × 10 6 M , yields ∆T • gvm 9.13 × 10 2 s.Also the analytical calculation of t, obtained by following the reasoning by Mashhoon et al. (2001) extended to the approach by Iorio (2016), agrees with Equation ( 12).One has, first, to compute with the theoretical model used the instant of times t ± corresponding to the fixed angles ϕ and 2π − ϕ for the prograde and retrograde directions, respectively.Then, by imposing the condition t + = t − = t, one gets t and ϕ.By using the calculational scheme in Iorio (2016), it can be obtained  1) for a circular and equatorial orbit arbitrarily oriented in space starting from the ascending node with Ĵ • ĥ = ±1.The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the two counter-orbiting test particles meet at t = P Kep /2 1 − 2 gvm , with gvm given by Equation ( 13), represented by the vertical dashed green line.
see below for further details.By requiring that the right-hand sides of Equations ( 15)-( 16) are equal and solving for u, one gets Inserting Equation ( 17) in any of Equations ( 15)-( 16) yields in agreement with Equation ( 12).
On the other hand, as far as ∆T gvm is concerned, a discrepancy by a factor of 4 occurs with respect to the so-far accepted result3 (Vladimirov et al. 1987;Cohen and Mashhoon 1993;Mashhoon 1997;You 1998;Tartaglia 2000c;Mashhoon et al. 1999;Mashhoon and Santos 2000;Mashhoon et al. 2001) which can be obtained by equating the centripetal acceleration ω 2 r 0 to the sum of the Newtonian monopole plus Equation (10).Indeed, from By integrating Equation ( 22) with respect to ϕ from 0 to +2π for the prograde motion and from 0 to −2π for the retrograde one, the LT orbital period is obtained.Note that ϕ is a polar angle counted from some fixed reference polar axis in the orbital plane aimed to instantaneously locate the test particle along its circular orbit; thus, for an equatorial orbit, it is straightforward to identify the fixed line of the nodes with the reference direction and ϕ with the argument of latitude u.
The explanation in the aforementioned discrepancy likely resides in the fact that the more general calculation in Iorio ( 2016), made by using the non-singular elements q and k, accounts for the fact that, during two consecutive crossings of the line of the nodes, the orbital elements in terms of which dt/du is parameterized, i.e. p, q and k, do actually change instantaneously.In the general case, also the line of the nodes does not stay fixed; such a feature is captured by the calculation in Iorio (2016) as well.It turns out that such an effect does not vanish even in the limit q, k → 0 corresponding to a circular orbit.Indeed, a step-by-step analysis of the calculation in Iorio (2016)  which, for q, k → 0, do not vanish yielding Instead, it turns out that ∂ (dt/du) /∂p ∆p (u) = 0 since, in the limit q, k → 0, the instantaneous variation ∆p (u) of the semilatus rectum p vanishes.Also the term due to the change of the line of the nodes containing dΩ/dt (Iorio 2016) is zero for an equatorial orbit since, according to Equation (4) and Equations ( 8)-( 9), A LT N = 0.The opposite sign in Equations ( 26)-( 27) is obtained for the retrograde motion.It can, now, be noted that Equation ( 15) comes from the sum of L. Iorio calculated with Ĵ • ĥ = +1.On the other hand, Equation ( 16) is the sum of calculated with Ĵ • ĥ = −1.
Instead, the integration based on Equation ( 20) is performed by considering only ϕ as variable during an orbital revolution, all the rest being kept fixed.

The apsidal period
Our numerical tests confirm also another result of Iorio (2016); in the case of an elliptical orbit, there is no LT correction to the apsidal period, i.e. the time interval between two consecutive crossings at the (moving) periapsis, for both the senses of motion.Figure 6 plots the numerically produced time series of the cosine of the angle between the position vector and the Laplace-Runge-Lenz vector A over an orbital revolution with and without the LT acceleration of Equation ( 1) for a highly elliptical and equatorial orbit arbitrarily oriented in space: both the runs share the same initial conditions corresponding to r0 • A 0 = +1 and Ĵ • ĥ = +1 (prograde motion along an equatorial orbit).It turns out that the apsidal periods are the same in both the Keplerian and LT cases.The same occurs also for Ĵ • ĥ = −1 (retrograde motion along an equatorial orbit), as shown by Figure 7.Further numerical integrations show that the anomalistic period is not impacted by the gravitomagnetic field of the primary also for non-equatorial orbits, in agreement with Iorio (2016).
The fact that the gravitomagnetic apsidal period is identical to the Keplerian one can be intuitively justified since, as per Equations ( 5)-( 6), there is not net shift per orbit of the mean anomaly at epoch η.Indeed, from the definition of the mean anomaly Thus, since n Kep stays constant because the semimajor axis is not secularly affected by the gravitomagnetic field, the rate of change of the mean anomaly at epoch is proportional to the opposite of the pace of variation of the time of passage at pericenter according to Should η increase, the crossing of the pericenter would be anticipated with respect to the Keplerian case since t p would decrease, and vice versa.In this case, the variation of η would result in an orbit-by-orbit advance or delay of the passages at the pericenter, which does not occur in the present case because, in fact, dη/dt = 0 for Equation (1).Instead, the 1pN gravitoelectric acceleration due to the mass monopole of the source affects η with a negative rate corresponding to an increase of t p (Iorio 2016).
7. A gravitomagnetic clock effect for LAGEOS and LARES 2; is it possible to measure it?
Strengthened by the numerical confirmations obtained in the previous Sections of our recent analytical results for the pK orbital periods (Iorio 2016), we can trustworthily apply them to some specific situations of interest.
It turns out that a form of gravitomagnetic clock effect can be derived for the laser-ranged geodetic satellites LAGEOS (L) (Cohen and Smith 1985) and LARES 2 (LR 2) (Paolozzi et al. 2019) in view of their peculiar orbital configuration characterized by essentially identical circular orbits lying in orbital planes inclined to the Earth's equator by  1) for an elliptical (e = 0.8) and equatorial orbit arbitrarily oriented in space starting from the periapsis with Ĵ • ĥ = +1 (prograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the apsidal line, which varies when Equation ( 1) is present, are identical in both the Keplerian and LT cases.
respectively.Indeed, although none of them move in the Earth's equatorial plane, their supplementary inclinations make them feasible, at least in principle, to detect the difference of the gravitomagnetic corrections to their sidereal periods meant as the time intervals between two consecutive crossings of a fixed reference direction in space.By aligning J along the reference z axis of a coordinate system whose fundamental plane is parallel to the Earth's equator, Equation (76) and Equation (78) of Iorio (2016), valid to zero order in e, yield  1) for an elliptical (e = 0.8) and equatorial orbit arbitrarily oriented in space starting from the periapsis with Ĵ • ĥ = −1 (retrograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the apsidal line, which varies when Equation ( 1) is present, are identical in both the Keplerian and LT cases.
By calculating Equation (37) with Equations ( 35 Unfortunately, the mutual cancelation of the Keplerian orbital periods P Kep , which should be ideally identical, is far from being perfect.Indeed, for (Ciufolini et al. 2023) a L = 12270.020705km, (39) it is What is worse, ∆P Kep cannot be even known with good enough accuracy to allow it to be safely subtracted from the measured difference of the satellites' orbital periods.Indeed, since the currently claimed uncertainties in the spacecraft's semimajor axes are (Ciufolini et al. 2023) the uncertainty in the difference of the Keplerian periods is as large as which is about six times larger than the gravitomagnetic clock effect of Equation ( 38).Incidentally, the Earth's gravitational parameter µ ⊕ would not be a limiting factor since its error (Petit and Luzum 2010) σ µ ⊕ = 8 × 10 5 m 3 s −2 contributes to σ ∆P Kep at the 6 × 10 −9 s level.Only an improvement of the already unrealistically small errors of Equations ( 42)-( 43) by three orders of magnitude would allow to bring down σ ∆P Kep to the 10 −9 s level.

Summary and conclusions
A recently published analytical calculation of the 1pN gravitomagnetic correction to the draconitic period of a test particle moving along a circular orbit arbitrarily oriented with respect to the equatorial plane of a massive, spinning primary can be specialized to the case in which the angular momentum J of the latter, on whose orientation in space no constraints are assumed as well, is (anti)parallel to the orbital angular momentum h.In such a scenario, in which the orbital plane lies in the equatorial one of the source, the line of the nodes is left unaffected by the pN gravitomagnetic field of the primary, and can be naturally used as reference axis in the orbital plane from which the polar angle ϕ, represented in this case by the argument of latitude u, is reckoned.It turns out that the corrections δT ± gvm , to be added to the Keplerian periods P Kep in order to have the time intervals between two consecutive crossings of the line of the nodes, are larger than the ones present so far in the literature, equal to ±2π J/ M c 2 depending on the prograde or retrograde directions, respectively, by a factor of 4 in both the senses of motion.Also the amount by which the angle ϕ marking the encounter of the two counter-rotating particles is reduced with respect to the Keplerian case is 4 times larger than the value quoted in the literature.
All such features are fully confirmed by numerical integrations of the equations of motion of a fictitious test particle orbiting a putative massive, spinning object and experiencing, to the 1pN order, the gravitomagnetic acceleration A LT induced by J in addition to the Newtonian inverse-square one A N .In particular, the time required by the cosine of the angle between the satellite's position vector r and the unit vector l of the line of the nodes to take its initial value again, assumed equal to unity, turns out to be given by the Keplerian orbital period increased or decreased by 8π J/ M c 2 for the prograde or retrograde senses of motion, respectively.Such a discrepancy with respect to the commonly known scenario may be due to the fact that the latter is based on the strict validity of the circularity condition throughout the integration over one orbital revolution, not allowing any instantaneous variations of the orbital parameters entering the analytical expression of the time derivative of the polar angle.Such a feature is, instead, accounted for in the aforementioned recent calculation which takes into account also the possible motion of the line of the nodes itself occurring, in general, for a non-equatorial orbits.As a result, the revised gravitomagnetic clock effect for identical circular and equatorial orbits traveled in opposite directions turns out to be 16π J/ M c 2 .Such a result has potential relevance to astrophysical observations of the environment of compact objects around which the gravitomagnetic field plays a major roles.
The numerical experiments performed in the present work confirmed also another analytically inferred distinctive feature of the temporal structure around a rotating body: there are no gravitomagnetic corrections to the anomalistic periods for both the senses of motion.It is so because the time of passage at pericenter is left unaffected by the gravitomagnetic acceleration, contrary to the 1pN gravitoelectric one.
Although they do not share the same orbital plane nor do any of them move in the Earth's equatorial plane, the artificial satellites LAGEOS and LARES 2 may be used, at least in principle, to detect a form of clock effect involving the difference of their sidereal orbital periods which is, ideally, entirely gravitomagnetic, amounting to 3 × 10 −7 s.Actually, despite the close values of their semimajor axes, their Keplerian orbital periods do not cancel each other, being their difference uncertain at a 10 −6 s level due to the current claimed errors in the aforementioned orbital elements.

Figure 1 .
Figure1.Prograde circular equatorial orbit arbitrarily oriented in space with, say, I = 30 • , Ω = 45 • .The orbital plane is aligned with the equator of the central body, and the test particle moves along the same sense of rotation of the latter, so that Ĵ • ĥ = +1.

Figure 2 .
Figure 2. Retrograde circular equatorial orbit arbitrarily oriented in space with, say, I = 150 • , Ω = 225 • .The orbital plane is aligned with the equator of the central body, and the test particle moves along the opposite sense of rotation of the latter, so that Ĵ • ĥ = −1.

Figure 3 .
Figure3.Plot of the numerically produced time series of r • l versus time t with and without the LT acceleration of Equation (1) for a circular and equatorial orbit arbitrarily oriented in space starting from the ascending node with Ĵ • ĥ = +1 (prograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the line of the nodes, which is fixed both in the Keplerian and LT cases, differ by +8π J/ M c 2 which, in units of P Kep , amounts to 0.01.The shaded area represents such a retardation with respect to the Keplerian period.

Figure 4 .
Figure4.Plot of the numerically produced time series of r • l versus time t with and without the LT acceleration of Equation (1) for a circular and equatorial orbit arbitrarily oriented in space starting from the ascending node with Ĵ • ĥ = −1 (retrograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the line of the nodes, which is fixed both in the Keplerian and LT cases, differ by −8π J/ M c 2 which, in units of P Kep , amounts to −0.01.The shaded area represents such an advance with respect to the Keplerian period.

Figure 5 .
Figure5.Plots of the numerically produced time series of r+ • l and r− • l versus time t with the LT acceleration of Equation (1) for a circular and equatorial orbit arbitrarily oriented in space starting from the ascending node with Ĵ • ĥ = ±1.The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the two counter-orbiting test particles meet at t = P Kep /2 1 − 2 gvm , with gvm given by Equation (13), represented by the vertical dashed green line.
made with the LT acceleration of Equation (1) shows that Equation (11) comes from the

Figure 6 .
Figure 6.Plot of the numerically produced time series of r • Â versus time t with and without the LT acceleration of Equation (1) for an elliptical (e = 0.8) and equatorial orbit arbitrarily oriented in space starting from the periapsis with Ĵ • ĥ = +1 (prograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the apsidal line, which varies when Equation (1) is present, are identical in both the Keplerian and LT cases.

Figure 7 .
Figure 7. Plot of the numerically produced time series of r • Â versus time t with and without the LT acceleration of Equation (1) for an elliptical (e = 0.8) and equatorial orbit arbitrarily oriented in space starting from the periapsis with Ĵ • ĥ = −1 (retrograde motion).The relevant physical parameters of the putative primary are in Table1.For the sake of clarity, a value of J such that A LT /A N = 0.005 was adopted.It turns out that the time intervals between two consecutive crossings of the apsidal line, which varies when Equation (1) is present, are identical in both the Keplerian and LT cases.

c
: speed of light in vacuum G : Newtonian constant of gravitation M : mass of the central body µ := GM : gravitational parameter of the central body J : angular momentum of the central body J : magnitude of the angular momentum of central body α J : right ascension (RA) of the north pole of rotation of the central body δ J : declination (DEC) of the north pole of rotation of the central body Ĵ = {cos α J cos δ J , sin α J cos δ J , sin δ J } : spin axis of the central body A : perturbing acceleration experienced by the test particle A N : Newtonian inverse-square acceleration ϕ : azimuthal angle reckoning the instantaneous position of the test particle in its orbital plane ω := dϕ/dt : azimuthal angular speed of the test particle in its orbital plane T gvm : orbital period in presence of the LT acceleration a : semimajor axis of the test particle n Kep := µ/a 3 : Keplerian mean motion of the test particle P Kep := 2π/n Kep : orbital period of the test particle e : eccentricity of the test particle p := a 1 − e 2 : semilatus rectum of the orbit of the test particle I : inclination of the orbital plane of the test particle to the reference plane {x, y} Ω : longitude of the ascending node of the test particle : ascending node : descending node ω : argument of pericenter of the test particle q := e cos ω : non-singular orbital element q k := e sin ω : non-singular orbital element k f (t) : true anomaly of the test particle u := ω + f argument of latitude of the test particle η : mean anomaly at epoch t 0 : initial instant of time t p : time of passage at pericenter M (t): mean anomaly r : position vector of the test particle with respect to the central body r : distance of the test particle from the central body r 0 : radius of a circular orbit r := r/r = {cos Ω cos u − cos I sin Ω sin u, sin Ω cos u + cos I cos Ω sin u, sin I sin u} : radial unit vector v : velocity vector of the test particle h = r × v : orbital angular momentum per unit mass ĥ := {sin I sin Ω, − sin I cos Ω, cos I} : unit vector of the orbital angular momentum such that l × m = ĥ A= v × h − µ r: Laplace-Runge-Lenz vector per unit mass l := {cos Ω, sin Ω, 0} : unit vector directed along the line of the nodes toward the ascending node m := {− cos I sin Ω, cos I cos Ω, sin I} : unit vector directed transversely to the line of the nodes in the orbital plane ŝ := ĥ × r = {− sin u cos Ω − cos I sin Ω cos u, − sin Ω sin u + cos I cos Ω cos u, sin I cos u} : transverse unit vector A R := A • r : radial component of the perturbing acceleration A A T := A • ŝ : transverse component of the perturbing acceleration A A N := A • ĥ : normal component of the perturbing acceleration A