Astronomy & Celestial Mechanics
the sky, computed
The collection

Every essay — page 26

One idea per essay, ordered so that the earlier ones set up the later ones — but nothing here depends on being read in sequence. Essays 501–520 of 534.

Exoplanets

Planets nobody has seen, weighed and measured from a dip, a wobble and a delay.

An eccentricity of 0.4 swings the surface by 126 K or by 0.1, depending on the length of the year. The peak-to-trough swing in surface temperature over one orbit, against orbital period, for a planet with eccentricity 0.4 receiving on average the flux the Earth does, for surface layers of 1, 10, 50 metres of water. The dashed line is the 126 K the surface would swing through with no heat capacity. Every curve rises from near zero at short periods, where the orbit is over before the layer can respond and the planet feels only the average flux, towards the full swing at long periods, where every part of the orbit lasts long enough to be felt in full. The crossover is where the orbital period is comparable with the layer's thermal time. The vertical marks are the orbital periods of the Earth-flux orbit round stars of 0.1 M☉ (7 days), 0.5 M☉ (0.2 years), 1 M☉ (0.9 years). At a fixed eccentricity a planet in the habitable zone of a small star is thermally averaging almost regardless of how much water it has, and one round a Sun-like star is not unless it has an ocean — the ordering is set by the star through the period, which is the one quantity the flux-averaged habitable zone discards.

A year too short to feel its own eccentricity

A planet on an eccentric orbit can have a comfortable average and murderous extremes, and the habitable zone is drawn from the average. Whether the surface lives on the average or on the extremes is not decided by the flux at all — it is the ratio of how long the surface takes to change temperature to how long the year lasts, and the star sets the year.

6 figures · Habitable zone A planet in the middle of a 0.08 M☉ star's zone spends 94 Myr too hot to keep an ocean. How long a planet spends receiving more flux than the runaway-greenhouse limit while its star contracts onto the main sequence, against stellar mass, for planets at the inner edge, in the middle, at the outer edge of the zone the star will have once it settles. The luminosity is the contraction law of a fully convective star, falling as t^(−2/3) from an age of 1 Myr until arrival; the limit is the runaway flux for the star's temperature. A planet at the inner edge is too hot for 313 Myr round a 0.08 M☉ star and 17.8 Myr round a 0.6 M☉ one; a planet in the middle is too hot for 94 Myr round a 0.08 M☉ star and 5.5 Myr round a 0.6 M☉ one; a planet at the outer edge is too hot for 39 Myr round a 0.08 M☉ star and 2.0 Myr round a 0.6 M☉ one. A runaway greenhouse is not a hot climate; it is a state in which the ocean is entirely in the atmosphere as steam, where ultraviolet light splits it and hydrogen escapes. A planet in the eventual zone of the commonest stars in the galaxy begins its life in that state for tens to hundreds of millions of years — a span comparable with the whole assembly of the Earth. The durations are measured from 1 Myr; a rocky planet may take tens of millions of years to finish forming, and one that formed later misses the start of its exposure, while the smallest stars' arrival times are somewhat short in this model, which lengthens the end of it.

Steam before the zone existed

The smallest stars take hundreds of millions of years to contract onto the main sequence, shining at many times the luminosity they will settle at. A planet in the habitable zone such a star will eventually have spends that time with its ocean in the air as steam, while starlight splits the water and the hydrogen leaves — so the zone of the commonest star in the galaxy is a place that had to survive being too hot first.

5 figures · Habitable zone

Galaxies

Where the unknown stops being a number and becomes a profile — and most of it is not light.

A slope of 0.57 at the low-mass end, against the 0.59 and 0.24 two winds give there. Gas-phase oxygen abundance against stellar mass for star-forming galaxies. The solid measured curve is the relation found from electron-temperature abundances in stacked spectra, which rises as a power of mass below a turnover near 10^8.9 solar masses and saturates above it at 12 + log(O/H) = 8.798. The two model curves are a galaxy in equilibrium with its gas supply, whose metallicity is the yield divided by one plus the mass it ejects per unit mass of stars and one plus the dilution by the gas it must keep accreting. Only the ejection changes with mass. A wind driven by the energy of supernovae must lift gas out of a potential whose depth goes as v², so its loading goes as v⁻² and, with v ∝ M^(1/3), as M^(−2/3); a wind driven by momentum goes as v⁻¹ and M^(−1/3). Those set the limiting slopes far below the turnover, 0.67 and 0.33; at log M = 7.5, where dilution still matters, the drawn curves have slopes of 0.59 and 0.24. The measured slope there is 0.57. The dashed measured curve uses a strong-line calibration of the same galaxies' spectra; it sits higher, flattens sooner and has a slope of only 0.40 at log M = 9 — so which wind the data prefer is decided as much by the choice of abundance calibration as by the galaxies. Above the turnover every curve saturates, because a galaxy that ejects almost nothing keeps what it makes and its abundance is the yield, diluted.

