arxiv: v1 [gr-qc] 30 Sep 2018

Transcription

1 A general relativistic gravitomagnetic effect on the orbital motion of a test particle around its primary induced by the spin of a distant third body Lorenzo Iorio 1 Ministero dell Istruzione, dell Università e della Ricerca M.I.U.R.-Istruzione Permanent address for correspondence: Viale Unità di Italia 68, 70125, Bari BA, Italy arxiv: v1 [gr-qc] 30 Sep 2018 lorenzo.iorio@libero.it Received ; accepted

2 2 Abstract We study a general relativistic gravitomagnetic 3rd-body effect induced by the spin angular momentum S X of a rotating mass M X orbited at distance r X by a local gravitationally bound restricted two-body systemsof size r r X consisting of a test particle revolving around a massive body M. We analytically work out the doubly averaged rates of change of the Keplerian orbital elements of the test particle by finding non-vanishing long-term effects for the inclination I, the node Ω and the pericenter ω. We numerically calculate their magnitudes for some astronomical scenarios in our solar system. For putative man-made orbiters of the natural moons Enceladus and Europa in the external fields of Saturn and Jupiter, the relativistic precessions due to the angular momenta of the gaseous giant planets can be as large as milliarcseconds per year mas yr 1. The effects induced by the Sun s angular momentum on artificial probes of Mercury and the Earth are at the level of microarcseconds per year µas yr 1. keywords Experimental studies of gravity; Experimental tests of gravitational theories; Satellite orbits 1. Introduction Let us consider a local gravitationally bound restricted two-body system S composed by a test particle orbiting a planet 1 of mass M at distance r, and a distant 3rd body X mass M X M and proper spin S X around whichsrevolves at distance r X r. In a kinematically rotating and dynamically non-rotating coordinate systemk Brumberg & Kopeikin 1989; Damour, Soffel & Xu 1994 attached to M in geodesic motion through the external spacetime deformed by the mass-energy currents of X, assumed stationary in a kinematically and dynamically non-rotating coordinate systemk X whose axes point towards the distant quasars Brumberg & Kopeikin 1989; Damour, Soffel & Xu 1994, the planetocentric motion of the test particle is affected two peculiar post-newtonian 3rd-body effects: the time-honored De Sitter precession due to solely the mass M X de Sitter 1916; Schouten 1918; Fokker 1920, and a further gravitomagnetic shift due to S X. Our purpose is to analytically work out the latter effect in full generality, out any a-priori simplifying assumptions concerning both the orbital configurations of the planetocentric satellite s motion and the astrocentric motion of the planet-satellite system, and for an arbitrary orientation of S X in space. For previous, approximate calculation restricted to the orbital angular momentum of the Moon orbiting the Earth in the field of the rotating Sun, see Gill et al. 1992, Sec In general, it may be endowed its own gravitational multipole moments and spin affecting the satellite s motion the standard Newtonian and post-newtonian orbital effects.

3 3 The plan of the paper is as follows. In Section 2, we analytically work out the long-term rates of change of the Keplerian orbital elements of the test particle. Section 3 is devoted to the application of the obtained results to some astronomical scenarios in our solar system. We summarize our results and offer our conclusions in Section The doubly averaged satellite s orbital precessions The 3rd-body gravitomagnetic potential induced by the angular momentum S X of the external spinning object X on the planetary satellite is U GM = G S X [ L+3 L ˆr c 2 r 3 X ˆr X ]. 1 X It comes from Eq of Barker & O Connell 1979, p. 155 for the interaction potential energy V S 1,S 2 of two spins S 1, S 2 of masses m 1, m 2 separated by a distance r and moving relative speed v in the limit m 2 M X m 1 M, and by assuming that the spin S 1 is the orbital angular momentum of the planetocentric satellite s motion while S 2 is the spin angular momentum S X of the distant 3rd body X. Thus, r in Eq of Barker & O Connell 1979, p. 155 has to be identified r X, and r P is nothing but the orbital angular momentum r X Mv X of the motion ofsaround M X. It is interesting to note that, the same identifications, V S 1 and V S 2 of Eqs in Barker & O Connell 1979, p. 155 yield the gravitoelectric De Sitter orbital precession for the planetocentric motion of the satellite and the gravitomagnetic Lense-Thirring effect for the astrocentric orbit of M, respectively. The perturbing potential U pert to be inserted into the Lagrange planetary equations for the rates of change of the osculating Keplerian orbital elements of the test particle Bertotti, Farinella & Vokrouhlický 2003; Kopeikin, Efroimsky & Kaplan 2011, obtained by doubly averaging Equation 1 respect to P b, P X b for arbitrary orbital configurations of both the external body X and the test particle and for a generic orientation of S X in space, is U pert = U GM = GS Xn b a 2 1 e 2 3/2 U, 2 2c 2 a 3 X 1 e 2 X U= cos I { 2Ŝ z 3 sin I X [Ŝ z sin I X + cos I X Ŝ y cosω X Ŝ x sinω X ]} + + sin I 2 { 2Ŝ y cosω 2Ŝ x sinω+3 cos Ω Ω X [ Ŝ z sin 2I X sin 2 I X Ŝ y cosω X + Ŝ x sinω X ]}. 3

