Astronomy & Celestial Mechanics
the sky, computed
Concept

Resonance — where it appears

The condition in which two orbital periods stand in a ratio of small integers, so that perturbations accumulate rather than average away. It is what clears the Kirkwood gaps and what locks the Galilean moons, and the same mechanism can be protective or destructive depending on the phase it holds.

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

A separation of 10⁻¹⁰ becomes 0.0285. The distance between two copies of Burrau's problem, started 10⁻¹⁰ apart in one coordinate of one body and then stepped in lockstep — one loop, one step size, both states advanced by it — for 60 time units. The vertical axis is logarithmic, so the straight stretch is exponential growth, and its slope is the Lyapunov exponent: 0.3413 per time unit, fitted over the 604 samples lying between ten times the starting offset and a tenth of the system's own size. That is an e-folding every 2.93 time units, so the separation multiplies by ten every 6.75 and by 7.8·10⁸ across the whole run. It has not saturated within the drawn interval: the growth rate over the last tenth of the run is still 0.6002 per time unit, 176% of the fitted exponent, and the separation has reached only 0.90% of the system's own size. The exponent is what sets a prediction horizon, and this run is drawn short of it on purpose: the straight line is the measurement, and the flat part that follows is arithmetic about how far apart two bounded systems can get.

A prediction with an expiry date

The inner solar system's Lyapunov time is about five million years, so a centimetre of error becomes an orbit in a hundred million. The ephemeris dies while the system survives — because the elements stay bounded when the phase does not, and only one of those is what stability means.

gravitation · Chaos A transit that will not keep time. Transit times of a planet of 8 Earth masses on a 3-day orbit, minus the best straight line through them, from a three-body integration in which the second planet at 4.62 days is the only thing perturbing it. The residual swings by ±1.7 minutes and repeats over 58 days — the super-period 1/|3/P₂ − 2/P₁|, which is long precisely because the pair is near the 3:2 resonance and gets longer the nearer it is. The amplitude is proportional to the perturbing planet's mass, so fitting this curve weighs a planet that may never transit at all. The short jagged component is not noise: it is the chopping signal, one kick per conjunction at the 8.6-day synodic period, sampled once every 3 days by the transits themselves and so barely above the rate at which it can be followed. The integration conserved energy to 1.6e-10.

Planets found by a transit running late

A planet on a fixed orbit transits like a clock. A second planet pulling on it makes the clock run fast and slow by minutes — and fitting that wander weighs a planet that may never cross the star at all.

exoplanets · Transit timing Where a 2-armed pattern is allowed to exist in the Milky Way. Angular speed against radius, with the two combinations that matter for a pattern of 2 arms. A wave turning at 25 km/s per kpc corotates with the stars at 8.7 kpc, has an inner Lindblad resonance at 1.9 kpc and an outer one at 13.8 kpc. Between the Lindblad radii the wave propagates; at them it is absorbed. Every one of these radii is found by bisection on the curves as drawn, not quoted — and none of them is a property of the pattern alone, because κ comes from the rotation curve, which comes from the mass.

Where a pattern is allowed to turn

A spiral wave has one pattern speed, and that single number decides where in the disc it can exist at all. The boundaries are set by the frequency at which a star rocks radially about its circular orbit — a frequency the rotation curve already contains.

galaxies · Spiral structure An edge where two torques balance, 21 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 12 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 20.8 kilometres out against a measured half-width of 21, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.

An edge is a balance, not a boundary

A ring of particles spreads, because collisions move angular momentum outward. Something has to push back, and the something is a small moon whose torque falls as the inverse cube of the gap. A flat curve and an inverse cube cross once, and the crossing is the sharp edge — which is why an edge exists at all rather than a gradient.

gravitation · Planetary rings 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 An edge where two torques balance, 21 kilometres from the shepherd. Two torques on the edge of a ring, against distance from a shepherding moon, both axes logarithmic and both scaled by the same combination of surface density, radius and orbital rate so that only their shapes are being compared. The flat line is the viscous torque, which comes from collisions between ring particles and does not care how far away anything is; it is drawn for a kinematic viscosity of 12 square centimetres a second, within the range ring seismology gives. The falling line is the moon's, summed over the first-order resonances that crowd together as the gap narrows, which makes it an inverse cube. A flat curve and an inverse cube cross once, and the crossing is where an edge can sit: closer in the moon wins and pushes the material back, further out viscosity wins and the ring spreads. For Daphnis and the Keeler gap the balance lands 20.8 kilometres out against a measured half-width of 21, which is agreement to well inside the uncertainty on the viscosity — and it is the only handle anybody has on that viscosity, since the quantity being inferred is the collision rate among particles a metre across, a billion kilometres away.

A torque that nearly cancels

A planet embedded in a gas disc pulls on the material inside its orbit and outside it, and the two torques are almost equal and opposite. What survives the subtraction is a per cent of either, and it is still enough to carry a planet from where it formed to its star in less time than the disc lasts.

exoplanets · Planet migration A resonance with the planet's rotation, and it sorts by size. Above: the radii at which a charged grain's orbital frequency is commensurate with the planet's spin, in planetary radii, with the measured edges of Jupiter's halo ring drawn over them. The synchronous radius is 2.24; the 3:2 resonance falls at 1.71 and the 2:1 at 1.41. The halo's outer boundary is at 1.71 and its inner extent near 1.4 — the ring ends where the resonances are, and both numbers were measured by a spacecraft camera with no reference to this arithmetic. Below: which grains care. The charge on a grain is proportional to its radius and the mass to the cube, so the force per unit mass goes exactly as the inverse square of the size, and the resonance grips sub-micron dust while leaving anything larger on a Keplerian orbit. That is why the halo is a cloud of fine dust puffed a thousand kilometres out of the ring plane while the coarse material stays flat: the same field acting on the same orbit sorts the material by size, which no gravitational resonance can do.

A resonance with the planet itself

Every other resonance in celestial mechanics is a commensurability between two orbits. A charged dust grain has a third clock available — the planet's rotation, which sweeps its magnetic field past the grain — and the commensurability with that selects by charge-to-mass ratio, which means by grain size.

orbits · Planetary rings Two longitudes a satellite falls towards, and 0.5 m/s a year to stay elsewhere. The along-track potential a geostationary satellite feels, against longitude. The Earth's equator is slightly elliptical — about seventy metres between its long and short axes — and that one harmonic of the gravity field gives the geostationary ring two minima and two maxima. A satellite left unattended drifts towards the nearer minimum at 75° or 255° east, overshoots, and librates about it with a period of a couple of years. The peak acceleration accumulates 0.5 metres a second of velocity change a year, and an operational east–west budget is a small multiple of that once the correction cycle is accounted for — a fixed cost of every commercial slot in the ring, paid forever. The two minima are the graveyard of the geostationary population: uncontrolled satellites accumulate there, which is why they are the most crowded longitudes in the belt and why an uncontrolled object is most likely to be found near one.

A satellite that drifts to one of two longitudes

The Earth's equator is elliptical by about seventy metres. That one harmonic of the gravity field gives geostationary orbit a potential with two minima, so an unattended satellite slides towards the nearer one and stays — and every operational slot in the ring is paid for with fuel, every year, forever.

spaceflight · Station-keeping One eccentric planet, drawn as two circular ones. Top: the velocity of a star with one companion on an orbit of eccentricity 0.25 (argument of periastron 0.9 radians), in units of its semi-amplitude K, over two orbits — and on top of it the sum of two circular orbits, one at the same period and one at half of it, fitted by taking the first two terms of the exact curve's Fourier series. The second orbit's amplitude comes out at 0.246 of the first, which is the eccentricity to first order. Bottom: the difference between the two, on a scale 12 times magnified; its largest excursion is 0.0889 K and its rms 0.0473 K, second order in the eccentricity. For a planet with K = 30 m/s that is 1.42 m/s rms — below the noise of most instruments — so an eccentric single planet and a pair of planets in a 2:1 period ratio are the same measurement unless the data are good enough to see a term of order e².

One eccentric planet, or two circular ones

The velocity curve of a star with one planet on an eccentric orbit is, to first order in the eccentricity, exactly the curve of two planets on circular orbits with periods in the ratio two to one. The difference between the two readings is a term in the eccentricity squared, and for the modest eccentricities most planets have, it is smaller than the noise.

exoplanets · Reflex velocity A doubling rule, and the planet that broke it. Each planet's real semi-major axis against the value the Titius–Bode rule 0.4 + 0.3 × 2ⁿ gives it, both in astronomical units on logarithmic axes; points on the diagonal are fitted exactly. The rule gives the Earth 1.0 by construction and fits Venus, Mars, Jupiter and Saturn within a few per cent. It needed the value n = 3 left empty until Ceres was found in its slot in 1801, and it predicted Uranus, discovered in 1781, within 2%. Neptune, discovered in 1846 by a calculation that assumed the rule, sits 22% inside its predicted place, and Pluto occupies the place Neptune should have had. Mercury is given n = −∞, which is a free choice made to fit it. The labels give each body's departure from the rule.

A spacing law that random orbits obey

The planets' distances roughly double from one to the next, and a rule written down in 1766 predicted the asteroid belt and Uranus before either was found. It failed on Neptune, and the reason it worked at all turns out to be that almost any set of orbits forbidden to crowd each other looks geometric on a logarithmic scale — as a test against thousands of random systems shows.

orbits · Harmonic law A chaotic orbit that waits, and then moves. The momentum of one trajectory of the standard map at K = 1.2, unwrapped and counted in whole turns of 2π, against the iteration number, over 39,892 iterations. The motion inside each band is chaotic from the first step — the trajectory's Lyapunov time is 3.7 iterations — but it lingers near one integer resonance for thousands of iterations before reaching the next: the waits were 3,098, 1,464, 9,461, 15,957, 6,058, 3,854 iterations. Each arrival at the next integer resonance, drawn dashed, is marked. Between the integer resonances lie the remnants of the broken invariant curves, which no longer block the trajectory but still hold it, and each passage is quick compared with the wait before it.

A barrier that leaks at a rate that can be computed

When the last invariant curve between two resonances breaks, it does not vanish. It becomes a cantorus — a curve full of gaps — and a chaotic orbit has to find the gaps to get through. The time that takes grows as the cube of how close the system is to the threshold, and it has nothing to do with the Lyapunov time, which is why an orbit can be hopelessly unpredictable and still stay put for the age of the solar system.

gravitation · Chaos The nutation lines and the core's resonance. The resonance factor |σ₀ / (σ − σ₀)|, σ₀ being the free core nutation's frequency, on a logarithmic scale, against frequency in the celestial frame in cycles per sidereal day, negative for retrograde motion, with the free core nutation at a period of 430.2 days retrograde. The vertical lines are the largest forced nutations the Sun and Moon drive: the 18.6-year, annual, semiannual and fortnightly terms, each prograde and retrograde, and precession at zero. Any part of the Earth's response that involves the core is multiplied by this factor. The retrograde annual term, 365 days, is 5.6 on this scale; the prograde annual term 0.46; the retrograde 18.6-year term 1.07. The annual retrograde nutation is the one the core's resonance distorts most, because a year is the forced period nearest the core's own.

A core that wobbles at its own period

The Earth's liquid core can turn about an axis slightly tilted from the mantle's, and the flattened boundary between them pushes it back. That gives the planet a second free mode of rotation: in space a slow retrograde nutation of about 430 days, on the ground a wobble just faster than once a day. Its period is a measurement of the shape of a surface under 2,900 kilometres of rock, and it is read twice — in the nodding of the Earth's axis against the quasars, and in the rise and fall of the solid ground under the daily tides.

sky · Precession Heliocentric latitude bought three ways, for the same launch. The highest inclination to the ecliptic a spacecraft can reach, against the hyperbolic excess speed it leaves the Earth with. Launched directly into an inclined orbit at 1 AU (dashed), it pays 2v⊕ sin(i/2) — 15.4° for 8 km/s and 23.2° for 12. Launched on a transfer to Mars whose plane is tilted — the wall across the porkchop plot, where transfers near half a turn are forced steeply out of the plane — it reaches latitude i half-way, for 13.7° at 8 km/s and 21.5° at 12: the same list price, starting from the cost of getting to Mars at all. Launched to Jupiter, it can turn there: a flyby rotates its Jupiter-relative velocity, and the steepest heliocentric orbit a flyby at 6 Jupiter radii can leave it on climbs to 90° for a 12 km/s launch, reaching the 79° of Ulysses's orbit (dotted) from about 10.6 km/s — close to the launch speed that mission actually used. Inclination bought near the Earth costs the Earth's orbital speed; bought at Jupiter it costs Jupiter's, less than half as much, and the flyby pays.

The wall sells inclination at the dearest price

Transfers to Mars that go almost exactly half-way round the Sun are forced into planes tilted steeply out of the ecliptic, which is why they form a wall of expensive dates across every launch chart. A mission that wanted to leave the ecliptic could use the wall on purpose. It would be buying inclination at the Earth's orbital speed — the dearest place to buy it — and the routes that reach the Sun's poles have always gone somewhere slower to make the turn.

spaceflight · Launch windows The tide's spectrum, with every frequency set by the sky. The equilibrium tide-raising potential broken into its largest constituents, each a sinusoid whose frequency is a sum of whole-number multiples of six astronomical rates — the Earth's turning relative to the Moon, the Moon's and the Sun's motions round the sky, and the slow turning of the lunar perigee and node — and whose height is its equilibrium amplitude relative to M₂. They fall in three bands: near two cycles a day (M₂ at 1.9323, S₂ at exactly 2, N₂, K₂), near one (K₁, O₁, P₁, Q₁), and near zero (fortnightly, monthly, semi-annual). Every frequency is known to many decimal places from astronomy; what a harbour does to each, its own amplitude and phase, is not, and is the only thing a tide table has to measure.

A tide table is astronomy with the harbour fitted

Nobody can compute the tide at a harbour from the Moon and the Sun, because the ocean's response depends on coastlines and depths nobody models well enough. Tides are predicted to centimetres years ahead anyway, because the frequencies are pure astronomy and only the response to each one has to be measured. The method needs a record long enough to tell neighbouring frequencies apart, and a correction for a cycle of 18.6 years.

gravitation · Tides A cluster sinking by Chandrasekhar's formula alone. The orbital radius of a 2.0 × 10⁵ solar-mass globular cluster released on a circular orbit at 1 kpc in the two Fornax halo models, decaying under Chandrasekhar's dynamical friction with a velocity dispersion of 11.7 km/s and a Coulomb logarithm of 3. In the cusp it reaches the centre in 6.7 billion years. In the core the formula gives no arrival within 12 billion years, crossing into the core at 2.8. The formula already slows the cluster inside the core, because the circular speed there drops below the stars' velocity dispersion and few stars move slower than the cluster; but it still lets the cluster creep inward without limit. N-body experiments of clusters sinking into cored haloes find something stronger: the decay stops altogether near the core, which the local formula cannot produce at any setting of its logarithm.

A drag that stops where every orbit keeps time

Chandrasekhar's formula reads the density and the velocities where a sinking body is, and nothing about the orbits those stars are on. Inside a constant-density core every orbit has the same frequency, the differential motion that drags a body inward disappears, and the sinking stops — which is why the five globular clusters of the Fornax dwarf galaxy are still where they are.

gravitation · Dynamical friction

Named alongside it

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

Angular momentumChaosFourier seriesGravity assistLindblad resonanceLyapunov timeOrbital elementsSurface densityThree-body problemTorqueViscosityAngular momentum transport

All concepts