Astronomy & Celestial Mechanics
the sky, computed
Generator

The three-body generator

Three bodies, integrated
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 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.

8 essays call three-body. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

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. Gravitation

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.

 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. Gravitation

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.

 An average that arrives, and a snapshot that never does. 2T/|U| against time for Burrau's problem — three masses released from rest, and never repeating, in two forms: the instantaneous ratio, and the ratio of the running time-averages. The instantaneous one runs between 0.03 and 1.97 — the system is a long way from equilibrium at almost every moment, because the bodies are alternately falling together and flying apart. The averaged one settles: after 50 time units it is 1.0142, against the exact 1 the virial theorem requires of any bound system. This system is not periodic and is not even permanently bound — Burrau's problem ejects its lightest body — so the average is drifting rather than converged, and that is the theorem's own condition made visible: it is a statement about bound systems and says nothing about anything that is leaving. What the figure cannot show is the error: a cluster observed once gives 2T/|U| with a scatter of this size, and what makes its mass believable is not the measurement but the relaxation. Gravitation

An average that weighs what cannot be watched

A cluster's mass can be had from its speeds alone — no orbit followed, no period observed, no distance to any single star. The theorem that allows it is an average over time, which is exactly what a photograph is not.

 A cold sphere collapsing onto the virial theorem. The virial ratio 2K/|W| — twice the kinetic energy over the magnitude of the potential energy — counted over the bound shells, of 400 concentric shells released almost at rest (2K/|W| = 0.1) from a sphere whose density falls as one over the radius, against time in units of the sphere's free-fall time. It starts near zero; overshoots to 1.7 as the collapse reaches the centre, the inner shells first; rebounds and oscillates; and settles to 1.03, the virial value, staying within ten per cent of it from about 4.9 free-fall times on. The shells never collide with anything; each feels only the mass inside it. What damps the oscillation is phase mixing — shells of different energy oscillate with different periods and drift out of step — and what redistributes the energy is the potential changing while they move through it. A system of stars reaches equilibrium this way in a few dynamical times, where encounters between the stars would take far longer than the age of the universe. Gravitation

An equilibrium reached without a single encounter

Every dynamical mass of a galaxy assumes its stars obey the virial theorem, and no two of those stars have ever exchanged energy. A collapse computed shell by shell shows how a system gets there anyway — in a few free-fall times, at half the size it fell from, by a potential that changes while its orbits are still moving through it.

 Which hierarchical triples survive, against the outer orbit's pericentre and eccentricity. Each cell is one integration of a triple: an equal-mass inner binary on a circular orbit of unit separation, and a third body of 1 times the binary's mass on an outer orbit of the given eccentricity with its pericentre at the given distance from the binary's centre of mass, all in one plane and going the same way round. Cells coloured as survivors kept the binary intact and ejected nothing for 30 outer orbits — 37 of 50; the others broke up, shaded darker the longer they lasted. The curve is the Mardling–Aarseth stability criterion, a pericentre of 2.8 inner separations times ((1 + q)(1 + e)/√(1 − e))^0.4, derived from where the resonances between the outer orbit and the inner one begin to overlap, and meant for survival over very many orbits. Every triple outside it survived. Inside it, 17 of 30 also survived this run; the outermost break-up is at a pericentre of 3.0 for a circular outer orbit and 2.6 at e = 0.6. The boundary drawn from a finite integration is therefore a lifetime rather than a line — triples between the break-ups and the curve survived a run of this length and are not thereby shown to be stable; the criterion is a statement about how the resonances are arranged, and a finite integration can only fail to refute it. Gravitation

A third body kept at a safe distance

Most triple stars last for billions of years because they are hierarchical — a close pair and a distant third. How distant is distant enough has a formula, the formula comes from overlapping resonances, and an integration of any finite length can only fail to refute it.

 Which of three bodies is thrown out. Of 60 integrations of three bodies of mass 1, 2, 3 released from rest in random triangles, 57 broke up within the time integrated, each leaving one body unbound and a binary behind. The escaper was the body of mass 1 in 46 runs (81 per cent), 2 in 6 runs (11 per cent), 3 in 5 runs (9 per cent). The bars beside them are the shares the statistical theory of break-up predicts to leading order, in proportion to the escaper's mass to the power −3: 86, 11, 3 per cent. Every run is deterministic and none of them could be predicted without integrating it, but the ensemble is: the lightest body leaves most often when the masses differ, because the energy of the binary left behind is shared between the escaper and the binary's recoil in inverse proportion to their masses, and the lightest body gets the most speed. Energy is conserved in every run to better than 3·10⁻⁴. Gravitation

The lightest body leaves

Three stars close together, with no hierarchy between them, cannot stay together. One of them is thrown out, which one is unpredictable in any single case, and across many cases it is almost always the lightest — a statistical law emerging from a deterministic chaos that cannot be integrated far enough to forecast.

 The figure-eight with equal masses, and with one mass changed. The paths of three bodies in the figure-eight choreography over 100 periods, with equal masses (left) and with the third body's mass changed by 0.0001 (right). With equal masses the three paths lie on one another: each body runs round the same figure-of-eight, a third of a period behind the next. With the mass changed the three bodies stay bound but leave the single curve, filling a broad band round where it was. The choreography is an exact solution only for exactly equal masses, and the island of stable orbits near it is narrow. Gravitation

A periodic orbit that the sky cannot keep

Three equal masses can chase one another for ever round a single figure of eight. The orbit exists, it is stable, and it will almost certainly never be seen — because its island of stability is a hundredth of a per cent wide in the masses, and nothing in nature makes three stars equal to that precision.

 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. Gravitation

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.

The whole library · All essays