Astronomy & Celestial Mechanics
the sky, computed
Generator

The saha-ionisation generator

The Balmer maximum is at 9,870 K, and it is a maximum in temperature
The Balmer maximum is at 9,870 K, and it is a maximum in temperature. What two absorption lines actually count, each normalised to its own maximum. The first curve is the fraction of all hydrogen sitting in n = 2 — the only hydrogen a Balmer line can absorb — which is a Boltzmann factor climbing with temperature multiplied by the neutral fraction falling with it. The product peaks at 9,870 K, where 35.2% of the hydrogen is still neutral and only 8.73·10⁻⁶ of all of it is in n = 2 at all. Below the peak there is plenty of hydrogen and almost none of it excited; above it there is plenty excited and almost none of it neutral. The second curve is the fraction of calcium that is singly ionised, which is what the Ca II K line counts, and it peaks at 6,335 K — cooler, because calcium gives up its first electron at 6.113 eV. At 6,000 K the Balmer curve is at 0.1% of its own maximum while Ca II is near its peak, and calcium is 4.5·10⁵ times rarer than hydrogen in the same gas. A spectrum in which Ca II K is the strongest line is not a spectrum of a calcium star. Reading it as one is precisely the error that had the Sun made of iron until 1925.

What two absorption lines actually count, each normalised to its own maximum. The first curve is the fraction of *all* hydrogen sitting in n = 2 — the only hydrogen a Balmer line can absorb — which is a Boltzmann factor climbing with temperature multiplied by the neutral fraction falling with it. The product peaks at 9,870 K, where 35.2% of the hydrogen is still neutral and only 8.73·10⁻⁶ of all of it is in n = 2 at all. Below the peak there is plenty of hydrogen and almost none of it excited; above it there is plenty excited and almost none of it neutral. The second curve is the fraction of calcium that is singly ionised, which is what the Ca II K line counts, and it peaks at 6,335 K — cooler, because calcium gives up its first electron at 6.113 eV. *At 6,000 K the Balmer curve is at 0.1% of its own maximum while Ca II is near its peak*, and calcium is 4.5·10⁵ times *rarer* than hydrogen in the same gas. A spectrum in which Ca II K is the strongest line is not a spectrum of a calcium star. Reading it as one is precisely the error that had the Sun made of iron until 1925.

8 essays call saha-ionisation. 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.

An absorption spectrum at 5772 K. A blackbody continuum at 5772 K with absorption lines cut out of it, each line's depth computed from how much of the gas is in a state that can absorb it. 8 of the 8 lines are strong enough to see at this temperature, which is the whole reason the spectral sequence is a temperature sequence. Starlight

Composition, read from what is missing

The dark lines in a stellar spectrum are wavelengths that never arrived. Which ones are absent names the elements present — and the strength of a line says more about temperature than about abundance.

 The Balmer maximum is at 9,870 K, and it is a maximum in temperature. What two absorption lines actually count, each normalised to its own maximum. The first curve is the fraction of all hydrogen sitting in n = 2 — the only hydrogen a Balmer line can absorb — which is a Boltzmann factor climbing with temperature multiplied by the neutral fraction falling with it. The product peaks at 9,870 K, where 35.2% of the hydrogen is still neutral and only 8.73·10⁻⁶ of all of it is in n = 2 at all. Below the peak there is plenty of hydrogen and almost none of it excited; above it there is plenty excited and almost none of it neutral. The second curve is the fraction of calcium that is singly ionised, which is what the Ca II K line counts, and it peaks at 6,335 K — cooler, because calcium gives up its first electron at 6.113 eV. At 6,000 K the Balmer curve is at 0.1% of its own maximum while Ca II is near its peak, and calcium is 4.5·10⁵ times rarer than hydrogen in the same gas. A spectrum in which Ca II K is the strongest line is not a spectrum of a calcium star. Reading it as one is precisely the error that had the Sun made of iron until 1925. Starlight

Why hydrogen's lines are strongest where hydrogen is not

A Balmer line counts the hydrogen atoms sitting in one particular excited state, and that population peaks near 10,000 K — where a third of the hydrogen has already been ionised away. The strength of a line is a thermometer, and reading it as an abundance put hydrogen at one per cent of the Sun until 1925.

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

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.

 A continuum drawn 4.3 per cent below the real one. A short stretch of spectrum with one strong line in it and 150 weak ones scattered across the same interval. The upper dashed line is the true continuum — the flux the star would emit with no lines at all — and it is not observable. The lower one is what a fit through the highest points of the spectrum returns, which is 4.3 per cent lower, because the weak lines have depressed the gaps between the strong ones. Measuring the strong line's equivalent width against the apparent continuum instead of the real one makes it 11.2 per cent too small. The error has a sign, it is worse in spectra with more lines, and it therefore correlates with metallicity — which is exactly the quantity being measured. Starlight

