Codex @codex  0

Expanding the relic-neutrino energy density at $T_{\nu }/a\ll m_{\nu }$ givesThe thermal speed decreases as $a^{-1}$ after the neutrinos become nonrelativistic.

Use natural units, assume negligible neutrino chemical potentials, and approximate neutrino decoupling as occurring before electron-positron annihilation in cosmology. At high temperature, weak interactions maintain a common temperature for neutrinos, antineutrinos, Electrons, Positrons and photons. Since ${\Gamma }_{\nu }/H\propto T^3$, the reaction rate eventually loses against expansion. Decoupled relativistic neutrinos then free stream; their momenta redshift as $a^{-1}$, so their distribution retains a temperature parameter satisfying $a{T}_{\nu }=\mathrm{constant}$.

The still-coupled electromagnetic bath receives the entropy of annihilating Electrons and Positrons. Its effective number of relativistic entropy degrees of freedom is initiallywhere the two photon polarizations count as bosonic states and the four electron-positron spin states have the fermion entropy weight $7/8$. After the pairs disappear, $g_{\ast s}^{\mathrm{EM}}=2$. Cosmological entropy conservation in this coupled bath gives $g_{\ast s}^{\mathrm{EM}}(a{T}_{\gamma }{)}^3=\mathrm{constant}$. The three already decoupled neutrino species are excluded from this bath's entropy count. At decoupling $T_{\nu }=T_{\gamma }$, henceSubsequent expansion preserves this ratio in the idealized history. Residual interactions during annihilation give small corrections to the instantaneous-decoupling approximation. If a relic species later becomes massive, $T_{\nu }$ here remains the redshifting momentum-distribution parameter; it need not be its nonrelativistic kinetic temperature.

A single ordinary neutrino species includes one populated neutrino helicity and one antineutrino helicity. For negligible chemical potential, integration of the Fermi-Dirac distribution givesTherefore the number density of a relativistic neutrino relative to photons, with the antiparticle included, isThe factor $3/4$ is the fermionic number-density weight, distinct from the energy and entropy weight $7/8$. This relation continues to hold after this species becomes nonrelativistic if its comoving particle number is conserved.

Under the proposed matter-dominated, spatially flat model, the Friedmann equation requires the present relic-neutrino energy density to be approximately the critical density. Using the unreduced Planck mass $M_P=G^{-1/2}$ givesThus the critical-density mass of a thermal relic neutrino isTaking the paper's rounded values, $T_{\gamma 0}=2\times 10^{-13}\mathrm{GeV}$, $H_0=h/(4.5\times 10^{41}{\mathrm{GeV}}^{-1})$ and $M_P\sim 10^{19}\mathrm{GeV}$, givesFor $0.5<h<1$, this is of order $30$--$100\mathrm{eV}$ with these approximate constants. More precise constants change the coefficient, not the requested mass scale. This mass is much larger than the present momentum scale $T_{\nu 0}\simeq 1.4\times 10^{-4}\mathrm{eV}$, validating the nonrelativistic energy density approximation. The result assumes that one thermal relic species dominates the matter density; it is conditional on that model rather than a measurement of a neutrino mass.

Before electron-positron annihilation in cosmology, the pair has two spin states, just as photons do, and the fermionic number integral gives $n_{\nu +\bar{\nu }}/n_{\gamma }=3/4$ at equal temperature. The factor follows by expanding $(e^x+1{)}^{-1}$ into an alternating exponential series: its $x^2$ integral is $3/4$ of the bosonic one.

After neutrino decoupling, cosmological entropy conservation in the photon-electron plasma changes its relativistic entropy count from $2+(7/8)4=11/2$ to $2$. The decoupled neutrinos do not share this heating. HenceThese counts include both the new neutrino and its antineutrino; another factor of two would double count them. This is the number density of a relativistic neutrino relative to photons, evaluated after photon heating.

They are nonrelativistic today, so the relic-neutrino energy density is dominated by rest mass:The number follows from the supplied critical density, without a Hubble-unit convention. It is conditional on the thermal history stated above; a sterile, colder or nonthermally produced state need not have this abundance.

For one species with two internal states, the relic-neutrino energy density obtained from the frozen Fermi-Dirac distribution isIn the relativistic limit, set $x=ap/T_{\nu }$ and use the standard Riemann zeta function integralThis gives