The metals a galaxy keeps measure what it threw away

Small star-forming galaxies are metal-poor and large ones are not, in a relation tight enough to be a law. Read as an equilibrium between inflow, star formation and wind, its slope says how the wind is driven, its turnover says where galaxies stop losing what they make — and its redshift evolution says galaxies were poorer because they were still being filled.

5 figures · Chemical evolution A disc built from the inside out, with a gradient of −0.066 dex per kiloparsec at 12 Gyr. The metallicity of the gas, in solar units on a logarithmic scale, against galactocentric radius, at ages of 2, 6, 12 Gyr, for a disc in which every ring is its own box with infall: pristine gas arrives on a timescale that grows with radius, from 1 Gyr in the centre to 7 Gyr at 8 kpc, turns into stars on the depletion time of a Kennicutt law, which is shorter where the gas is denser and much longer below a threshold of 7 M☉ pc⁻², and keeps everything it makes. The slopes fitted between 4 and 14 kpc: −0.259 dex/kpc at 2 Gyr, −0.124 dex/kpc at 6 Gyr, −0.066 dex/kpc at 12 Gyr. The inner disc has had its gas early and turned it over many times, so it is near the yield; the outer disc is still accreting and forming stars slowly, so its gas is diluted and young in the chemical sense. At 8 kpc the model's present abundance is 1.16 of the yield. Every ring is independent here: no gas flows between them and no star moves, which are the two processes that real discs add and which both act to flatten what is drawn.

A gradient the old stars have walked away from

The gas in a disc galaxy is richer in metals near the centre than at the edge, by about six-hundredths of a dex per kiloparsec in the Milky Way. Two ingredients of disc growth make that slope, a disc that grows from the inside out makes it flatten with time — and the old stars that should carry the steeper history have moved several kiloparsecs from where they were born.

6 figures · Chemical evolution A power law has no timescale: half the Type Ia supernovae by 740 Myr and a tail to 13.7 Gyr. The fraction of all the Type Ia supernovae a single burst of star formation will ever produce that have exploded by a given delay, on a logarithmic time axis, for t⁻¹ from 40 Myr; t⁻¹·⁴ from 40 Myr; single delay of 1 Gyr; Gaussian, 3 ± 1 Gyr. The power law is what rates measured against host-galaxy ages and against the cosmic star-formation history both favour, and its cumulative fraction rises as the logarithm of the delay — equal numbers per decade of time. Half have exploded by the geometric mean of its limits, 740 Myr, 55 per cent by 1 Gyr, and the last are still exploding after a Hubble time. A single delay turns the whole population on at once; a Gaussian concentrates it at a characteristic age. The power law's shape has a physical reading: if white dwarfs explode when a pair of them merges by emitting gravitational waves, the merger time goes as the fourth power of their separation, and a broad distribution of separations becomes a distribution of delays with no preferred scale. What the drawing cannot say is which progenitors are involved — the measured rates constrain the shape and the normalisation, about one Ia per thousand solar masses of stars formed, and not the mechanism.

The iron clock has no single delay

The α-element knee is drawn as though Type Ia supernovae switched on a billion years after the stars that made them. Measured rates say otherwise — the delays are spread evenly over every decade from forty million years to a Hubble time, as a power law with no timescale in it — and a clock with no timescale bends where a clock with one would break.

6 figures · Chemical evolution

The observed sky

Coordinates, seasons, phases and shadows — geometry seen from inside it.

