Lecture 04: HI Column Density: From Spectra to the Atomic Gas Budget

ASTR 360 The Interstellar Medium

Context and Why This Matters

Lecture 03 derived the radiative-transfer physics of the 21 cm line; this lecture turns that physics into the single most-used quantitative tool of Galactic HI astronomy -- the column-density formula -- and applies it to a real, published measurement toward a well-studied extragalactic sightline.

Learning Goals

  1. Derive the optically thin HI column-density formula N(HI) = 1.8224e18 integral(T_B dv) [cm^-2] from the 21 cm radiative-transfer solution.
  2. Apply the formula to a representative two-component (CNM+WNM) 21 cm spectrum and to a real, published all-sky-survey column density.
  3. Explain the physical origin of the numerical coefficient 1.8224e18 in terms of the 21 cm transition’s atomic constants.

Opening Phenomenon

The HI4PI all-sky 21 cm survey measures a Galactic-foreground atomic hydrogen column density of N(HI)=1.79e+20 cm^-2 toward the bright, well-studied quasar 3C 273 -- a single number, derived entirely from a radio spectrum, that astronomers routinely use to correct X-ray and ultraviolet observations of that same distant quasar for foreground Galactic absorption.

Terms Developed in Context

column density

The number of particles per unit projected area integrated along the line of sight, commonly expressed in cm^-2. In this lecture, the term is not an isolated label: This concept is needed to reason about hi column density: from spectra to the atomic gas budget. Starting from Lecture 03’s optically thin solution T_B(v) ~= T_spin*tau(v), and the definition of 21 cm optical depth in terms of the lower-level column density and the Einstein A coefficient, the total (velocity-integrated) HI column density reduces to a simple integral of the observed brightness temperature over velocity, independent of the (unknown, but nearly canceling) spin temperature: N(HI) = C * integral(T_B(v) dv), where C = 1.8224e18 cm^-2 (K km/s)^-1 is a fixed combination of fundamental constants and the 21 cm transition’s atomic properties (the Einstein A coefficient, statistical weights, and hnu/k).

It is tied to the lecture’s observables or evidence in this way: The total number of hydrogen atoms per unit area along a line of sight (the column density N(HI), typically reported in cm^-2) is the single most useful summary statistic of an HI spectrum for nearly every downstream application, from Galactic structure studies to correcting extragalactic X-ray/UV absorption. A second connection is important for using the term in physical reasoning: For a Gaussian line profile with peak brightness temperature T_peak and full width at half maximum (FWHM) Delta_v, the velocity integral has the closed form integral(T_B dv) = T_peak * Delta_v * sqrt(pi/(4 ln2)), letting the column-density formula be evaluated directly from two observable quantities (peak and width) rather than requiring numerical integration of every spectrum. The observational interpretation is constrained by All-sky 21 cm surveys such as the Leiden/Argentine/Bonn (LAB) survey and its modern successor, HI4PI (HI4PI Collaboration 2016), provide column densities toward essentially any Galactic direction, including toward extragalactic background sources used for X-ray absorption studies (e.g., 3C 273, Cygnus X-1).

Boundary and common confusion. Applying the optically thin approximation uncritically to every CNM component. Narrow, cold, high-column CNM clumps can have tau approaching or exceeding order unity, in which case the true column density (from the full T_B=T_spin(1-e^-tau) solution) is larger than the optically thin formula would suggest; the assumption should always be checked, not assumed, particularly for narrow, high-peak-brightness features. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

optically thin approximation

The first-order limit in which optical depth is much less than one, allowing 1-e^-tau to be replaced by tau. In this lecture, the term is not an isolated label: This concept is needed to reason about hi column density: from spectra to the atomic gas budget. For a Gaussian line profile with peak brightness temperature T_peak and full width at half maximum (FWHM) Delta_v, the velocity integral has the closed form integral(T_B dv) = T_peak * Delta_v * sqrt(pi/(4 ln2)), letting the column-density formula be evaluated directly from two observable quantities (peak and width) rather than requiring numerical integration of every spectrum.

It is tied to the lecture’s observables or evidence in this way: All-sky 21 cm surveys such as the Leiden/Argentine/Bonn (LAB) survey and its modern successor, HI4PI (HI4PI Collaboration 2016), provide column densities toward essentially any Galactic direction, including toward extragalactic background sources used for X-ray absorption studies (e.g., 3C 273, Cygnus X-1). A second connection is important for using the term in physical reasoning: A real sightline’s spectrum is typically the sum of several distinct velocity components (CNM clouds embedded in a broader WNM background), so the total column density is the sum of each component’s individually computed contribution -- exactly the two-component decomposition this course’s representative CNM+WNM spectrum illustrates. The observational interpretation is constrained by The optically thin approximation, while not universally valid (it can underestimate N(HI) in unusually cold, dense, or high-column CNM structures where self-absorption becomes significant), is an excellent approximation for the majority of high-latitude Galactic sightlines used in this way.

