Matter power spectrum

Matter power spectrum inferred from various cosmological probes.

The matter power spectrum describes the density contrast of the universe (the difference between the local density and the mean density) as a function of scale. It is the Fourier transform of the matter correlation function. On large scales, gravity competes with cosmic expansion, and structures grow according to linear theory. In this regime, the density contrast field is Gaussian, Fourier modes evolve independently, and the power spectrum is sufficient to completely describe the density field. On small scales, gravitational collapse is non-linear, and can only be computed accurately using N-body simulations. Higher-order statistics are necessary to describe the full field at small scales.

Definition

Let ${\displaystyle \delta (\mathbf {x} )}$ represent the matter overdensity, a dimensionless quantity defined as:

$${\displaystyle \delta (\mathbf {x} )={\frac {\rho (\mathbf {x} )-{\bar {\rho }}}{\bar {\rho }}},}$$

where ${\displaystyle {\bar {\rho }}}$ is the average matter density over all space.

The power spectrum is most commonly understood as the Fourier transform of the autocorrelation function, ${\displaystyle \xi }$, mathematically defined as:

$${\displaystyle \xi (r)=\langle \delta (\mathbf {x} )\delta (\mathbf {x} ')\rangle ={\frac {1}{V}}\int d^{3}\mathbf {x} \,\delta (\mathbf {x} )\delta (\mathbf {x} -\mathbf {r} ),}$$

for ${\displaystyle \mathbf {r} =\mathbf {x} -\mathbf {x} '}$. This then determines the easily derived relationship to the power spectrum, ${\displaystyle P(\mathbf {k} )}$, that is ${\displaystyle \xi (r)=\int {\frac {d^{3}k}{(2\pi )^{3}}}P(k)e^{i\mathbf {k} \cdot (\mathbf {x} -\mathbf {x} ')}.}$

Equivalently, letting ${\displaystyle {\tilde {\delta }}(\mathbf {k} )}$ denote the Fourier transform of the overdensity ${\displaystyle \delta (\mathbf {x} )}$, the power spectrum is given by the following average over Fourier space: ^[1]

$${\displaystyle \langle {\tilde {\delta }}(\mathbf {k} ){\tilde {\delta }}^{*}(\mathbf {k} ')\rangle =(2\pi )^{3}P(k)\delta ^{3}(\mathbf {k} -\mathbf {k} ')}$$

(note that ${\displaystyle \delta ^{3}}$ is not an overdensity but the Dirac delta function).

Since ${\displaystyle P(k)}$ has dimensions of (length)³, the power spectrum is also sometimes given in terms of the dimensionless function: ^[1]

$${\displaystyle \Delta ^{2}(k)={\frac {k^{3}P(k)}{2\pi ^{2}}}.}$$

Development according to gravitational expansion

If the autocorrelation function describes the probability of a galaxy at a distance ${\displaystyle r}$ from another galaxy, the matter power spectrum decomposes this probability into characteristic lengths, ${\displaystyle k\approx 2\pi /L}$, and its amplitude describes the degree to which each characteristic length contributes to the total over-probability.

The overall shape of the matter power spectrum is best understood in terms of the linear perturbation theory analysis of the growth of structure, which predicts to first order that the power spectrum grows according to:

$${\displaystyle P(\mathbf {k} ,t)=D_{+}^{2}(t)P(\mathbf {k} ,t_{0})=D_{+}^{2}(t)P_{0}(\mathbf {k} )}$$

Where ${\displaystyle D_{+}(t)}$ is the linear growth factor in the density, that is to first order ${\displaystyle \delta (r,t)=D_{+}(t)\delta _{0}(r)}$, and ${\displaystyle P_{0}(\mathbf {k} )}$ is commonly referred to as the primordial matter power spectrum. Determining the primordial ${\displaystyle P_{0}(\mathbf {k} )}$ is a question that relates to the physics of inflation.

