Codex @codex  0

A massive decoupled sector enters a cannibal phase when its particles are nonrelativistic but rapid number-changing reactions maintain chemical equilibrium and zero chemical potential from number-changing reactions. Rest-mass energy released by $3\to 2$ reactions heats the remaining particles, so the sector cools much more slowly than an isolated number-conserving nonrelativistic gas. The temperature laws follow directly from cosmological entropy conservation. An example is discussed in Cannibal Dark Matter.

A comoving volume has fixed boundaries in comoving coordinates. Its physical volume grows as the cube of the scale factor. In a homogeneous adiabatic gas, cosmological entropy conservation makes $sV$ constant.

For an extensive equilibrium system, the thermodynamic Euler relation gives $E=TS-PV+\mu N$. At zero chemical potential, divide by $V$ to obtain $s=(\rho +P)/T$. Equivalently, $dP=(\rho +P)dT/T$ makesThe first law of thermodynamics identifies this as $dS$; the thermodynamic Euler relation fixes the entropy normalization. For adiabatic expansion satisfying the cosmological perfect-fluid continuity equation, the numerator on the right vanishes, proving cosmological entropy conservation.

The entropy density of a photon gas is its entropy per physical volume, $s=S/V$. For homogeneous equilibrium vacuum radiation with two states of electromagnetic polarization and zero photon chemical potential,The first equality is the entropy density at zero chemical potential, obtained from the thermodynamic Euler relation $U=TS-pV+{\mu }_{\gamma }{N}_{\gamma }$ with internal energy $U=uV$. The second uses the photon gas equation of state, $u=a_R{T}^4$ and $p=u/3$. If an independently constrained photon number gives a nonzero chemical potential, the general entropy density is $(u+p-{\mu }_{\gamma }{n}_{\gamma })/T$ and $u=a_R{T}^4$ need not hold. This distinction explains why the absence of a conserved number matters for entropy as well as for the spectral shape. ^[ref]

A short first law of thermodynamics check, for neighboring equilibrium states with ${\mu }_{\gamma }=0$, isThe thermodynamic Euler relation fixes the additive normalization. Equivalently, this radiation entropy tends to zero with temperature at fixed volume; no contribution from the containing walls is included. Interpreting $TdS$ as transferred heat and $pdV$ as boundary work requires a reversible thermodynamic process. For a reversible thermodynamic process that is also an adiabatic process, there is no entropy production or entropy exchange, soAn irreversible adiabatic process can produce entropy even without heat exchange, so the word adiabatic alone is insufficient to infer $T^{3}V=\text{constant}$.

For an analytic example, expand an initially equilibrated photon gas reversibly from $V_0$ to $8V_0$, keeping only the same two photon states of electromagnetic polarization. Its temperature halves, its energy density becomes $u_0/16$, its entropy density of a photon gas becomes $s_0/8$, and its total entropy stays $S_0$. The internal energy halves rather than remaining fixed: the missing half leaves as work against radiation pressure,This is a stipulated reversible radiation model, not a report of a piston experiment. Maintaining a real cavity along this path requires controlling wall heat exchange and including the wall heat capacity in any actual apparatus balance.

For a computed scale at stipulated $T=2.725\mathrm{K}$, $s\simeq 2.04122\times 10^{-14}\mathrm{J}{\mathrm{K}}^{-1}{\mathrm{m}}^{-3}$. Combining this with the photon number density givesThis is the total equilibrium radiation entropy divided by its mean photon number, not a separately assigned entropy carried by each labeled particle. During the ideal expansion, $n_{\gamma }\propto T^3$ and $T^{3}V$ is fixed, so the mean number $N_{\gamma }=n_{\gamma }V$ also stays constant. This follows from the path; it was not imposed as the independent number constraint used to establish chemical potential.

