Past exam of the mathematics course of the University of Cambridge / 2006 / iii / Paper 55 (Advanced Cosmology)

Write $\mathcal{H}=a^{\prime }/a$ for the conformal Hubble parameter and ${\theta }_N=i\mathbf{k}\cdot {\mathbf{v}}_N$ for the velocity divergence of a Fourier mode. Work in the synchronous gauge comoving with the cold dark matter, so ${\mathbf{v}}_C=0$. This remains consistent with its pressureless momentum equation ${\mathbf{v}}_C^{\prime }+\mathcal{H}{\mathbf{v}}_C=0$. Its density equation then gives ${\delta }_C^{\prime }=-h^{\prime }/2$.

For the string component, $1+3w_S=0$. Thus it makes no contribution to the trace gravitational source in the supplied perfect-fluid model. Substitute $h^{\prime }=-2{\delta }_C^{\prime }$ and $h^{\prime \prime }=-2{\delta }_C^{\prime \prime }$ in that trace equation to obtainThe absence of a pressure-gradient term reflects the zero cosmological sound speed of the cold matter.

For the $w_S=-1/3$ coasting fluid, the continuity and momentum equations becomeThe second equation follows by taking $i\mathbf{k}\cdot$ the vector equation; in particular, $i\mathbf{k}\cdot i\mathbf{k}=-k^2$, which fixes the pressure sign. Differentiate the first equation and use ${\theta }_S={\delta }_C^{\prime }-3{\delta }_S^{\prime }/2$:Eliminating ${\delta }_C^{\prime \prime }$ with the cold-matter equation givesThese are equations for the specified idealized barotropic string fluid. A cosmic string network with additional anisotropic stress need not obey the same pressure response.

Separate background stress-energy conservation gives ${\rho }_C^{\prime }=-3\mathcal{H}{\rho }_C$ and ${\rho }_S^{\prime }=-2\mathcal{H}{\rho }_S$. ThereforeFor consistency, compute the expansion derivative from the Friedmann equation, rather than assume the printed second hint. Differentiating ${\mathcal{H}}^2=(8\pi G/3)a^2{\rho }_{\mathrm{tot}}$ and using total stress-energy conservation yieldsThe PDF's hint has ${\rho }_{\mathrm{tot}}+P_{\mathrm{tot}}$ where ${\rho }_{\mathrm{tot}}+3P_{\mathrm{tot}}$ is required. It is incompatible with the first hint and conservation, and would not give the requested equation. The expression above is the corrected conformal-time acceleration identity.

For a function $D(\eta )={\delta }_C$, the chain rule givesSubstituting in the cold-matter equation and dividing by ${\mathcal{H}}^2{\eta }^2$ producesorThis is cold-matter growth in a matter-coasting-fluid universe; its coefficients depend on the background density ratio, not on $k$.

Multiply the growth equation by $2{\eta }^2(1+\eta )$:Use the Frobenius method, writing $D={\sum }_{n\ge 0}{b}_n{\eta }^{\alpha +n}$ with $b_0\ne 0$. The lowest power gives the indicial equationThus the leading powers are $\alpha =1$ and $\alpha =-3/2$. For $n\ge 1$, equating the coefficient of ${\eta }^{\alpha +n}$ givesFor the growing branch $\alpha =1$, $b_n=-2(n+1)b_{n-1}/(2n+5)$, so $b_1=-4b_0/7$ and $b_2=8b_0/21$. Taking $b_0=A_k$,Since $\eta \propto a$, the leading behavior is precisely the usual growing density contrast ${\delta }_C\propto a\propto t^{2/3}$ during matter domination. The negative first correction shows the beginning of suppressed growth as the string fraction increases. The other branch starts as $a^{-3/2}$, the familiar decaying matter mode; its exact form is proportional to $\sqrt{1+\eta }/{\eta }^{3/2}$.

The Jeans length is the physical wavelength at which the restoring pressure term and self-gravity balance. In the elementary static-fluid approximation,Wavelengths longer than ${\lambda }_J$ are unstable to self-gravitating collapse, whereas shorter ones have restoring sound oscillations. Expansion and relativistic effects modify this elementary criterion. Ideal cold dark matter has $c_s=0$, hence no positive pressure-supported Jeans length in these equations. For the specified negative-pressure barotropic string fluid, $c_s^2=-1/3$; there is no real restoring-pressure Jeans threshold.

