Astronomy & Celestial Mechanics
the sky, computed
Concept

Effective potential — where it appears

The radial potential obtained by folding a body's conserved angular momentum into the central potential, so that the motion becomes one-dimensional. The centrifugal barrier it produces is why a body with any angular momentum at all cannot reach the centre, however much energy it has.

Named by 8 essays across 3 fields — each of them below, with the objects they name alongside it.

A vector that does not change, drawn five times. The Laplace–Runge–Lenz vector A = v × L − GM r̂, constructed at five points of one orbit at e = 0.45 under a force ∝ r^−2, each drawn from the body rather than from the focus so that its length and direction can be compared point by point. Under the inverse square every one of them is the same vector: identical to machine precision in length and in direction. Its length is 0.4500, which is the eccentricity, and it points at pericentre — so the orbit's orientation is a conserved quantity and not a constant of integration, which is why the ellipse does not turn.

The third thing that is conserved

Energy fixes an orbit's size and angular momentum fixes its shape. Neither fixes which way it points — so a curve that closes needs a third conserved quantity, and only two force laws in the universe supply one.

orbits · Conic sections The effective potential, for three angular momenta. The radial motion of an orbiting body is one-dimensional motion in an effective potential: the attraction −GM/r plus the centrifugal term L²/2r² that the angular momentum contributes. The barrier at small radius is what stops a body with any angular momentum at all from reaching the centre, and the bottom of each well is the circular orbit.

The wall that angular momentum builds

A body falling toward a star almost never arrives. Sideways motion, which looks like a detail of the initial conditions, turns the attraction into a well with a wall around the middle of it — and the wall is why it costs more fuel to hit the Sun than to leave the solar system.

orbits · Effective potential The five Lagrange points at mass fraction 0.12. The five points at which a small body can keep station with two larger ones. The three on the line of centres are roots of a quintic and are unstable; the two forming equilateral triangles are stable for a sufficiently lopsided mass ratio.

Five places that keep station, in a problem with no solution

Three bodies under gravity cannot be solved. Restrict the problem slightly and five exact answers fall out anyway — three of them roots of a quintic, two of them perfect equilateral triangles.

gravitation · Lagrange points Three bodies, integrated. Three equal masses integrated forward under mutual gravity: the figure-eight choreography — all three bodies on one closed curve. Every point is a step of the equations of motion, and the total energy is conserved to 2.5e-11 across the run.

Three bodies, and what "no solution" actually means

The three-body problem is routinely called unsolvable. Trajectories are computed for it every day, exact periodic solutions are known, and both statements are true — the word is doing more work than it looks.

gravitation · The three-body problem Mercury's perihelion, term by term. The observed advance of Mercury's perihelion is 5600 arcseconds per century against the equinox. Almost all of it is the equinox: the coordinate frame itself turns, and removing it leaves 574. Subtracting the perturbations of the other planets — computed by Le Verrier in 1859 and refined many times since — leaves 42.98 arcseconds a century that nothing in Newtonian gravitation accounts for. The bars are logarithmic in nothing; they are the real proportions, which is why the residual is barely visible beside the frame term.

Forty-three arcseconds, after everything else

Mercury's perihelion moves through 5,600 arcseconds a century. Nearly all of that is the coordinate system turning, and almost all of the rest is the other planets. What was left over was 43 — under one per cent of the raw number, and the most consequential residual in the history of the subject.

gravitation · Relativistic orbits The horizon and the orbit that cannot come back, against the hole's spin. Two radii round a rotating black hole, in gravitational radii GM/c², against its spin a from −1 (an orbit against the rotation) to +1 (with it), for orbits in the equatorial plane. The lower curve is the horizon, 1 + √(1 − a²), which is the same whichever way the orbit goes. The upper one is the marginally bound orbit, 2 − a + 2√(1 − a): the closest a body falling in from far away can pass and still escape back out. For a hole with no spin the horizon is at 2 and the marginally bound orbit at 4 — twice as far out. With maximal spin the marginally bound orbit comes in to 1.09 for a prograde orbit (a = 0.998) and moves out to 5.83 for a retrograde one. A star whose tidal radius lies inside this curve is swallowed whole, so the curve, not the horizon, is the line a disruption flare is measured against — and it moves by a factor of 5.3 with the spin.

The line a star is swallowed at is not the horizon

A star torn apart by a black hole makes a flare, and a star swallowed whole makes nothing, so the heaviest hole that can produce a flare is a measurement. That ceiling is set not at the horizon but at the closest orbit from which infalling matter can still come back out — twice as far out for a hole that does not spin, and moved by a factor of five by the hole's rotation.

galaxies · Tidal disruption The two force laws whose orbits close. The apsidal angle — the angle swept from periapsis to the next apoapsis — against the exponent of the force law, for F ∝ r^p. The dashed curve is the near-circular limit π/√(3+p), which has a closed form; the solid curve is the same angle for an orbit of eccentricity 0.4, computed by quadrature of ∫(L/r²)dr/√(2(E−U)) between its two turning points, with U the effective potential. An orbit closes when the apsidal angle is a rational multiple of π, and an orbit closes at every eccentricity only where the two curves meet: p = −2 at exactly 180° and p = +1 at exactly 90°, which is Bertrand's theorem. The quadrature returns 180.0000° and 90.0000° at those two exponents and departs from the near-circular curve by 1.4° at p = 0. The angle diverges as p approaches −3, where the circular orbit stops being stable and there is no well left to oscillate in.

Only two force laws let an orbit come back

That a planet returns to the same point of its own path after one lap is not a fact about orbits. It is a fact about the exponent in the force, and out of the whole continuum of attractions only two — the inverse square, and a spring — bring every bound orbit back to where it started.

orbits · Effective potential The same well after the primary has lost 45 per cent of its mass. Two effective potentials for one body: the solid curve before the primary loses mass and the faint one after, both at the same angular momentum, because a central force of any strength exerts no torque. The well shallows and its floor moves out from r = 1.00 to 1.82. The body's own level moves with it — from E = -0.420 to -0.127 — and the two horizontal lines are drawn where the radial action is conserved, which puts the turning points at 0.71–1.67 before and 1.30–3.03 after. The ratio between them is 2.3333 in both, so the orbit is the same shape at a larger size: everything about the body's path has scaled and its eccentricity of 0.4 has not moved. That is what a slow change leaves behind, and it is not what a sudden one leaves.

The well moves, and the body keeps its share of it

When the Sun becomes a white dwarf it will throw away half its mass, and every planet's orbit will swell by the same factor. Their eccentricities will not change at all — provided the loss is slow, and the only meaning "slow" has here is slow compared with one orbital period.

orbits · Effective potential

Named alongside it

The objects these essays reach for when they reach for this one.

Angular momentumEccentricityApsidal precessionCentrifugal barrierInverse square lawLagrange pointsMass ratioOrbital stabilitySpecific orbital energyThree-body problemTurning pointAdiabatic invariant

All concepts