Astronomy & Celestial Mechanics
the sky, computed
Field

Stars

The diagram that sorted them, and the one quantity that decides a star's whole life.
The Hertzsprung–Russell diagram. Luminosity against surface temperature, both in solar units and on logarithmic axes, with temperature increasing to the left. The main sequence is computed from the mass–luminosity and mass–radius relations; the dashed diagonals are lines of constant radius.

The diagram that sorted the stars, by plotting two things against each other

Plot brightness against colour for a few thousand stars and they do not scatter. They fall on a narrow band with two islands off it, and explaining that structure is most of stellar astronomy.

 Luminosity against mass, against a slope of 3.5. Main-sequence luminosity against mass, both in solar units, on logarithmic axes, over the range 0.079 to 63 solar masses. The measured curve comes from the eclipsing binaries and is the same in every drawing of it; what changes here is what it is compared against. The dashed line is a pure power law of exponent 3.5, and the curve crosses it rather than following it — the local slope runs from about 2.3 at the bottom of the range, where the interiors are convective, through nearly 4 near a solar mass where bound-free opacity dominates, to about 3 among the massive stars where electron scattering does. Quoting one exponent across the whole sequence is a convenience and the places it fails are the places the interior physics changes. Because the slope is between three and four across most of the range, a small spread in mass becomes an enormous spread in output: the 63-solar-mass end is 3.0e+9 times brighter than the 0.079-solar-mass end.

Mass decides everything, by a power of three and a half

Two stars of the same mass are almost the same star. Double the mass and the output multiplies by eleven — which is why a modest range of masses produces a colossal range of stars.

 How long a star lasts, against its mass. Main-sequence lifetime against mass, on logarithmic axes. Fuel grows in proportion to mass and consumption grows as its three-and-a-half power, so the lifetime falls steeply — a star of thirty solar masses lives for a few million years.

Why the biggest stars die first, and take the galaxy with them

A star thirty times the Sun's mass has thirty times the fuel and burns it forty thousand times faster. It lasts a few million years, and everything heavier than iron exists because of it.

 Energy generation against core temperature. The proton–proton chain and the CNO cycle, in solar units, against core temperature on logarithmic axes. The CNO curve is far steeper, so the two cross at 18.8 million kelvin — above that temperature a star runs mostly on CNO, and below it mostly on pp.

The furnace that runs cooler than a compost heap

The Sun's core produces 276 watts per cubic metre. A human body produces more. The Sun is bright because it is enormous, and its fusion is astonishingly slow — which is exactly why it lasts.

 The Hertzsprung–Russell diagram. Luminosity against surface temperature, both in solar units and on logarithmic axes, with temperature increasing to the left. The main sequence is computed from the mass–luminosity and mass–radius relations; the dashed diagonals are lines of constant radius.

The main sequence is a place stars sit, not a track they travel

The commonest misreading of the Hertzsprung–Russell diagram is that stars slide down the band as they age. They do not. They sit at one point for ninety percent of their lives and then leave sideways.

 Radius against mass for a degenerate star. The mass–radius relation for electron-degenerate matter. More mass gives a smaller star, and the radius reaches zero at 1.46 solar masses — the Chandrasekhar limit, solved from the same expression that draws the curve rather than quoted alongside it.

A star held up by a rule about counting

A white dwarf makes no energy and does not collapse. What holds it up is not heat or pressure in any ordinary sense — it is a quantum rule forbidding two electrons from occupying the same state.

 Inside a star that is holding itself up (polytrope n = 3). Temperature, density and pressure through a star, as fractions of their central values, against fractional radius. All three come from one numerical integration of the Lane–Emden equation at index 3, which is the hydrostatic balance written for a gas whose pressure is a power of its density. The inner half of the radius holds 90% of the mass, and the outer half is nearly weightless — which is why the load, and so the temperature, is concentrated where the burning is.

A star is held up by its own weight

A star has a central temperature because it has a central pressure, and it has a central pressure because everything above is pressing down. The nuclear reactions do not set that temperature — they obey it.

 Leavitt's law, through the Milky Way calibrators. Mean absolute magnitude against the logarithm of the pulsation period, for classical Cepheids with independently known distances, and RR Lyrae itself for comparison. The line is M = −2.81 log P − 1.43: a tenfold longer period is 2.81 magnitudes brighter, a factor of 13 in luminosity. The vertical axis runs the astronomers' way, with brighter upward.

A star that tells its distance by how slowly it blinks

Some stars pulsate, and the slow ones are the bright ones. That single correlation turns a clock into a ruler, and it is how the size of the universe was first measured.

 Two radial-velocity curves, and one mass ratio. The line-of-sight velocity of each star through one orbit of AI Phoenicis. Both curves are computed from the two masses and the period; what a spectrograph delivers is the reverse. The ratio of the amplitudes is the inverse ratio of the masses — 48.2 to 50.3 kilometres a second, so the heavier star moves more slowly — and the sum of the amplitudes with the period gives the mass sum, 2.437 solar masses, once the inclination is known from the eclipses.

The only stars whose masses are known

A star's mass cannot be measured by looking at it. It can be measured by watching two stars pull on each other, and if the pair also eclipses, the same observations give both radii as well — with no stellar model anywhere in the chain. A few hundred such systems calibrate everything else.

 The main sequence against its own ceiling. Luminosity against mass, in solar units, with the Eddington limit drawn on the same axes. The limit is exactly linear in mass. The main sequence is not — until it is: above about 55 solar masses the relation flattens to a straight proportionality as well, and the two lines then run parallel at a fixed ratio of 83 per cent. They never cross on this range, which is the more interesting outcome: the most massive stars do not approach their limit gradually, they are built at a fixed fraction of it. The Sun, by contrast, radiates 0.0026 per cent of its own ceiling.

A brightness that would blow the star apart

Light pushes. For a star of a given mass there is a luminosity at which the outward push on its own outer layers equals the inward pull of gravity, and it depends on nothing but the mass and the opacity — not on the star's structure, its composition or its age.

 A furnace with the thermostat taken out. Left, the two pressures in a red giant's helium core at 10⁶ g/cm³. The ideal-gas pressure rises with temperature, as the whole regulation of a star depends on it doing; the degenerate electron pressure is a horizontal line, because the exclusion principle does not know the temperature. At the ignition point near 10⁸ K the degenerate term is 5.1 times the gas term, so a temperature rise there raises the total pressure by almost nothing and the core does not expand. Right, what that costs. Both curves start at the same temperature with the same triple-alpha rate, whose logarithmic slope is 41 — the famous T⁴⁰, differentiated here rather than quoted. The regulated core settles; the degenerate one, with nowhere to put the extra energy, goes vertical. The flash reaches 10¹⁰ solar luminosities for a few seconds and not one photon of it is seen: every erg goes into lifting the degeneracy, and the star's visible response is to become fainter and settle on the horizontal branch.

The resonance that had to exist

Three alpha particles cannot meet at once, and the two-body intermediate falls apart in 8×10⁻¹⁷ seconds. That is a hundred thousand times longer than a crossing time, which is just enough — provided a nuclear level sits at 7.65 MeV. Hoyle argued from the existence of carbon to the existence of the level, and it was found where he said.

 From 1.667 to 1.343, and the limit is in the second number. The local slope Γ = d ln P / d ln ρ of the exact degenerate electron equation of state, against density, differentiated numerically from the drawn pressure rather than quoted. It is 1.6666 at 10⁴ kg/m³ — the 5/3 of a non-relativistic Fermi gas — and 1.3430 at 10¹² — the 4/3 of a relativistic one — falling monotonically between, and passing 1.5 at 3.63·10⁹ kg/m³. Nothing gives way at either end; the electrons simply run out of room to go faster, because they are already near c. The mass limit is in the exponent and not in the strength. Hydrostatic balance for a polytrope of index Γ gives M ∝ ρ_c^(3Γ−4)/2, so at Γ = 4/3 the exponent is zero: the mass no longer depends on the central density at all, and there is exactly one mass a fully relativistic degenerate star can have — 1.459 M☉ at μ_e = 2. Above it, squeezing harder buys pressure that rises more slowly than gravity does, and there is nothing left to stop at.

The mass a cold star cannot exceed

