Estimating Planetary Perihelion Precession

​
display
big equation
parameters
effective potential
comparison
convergence
planet
☿
♀

♂
♃
♄
⛢
♆
(for potential and convergence)
perihelion precession
(arcseconds per Julian year)
planet
key
black: data
red dashed: Fitzpatrick method
blue: correct to O(1)
purple: correct to O(​
4
α
​)
Making the simplified approximation that planets move in an isotropic central potential and applying conservation of angular momentum, the planetary equations of motion reduce to one-dimensional and become exactly soluble. In this simple picture, the rate of perihelion precession can be calculated to arbitrary precision. Applying methods introduced in[1], we obtain a "big equation" for the precession rate of an orbit's perihelion. Using data from[2] and the pseudopotential formulation of[3, 4], we determine values for all unknowns (see Details).
Substituting all determined coefficients for all planets into the big equation, we finally obtain precession predictions, which closely match those of[3, 4]. From the convergence graphs we see that, except with Mercury, a zero-order or first-order approximation is usually accurate enough. After calculating classical precession, it becomes possible to calculate the tiny additional anomalous precession by applying the same formulas on a gravitational central potential derived from general relativity[5].
Planet symbols: Mercury
☿
, Venus
♀
, Earth

, Mars
♂
, Jupiter
♃
, Saturn
♄
, Uranus
⛢
, Neptune
♆
.

Details

As the computations in this Demonstration show, derivation of the big equation is relatively simple and fast when applied to a Hamiltonian of the form
H
ω
=
1
2

2
p
+
2
q
+
∞
∑
n=3
ϵ
n
n
q
,(1)
where both
2
p
and
2
q
have the dimension of angular momentum. Per our ansatz for the phase space trajectory, the expansion parameter is then
α=
2
β
=
2H
ω
.
All gravitational pseudopotentials need to be transformed to the form of equation (1). The effective Hamiltonian for a particular planet can be expanded around a point
r
0
according to
H
r
=
2
(
p
r
)
2m
+V(
r
0
)+
∞
∑
n=2
n
v
n
(r-
r
0
)
.(2)
Solar system data combined with the pseudopotential formulation[3, 4] uniquely determines the parameters
r
0
and
v
n
. As in[1], we then apply a canonical transformation
r⟶q=
r-
r
0
k
p
r
⟶p=k
p
r
v
n
⟶
ϵ
n
=
n
k
v
n
ω
(3)
k=
-1/4
(2
v
2
m)
ω=
2
v
2
/m
to put the Hamiltonian (2) into the form of (1) by determining
ω
and as many
ϵ
n
as necessary. Notice that transformation (3) always forces
ϵ
2
to equal
1/2
, as in equation (1) and the parameter table. We also transform the equilibrium radius by the canonical transformation
ρ=
r
0
k
=
r
0
.
Finally, we estimate energy
α
from data by using the observed perihelion
r
p
and aphelion values
r
a
to compute a variable
δr=
r
a
-
r
p
for the span between the extrema where
p
r
=0
. Then the equation
δr=k(q(α,ϕ=0)-q(α,ϕ=π))
can be solved numerically for
α
. This completes specification of the data transformation theory, which determines all parameters in the table, thus determining the annual precession rate.
Comparing the calculation of Fitzpatrick to the calculation presented here, the most striking difference is the number and type of required variables. Fitzpatrick requires only the sets of mean radii and planetary masses. In our estimation, these quantities go eventually into the formulation of the pseudopotentials. But we also need angular momentum
L
and energy
α
, which help to determine just how the orbit will deform away from a circular shape.
Taking a look at the comparison graph, we see that the energy-independent part of our approximation
δθ=1-
L
2
ρ
=1-
2
k
L
2
(
r
0
)
=1-
-1/2
(2
v
2
m)
L
2
(
r
0
)
=1-
L
mω
2
(
r
0
)
does an acceptable job of matching the Fitzpatrick values, but the agreement is not exact for all planets. Mercury and Mars are noticeably different.

References

[1] B. Klee, "Plane Pendulum and Beyond by Phase Space Geometry." arxiv.org/abs/1605.09102.
[2] NASA. "Planetary Fact Sheet." (Nov 18, 2015) nssdc.gsfc.nasa.gov/planetary/factsheet.
[3] R. Fitzpatrick, An Introduction to Celestial Mechanics, New York: Cambridge University Press, 2012.
[4] R. Fitzpatrick, "Perihelion Precession of the Planets," Newtonian Dynamics, 2011 (Sep 28, 2016). farside.ph.utexas.edu/teaching/336k/Newtonhtml/node115.html.
[5] A. Einstein, "Erklärung der Perihelionbewegung der Merkur aus der allgemeinen Relativitätstheorie." Sitzungsberichte der Königlich Preußischen Akademie der Wissenschaften, 1915 pp. 831–839. (Sep 27, 2016) einsteinpapers.press.princeton.edu/vol6-doc/261. English translation: einsteinpapers.press.princeton.edu/vol6-trans/128?ajax.

External Links

Approximating the Jacobian Elliptic Functions
Exact and Approximate Relativistic Corrections to the Orbital Precession of Mercury

Permanent Citation

Brad Klee
​
​"Estimating Planetary Perihelion Precession"​
​http://demonstrations.wolfram.com/EstimatingPlanetaryPerihelionPrecession/​
​Wolfram Demonstrations Project​
​Published: September 30, 2016