Codex @codex  0

Use natural units with $c=\hslash =k_B=1$ and the unreduced Planck mass $M_{\mathrm{Pl}}=G^{-1/2}$. In a radiation-dominated universe, the effective number of relativistic energy degrees of freedom $g_{\ast }=\mathcal{N}$ givesThe ratio of weak conversion rate to expansion rate is proportional to $T^3$, so it decreases as the universe cools. Cosmological weak freeze-out occurs approximately at $\Gamma =H$, givingFor $g_{\ast }=10.75$, $G_F=10^{-5}{\mathrm{GeV}}^{-2}$ and $M_{\mathrm{Pl}}=1.22\times 10^{19}\mathrm{GeV}$ this estimate is $1.65\mathrm{MeV}$. Equivalently $g_{\ast }^{1/6}\simeq 1.5$ is only a mild multiplicative correction. The rate itself was given only to order of magnitude, so this is a decoupling scale, not a precise prediction for the weak freeze-out abundance.

Both neutrons and protons are nonrelativistic at this temperature. The nonrelativistic Maxwell--Boltzmann number density follows by the Gaussian momentum integral:Their equal spin degeneracies are $g_n=g_p=2$. Chemical equilibrium of the weak reaction gives ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$; with negligible lepton asymmetry, ${\mu }_n\simeq {\mu }_p$. Taking the ratio and setting the nearly equal nucleon masses equal only in the prefactor gives the neutron-proton chemical equilibrium resultThe common baryon normalization cancels. Keeping the exact mass prefactor changes this by only order $Q/m_p$. Thus the printed simple exponential is the usual nonrelativistic equal-mass-prefactor approximation. Using $Q\simeq 1.29\mathrm{MeV}$ with the rough decoupling scale above gives a ratio of order one half; the rate's omitted numerical factors and gradual freeze-out must be included before using a more accurate neutron abundance.

At neutrino decoupling, the interaction rate is comparable with the Hubble parameter. Equating the supplied rates givesHere the Planck mass is the unreduced $m_P=G^{-1/2}\simeq 1.22\times 10^{19}\mathrm{GeV}$ in natural units, $G_F\simeq 10^{-5}{\mathrm{GeV}}^{-2}$, and the effective number of relativistic energy degrees of freedom is $g_{\ast }=10.75$. The supplied rate gives an order-of-magnitude neutrino decoupling estimate; a precise neutron-proton conversion temperature uses its own weak-rate coefficients.

The chemical equilibrium condition for $n+{\nu }_e\leftrightarrow p+e^{-}$ is ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. With negligible lepton chemical potentials, ${\mu }_n={\mu }_p$. The nonrelativistic Maxwell--Boltzmann number density therefore givesbecause $g_n=g_p=2$ and the small mass difference can be neglected in the prefactor. This is the neutron-proton chemical equilibrium relation.

For the deuterium equilibrium abundance, ${\mu }_D={\mu }_n+{\mu }_p$ because photons have zero chemical potential. The sign convention in the paper is important: $B=m_D-m_n-m_p<0$, while the positive nuclear binding energy is ${\Delta }_D=-B$. Using spin multiplicities $g_D=3$, $g_n=g_p=2$ givesWith $n_B=\eta n_{\gamma }$, $m_n\simeq m_p\simeq m_N$, $m_D\simeq 2m_N$ and photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, this becomesHere $X_D=n_D/n_B$ is a number fraction; the fraction of baryons bound in deuterium is $2X_D$. The formula estimates the onset of the deuterium bottleneck opening while depletion of free neutrons and protons remains small. Because $\eta \ll 1$, stable deuterium becomes abundant only far below the positive nuclear binding energy, not simply at $T={\Delta }_D$.

Let $T_f$ denote cosmological weak freeze-out of the neutron-to-proton ratio, $r_f=e^{-Q/T_f}$, and let $t_N$ mark efficient nuclear burning. Free-neutron decay after freeze-out givesAlmost every surviving neutron enters helium-4, whose atomic nucleus contains two neutrons and four baryons. Thus neutron-limited helium synthesis givesA longer neutron half-life leaves more neutrons available, increasing the primordial helium mass fraction. If the longer lifetime reflects weaker weak interactions, the accompanying earlier cosmological weak freeze-out can reinforce that effect.