For $\eta \gg 1$, the matter growth equation has leading termsThe gravitational term is one inverse power smaller than the derivative coefficients. Dropping it at leading order gives $({\eta }^2{D}_{\eta }{)}_{\eta }=0$, hence a constant mode and a decaying ${\eta }^{-1}$ mode. A more controlled check is possible: direct substitution gives $D_{-}=\sqrt{1+\eta }/{\eta }^{3/2}$ as an exact solution. Reduction of order then gives the growing mode, normalized by $D_{+}/\eta \to 1$ at early times,The integral grows as $\eta +O(\log \eta )$, while the prefactor is ${\eta }^{-1}[1+O({\eta }^{-1})]$. Thuswith corrections of order $\log \eta /\eta$. The perturbation freezes because the matter fraction tends to zero while the coasting fluid controls the expansion.

To describe the requested cold-dark-matter transfer function, its normalization must be specified. The early matter-era Hamiltonian constraint, in the negative-expansion convention for $K$, gives $-4k^2\Psi /a^2-4H\kappa =16\pi G\delta \rho$. For the growing comoving dust mode, the continuity equation gives $\kappa ={\dot{\delta }}_C=H{\delta }_C$ in proper time. Since $16\pi G\bar{\rho }=6H^2$, this yields ${\delta }_C=-2k^2\Psi /[5(aH{)}^2]$. On early superhorizon scales $\Psi$ is the conserved curvature amplitude, up to the chosen overall sign, and $(aH{)}^2\propto a^{-1}$. This establishes the $k^{2}a$ initial growing behavior rather than assuming identical horizon-entry amplitudes at all epochs. For regular adiabatic modes outside the horizon during early matter domination, $A_k=\mathcal{A}{k}^2{\mathcal{R}}_{\ast }(k)$, with a $k$-independent constant $\mathcal{A}$ and primordial curvature amplitude ${\mathcal{R}}_{\ast }$. The matter equation above has no $k$, so at fixed late timeConsequently the density-per-curvature transfer is proportional to $k^2{D}_{+}$ on scales entering after equality. The conventional shape transfer, which factors out $k^2$ and the common growth factor, is a plateau, $T(k)\simeq 1$ on these scales. Comparing with a matter-only universe that continues growing as $\eta$ instead gives a common amplitude suppression $D_{+}/\eta \simeq 5/(2\eta )$, still independent of $k$. A decline caused solely by later horizon entry would contradict the supplied scale-independent matter equation: the regular synchronous growing solution already evolves outside the Hubble radius. This is the scale-independent matter suppression by a coasting fluid.

There is also a limiting entry scale. Since $a^2{\rho }_S$ is constant, define $k_{\ast }^2=(8\pi G/3)a^2{\rho }_S$. ThenA mode enters the Hubble radius only if $k>k_{\ast }$; $k\gg k_{\ast }$ enters during matter domination, whereas $k\downarrow k_{\ast }$ enters arbitrarily late. Modes $k\le k_{\ast }$ never enter in this ideal coasting future. The requested post-equality entrants occupy $k_{\ast }<k<k_{\mathrm{eq}}$. The following original diagram distinguishes the two transfer conventions and marks this restricted interval; it does not attempt to model modes entering during radiation domination.

The adiabatic initial conditions for separately conserved components have equal ${\delta }_N/(1+w_N)$. Therefore the radiation relation is replaced here byTo check dynamically, put $S={\delta }_S-2{\delta }_C/3$. At $k=0$, subtracting two-thirds of the matter equation from the string equation gives $S^{\prime \prime }+2\mathcal{H}{S}^{\prime }=0$. Selecting $S=0$ and no independent relative-velocity mode preserves the relation. Thus the growing adiabatic string density contrast follows the matter mode and is also approximately frozen at late times to leading order in gradients. This is the adiabatic matter-string density relation.

At late times $\mathcal{H}\to k_{\ast }$, while ${\Omega }_C$ and ${\delta }_C^{\prime }$ vanish. The homogeneous string equation becomesFor $k\gg k_{\ast }$, one solution grows rapidly with exponent approximately $k/\sqrt{3}-k_{\ast }$, instead of undergoing radiation-like acoustic oscillations. This is a gradient instability of a negative-pressure perfect fluid. For $k\ll k_{\ast }$, the growing exponent is only $k^2/(6k_{\ast })$, while the other exponent is approximately $-2k_{\ast }$. The leading zero-gradient approximation is frozen, but at fixed nonzero superhorizon $k$ a slow instability can accumulate over sufficiently long conformal time. The strict gradient limit and infinite-future limit therefore need not commute.