A right angle short by 0.147°, and a ratio of 389 hanging on it. The ratio of the Sun's distance to the Moon's implied by the angle between them at the moment the Moon is exactly half lit, on a logarithmic scale, for angles from 80° to 89.95°. At that moment the angle at the Moon between the directions to the Sun and the Earth is a right angle, so the ratio is the secant of the observed angle — a construction with no distance in it, the same right triangle that gives an inferior planet's orbit from its greatest elongation. Aristarchus measured 87°, which gives 19.1. The mean distances give 389, which corresponds to 89.853°. The curve's steepness is the whole story: across a tenth of a degree centred on each marked angle the ratio changes by 3.4 per cent at 87°, 10.5 per cent at 89° and 103 per cent at the true angle. The method was exact and it asked for a right angle measured to a hundredth of a degree, at an instant the eye can judge only to within hours, on a terminator that is never quite straight.

A right angle short by a seventh of a degree

When the Moon is exactly half lit, the angle at the Moon between the Sun and the Earth is a right angle, so the angle seen from the Earth gives the Sun's distance in units of the Moon's. Aristarchus measured 87° and concluded the Sun was nineteen times further away. The construction was exact; the angle he needed was 89.85°, and at that angle a tenth of a degree is the whole answer.

5 figures · Apparent motion Jupiter and Saturn meet every 19.86 years, tracing a three-cornered figure that turns 8.5° each round. The heliocentric longitudes at which Jupiter and Saturn are in conjunction — the same longitude seen from the Sun — for 21 successive conjunctions from 1800 to 2200, computed from Keplerian elements and dotted in three colours for the first, middle and last thirds of the span. The mean interval is 19.857 years, the synodic period the two mean motions give. Each conjunction falls 242.8° further round the orbit of Saturn than the one before, so 3 of them come back within 8.5° of where they started: the conjunctions sit near the corners of a 3-sided figure, and the figure itself rotates by 8.5° every 59.6 years. At that rate it returns to its starting orientation — a figure with 3 identical corners only needs to turn by a third of a turn — after about 838 years. The drawn corners are not exactly repeated because the orbits are ellipses: the planets move faster near perihelion, and the conjunction longitudes cluster where both are slow. That near-return is not a coincidence of dates. It is the statement that 3 synodic periods are close to a whole number of each planet's years, which is a near-commensurability of the two mean motions — and near-commensurabilities are where planets perturb one another most. The elements are a fit valid between 1800 and 2050; outside those years they are carried as fixed ellipses turning at their mean rates, which is right for the pattern and not for any individual date.

A triangle of meetings that turns in eight centuries

Jupiter and Saturn meet every twenty years, and each meeting falls about two-thirds of the way round the sky from the last, so the meetings trace a triangle. The triangle turns a third of a turn in 838 years because five of Jupiter's years almost equal two of Saturn's — and that same near-fit is the largest perturbation in the solar system, the one that made Saturn appear to be slowing down.

6 figures · Apparent motion Venus is brightest 51 days from inferior conjunction, 38 per cent lit. The brightness of Venus through one synodic period of 584 days, in magnitudes below its brightest, against days from superior conjunction — inferior conjunction at the two ends of the axis. The planet is treated as a matte, Lambert-scattering sphere on a circular orbit of 0.723 AU, so its flux is its phase function divided by the square of its distance from the Earth. The two factors fight: near inferior conjunction the planet is closest but shows only a thin crescent, and near superior conjunction it is fully lit but 1.72 AU away. For this orbit the contest has an interior winner. The brightest moment is 51 days either side of inferior conjunction, at an elongation of 44.6° from the Sun, with 38 per cent of the disc lit and the planet 0.539 AU away; greatest elongation, at 46.3°, comes 71 days from inferior conjunction, after the brightness peak on the way out from inferior conjunction. At superior conjunction the planet is 0.82 magnitudes fainter than its best. The marked point is the observed greatest brilliancy, about 36 days from inferior conjunction at an elongation near 39° — closer to conjunction and to a thinner crescent than any matte sphere on this orbit can be brightest at. A surface that sends more light forward, towards large phase angles, would move the peak exactly that way. Very near either conjunction the planet is lost in the Sun's glare, and the curve there describes light nobody sees.

Brightest as a crescent, and not as a disc

Venus is fully lit when it is furthest away and nearest when it is barely lit, and it is brightest in between, as a crescent weeks from inferior conjunction. A matte planet only has such a peak if its orbit is wider than about 0.46 of the Earth's — Mercury's is not, and Mercury is brightest full — and Venus's real peak sits closer to conjunction than any matte sphere allows, which is its clouds throwing light forward.