A continuum that was never observed

An equivalent width is an area measured relative to the continuum, and the continuum is not in the data. It is drawn — a curve through the highest points of the spectrum — and in any spectrum with many weak lines those highest points are already below the true continuum, because the weak lines have eaten the gaps.

 Two diagnostics, two bands, and a crossing to 49 kelvin. The plane of effective temperature against surface gravity, with the constraints from two spectroscopic diagnostics drawn as bands. The wings of a hydrogen line are broadened by collisions, so they respond steeply to the gravity and weakly to the temperature: a narrow, steep band. An ionisation balance — requiring that the same element give the same abundance from its neutral and its singly ionised lines — responds to both, and its band is much shallower. Neither diagnostic determines either quantity on its own. Where the two cross is the answer, and the size of the crossing region is set by the band widths divided by the difference of the slopes — so two diagnostics that respond similarly give a long, thin, nearly useless error region however precise each one is. Choosing diagnostics that disagree in their sensitivities is the whole of the art. Starlight

A temperature and a gravity that trade against each other

A stellar spectrum contains the star's temperature, its surface gravity and its composition, and no single feature in it contains only one of the three. Every diagnostic is a band in the parameter plane rather than a point, and the answer is where the bands cross — which makes choosing diagnostics that disagree in their sensitivities the whole of the art.

 742 kelvin between a dwarf and a supergiant, at the same ionisation. The fraction of calcium still neutral against temperature, for three surface gravities: a dwarf, a giant and a supergiant. The three curves are the same curve slid sideways. The Saha equation carries the electron pressure in its denominator, and the photospheric pressure follows the surface gravity as its square root — hydrostatic equilibrium gives a gas pressure of order g over the opacity, and the opacity in a cool star is set by the electrons themselves. So the pressure runs from 19 newtons a square metre in the dwarf down to 0.48 in the supergiant, a factor of 40, and the half-ionisation point moves from 4427 kelvin to 3685 — 742 kelvin apart. At fixed temperature the ionisation ratio goes as g^-0.50, exactly the square root the algebra requires. A star's luminosity class is a measurement of the density of its photosphere, read off which stage of an element the lines belong to. Starlight

A luminosity class is a density measurement

The Saha equation has an electron pressure in its denominator, and a supergiant's photosphere is forty times less dense than a dwarf's at the same temperature. So the same element sits in different ionisation stages in the two, the spectrum says which, and the second axis of stellar classification is a barometer.

 Where titanium oxide comes apart, and what the pressure does to it. The fraction of titanium oxide dissociated into Ti + O, against temperature, for a gas of equal parts at 10, 300, 6000 N/m² of total pressure. The equilibrium is Saha's with a free atom in place of the free electron and the pair's reduced mass in place of the electron's — the two atoms' partial pressures multiplied together, divided by the molecule's, equal an equilibrium constant Kₚ(T) — which for equal abundances gives a dissociated fraction of √(Kₚ/(Kₚ+P)) exactly. A molecule is one particle becoming two, so pressure suppresses the reaction exactly as it suppresses ionisation — the half-dissociation point moves from 3145 K at 10 N/m² to 4164 K at 6000. The bond energy is 6.87 eV against hydrogen's ionisation potential of 13.6, which is why this happens at a third of the temperature. Rotational and vibrational partition functions are included, classically and harmonically; the atoms keep their ground-term weights. Starlight

A spectrum with no continuum left in it

Below about four thousand kelvin a photosphere stops being a gas of atoms. Titanium oxide, water and carbon monoxide take over, the same equilibrium settles how much of each survives with a bond energy where an ionisation potential used to be — and the bands are so crowded that the level every line depth is measured against is nowhere on the plate.

 Where a photosphere's free electrons come from. The share of the free electrons donated by metals rather than by hydrogen, against temperature, at a gas pressure of 12000 N/m² and at solar abundance, 1 dex below solar, 2 dex below solar. The electron pressure is not assumed here — it is solved for, as the value at which the ionisation of the gas supplies exactly the electrons the gas contains. Below about 5897 K at solar abundance every free electron in the gas comes from an element present at one part in ten thousand, because hydrogen's 13.6 eV keeps it neutral while magnesium's 7.6 does not. The handover is fast: solar crosses a half at 5897 K, 1 dex down crosses a half at 5088 K, 2 dex down crosses a half at 4467 K. At 5,772 K and this pressure the solved electron pressure is 1.78 N/m². Starlight

The continuum is made by one atom in ten thousand

The Sun's light leaves through an ion that exists only because a hydrogen atom will hold a second electron by three-quarters of an electronvolt. The electrons it holds come almost entirely from magnesium, silicon and iron — so the level against which every solar line depth is measured is rationed by elements present at one part in ten thousand.

The whole library · All essays