4 4 The resulting doubly averaged rates of change of the relevant Keplerian orbital elements turn out to be ȧ=0, 4 ė=0, 5 GS X İ= 2a 3 X c2 3/2 I, 6 1 e 2 X GS X Ω= 2a 3 X c2 3/2 O, 7 1 e 2 X GS X csc I ω= 8a 3 X c2 3/2 P. 8 1 e 2 X I=sinΩ { [ Ŝ ]} Ŝ y + 3 sin I X cosω X Ŝ z cos I X + sin I X y cosω X Ŝ x sinω X + + cosω { Ŝ x + 3 sin I X sinω X [Ŝ z cos I X + sin I X Ŝ y cosω X + Ŝ x sinω X ]}, 9 O=2Ŝ z + Ŝ x cot I sinω 3 cos I X sin I X Ŝ y cosω X Ŝ x sinω X + Ŝ z cot I sinωsinω X 3 sin 2 I X [Ŝ z + cot I sinωsinω X Ŝ y cosω X + Ŝ x sinω X ] + + cosωcot I { Ŝ y + 3 sin I X cosω X [ Ŝ z cos I X + sin I X Ŝ y cosω X Ŝ x sinω X ]}, 10 P=Ŝ y [cosω 3 cos Ω 2Ω X ] Ŝ x [sinω+3 sin Ω 2Ω X ]+ + 6 cos Ω Ω X [ Ŝ z sin 2I X + cos 2I X Ŝ y cosω X Ŝ x sinω X ]. 11 We remark that Equations 4 to 11 are exact in both e and e X in the sense that the low-eccentricity approximation was not adopted in the calculation. A more computationally cumbersome approach to obtain the same long-term rates of change of Equations 4 to 11 consists, first of all, in deriving a perturbing acceleration from

5 5 Equation 1. By writing the Lagrangian per unit mass of a gravitationally bound restricted two-body system affected by a generic perturbing potential as the conjugate momentum per unit mass is, by definition, L=L 0 +L pert = v2 2 +µ r +L pert, 12 Thus, The Hamiltonian per unit mass is p L v = v+ L pert v. 13 ṗ= v+ d dt Lpert v. 14 From the Hamilton equations of motion, it is H=H 0 +H pert = v2 2 µ r +H pert. 15 ṗ= H r = µ r r H pert r SinceL pert = H pert, by comparing Equation 14 and Equation 16, it turns out that the perturbing acceleration is just A pert = d Hpert H pert. 17 dt v r In our specific case, sinceh pert = U GM, we have A GM = d UGM U GM = 2G v [S dt v r c 2 r 3 X 3 S X ˆr X ˆr X ]. 18 X Then, Equation 18 must be decomposed into its radial R, transverse T and out-of-plane N