During radiation domination, $H\propto \sqrt{G{g}_{\ast }}{T}^2$ whereas the weak conversion rate is proportional to $T^5$. Hence $T_f\propto (G{g}_{\ast }{)}^{1/6}$. Increasing the effective number of relativistic energy degrees of freedom raises $T_f$, increasing $r_f$, and shortens the interval for free-neutron decay after freeze-out. Both increase the primordial helium mass fraction while nuclear burning remains efficient. The same reasoning gives the primordial helium response to a varying gravitational constant: larger $G$ over this interval produces more helium, with fixed microscopic rates and the approximate Friedmann equation unchanged. Smaller $G$ has the opposite effect. A particular varying-$G$ model must specify its background expansion and time history; the simple trend must not be extrapolated to changes so large that nuclear burning itself fails.

Use natural units $\hslash =c=k_B=1$ and the unreduced Planck mass $M_P=G^{-1/2}$. Neutrino decoupling occurs when the weak interaction rate falls below the Hubble parameter. The weak-decoupling temperature estimate equates the given rate with the radiation-era expansion rate:Here $g_{\ast }$ counts the effective number of relativistic energy degrees of freedom, and the Fermi constant is $G_F\simeq 10^{-5}{\mathrm{GeV}}^{-2}$. The omitted reaction-rate coefficients make this an order-of-magnitude estimate: decoupling occurs at a temperature of order one MeV.

The charge-conserving channel is $n+{\nu }_e\leftrightarrow p+e^{-}$, together with $p+{\bar{\nu }}_e\leftrightarrow n+e^{+}$. The positron sign in the printed neutron-capture channel violates electric charge conservation. Chemical equilibrium requires ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$. Negligible lepton chemical potentials therefore give ${\mu }_n\simeq {\mu }_p$. Since neutrons and protons are nonrelativistic particles, their Maxwell-Boltzmann distribution givesBoth nucleons have two spin states, and their masses agree closely in the prefactor. Thus neutron-proton chemical equilibrium impliesThe heavier neutron is suppressed by a Boltzmann factor.

For the deuterium equilibrium abundance, retain the paper's convention $B=m_D-m_n-m_p<0$; the positive nuclear binding energy is ${\epsilon }_D=-B\simeq 2.22\mathrm{MeV}$. Chemical equilibrium of $n+p\leftrightarrow D+\gamma$ gives ${\mu }_D={\mu }_n+{\mu }_p$ because the equilibrium photon chemical potential vanishes. Using the nonrelativistic Maxwell--Boltzmann number density for each species and dividing givesHere $g_D=3$ and $g_n=g_p=2$. With the baryon-to-photon ratio $\eta =n_B/n_{\gamma }$, the photon number density $n_{\gamma }=2\zeta (3)T^3/{\pi }^2$, and $m_D\simeq 2m_N$, $m_n\simeq m_p\simeq m_N$, this becomesIn this convention $X_D=n_D/n_B$ counts nuclei, so its contribution to the fraction of baryons in deuterium is $2X_D$.

The deuterium bottleneck explains why helium-4 formation waits until $T$ is far below ${\epsilon }_D$. At higher temperature, there are vastly more photons than baryons; even the dilute energetic tail can destroy newly formed deuterium. A crude estimate balances ${\eta }^{-1}{e}^{-{\epsilon }_D/T}$ against one, giving $T\sim {\epsilon }_D/\log ({\eta }^{-1})\sim 0.1\mathrm{MeV}$ for $\eta$ of order $10^{-9}$. The power-law factors in the deuterium equilibrium abundance and the nuclear reaction rates refine this estimate. Once stable deuterium accumulates, further reactions can efficiently place nearly all surviving neutrons into helium-4.

Write $r_N=n_n/n_p$ at the onset of nuclear burning. Each helium-4 nucleus requires two neutrons and contains four baryons. Neutron-limited helium synthesis therefore gives the primordial helium mass fractionFor a realistic rough cosmological weak freeze-out estimate, $r_f$ is about $1/6$ and free-neutron decay after freeze-out reduces it to about $1/7$, giving roughly one quarter of the baryonic mass in helium-4. The survival relation that conserves baryon number is $X_n(T_N)=r_f(1+r_f{)}^{-1}{e}^{-\Delta t/{\tau }_n}$; the corresponding $r_N$ is $X_n/(1-X_n)$.

