Astronomy & Celestial Mechanics
the sky, computed
Generator

The curve-of-growth generator

The curve of growth
The curve of growth. The equivalent width of an absorption line against the number of absorbers along the sight line, both logarithmic, computed by integrating a Voigt profile with damping parameter 0.005. Three regimes: the width grows in proportion to the abundance while the line is weak, then almost not at all for two decades once the core saturates, then as the square root once the damping wings dominate. A line measured in the middle stretch carries almost no information about the abundance, and most strong lines in a stellar spectrum are there.

The equivalent width of an absorption line against the number of absorbers along the sight line, both logarithmic, computed by integrating a Voigt profile with damping parameter 0.005. Three regimes: the width grows in proportion to the abundance while the line is weak, then almost not at all for two decades once the core saturates, then as the square root once the damping wings dominate. A line measured in the middle stretch carries almost no information about the abundance, and most strong lines in a stellar spectrum are there.

5 essays call curve-of-growth. 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.

The curve of growth. The equivalent width of an absorption line against the number of absorbers along the sight line, both logarithmic, computed by integrating a Voigt profile with damping parameter 0.005. Three regimes: the width grows in proportion to the abundance while the line is weak, then almost not at all for two decades once the core saturates, then as the square root once the damping wings dominate. A line measured in the middle stretch carries almost no information about the abundance, and most strong lines in a stellar spectrum are there. Starlight

How much of a line is not in its depth

An absorption line stops getting deeper long before it stops getting stronger. What keeps growing is its area, and the way the area depends on the number of absorbers has three distinct regimes over five decades — one of which carries almost no information at all.

 Three profiles of equal equivalent width, 28.6 mÅ. Three absorption profiles with the same equivalent width — 28.6 milliångström, 1.72 km/s at 500 nm, matched to better than 0.1% by root-finding over the quadrature — differing in nothing but shape. Left: the cores, on a common velocity axis. Right: the same three normalised to their own half widths, on a logarithmic depth scale. The thermal profile is a Gaussian set by Fe's mass at 6000 K, 1.34 km/s; the collisional one a Lorentzian of γ = 2.28e-3 nm; the rotational one the classical kernel of a disc turning at v sin i = 1.73 km/s, which is exactly zero beyond 1.28 half widths and is the only one of the three with an edge. Matching the areas does not match the widths: the half widths are 1.11 km/s (thermal), 1.35 km/s (rotational), 0.68 km/s (collisional), a factor of 1.98 between the widest and the narrowest. At three half widths the collisional wing is 49 times the thermal one and at five it is 1.27·10⁶ times — six decades, which is why a line's shape stays diagnostic long after its width has stopped being so. 11% of the Lorentzian's own equivalent width lies beyond the right-hand panel's edge and is not drawn anywhere. Starlight

The same width for three different reasons

Thermal motion, rotation and collisions each widen an absorption line, and they can be tuned to areas that agree to a part in a thousand. What separates them is the shape, and the shape carries a rotation speed from one profile and a surface gravity from another.

 A wind that writes its own terminal speed on the blue edge, at 2017 km/s. A P Cygni profile, computed by integrating a spherical wind on a grid rather than drawn. The wind accelerates outward as one minus the inverse radius to the power 0.8, reaching 2000 kilometres a second, and its density falls as the inverse square of the radius over its own speed; the horizontal axis is the Doppler shift in units of that terminal speed, with blue to the left. The column directly in front of the stellar disc is moving toward the observer at every speed from nearly nothing near the surface up to the terminal value far out, so it takes light out of the beam across that whole range and the absorption trough runs from the rest wavelength to the blue edge at -1.01 of the terminal speed. That edge is the measurement: no model of the star, its distance or its mass enters it, only the geometry of a column seen end-on. The emission comes from everywhere else, where the wind moves across the line of sight and can only add photons; it is Infinity times stronger to the red than to the blue, because the part of the shell receding directly away is hidden behind the star and its blueshifted counterpart is not. Two features, one shell, and between them a velocity and a mass-loss rate. Starlight

A speed read off an edge

A hot star drives material off itself at thousands of kilometres a second, and the profile that wind prints on the star's own spectrum has a sharp blue edge. That edge is the terminal velocity, measured with no model of the star, no distance, and no calibration — one of the few numbers in stellar astrophysics obtained from geometry alone.

 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.

 A synthetic spectrum fitted with a line list that is missing the weak lines. Four ångströms of a simulated metal-line spectrum near 5000 Å at resolving power 20,000 and a signal-to-noise ratio of 200 (thin, noisy line), and the best-fitting synthesis (heavy line) computed from a line list that omits every line with a central optical depth below 0.15 — 137 of the 160 lines in the ten-ångström window, carrying 32 per cent of its total absorption. The fit adjusts one abundance for all the listed lines, with the continuum held at its true level; it returns an abundance 0.188 dex too high: the missing lines' absorption has nowhere to go but into the parameters the synthesis has. The residuals (below, offset) show the missing lines as dips that no parameter of the synthesis can fill — which is how an incomplete list announces itself at high resolution, and why it does not at low. Starlight

The lines a synthesis keeps pay for the ones it lacks

Modern abundance work fits a computed spectrum to a whole stretch of the observed one rather than measuring lines one at a time, and it averages away the errors in each line's laboratory data as it does so. It cannot average away a line that is not in its list. The absorption of the missing lines has to be put somewhere, and the fit puts it into the parameters it has — the abundance of the lines it knows, or the level of the continuum — in a direction set by which of the two it was allowed to move.

The whole library · All essays