Next Article in Journal / Special Issue
The Gravitational Magnetoelectric Effect
Previous Article in Journal
Progress in Solving the Nonperturbative Renormalization Group for Tensorial Group Field Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Post-Newtonian Gravitomagnetic Effect on the Orbital Motion of a Test Particle around Its Primary Induced by the Spin of a Distant Third Body

Ministero dell’Istruzione, dell’Università e della Ricerca (M.I.U.R.)-Istruzione, Viale Unità di Italia 68, 70125 Bari (BA), Italy
Universe 2019, 5(4), 87; https://doi.org/10.3390/universe5040087
Submission received: 4 March 2019 / Revised: 20 March 2019 / Accepted: 28 March 2019 / Published: 31 March 2019
(This article belongs to the Special Issue Universe: Feature Papers 2019 - Gravitational Physics)

Abstract

:
We study a general relativistic gravitomagnetic 3-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 system S of size r r X consisting of a test particle revolving around a massive body M. At the lowest post-Newtonian order, 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 ω . Such theoretical results are confirmed by a numerical integration of the equations of motion for a fictitious 3-body system. We numerically calculate the magnitudes of the post-Newtonian gravitomagnetic 3-body precessions 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 ≃10 − 50 milliarcseconds per year (mas year−1). A preliminary numerical simulation shows that, for certain orbital configurations of a hypothetical Europa orbiter, its range-rate signal Δ ρ ˙ can become larger than the current Doppler accuracy of the existing spacecraft Juno at Jupiter, i.e., σ ρ ˙ = 0.015 mm s−1, after 1 d. The effects induced by the Sun’s angular momentum on artificial probes of Mercury and the Earth are at the level of ≃1 − 0.1 microarcseconds per year (μas year−1).

1. Introduction

