Past exam of the mathematics course of the University of Cambridge / 2017 / iii / Paper 310 / 2

Here primes mean conformal time derivatives, $d\tau =dt/a$, and spatial derivatives use comoving coordinates in Newtonian gauge in cosmology. Set ${\theta }_m=\nabla \cdot {\mathbf{v}}_m$. The linear cosmological Euler equation gives ${\theta }_m^{\prime }=-\mathcal{H}{\theta }_m-{\nabla }^2\Phi$, while the continuity equation gives ${\theta }_m=-{\delta }_m^{\prime }+3{\Phi }^{\prime }$.

Differentiate the continuity equation and substitute:Thus the linear matter density equation in Newtonian gauge isThe common potential in these supplied equations assumes that the relevant scalar anisotropic stress is neglected. No Fourier-sign convention is needed for the derivation in position space.

Normalize the scale factor so that $a_0=1$. With separately conserved pressureless matter and radiation in cosmology, their energy densities scale as $a^{-3}$ and $a^{-4}$. Hence the flat Friedmann equation isEquality of the two densities defines matter-radiation equality, with $a_{\mathrm{eq}}={\Omega }_{r,0}/{\Omega }_{m,0}$. Define the dimensionless variableSubstituting $a=a_{\mathrm{eq}}y$ givesBoth reference density fractions are positive. In the two-component flat model their sum at the reference epoch is one; no dark-energy component is being included in this background calculation.

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.

In the smooth-radiation, subhorizon approximation of part (c), the cosmological Poisson equation gives

From $y^{\prime }=y\mathcal{H}$ and part (b),The expanding branch has $y^{\prime }>0$. Applying the chain rule yieldsSubstitute these expressions and divide by $C(1+y)>0$:This is the Mészáros equation for $y>0$; it keeps radiation in the expansion rate and matter in the slowly varying self-gravity source.

The Mészáros effect suppresses the growth of modes which enter the horizon during radiation domination. Modes with comoving wavenumber $k\gg k_{\mathrm{eq}}$ enter earlier, but gain only a logarithmic growth factor until matter-radiation equality; modes entering during matter domination grow as $a$ after entry. Consequently small-scale amplitudes are reduced relative to an extrapolation of the unsuppressed large-scale cosmological density power spectrum, creating a turnover near the matter-radiation equality scale $k_{\mathrm{eq}}={\mathcal{H}}_{\mathrm{eq}}$ in units $c=1$.

More quantitatively, for adiabatic primordial fluctuations with scalar spectral index $n_s$, the late-time linear spectrum can be written, up to a time-dependent normalization,The cold-dark-matter transfer function approaches a constant at $k\ll k_{\mathrm{eq}}$, whereas logarithmic small-scale matter transfer gives $T(k)\propto (k_{\mathrm{eq}}/k{)}^2\log (k/k_{\mathrm{eq}})$ at large $k$. ThusFor approximately scale-invariant initial curvature fluctuations, $n_s\simeq 1$, the large-$k$ scaling is $k^{-3}{\log }^{2}k$, rather than $k$. This discussion ignores baryonic acoustic structure and nonlinear clustering. The Mészáros equation solution basis includes both a constant and a logarithmic branch during radiation domination; the physical logarithmic growth is selected by matching the earlier radiation forcing, so not every mathematical solution grows logarithmically.