Kepler Orbits and the Runge-Lenz Vector
~6 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)
You should be comfortable with vector calculus identities, specifically the cross product rule for derivatives, and have reviewed the derivation of the orbit equation using Binet's substitution in the previous lecture.
Why this matters
The inverse-square law is unique in classical mechanics because it possesses a hidden symmetry that forces orbits to close perfectly.
While Newton's laws describe the instantaneous force, the Lagrangian and Hamiltonian frameworks reveal deeper structural properties. In the case of the Kepler problem, there exists an additional conserved quantity beyond energy and angular momentum. This vector, known as the Laplace–Runge–Lenz vector, points directly at the perihelion (closest approach) and fixes the orientation of the orbit. Its conservation is the mathematical reason why planetary orbits are stable ellipses rather than precessing rosettes, a feature that distinguishes the inverse-square potential from almost all other central forces.
The Geometry of Closed Orbits
Imagine a particle moving in a plane under a central force. We already know from the previous lecture that the shape of the path is a conic section. However, knowing the shape is not enough; we need to know where the "focus" of that conic section lies relative to the particle's motion. In a generic central force field, the orbit might slowly rotate or "precess," meaning the point of closest approach moves around the center over time. This would look like a flower petal pattern.
For the specific case of gravity or electrostatics, where the force drops off exactly as the inverse square of the distance, something special happens. The orbit does not precess. The ellipse (or hyperbola) stays fixed in space. This implies that there is a vector quantity that always points from the center of force to the perihelion, and its magnitude remains constant. This vector is the Laplace–Runge–Lenz vector. It acts as a "compass needle" for the orbit, locking the major axis in place. This is a rare symmetry in classical mechanics, and it is the key to understanding why Kepler's laws hold so precisely for ideal two-body systems.
Defining the Conserved Vector
Let us define the Laplace–Runge–Lenz vector \(\mathbf{A}\) for a particle of mass \(m\) moving under the potential \(V(r) = -k/r\). The vector is constructed from the linear momentum \(\mathbf{p}\), the angular momentum \(\mathbf{L}\), and the position vector \(\mathbf{r}\):
Here, \(k\) is the force constant (e.g., \(G M m\) for gravity), and \(\hat{\mathbf{r}}\) is the unit vector pointing from the center of force to the particle. Both terms lie in the plane of motion: \(\mathbf{p} \times \mathbf{L}\) is perpendicular to \(\mathbf{L}\), which is itself perpendicular to that plane, and \(\hat{\mathbf{r}}\) lies in it. The combination is carefully chosen so that its time derivative vanishes when the force is inverse-square.
It is crucial to note that \(\mathbf{A}\) is perpendicular to the angular momentum \(\mathbf{L}\). We can verify this by taking the dot product:
Since \(\mathbf{A}\) is perpendicular to \(\mathbf{L}\), it must lie entirely within the plane of the orbit. This confirms that \(\mathbf{A}\) points in a specific direction within the orbital plane, which we will identify as the direction of the perihelion.
Direction and the Orbit Equation
Take the dot product of \(\mathbf{A}\) with \(\mathbf{r}\). By the scalar triple product, \((\mathbf{p} \times \mathbf{L}) \cdot \mathbf{r} = \mathbf{p} \cdot (\mathbf{L} \times \mathbf{r})\), and \(\mathbf{L} \times \mathbf{r} = (\mathbf{r} \times \mathbf{p}) \times \mathbf{r} = \mathbf{p}\, r^2 - \mathbf{r}\,(\mathbf{p} \cdot \mathbf{r})\), so
Write \(\mathbf{A} \cdot \mathbf{r} = |\mathbf{A}|\, r \cos\theta\), with \(\theta\) measured from the direction of \(\mathbf{A}\), and solve for \(r\):
This is the conic of the previous lecture with eccentricity \(e = |\mathbf{A}|/(mk)\), that is \(|\mathbf{A}| = mke\). Because \(r\) is smallest at \(\theta = 0\), \(\mathbf{A}\) points from the centre of force to the perihelion. For a circular orbit \(e = 0\) and \(\mathbf{A} = 0\): a circle has no distinguished perihelion.
Magnitude and Energy
Squaring the definition, with \((\mathbf{p} \times \mathbf{L})^2 = p^2 L^2\) (because \(\mathbf{p} \perp \mathbf{L}\)) and \(\hat{\mathbf{r}} \cdot (\mathbf{p} \times \mathbf{L}) = L^2/r\) from above:
so \(|\mathbf{A}| = mk\sqrt{1 + 2EL^2/(mk^2)} = mke\). Bound orbits (\(E < 0\)) have \(e < 1\) and are ellipses; \(E = 0\) gives the parabola \(e = 1\), and \(E > 0\) the hyperbolae \(e > 1\).
The derivation
Kepler Orbits and the Laplace-Runge-Lenz Vector
For a particle of (reduced) mass m moving under the attractive inverse-square central force F(r) = −k/r2 r̂ (k > 0), the bound and unbound orbits are conic sections r(θ) = p/(1 + e cosθ) with semi-latus rectum p = L2/(mk) and eccentricity e = √(1 + 2EL2/mk2); and the Laplace–Runge–Lenz vector A = p × L − mk r̂ is a constant of the motion whose conservation forces the orbit to close and its major axis to stay fixed (no precession).
The Laplace–Runge–Lenz vector is a constant of motion only for the pure inverse-square potential \(V(r) = -k/r\). If you add any small perturbation, such as a \(1/r^2\) term (from general relativity or quadrupole moments), \(\mathbf{A}\) is no longer conserved. Its direction will slowly rotate, leading to perihelion precession. This is why the conservation of \(\mathbf{A}\) is a diagnostic tool: if you observe precession, you know the force is not exactly inverse-square. C
Worked example
Problem. A particle of mass \(m = 1.0\ \mathrm{kg}\) moves in the attractive potential \(V = -k/r\) with \(k = 4.0\ \mathrm{N\,m^2}\). Its angular momentum is \(L = 2.0\ \mathrm{kg\,m^2/s}\) and its energy \(E = -1.5\ \mathrm{J}\). Find the eccentricity, \(|\mathbf{A}|\), and the perihelion and aphelion distances.
- \(\dfrac{2EL^2}{mk^2} = \dfrac{2(-1.5)(2.0)^2}{(1.0)(4.0)^2} = -0.75\), so \(e = \sqrt{1 - 0.75} = 0.50\): an ellipse.
- \(|\mathbf{A}| = mke = (1.0)(4.0)(0.50) = 2.0\ \mathrm{kg^2\,m^3/s^2}\), pointing at the perihelion.
- Semi-latus rectum \(L^2/(mk) = 4.0/4.0 = 1.0\ \mathrm{m}\); perihelion \(r_{\min} = 1.0/(1 + e) \approx 0.67\ \mathrm{m}\) at \(\theta = 0\); aphelion \(r_{\max} = 1.0/(1 - e) = 2.0\ \mathrm{m}\) at \(\theta = \pi\).
- Check against the energy: at the turning points \(E = V_{\mathrm{eff}} = -k/r + L^2/(2mr^2)\); at \(r = 2.0\ \mathrm{m}\): \(-2.0 + 0.5 = -1.5\ \mathrm{J}\).
Units check. \(EL^2/(mk^2)\): \(\mathrm{J\,(kg\,m^2/s)^2 / (kg\,(N\,m^2)^2)}\) is dimensionless; \(mk\) is in \(\mathrm{kg \cdot kg\,m^3/s^2}\), the same units as \(\mathbf{p} \times \mathbf{L}\); \(L^2/(mk)\) is in metres.
Three quick tests
- Is the Laplace–Runge–Lenz vector conserved for a particle in a harmonic oscillator potential \(V(r) = \frac{1}{2}kr^2\)?
Answer
No. The harmonic oscillator has its own hidden symmetry (a conserved tensor, not a vector), and the Runge–Lenz vector is specific to the \(1/r\) potential. The orbits are ellipses centered on the force center, not focused on it. - If the eccentricity \(e\) is zero, what is the value of the Laplace–Runge–Lenz vector?
Answer
\(\mathbf{A} = 0\). A circular orbit has no unique perihelion, so the vector that points to the perihelion must vanish. - Why does the conservation of \(\mathbf{A}\) imply that the orbit does not precess?
Answer
\(\mathbf{A}\) points to the perihelion. If \(\mathbf{A}\) is constant in both magnitude and direction, the perihelion direction is fixed in space. Therefore, the major axis of the ellipse does not rotate, and the orbit closes perfectly after one revolution.
Problem set for PU-201 →
Each derivation in this unit ends with five graded problems with full solutions.
Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.