Codex @codex  0

Use $\mathcal{H}=a^{\prime }/a$ and the Fourier convention $\nabla \to i\mathbf{k}$. A valid linearization requires all components of $h_{ij}$ to be small: the printed small-determinant condition alone is insufficient, as the matrix $\mathrm{diag}(1,1,0)$ has zero determinant without being small. This is small determinant does not imply a small metric perturbation.

Define the density contrast $\delta =(\rho -\bar{\rho })/\bar{\rho }$ and the scalar velocity potential by $v^i={\partial }^i\theta$, or $v_{\mathbf{k}}^i=i{k}^i{\theta }_{\mathbf{k}}$. Hence $\nabla \cdot \mathbf{v}=-k^2\theta$. To first order, $u^0=a^{-1}$, $u^i=a^{-1}{v}^i$ and $g^{ij}=a^{-2}({\delta }^{ij}-h^{ij})$. The product $u^i{u}^j$ is second order. Substitution into the perfect fluid stress-energy tensor givesHere the equation of state is used locally, so the pressure perturbation is $w\bar{\rho }\delta$.

The contraction in ${\nabla }_{\mu }{T}^{0\mu }=0$ isUsing the displayed connections gives, through first order,In particular ${\Gamma }_{\mu 0}^{\mu }=4\mathcal{H}+h^{\prime }/2$; the spatial connection contraction adds $3\mathcal{H}P+\bar{P}{h}^{\prime }/2$, while differentiating $a^{-2}\rho$ supplies $-2\mathcal{H}\rho$. Subtract the background relation ${\bar{\rho }}^{\prime }+3\mathcal{H}(1+w)\bar{\rho }=0$. The background terms multiplying $\delta$ cancel, leavingFor comoving cold dark matter, $w=0$ and $\theta =0$, so ${\delta }_c^{\prime }=-h^{\prime }/2$ and ${\delta }_c+h/2$ is constant at each spatial point. The proper-volume element of a comoving cell is $a^3\sqrt{\det (I+h)}{d}^{3}x=a^3(1+h/2)d^{3}x$. Conserved mass in the cell then gives the density change opposite to the perturbation of its volume. In comoving synchronous coordinates the linear CDM density change is a volume change, not transport across cell boundaries.

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.

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 numerical fractions describe physically very different components. Use $c=1$ so that pressure and energy density have the same units. The stress-energy conservation equation for a homogeneous component, together with its equation of state $p_i=w_i{\rho }_i$, givesNonrelativistic baryons and cold dark matter have $w\simeq 0$ and dilute as $a^{-3}$; radiation has $w=1/3$ and dilutes as $a^{-4}$; vacuum energy has $w=-1$ and remains constant. The curvature term scales as $a^{-2}$ in the Friedmann equation but is a geometric contribution, not an additional material fluid.

There is consequently a coincidence in the present mixture. Normalize $a_0=1$. ThenAt early times vacuum is negligible, while at late times it dominates overwhelmingly. Observing order-unity matter and vacuum fractions during their comparatively recent crossover motivates the cosmic coincidence problem: why should this crossover occur near our observing epoch?

The absolute scale gives a second puzzle. Restoring SI units, $H_0\simeq 2.27\times 10^{-18}{\mathrm{s}}^{-1}$ yieldsIts energy density is ${\rho }_{\Lambda ,0}{c}^2\simeq 5.8\times 10^{-10}\mathrm{J}{\mathrm{m}}^{-3}$. Quantum vacuum contributions associated with ordinary particle-physics scales can be far larger than this observed value. Explaining the tiny net value, including cancellations among contributions, is the cosmological constant problem. The required negative pressure also differs qualitatively from ordinary dilute matter.

The matter budget supplies another explanatory gap: only $0.045$ is baryonic, while $0.255$ must be nonbaryonic matter in this benchmark. Together with $0.7$ in vacuum energy, most of the cosmic density is in components whose microscopic explanation is not supplied by fitting their fractions. Finally, near-zero curvature raises the flatness problem. For nonzero $k$,During radiation domination $H\propto a^{-2}$, so $|{\Omega }_k|\propto a^2$; during matter domination $H\propto a^{-3/2}$, so $|{\Omega }_k|\propto a$. Near-flatness today therefore requires extremely small early curvature in an unmodified decelerating history. Cosmic inflation offers a mechanism that reduces this curvature fraction through accelerated expansion.

The apparent contrivances are the recent matter-vacuum coincidence, the tiny vacuum-energy scale, the near-flat initial geometry and the unexplained microscopic dark components. Their observationally inferred values impose concrete requirements on a deeper physical account.

