Astronomy & Celestial Mechanics
the sky, computed
Series

Stellar evolution — the series

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

    part 1 · stars
  2. 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.

    part 2 · stars
  3. 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.

    part 3 · stars
  4. 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.

    part 4 · stars
  5. 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.

    part 5 · stars

All series