Kepler’s Laws from Comoment Maps

3.1 The Vector Field

Definition 4 Kepler vector field
#

The Kepler equations of motion are

\[ \dot q = \frac{1}{m}p, \qquad \dot p = -\mu \frac{q}{\| q\| ^3}. \]

The first equation says velocity equals momentum divided by mass. The second is Newton’s law for the inverse-square gravitational force. In Lean, the second component is defined by cases so the expression makes sense on all of \(\mathbb {R}^3 \times \mathbb {R}^3\). The case \(q = 0\) is never reached by a non-collision solution.

Definition 5 Solutions
#

A curve \(\gamma \colon \mathbb {R}\to \mathbb {R}^3 \times \mathbb {R}^3\) is a solution if:

  1. \(\gamma '(t)\) equals the Kepler vector field at \(\gamma (t)\) for every \(t\),

  2. \(\| q(t)\| \neq 0\) for every \(t\) (no collision).

Proposition 6 Position equation
#

For every solution, \(\dot q = p/m\).

Proof

By Definition 5, the derivative \(\gamma '(t)\) agrees with the Kepler vector field at \(\gamma (t)\). Taking the first component of that vector field gives exactly

\[ \dot q(t) = \frac{1}{m}p(t). \]
Proposition 7 Momentum equation

For every solution, \(\dot p = -\mu q/\| q\| ^3\).

Proof

Again \(\gamma '(t)\) equals the Kepler vector field at \(\gamma (t)\), so taking the second component gives the momentum equation. In the formalization the vector field is defined by cases so that it makes sense on all of \(\mathbb {R}^3 \times \mathbb {R}^3\). However, a non-collision solution never reaches \(q = 0\), so along the trajectory we are always in the branch

\[ \dot p(t) = -\mu \frac{q(t)}{\| q(t)\| ^3}. \]
Proposition 8 The Kepler force is radial

For every solution, \(q \times \dot p = 0\).

Proof

By the momentum equation, \(\dot p = (-\mu /\| q\| ^3)q\) is a scalar multiple of \(q\). The cross product of any vector with a scalar multiple of itself is zero:

\[ q \times \dot p = q \times \! \left(-\frac{\mu }{\| q\| ^3}q\right) = -\frac{\mu }{\| q\| ^3}(q \times q) = 0. \]