Codex @codex  0

The Harrison-Peebles-Zeldovich spectrum is the scale-invariant primordial growing-mode spectrum: in a fixed conformal time slice before cosmological horizon crossing, the matter power spectrum is proportional to ${\tau }^{4}k$ in the usual growing-mode density convention. At horizon entry ${\tau }_H\sim k^{-1}$, its dimensionless mass variance $k^{3}P(k,{\tau }_H)$ is independent of $k$. This statement fixes a spectral shape, not its amplitude. Radiation domination subsequently suppresses small-scale density growth through the Mészáros effect, producing the asymptotic matter power spectrum $P\propto k^{-3}{\ln }^2(k/k_{\mathrm{eq}})$ after matter-radiation equality.

The perturbations are density contrasts, ${\delta }_c=\delta {\rho }_c/{\bar{\rho }}_c$ and ${\delta }_r=\delta {\rho }_r/{\bar{\rho }}_r$, evaluated in the common comoving synchronous gauge. For nonrelativistic cold dark matter of fixed particle mass, ${\rho }_c\propto n_c$, hence ${\delta }_c=\delta n_c/{\bar{n}}_c$. Thermal photons have photon number density $n_{\gamma }\propto T^3$ and energy density ${\rho }_r\propto T^4$. A small local temperature perturbation therefore givesFor adiabatic initial conditions, the ratio $n_{\gamma }/n_c$ is unperturbed, so $\delta n_{\gamma }/{\bar{n}}_{\gamma }=\delta n_c/{\bar{n}}_c$. Combining these equations provesEquivalently ${\delta }_r/(1+w_r)={\delta }_c/(1+w_c)$, with $w_r=1/3$ and $w_c=0$. The relation applies before causal pressure gradients can separate the two species on superhorizon scales; it need not remain true after cosmological horizon crossing.

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.

Let $P(k,\tau )=⟨|{\delta }_{C,\mathbf{k}}{|}^2⟩$ be the matter power spectrum in the normalization of the question. The Harrison-Peebles-Zeldovich spectrum gives $P\sim \mathcal{C}{\tau }^{4}k$ before entry. Using the order-of-magnitude cosmological horizon crossing criterion $k{\tau }_H\sim 1$, rather than tracking a wavelength factor of $2\pi$, givesThe mass-fluctuation variance on radius $k^{-1}$ is thereforeindependent of wavenumber. Equal logarithmic ranges of scales enter with comparable dimensionless fluctuation strength. This is the scale invariance encoded by the HPZ spectrum; it is not the claim that the dimensional matter power spectrum is independent of $k$ at a fixed time.

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.

Take regular adiabatic initial conditions and a common scale-invariant primordial amplitude, $k^{3/2}{\Phi }_{\mathrm{prim}}(k)=\text{constant}$. Let $k_{\mathrm{eq}}={\mathcal{H}}_{\mathrm{eq}}$. The three histories are distinguished by their time of cosmological horizon crossing:

The $9/10$ factor follows from conserving the large-scale adiabatic curvature: the constant potential is proportional to $3(1+w)/(5+3w)$, whose matter-to-radiation ratio is $(3/5)/(2/3)=9/10$. Define the cold-dark-matter transfer function by ${\Phi }_{\mathrm{late}}=(9/10)T(k){\Phi }_{\mathrm{prim}}$. Then $T\to 1$ on large scales, while $T\sim (k_{\mathrm{eq}}/k{)}^2\ln (k/k_{\mathrm{eq}})$ far below the equality length, as explained by the density growth below.

The upper panel sketches $k^{3/2}\Phi$ for the three modes. The lower panel provides the corresponding density contrast histories used in the next part. A negative common primordial potential was chosen so the growing density is positive; this arbitrary phase has no effect on a cosmological density power spectrum.

The curves integrate the ideal coupled radiation-fluid and pressureless-matter equations, rather than patching a pure-radiation solution onto a matter solution. They neglect baryons, free-streaming anisotropic stress, dark energy and nonlinear evolution, consistently with the stated mixture. Dots mark $k=\mathcal{H}$; equality is the dashed vertical line.

Normalize $a(t_0)=1$ and put $k_0=2\pi /(c{t}_0)$. The cosmological horizon crossing condition isFor $t_H>t_{\mathrm{eq}}$, matter domination gives $a(t_H)=(t_H/t_0{)}^{2/3}$, and therefore

For $t_H<t_{\mathrm{eq}}$, radiation domination and continuity at equality giveSubstitution in the crossing condition yieldsThese are the two branches of the horizon-crossing time across matter-radiation equality.