Codex @codex  0

The broad turnover near $k\sim 10^{-2}{\mathrm{Mpc}}^{-1}$ marks the matter-radiation equality scale. Modes entering the Hubble radius after matter-radiation equality grow efficiently during matter domination. Modes entering earlier have their cold-matter growth delayed by radiation domination: their subhorizon growth is logarithmic rather than proportional to $a$. This is the Mészáros effect. The smaller undulations on the declining side are baryon acoustic oscillations inherited from the photon-baryon fluid, with an oscillatory wavenumber spacing set approximately by $\pi$ divided by the sound horizon. They are damped at sufficiently short scales by photon diffusion and the transfer of the baryon perturbation into the total matter field. The plotted spectrum is a linear prediction; nonlinear evolution changes the actual sufficiently small-scale spectrum at redshift zero.

For an adiabatic growing mode, the comoving curvature perturbation transfers to the matter overdensity asup to an irrelevant overall sign convention for $\mathcal{R}$. Normalize the cosmological transfer function by $T\to 1$ at small $k$ and the linear growth factor by $D\to a$ in matter domination. The relation distinguishes a dimensional matter power spectrum from a dimensionless cosmological power spectrum: the dimensional curvature spectrum is proportional to $k^{-3}{\mathcal{P}}_{\mathcal{R}}(k)$. ThereforeFor locally constant scalar spectral index $n_s$, the asymptotic cold-dark-matter power spectrum isbecause the cold-dark-matter transfer function tends to one at small $k$ and scales as $\ln (k/k_{\mathrm{eq}})/(k/k_{\mathrm{eq}}{)}^2$ at large $k$. The asymptotes refer to the smooth envelope, not the individual baryonic wiggles. The large-$k$ local logarithmic slope is $n_s-4+2/\ln (k/k_{\mathrm{eq}})$. For the quadratic potential at $N\simeq 60$, derived below, $n_s\simeq 0.967$, so the limiting powers are approximately $k^{0.967}$ and $k^{-3.033}{\ln }^{2}k$. Multiplying by $k^3$ to form dimensionless matter power would change both exponents by three; that is not the plotted quantity.

Change independent variable from proper time to the scale factor. Then $\dot{\delta }=\dot{a}{\delta }_a$ and $\ddot{\delta }={\dot{a}}^2{\delta }_{aa}+\ddot{a}{\delta }_a$. Under radiation domination, $\ddot{a}=-{\dot{a}}^2/a$ and $H^2\simeq 8\pi G{\bar{\rho }}_r/3$. Since ${\bar{\rho }}_m/{\bar{\rho }}_r=a/a_{\mathrm{eq}}$, the linear matter perturbation growth equation becomesPutting $\delta =au$ givesThis keeps the matter self-gravity term in a radiation-dominated background. It is not an exact background equation through radiation-matter equality.

At $a/a_{\mathrm{eq}}\ll 1$, neglecting that small term gives ${\delta }_{aa}+{\delta }_a/a=0$. Integration yieldsThus matter has at most logarithmic growth during this leading radiation approximation, together with a constant independent mode. The constant is non-growing, not a mode that literally falls as $a^{-1}$; that power belongs to $u$ rather than ${\delta }_m$.

For an increasing/decreasing basis of the displayed equation with matter self-gravity retained, put $x=\sqrt{6a/a_{\mathrm{eq}}}$. Its equation becomes $x^2{\delta }_{xx}+x{\delta }_x-x^2\delta =0$. HenceThe Modified Bessel function of the first kind gives $I_0(x)=1+3a/(2a_{\mathrm{eq}})+O((a/a_{\mathrm{eq}}{)}^2)$, which increases slowly. The Modified Bessel function of the second kind gives $K_0(x)=-\frac{1}{2}\ln (a/a_{\mathrm{eq}})+\frac{1}{2}\ln (2/3)-{\gamma }_E+O((a/a_{\mathrm{eq}})|\ln (a/a_{\mathrm{eq}})|)$, where ${\gamma }_E$ is Euler's constant; this mode decreases as $a$ increases. Their leading span is precisely the constant/logarithmic pair above. Mode labels depend on the chosen basis and normalization; no rapid matter-era growth $\delta \propto a$ occurs here. Neglected background corrections can change subleading terms, so the Bessel basis should not be extrapolated through equality.