These are the apparent behaviors of the equations given. They do not establish rapid collapse of a realistic cosmic string network: tension, anisotropic stress and the constitutive sound response were excluded in replacing the network by this barotropic perfect fluid. Once the modeled string perturbation becomes nonlinear, the linear calculation itself ceases to apply.

The comoving observer's four-velocity has $u_0=-a$, so its measured photon energy is $E=-u_{\mu }{p}^{\mu }=a{p}^0$, and the comoving momentum is $q=aE=a^2{p}^0$. Let $P^i=p^i/p^0=d{x}^i/d\tau$. The null geodesic condition isDivide the temporal geodesic equation by $p^0=d\tau /d\lambda$, and use the supplied Christoffel symbols:Differentiating $q=a^2{p}^0$ cancels the homogeneous redshift terms. In the remaining first-order term, replace $P^i$ by the unperturbed unit propagation direction $n^i$:This is the synchronous photon momentum redshift. Background physical momentum decreases as $a^{-1}$, while the background comoving momentum stays constant.

The spatial geodesic equation gives $d{p}^i/d\tau =-2\mathcal{H}{p}^i$ in the unperturbed universe. Since the temporal equation gives the identical factor for $p^0$, their ratio $P^i$ is constant at zeroth order. More explicitly, at first order,Every displayed term contains a metric perturbation or its derivative. The orthonormal-frame unit direction differs from $P^i$ by a first-order metric correction, $n^i=P^i+h^i{}_j{P}^j/2+O(h^2)$, whose derivative is also first order. Hence $d{n}^i/d\tau =O(h)$, with no background change of direction. Direction corrections multiplied by an already first-order source contribute only at second order.

Express the Collisionless Boltzmann equation in variables $(\mathbf{x},q,\mathbf{n},\tau )$ and divide the derivative along the photon ray by $d\tau /d\lambda$. Since $f_0(q)$ has no explicit position, direction or conformal-time dependence,This is the chain rule in momentum magnitude and direction. At linear order, $P^i=n^i+O(h)$, while the last two terms are products of first-order perturbations and can be dropped. The synchronous photon momentum redshift from part i then gives, for spatial Fourier convention $e^{i\mathbf{k}\cdot \mathbf{x}}$,The factor $q$ is essential. It is missing in the PDF's displayed intermediate equation, but follows directly from part i and is needed to obtain its final brightness equation.

Let $I_0={\int }_0^{\infty }{q}^3{f}_0(q)dq$. Absorbing the common phase-space normalization into $f$, the homogeneous photon energy density is ${\rho }_{\gamma }=4\pi I_0/a^4$, so the photon brightness perturbation isA small directional temperature change $\Theta =\Delta T/T$ in a Planck distribution gives $f_1=-q{f}_0^{\prime }(q)\Theta$. Integration by parts yieldsbecause the Planck distribution behaves as $q^{-1}$ at low momentum and decays exponentially at high momentum. Thus $\Delta =4\Theta$. Equivalently, frequency-integrated blackbody radiation energy scales as $T^4$. Identifying brightness with a single temperature perturbation assumes this blackbody spectral form; the brightness integral itself is defined more generally.

Finally, multiply the corrected transport equation by $q^3$, integrate and divide by $I_0$. The spatial streaming factor is independent of $q$, while its gravitational source integrates to $(1/2)(-4)h_{ij}^{\prime }{n}^i{n}^j$. ThereforeThis derives the synchronous photon brightness equation with the normalization and sign fixed, rather than inferring it from the dimensionally inconsistent intermediate formula.

There is first an angular-scale typo in the PDF: smaller angles correspond to larger multipoles, approximately $\theta \sim \pi /\ell$. Small scales relative to the first acoustic peak therefore have $\ell \gtrsim 200$, not the printed inequality in the other direction.

The photon-baryon acoustic oscillator produces rapidly varying density and velocity sources around and below the sound horizon. Treating the emission density as a simple long-wavelength gravitational source misses these acoustic peaks and their phases. Silk damping suppresses short-scale perturbations through diffusion before decoupling. The cosmological visibility function has finite width, so different emission times and acoustic phases are averaged rather than sampled on one sharp surface. As tight coupling breaks down, the photon quadrupole and higher moments are generated; Thomson scattering of the quadrupole produces polarization. These require the full photon Boltzmann hierarchy, rather than exactly vanishing emission moments above the dipole. They are the angular-scale limits of instantaneous photon decoupling.