A degenerate star's radius falls as its mass rises, and nothing in that relation suggests a limit. Making the electrons relativistic softens the pressure law until pressure and gravity scale the same way with radius — and the radius drops out of the balance, leaving one mass and no room to argue.

 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.

 Three clusters, three ages, one diagram. Isochrones for populations of 100 Myr, 1 Gyr, 12 Gyr, each drawn as the main sequence up to its own turnoff and then the post-main-sequence track of the turnoff mass. The turnoffs are at 5.52 M☉, 2.20 M☉, 0.94 M☉, from t = 10¹⁰ M/L with this file's own mass–luminosity relation. Nothing here is a track along which a star moves. Every point is a different star of a different mass, all the same age, and the bend is simply where the population runs out of stars that have had time to leave. That is why a cluster has an age and a field star does not: the bend needs a population, and one star is not one. The oldest globular clusters sit near the 12 Gyr line, and in the 1990s the same construction gave them ages of 16 to 18 Gyr against a universe measured at 10 — a two-standard-deviation contradiction that was resolved from the distance side, by Hipparcos, and not from this one.

An age read off a bend

A cluster is a population of one age and many masses, so the place where its stars leave the main sequence is a clock. The bend needs a population — which is why a cluster has an age and a single star does not.

 The mirror: an envelope 62 times larger around a core 1.48 times smaller. Envelope radius and core radius against helium-core mass for a 1 solar-mass star, both in solar radii on one logarithmic axis. Two curves with opposite slopes at every point: the envelope grows from 2.0 to 122 solar radii — 0.57 astronomical units — while the degenerate core contracts from 13,848 to 9,368 kilometres, and the two are separated by a factor of 364 in the middle of the range. The figure is drawn from two stated relations rather than from a stellar model: a core-mass–luminosity law of the sixth power, calibrated to 2,500 solar luminosities at 0.48 solar masses of core, and a Hayashi temperature falling from 5000 to 3700 kelvin across the range. Everything else follows exactly — the radius by the Stefan–Boltzmann law, the core radius by the cold degenerate relation R ∝ M⁻¹ᐟ³. Note what the tip radius is and is not: it is the red-giant branch, and the later asymptotic-giant tip is a different and larger number.

The star that swells because its centre shrank

A helium core contracts, and the envelope around it expands by a factor of sixty. The central density rises at the same time as the radius, which is the opposite of what a self-gravitating body is expected to do, and a burning shell between the two is the whole reason it happens.

 14 radial orders of the Sun, at 135.1 μHz apart. The p-mode spectrum of the Sun — 1 solar mass in 1 solar radius — from the asymptotic relation with its second-order term, drawn as 14 radial orders of ℓ = 0, 1 and 2 under a Gaussian envelope centred on ν_max = 3,090 μHz. Two numbers are marked and they do very different work. The large separation, 135.1 μHz, is the spacing between consecutive ℓ = 0 modes and fixes the mean density. The small separation, 9.00 μHz, is 2.7 pixels on this axis — it fixes the age, and it is why the échelle diagram exists rather than being a convenience. The vertical axis is the measurement: each mode moves the surface by about 20.0 cm s⁻¹ at the peak, and brightens it by a few parts per million, which is why this was impossible before a decade-long velocity series. Each mode is one line: its true width is set by its lifetime and is far below a pixel here.

The interior read from a comb of frequencies

A star's surface moves by about twenty centimetres a second, in thousands of overlapping sound modes at once. Two numbers off that spectrum give a mass and a radius with almost no stellar model in the chain, and a third gives an age.

 One population, three integrals, two ends. The same mass function weighted three ways, each normalised so that the area under it is one, against mass on a logarithmic axis — so a share of the page is a share of the total. The light is a zero-age population's, every star still on the main sequence, which is the only age at which the top of the range is present at all. The steep curve on the left is the number of stars, which peaks at 0.08 M☉ and falls away because the function is steeper than m⁻¹; the middle one is the mass they carry; the one on the right is the light they emit, computed from this file's own main-sequence relation and peaking at 55.0 M☉ — 682 times further up the axis. The population that is counted and the population that is seen are two different populations. Above 8 M☉ there are 0.63% of the stars, 21% of the mass and 99% of the light. The slopes are α = 1.3 above 0.08 M☉, α = 2.3 above 0.5 M☉, and the low-mass end is where the honesty runs out: it has never been measured in another galaxy, and every mass inferred from a luminosity assumes it.

Most stars are small, and most of the light is not

The mass function is steep and the mass–luminosity relation is steeper, so counting stars and measuring their output are integrals of the same function dominated by opposite ends of it. Every mass inferred from a brightness passes through that mismatch.

 A light curve that is a decay chain. Bolometric luminosity against days since explosion, for three Type Ia models synthesising 0.3 M☉, 0.6 M☉, 0.9 M☉ of ⁵⁶Ni, from Arnett's one-zone diffusion model. Nothing in the shape is fitted: the two timescales are the laboratory half-lives of ⁵⁶Ni and ⁵⁶Co, 8.8 and 111.3 days, and the rise is those decays seen through an envelope that takes 11.3, 15.5, 18.6 days to leak. Each peaks at 14, 18, 21 days, and at maximum the luminosity equals the instantaneous deposited power to within 0.6% — Arnett's rule, and the only reason a peak brightness can be read as a mass of nickel. After maximum the curves settle to nearly straight lines on this logarithmic axis, declining at 0.0360 magnitudes a day against ⁵⁶Co's own 0.0098. The difference is gamma-ray escape: the ejecta become transparent to the very photons that are supposed to be heating them, and the tail is therefore steeper than the isotope. The one thing the tail is not is a property of the star — it is a half-life, plus a column density that is falling as t⁻².

The candle that has to be standardised

A Type Ia's light is the decay of half a solar mass of nickel-56 seen through an expanding envelope. Their peak brightnesses span six-tenths of a magnitude — and it is the width of the light curve, not anything else, that says which.

 A furnace with the thermostat taken out. Left, the two pressures in a red giant's helium core at 10⁶ g/cm³. The ideal-gas pressure rises with temperature, as the whole regulation of a star depends on it doing; the degenerate electron pressure is a horizontal line, because the exclusion principle does not know the temperature. At the ignition point near 10⁸ K the degenerate term is 5.1 times the gas term, so a temperature rise there raises the total pressure by almost nothing and the core does not expand. Right, what that costs. Both curves start at the same temperature with the same triple-alpha rate, whose logarithmic slope is 41 — the famous T⁴⁰, differentiated here rather than quoted. The regulated core settles; the degenerate one, with nowhere to put the extra energy, goes vertical. The flash reaches 10¹⁰ solar luminosities for a few seconds and not one photon of it is seen: every erg goes into lifting the degeneracy, and the star's visible response is to become fainter and settle on the horizontal branch.

The explosion that never reaches the surface

When helium ignites in a degenerate core the thermostat every other burning stage runs on is disconnected. The runaway reaches ten billion solar luminosities for a few seconds, and the star's visible response is to get fainter.

 The temperature of a disc around a stellar black hole. Effective temperature against radius, in units of the inner edge, for a stellar black hole of 10 solar masses accreting 10⁻⁸ solar masses a year. Two features are structural. The profile turns over rather than rising all the way in: the factor (1 − √(r_in/r)) is the statement that no torque acts across the inner edge, so nothing is dissipated there and the peak sits at 49/36 of it, measured here at 1.361. And outside a few inner radii the run is exactly r^−3/4, drawn as the dashed line, which is what makes a disc's spectrum broad: every decade of radius contributes at a temperature a factor of 5.6 lower. The peak is 3.46·10⁶ K here, so the disc radiates in X-rays, and integrating the whole profile gives 4.72·10³⁰ W — which is GMṀ/2r_in to a per cent, half the binding energy released and no more, because the other half is still going round.

The disc that has to throw angular momentum away