6 figures · Apparent motion

Spaceflight

Celestial mechanics used forwards — where to burn, and what it costs.

At the Sun–Earth L₂ the cheapest correction is every 23 days, and it costs e σ per e-folding. The annual station-keeping cost at the Sun–Earth L₂ point, against the interval between corrections, on logarithmic axes, for velocity errors of 0.5 cm/s, 2.0 cm/s, 5.0 cm/s along the unstable direction at each correction. Correcting often costs a lot because every correction carries its own error σ; correcting rarely costs a lot because the error has grown by e^(T/τ) in between, with an e-folding time τ = 23.4 days set by the point's growth rate of 2.484 times the orbital mean motion. The product (365.25/T) σ e^(T/τ) has its minimum at exactly T = τ, where the annual cost is 365.25 e σ/τ: 0.21 m/s a year for σ = 0.5 cm/s, 0.85 m/s a year for σ = 2.0 cm/s, 2.12 m/s a year for σ = 5.0 cm/s. The minimum is broad, so an operator can correct at a convenient interval near the e-folding time for little penalty, and the cost scales linearly with how well the spacecraft's velocity is known and executed. This is a one-dimensional caricature: a real halo orbit's correction also removes a stable component it need not, and solar radiation pressure on a large sunshield is a steady error source of its own. The figure's claim is the structure — an unstable equilibrium is cheap to hold if the instability is caught while it is still small, and its cost is a navigation budget rather than a force budget.

An unstable point that costs less to hold than a stable orbit

A spacecraft at the Sun–Earth L₂ point sits on an equilibrium that throws it away, doubling any error every sixteen days. It holds station for a few metres per second a year — a twentieth of what a geostationary satellite pays to stay on an orbit that is stable. The difference is what is being paid for — an instability caught small costs a navigation budget, and a steady torque costs a force budget.

6 figures · Station-keeping

Starlight

Distance, brightness and colour: what can be measured when nothing can be visited.

A bright star wants a wide aperture and a faint one wants 0.68 of the seeing. Signal-to-noise of simple aperture photometry against the aperture radius, in units of the seeing's full width at half maximum (1″), each divided by what optimal pixel weighting achieves for the same star, for stars of V = 12, 17, 20, 23 observed for 60 s through a 1 m telescope under a sky of 21 mag/arcsec². A small aperture loses starlight; a large one admits sky, and the balance depends on which dominates. For a bright star its own photons are most of the noise, so a wider aperture keeps gaining light almost for free and the best radius is large — 1.63 FWHM at V = 12, reaching 100.0 per cent of the optimum. For a star fainter than its sky the best radius shrinks to 0.680 FWHM and the best aperture reaches only 90.5 per cent of what weighting each pixel by its share of starlight divided by its variance achieves. That residual is exact in the background-limited limit: the best aperture captures 71.5 per cent of the light and 0.902 of the optimal signal-to-noise, so optimal weighting is worth 11 per cent in signal-to-noise, or 23 per cent in exposure time, and no more. The image is taken to be Gaussian; a real point-spread function has broader wings, which makes a fixed aperture a little worse and the optimal weights harder to know.

The best aperture throws away a tenth

Aperture photometry counts every pixel inside a circle equally and every pixel outside it not at all. For a faint star against its sky the best circle is two-thirds of the seeing wide, catches 71.5 per cent of the light, and reaches 90.2 per cent of the signal-to-noise that weighting each pixel by what it is worth achieves — a loss of 23 per cent in exposure time that no algorithm can beat by more.

6 figures · Photon noise Ten comparison stars as bright as the target cost 5 per cent in precision; ten 2 magnitudes fainter cost 23. The precision of a V = 12 target measured relative to an ensemble of comparison stars on the same 60-second frames, against the number of comparison stars, on a logarithmic precision axis, for comparisons 1 mag brighter, as bright as the target, 1 mag fainter, 2 mag fainter. A change in the atmosphere's transparency of 2.0 per cent — which would put the target's raw brightness out by 20.0 mmag — multiplies every star by the same factor and vanishes from the ratio. What is left is the target's own noise, 0.89 mmag, plus the ensemble's, which falls as the inverse square root of the number of stars in it. With comparisons as bright as the target the result is σ√(1 + 1/N): one comparison costs 41 per cent, ten cost 5. Fainter comparisons are noisier and need many more to reach the same point; brighter ones help, but the target's own noise is a floor the ensemble can only approach. Scintillation is treated as independent from star to star, which is right for stars more than a few arcseconds apart on a large telescope and makes it part of the noise that does not cancel. The figure also cannot show the defining weakness: every comparison star is assumed constant, and a variable among them injects its variability into every measurement made against the ensemble.