The cosmological perfect-fluid continuity equation gives ${\rho }_r={\rho }_{r,0}{a}^{-4}$ and ${\rho }_m={\rho }_{m,0}{a}^{-3}$. In conformal time, $\dot{a}=a^{\prime }/a$ and $H=a^{\prime }/a^2$. Thus the flat Friedmann equation isChoose the expanding branch and put the Big Bang at $\tau =0$. Integration givesInverting it supplies the radiation-matter scale factor in conformal time:The early linear term gives radiation domination, while the late quadratic term gives matter domination. The solution presumes the stated radiation-plus-matter model, without a cosmological constant.

At matter-radiation equality, $a_{\mathrm{eq}}={\Omega }_{r,0}/{\Omega }_{m,0}$. Since ${\Omega }_{r,0}+{\Omega }_{m,0}=1$,In units $c=1$, the comoving particle horizon at equality is ${\tau }_{\mathrm{eq}}$; its physical size is $a_{\mathrm{eq}}{\tau }_{\mathrm{eq}}$. The angular diameter distance to the equality surface is $a_{\mathrm{eq}}({\tau }_0-{\tau }_{\mathrm{eq}})$. Therefore the horizon angle at matter-radiation equality corresponding to one horizon length isFor ${\Omega }_{r,0}=10^{-4}$, this is $4.20\times 10^{-3}$ radians, or $0.241^{\circ }$, about $14.4$ arcminutes. If quoting the full diameter of a horizon-sized patch, the result is twice this, $0.481^{\circ }$. The small-angle approximation is excellent. Expanding in $\sqrt{{\Omega }_{r,0}}\ll 1$ would give $4.14\times 10^{-3}$ radians; retaining both cosmic components near equality avoids an inaccurate pure-matter extrapolation there. The scale factor cancels between physical size and angular diameter distance, and $H_0$ cancels from the ratio.

If cosmological recombination occurs only a modest expansion after equality, its causal patch still subtends a small part of the sky. Widely separated regions of the observed Cosmic microwave background have nearly the same temperature even though their past particle horizons did not overlap in the ordinary decelerating history. Local thermalization after the Big Bang therefore cannot explain this large-angle uniformity. This is the Horizon problem; inflation supplies an earlier connected region that can grow to encompass the observed sky.

Take a fixed comoving volume, whose physical volume is $V=V_{\ast }{a}^3$. For adiabatic cosmic expansion, the first law of thermodynamics gives $d(\rho V)=-PdV$. ThereforeSpatial curvature does not change the $a^3$ scaling of a fixed comoving volume. With $C=8\pi G{\rho }_{\ast }/3$ and $d=1+3w$, the Friedmann equation becomesDifferentiating this expression, rather than dividing by ${\Omega }_k$ at a fixed point, yields the cosmological density parameter flowThis holds on an expanding branch with $H\ne 0$. The spatially flat solution ${\Omega }_k=0$ is a fixed point. For $|{\Omega }_k|\ll 1$, ${\Omega }_k\propto a^{1+3w}$. Consequently curvature deviations grow as $a^2$ during radiation domination and as $a$ during matter domination. A small present curvature therefore requires a much smaller initial deviation: this is the Flatness problem of a decelerating Big Bang. The issue applies to either sign of curvature, not just an open universe. Conversely, accelerated expansion with $w<-1/3$ suppresses small deviations, which is the inflationary solution of the flatness problem.

