Astronomy & Celestial Mechanics
the sky, computed
Series

Line formation — the series

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

    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.

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

    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.

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

    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.

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

    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.

    part 4 · starlight
  5. The same two lines, seen as frequencies. The amplitude of the Fourier transform of the two line profiles of the previous comparison, on a logarithmic scale, against frequency in cycles per kilometre per second. At low frequency — the overall width — they agree, as the profiles did. At higher frequency they separate completely: the first star's transform falls to zero at 0.110 cycles per km/s and the second's at 0.132, because the rotation kernel's transform vanishes at 0.660/v sin i and no Gaussian-like broadening has a zero anywhere. The horizontal dashed line marks a noise level of 0.0005 in the transform, which a high signal-to-noise spectrum averaged over several lines reaches; the dip and the side lobe beyond it rise above it for both stars, so the rotation of each is measured by where its transform dips, independent of the macroturbulence. In a noisier spectrum the second star's weaker side lobe would be lost first.

    A rotation found where the transform goes dark

    A star's rotation and the turbulent flows on its surface both widen its spectral lines, and a profile fitted for one can be refitted for the other with almost no change. Convolution becomes multiplication in Fourier space, though, and the transform of a rotating disc has zeros — frequencies where it is exactly nothing — while the transforms of turbulence and of the spectrograph have none. Where the transform of a line goes dark is a measurement of rotation that turbulence cannot fake.

    part 5 · starlight

All series