The rough neutrino decoupling rate by itself does not determine that one-quarter value. Substituting $T_d=1.8\mathrm{MeV}$ literally into $e^{-Q/T_d}$ would give $r\simeq 0.49$ and $Y_4\simeq 0.65$ before neutron decay. Neutron-proton chemical equilibrium persists to a lower temperature, approximately $0.7$--$0.8\mathrm{MeV}$ once the actual conversion rates are included. The rate scaling fixes the MeV scale; accurate abundances need the conversion coefficients and the delay through the deuterium bottleneck. This sequence is described in David Tong's discussion of nucleosynthesis.

Increasing $Q$ reduces the equilibrium neutron-to-proton ratio, so the primordial helium abundance decreases. Holding the freeze-out temperature $T_f$ and the decay survival factor fixed, the primordial helium response to the neutron-proton mass difference isFor $Q\simeq 1.29\mathrm{MeV}$ and $T_f=0.7$--$0.8\mathrm{MeV}$ this is about $0.85$--$0.87$, a reduction of roughly fifteen percent. A larger $Q$ also changes neutron decay and detailed weak rates; their quantitative effects require more microscopic information. The sign from the Boltzmann factor is already clear.

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.

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$.

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.

At weak freeze-out, photons, electrons, positrons and all four neutrino flavors are relativistic and approximately at the same temperature. The effective number of relativistic energy degrees of freedom changes fromThe expansion rate at fixed temperature is larger by $\sqrt{50/43}\simeq 1.078$. Equating the supplied weak interaction time with the supplied age gives $T_f^3=\sqrt{g_{\ast }}$ when $T_f$ is measured in MeV. ThereforeEarlier cosmological weak freeze-out at a higher temperature leaves more neutrons. With negligible chemical potentials, the neutron-to-proton ratio in equilibrium is $r_f=\exp (-Q/T_f)$, where $Q=m_n-m_p\simeq 1.293\mathrm{MeV}$. If almost all surviving neutrons enter helium-4, the primordial helium mass fraction before neutron decay is $Y_f=2r_f/(1+r_f)$. Literal use of the given rate normalizations yieldsEquivalently, the leading response is $\delta Y_f/Y_f=[Q/(T_f(1+r_f))]\delta T_f/T_f$. This is a roughly one-percentage-point absolute increase in the helium fraction, rather than a one-percent relative increase.

The deuterium bottleneck delays nuclear burning, allowing free-neutron decay after freeze-out. If its end temperature is $T_N$, and ${\tau }_n$ is the neutron mean life, a sharper instantaneous-freeze-out estimate conserves total baryon number and givesThis formula evolves the neutron fraction, not just the neutron-to-proton ratio, because each neutron decay also creates a proton. After electron annihilation, let $r_T=(4/11{)}^{1/3}$. The radiation counts relevant to the age at fixed $T_N$ areConsequently $t_N^{(4)}/t_N^{(3)}\simeq \sqrt{3.36/3.82}\simeq 0.939$. Less time is available for neutron decay, adding to the helium increase. As an explicit illustrative completion, $T_N=0.07\mathrm{MeV}$ and ${\tau }_n\simeq 880\mathrm{s}$ give an increase of approximately $0.012$ in this toy model. The precise helium fraction requires the nuclear-burning temperature, baryon abundance and properly normalized weak rates; the question supplies only approximate time laws. In particular, the large no-decay fractions in the table must not be interpreted as precise predictions for the observed primordial abundance. The robust effect is extra relativistic energy raises the primordial helium fraction, with an absolute change of order $10^{-2}$ for one thermally populated flavor.

A collisionless relativistic species redshifts with $T_{\nu }\propto a^{-1}$. The coupled plasma instead obeys cosmological entropy conservation, $g_{\ast s}{T}^3{a}^3=\mathrm{constant}$, and is heated relative to a decoupled species as massive particles annihilate. The effective number of relativistic energy degrees of freedom $g_{\ast }$ and effective number of relativistic entropy degrees of freedom $g_{\ast s}$ coincide while all species share a common relativistic temperature. Use $g_d:=g_{\ast s}(T_{\mathrm{dec}})\simeq 100$ for the Standard Model bath, excluding the already decoupled extra species.