The distinction is important: acoustic evolution and diffusion can already be included in correctly calculated emission ${\delta }_{\gamma }$ and $\mathbf{v}$. They do not themselves invalidate collisionless propagation after decoupling. The two-moment, single-time initial condition and omission of scattering after that surface are the approximations used in the requested line-of-sight formula.

For the stated assumptions, multiply the synchronous photon brightness equation by $e^{ik\mu \tau }$ and integrate:Here ${\tau }_d$ is the decoupling time. Inverse Fourier transformation shifts each source to the unperturbed rayInsert the initial brightness monopole and Doppler dipole, and divide by four:This is the synchronous Sachs-Wolfe line-of-sight formula. The first term is intrinsic emission temperature, the second is the Doppler effect, and the last is the accumulated gravitational redshift in this gauge. The derivative $h_{ij}^{\prime }$ is the partial conformal-time derivative evaluated on the ray, not its total derivative along the ray. Its integral contains endpoint gravitational effects as well as evolving-potential effects when rewritten in Newtonian gauge; it should not be identified solely with the Integrated Sachs-Wolfe effect.

The direction $\mathbf{n}$ is photon propagation toward the observer. The observer-to-source sky direction is its negative, so changing that convention also changes the displayed Doppler sign and the ray parametrization. Correct emitter positions and this direction convention are needed to interpret the formula.

Let ${\gamma }_{ij}$ denote the spatial induced metric and $\gamma =\det {\gamma }_{ij}$. Its conformally rescaled version ${\tilde{\gamma }}_{ij}={\gamma }^{-1/3}{\gamma }_{ij}$ has determinant one. The Jacobi determinant derivative formula consequently givesThus the determinant measures the volume-changing part of the spatial metric, while the unit-determinant part contains shape changes. Metric compatibility for the spatial covariant derivative also gives ${\gamma }^{ij}{D}_i{N}_j=D_i{N}^i$.

There is a sign inconsistency in the PDF. Contracting its defining expression for $K_{ij}$, with its threading convention $d{x}^i-N^{i}dt$ and its future normal $n^{\mu }=N^{-1}(1,N^i)$, givesnot the displayed trace with a minus sign before the divergence. This is the trace of extrinsic curvature with a negative shift. For a concrete check, take ${\gamma }_{ij}={\delta }_{ij}$, $N=1$ and $N^i=b{x}^i$. The defining formula gives $K_{ij}=-b{\delta }_{ij}$ and $K=-3b$, whereas the printed trace would give $+3b$. Reversing the shift convention can produce the other trace formula, but also changes the shift terms in $K_{ij}$ and the normal; the two conventions cannot be combined.

For the requested expansion interpretation, set the shift vector to zero. The sign discrepancy then disappears. A small coordinate volume has proper volume proportional to $\sqrt{\gamma }$, and proper time along its normal is $ds=Ndt$. Its expansion rate isThereforeThis is the local volume Hubble parameter, the mean of the three local directional expansion rates. In a homogeneous isotropic background $\gamma =a^6$, it reduces to $\dot{a}/(Na)$. The interpretation does not require the local expansion to be isotropic: the trace-free part of the extrinsic curvature can still describe shear.

Use $T={\xi }^0$ and ${\xi }^i={\partial }^i\lambda$, and define the physical background Hubble parameter by $H=\dot{a}/(\bar{N}a)$. At first order, the spacetime metric components are $\delta g_{00}=-2{\bar{N}}^2\Phi$, $\delta g_{0i}=a^2{B}_{,i}$ and $\delta g_{ij}=-2a^2(\Psi {\delta }_{ij}+E_{,ij})$. Apply the given passive cosmological gauge transformation component by component. The temporal component changes by $2\bar{N}\dot{\bar{N}}T+2{\bar{N}}^2\dot{T}$; the mixed component changes by ${\bar{N}}^2{T}_{,i}-a^2{\dot{\lambda }}_{,i}$; and the spatial component changes by $-2a\dot{a}T{\delta }_{ij}-2a^2{\lambda }_{,ij}$. HenceThese are the scalar gauge transformations with a background lapse in the paper's sign convention.

