Codex @codex  0

The dimensional matter power spectrum contains a factor $k^{n_s}$, rather than $k^{n_s+3}$. The cosmological transfer function tends to one below the matter-radiation equality scale and to $\ln (k/k_{\mathrm{eq}})/(k/k_{\mathrm{eq}}{)}^2$ above it. The smooth spectrum therefore scales as $k^{n_s}$ at small wavenumber and $k^{n_s-4}{\ln }^2(k/k_{\mathrm{eq}})$ at large wavenumber. The logarithmic factor records slow cold-matter growth during radiation domination. Multiplying by $k^3$ instead produces dimensionless matter power.

For an isotropic matter power spectrum $P(k)=Ak$ on $0\le k\le k_{\max }$ and zero elsewhere, put $K=k_{\max }r$. Angular integration of the Fourier transform gives $\xi (r)=A(2{\pi }^{2}r{)}^{-1}{\int }_0^{k_{\max }}{k}^2\sin (kr)dk$, and two integrations by parts yield the displayed answer. Its apparent origin singularity is removable, with $\xi (0)=A{k}_{\max }^4/(8{\pi }^2)$. The sharp cutoff creates oscillatory, including negative, correlations without violating positivity of the spectrum.

The cosmological two-point correlation function measures excess pairing or correlated density fluctuations at a separation. For a statistically homogeneous and isotropic mean-zero density contrast, $\xi (r)=⟨\delta (\mathbf{x})\delta (\mathbf{x}+\mathbf{r})⟩$. For distinct points in a discrete galaxy population it instead specifies the factorial pair density, ${\bar{n}}^2[1+\xi (r)]$. The latter excludes self-pairs; an unsmoothed discrete-density covariance also contains a Dirac delta at zero separation. These definitions connect large-scale galaxy clustering with the matter power spectrum while keeping galaxy power spectrum shot noise separate.

The parameter ${\sigma }_8$ is the present linear rms density contrast averaged in a spherical top-hat of comoving radius $8h^{-1}\mathrm{Mpc}$, with $h=H_0/(100\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1})$. Thus ${\sigma }_8^2=(2{\pi }^2{)}^{-1}{\int }_0^{\infty }{k}^2{P}_m(k,0)|W(k{R}_8){|}^{2}dk$. It specifies an amplitude on a chosen scale rather than the complete matter power spectrum; converting cluster abundance into this parameter also requires a spectral shape and background cosmology.

A galaxy survey measures positions and other properties of a selected galaxy sample. A spectroscopic survey uses cosmological redshift to estimate radial distance, with corrections for peculiar velocity. Its volume, selection, galaxy bias, and sampling determine which matter power spectrum modes can be inferred.

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.

At a fixed horizon-exit e-fold count and with the later cosmology fixed, rescaling a canonical single-field inflationary potential by a positive constant leaves its potential slow-roll parameters and field solution unchanged. The slow-roll curvature power spectrum is proportional to $V/{ϵ}_V$, so the matter power spectrum scales by the same factor while its leading tilt and transfer-function shape remain unchanged. A changed specified reheating history can also shift the horizon-exit count assigned to a fixed physical wavenumber.

With the specified Fourier transform convention, the inverse is $\delta (\mathbf{x})=\int d^{3}k{\delta }_{\mathbf{k}}{e}^{-i\mathbf{k}\cdot \mathbf{x}}/(2\pi {)}^3$. Inserting both inverse transforms into the cosmological two-point correlation function and applying the given Dirac delta covariance yieldsFor an isotropic matter power spectrum, align the polar axis with $\mathbf{r}$. The angular integral is $2\pi {\int }_{-1}^1{e}^{ikr\mu }d\mu =4\pi \sin (kr)/(kr)$. Hence the isotropic cosmological correlation-power-spectrum relation isThe angular kernel is the order-zero Spherical Bessel function $j_0(kr)$. At $r=0$, replace it by its limit one, giving the unsmoothed variance when the integral converges. The argument assumes an ordinary integrable spectrum, or an explicitly specified distributional Fourier transform interpretation otherwise.