In the nonrelativistic dilute regime, $E(p)=m{c}^2+p^2/(2m)+\cdots$ and the Fermi-Dirac distribution reduces to the Maxwell-Boltzmann distribution. HenceDifferentiating the Gaussian integral with respect to its coefficient gives ${\int }_0^{\infty }{p}^2{e}^{-a{p}^2}dp=\sqrt{\pi }/(4a^{3/2})$. ThereforeThe degeneracy factor here is the one in the supplied integral, namely one. Substituting $kT=m{c}^2/\alpha$, $\mu =0$ and dividing by the supplied photon density givesAt $\alpha =20$ this is about $4.8\times 10^{-8}$. The relic is already nonrelativistic at decoupling and becomes still colder as its momenta redshift, making it a possible cold dark matter component if it is stable and interacts sufficiently weakly with light.

For orientation, write ${\eta }_b=n_b/n_{\gamma }$. Ignoring subsequent entropy transfer, its density relative to baryons would be ${\rho }_X/{\rho }_b=(1/20)(n/n_{\gamma })/{\eta }_b$. With the cosmological baryon-to-photon ratio of order $6\times 10^{-10}$, this is of order four, demonstrating the right cosmological scale. The measured dark-to-baryonic density ratio is likewise of order five; an observational reference is arxiv.org/abs/1807.06209 . This is a candidate abundance argument, not an exact relic-density calculation: photon heating after decoupling changes $n/n_{\gamma }$, and stability and interaction properties were not specified. The supplied mass and freeze-out parameter alone therefore do not establish a precise fit.

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.

Use units with $c=1$. For the flat FRW metric, the Friedmann equation is $H^2=8\pi G\rho /3+\Lambda /3$, where $H=\dot{a}/a$ is the Hubble parameter. The cosmological continuity equation for a separately conserved component with constant equation of state $P=w\rho$ gives $\rho \propto a^{-3(1+w)}$. Consequently radiation in cosmology scales as $a^{-4}$, pressureless matter as $a^{-3}$, and cosmological constant energy is constant. With $a_0=1$, define the present critical density and cosmological density parameters byThe Hubble constant $H_0$ measures the expansion rate today. Here ${\Omega }_r$ includes photons and relativistic neutrinos, ${\Omega }_m$ includes baryons and cold dark matter, and ${\Omega }_{\Lambda }$ measures the vacuum contribution. Spatial flatness implies ${\Omega }_r+{\Omega }_m+{\Omega }_{\Lambda }=1$. ThusSince conformal time satisfies $d\tau =dt/a$, we have $a^{\prime }=a\dot{a}=a^{2}H$. Multiplication by $a^4$ provesFor an exam-level estimate, $H_0=100h\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$ with $h\simeq 0.7$, ${\Omega }_m\simeq 0.3$, ${\Omega }_{\Lambda }\simeq 0.7$, and the specified radiation in cosmology gives ${\Omega }_r=4.2\times 10^{-5}{h}^{-2}\simeq 8.6\times 10^{-5}$. These are rounded illustrative flat-universe values; precision cosmological density parameters depend on the fitted model. For comparison, the base flat-model fit in Planck 2018 cosmological parameters gives $H_0\simeq 67.4\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1}$ and ${\Omega }_m\simeq 0.315$.

After cosmological horizon crossing, the rapidly oscillating radiation source can be averaged away as specified. Well inside radiation domination, the remaining cold dark matter self-gravity is also subleading because ${\rho }_C/{\rho }_R\ll 1$. To leading order the linear cosmological density perturbation equation is thereforeTwo integrations giveThe reference time makes the logarithm dimensionless. A growing superhorizon mode reaches entry with a nonzero derivative, so matching generally excites $B\ne 0$. For example, an abrupt leading-order match from $\delta \propto {\tau }^2$ gives $B={\tau }_H{\delta }_H^{\prime }\simeq 2{\delta }_H$, whereas its exact numerical coefficient depends on the transition and the radiation forcing.

This logarithmic growth is the Mészáros effect. It is an asymptotic radiation-era result: after dropping the oscillating radiation perturbation, the original equation still contains $4\pi G{a}^2{\rho }_C\delta$. Retaining that term and the changing background leads to the Mészáros equation, which provides the interpolation near matter-radiation equality. It would be incorrect to claim that a pure logarithm solves the complete mixed-background equation exactly.

