Confirming Einstein’s Theory of General Relativity With Calculus, Part 7d: Predicting Precession IV

In this series, I’m discussing how ideas from calculus and precalculus (with a touch of differential equations) can predict the precession in Mercury’s orbit and thus confirm Einstein’s theory of general relativity. The origins of this series came from a class project that I assigned to my Differential Equations students maybe 20 years ago.

We have shown that the motion of a planet around the Sun, expressed in polar coordinates $(r,theta)$ with the Sun at the origin, under general relativity is

$u(\theta) \approx \displaystyle \frac{1}{\alpha} \left[ 1 + \epsilon \cos \left( \theta - \frac{\delta \theta}{\alpha} \right) \right]$ ,

where $u = \displaystyle \frac{1}{r}$ , $\alpha = a(1-\epsilon^2)$ , $a$ is the semi-major axis of the planet’s orbit, $\epsilon$ is the orbit’s eccentricity, $\delta = \displaystyle \frac{3GM}{c^2}$ , $G$ is the gravitational constant of the universe, $m$ is the mass of the planet, $M$ is the mass of the Sun, $P$ is the planet’s perihelion, $\ell$ is the constant angular momentum of the planet, and $c$ is the speed of light.

The above function $u(\theta)$ is maximized (i.e., the distance from the Sun $r(\theta)$ is minimized) when $\displaystyle \cos \left( \theta - \frac{\delta \theta}{\alpha} \right)$ is as large as possible. This occurs when $\theta - \displaystyle \frac{\delta \theta}{\alpha}$ is a multiple of $2\pi$ .

Said another way, the planet is at its closest point to the Sun when $\theta = 0$ . One orbit later, the planet returns to its closest point to the Sun when

$\theta - \displaystyle \frac{\delta \theta}{\alpha} = 2\pi$

$\theta \displaystyle\left(1 - \frac{\delta}{\alpha} \right) = 2\pi$

$\theta = 2\pi \displaystyle\frac{1}{1 - (\delta/\alpha)}$

We now use the approximation

$\displaystyle \frac{1}{1-x} \approx 1 + x \qquad \hbox{if} \qquad x \approx 0$ ;

this can be demonstrated by linearization, Taylor series, or using the first two terms of the geometric series $1 + x + x^2 + x^3 + \dots$ . With this approximation, the closest approach to the Sun in the next orbit occurs when

$\theta = 2\pi \displaystyle\left(1 + \frac{\delta}{\alpha} \right) = 2\pi + \frac{2\pi \delta}{\alpha}$ ,

which is coterminal with the angle

$\phi = \displaystyle \frac{2\pi \delta}{\alpha}$ .

Substituting $\alpha = a(1-\epsilon^2)$ and $\delta = \displaystyle \frac{3GM}{c^2}$ , we see that the amount of precession per orbit is

$\phi = \displaystyle 2 \pi \frac{3GM}{c^2} \frac{1}{a(1-\epsilon^2)} = \frac{6\pi G M}{ac^2(1-\epsilon^2)}$ .

The units of $\phi$ are radians per orbit. In the next post, we will use Mercury’s data to find $\phi$ in seconds of arc per century.

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 7d: Predicting Precession IV”

Confirming Einstein’s Theory of General Relativity With Calculus, Part 7d: Predicting Precession IV

Published by John Quintanilla

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 7d: Predicting Precession IV”

Leave a comment Cancel reply

Share this:

Related

Published by John Quintanilla

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 7d: Predicting Precession IV”

Leave a comment Cancel reply