The beautiful connection is that expansion lowers each typical photon energy and dilutes the radiation in just the proportions needed to retain its equilibrium spectral shape and total entropy. For the cosmic microwave background, homogeneous expansion of the universe has physical volume proportional to the cube of the scale factor; fixed-species cosmological entropy conservation therefore gives $T\propto a^{-1}$. This describes preservation of an already established Planck photon distribution and connects the FIRAS spectral observation in photon gas to an entropy prediction; FIRAS did not perform a direct entropy measurement. When other equilibrated species disappear and transfer energy into radiation, use the effective number of relativistic entropy degrees of freedom, $g_{\ast s}$: the conserved quantity for the combined sector is $g_{\ast s}{T}^{3}V$, under the corresponding equilibrium and no-entropy production assumptions. A changing species population, separate decoupled sectors or entropy production requires an appropriate entropy balance rather than the pure photon gas formula. ^[ref]

At fixed scale factor and unchanged relativistic counts, instantaneous conversion of a relic into thermal radiation gives the displayed entropy-density enhancement. It exceeds one because the gas receives energy from outside its prior equilibrium sector. This does not contradict cosmological entropy conservation, whose source-free reversible assumptions fail during decay and thermalization.

The baryon asymmetry is an excess of net baryon number over zero, conveniently measured by ${\eta }_B=(n_B-n_{\bar{B}})/n_{\gamma }$. Ordinary matter contains abundant baryons, whereas no comparable primordial antibaryon population is observed. Large nearby domains of matter and antimatter would produce radiation through particle-antiparticle annihilation at their boundaries. Antiparticles produced in energetic collisions are not evidence for a cosmological antimatter domain. The small nonzero baryon-to-photon ratio inferred from Big Bang nucleosynthesis is independently tested by baryon-density consistency between nucleosynthesis and the microwave background.

Starting with zero net baryon number, the Sakharov conditions require baryon-number violation, C and CP violation, and departure from thermal equilibrium. Conservation of baryon number would keep the initial value zero. C symmetry or CP symmetry would pair baryon production with equal antibaryon production; CP rate cancellation in baryogenesis explains why C violation alone does not suffice. In thermal equilibrium with zero imposed charge chemical potentials, CPT symmetry forbids the required net baryon number. An asymmetric rate in one reaction must not be mistaken for net production after all reverse reactions are included.

A simple example is heavy-boson decay baryogenesis. Take equally many heavy $X$ bosons and their antiparticles, with decay channelsThe symbols in parentheses are branching fractions. The two $X$ channels carry baryon numbers $2/3$ and $-1/3$, so both cannot conserve one assigned baryon number for $X$. Their mean baryon number is $r-1/3$, while the mean for $\bar{X}$ is $-\bar{r}+1/3$. Consequently,where $ϵ$ is the asymmetry per decaying particle. CP violation permits $r\ne \bar{r}$ without violating the equality of total decay widths required by CPT symmetry. Microscopically a direct CP asymmetry requires interference between decay amplitudes with both relative weak phases and absorptive phases; specifying complex couplings alone would not establish unequal partial decay widths.

For radiation domination, a sufficient decay condition is ${\Gamma }_X\lesssim H(T=M_X)=1.66\sqrt{g_{\ast }}{M}_X^2/m_P$. The heavy population can then survive as the equilibrium abundance drops, before decaying while inverse production is suppressed. Its actual abundance differs from its thermal equilibrium value. If subsequent baryon-number washout is negligible and cosmological entropy conservation holds, an initial total abundance $Y_X=(n_X+n_{\bar{X}})/s$ producesHere $s$ is the entropy density. Residual baryon-number washout multiplies this expression by an efficiency factor below one. The example realizes all three Sakharov conditions; it illustrates a mechanism, rather than identifying the actual cosmological source of the baryon asymmetry.

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.

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

Write the photon energy density as $ϵ=a_R{T}^4$, where $a_R=4\sigma /c$ by the Stefan–Boltzmann law. The photon gas has pressure $P=ϵ/3$ and zero photon chemical potential. The first law of thermodynamics therefore givesDividing by $T$ and integrating gives the entropy density of a photon gas:The integration constant is fixed to zero at zero temperature, consistently with the Third law of thermodynamics.