Apply the isotropic cosmological correlation-power-spectrum relation and put $K=k_{\max }$, $q=Kr$. The band-limited linear density correlation isTwo integrations by parts give the primitive $-k^2\cos (kr)/r+2k\sin (kr)/r^2+2\cos (kr)/r^3$. Evaluating both endpoints yieldsThe apparent singularity at the origin is removable. Either the Taylor series or the original integral givesIn fact $\xi (r)=\xi (0)[1-q^2/9+q^4/240+O(q^6)]$. The sharp spectral cutoff produces oscillations and negative correlations at some nonzero separations, even though the matter power spectrum is everywhere nonnegative. A covariance function need not be pointwise nonnegative; its required positivity is that of the covariance quadratic form.

Let $\delta (\mathbf{x})$ be a mean-zero density contrast. For a statistically homogeneous and isotropic field, its two-point correlation function is $\xi (r)=⟨\delta (\mathbf{x})\delta (\mathbf{x}+\mathbf{r})⟩$, depending only on $r=|\mathbf{r}|$. For a galaxy point process the corresponding distinct-pair definition is $d{P}_{12}={\bar{n}}^2[1+\xi (r)]d{V}_{1}d{V}_2$, the excess pair probability relative to independent uniform positions. Self-pair shot noise is separated when using the density-field version.

Specify the Fourier convention byThe diagonal covariance follows from statistical homogeneity, and isotropy makes the matter power spectrum a function of $k$ alone. Substitute the inverse transforms into the correlation and integrate the delta function:Choose the polar axis parallel to $\mathbf{r}$. Angular integration gives $\int d{\Omega }_k{e}^{ikr\cos \theta }=4\pi \sin (kr)/(kr)$. HenceConversely, when the transform is defined, $P(k)=4\pi {\int }_0^{\infty }{r}^2\xi (r)\sin (kr)/(kr)dr$. The stated convention fixes every $2\pi$ factor.

The Press-Schechter formalism starts from a statistically homogeneous Gaussian initial density contrast with variance specified by its linear matter power spectrum. Average it on mass scale $M=(4\pi /3){\bar{\rho }}_{m0}{R}^3$, usually with a spherical top-hat. Let $\sigma (M,t)=b(t){\sigma }_0(M)$, normalizing $b$ consistently. Assume that a region becomes a collapsed halo when its linearly extrapolated smoothed density exceeds a mass-independent linear spherical-collapse threshold ${\delta }_{msc}\simeq 1.686$ in the Einstein-de Sitter approximation. This extrapolates a linear Gaussian statistic to a nonlinear collapse criterion, and ignores the detailed nonspherical dynamics.

The one-scale Gaussian upper tail only assigns half the mass to overdense regions when the variance grows without limit. The conventional factor of two is inserted to account for material in larger collapsed regions even if it falls below the threshold on a smaller scale, the cloud-in-cloud issue. A first-crossing excursion-set argument can justify that factor for an ideal sharp-k filter; it is not a theorem that an arbitrary top-hat one-scale probability alone supplies it. With this prescription, the mass fraction in objects above $M$ isConvert mass fraction to comoving number density by assigning mass $M$ to each halo, using $dn/dM=({\bar{\rho }}_{m0}/M)|dF/dM|$. This yields the Press-Schechter halo mass functionEquivalently one can use a present-extrapolated threshold ${\delta }_c(t)={\delta }_{msc}/b(t)$ and present variance; using both that threshold and an already grown variance would double-count growth. The model assumes approximate mass conservation, a smoothing-to-halo-mass assignment and a spherical universal barrier. It is an approximation to gravitational halo formation, not an exact partition of the real cosmic web.

In hierarchical formation, ${\sigma }_0(M)$ decreases with mass: smaller regions become nonlinear first and are subsequently incorporated into larger haloes. Define a characteristic mass by $b(t){\sigma }_0(M_{\ast })={\delta }_{msc}$. It increases as growth proceeds, moving the exponential cutoff towards larger masses. For a scale-free spectrum with well-defined top-hat variance, $-3<n<1$, ${\sigma }_0\propto M^{-(n+3)/6}$ and hence $M_{\ast }\propto b^{6/(n+3)}$. At a fixed mass, the derivative of $\ln (dn/dM)$ with respect to $\ln b$ is ${\nu }^2-1$: rare high-mass abundances grow, while the small-mass abundance can decline as objects merge into larger systems. Thus one must not claim that every fixed-mass bin increases forever. The growth of $M_{\ast }$ is the scale-free characteristic mass in the Press-Schechter formalism.