The comparison stars are part of the measurement

Measuring a star against others on the same frame cancels everything the atmosphere and the instrument do to all of them at once — a two per cent change in transparency vanishes completely. What does not vanish is the comparison stars' own noise, which the target inherits, and the variability of any comparison that is not constant, which the target reports as its own.

6 figures · Photon noise

Orbits

Kepler's three laws, and the family of curves a single force allows.

A wall at 2.33 hours that bends into a slope below 598 metres. The fastest rotation period a body of bulk density 2 g/cm³ can hold against its own spin, against its diameter, on logarithmic axes, for cohesive strengths of 0, 10, 100, 1000 pascals. With no cohesion the limit is the density-only barrier √(3π/Gρ) = 2.33 hours at every size. A cohesion C adds a stress that does not depend on size to a gravitational stress that goes as the square of the radius, so small bodies are held mostly by cohesion and can spin faster in proportion to their smallness: below the corner the limiting period goes as the diameter. The corner is where the two stresses are equal, at a diameter of 189 m for 10 Pa, 598 m for 100 Pa, 1.9 km for 1000 Pa. The scaling is the strength-regime form with a single coefficient of order one; detailed limits depend on the angle of friction and the shape. What the figure makes plain is why the observed spin barrier is sharp for kilometre-sized asteroids and fades below a few hundred metres, and why a handful of fast rotators a few hundred metres across can be rubble piles with a few tens of pascals of cohesion — a strength far below that of any rock — rather than monoliths.

A spin barrier with a corner in it

A rubble pile cannot spin faster than a period set by its density alone — 2.3 hours for most asteroids — at any size. Give it a cohesion of a few tens of pascals, weaker than any rock, and the barrier bends into a slope below a corner at a few hundred metres. The thermal torque that drives bodies to the barrier doubles their spin in a time that grows as the square of their size, so the bodies it pushes hardest are the ones that can go past.

6 figures · Rubble piles A rubble pile that splits below a mass ratio of 0.204 can lose its piece; above it, the piece stays. The total energy of two spherical components of equal density in contact, spinning together at the rate at which their mutual gravity just holds them against the spin, against the mass ratio of the smaller to the larger, in units of G m₁²/R₁. The energy is the kinetic energy of the rotating pair minus their mutual gravitational binding. When the smaller piece is a small fraction of the whole, the spin carries more energy than the binding and the total is positive: a body spun to breakup that sheds a fragment of that size has enough energy for the fragment to escape entirely, becoming a separate asteroid on a nearly identical orbit. The total changes sign at q = 0.204. Above that ratio the pair cannot separate without an energy source; it stays as a binary, orbiting and eventually synchronising, or re-accretes. As q goes to zero the energy tends to 0.2 G m₁²/R₁, the rotational energy of the primary alone at its breakup rate. Nothing in the threshold depends on the size or the density of the body — it is a pure number from the geometry of two touching spheres — and asteroid pairs sharing an orbit have been found overwhelmingly with estimated mass ratios below it.

A split that decides whether the piece can leave

A rubble pile spun past its limit splits in two, and whether the smaller piece escapes or stays in orbit is not a matter of luck. Two touching spheres spinning at their shared limit have positive total energy only when the smaller is less than 0.204 of the larger's mass — a number with no size and no density in it. Below it the pieces can become a pair of asteroids on nearly identical orbits; above it, a binary. And the larger the piece that leaves, the slower the body left behind.

6 figures · Rubble piles

Spaceflight

Celestial mechanics used forwards — where to burn, and what it costs.

