Comoving entropy 2026-10-14
Comoving entropy is the entropy density multiplied by the physical volume factor for a fixed comoving region in an expanding homogeneous universe. In adiabatic expansion it is conserved. Separate decoupled radiation sectors conserve their own comoving entropy, so disappearance of massive species can heat the still-equilibrated sector relative to a decoupled species.
Cosmological entropy conservation Created 2026-09-24 Updated 2026-09-24
For an adiabatically expanding homogeneous universe, the entropy density satisfies
so the entropy in a comoving volume is constant.
The entropy density of a photon gas is its entropy per physical volume, . 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 with internal energy . The second uses the photon gas equation of state, and . If an independently constrained photon number gives a nonzero chemical potential, the general entropy density is and 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 , is
The 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 as transferred heat and 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, so
An irreversible adiabatic process can produce entropy even without heat exchange, so the word adiabatic alone is insufficient to infer .
For an analytic example, expand an initially equilibrated photon gas reversibly from to , keeping only the same two photon states of electromagnetic polarization. Its temperature halves, its energy density becomes , its entropy density of a photon gas becomes , and its total entropy stays . 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.
Figure 1.
Analytic reversible photon gas expansion: entropy conservation and pressure work
. Original fixed-species model with from 1 to 8 and no entropy production. The left panel shows , and on a logarithmic vertical scale. The right panel plots ; the shaded integral is . Total entropy remains constant even though both entropy density and internal energy decrease. Details.
For a computed scale at stipulated , . Combining this with the photon number density gives
This 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, and is fixed, so the mean number 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 . 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, : the conserved quantity for the combined sector is , 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]
For a single-component equilibrium system with extensive internal energy , scaling its entropy, volume and particle number together gives . The Euler theorem for homogeneous functions therefore gives . The first law of thermodynamics identifies these derivatives as , respectively, proving . Dividing by volume yields the entropy density formula . This step fixes the entropy normalization rather than merely identifying its differential. Differentiating the Euler relation and subtracting gives the Gibbs-Duhem equation . The argument requires extensivity; it is not an assumption about arbitrary nonadditive systems.
The baryon asymmetry is an excess of net baryon number over zero, conveniently measured by . 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 bosons and their antiparticles, with decay channels
The symbols in parentheses are branching fractions. The two channels carry baryon numbers and , so both cannot conserve one assigned baryon number for . Their mean baryon number is , while the mean for is . Consequently,
where is the asymmetry per decaying particle. CP violation permits 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 . 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 produces
Here 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.
Work in natural units , initially neglecting helium and excited hydrogen states. The reaction is in chemical equilibrium when , since equilibrium photons have zero chemical potential. Substitution of the nonrelativistic Maxwell-Boltzmann distribution for the three densities gives
For ground-state hydrogen use , , and in the translational prefactor. Writing its Hydrogen binding energy as , this becomes . Charge neutrality gives , while baryon conservation gives . Hence the hydrogen-only Saha equation is
The photon Bose-Einstein distribution supplies , where is the usual Riemann zeta function value in the thermal number integral. With baryon-to-photon ratio , substitution yields
This reproduces the requested numerical coefficient, with the printed interpreted as the thermal value. There is a normalization ambiguity in the PDF: it writes without defining . To reproduce this coefficient, must mean photon number density. If has its usual entropy density meaning, the Saha baryon-abundance normalization instead gives
The two forms agree through , with the entropy sector specified. A baryon-to-entropy ratio cannot simply be inserted into the photon-normalized expression.
As the universe cools, the Boltzmann factor favors neutral atoms and cosmological recombination sharply lowers . For an order-one ionized fraction the photon-normalized equation gives the balance
up to an order-one ionization-fraction factor. The small baryon abundance means that there are very many photons per baryon: even a small high-energy photon tail can ionize most atoms. Thus the recombination temperature is far below the microscopic binding scale , typically a few tenths of an eV for the standard abundance, rather than of order .
Photon decoupling is a related but distinct kinetic event. The Thomson scattering rate falls relative to as free electrons disappear, and photons begin free streaming. The precise last-scattering epoch depends on integrated optical depth, not only the local equality of rates. At still lower density the effective recombination rate per electron also drops below the expansion rate. Residual electron freeze-out leaves a small nonzero ionized fraction, larger than the vanishingly small equilibrium Saha prediction; radiative atomic bottlenecks mean the entire history cannot be obtained from equilibrium alone.
The temperature comparison in the last request needs an interpretation. If “reionisation temperature” means the microscopic thermal ionization scale, the large photon-to-baryon ratio explains why cosmological recombination requires . If it means later cosmic reionization, there is no universal thermal threshold: stars or other sources supply ionizing photons and can reheat gas. The cosmological photon background is actually colder at that later epoch, because . Its temperature therefore cannot literally be higher than the background temperature at recombination merely because hydrogen is reionized.
An adiabatic change of box size preserves the mode labels, whose momenta scale as . Since , differentiation gives . Occupation numbers are fixed under this slow mechanical change, so summing the changes in the individual particle energies gives
Comparison with at fixed entropy and particle number yields
The integration notation means summing over the original momentum-mode labels; it does not assume that the density in momentum space is unchanged after rescaling. In the ultrarelativistic limit , so and .
Adiabatic radiation obeys , giving energy density proportional to . Equilibrium photon energy density is proportional to , so . Therefore the photon number in a comoving volume, , is constant. Its entropy is also constant: photon entropy density is . Photon number is not a microscopic conserved charge, but equilibrium processes maintaining the spectrum do not alter this comoving-number scaling in adiabatic expansion; entropy-producing processes would invalidate the assumptions.
For a relativistic species at zero chemical potential, and , so its entropy density is proportional to . The Bose-Einstein distribution gives the bosonic normalization, while the fermionic thermal integral is of it. Species still in thermal equilibrium have ; freely streaming decoupled species can have another temperature and therefore contribute an extra factor when the total is expressed in units of . This gives the two terms and and their stated effective-degree-of-freedom weights.
After neutrino decoupling, the comoving entropy of the photon-electron-positron plasma is conserved separately from the neutrinos. Before electron-positron annihilation its effective entropy count is ; afterwards it is . Consequently
whereas stays equal to its pre-annihilation value. Hence
This is the instantaneous-decoupling, negligible residual-electron approximation used in the question. Heating of the still-coupled plasma, not cooling peculiar to neutrinos, creates the temperature difference.
The photon number density is obtained by integrating the Planck photon distribution. With ,
The integral converges, since its integrand is of order near zero and decays exponentially at infinity. Hence , with independent of temperature.
For equilibrium radiation undergoing adiabatic expansion, the total entropy in a comoving volume is constant. Since the entropy density is , this means is constant. A comoving volume scales as the cube of the scale factor, so
This agrees with the cosmological redshift of each photon frequency. The conclusion uses entropy conservation during the expansion.
When the number-changing reaction is in chemical equilibrium, its chemical balance condition is . Therefore zero chemical potential from number-changing reactions gives . Elastic scattering maintains kinetic equilibrium but imposes no independent condition on . The conclusion requires the number-changing rate to remain large compared with the Hubble parameter; merely listing an allowed reaction does not guarantee this after freeze-out.
Before thermal decoupling in cosmology, the sectors have a common temperature . For relativistic species at that temperature their entropy densities are
so
These are the effective number of relativistic degrees of freedom for entropy. The paper's agrees with for relativistic components at a common temperature. For one real scalar field, . At temperatures where all Standard Model species are relativistic, both counts are constant, and so is before decoupling; thresholds can change this equilibrium ratio if they occur while the sectors remain coupled.
After decoupling, adiabatic isolation gives separately conserved visible and dark entropy: and are individually constant. Thus even when either sector's entropy degrees of freedom change. Changes in particle content then alter their temperature ratio rather than the ratio of their comoving entropies. This assumes no late inter-sector transfer or entropy injection.
In natural units, an energy density and a pressure have mass dimension four, while a number density and an entropy density have mass dimension three. Since and each have mass dimension one,
For a relativistic species the supplied relation makes its entropy density proportional to . The fermion-to-boson thermal integral ratio gives the factor for each fermionic degree of freedom, so the electromagnetic plasma immediately before electron-positron annihilation in cosmology has
After annihilation only the two photon polarizations remain, so .
The neutrinos have already undergone neutrino decoupling. They therefore receive none of the electron-positron entropy and cool as . In the still-coupled electromagnetic plasma, cosmological entropy conservation gives . If and denote the common temperature and scale factor just before annihilation, then at a later scale factor ,
Dividing gives the Cosmic neutrino background temperature ratio
The entropy density is
For the homogeneous cosmological plasma, adiabaticity gives the cosmological entropy conservation equation
or equivalently in a comoving volume.
For a hydrogen-only gas with , and , the Saha equation is . Defining and using gives coefficient multiplying . If instead , with entropy density , that coefficient is . The two expressions agree because . The reference entropy sector must be specified, especially after particles have decoupled.
Separately conserved visible and dark entropy Created 2026-10-06 Updated 2026-10-14
After thermal decoupling in cosmology, two homogeneous sectors that evolve adiabatically without entropy injection satisfy and , hence their entropy density ratio is fixed. At a common relativistic temperature before decoupling, . Internal annihilations can change each sector's temperature and effective entropy degrees of freedom without changing its total comoving entropy.