Matter cannot simply fall onto a compact object. At the energy it arrives with, it has far too much angular momentum, and the only way in is for some of it to be carried outwards — which is what a disc is for.

 Two stars at the same point of every other diagram, 4.1 times apart in one. Above: the gravity-mode period spacing against the large frequency separation, for 90 shell-burning giants and 55 core-burning ones. These are the same stars in every other measurement. They have the same luminosity, the same temperature, the same colour, the same surface gravity and — inside the band drawn — the same Δν, which means the same mean density; the scaling relations of the rung below return the same mass and the same radius for both. What separates them is a quantity that comes from nowhere near the surface: ΔΠ₁ is set by the buoyancy frequency integrated across the core, and a core that has ignited helium is expanded and convective, so its integral is smaller and its period spacing larger. The two sequences do not touch — 61 seconds against 248 — and the gap sorts a catalogue of tens of thousands of giants into stars burning hydrogen in a shell and stars burning helium in a core, by a Fourier transform of a light curve. Below: what that looks like in the spectrum itself, drawn over two radial orders. The number of mixed ℓ = 1 modes between consecutive radial modes is Δν/(ΔΠ₁ν²), so the shell burner has about 12 of them and the core burner about 3: the star with the larger period spacing has the sparser spectrum. The picture cannot show what the same modes are also used for and cannot settle — the splitting of each mixed mode gives the rotation rate of the core separately from the envelope, and the cores come out spinning some ten times faster than the surface and a hundred times slower than any model of angular-momentum transport predicts.

Two stars only a Fourier transform can tell apart

A giant burning hydrogen in a shell and one burning helium in its core sit at the same luminosity, the same temperature and the same mean density. Every scaling relation returns the same mass and radius for both. The gravity-mode period spacing is fifty seconds for one and three hundred for the other.

 Two measured numbers, and everything else on the page derived from them. 25 pulsars in the plane of period against period derivative, at their catalogued values. Only the two axes are measurements; the three families of contour are models. Constant surface field runs at slope −1 because B ∝ √(PṖ), constant characteristic age at slope +1 because τ = P/2Ṗ, and the two families cross the population at right angles — which is why a single dot fixes both. The Crab sits at 3.8·10¹² G and 1257 years, and its true age is 972; the millisecond pulsars at the lower left have fields ten thousand times weaker and characteristic ages of billions of years, because they were spun back up by a companion long after they died. The line at the lower right is the death line, B/P² below which the model says no pair production and therefore no radio emission — and J2144−3933 is drawn below it, an 8.5-second pulsar that is radiating anyway.

A clock that runs down and says what it is

A pulsar hands over two measured numbers, a period and its derivative. Everything else usually quoted about it — an age, a magnetic field, a luminosity — is derived from those two through a model, and the one measurement that tests the model refutes it.

 The period spans 492× and the pulsation constant 2.1×. The pulsation constant Q = P√(ρ̄/ρ̄_⊙) for 7 radial pulsators, against their periods. A radial pulsation is a standing sound wave across a star and the sound speed in a self-gravitating body is set by its own gravity, so the period should go as (Gρ̄)^−1/2 and Q should be the same for every star pulsating in the same mode. Across a factor of 492 in period — from 1.3 hours to 27 days — Q varies by 2.07, with a mean of 0.0393 days. That is the whole reason a period–luminosity relation can exist: a period is a density, a density is a mass and a radius, and a radius with a temperature is a luminosity. The masses here are the weak link and the caption should say so — for the Cepheids they come from evolutionary models rather than from a dynamical measurement, and the drift of Q with period is where that assumption is showing.

The valve that has to sit at the right depth

A period–luminosity relation only exists because a period is a density in disguise. Which stars have a period at all is decided by whether a partial ionisation zone sits deep enough to have mass behind it and shallow enough that convection has not taken it over.

 A distance of 52.0 parsecs with nothing underneath it. Two ways to a distance for the same pair. The orbital parallax needs no iteration and no assumption: a double-lined spectroscopic orbit gives the relative orbit's linear size as (K₁+K₂)P√(1−e²)/2π sin i = 0.2268 AU, an astrometric orbit gives its angular size as 4.36 milliarcseconds, and the ratio is 52.0 parsecs — a length divided by an angle, with no rung of the distance ladder below it and no property of the stars assumed. The curves show the dynamical parallax, the version available when only one spectrum can be measured: guess the mass sum, take the linear size from the harmonic law, divide by the angular size, convert the apparent magnitude to an absolute one and read a new mass sum off a mass–luminosity relation. Three starting guesses spanning a factor of 10 in mass converge to the same distance in 8 passes and agree to 0.001 per cent. It converges because the distance depends on the assumed mass only as its cube root — the measured exponent here is 0.3333 — so a factor of two in the mass is 26 per cent in the distance, and one pass removes most of that. What it converges to is not the orbital parallax: the iteration settles at 54.2 pc against 52.0, 4.2 per cent away, because the fixed point is set by the mass–luminosity relation and the apparent magnitude rather than by anything measured about this orbit. The same insensitivity that makes it converge is why it is never better than the relation it leans on.

Two orbits of one pair, and a distance falls out

Measure the same binary spectroscopically and astrometrically and the orbit comes back twice — once as a length in kilometres and once as an angle on the sky. The ratio is a distance that owes nothing to parallax, nothing to a standard candle, and nothing to any assumption about the stars.

 Three ways for a timing model to be wrong, and three shapes that say which. Timing residuals over 6 years for the Crab pulsar, one curve per kind of error in the model, in microseconds. A position error of 1.2 mas leaves a sinusoid of period exactly one year — measured off the drawn curve as 1.000 — with amplitude (a/c)·δθ·cos β = 2.7 μs, because the error is being projected onto a baseline that is the Earth's own orbit and nothing about the pulsar. An unmodelled proper motion of 0.9 mas/yr leaves the same sinusoid with an envelope growing linearly: twice as large at 6 years as at 3. An error of one part in 10⁹ in Ṗ leaves a parabola, the second integral of a frequency drift, whose second derivative is constant to 4e-12 across the span — and that fractional error is deliberately minute, because anything larger produces a residual thousands of times the other two and draws them as flat lines. The shapes do not resemble each other, which is the whole reason a pulsar is an instrument rather than a clock: fitting them simultaneously delivers a position, a proper motion and — from the annual curvature term, not drawn here — a parallax, all from the arrival times of pulses and no image of anything. What is left when every known shape has been removed is the science: glitches, red noise, and the correlated residual between pairs of pulsars that a timing array exists to find.

A clock read against a model of everything in between

A pulsar delivers arrival times and nothing else. Everything else is a model, and what is measured is the difference between the model and the arrivals — a residual whose shape says which term is wrong, and whose annual sinusoid is a position measured from pulses rather than from an image.

 4 cycles of wings, and the polarity reverses at every boundary. Sunspot latitude against date, one mark per spot group, over 4 cycles of 11 years. The pattern is the reason the plot is called a butterfly diagram, and both of its features are laws with names. Spörer's law is the downward slope: spots emerge near ±28° at the start of a cycle and near ±7° at the end, so each wing narrows towards the equator and never crosses it. Hale's law is what the two mark shapes say: the leading spot of a pair has one magnetic polarity in the north and the other in the south, and both reverse when a new cycle starts — so a diagram that repeats every 11 years in appearance repeats only every 22 in magnetism. The wings overlap: the first high-latitude spots of a cycle appear about 1.6 years before the last low-latitude spots of the one before, which is why counting spots gives a cycle length slightly different from measuring one between polarity reversals. The dot density in time is the sunspot number itself, drawn from the standard skewed fitting function — the rise to maximum takes about four years and the decline about seven, in every cycle ever recorded.

A magnetic clock read off a butterfly

Plot sunspot latitude against date and the marks form wings that open at thirty degrees and march to the equator. The polarities reverse between wings, so the magnetic period is twenty-two years and the famous eleven is an artefact of counting spots rather than fields.

 7 solar masses of progenitor become 0.56 of white dwarf. Above: the initial–final mass relation. The steep line is what a star would leave if it kept everything; the shallow one is what it actually leaves, M_f = 0.08M_i + 0.489, fitted to white dwarfs in open clusters. Everything between 1 and 8 solar masses ends up between 0.57 and 1.13 — a range of 7 compressed by a factor of 12 — and the fraction thrown away rises from 43% at the bottom to 86% at the top. A star spends most of its mass in the last per cent of its life, at rates no theory predicts from first principles, and this relation is measured by running stellar evolution backwards: a cluster's turn-off gives the age, a white dwarf's temperature gives its cooling time, the difference gives how long its progenitor lived, and a model turns that into the mass it was born with. Below: the white dwarf mass distribution that follows, from a Salpeter initial mass function pushed through the relation above. It peaks at 0.61 solar masses, which is what every survey of white dwarfs measures. Two things make that peak and the flatness of the relation is only one of them: the unsmoothed distribution (the stepped curve) falls monotonically and is largest at its low-mass edge, because the initial mass function is steep. The edge is where it is because no star lighter than about 1.0 solar masses has had time to die yet, and the scatter of the measurement rounds that cliff into a maximum a little above it. A flat relation produces a narrow range; it takes the age of the Galaxy to produce a peak.