Boundary and common confusion. Applying the optically thin approximation uncritically to every CNM component. Narrow, cold, high-column CNM clumps can have tau approaching or exceeding order unity, in which case the true column density (from the full T_B=T_spin(1-e^-tau) solution) is larger than the optically thin formula would suggest; the assumption should always be checked, not assumed, particularly for narrow, high-peak-brightness features. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

brightness-temperature integral

The area under a radio-line brightness-temperature profile as a function of velocity, measured in K km/s. In this lecture, the term is not an isolated label: This concept is needed to reason about hi column density: from spectra to the atomic gas budget. A real sightline’s spectrum is typically the sum of several distinct velocity components (CNM clouds embedded in a broader WNM background), so the total column density is the sum of each component’s individually computed contribution -- exactly the two-component decomposition this course’s representative CNM+WNM spectrum illustrates.

It is tied to the lecture’s observables or evidence in this way: The optically thin approximation, while not universally valid (it can underestimate N(HI) in unusually cold, dense, or high-column CNM structures where self-absorption becomes significant), is an excellent approximation for the majority of high-latitude Galactic sightlines used in this way. A second connection is important for using the term in physical reasoning: Starting from Lecture 03’s optically thin solution T_B(v) ~= T_spin*tau(v), and the definition of 21 cm optical depth in terms of the lower-level column density and the Einstein A coefficient, the total (velocity-integrated) HI column density reduces to a simple integral of the observed brightness temperature over velocity, independent of the (unknown, but nearly canceling) spin temperature: N(HI) = C * integral(T_B(v) dv), where C = 1.8224e18 cm^-2 (K km/s)^-1 is a fixed combination of fundamental constants and the 21 cm transition’s atomic properties (the Einstein A coefficient, statistical weights, and hnu/k). The observational interpretation is constrained by The total number of hydrogen atoms per unit area along a line of sight (the column density N(HI), typically reported in cm^-2) is the single most useful summary statistic of an HI spectrum for nearly every downstream application, from Galactic structure studies to correcting extragalactic X-ray/UV absorption.

Boundary and common confusion. Applying the optically thin approximation uncritically to every CNM component. Narrow, cold, high-column CNM clumps can have tau approaching or exceeding order unity, in which case the true column density (from the full T_B=T_spin(1-e^-tau) solution) is larger than the optically thin formula would suggest; the assumption should always be checked, not assumed, particularly for narrow, high-peak-brightness features. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

velocity-integrated intensity

Spectral-line intensity summed over Doppler velocity, combining emission from all channels assigned to the component of interest. In this lecture, the term is not an isolated label: This concept is needed to reason about hi column density: from spectra to the atomic gas budget. Starting from Lecture 03’s optically thin solution T_B(v) ~= T_spin*tau(v), and the definition of 21 cm optical depth in terms of the lower-level column density and the Einstein A coefficient, the total (velocity-integrated) HI column density reduces to a simple integral of the observed brightness temperature over velocity, independent of the (unknown, but nearly canceling) spin temperature: N(HI) = C * integral(T_B(v) dv), where C = 1.8224e18 cm^-2 (K km/s)^-1 is a fixed combination of fundamental constants and the 21 cm transition’s atomic properties (the Einstein A coefficient, statistical weights, and hnu/k).

It is tied to the lecture’s observables or evidence in this way: The total number of hydrogen atoms per unit area along a line of sight (the column density N(HI), typically reported in cm^-2) is the single most useful summary statistic of an HI spectrum for nearly every downstream application, from Galactic structure studies to correcting extragalactic X-ray/UV absorption. A second connection is important for using the term in physical reasoning: For a Gaussian line profile with peak brightness temperature T_peak and full width at half maximum (FWHM) Delta_v, the velocity integral has the closed form integral(T_B dv) = T_peak * Delta_v * sqrt(pi/(4 ln2)), letting the column-density formula be evaluated directly from two observable quantities (peak and width) rather than requiring numerical integration of every spectrum. The observational interpretation is constrained by All-sky 21 cm surveys such as the Leiden/Argentine/Bonn (LAB) survey and its modern successor, HI4PI (HI4PI Collaboration 2016), provide column densities toward essentially any Galactic direction, including toward extragalactic background sources used for X-ray absorption studies (e.g., 3C 273, Cygnus X-1).

