09.01.04 · classical-mech / newtonian

Two-body central-force problem, Kepler orbits, and Rutherford scattering

shipped3 tiersLean: none

Anchor (Master): Goldstein, Poole & Safko, *Classical Mechanics*, 3e (Addison-Wesley, 2002), §3 (full Kepler treatment, Laplace-Runge-Lenz vector); Landau & Lifshitz, *Mechanics*, 3e (Pergamon, 1976), §11-19; Arnold, *Mathematical Methods of Classical Mechanics*, 2e (Springer GTM 60, 1989), §8

Intuition Beginner

When two objects pull on each other only along the line that joins them — like the Sun and a planet bound by gravity, or an electron and a proton bound by electrical attraction — the motion has a special structure. The pull depends only on how far apart they are, never on the direction the line is pointing. Forces of this kind are called central forces, and the resulting motion is one of the oldest solved problems in physics. Kepler worked out the planetary version from naked-eye data between 1605 and 1619. Newton derived it from his law of gravity in 1687. Rutherford repeated the analysis for repulsive electrical forces in 1911 and discovered the atomic nucleus.

The first simplification is that a two-body problem is really a one-body problem in disguise. The centre of mass of the pair drifts in a straight line at constant speed (nothing pushes the pair from outside), so it can be ignored. What matters is the relative position — how the planet sits with respect to the Sun, or the electron with respect to the proton. That relative position behaves as if a single particle of a certain "reduced mass" were orbiting a fixed centre.

The second simplification is that the orbit lies in a plane. Because the force always points along the line between the two bodies, it can never twist the orbit out of the plane it started in. The plane is fixed by the initial velocity and the initial separation. Inside that plane, the orbit is a curve of one of three shapes: an ellipse if the particle is bound and cannot escape, a parabola if it has the exact escape energy, or a hyperbola if it has more than enough energy to leave forever.

Kepler's three laws are the empirical face of this geometry. First law (1609): planetary orbits are ellipses with the Sun at one focus. Second law (1609): the line from the Sun to the planet sweeps equal areas in equal times — the planet moves faster when it is closer to the Sun. Third law (1619): the square of the orbital period is proportional to the cube of the semi-major axis, with the same constant for every planet in the solar system. Earth: semi-major axis 1 AU, period 1 year. Halley's comet: semi-major axis 17.8 AU, period 75.3 years. The constant matches.

Rutherford's scattering experiment is the same problem in a different costume. Fire an alpha particle (positive charge $+ 2 e$ , four times the proton mass) at a thin gold foil. Most pass through with only a small deflection, but a few bounce nearly straight back. The repulsive electrical force between the alpha particle and the gold nucleus follows the same inverse-square law as gravity, only with the opposite sign. The resulting orbit is a hyperbola: the particle comes in along one branch and leaves along the other. From the rate of back-scattering, Rutherford inferred that the gold atom's positive charge is concentrated in a tiny nucleus, not spread through the atom as J.J. Thomson's plum-pudding model had assumed.

The takeaway: one piece of mathematics — the inverse-square central-force orbit — covers the motion of the planets, the bouncing of charged particles off nuclei, and (once the rules of quantum mechanics are imposed on top) the energy levels of the hydrogen atom.

Visual Beginner

Figure: Three orbits sharing the same focus (the Sun, marked with a yellow dot at the origin). An ellipse (blue) for a bound planet, a parabola (green) for an object with exactly the escape energy, and a hyperbola (red) for an unbound flyby trajectory. The ellipse has its second focus marked with a small grey dot; the line from the Sun to a moving planet sweeps equal grey-shaded areas in equal times, illustrating Kepler's second law.

The same figure works for Rutherford scattering: the red hyperbolic branch becomes the trajectory of an alpha particle being deflected by a gold nucleus at the focus. The angle between the asymptotes is the scattering angle.

Worked example Beginner

A satellite of mass 1000 kg is placed in a circular orbit 400 km above Earth's surface (the altitude of the International Space Station). Use Earth's mass $M = 5.972 \times 1 0^{24}$ kg, Earth's radius $R = 6371$ km, and the gravitational constant $G = 6.674 \times 1 0^{- 11}$ N m²/kg². Find the orbital speed and the period.

Step 1. Orbital radius from the centre of the Earth: $r = R + h = 6371 + 400 = 6771$ km $= 6.771 \times 1 0^{6}$ m.

Step 2. For a circular orbit, the gravitational pull provides the centripetal force needed to bend the satellite around the curve. The speed comes from $GM m / r^{2} = m v^{2} / r$ , which rearranges to $v = GM / r$ .

Step 3. Plug in numbers: $GM = (6.674 \times 1 0^{- 11}) \times (5.972 \times 1 0^{24}) = 3.986 \times 1 0^{14}$ m³/s². Then $v = 3.986 \times 1 0^{14} /6.771 \times 1 0^{6} = 5.887 \times 1 0^{7} = 7672$ m/s, or about 7.67 km/s.

Step 4. Period $T = 2 π r / v = 2 π (6.771 \times 1 0^{6}) /7672 = 5547$ s, or about 92.4 minutes.

Step 5. Escape velocity at the same radius is $v_{esc} = 2 GM / r = 2 \times 7672 = 10848$ m/s, about 10.85 km/s. From Earth's surface the corresponding number is $2 GM / R = 11.18$ km/s.

What this tells us: the ISS travels at almost 28000 km/h and circles the Earth in just over an hour and a half. The orbital speed depends only on the radius and the central mass, not on the satellite's mass — a tennis ball and a space station at the same altitude have the same orbital speed.

Check your understanding Beginner

Exercise (easy, multiple choice).

Kepler's second law states that a planet sweeps out equal areas in equal times. Which conserved quantity is this statement equivalent to?

A. Total energy of the orbit. B. Angular momentum about the Sun. C. Linear momentum of the planet. D. Eccentricity of the orbit.

Hint

The area swept per unit time is half the magnitude of the cross product of the radius vector with the velocity, times the mass. What quantity is that?

Answer

B — Angular momentum about the Sun.

The area swept per unit time is $d A / d t = \frac{1}{2} ∣ r \times v ∣$ , which equals $L / (2 m)$ where $L$ is the magnitude of the angular momentum. Because the gravitational force is central (always parallel to $r$ ), it exerts no torque about the Sun, so $L$ is constant — and therefore the area-sweep rate is constant. Energy is also conserved, but it is not the quantity that sets the area rate. Feedback for option A: energy conservation is a separate law, leading to the shape (conic section) of the orbit but not directly to the area rate. Feedback for D: eccentricity is conserved for a Kepler orbit, but it does not control the area rate.

Exercise (easy, true/false).

In a Rutherford scattering experiment, an alpha particle aimed straight at a gold nucleus (zero impact parameter) is scattered through 180° — it bounces straight back.

Hint

Think about what happens for a head-on approach: where does the alpha particle stop, and which way does it go after it stops?

Answer

True.

For a head-on collision (impact parameter zero), the alpha particle decelerates along the line of approach, comes to a momentary stop at the distance where its kinetic energy equals the Coulomb potential energy, then is pushed back along the same line. The scattering angle is 180°. Larger impact parameters give smaller deflections. The smooth dependence of scattering angle on impact parameter is what Rutherford used to deduce the structure of the atom.

