Past exam of the mathematics course of the University of Cambridge / 2010 / iii / Paper 53 / 2 / i / Solution

At temperatures above a few MeV, weak interactions such as $n+{\nu }_e\leftrightarrow p+e^{-}$ and $n+e^{+}\leftrightarrow p+{\bar{\nu }}_e$ maintain neutron-proton chemical equilibrium. Negligible lepton chemical potentials give the neutron-to-proton ratio $r_{\mathrm{eq}}=n_n/n_p=\exp [-\Delta /(k_{B}T)]$, where $\Delta =(m_n-m_p)c^2\simeq 1.293\mathrm{MeV}$. At high temperature this ratio is nearly one.

During radiation domination, the interaction time grows as $T^{-5}$, faster than the expansion time $T^{-2}$. Their crossing therefore produces cosmological weak freeze-out. If $\Theta =T/(10^{10}\mathrm{K})$, equating the displayed order-of-magnitude times gives ${\Theta }_f^{-5}\sim g_{\ast }^{-1/2}{\Theta }_f^{-2}$, soUsing $H^{-1}=2t$ instead of $t$ changes only the order-one coefficient. Before electron-positron annihilation, the usual effective relativistic count is $g_{\ast }=2+(7/8)(4+6)=10.75$. The factor $7/8$ weights fermionic relative to bosonic energy density.

After freeze-out, free neutrons decay while protons increase correspondingly. Define the neutron baryon fraction $X_n=n_n/(n_n+n_p)$. In the instantaneous-freeze-out estimate,where ${\tau }_n=t_{1/2}/\ln 2$ is the neutron mean lifetime. Hence the ratio at the onset of nuclear synthesis isThis expression accounts for the protons produced by neutron decay.

Nuclear burning initially waits at the deuterium bottleneck: although deuterium is bound by about $2.2\mathrm{MeV}$, the many photons per baryon can still photodissociate it. Roughly, ${\eta }_b^{-1}{e}^{-B_D/(k_B{T}_N)}\sim 1$ sets the bottleneck temperature, with refinements from reaction rates and phase space. Around $k_B{T}_N\sim 0.1\mathrm{MeV}$, deuterium survives, reactions through deuterium, tritium and helium-3 rapidly build tightly bound helium-4, and nearly all surviving neutrons are captured.

Thus $n_{\mathrm{He}}\simeq n_n/2$ and the primordial helium mass fraction isFor the standard estimate $r_N\simeq 1/7$, neutron and proton constituent fractions are $1/8$ and $7/8$, respectively, and $Y_{\mathrm{He}}\simeq 1/4$. Free neutrons are negligible after capture. The remaining hydrogen mass fraction is about $3/4$, and the helium-to-hydrogen number ratio is about $1/12$.

The printed numerical time coefficients are schematic, not a precision prediction of this quarter fraction: their literal crossing with $g_{\ast }=10.75$ gives $k_B{T}_f\simeq 1.28\mathrm{MeV}$. Keeping physical weak-rate coefficients and the gradual freeze-out instead gives an effective freeze-out around $0.7$–$0.8\mathrm{MeV}$, followed by decay over a few minutes. The abundance formulas above explain the standard quarter fraction and its parameter dependence without treating the rounded time formulas as exact nuclear-rate data.