Past exam of the mathematics course of the University of Cambridge / 2008 / iii / Paper 62 / 3

Work in natural units $\hslash =c=k_B=1$, at temperatures far below the $W$-boson mass but of order or above the electron mass. The Fermi interaction then supplies a local four-fermion amplitude. The Fermi constant has dimension energy${}^{-2}$. In scattering of a thermal lepton from a nonrelativistic nucleon, the nucleon current matrix element has size $m_N$, while the lepton current has size $T$. Thus $\mathcal{M}\sim G_F{m}_{N}T$. The two-body relativistic scattering cross-section contains the kinematic factor $s^{-1}\sim m_N^{-2}$; when the final and initial lepton energies are comparable, their momentum ratio is of order one. Consequently $\sigma \sim |\mathcal{M}{|}^2/m_N^2\sim G_F^2{T}^2$. This estimates temperatures above the electron mass and mass-splitting thresholds; those scales only affect the dimensionless rate function near freeze-out.

The light leptons in the thermal bath have number density $n\sim T^3$ and speed of order unity. Their conversion rate per nucleon is thereforeFor example, $n+{\nu }_e\leftrightarrow p+e^{-}$ and $n+e^{+}\leftrightarrow p+{\bar{\nu }}_e$ change the neutron-to-proton ratio. The bath density here is the lepton density, not the much smaller baryon density. Dimensionless weak couplings, phase-space integrals and the several reaction channels set a numerical coefficient. Electron masses and the neutron-proton splitting also become important for a precision rate near cosmological weak freeze-out.

During radiation domination, the effective number of relativistic energy degrees of freedom, denoted $g_{\ast }$, gives the energy densityThe paper's ${\mathcal{N}}_{\mathrm{eff}}$ is this $g_{\ast }$; it is not the conventionally normalized effective number of neutrino species. With the specified reduced Planck mass, the Friedmann equation reads $H^2={\rho }_R/(3M_{\mathrm{Pl}}^2)$, whenceWrite the weak interaction conversion rate as $\Gamma =A{G}_F^2{T}^5$. The ratio $\Gamma /H$ scales as $T^3$, so it falls during cooling: conversion eventually becomes too slow to maintain neutron-proton chemical equilibrium. Equating the conversion and expansion rates givesOn suppressing the numerical coefficient, the cosmological weak freeze-out temperature is consequentlyThe one-sixth power reflects a balance between a conversion rate proportional to $T^5$ and a Hubble parameter proportional to $\sqrt{g_{\ast }}{T}^2$.

Before electron-positron annihilation in cosmology, photons contribute two helicity states; electrons and positrons contribute two spin states each. Each of the three light neutrino flavours contributes one active neutrino helicity and one antineutrino helicity. At a common temperature the effective number of relativistic energy degrees of freedom is thereforeUsing $G_F=1.17\times 10^{-5}{\mathrm{GeV}}^{-2}$ and $M_{\mathrm{Pl}}\simeq 2.435\times 10^{18}\mathrm{GeV}$ in the scaling formula alone gives $T_F\sim 2\mathrm{MeV}$. Including the explicit expansion coefficient from (ii), but setting the unknown rate coefficient $A$ to one, gives $T_F\simeq 1.5\mathrm{MeV}$. Thus the warranted conclusion is a freeze-out temperature of order $1\mathrm{MeV}$, above the $0.511\mathrm{MeV}$ electron mass. The omitted weak-rate coefficient prevents a more precise prediction from this estimate.

Neutrino decoupling occurs before most electron-positron pairs annihilate. Afterwards neutrinos redshift freely, $T_{\nu }a=\mathrm{constant}$, while the photons receive the pair entropy. For the still-coupled photon-pair gas, cosmological entropy conservation gives $g_{\ast s}{T}_{\gamma }^3{a}^3=\mathrm{constant}$. Its entropy degrees of freedom drop from $2+(7/8)4=11/2$ to $2$. Comparing this with the freely cooling neutrinos givesAfter annihilation, photon and neutrino temperatures differ. To write the total energy density as $({\pi }^2/30)g_{\ast }{T}_{\gamma }^4$, the neutrino contribution must include that temperature ratio:This is the paper's post-annihilation ${\mathcal{N}}_{\mathrm{eff}}$. Merely deleting the electron term while keeping all species at the same temperature would be incorrect. For comparison the entropy count, normalized to $T_{\gamma }^3$, is $g_{\ast s}=2+(7/8)6(4/11)=43/11\simeq 3.91$. These results assume instantaneous neutrino decoupling and relativistic neutrinos; the two counts represent different thermal moments.