A gauge good to 7 per cent with a tenth of the load left, and to 23 per cent with three hundredths. The uncertainty in the propellant remaining in a spacecraft tank, as a percentage of what remains, against the fraction of the 450-kilogram load still in the tank, on a logarithmic uncertainty axis with the tank emptying to the right. Bookkeeping — summing every thruster firing through a flow-rate model — carries an error common to all burns of 2 per cent of the mass used, plus an independent 5 per cent per burn that averages down over 2000 firings; its absolute error grows with the mass used. Gauging by pressure and temperature infers the empty volume of the tank from the gas law applied to a known mass of pressurant, with a combined 0.66 per cent uncertainty in n R T / P and a 0.2 per cent uncertainty in the tank's volume; its absolute error grows as the gas fills the tank. The two methods are independent and are combined by inverse variance. With a tenth of the load left the combined estimate is uncertain by 2.9 kg, 7 per cent of what remains; with three per cent left, by 3.1 kg, 23 per cent. Near empty the absolute error barely changes, so halving what is left doubles the relative error — the gauge is at its worst exactly when the last manoeuvre has to be planned from it.

A fuel gauge that is worst when it is needed

A spacecraft's tank has no float and no dial. The propellant left is estimated by adding up every burn or by reading the pressure and temperature of the gas above the liquid, and both methods' errors grow with the propellant used. Relative to what remains, the error doubles every time what remains halves — so a geostationary satellite has to hold back months of station-keeping as a margin against a gauge that cannot see the last few kilograms.

6 figures · Station-keeping

Gravitation

Two bodies pulling on each other, and everything that goes wrong at three.

A hole's shadow at 17° to its spin axis. The critical curve — the outline of a black hole's shadow on a distant observer's sky — for spins from 0 to 0.998, seen 17 degrees from the spin axis, in gravitational radii. Every curve is the projection of the unstable spherical photon orbits, and nothing about an accretion flow enters it. The sizes are 5.20, 5.12, 4.92, 4.84: a spread of 6.9 per cent across the entire range of spin, which is what makes an image of one of these a mass measurement with almost no model in it. What spin does instead is push the curve sideways and flatten one edge — the fast hole's outline sits 0.70 gravitational radii off centre, because photons co-rotating with the hole are dragged round and escape from closer in than counter-rotating ones can. That displacement is the spin signal, and it is measured against a centre nothing else marks.

A ring whose size is almost only a mass

The dark patch at the centre of an image of a black hole is the projection of the orbits light cannot escape from, and its size changes by seven per cent across the whole range of spin and inclination. That insensitivity is what makes an image a mass measurement with almost no model in it — and what leaves the spin in the part of the picture that barely moves.

7 figures · Black hole spin Two assembly histories, and the sign that separates them. The distribution of effective spin produced by two ways of making a black-hole binary, from 40,000 draws each. A pair that evolved together as two stars keeps its spins near the orbital axis — drawn here aligned to within 25 degrees — and cannot produce a negative projection at all: 0.0 per cent of its mergers fall below zero, and its mean is 0.45. A pair assembled by chance encounters in a dense cluster has no memory of any plane, so its tilts are isotropic, its distribution is symmetric, and exactly half of its mergers have an effective spin below zero. The sign of one number distinguishes two histories — and it does so without any of the individual spins being measured, which is the only reason the question is answerable at all from a catalogue in which every event's spins are individually uncertain. The third curve is a mixture with 50 per cent of its mergers drawn isotropically: its mean is 0.22 and 25 per cent of it lies below zero, which is the shape an observed distribution with a small positive mean and a real negative tail requires. Neither channel alone produces it.

A spin that is measured as an angular momentum

Every other method of measuring a black hole's spin reads it from what the hole does to matter. A merging pair's waveform reads it from the orbit instead — and what it recovers is one mass-weighted projection onto the orbital axis, which is a number that two non-rotating holes and two extremal ones can share.

7 figures · Black hole spin The quantity a capture has to change. The Jacobi constant along three trajectories through the same encounter with the secondary, integrated in the rotating frame at mass fraction 0.0009543. The undissipated one holds its value to 3e-11 over the whole passage — that is not an approximation but an identity of the equations of motion, and the residual is the integrator's. A body that enters through the neck at L₁ with the neck open leaves through it, because C has not changed and the neck is still open; a purely gravitational two-body encounter is time-reversible and cannot end in a bound orbit. The other two runs remove energy in the rotating frame. A drag proportional to velocity raises C continuously at a rate 2k v², and a single impulse against the motion raises it in a step. Both finish above C(L₁) = 3.039, which is the line at which the neck shuts and the body is trapped in the secondary's lobe. Which mechanism operated for a given satellite is not recoverable from its orbit today — but it is partly recoverable from the population, because the three leave different distributions behind.

