Astrophysics

The orbit special relativity cannot close

An inverse-square orbit closes on itself because of a conserved vector that nothing else has. Give the orbiting body an inertia that grows with its speed, as special relativity requires, and the vector turns — the orbit becomes a rosette, advancing by a sixth of Mercury's famous 43 arcseconds a century. A different theory that respects special relativity just as well turns the rosette backwards by the same amount. The 43 is not special relativity plus a correction; it is a measurement of what gravity pulls on.

Assumes: The arrow that says which way the orbit points · Mass is a form of energy, which is not the same as a source of it

An orbit under an exact inverse-square attraction comes back to the point it started from, and it does so because of a conserved vector pointing at the perihelion that no other force law has. Any change to the force turns that vector, and the orbit becomes a rosette. Mercury’s perihelion turns by 43 arcseconds a century more than the other planets explain, and the usual summary is that relativity is responsible. Which relativity is a question worth asking precisely, because the answer is not the one the summary suggests. Keep Newton’s attraction exactly as it is and change only one thing — let the orbiting body’s inertia grow with its energy, as special relativity insists it must — and the orbit already fails to close.

A rosette: the orbit of a body whose inertia is its energy. An orbit of eccentricity 0.45 under inertia from energy, Newton's pull, integrated for 7 revolutions from the closest approach of a body with k/Lc = 0.35. The dashed curve is the Newtonian ellipse from the same starting position and momentum, which closes on itself. The integrated orbit does not: its closest approach moves forward by 24.3° each revolution, against the exact 24.3° of 2π(1/Γ − 1) with Γ = √(1 − k²/L²c²). The dots are the successive closest approaches.
Fig. 1 Seven revolutions of a body in an exact inverse-square attraction whose inertia is its energy, integrated from the closest approach. The strength is exaggerated — the ratio k/Lck/Lc of the attraction to angular momentum times cc is 0.35, where Mercury’s is 1.6×1041.6 \times 10^{-4} — so that the effect is visible in a single turn. The dashed ellipse is Newton’s orbit from the same position and momentum, and it closes; the integrated orbit advances its closest approach by 24.3 degrees every revolution, which is what the exact solution gives.

Heavier where it turns hardest

The mechanism takes one sentence once the right quantity is watched. A body on an eccentric orbit is fastest at its closest approach, so that is where its inertia is largest.

Where on the orbit the body is heaviest. The factor γ by which a body's inertia exceeds its rest mass, through one radial period of an orbit of eccentricity 0.45 at k/Lc = 0.35, against the angle swept from the first closest approach. It peaks at 1.141 at closest approach, where the body is fastest, and falls to 1.021 at the far end. A push across the path turns a body by an amount inversely proportional to its inertia, so near closest approach this body is turned less sharply than a Newtonian one at the same speed would be, and it overshoots. The next closest approach arrives at 384.3° rather than 360°: the 24.3° of overshoot is the step the rosette takes each revolution. The shading marks where the inertia is more than halfway to its peak.
Fig. 2 The factor γ\gamma by which the body’s inertia exceeds its rest mass, through one radial period of the orbit in the first figure, against the angle swept from the closest approach. The body is 14 per cent heavier at closest approach than at rest and 2 per cent heavier at the far end. The next closest approach arrives at 384.3 degrees rather than 360: that overshoot of 24.3 degrees is the rosette’s step.

A push across a body’s path turns it by an amount inversely proportional to its inertia. That is the content of the relativistic second law for a force at right angles to the motion: the rate at which the direction changes is the force divided by γmv\gamma m v, not by mvmv. So at closest approach, where the attraction is strongest and the turn is sharpest, the relativistic body is turned less than a Newtonian one would be at the same speed. It swings wide, and the swing carries the whole orbit round. By the time the body comes back to its closest approach it has swept more than a full turn — 384.3 degrees in the figure — and the perihelion has advanced by the excess.

Nothing here requires curved spacetime, a field with energy of its own, or any modification of the inverse square. The attraction is Newton’s, the geometry is flat, and the only relativistic ingredient is that a moving body’s inertia is its energy. That is enough to break the one symmetry that made the Kepler orbit close.