Boundary and common confusion. Applying the optically thin approximation uncritically to every CNM component. Narrow, cold, high-column CNM clumps can have tau approaching or exceeding order unity, in which case the true column density (from the full T_B=T_spin(1-e^-tau) solution) is larger than the optically thin formula would suggest; the assumption should always be checked, not assumed, particularly for narrow, high-peak-brightness features. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

HI4PI survey

A full-sky 21 cm neutral-hydrogen survey combining northern and southern single-dish data into a uniform spectral data set. In this lecture, the term is not an isolated label: This concept is needed to reason about hi column density: from spectra to the atomic gas budget. For a Gaussian line profile with peak brightness temperature T_peak and full width at half maximum (FWHM) Delta_v, the velocity integral has the closed form integral(T_B dv) = T_peak * Delta_v * sqrt(pi/(4 ln2)), letting the column-density formula be evaluated directly from two observable quantities (peak and width) rather than requiring numerical integration of every spectrum.

It is tied to the lecture’s observables or evidence in this way: All-sky 21 cm surveys such as the Leiden/Argentine/Bonn (LAB) survey and its modern successor, HI4PI (HI4PI Collaboration 2016), provide column densities toward essentially any Galactic direction, including toward extragalactic background sources used for X-ray absorption studies (e.g., 3C 273, Cygnus X-1). A second connection is important for using the term in physical reasoning: A real sightline’s spectrum is typically the sum of several distinct velocity components (CNM clouds embedded in a broader WNM background), so the total column density is the sum of each component’s individually computed contribution -- exactly the two-component decomposition this course’s representative CNM+WNM spectrum illustrates. The observational interpretation is constrained by The optically thin approximation, while not universally valid (it can underestimate N(HI) in unusually cold, dense, or high-column CNM structures where self-absorption becomes significant), is an excellent approximation for the majority of high-latitude Galactic sightlines used in this way.

Boundary and common confusion. Applying the optically thin approximation uncritically to every CNM component. Narrow, cold, high-column CNM clumps can have tau approaching or exceeding order unity, in which case the true column density (from the full T_B=T_spin(1-e^-tau) solution) is larger than the optically thin formula would suggest; the assumption should always be checked, not assumed, particularly for narrow, high-peak-brightness features. The term should therefore be applied only after identifying the measured quantity, the model assumption that links it to the desired physical quantity, and the regime in which that link remains valid.

Evidence

Derivation and Model

Working Equation

\[ N(\mathrm{HI}) = 1.8224\times10^{18}\,\mathrm{cm^{-2}}\int T_B(v)\,dv, \qquad \int T_B\,dv = T_{\rm peak}\,\Delta v\sqrt{\pi/(4\ln 2)}\ \text{(Gaussian)} \]

Worked Example

  1. This course’s representative CNM component (peak T_B=25 K, FWHM=3.0 km/s, illustrative synthetic profile shape but real, standard CNM linewidth/T_spin values) gives N(HI)_CNM=1.455e+20 cm^-2 by direct Gaussian integration.
  2. The representative WNM component (peak T_B=4 K, FWHM=25 km/s) gives N(HI)_WNM=1.940e+20 cm^-2, and the total two-component column density is N(HI)_total=3.395e+20 cm^-2 -- of the same general order of magnitude as the real, published 3C 273 sightline value of 1.79e+20 cm^-2, though not an attempt to reproduce that specific sightline’s exact spectrum (this course’s synthetic profile is explicitly a representative illustration, not a claimed measurement of 3C 273 itself; see reference-log.md).
  3. By contrast, the real published Cygnus X-1 sightline column density (N_H=6.0e+21 cm^-2, from X-ray spectral absorption studies) is roughly 34 times larger than toward 3C 273 -- a direct, real illustration of how strongly Galactic HI column density depends on sightline (Cygnus X-1 lies much closer to the Galactic plane, b=3.07 degrees, than 3C 273, b=64.36 degrees).

Common Misconception

Applying the optically thin approximation uncritically to every CNM component. Narrow, cold, high-column CNM clumps can have tau approaching or exceeding order unity, in which case the true column density (from the full T_B=T_spin(1-e^-tau) solution) is larger than the optically thin formula would suggest; the assumption should always be checked, not assumed, particularly for narrow, high-peak-brightness features.

Active-Learning Guidance

Using the Gaussian-integral formula, explain why a broad, low-peak WNM component and a narrow, high-peak CNM component can contribute comparable total column densities even though their peak brightness temperatures differ by a large factor -- and identify which observable (peak temperature or linewidth) matters more for the total column density of a broad component.

Lab Connection

Lab 02 performs this lecture’s full worked HI column-density calculation for a two-component spectrum and compares the result to real, published all-sky-survey values toward two different real extragalactic sightlines.

Synthesis Questions

References