Perturbations and Perihelion Precession

In 1859 Urbain Le Verrier — the man who had discovered Neptune with a pen, by reading its gravity in the perturbations of Uranus — announced a scandal at the other end of the Solar System. After accounting for the tug of every known planet on Mercury, its perihelion was still creeping forward by a stubborn 38 seconds of arc per century (later refined to 43). The machinery that had conjured a new planet out of Uranus's wobbles now demanded one inside Mercury's orbit. Astronomers hunted the phantom planet "Vulcan" through eclipses for half a century. It was never found, because it does not exist: the missing 43 arcseconds were the first visible crack in Newton's law of gravity itself. This article builds the classical machinery behind that story. A pure $1/r$ potential is special: its bound orbits are perfect closed ellipses whose major axis never moves, a degeneracy protected by the conserved Laplace–Runge–Lenz vector. Add any small correction $\delta V(r)$ — another planet's pull, the flattening of the central body, the first whisper of general relativity — and the closure breaks. The orbit becomes a slowly turning ellipse, a rosette, and the turning rate is exquisitely measurable because it accumulates forever. You already own every tool this needs: the effective potential, the apsidal-angle integral, Binet's equation. Here they combine into a first-order formula converting any $\delta V(r)$ into a precession rate — the calculation that failed to explain Mercury with matter and then, in Einstein's hands in November 1915, explained it exactly with geometry.

Key takeaways