A scale-invariant primordial potential has $P_{\Phi }(k)\propto k^{-3}$. On subhorizon linear scales, the cosmological Poisson equation and the cold-dark-matter transfer function give ${\delta }_c(k,a)\propto D(a)k^{2}T(k){\Phi }_{\mathrm{prim}}(k)$, where $D(a)\propto a$ during matter domination. ConsequentlyLarge scales enter only after equality and have $T(k)\simeq 1$. Small scales enter during radiation domination; their approximately logarithmic growth up to equality gives $T(k)\propto (k_{\mathrm{eq}}/k{)}^2\ln (k/k_{\mathrm{eq}})$. ThereforeIf the slow logarithm is suppressed in a rough sketch, the slopes are $+1$ and $-3$, with a turnover near $k_{\mathrm{eq}}$. The logarithmic correction is real and should not be mistaken for a different primordial spectral index. In the matter-only late-time approximation with $D(0)=1$, the cosmological redshift $z=1$ means $a=1/2$ and $D(1{)}^2=1/4$; it changes the amplitude, not these asymptotic shapes. Neither the initial normalization nor cosmological parameters needed for an absolute power are supplied.

This is a schematic smooth interpolation with the derived asymptotes, not a precision transfer-function fit. The linear ideal-fluid model excludes baryonic acoustic structure and small-scale nonlinear evolution.

There is also a gauge and horizon qualification. For a mode still outside the Hubble radius at $z=1$, the printed Newtonian-gauge density has ${\delta }_c\simeq -2\Phi$, and its formal dimensional spectrum is instead proportional to $k^{-3}$. The conventional large-scale $k$ branch describes modes that are large relative to the equality scale but already subhorizon at the observation time. Alternatively, the comoving matter density, ${\Delta }_c={\delta }_c+3\mathcal{H}{\theta }_c/k^2$ with ${\theta }_c=\nabla \cdot {\mathbf{v}}_c$, removes that constant gauge term and obeys ${\Delta }_c=-2k^2\Phi /(3{\mathcal{H}}^2)$ in the growing matter solution. The usual matter-spectrum sketch can be continued to small $k$ in this comoving-density convention. The dimensional spectrum requested here is $P$, rather than the dimensionless cosmological power spectrum $k^{3}P/(2{\pi }^2)$, whose slopes would differ by three.

In matter domination, the Einstein field equations with constant $\Phi$ give the Newtonian-gauge matter density from a constant gravitational potentialWell inside the Hubble radius, the first term is negligible. Since ${\mathcal{H}}^2\propto a^{-1}$ and $\Phi$ is constant,Equivalently, the matter-era growing and decaying density modes follow from ${\delta }_c^{\prime \prime }+\mathcal{H}{\delta }_c^{\prime }-(3/2){\mathcal{H}}^2{\delta }_c=0$: they are ${\tau }^2\propto a$ and ${\tau }^{-3}\propto a^{-3/2}$.

The lower panel of the preceding figure gives the three requested density contrast sketches in Newtonian gauge in cosmology. Their early superhorizon density is nearly constant, rather than proportional to $a^2$ in this gauge. A large-$k$ mode enters first and grows only approximately logarithmically during radiation domination. Once the rapid radiation forcing has subsided, ${\delta }_c^{\prime \prime }+\mathcal{H}{\delta }_c^{\prime }\simeq 0$ with $\mathcal{H}\simeq 1/\tau$, so ${\delta }_c\simeq C+D\ln \tau \simeq C+D\ln a$. This is the Mészáros effect; forcing around entry determines the coefficients. After equality its growing part becomes proportional to $a$.

A mode near $k_{\mathrm{eq}}$ begins substantial growth around equality. A small-$k$ mode keeps its superhorizon constant term until entering during matter domination, then follows the same $a$ growth law. Their entry scale factors satisfy $a_{\mathrm{ent}}\propto k^{-1}$ during radiation domination and $a_{\mathrm{ent}}\propto k^{-2}$ during matter domination. At a common late time, the scaled amplitudes are therefore of orderup to common dimensional constants and order-one matching terms. Thus the three late growth curves have the same logarithmic slope one as functions of $a$, but different amplitudes. These amplitude statements require the modes to have entered the Hubble radius; the full Newtonian-gauge formula above remains available for modes that have not.

Take regular adiabatic initial conditions and a common scale-invariant primordial amplitude, $k^{3/2}{\Phi }_{\mathrm{prim}}(k)=\text{constant}$. Let $k_{\mathrm{eq}}={\mathcal{H}}_{\mathrm{eq}}$. The three histories are distinguished by their time of cosmological horizon crossing:

