Astronomy & Celestial Mechanics
the sky, computed
Series

Chaos — the series

2 essays on one idea, from the one that introduces it to the one that assumes the rest.
  1. 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.

    part 1 · gravitation
  2. The last curve across the cylinder, before and after it breaks. Surfaces of section for the standard map at K = 0.5, 0.9716, 1.4, each 22 trajectories iterated 41 times from a column of starting values. Nothing here is placed: every dot is an iterate. At K = 0.5 the islands are separate and the space between them is filled with curves that run all the way round in θ — a trajectory cannot get from one island to the next, and one started beside the hyperbolic fixed point wanders 2.91 in p — 0.46 of a cylinder — and stops. At K = 0.9716, which is Greene's threshold to four figures, the last of those curves is on the point of going and the same trajectory still only reaches 4.84, which is 0.77 of a cylinder. At K = 1.4 it covers 6.1 cylinders: the barrier is gone and there is nothing left to stop it. That transition is what a chaotic zone in the asteroid belt is, drawn without any asteroids.

    Where the chaos comes from

    A Lyapunov time says how long a prediction lasts. It does not say what destroyed it. The mechanism is two resonances whose libration widths overlap, and the transition can be watched happening on a surface of section as one number is turned up.

    part 2 · gravitation

All series