UDC 530.1; 539.1
Post-Newtonian effects in the motion of the nearest satellites of Jupiter
T. S. Boronenko
Department of Theoretical Physics, Tomsk State Pedagogical University, Tomsk, 634061, Russia.
E-mail: boron@tspu.edu.ru
In this article the possibility of measurement of post-Newtonian effects in the motion of close satellites of Jupiter is discussed. On an example of Jupiter’s fifth satellite Amalthea we study a question of, whether can be isolated the PN component of orbital precession of the satellite from much bigger Newtonian components. Results of researches have shown that all larger contributions of Newtonian perturbations can be modeled and subtracted out.
Keywords: celestial mechanics, perturbation theory, satellites dynamics, relativistic effects, Jupiter’s satellites.
1 Introduction
It is known that gravitation fields of the gas giant planets are the best laboratories for measuring the effects of Post-Newtonian (PN) gravity on trajectories of natural satellites and spacecrafts [1]. In this work the problem connected with the PN shift of pericenters of orbits of close satellites of Jupiter is considered. Theoretically predicted the PN shift of pericenters of orbits of the close satellites of Jupiter much more, than for the Mercury orbit (b = 43"/per century) [1], [2].
Table 1. PN precession rates b (arcsec/per century)
MetisJ16 AdrasteaJ15 AmaltheaJ5 ThebeJ14
5286".64 5283.65 2212.85 1335.57
Although the post-Newtonian effects are very large for these satellites (Table 1), there are serious problems to measure and separate out such effects. In article J.D. Anderson, et al. "Gravitation and Celestial Mechanics Investigations with Galileo" [3] it is possible to read: "It is unknown whether the relativistic components of orbital precession for the inner satellites can be isolated from the far larger Newtonian precessions". In this work the attempt of studying of this question is undertaken.
2 Perturbing factors in movement of of the Jupiter’s close satellites
The Hamiltonian function of a considered problem within the PN formalism can be written down in a form:
H = Ho + Hn + Hpn , (1)
where H0 is the Hamiltonian function of a separable system (two-body problem), HN is the disturbing Newtonian potential function and HPN is the relativistic disturbing part of H. Function HN is
Here Hoi,Ho2,Ho3 are perturbations caused respectively by the Jupiter’s oblateness, Galilean satellites and the Sun.
The orbit of the fifth satellite of Jupiter Amalthea is the most studied. Therefore we use it as the test.
The first term of Hamiltonian function can be introduced in form:
TT MJ
Ho = ^-, Mj = Gmj, 2a
(3)
where a is the semi-magor axis of an orbit of the satellite, G is the Newtonian gravitational constant, mj is mass of Jupiter. This term establishes the zero order level for the Hamiltonian function.
The most considerable disturbing influence is defined by the second zonal harmonic of the Jupiter’s potential. Considering only a leading factor in expression for a zonal harmonic, we receive
MjRj J2 HorA :
J2 = 1473.6 • 10-
(4)
Here RJ is the equatorial radius of Jupiter, rA is the mean orbit radius of a satellite, J2 is the dimensionless coefficients of the second zonal harmonic. The ratio (4) defines the small parameter e of our problem. The same way the valuations of perturbations from Galilean satellites and from the Sun were made.
The term HPN is treated as the additional potential to the Newtonian disturbing function which describes the geodesic motion of a test particle (satellite) in the gravitational field of a spherically symmetric body (Jupiter) of mass mJ. In the framework of General Relativity (GR) the dominant GR effect can be described by the following expression
TPN
Mh2
(5)
Hn
Roi + H02 + R
¿03.
(2)
For a case which is considered in this work , the right part of (5) is interpreted as follows: ^ = ^J,r =
5
c2r3
TSPU Bulletin. 2012. 13 (128)
rA- h = r29 = \JMa • (1 — e2) is the orbital angular momentum of the satellite per unit mass, 9 is the poe
contributions from Hpn is tpn/R0.
The numerical values of contributions from the disturbing potential terms are given in Table 2:
Table 2. The disturbing function of the problem.
Disturbing action (Leading factor) /H0 Order
Zonal harmonic (J2) 4.55866 • 10-3 e
Galilean satellites) 9.43476 • 10-6 e2
Zonal harmonic (J4) 2.81962 • 10-5 e2
Zonal harmonic(J6) 2.30768 • 10-7 e3
Sun 2.64792 • 10-8 e3
Hpn 1.55476 • 10-8 e3
Analysis of results shows that the influence of the PN effects of the gravity are rather big and it is necessary to take them into consideration in the process modeling of movement of the inner satellites of Jupiter. The contribution from HPN is commensurable with Solar perturbations. And, if it is possible to separate Solar perturbations, it can be possible to separate and PN effects also. Influence of other bodies of the Solar system on movement of the close satellites of Jupiter is negligible.
3 Secular perturbations
Usually the disturbing function is represented in the form of the truncated series at which there are century and periodic terms. Our problem (1) is nonintegrable. However, by means of suitable approximations it is possible to find an analytical solution for (1) with sufficient accuracy. In this work we concentrate only on studying of the secular motions of satellites.
With help of canonical transformations and the averaging method the function H0i ( oblateness part) was presented in a form:
1
(#01) = 2 n2a2
1 2 2 ----n a
15 {R\4 1 T (R
~T J4 ( — ) + 15J6 ( —
4 V a J \ a
e2-
e2-
1 2 2 —n a 2
3 jJ—)2 - 27 я(—'4
2 2 V W 8 2 I a
2
sin2 г +
, 1 2 2 + 2 n a
15 R 4 R 6
~YJ4 ( — ) + 15J6 ( —
4 a a
sin2 г
(6)
Here n = 2n/Ts is the mean motion of a satellite, where Ts is its orbital period, i is the orbital inclination of a satellite to the equatorial plane of Jupiter. Expression (6) by accuracy to J4 can be found in work [4].
Then, using the equations of motion of the satellite in the orbital elements, we received expressions for the
rate of change of pericentre W01 (Newtonian part from oblateness):
W 01 — n
T J‘(R
4a
4
R
+ 15J6 —
(7)
The full velue of the Newtonian shift of a pericenter is determined by a formula:
Wn — Woi + W 02 + Wo3,
(8)
Calculation of perturbations in zu02 (from Galilean satellites) and in b03 (from the Sun)was carried out on known analytical expressions [4].
4 Numerical experiment
Observable value of shift of a pericenter of Amalthea was received by P.V. Sudbury from the analysis of a large number of observations of the satellite [5]
Ws — 3". 30264001 • 108per century.
(9)
At first we received total value of the Newtonian shift of a pericenter without having included the sixth harmonica:
wff — 3". 30426244 • 108 per century.
(10)
The similar calculations made by S. Breiter [6] , give quite close value of the shift of a pericenter:
zbff = 3". 30409477 • 108per century (Breiter). (11)
The difference between the observable value and the calculated value of the shift of a pericenter is:
AW(4) — W s — wff
-162242".
(12)
The difference has negative value, and the problem becomes uncertain.
Then the sixth harmonica was included. The following result was received:
W^ — 3". 30261439 • 108 per century. and
AW(6) — Ws — W^ — 2566". 5.
(13)
(14)
Predicted relativistic shift of pericenter of an orbit of Amalthea has the following value (Table 1):
шpn = 2212". 85per century. (15)
Thus the value Aw(6) is commensurable with a value of the predicted relativistic shift of a pericenter of the satellite.
6
n
a
6
5 Conclusion ter is commensurable with the value of the predicted
relativistic shift of the pericenter of a satellite. The The presented analytical solution for the secular received result says that, apparently, possibility of isoshift of the pericenter of a satellite includes the sixth iation 0f a relativistic part of pertutbations from New-
harmonic of the gravitational potential of Jupiter. The tonian exists. But this problem is very difficult and it
difference between the observable value and the calcu- is necessary to consider the received result preliminary,
lated value of the Newtonian part of shift of a pericen-
References
[1] Hiscock, W. A., Lindblom, L. 1979. The Astrophysical Journal 231, 224-228.
[2] Boronenko T. S. 2012. Tomsk. TSPU Bulletin 7(122), 70-75.
[3] Anderson J. D. et al. 1992. Space Science Reviews 60, 591-610.
[4] Murray C. D., Dermott S. F. 1999. Solar System Dynamics (Cambridge university press) 592 p.
[5] Sudbury P. V. 1968. Icarus 110 , 116-143.
[6] Breiter S. 1996. Astron. Astrophis 314, 966-976.
Received 01.10.2012
Т. С. Вороненко
ПОСТНЬЮТОНОВСКИЕ ЭФФЕКТЫ В ДВИЖЕНИИ БЛИЗКИХ СПУТНИКОВ
ЮПИТЕРА
В статье обсуждается проблема измерения постныотоновских (ПН) эффектов в движении близких спутников Юпитера. На примере пятого спутника Юпитера Амальтеи решается вопрос об отделении ПН компоненты орбитальной прецессии спутника от сравнительно большой по величине ньютоновской части. Результаты исследования показали, что все нерелятивистские возмущения в вековом движении перицентра спутника могут быть получены с достаточно высокой степенью точности, позволяющей отделить их от релятивистской компоненты.
Ключевые слова: небесная механика, теория возмущений, релятивистские эффекты, спутники Юпитера.
Вороненко Т. С., кандидат физико-математических наук, доцент. Томский государственный педагогический университет.
Ул. Киевская, 60, Томск, Россия, 634061.
E-mail: boron@tspu.edu.ru