During matter domination, the radiation source can be neglected irrespective of whether the mode is inside or outside the horizon in the supplied linear cosmological density perturbation equation. The Friedmann equation gives $\mathcal{H}=2/\tau$ and $4\pi G{a}^2{\rho }_C=6/{\tau }^2$, just as in part (a). ThereforeSubstituting ${\tau }^q$ gives $(q-2)(q+3)=0$, so for every Fourier modeThus the matter-era growing and decaying density modes have growing behaviour ${\delta }_{C,\mathbf{k}}\propto {\tau }^2\propto a$ for all $k$. There is no pressure-gradient term introducing a $k$-dependent growth rate for cold dark matter. The coefficients retain their different earlier histories; common growth does not erase the spectral shape. A specially chosen pure decaying solution is of course possible, so the growing-mode statement does not assert that every initial condition grows.

Put $\theta =i\mathbf{k}\cdot \mathbf{v}$. The synchronous perfect-fluid density equation and Cosmological Euler equation in synchronous gauge giveDifferentiate the first equation and use the second and the scalar trace metric equation. Before making a nonrelativistic approximation, the result isFor $w=w_m=c_s^2\ll 1$, discard relative pressure corrections to the background gravitational and damping terms, while retaining $c_s^2{k}^2$ because a large wavenumber can compensate a small sound speed. Matter domination and the flat Friedmann equation give $\frac{3}{2}{\mathcal{H}}^2=4\pi G{\bar{\rho }}_m{a}^2$. Thus the leading nonrelativistic density equation in synchronous gauge isThe approximation sign records the small finite-$w$ terms just displayed; the printed Newtonian-form equation is not an exact relativistic equation at nonzero $w$.

The comoving Jeans wavenumber balances pressure and gravity: $c_s^2{k}_J^2=4\pi G{\bar{\rho }}_m{a}^2$. Defining physical wavelength as $2\pi a/k$, the physical Jeans length isFor wavelengths longer than this scale, gravity can overcome pressure and density modes grow; shorter modes are pressure-supported acoustic oscillations, with expansion modifying both growth and damping. Before cosmological recombination, the tightly coupled photon-baryon fluid has substantial photon pressure, $c_{s,\gamma b}^2=1/[3(1+R_b)]$, where $R_b=3{\bar{\rho }}_b/(4{\bar{\rho }}_{\gamma })$. The baryonic Jeans scale is correspondingly large, so baryons oscillate instead of freely collapsing on small scales. Cold dark matter has a much smaller pressure scale and can supply gravitational potential wells. After photon decoupling, baryons lose photon pressure support and their sound speed falls to its thermal-gas value; the baryon Jeans length across recombination drops, allowing baryons to fall into those wells. The earlier photon-coupled fluid is not itself a constant-small-$w$ fluid; its role is a physical comparison of pressure support, not an extension of the preceding approximation through recombination.

Define a normalized smoothing window on physical radius $R$, with comoving radius $r=R/a$. If $⟨{\delta }_{\mathbf{k}}{\delta }_{{\mathbf{k}}^{\prime }}^{\ast }⟩=(2\pi {)}^3{\delta }_D(\mathbf{k}-{\mathbf{k}}^{\prime })P(k,\tau )$, the smoothed matter density variance isThus ${\sigma }_R$ is the root-mean-square amplitude. For a Gaussian window, $W(u)=e^{-u^2/2}$; this choice makes the scale-free integral ultraviolet convergent.

On scales where pressure is negligible, a flat matter-dominated universe has $a\propto {\tau }^2$, $\mathcal{H}=2/\tau$, and $4\pi G{\bar{\rho }}_m{a}^2=6/{\tau }^2$. The density equation has powers satisfying $s(s-1)+2s-6=0$, soSelect the growing mode and exclude residual synchronous gauge modes. Starting with $P(k,{\tau }_{\mathrm{eq}})=Ak$ gives $P(k,\tau )=Ak(\tau /{\tau }_{\mathrm{eq}}{)}^4$. The dimensionless cosmological power spectrum is thenAt cosmological horizon crossing, $k\sim aH=\mathcal{H}$, hence $k\tau \sim 2$, andFor the smoothed variance, change variables to $u=kr$:The Gaussian integral equals $1/2$. At a physical Hubble-radius window, $R=H^{-1}$ and $r=\tau /2$, giving the Gaussian horizon-crossing variance for a Harrison-Zeldovich spectrumThe window-dependent numerical constant is secondary; cancellation of the scale dependence is the Harrison-Zeldovich spectrum property. This demonstration uses the growing pressureless matter-era idealization and therefore applies to modes crossing in that era. Modes which crossed before equality require radiation-era transfer and cannot be evolved through that period with the matter equation. A pure $P\propto k$ spectrum with a real-space top-hat window also has a logarithmically divergent ultraviolet variance unless a physical high-$k$ cutoff or transfer function is supplied; the Gaussian choice makes the stated idealized variance well defined.