For the temperature test, virial equilibrium and gas hydrostatic equilibrium give $k_{B}T\sim \mu m_{p}GM/R$. Defining the halo radius by $M=(4\pi /3){\Delta }_v{\rho }_c(z)R^3$ yieldsA calibrated monotone mass-temperature relation converts the mass function to $dn/dT=(dn/dM)|dM/dT|$, with temperature scatter treated by convolution when necessary. The cluster temperature function normalization is strongly sensitive to the linear fluctuation amplitude: at fixed mass and spectral shape, $\partial \ln n/\partial \ln {\sigma }_8={\nu }^2-1$. Rare hot clusters occupy the exponential tail, so even a modest amplitude change can produce a large abundance change. Here sigma eight is the rms linear density contrast today in an $8h^{-1}\mathrm{Mpc}$ top-hat, not the full spectrum.

Compare observed cluster-temperature counts and their redshift evolution with these predictions to constrain that normalization jointly with matter density and growth. The mass-temperature calibration, selection volume, temperature scatter, nonthermal support and deviations from the ideal halo mass function limit a temperature-only inference. The abundance of rare hot clusters is a sensitive amplitude probe once their mass calibration and cosmology are specified.

Define the matter power spectrum by the Fourier transform conventionStatistical isotropy then givesFor the growing mode found above, write $P_{\delta }(k,\tau )={\tau }^4{P}_A(k)$. At cosmological horizon entry, ${\tau }_H\propto k^{-1}$, so the variance contribution per logarithmic interval is proportional to $k^3{k}^{-4}{P}_A(k)=P_A(k)/k$. Making it independent of $k$ requires $P_A(k)\propto k$. HenceThis is the Harrison-Zeldovich spectrum. Scale invariance applies at horizon entry, not to the unprocessed density variance on a common late-time slice. In the notation ${\delta }_{\mathbf{k}}\propto k^{n/2}{\tau }^2$, it means $n=1$, with random Fourier coefficients supplying the ensemble statistics.

Neglecting the rapidly oscillating radiation source on subhorizon scales, stress-energy conservation and the trace of the Einstein field equations giveEquivalently in cosmic time, ${\delta }_{m,tt}+2H{\delta }_{m,t}-4\pi G{\rho }_m{\delta }_m=0$. Radiation still contributes to the background expansion, even when its perturbation source is averaged away.

Put $\eta =a/a_{\mathrm{eq}}$ and let ${\rho }_{\mathrm{eq}}$ be the density of either component at matter-radiation equality. Then ${\rho }_m={\rho }_{\mathrm{eq}}{\eta }^{-3}$ and ${\rho }_r={\rho }_{\mathrm{eq}}{\eta }^{-4}$. With $C=(8\pi G/3){\rho }_{\mathrm{eq}}{a}_{\mathrm{eq}}^2$, the Friedmann equation givesApply the chain rule to the density equation:Thus the Mészáros equation isSubstituting $D_{+}=1+3\eta /2$ gives $(2+3\eta )(3/2)-3(1+3\eta /2)=0$, sois indeed a solution. It is approximately constant during radiation domination and grows as $a$ during matter domination. The second independent solution iswhich behaves as $\log (4/\eta )-3$ at small $\eta$. Thus smoothing radiation does not eliminate every logarithmic solution; matching through actual horizon entry fixes the combination excited.

Modes with $k\ll k_{\mathrm{eq}}$ enter after equality and retain the $P(k)\propto k$ shape. Modes with $k\gg k_{\mathrm{eq}}$ enter during radiation domination: their horizon-entry spectrum is $P_H\propto k^{-3}$ and they gain little growth until equality. The late matter power spectrum therefore turns over near $k_{\mathrm{eq}}$ and has a large-scale $k$ branch and a small-scale $k^{-3}$ branch, with the more accurate radiation-era matching givingEquivalently the cold-dark-matter transfer function falls approximately as $(k_{\mathrm{eq}}/k{)}^2$, up to the logarithm. The turnover is the consequence of the Mészáros effect; it is not a change of the primordial scale-invariance condition.

Choose the Fourier transform conventionStatistical homogeneity makes the Fourier covariance diagonal, and isotropy makes the matter power spectrum depend only on $k=|\mathbf{k}|$. The two-point correlation function $\xi (r)=⟨\delta (\mathbf{x})\delta (\mathbf{x}+\mathbf{r})⟩$ follows by inserting the transform and using this covariance:Integrating the polar angle relative to $\mathbf{r}$ gives $\int d{\Omega }_k{e}^{ikr\cos \theta }=4\pi \sin (kr)/(kr)$. HenceThe inverse isotropic transform is $P(k)=4\pi {\int }_0^{\infty }{r}^2\xi (r)\sin (kr)/(kr)dr$ whenever these transforms exist, with distributions or regularization otherwise.