The mass a star does not keep

The initial–final mass relation is nearly flat, so everything from one to eight solar masses ends as a white dwarf between 0.57 and 1.13 — and the observed distribution piles up at 0.6 almost regardless of what went in. The relation is measured by running stellar evolution backwards, because nothing predicts the loss.

 Seven sources over twelve decades of flux, and one threshold below the one that matters. The Sun's neutrino spectrum at the Earth, with the continua drawn per unit energy and the two monoenergetic lines as spikes at their own energies. The pp reaction supplies 91 per cent of all of them and its endpoint is at 0.4233 MeV. The vertical lines are experimental thresholds, and they are the figure's argument: only gallium sits below that endpoint. Chlorine, which produced the deficit and held it for twenty years, could not see a single pp neutrino — it counted ⁷Be and ⁸B, which are a rare branch of a rare branch, together under a per cent of the total — and the water detectors that followed were higher still. So the discrepancy that eventually turned out to be a property of the neutrino was measured, for two decades, using the least representative one per cent of the flux available. The total drawn here is 6.54·10¹⁰ cm⁻² s⁻¹, and it is checkable without any stellar model at all: every completed chain turns four protons into helium, releases 26.73 MeV of which 0.59 leaves as neutrinos, and emits two neutrinos — so the solar constant of 1361 W m⁻² fixes the number at 6.5·10¹⁰ cm⁻² s⁻¹, within 0.6 per cent of the sum of the model's own branches. The dominant flux is a consequence of the Sun shining and of nothing else.

The only thing that leaves the centre

A photon made in the Sun's core takes a hundred thousand years to get out and arrives thermalised past recognition. A neutrino takes 2.3 seconds and arrives unchanged, so its flux is the fusion rate now — which is why a factor of three could not be absorbed by any adjustment to the Sun.

 The light curve of AI Phoenicis, computed from its elements. Total light against orbital phase, computed by overlapping two discs of radius 1.805 and 2.9303 solar radii at an inclination of 88.5°, each weighted by its own surface brightness. The two eclipses hide the same area of sky and have different depths — 48.0 and 19.1 per cent — because what is lost is the light of whichever star is behind, and the ratio of the depths is therefore the ratio of the two surface brightnesses. Two things are left out and both matter to a real solution: this is the bolometric light rather than the light in a filter, and the discs are uniform, where a real one is limb-darkened and so has a deeper, rounder eclipse than the flat-bottomed one drawn here.

Two radii, from a light curve alone

The radius of a star is not measured. It is inferred, from a temperature and a luminosity, through a model. There is one exception — a pair of stars that eclipse each other, whose light curve and velocity curves between them give both radii, both masses and the ratio of temperatures with no model of a stellar interior anywhere in the chain.

 A spectrum belonging to no temperature at all. The disc's summed emission, with the individual annuli drawn faintly beneath it. Each ring is a blackbody at its own temperature, and each is drawn at the area it actually has — the outer rings are cool and enormous, the inner ones hot and small. The sum has three parts and only the two ends belong to a temperature: a Rayleigh–Jeans rise of slope 2 from the outermost ring, a Wien cutoff at the hottest, and between them a stretch of slope 0.316, against the 1/3 that comes out of integrating ν²T(r)r dr with T ∝ r^−3/4. That middle section is the observational signature of a disc: no single blackbody produces it, no photosphere produces it, and its width rather than its peak is what says how far in the disc goes. What the figure cannot show is that a real disc's innermost rings are neither thin nor blackbodies, which is where the model's clean edges stop.

A spectrum that is a stack of temperatures

An accretion disc is not hot. Its inner edge is, its outer edge is not, and the temperature runs continuously between them as a power of radius — so what leaves the disc is the sum of a great many blackbodies at different temperatures, which is a spectrum with no temperature in it and a slope no single body can produce.

 A measured mass of 2.08 deletes an equation of state. Mass against radius for three neutron-star equations of state, each a polytrope P = Kρ² integrated through the Tolman–Oppenheimer–Volkoff equation from the centre outwards until the pressure reaches zero. Every sequence rises, turns over and falls; only the rising part is stable, because past the maximum adding mass makes the star smaller and the smaller star cannot hold itself up. The maxima here are 1.60, 2.10, 2.57 solar masses, in the order of increasing stiffness — the same nuclear matter with a slightly harder response to compression supports a heavier star, and nothing else in the calculation changes. The horizontal band is PSR J0740+6620, whose mass of 2.08 ± 0.07 solar masses comes from the Shapiro delay of its own pulses passing its companion, which is a timing measurement and involves no model of the star at all. It sits above the maximum of one of the three, and those are not disfavoured but excluded: an equation of state that cannot hold up two solar masses is wrong, whatever else recommends it. The two lines at the left are exact and no star may cross them — the Schwarzschild radius, and the bound above it inside which the speed of sound in the matter would exceed the speed of light. A caution about the curves themselves: a Γ = 2 polytrope is a stand-in for nuclear matter and runs a kilometre or two large in radius at fixed mass, so read the ordering and the maxima rather than the radii.

A radius that decides what matter can be

Nobody can make matter at four times the density of an atomic nucleus, and no calculation settles what it does there. What can be done is to weigh a neutron star — every candidate description of that matter predicts a heaviest star it could hold up, and a single measured mass above that value deletes the description permanently.

 Four kinds of explosion, told apart by which line is missing. Spectra of 4 supernova types — Ia, II, Ib, Ic — near maximum light, stacked, over the same wavelength range and drawn with the same photospheric expansion speed of 10,000 kilometres per second. Every feature is a P Cygni profile: a trough blueshifted by the material expanding towards the observer and a peak at rest wavelength from the material moving across the line of sight, and the width of both is the expansion speed. Reading the silicon trough of the type Ia off the drawing puts it at 614.0 nanometres, which recovers 10,131 km/s from a rest wavelength of 635.5. The classification is a decision tree on absences. Hydrogen present makes it a type II; no hydrogen but silicon makes it a Ia; no hydrogen and no silicon but helium makes it a Ib; none of the three makes it a Ic. That tree was written down before anybody knew what these objects were, and it separates a detonating white dwarf from a collapsing massive core almost perfectly — because what the letters are actually reading is how much of the star's outer envelope was still there when it exploded. The type II kept its hydrogen, the Ib had lost it and shows the helium beneath, the Ic had lost that too.

Two explosions told apart by a missing line

The classification of supernovae is a decision tree on absences — is hydrogen there, is silicon there, is helium there — and it was drawn up decades before anybody knew what any of these events were. It nevertheless separates a detonating white dwarf from a collapsing stellar core almost perfectly, and the reason it does is worth the essay.

 Two responses to the same transfer, turning round at q = 0.79 and q = 1. What conservative mass transfer does to the orbit and to the lobe, plotted against the mass ratio of donor to accretor on a logarithmic axis. Both curves are logarithmic derivatives with respect to the donor's mass, so a positive value means the quantity shrinks as the donor loses mass and a negative one means it grows. The orbit's response is exactly twice the mass ratio less one, which follows from holding the total mass and the total angular momentum fixed and nothing else, and it crosses zero at equal masses: transfer from the heavier star draws the orbit in, transfer from the lighter one pushes it out. The lobe's response adds to that the change in the lobe's shape, and it crosses zero earlier, at a mass ratio of 0.788. Between those two crossings the orbit is still widening while the lobe is already closing. To the right of both, a donor that loses mass finds its lobe shrinking around it, which is the runaway the essay is about: the transfer narrows the valve it is flowing through.

The flow that narrows its own channel

Two stars close enough share a surface, and the point where that surface pinches is a valve. What comes through it changes the orbit, and the orbit changes the valve — with a sign that reverses at equal masses, which is why some binaries transfer quietly for a hundred million years and others tear themselves apart in a thousand.

 Accretion from a medium a body is moving through, against how fast it moves. The rate at which a gravitating body captures gas out of a medium it is ploughing through, divided by the rate it would capture at rest, against its Mach number, both axes logarithmic. The accretion radius is set by where the body's escape speed matches the speed of the gas relative to it, and that relative speed combines the sound speed with the motion in quadrature; the rate carries the square of that radius times the speed, which leaves one plus the square of the Mach number to the power minus three halves. Subsonic motion therefore costs almost nothing — the curve is flat below Mach one third — while supersonic motion costs the inverse cube of the speed, drawn here with a measured logarithmic slope of -2.99. The same focusing produces the wake and the drag: the body pulls a denser column behind it, that column pulls back, and the material closest to the axis is captured. Accretion and dynamical friction are not two processes but two accounts of one, which is why a body that grows by this mechanism is also being slowed by it.