Formal definition Intermediate+

Let two point particles of masses $m_{1}$ and $m_{2}$ at positions $r_{1}, r_{2} \in R^{3}$ interact via a central potential $V (r)$ where $r = ∣ r_{1} - r_{2} ∣$ . The Lagrangian is

L = \frac{1}{2} m_{1} ∣ \dot{r}_{1} ∣^{2} + \frac{1}{2} m_{2} ∣ \dot{r}_{2} ∣^{2} - V (r) .

Centre-of-mass / relative coordinates. Set $M = m_{1} + m_{2}$ , $μ = m_{1} m_{2} / M$ (the reduced mass), $R = (m_{1} r_{1} + m_{2} r_{2}) / M$ (the centre of mass), and $r = r_{1} - r_{2}$ (the relative coordinate). The Lagrangian separates:

L = \frac{1}{2} M ∣ \dot{R} ∣^{2} + \frac{1}{2} μ ∣ \dot{r} ∣^{2} - V (r) .

The Euler-Lagrange equation for $R$ gives $M \ddot{R} = 0$ — the centre of mass drifts in a straight line. The remaining problem is a single particle of mass $μ$ in the central potential $V (r)$ , with Lagrangian

L_{rel} = \frac{1}{2} μ ∣ \dot{r} ∣^{2} - V (r) .

