Astronomy & Celestial Mechanics
the sky, computed
Concept

Binomial statistics — where it appears

The statistics of counting how many of a fixed number of independent trials succeed. A probability estimated this way has a relative error of about the square root of (1 − p)/(pN), so rare events need proportionally more trials.

Named by 2 essays across one field — each of them below, with the objects they name alongside it.

Also named here as ensemble — the same set of essays touches all of them, so they are one junction rather than several.

Sixty futures that began a trillionth apart. The momentum of the standard map at K = 3, a kicked rotor whose phase space is almost entirely chaotic, for 60 members started 10⁻¹² apart in angle, against iteration. For 41 iterations they are one trajectory; then the separation, growing by a factor of about 2.0 each step, reaches the size of the system and they become independent random walks. Three of the 60 cross the threshold at |p| = 90 within 400 iterations (dark), the stand-in for an orbit that goes unstable. Which members cross is decided by differences far below anything that could be measured; how many cross is a property of the dynamics, and it is the only thing the ensemble can be asked for.

A probability for a system with one future

The Solar System will do one thing, and an integration of it is wrong about which thing within a hundred million years. What long integrations report instead is a probability — one per cent that Mercury's orbit goes unstable before the Sun dies — and a probability of a single event is a strange object. It is a count over an ensemble of futures that differ by less than anything measurable, which makes sense only because chaos makes those futures independent; and how many of them are needed is a binomial question with an unforgiving answer for rare events.

gravitation · Numerical integration A rare crossing reached in stages. Multilevel splitting on a kicked rotor (the standard map at K = 10), aimed at the event |p| > 450 within 400 iterations. 16 members are run until each reaches the first level or runs out of time; those that reach it are copied back up to 16, each copy restarted from its parent's crossing point with its angle nudged by a thousandth of a radian, and run on toward the next. The levels (dashed) are not fixed values of |p| but curves that rise as time runs out, so a member counts as closer to the event only if it has the time left to finish. 4 levels; the fractions passing them are 0.13, 0.06, 0.63, 0.88, and their product, 0.0043, is the estimate. Each stage is a common event; the rare one is never waited for.

A rare event reached by cloning the near misses

Counting how often an ensemble of integrations goes unstable works for an event of one in a hundred and becomes hopeless for one in a million, because a ten per cent answer needs a hundred events and so a hundred million futures. Splitting gets there by a different route: wait for the members that come closest, copy them, and run the copies on. Each stage measures a common event and the rare one is a product of common ones — provided the copies forget their parents and the stages know how much time is left.

gravitation · Numerical integration

Named alongside it

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

ChaosEnsembleLyapunov exponentNumerical integrationSolar system stabilityStandard mapRare event samplingRoundoff error

All concepts