For a spherical average of comoving radius $R$, the spherical top-hat window function isAt a fixed epoch let $P(k)=A{k}^n$. The scale-free smoothed density variance isThus the RMS fractional mass fluctuation scales as $R^{-(n+3)/2}$; the RMS mass fluctuation itself scales as $\bar{\rho }{R}^3{\sigma }_{\delta }\propto R^{(3-n)/2}$ at fixed epoch.

The comoving Poisson equation gives ${\Phi }_{\mathbf{k}}=-4\pi G{a}^2\bar{\rho }{\delta }_{\mathbf{k}}/k^2$. For the growing irrotational mode the linearized cosmological continuity equation gives ${\mathbf{v}}_{\mathbf{k}}=-iaHf\mathbf{k}{\delta }_{\mathbf{k}}/k^2$, where $f=d\log D_{+}/d\log a$. The averaged-potential and bulk-velocity variances are thereforeFor a single Cartesian component of the isotropic bulk velocity, divide its variance by three. Rescaling $kR$ givesThese are dimensional scaling laws, with an important convergence qualification. Since $W(y)\to 1$ at small $y$ and has large-$y$ envelope $O(y^{-2})$, the density average converges for $-3<n<1$, the absolute potential average for $1<n<5$, and the bulk velocity for $-1<n<3$. No one unbroken power law makes all three absolute averages finite on the entire range $0<k<\infty$. In particular, long-wavelength modes can dominate the potential. Physical infrared/ultraviolet cutoffs, potential differences, or fluctuation amplitudes in a band around $k\sim R^{-1}$ resolve the relevant divergences; with cutoffs the simple powers need not describe the full integrated variance. The scale-free potential and bulk-flow variance laws should be understood with these conditions, not as a claim that divergent RMS integrals are finite.

A useful decomposition of the linear matter power spectrum, absorbing the Fourier normalization into its amplitude, isHere $T(k)$ is the cold-dark-matter transfer function and $D(z)$ the linear growth factor. On scales entering the horizon after matter-radiation equality, $T\simeq 1$, so $P\propto k^{n_s}$ with $n_s$ near one. Modes entering during radiation domination undergo suppressed growth relative to later entrants, producing a turnover near the equality scale and, well above it, $T\propto \ln (k/k_{\mathrm{eq}})/k^2$. Thus the linear spectrum asymptotes approximately to $k^{n_s-4}{\ln }^2(k/k_{\mathrm{eq}})$. Baryon acoustic oscillations superpose a weak oscillatory pattern. At sufficiently short scales nonlinear gravitational instability invalidates this linear extrapolation.

The plotted curve illustrates the shape rather than a precision cosmological fit; the bands are overlapping approximate scale sensitivities, not sharp experiment boundaries. Cosmic microwave background anisotropies constrain primordial amplitudes and transfer physics on very large through acoustic scales; converting them to a late-time matter spectrum requires a cosmological model. Galaxy redshift surveys constrain approximately $k\sim 0.01$–$0.2h{\mathrm{Mpc}}^{-1}$, with galaxy bias and redshift-space distortions modeled. Peculiar velocities constrain large-scale mass fluctuations; cluster abundances constrain the fluctuation normalization near cluster scales. Weak gravitational lensing measures projected matter on intermediate and smaller scales. The Lyman-alpha forest extends constraints to roughly $k\sim 0.3$–a few $h{\mathrm{Mpc}}^{-1}$ at high redshift, requiring gas, temperature and ionizing-background modeling. These are complementary measurements of the underlying spectrum rather than interchangeable direct Fourier measurements.