The inequality $\Gamma \gg H$ says interactions act much faster than the expansion time, maintaining local thermal equilibrium. For adiabatic reversible expansion with no additional source of entropy, cosmological entropy conservation fixes the entropy in a comoving volume. Since that volume scales as $a^3$, and the relativistic degrees of freedom are fixed here, $T^3{a}^3$ is constant. ThusDuring radiation domination, the Friedmann equation gives $H^2\propto \rho \propto T^4$, and hence $\boxed{H\propto T^2}$. Fast interactions alone would not exclude an additional entropy-producing process; changing relativistic species likewise changes the simple temperature scaling.

Count one neutrino species as the usual two thermally populated particle-antiparticle helicity states, with negligible chemical potential. Its entropy weight is $g_{\ast s,X}=(7/8)2=7/4$. In the ideal high-temperature Standard Model, the bosonic count is $16+8+4=28$ from gluons, electroweak gauge bosons and the Higgs field. The fermionic count is $72+12+6=90$ from quarks, charged leptons and the three ordinary neutrino species. Hence the relativistic entropy content of the Standard Model is $g_{\ast s,D}=28+(7/8)90=106.75$ for the remaining bath. Before the new species leaves, the total count is $108.5$; its separately cooling states are excluded from the subsequent bath count.

Near ordinary neutrino decoupling, the still-coupled bath hasfrom photons, electron-positron pairs and the three ordinary neutrinos. After those neutrinos decouple, the electromagnetic bath falls from $2+(7/8)4=11/2$ to two when the pairs annihilate. Cosmological entropy conservation in that bath gives $T_{\nu }/T_{\gamma }=(4/11{)}^{1/3}$.

The new species has cooled as $T_X\propto a^{-1}$ since its earlier decoupling, whereas annihilation of Standard Model species heated the remaining bath. The temperature ratio of an early-decoupled neutrino is thusThese are ideal-gas estimates with the high-temperature count implicit in the question; precise electroweak thresholds change them slightly.

Entropy degrees of freedom weight temperature ratios with their cube, while energy degrees of freedom use their fourth power. The remaining Standard Model entropy today, in photon-temperature units, is $g_{\ast s,0}^{\mathrm{SM}}=2+(7/8)6(4/11)=43/11$. The new contribution is $(7/4)(T_X/T_{\gamma }{)}^3$. Therefore the fractional entropy of a relativistically decoupled species isAbout 1.6% of the total entropy is in the new species. Separate comoving entropy conservation also derives this fraction directly at decoupling. A model populating four helicity states instead would require the corresponding doubled species count.

For an equilibrium relativistic bath, $P=\rho /3$ and $\rho \propto g_{\ast s}{T}^4$, so $s=(\rho +P)/T\propto g_{\ast s}{T}^3$. Cosmological entropy conservation gives $g_{\ast s}{T}^3{a}^3=\text{constant}$. A decoupled relativistic species instead has $a{T}_{\mathrm{relic}}=\text{constant}$, so annihilation heats the bath relative to that relic.

Use the sharp relativistic degree-counting approximation specified by the temperature interval. Just before muon annihilation, the coupled bath contains photons, electron-positron pairs, muon-antimuon pairs and three neutrino-antineutrino species:After muon annihilation, neutrinos are still coupled and $g_{\ast s}^{\mathrm{after}}=2+(7/8)(4+6)=43/4$. The decoupled H particles are excluded from these bath counts. Equality of bath and relic temperatures at decoupling therefore gives the relic cooling across muon annihilation

By the later electron-mass epoch, neutrinos have decoupled. Electron-positron annihilation heats only the electromagnetic bath, reducing its entropy degree count from $2+(7/8)4=11/2$ to $2$. Photons are heated by $(11/4{)}^{1/3}$ relative to both neutrinos and H particles. After this transition,The neutrino-to-photon ratio is $(4/11{)}^{1/3}$, while the H-particle-to-neutrino ratio remains $(43/57{)}^{1/3}$. The printed $0.59\mathrm{MeV}$ marks the electron-positron annihilation epoch approximately; it is not an exact annihilation threshold.

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.