It is worth being precise about which symmetry, because the answer says what kind of effect this is. Of all the central forces, only the inverse square and the linear spring close every bound orbit, and the inverse square does it because its orbits have a conserved vector over and above energy and angular momentum. Conserved quantities come from symmetries — each one is handed over by an invariance of the equations — and the one behind this vector is not a symmetry of space at all. It is a hidden symmetry of the Kepler problem’s equations in a four-dimensional space of positions and momenta, and it holds only because the kinetic energy is exactly p2/2mp^2/2m and the potential exactly k/r-k/r. Replace the kinetic energy by p2c2+m2c4\sqrt{p^2c^2 + m^2c^4} and the hidden symmetry is gone, whatever the potential is. The rotation of space survives — angular momentum is still conserved, the orbit still lies in a plane — but the orientation of the ellipse within the plane is no longer fixed by anything, and it drifts.

An exact rosette

The special-relativistic Kepler problem has an exact solution, and it is worth having, because it turns the qualitative story into a number and shows where the number stops making sense.

The body’s energy is E=p2c2+m2c4k/rE = \sqrt{p^2c^2 + m^2c^4} - k/r, with k=GMmk = GMm for gravity or k=Ze2/4πε0k = Ze^2/4\pi\varepsilon_0 for an electron round a nucleus. Squaring and separating the radial motion gives

pr2c2=(E+kr)2m2c4L2c2k2r2,p_r^2 c^2 = \left(E + \frac{k}{r}\right)^2 - m^2c^4 - \frac{L^2c^2 - k^2}{r^2},

and the last term is the telling one. In Newtonian mechanics the centrifugal barrier is L2/r2L^2/r^2; here it is (L2k2/c2)/r2(L^2 - k^2/c^2)/r^2. The attraction has eaten into the angular momentum’s barrier, because the attraction’s own energy k/rk/r is squared along with everything else and one piece of the square has the same 1/r21/r^2 shape as the barrier. The orbit that results is

1r=A[1+ecos(Γφ)],Γ=1k2L2c2,\frac{1}{r} = A\left[1 + e \cos(\Gamma \varphi)\right], \qquad \Gamma = \sqrt{1 - \frac{k^2}{L^2c^2}},

an ellipse whose angle has been rescaled by Γ\Gamma. The body returns to its closest approach after φ=2π/Γ\varphi = 2\pi/\Gamma rather than 2π2\pi, so the perihelion advances by 2π(1/Γ1)2\pi(1/\Gamma - 1) every revolution, which for a weak attraction is πk2/L2c2\pi k^2/L^2c^2. For gravity, where L2=GMm2a(1e2)L^2 = GMm^2 a(1 - e^2) for an orbit of semi-major axis aa, that is

Δφ=πGMc2a(1e2)per revolution.\Delta\varphi = \frac{\pi GM}{c^2 a(1-e^2)} \quad \text{per revolution}.

General relativity’s answer, the one the Laplace–Runge–Lenz argument reconciles with Mercury, is 6πGM/c2a(1e2)6\pi GM/c^2a(1-e^2). Special relativity alone supplies exactly one sixth of it.

Two properties of the exact rosette are not obvious from the approximate one. The advance depends on the angular momentum and on nothing else — not on the eccentricity, not on the energy — so a nearly circular orbit and a very eccentric one with the same LL precess at the same rate. And Γ\Gamma reaches zero at k=Lck = Lc, which is not a small correction becoming large but the end of the solution.

One curve, from a planet to the inside of a heavy atom

The ratio k/Lck/Lc is the only number the rosette depends on, and it runs over an enormous range across systems that obey the same equation.

One curve from a planet to a heavy atom, and where it ends. The advance of the closest approach per orbit, in radians, of a body whose inertia is its energy in an inverse-square attraction, against the ratio k/Lc of the attraction's strength to angular momentum times the speed of light, on logarithmic axes. The curve is 2π(1/√(1 − (k/Lc)²) − 1) and the open points are integrated orbits. Mercury's orbit round the Sun sits at k/Lc = 1.63·10⁻⁴ and advances 8.37·10⁻⁸ rad an orbit. The hydrogen atom's lowest orbit sits at k/Lc = 0.0073 and advances 1.67·10⁻⁴ rad an orbit. The innermost orbit in a mercury atom sits at k/Lc = 0.584 and advances 1.46 rad an orbit. The curve runs off the top at k/Lc = 1, where the orbit stops closing at all and the body falls in — the classical form of the limit near Z = 137 beyond which a point nucleus has no stable innermost orbit.
Fig. 3 The advance of the closest approach per revolution against k/Lck/Lc, on logarithmic axes. The curve is 2π(1/Γ1)2\pi(1/\Gamma - 1) and the open circles are integrated orbits. Mercury’s orbit round the Sun sits at k/Lc=1.6×104k/Lc = 1.6 \times 10^{-4} and advances 8.4×1088.4 \times 10^{-8} radians a revolution; the lowest orbit in hydrogen at k/Lc=αk/Lc = \alpha advances 1.7×1041.7 \times 10^{-4}; the innermost orbit in an atom of mercury, where the nuclear charge is 80, advances by 1.46 radians — 83 degrees — every turn. The dashed line at k/Lc=1k/Lc = 1 is where the curve ends.

