Astronomy & Celestial Mechanics
the sky, computed
Concept

Random walk — where it appears

Accumulation of independent steps of random sign, so that the total grows as the square root of their number. An unbiased numerical error grows this way; a biased one grows linearly in the number of steps and quadratically in time.

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

Which gradient is steeper, and where. The temperature gradient radiation would need in order to carry the whole luminosity, against the gradient a rising blob of gas follows, through an n = 3 polytrope with a Kramers opacity and an energy generation going as ρT^4. Wherever the first exceeds the second the layer is unstable and convects. This model puts the outer boundary at 89 per cent of the radius; helioseismology measures the base of the Sun's convection zone at 71.3 per cent, and the gap between the two numbers is the price of a polytrope with one opacity law rather than a solar model. The curve also rises above the adiabatic value inside 4 per cent of the radius, which is a convective core — it appears when the energy generation is concentrated enough, and it is the structure a star burning by the CNO cycle has.

The part of a star that boils

Energy leaves a star's interior either by radiation or by bulk motion, and which one happens is settled by comparing two temperature gradients. The comparison decides that the Sun's outer third churns and its core does not, and that a massive star does exactly the opposite.

stars · Energy transport A step size below which a smaller step is worse. The error left in a long integration against the step size, for methods of three different orders, with both contributions drawn. The falling lines are truncation error, whose slope on these axes is exactly the order of the method. The rising line is round-off, identical for all three because it is a property of the arithmetic and not of the algorithm: every operation loses a few bits, the losses are independent, and they accumulate as the square root of the number of steps — which is why its slope is −1/2 and why it rises as the step shrinks. Each method's total has a minimum, at a step of 1.0e-6, 5.0e-6, 9.7e-4 for orders 1, 2, 4. Below that minimum every halving of the step costs time and makes the answer worse. That is the practical reason a solar-system integration is not run at an arbitrarily fine step, and it is a reason with nothing to do with computer time.

An error that grows like a random walk

A long integration accumulates two errors with opposite habits. One falls when the step is made smaller and grows in proportion to the time; the other grows when the step is made smaller and accumulates as a square root. Which of the two dominates decides whether a billion-year integration means anything.

gravitation · Numerical integration The spike, and where one passage through the planets sends it. The distribution of the reciprocal semi-major axis 1/a of long-period comets, in units of 10⁻⁴ AU⁻¹ and each curve scaled to its own peak. The narrow shaded peak is the Oort spike: comets arriving from the Oort cloud for the first time, with 1/a of about 4·10⁻⁵ AU⁻¹ and a spread of a few times 10⁻⁵. The broad curve is where those same comets go after one passage through the planetary region, whose perturbations change 1/a by a typical 5·10⁻⁴ AU⁻¹ — ten times the spike's width. 47 per cent end with negative 1/a, on hyperbolic orbits that leave the Solar System; only 4 per cent remain inside the spike's width. A comet in the spike is therefore on its first passage, almost by definition: the planets cannot put a comet back into the spike once it has been through them.

A spike only first-time visitors can be in

Long-period comets arriving from the Oort cloud have orbital energies packed into a spike a few hundred-thousandths of an inverse astronomical unit wide. One passage through the planets changes that energy by ten times as much, so no comet can be in the spike twice — and the spike's position, near thirty thousand astronomical units, is where the Galaxy's own tide is just strong enough to lift a comet's perihelion over Jupiter and Saturn in a single orbit.

orbits · Hyperbolic orbits One history of mass loss, three orbits, three outcomes. The eccentricity of three planets' orbits against the fraction of the star's mass lost so far, as the star sheds 45 per cent of its mass in 6 equal pulses at irregular intervals. Each planet starts with e = 0.05. The pulses are the same events in the star's own time, and the planets differ only in their periods, so each pulse lasts 3, 0.3, 0.03 orbital periods respectively. For the innermost planet, to which every pulse is slow, the eccentricity ends at 0.050: the orbit swells and keeps its shape. For the outermost, to which every pulse is sudden, each pulse kicks the eccentricity by an amount set by where the planet happened to be, and it ends at 0.298. The middle planet ends at 0.090. The same total loss leaves three different shapes, and only one of them is predicted by either limit alone.

The size forgets how the mass was lost, and the shape does not

