Past exam of the mathematics course of the University of Cambridge / 2014 / iii / Paper 49 / 4

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.

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.

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.

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.