Codex @codex  0

Combining the degenerate neutrino and antineutrino energy density with a critical density bound constrains the neutrino degeneracy parameter through the fourth root of the available energy budget. The numerical bound depends on $T_{\nu ,0}/T_{\gamma ,0}$ and the number of degenerate flavors. The standard ratio $(4/11{)}^{1/3}$ assumes neutrino decoupling before electron-positron annihilation in cosmology; it should be stated as a thermal-history assumption if the degeneracy is very large. The massless formula is inappropriate for a species that is nonrelativistic today.

Use $\hslash =c=k_B=1$ and count one populated helicity per neutrino or antineutrino. For massless particles, the Fermi-Dirac distribution givesWhen ${\mu }_{\nu }/T_{\nu }\gg 1$, the occupied states form an almost sharp Fermi sea. Replacing the distribution by $\Theta ({\mu }_{\nu }-p)$ givesThe relative finite-temperature correction is of order $(T_{\nu }/{\mu }_{\nu }{)}^2$.

Chemical equilibrium gives ${\mu }_{\bar{\nu }}=-{\mu }_{\nu }$. For positive large ${\mu }_{\nu }$, the corresponding antineutrino distribution is dilute, soIt is exponentially negligible. The neutrino degeneracy parameter ${\xi }_{\nu }={\mu }_{\nu }/T_{\nu }$ is conserved in the assumed adiabatic massless evolution: redshifting preserves the distribution with both $T_{\nu }$ and ${\mu }_{\nu }$ proportional to $a^{-1}$.

There is a sign qualification in the printed absolute-value formula. If ${\mu }_{\nu }<0$, the distribution of neutrinos themselves is exponentially suppressed, and the antineutrinos, with positive chemical potential, form the degenerate sea. In either case the degenerate neutrino and antineutrino energy density isThus the formula with $|{\mu }_{\nu }|$ describes the dominant member of the pair, or their leading total, rather than the named neutrino distribution for both signs. No extra factor of two is present: only one member is densely occupied. Additional populated internal states multiply the result by their degeneracy.