Let us consider a local gravitationally bound restricted two-body system S composed by a test particle completing a full orbital revolution around a planet of mass M at distance r in a time interval P b , and a distant 3rd body X with mass M X M and proper spin S X around which S revolves at distance r X r with orbital period P b X . In general, M may be endowed with its own Newtonian and post-Newtonian (pN) mass and spin multipole moments [1,2] affecting the satellite’s motion with known [3,4] and less known [5,6,7] Newtonian and pN orbital effects like the classical oblateness-driven orbital precessions, the gravitoelectric Einstein pericentre shift, the gravitomagnetic Lense-Thirring effect, etc. Let us consider a kinematically rotating and dynamically non-rotating coordinate system K [8,9,10] 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 system K X whose axes point towards the distant quasars [8,9,10]. The planetocentric motion of the test particle referred to K is further affected by two peculiar pN 3-body effects: the time-honored De Sitter precession due to solely the mass M X [11,12,13], and a gravitomagnetic shift due to S X which, to our knowledge, has never been explicitly and clearly calculated in the literature, if it had ever been. Our purpose is to analytically work out the latter effect at the lowest pN order without any a-priori simplifying assumptions concerning both the orbital configurations of the planetocentric satellite’s motion and the trajectory of the planet-satellite system S in the external field of X, and for an arbitrary orientation of S X in space. For the 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. [14] (Section 3.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 4. For the benefit of the reader, Appendix A contains a list of the definitions of the symbols used in the paper, while their numerical values and tables are collected in Appendix B.

2. The Doubly Averaged Satellite’s Orbital Precessions

In the weak-field and slow-motion approximation of general relativity, the gravitomagnetic 3-body potential induced by the angular momentum S X of the external spinning object X on the planetary satellite is
U GM = G c 2 r X 3 S X · L + 3 L · r ^ X r ^ X .
In Equation (1), G , c are the Newtonian constant of gravitation and the speed of light in vacuum, respectively, while L being the angular momentum per unit mass of the test particle’s orbital motion around M. Equation (1) comes from Equation (2.19) of Barker & O’Connell [15] (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 with 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 Equation (2.19) of Barker & O’Connell [15] (p. 155) has to be identified with r X , and r × P is nothing but the orbital angular momentum r X × M v X of the motion of S around M X . It is interesting to note that, with the same identifications, V S 1 and V S 2 of Equations (2.17) and (2.18) in Barker & O’Connell [15] (p. 155) yield the gravitoelectric De Sitter orbital precession for the planetocentric motion of the satellite and the gravitomagnetic Lense-Thirring effect for the X-centric orbit of M, respectively.
The velocity-dependent perturbing potential U pert to be inserted into the Lagrange planetary equations [10,16] for the rates of change of the nonosculating [10,17] Keplerian orbital elements of the test particle, obtained by doubly averaging Equation (1) with respect to P b , P b X 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 = G S X n b a 2 1 e 2 2 c 2 a X 3 1 e X 2 3 / 2 U ,
with
U = cos I 2 S ^ z 3 sin I X S ^ z sin I X + cos I X S ^ y cos Ω X S ^ x sin Ω X + + sin I 2 2 S ^ y cos Ω 2 S ^ x sin Ω + 3 cos Ω Ω X S ^ z sin 2 I X + + 2 sin 2 I X S ^ y cos Ω X + S ^ x sin Ω X .
In Equations (2) and (3), a , e , I , a X , e X , I X are the semimajor axes, the eccentricities and the inclinations of the orbits of the test particle and of S , respectively, while n b is the Keplerian orbital motion of the satellite’s planetary motion about M. Equations (2) and (3) were obtained in two steps. First, U GM of Equation (1) was evaluated onto the unperturbed ellipse of the planetocentric satellite motion through the standard Keplerian formulas of the restricted two-body problem (see, e.g., Equations (3.40a) to (3.41c) of Poisson & Will [4]). Then, it was averaged over one orbital period P b to the first order in the disturbing potential by using just the Keplerian part of Equation (3.66) of Poisson & Will [4] for d f / d t , where f is the true anomaly. Then, the resulting averaged potential U ¯ GM was, in turn, calculated onto the unperturbed Keplerian trajectory of S about X and averaged over P b X to the first order in the perturbation under consideration, thus finally obtaining the double average of Equations (2) and (3).
Inserting Equations (2) and (3) into the right-hand-sides of the Lagrange planetary equations allows to calculate the doubly averaged rates of change of the Keplerian orbital elements. They turn out to be
a ˙ = 0 ,
e ˙ = 0 ,
I ˙ = G S X 2 a X 3 c 2 1 e X 2 3 / 2 I ,
Ω ˙ = G S X 2 a X 3 c 2 1 e X 2 3 / 2 O ,
ω ˙ = G S X csc I 8 a X 3 c 2 1 e X 2 3 / 2 P .
with
I = sin Ω S ^ y + 3 sin I X cos Ω X S ^ z cos I X + sin I X S ^ y cos Ω X S ^ x sin Ω X + + cos Ω S ^ x + 3 sin I X sin Ω X S ^ z cos I X + sin I X S ^ y cos Ω X + S ^ x sin Ω X ,
O = 2 S ^ z + S ^ x cot I sin Ω 3 cos I X sin I X S ^ y cos Ω X S ^ x sin Ω X + S ^ z cot I sin Ω sin Ω X 3 sin 2 I X S ^ z + cot I sin Ω sin Ω X S ^ y cos Ω X + S ^ x sin Ω X + + cos Ω cot I S ^ y + 3 sin I X cos Ω X S ^ z cos I X + sin I X S ^ y cos Ω X S ^ x sin Ω X ,
P = S ^ y cos Ω 3 cos Ω 2 Ω X S ^ x sin Ω + 3 sin Ω 2 Ω X   + +   6 cos Ω Ω X S ^ z sin 2 I X + cos 2 I X S ^ y cos Ω X S ^ x sin Ω X .
We remark that Equations (4) and (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) and (11) consists, first of all, in deriving a perturbing acceleration from Equation (1). By writing the Lagrangian per unit mass of a gravitationally bound restricted two-body system affected by a generic perturbing potential as
L = L 0 + L pert = v 2 2 + μ r + L pert ,
the conjugate momentum per unit mass is, by definition,
p L v = v + L pert v .
Thus,
p ˙ = v ˙ + d d t L pert v .
The Hamiltonian per unit mass is
H = H 0 + H pert = v 2 2 μ r + H pert .
From the Hamilton equations of motion, it is
p ˙ = H r = μ r 3 r H pert r .
Since L pert = H pert , by comparing Equations (14) and (16), it turns out that the perturbing acceleration is just
A pert = d d t H pert v H pert r .
In our specific case, since H pert = U GM , we have
A GM = d d t U GM v U GM r = 2 G c 2 r X 3 v × S X 3 S X · r ^ X r ^ X .
Then, Equation (18) must be decomposed into its radial ( ρ ), transverse ( τ ) and out-of-plane ( ν ) components, which are
A ρ GM = G S X c 2 r X 5 S ^ x r X 2 + 3 S ^ x x X 2 + 3 S ^ y x X y X + 3 S ^ z x X z X cos Ω + + S ^ y r X 2 + 3 S ^ x x X y X + 3 S ^ y y X 2 + 3 S ^ z y X z X sin Ω ,
A τ GM = G S X csc I c 2 r X 5 sin I 3 S ^ x x X + S ^ y y X z X + S ^ z r X 2 3 z X 2 + + cos I S ^ y r X 2 3 y X 2 3 S ^ x x X + S ^ z z X y X cos Ω + + cos I S ^ x r X 2 + 3 x X 2 + 3 S ^ y x X y X + 3 S ^ z x X z X sin Ω ,
A ν GM = G S X csc I c 2 r X 5 S ^ y r X 2 3 y X 2 3 y X S ^ x x X + S ^ z z X cos Ω + + S ^ x r X 2 + 3 x X 2 + 3 S ^ y x X y X + 3 S ^ z x X z X sin Ω .
They have to be inserted into the right-hand-sides of the standard Gauss equations for the variation of the orbital elements [4,10,16] which, finally, are doubly averaged with respect to P b , P b X in the same way as previously described for the disturbing potential of Equation (1).
Even putting aside the post-Keplerian effects of classical and pN nature depending on the specific characteristics of the field of M itself, the motion of the test particle is perturbed also by the Newtonian 3-body acceleration due to the gravitational pull of X. By coupling with Equation (18), it would give rise to mixed orbital perturbations which can be worked out, in principle, as in Iorio [18], Will [19]. A rigorous calculation of such effects is outside the scope of the present paper. An approximate evaluation of their order of magnitude can be made as follows. For the sake of simplicity, let the motion of S occur in a plane perpendicular to S X , assumed as reference x , y plane so that I X = 0 , S ^ X · r ^ X = 0 , S ^ x = S ^ y = 0 , S ^ z = 1 , along a circular orbit with r / r X 1 . Thus, the Newtonian 3-body acceleration felt by the test particle, obtainable from the tidal-type potential by Hogg, Quinlan & Tremaine [20], can be approximately posed equal to A X G M X / r X 2 r ^ X by neglecting the other two terms proportional to r / r X r ^ , r / r X r ^ · r ^ X r ^ X . Then, following Iorio [18], it can be shown that the resulting mixed precessions are proportional to ξ X v X / c 2 M X / M R X / r X 2 r / r X 2 Ψ X , where v X G M X / r X is the orbital velocity of S about X, and ξ X < 1 , R X , Ψ X are the normalized moment of inertia, the equatorial radius and the angular speed of X, respectively. This implies that the mixed effects are smaller than Equations (4)–(8) by a scaling factor of the order of M X / M r / r X 2 . In the case of, say, Jupiter and Europa, by assuming r approximately equal to the moon’s radius R, it amounts to ≃0.2.
We successfully checked our analytical results of Equations (4)–(11) as follows. We considered a fictitious system S orbiting a Jupiter-like body X along the same orbit of, say, Callisto, whose mass was assumed for the particle’s primary M, and numerically integrated its equations of motion over a time span much longer than P b , P b X with and without the pN gravitomagnetic acceleration of Equation (18) affecting the test particle; both the integrations, which assumed a purely Keplerian motion of S about the fictitious body X, shared the same initial conditions for the test particle and its primary. For X, the same physical properties of Jupiter were assumed, including the size and the orientation of its angular momentum S . As a result, numerically produced times series of the orbital elements of the putative probe were produced by subtracting the purely classical ones from those obtained by including also Equation (18) in the equations of motion; they are displayed in Figure 1. It turned out that the resulting numerically calculated pN gravitomagnetic 3-body orbital shifts agree with those computed by means of the analytical formulas of Equations (4)–(11).

3. Some Potentially Interesting Astronomical Scenarios

In the case of a hypothetical orbiter of the Kronian natural satellite Enceladus in the external field of Saturn, Equations (6), (7), (9) and (10), referred to the mean Earth’s equator at the reference epoch J2000.0 as reference x , y plane, yield for a circular orbit ( e = 0 ), assumed just for the sake of simplicity,
I ˙ = A eq sin Ω + φ eq ,
Ω ˙ = 49.9 mas year 1 + cot I A eq cos Ω + φ eq ,
with
A eq = 5.7 mas year 1 ,
φ eq = 49.4 deg .
Instead, if the mean ecliptic at the reference epoch J2000.0 is adopted as reference x , y plane, we have
I ˙ = A ecl sin Ω + φ ecl ,
Ω ˙ = 34.0 mas year 1 + cot I A ecl cos Ω + φ ecl ,
with
A ecl = 23.9 mas year 1 ,
φ ecl = 10.4 deg .
By looking at a putative orbiter of the Jovian natural satellite Europa in the external field of Jupiter, we have
I ˙ = A eq sin Ω + φ eq ,
Ω ˙ = 9.9 mas year 1 + cot I A eq cos Ω + φ eq ,
with
A eq = 4.8 mas year 1 ,
φ eq = 2.9 deg ,
and
I ˙ = A ecl sin Ω + φ ecl ,
Ω ˙ = 11.0 mas year 1 + cot I A ecl cos Ω + φ ecl ,
with
A ecl = 0.3 mas year 1 ,
φ ecl = 31.0 deg .
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 [21]. As a consequence, they are the natural targets of several concept studies and proposals for dedicated missions to them, including also orbiters [22,23,24,25,26]. Since, at present, sending a spacecraft to Europa seems more likely than to Enceladus, as it can be learnt at https://europa.nasa.gov/about-clipper/overview/ and http://sci.esa.int/juice/ on the Internet, we investigated in a little more detail this potentially appealing Jovian scenario, even if it is not said that the actually approved missions will finally involve the use of an orbiter. In such kind of endeavours, the observable quantity is typically the Earth-probe range-rate ρ ˙ , whose accuracy for, e.g., the ongoing mission Juno [27] around Jupiter is σ ρ ˙ 0.015 mm s 1 [28]. Figure 2 shows the numerically simulated Earth-spacecraft range-rate signature due to the pN gravitomagnetic 3-body acceleration of Equation (18) for a generic orbital configuration of the hypothesized orbiter. To produce it, we numerically integrated the equations of motion in Cartesian rectangular coordinates of the Earth, Jupiter, its Galilean moons and a fictitious test particle orbiting Europa over 1 d. In both runs, sharing the same initial conditions retrieved from the database JPL HORIZONS (https://ssd.jpl.nasa.gov/?horizons) at the arbitrary epoch of midnight on 1 January 2030, we modeled the mutual attractions among all the bodies involved to the Newtonian level, with the exception of Equation (18) which was added to the other classical gravitational pulls felt by the probe in one of the runs. Then, we numerically calculated two range-rate time series, with and without Equation (18) ceteris paribus, and subtracted the purely Newtonian one from that including also the pN gravitomagnetic acceleration. It can be noted that, for the orbital configuration chosen, the range-rate relativistic signature Δ ρ ˙ reaches the 0.05 mm s 1 level after just 1 d. Thus, the scenario considered seems worthy of further, dedicated analyses investigating the actual measurability of Equation (18) in a realistic error budget analysis and mission proposal. It should take into account several concurring perturbations of gravitational and non-gravitational nature, and also several technological and engineering issues.
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 in a circular orbit, we have
I ˙ = A eq sin Ω + φ eq ,
Ω ˙ = 0.2 μ as year 1 + cot I A eq cos Ω + φ eq ,
with
A eq = 0.1 μ as year 1 ,
φ eq = 9.13 deg ,
and
I ˙ = A ecl sin Ω + φ ecl ,
Ω ˙ = 0.3 μ as year 1 + cot I A ecl cos Ω + φ ecl ,
with
A ecl = 0.02 μ as year 1 ,
φ ecl = 104.2 deg .
For a probe orbiting Mercury with e = 0 , one gets
I ˙ = A eq sin Ω + φ eq ,
Ω ˙ = 4.3 μ as year 1 + cot I A eq cos Ω + φ eq ,
with
A eq = 2.5 μ as year 1 ,
φ eq = 171.3 deg ,
and
I ˙ = A ecl sin Ω + φ ecl ,
Ω ˙ = 5 μ as year 1 + cot I A ecl cos Ω + φ ecl ,
with
A ecl = 0.6 μ as year 1 ,
φ ecl = 144.6 deg .

4. Summary and Overview

In the weak-field and slow-motion approximation of general relativity, we analytically worked out the pN 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 mass-energy currents of the spin angular momentum S X of a distant ( r X r ) 3rd body X with mass M X M . We did not assume any preferred orientation for the spin axis S ^ 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 has a general validity, being applicable to arbitrary astronomical systems of potential interest. It turns out that, by doubly averaging the perturbing potential employed in the calculation with respect to the orbital periods P b , P b X 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 I ˙ 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 with the frequency of the node itself. A numerical integration of the equations of motion of a fictitious 3-body system made of a distant spinning body with the same physical properties of Jupiter, a primary with the same orbital and physical characteristics of Callisto and a test particle orbiting it confirms our analytical results.
The Sun’s angular momentum exerts very small effects on spacecraft orbiting Mercury ( 1 μ as year 1 ) and the Earth ( 0.1 μ as year 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 year 1 ) and Enceladus (≲ 50 mas year 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. In the case of Europa, for whose exploration there are already approved missions by NASA and ESA, a preliminary numerical simulation of the signature induced by the pN gravitomagnetic 3-body effect of interest on the range-rate of a putative orbiter shows that, for certain orbital configurations, its magnitude can become larger than the present-day accuracy σ ρ ˙ = 0.015 mm s 1 of the current Juno mission around Jupiter after 1 d.

Acknowledgments

I thank M. Efroimsky for his useful remarks which contributed to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

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 Table A1 and Table A2.
  • 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 G M X : gravitational parameter of the 3rd body X
  • R X : equatorial radius of the 3rd body X
  • S X : magnitude of the angular momentum of the 3rd body X
  • ξ X : normalized moment of inertia of the 3rd body X
  • Ψ X : angular speed of the 3rd body X
  • S ^ X = S ^ x , S ^ y , S ^ 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
  • S ^ x eq = cos δ X cos α X : component of the 3rd body’s spin axis w.r.t. the reference x axis of an equatorial coordinate system
  • S ^ y eq = cos δ X sin α X : component of the 3rd body’s spin axis w.r.t. the reference y axis of an equatorial coordinate system
  • S ^ z eq = 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 of S to 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 b X μ X / a X 3 : mean motion of the orbit about the 3rd body X
  • P b X 2 π / n b X : orbital period of the orbit about the 3rd body X
  • e X : eccentricity of the orbit about the 3rd body X
  • 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
  • μ G M : gravitational parameter of the primary orbited by the test particle and moving in the external field of the 3rd body X
  • R : radius of the primary (planet or planetary natural satellite) orbited by the test particle and moving in the external field of the 3rd body X
  • S : angular momentum of the primary
  • r : position vector of the test particle with 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
  • f : true anomaly 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

Table A1. Relevant physical and orbital parameters for Saturn, Jupiter, Enceladus and Europa. Most of the reported values come from Petit, Luzum & et al. [29], Seidelmann et al. [30], Soffel et al. [31] and references therein. The source for the orbital elements referred to either the mean ecliptic (ecl) at the reference epoch J2000.0 or the mean Earth’s equator (eq) at the same epoch is the freely consultable database JPL HORIZONS on the Internet at https://ssd.jpl.nasa.gov/?horizons from which they were retrieved by choosing the time of writing this paper as input epoch.
Table A1. Relevant physical and orbital parameters for Saturn, Jupiter, Enceladus and Europa. Most of the reported values come from Petit, Luzum & et al. [29], Seidelmann et al. [30], Soffel et al. [31] and references therein. The source for the orbital elements referred to either the mean ecliptic (ecl) at the reference epoch J2000.0 or the mean Earth’s equator (eq) at the same epoch is the freely consultable database JPL HORIZONS on the Internet at https://ssd.jpl.nasa.gov/?horizons from which they were retrieved by choosing the time of writing this paper as input epoch.
ParameterUnitsNumerical Value
G kg 1 kgm 3 s 2 6.67259 × 10 11
c m s 1 2.99792458 × 10 8
S kg m 2 s 1 1.4 × 10 38
α deg 40.59
δ deg 83.54
a Enc km 237,948
e Enc 0.0047
I Enc eq deg 6.475336858877378
I Enc ecl deg 28.06170970578348
Ω Enc eq deg 130.5900992493321
Ω Enc ecl deg 169.5108697290241
P b Enc d 1.370218
S kg m 2 s 1 6.9 × 10 38
α deg 268.05
δ deg 64.49
a Eur km 671,034
e Eur 0.0094
I Eur eq deg 25.88280598312641
I Eur ecl deg 1.790876103183550
Ω Eur eq deg 357.4169659423443
Ω Eur ecl deg 332.6268549691798
P b Eur d 3.551810
Table A2. Relevant physical and orbital parameters used in the text for the Sun, Mercury and the Earth. Most of the reported values come from Petit, Luzum & et al. [29], Seidelmann et al. [30], Soffel et al. [31] and references therein. The source for the orbital elements referred to either the mean ecliptic (ecl) at the reference epoch J2000.0 or the mean Earth’s equator (eq) at the same epoch is the freely consultable database JPL HORIZONS on the Internet at https://ssd.jpl.nasa.gov/?horizons from which they were retrieved by choosing the time of writing this paper as input epoch.
Table A2. Relevant physical and orbital parameters used in the text for the Sun, Mercury and the Earth. Most of the reported values come from Petit, Luzum & et al. [29], Seidelmann et al. [30], Soffel et al. [31] and references therein. The source for the orbital elements referred to either the mean ecliptic (ecl) at the reference epoch J2000.0 or the mean Earth’s equator (eq) at the same epoch is the freely consultable database JPL HORIZONS on the Internet at https://ssd.jpl.nasa.gov/?horizons from which they were retrieved by choosing the time of writing this paper as input epoch.
ParameterUnitsNumerical Value
G kg 1 m 3 s 2 6.67259 × 10 11
c m s 1 2.99792458 × 10 8
S kg m 2 s 1 1.90 × 10 41
α deg 286.13
δ deg 63.87
a au 0.3870982252717257
e 0.2056302512089075
I eq deg 28.55225598038233
I ecl deg 7.005014199657344
Ω eq deg 10.98794759075666
Ω ecl deg 48.33053756455964
P b yr 0.2408467
a au 0.9992521882390240
e 0.01731885059206812
I eq deg 23.43903457134406
I ecl deg 2.669113820737183 × 10 4
Ω eq deg 1.852352676284691 × 10 4
Ω ecl deg 163.9752443600624
P b yr 1.0000174

References

  1. Frutos-Alfaro, F.; Soffel, M. On relativistic multipole moments of stationary space-times. R. Soc. Open Sci. 2018, 5, 180640. [Google Scholar] [CrossRef] [PubMed]
  2. Soffel, M.; Frutos, F. On the Usefulness of Relativistic Space-Times for the Description of the Earth’s Gravitational Field. J. Geod. 2016, 90, 1345–1357. [Google Scholar] [CrossRef]
  3. Capderou, M. Satellites: Orbits and Missions; Springer: Berlin, Germany, 2005. [Google Scholar]
  4. Poisson, E.; Will, C.M. Gravity; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  5. Angélil, R.; Saha, P.; Bondarescu, R.; Jetzer, P.; Schärer, A.; Lundgren, A. Spacecraft clocks and relativity: Prospects for future satellite missions. Phys. Rev. D 2014, 89, 064067. [Google Scholar] [CrossRef]
  6. Schanner, M.; Soffel, M. Relativistic satellite orbits: Central body with higher zonal harmonics. Celest. Mech. Dyn. Astr. 2018, 130, 40. [Google Scholar] [CrossRef]
  7. Schärer, A.; Bondarescu, R.; Saha, P.; Angélil, R.; Helled, R.; Jetzer, P. Prospects for Measuring Planetary Spin and Frame-Dragging in Spacecraft Timing Signals. Front. Astron. Space Sci. 2017, 4, 11. [Google Scholar] [CrossRef]
  8. Brumberg, V.A.; Kopeikin, S.M. Relativistic reference systems and motion of test bodies in the vicinity of the Earth. Nuovo Cimento B 1989, 103, 63–98. [Google Scholar] [CrossRef]
  9. Damour, T.; Soffel, M.; Xu, C. General-relativistic celestial mechanics. IV. Theory of satellite motion. Phys. Rev. D 1994, 49, 618. [Google Scholar] [CrossRef]
  10. Kopeikin, S.; Efroimsky, M.; Kaplan, G. Relativistic Celestial Mechanics of the Solar System; Wiley-VCH: Weinheim, Germany, 2011. [Google Scholar]
  11. de Sitter, W. Einstein’s theory of gravitation and its astronomical consequences, Second Paper. Mon. Not. R. Astron. Soc. 1916, 77, 155–184. [Google Scholar] [CrossRef]
  12. Fokker, A.D. De geodetische precessie; een uitvloeisel van Einstein’s gravitatie- theorie. Versl. Kon. Ak. Wet. 1920, 29, 611–621. [Google Scholar]
  13. Schouten, W.J.A. Over het ontstaan eener praecessiebeweging tengevolge van het niet euklidisch zijn der ruimte in de nabijheiud van zon. Versl. Kon. Ak. Wet. 1918, 27, 214. [Google Scholar]
  14. Gill, E.; Soffel, M.; Ruder, H.; Schneider, M. Relativistic Motion of Gyroscopes and Space Gradiometry; Deutsche Geodätische Kommission: München, Germany, 1992. [Google Scholar]
  15. Barker, B.M.; O’Connell, R.F. The gravitational interaction: Spin, rotation, and quantum effects-A review. Gen. Relat. Gravit. 1979, 11, 149–175. [Google Scholar] [CrossRef]
  16. Bertotti, B.; Farinella, P.; Vokrouhlický, D. Physics of the Solar System; Kluwer Academic Press: Dordrecht, The Netherlands, 2003. [Google Scholar]
  17. Efroimsky, M. Gauge freedom in orbital mechanics. Ann. N. Y. Acad. Sci. 2005, 1065, 346–374. [Google Scholar] [CrossRef]
  18. Iorio, L. Post-Newtonian direct and mixed orbital effects due to the oblateness of the central body. Int. J. Mod. Phys. D 2015, 24, 1550067. [Google Scholar] [CrossRef]
  19. Will, C.M. Incorporating post-Newtonian effects in N-body dynamics. Phys. Rev. D 2014, 89, 044043. [Google Scholar] [CrossRef]
  20. Hogg, D.W.; Quinlan, G.D.; Tremaine, S. Dynamical limits on dark mass in the outer solar system. Astron. J. 1991, 101, 2274–2286. [Google Scholar] [CrossRef]
  21. Lunine, J.I. Ocean worlds exploration. Acta Astronaut. 2017, 131, 123–130. [Google Scholar] [CrossRef]
  22. MacKenzie, S.M.; Caswell, T.E.; Phillips-Lander, C.M.; Stavros, E.N.; Hofgartner, J.D.; Sun, V.Z.; Powell, K.E.; Steuer, C.J.; O’Rourke, J.G.; Dhaliwal, J.K.; et al. THEO concept mission: Testing the habitability of Enceladus’s Ocean. Adv. Space Res. 2016, 58, 1117–1137. [Google Scholar] [CrossRef]
  23. Razzaghi, A.I.; di Pietro, D.A.; Quinn, D.A.; Simon-Miller, A.A.; Tompkins, S.D. Space Technology and Applications International Forum-STAIF 2008; American Institute of Physics Conference Series; El-Genk, M.S., Ed.; American Institute of Physics: College Park, MD, USA, 2008; Volume 969, pp. 388–395. [Google Scholar]
  24. Sherwood, B.; Lunine, J.; Sotin, C.; Cwik, T.; Naderi, F. Program options to explore ocean worlds. Acta Astronaut. 2018, 143, 285–296. [Google Scholar] [CrossRef]
  25. Spencer, J.; Niebur, C. Planetary Science Decadal Survey. Enceladus Orbiter; National Aeronautics and Space Administration: Washington, DC, USA, 2010.
  26. Verma, A.K.; Margot, J.-L. Expected precision of Europa Clipper gravity measurements. Icarus 2018, 314, 35–49. [Google Scholar] [CrossRef]
  27. Bolton, S.J.; Lunine, J.; Stevenson, D.; Connerney, J.E.P.; Levin, S.; Owen, T.C.; Bagenal, F.; Gautier, D.; Ingersoll, A.P.; Orton, G.S.; et al. The Juno mission. Space Sci. Rev. 2017, 213, 5–37. [Google Scholar] [CrossRef]
  28. Iess, L.; Folkner, W.M.; Durante, D.; Parisi, M.; Kaspi, Y.; Galanti, E.; Guillot, T.; Hubbard, W.B.; Stevenson, D.J.; Anderson, J.D.; et al. Measurement of Jupiter’s asymmetric gravity field. Nature 2018, 555, 220–222. [Google Scholar] [CrossRef]
  29. Petit, G.; Luzum, B. IERS Conventions. IERS Tech. Note 2010, 36, 1. [Google Scholar]
  30. Seidelmann, P.K.; Archinal, B.A.; A’hearn, M.F.; Conrad, A.; Consolmagno, G.J.; Hestroffer, D.; Hilton, J.L.; Krasinsky, G.A.; Neumann, G.; Oberst, J.; et al. Report of the IAU/IAG Working Group on cartographic coordinates and rotational elements: 2006. Celest. Mech. Dyn. Astr. 2007, 98, 155–180. [Google Scholar] [CrossRef]
  31. Soffel, M.; Klioner, S.A.; Petit, G.; Wolf, P.; Kopeikin, S.M.; Bretagnon, P.; Brumberg, V.A.; Capitaine, N.; Damour, T.; Fukushima, T.; et al. The IAU 2000 resolutions for astrometry, celestial mechanics, and metrology in the relativistic framework: Explanatory supplement. Astron. J. 2003, 126, 2687. [Google Scholar] [CrossRef]
Figure 1. Numerical time series of the pN 3-body shifts of I , Ω , ω of a fictitious test particle moving around a Callisto-like primary M which orbits a Jupiter-type 3rd body X. They were obtained by integrating the equations of motion of the orbiter about M and of M about X in Cartesian rectangular coordinates referred to the Earth’s mean equator at the epoch J2000.0 with and without Equation (18) acting on the test particle. Both runs shared the same set of arbitrary initial conditions for the probe P b = 10.07 d , e 0 = 0.3 , I 0 = 80 deg , Ω 0 = 230 deg , ω 0 = 40 deg , f 0 = 50 deg and the primary; the initial state vector of the Callisto-Jupiter relative motion was adopted from the database JPL HORIZONS. For each Keplerian orbital element, its time series calculated from the purely Newtonian run was subtracted from that obtained from the pN integration in order to obtain the signatures displayed here. The resulting rates, in mas year 1 , agree with those computed in Equations (6)–(8).
Figure 1. Numerical time series of the pN 3-body shifts of I , Ω , ω of a fictitious test particle moving around a Callisto-like primary M which orbits a Jupiter-type 3rd body X. They were obtained by integrating the equations of motion of the orbiter about M and of M about X in Cartesian rectangular coordinates referred to the Earth’s mean equator at the epoch J2000.0 with and without Equation (18) acting on the test particle. Both runs shared the same set of arbitrary initial conditions for the probe P b = 10.07 d , e 0 = 0.3 , I 0 = 80 deg , Ω 0 = 230 deg , ω 0 = 40 deg , f 0 = 50 deg and the primary; the initial state vector of the Callisto-Jupiter relative motion was adopted from the database JPL HORIZONS. For each Keplerian orbital element, its time series calculated from the purely Newtonian run was subtracted from that obtained from the pN integration in order to obtain the signatures displayed here. The resulting rates, in mas year 1 , agree with those computed in Equations (6)–(8).
Universe 05 00087 g001
Figure 2. Numerically produced Earth-probe range-rate shift Δ ρ ˙ t due to the pN gravitomagnetic 3-body acceleration of Equation (18). We numerically integrated the solar system barycentric equations of motion in Cartesian rectangular coordinates of the Earth, Jupiter, its Galilean moons and a fictitious test particle orbiting Europa over 1 d. In both runs, sharing the same initial conditions for all the existing natural bodies retrieved from the database JPL HORIZONS (https://ssd.jpl.nasa.gov/?horizons) at the arbitrary epoch of midnight of 1 January 2030, we modeled the mutual attractions among all the planets and the satellites involved to the Newtonian level, with the exception of Equation (18) which was added to the other classical gravitational pulls felt by the orbiter in one of the runs. Then, we numerically calculated two range-rate time series, and subtracted the purely Newtonian one from that including also Equation (18). The orbital configuration adopted for the spacecraft, referred to Europa, was a 0 = 3.55 R , e 0 = 0.69 , I 0 = 100 deg , Ω 0 = 90 deg , ω 0 = 40 deg , f 0 = 50 deg , where R is the radius of Jovian moon.
Figure 2. Numerically produced Earth-probe range-rate shift Δ ρ ˙ t due to the pN gravitomagnetic 3-body acceleration of Equation (18). We numerically integrated the solar system barycentric equations of motion in Cartesian rectangular coordinates of the Earth, Jupiter, its Galilean moons and a fictitious test particle orbiting Europa over 1 d. In both runs, sharing the same initial conditions for all the existing natural bodies retrieved from the database JPL HORIZONS (https://ssd.jpl.nasa.gov/?horizons) at the arbitrary epoch of midnight of 1 January 2030, we modeled the mutual attractions among all the planets and the satellites involved to the Newtonian level, with the exception of Equation (18) which was added to the other classical gravitational pulls felt by the orbiter in one of the runs. Then, we numerically calculated two range-rate time series, and subtracted the purely Newtonian one from that including also Equation (18). The orbital configuration adopted for the spacecraft, referred to Europa, was a 0 = 3.55 R , e 0 = 0.69 , I 0 = 100 deg , Ω 0 = 90 deg , ω 0 = 40 deg , f 0 = 50 deg , where R is the radius of Jovian moon.
Universe 05 00087 g002

Share and Cite

MDPI and ACS Style

Iorio, L. A Post-Newtonian Gravitomagnetic Effect on the Orbital Motion of a Test Particle around Its Primary Induced by the Spin of a Distant Third Body. Universe 2019, 5, 87. https://doi.org/10.3390/universe5040087

AMA Style

Iorio L. A Post-Newtonian Gravitomagnetic Effect on the Orbital Motion of a Test Particle around Its Primary Induced by the Spin of a Distant Third Body. Universe. 2019; 5(4):87. https://doi.org/10.3390/universe5040087

Chicago/Turabian Style

Iorio, Lorenzo. 2019. "A Post-Newtonian Gravitomagnetic Effect on the Orbital Motion of a Test Particle around Its Primary Induced by the Spin of a Distant Third Body" Universe 5, no. 4: 87. https://doi.org/10.3390/universe5040087

APA Style

Iorio, L. (2019). A Post-Newtonian Gravitomagnetic Effect on the Orbital Motion of a Test Particle around Its Primary Induced by the Spin of a Distant Third Body. Universe, 5(4), 87. https://doi.org/10.3390/universe5040087

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop