 Solution

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.