For Mercury the advance is 8.4×1088.4 \times 10^{-8} radians a revolution. The planet goes round 415 times a century, so special relativity alone gives 7.2 arcseconds a century — a sixth of the 43 that are observed after the other planets are subtracted, which is exactly what the formulas above say it must be.

For hydrogen the ratio k/Lck/Lc in the lowest orbit is the fine-structure constant, α1/137\alpha \approx 1/137, and the advance is πα2\pi\alpha^2, a part in six thousand of a turn every revolution. This is the rosette Arnold Sommerfeld computed in 1916, three years after Bohr had put the electron on a circle whose radius survives in the modern atom as the most probable distance and nothing else. Quantising it — requiring the rosette’s radial and angular motions each to hold a whole number of quanta — split every level of the old Bohr atom by an amount of order α2\alpha^2, and gave a formula for hydrogen’s fine structure that agrees with measurement. It also agrees, exactly, with the formula Dirac derived twelve years later from a relativistic wave equation with the electron’s spin built in — although Sommerfeld’s electron had no spin at all. The agreement is not a sign that the rosette was right. Two omissions in the older calculation, the missing spin and the crude quantisation rule, cancel term by term, and the cancellation has since been traced to a hidden symmetry the two calculations share. It remains one of the strangest coincidences in the history of the subject, and it is the reason the fine structure constant has its name.

For the innermost electron of a heavy atom the ratio is ZαZ\alpha, and at Z=80Z = 80 it is 0.58. The rosette advances 83 degrees a turn; relativity is not a correction here but the dominant fact about the orbit. That is where the heavy elements start behaving in ways their lighter relatives do not, with inner orbits pulled in and everything outside them rearranged.

Past k=Lck = Lc there is no orbit

At k/Lc=1k/Lc = 1 the effective centrifugal barrier, (L2k2/c2)/r2(L^2 - k^2/c^2)/r^2, vanishes, and beyond it the term changes sign and joins the attraction.

An orbit past the point where special relativity lets one close at all. A body with a relativistic inertia moving in an inverse-square attraction with k/Lc = 1.05, integrated from its closest approach. At this strength the attraction wins against the centrifugal term for every radius: the body does not reach an outer turning point but spirals inward and falls into the centre within a single revolution. Below k/Lc = 1 every bound orbit is a rosette; at and above it there are no orbits at all.
Fig. 4 A body started from what would be a closest approach at k/Lc=1.05k/Lc = 1.05. There is no outer turning point and no second closest approach: the 1/r21/r^2 term that holds every other orbit away from the centre has changed sign, the body spirals inward, and it reaches the centre in less than a revolution. The dashed curve is the Newtonian ellipse the same start would give.

The Newtonian barrier is invincible: at small enough radius the L2/r2L^2/r^2 term beats any 1/r1/r attraction, so an orbit with any angular momentum at all has a closest approach. The relativistic barrier is not, because the attraction’s square has the same 1/r21/r^2 shape and can outweigh it. When it does, nothing stops the fall. A classical body with k>Lck > Lc is not in an orbit with a large precession; it is not in an orbit.

In an atom the lowest possible LL is fixed by quantum mechanics, and the classical condition k>Lck > Lc becomes Zα>1Z\alpha > 1 for the innermost state — a nuclear charge near 137. The full quantum treatment moves the threshold and softens it, and a nucleus of finite size moves it again, to a charge of about 170, where the innermost level would plunge into the negative-energy continuum and the vacuum near the nucleus would be expected to produce electron–positron pairs spontaneously. The rosette’s collapse at k=Lck = Lc is the classical shadow of that limit, and it appears from nothing more than an inertia that grows with energy and an inverse-square attraction.

Three relativistic theories, three answers