Conserved quantities. Angular momentum $L = μ r \times \dot{r}$ is constant because the force $- \nabla V (r) = - V^{'} (r) \hat{r}$ is parallel to $r$ . The motion is therefore confined to the plane perpendicular to $L$ . In polar coordinates $(r, θ)$ on that plane,

L = μ r^{2} \dot{θ}, E = \frac{1}{2} μ \overset{r}{˙}^{2} + \frac{L ^{2}}{2 μ r ^{2}} + V (r),

with the second term in $E$ the centrifugal barrier. The function $V_{eff} (r) = L^{2} / (2 μ r^{2}) + V (r)$ is the effective potential controlling the one-dimensional radial motion.

Orbit equation. Substitute $u = 1/ r$ and use $\dot{θ} = L / (μ r^{2})$ to eliminate time. A short computation gives the Binet equation:

\frac{d ^{2} u}{d θ ^{2}} + u = - \frac{μ}{L ^{2}} \frac{d V}{d u}_{r = 1/ u} .

For the Kepler potential $V (r) = - k / r$ with $k > 0$ , $d V / d u = - k$ , so the equation reduces to $u^{''} + u = μ k / L^{2}$ , a linear inhomogeneous second-order ODE with the general solution

u (θ) = \frac{μ k}{L ^{2}} (1 + e cos (θ - θ_{0})) .

The integration constants are the eccentricity $e$ and the orientation angle $θ_{0}$ (orient axes so $θ_{0} = 0$ ). Reading off, the orbit $r (θ) = p / (1 + e cos θ)$ is a conic section with semi-latus rectum $p = L^{2} / (μ k)$ and focus at the origin.

Counterexamples to common slips Intermediate+

"The orbit equation is exact for any central potential." The closed-form $u (θ) = (μ k / L^{2}) (1 + e cos θ)$ is special to the inverse-square law (and to the harmonic oscillator $V = \frac{1}{2} k r^{2}$ , by Bertrand's theorem below). A potential like $V = - k / r^{1.001}$ produces precessing rosettes, not closed conics.
"Reduced mass is always close to the smaller mass." Only when one body dominates. For two equal masses, $μ = m /2$ , not $m$ . The earth-moon system has $μ = 0.987 m_{moon}$ since $m_{earth} \approx 81 m_{moon}$ ; the proton-electron system has $μ = 0.99946 m_{e}$ . For a binary star with equal components, $μ = m /2$ , and the period formula changes accordingly.
"Kepler's third law $T^{2} \propto a^{3}$ has a universal constant." The constant is $4 π^{2} / (G M_{total})$ , which depends on the combined mass of the two bodies. Solar-system planets share a common constant because $M_{total} \approx M_{⊙}$ for each one. For a binary star or a binary asteroid the constant is different.

Key theorem with proof Intermediate+

Theorem (Conic-section orbits and Kepler's three laws). Consider a particle of mass $μ$ in the inverse-square central potential $V (r) = - k / r$ with $k > 0$ . Let $E$ and $L$ denote the conserved energy and angular-momentum magnitude. Then:

(i) The orbit lies in a fixed plane and has the form $r (θ) = p / (1 + e cos θ)$ with $p = L^{2} / (μ k)$ and $e = 1 + 2 E L^{2} / (μ k^{2})$ — a conic with focus at the origin.

(ii) Kepler's first law (1609): bound orbits ( $E < 0$ , equivalently $0 \leq e < 1$ ) are ellipses with semi-major axis $a = - k / (2 E)$ . Parabolic orbits have $E = 0$ ( $e = 1$ ); hyperbolic orbits have $E > 0$ ( $e > 1$ ).

(iii) Kepler's second law (1609): the radius vector sweeps area at the constant rate $d A / d t = L / (2 μ)$ .

(iv) Kepler's third law (1619): for bound orbits the period is $T = 2 π μ a^{3} / k$ , equivalently $T^{2} = (4 π^{2} μ / k) a^{3}$ . For the gravitational case $k = G m_{1} m_{2}$ and $μ \approx m_{2}$ when $m_{1} ≫ m_{2}$ , this becomes $T^{2} = 4 π^{2} a^{3} / (G m_{1})$ .

Proof. The conservation of $L$ from the central-force condition was established in 09.01.03; it confines the motion to a fixed plane. In polar coordinates on that plane, $L = μ r^{2} \dot{θ}$ is constant and the area element is $d A = \frac{1}{2} r^{2} d θ$ , so $d A / d t = \frac{1}{2} r^{2} \dot{θ} = L / (2 μ)$ , proving (iii) without any property of $V$ beyond centrality.

For (i), substitute $u = 1/ r$ . Then $\overset{r}{˙} = - L u^{'} / μ$ where $u^{'} = d u / d θ$ , and the energy becomes

E = \frac{L ^{2}}{2 μ} ((u^{'})^{2} + u^{2}) - k u .

Differentiating with respect to $θ$ and dividing by $u^{'}$ gives $u^{''} + u = μ k / L^{2}$ . The general solution is $u (θ) = μ k / L^{2} + C cos (θ - θ_{0})$ . Setting $p = L^{2} / (μ k)$ and $e = C p$ and orienting axes so $θ_{0} = 0$ , the orbit is $r = p / (1 + e cos θ)$ .

The eccentricity formula comes from substituting back into the energy equation. At $θ = 0$ (periapsis, the point of closest approach) $r_{min} = p / (1 + e)$ and $\overset{r}{˙} = 0$ , so $E = L^{2} / (2 μ r_{min}^{2}) - k / r_{min}$ . A short algebraic rearrangement yields $e^{2} = 1 + 2 E L^{2} / (μ k^{2})$ .

The geometry of conic sections classifies the orbits by $e$ : $e = 0$ is a circle, $0 < e < 1$ an ellipse, $e = 1$ a parabola, $e > 1$ a hyperbola. The sign of $E$ correlates: $E < 0$ gives $e < 1$ (bound), $E = 0$ gives $e = 1$ (marginal), $E > 0$ gives $e > 1$ (unbound). For ellipses the major and minor radii are $r_{min} = p / (1 + e)$ and $r_{max} = p / (1 - e)$ , with semi-major axis $a = (r_{min} + r_{max}) /2 = p / (1 - e^{2}) = - k / (2 E)$ , completing (ii).

For (iv), integrate Kepler's second law over a full revolution. The total area of the ellipse is $π ab$ where $b = a 1 - e^{2} = a p$ is the semi-minor axis. Then $T = π ab / (L / (2 μ)) = 2 π μ ab / L = 2 π μ a^{3} p / (μ k /1) \cdot 1/ L$ . Substituting $p = L^{2} / (μ k)$ collapses this to $T = 2 π μ a^{3} / k$ . Squaring gives $T^{2} = (4 π^{2} μ / k) a^{3}$ .

$□$

Bridge. The three laws of Kepler, which appeared in 1609 and 1619 as inductive summaries of Tycho Brahe's planetary observations, are recovered here as direct consequences of two conservation laws (energy, angular momentum) and the inverse-square character of the gravitational force; this is exactly the unification that Newton accomplished in the Principia. The conic-section solution builds toward 09.02.01 (action principle, where the orbit equation re-emerges from a variational integral) and appears again in 09.06.01 (action-angle variables, where the Kepler problem is the canonical integrable system and its action variables generate the periodic motion). The bridge is between the empirical inverse-square law and the geometric closure of bounded orbits, with the foundational reason being the algebraic miracle that the Binet equation linearises for $V \propto 1/ r$ — a property shared only with the harmonic oscillator and codified in Bertrand's theorem below. The result generalises in two directions: relativistically (perihelion precession from a $1/ r^{3}$ correction to $V_{eff}$ , appearing in 13.05.03) and quantum-mechanically (discrete bound-state spectra of the Coulomb potential in 12.06.01).

Exercises Intermediate+

Exercise 3 (medium, symbolic).

Derive the radius of the stable circular orbit in the Kepler potential $V (r) = - k / r$ in terms of the angular momentum $L$ and reduced mass $μ$ . Then express the orbital speed in that circular orbit in terms of $k$ , $μ$ , and the radius $r_{0}$ .

Hint

A circular orbit corresponds to the minimum of the effective potential $V_{eff} (r) = L^{2} / (2 μ r^{2}) - k / r$ .

Answer

Setting $V_{eff}^{'} (r_{0}) = - L^{2} / (μ r_{0}^{3}) + k / r_{0}^{2} = 0$ gives $r_{0} = L^{2} / (μ k)$ — which is also the semi-latus rectum $p$ for $e = 0$ . The orbital speed satisfies $μ v^{2} / r_{0} = k / r_{0}^{2}$ (centripetal equation), so $v = k / (μ r_{0})$ . For the Kepler-gravity case $k = GM μ$ with $M$ the central mass, this reduces to $v = GM / r_{0}$ , the familiar formula. Stability follows from $V_{eff}^{''} (r_{0}) = 3 L^{2} / (μ r_{0}^{4}) - 2 k / r_{0}^{3} = k / r_{0}^{3} > 0$ .

Exercise 6 (medium, symbolic).

Derive the Rutherford differential cross section. Consider a particle of charge $z e$ and energy $E$ scattering off a fixed charge $Z e$ via the repulsive Coulomb potential $V (r) = k / r$ with $k = z Z e^{2} / (4 π ϵ_{0}) > 0$ . Show that the impact parameter $b$ and the scattering angle $θ$ are related by $b = (k / (2 E)) cot (θ /2)$ , and conclude that

\frac{d σ}{d Ω} = (\frac{k}{4 E})^{2} \frac{1}{sin ^{4} ( θ /2 )} .

Hint

The hyperbolic-orbit angle between the incoming asymptote and the perihelion radius is $α$ with $cos α = 1/ e$ . The scattering angle is $θ = π - 2 α$ , equivalently $sin (θ /2) = 1/ e$ . Use the eccentricity formula $e = 1 + (2 E b / k)^{2}$ for a repulsive Coulomb orbit with $E > 0$ .

Answer

For a hyperbolic orbit the asymptotic angle satisfies $cos α = - 1/ e$ for an attractive force (Kepler) or $cos α = 1/ e$ for a repulsive force (Rutherford), with the scattering angle $θ = π - 2 α$ . Then $sin (θ /2) = cos α = 1/ e$ , so $cot (θ /2) = e^{2} - 1$ . The eccentricity in terms of $E$ and the impact parameter is $e^{2} = 1 + 2 E L^{2} / (μ k^{2}) = 1 + (2 E b / k)^{2}$ using $L = μ v_{\infty} b$ and $E = \frac{1}{2} μ v_{\infty}^{2}$ . Thus $cot (θ /2) = 2 E b / k$ , giving $b = (k / (2 E)) cot (θ /2)$ .

The differential cross section is $d σ / d Ω = (b / sin θ) ∣ d b / d θ ∣$ (the Jacobian from impact parameter to scattering angle in axial symmetry). Differentiating $b (θ)$ gives $d b / d θ = - (k / (4 E)) csc^{2} (θ /2)$ . Putting it together: $d σ / d Ω = b \cdot (k / (4 E)) csc^{2} (θ /2) / sin θ$ . Using $sin θ = 2 sin (θ /2) cos (θ /2)$ and $b = (k / (2 E)) cos (θ /2) / sin (θ /2)$ , the result simplifies to $d σ / d Ω = (k / (4 E))^{2} / sin^{4} (θ /2)$ — the Rutherford formula of 1911.

Exercise 7 (hard, symbolic).

The Laplace-Runge-Lenz vector for the Kepler potential $V = - k / r$ is $A = p \times L - μ k \hat{r}$ . Show that $\dot{A} = 0$ along orbits, and prove $A \cdot L = 0$ and $∣ A ∣ = μ k e$ .

Hint

Use Newton's equation $\dot{p} = - k \hat{r} / r^{2}$ and the identity $\hat{r} \times L = μ r (\dot{r} - \hat{r} \overset{r}{˙})$ . For the orientation $A \cdot L$ , exploit that $L$ is perpendicular to the orbital plane while $A$ lies in it.

Answer

Since $\dot{L} = 0$ , $\dot{A} = \dot{p} \times L - μ k d \hat{r} / d t$ . Use $\dot{p} = - k \hat{r} / r^{2}$ and $L = μ r \times \dot{r}$ to expand the first term: $\dot{p} \times L = (- k / r^{2}) \hat{r} \times (μ r \times \dot{r}) = - μ k / r^{2} [r (\hat{r} \cdot \dot{r}) - \dot{r} (\hat{r} \cdot r)]$ via $a \times (b \times c) = b (a \cdot c) - c (a \cdot b)$ . This collapses to $- μ k [\hat{r} \overset{r}{˙} / r - \dot{r} / r] = μ k (\dot{r} / r - \hat{r} \overset{r}{˙} / r) = μ k d \hat{r} / d t$ (since $d \hat{r} / d t = (\dot{r} r - r \overset{r}{˙}) / r^{2} = \dot{r} / r - \hat{r} \overset{r}{˙} / r$ ). Thus $\dot{A} = μ k d \hat{r} / d t - μ k d \hat{r} / d t = 0$ .

For the perpendicularity: $L \cdot (p \times L) = p \cdot (L \times L) = 0$ , and $L \cdot \hat{r} = 0$ because $L$ is perpendicular to the orbital plane while $\hat{r}$ lies in it. Hence $A \cdot L = 0$ .

For the magnitude, evaluate at perihelion ( $r = a (1 - e)$ , $p$ tangential, $p \times L$ radial outward): $∣ p \times L ∣ = p_{peri} L = (μ v_{peri}) L$ . Using $v_{peri} = k (1 + e) / (μ a (1 - e))$ from Exercise 5 and $L = μ k a (1 - e^{2})$ , one gets $∣ p \times L ∣ = μ k (1 + e)$ . The radial $- μ k \hat{r}$ contributes $- μ k$ in the same direction, leaving $∣ A ∣ = μ k (1 + e) - μ k = μ k e$ .

Exercise 8 (hard, symbolic).

State the Kepler equation relating the eccentric anomaly $E$ to the mean anomaly $M$ for an elliptical orbit, and explain how it parameterises the time-dependence of the orbital position. Sketch one Newton-iteration step for solving it numerically.

Hint

Parameterise the elliptic orbit by the eccentric anomaly $E$ via $r = a (1 - e cos E)$ . Integrate Kepler's area-sweep rate to relate $E$ and $t$ .

Answer

The eccentric anomaly $E$ parameterises position on the elliptic orbit via $r = a (1 - e cos E)$ , with the true anomaly $ν$ (angle from the focus) related by $tan (ν /2) = (1 + e) / (1 - e) tan (E /2)$ . Integrating Kepler's areal rule from perihelion (where $E = 0$ and $t = 0$ ): $\frac{1}{2} (L / μ) t = \int_{0}^{E} \frac{1}{2} r^{2} d ν$ . Carrying out the integral (using $r^{2} d ν = a^{2} 1 - e^{2} (1 - e cos E) d E$ after parametric substitution) yields

M = E - e sin E, M := \frac{2 π}{T} t .

This is Kepler's equation. Given the mean anomaly $M$ (linear in time), the transcendental equation $E - e sin E = M$ determines $E$ , and the orbit's position follows from $r$ and $ν$ . Newton-Raphson iteration: $E_{n + 1} = E_{n} - (E_{n} - e sin E_{n} - M) / (1 - e cos E_{n})$ . Start with $E_{0} = M$ ; convergence is quadratic. This algorithm is the basis of every modern ephemeris computation, from satellite ground tracks to NASA's Deep Space Network pointing.

Lean formalization Intermediate+

Mathlib has the topology + ODE infrastructure (Picard-Lindelöf existence, smoothness of integral curves) that would support a Lean proof of the conic-section orbit theorem, but it does not yet contain the Kepler-problem framework. A formalisation would need (i) a CentralForceProblem structure carrying the potential $V : R_{> 0} \to R$ , the reduced mass $μ$ , and the conserved energy / angular-momentum maps; (ii) the Binet substitution as a Lean tactic relating $u$ and $θ$ ; (iii) a proof that the linear ODE $u^{''} + u = μ k / L^{2}$ has the closed-form solution given above; (iv) the Laplace-Runge-Lenz vector with its conservation lemma. Status: none.

Bertrand's theorem and superintegrability of Kepler and harmonic oscillator Master

Theorem 1 (Bertrand 1873). Among all attractive central potentials $V (r)$ for which the equation of motion admits at least one stable circular orbit, the only ones for which every bound orbit is closed are the Kepler potential $V (r) = - k / r$ and the isotropic three-dimensional harmonic oscillator $V (r) = \frac{1}{2} k r^{2}$ . Every other potential produces precessing (non-closed) rosette orbits for generic initial conditions.

The proof, given in Bertrand's three-page note in the Comptes Rendus ^{[Bertrand 1873 *Comptes Rendus* 77, 849-853]}, proceeds by perturbing a circular orbit and computing the apsidal angle (angle from perihelion to aphelion). For a stable circular orbit at radius $r_{0}$ in a potential $V (r)$ , the apsidal angle for small radial oscillations is $Δ θ = π V^{'} (r_{0}) / (r_{0} V^{''} (r_{0}) + 3 V^{'} (r_{0}))$ , which evaluates to $π$ for $V \propto r^{- 1}$ (Kepler) and $π /2$ for $V \propto r^{2}$ (oscillator). For closure of all bounded orbits, the apsidal angle must be a rational multiple of $π$ for every value of the energy and angular momentum; expanding to higher order in the perturbation, this rationality constraint forces $V \propto r^{- 1}$ or $V \propto r^{2}$ . The two exceptional potentials are precisely the ones admitting an additional vector-valued conserved quantity beyond energy and angular momentum: the Laplace-Runge-Lenz vector $A$ for Kepler, the Fradkin (or Demkov) symmetric tensor $T_{ij} = p_{i} p_{j} / (2 μ) + (k /2) r_{i} r_{j}$ for the oscillator. Bertrand's theorem can be regarded as the converse to the orbit-closure half of the Kepler-problem analysis: closure is not generic in central-force dynamics but is a delicate symmetry phenomenon, occurring precisely when the equations of motion possess a Lie algebra of conserved quantities larger than $R \oplus so (3)$ (the energy plus angular-momentum generators).

Theorem 2 (LRL conservation and the SO(4) algebra). For the bound Kepler problem ( $E < 0$ ), define $A = p \times L - μ k \hat{r}$ and the rescaled vector $D = A / - 2 μ E$ . Then the six quantities $L$ and $D$ generate the Lie algebra $so (4)$ under the Poisson bracket:

{L_{i}, L_{j}} = ϵ_{ij k} L_{k}, {L_{i}, D_{j}} = ϵ_{ij k} D_{k}, {D_{i}, D_{j}} = ϵ_{ij k} L_{k} .

For unbound orbits ( $E > 0$ ), the algebra becomes $so (3, 1)$ (using $A / 2 μ E$ ); for $E = 0$ it contracts to $iso (3)$ , the Euclidean algebra in three dimensions.

This is the algebraic content of the celebrated "hidden symmetry" of the Kepler problem. The closure of orbits is the dynamical face of the fact that the symmetry group acts transitively on the level surfaces of the Hamiltonian within each energy stratum. Bargmann 1936 ^{[Bargmann 1936 *Z. Phys.* 99, 576-582]} gave the modern statement of $so (4)$ as the bound-state symmetry; Pauli 1926 ^{[Pauli 1926 *Z. Phys.* 36, 336-363]} had already used the quantum version of this algebra — six months before Schrödinger's equation appeared — to derive the entire non-relativistic hydrogen spectrum $E_{n} = - m_{e} e^{4} / (2 ℏ^{2} n^{2} (4 π ϵ_{0})^{2})$ purely from the commutation relations. The $(n + 1)^{2}$ -fold degeneracy of the hydrogen levels (states with different orbital angular momentum but the same principal quantum number share an energy) is the representation-theoretic statement that $so (4) ≅ so (3) \oplus so (3)$ acts on hydrogen states with both copies carrying the same spin $j = (n - 1) /2$ , giving dimension $(2 j + 1)^{2} = n^{2}$ .

Theorem 3 (Kepler equation). Let $E$ denote the eccentric anomaly (a parametric angle on the auxiliary circle of radius $a$ circumscribing the ellipse) and $M = (2 π / T) t$ the mean anomaly (linear in time). They satisfy the transcendental relation $M = E - e sin E$ .

This equation, present in Kepler's Astronomia Nova ^{[Kepler 1609]} in geometric form, is the workhorse of every ephemeris computation: given the orbital elements $(a, e, ω, Ω, i, T_{peri})$ and an observation epoch $t$ , one computes $M$ from elapsed time, solves the transcendental equation for $E$ by Newton-Raphson iteration (quadratic convergence in $\leq 5$ iterations to double precision for any orbit), and reads off the orbit position from $r = a (1 - e cos E)$ and $ν$ via $tan (ν /2) = (1 + e) / (1 - e) tan (E /2)$ . Modern variants (Markley 1995, Mikkola series-expansions for high eccentricity) reach 12 decimal places in two iterations and are used by JPL HORIZONS, NASA Deep Space Network tracking, and every spacecraft trajectory simulator.

A concrete iteration trace illustrates the convergence rate. For Mars's orbit ( $e = 0.0934$ ) at $M = 1.0$ radians, Newton-Raphson with seed $E_{0} = M = 1.0$ gives $E_{1} = 1.0 - (1.0 - 0.0934 sin 1.0 - 1.0) / (1 - 0.0934 cos 1.0) = 1.0 + 0.0786/ (0.9495) = 1.0828$ , then $E_{2} = 1.0828 - (1.0828 - 0.0934 sin 1.0828 - 1.0) / (1 - 0.0934 cos 1.0828) = 1.0827$ , with the third iterate matching to twelve decimals. For Halley's comet ( $e = 0.967$ ) the same seed converges in $\sim 5$ iterations; for orbits with $e > 0.99$ (long-period comets, hyperbolic comets) the seed $E_{0} = M$ can land in a region where the iteration cycles, and the standard fix is to seed with $E_{0} = π$ if $M > π - e$ or to use a Markley-style higher-order opener. Brouwer 1959's Methods of Celestial Mechanics treats the convergence theory in detail; the orbital-mechanics community uses Battin's 1999 Introduction to the Mathematics and Methods of Astrodynamics as the canonical operational reference.

Theorem 4 (Restricted three-body problem and Lagrange points). In the planar circular restricted three-body problem (two primaries on a circular orbit, test mass much lighter than either), the rotating-frame Hamiltonian admits five equilibrium points: three collinear (L1, L2, L3, unstable saddle points discovered by Euler 1767) and two equilateral (L4, L5, stable for mass ratio $m_{2} / m_{1} < 0.0385$ , discovered by Lagrange 1772).

The collinear points sit on the line through the two primaries; for the Sun-Earth system, L1 (sunward, 1.5 million km in) hosts the SOHO and DSCOVR space observatories, and L2 (anti-sunward, 1.5 million km out) hosts JWST and Euclid. L4 and L5, leading and trailing Jupiter by 60° on its orbit, contain the Trojan asteroids — over 12 000 catalogued as of 2024. The stability condition $m_{2} / m_{1} < (1/2) (1 - 1 - 27/25) \approx 0.0385$ is satisfied by Jupiter ( $m_{J} / m_{⊙} \approx 1 0^{- 3}$ ), Earth, and Neptune, all of which host known Trojans. The Hill region $∣ H - H_{L} ∣ < 0$ around each primary defines its gravitational sphere of influence; the Hill radius $r_{H} \approx a (m /3 M)^{1/3}$ for a small body of mass $m$ orbiting a primary of mass $M$ at semi-major axis $a$ sets the size of its retinue (e.g., the Earth's Hill radius is 1.5 million km, beyond which solar gravity dominates).

Poincaré 1890s ^{[Poincare1890]} established the non-integrability of the general three-body problem in his Méthodes Nouvelles de la Mécanique Céleste, exhibiting an exponentially divergent stable / unstable-manifold tangle at the L1 saddle that violates analytic integrability and forecasts deterministic chaos. The patched-conic approximation used for spacecraft trajectories (Voyager, Cassini, New Horizons, Parker Solar Probe) chains successive two-body Kepler arcs joined at sphere-of-influence crossings: a heliocentric arc, a Jupiter-centric hyperbolic flyby, a heliocentric departure arc with new orbital elements. The Voyager 2 Grand Tour (Jupiter 1979, Saturn 1981, Uranus 1986, Neptune 1989) exploits a $176$ -year alignment of the outer planets that permits four sequential gravity-assist Kepler manoeuvres, increasing the spacecraft's heliocentric energy at each flyby by twice the planet's orbital velocity projected onto the relative-velocity vector. The Interplanetary Transport Network (Koon-Lo-Marsden-Ross 2000) exploits the stable / unstable manifold structure around L1 and L2 to construct very-low-energy trajectories between Earth orbit, the lunar L1/L2 region, and the inner Lagrange points of other planets — the technological output of the dynamical-systems generalisation of Lagrange's 1772 result.

Theorem 5 (Roche limit and tidal disruption). A self-gravitating fluid satellite of density $ρ_{s}$ orbiting a primary of density $ρ_{p}$ and radius $R_{p}$ is tidally disrupted if its orbital radius falls below the Roche limit $d_{Roche} \approx 2.44 R_{p} (ρ_{p} / ρ_{s})^{1/3}$ .

Roche's 1849 derivation ^{[Roche 1849 *Mém. Acad. Sci. Lett. Montpellier* 1, 243-262]} balances the tidal stretching across the satellite (proportional to $G M_{p} R_{s} / d^{3}$ ) against the satellite's self-gravity ( $G M_{s} / R_{s}^{2} \propto G ρ_{s} R_{s}$ ). For a rigid satellite the numerical prefactor is $1.26$ ; the $2.44$ figure is for a fluid body (Saturn's rings, located at 1.10-2.27 Saturn radii, are inside the Roche limit and cannot coalesce into a moon). The same calculation, applied to a comet (Shoemaker-Levy 9, broken into 21 fragments by Jupiter's tides in 1992 before impacting Jupiter in 1994), to neutron-star inspirals (the matter is sheared at $\sim 20$ km separation, gauging the equation of state at supranuclear density), and to the gradual outward migration of the Moon (3.8 cm per year, tracked by lunar laser ranging since the Apollo retroreflectors of 1969-1972) furnishes some of the most precise tests of the central-force framework.

Theorem 6 (Schwarzschild perihelion precession). In the post-Newtonian expansion of general relativity, a bound orbit around a central mass $M$ acquires a perihelion advance of $Δ ω = 6 π GM / (c^{2} a (1 - e^{2}))$ per orbit, beyond the Newtonian closed-ellipse result.

For Mercury ( $a = 0.387$ AU, $e = 0.206$ , $T = 87.97$ d), this gives 42.98 arcseconds per century, matching the observed anomaly of $43.11 \pm 0.21$ ″/cy that remained after subtracting the Newtonian planetary perturbations (Le Verrier 1859, Newcomb 1882). The relativistic prediction was Einstein's 1915 calculation ^{[Einstein 1915]} in his definitive paper on the field equations and remains the first quantitative test of general relativity. The same formula yields 8.6 ″/cy for Venus, 3.8 ″/cy for Earth, and 100s/cy for the Hulse-Taylor binary pulsar PSR B1913+16, which (in the strong-field regime) tests GR to one part in $1 0^{14}$ . The Newtonian-Kepler problem is the leading-order term in a systematic post-Newtonian expansion treated in 13.05.03.

Theorem 7 (Rutherford-Bohr atomic transition). Combined with quantization of angular momentum $L = n ℏ$ ( $n = 1, 2, 3, \dots$ ), the Kepler-problem solution for the Coulomb potential $V = - e^{2} / (4 π ϵ_{0} r)$ yields the bound-state energies $E_{n} = - m_{e} e^{4} / (2 (4 π ϵ_{0})^{2} ℏ^{2} n^{2}) = - 13.6 eV / n^{2}$ of the hydrogen atom.

This is the Bohr 1913 atom ^{[Bohr 1913 *Phil. Mag.* (6) 26, 1-25]}, the first quantum-mechanical theory of an atom. Rutherford's 1911 scattering experiments ^{[Rutherford 1911 *Phil. Mag.* (6) 21, 669-688]} had established the nuclear-atom model; Bohr's hand-tooled quantization rule reproduced the Rydberg-Ritz spectroscopic combination principle and the Balmer-Lyman-Paschen series of hydrogen emission lines to part-per-thousand precision before any wave-mechanical or matrix-mechanical formalism existed. The full quantum treatment in 12.06.01 replaces Bohr's quantization with the eigenvalue problem $\hat{H} ψ = E ψ$ for the hydrogenic Hamiltonian; remarkably, the energy spectrum is unchanged — Bohr's ansatz hit the right answer for the wrong reason (he assumed circular orbits, but the spectrum is determined by the $so (4)$ symmetry, which is independent of the orbital shape).

Theorem 8 (Yarkovsky and YORP effects). The asymmetric thermal re-emission of absorbed sunlight by a rotating asteroid produces a recoil acceleration of order $1 0^{- 10}$ m/s $^{2}$ for a 100-m near-Earth asteroid, perturbing its orbit by tens of kilometres per century and altering its spin state.

The Yarkovsky effect (named after the 19th-century Russian engineer who described it in 1900) is a small non-gravitational perturbation that has been measured directly for several near-Earth asteroids: most famously 6489 Golevka (Chesley et al. 2003, Science 302, 1739) and 101955 Bennu (target of NASA's OSIRIS-REx mission, which returned samples in 2023). The YORP (Yarkovsky-O'Keefe-Radzievskii-Paddack) effect is the rotational analogue: thermal photons exert a net torque on irregularly shaped bodies, spinning them up to fission speeds in some cases. Both effects are inverse-square-Coulomb (thermal radiation pressure) corrections to the central-force Kepler picture and are essential for predicting impact probabilities of potentially hazardous asteroids.

Theorem 9 (Hulse-Taylor binary and gravitational radiation). The binary pulsar PSR B1913+16 (Hulse-Taylor 1974) is a Kepler-like orbit ( $a \sim 1.95 \times 1 0^{9}$ m, $e = 0.617$ , $T = 7.75$ hr) whose period decreases by $\dot{T} / T = - 2.40 \times 1 0^{- 12}$ s/s in agreement with general-relativistic gravitational-wave emission to the 0.1% level.

The quadrupole-radiation formula of Peters 1964 gives the orbit-averaged energy loss $d E / d t = - (32/5) G^{4} μ^{2} M^{3} / (c^{5} a^{5}) f (e)$ with $f (e) = (1 + (73/24) e^{2} + (37/96) e^{4}) / (1 - e^{2})^{7/2}$ for elliptical orbits, leading to a slow secular shrinkage of the semi-major axis $\overset{a}{˙} = - (64/5) G^{3} μ M^{2} / (c^{5} a^{3}) f (e)$ and corresponding period decrease via Kepler's third law. Hulse and Taylor (1993 Nobel Prize) measured this orbital decay over decades via pulse timing, demonstrating quantitative agreement with GR's gravitational-radiation prediction long before the 2015 direct detection of GW150914. The same Kepler-orbit-plus-radiation-reaction framework applies to the merger waveforms detected by LIGO/Virgo: GW150914's signal exhibits the characteristic chirp from a Kepler-like inspiral, with the orbital frequency sweeping through the LIGO band from $\sim 35$ to $\sim 250$ Hz in the final 0.2 seconds as two black holes spiral together.

Synthesis. The Kepler-Rutherford framework is the foundational reason that classical mechanics is more than an empirical compendium of force laws: it is a tightly woven analytical structure in which a single set of conservation laws — energy, angular momentum, and (for $V \propto 1/ r$ alone) the Laplace-Runge-Lenz vector — completely determines the dynamics. The central insight is that the inverse-square law is doubly distinguished: it is the only attractive central potential (alongside the isotropic oscillator) for which all bound orbits close (Bertrand's theorem), and the LRL conservation that enforces this closure builds toward the $so (4)$ symmetry that organises the hydrogen-atom spectrum 12.06.01. Putting these together identifies the classical and quantum versions of the Coulomb problem as instances of the same algebraic structure: in classical mechanics the LRL vector orients the elliptical orbit, in quantum mechanics it organises the principal-quantum-number degeneracy. The bridge is between Newton's 1687 derivation of Kepler's laws and Pauli's 1926 derivation of the hydrogen spectrum — both rest on the same hidden symmetry, separated by 239 years of mathematical development and by the introduction of the quantum-mechanical formalism in 12.01.01.

The pattern recurs across physics. Gravitational wave inspirals (LIGO observations from GW150914 onward) are post-Newtonian expansions whose leading order is the Kepler problem proved here, with relativistic corrections coming from a perihelion-precession term 13.05.03 generalised to the strong-field regime. Spacecraft slingshot trajectories (Voyager, Cassini, New Horizons) solve the restricted three-body problem in successive segments of two-body Kepler arcs joined at the spheres of influence. Rutherford-style scattering generalises to all of nuclear and particle physics: the same $1/ sin^{4} (θ /2)$ singularity at small scattering angle persists in the Mott formula for electron scattering off nuclei, in the Born approximation for slow neutron scattering off nuclei, and in the Rosenbluth formula for electron-proton scattering from form factors. The Coulomb logarithm that arises when one integrates the Rutherford cross section over impact parameters cuts off at the Debye length in a plasma, producing the classical Coulomb-collision frequency that controls transport in fusion reactors and in the solar wind. Every one of these branches grew from the two-body problem solved in this unit.

Connections Master

Conservation laws 09.01.03. The Kepler-problem solution rests entirely on the conserved energy and angular momentum proved in the preceding unit. Energy conservation reduces the radial motion to a one-dimensional problem in the effective potential; angular-momentum conservation confines the motion to a fixed plane and yields Kepler's second law as an immediate consequence. The Laplace-Runge-Lenz vector is the extra (Kepler-specific) conserved quantity that enforces closure of bound orbits.
Action principle 09.02.01. The two-body central-force Lagrangian reduces to a single-particle problem in a central potential, and the Euler-Lagrange equations yield the Binet equation derived here. The action principle re-derives Kepler's laws from a variational integral and feeds into the Hamilton-Jacobi treatment that organises the orbit in terms of action variables.
Action-angle variables 09.06.01. The Kepler problem is the canonical example of a completely integrable Hamiltonian system; its three action variables $I_{r}, I_{θ}, I_{ϕ}$ generate the periodic motion and the energy depends only on the combination $J = I_{r} + I_{θ} + I_{ϕ}$ , which is the algebraic statement of the orbital closure proved in this unit. The Liouville-Arnold theorem applies in textbook form; the LRL super-integrability is what reduces the three frequencies to one.
Solar-system tests of GR 13.05.03. The Newtonian Kepler problem treated here is the leading-order term in a post-Newtonian expansion; the next-order term produces the Schwarzschild perihelion precession of $6 π GM / (c^{2} a (1 - e^{2}))$ per orbit, the classical observable that distinguishes general relativity from Newtonian gravity. Mercury's 43″/century anomaly, lunar laser ranging, and Gravity Probe B all generalise the framework of this unit to the relativistic regime.
Hydrogen atom 12.06.01. The quantum-mechanical Coulomb problem inherits the orbit equation and the LRL vector from this classical treatment, with the conservation of $A$ becoming the operator identity that fixes the hydrogen spectrum's accidental degeneracy. Pauli's 1926 derivation uses precisely the $so (4)$ algebra introduced here.

Historical & philosophical context Master

Kepler 1609 ^{[Kepler 1609]} derived his first two laws from Tycho Brahe's naked-eye observations of Mars's position over twenty years. He found that no combination of circular epicycles could match the Martian data to better than eight arcminutes — the limit of the naked eye — and was forced to abandon two thousand years of Aristotelian-Ptolemaic-Copernican circular dogma in favour of a non-circular curve. After several false starts (an egg-shaped orbit, an oval) he identified the ellipse as the correct shape, established the area-sweep rule, and published both in Astronomia Nova. The third law (the $T^{2} \propto a^{3}$ harmonic relation) followed ten years later in Harmonices Mundi ^{[Kepler 1619]}, where Kepler also speculated on a "harmony of the spheres" reading musical intervals into the orbital ratios.

Newton 1687 ^{[Newton 1687]} derived all three laws from his law of universal gravitation in Principia Book I, working backwards: he proved that an inverse-square central force produces conic-section orbits (Propositions 11-13), that the centre of force coincides with one focus, and that the area sweep rate and the period-axis relation follow as corollaries. The geometric proofs in Principia use a synthetic style that obscures the underlying ODE structure; the analytic formulation in terms of the Binet equation and the polar-coordinate ODE is due to the eighteenth-century mathematicians (Bernoulli, Euler, Clairaut), culminating in Laplace's Mécanique Céleste ^{[Laplace 1799]}, where the eccentricity vector that came to be known as the Laplace-Runge-Lenz vector first appears in the context of perturbation theory for planetary orbits. Runge 1919 ^{[Runge 1919]} and Lenz 1924 ^{[Lenz 1924]} gave the modern derivations, with Lenz applying the vector to the perturbed Kepler problem of the old quantum theory.

Bertrand 1873 ^{[Bertrand 1873]} proved the closed-orbits theorem that restricts the dynamical exception to Kepler and the harmonic oscillator. Pauli 1926 ^{[Pauli 1926]} discovered that the LRL vector becomes a conserved quantum operator on the hydrogen-atom Hilbert space and used the resulting $so (4)$ algebra to compute the energy spectrum a year before Schrödinger's wave-mechanical derivation. Bargmann 1936 ^{[Bargmann 1936]} gave the modern statement of the group-theoretic content. The Kepler problem's role as the canonical "hidden-symmetry" example in mathematical physics — visible in classical mechanics, in non-relativistic quantum mechanics, in the conformal compactification of momentum space (Fock 1935), in symplectic geometry as the affine motion on the cotangent bundle of $S^{3}$ — places it among the most-revisited problems in the discipline.

Rutherford 1911 ^{[Rutherford 1911]} turned the central-force framework into an experimental probe of atomic structure. His students Geiger and Marsden had measured the angular distribution of alpha particles scattered off gold foil and observed that one in eight thousand were deflected by more than 90° — a result Rutherford famously described as "almost as incredible as if you fired a 15-inch shell at a piece of tissue paper and it came back and hit you." Working out the differential cross section from the hyperbolic-orbit solution of the repulsive Coulomb problem, Rutherford deduced that the positive charge must be concentrated in a region smaller than $\sim 1 0^{- 14}$ m — the nucleus. Bohr 1913 ^{[Bohr 1913]} added the quantization rule $L = n ℏ$ to the Rutherford atom and produced the first correct theoretical value of the Rydberg constant, opening the path to modern atomic physics.

Bibliography Master

@book{Kepler1609,
  author = {Kepler, Johannes},
  title = {Astronomia Nova},
  year = {1609},
  publisher = {Heitzhofer},
  address = {Prague},
}

@book{Kepler1619,
  author = {Kepler, Johannes},
  title = {Harmonices Mundi},
  year = {1619},
  publisher = {Plancus},
  address = {Linz},
}

@book{Newton1687,
  author = {Newton, Isaac},
  title = {Philosophi{\ae} Naturalis Principia Mathematica},
  year = {1687},
  publisher = {Royal Society},
  address = {London},
}

@article{Rutherford1911,
  author = {Rutherford, Ernest},
  title = {The scattering of {$\alpha$} and {$\beta$} particles by matter and the structure of the atom},
  journal = {Philosophical Magazine (6)},
  volume = {21},
  year = {1911},
  pages = {669--688},
}

@article{Bohr1913,
  author = {Bohr, Niels},
  title = {On the constitution of atoms and molecules},
  journal = {Philosophical Magazine (6)},
  volume = {26},
  year = {1913},
  pages = {1--25},
}

@article{Bertrand1873,
  author = {Bertrand, Joseph},
  title = {Th{\'e}or{\`e}me relatif au mouvement d'un point attir{\'e} vers un centre fixe},
  journal = {Comptes Rendus},
  volume = {77},
  year = {1873},
  pages = {849--853},
}

@book{Laplace1799,
  author = {Laplace, Pierre-Simon},
  title = {Trait{\'e} de M{\'e}canique C{\'e}leste, Vol. 1},
  year = {1799},
  publisher = {Crapelet},
  address = {Paris},
}

@book{Runge1919,
  author = {Runge, Carl},
  title = {Vektoranalysis},
  year = {1919},
  publisher = {Hirzel},
  address = {Leipzig},
}

@article{Lenz1924,
  author = {Lenz, Wilhelm},
  title = {{\"U}ber den Bewegungsverlauf und die Quantenzust{\"a}nde der gest{\"o}rten Keplerbewegung},
  journal = {Z. Phys.},
  volume = {24},
  year = {1924},
  pages = {197--207},
}

@article{Pauli1926,
  author = {Pauli, Wolfgang},
  title = {{\"U}ber das Wasserstoffspektrum vom Standpunkt der neuen Quantenmechanik},
  journal = {Z. Phys.},
  volume = {36},
  year = {1926},
  pages = {336--363},
}

@article{Bargmann1936,
  author = {Bargmann, Valentine},
  title = {Zur Theorie des Wasserstoffatoms},
  journal = {Z. Phys.},
  volume = {99},
  year = {1936},
  pages = {576--582},
}

@article{Roche1849,
  author = {Roche, {\'E}douard Albert},
  title = {La figure d'une masse fluide soumise {\`a} l'attraction d'un point {\'e}loign{\'e}},
  journal = {M{\'e}moires de la Section des Sciences de l'Acad{\'e}mie des Sciences et des Lettres de Montpellier},
  volume = {1},
  year = {1849},
  pages = {243--262},
}

@article{Einstein1915,
  author = {Einstein, Albert},
  title = {Erkl{\"a}rung der Perihelbewegung des Merkur aus der allgemeinen Relativit{\"a}tstheorie},
  journal = {Sitzungsber. Preuss. Akad. Wiss.},
  year = {1915},
  pages = {831--839},
}

@book{Goldstein2002,
  author = {Goldstein, Herbert and Poole, Charles P. and Safko, John},
  title = {Classical Mechanics},
  edition = {3rd},
  year = {2002},
  publisher = {Addison-Wesley},
}

@book{LandauLifshitz1976,
  author = {Landau, Lev D. and Lifshitz, Evgeny M.},
  title = {Mechanics},
  edition = {3rd},
  year = {1976},
  publisher = {Pergamon},
  series = {Course of Theoretical Physics Vol. 1},
}

@book{Arnold1989,
  author = {Arnold, Vladimir I.},
  title = {Mathematical Methods of Classical Mechanics},
  edition = {2nd},
  year = {1989},
  publisher = {Springer},
  series = {GTM 60},
}

@book{Taylor2005,
  author = {Taylor, John R.},
  title = {Classical Mechanics},
  year = {2005},
  publisher = {University Science Books},
}

@book{Poincare1890,
  author = {Poincar{\'e}, Henri},
  title = {Les M{\'e}thodes Nouvelles de la M{\'e}canique C{\'e}leste},
  year = {1892--1899},
  publisher = {Gauthier-Villars},
  address = {Paris},
  note = {3 volumes},
}

@article{Peters1964,
  author = {Peters, Philip C.},
  title = {Gravitational radiation and the motion of two point masses},
  journal = {Phys. Rev.},
  volume = {136},
  year = {1964},
  pages = {B1224--B1232},
}

@article{HulseTaylor1975,
  author = {Hulse, Russell A. and Taylor, Joseph H.},
  title = {Discovery of a pulsar in a binary system},
  journal = {Astrophys. J.},
  volume = {195},
  year = {1975},
  pages = {L51--L53},
}

@book{HallidayResnickWalker2014,
  author = {Halliday, David and Resnick, Robert and Walker, Jearl},
  title = {Fundamentals of Physics},
  edition = {10th},
  year = {2014},
  publisher = {Wiley},
}

Prerequisites

09.01.01
09.01.02
09.01.03

Tier anchors

beginner: Halliday, Resnick & Walker, *Fundamentals of Physics*, 10e (Wiley, 2014), Ch. 13; Khan Academy lessons on Kepler's laws and gravitation
intermediate: Taylor, *Classical Mechanics* (University Science Books, 2005), §8; Goldstein, Poole & Safko, *Classical Mechanics*, 3e (Addison-Wesley, 2002), §3
master: Goldstein, Poole & Safko, *Classical Mechanics*, 3e (Addison-Wesley, 2002), §3 (full Kepler treatment, Laplace-Runge-Lenz vector); Landau & Lifshitz, *Mechanics*, 3e (Pergamon, 1976), §11-19; Arnold, *Mathematical Methods of Classical Mechanics*, 2e (Springer GTM 60, 1989), §8

References

Kepler — Astronomia Nova (Prague, 1609) · Part III (Mars orbit, areas law); Part IV (ellipse law)
Kepler — Harmonices Mundi (Linz, 1619) · Book V, Ch. 3 (harmonic law T^2 ∝ a^3)
Newton — Philosophiæ Naturalis Principia Mathematica (Royal Society, 1687) · Book I, Propositions 1-17, 30-41 (Kepler's laws from inverse-square gravitation)
Rutherford — The scattering of α and β particles by matter and the structure of the atom (Phil. Mag. 6 21, 1911) · Pages 669-688 (Coulomb-scattering cross section; nuclear-atom model)
Bohr — On the constitution of atoms and molecules (Phil. Mag. 6 26, 1913) · Pages 1-25 (quantum quantization of Kepler orbits in the atom)
Bertrand — Théorème relatif au mouvement d'un point attiré vers un centre fixe (Comptes Rendus 77, 1873) · Pages 849-853 (closed-orbits classification)
Laplace — Traité de Mécanique Céleste, Vol. 1 (Paris, 1799) · Livre II, Ch. III (eccentricity vector of planetary orbits)
Runge — Vektoranalysis (Hirzel, Leipzig, 1919) · §§ on conserved vectors in central-force motion
Lenz — Über den Bewegungsverlauf und die Quantenzustände der gestörten Keplerbewegung (Z. Phys. 24, 1924) · Pages 197-207 (perturbed-Kepler conserved vector)
Pauli — Über das Wasserstoffspektrum vom Standpunkt der neuen Quantenmechanik (Z. Phys. 36, 1926) · Pages 336-363 (LRL algebra and hydrogen spectrum)
Bargmann — Zur Theorie des Wasserstoffatoms (Z. Phys. 99, 1936) · Pages 576-582 (SO(4) symmetry of bound Kepler problem)
Roche — La figure d'une masse fluide soumise à l'attraction d'un point éloigné (Mém. Acad. Sci. Lett. Montpellier 1, 1849) · Pages 243-262 (tidal-disruption radius for a fluid satellite)
Poincaré — Les Méthodes Nouvelles de la Mécanique Céleste (Gauthier-Villars, 1892-1899), 3 vols. · Vol. III, §§391-401 (non-integrability of the three-body problem; homoclinic-orbit tangle at L1)
Goldstein, Poole & Safko — Classical Mechanics, 3e (Addison-Wesley, 2002) · §3 (central force problem, Kepler orbits, scattering)
Landau & Lifshitz — Mechanics, 3e (Pergamon, 1976) · §11-19 (Kepler problem; scattering; Coulomb collisions)
Arnold — Mathematical Methods of Classical Mechanics, 2e (Springer GTM 60, 1989) · §8 (central-force motion as integrable system; Bertrand's theorem)

Estimated time

beginner: 18m
intermediate: 40m
master: 80m