Starting in synchronous gauge, preserving both $\Phi =0$ and $B=0$ requiresThe first equation integrates to $\bar{N}T=C(\mathbf{x})$; substituting in the second and integrating givesThe arbitrary lower integration limit is absorbed into $D$. Thus specifying synchronous lapse and shift does not fix the origins of the freely falling clocks or all spatial coordinate labels. This is the residual gauge freedom in synchronous gauge.

During cosmic inflation, choosing these clocks remains arbitrary. For example, the $C$ mode shifts the spatial potential by $HC$ and the density perturbation by $3HC(\bar{\rho }+\bar{P})$. In the de Sitter approximation, $H$ varies little, so a nearly constant metric perturbation can contain a pure coordinate contribution. An inflaton fluctuation or density contrast in synchronous gauge therefore cannot be declared physical solely from its time behavior. Use a uniform-density curvature perturbation or an appropriate field-based gauge-invariant variable, or fix a physical initial slicing. The exact de Sitter case has no evolving homogeneous density clock, which makes density-defined slicing degenerate.

During the standard hot Big Bang, a cold-matter rest frame supplies a convenient additional condition: imposing zero cold-matter velocity removes the nontrivial spatially varying time-shift mode. Before that condition is imposed, the pure-gauge density contrast for constant $w$ is $3H(1+w)C$, which decays as $t^{-1}$ in either radiation or matter domination. In conformal time this is proportional to ${\tau }^{-2}$ during radiation domination and ${\tau }^{-3}$ during matter domination. It can be confused with an independent decaying physical mode unless gauge freedom is fixed. The $D$ mode is a time-independent spatial relabeling and must likewise be fixed by a coordinate or initial-metric convention.

Finally, a scalar density transforms by evaluating the homogeneous scalar at the shifted time:To reach Newtonian gauge from synchronous quantities, impose $\tilde{E}=\tilde{B}=0$. For the nonhomogeneous scalar modes, chooseUsing background stress-energy conservation, $\dot{\bar{\rho }}=-3\bar{N}H(\bar{\rho }+\bar{P})$, gives the synchronous-to-Newtonian density transformation with a lapse:The expression includes the lapse factors because dots refer to the arbitrary time coordinate, not automatically to proper time.

Let $A=\bar{\rho }+\bar{P}$, and assume $A\ne 0$ so that the density-defined curvature is meaningful. Under the cosmological gauge transformation in part ii,Therefore the two time-shift terms cancel:This is the uniform-density curvature perturbation with the paper's overall sign. On a uniform-density slice it is the spatial curvature potential $\Psi$. Its opposite overall sign is also a common convention, so the sign must be kept consistent rather than imported from another metric convention.

To establish time conservation, write $D_t={\bar{N}}^{-1}{\partial }_t$ and $c_a^2=\dot{\bar{P}}/\dot{\bar{\rho }}$. Background stress-energy conservation impliesFor constant $w$ in the given barotropic equation of state, $c_a^2=w$. Differentiate $\zeta =\Psi -\delta \rho /(3A)$:Insert both perturbed evolution equations from the hint, keeping the spatial Laplacian $\Delta ={\nabla }^2/a^2$:Thus the lapse perturbation $\Phi$ and the trace perturbation $\kappa$ cancel exactly. The remaining pressure term is the non-adiabatic pressure perturbation $\delta P_{\mathrm{nad}}=\delta P-c_a^2\delta \rho$. In the assumed barotropic fluid, the perturbed equation of state gives $\delta P=w\delta \rho$ and $c_a^2=w$, so it vanishes. For a Fourier mode,For regular long-wavelength perturbations, the velocity and shear potentials in this combination do not grow as $k^{-2}$, so the right-hand side is gradient-suppressed when $k\ll aH$. Therefore $\zeta$ is conserved to leading superhorizon order, which is the superhorizon conservation of uniform-density curvature requested.

The statement is not exact for arbitrary finite $k$: the displayed gradient term gives its correction. It also requires adiabatic cosmological perturbations; entropy perturbations in a multicomponent system or a nonbarotropic fluid can give $\delta P_{\mathrm{nad}}\ne 0$ and evolve $\zeta$ even at zero gradient order. At exact $w=-1$, $A=0$ and the question's variable is undefined, rather than a well-defined conserved density-slicing perturbation.