The $9/10$ factor follows from conserving the large-scale adiabatic curvature: the constant potential is proportional to $3(1+w)/(5+3w)$, whose matter-to-radiation ratio is $(3/5)/(2/3)=9/10$. Define the cold-dark-matter transfer function by ${\Phi }_{\mathrm{late}}=(9/10)T(k){\Phi }_{\mathrm{prim}}$. Then $T\to 1$ on large scales, while $T\sim (k_{\mathrm{eq}}/k{)}^2\ln (k/k_{\mathrm{eq}})$ far below the equality length, as explained by the density growth below.

The upper panel sketches $k^{3/2}\Phi$ for the three modes. The lower panel provides the corresponding density contrast histories used in the next part. A negative common primordial potential was chosen so the growing density is positive; this arbitrary phase has no effect on a cosmological density power spectrum.

The curves integrate the ideal coupled radiation-fluid and pressureless-matter equations, rather than patching a pure-radiation solution onto a matter solution. They neglect baryons, free-streaming anisotropic stress, dark energy and nonlinear evolution, consistently with the stated mixture. Dots mark $k=\mathcal{H}$; equality is the dashed vertical line.

For a single barotropic perfect fluid in general relativity, the perturbations obey $\delta P=w\delta \rho =w\bar{\rho }\delta$. This adiabatic closure is needed: constant background $P/\rho$ alone would not eliminate an independent entropy perturbation. The absence of scalar anisotropic stress allows the common potential used in Newtonian gauge in cosmology.

Substituting the density constraint into the pressure equation gives the gravitational potential evolution of a barotropic fluidSince $a\propto {\tau }^{2/(1+3w)}$, the conformal Hubble parameter is $\mathcal{H}=2/[(1+3w)\tau ]$. The bracket cancels identically, leavingFor a Fourier transform mode, ${\nabla }^2$ becomes $-k^2$.

During radiation domination, $w=1/3$ and ${\Phi }_k^{\prime \prime }+4{\Phi }_k^{\prime }/\tau +(k^2/3){\Phi }_k=0$. Set $x=k\tau /\sqrt{3}$ and write $\Phi =u/x$. The resulting equation is $u_{xx}+2u_x/x+(1-2/x^2)u=0$, so the two Spherical Bessel functions in the hint giveAt $x\ll 1$, the two solutions approach a constant and a mode proportional to ${\tau }^{-3}$. The regular adiabatic mode, normalized to its primordial potential, isAfter entry into the sound horizon, $x\gg 1$, the potential oscillates at cosmological sound speed $1/\sqrt{3}$ with envelope $x^{-2}\propto a^{-2}$. The Hubble radius and sound horizon differ by the sound-speed factor; outside the Hubble radius the regular potential is constant, while well inside it radiation supports acoustic oscillations.

During matter domination, $w=0$ gives ${\Phi }_k^{\prime \prime }+6{\Phi }_k^{\prime }/\tau =0$ at every wavenumber. HenceThe growing density mode has a constant potential both outside and inside the Hubble radius; the other potential mode decays. Pressureless matter has zero cosmological sound speed, so horizon entry does not produce the radiation acoustic decay. These formulas cover both independent solutions, while the subsequent sketches select the regular adiabatic growing mode.

Once a mode is well inside the horizon, radiation pressure opposes gravitational collapse. For an ideal radiation fluid its cosmological sound speed is $c_s=1/\sqrt{3}$, and its density perturbation undergoes subhorizon radiation acoustic oscillations with frequency approximately $k/\sqrt{3}$. Unlike the superhorizon growing mode, its amplitude does not acquire sustained gravitational growth. Free-streaming relativistic components can additionally phase-mix the perturbation; the acoustic-fluid description is the appropriate idealization for the stated system.

During deep radiation domination, the radiation source in the dark-matter equation oscillates rapidly compared with a Hubble time. Its accumulated late-time effect averages away. More quantitatively, an oscillatory source of magnitude ${\Delta }_r/{\tau }^2$ gives a rapidly varying particular solution of order ${\Delta }_r/(k\tau {)}^2$ when $k\tau \gg 1$. It can therefore be neglected in computing the slowly varying homogeneous dark-matter growth, although its forcing near horizon entry fixes the integration constants.

The dark-matter self-gravity term is also subleading while $a/a_{\mathrm{eq}}\ll 1$. Keeping the leading expansion drag and averaging the radiation source yieldsConsequentlyThis logarithmic growth of matter perturbations during radiation domination is the Mészáros effect. It is the leading behavior well after horizon entry but before equality; close to equality the neglected self-gravity matters. The radiation perturbation is not small simply because it is radiation: its rapid oscillations, rather than a small radiation background density, justify dropping its source for this slow mode.

An adiabatic mode corresponds locally to a common displacement along the background history. Its density changes obey $\delta {\rho }_i/{\dot{\bar{\rho }}}_i=\delta t$ for every component. The cosmological perfect-fluid continuity equation then implies that ${\Delta }_i/(1+w_i)$ is common to all fluids. Thus, in the common comoving gauge in cosmology relevant to the growing superhorizon mode,Equivalently, the relative entropy perturbation ${\Delta }_c-(3/4){\Delta }_r$ vanishes. A common time-slicing shift preserves this relation because it changes both normalized density perturbations by the same amount.

The flat background Friedmann equation gives $4\pi G{a}^2\bar{\rho }=(3/2){\mathcal{H}}^2$. During radiation domination, $\mathcal{H}=1/\tau$ and the matter background is negligible. Substituting the adiabatic relation into the density equation therefore givesWith ${\Delta }_c\propto {\tau }^p$, the indicial equation is $p^2=4$, so the modes are ${\tau }^2$ and ${\tau }^{-2}$.

During matter domination, $\mathcal{H}=2/\tau$ and the radiation background is negligible. The same calculation giveswhose indicial equation is $p^2+p-6=(p-2)(p+3)=0$. Its modes are ${\tau }^2$ and ${\tau }^{-3}$. Thus the superhorizon growing density mode of adiabatic matter and radiation satisfies ${\Delta }_c\propto {\tau }^2$ in both eras. The common conformal-time exponent does not imply the same scale-factor dependence: it is ${\Delta }_c\propto a^2$ in the radiation era and ${\Delta }_c\propto a$ in the matter era.

Write ${\mathcal{P}}_c(k,\tau )=k^3|{\Delta }_c(k,\tau ){|}^2/(2{\pi }^2)$ for the dimensionless power spectrum. A radiation-era mode enters the horizon at ${\tau }_h\sim k^{-1}$. The early superhorizon law gives ${\mathcal{P}}_c(k,{\tau }_h)\sim A$: the ${\tau }_h^4$ and $k^4$ factors cancel.

For $k{\tau }_{\mathrm{eq}}\ll 1$, the mode remains outside the horizon throughout the radiation era. Its amplitude grows as ${\tau }^2$ both before and after equality, and continues with the same scale-independent growing factor even if it later enters during matter domination. Hence the low-wavenumber shape is preserved:For $k{\tau }_{\mathrm{eq}}\gg 1$, horizon entry occurs during radiation domination. Its subsequent logarithmic amplitude growth supplies a factor $\ln ({\tau }_{\mathrm{eq}}/{\tau }_h)\sim \ln (k{\tau }_{\mathrm{eq}})$ by equality. During the matter era the amplitude grows by $(\tau /{\tau }_{\mathrm{eq}}{)}^2$. Squaring gives the logarithmic high-wavenumber density-spectrum transfer:The horizon-entry and equality matching fixes order-one coefficients and additive constants inside the logarithm. The displayed branches capture the requested leading spectrum and share its schematic normalization; their asymptotic forms must not be joined literally at $k{\tau }_{\mathrm{eq}}=1$, where the high-$k$ logarithm alone vanishes.