6 6 components, which are A GM R = GS X c 2 r 5 X [ Ŝ x r 2 X + 3Ŝ x x 2 X + 3Ŝ y x X y X + 3Ŝ z x X z X cosω+ + Ŝ y r 2 X + 3Ŝ x x X y X + 3Ŝ y y X 2 + 3Ŝ z y X z X sinω ], 19 A GM T = GS X csc I c 2 r 5 X { sin I [ 3 Ŝ x x X + Ŝ y y X zx + Ŝ z r 2 X 3z X 2] + + cos I [ Ŝ y r 2 X 3y X 2 3 Ŝ x x X + Ŝ z z X yx ] cosω+ + cos I [ Ŝ x r 2 X + 3x2 X + 3Ŝ y x X y X + 3Ŝ z x X z X ] sinω }, 20 A GM N = GS X csc I c 2 r 5 X {[Ŝ y r 2 X 3y X 2 3y X Ŝ x x X + Ŝ z z X ] cosω+ + [ Ŝ x r 2 X + 3x2 X + 3Ŝ y x X y X + 3Ŝ z x X z X ] sinω }. 21 They have to be inserted into the right-hand-sides of the standard Gauss equations for the variation of the orbital elements Poisson & Will 2014 which, finally, are doubly averaged respect to P b, P X b. It should be noted that the averaging procedures previously outlined for both the Lagrange and Gauss equations hold if the orbital elements{ψ} of the test particle can be considered as approximately constant during one orbital period P X b. Indeed, the Keplerian ellipses were used as unperturbed, reference orbits for both the test particle and S itself. This is usually the case in several scenarios of interests, some of which are treated in Section 3. Indeed, the frequencies Ψ of the long-term variations of the test particle s orbital elements{ψ} induced by the usual internal dynamics ofs, determined by, e.g., the quadrupole mass moments J l,l=2, 3, 4... of M and the standard post-newtonian effects like the gravitoelectric Einstein pericenter precession and the Lense-Thirring effect due to the planetary angular momentum S, are often much smaller than the mean motion n X b characterizing the orbital revolution ofsitself around X. However, there may be cases in which the average respect to P X b must take into account also for the relatively fast precessing ellipse of the test particle. For example, some low-altitude Earth s satellites are impacted by the Earth s oblateness J 2 in such a way that the frequency of their node is comparable to, or even larger than the mean motion n b of the heliocentric terrestrial revolution.

7 7 3. Some potentially interesting astronomical scenarios For the sake of simplicity, we will consider a circular orbit e=0 for the test particle motion around M. In the case of a hypothetical orbiter of the Kronian natural satellite Enceladus in the external field of Saturn, Equations 6 to 7 and Equations 9 to 10, referred to the mean Earth s equator at the reference epoch J as reference{x, y} plane, yield İ=A eq sin Ω+ϕ eq, 22 Ω= 49.9 mas yr 1 + cot I A eq cos Ω+ϕ eq, 23 A eq = 5.7 mas yr 1, 24 ϕ eq = 49.4 deg. 25 Instead, if the mean ecliptic at the reference epoch J is adopted as reference{x, y} plane, we have İ= A ecl sin Ω+ϕ ecl, 26 Ω= 34.0 mas yr 1 + cot I A ecl cos Ω+ϕ ecl, 27 A ecl = 23.9 mas yr 1, 28 ϕ ecl = 10.4 deg. 29 By looking at a putative orbiter of the Jovian natural satellite Europa in the external field of Jupiter, we have İ=A eq sin Ω+ϕ eq, 30 Ω= 9.9 mas yr 1 + cot I A eq cos Ω+ϕ eq, 31

8 8 A eq = 4.8 mas yr 1, 32 ϕ eq = 2.9 deg, 33 and 34 İ= A ecl sin Ω+ϕ ecl, 35 Ω= 11.0 mas yr 1 + cot I A ecl cos Ω+ϕ ecl, 36 A ecl = 0.3 mas yr 1, 37 ϕ ecl = 31.0 deg. 38 We considered just Enceladus and Europa because they are of great planetological interest in view of the possible habitability of their oceans beneath their icy crusts Lunine As a consequence, they are the natural targets of several concept studies and proposals for dedicated missions to them, including also orbiters Razzaghi et al. 2008; Spencer & Niebur 2010; MacKenzie et al. 2016; Verma & Margot 2018; Sherwood et al In the case of an artificial satellite orbiting a planet in the field of the Sun, the effects are much smaller. For an Earth s spacecraft, we have İ=A eq sin Ω+ϕ eq, 39 Ω= 0.2 µas yr 1 + cot I A eq cos Ω+ϕ eq, 40 A eq = 0.1 µas yr 1, 41 and ϕ eq = 9.13 deg, 42 İ= A ecl sin Ω+ϕ ecl, 43 Ω= 0.3 µas yr 1 + cot I A ecl cos Ω+ϕ ecl, 44

