Lecture 04: Why Spectral Lines Look the Way They Do: The Boltzmann and Saha Equations

Context and Why This Matters

ASTR 230 explained the non-monotonic strength of Balmer lines across spectral type qualitatively, as a competition between excitation and ionization. This lecture makes that competition quantitative: the Boltzmann and Saha equations are the actual statistical-mechanics tools astronomers use to convert a measured line strength into a precise temperature, and they underlie every stellar-atmosphere model used in modern astrophysics.

Learning Goals

1. Apply the Boltzmann equation to compute the relative population of two atomic energy levels at a given temperature.
2. Explain qualitatively how the Saha equation extends this to ionization states, and why hydrogen Balmer lines peak near spectral type A.
3. Connect these microphysical population equations to the macroscopic spectral classification introduced in ASTR 230.

Opening Phenomenon

Two stars can have almost the same total hydrogen abundance (essentially all stars do) yet dramatically different Balmer line strengths, because the fraction of hydrogen atoms with an electron in the first excited state (n=2, required to produce a Balmer absorption line) depends exponentially on temperature through the Boltzmann factor.

Terms Developed in Context

These terms form the reasoning chain from observation to model, derivation, and interpretation. Each entry expands the corresponding term-explorer panel and should be read alongside the worked example.

Boltzmann equation

The Boltzmann equation gives the ratio of the number of atoms in excited state 2 to state 1 as n2/n1 = (g2/g1) e^{-(E2-E1)/kT}, where g is the statistical weight (degeneracy) of each level; higher T always increases the population of the higher level relative to the lower one.

Use and relationship: The Saha equation gives the ratio of ionized to neutral atom number densities as a function of T and electron density, with an exponential dependence on the ionization energy divided by kT, analogous in form to the Boltzmann equation but for a continuum (free-electron) final state rather than a discrete bound level. The lecture's quantitative connection is \( \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \), which shows how this term participates in the model rather than functioning as an isolated label.

Boundary or confusion: Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

excitation temperature

At still higher temperatures (O and B stars), the Saha equation predicts the neutral-hydrogen fraction itself collapses toward zero as hydrogen ionizes, so even though the Boltzmann n2/n1 ratio keeps rising, there are too few neutral atoms left to produce strong Balmer absorption — reproducing the observed turnover.

Use and relationship: This lecture makes that competition quantitative: the Boltzmann and Saha equations are the actual statistical-mechanics tools astronomers use to convert a measured line strength into a precise temperature, and they underlie every stellar-atmosphere model used in modern astrophysics. The lecture's quantitative connection is \( \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \), which shows how this term participates in the model rather than functioning as an isolated label.

Boundary or confusion: Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

statistical weight

The Boltzmann equation gives the ratio of the number of atoms in excited state 2 to state 1 as n2/n1 = (g2/g1) e^{-(E2-E1)/kT}, where g is the statistical weight (degeneracy) of each level; higher T always increases the population of the higher level relative to the lower one.

Use and relationship: This lecture makes that competition quantitative: the Boltzmann and Saha equations are the actual statistical-mechanics tools astronomers use to convert a measured line strength into a precise temperature, and they underlie every stellar-atmosphere model used in modern astrophysics. The lecture's quantitative connection is \( \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \), which shows how this term participates in the model rather than functioning as an isolated label.

Boundary or confusion: Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

Saha equation

The Saha equation gives the ratio of ionized to neutral atom number densities as a function of T and electron density, with an exponential dependence on the ionization energy divided by kT, analogous in form to the Boltzmann equation but for a continuum (free-electron) final state rather than a discrete bound level.

Use and relationship: At still higher temperatures (O and B stars), the Saha equation predicts the neutral-hydrogen fraction itself collapses toward zero as hydrogen ionizes, so even though the Boltzmann n2/n1 ratio keeps rising, there are too few neutral atoms left to produce strong Balmer absorption — reproducing the observed turnover. The lecture's quantitative connection is \( \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \), which shows how this term participates in the model rather than functioning as an isolated label.

Boundary or confusion: Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

ionization fraction

The Saha equation gives the ratio of ionized to neutral atom number densities as a function of T and electron density, with an exponential dependence on the ionization energy divided by kT, analogous in form to the Boltzmann equation but for a continuum (free-electron) final state rather than a discrete bound level.

Use and relationship: Balmer-line strength is proportional to the number of neutral hydrogen atoms in the n=2 state, which is the product of (fraction neutral, from Saha) and (fraction of neutral atoms in n=2, from Boltzmann); this product is non-monotonic in T because the two factors move in opposite directions as T rises. The lecture's quantitative connection is \( \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \), which shows how this term participates in the model rather than functioning as an isolated label.

Boundary or confusion: Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

partition function