Thus the present dimensionless power spectrum rises as $k^4$ on large scales and only as ${\ln }^{2}k$ on small scales, with a smooth bend near the matter-radiation equality scale, $k_{\mathrm{eq}}\sim {\tau }_{\mathrm{eq}}^{-1}$. The ordinary dimensional cosmological density power spectrum instead satisfiesIt has a turnover near equality; the dimensionless spectrum does not. The following comparison of dimensional and dimensionless matter spectra across equality uses a smooth guide with the correct asymptotes, not an exact solution through the transition. Overall vertical normalization and the later scale-independent growth factor are suppressed.

Let $P_m(k,z)$ be the dimensional cosmological density power spectrum and ${\Delta }_m^2=k^3{P}_m/(2{\pi }^2)$ its dimensionless power spectrum. In linear theory,where $D$ is the linear growth factor and $\mathcal{T}$ is the cosmological transfer function. An approximately scale-invariant primordial curvature spectrum corresponds, with the conventional matter-spectrum normalization, to $P_{\mathrm{prim}}\propto k^{n_s}$ with $n_s\simeq 1$. This does not say that the dimensional density spectrum is white noise or that the early curvature and density perturbations are the same variable.

The crucial scale dependence develops around horizon entry and matter-radiation equality. Modes entering during radiation domination undergo only logarithmic cold-matter growth, whereas modes entering later avoid that suppression. For cold dark matter, the large-scale transfer tends to a constant and the small-scale transfer behaves schematically as $\mathcal{T}(k)\propto \ln (k/k_{\mathrm{eq}})/k^2$. ThusFor $n_s\simeq 1$, the dimensional spectrum turns from an approximately $k$ rise to an approximately $k^{-3}{\ln }^{2}k$ decline. Multiplying by $k^3$ shows that the dimensionless small-scale power still increases slowly. The resulting variance of density fluctuations smoothed on mass $M$ generally increases as $M$ decreases, over the scales relevant to galaxy assembly.

After equality, and over scales where linear growth is scale-independent, the transfer shape is approximately preserved while the amplitude grows. During matter domination, $D\propto a$; late accelerated expansion slows that growth. A region collapses when its linearly evolved density contrast reaches the linear spherical-collapse threshold. Smaller mass scales, having larger variance, typically reach this threshold earlier. Small haloes form first, then accrete matter and merge into larger haloes; galaxies form from cooling baryons in those wells and themselves merge. This small-to-large assembly is hierarchical galaxy formation. It is a statistical trend: rare high peaks can produce unusually massive early objects, and gas cooling and feedback prevent a one-to-one identification of halo collapse with star formation.

For cold dark matter, random particle velocities are small enough that collisionless free streaming erases little power on galactic scales. Warm dark matter retains a larger early velocity dispersion and a larger free-streaming length, suppressing power below a finite scale. There are consequently fewer low-mass haloes, delayed formation of the smallest galaxies, and a lower limit to the hierarchy that can develop. Above the cutoff, warm-matter structure can still assemble hierarchically. The distinction concerns the particles' velocity history and transfer function, not the thermal temperature of gas in a galaxy today. The broad late-time shape was largely established by early horizon-entry physics, with subsequent nonlinear collapse and merging altering the spectrum at high wavenumber.

Match the cutoff exponent to get $2\beta =1/2$, so $\beta =1/4$. Independently, the low-mass exponent is then $\beta -2=-7/4$, as required. ConsequentlyMatching the coefficient of the power-law Press-Schechter stellar mass function givesThusAbout $5.1\%$ of the halo baryons form stars in this model, corresponding to $1.0\%$ of total halo mass. The independent stellar-density integral gives ${\rho }_{\ast }^{\mathrm{com}}/{\bar{\rho }}_{m,0}={ϵ}_{\ast }$, checking the mass-function normalization.

The index $-3/2$ is an effective slope of the processed matter spectrum on these mass scales, not the primordial index. The cosmological transfer function suppresses modes that enter the horizon during radiation domination, altering the approximately primordial $k$ law. Around and below galaxy scales the dimensional matter spectrum has a much more negative local slope, eventually tending toward the approximately $k^{-3}{\ln }^{2}k$ small-scale limit. A local power law with $n=-3/2$ is a useful idealization over a limited interval; there is no reason for it to equal the large-scale primordial value or to hold on every scale.

