Astronomy & Celestial Mechanics
the sky, computed
Series

Chemical evolution — the series

4 essays on one idea, from the one that introduces it to the one that assumes the rest.
  1. A one-parameter model, and 9 times too many metal-poor stars. The metallicity distribution of long-lived stars near the Sun — the bars — against two models of how a galaxy enriches itself. The closed box is as simple as a model of a galaxy can be: gas turns into stars, stars make metals and return them, nothing enters and nothing leaves. It has exactly one parameter, the yield, and it predicts the whole curve. It gets the peak roughly right and the tail catastrophically wrong: 13.7 per cent of its stars fall below [Fe/H] = -1 against an observed 1.6 per cent, a factor of 9. This is the G-dwarf problem, and the reason it is an argument rather than a discrepancy is that the failure cannot be fixed by changing the yield: the yield sets where the peak is, and moving the peak to fix the tail moves it away from the data. What is wrong is a boundary condition. Let gas keep arriving — pristine, at roughly the rate it is being consumed — and the gas is never both abundant and metal-poor for long, so few stars form while it is. That is the second curve, with 3.2 per cent below -1, and it needed no new nucleosynthesis and no new parameter beyond the fact of accretion. The picture cannot show the thing that would settle it directly, which is the infall itself: the gas arriving on the disc now is a few solar masses a year spread over twenty kiloparsecs, and it has never been securely detected.

    A histogram that says the box was not closed

    The simplest model of a galaxy enriching itself has exactly one parameter and predicts the whole metallicity distribution of its surviving stars. The solar neighbourhood's disagrees, in a specific direction — and the failure is not in the nucleosynthesis but in a boundary condition.

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

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

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

    part 4 · galaxies

All series