Perturbations — the series · Astronomy & Celestial Mechanics

The osculating semi-major axis of a perturbed orbit. The semi-major axis a test particle would have if the perturber vanished, computed from its position and velocity at every step of an integration over 26 orbits of the perturber. It is not constant: a short-period ripple rides on a slow trend, and only the trend accumulates.

Elements that do not stay constant

Six numbers fix an orbit for all time, and the phrase is only true in a universe containing two bodies. Add a third and the six start moving — some of them wandering and returning, one or two of them drifting in one direction forever, and the difference between those two behaviours is the whole of celestial mechanics after Newton.

part 1 · orbits

The same thrust is worth 2.08 times more at perigee, and out of plane it is worth nothing at all. The rate of change of the semi-major axis under a unit acceleration in each of the three directions, against position around an orbit of eccentricity 0.35, from Gauss's variational equations. The along-track curve carries the factor p/r = 1 + e cos f and therefore peaks at perigee, where the same impulse is worth 2.08 times what it is worth at apogee — the whole of the Oberth effect, arriving as a term in a differential equation rather than as an argument about kinetic energy. The radial curve is antisymmetric about apoapsis and integrates to exactly zero over a revolution: pushing outwards for half an orbit and being pushed back for the other half changes the energy by nothing, which is checked here by quadrature and comes out at -9.0e-18. And the out-of-plane response is identically zero at every point of the orbit, because W is perpendicular to the velocity and does no work. An orbital plane can be rotated without touching the energy, and that is why a plane change is so expensive: none of what is spent goes anywhere useful.

Which direction moves which element

Resolve a small force into three components and Gauss's equations say exactly what each one does. The out-of-plane component can rotate an orbit and can never change its energy; the along-track component owns the semi-major axis outright and is worth more at perigee than at apogee by a factor that is pure geometry.

part 2 · orbits

The two halves of what is left over. The disturbing potential of a perturber on a test particle, against the difference in longitude between them, at a semi-major axis ratio of 0.62. The upper curve is the direct term — the perturber's own attraction, which peaks at conjunction where the separation is smallest and falls to 1/(1+α) half a turn later. The lower one is the indirect term, which exists only because the coordinates are centred on a primary that is itself being accelerated, and which is a pure cosine of the longitude difference. The indirect term is the larger of the two over 3 per cent of the circle, and it averages to exactly zero while the direct term averages to something positive. Everything that happens slowly in a planetary system comes from that asymmetry: the part that survives averaging is not the part that dominates the instantaneous force.

The series that is subtracted

The two-body problem is solved, so nobody solves it twice. Every planetary theory since Newton begins by taking that solution away and asking what is left — and what is left is an infinite series whose terms are stacked in a hierarchy that makes the first half-dozen of them enough.

part 3 · orbits

Jupiter and Saturn trading a degree of longitude over nine centuries. Mean longitudes of Jupiter and Saturn in arcminutes, with their straight trends taken out and a 40-year smoothing applied, from a planar integration of the Sun, Jupiter and Saturn from their elements at the year 2000. Jupiter swings ±20′ and Saturn ±50′, in opposite senses, with a period of about 870 years. Comparing Ptolemy's positions with those of his own day, Kepler could not fit one mean motion to both, and by 1676 Halley reported Jupiter gaining and Saturn losing. Before 1785 that looked as if the Solar System were unstable; the curve shows that it is a very long period.

A term too small to keep that a denominator kept

Every planetary theory is an expansion, and the order of its terms is meant to say which can be dropped. In the 1770s the largest unexplained motion in the Solar System, Jupiter speeding up and Saturn slowing down by a degree over centuries, came from a term that the ordering called negligible. Five of Saturn's orbits nearly match two of Jupiter's, and that one near-match divides a tiny third-order coefficient twice by a tiny frequency. An integration of the three bodies shows the swap directly, and shows what Laplace saw: it reverses every nine centuries.

part 4 · orbits

The first-order perigee rate is right for a slower Sun, and half wrong for ours. The rate at which an orbit's perigee turns, in turns per orbit, against the ratio m of the Sun's motion to the orbit's, from separate integrations of a satellite with the Moon's eccentricity under a solar tide on a circular path. The first-order rate (3/4)m² gives 0.0042 at the Moon's m of 0.0748, a perigee period of 18 years. The integrations give 0.0085, 8.9 years; the Moon's perigee takes 8.85. At m = 0.01 the integration gives 8.3e-5 against a first-order 7.5e-5. So the disagreement is not a matter of neglected forces: it depends on how large m is.

A perigee that moved twice as fast as Newton

The Sun's tide on the Moon should turn the line of its perigee once in eighteen years. It takes less than nine. For sixty years that was the outstanding failure of the law of gravitation, and in 1747 Clairaut told the Paris Academy that the inverse-square law would have to be amended. Two years later he withdrew the remark. The law was right and the calculation was unfinished: the series for the perigee has a small parameter that is not small, and its second term is seventy per cent of its first.

part 5 · orbits