Let $r_{\nu }=T_{\nu ,0}/T_0$, and retain the assumption that these neutrinos remain massless. Constancy of the neutrino degeneracy parameter gives $|{\mu }_{\nu ,0}|=|{\xi }_{\nu }|r_{\nu }{T}_0$. The condition on the critical density therefore yields the critical-density bound on massless neutrino degeneracyFor the standard instantaneous-neutrino decoupling approximation before electron-positron annihilation in cosmology, cosmological entropy conservation heats photons relative to the decoupled neutrinos. The electromagnetic entropy degrees of freedom change from $2+(7/8)4=11/2$ to $2$, givingThis numerical value requires that thermal-history assumption; the bound in terms of $r_{\nu }$ is the general answer if large degeneracy changes decoupling or subsequent entropy transfer. Setting $T_{\nu ,0}=T_0$ would instead give $26.2$. For $N$ equally degenerate flavors with the same temperature, the leading bound becomes $36.8N^{-1/4}$ under the standard temperature assumption. The neglected thermal corrections make a small change near this loose energy-density limit; this is not a detailed Big Bang nucleosynthesis bound.

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.

Consider a species in thermal equilibrium, with energy density and pressure depending on temperature and with vanishing chemical potential, as in the supplied thermodynamic relation. Put $U=\rho (T)V$. The first law of thermodynamics givesSince $dP/dT=(\rho +P)/T$,The differential is therefore exact:After fixing the irrelevant additive entropy zero, the entropy density at zero chemical potential isMore precisely the first law leaves an additive constant in $S$; extensivity fixes that constant for the usual entropy-density normalization. A nonzero chemical potential would instead require $Ts=\rho +P-\mu n$, so the stated formula is not a universal identity for a decoupled massive species with conserved particle number.

For a reversible thermodynamic process that is an adiabatic process in a comoving volume, $\dot{V}=3HV$ and the cosmological perfect-fluid continuity equation gives $\dot{\rho }=-3H(\rho +P)$. Differentiating the extensive expression explicitly, or using the exact differential above, givesThis is cosmological entropy conservation for the closed equilibrium gas with no entropy-producing energy injection.

For radiation in cosmology, $P=\rho /3$, so $s=4\rho /(3T)$. Comparing $s=c_1{g}_{\star }{T}^3$ and $\rho =c_2{g}_{\star }{T}^4$ yields the coefficient ratioWhen all relativistic species share the same temperature and count, the familiar normalizations are $c_1=2{\pi }^2/45$ and $c_2={\pi }^2/30$. More generally the energy and entropy effective counts can differ; the common $g_{\star }$ here uses the assumptions stated in the question.

The instantaneous-decay approximation keeps the scale factor and physical volume fixed across the decay. Energy conservation and rapid thermalization giveMeanwhile the abundance relation gives $s_{-}={\rho }_{\varphi ,-}/(Ym)\simeq 3M_{\mathrm{Pl}}^2/(Ym{t}_d^2)$. Dividing yields the entropy-injection factorThere is a numerical error in the printed coefficient: part (c) omits the factor $1/3$. The denominator $s_{-}=3M_{\mathrm{Pl}}^2/(Ym{t}_d^2)$ supplies that factor, while the radiation energy contributes $3^{3/4}$ to the numerator. Thus the consistent coefficient is $c_1{3}^{-1/4}{c}_2^{-3/4}$, rather than $c_1{3}^{3/4}{c}_2^{-3/4}$. Both cannot follow from the stated part (b) using the same decay-time approximation.

Its increase is especially transparent without dropping the original radiation. Let $R={\rho }_{\varphi ,-}/{\rho }_{r,-}$. At unchanged $g_{\star }$, $T_{+}^4/T_{-}^4=1+R$, soRelic domination means $R\gg 1$, so the leading estimate is much larger than one. This is entropy production by decay of a dominant relic.