The wake and the meal are one calculation

A body moving through gas gathers what passes inside the radius at which its escape speed matches the flow. That single radius sets both the drag it feels and the rate at which it grows, so accretion and friction are not two processes but two readings of one — and the rate falls as the inverse cube of the speed.

 A line nothing spun up by accretion can lie above, and the millisecond pulsars beneath it. The period–period-derivative diagram with the spin-up line drawn on it. An accreting neutron star is torqued by the disc until its magnetosphere turns at the same rate as the material arriving there, which fixes an equilibrium period as a function of the magnetic field and the accretion rate. Eliminating the field between that relation and the dipole formula that every point in this diagram is already read through leaves a straight line of slope 1.33, drawn here for accretion at the Eddington rate — the fastest a star can be pushed. The 7 recycled pulsars all sit below it, which is what the figure is for: none of them was spun up faster than the limit allows, and their positions are a record of how much mass each one received rather than of how old it is. The young pulsars are in the opposite corner, above the line and to the right, spinning down from birth. The two populations are not two stages of one life. A star that reaches the bottom left has been fed by a companion for a hundred million years, which is why almost every millisecond pulsar has one and almost no young pulsar does.

A corner of the diagram that has to be earned

A pulsar spinning a thousand times a second cannot have been born that way and stayed that way, because its own radiation would have slowed it in a few million years. It got there by being fed, and the line it cannot lie above is where the accretion torque balances the magnetic one.

 Two bodies at a mass ratio of 2.6 to 1. Both bodies orbit their common centre of mass, on similar ellipses whose sizes are in inverse proportion to the masses — here 2.6 to 1, so the heavier body's path is 2.6 times smaller.

The companion survives, and it is moving

When one star of a close pair explodes, the other is left holding an orbital velocity and no orbit. Whether the pair stays bound turns on a single question — did more than half the total mass leave? — and the stars that answer no are running across the Galaxy at a hundred kilometres a second with nothing visible behind them.

 Lane–Emden solutions for n = 0, 1, 1.5, 3, 4.5, 5, and the one that has no surface. The dimensionless density θ against the dimensionless radius ξ, for polytropic indices 0, 1, 1.5, 3, 4.5, 5. Each curve is the whole structure of a star whose pressure is K times its density to the power 1 + 1/n: the equation of state and hydrostatic equilibrium leave one second-order differential equation, and this is its solution. Every curve starts at θ = 1 with zero slope, because the density is greatest at the centre and has no cusp there. What separates them is where they end. At n = 0 the density is uniform and the surface is at ξ₁ = 2.4495; by n = 3 it has moved out to 6.8968 and the central density is 54.2 times the mean. At n = 5 the curve reaches zero only at infinity — a configuration of infinite radius and, remarkably, finite mass — and every index above it has neither. The three curves that have closed forms, n = 0, 1 and 5, are drawn from the same numerical integration as the rest and agree with those forms to better than two parts in a million, which is what licenses reading the others off the picture. What the figure cannot show is the scale: ξ is radius divided by a length that depends on the central density and on K, so two stars of the same index and wildly different sizes have the same curve here.

An equation of state is already a star

Write down how a gas's pressure depends on its density, insist that the pressure hold the weight up, and everything else follows — the run of density, the fraction of the mass inside each radius, and, at one particular index, a mass that does not care what the radius is.

 Mestel cooling for white dwarfs of 0.4 to 1 solar masses. Luminosity against age for white dwarfs of 0.4, 0.6, 0.8, 1 solar masses, both axes logarithmic. Each line is the whole of the star's later life: there is no nuclear source, so what is radiated is the thermal energy the ions had when the star was made, leaking out through a thin envelope whose opacity is Kramers'. That one sentence gives L ∝ t^(−7/5), and the lines are straight here because a power law is straight on these axes — the slope is the exponent and can be measured off the picture. Two features are worth more than the numbers. The first is how flat the sequence becomes: a factor of ten in luminosity costs a factor of 5.2 in age, so a white dwarf spends most of its existence faint, and any population of them piles up at the bottom. The second is that the heavier tracks lie above the lighter ones. A heavier white dwarf is smaller, so it has more ions to cool and less surface to lose them through, and at a given age it is the brighter object. The picture does not show what happens at the faint end, where the ions crystallise and the heat capacity stops being constant; that is a separate figure, and it is where this law stops being enough.

A clock with no fuel in it

A white dwarf has nothing left to burn, so its brightness is a record of how long it has been cooling. The faintest ones in the Galaxy are not faint because they are small — they are as faint as anything has had time to become, and the place where they stop is a date.

 The Sun's internal rotation: differential above 0.693 R, rigid below, and the 15 nHz between two modes. Angular velocity against fractional radius, in nanohertz, from the inversion of hundreds of thousands of measured rotational splittings. Two features had to be discovered rather than deduced. The convection zone rotates differentially in latitude — 473 nanohertz at the equator, 330 at sixty degrees — and that latitude dependence persists all the way down through it, in surfaces that are very nearly radial rather than the cylinders a rotating convective fluid was expected to produce. And below about 0.693 of the radius the latitude dependence stops: the radiative interior turns as a single rigid body at 430 nanohertz, which is a remarkable thing for a fluid with no strength to do, and requires something — most likely a weak internal magnetic field — to be enforcing it. The two regimes are joined by a shear layer a few per cent of the radius thick, the tachocline, and it is where the solar magnetic cycle is generally thought to be generated, because it is the only place in the Sun with the shear a dynamo of that strength needs. The horizontal marks are what individual modes would report: the rotation averaged over the cavity each one occupies, weighted by the time the wave spends at each radius. A mode of degree 100 is trapped near the surface and reports 473 nanohertz; one of degree 1 passes through the deep interior and reports 458. The 15-nanohertz difference between them is the measurement, and splittings are measured to about a nanohertz — which is why the profile can be resolved at all. The kernels used here have the right support and the right sense; a real inversion uses computed eigenfunctions, and its resolution below 0.2 R is poor for the reason the previous figure gives.

A fluid that turns as one piece

The Sun's surface turns faster at its equator than at its poles, and everyone expected the inside to do something similar. Below seven-tenths of the way down it does not — the radiative interior rotates as rigidly as a bell, and nothing in hydrodynamics makes a fluid do that.

 Central pressure bracketed without a model: 6 bodies, 23 decades apart. What can be said about the middle of a body from its mass and its radius alone. The lower end of each bar is GM²/8πR⁴, which follows from hydrostatic equilibrium and nothing else — no equation of state, no composition, no temperature, no assumption whatever about how the density is arranged inside. The upper end costs one more assumption, that the density does not increase outward, and it needs a central density, which is a model output rather than an observation and is why that edge is drawn as the softer one. The dot is what a full structural model gives. For the first five bodies every dot lies inside its bar, and what is worth noticing is how wide the bar is: Sun's rigorous floor is 4.48e+13 pascals against a modelled 2.34e+16, a factor of 522. The bound is true and nearly useless there, because most of a centrally condensed body's pressure comes from the concentration and the derivation deliberately knows nothing about it. The relativistic entry is the exception, and the reason to draw the figure at all. neutron star's modelled central pressure is 8.5 times the Newtonian ceiling — a body no Newtonian arrangement of matter with density falling outward can produce. The floor still holds, and holds for a statable reason: relativity makes the pressure gradient steeper than Newtonian gravity does, so the true central pressure can only exceed what the Newtonian derivation demands. The bracket therefore does more than constrain an interior. Applied at a small enough radius it breaks, and where it breaks is where Newtonian hydrostatics has stopped being the right equation.

A floor under the centre that assumes nothing