For photon-dominated inertia, $R\gg 1$ and $c_s^2\simeq 1/3$. The baryon expansion and pressure terms are suppressed by $1/R$. Ignoring metric driving as stipulated, the common-velocity Euler equation is ${\theta }^{\prime }=-{\delta }_{\gamma }/4+{\sigma }_{\gamma }$. Differentiating the photon continuity equation givesThis limit is compatible with a matter-dominated background: cold dark matter can dominate the expansion while photons dominate the inertia of the photon-baryon subsystem.

In the first nonzero quadrupole approximation, neglect ${\sigma }_{\gamma }^{\prime }$ relative to its collision damping. The shear equation givesSubstitution yields the shear-only photon diffusion damping equationThe sign of the first-derivative term is positive, so finite photon mean free time damps acoustic energy.

Put $\omega =k/\sqrt{3}$ and $\gamma (\tau )=4k^2{\tau }_c(\tau )/15$. For constant ${\tau }_c$, the characteristic roots are $-\gamma /2\pm i\sqrt{{\omega }^2-{\gamma }^2/4}$. The frequency shift begins at second order in the mean free time. To the requested order the general solution iswhere the error describes the frequency expansion at fixed $k$ and time interval. Keeping the exponential retains the accumulated damping rather than expanding away a possibly important effect.

For a general time-dependent coefficient, a first-order solution with prescribed initial conditions can also be written aswhere ${\delta }_0=A\cos [\omega (\tau -{\tau }_i)]+B\sin [\omega (\tau -{\tau }_i)]$. Substitution verifies that the integral solves the oscillator forced by $-\gamma {\delta }_0^{\prime }$, and its value and first derivative vanish initially. This gives the general first-order result without assuming arbitrary time dependence may be absorbed exactly into an exponential.

To obtain the slowly varying subhorizon form, set ${\delta }_{\gamma }=e^{-G}y$ with $G^{\prime }=\gamma /2$. The exact transformed equation isDropping the quadratic term alone leaves ${\omega }^2-{\gamma }^{\prime }/2$, so slow variation is also needed. When $k\gg \mathcal{H}$, $k{\tau }_c\ll 1$, and the mean free time changes on an expansion timescale, $|{\gamma }^{\prime }|/{\omega }^2\ll 1$. The WKB method solution has local frequency $\Omega =({\omega }^2-{\gamma }^{\prime }/2{)}^{1/2}$ and prefactor ${\Omega }^{-1/2}$. At leading acoustic order both reduce to constants; the small derivative corrections to phase and prefactor can be neglected. This slowly varying shear-damped photon oscillator therefore givesA change in the phase origin is absorbed into $A,B$. At zero mean free time the damping vanishes; as time passes the diffusion length grows and $k_D$ decreases. The coefficient is that of the question's isotropic-scattering, shear-only truncation, not the polarization-inclusive full Silk damping formula.

If ionization remains approximately constant during matter domination, $n_e\propto a^{-3}$, so ${\tau }_c\propto a^2\propto {\tau }^4$. Away from the lower limit this givesNear recombination the rapidly changing ionization fraction and eventual failure of tight coupling require the full hierarchy; the integral formula specifies the damping accumulated during the regime used here.

Let $\delta =(\rho -\bar{\rho })/\bar{\rho }$ and define the cosmological density power spectrum by $⟨{\delta }_{\mathbf{k}}{\delta }_{{\mathbf{k}}^{\prime }}^{\ast }⟩=(2\pi {)}^3{\delta }_{\mathrm{D}}(\mathbf{k}-{\mathbf{k}}^{\prime })P(k)$. The contribution per logarithmic wave-number interval is the dimensionless cosmological power spectrum ${\Delta }^2(k)=k^{3}P(k)/(2{\pi }^2)$. The linear cosmological mass variance on mass scale $M$ iswith $R$ a comoving smoothing radius and $W$ the top-hat filter in Fourier space. On halo scales the observed spectrum is consistent with greater variance at smaller masses. Growth therefore brings smaller objects to the collapse threshold earlier, statistically; larger structures assemble by accretion and dark-matter halo mergers. This statistical growth from smaller bound systems to larger ones is hierarchical galaxy formation. It is an ordering of dark-matter assembly, not a rule that every small visible galaxy must precede every large visible galaxy.

