Codex @codex  0

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.

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.

Use primes for conformal-time derivatives. The cold-dark-matter Cosmological Euler equation in synchronous gauge is $v_C^{\prime }+\mathcal{H}{v}_C=0$, so multiplying by $a$ gives $(a{v}_C{)}^{\prime }=0$. ThereforeIn an expanding universe $a$ increases, damping any initial peculiar velocity as $a^{-1}$. During radiation domination $a\propto \tau$, so the velocity scales as ${\tau }^{-1}$; during matter domination $a\propto {\tau }^2$, it scales as ${\tau }^{-2}$. This is the nonrelativistic limit of the physical-momentum redshift derived in Q1. It concerns the bulk velocity in the synchronous coordinate frame, not a growing peculiar velocity sourced in another gauge. A cold-matter rest-frame choice also removes the residual synchronous time displacement relevant to that component. The decay statement assumes expanding evolution and the pressureless, linear approximation used by the equation.

The continuity and Cosmological Euler equation in synchronous gauge combine particularly simply. Differentiate ${\delta }_C^{\prime }=-\nabla \cdot v_C-h^{\prime }/2$ and add $\mathcal{H}{\delta }_C^{\prime }$:The velocity terms cancel by the Cosmological Euler equation in synchronous gauge, so this step does not require setting the initial velocity exactly to zero.

For subhorizon modes deep in radiation domination, take the rapidly oscillating radiation perturbation to be negligible in the averaged gravitational forcing. Also neglect the subdominant matter self-gravity to leading order in ${\rho }_C/{\rho }_r$. With $a\propto \tau$ and $\mathcal{H}=1/\tau$, the resulting equation isIntegrating twice giveswhere the arbitrary reference time ${\tau }_{\ast }$ makes the logarithm dimensionless and can be absorbed into $B$. This is logarithmic growth of matter perturbations during radiation domination.

The neglect of matter self-gravity is needed for the displayed form to be exact within the leading radiation-background approximation. Uniform radiation alone leaves a nonzero term $4\pi G{a}^2{\rho }_C{\delta }_C$ if the cold-matter density is retained. Relative to ${\mathcal{H}}^2$, its coefficient is $3{\Omega }_C(a)/2$, small deep in the radiation era but not near equality. The Mészáros equation retains that effect in a matter-plus-radiation background; its solution basis tends to a constant and a logarithm at $a/a_{\mathrm{eq}}\ll 1$, recovering the present leading result. It connects smoothly to matter-era growth rather than allowing the logarithmic approximation to be extrapolated indefinitely.

Galaxy-scale modes enter the Hubble radius before matter-radiation equality. The earlier radiation forcing and horizon-entry matching ordinarily generate a nonzero logarithmic coefficient, even though subsequent oscillatory radiation forcing is small. A Fourier mode with entry time ${\tau }_h\sim k^{-1}$ accumulates $\log ({\tau }_{\mathrm{eq}}/{\tau }_h)\sim \log (k{\tau }_{\mathrm{eq}})$ before equality. Across a large scale-factor range this factor need not be small, so treating the density contrast as exactly frozen would lose an important part of the galaxy-scale seed amplitude and its scale dependence. Growth is nevertheless much slower than ${\delta }_C\propto a$ in matter domination. This is the Mészáros effect, and the logarithmic matching contributes to the small-scale matter cold-dark-matter transfer function.

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.

During matter domination, ${\rho }_C\propto a^{-3}$. The flat Friedmann equation in conformal time implies $a^{\prime 2}\propto a$, and henceWith radiation neglected, the CDM density equation in a matter-radiation universe reduces toFor this Euler-Cauchy equation, ${\delta }_C\propto {\tau }^p$ requires $p(p-1)+2p-6=(p-2)(p+3)=0$. Thus the matter-era growing and decaying density modes are ${\tau }^2$ and ${\tau }^{-3}$, respectively, andThere is no explicit $k$ in this pressureless equation, so the result holds on all scales within the stated CDM-comoving synchronous gauge description. The amplitude can still depend on $k$ through its earlier evolution; the result does not equate density contrasts on different time slicings.

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.

