Past exam of the mathematics course of the University of Cambridge / 2007 / iii / Paper 64 / 2

Let $\mathcal{H}=a^{\prime }/a$ and ${\theta }_N=i\mathbf{k}\cdot {\mathbf{v}}_N$. Choose the synchronous gauge rest frame of cold dark matter, so ${\theta }_C=0$. Its continuity equation then gives $h^{\prime }=-2{\delta }_C^{\prime }$. Substitution in the trace Einstein equation givesFor radiation in cosmology the continuity and Euler equations are ${\delta }_R^{\prime }+4{\theta }_R/3+2h^{\prime }/3=0$ and ${\theta }_R^{\prime }=k^2{\delta }_R/4$. Differentiate the first, eliminate ${\theta }_R^{\prime }$ and $h^{\prime \prime }$, and obtainradiation in cosmology pressure supplies the acoustic restoring term while its energy density contributes twice as strongly to the trace gravitational source through $1+3w_R=2$.

Deep in radiation domination, $a\propto \tau$, $\mathcal{H}=1/\tau$, ${\Omega }_R\simeq 1$ and ${\Omega }_C\simeq 0$. At long wavelengths $k\tau \ll 1$, the adiabatic initial conditions require ${\delta }_C=3{\delta }_R/4$. The cold-matter equation reduces toA power ${\tau }^p$ gives $p^2-4=0$. The regular growing solution therefore hasThe other adiabatic solution diverges as ${\tau }^{-2}$ and is discarded for the regular primordial growing mode. The relation between the density contrasts fixes their relative entropy to zero, not their subsequent equality inside the acoustic horizon.

One can also solve the leading radiation-era equations explicitly, which proves that the same regular mode develops logarithmic cold-matter growth after horizon entry. Put $x=k\tau /\sqrt{3}$. The first equation gives ${\delta }_R=(x^2{\delta }_{C,xx}+x{\delta }_{C,x})/3$. Substituting in the second givesDefineIts useful derivative is $F^{\prime }=1/x-2\sin x/x^2+2(1-\cos x)/x^3$. Direct differentiation verifies the fourth-order equation and gives the associated radiation in cosmology solution:At zero, $F=x^2/8+O(x^4)$ and ${\delta }_R=K{x}^2/6+O(x^4)$, so this is precisely the regular adiabatic radiation-era cold-dark-matter transfer solution with $K=24A/k^2$. The other indicial behaviors of the fourth-order equation are $1,x,x^{-2}$; they are either nonadiabatic at leading order or singular and do not add another regular growing adiabatic mode.

For $x\gg 1$, the integral is $\gamma +\ln x-\mathrm{Ci}(x)$, where $\gamma$ is the Euler's constant and $\mathrm{Ci}$ the cosine integral. Its oscillatory tail cancels the explicit leading $\sin x/x$ term, givingThus the requested late sound-horizon behavior isAn additive constant in ${\delta }_C$ is absorbed into the reference ${\tau }_{\ast }$; the regular adiabatic solution fixes the acoustic phase, with $D=0$ in the explicit time-origin convention above. A general acoustic mixture has both terms. The physical expansion parameter is $k\tau /\sqrt{3}$; comparisons with $2\pi /k$ describe the same long- and short-wavelength regimes up to fixed numerical factors.

The contrast grows much more slowly after entry during radiation domination: radiation in cosmology undergoes pressure-supported acoustic oscillations and cold matter only the logarithmic Mészáros effect. Modes entering well before matter-radiation equality consequently have suppressed late amplitudes relative to modes that stay outside the horizon until near equality. After matter domination begins, pressureless growing modes can grow as $a$. This scale-dependent history produces the turnover and small-scale suppression of the matter transfer function and affects the initial conditions for large-scale structure of the universe. The early subhorizon limit requires both $k\tau \gg 1$ and $\tau \ll {\tau }_{\mathrm{eq}}$; it is not applicable to a mode that enters only in the matter domination.

Let $w=c_s^2$ be constant and small, and write $\theta =i\mathbf{k}\cdot \mathbf{v}$. The supplied synchronous gauge equations implyDifferentiating the first and eliminating $\theta$ givesFor the dominant single component, the Einstein equation gives $h^{\prime \prime }+\mathcal{H}{h}^{\prime }=-3{\mathcal{H}}^2(1+3w)\delta$. Therefore the exact algebraic intermediate relation isThe requested nonrelativistic density equation in synchronous gauge follows to leading nonrelativistic order: discard pressure corrections to the background, Hubble damping, and metric tensor terms while retaining $c_s^2{k}^2\delta$, which can be comparable to gravity at the Jeans wavenumber. Since $w\ll 1$, that scale lies well inside the Hubble radius, so keeping its enhanced gradient term is consistent. The resulting equation iswhere the Friedmann equation gave $3{\mathcal{H}}^2/2=4\pi G{a}^2{\bar{\rho }}_m$. Now ${\delta }^{\prime }=a\dot{\delta }$ and ${\delta }^{\prime \prime }=a^2(\ddot{\delta }+H\dot{\delta })$, soThis is the specified approximation, not an exact finite-$w$ relativistic identity. The pressure-gradient term stabilizes short wavelengths, while gravity drives the long-wavelength Jeans instability.

Now apply the stipulated equation to $P=K{\rho }^{4/3}$ with $K>0$. In the pressure-negligible matter background,Thus $c_s^2{k}^2/a^2$ is proportional to $t^{-2}$, exactly like the gravitational term. DefineThe comoving Jeans wavenumber and $q$ are constant in this era. The equation becomes the Euler-Cauchy equationSetting $\delta \propto t^p$ gives $p^2+p/3+2(q-1)/3=0$. For distinct real roots the full pair of modes isThese are the Jeans modes of a four-thirds polytropic cosmological fluid. The physical Jeans length isHere ${\lambda }_J\propto a$, so its comoving value is constant. For $\lambda >{\lambda }_J$, $q<1$ and $p_{+}>0$, giving genuine growth; the other solution decays. In the long-wavelength pressure-free limit, $p_{+}=2/3$ and $p_{-}=-1$, the familiar matter-era modes. At $\lambda ={\lambda }_J$, the two powers are $0$ and $-1/3$, so the leading mode is constant rather than growing.

For $1<q<25/24$ both powers are negative: pressure prevents growth, although the modes still decay without oscillation. At $q=25/24$ they coincide and the complete solution isFor $q>25/24$ they become a complex-conjugate pair; a real basis isThese short-wavelength acoustic oscillations have a decaying envelope. The onset of instability is $q=1$, while the threshold for oscillatory time dependence is slightly larger, $q=25/24$; expansion damping explains why these two thresholds are not identical.