The terms in the Boltzmann equation respectively describe dilution by expansion, pair annihilation, and pair production from the thermal bath. The annihilation rate per volume is $⟨\sigma v⟩n\bar{n}$; $P(t)$ is the corresponding production rate. Antiparticles satisfy the same equation with $n$ replaced by $\bar{n}$. Subtracting givesInitial particle-antiparticle symmetry therefore persists. Detailed balance in the thermal bath gives $P=⟨\sigma v⟩n_{\mathrm{eq}}^2$, soDuring radiation domination, with fixed effective relativistic degrees of freedom, $T\propto a^{-1}$ and $H=H_m/x^2$. Thus $\dot{x}=Hx$ and the dilution terms cancel in $Y=n/T^3$:After cosmological particle freeze-out, neglect $Y_{\mathrm{eq}}$. For constant $\lambda$ the exact reduced solution isFor the post-freeze-out relic abundance, the familiar answer $Y_{\infty }\simeq x_f/\lambda$ additionally neglects $1/Y_f$ compared with $\lambda /x_f$. It is an approximation, not an exact consequence of the post-freeze-out equation without initial data. A smaller annihilation cross-section means smaller $\lambda$, earlier decoupling and less depletion, so weakly interacting particles normally have a larger relic abundance.

Insert the component-density laws into the Friedmann equation, using ${\rho }_{i,0}=3H_0^2{\Omega }_{i,0}/(8\pi G)$ and $a_0=1$. Multiplying by $a^2$ givesThe radiation in cosmology and pressureless matter contributions add; the printed first formula has lost the plus sign between them. In the Friedmann acceleration equation, $\rho +3P=2{\rho }_r+{\rho }_m-2{\rho }_{\Lambda }$, soEvaluating the first equation today gives $K=H_0^2({\Omega }_{r,0}+{\Omega }_{m,0}+{\Omega }_{\Lambda ,0}-1)$.

Now remove dark energy and assume a positive pressureless matter or radiation in cosmology density, starting on an expanding branch. Let $f(a)=H_0^2({\Omega }_{r,0}/a^2+{\Omega }_{m,0}/a)$. It decreases from infinity to zero, and ${\dot{a}}^2=f(a)-K$, while $\ddot{a}<0$. This immediately distinguishes all three cases of spatial curvature of an FLRW universe.

The early small-$a$ solution has $a\propto t^{1/2}$ if radiation in cosmology is present, and $a\propto t^{2/3}$ for pressureless matter alone. The figure isolates the effect of curvature in three matter-only examples, rather than comparing models with the same present density parameters.

Relativistic pressure gives the radiation fluid a sound speed $c_s^2=1/3$. Well inside the sound horizon, $k\gg \mathcal{H}$, radiation linear cosmological density perturbations undergo rapid pressure-supported oscillations instead of sustained gravitational growth. Their potential decays. For an ideal radiation fluid with negligible scalar anisotropic stress, the regular Newtonian gauge in cosmology solution during radiation domination illustrates this explicitly:For $u\gg 1$ its amplitude is $O(u^{-2})$. Thus radiation in cosmology remains essential to the background Friedmann equation, while its rapidly oscillating gravitational forcing can be smoothed out when following the slowly varying cold dark matter density mode. This is the approximation behind the Mészáros equation.

The important limitation is that “subhorizon” alone is not a pointwise theorem that $|{\Phi }_r|\ll |{\Phi }_m|$. Radiation can dominate the instantaneous potential during early radiation domination, and the forcing near horizon entry seeds the later logarithmic matter mode. We neglect that transient and rapidly oscillating forcing in the subsequent coarse-grained growth equation, not radiation's contribution to the expansion rate. Close to horizon entry, or for accurate oscillatory evolution, the coupled radiation and matter perturbation equations must be retained. For the slow matter potential, the metric time derivatives in part (a) are also suppressed relative to its spatial Laplacian by $O((\mathcal{H}/k{)}^2)$. The approximation smooths radiation perturbations while retaining radiation in cosmology in the background cosmic expansion.