The claim that special relativity “supplies one sixth” needs its qualification stated plainly, because it is where the argument actually lands. The rosette above assumed one particular coupling: the attraction is exactly Newton’s k/rk/r, and only the body’s inertia changes with speed. That is the natural model for an electron round a nucleus. It is not the only relativistically consistent model of gravity, and the obvious alternatives give different answers.

Three relativistic theories of one orbit, and three answers. The advance of the closest approach per orbit, in units of πGM/c²ℓ where ℓ is the orbit's semi-latus rectum, measured by integrating one weak-field orbit under each of three couplings that all respect special relativity. Inertia from energy, Newton's pull: 0.99. A scalar field (Nordström): −1.01. Curved spacetime (general relativity): 6.11. For Mercury the unit is 7.2 arcseconds a century, so the three theories predict +7.2, −7.2 and +43.0. Special relativity fixes none of the three by itself; what fixes the answer is what gravity couples to.
Fig. 5 The advance per revolution in units of πGM/c2\pi GM/c^2\ell, where \ell is the orbit’s semi-latus rectum, measured by integrating one weak-field orbit under each of three couplings. A body whose inertia is its energy in Newton’s attraction: +0.99. A body in Nordström’s scalar theory of gravity, in which the field couples to the body’s rest-frame energy and the pull weakens as the body speeds up: −1.01. The Schwarzschild geodesic of general relativity: +6.11. For Mercury one unit is 7.2 arcseconds a century.

Gunnar Nordström’s theory of 1913 is the instructive case. It is a perfectly Lorentz-invariant field theory of gravity with a single scalar field, it reproduces Newton in the weak static limit, and it was taken seriously by Einstein himself for a year. In it, a moving body’s coupling to the field is its energy measured in its own rest frame, which falls as 1/γ1/\gamma relative to its total energy. At closest approach, where the body is fastest, the pull is weakened more than the inertia is increased. The body is turned more sharply than a Newtonian one, not less, and the rosette runs backwards.

The same start under a scalar field: the rosette turns backwards. An orbit of eccentricity 0.45 under a scalar field (Nordström), integrated for 7 revolutions from the closest approach of a body with k/Lc = 0.35. The dashed curve is the Newtonian ellipse from the same starting position and momentum, which closes on itself. The integrated orbit does not: its closest approach moves backward by 29.7° each revolution. The dots are the successive closest approaches.
Fig. 6 The same start as the first figure, under Nordström’s scalar theory. The closest approach retreats by 29.7 degrees a revolution instead of advancing, and the rosette turns the other way. In the weak field the backward advance is one sixth of general relativity’s forward one — minus 7.2 arcseconds a century for Mercury.

Nordström’s theory failed on exactly this count, among others. It predicts that Mercury’s perihelion should regress by 7 arcseconds a century, where the measured residual is an advance of 43. It also predicts no bending of light at all, since its geometry differs from flat spacetime only by an overall scale factor, and a light ray’s path does not care about the scale of the geometry it crosses. The 1919 eclipse and the perihelion together ruled it out.

Which of the three is the honest special-relativistic theory of gravity is the wrong question, and the reason is instructive. The Sommerfeld coupling, natural for an electron round a nucleus, is unnatural for gravity: its attraction stays fixed while the body’s inertia grows, so a fast body is pulled with a force that does not grow with its energy, and it accelerates more slowly in the field than a slow one. That violates the universality of free fall — the principle, tested to parts in 101510^{15}, that everything falls alike whatever it is made of and however much of its mass is kinetic energy. For an electric charge there is no such principle, since charge and inertia are different quantities. For gravity there is, and it says the gravitational charge is the energy. Nordström’s theory respects it — it is a theory in which every body follows the geodesics of one geometry — and fails instead because its geometry is the wrong shape, a flat spacetime rescaled point by point, which bends no light and turns the perihelion backwards. General relativity respects it too, and is the only one of the three that also bends light by the measured amount.

So the famous 43 arcseconds are not special relativity with a correction added. Special relativity is a constraint every candidate satisfies, and the candidates it admits give answers ranging from minus one to plus six in the natural unit. What the measurement actually determines is how gravity couples to a moving body and how the field behaves at second order in its strength — the parameters that the post-Newtonian framework calls γ\gamma and β\beta, for which general relativity’s values combine to give six. The measured advance is one of the two or three numbers that pin them down.