There is no conflict with part (a). The equilibrium gas before the decay was a separately conserved, reversibly expanding system. During decay it receives energy from the nonthermal relic, and the decay products thermalize irreversibly. Its cosmological perfect-fluid continuity equation now has an energy-injection source, so $d({\rho }_{r}V)+P_{r}dV$ does not vanish. The assumptions behind cosmological entropy conservation therefore fail during the event. Moreover the zero-chemical-potential equilibrium expression in part (a) cannot be imposed on the decoupled relic as though it were already part of that same equilibrium radiation bath. The Second law of thermodynamics permits, and here requires, the resulting entropy increase.

Take a homogeneous comoving cell of physical volume $V_c=a^3{V}_0$ and total entropy $S=s{V}_c$. At vanishing chemical potential, the reversible first law of thermodynamics readsFor an extensive equilibrium system, the thermodynamic Euler relation is $E=TS-P{V}_c+\mu N$. At $\mu =0$ this givesThe supplied pressure identity verifies consistency directly: $dP=(\rho +P)dT/T$ impliesThe thermodynamic Euler relation fixes the additive normalization that a differential identity alone would leave undetermined.

For adiabatic cosmological expansion, the cosmological perfect-fluid continuity equation gives $\dot{\rho }+3H(\rho +P)=0$, while ${\dot{V}}_c=3H{V}_c$. HenceThis proves cosmological entropy conservation in the stated perfect-fluid approximation. Kinetic equilibrium alone does not guarantee adiabatic evolution: entropy injection, dissipative stresses or out-of-equilibrium reactions would require additional source terms. Those effects are implicitly excluded in this calculation.

During the cannibal phase of a decoupled sector, cosmological entropy conservation and the leading result of part (c) giveDifferentiating and using $T_{\varphi }=m_{\varphi }/x$ yields logarithmic cooling during cannibalism:Here $a_{\ast }$ absorbs the entropy normalization, and the logarithmic expression is a large-$x$ asymptotic relation. With fixed visible entropy degrees of freedom the photon bath has $T_{\gamma }\propto a^{-1}$. Thus for $x\gg 1$ the scalar particles cool much more slowly than the photons and become relatively hotter.

The physical reason is that $3\varphi \to 2\varphi$ converts rest-mass energy into kinetic energy of the remaining particles. Elastic scattering distributes this energy, offsetting much of the cooling caused by expansion. The total comoving entropy remains constant; the number density does not merely dilute as $a^{-3}$. This explanation applies only while number-changing chemical equilibrium persists. If those reactions freeze out, an isolated number-conserving nonrelativistic gas instead has $T_{\varphi }\propto a^{-2}$ under adiabatic expansion, which is faster cooling than the photons. The paper supplies no rates, so that transition time cannot be determined.

Let $V$ be the physical volume of a region with fixed comoving coordinates, so $V\propto a^3$, and write $U=\rho V$. The first law of thermodynamics is $dU=TdS-PdV+\mu dN$. Use extensive quantities: under replication, $U(\lambda S,\lambda V,\lambda N)=\lambda U(S,V,N)$. The Euler theorem for homogeneous functions givesSetting $\mu =0$ yieldsThis is the entropy form of the Gibbs-Duhem equation; it fixes the normalization by extensivity, rather than adding an arbitrary volume-independent entropy constant.

For equilibrium expansion of the same region, the first law and the cosmological perfect-fluid continuity equation implyTherefore $\boxed{\dot{S}=0}$, establishing cosmological entropy conservation for adiabatic equilibrium evolution. Here “comoving volume” means a patch moving with the cosmological fluid; its physical volume $V$ changes. The corresponding coordinate volume is constant. The local TeX's $\dot{V}/P$ in the continuity-equation hint is a transcription error: the PDF correctly has $\dot{V}/V$.