Recombination and photon decoupling
Distinguish equilibrium ionization, radiative recombination kinetics, and the probability of last scattering.
On this page
The question for today
What you’ll understand
- Apply the hydrogen-only Saha equation with a defined ionization fraction.
- Explain the radiative bottlenecks that invalidate instantaneous equilibrium.
- Separate an instantaneous scattering rate from integrated optical depth.
Hydrogen recombination reduces the free-electron density, enabling cosmic photons to travel much farther between scatterings. Atomic neutralization and photon decoupling are related but distinct processes. One concerns populations of atomic and ionized states; the other concerns the probability that a photon scatters along its remaining path. A thermal equilibrium formula provides an instructive baseline, but accurate recombination requires radiative kinetics, and accurate last-scattering predictions require an optical-depth history. We first use a hydrogen-only, electrically neutral gas, then identify the missing physics needed before comparing with a precision cosmic microwave background calculation.
Saha predicts an equilibrium population, not a reaction speed
Let denote the total number density of hydrogen nuclei, whether bound or ionized. Define the free-electron fraction by ; electrical neutrality gives the same free-proton fraction. In chemical equilibrium, the reaction linking an electron, proton, hydrogen atom, and photon fixes the chemical-potential relation. Nondegenerate nonrelativistic matter distributions and a thermal photon bath then give the Saha equation. The ground-state spin multiplicities are included in the conventional prefactor displayed here.
The hydrogen binding energy alone does not set the recombination temperature. Because photons outnumber hydrogen nuclei, the exponential photon tail can keep the gas ionized at temperatures far below the binding energy. The explicit inverse density in the Saha expression matters: at fixed temperature, a denser hydrogen gas favors neutralization. In a real expanding history, both temperature and density change, so one must specify their relation before solving for an epoch.
Atomic radiative bottlenecks leave a nonequilibrium history
Direct capture into the ground state emits a photon energetic enough to ionize another atom, so it need not yield a lasting increase in neutral hydrogen. Population must instead pass through excited states and lose radiation without immediately reversing the net recombination. Two important escape routes are the metastable two-photon transition from the second s state and redshifting of resonant Lyman-alpha photons out of the line. These processes introduce timescales absent from an equilibrium population formula.
As expansion outpaces effective atomic relaxation, the electron fraction departs from the Saha value. A simplified effective-level treatment captures the bottleneck qualitatively; precision work follows more levels and radiative transfer. Helium recombination and later astrophysical reionization also alter the electron history. Thus using Saha alone all the way through last scattering cannot justify a precise visibility peak or residual ionization fraction, even if it gives a reasonable initial scale for when atoms begin forming.
Transparency depends on accumulated scattering probability
At low photon energy, free electrons scatter photons with the Thomson cross section. The rate per photon is free-electron density times cross section times light speed. Compare it with expansion for a local diagnostic, but integrate the rate along a trajectory to obtain optical depth. The probability of no further scattering is an exponential of minus that integrated depth. A particular ionized fraction therefore cannot universally define decoupling independently of density and expansion.
For optical depth defined from an emission epoch to the present, the last-scattering visibility combines the chance to scatter near that epoch with the chance of no scattering afterward. In conformal time it is , since the remaining optical depth decreases toward today. Recombination, the peak of this visibility, and baryon drag decoupling answer different questions. A later reionized population can add optical depth without undoing the earlier atomic history.
Work it through
Half neutral need not mean transparent
At one stipulated instant, a hydrogen-only equilibrium toy model has Saha right-hand side . Take , , and . Find the physical electron fraction and scattering-to-expansion ratio. These rounded inputs do not specify a fitted recombination epoch.
Solve , giving . The roots are and ; only the first lies in the allowed interval from zero to one. The toy gas is half ionized.
Multiply the total hydrogen-nucleus density by the free fraction to obtain . Bound electrons do not enter the free-electron Thomson rate used here.
Compute . Division by the stipulated expansion rate gives . Many scatterings fit into one expansion time despite the substantial neutral fraction.
The equilibrium fraction is , but the instantaneous scattering-to-expansion ratio is . Determining a last-scattering probability still requires the future density and ionization history, rather than this single instant.
The living notebook
Explore the model
Interactive lab · heat
An exponential attenuation analogy
Reset and set both controls to mode , diffusivity , then Step once: the amplitude is . Reset and set both controls to , , then Step once: it is . This is decay of a periodic diffusion mode; its exponential illustrates accumulated attenuation only. It does not model photon paths, atomic populations, Thomson scattering, or cosmological optical depth.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forAn exponential depends on the accumulated decay parameter; a real transparency calculation must supply the correct physical integral.
Pause & reason
Check your understanding
Your turn
Try it yourself
In a separate constant-rate segment along a photon path, let for . Find the optical depth and probability of zero scatterings in that segment. Repeat for twice the duration, with the same rate.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Integrate the rate over time before exponentiating; this segment model is independent of Saha equilibrium.
Compare with the solution
The first optical depth is , so the no-scattering probability is . Doubling the duration gives and probability . The rate is unchanged while the survival probability changes substantially. A real cosmological calculation integrates the time-dependent free-electron density and does not assume a constant rate over the whole recombination history.
Keep the idea
Take it with you
- Saha supplies equilibrium populations under specified thermodynamic assumptions.
- Radiative bottlenecks make recombination a kinetic problem.
- Photon decoupling is governed by an optical-depth history, not a universal ionized fraction.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.