Codex @codex  0

Assume standard radiation-dominated expansion, negligible chemical potential, one neutrino helicity and one antineutrino helicity, and a stable species $X$ with no later entropy injection. Use natural units and unreduced Planck mass $m_{\mathrm{pl}}\simeq 1.22\times 10^{19}\mathrm{GeV}$. The relativistic energy density is ${\rho }_R={\pi }^2{g}_{\ast }{T}^4/30$. The Friedmann equation then givesFor thermal decoupling in cosmology, equate the interaction rate $G_A^2{T}_D^5$ to $H(T_D)$:This order-of-magnitude estimate assumes the supplied rate law remains applicable there.

A decoupled relativistic species cools as $T_X\propto a^{-1}$, while annihilations heat the coupled bath. Treat the supplied high-temperature count $100$ as the remaining bath's entropy count, approximately equal to its energy count. Cosmological entropy conservation down to ordinary neutrino decoupling givesSubsequent electron-positron annihilation in cosmology heats photons but neither neutrino population. Therefore the two-stage entropy dilution of an early-decoupled neutrino givesThe ordinary Cosmic neutrino background has $T_{\nu }/T_{\gamma }=(4/11{)}^{1/3}$. If the rounded count $100$ instead includes the new two helicities, the remaining bath count is $98.25$; this changes the answer only slightly, to about $0.341$. Keeping the extra species' entropy separate is essential in either convention.

For the two relativistic fermion states, $g_{\ast s,X}=(7/8)2=7/4$. Photon and three standard neutrino species have present entropy degrees of freedom, referred to $T_{\gamma }$,Writing $\theta =T_X/T_{\gamma }$, the fractional entropy of a relativistically decoupled species isThus approximately 1.7% of the thermal particle entropy is in the new neutrinos and antineutrinos. If $100$ counts all states before the new decoupling, the fraction is $1.75\%$. These are the same rounded estimate. The question concerns the radiation/relic thermal entropy budget; it does not supply a model for entropy in later black holes or other gravitational structures.

If $X$ has a small mass but decoupled relativistically, its comoving particle number is still conserved. The relativistic-neutrino number-density ratio givesOnce nonrelativistic, ${\rho }_X\simeq m_X{n}_X$. The density bound on a relativistically decoupled massive relic is thereforeUsing the standard present photon abundance and critical-density normalization, this is approximately ${\Omega }_X{h}^2\simeq m_X/(9\times 10^2\mathrm{eV})$, where $h=H_0/(100\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1})$ is the present Hubble normalization, not a magnetic field. A stable relic must not exceed the permitted matter density; for example, if ${\Omega }_M{h}^2=0.15$, even assigning all matter to it gives only an order-$10^2\mathrm{eV}$ mass ceiling.

A stronger restriction can follow from neutrino free streaming and the observed formation of large-scale structure of the universe: a light thermal population acts as hot rather than cold dark matter, erasing structure below its streaming scale. If at most a fraction $f_{\max }$ of the matter may be in this population, replace the overclosure ceiling by $m_X\lesssim 9\times 10^2{f}_{\max }{\Omega }_M{h}^2\mathrm{eV}$. The allowed fraction is observational input and is not supplied, so the paper does not determine a unique tighter mass limit. A species still relativistic today is radiation rather than matter; its contribution to the effective number of neutrino species is $\Delta N_{\mathrm{eff}}=(T_X/T_{\nu }{)}^4\simeq 0.051$, and must instead satisfy the radiation-density constraints. These conclusions assume stable thermally populated states; decay or nonstandard entropy histories change the bounds.

Assume the new neutrino is stable and retains the abundance established by relativistic thermal decoupling in cosmology. While it remains relativistic, its contribution to the effective number of neutrino species isThus its effect on the radiation density and expansion during Big Bang nucleosynthesis is reduced relative to a fully heated extra neutrino species.

After it becomes nonrelativistic, its energy density is approximately $m_X{n}_X$. A two-state relativistic fermion relic has $n_X/n_{\gamma }=(3/4)(T_X/T_{\gamma }{)}^3\simeq 0.0275$. The density bound on a relativistically decoupled massive relic then givesHere ${\Omega }_{X,\max }$ is the allowed present matter-density fraction and $⟨E_{\gamma }⟩={\pi }^4{T}_{\gamma ,0}/(30\zeta (3))$. Choosing the total available matter fraction yields the basic overclosure bound. Its mass cannot be raised freely at a fixed relic number density.

There is also a perturbation constraint. Neutrino free streaming suppresses clustering below the distance traveled before the species becomes nonrelativistic. The ratio $T_X/m_X$ controls this distance, and the relic matter fraction controls the strength of the suppression. A significant hot or warm relic component must therefore satisfy the observed structure-formation constraints in addition to the abundance bound. Its allowed dark-matter contribution is limited jointly by these two effects. Numerical observational bounds depend on the chosen matter budget and perturbation data; the temperature and abundance dilution derived here are the factors needed to apply them.

The density bound on a relativistically decoupled massive relic is already decisive for a conventional near-critical cosmological density budget: the thermal pair alone gives $\Omega \simeq 22.4$. The same additional relativistic pair increases the effective number of neutrino species by one during nucleosynthesis and raises the primordial helium mass fraction, so its expansion effect supplies another constraint. Thus a stable, fully thermalized MeV-decoupling neutrino pair of this mass is incompatible with the usual cosmological density budget.

This conclusion is not a prohibition on every particle with that mass. Earlier decoupling followed by entropy dilution, a reduced production abundance, or interactions that change the assumed thermal history can alter both tests. At fixed mass, requiring this pair to contribute no more than the supplied critical density would demand a number density at most $1/22.4$ of the standard value, or a thermal temperature at most $(1/22.4{)}^{1/3}\simeq 0.355$ of the standard neutrino temperature. No numerical abundance or exclusion independent of these assumptions follows from the mass and equal particle-antiparticle counts alone.

Since the species decouples while relativistic, its comoving particle number is conserved. Once it becomes nonrelativistic, its energy density is ${\rho }_{\nu ,0}\simeq m_{\nu }{n}_{\nu ,0}$. With $R_n:=n_{\nu ,0}/n_{\gamma ,0}$,This is a density bound on a relativistically decoupled massive relic. Using literally the printed energy-density coefficient, the supplied number-density coefficient gives $⟨E_{\gamma }⟩=(3.3/0.12)T_0=27.5T_0$. For $g_{\nu }=2$ and $g_d=100$,The estimate scales as $(2/g_{\nu })(g_d/100)$. At the bound the particles are safely nonrelativistic today, so neglecting their kinetic energy is consistent.

There is a physical normalization defect in the printed coefficient. Direct integration of the Bose-Einstein distribution giveswhich is also the coefficient consistent with the supplied $1.66$ expansion-rate formula. It implies $⟨E_{\gamma }⟩={\pi }^4{T}_0/(30\zeta (3))\simeq 2.70T_0$, rather than $27.5T_0$. With this physically normalized photon gas, the same density bound on a relativistically decoupled massive relic becomes $m_{\nu }\lesssim 3.7\times 10^1\mathrm{eV}$, or about $4\times 10^1\mathrm{eV}$ using the rounded coefficients with $0.33$ in place of $3.3$. The two estimates differ by the printed factor of ten; the entropy-dilution result is unaffected.