Write $\delta ={\delta }_{C,\mathbf{k}}$ and $\mathcal{H}=a^{\prime }/a$, the conformal Hubble parameter. The flat Friedmann equation givesOutside the horizon, the stated adiabatic cosmological perturbation relation is ${\delta }_R=(4/3)\delta$. It expresses equality of ${\delta }_i/(1+w_i)$ for dust and radiation in the density convention of the supplied equations. During radiation domination, $a\propto \tau$, so $\mathcal{H}=1/\tau$ and ${\rho }_C$ is negligible. Consequently $8\pi G{a}^2{\rho }_R=3/{\tau }^2$, and the linear cosmological density perturbation equation becomesThe trial solution $\delta \propto {\tau }^q$ gives $q(q-1)+q=4$, hence $q=2,-2$.

During matter domination, $a\propto {\tau }^2$, so $\mathcal{H}=2/\tau$. Now the radiation source is negligible and $4\pi G{a}^2{\rho }_C=(3/2){\mathcal{H}}^2=6/{\tau }^2$. The equation becomeswhose exponents satisfy $q(q-1)+2q=6$, or $(q-2)(q+3)=0$. Thus the growing solutions in both eras areThe other modes are ${\tau }^{-2}$ in the radiation limit and ${\tau }^{-3}$ in the matter limit. These are statements in the density and gauge convention of the supplied evolution equation; a density contrast outside the horizon is not a gauge-independent observable by itself.

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.

First consider $k{\tau }_{\mathrm{eq}}<1$. Such a mode remains outside the horizon through matter-radiation equality, grows as ${\tau }^2$ on both sides, and keeps the same matter-era growth after entry. Squaring this growth preserves the shape of the HPZ spectrum:up to a scale-independent matching coefficient.

For $k{\tau }_{\mathrm{eq}}\gg 1$, entry occurs in radiation domination at ${\tau }_H\sim k^{-1}$, with $P_H\sim \mathcal{C}{k}^{-3}$. The density amplitude then grows only logarithmically until equality. ThusDuring matter domination, the matter-era linear growth factor of the density is $(\tau /{\tau }_{\mathrm{eq}}{)}^2$, so its square multiplies the matter power spectrum:The turnover at $k\sim {\tau }_{\mathrm{eq}}^{-1}$ records whether a mode entered during radiation domination. Earlier-entering small-scale modes miss the rapid ${\tau }^2$ growth and receive only the logarithmic Mészáros effect. This is the origin of the logarithmic small-scale matter transfer. The formula is a large-$k$ asymptote: the additive constant in the logarithmic solution and smooth equality matching must be retained near $k{\tau }_{\mathrm{eq}}=1$, where the spectrum does not vanish.

Multiplying by $k^3$ gives the small-scale mass-fluctuation varianceThis gives the requested estimate for modes entering before equality, with the logarithm replaced by a smooth order-one factor in the transition region. On comoving scale $L\sim k^{-1}$, unit variance today requiresFor a numerical illustration faithful to this question's flat, matter-plus-radiation background, take $h=0.7$ and ${\Omega }_r=4.2\times 10^{-5}{h}^{-2}$, hence ${\Omega }_C=1-{\Omega }_r\simeq 1$ rather than the $0.3$ appropriate to a vacuum-dominated model. Part 1(b), now with ${\Omega }_m={\Omega }_C$, gives ${\tau }_{\mathrm{eq}}\simeq 33\mathrm{Mpc}$. In this two-component model the same integral is valid all the way to today, givingWith $L=3\mathrm{Mpc}$, ${\tau }_0/{\tau }_{\mathrm{eq}}\simeq 258$ and $\ln ({\tau }_{\mathrm{eq}}/L)\simeq 2.39$. Therefore $\mathcal{C}\sim [258^4(2.39{)}^2{]}^{-1}\simeq 4\times 10^{-11}$, corresponding to a horizon-entry root-mean-square fluctuation of roughly $6\times 10^{-6}$. A very small primordial amplitude, of order $10^{-10}$ in this schematic normalization, can reach unit small-scale variance by today. Factors from the Fourier convention, smoothing window and equality matching preclude a precise normalization here. The calculation marks the onset of nonlinear structure formation, where extrapolating the linear solution ceases to be valid. For a universe with late vacuum domination, one must replace the matter-only growth law by its actual growth factor before making a numerical inference.