The dependence on second-order structure is sharper than it sounds. Nordström’s field for a point mass is φ=1GM/rc2\varphi = 1 - GM/rc^2. Replace it with eGM/rc2e^{-GM/rc^2}, which agrees with it to first order and therefore gives exactly the same Newtonian gravity, the same light bending and the same everything else a first-order experiment can see — and the backward precession doubles, from −1 unit to −2. A precession is a second-order effect, so the (GM/rc2)2(GM/rc^2)^2 term in the field matters as much as the leading one. A theory of gravity is not fixed by its Newtonian limit, and the perihelion is one of the few measurements that looks past it.

Where the rosettes stop being the answer

The attraction is taken to be instantaneous. Every coupling here uses a static field centred on a fixed source. A source that moves, or a field that propagates at a finite speed, adds velocity-dependent terms of the same order as the effects computed, and for two bodies of comparable mass those terms are not optional.

The body does not radiate. A charge on the rosettes of the fourth figure would radiate electromagnetically and spiral inward in a fraction of a nanosecond; a classical atom is not stable at all, and the Sommerfeld orbits are a calculation performed on an object that cannot exist classically, redeemed only by the quantum rules applied to it. A planet radiates gravitationally, far too weakly to matter.

Spin is left out. An electron’s magnetic moment couples to the field it sees in its own frame, and that coupling is as large as the rosette’s own effect in hydrogen — it is why the missing spin and the crude quantisation happened to cancel. The frame the electron sees is itself rotating, because two boosts in different directions leave a rotation behind, and that kinematic rotation halves the naive magnetic coupling. A spinning planet’s spin couples to the curvature, negligibly for any real planet.

The general-relativistic bar is the test-particle geodesic. It assumes the orbiting body has no gravity of its own and that the central mass does not rotate. Both are excellent for Mercury and the Sun and neither is excellent for two neutron stars. General relativity has its own version of the collapse at k=Lck = Lc, and it arrives sooner: below six gravitational radii no circular orbit is stable at any angular momentum, because the curvature adds an attractive 1/r31/r^3 term that eventually beats any barrier. The special-relativistic collapse needs the attraction to outweigh LcLc; the general-relativistic one needs only that the body come close enough.

The clock carried round the rosette

The drawings show one plane and one orbit. What they cannot show is the time along it. A body on the special-relativistic rosette carries a clock that runs slow in proportion to its speed, most at closest approach, and a body on the general-relativistic orbit carries one slowed by both its speed and its depth in the field. The rosettes are the same shape as they would be on a distant observer’s clock, but a traveller on each counts a different number of its own seconds per revolution, and that difference is a second observable that none of these pictures carries.

Nor can the pictures show that the exaggeration is an exaggeration. Every rosette is drawn at a strength where the orbit is relativistic in every sense; Mercury’s special-relativistic rosette, drawn at the same scale, would be an ellipse whose closest approach takes seventy-five million revolutions to go once round — and the full general-relativistic one, twelve million.

Still open: how well the second-order term is known

The advance of Mercury’s perihelion is now measured by radio ranging to spacecraft in orbit around it, and the post-Newtonian parameters it constrains are known to parts in 10510^5: the space-curvature parameter from the deflection and delay of radio signals past the Sun, the second-order parameter from the perihelion and from the Moon’s orbit measured by laser ranging. Both agree with general relativity’s value of exactly one. Missions now in progress aim to improve the perihelion-based measurement by a further order of magnitude, and one of the obstacles is the Sun’s own flattening, whose contribution to the precession has a similar form and has to be measured separately. Whether the second-order coupling departs from one at a part in a million is not known, and theories that modify gravity at large distances generally predict that it should somewhere.

The natural question beyond this one is what happens when the orbiting body is not a test particle. Two bodies of comparable mass each move in the other’s field while radiating, and the orbit neither closes nor keeps its size. There is no rosette with a fixed step there, and computing what there is took most of a century. The habit worth carrying from here is to ask what a famous residual is a residual of. The 43 arcseconds are not relativity versus Newton; they are one theory’s coupling against the alternatives special relativity leaves open, and a sixth of them would have been there in a universe where space did not curve at all.

Part 4 of 4

This essay is one argument about Orbit stability. The others:

What links here

Essays that reach for this one mid-argument — the half of a link its own author cannot write down.

The objects named here

The third axis, after the field and the reading path: the things themselves, and every essay that touches each one.

Central forceConserved quantityEffective potentialFine structure constantGeneral relativityThe Lorentz factorOrbitPrecessionRelativistic dynamics