Astronomy & Celestial Mechanics
the sky, computed
Concept

Osculating elements — where it appears

The elements of the unperturbed orbit a body would follow if all perturbations ceased at that instant. They change continuously under perturbation, so an element quoted for a real body is a snapshot with a date rather than a property.

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

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.

orbits · Perturbations The angle that has nowhere to be measured from. An eccentricity vector carried round a circle of radius 0.034 centred at 0.031 — which is what a secular perturbation does to one, a forced eccentricity with a free one turning about it. Below, the two components e cos ϖ and e sin ϖ, which are smooth, bounded and perfectly ordinary throughout. Above, the longitude of pericentre read off them, which is not: as the eccentricity passes its minimum of 0.0030 the pericentre sweeps through most of a circle, at up to 4080° per unit time against the 12° the free vector itself turns in the same interval. Nothing has happened to the orbit. The pericentre is a place on the orbit, and a nearly circular orbit does not have one — so ω, and Ω with it at zero inclination, are angles measured from a feature that is not there. The equinoctial elements are the pair drawn below, and a propagator written in them steps through this instant without noticing it.

The elements that stop existing

An orbit needs six numbers, and three of the usual six are angles measured from features a perfectly ordinary orbit may not have. At zero eccentricity there is no pericentre to measure from, and the arithmetic knows it.

orbits · Orbital elements 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.

orbits · Perturbations Earth's eccentricity is a sum of 8 sinusoids. The eccentricity of Earth over 800 thousand years, from the Laplace–Lagrange solution for all eight planets — the secular matrix built from the JPL masses and semi-major axes, symmetrised, and diagonalised by Jacobi rotations. It runs between 0.0035 and 0.0436, and it has no period, because it is a sum of 8 incommensurable frequencies. The two largest contributions to this planet are the modes at 3.73 and 7.33 arcseconds per year, drawn as the flat lines: those are constants, and everything moving in the figure is their beat. The check is the quadratic form ½ΣΛe², which a symmetric secular matrix conserves exactly and which drifts by 6.7e-16 across the whole interval — computed from the same curves the figure draws and from nothing the eigenvalues were fitted to. The true angular momentum deficit, Σ Λ(1 − √(1−e²)), drifts by 7.1e-4, and that difference is not an error either: the two agree only to fourth order in e, and Mercury at 0.206 supplies almost all of the gap.

No planet has an eccentricity of its own

Strip the short-period terms out of the planetary equations and what is left is a linear system. Its eigenvectors are modes of the whole solar system, and the number a catalogue quotes for a planet's eccentricity turns out to be a reading of a clock rather than a property of the planet.

orbits · Secular theory The contours a planetary encounter cannot cross. Contours of the Tisserand parameter with respect to Jupiter in the plane of semi-major axis and eccentricity, drawn for a coplanar orbit. An encounter with Jupiter moves a comet along one of these curves and never across one, because T is what the encounter conserves. The heavy contours are at T = 3 and T = 2, and they are the boundaries the comet families are defined by: T > 3 means no encounter is possible at all, since v∞²/v_J² = 3 − T and a negative squared speed is not a trajectory; 2 < T < 3 is the Jupiter family; below 2 the approach speed exceeds Jupiter's own orbital speed and the orbits are the nearly isotropic ones. The shaded boundary is the crossing condition — an orbit whose pericentre is outside Jupiter's, or whose apocentre is inside it, never meets the planet whatever its T. The five comets are placed at their JPL elements and labelled with the T the literature quotes; all five with inclination sit off the coplanar contours by exactly the cos i in the definition, which is why 1P/Halley's is negative — a retrograde orbit meets Jupiter at nearly twice Jupiter's speed.

The number that survives the encounter

A comet that passes Jupiter comes away with every orbital element changed. One combination of them is not changed, and it is enough to recognise the comet afterwards, to sort the comet families, and to say where a spacecraft can and cannot go.

gravitation · Tisserand parameter An invariant that moves by 8.0e-3 once the planet's orbit is real. The Tisserand parameter of a comet on an orbit of semi-major axis 5 and eccentricity 0.8, followed through a close passage of Jupiter, integrated twice. The lower trace has Jupiter on a perfect circle, which is the problem the parameter is an exact constant of: it survives the encounter having moved by 1.0e-4, which is the integrator's own error and not a physical change, and the spike at the moment of closest approach is the osculating elements being briefly meaningless while the comet is inside Jupiter's sphere of influence rather than the constant failing. The upper trace is the same encounter with Jupiter on its real orbit, eccentricity 0.0489. The parameter comes out changed by 8.0e-3, 79 times as much, because the Jacobi constant exists only when the rotating frame is uniformly rotating and a planet on an ellipse does not provide one. That number is small and it is not negligible: comet families are separated by boundaries in this parameter placed to two decimal places, and a comet that drifts across one over several encounters has changed class without anything having happened to it that a single encounter could account for.

An invariant that is only almost one

The Tisserand parameter survives a close encounter with Jupiter exactly, and comet families are separated by boundaries in it drawn to two decimal places. The exactness holds for a Jupiter on a circle. Jupiter's eccentricity is 0.0489, and integrating the same encounter twice shows what that costs.

orbits · Tisserand parameter 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.

orbits · Perturbations An ephemeris fitted to 120 days, 43 minutes wrong within a year against a band of ±5.1. Transit times of a 6 Earth-mass planet on a 10-day orbit, perturbed by a 14 Earth-mass planet at 15.24 days, integrated for 1460 days and compared with a straight-line ephemeris fitted only to the transits in the first 120 days — the shaded window. Inside the window the line fits to 2.1 minutes. Outside it the pair's 317-day super-period carries the transits away from the line, and within a year of the window closing the prediction is 42.8 minutes early of the observed transit, 347 days after the last fitted one. The narrow band is the formal three-sigma uncertainty of the same line for a timing precision of 0.5 minutes per transit, which at that date is ±5.1 minutes: the error is 8.4 times the band. A statistical uncertainty assumes the residuals are noise, and these are a signal, so the band describes a planet that does not exist.

A forecast that fails on a schedule

A transiting planet perturbed near a resonance keeps a clock that wanders, and a straight-line ephemeris fitted to part of the wander predicts the next transit with a confidence the wander does not deserve. The error is not noise and does not average down; it grows on the pair's super-period, it is many times the formal uncertainty within a year, and how soon it appears depends on which stretch of the wander happened to be observed. When a model that includes the known perturber still fails, the failure has a period, and the period is a planet.

exoplanets · Transit-timing

Named alongside it

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

Disturbing functionSecular perturbationForced eccentricityGravity assistJacobi constantLaplace coefficientNode regressionOrbital elementsRestricted three-body problemSecular variationTisserand parameterAngular momentum deficit

All concepts