Before cosmological recombination, photon-baryon sound speed controls the baryonic comoving Jeans length. The tightly coupled photon-baryon fluid has , where . Using total background density in the instantaneous dynamical Jeans estimate, the scale grows as during radiation domination, approaches a plateau in strongly baryon-loaded matter domination, and drops when photon decoupling removes photon pressure support. Subsequently, adiabatic monatomic-gas cooling gives and a comoving Jeans length proportional to . Residual Compton scattering can delay that cooling law.
Cosmic microwave background last-scattering surface Created 2026-09-28 Updated 2026-10-14
The cosmic microwave background last-scattering surface is the shell from which the photons observed in the cosmic microwave background last scattered around photon decoupling.
Cosmic neutrino background last-scattering surface Created 2026-09-28 Updated 2026-10-14
The cosmic neutrino background last-scattering surface is the locus from which relic neutrinos that reach an observer underwent their final appreciable scattering. Although neutrino decoupling occurs before photon decoupling, a massive neutrino travels slower than light after becoming nonrelativistic, so its last-scattering surface need not be farther away than the cosmic microwave background last-scattering surface.
For small baryon fraction the CDM growing mode is , and its gravity forces
For an adiabatic cosmological perturbation to follow the unsuppressed growing mode, pressure must be below the gravitational forcing frequency: . With physical wavelength , this becomes
This is Jeans support versus forced baryon response. The wording “only grow” must refer to unsuppressed gravitational structure growth; it is not a theorem that all shorter-wavelength baryon perturbations are absent or constant. In the pressure-dominated regime a forced response remains. In particular, a falling sound speed can make this small response grow even before it catches up to CDM.
Before photon decoupling, the photon-baryon fluid has sound speed of order . Its Jeans scale is roughly a horizon scale, so smaller perturbations oscillate acoustically instead of freely following the growing CDM structures. After decoupling the baryon sound speed drops sharply. Above the new Jeans scale, pressure is negligible and subtraction of the CDM equation gives . For , the difference is a constant plus a decaying term; its ratio to the growing CDM contrast tends to zero. Baryons therefore catch up. Under the stipulated post-decoupling , the physical Jeans length scales as , the comoving scale as and the corresponding Jeans mass as .
For the requested instantaneous mass estimates, define the baryonic Jeans mass as the baryon content of a sphere of radius :
Using the supplied matter-era age gives s and . Then
The drop is by . A definition using the total gravitating mass instead of baryon content gives values about ten times larger. The pre-decoupling number is an acoustic Jeans order estimate, not an accurate application of the nonrelativistic equations from part a when .
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.
There are about photons per baryon. Even when the typical photon energy is below the Hydrogen binding energy, the high-energy tail can provide abundant ionizing photons. The nonrelativistic Electron's translational phase space also favours the ionized state. Quantitatively, the Saha equation contains both suppressing factors and . Recombination needs their tiny product to be overcome by , so must be several tens, not of order one.
For an illustrative , the supplied baryon-to-photon relation gives . With , direct evaluation gives
Thus the equilibrium ionization fraction reaches a few per cent just below , while at it is essentially unity. The exact temperature depends weakly on the baryon abundance and on the fraction used to define cosmological recombination. This is small baryon abundance delays hydrogen recombination. Recombination, photon decoupling and residual electron freeze-out are related but distinct: the last two involve rates and cannot be determined by chemical equilibrium alone.
Figure 1.
Hydrogen-only Saha ionization fraction during cooling, with the three-per-cent crossing near 0.28 electronvolts for baryon-to-photon ratio 5.36 times ten to the minus ten
. Details.
Put . The synchronous perfect-fluid density equation and Cosmological Euler equation in synchronous gauge give
Differentiate the first equation and use the second and the scalar trace metric equation. Before making a nonrelativistic approximation, the result is
For , discard relative pressure corrections to the background gravitational and damping terms, while retaining because a large wavenumber can compensate a small sound speed. Matter domination and the flat Friedmann equation give . Thus the leading nonrelativistic density equation in synchronous gauge is
The approximation sign records the small finite- terms just displayed; the printed Newtonian-form equation is not an exact relativistic equation at nonzero .
The comoving Jeans wavenumber balances pressure and gravity: . Defining physical wavelength as , the physical Jeans length is
For wavelengths longer than this scale, gravity can overcome pressure and density modes grow; shorter modes are pressure-supported acoustic oscillations, with expansion modifying both growth and damping. Before cosmological recombination, the tightly coupled photon-baryon fluid has substantial photon pressure, , where . The baryonic Jeans scale is correspondingly large, so baryons oscillate instead of freely collapsing on small scales. Cold dark matter has a much smaller pressure scale and can supply gravitational potential wells. After photon decoupling, baryons lose photon pressure support and their sound speed falls to its thermal-gas value; the baryon Jeans length across recombination drops, allowing baryons to fall into those wells. The earlier photon-coupled fluid is not itself a constant-small- fluid; its role is a physical comparison of pressure support, not an extension of the preceding approximation through recombination.
Define a normalized smoothing window on physical radius , with comoving radius . If , the smoothed matter density variance is
Thus is the root-mean-square amplitude. For a Gaussian window, ; this choice makes the scale-free integral ultraviolet convergent.
On scales where pressure is negligible, a flat matter-dominated universe has , , and . The density equation has powers satisfying , so
Select the growing mode and exclude residual synchronous gauge modes. Starting with gives . The dimensionless cosmological power spectrum is then
At cosmological horizon crossing, , hence , and
For the smoothed variance, change variables to :
The Gaussian integral equals . At a physical Hubble-radius window, and , giving the Gaussian horizon-crossing variance for a Harrison-Zeldovich spectrum
The window-dependent numerical constant is secondary; cancellation of the scale dependence is the Harrison-Zeldovich spectrum property. This demonstration uses the growing pressureless matter-era idealization and therefore applies to modes crossing in that era. Modes which crossed before equality require radiation-era transfer and cannot be evolved through that period with the matter equation. A pure spectrum with a real-space top-hat window also has a logarithmically divergent ultraviolet variance unless a physical high- cutoff or transfer function is supplied; the Gaussian choice makes the stated idealized variance well defined.
Photon decoupling is a dynamical loss of frequent scattering, distinct from the chemical conversion of charged particles into neutral atoms. Before recombination, Thomson scattering on free Electrons keeps photons coupled to the baryon plasma. The physical scattering rate is
in units with . Define the approximate decoupling temperature by
Equivalently the photon mean scattering time becomes comparable to an expansion time. A refined last-scattering definition uses the optical depth and the maximum of the visibility function; it need not coincide exactly with the local rate criterion. As falls, scattering becomes inefficient, ordinarily after the half-ionized recombination stage, so .
Residual electron freeze-out concerns the remaining ionized fraction. The Saha equation presumes sufficiently rapid forward and reverse reactions to maintain chemical equilibrium. Real recombination is slowed by radiative bottlenecks and expansion; the free-electron fraction can therefore exceed its equilibrium value. Eventually the effective recombination rate per free Electron, roughly , falls below . The remaining Electrons cannot all find and recombine with protons within an expansion time. A small nonzero residual fraction survives while the Saha prediction decreases exponentially toward zero.
The sketch shows this lag and residual floor. Its kinetic curve and the position of the decoupling marker are schematic: the data supplied determine the equilibrium curve, but not a quantitative recombination history or an exact . Determining those needs the expansion history and atomic transition rates. The horizontal axis decreases toward the right to follow cosmological cooling.
Figure 1.
Saha equilibrium and a schematic delayed recombination curve, with recombination and photon-decoupling temperatures marked during cooling
. Details.
Put , the Thomson scattering rate per unit conformal time. The cosmological optical depth is the integrated scattering rate from emission time to observation, so a Poisson scattering process has survival probability
Since , the cosmological visibility function is
A photon scatters in an interval with rate factor and then reaches us without another scattering with probability , so is its probability of last scattering in that interval. For an optically thick initial epoch,
If the initial optical depth is finite, the missing weight represents photons that have already stopped scattering before .
Before cosmological recombination, the optical depth to the present is enormous: , and is also small because almost every scattering is followed by another. Through photon decoupling, rises rapidly and has its dominant positive peak. After decoupling, without reionization, is close to one and reaches exactly one at , while is small because the remaining electron density and scattering rate are small. Residual ionization can give a weak tail; there is no second reionization peak. The peak marks the Cosmic microwave background last-scattering surface.
Figure 1.
The photon visibility function during recombination
. Normalized schematic without reionization. The right curve is the derivative of the left curve and has unit area; the horizontal scale and width are illustrative, not a numerical recombination calculation. Details.
Write for a comoving coordinate, for the peculiar velocity, for the Hubble parameter, and for the density contrast. For a homogeneous self-gravitating fluid, the linearized continuity equation, Euler equations for an inviscid fluid, and cosmological Poisson equation are
Here is the physical adiabatic sound speed, gradients refer to comoving coordinates, and the homogeneous background force has already been subtracted. Differentiating the first equation and eliminating and gives
A Fourier mode with comoving wavenumber therefore obeys
The Jeans wavenumber occurs where the coefficient in parentheses vanishes. Since a comoving wavelength is , the result is
Below this comoving Jeans length, pressure restores a displaced fluid element: there are acoustic waves rather than growing Jeans instability. In a static background they oscillate; expansion changes their amplitude and frequency, and dissipation can damp them. Pressure support alone does not erase them. In the pre-recombination photon-baryon fluid, Silk damping supplies a separate erasure mechanism.
For collisionless matter, replace the sound speed by a characteristic velocity dispersion. For an isotropic Maxwell-Boltzmann velocity distribution, defining , the static collisionless marginal-stability calculation gives
The coefficient depends on the distribution and the convention for velocity dispersion. Small-scale suppression is now collisionless free streaming and phase mixing, rather than collisional acoustic waves. An instantaneous collisionless Jeans length must be distinguished from the accumulated comoving distance traveled since decoupling.
For the requested thermal histories, an Einstein-de Sitter universe describes the matter-era limit: it cannot literally also have a radiation-dominated epoch. Interpret the question as a spatially flat universe with negligible cosmological constant, passing from radiation domination, , to matter domination, . For the usual schematic dynamical Jeans scale use in the estimate, so
Near the particle horizon a relativistic perturbation calculation replaces the Newtonian derivation; the scaling estimate still identifies the relevant sound-crossing or streaming scale.
Before cosmological recombination, adiabatic initial conditions keep the entropy per baryon fixed in the tightly coupled photon-baryon fluid. With , its photon-baryon sound speed is
For example, this follows from and . In the radiation-dominated limit , ; in the strongly baryon-loaded matter-era limit , . Thus the baryon Jeans length across recombination has the asymptotic history
The transition to the plateau is smooth: retaining gives in matter domination. If baryon loading is still small, this intermediate segment instead rises as until loading becomes important.
At photon decoupling, radiation ceases to provide pressure support to the baryons, so the sound speed drops from the coupled-fluid value to that of a nonrelativistic ideal gas, , with . The drop in comoving Jeans length at fixed density is the same ratio of sound speeds, typically several orders of magnitude. Subsequently an isentropic process gives , so and . Residual energy exchange through Compton scattering can initially keep , producing a short approximately flat gas segment before the decreasing segment. The idealized sketch assumes immediate adiabatic cooling:
log(comoving baryon Jeans scale)
  ^
  |                   __________________
  |                 /                   |
  |               /                     |  loss of photon pressure
  |             /                       |
  |           /                         +---\
  |         /                                \
  +-------------------|-----------------|----------> log(time)
                    t_eq              t_rec
        slope 1/2          slope ~0        slope -1/3