Preserving synchronous gauge means keeping both $\delta g_{00}$ and $\delta g_{0i}$ zero. Applying the metric transformation to the background givesThe first equation requires $(a{\xi }^0{)}^{\prime }=0$, so ${\xi }^0=C(\mathbf{x})/a$. For the scalar displacement ${\xi }^i={\partial }^i\lambda$, the second requires ${\lambda }^{\prime }={\xi }^0$, up to a spatially homogeneous coordinate displacement that does not affect scalar perturbations. Integrating givesThis is the residual gauge freedom in synchronous gauge. The lapse and shift conditions do not uniquely choose a slicing and spatial labeling.

A background density is a scalar, so under this remaining time displacementFor a fluid with constant equation of state, background conservation gives ${\bar{\rho }}^{\prime }/\bar{\rho }=-3(1+w)\mathcal{H}$, henceThe synchronous density gauge mode scales as ${\tau }^{-2}$ in radiation domination and as ${\tau }^{-3}$ in matter domination. It can mimic a changing density perturbation without changing physical observables. A density time dependence alone therefore does not identify a physical growing or decaying mode. Choosing coordinates comoving with cold matter fixes the residual time freedom, and a remaining spatial relabeling fixes the initial metric-trace integration convention. Gauge-invariant curvature perturbations provide another way to specify physical initial conditions.

For consistency, linear metric theory requires all background-orthonormal components of $h_{ij}$ to be small. The PDF's small-determinant condition by itself is insufficient: $h_{ij}=\mathrm{diag}(L,0,0)$ has zero determinant even for large $L$. This is small determinant does not imply a small metric perturbation.

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.

Assume both densities are nonnegative and ${\Omega }_{\mathrm{DE},0}<1$, so a nonzero pressureless matter component exists. Their ratio obeysThus matter domination givesAt fixed $H_0$, the early expansion rate does depend on the present dark energy fraction: increasing it reduces ${\rho }_{M,0}=(1-{\Omega }_{\mathrm{DE},0})3H_0^2/(8\pi G)$. There is no contradiction with negligible early dark energy, since the leading expression is simply $H^2\sim 8\pi G{\rho }_M/3$. At fixed physical matter density, that leading rate is independent of the negligible component.

The endpoint ${\Omega }_{\mathrm{DE},0}=1$ has no matter, so the ratio argument does not apply. Its Hubble parameter is $H=H_0{a}^{-3(1+w)/2}$: it diverges for $-1<w<0$, is constant for $w=-1$, and tends to zero for $w<-1$.

Well before matter-radiation equality, radiation domination makes the integrand approximately constant:At equality, ${\Omega }_{r,0}{a}_{\mathrm{eq}}^{-4}={\Omega }_{m,0}{a}_{\mathrm{eq}}^{-3}$, so $a_{\mathrm{eq}}\simeq 1/3000$ and ${\Omega }_{r,0}\simeq 10^{-4}$. Extrapolating the radiation-only result gives $r_p(a_{\mathrm{eq}})\simeq 133\mathrm{Mpc}$. This overestimates the horizon because matter is no longer negligible at equality and increases the expansion rate.

Retain both matter and radiation, still neglecting the insignificant early cosmological constant. The particle horizon through radiation-matter equality isThereforeThis is a comoving distance; the proper horizon at equality is approximately $110/3000\mathrm{Mpc}\simeq 37$ kiloparsecs.

For today's horizon, extrapolating matter domination givesThis is an order-of-magnitude estimate: dark energy is significant late, and radiation matters early. Numerical integration of the complete expression with the supplied density parameters and ${\Omega }_{r,0}\simeq 10^{-4}$ givesBoth neglected positive-density terms increase $H$ at fixed $a$ and reduce the horizon relative to the matter-only extrapolation.