There is a lower bound on the pressure at the centre of any body in hydrostatic equilibrium, and it needs no equation of state, no composition and no temperature — only a mass and a radius. For the Sun it is nearly useless. For a neutron star it says which theory of gravity the interior needs.

 Where each mode turns back: ℓ = 0 through the centre, ℓ = 300 in the outer 2 per cent. Why a set of frequencies is a depth profile and a single frequency is not. An acoustic wave travelling into a star meets a rising sound speed and is refracted back; it turns where its horizontal phase speed matches the local sound speed, which happens at c(r)/r = 2πν/√(ℓ(ℓ+1)). The horizontal axis is the angular degree on a logarithmic scale and the vertical axis is the fractional radius of that turning point, drawn at 2000, 3090, 4000 microhertz. The ordering is the content. A radial mode, ℓ = 0, has no horizontal phase speed at all and passes straight through the centre. Degrees one and two turn deep in the core. By ℓ = 300 the mode is trapped in the outer 2 per cent and knows nothing about anything below. So a frequency measured to a part in ten thousand constrains an average of the interior weighted in a way the mode itself decides, and measuring thousands of modes of different degree gives thousands of differently weighted averages — which is a solvable inverse problem, and is how the base of the convection zone was located at 0.713 of the radius rather than assumed. The sound speed here is a polytrope's rather than a tabulated solar model's, so the curve is the right shape and the wrong star in its outer tenth, where the real Sun is convective and this one is not. What the picture cannot show is the frequency dependence at fixed degree, which is weaker but not negligible: a higher-frequency mode of the same degree turns slightly deeper, and the three curves separating toward the right is that effect.

Every note turns back at its own depth

A sound wave heading into a star is refracted by the rising sound speed and turns around before it reaches the middle. Where it turns depends on which mode it is — so a list of frequencies is not one average of the interior but thousands of differently weighted ones, and that is a solvable problem.

 14 glitches, and the 1.5 per cent of Vela that is not slowing down. Accumulated fractional spin-up against time for a Vela-like pulsar over 40 years, in parts per million. The underlying spin-down has been removed, so a perfectly braking pulsar would be a flat line at zero. What is drawn instead is a staircase: 14 sudden jumps of a few parts per million, each rising in less than a minute and then relaxing partway back over a couple of months, leaving a permanent step behind. Nothing outside the star can deliver angular momentum on a timescale of seconds, so the source is internal, and the only internal component that could have any to give is one that has not been slowing down with the rest. That is a neutron superfluid: it carries its rotation in quantised vortices, the vortices pin to the crustal lattice and cannot migrate outward, and so the superfluid keeps the spin it had while the crust brakes past it. The reservoir grows until the pinning fails somewhere, and a glitch is the unpinning. The straight line through the staircase is the glitch activity, 0.68 parts per million per year, and it converts directly into an interior measurement: the crust cannot on average take more than the superfluid stores, so the decoupled component must hold at least 2τ_c times the activity of the star's moment of inertia, which here is 1.5 per cent. It is a lower bound rather than a value, and it is one of the very few quantitative statements about the inside of a neutron star that needs no equation of state at all. The individual glitch times and sizes here are drawn from a seeded generator rather than a catalogue; what is real is the staircase's shape, the partial healing, and the arithmetic that turns a slope into a fraction.

The part of a star that never slowed down

A pulsar's spin decays smoothly for years and then jumps upward inside a minute. Nothing outside can deliver angular momentum that fast, so something inside has been storing it — and the rate at which the jumps accumulate is a lower bound on how much of the star is not braking with the rest.

 A lever 12 radii long, and the spin it removes. A magnetised stellar wind, drawn with the Alfvén surface at 12 stellar radii — the schematic distance at which the wind's inertia finally beats the field. Inside it the gas is forced to turn with the star, so every gram that leaves carries the specific angular momentum of the radius at which it broke free rather than of the surface it came from, and the lever squares: J̇ = (2/3) Ṁ Ω r_A². Beyond the surface the streamlines curve backwards, because angular momentum conservation makes the azimuthal speed fall as 1/r while the radial speed does not. With a moment of inertia coefficient of 0.073 and a mass loss of 2.3·10⁻¹⁴ solar masses a year, the star loses a fraction 2.3·10⁻¹⁴ of its mass and a fraction 3·10⁻¹¹ of its angular momentum in the same year — a ratio of 1,315, which is (2/3)(r_A/R)²/k² and nothing else. The e-folding time for the spin is 3.3·10¹⁰ years against 4.3·10¹³ years for the mass. Nothing here is to scale in one respect that matters: the wind's density falls by more than ten orders of magnitude across the drawn region, so the streamlines are drawn as though the flow were visible when almost none of it is.

A wind that takes no mass and all the spin

The Sun loses about a ten-thousandth of itself to its own wind over its whole life. It loses most of its rotation to the same wind, and the whole of that asymmetry is one geometric factor — the gas is forced to keep turning with the star until it is a dozen or more radii out from the surface it left.

 A period, a colour, and an age. Rotation period against colour for stars of 125, 625, 1000, 2500, 4570 million years, under the empirical relation P = t^0.5189 × 0.7725(B−V − 0.4)^0.601. The isochrones do not cross and are separated at every colour by exactly the age ratio raised to 0.5189, which is what allows a single measured period to be inverted for an age once the colour is known. The Sun, at B−V = 0.653 and 4570 million years, is placed by the relation at 26.8 days against the 25.4 days it is observed to have. Three clusters are marked at a common colour to show the spacing directly. The dashed boundary at the left is where the Rossby number — the period divided by the convective turnover time — passes 2 on the oldest isochrone, at B−V = 1.35: past that point the braking weakens and the relation is known to over-predict the age, which is the one place a rotation period stops being a clock. What the figure cannot show is the scatter, which is a few days at fixed colour and age and is the real error bar on any single star.

The clock that starts by forgetting

An ordinary star tells nothing about its age. It sits on the main sequence for billions of years at almost fixed brightness and colour, and the one property that changes monotonically is how fast it turns — but only because the braking law destroys the initial condition first, and only until it stops.

 The Sun's internal rotation: differential above 0.693 R, rigid below, and the 15 nHz between two modes. Angular velocity against fractional radius, in nanohertz, from the inversion of hundreds of thousands of measured rotational splittings. Two features had to be discovered rather than deduced. The convection zone rotates differentially in latitude — 473 nanohertz at the equator, 330 at sixty degrees — and that latitude dependence persists all the way down through it, in surfaces that are very nearly radial rather than the cylinders a rotating convective fluid was expected to produce. And below about 0.693 of the radius the latitude dependence stops: the radiative interior turns as a single rigid body at 430 nanohertz, which is a remarkable thing for a fluid with no strength to do, and requires something — most likely a weak internal magnetic field — to be enforcing it. The two regimes are joined by a shear layer a few per cent of the radius thick, the tachocline, and it is where the solar magnetic cycle is generally thought to be generated, because it is the only place in the Sun with the shear a dynamo of that strength needs. The horizontal marks are what individual modes would report: the rotation averaged over the cavity each one occupies, weighted by the time the wave spends at each radius. A mode of degree 100 is trapped near the surface and reports 473 nanohertz; one of degree 1 passes through the deep interior and reports 458. The 15-nanohertz difference between them is the measurement, and splittings are measured to about a nanohertz — which is why the profile can be resolved at all. The kernels used here have the right support and the right sense; a real inversion uses computed eigenfunctions, and its resolution below 0.2 R is poor for the reason the previous figure gives.

A shear layer that should have spread

The Sun's convection zone turns differentially — its equator laps its poles about once every three months — and the radiative interior below turns as one rigid piece. Between them is a transition four per cent of the radius thick. Nothing in hydrodynamics keeps a velocity discontinuity that thin for four and a half billion years.

 Two measured numbers, and everything else on the page derived from them. 25 pulsars in the plane of period against period derivative, at their catalogued values. Only the two axes are measurements; the three families of contour are models. Constant surface field runs at slope −1 because B ∝ √(PṖ), constant characteristic age at slope +1 because τ = P/2Ṗ, and the two families cross the population at right angles — which is why a single dot fixes both. The Crab sits at 3.8·10¹² G and 1257 years, and its true age is 972; the millisecond pulsars at the lower left have fields ten thousand times weaker and characteristic ages of billions of years, because they were spun back up by a companion long after they died. The line at the lower right is the death line, B/P² below which the model says no pair production and therefore no radio emission — and J2144−3933 is drawn below it, an 8.5-second pulsar that is radiating anyway.

A neutron star born turning too slowly