The simplest ${\displaystyle P_{0}(\mathbf {k} )}$ is the Harrison–Zeldovich spectrum (named after Edward R. Harrison and Yakov Zeldovich), ^{[2] [3]} which characterizes ${\displaystyle P_{0}(\mathbf {k} )}$ according to a power law, ${\displaystyle P_{0}(\mathbf {k} )=Ak}$. More advanced primordial spectra include the use of a transfer function which mediates the transition from the universe being radiation dominated to being matter dominated.

The broad shape of the matter power spectrum is determined by the growth of large-scale structure, with the turnover (the point where the spectrum goes from increasing with k to decreasing with k) at ${\displaystyle k_{eq}\approx 2\times 10^{-2}~h~{\text{Mpc}}^{-1}}$, corresponding to ${\displaystyle \lambda _{eq}=350~h^{-1}~{\text{Mpc}}}$ (where h is the dimensionless Hubble constant). ^[4] The co-moving wavenumber corresponding to the maximum power in the mass power spectrum is determined by the size of the cosmic particle horizon at the time of matter-radiation equality, and therefore depends on the mean density of matter and to a lesser extent on the number of neutrino families (${\displaystyle N\geq 3}$), ${\displaystyle k_{\text{eq}}=(2\Omega _{M}H_{0}^{2}z_{\text{eq}}/c^{2})=7.3\times 10^{-2}\Omega _{M}h^{2}{\text{Mpc}}^{-1}}$, for ${\displaystyle T_{\text{cmb}}=2.728~\mathrm {K} }$. The ${\displaystyle P_{0}(k)}$ at smaller k (equivalently, larger scales) corresponds to scales which were larger than the particle horizon at the time of the transition from the regime of radiation dominance to that of matter dominance. ^{[5] [6]} At linear order in perturbations ${\displaystyle \delta }$, the power spectrum's broad shape follows

$${\displaystyle P({\bf {k,t)\propto {\begin{cases}k^{n_{s}}&k\ll k_{eq}\\k^{n_{s}-1}&k\gg k_{eq}\end{cases}}~,}}}$$

where ${\displaystyle n_{s}\simeq 1}$ is the scalar spectral index. ^[7]

References

1. Dodelson, Scott; Schmidt, Fabian (2020). Modern Cosmology - 2nd Edition. Academic Press. ISBN   978-0128159491.
2. Harrison, E. (1970). "Fluctuations at the threshold of classical cosmology". Physical Review. D1 (10): 2726–2730. Bibcode:1970PhRvD...1.2726H. doi:10.1103/PhysRevD.1.2726.
3. Zeldovich, Y. (1972). "A hypothesis, unifying the structure and the entropy of the universe". MNRAS. 160: 1P-3P. doi: 10.1093/mnras/160.1.1P .
4. Michael, Norman (2010). "Simulating Galaxy Clusters, 2. COSMOLOGICAL FRAMEWORK AND PERTURBATION GROWTH IN THE LINEAR REGIME".
5. Hu, Wayne; Sugiyama, Naoshi; Silk, Joseph (1997). "The physics of microwave background anisotropies". Nature. 386 (6620). Springer Science and Business Media LLC: 37–43. arXiv: astro-ph/9604166 . doi:10.1038/386037a0. ISSN   0028-0836.
6. Eisenstein, Daniel (1998). "Baryonic Features in the Matter Transfer Function". The Astrophysical Journal. 496 (2): 605. arXiv: astro-ph/9709112 . Bibcode:1998ApJ...496..605E. doi:10.1086/305424. S2CID   6505927.
7. Huterer, Dragan (2023). A course in cosmology: from theory to practice. Cambridge, United Kingdom New York, NY, USA: Cambridge University Press. ISBN   978-1-316-51359-0.

• Dodelson, Scott (2003). Modern Cosmology. Academic Press. ISBN   978-0-12-219141-1.
• Theuns, Physical Cosmology
• Michael L. Norman, Simulating Galaxy Clusters