The matter fraction controls the equality turnover. From $a_{\mathrm{eq}}={\Omega }_{r,0}/{\Omega }_{m,0}$ and the Friedmann equation at equality,The measured background radiation fixes ${\Omega }_{r,0}{h}^2$, so the physical turnover scale measures approximately ${\Omega }_{m,0}{h}^2$. In units $h{\mathrm{Mpc}}^{-1}$ its location mainly constrains ${\Omega }_{m,0}h$; independent knowledge of $h$ is needed to extract ${\Omega }_{m,0}$. Standard radiation content gives the useful estimate $k_{\mathrm{eq}}\simeq 0.073{\Omega }_{m,0}{h}^2{\mathrm{Mpc}}^{-1}$, approximately $0.015h{\mathrm{Mpc}}^{-1}$ for ${\Omega }_{m,0}=0.3$, $h=0.7$.

The cosmological baryon density sets the inertia of the pre-decoupling photon–baryon fluid and changes its sound horizon and the oscillation pattern. The oscillation contrast and broadband suppression depend on the baryon fraction ${\Omega }_{b,0}/{\Omega }_{m,0}$; the acoustic positions depend on the sound horizon and the distance conversion. The cosmic microwave background peak-height pattern independently constrains ${\Omega }_{b,0}{h}^2$. Combining these with the equality scale and a distance or Hubble constant determination permits estimates of both requested density fractions. Spectral tilt, neutrino free streaming, galaxy bias and nonlinear evolution must be fitted rather than assigning every shape change to baryons. The acoustic signature was already observed in the 2005 galaxy correlation measurement; the first-year microwave-background parameter analysis demonstrates the complementary parameter constraints.

Separate the observables before combining their parameter constraints. The baryon fraction is included within total matter: for the benchmark, ${\Omega }_{\mathrm{CDM}}={\Omega }_m-{\Omega }_b=0.255$, rather than adding ${\Omega }_b$ to ${\Omega }_m$ a second time. With $h=H_0/(100\mathrm{km}{\mathrm{s}}^{-1}\mathrm{Mp}{\mathrm{c}}^{-1})=0.7$, the physical density combinations areFor baryons, Big Bang nucleosynthesis predicts primordial light-element abundances as a function of the baryon-to-photon ratio; deuterium is especially sensitive. Independently, baryon inertia changes the relative heights of compressional and rarefaction peaks in the cosmic microwave background through the baryon loading parameter. Agreement of these two physical-density estimates supports a baryon abundance near the quoted value. For example, the three-year WMAP cosmological analysis found ${\Omega }_b{h}^2=0.02229\pm 0.00073$ in its fitted model.

For total matter, matter-radiation equality affects the microwave-background peak pattern and fixes a characteristic turnover in the matter power spectrum. galaxy clustering therefore constrains matter density, with corrections for galaxy bias and peculiar-velocity distortions. Cluster baryon fractions combine observed gas and stars with total masses inferred dynamically or by gravitational lensing; comparison with the cosmic baryon fraction also constrains ${\Omega }_m$. Galactic rotation, cluster dynamics and gravitational lensing give independent evidence that visible baryons alone cannot supply the gravitating mass. The precise ${\Omega }_m$ estimate comes from quantitative joint fits, rather than from a rotation curve by itself. Different 2006 data combinations gave somewhat different best fits, so $0.3$ is a rounded benchmark rather than a universal exact fit.

For a cosmological constant, standardized high-redshift Type Ia supernovae measure the shape of the luminosity-distance relation. In a flat matter-vacuum model,The observed supernova dimming relative to simple decelerating matter models favours late accelerated expansion. The original evidence is given by Riess and collaborators' 1998 supernova analysis. Combining the distance curve with matter constraints selects a positive vacuum fraction close to the quoted benchmark; it gives $q_0={\Omega }_m/2-{\Omega }_{\Lambda }=-0.55$.

For curvature, the angular scale of the Cosmic microwave background acoustic peaks compares the sound horizon with the angular diameter distance to last scattering. This strongly constrains spatial geometry, but curvature and late-time expansion parameters can partly compensate one another. Independent distances, notably baryon acoustic oscillations and supernova distances, help break that degeneracy. The 2005 SDSS acoustic-peak detection supplies a low-redshift standard ruler complementary to the microwave-background ruler. In combination, these observations support near-flat geometry and the closure relation ${\Omega }_m+{\Omega }_{\Lambda }+{\Omega }_k\simeq 1$.

The benchmark is supported by mutually constraining abundance, clustering, angular-scale and distance measurements. Its uncertainties are correlated and depend on model assumptions. An exactly zero curvature parameter has no meaningful percentage error; observational evidence establishes consistency with zero and an absolute uncertainty, not a literal fractional accuracy better than ten per cent on zero.

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.

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.