Past exam of the mathematics course of the University of Cambridge / 2013 / iii / Paper 48 / 2

For the reaction $p+e\leftrightarrow H+\gamma$, chemical equilibrium requires ${\mu }_p+{\mu }_e={\mu }_H$ because thermal photons have zero chemical potential. Insert the nonrelativistic Maxwell-Boltzmann distribution number densities and use $m_H=m_p+m_e-\mathcal{B}$:For ground-state hydrogen, including its spin states, $g_e=g_p=2$, $g_H=4$; hence the degeneracy factor is one. Since $m_p/m_H\simeq 1$, the right side is $(m_{e}T/2\pi {)}^{3/2}{e}^{-\mathcal{B}/T}$.

Charge neutrality gives $n_p=n_e$, while baryon conservation gives $n_b=n_p+n_H$. Consequently $n_e=n_p=X_e{n}_b$ and $n_H=(1-X_e)n_b$. Inverting the preceding ratio and inserting the photon number density with $n_b=\eta n_{\gamma }$ gives the hydrogen-only Saha equationWriting $y=\mathcal{B}/T$, its coefficient is $(2\zeta (3)/{\pi }^2)\eta (2\pi \mathcal{B}/m_e{)}^{3/2}\simeq 3.2\times 10^{-16}$, consistent with the rounded $3\times 10^{-16}$. The negligible mass-ratio correction and ground-state approximation are the assumptions behind this form.

At $X_e=1/2$, the left side of the Saha equation is two. With its rounded coefficient, the recombination temperature therefore satisfiesNeglecting the logarithmic prefactor gives $y\sim 36$ andThis is the requested rough scale. Retaining the prefactor and the unrounded logarithm gives $y\simeq 42.0$ and $T_{\mathrm{rec}}\simeq 0.32$ eV; the displayed $0.4$ eV is not a high-accuracy numerical root of the Saha equation.

The small baryon-to-photon ratio means that radiation and the entropy of the free charged particles strongly favor ionization. A low mean photon energy does not eliminate the energetic tail, and ionized particles have a large phase-space advantage. The exponential binding factor must become very large before neutral atoms dominate this dilute gas. Accordingly $\mathcal{B}/T$ is of order forty rather than order one, so recombination occurs well below the hydrogen binding energy.

Photon decoupling is a dynamical loss of frequent scattering, distinct from the chemical conversion of charged particles into neutral atoms. Before recombination, Thomson scattering on free Electrons keeps photons coupled to the baryon plasma. The physical scattering rate isin units with $c=1$. Define the approximate decoupling temperature byEquivalently the photon mean scattering time becomes comparable to an expansion time. A refined last-scattering definition uses the optical depth and the maximum of the visibility function; it need not coincide exactly with the local rate criterion. As $X_e$ falls, scattering becomes inefficient, ordinarily after the half-ionized recombination stage, so $T_{\mathrm{dec}}<T_{\mathrm{rec}}$.

Residual electron freeze-out concerns the remaining ionized fraction. The Saha equation presumes sufficiently rapid forward and reverse reactions to maintain chemical equilibrium. Real recombination is slowed by radiative bottlenecks and expansion; the free-electron fraction can therefore exceed its equilibrium value. Eventually the effective recombination rate per free Electron, roughly ${\alpha }_{\mathrm{eff}}(T)n_e$, falls below $H$. The remaining Electrons cannot all find and recombine with protons within an expansion time. A small nonzero residual fraction survives while the Saha prediction decreases exponentially toward zero.

The sketch shows this lag and residual floor. Its kinetic curve and the position of the decoupling marker are schematic: the data supplied determine the equilibrium curve, but not a quantitative recombination history or an exact $T_{\mathrm{dec}}$. Determining those needs the expansion history and atomic transition rates. The horizontal axis decreases toward the right to follow cosmological cooling.