For linear modes the linear growth factor gives $P(k,z)=D^2(z)P(k,0)$. During matter domination $D\propto a$, so $P\propto a^2$ and the characteristic nonlinear mass increases. Once modes become nonlinear, coupling between scales and halo formation change the shape, so the same multiplication cannot be used for the entire late-time spectrum. At low redshift accelerated expansion suppresses linear growth. Galaxy clustering traces the matter spectrum with galaxy bias; it should not be equated directly with an unbiased matter measurement.

The broad linear shape was set by the cosmological transfer function before and around matter-radiation equality, with baryonic acoustic structure also imprinted before recombination. A nearly scale-invariant primordial curvature spectrum has ${\mathcal{P}}_{\mathcal{R}}(k)\propto k^{n_s-1}$, with $n_s$ close to one. After converting curvature perturbations to matter-density perturbations,Thus nearly scale-invariant primordial curvature does not mean a constant density $P(k)$. Modes with $k\ll k_{\mathrm{eq}}$ enter the horizon after equality and have $T\simeq 1$. Modes with $k\gg k_{\mathrm{eq}}$ enter during radiation domination, when radiation controls the expansion and cold-matter perturbations grow only slowly. The cold-dark-matter transfer function behaves approximately as $T\propto \ln (k/k_{\mathrm{eq}})/(k/k_{\mathrm{eq}}{)}^2$ at large $k$. Hence the density spectrum turns over near $k_{\mathrm{eq}}=a_{\mathrm{eq}}{H}_{\mathrm{eq}}/c$:The turnover records the equality horizon, while the late nonlinear excess records gravitational clustering. On galactic scales the effective slope is greater than $-3$, giving the growing small-scale variance needed for hierarchical galaxy formation.

Cold dark matter has negligible primordial thermal velocities and a very short collisionless free streaming length on galactic scales. Warm dark matter has appreciable residual velocities while structure is being seeded; particles stream across small fluctuations and reduce their contrast. Its cosmological transfer function is consequently cut off below a characteristic length, suppressing low-mass halos and delaying their formation. Above that cutoff its assembly can still be hierarchical. The important distinction is the free-streaming scale, rather than the present temperature or an arbitrary particle-mass label.

Linear evolution assumes $|\delta |\ll 1$. When $\delta$ becomes of order unity, overdense regions depart strongly from the Hubble flow and can turn around and collapse. Collisionless dark matter develops multistream motion after trajectories cross; gravitational mixing redistributes energy and produces a bound dark matter halo. The formal infinite-density collapse of an ideal spherical pressureless solution is not the physical endpoint. A roughly virialized halo has $2K+W\simeq 0$ and a characteristic virial velocity $V_{\mathrm{vir}}^2=G{M}_h/r_{\mathrm{vir}}$. Its gas virial temperature is conventionallyIt measures the thermal energy associated with the gravitational potential; the numerical factor depends on the velocity-dispersion convention. It does not imply that the collisionless dark matter has a thermodynamic gas temperature.

The unheaded request about baryon conversion efficiency of a halo is also answered here. Define $f_{\ast }=M_{\mathrm{stars}}/(f_b{M}_h)$ using the cosmic baryon fraction $f_b$. In small halos, shallow potentials let stellar feedback drive outflows or repeatedly heat star-forming gas; supernova energy per stellar mass is roughly fixed while binding energy per gas mass scales as $V_{\mathrm{vir}}^2$. Photoheating during reionization also prevents very small halos from retaining or accreting cool gas. Molecular/atomic cooling thresholds further reduce star formation in the smallest systems. These effects make $f_{\ast }$ fall toward low mass.

Near $M_h\sim 10^{12}{M}_{\odot }$, gas can cool efficiently and the potential is deep enough to retain more of it, while a long-lived hot atmosphere and maintenance heating are less effective than in larger systems. At high mass, higher virial temperature and lower cooling efficiency let a substantial hot atmosphere persist. Active-galactic-nucleus feedback can prevent that atmosphere from supplying cold gas and can expel some gas; the cooling-time bottleneck alone is not an adequate explanation for the low stellar fractions of massive groups and clusters. The peak reflects a competition between gas supply/cooling and feedback, rather than complete conversion of all baryons at a sharply universal mass. Its exact location and height depend on epoch, metallicity, gas history and the stellar population included.

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: