Astronomy & Celestial Mechanics
the sky, computed
Concept

Phase-space — where it appears

The space whose points are complete states — every position and every momentum at once — in which a trajectory never crosses itself. A surface of section through it turns a trajectory into a map, and whether the resulting points fall on a curve or scatter is what distinguishes regular motion from chaos.

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

Energy error over 240 revolutions, at one step size. The relative error in total energy against revolution number, for three integrators run on the same Kepler orbit at e = 0.5 with the same step of 200 per revolution. The exact energy is a constant, so every curve here is the method rather than the problem. Euler climbs steadily: its energy at the end is 106.5% wrong, and the orbit it draws has spiralled outwards. Runge–Kutta 4 begins 4.7e+2 times more accurate than leapfrog and ends at 2.22e-4, having grown by a factor of 10 across the run: the error is SECULAR. Leapfrog oscillates inside a band and stays there — worst error 2.62e-3, and the second half of the run is no worse than the first, which is measured here rather than claimed. That is the property that decides whether a five-billion-year integration means anything, and it is not accuracy: a symplectic method is the exact solution of a Hamiltonian a step-size away from the intended one, so its energy cannot wander, while a more accurate non-symplectic method has no such constraint and eventually wanders further.

Wrong about where, and right about how much

Runge–Kutta is the more accurate method and loses energy steadily; leapfrog is cruder and its energy error never leaves a band. Over five billion years only one of those properties survives — and neither method knows where the planet is.

gravitation · Numerical integration 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 stream that is not the orbit it came from. 456 stars released in pairs from the two saddles of a 10⁵ solar-mass cluster over 4.0 billion years, integrated in a halo whose circular speed is 220 kilometres a second, drawn with the progenitor's own orbit. The cluster runs between 10 and 25 kiloparsecs and the orbit is the thin closed-looking curve; the stars are everything else. The point of the figure is the discrepancy. Stars leaving through the inner saddle are on slightly smaller orbits, turning round at a median of 24.77 kiloparsecs rather than the progenitor's 25.00, and therefore running ahead; stars leaving through the outer saddle reach 25.22 and fall behind. The whole spread is 1.8 per cent of the apocentre, which is the number worth carrying: an offset far too small to see in this drawing builds the entire stream, because it acts for four billion years. The two arms are therefore not merely displaced along the orbit, they are on different orbits, and the track a survey measures is a family of them rather than any single one. Fitting a Galactic potential by demanding that a stream lie along an orbit is wrong by exactly this much, and the size of the error grows with the mass of the progenitor, because the mass is what sets the distance between the two doors.

A stream is not the orbit it came from

A cluster torn apart by a galaxy leaves a thin trail of stars across the sky, and the obvious thing to do with it is fit an orbit. That is wrong by a knowable amount, because stars leave through two doors with a small energy offset and end up on a family of orbits rather than on one.

galaxies · Stellar streams A ceiling the outer orbit imposes. The greatest eccentricity a cycle can drive an inner orbit to, against the mutual inclination it started at, for four ratios of the inner orbit's angular momentum to the outer one's. The bottom curve is the test-particle case and it is the familiar closed form, √(1 − (5/3)cos²i₀), reaching one at ninety degrees — recovered here by a solver that knows nothing about it, which is the check. Every other curve is lower. The reason is conservation: the cycle works by moving angular momentum out of the inner orbit, and it has nowhere to put it but the outer one. A test particle can give up all of its angular momentum because the outer orbit is infinite by comparison; an inner binary carrying 1 times as much as its companion cannot, because tilting the outer orbit by the required amount changes the mutual inclination and shuts the cycle down. At ninety degrees the ceiling falls from 1.000 to 0.866. The pericentre a mechanism can deliver is bounded by what the perturber can absorb, and the test-particle formula is an upper bound rather than a prediction.

A ceiling the outer orbit imposes

The closed form for how far a cycle can drive an eccentricity assumes the inner orbit is a test particle, so the angular momentum it gives up goes nowhere. A real inner binary of comparable mass has to be paid for: the outer orbit absorbs what the inner one loses, tilts, and changes the very inclination the cycle depends on.

orbits · Kozai–Lidov A cold slab collapsing into a spiral. 800 self-gravitating sheets in phase space — position across, velocity up — at times 0, 3, 8, 120 in units where a sheet released from the edge reaches the centre in about 1.4. They start as a cold rectangle, uniform in position, with 10 per cent of the kinetic energy equilibrium needs. The slab falls in, overshoots, and the sheets, each orbiting at a period set by its own energy, wind into a spiral whose arms grow thinner and closer until they can no longer be told apart. By the last panel the system is in a steady state: a dense core and a halo of sheets that were flung out in the first collapse. Nothing in the dynamics is an encounter; every change of a sheet's energy comes from the potential changing in time.

A relaxation that stops before it is finished

A cold, collapsing system of stars reaches virial equilibrium in a few crossing times without any star meeting another, because the potential changing in time reshuffles their energies. Lynden-Bell's theory predicted the state it should reach — the most probable one. Computed exactly, in a slab of self-gravitating sheets where every force is known, the collapse reaches equilibrium and not that state: the fluctuations that drive the mixing are what the mixing destroys, and they die before the job is done, leaving a profile and an ordering of energies that remember how the collapse began.

gravitation · Virial theorem

Named alongside it

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

Dark matter haloNumerical integrationOrbital energyPerturbationsThree-body problemAngular momentumChaosCollisionless systemConditioningConservation lawsConserved quantityCritical inclination

All concepts