9 9 A ecl = 0.02 µas yr 1, 45 ϕ ecl = deg. 46 For a probe orbiting Mercury one gets İ=A eq sin Ω+ϕ eq, 47 Ω= 4.3 µas yr 1 + cot I A eq cos Ω+ϕ eq, 48 A eq = 2.5 µas yr 1, 49 and ϕ eq = deg, 50 İ=A ecl sin Ω+ϕ ecl, 51 Ω= 5 µas yr 1 + cot I A ecl cos Ω+ϕ ecl, 52 A ecl = 0.6 µas yr 1, 53 ϕ ecl = deg Summary and overview We analytically worked out the general relativistic gravitomagnetic long-term rates of change of the relevant Keplerian orbital elements of a test particle orbiting a primary M at distance r from it which, in turn, moves in the external spacetime deformed by the spin angular momentum S X of a distant r X r 3rd body X mass M X M. We did not assume any preferred orientation for the spin axis Ŝ X of the external body; moreover, we did not make simplifying assumptions pertaining the orbital configurations of both the M s satellite and of M itself in its motion around M X. Thus, our calculation have a general validity, being applicable to arbitrary astronomical systems of potential interest. It turns out that, by doubly averaging the perturbing

10 10 potential employed in the calculation respect to the orbital periods P b, P X b of both M and M X, the semimajor axis a and the eccentricity e do not experience long-term variations, contrary to the inclination I of the orbital plane, the longitude of the ascending nodeωand the argument of pericenter ω. While the gravitomagnetic rates İ and ω are harmonic signatures characterized by the frequency of the possible variation of the nodeω, induced by other dominant perturbations like, e.g., the Newtonian quadrupole mass moment of the satellite s primary M, the gravitomagnetic node rate Ω exhibits also a secular trend in addition to a harmonic component the frequency of the node itself. The Sun s angular momentum exerts very small effects on spacecraft orbiting Mercury 1 µas yr 1 and the Earth 0.1 µas yr 1. Instead, the angular momenta of the gaseous giant planets like Jupiter and Saturn may induce much larger perturbations of the orbital motions of hypothetical anthropogenic orbiters of some of their major natural moons like, e.g., Europa 10 mas yr 1 and Enceladus 50 mas yr 1. Such natural satellites have preeminent interest in planetology making them ideal targets for future, dedicated spacecraft-based missions which may be opportunistically exploited to attempt to measure such relativistic effects as well.

11 11 REFERENCES Barker B. M., O Connell R. F., 1979, Gen. Relat. Gravit., 11, 149 Bertotti B., Farinella P., Vokrouhlický D., 2003, Physics of the Solar System. Kluwer Academic Press, Dordrecht Brumberg V. A., Kopeikin S. M., 1989, Nuovo Cimento B, 103, 63 Damour T., Soffel M., Xu C., 1994, Phys. Rev. D, 49, 618 de Sitter W., 1916, Mon. Not. Roy. Astron. Soc., 77, 155 Fokker A. D., 1920, Versl. Kon. Ak. Wet., 29, 611 Gill E., Soffel M., Ruder H., Schneider M., 1992, Relativistic Motion of Gyroscopes and Space Gradiometry. Deutsche Geodätische Kommission, München Kopeikin S., Efroimsky M., Kaplan G., 2011, Relativistic Celestial Mechanics of the Solar System. Wiley-VCH, Weinheim Lunine J. I., 2017, Acta Astronaut., 131, 123 MacKenzie S. M. et al., 2016, Adv. Space Res., 58, 1117 Petit G., Luzum B., et al., 2010, IERS Technical Note, 36, 1 Poisson E., Will C. M., 2014, Gravity. Cambridge University Press, Cambridge Razzaghi A. I., di Pietro D. A., Quinn D. A., Simon-Miller A. A., Tompkins S. D., 2008, in American Institute of Physics Conference Series, Vol. 969, Space Technology and Applications International Forum-STAIF 2008, El-Genk M. S., ed., pp Schouten W. J. A., 1918, Versl. Kon. Ak. Wet., 27, 214 Seidelmann P. K. et al., 2007, Celestial Mechanics and Dynamical Astronomy, 98, 155 Sherwood B., Lunine J., Sotin C., Cwik T., Naderi F., 2018, Acta Astronaut., 143, 285 Soffel M. et al., 2003, AJ, 126, 2687 Spencer J., Niebur C., 2010, Planetary Science Decadal Survey. Enceladus Orbiter. National Aeronautics and Space Administration Verma A. K., Margot J.-L., 2018, Icarus, 314, 35 This manuscript was prepared the AAS L A TEX macros v5.2.