The nucleons are nonrelativistic, so the Maxwell-Boltzmann distribution givesBoth spin degeneracies are two. Taking the ratio and using the specified negligible chemical potentials yields the neutron-to-proton ratioThe mass prefactor differs from unity only at the per-mille level. The dependence on $T_F$ is much stronger: $d\ln r_F/d\ln T_F=Q/T_F>0$. Thus earlier, hotter weak freeze-out leaves more neutrons relative to protons. For illustration, $T_F=1\mathrm{MeV}$ gives $r_F\simeq 0.28$; $T_F=0.8\mathrm{MeV}$ gives $r_F\simeq 0.20$. Free-neutron decay after freeze-out lowers the ratio further before most Big Bang nucleosynthesis takes place.

The small baryon-to-photon ratio suppresses nuclear equilibrium abundances through the factor ${\eta }^{A-1}$. Cooling eventually makes the binding-energy exponential favourable, but the deuterium bottleneck delays substantial nuclear burning until well after cosmological weak freeze-out. In ordinary Big Bang nucleosynthesis, nearly every neutron surviving that delay ends up in helium-4. This makes the final primordial helium mass fraction chiefly a neutron budget rather than an unrestricted equilibrium abundance.

Let $r_F$ be the neutron-to-proton ratio at freeze-out, let $\Delta t$ be the delay until effective nuclear burning, and let ${\tau }_n$ be the neutron mean lifetime. The initial fraction of baryons which are neutrons is $x_{n,F}=r_F/(1+r_F)$. Because neutron decay converts a neutron into a proton, total baryon number remains constant, while the neutron fraction becomes $x_{n,B}=x_{n,F}{e}^{-\Delta t/{\tau }_n}$. Each helium nucleus uses two neutrons and four baryons. Hence neutron-limited helium synthesis givesEquivalently, $Y_4=2r_B/(1+r_B)$ with $r_B=x_{n,B}/(1-x_{n,B})$. This accounts for the extra protons produced by neutron decay; multiplying $r_F$ alone by the survival factor would miss their contribution. The following responses refer to ordinary neutron-limited burning, without exotic particle injection or a change to the nuclear network.

Extra relativistic neutrino species increase the effective number of relativistic energy degrees of freedom. At fixed temperature the Hubble parameter rises as $g_{\ast }^{1/2}$, while the standard conversion rate is unchanged. Consequently cosmological weak freeze-out moves to higher temperature, $T_F\propto g_{\ast }^{1/6}$, and the neutron-to-proton ratio $r_F=e^{-Q/T_F}$ increases. The initial neutron fraction $r_F/(1+r_F)$ therefore increases as well.

The quicker expansion also reduces the cosmic time required to cool to the deuterium bottleneck temperature: during radiation domination, $t\propto g_{\ast }^{-1/2}{T}^{-2}$ while the thermal count is roughly constant. Fewer neutrons decay before burning. Both factors in the neutron budget increase, giving a larger final helium-4 mass fraction. This is the primordial helium response to extra relativistic species, under the assumption that the additional species contribute mainly to expansion rather than independently changing the conversion rates.

A longer neutron half-life means a larger mean lifetime ${\tau }_n={\tau }_{1/2}/\ln 2$. Holding the initial neutron-to-proton ratio and the nuclear-burning delay fixed, the derivative of the neutron survival factor is positive:Thus fewer neutrons disappear during the deuterium bottleneck, and the final helium-4 mass fraction increases.

There can be a second effect if the increased lifetime represents a uniform reduction of the relevant weak interaction strength. With fixed masses and other couplings, neutron decay and the conversion reactions both scale as $G_F^2$. Thus ${\tau }_n\propto G_F^{-2}$, whereas $T_F\propto G_F^{-2/3}\propto {\tau }_n^{1/3}$. Weaker conversion then freezes out earlier and increases $r_F$, reinforcing the survival effect. The survival argument alone already establishes the requested neutron half-life effect on primordial helium; linking lifetime to all conversion rates additionally requires this microscopic assumption.

A larger baryon-to-photon ratio favours bound nuclei at a fixed temperature through the equilibrium factor ${\eta }^{A-1}$. In particular it helps deuterium survive photodissociation sooner: there are fewer energetic photons per baryon, and nuclear reaction rates per baryon also rise. A schematic deuterium bottleneck condition is $\eta e^{B_D/T_B}\sim \mathrm{constant}$, suppressing slowly varying powers of temperature. It givesThus higher $\eta$ permits burning at a higher temperature. During radiation domination that means earlier cosmic time, $t_B\propto T_B^{-2}$, leaving less time for free-neutron decay after freeze-out. More neutrons enter helium, so the final helium-4 mass fraction increases, ordinarily rather weakly. This is the baryon-density effect on primordial helium.

The equilibrium helium factor ${\eta }^3$ indicates the direction of the change but is not a prediction that the final helium abundance grows without bound as ${\eta }^3$. Once almost all surviving neutrons are captured, the final primordial helium mass fraction is limited by their supply, as the shared neutron budget shows.