Codex @codex  0

After the decay products dominate the background, radiation domination gives ${\rho }_r\propto a^{-4}$. Neglecting curvature and vacuum terms, the Friedmann equation implies $\dot{a}/a=C/a^2$, so $a\propto (t-t_B{)}^{1/2}$ andThe displayed $1/(2t)$ form uses a cosmic-time origin with $t_B=0$. After a late transition the extrapolated radiation-era origin need not be the original Big Bang; this time shift does not change the growth law.

For the remaining pressureless matter, neglect its self-gravity to leading order in ${\rho }_m/{\rho }_r$. Its density contrast obeys $\ddot{\delta }+(t-t_B{)}^{-1}\dot{\delta }=0$. Using the conserved expansion-weighted derivative of a passive density perturbation givesThus logarithmic growth of matter perturbations during radiation domination replaces the static exponential growth and the matter-era $t^{2/3}$ growth. The expansion persists while its dominant radiation component supplies little sustained clustering in this pressureless test-component approximation. The residual matter can eventually become important again because ${\rho }_m/{\rho }_r\propto a$; the approximation is for the radiation-dominated interval. The equation is not a pressureless growth equation for the relativistic decay products themselves, whose pressure and, when a fluid description applies, radiation acoustic oscillations must be retained.

The continuity and Cosmological Euler equation in synchronous gauge combine particularly simply. Differentiate ${\delta }_C^{\prime }=-\nabla \cdot v_C-h^{\prime }/2$ and add $\mathcal{H}{\delta }_C^{\prime }$:The velocity terms cancel by the Cosmological Euler equation in synchronous gauge, so this step does not require setting the initial velocity exactly to zero.

For subhorizon modes deep in radiation domination, take the rapidly oscillating radiation perturbation to be negligible in the averaged gravitational forcing. Also neglect the subdominant matter self-gravity to leading order in ${\rho }_C/{\rho }_r$. With $a\propto \tau$ and $\mathcal{H}=1/\tau$, the resulting equation isIntegrating twice giveswhere the arbitrary reference time ${\tau }_{\ast }$ makes the logarithm dimensionless and can be absorbed into $B$. This is logarithmic growth of matter perturbations during radiation domination.

The neglect of matter self-gravity is needed for the displayed form to be exact within the leading radiation-background approximation. Uniform radiation alone leaves a nonzero term $4\pi G{a}^2{\rho }_C{\delta }_C$ if the cold-matter density is retained. Relative to ${\mathcal{H}}^2$, its coefficient is $3{\Omega }_C(a)/2$, small deep in the radiation era but not near equality. The Mészáros equation retains that effect in a matter-plus-radiation background; its solution basis tends to a constant and a logarithm at $a/a_{\mathrm{eq}}\ll 1$, recovering the present leading result. It connects smoothly to matter-era growth rather than allowing the logarithmic approximation to be extrapolated indefinitely.

Galaxy-scale modes enter the Hubble radius before matter-radiation equality. The earlier radiation forcing and horizon-entry matching ordinarily generate a nonzero logarithmic coefficient, even though subsequent oscillatory radiation forcing is small. A Fourier mode with entry time ${\tau }_h\sim k^{-1}$ accumulates $\log ({\tau }_{\mathrm{eq}}/{\tau }_h)\sim \log (k{\tau }_{\mathrm{eq}})$ before equality. Across a large scale-factor range this factor need not be small, so treating the density contrast as exactly frozen would lose an important part of the galaxy-scale seed amplitude and its scale dependence. Growth is nevertheless much slower than ${\delta }_C\propto a$ in matter domination. This is the Mészáros effect, and the logarithmic matching contributes to the small-scale matter cold-dark-matter transfer function.

In a flat matter-dominated universe, $a\propto t^{2/3}$, $H=2/(3t)$ and $4\pi G{\bar{\rho }}_c=2/(3t^2)$. Substitution of ${\delta }_k\propto t^{\beta }$ gives $\beta (\beta -1)+(4/3)\beta -2/3=0$, whose roots are $2/3$ and $-1$. HenceDuring radiation domination, $a\propto t^{1/2}$ while ${\bar{\rho }}_c/{\bar{\rho }}_{\mathrm{total}}\propto a\ll 1$. The matter self-gravity term is suppressed relative to $H^2$, leaving approximately ${\ddot{\delta }}_k+t^{-1}{\dot{\delta }}_k=0$, with solutions constant and $\log t$. Thus there is no growing power-law mode comparable to the later $t^{2/3}$ growth. The elementary estimate in part (ii) neglects this slow logarithmic growth of matter perturbations during radiation domination.

The plot uses comoving wavenumber in $\mathrm{Mp}{\mathrm{c}}^{-1}$. Using the inverse comoving particle horizon from part (ii) gives $1/r_{p,\mathrm{eq}}\simeq 9\times 10^{-3}\mathrm{Mp}{\mathrm{c}}^{-1}$. The conventional matter-radiation equality scale is insteadThese are consistent order-one conventions for locating the turnover near $10^{-2}\mathrm{Mp}{\mathrm{c}}^{-1}$. The horizon-entry condition is $k/(aH/c)\sim 1$; equating an entire Fourier wavelength $2\pi /k$ to the particle horizon would attach a different numerical factor.

For adiabatic initial conditions, the asymptotic cold-dark-matter power spectrum has $P(k)\propto k^{n_s}T(k{)}^2$, where the scalar spectral index is $n_s\simeq 1$ and $T$ is the cosmological transfer function. Modes below $k_{\mathrm{eq}}$ enter during matter domination and have little relative suppression, giving approximately $P\propto k$. Modes above $k_{\mathrm{eq}}$ enter during radiation domination, when logarithmic growth of matter perturbations during radiation domination is much slower. Their transfer is approximately $T\propto \ln (k/k_{\mathrm{eq}})/k^2$, giving $P\propto k^{-3}{\ln }^2(k/k_{\mathrm{eq}})$ well above equality. This explains the broad turnover; baryon acoustic oscillations add structure, and nonlinear evolution modifies sufficiently small scales.

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.

Under the stated approximation the perturbation equation isDuring radiation domination, $a\propto t^{1/2}$, so $\dot{\delta }\propto t^{-1}$ andMatter perturbations therefore have only logarithmic growth of matter perturbations during radiation domination.

In a curvature-dominated universe, $a\propto t$, so $\dot{\delta }\propto t^{-2}$ andThe nondecaying mode is constant: curvature domination gives freezing of matter perturbations during curvature domination, rather than continued growth.