ASTR 360 · Lecture 04

HI Column Density: From Spectra to the Atomic Gas Budget

Turning a radio spectrum into a real number of hydrogen atoms per square centimeter

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.

Reading anchor: OpenStax Astronomy 2e (local extracted reference copy: references/openstax-astronomy-2e-extracted.txt) covers the interstellar medium in Chapter 20 ("Between the Stars: Gas and Dust in Space"): 20.1 The Interstellar Medium, 20.2 Interstellar Gas (including the 21 cm line and H II regions), 20.3 Cosmic Dust (extinction and reddening), 20.4 Cosmic Rays, 20.5 The Life Cycle of Cosmic Material, and 20.6 Interstellar Matter around the Sun. Chapter 20’s section numbering was cross-checked directly against its own embedded figure numbering this session (Figures 20.1 through 20.20 all appear under their matching section headers with no offset), so this citation is verified against the local extracted text rather than extrapolated by analogy. OpenStax’s treatment throughout Chapter 20 is descriptive and qualitative; nearly every quantitative derivation in this course (the 21 cm hyperfine radiative-transfer solution and HI column-density formula, the Jeans mass/length instability criterion, the free-fall collapse time, the Stromgren-sphere ionization-balance derivation, the Rankine-Hugoniot shock-jump conditions, and the cosmic-ray power-law spectrum) goes well beyond OpenStax Astronomy 2e’s introductory level; each lecture states explicitly where this course extends beyond the assigned reading.

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.

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.

First question: what here is directly observed or a published/live-verified measurement, and what follows only once a physical law is derived and applied?

Terms for This Lecture

column density

Definition: The number of particles per unit projected area integrated along the line of sight, commonly expressed in cm^-2.

Why it matters here: 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).

Relationship: 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.

Limitation or 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.

Evidence We Need to Explain

Derivation and Model

Quantitative Tool

\[ 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).

Visual Reasoning

Trace the reasoning chain from raw observation or live-verified/published data to derived physical law to inferred quantity in this lecture’s figure.

  • Which stage is directly observed or a real published measurement?
  • Which stage is the physical law or derivation step?
  • What would change if an assumption in the derivation failed?
Lecture 04: HI Column Density: From Spectra to the Atomic Gas BudgetLSR velocity (km/s)Brightness temperature T_B (K)Representative CNM+WNM spectrum (illustrative shape): integrating gives N(HI)=3.39e+20 cm^-2 (CNM 1.45e+20 + WNM 1.94e+20).
HI Column Density: From Spectra to the Atomic Gas Budget: from observation to inferred quantity.

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 Segment

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

The 21 cm column-density formula converts an entirely passive radio observation into a direct census of atomic hydrogen along any line of sight in the Galaxy -- closing out this course’s atomic-ISM unit and setting up Unit 3’s shift to the interstellar medium’s solid component: dust.

References