For the cold-dark-matter kinetic Jeans scale, take the stated toy history to mean that collisions maintain a particle temperature equal to the photon temperature until kinetic decoupling. This is an assumption about thermal coupling; a negligible collision rate after decoupling makes the subsequent treatment collisionless. Assume the ordering . While relativistic, the characteristic particle speed is constant; while nonrelativistic but thermally coupled, ; after kinetic decoupling, particle momentum redshifts as and . Combining these with the background dynamical density gives
There is no sudden change at cosmological recombination: cold dark matter has already decoupled. Its small initial velocity dispersion makes its suppression scale much smaller than the pre-recombination baryon scale.
log(comoving CDM dynamical Jeans scale)
  ^
  |                        _______________
  |                    ___/               \
  |                ___/                    \       no jump
  |              /                          \      at t_rec
  |            /                             \
  +-----------|------------|--------------|----|------> log(time)
            t_NR         t_dec          t_eq t_rec
   slope 1/2    slope 1/4       slope 0       slope -1/3
The density convention matters. If one instead defines a species' self-gravitating collisionless Jeans length using only its own nonrelativistic density , its thermally coupled segment is constant and its decoupled segment decreases as even during radiation domination. The constant radiation-era segment in the second sketch is the usual instantaneous background dynamical streaming estimate, not an exact single-species self-gravity threshold. A complete multicomponent calculation uses the separate perturbation equations and their combined cosmological Poisson equation.
Photon decoupling occurs when the Thomson interaction rate falls below the expansion rate,
At decoupling the universe is approximately matter dominated, so
Using in gives
For , the Saha equation makes exponentially steep through . Power-law changes of or therefore cause only a small shift in , explaining why is a good approximation.
Residual electron freeze-out Created 2026-10-07 Updated 2026-10-14
When the effective recombination rate per Electron becomes smaller than the expansion rate, the small remaining ionized fraction can no longer track chemical equilibrium. A nonzero Electron fraction remains, exceeding the exponentially small Saha equation prediction. The effective rate includes atomic radiative bottlenecks; chemical freeze-out is distinct from photon decoupling by scattering.
Saha half-ionization temperature Created 2026-10-07 Updated 2026-10-14
With the rounded hydrogen Saha equation coefficient, the condition gives the displayed equation for . Its low-temperature root is about , hence eV for hydrogen binding energy eV. Dropping the slowly varying logarithm gives the rough leading-log estimate eV. This equilibrium convention is distinct from photon decoupling and from a full kinetic recombination calculation.