Astronomy & Celestial Mechanics
the sky, computed
Generator

The hr-diagram generator

The Hertzsprung–Russell diagram
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.

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.

12 essays call hr-diagram. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

Blackbody curves at 3000, 5800, 10000 K. Thermal emission against wavelength, each curve scaled to its own peak so the shift can be seen on one plot. The peak moves to shorter wavelengths as the temperature rises, which is why colour is a thermometer. Starlight

Colour is a thermometer, and it reads across the galaxy

A star's colour gives its surface temperature, from two brightness measurements and no other information. It is the cheapest useful measurement in astronomy.

 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. Stars

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.

 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. Stars

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.

 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. Stars

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.

 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. Stars

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. Stars

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.

 Hard below 8.3 astronomical units, soft above, and nothing settles at the line. The binding energy of a binary of two 0.7 solar-mass stars against its separation, both axes logarithmic, with the mean kinetic energy of a single cluster star at a velocity dispersion of 5 kilometres a second drawn as a level. Where the curve is above the level the binary is bound more tightly than a passing star's motion, and encounters on average take energy out of the field and put it into the pair; where it is below, they do the reverse. The crossing at 8.3 astronomical units is the hard–soft boundary, and it is the entire content of Heggie's law: hard binaries harden and soft binaries soften. The arrows are the direction each side moves, and they point away from the crossing in both directions rather than toward it. That is not a coincidence but a negative heat capacity, the same property that makes a star contract when it radiates: taking energy out of a bound pair moves it closer together and speeds it up, so a hard binary that gives energy to the cluster becomes harder still and gives more. A cluster with binaries in it therefore has a heat source that turns itself up, and the boundary drawn here is a watershed rather than an equilibrium. Gravitation

A pair that heats what is trying to cool

A star cluster has a negative heat capacity, so it cannot reach equilibrium — its core contracts and gets hotter without limit. What stops it is a binary, and which way a binary exchanges energy with the stars around it is settled by one comparison of two energies.

 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. Stars

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.

 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. Stars

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 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. Stars

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.

 Isochrones that crowd on the main sequence and fan out above it. Isochrones for stars of solar metallicity at 1, 2, 4, 6, 8, 10, 12 billion years, luminosity against temperature with temperature falling to the right. On the main sequence they lie almost on top of one another — a star brightens and barely changes colour as it ages — and only where they turn off and cross towards the giants do they separate. Two boxes show the measurement errors of a star with a good parallax and spectrum (twice 70 K and 3 per cent): one around a Sun-like star of 4.6 billion years, whose age they allow anywhere in 4.53 (3.38–5.63) billion years; one around a subgiant of the same mass at 10.5, allowed only 10.5 (10.3–10.7). Stars

A lone star is dated only as it leaves

A star in a cluster is dated by the bend its neighbours trace. A star alone has one position on the Hertzsprung–Russell diagram, and an age for it means finding which isochrone passes through that point. On the main sequence they almost all do: a Sun-like star measured as well as a parallax and a spectrum allow is anywhere from three and a half to five and a half billion years old. A star the same mass on its way off the main sequence is dated to two per cent.

 Four clocks for a 1 M☉ dwarf, and the best of them. The fractional age uncertainty each method reaches for an isolated 1 solar-mass main-sequence star, against its age, logarithmic, for isochrones, rotation, oscillations and chemical clock, with representative published precisions. Rotation is best while the star is young, once its spin has forgotten its start and before its braking weakens; isochrones and oscillations improve with age as the star moves further on the diagram. The combination (heavy line) is 14 per cent at 1 Gyr, 9 at 3 and 4 at 8. It is better than ten per cent (dashed) from 2.6 to 9.4 Gyr. Stars

Four clocks and the stars none of them reads

An isolated dwarf star can be dated four ways — by its position on the Hertzsprung–Russell diagram, by how fast it spins, by its oscillations, and by the ratio of two elements its galaxy made at different rates. Each clock has an age range where it works and masses it cannot read. Laid side by side for a Sun-like star they cover most of its life to ten per cent or so. For a star of eight-tenths of the Sun's mass only one clock is left, it works only while the star is young, and the commonest stars in the Galaxy's disc are dated one at a time no better than fifteen per cent at any age.

The whole library · All essays