Collapse an iron core a few thousand kilometres across down to twelve, and conservation of angular momentum multiplies its rotation rate by about eighteen thousand. A model with no transport in it delivers a newborn pulsar at the break-up limit; the ones that are observed turn twenty times slower, which is a measurement of the core before it fell.

 A hundredfold expansion is a four-order slowdown. Equatorial rotation speed against radius for three stars leaving the main sequence at 2, 10, 100 km/s, on the single assumption that nothing exerts a torque. Both axes are logarithmic. The specific angular momentum is held fixed along each track and checked at four radii rather than asserted, so the speed falls as the reciprocal of the radius and then faster, as the moment of inertia coefficient slides from 0.073 on the main sequence to 0.02 in a centrally condensed giant envelope. A star at 100 km/s crosses the 8 km/s line — the usual boundary for calling a giant a rapid rotator — at 45.6 solar radii, and everything larger is slower. That is why the observed giants are almost all under two kilometres a second, and why the one or two per cent that are not cannot be explained by anything the star did on its own: the angular momentum has to have been delivered, by a swallowed companion or a merger. The figure assumes no mass loss, which for the largest radii drawn is the weakest of its assumptions.

A surface that slowed because the star grew

A star leaving the main sequence expands by a factor of a hundred. Conserve its angular momentum and its equator slows by four orders of magnitude, which is why every red giant is a slow rotator. One or two per cent are not, and a giant turning at fifty kilometres a second has to have been given the angular momentum by something else.

 The field a collapse would arrive with, and the one it has. Field strength against hydrogen density, both logarithmic, over the eight decades between diffuse gas and a protostellar core. The steeper line is what perfect flux freezing demands: a sphere collapsing conserves both mass and flux, so B goes as R⁻² while ρ goes as R⁻³, and therefore B goes as the two-thirds power of the density exactly — an exponent with no free parameter in it. The shallower locus is what Zeeman measurements find: a flat branch at about 10 µG up to 300 per cubic centimetre, where the density is rising and the field is not, and a rise as the 0.65 power of the density above it. By the density of a core the two differ by a factor of 1, and a star built at the frozen-flux value would carry a field four orders of magnitude beyond anything measured on one. The flat branch is the important half: it says the gas is moving along field lines without dragging them, which is what gravity does to a cloud that is still magnetically supported, and it locates where the freezing has to break.

A field that would have arrived ten thousand times too strong

A collapsing cloud carries its magnetic field with it, and the arithmetic of that is not negotiable — the field grows as the two-thirds power of the density. Run it from a diffuse cloud to a protostellar core and the answer is four orders of magnitude above anything measured on a young star.

 A spot is dark because it is squeezed. Pressure against depth below the quiet photosphere's optical surface, drawn as a ratio to the pressure there. The rising curve is the surrounding gas, which grows exponentially with a 140-kilometre scale height because that is what hydrostatic equilibrium in an ideal gas produces. The flat pair of bands is the spot's own budget: a magnetic pressure of 3.58·10⁵ dyn/cm² from a 3000-gauss field, which is 25 per cent of the total, plus the gas pressure left over. Horizontal balance requires the two columns to reach the same total at the same geometric level, and the level at which they do is 348 kilometres below the quiet surface — so the spot's own optical surface sits in a hollow. Measured Wilson depressions, obtained from the foreshortening of a spot near the limb, are four to six hundred kilometres. Nothing about the darkness was assumed: a field strength read off a Zeeman splitting fixes the magnetic share, the share fixes the depression, and the depression fixes the temperature the deeper layer must have to carry the reduced flux.

The darkness a field pays for

A sunspot is not cool because something is missing. It is cool because a three-thousand-gauss field supplies part of the pressure that holds the column up, the gas therefore supplies less, and the level at which that gas becomes opaque sits some hundreds of kilometres deeper than the surface around it.

 One step of memory, and it is kept at the poles. Two predictors of a solar cycle's amplitude, each tested against the seven cycles for which both quantities exist, with amplitudes scaled so that cycle 21 is one. Left: the polar field measured at the minimum before the cycle begins, which correlates with what follows at r = 1.00. Right: the amplitude of the previous maximum, which correlates at r = 0.34 — that is, not at all. The cycle is therefore not a pendulum with momentum; it is a process with exactly one state variable, and the variable is the poloidal field that the decay of the previous cycle's spots leaves behind at high latitude. That is the observational core of the flux-transport picture: spots emerge tilted, their trailing polarity drifts poleward, and what accumulates there is the seed the next cycle's shear will wind up. It is also the only prediction in solar physics with a lead time of years, and it was what said, correctly and unpopularly, that cycle 24 would be the weakest in a century.

One step of memory kept at the poles

The amplitude of a solar maximum tells almost nothing about the next one. The polar field measured at the minimum in between tells a great deal — so the cycle is not a pendulum with momentum but a process with exactly one state variable, and the variable is a field of a few gauss at latitudes nobody can see well.

 A formula everyone uses, and the number no pulsar has. Above: the braking index measured for the 4 pulsars whose spin-down has been followed long enough to give a second derivative, against the value a magnetic dipole rotating in vacuum requires. That value is exactly 3, and it is what every catalogued field strength and every characteristic age assumes. Not one measurement reaches it: they run from 1.4 to 2.839, and all of them fall short in the same direction, which is the signature of a systematic rather than of noise. Below: what that costs. The age a spin-down history gives is the period divided by (n − 1) times its derivative, so the ratio to the quoted characteristic age is 2/(n − 1) — 5.00 for Vela. The numbers are not thereby useless: an exponent recovered from the data is exactly the kind of correction a measurement can absorb. What has gone is the claim that the field strength printed beside a pulsar is a measurement of a field. It is a measurement of a spin-down rate, read through a model the same pulsar refutes.

The exponent no pulsar has

Every field strength in the pulsar catalogue comes from one formula, which assumes the star is a magnetic dipole rotating in a vacuum and therefore that its spin-down obeys an exponent of exactly three. Where that exponent has been measured it is 2.51, 2.84, 1.4 — never three, and always short.

 One dimensionless number decides whether a cloud may collapse. Field strength against hydrogen column density, with three loci of constant mass-to-flux ratio. The ratio of mass to magnetic flux is conserved under flux freezing, so it cannot be changed by anything that happens during a collapse — which is what makes it a criterion rather than a description. Its critical value is the pure constant 1/(2π√G), and dividing by that gives the dimensionless λ drawn here: below λ = 1 the field can hold the cloud up for ever, however cold it gets, because gravity and the magnetic force scale the same way with radius; above it no field strength suffices. Each locus is a straight line of slope exactly one, because at fixed λ the required field is exactly proportional to the column. Zeeman measurements of dense cores put them a little above the line and their envelopes a little below it, which is the arrangement a slow leak of flux out of the centre would produce and is the observational case for ambipolar diffusion.

A threshold with no free parameter in it

Divide a cloud's mass by the magnetic flux threading it. Gravity and the magnetic force both fall as the inverse square of the radius, so the ratio cannot change during a collapse — and the critical value that separates a cloud which must collapse from one that never can is one over two pi root G, a pure constant.

 Three causes, three shapes, one residual. An eclipse-timing residual: the observed minus the computed time of each eclipse, against cycle number, in days. Three effects are superposed and each has its own functional form. A slow change in the orbital period — from mass transfer or from magnetic braking — integrates to a parabola. A third body in a wide orbit moves the whole binary towards and away from the observer, so its light-travel time adds a sinusoid at the third body's period. And an eccentric orbit whose apsides are precessing moves the two eclipses in opposite directions, which is a sinusoid that changes sign between primary and secondary minima. The last is why both eclipses have to be timed: a third body moves them together and apsidal motion moves them apart, and a series of primary minima alone cannot tell the two apart at all.

Three causes with three shapes in one curve

The times of eclipse in a binary star are a clock, and the clock runs late and early. Three completely different things make it do so — a third body, a changing period, and a slowly turning orbit — and they are separable only because each imposes a different shape on the residual.

 A high end that is 2.5 in slope steeper than the one stars were born with. The mass function a star count measures against the one stars were born with, for a population that has been forming stars steadily for 10 billion years. Below about a solar mass nothing has had time to die, so the two coincide exactly. Above it the fraction still alive is the main-sequence lifetime divided by the age, and since the lifetime falls as the two-and-a-half power of the mass, the present-day function is steeper than the initial one by exactly that exponent. A count of massive stars in an old population therefore under-represents them by orders of magnitude, and reading it as an initial mass function gives a slope far too steep. The correction is large, it is calculable, and it depends on the star-formation history — which is usually the thing the mass function was going to be used to constrain.

