Kepler's laws

  1. The orbit of a planet is an ellipse with the sun at one of the two foci.

  2. A line segment joining a planet and the sun sweeps out equal areas

during equal time intervals.

  1. The planet’s orbital period TT is related to the semi-major axis of its

orbit aa by T2a3T^2 \propto a^3.

Kepler’s three laws follow from Newton’s law of gravitation: that the force between two bodies of mass m1m_1 and m2m_2 is

FG=Gm1m2r2r^. \begin{align*} \mathbf F_G = -\frac{Gm_1m_2}{r^2} \mathbf{\hat r}. \end{align*}

We consider the motion of the two bodies in the center of mass frame. Let rr be their separation and MM be their total mass. We see that the distances from the COM are rm2/Mrm_2/M and rm1/Mrm_1/M. The Lagrangian is then

L=12m1[(m2r˙M)2+(m2rθ˙M)2]12m2[(m1r˙M)2+(m1rθ˙M)2]+Gm1m2r. \begin{align*} L = \frac12 m_1 \left[ \left(\frac{m_2\dot r}{M}\right)^2 + \left(\frac{m_2r\dot\theta}{M}\right)^2 \right] \frac12 m_2 \left[ \left(\frac{m_1\dot r}{M}\right)^2 + \left(\frac{m_1r\dot\theta}{M}\right)^2 \right] + \frac{Gm_1m_2}{r}. \end{align*}

Following some algebra, we define the reduced mass μ=m1m2/M\mu = m1_m2/M to write

L=12μ(r˙2+r2θ˙2)+Gm1m2r. \begin{align*} L = \frac12 \mu(\dot r^2 + r^2\dot\theta^2) + \frac{Gm_1m_2}{r}. \end{align*}

We see there is no θ\theta dependence so L/θ˙\partial L/\partial\dot\theta is conserved. In this case that is the angular momentum, l=μr2θ˙l = \mu r^2 \dot\theta. This is enough to prove Kepler’s 2nd law.

The area swept out in some slice of the orbit is approximately  ⁣dA1/2(r+ ⁣dr)r ⁣dθ1/2r2 ⁣dθ\d A \approx 1/2\, (r+\d r)r \d\theta \approx 1/2\,r^2 \d\theta.

Dividing by  ⁣dt\dt we get A˙=1/2r2θ˙\dot A = 1/2\,r^2 \dot\theta. We substitute in the angular momentum ll to show A˙=l/2μ\dot A = l / 2\mu is constant.