A dying star loses half its mass, and a planet round it keeps its orbit's shape if the loss is slow and has it scrambled if the loss is sudden. Real stars lose mass in pulses — each one sudden or slow depending on the planet — and what survives a train of pulses is neither answer: the orbit's size settles to the slow value quickly, while its eccentricity random-walks and keeps a record of the history.

orbits · Effective potential Spun apart, or run out of ice: two clocks against size. Numbers of orbits against nucleus radius, both axes logarithmic, for a comet on 41P's orbit losing mass at 41P's rate per unit area. Solid: the orbits until the nucleus has lost all its mass to sublimation, growing in proportion to the radius. Dashed: the orbits until outgassing torques, changing sense at random from one orbit to the next, spin it past the rate at which a strengthless body of density 500 kg/m³ comes apart (a period of 4.7 hours) — growing as the fourth power of the radius, because each orbit's change of spin goes as the inverse square and a random walk needs the square of the number of steps. With a torque efficiency of 0.01, about what other comets' measured spin changes suggest, the two clocks cross at a radius of 1.2 km; with 41P's own 0.044, at 3.3 km. Smaller nuclei are spun apart before they run out of ice; larger ones run out first. The crossing moves as the two-thirds power of the efficiency.

The jet that pushes a comet also turns it

The gas a comet sheds near the Sun pushes its orbit, and because the vents are not arranged symmetrically about its spin axis, it also twists the nucleus. In 2017 one small comet's day more than doubled in two months. The twist per orbit grows as the inverse square of the nucleus's radius while its supply of ice grows only as the radius, so below about a kilometre a comet is spun apart before it runs out of ice — and the shortage of small comet nuclei may be the result.

orbits · Non-gravitational forces How much a dying star stirs its companion's orbit, against their separation. The rms change in orbital eccentricity of a 0.8 M☉ companion when a 2 M☉ star becomes a 0.6 M☉ white dwarf, against the binary's separation before the loss, both logarithmic. Blue: the mass loss alone — twelve thermal-pulse episodes 500 years wide and a superwind 10,000 years wide, each adding a kick of random sign whose size falls off exponentially once it lasts longer than about a sixth of an orbit. Below 100 au the orbit is slow to every episode and the eccentricity is untouched. The pulses switch on near 229 au, and by 3,000 au the change is 0.26 and the superwind too has become sudden. Dashed: the change a velocity kick to the white dwarf of 0.05 or 0.2 km/s would make on its own, which grows as the square root of the separation. A 0.05 km/s kick outweighs the pulses out to 251 au; a 0.2 km/s kick outweighs the pulses at every separation drawn.

A kick that writes over the record of the pulses

A star that becomes a white dwarf loses half its mass in pulses, and a companion far enough away feels each pulse as a sudden jolt to its orbit's shape. The wide binaries of white dwarfs should therefore carry the pulses' record in their eccentricities. They would, except that the white dwarf itself is born moving at a few tenths of a kilometre a second, and a kick that small writes more eccentricity than all the pulses together.

orbits · Effective potential A tilt that wanders up, one passing body at a time. The inclination of the Moon's orbit to the plane it started in, against the number of lunar-mass bodies that have passed the Earth–Moon pair at 5 km/s within 300 Earth radii, with the Moon at 30 Earth radii, for four different sequences of passages. Each passage tilts the orbit by a fraction of a degree in a random direction, so the inclination performs a random walk and grows as the square root of the number of passages: after 1,000 the four runs stand at 3.8, 1.8, 3.9, 0.2 degrees. The dashed line is the Moon's present inclination, 5.145 degrees. Whether a given history reaches it is a matter of chance; how far a typical one gets is set by the passages' number, mass and speed.

A tilt built from a thousand near misses

The Moon's orbit is tilted five degrees to the ecliptic, and slow tidal recession cannot have done it. One candidate is the leftover debris of planet formation: bodies the size of the Moon, passing the young Earth–Moon pair, each giving its orbit a kick of a fifth of a degree in a random direction. A random walk of a few thousand such kicks reaches five degrees — and the bodies that would supply them are the same population whose impacts left the Earth's mantle its gold.

orbits · Angular momentum

Named alongside it

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

Adiabatic invariantAsymptotic giant branchEccentricityImpulsive approximationMass lossOsculating elementsWhite dwarfAdiabatic gradientAngular momentumBrouwers lawCNO cycleComet nucleus

All concepts