12 12 Appendix A Notations and definitions Here, some basic notations and definitions pertaining the restricted two-body system S moving in the external field of the distant 3rd body X considered in the text are presented. For the numerical values of some of them, see Tables 1 to 2. G : Newtonian constant of gravitation c : speed of light in vacuum ǫ : mean obliquity M X : mass of the distant 3rd body X a star like the Sun or a planet like, e.g., Jupiter or Saturn µ X GM X : gravitational parameter of the 3rd body X S X : magnitude of the angular momentum of the 3rd body X Ŝ X = { Ŝ x, Ŝ y, Ŝ z } : spin axis of the 3rd body X in some coordinate system α X : right ascension RA of the 3rd body s spin axis δ X : declination DEC of the 3rd body s spin axis Ŝ eq x Ŝ eq y Ŝ eq z = cosδ X cosα X : component of the 3rd body s spin axis w.r.t. the reference x axis of an equatorial coordinate system = cosδ X sinα X : component of the 3rd body s spin axis w.r.t. the reference y axis of an equatorial coordinate system = sinδ X : component of the 3rd body s spin axis w.r.t. the reference z axis of an equatorial coordinate system r X : position vector towards the 3rd body X r X : distance ofsto the 3rd body X ˆr X r X /r X : versor of the position vector towards the 3rd body X a X : semimajor axis of the orbit about the 3rd body X n X b µ X /a 3 X : mean motion of the orbit about the 3rd body X P X b 2π/nX b : orbital period of the orbit about the 3rd body X e X : eccentricity of the orbit about the 3rd body X

13 13 I X : inclination of the orbital plane of orbit about the 3rd body X to the reference{x, y} plane of some coordinate system Ω X : longitude of the ascending node of the orbit about the 3rd body X referred to the reference {x, y} plane of some coordinate system M : mass of the primary planet or planetary natural satellite orbited by the test particle and moving in the external field of the 3rd body X µ GM : gravitational parameter of the primary orbited by the test particle and moving in the external field of the 3rd body X S : angular momentum of the primary J l,l=2, 3, 4 : zonal multipole moments of the classical gravitational potential of the primary r : position vector of the test particle respect to its primary r : magnitude of the position vector of the test particle v : velocity vector of the test particle L r v : orbital angular momentum per unit mass of the test particle a : semimajor axis of the test particle s orbit n b µ/a 3 : Keplerian mean motion of the test particle s orbit P b 2π/n b : orbital period of the test particle s orbit e : eccentricity of the test particle s orbit I : inclination of the orbital plane of the test particle s orbit to the reference{x, y} plane of some coordinate system Ω : longitude of the ascending node of the test particle s orbit referred to the reference{x, y} plane of some coordinate system Appendix B Tables

14 14 Table 1: Relevant physical and orbital parameters for Saturn, Jupiter, Enceladus and Europa. Most of the reported values come from Soffel et al. 2003; Seidelmann et al. 2007; Petit, Luzum & et al and references therein. The source for the orbital elements referred to either the mean ecliptic ecl at the reference epoch J or the mean Earth s equator eq at the same epoch is the freely consultable database JPL HORIZONS on the Internet at from which they were retrieved by choosing the time of writing this paper as input epoch. Parameter Units Numerical value G kg m 3 s c m s S kg m 2 s α deg δ deg a Enc km 237, 948 e Enc I eq Enc deg I ecl Enc deg Ω eq Enc deg Ω ecl Enc deg P Enc b d S kg m 2 s α deg δ deg a Eur km 671, 034 e Eur I eq Eur deg I ecl Eur deg Ω eq Eur deg Ω ecl Eur deg P Eur b d

15 15 Table 2: Relevant physical and orbital parameters used in the text for the Sun, Mercury and the Earth. Most of the reported values come from Soffel et al. 2003; Seidelmann et al. 2007; Petit, Luzum & et al and references therein. The source for the orbital elements referred to either the mean ecliptic ecl at the reference epoch J or the mean Earth s equator eq at the same epoch is the freely consultable database JPL HORIZONS on the Internet at from which they were retrieved by choosing the time of writing this paper as input epoch. Parameter Units Numerical value G kg m 3 s c m s S kg m 2 s α deg δ deg a au e I eq I ecl Ω eq Ω ecl deg deg deg deg P b yr a au e I eq deg I ecl deg Ω eq deg Ω ecl deg P b yr