The Saha equation gives the ratio of ionized to neutral atom number densities as a function of T and electron density, with an exponential dependence on the ionization energy divided by kT, analogous in form to the Boltzmann equation but for a continuum (free-electron) final state rather than a discrete bound level.

Use and relationship: The Boltzmann equation gives the ratio of the number of atoms in excited state 2 to state 1 as n2/n1 = (g2/g1) e^{-(E2-E1)/kT}, where g is the statistical weight (degeneracy) of each level; higher T always increases the population of the higher level relative to the lower one. The lecture's quantitative connection is \( \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \), which shows how this term participates in the model rather than functioning as an isolated label.

Boundary or confusion: Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

Evidence

• The Boltzmann equation predicts that the ratio of atoms in an excited state to the ground state rises steeply (exponentially) with temperature, exactly matching the observed steep rise in Balmer-line strength from cool M stars toward intermediate A stars.
• The Saha equation predicts that above roughly 10,000 K, hydrogen becomes predominantly ionized, removing the bound electron needed for any Balmer transition — exactly matching the observed decline in Balmer-line strength from A stars toward the hottest O and B stars.
• Combining both equations correctly predicts the observed peak of Balmer-line strength near spectral type A0 (T_eff ≈ 9,000-10,000 K), a result that neither equation alone can explain.

Derivation and Model

• The Boltzmann equation gives the ratio of the number of atoms in excited state 2 to state 1 as n2/n1 = (g2/g1) e^{-(E2-E1)/kT}, where g is the statistical weight (degeneracy) of each level; higher T always increases the population of the higher level relative to the lower one.
• The Saha equation gives the ratio of ionized to neutral atom number densities as a function of T and electron density, with an exponential dependence on the ionization energy divided by kT, analogous in form to the Boltzmann equation but for a continuum (free-electron) final state rather than a discrete bound level.
• Balmer-line strength is proportional to the number of neutral hydrogen atoms in the n=2 state, which is the product of (fraction neutral, from Saha) and (fraction of neutral atoms in n=2, from Boltzmann); this product is non-monotonic in T because the two factors move in opposite directions as T rises.

Working Equation

\[ \dfrac{n_2}{n_1} = \dfrac{g_2}{g_1}\,e^{-(E_2-E_1)/kT}, \qquad \dfrac{n_{i+1}}{n_i} \propto T^{3/2} e^{-\chi_i/kT}\ \ (\text{Saha, schematic form}) \]

Worked Example

1. Boltzmann ratio for hydrogen n=2 relative to n=1 (E2 - E1 = 10.2 eV = 1.634e-18 J, g2/g1 = 4): at T = 5,750 K (near the Sun and Alpha Centauri A/B), kT = 7.941e-20 J, giving n2/n1 = 4 × exp(-20.6) = 4.628e-09 — an extremely small fraction of atoms are excited, so Balmer lines are intrinsically weak in G/K stars like the Sun and Alpha Centauri A and B.
2. At T = 9,900 K (near Vega/Sirius A, spectral type A): kT = 1.367e-19 J, giving n2/n1 = 4 × exp(-11.95) = 2.580e-05 — roughly 5574 times larger than at 5,750 K, consistent with the much stronger observed Balmer lines in A-type stars.
3. At still higher temperatures (O and B stars), the Saha equation predicts the neutral-hydrogen fraction itself collapses toward zero as hydrogen ionizes, so even though the Boltzmann n2/n1 ratio keeps rising, there are too few neutral atoms left to produce strong Balmer absorption — reproducing the observed turnover.

Common Misconception

Using the Boltzmann equation alone to explain the full OBAFGKM Balmer-strength sequence. The Boltzmann equation alone predicts monotonically increasing line strength with temperature; only combining it with the Saha equation’s ionization turnover reproduces the observed peak at spectral type A.

Active-Learning Guidance

Using only the sign of the exponent in each equation, explain in words (no calculation needed) why raising temperature always increases excitation (Boltzmann) but can either increase or asymptotically saturate ionization (Saha) depending on how close the gas already is to fully ionized.

Lab Connection

Lab 02 computes the Boltzmann excitation ratio for hydrogen across a range of stellar temperatures and reproduces the qualitative OBAFGKM Balmer-strength curve from ASTR 230 quantitatively.

Synthesis Questions

• What was measured directly in this lecture’s worked example, and what was derived from a physical law?
• Which assumption in the derivation would most change the interpretation if it were wrong?
• How does this lecture’s technique connect to the lab and problem set that follow it?

References

• OpenStax Astronomy 2e, Chapter 5.5 (formation of spectral lines) and Chapter 17.3 (spectral types); this course adds the Boltzmann/Saha formalism omitted at the introductory level.
• Every numeric result above is computed programmatically in materials/ASTR310/src/generate_astr310_content.py, not hand-typed.