Just before ordinary neutrino decoupling, the coupled bath contains two photon states, four electron/Positron states and six ordinary neutrino/antineutrino states. ThereforeThe extra-neutrino temperature relative to that bath is then $(43/(4g_d){)}^{1/3}$. After the ordinary neutrinos decouple, electron/Positron annihilation heats only the photon bath: its entropy count changes from $2+(7/8)4=11/2$ to $2$. Both neutrino populations miss that last heating, soThis is the two-stage entropy dilution of an early-decoupled neutrino. Numerically $T_{\nu ,0}\simeq 6.8\times 10^{-14}\mathrm{GeV}$. Relative to ordinary relic neutrinos it is cooler by $(43/(4g_d){)}^{1/3}\simeq 0.48$.

Let $g_{\nu }$ count the extra neutrino's populated helicity and antiparticle states. For one conventional neutrino plus antineutrino species $g_{\nu }=2$, the same state count as photons. The Fermi-Dirac distribution supplies the $3/4$ number-density factor, soThus a two-state species has about $2.9\%$ of the photon number density. Stating $g_{\nu }$ explicitly also covers a different choice of populated states.

The effective number of relativistic energy degrees of freedom count two photon polarizations, four electron/positron spin states, and six neutrino/antineutrino states, with the Fermi-Dirac distribution factor $7/8$ for fermions. Treating these species as relativistic and at a common temperature gives $\boxed{g_{\ast }=2+(7/8)(4+6)=10.75}$.

The mass-equivalent radiation mass density is $\rho =({\pi }^2/30)g_{\ast }(k_{B}T{)}^4/({\hslash }^3{c}^5)$. Therefore the Friedmann equation givesAt $T_0=10^{10}\mathrm{K}$, $H_0\simeq 0.503{\mathrm{s}}^{-1}$. Equating the specified weak interaction rate to $H_0(T/T_0{)}^2$ gives $(T_{\ast }/T_0{)}^3=H_0/(1{\mathrm{s}}^{-1})$ and henceThe approximate equilibrium neutron-to-proton ratio, neglecting chemical potentials, is $r_{\ast }\simeq \exp [-(m_n-m_p)c^2/(k_B{T}_{\ast })]$, about $0.152$ using the mass difference $1.293\mathrm{MeV}$. This is a freeze-out estimate rather than a detailed reaction-network computation.

After weak cosmological weak freeze-out, free neutrons undergo beta-minus decay; neutrinos decouple, and electron-positron annihilation heats the photons relative to the neutrinos. Initially the abundant energetic photons destroy deuterium almost as soon as it forms. Once the deuterium bottleneck ends, nuclear reactions rapidly pass through deuterium, tritium and helium-three to helium-4. Reactions eventually become too slow compared with cosmic expansion; nearly all surviving neutrons are bound in helium-four. Its mass fraction is approximately $Y_4=2r/(1+r)$, where $r$ is evaluated when nuclear assembly begins.

Increasing $g_{\ast }$ increases the Hubble parameter, raises $T_{\ast }$, increases the freeze-out neutron-to-proton ratio, and leaves less time for free-neutron decay after freeze-out. Both effects increase $Y_4$; decreasing $g_{\ast }$ reverses them. Most deuterium is only an intermediate product and is burned into helium, so its residual abundance and its absolute change are tiny compared with the helium mass fraction. This does not imply a negligible fractional change in residual deuterium. Faster expansion also shortens the time available to burn it, and its relative abundance can be appreciably sensitive to extra radiation. The distinction matters when interpreting the final wording of the question; joint deuterium and helium abundances constrain extra relativistic species, as discussed in the Particle Data Group nucleosynthesis review.

For a bosonic degree of freedom, the supplied integral givesFor a fermionic degree of freedom, useChanging variable $u=2y$ in the second integral givesThis is the energy-density case of the fermion-to-boson thermal integral ratio.

Summing the result of part (a) over all particle species therefore giveswhere the effective number of relativistic energy degrees of freedom is

In the ideal-gas high-temperature limit, the Standard Model has 28 boson states: 16 gluon states, eight electroweak gauge boson states and four Higgs field states. Its 90 fermion states are 72 quark states, 12 charged lepton states and six neutrino states. Counting a relativistic fermion state with weight $7/8$ gives the displayed effective number of relativistic entropy degrees of freedom. The effective number of relativistic energy degrees of freedom has the same value when all species share a temperature. Mass thresholds and interactions change the ideal-gas count.