A moon that two bodies could not have caught

A body that falls into a planet's Hill sphere under gravity alone leaves it again, because the quantity the rotating frame conserves is the same on the way out as on the way in. Every captured moon is therefore a record of something that removed energy — and which of the three candidates did it is written in the inclinations rather than in any one orbit.

7 figures · Hill sphere A limit with a notch cut in the middle of it. The stability limit for a satellite, in units of its planet's Hill radius, against the inclination of its orbit to the planet's. The two flat stretches are the coplanar answers — 0.5 prograde and 0.7 retrograde, the difference between them being the Coriolis asymmetry that a rotating frame imposes and no potential can hold. Between them the limit is not an interpolation. Inside the shaded band, from 39.23° to 140.77°, the Sun's averaged perturbation exchanges the orbit's inclination for its eccentricity, and a satellite whose semi-major axis is comfortably inside the lobe has an apocentre that is not — so the boundary is set by a(1 + eₘₐₓ) rather than by a, and it falls to exactly half the coplanar value at ninety degrees. The band edges are cos²i = 3/5 exactly, with no fitted quantity anywhere in them, and the notch is the reason a quoted stability limit in Hill radii is a statement about a coplanar orbit and not about a satellite.

A stability limit with a notch cut in it

The two numbers usually quoted for how far out a moon may orbit are coplanar numbers. A tilted orbit is not somewhere between them: inside a band of inclination fifty-six degrees wide the Sun converts the tilt into an eccentricity, the apocentre rather than the semi-major axis has to stay inside the lobe, and the limit falls to half. The band is empty in the sky.

6 figures · Hill sphere

Spaceflight

Celestial mechanics used forwards — where to burn, and what it costs.

The ratio that is √2 only for a point. Escape speed divided by circular speed, against distance in units of the Sun's, for a point mass and for flat-rotation-curve haloes truncated at four different radii. The point-mass value is √2 at every distance and it is a floor rather than a rule. Mass outside an orbit contributes nothing to the circular speed — a spherical shell exerts no force inside itself — and everything to the depth of the well, so any extended distribution has a ratio above √2, by an amount that measures how much lies beyond. For a flat rotation curve the escape speed is not even defined without an outer edge: the potential grows as the logarithm of radius and diverges, so the ratio is √(2 ln(R/r)) with R the truncation radius, and the answer depends on a boundary nobody can see. The measured value near the Sun is about 2.27, which on this curve puts the edge at 13 times the Sun's distance — and the strength of that inference is exactly the weakness of the method, because a ratio of 2.50 would put it at 23.

The ratio that is not √2

Escape speed is exactly √2 times circular speed at every distance from a point mass, and the Sun's own escape speed from the Galaxy is 2.3 times its orbital speed. The excess is not an error in either number. It is a measurement of the mass that lies outside the Sun's orbit, which pulls on nothing and holds everything.

6 figures · Escape What is left at infinity, and what it costs from 200 km. The velocity change needed to leave Earth from a circular orbit 200 kilometres up, against the speed the vehicle still has when the pull has run out. Escaping with nothing left over costs 3.22 km/s. Arriving at infinity with three kilometres a second in hand costs 3.63 — an extra 0.40 km/s to buy three, a discount of 87 per cent, because the two speeds do not add and their squares do — the departure speed is √(2μ/r + v∞²) — and adding a small square to a large one barely moves the root. The curve is therefore almost flat where the planetary missions are and straightens to slope 0.79 at the right-hand edge, where v∞ has overtaken escape speed and each metre a second bought costs a metre a second spent. The practical reading is that the expensive part of an interplanetary departure is getting to escape, and the destination is a surcharge on it.

The speed that is left over

No mission is designed to sit on the escape threshold. What an interplanetary departure is specified by is the speed still remaining when the planet's pull has run out — and because energies add while speeds do not, the first three kilometres a second of it cost four hundred metres a second to buy.

7 figures · Escape

Every sequence · Every named object · Search