To the first postNewtonian order the gravitational action of massenergy currents is encoded by the offdiagonal gravitomagnetic components of the spacetime metric tensor. If they are timedependent a further acceleration enters the equations of motion of a moving test particle. Let the source of the gravitational field be an isolated massive body rigidly rotating whose spin angular momentum experiences a slow precessional motion. The impact of the aforementioned acceleration on the orbital motion of a test particle is analytically worked out in full generality. The resulting averaged rates of change are valid for any orbital configuration of the satellite furthermore they hold for an arbitrary orientation of the precessional velocity vector of the spin of the central object. In general all the orbital elements with the exception of the mean anomaly at epoch undergo nonvanishing longterm variations which in the case of the Juno spacecraft currently orbiting Jupiter and the double pulsar PSR J07373039 AB turn out to be quite small. Such effects might become much more relevant in a starsupermassive black hole scenario as an example the relative change of the semimajor axis of a putative test particle orbiting a Kerr black hole as massive as the one at the Galactic Centre at say 100 Schwarzschild radii may amount up to about inlineformulammlmath idmm1mmlsemanticsmmlmrowmmlmn7mmlmnmmlmommlmommlmrowmmlsemanticsmmlmathinlineformula per year if the holes spin precessional frequency is inlineformulammlmath idmm2mmlsemanticsmmlmrowmmlmn10mmlmnmmlmommlmommlmrowmmlsemanticsmmlmathinlineformula of the particles orbital one.
['classical general relativity', 'experimental studies of gravity', 'experimental tests of gravitational theories', 'General Relativity and Quantum Cosmology', 'Astrophysics - Earth and Planetary Astrophysics', 'Physics - Space Physics']
On the EulerType Gravitomagnetic Orbital Effects in the Field of a Precessing Body
{'On the Euler-type gravitomagnetic orbital effects in the field of a precessing body': "Lorenzo Iorio 1  \n1 Ministero dell' Istruzione e del Merito. Viale Unit'a di Italia 68, I-70125, Bari (BA), Italy", 'Abstract': "To the first post-Newtonian order, the gravitational action of mass-energy currents is encoded by the off-diagonal gravitomagnetic components of the spacetime metric tensor. If they are time-dependent, a further acceleration enters the equations of motion of a moving test particle. Let the source of the gravitational field be an isolated, massive body rigidly rotating whose spin angular momentum experiences a slow precessional motion. The impact of the aforementioned acceleration on the orbital motion of a test particle is analytically worked out in full generality. The resulting averaged rates of change are valid for any orbital configuration of the satellite; furthermore, they hold for an arbitrary orientation of the precessional velocity vector of the spin of the central object. In general, all the orbital elements, with the exception of the mean anomaly at epoch, undergo nonvanishing long-term variations which, in the case of the Juno spacecraft currently orbiting Jupiter and the double pulsar PSR J0737-3039 A/B turn out to be quite small. Such effects might become much more relevant in a star-supermassive black hole scenario; as an example, the relative change of the semimajor axis of a putative test particle orbiting a Kerr black hole as massive as the one at the Galactic Centre at, say, 100 Schwarzschild radii may amount up to about 7% per year if the hole's spin precessional frequency is 10% of the particle's orbital one.  \nKeywords: Classical general relativity; Experimental studies of gravity; Experimental tests of gravitational theories", '1. INTRODUCTION': "In an accelerated reference frame rigidly rotating with time-dependent angular velocity Ω L ( t ), a particle located at position r experiences, among other things, also the fictitious Euler acceleration (Marsden and Ratiu 1999; Morin 2008; Fowles and Cassiday 1999; Battin 1999)  \nA E = -d Ω L d t × r , (1)  \nwhich is often neglected. On the basis of the equivalence principle, at the foundation of the General Theory of Relativity (GTR), one may expect that, within certain limits, an analogous acceleration of gravitational origin should act on a test particle orbiting a rotating body as seen in a local 1 inertial frame attached to the latter. As it will be shown explicitly in the following, it is just the case. After all, the Lense-Thirring (LT) acceleration (Lense and Thirring 1918; Mashhoon et al. 1984) is the general relativistic counterpart, to the first post-Newtonian order (1pN), of the largely known fictitious Coriolis acceleration affecting the motion of an object referred to a rigidly rotating frame.  \nThe spin angular momentum J of an extended, rigidly rotating body of mass M is often displaced, to the Newtonian level, by the differential gravitational tugs exerted on different parts of its centrifugal bulge by other distant masses. It is just the case of the precession and nutation of the Earth mainly due to the torques exerted by the Moon and the Sun (Souchay and Capitaine 2013). Precessional motions of J occur also to the 1pN order when the body of interest moves in the deformed spacetime of other objects (Damour and Ruffini 1974; Barker and O'Connell 1975).", 'lorenzo.iorio@libero.it': "Does the temporal variation of J directly affect the orbital motion of a test particle revolving about its spinning primary? While the answer is negative to the Newtonian level, it is generally affirmative to the 1pN order. If, on the one hand, such general relativistic effects are usually quite small for ordinary restricted two-body systems, on the other hand, they might become relevant when supermassive black holes are involved; be that as it may, it is worthwhile calculating them in full generality.  \nAccording to Brumberg (1991, Eq. (2.2.49), pag. 56), the 1pN equations of motion of a test particle contain, among other things, a time-varying gravitomagnetic acceleration which reads  \nA = c 2 ∂ h ∂x 0 = c ∂ h ∂t . (2)  \nIn Equation (2), c is the speed of light in vacuum, x 0 := ct is the time-like coordinate, and  \nh := { h 01 , h 02 , h 03 } (3)  \nis made of the pN corrections h 0 i , i = 1 , 2 , 3 to the otherwise vanishing off-diagonal components η 0 i = 0 , i = 1 , 2 , 3 of the Minkowskian spacetime metric tensor η µν , µ, ν = 0 , 1 , 2 , 3. They are due to the mass-energy currents of the source of the gravitational field, and are conventionally dubbed as gravitomagnetic. In general, their effects can formally be described in terms of a gravitomagnetic field B g which can be obtained from a gravitomagnetic potential vector (Mashhoon 2001, 2007)  \nA g := c 2 h 2 (4)  \nas  \nand the resulting gravitomagnetic field is  \nB g = 2 G [ J -( J · ˆ r ) ˆ r ] cr 3 . (8)  \nIf J varies over time, Equation (2) and Equation (6) yield  \nA = 2 G c 2 r 3 ∂ J ∂t × r . (9)  \nIn a body-fixed rotating frame with angular velocity ω 0 , the dynamical Euler equations  \nd J d t = ∂ J ∂t + ω 0 × J (10)  \nhold; since an inertial frame is assumed, Equation (10) reduces to  \nd J d t = ∂ J ∂t . (11)  \nAccording to Equation (11), Equation (9) becomes  \nA = 2 G c 2 r 2 d J d t × ˆ r . (12)  \n∇× A g = -B g 2 . (5)  \nIn the case of an isolated, rigidly rotating body, Equation (3) turn out to be (Rindler 2001; Ruggiero and Tartaglia 2002; Poisson and Will 2014)  \nh = 2 G J × r c 3 r 3 , (6)  \nwhere G is the Newtonian constant of gravitation, and r is the position vector of a test particle revolving about the central object, being r their mutual distance. Thus, for such a source, the gravitomagnetic potential can be cast into the form  \nA g = G J × r cr 3 , (7)  \nCan the existence of Equation (12) be guessed on the basis of the equivalence principle, in analogy with the fictitious Euler acceleration of Equation (1)? The gravitational analogue of the Larmor theorem (Mashhoon 1993) tells that the geodesic motion of a test particle in a spatially uniform gravitomagnetic field B g , characterized by 2  \n-4 A g = B g × r , (13)  \noccurs as if it were affected by the Coriolis acceleration 3  \nA C = 2 v × Ω L (14)  \nexperienced in a noninertial frame rigidly rotating with angular velocity  \nΩ L := B g 2 c . (15)  \nBy means of Equation (4), Equation (13) and Equation (15), Equation (2) can be cast just into the form 4  \nA = c ∂ h ∂t = 2 c ∂ A g ∂t = -1 2 c ∂ B g ∂t × r = -1 2 c d B g d t × r = -d Ω L d t × r . (16)  \nThus, the gravitational analogue of the Larmor's theorem holds exactly also when the gravitomagnetic field of the source, assumed spatially uniform, is explicitly time-dependent (Iorio 2002).  \nAim of the paper is the calculation of the impact of Equation (12), viewed as a small correction to the dominant Newtonian inverse-square term, on the orbital motion of a satellite of the spinning body by assuming a purely precessional 5 motion of J ; such a task is accomplished in full generality in Section 2. Numerical evaluations of the size of the resulting effects on the motion of the spacecraft Juno currently orbiting Jupiter and of the double pulsar PSR J0737-3039 A/B are given in Section 3 where also the case of a star orbiting a supermassive Kerr black hole is treated. Section 4 summarizes the findings and offers conclusions.", '2. THE AVERAGED RATES OF CHANGE OF THE KEPLERIAN ORBITAL ELEMENTS': "The net rates of change of the Keplerian orbital elements of a test particle revolving about the spinning primary can be analytically worked out by means of the Gauss equations, valid for any perturbing acceleration A  \nd a d t = 2 n b √ 1 -e 2 [ eA r sin f + ( p r ) A τ ] , (17)  \nd e d t = √ 1 -e 2 n b a { A r sin f + A τ [ cos f + 1 e ( 1 -r a ) ]} , (18)  \nd I d t = 1 n b a √ 1 -e 2 A h ( r a ) cos u, (19)  \ndΩ d t = 1 n b a sin I √ 1 -e 2 A h ( r a ) sin u, (20)  \nd ω d t = √ 1 -e 2 n b ae [ -A r cos f + A τ ( 1 + r p ) sin f ] -cos I dΩ d t , (21)  \nL. Iorio  \nd η d t = -2 n b a A r ( r a ) -( 1 -e 2 ) n b ae [ -A r cos f + A τ ( 1 + r p ) sin f ] . (22)  \nIn Equations (17)-(22), a is the semimajor axis, e is the eccentricity, n b := √ GM/a 3 is the Keplerian mean motion, p := a ( 1 -e 2 ) is the semilatus rectum, I is the inclination of the orbital plane to the reference plane adopted, Ω is the longitude of the ascending node, ω is the argument of pericentre, η is the mean anomaly at epoch, f is the true anomaly, and u := ω + f is the argument of latitude. Furthermore, A r , A τ , A h are the projections of the perturbing acceleration at hand onto the radial, transverse and normal directions, respectively determined by the mutually orthonormal vectors  \nˆ r = { cos Ω cos u -cos I sin Ω sin u, sin Ω cos u +cos I cos Ω sin u, sin I sin u } , (23)  \nˆ τ = {-cos Ω sin u -cos I sin Ω cos u, -sin Ω sin u +cos I cos Ω cos u, sin I cos u } , (24)  \nˆ h = { sin I sin Ω , -sin I cos Ω , cos I } ; (25)  \nˆ h is normal to the orbital plane, being directed along the orbital angular momentum. Equation (23) can be also expressed as  \nˆ r = ˆ l cos u + ˆ m sin u. (26)  \nIn Equation (26),  \nˆ l := { cos Ω , sin Ω , 0 } , (27)  \nˆ m := {-cos I sin Ω , cos I cos Ω , sin I } (28)  \nare two unit vectors lying in the orbital plane; ˆ l is directed along the line of nodes, while ˆ m is perpendicular to ˆ l in such a way that  \nˆ l × ˆ m = ˆ h . (29)  \nThe right-hand-sides of Equations (17)-(22), calculated for the disturbing acceleration under consideration, are to be first evaluated onto the Keplerian ellipse  \nr = p 1 + e cos f , (30)  \nassumed as unperturbed reference trajectory, and then integrated over one full orbital revolution of the test particle by means of  \nd t d f = r 2 √ GMp ; (31)  \nthe averaged rates are finally obtained by dividing the resulting expressions for the aforementioned integrations by the Keplerian orbital period P b = 2 p /n b .  \nBy assuming a purely precessional motion for the spin angular momentum of the primary, i.e. for  \nd J d t = Ω p × J , (32)  \nwhere Ω p is the precession velocity vector of J , and by means of the Binet-Cauchy identity (Gibbs and Wilson 1901, Eq. (25), pag. 76)  \n( C × D ) · ( E × F ) = ( C · E ) ( D · F ) -( C · F ) ( D · E ) , (33)  \nthe radial, transverse and normal components of Equation (12) can be finally cast into the form  \nA r = 0 , (34)  \nA τ = 2 GJK 1 c 2 r 2 , (35)  \nwhere  \nA h = -2 GJ ( K 2 cos u + K 3 sin u ) c 2 r 2 , (36)  \nK 1 := ( Ω p × ˆ J ) · ˆ h , (37)  \nK 2 := ( Ω p · ˆ h )( ˆ J · ˆ l ) -( Ω p · ˆ l )( ˆ J · ˆ h ) , (38)  \nK 3 := ( Ω p · ˆ h )( ˆ J · ˆ m ) -( Ω p · ˆ m ) ( ˆ J · ˆ h ) . (39)  \nBy assuming that ˆ J stays approximately constant during an orbital revolution, the integration of the right-handsides of Equations (17)-(22), calculated with Equations (34)-(36), straightforwardly yields  \nd a d t = 4 GJK 1 c 2 n b a 2 (1 -e 2 ) , (40)  \nd e d t = 2 GJ ( 1 -√ 1 -e 2 ) K 1 c 2 n b a 3 e , (41)  \nd I d t = GJ { K 2 [ -e 2 + ( -2 + e 2 +2 √ 1 -e 2 ) cos 2 ω ] + K 3 ( -2 + e 2 +2 √ 1 -e 2 ) sin 2 ω } c 2 n b a 3 e 2 √ 1 -e 2 , (42)  \ndΩ d t = -GJ csc I { K 3 [ e 2 + ( -2 + e 2 +2 √ 1 -e 2 ) cos 2 ω ] -K 2 ( -2 + e 2 +2 √ 1 -e 2 ) sin 2 ω } c 2 n b a 3 e 2 √ 1 -e 2 , (43)  \nd ω d t = GJ cot I { K 3 [ e 2 + ( -2 + e 2 +2 √ 1 -e 2 ) cos 2 ω ] -K 2 ( -2 + e 2 +2 √ 1 -e 2 ) sin 2 ω } c 2 n b a 3 e 2 √ 1 -e 2 , (44)  \nd η d t = 0 . (45)  \nAccording to Equations (40)-(45), all the orbital elements experience generally nonvanishing secular variations, apart from the mean anomaly at epoch whose precession is identically zero. Furthermore, Equations (40)-(45) have a general validity since they hold for any orbital configuration of the satellite, and for an arbitrary orientation of the primary's spin axis as well. From Equation (37) and Equations (17)-(18) it turns out that the rates of change of the semimajor axis and the eccentricity vanish for equatorial orbits, i.e. if ˆ J = ± ˆ h . On the contrary, the precessions of the inclination, the node and the pericentre do not vanish in such a scenario, as per Equation (32), Equations (38)-(39) and Equations (42)-(44). If the orbit is polar, i.e., if ˆ J · ˆ h = 0, all the rates of the orbital elements do not vanish provided that Ω p does not lie in the orbital plane; in this case, the inclination, the node and the pericentre stay constant, as per Equations (38)-(39) and Equations (43)-(44).  \nIn the limit of small eccentricities, Equations (40)-(45) reduce to  \nd a d t = 4 GJK 1 c 2 n b a 2 + O ( e 2 ) , (46)  \nd e d t = eGJK 1 c 2 n b a 3 + O ( e 2 ) , (47)  \nd I d t = -GJK 2 c 2 n b a 3 + O ( e 2 ) , (48)  \ndΩ d t = -GJ csc IK 3 c 2 n b a 3 + O ( e 2 ) , (49)  \nd ω d t = GJ cot IK 3 c 2 n b a 3 + O ( e 2 ) . (50)  \nd η d t = 0 . (51)  \nIt turns out that, to the lowest order in e , the semimajor axis, the inclination, the node and the pericentre formally undergo secular variations even for circular orbits, while the rate of the eccentricity is proportional to e itself.  \nIn a previous work (Iorio 2002), only the rates of change of the semimajor axis, the eccentricity, the inclination and the node were calculated, to the zeroth order in e , in the particular case of the LAGEOS satellite (Cohen and Smith 1985) orbiting the Earth whose spin axis undergoes the lunisolar precession.", '3.1. Juno and Jupiter': "According to Le Maistre et al. (2016), the Jupiter's pole is precessing due to the gravitational tugs of the Sun, its satellites and the other bodies of the solar system about the normal ˆ w 0 to the Sun-Jupiter invariable plane which can be approximately assumed equal to the solar system's invariable plane (Souami and Souchay 2012). The unit vector ˆ w 0 can be expressed as (Souami and Souchay 2012)  \nˆ w 0 = { sin i p sin θ p , -sin i p cos θ p , cos i p } , (52)  \nwhere (Souami and Souchay 2012)  \ni p glyph[similarequal] 23 · , (53)  \nθ p glyph[similarequal] 3 . 8 · ; (54)  \nsuch figures are referred to the International Celestial Reference Frame (ICRF) having the mean Earth's equator at epoch as reference plane.  \nThe Jovian spin axis can be parameterized as  \nˆ J = { cos α cos δ, sin α cos δ, sin δ } , (55)  \nwhere the nominal values of the R.A. α and decl. δ of the Jupiter's pole are approximately (Le Maistre et al. 2016)  \nα glyph[similarequal] 268 · , (56)  \nδ glyph[similarequal] 64 · . (57)  \nThe spin precession rate of Jupiter is about (Le Maistre et al. 2016)  \nΩ p glyph[similarequal] 3700 mas yr -1 , (58)  \nwhere masyr -1 stands for milliarcseconds per year.  \nBy calculating Equation (25) with the values  \nI = 92 . 99 · , (59)  \nΩ = 267 . 52 · (60)  \nfor the inclination and the node of the Juno spacecraft (Bolton et al. 2017), referred to the ICRF, retrieved from the web interface HORIZONS, maintained by the NASA Jet Propulsion Laboratory (JPL), and inserting Equations (52)-(58) in Equation (37), one obtains  \nK 1 glyph[similarequal] -2 . 5 × 10 -14 s -1 = -0 . 00005 · yr -1 . (61)  \nThus, the rate of the semimajor axis of Juno, obtained from Equation (40) calculated with Equation (61), the value of the Jovian angular momentum  \nJ glyph[similarequal] 6 . 9 × 10 38 kg m 2 s -1 (62)  \nreported in Soffel et al. (2003), and the figures  \na = 4 . 06 × 10 6 km , (63)  \ne = 0 . 981 (64)  \nfor the semimajor axis and the eccentricity of the spacecraft retrieved from HORIZONS, turns out to be  \nd a d t glyph[similarequal] -2 µ myr -1 . (65)  \nThe rate of change of the eccentricity can be calculated with Equation (41) in the same way as just done for the semimajor axis getting  \nd e d t glyph[similarequal] -1 . 5 pas yr -1 , (66)  \nwhere pas yr -1 stands for picoarcseconds per year.  \nSince the calculated values of K 2 and K 3 are similar to that of Equation (61), the precessions of the inclination, the node and the perijove of Juno are of the same order of magnitude of Equation (66).", '3.2. The double pulsar': 'The spin angular momentum J B of the member B of the double pulsar PSR J0737-3039 A/B (Burgay et al. 2003; Lyne et al. 2004) undergoes a precession 6 of 4 . 77 · yr -1 (Breton et al. 2008). The orbital period of the binary is as short as P b = 2 . 45 hr (Kramer et al. 2006), while the spin period of B amounts to T B = 2 . 77 s (Kramer et al. 2006); thus, by assuming for its moment of inertia the standard value  \nI glyph[similarequal] 1 × 10 38 kg m 2 (67)  \nfor neutron stars (Lorimer and Kramer 2005), the size of its angular momentum is of the order of  \nJ B glyph[similarequal] 2 . 2 × 10 38 kg m 2 s -1 . (68)  \nIt turns out that the semimajor axis rate amounts to, at most,  \n∣ ∣ ∣ ∣ d a d t ∣ ∣ ∣ ∣ glyph[lessorsimilar] 0 . 1 mmyr -1 , (69)  \nwhile the orbital precessions are at the glyph[similarequal] 0 . 1 -10 nas yr -1 (nanoarcseconds per year) level.', '3.3. A supermassive black hole-star scenario': "The effects under examination might become relevant in the case of, say, a star orbiting a spinning supermassive black hole (SMBH) whose spin axis ˆ J · undergoes, for some reasons 7 , a relatively fast precession.  \nLet a hypothetical test particle orbits a SMBH at a distance of, say, 100 Schwarzschild radii 8 along an almost circular orbit. It should be recalled that the magnitude of the angular momentum a Kerr black hole (Kerr 1963; Teukolsky 2015) is (Shapiro and Teukolsky 1986)  \nJ · = χ M 2 · G c , (70)  \nwith | χ | ≤ 1. If one assumes, say,  \nM · = 4 . 5 × 10 6 M glyph[circledot] , (71)  \nχ = 1 , (72)  \nwhere M glyph[circledot] is the Sun's mass, it turns out that, according to Equation (46), the particle's semimajor axis would be changed by at most  \n1 a ∣ ∣ ∣ ∣ d a d t ∣ ∣ ∣ ∣ glyph[lessorsimilar] 7%yr -1 , (73)  \nprovided that the hole's spin axis precessional frequency Ω p is 10% of the star's orbital one n b , i.e., if  \nΩ p n b = 0 . 1 . (74)  \nBy relying upon the same assumptions, Equations (48)-(50) tell that the precessions of the other orbital elements would be  \n∣ ∣ ∣ ∣ d κ d t ∣ ∣ ∣ ∣ glyph[lessorsimilar] 1 · yr -1 , κ = I, Ω , ω (75)  \nwhile the rate of change of the eccentricity would be of the order of  \n∣ ∣ ∣ ∣ d e d t ∣ ∣ ∣ ∣ glyph[lessorsimilar] 10 -6 yr -1 , (76)  \nas per Equation (47). Since Equations (46)-(51) are proportional to Ω p , the previous values can be straightforwardly rescaled for different values of Equation (74). As far as the distance from the black hole is concerned, both the relative rate of change of a and Equations (47)-(50) fall as a -3 / 2 .", '4. SUMMARY AND CONCLUSIONS': "The impact of the general relativistic Euler-type gravitomagnetic acceleration induced by the temporal variation of mass-energy currents was analytically calculated, to the first post-Newtonian order, for a restricted two-body system in the hypothesis that the source of the gravitational field is an isolated, massive body rigidly rotating whose spin angular momentum undergoes a purely precessional motion.  \nThe calculation has a full generality since it holds for any orbital configuration of the test particle, and for an arbitrary orientation of the precession velocity vector of the central object as well.  \nIt turns out that, in general, all the orbital elements, apart from the mean anomaly at epoch, undergo long-term variations calculated by assuming that the primary's spin axis precession is much slower than the satellite's orbital revolution.  \nThe resulting effects are usually quite small; for the Juno spacecraft currently orbiting Jupiter, the semimajor axis changes at a rate as small as a few microns per year, while the shifts of the other orbital elements are at the picoarcseconds per year level. For the double pulsar, the resulting figures are larger by just a few orders of magnitude. On the other hand, they may become relevant around a supermassive black hole. Indeed, by assuming a precessional frequency of the hole's spin axis equal to 10% of the mean motion of a putative test particle orbiting it at 100 Schwarzschild radii, it turns out that its semimajor axis changes by up to 7% per year if the mass of the hole at the Galactic Centre is assumed.", 'DATA AVAILABILITY': 'No new data were generated or analysed in support of this research.', 'CONFLICT OF INTEREST STATEMENT': 'I declare no conflicts of interest.', 'REFERENCES': "- B. M. Barker and R. F. O'Connell. Gravitational two-body problem with arbitrary masses, spins, and quadrupole moments. Phys. Rev. D , 12:329-335, 1975. https://doi.org/10.1103/PhysRevD.12.329.\n- R. H. Battin. An Introduction to the Mathematics and Methods of Astrodynamics . American Institute of Aeronautics and Astronautics, 1999.\n- S. J. Bolton, J. Lunine, D. Stevenson, et al. The Juno Mission. Space Sci. Rev. , 213:5-37, 2017. https://doi.org/10.1007/s11214-017-0429-6.\n- R. P. Breton, V. M. Kaspi, M. Kramer, et al. Relativistic Spin Precession in the Double Pulsar. Science , 321:104, 2008. https://doi.org/10.1126/science.1159295.\n- V. A. Brumberg. Essential Relativistic Celestial Mechanics . Adam Hilger, 1991.\n- M. Burgay, N. D'Amico, A. Possenti, et al. An increased estimate of the merger rate of double neutron stars from observations of a highly relativistic system. Nature , 426: 531-533, 2003. https://doi.org/10.1038/nature02124.\n- S. C. Cohen and D. E. Smith. LAGEOS Scientific results: Introduction. J. Geophys. Res. , 90:9217-9220, 1985. https://doi.org/10.1029/JB090iB11p09217.\n- T. Damour and R. Ruffini. Sur certaines v'erifications nouvelles de la Relativit'e g'en'erale rendues possibles par la d'ecouverte d'un pulsar membre d'un syst'eme binaire. C.R. Acad. Sc. Paris, S'erie A , 279:971-973, 1974.\n- G. R. Fowles and G. L. Cassiday. Analytical Mechanics. Sixth Edition . Harcourt, 1999.\n- J. W. Gibbs and E. B. Wilson. Vector analysis. A text-book for the use of students of mathematics and physics . Scribner, 1901.\n- L. Iorio. A Gravitomagnetic Effect on the Orbit of a Test Body Due to the Earth's Variable Angular Momentum. Int. J. Mod. Phys. D , 11:781-787, 2002.  \nhttps://doi.org/10.1142/S0218271802002086.  \n- R. P. Kerr. Gravitational Field of a Spinning Mass as an Example of Algebraically Special Metrics. Phys. Rev. Lett. , 11:237-238, 1963.  \nhttps://doi.org/10.1103/PhysRevLett.11.237.  \n- M. Kramer, I. H. Stairs, R. N. Manchester, et al. Tests of General Relativity from Timing the Double Pulsar. Science , 314:97-102, 2006.\n- https://doi.org/10.1126/science.1132305.\n- S. Le Maistre, W. M. Folkner, R. A. Jacobson, and D. Serra. Jupiter spin-pole precession rate and moment of inertia from Juno radio-science observations. Planet. Space Sci. , 126:78-92, 2016. https://doi.org/10.1016/j.pss.2016.03.006.\n- J. Lense and H. Thirring. Uber den Einfluß der Eigenrotation der Zentralkorper auf die Bewegung der Planeten und Monde nach der Einsteinschen Gravitationstheorie. Phys. Z , 19:156-163, 1918.\n- D. Lorimer and M. Kramer. Handbook of Pulsar Astronomy . Cambridge Observing Handbooks for Research Astronomers. Cambridge University Press, 2005.\n- A. G. Lyne, M. Burgay, M. Kramer, et al. A Double-Pulsar System: A Rare Laboratory for Relativistic Gravity and Plasma Physics. Science , 303: 1153-1157, 2004. https://doi.org/10.1126/science.1094645.\n- J. E. Marsden and T. S. Ratiu. Introduction to Mechanics and Symmetry. A Basic Exposition of Classical Mechanical Systems. Second edition , volume 117 of Texts in Applied Mathematics . Springer, 1999. https://doi.org/10.1007/978-0-387-21792-5.\n- B. Mashhoon. On the gravitational analogue of Larmor's theorem. Phys. Lett. A , 173:347-354, 1993. https://doi.org/10.1016/0375-9601(93)90248-X.\n- B. Mashhoon. Gravitoelectromagnetism. In J. F. Pascual-S'anchez, L. Flor'ıa, A. San Miguel, and F. Vicente, editors, Reference Frames and Gravitomagnetism . World Scientific, 2001.\n- B. Mashhoon. Gravitoelectromagnetism: A Brief Review. In L. Iorio, editor, The Measurement of Gravitomagnetism: A Challenging Enterprise , pages 29-39. Nova Science, 2007.\n- B. Mashhoon, F. W. Hehl, and D. S. Theiss. On the gravitational effects of rotating masses: The Thirring-Lense papers. Gen. Relativ. Gravit. , 16: 711-750, 1984. https://doi.org/10.1007/BF00762913.  \n- D. Morin. Introduction to Classical Mechanics. With Problems and Solutions . Cambridge University Press, 2008.\n- G. Petit and B. Luzum, editors. IERS Conventions (2010) , volume 36 of IERS Technical Note . Verlag des Bundesamts fur Kartographie und Geodasie, Frankfurt am Main, 2010.\n- E. Poisson and C. M. Will. Gravity. Newtonian, Post-Newtonian, Relativistic . Cambridge University Press, 2014.\n- https://doi.org/10.1017/CBO9781139507486.\n- ' E. Racine. Analysis of spin precession in binary black hole systems including quadrupole-monopole interaction. Phys. Rev. D , 78:044021, 2008. https://doi.org/10.1103/PhysRevD.78.044021.\n- W. Rindler. Relativity: special, general, and cosmological . Oxford University Press, 2001.\n- M. L. Ruggiero and L. Iorio. Gravitomagnetic time-varying effects on the motion of a test particle. Gen. Relativ. Gravit. , 42:2393-2402, 2010. https://doi.org/10.1007/s10714-010-0987-3.\n- M. L. Ruggiero and A. Tartaglia. Gravitomagnetic effects. Nuovo Cim. B , 117:743, 2002.\n- M. Sayeb, L. Blecha, L. Z. Kelley, et al. Massive black hole binary inspiral and spin evolution in a cosmological framework. Mon. Not. Roy. Astron. Soc. , 501:2531-2546, 2021. https://doi.org/10.1093/mnras/staa3826.\n- S. L. Shapiro and S. A. Teukolsky. Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects . Wiley, 1986.\n- M. H. Soffel, S. A. Klioner, G. Petit, et al. The IAU 2000 Resolutions for Astrometry, Celestial Mechanics, and Metrology in the Relativistic Framework: Explanatory Supplement. Astron J. , 126:2687-2706, 2003. https://doi.org/10.1086/378162.\n- D. Souami and J. Souchay. The solar system's invariable plane. Astron. Astrophys. , 543:A133, 2012. https://doi.org/10.1051/0004-6361/201219011.\n- J. Souchay and N. Capitaine. Precession and Nutation of the Earth. In J. Souchay, S. Mathis, and T. Tokieda, editors, Tides in Astronomy and Astrophysics , volume 861 of Lecture Notes in Physics , pages 115-166. Springer, 2013. https://doi.org/10.1007/978-3-642-32961-6 4.\n- S. A. Teukolsky. The Kerr metric. Class. Quantum Gravit. , 32:124006, 2015.  \nhttps://doi.org/10.1088/0264-9381/32/12/124006."}