A mass function corrected by an age

Counting stars by mass gives the stars that are alive. What every argument needs is the stars that were born, and the two differ by the fraction of each mass still on the main sequence — which is a lifetime divided by an age, and which for massive stars in an old population is a very small number.

 4 per cent in one observable is 48 per cent in an age. How an error in the calibration of the large frequency separation propagates into the quantities derived from it. The two scaling relations are exact in their exponents, so a fractional error in the separation appears as twice that in the radius, four times in the mass, and — because a main-sequence lifetime falls as roughly the two-and-a-half power of the mass — ten times in an age. At the 4 per cent level, which is about what the theoretical corrections to the relation amount to for a red giant, that is 15 per cent in mass and 48 per cent in age. Nothing about the seismology is uncertain at that level; the frequencies are measured to parts in a thousand. What is uncertain is the constant of proportionality, and it is uncertain because it was calibrated on one star.

Two scaling relations calibrated on one star

Asteroseismology gives a star's mass and radius from two numbers read off its oscillation spectrum. The two relations are exact in their exponents and approximate in their constants, and the constants were fixed by requiring that the Sun come out right — so an error of a few per cent in one observable is tens of per cent in a mass and nearly a factor in an age.

 A free parameter worth 88 kelvin across its plausible range. The effective temperature a stellar model predicts, against mass, for three values of the mixing-length parameter. The parameter has no derivation: it is the distance a convective blob is supposed to travel before dissolving, in units of the local pressure scale height, and it is fixed by requiring that a model of the Sun reproduce the Sun. The three curves span 88 kelvin, which at fixed luminosity is a radius difference of 1.5 per cent — comparable to the precision with which radii are now measured by interferometry and by eclipsing binaries. Every stellar age, every isochrone and every mass inferred from a position in the temperature–luminosity plane depends on the value chosen, and there is no reason beyond convenience to expect the solar value to apply to a red giant or to a metal-poor dwarf.

A length nobody derived, fitted to one star

Convection in a star is turbulent, three-dimensional and impossible to compute inside an evolution code. What is used instead is one number — how far a blob of gas travels before dissolving — fixed by requiring that a model of the Sun come out with the Sun's radius, and then applied to every star ever modelled.

 A free parameter worth 42 per cent of an age. The main-sequence lifetime of a star against its mass, for four values of the convective overshoot parameter. A convective core has a boundary where buoyancy vanishes, and a rising blob arriving there still has momentum, so it penetrates into the stable region above and mixes fresh hydrogen into the core. How far it penetrates is not computed — it is a parameter, quoted in pressure scale heights, and the values in current use range from zero to about a third. A larger core is a larger fuel supply, so the lifetime rises and the turn-off at a given age is more massive. Across the plausible range the lifetime of a two-solar-mass star changes by 42 per cent, and an age read off a cluster's turn-off changes by the same amount. That is larger than the quoted uncertainty on almost every cluster age in the literature.

A convective boundary with no theory to fix it

A convective core has an edge where buoyancy vanishes. A blob arriving there still has momentum, so it carries on, mixing fresh hydrogen into the core and extending the star's life. How far it carries on is a fitted parameter, and across its plausible range every stellar age changes by nearly half.

 A flux that measures a temperature to 0.12 per cent. Neutrino flux against central temperature, in units of the standard model's, for the three main solar channels. The exponents are not arbitrary: each reflects how far up the Gamow peak the reaction has to reach, so the channel with the largest Coulomb barrier is the steepest. The boron-8 flux goes as roughly the twenty-fourth power, which means a measurement good to 3 per cent constrains the Sun's central temperature to 0.12 per cent — better than any other technique by an order of magnitude. The same steepness is why the flux is useless as a check on anything else: a stellar model whose central temperature is uncertain at the half-per-cent level predicts this flux to within a factor, and the disagreement between two model families is far larger than the measurement.

A flux that is a thermometer to a tenth of a per cent

The boron-8 neutrino flux from the Sun's core rises as roughly the twenty-fourth power of the central temperature. That makes it the sharpest thermometer in astrophysics and simultaneously the most fragile prediction — a model uncertain in its central temperature by half a per cent predicts the flux to within a factor.

 An exponent between 2.3 and 4, not a number. The local logarithmic slope of the mass–luminosity relation against mass — the exponent that "L ∝ M^3.5" stands in for — read off the same fit to eclipsing binaries that the relation itself is drawn from. It is a step function with two corners, and the steps are 2.3 below 0.43 M☉, 4 0.43 to 2 M☉, 3.5 2 to 55 M☉. The horizontal lines are what homology predicts: a star in radiative equilibrium has L ∝ μ⁴M³/κ̄, which gives an exponent of 3 when the opacity is electron scattering and a constant, 5.5 when it is Kramers and depends on density and temperature, and 1 where radiation pressure holds the star up and the luminosity is pinned to the Eddington value. The measured steps sit between those predictions rather than on them, which is what a real star does — no star is homologous, the low-mass ones are convective throughout and the heavy ones have convective cores, and the fitted exponents are what is left after all of that. The useful statement is not that the exponent is 3.5 but that it is between 2 and 5 everywhere and 4 in the middle, and that the reason it moves is a change of which opacity is carrying the energy out.

An exponent that is a slope, not a law

The phrase “L goes as M to the three and a half” stands in for a curve with three straight pieces and two corners. The exponent is 2.3 below half a solar mass, 4 in the middle and falls towards 1 at the top — and each change of exponent is a change in which opacity is carrying the energy out.

 A 1.1 solar-mass white dwarf held up for 2.8 extra billion years. The time a white dwarf takes to reach the crystallisation luminosity, and the two delays that follow, against mass. The lower band is bare Mestel cooling — thermal energy of the ions leaking out through an envelope whose opacity is Kramers'. On top of it sits the latent heat of crystallisation, 0.85 kT per ion released when the liquid interior freezes into a lattice; and on top of that the gravitational energy of ²²Ne settling through what is left, taken here as 0.6 kT per ion. Neither is fuel: both are energy the star already had, released late and radiated at the low luminosity it has by then, which is why so little of it buys so much time. The delay rises from 1.88 billion years at 0.5 solar masses to 2.81 at 1.1 — 58 per cent of the cooling already done. A white dwarf age computed from the bare law is too young, and it is too young by more the heavier the star is, which is exactly the direction that matters, because the massive white dwarfs are the ones used to date the oldest populations.

A clock that stops while its interior freezes

A white dwarf has nothing left to burn, so its brightness is a record of how long it has been cooling — until the interior crystallises. The latent heat of that phase change, and the settling of a heavy isotope through what is left, hold a massive white dwarf up for nearly three extra billion years.

 The same equation, run away at q = 0.8 and settled at q = 0.5. The donor's overfill of its own Roche lobe against time, integrated for 3 mass ratios at a donor adiabatic response of −0.33 — a star with a deep convective envelope, which expands as it loses mass. The overfill sets the transfer rate and the transfer rate changes the overfill, and the sign of that feedback is the stability criterion: where the lobe shrinks faster than the star does, ζ_L > ζ_ad, the overfill grows and the growth is exponential. At q = 0.8 the lobe responds at 0.03 and the overfill runs away; At q = 1.4 the lobe responds at 1.32 and the overfill runs away. At q = 0.5 it responds at -0.62 and the transfer throttles itself back. The critical ratio for this donor is 0.63 — below equal masses, which is the result the whole subject turns on: a giant transferring to a lighter companion is already unstable before the mass ratio has reversed, and what follows is not accretion but a common envelope. Nothing here is a fitted rate. The vertical scale is logarithmic and the runaway is a straight line on it, which is what an exponential is.

Three clocks and a runaway

Whether mass transfer between two stars is stable is a comparison of two logarithmic derivatives. How fast it runs is a separate question with three possible answers fourteen orders of magnitude apart — and the answer decides whether the companion accretes, is buried, or is swallowed.

The ladders in this field

26 anchors · one idea each

The HR diagramThe mass–luminosity relationStellar lifetimesFusionDegeneracyHydrostatic equilibriumVariable starsBinary starsEddington limitEnergy transportStellar evolutionAsteroseismologyInitial mass functionSupernovaeAccretionPulsarsSolar cycleSolar neutrinosMass transferPolytropesWhite dwarf coolingInternal rotationMagnetic brakingFlux freezingPlasma betaStar formation

All fields · All essays