A pure inverse-square force produces closed ellipses with a fixed major axis — a fragile degeneracy protected by the conserved Laplace–Runge–Lenz vector, not a generic feature of central forces.
Any perturbation $\delta V(r)$ shifts the apsidal angle $\Phi$ away from the Kepler value $\pi$ , and the perihelion turns by $\Delta\phi = 2\Phi - 2\pi$ per radial period: advance if the perturbation strengthens the attraction at small $r$ , regression if it weakens it.
For near-circular orbits, $\Phi = \pi/\sqrt{3 + (r/F)(\mathrm{d}F/\mathrm{d}r)}$ ; for $V = kr^n$ this gives $\Phi = \pi/\sqrt{n+2}$ , closed only for $n=-1$ and $n=2$ (Bertrand's theorem).
First-order perturbation theory compresses everything into Landau's formula $\Delta\phi = \frac{\partial}{\partial L}\big(\frac{2m}{L}\int_0^{\pi} r^2\,\delta V\,\mathrm{d}\psi\big)$ , an integral of $\delta V$ over the unperturbed ellipse.
General relativity contributes an effective $-1/r^3$ correction, giving $\Delta\phi_{\mathrm{GR}} = 6\pi GM/[c^2 a(1-e^2)]$ — for Mercury $0.104''$ per orbit, $43.0''$ per century, exactly the observed residual.
Key idea: closed orbits are a degeneracy of the exact $1/r$ potential; any deviation makes the apsidal angle $\Phi \neq \pi$ , so the perihelion precesses by $\Delta\phi = 2\Phi - 2\pi$ per orbit — computable, to first order, as an integral of the perturbing potential over the unperturbed Kepler ellipse.

The core insight

Closed orbits are special; any deviation from $1/r^2$ makes the apsidal angle $\Phi\neq\pi$ , so the perihelion precesses by $\Delta\phi=2\Phi-2\pi$ per revolution — a quantity computable, to first order, as an integral of the perturbing potential over the unperturbed Kepler ellipse.

Where this lives today

Perihelion precession has been promoted from anomaly to precision instrument. The GRAVITY interferometer tracked the star S2 through its 16-year orbit around Sagittarius A* and in 2020 detected its Schwarzschild precession — about 12 arcminutes per orbit, formula (9) evaluated at a four-million-solar-mass focus — and the upgraded GRAVITY+ is pushing the same measurement to fainter S-stars deeper in the potential. In the Solar System, the BepiColombo mission, scheduled to settle into Mercury orbit in late 2026, will radio-track exactly this orbit to tighten the relativistic parameters and the Sun's quadrupole moment simultaneously. And for the double pulsar J0737-3039 the measured periastron advance is $16.9°$ per year — formula (9) as the leading term of a post-Keplerian expansion tested to parts in $10^4$ . The same physics runs satellite engineering. Earth's oblateness ( $J_2$ ) adds a short-range perturbation that rotates a low satellite's perigee by several degrees per day — a hundred thousand times the relativistic term — and mission designers exploit rather than fight it: Molniya-type communications orbits sit at the critical inclination $63.4°$ , where the $J_2$ apsidal precession vanishes ( $5\cos^2 i - 1 = 0$ ) and the apogee stays parked over the service hemisphere. Anything that orbits, precesses; the art is in reading which $\delta V$ is talking.

The physical picture

Strip the orbit problem down to its two clocks. The radial coordinate $r$ oscillates in the one-dimensional effective potential $V_{\mathrm{eff}}(r)$ between $r_{\min}$ and $r_{\max}$ — one clock, ticking with the radial period. Meanwhile the polar angle $\phi$ winds steadily forward under angular momentum conservation — the second clock. Closure is a question of resonance between the two clocks: the path repeats only if, during one full radial oscillation, the angle advances by a rational multiple of $2\pi$ . The Kepler potential achieves the most perfect resonance imaginable — the angle advances by exactly $2\pi$ per radial period, at every energy and every eccentricity. That is why the ellipse closes.

Figure

Two fates for an orbit — the Kepler ellipse closes on itself; a perturbed orbit becomes a rosette whose axis turns by a fixed angle each revolution

Here is the crucial reframe: closure is a degeneracy, not a robust property. The inverse-square law hides an extra conserved quantity beyond energy and angular momentum — the Laplace–Runge–Lenz vector, which points permanently from the focus toward perihelion. As long as it is conserved, the major axis is pinned to the stars. Perturb the potential by any amount, however small, and that conservation law fails: the LRL vector begins to rotate, and its rotation rate is the precession rate. Nothing dramatic happens to the orbit's shape or energy — each pass still looks like the same ellipse — but the axis drifts, and drift accumulates: an effect of one part in ten million per orbit becomes an unmistakable signal after a century of orbits. Which way does the axis drift? A perturbation that strengthens the attraction at small $r$ (an extra $-\alpha/r^2$ or $-\beta/r^3$ term, say) makes the radial clock lag the angular one: a radial cycle now completes only after the angle has swept more than $2\pi$ . The perihelion arrives late in angle — it has moved forward. So extra attraction at short range means advance (precession in the same sense as the orbit); anything that weakens the force relative to $1/r^2$ produces regression. This one sign rule instantly predicts that general relativity — effectively extra attraction at small $r$ — must advance Mercury's perihelion, not retard it.

Figure

The apsidal angle is the whole story — Kepler sweeps exactly 180 degrees from perihelion to aphelion; a perturbed orbit sweeps more or less, and the excess is the precession

The single number capturing all of this is the apsidal angle $\Phi$ : the polar angle swept while $r$ runs from perihelion to aphelion. Kepler orbits have $\Phi = \pi$ exactly, so a full radial period sweeps $2\pi$ and the figure closes. Everything that follows is a set of increasingly powerful ways to compute $\Phi - \pi$ .

The apsidal machinery

Start from what conservation laws alone provide. In any central potential, the two-body problem reduces to a single body of reduced mass $m$ with conserved $E$ and $L$ , and the radial motion is one-dimensional motion in the effective potential — the launching point for every formula in this article:

V_{\mathrm{eff}}(r) = V(r) + \frac{L^2}{2mr^2} \qquad (1)

where: $V(r)$ — the central potential, here the Kepler term $-k/r$ plus a small perturbation $\delta V(r)$ ; $L$ — the conserved angular momentum; $m$ — the reduced mass. The centrifugal term $L^2/2mr^2$ is what turns the radial problem into oscillation between two turning points $r_{\min}$ and $r_{\max}$ , the roots of $E = V_{\mathrm{eff}}(r)$ . To find how much angle is swept during that radial oscillation, divide the two rates the conservation laws give you. Angular momentum fixes $\dot\phi = L/mr^2$ , and energy conservation fixes $\dot r = \sqrt{(2/m)(E - V_{\mathrm{eff}})}$ ; their ratio $\mathrm{d}\phi/\mathrm{d}r$ integrates to the apsidal angle:

\Phi = \int_{r_{\min}}^{r_{\max}} \frac{(L/r^2)\,\mathrm{d}r}{\sqrt{2m\big(E - V_{\mathrm{eff}}(r)\big)}} \qquad (2)

where: $\Phi$ — the polar angle accumulated as $r$ climbs once from perihelion to aphelion. The integrand diverges like an inverse square root at both endpoints (where $\dot r \to 0$ ), but integrably so. For $V = -k/r$ this integral works out to exactly $\pi$ , independent of $E$ and $L$ : that stubborn, structureless answer is the closure degeneracy in analytic form. The observable is the mismatch between one radial period's angular sweep, $2\Phi$ , and a full turn:

\Delta\phi = 2\Phi - 2\pi \qquad (3)

where: $\Delta\phi$ — the net rotation of the major axis per radial period, positive for advance (rotation in the same sense as the orbital motion), negative for regression. For Mercury's relativistic anomaly $\Delta\phi \approx 5.0\times 10^{-7}\ \mathrm{rad}$ — a tenth of an arcsecond per orbit, yet decisive over a century of 415 orbits. For actually solving orbits, the workhorse is Binet's equation — the equation of motion rewritten with $u = 1/r$ as a function of angle $\phi$ instead of time, which converts orbit geometry questions into oscillator questions:

\frac{\mathrm{d}^2 u}{\mathrm{d}\phi^2} + u = -\frac{m}{L^2 u^2}\,F(1/u) \qquad (4)

where: $u = 1/r$ ; $F = -\mathrm{d}V/\mathrm{d}r$ — the radial force. For the Kepler force $F = -k/r^2$ the right side is the constant $mk/L^2$ , so $u$ oscillates about it with angular frequency exactly 1 in $\phi$ — perihelion recurs every $\Delta\phi_{\mathrm{recur}} = 2\pi$ of angle, the closed ellipse again. Any force that makes the right side $u$ -dependent shifts that frequency away from 1, and the orbit precesses.

Derivation AExact precession from an inverse-cube perturbation

We are deriving the precession caused by a perturbation $\delta V = -\alpha/r^2$ , whose force correction falls as $1/r^3$ — the one case Binet's equation solves exactly, with no small-perturbation expansion at all. The strategy: write the total force for $V = -k/r - \alpha/r^2$ and insert it into Binet's equation (4); because $\delta V \propto 1/r^2$ has the same $r$ -dependence as the centrifugal term, it only retunes the oscillator's frequency, leaving the equation linear. The thing to watch is that frequency: perihelion recurs when the phase $\beta\phi$ advances by $2\pi$ , so any $\beta \neq 1$ forces the angle $\phi$ to overshoot or undershoot a full turn. This exactly solvable case calibrates every approximate method that follows.

Step 1. Write the force and feed it to Binet. With $V = -k/r - \alpha/r^2$ , the force is $F = -k/r^2 - 2\alpha/r^3$ , i.e. $F(1/u) = -ku^2 - 2\alpha u^3$ . Then

\frac{\mathrm{d}^2u}{\mathrm{d}\phi^2} + u = -\frac{m}{L^2u^2}\,F = \frac{m}{L^2}\left(k + 2\alpha u\right)

Why this stays exactly solvable: the Kepler term lands as a constant, and the perturbation lands as a term linear in $u$ — no $u^2$ or worse appears, so this is still a driven harmonic oscillator in disguise.

Step 2. Collect the terms proportional to $u$ on the left:

\frac{\mathrm{d}^2u}{\mathrm{d}\phi^2} + \left(1 - \frac{2m\alpha}{L^2}\right)u = \frac{mk}{L^2}

Why the frequency shifted: the perturbation has exactly the radial profile of the centrifugal barrier, so it acts like a small change of the effective angular momentum — it retunes the restoring constant of the $u$ -oscillator from $1$ to $\beta^2 = 1 - 2m\alpha/L^2$ .

Step 3. Solve. The general solution is a displaced oscillation,

u(\phi) = \frac{mk}{L^2\beta^2} + A\cos(\beta\phi), \qquad \beta = \sqrt{1 - \frac{2m\alpha}{L^2}}

with $A$ set by the eccentricity. Perihelion (maximum $u$ ) recurs when $\beta\phi$ has advanced by $2\pi$ , i.e. after an angle $\phi = 2\pi/\beta$ .

Step 4. Read off the precession per orbit:

\Delta\phi = \frac{2\pi}{\beta} - 2\pi \;\approx\; \frac{2\pi m\alpha}{L^2} \quad (\alpha \text{ small}) \qquad (5)

where: the approximation expands $\beta^{-1} \approx 1 + m\alpha/L^2$ . Why the sign is right: $\alpha > 0$ means extra attraction, $\beta < 1$ , so the $u$ -oscillation is slow in $\phi$ and the angle overshoots $2\pi$ before perihelion returns — the perihelion advances, confirming the short-range-attraction sign rule. *As a check: $\alpha \to 0$ gives $\beta = 1$ and $\Delta\phi = 0$ , recovering the closed Kepler ellipse; and applying the first-order formula (7) below to $\delta V = -\alpha/r^2$ returns the same result.*

For nearly circular orbits there is a one-line diagnostic. Linearizing Binet's equation about the circular radius $r_0$ (write $u = u_0 + \epsilon$ , expand the right side of (4) to first order) gives a simple harmonic oscillator $\epsilon'' + \beta^2\epsilon = 0$ whose stiffness is the local slope of the force law:

\Phi = \frac{\pi}{\beta}, \qquad \beta^2 = 3 + \frac{r}{F}\frac{\mathrm{d}F}{\mathrm{d}r}\bigg|_{r_0} \qquad (6)

where: $F(r)$ — the attractive central force evaluated at the circular-orbit radius $r_0$ ; $\beta$ — the ratio of radial to angular frequency. For a power-law potential $V = kr^n$ (force $\propto r^{n-1}$ ) this collapses to $\Phi = \pi/\sqrt{n+2}$ : Kepler ( $n=-1$ ) gives $\Phi = \pi$ , the harmonic oscillator ( $n=2$ ) gives $\Phi = \pi/2$ — and Bertrand's theorem sharpens this into the statement that those two are the only potentials whose orbits close at every energy, not just near circularity. Note the built-in stability check: $\beta^2 < 0$ (force falling faster than $1/r^3$ ) means no oscillation at all — circular orbits are unstable and there is no apsidal angle to speak of.

Which formula answers which question

you want	use	notes
definition of the precession	$\Delta\phi = 2\Phi - 2\pi$ , eqs. (2)+(3)	exact, any central potential
quick sign and size, near-circular orbit	$\Phi = \pi/\sqrt{3 + rF'/F}$ , eq. (6)	one line; $V=kr^n$ gives $\pi/\sqrt{n+2}$
exact answer for a $1/r^2$ potential correction	Binet retuning, eq. (5)	the only exactly solvable perturbation
any small $\delta V(r)$ , any eccentricity	Landau's formula, eq. (7)	integrate over the unperturbed ellipse
the relativistic advance	$\Delta\phi_{\mathrm{GR}} = 6\pi GM/c^2 p$ , eq. (9)	$p = a(1-e^2)$ ; always an advance

Go deeper · Derivation B · 2 minThe first-order precession formula for any small perturbation

Derivation BThe first-order precession formula for any small perturbation

We are deriving Landau's compact master formula: for any small perturbation $\delta V(r)$ added to the Kepler potential, the precession per orbit is $\Delta\phi = \frac{\partial}{\partial L}\big(\frac{2m}{L}\int_0^{\pi} r^2\,\delta V\,\mathrm{d}\psi\big)$ , with the integral taken over the unperturbed ellipse. The strategy has three moves: express the swept angle per radial period as an $L$ -derivative of the radial action integral; expand that integral to first order in $\delta V$ ; then trade the radial integration variable for the true anomaly $\psi$ of the unperturbed orbit. The point to watch is the first-order logic — evaluating the correction on the unperturbed orbit is legitimate precisely because the correction is already first-order small; chasing the orbit's own distortion would be a second-order refinement.

Step 1. Write the angle swept in one radial period as an action derivative. With radial momentum $p_r = \sqrt{2m(E - V) - L^2/r^2}$ ,

\Phi_{\mathrm{full}} = -\frac{\partial}{\partial L}\oint p_r\,\mathrm{d}r

Why this identity holds: differentiating $p_r$ under the integral gives $\partial p_r/\partial L = -L/(r^2 p_r)$ , and $L\,\mathrm{d}r/(r^2 p_r) = (L/mr^2)(\mathrm{d}r/\dot r) = \dot\phi\,\mathrm{d}t = \mathrm{d}\phi$ — the derivative literally counts swept angle. This is the same identity that powers action-angle theory.

Step 2. Expand the action to first order in the perturbation. With $V = -k/r + \delta V$ and $p_r^{(0)}$ the unperturbed radial momentum,

\oint p_r\,\mathrm{d}r = \oint p_r^{(0)}\,\mathrm{d}r \;-\; \oint \frac{m\,\delta V}{p_r^{(0)}}\,\mathrm{d}r

Why the zeroth-order piece drops out of the physics: $-\partial_L\oint p_r^{(0)}\mathrm{d}r = 2\pi$ exactly — the closed Kepler orbit — so the entire precession lives in the correction term.

Step 3. Convert the correction to an integral over the true anomaly. Along the unperturbed orbit, $m\,\mathrm{d}r/p_r^{(0)} = \mathrm{d}r/\dot r = \mathrm{d}t$ , and angular momentum conservation gives $\mathrm{d}t = (mr^2/L)\,\mathrm{d}\psi$ ; one radial period corresponds to $\psi$ running from $0$ to $2\pi$ , i.e. twice the range $0$ to $\pi$ :

\oint \frac{m\,\delta V}{p_r^{(0)}}\,\mathrm{d}r = \frac{2m}{L}\int_0^{\pi} r^2\,\delta V(r)\,\mathrm{d}\psi

Step 4. Assemble the two pieces:

\Delta\phi = \Phi_{\mathrm{full}} - 2\pi = \frac{\partial}{\partial L}\left(\frac{2m}{L}\int_0^{\pi} r^2\,\delta V(r)\,\mathrm{d}\psi\right) \qquad (7)

where: $\psi$ — the true anomaly, the polar angle from perihelion of the unperturbed ellipse; $r(\psi) = p/(1+e\cos\psi)$ with $p = L^2/mk$ and $e = \sqrt{1 + 2EL^2/mk^2}$ . The recipe: substitute the unperturbed orbit, do the $\psi$ -integral, and only then differentiate with respect to $L$ at fixed $E$ — the $L$ -dependence hiding inside $p$ and $e$ is part of the derivative, but $\delta V$ itself must be held as a fixed function of $r$ (see the warning block below). As a consistency check: applying (7) to $\delta V = -\alpha/r^2$ gives $\int_0^\pi r^2\delta V\,\mathrm{d}\psi = -\pi\alpha$ , so $\Delta\phi = \partial_L(-2\pi m\alpha/L) = 2\pi m\alpha/L^2$ — exactly the Binet result (5).

The unperturbed orbit that feeds formula (7) is worth displaying once in full, since every application below leans on it. The bound Kepler orbit is the focus-centered conic

r(\psi) = \frac{p}{1 + e\cos\psi}, \qquad p = \frac{L^2}{mk}, \qquad e = \sqrt{1 + \frac{2EL^2}{mk^2}} \qquad (8)

where: $p$ — the semi-latus rectum, related to the semi-major axis by $p = a(1-e^2)$ ; $e$ — the eccentricity. The key property for precession integrals: $1/r$ is linear in $\cos\psi$ , so powers of $1/r$ integrate against $\mathrm{d}\psi$ in closed form — which is why short-range perturbations give such clean answers.

Go deeper · Derivation C · 2 minThe forty-three arcseconds of Mercury

Derivation CThe forty-three arcseconds of Mercury

We are deriving the general-relativistic perihelion advance, $\Delta\phi_{\mathrm{GR}} = 6\pi GM/[c^2a(1-e^2)]$ per orbit, and evaluating it for Mercury. The input from general relativity is a single fact: the Schwarzschild geometry corrects the Newtonian orbit equation by one extra term, turning Binet's equation (per unit mass, with $h = L/m$ the specific angular momentum) into $u'' + u = GM/h^2 + (3GM/c^2)u^2$ . The strategy is textbook perturbation theory on an oscillator: solve the Kepler part, feed that solution into the small quadratic term, and keep only the piece of the response that is resonant — oscillating at the orbit frequency itself — because only resonant forcing accumulates orbit after orbit. The non-resonant pieces merely wiggle the ellipse's shape by parts in $10^7$ ; the resonant piece masquerades as a slow frequency shift, which is precisely a precession.

Step 1. Write the relativistic orbit equation and identify the small parameter:

\frac{\mathrm{d}^2u}{\mathrm{d}\phi^2} + u = \frac{GM}{h^2} + \frac{3GM}{c^2}\,u^2

Why the correction is tiny: the new term over the Newtonian term is $3h^2u^2/c^2 \sim 3(v_\perp/c)^2$ , about $2.6\times10^{-8}$ for Mercury — comfortably in first-order-perturbation territory.

Step 2. Insert the unperturbed Kepler solution $u_0 = (GM/h^2)(1 + e\cos\phi)$ into the quadratic term. Expanding $u_0^2 \propto 1 + 2e\cos\phi + e^2\cos^2\phi$ , the forcing splits into a constant, a $\cos 2\phi$ piece (via $\cos^2\phi$ ), and the term $2e\cos\phi$ . Why only $\cos\phi$ matters: the homogeneous solution oscillates as $\cos\phi$ , so forcing at that same frequency is resonant; the constant and $\cos2\phi$ pieces produce only bounded, non-accumulating distortions.

Step 3. Solve for the resonant response. Forcing $A\cos\phi$ on the left-hand operator $u'' + u$ yields the secular particular solution $(A/2)\,\phi\sin\phi$ ; with $A = (3GM/c^2)\cdot 2e\,(GM/h^2)^2$ this gives

\delta u = \frac{3(GM)^3 e}{c^2 h^4}\,\phi\sin\phi

Step 4. Recognize the growing term as a frequency shift. Since $\cos\big((1-\delta)\phi\big) \approx \cos\phi + \delta\,\phi\sin\phi$ for small $\delta$ ,

u_0 + \delta u \approx \frac{GM}{h^2}\Big[1 + e\cos\big((1-\delta)\phi\big)\Big], \qquad \delta = \frac{3(GM)^2}{c^2h^2}

Why this reading is forced on us: a term growing linearly in $\phi$ cannot persist in a bounded orbit; summing the perturbation series to all orders repackages it as a slightly detuned cosine — the ellipse is fine, its clock runs slow.

Step 5. Extract the advance. Perihelion recurs when $(1-\delta)\phi$ gains $2\pi$ , i.e. at $\phi = 2\pi(1+\delta)$ to first order, so with $h^2 = GMp = GMa(1-e^2)$ :

\Delta\phi_{\mathrm{GR}} = 2\pi\delta = \frac{6\pi (GM)^2}{c^2h^2} = \frac{6\pi GM}{c^2\,a(1-e^2)} \qquad (9)

where: $M$ — the central mass; $a, e$ — the orbit's semi-major axis and eccentricity. Every factor is positive: general relativity always advances the perihelion, and it punishes small, eccentric orbits hardest — hence Mercury.

Step 6. Put in Mercury's numbers ( $GM_\odot = 1.327\times10^{20}\ \mathrm{m^3/s^2}$ , $p = 5.55\times10^{10}\ \mathrm{m}$ ): $\Delta\phi = 5.02\times10^{-7}$ rad $= 0.104''$ per orbit; at $415$ orbits per century, $43.0''$ per century. The full arithmetic is laid out in worked example 2 below. As a cross-check: the observed anomalous advance is $42.98'' \pm 0.04''$ per century; feeding $\delta V_{\mathrm{eff}} = -GML^2/(mc^2r^3)$ through formula (7) with its coefficient held fixed gives the identical $6\pi GM/c^2p$ .

Limiting cases

Case	Result	Meaning
pure Kepler $V = -k/r$	$\Phi = \pi$ , $\Delta\phi = 0$	closed ellipse at every $E, L$ ; the LRL vector pins the axis
isotropic oscillator $V = \tfrac12 m\omega^2r^2$	$\Phi = \pi/2$	also closed (center-based ellipse, two radial cycles per turn) — the other Bertrand potential
power law $V = kr^n$ , near-circular	$\Phi = \pi/\sqrt{n+2}$	precesses for every $n$ except $-1$ and $2$ ; circular orbits unstable for $n \le -2$
flat rotation curve $F \propto 1/r$	$\Phi = \pi/\sqrt2 \approx 127.3°$	rosette orbits of stars in galactic halos — apsides regress $105°$ per radial period
GR correction	$\Delta\phi = \dfrac{6\pi GM}{c^2a(1-e^2)} > 0$	always advance; strongest for small $a$ , large $e$ — Mercury, close binaries, S-stars
$\delta V = -\alpha/r^2$	$\Delta\phi = 2\pi m\alpha/L^2$ exactly	the calibration case: Binet solves it with no approximation

Common trapPlugging the relativistic correction into the first-order formula gives 2π, not 6π — so the formula must be wrong

The formula is fine; the trap is in the differentiation. The GR effective perturbation is $\delta V = -GML^2/(mc^2r^3)$ — it carries an explicit factor $L^2$ , inherited from the velocity-dependence of the relativistic correction. Differentiate formula (7) treating that $L^2$ as live and two thirds of the answer cancels: you get $2\pi GM/c^2p$ , too small by exactly 3. The rule: Landau's formula was derived for a perturbing potential, a fixed function of $r$ alone — freeze the coefficient ( $\beta \equiv GML^2/mc^2$ , evaluated on the actual orbit), differentiate only the $L$ -dependence of the unperturbed motion, and substitute back at the end. Done that way, $\delta V = -\beta/r^3$ gives $\Delta\phi = 6\pi\beta m^2k/L^4 = 6\pi GM/c^2p$ , agreeing with the resonance derivation. Two smaller traps: $\Delta\phi$ is the net axis rotation $2\Phi - 2\pi$ , not the apsidal angle $\Phi$ itself; and the radial period (perihelion to perihelion) is not the angular period (return to the same direction in space) — they differ precisely because the orbit precesses.

Worked examples

1 · A star's rosette orbit in a galactic halo (confidence builder)

Stars in the flat-rotation-curve region of a galaxy feel a central force $F \propto 1/r$ (this is what "flat rotation curve" means: $v_c = \sqrt{Fr/m}$ independent of $r$ ). For a star like the Sun on a nearly circular orbit at $R = 8\ \mathrm{kpc}$ with $v_c = 220\ \mathrm{km/s}$ , find the apsidal angle, the precession per radial period, and the time between perihelion passages.

Step 1. Apply the near-circular diagnostic (6). For $F = c_1/r$ : $(r/F)(\mathrm{d}F/\mathrm{d}r) = (r/F)(-c_1/r^2) = -1$ , so $\beta^2 = 3 - 1 = 2$ and $\Phi = \pi/\sqrt2 = \boxed{127.3°}$ .

Step 2. Precession per radial period from (3): $\Delta\phi = 2\Phi - 360° = 254.6° - 360° = \boxed{-105.4°}$ — the apsides regress by about $105°$ each radial cycle: a rosette that never closes, the flower-petal orbit of every galactic-dynamics simulation.

Step 3. Time scale. The angular frequency is $\Omega = v_c/R = 2.2\times10^5/(8\times3.09\times10^{19}) = 8.9\times10^{-16}\ \mathrm{rad/s}$ , an orbital period $2\pi/\Omega \approx 224\ \mathrm{Myr}$ — the galactic year. The radial (epicyclic) frequency is $\beta\Omega = \sqrt2\,\Omega$ , so perihelion passages come every $224/\sqrt2 = \boxed{\approx 158\ \mathrm{Myr}}$ .

Step 4. Sanity-check the sign: a $1/r$ force falls off more slowly than $1/r^2$ — relatively weaker attraction at small $r$ — so the perihelion should regress. It does. (Same physics, opposite sign to general relativity.)

2 · Mercury versus Earth — computing the relativistic advance (exam level)

Compute the general-relativistic perihelion advance per orbit and per century for Mercury ( $a = 5.79\times10^{10}\ \mathrm{m}$ , $e = 0.206$ , $T = 88.0$ d) and for Earth ( $a = 1.496\times10^{11}\ \mathrm{m}$ , $e = 0.0167$ , $T = 365.25$ d). Use $GM_\odot = 1.327\times10^{20}\ \mathrm{m^3/s^2}$ , $c = 2.998\times10^8\ \mathrm{m/s}$ .

Step 1. Mercury's semi-latus rectum: $p = a(1-e^2) = 5.79\times10^{10}\times(1 - 0.0424) = 5.55\times10^{10}\ \mathrm{m}$ .

Step 2. Advance per orbit from (9): the numerator is $6\pi GM_\odot = 18.85\times1.327\times10^{20} = 2.50\times10^{21}$ , the denominator $c^2p = 8.99\times10^{16}\times5.55\times10^{10} = 4.99\times10^{27}$ , so $\Delta\phi = 2.50\times10^{21}/4.99\times10^{27} = 5.02\times10^{-7}\ \mathrm{rad}$ . In arcseconds: $5.02\times10^{-7}\times206{,}265 = \boxed{0.104''\ \text{per orbit}}$ .

Step 3. Per century: Mercury completes $36{,}525/88.0 = 415$ orbits, so $415\times0.104'' = \boxed{43.0''\ \text{per century}}$ — precisely Le Verrier's unexplained residual.

Step 4. Earth: $p = 1.496\times10^{11}\times(1-0.00028) = 1.4956\times10^{11}\ \mathrm{m}$ , so $\Delta\phi = 2.50\times10^{21}/(8.99\times10^{16}\times1.4956\times10^{11}) = 1.86\times10^{-7}\ \mathrm{rad} = 0.038''$ per orbit; at 100 orbits per century, $\boxed{3.8''\ \text{per century}}$ .

Step 5. Read the comparison physically. Earth's per-orbit advance is only $2.7$ times smaller than Mercury's, but per century Mercury wins by a factor of $11$ : the $1/p$ scaling and the short period multiply. That double advantage is why the innermost planet, and no other, carried the first evidence against Newtonian gravity.

Origins · 1 minAlbert Einstein and how this was discovered

Albert Einstein

1879–1955 · Berlin, November 1915

In November 1915, midway through the four legendary Thursday lectures in which he wrestled general relativity into its final form, Einstein set his brand-new field equations the one test the Solar System had held open for fifty-six years: Mercury's unexplained 43 arcseconds. He solved the orbit equation perturbatively — essentially the resonance calculation on this page — and out came $43''$ per century, from geometry alone, with no adjustable constant anywhere in the theory. He later told friends he had palpitations of the heart and was "beside myself with joy" for days. The calculation also marked a turning point in method: an anomaly that perturbation theory had failed to explain with unseen matter (the phantom planet Vulcan) was resolved by changing the law of gravity itself — making precession the archetype for every "missing mass or new physics?" debate since, dark matter included.

History

1846

— Le Verrier predicts Neptune's position from unexplained perturbations of Uranus; Galle finds it within one degree — perturbation theory's greatest triumph.

1859

— Le Verrier turns the same machinery on Mercury and reports an unexplained perihelion advance; the hunt for the intramercurial planet Vulcan begins.

1882

— Simon Newcomb refines the anomalous residual to about 43 arcseconds per century; decades of Vulcan searches during eclipses have found nothing.

1915

— Einstein computes the perihelion advance from general relativity: 43 arcseconds per century, no free parameters. He reports palpitations of the heart.

1974

— Hulse and Taylor discover the binary pulsar PSR B1913+16; its periastron advances 4.2 degrees per year — thirty-five thousand Mercurys — opening strong-field tests of relativity.

2020

— The GRAVITY interferometer detects the Schwarzschild precession of the star S2 around the Milky Way's central black hole, the same formula at four million solar masses.

Figure

Mercury's precession budget — of the observed 574 arcseconds per century, Newtonian planetary perturbations supply 531; the last 43 belong to general relativity

Going further

Texts. Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., §3.6, §3.8 · Landau & Lifshitz, Mechanics, 3rd ed., §14–15 and the §15 problems · Morin, Introduction to Classical Mechanics, ch. 7 · Weinberg, Gravitation and Cosmology, §8.6 · Taylor, Classical Mechanics, ch. 8. The near-circular formula $\Phi=\pi/\sqrt{3+(r/F)F'}$ reduces in the power-law case to $\Phi=\pi/\sqrt{n+2}$ . The $1/r^3$ Binet result is $\Delta\phi=2\pi m\alpha/L^2$ . The GR advance $\Delta\phi=6\pi GM/c^2 a(1-e^2)$ evaluated with $a=5.79\times10^{10}$ m, $e=0.206$ , $M_\odot=1.989\times10^{30}$ kg, $G=6.674\times10^{-11}$ N m² kg $^{-2}$ , $c=2.998\times10^8$ m/s, and period $87.97$ d gives Mercury's $43''$ /century. The $574''$ / $531''$ / $43''$ precession budget and the Hulse–Taylor periastron advance of $4.2°$ /yr are standard observational values.

Builds on: The Effective Potential · The Orbit Equation and Binet's Equation · The Kepler Problem · Angular Momentum Conservation in Central Fields · Conic-Section Orbits. Leads to: Bertrand's Theorem · General Relativistic Orbits · Kepler's First Law: Elliptical Orbits.

Does the orbit lose energy or angular momentum as it precesses?

No — precession is geometry, not dissipation. Both $E$ and $L$ are exactly conserved in any static central potential, perturbed or not. What is lost is a third, accidental conservation law: the Laplace–Runge–Lenz vector, which only the exact $1/r$ potential protects. With the LRL vector no longer constant, nothing pins the major axis, and it rotates at a steady rate while the orbit's size, shape, and energy stay fixed. Contrast a genuinely dissipative effect like drag or gravitational-wave emission, which shrinks $a$ and circularizes the orbit. The Hulse–Taylor binary shows both at once: a huge conservative periastron advance ( $4.2°$ /yr) and a slow orbital decay from radiated waves — two different physical effects, cleanly separated in the timing data.

Why are the Kepler potential and the harmonic oscillator the only two with closed orbits?

That is Bertrand's theorem, and the near-circular formula (6) shows where the shortlist comes from: for $V = kr^n$ the apsidal angle is $\Phi = \pi/\sqrt{n+2}$ . Closure at all energies requires $\Phi$ to stay a rational multiple of $\pi$ as the orbit becomes eccentric — but away from circularity, $\Phi$ generally drifts with energy. Bertrand's argument tracks that drift and finds it vanishes only for $n = -1$ ( $\Phi = \pi$ always) and $n = 2$ ( $\Phi = \pi/2$ always). Both special cases hide an extra conserved quantity — the LRL vector for Kepler, a symmetric tensor invariant for the oscillator — and hidden symmetry, not coincidence, is what closure means. Any other force law closes at most at isolated, fine-tuned energies.

I applied Landau's formula to the GR correction and got 2π GM/c²p instead of 6π. What did I do wrong?

You differentiated the $L^2$ inside the perturbation. The relativistic effective correction $\delta V = -GML^2/(mc^2r^3)$ carries an explicit $L^2$ because it descends from a velocity-dependent term, not from a genuine static potential. Landau's formula (7) says: for a perturbation that is a fixed function of $r$ , the precession is $\partial_L$ of the orbit-averaged perturbation. The $\partial_L$ is supposed to probe how the unperturbed motion responds to a change of angular momentum — the perturbing function must be held fixed while you differentiate. So set $\beta = GML^2/(mc^2)$ , treat $\delta V = -\beta/r^3$ with $\beta$ constant, compute $\Delta\phi = 6\pi\beta m^2k/L^4$ , and only then substitute $\beta$ 's value: the result is $6\pi GM/c^2p$ , matching the resonance derivation exactly. The factor-3 discrepancy is a classic exam trap precisely because both procedures look plausible.

Mercury's perihelion actually advances 574 arcseconds per century — why do we only celebrate 43 of them?

Because 531 of the 574 are ordinary Newtonian physics, computed to high precision by Le Verrier's perturbation methods: Venus contributes about $277''$ (closest), Jupiter $154''$ (heaviest), Earth $90''$ , the rest about $10''$ — each piece is formula (7) with $\delta V$ the orbit-averaged potential of one perturbing planet. By 1859 this budget was solid enough that the remaining $43''$ constituted a crisis. The celebrated number is the residual — the part no arrangement of known matter could supply, hunted as the planet Vulcan for half a century, and finally delivered by general relativity with zero adjustable parameters. The discovery is only as good as the boring bookkeeping underneath it.

Which way does an oblate central body drive the perihelion?

An equatorial bulge adds a quadrupole term falling as $1/r^3$ ; for a near-equatorial orbit this is extra attraction at small $r$ — the same structure as the GR term — so it advances the perigee, and for a low satellite the effect is enormous by celestial standards: several degrees per day. The three-dimensional story has a twist: at inclination $i$ the apsidal rate scales as $(5\cos^2 i - 1)$ , which changes sign at $i = 63.43°$ . Below the critical inclination the perigee advances, above it regresses, and at it the perigee freezes — the trick behind Molniya and Tundra communications orbits, whose apogees must loiter over high latitudes for decades without expensive station-keeping. One formula, one sign flip, one billion-dollar orbit design.

Is the precession rate the same thing as the rotation rate of the Laplace-Runge-Lenz vector?

Yes, exactly — and it is often the slickest way to compute it. The LRL vector $\mathbf{A} = \mathbf{p}\times\mathbf{L} - mk\,\hat{\mathbf{r}}$ points from focus to perihelion with magnitude $mke$ ; under the pure Kepler force, $\dot{\mathbf{A}} = 0$ . Add a perturbing force and $\langle\dot{\mathbf{A}}\rangle$ , averaged over one unperturbed orbit, comes out perpendicular to $\mathbf{A}$ with magnitude $|\mathbf{A}|\,\Delta\phi/T$ — pure rotation at the precession rate, with no change of length to first order (no change of eccentricity). This orbit-averaged LRL method is standard in celestial mechanics and atomic physics alike, and it makes the geometry vivid: the perturbation grabs the one arrow that Kepler froze, and slowly swings it around the focus.

Perturbations and Perihelion Precession

Where this lives today

The physical picture

The apsidal machinery

Which formula answers which question

Limiting cases

Worked examples

1 · A star's rosette orbit in a galactic halo (confidence builder)

2 · Mercury versus Earth — computing the relativistic advance (exam level)

Going further

Before this

Where this leads