REVIEW
Power Spectra for Cold Dark Matter and its Variants
Astrophysics · Submitted 23 Oct 1997 · 75 references
Pith's verdict: A genuinely useful closed-form transfer function for mixed dark matter and its variants, with accuracy claims that are honest fit residuals rather than blind accuracy bounds.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The authors replace a Boltzmann integration with a single closed-form formula that reproduces the matter transfer function for cold dark matter and its common variants, including mixtures with massive neutrinos.
The reading
What carries the argument
The load-bearing device is the separation $T_{rm master}(k) = T_{rm sup}(k) B(k)$ combined with $T_{cb}(q,z) = T_{rm master}(q) D_{cb}(q,z)/D_1(z)$. $T_{rm sup}$ is a zero-baryon transfer function evaluated at an effective wavenumber $q_{rm eff} = k Theta_{2.7}^2 / (Gamma_{rm eff} {rm Mpc}^{-1})$, where $Gamma_{rm eff} = Omega_0 h^2 [sqrt{alpha_nu} + (1 - sqrt{alpha_nu})/(1 + (0.43 k s)^4)]$ rolls the curve from the pure-CDM shape at large scales to the analytically known small-scale suppression $alpha_nu$ at scales below the sound horizon. The growth functions $D_{cb}$ and $D_{cbnu}$ carry the free-streaming epoch $y_{fs}(q) = 17.2 f_nu (1 + 0.488 f_nu^{-7/6})(N_nu q / f_nu)^2$ and the altered logarithmic growth rates $p_i = [5 - sqrt{1 + 24 f_i}]/4$, so a single master curve plus a growth factor reproduces the whole time history. The correction $B(k) = 1 + 1.24 f_nu^{0.64} N_nu^{0.3 + 0.6 f_nu} / (q_nu^{-1.6} + q_nu^{0.8})$ patches the region around the scale where neutrinos finally became non-relativistic.
Extended reading notes
Core claim
The central claim is that the late-time transfer function of an adiabatic cosmology with cold dark matter, baryons, and massive neutrinos factors into two pieces: a time-independent master function fixed at the drag epoch, and a scale-dependent growth function that carries all the subsequent time evolution, cosmological constant, and curvature dependence. Because the master function depends only on $Omega_0 h^2$, the baryon fraction, the neutrino fraction, and $N_nu$, it can be fit once and reused at any redshift. The authors build that fit by anchoring two limits: on large scales the curve must match a pure-CDM form, and on small scales the amount of suppression from baryons and neutrinos is known analytically from their earlier work. They interpolate between the limits by rescaling the shape parameter $Gamma$ with an effective $Gamma_{rm eff}$ that switches on near the sound horizon $s$, then add a multiplicative correction $B(k)$ that repairs an overestimate of the maximal neutrino free-streaming scale, which arises because the fit assumes neutrino velocities keep scaling as $(1+z)$ when in fact only the momenta do. The result is a formula whose stated accuracy is better than 5 percent in the transfer function over the central parameter box, and better than 2 percent for $q > 0.25$.
Load-bearing premise
The accuracy claim is measured against a single Boltzmann code, CMBfast v2.3, and the paper itself notes in section 3.4 that different numerical codes disagree at the several-percent level for $q gg 1$; if the benchmark carries that much uncertainty on small scales, the quoted 5 percent agreement is bounded below by the benchmark's own error.
Editorial extensions
If this is right
- Cosmological parameter scans over baryon and neutrino fraction become cheap enough to map constraints in two dimensions, which the paper does for the power-spectrum shape ratio $sigma_{50}/sigma_8$, cluster abundance, damped Lyman-$alpha$ halos, and the Lyman-$alpha$ forest.
- Neutrino mass limits from small-scale power can be quoted as functions of $Omega_0$, $h$, $Omega_b$, and tilt rather than pinned to one fiducial cosmology, since the suppression depends on $Omega_nu/Omega_0$.
- Baryons and neutrinos both suppress small-scale power but at different scales, so the formula lets a fit separate the two contributions instead of treating them as one effective shape parameter.
- Any calculation that needs a linear matter power spectrum as input, from N-body initial conditions to weak-lensing predictions, can generate it without running a Boltzmann solver.
- The paper's worked examples show that a baryon fraction near 0.1 with a tilt of $n = 0.95$, or a second degenerate neutrino species, opens regions of parameter space that satisfy both large-scale structure and high-redshift constraints simultaneously.
Reading between the lines
- The real test of the accuracy claim is out-of-sample: comparing against a second, independently written Boltzmann solver at the edges of the stated box would turn a statement about agreement with one code into a statement about the true transfer function.
- Because the fit deliberately runs through the middle of the baryon acoustic oscillations rather than reproducing them, it is the wrong tool for the applications where those oscillations are the signal; a user wanting peak positions and amplitudes should use the oscillatory form instead.
- The separability of master function from scale-dependent growth is the reusable idea here; the same decomposition plausibly extends to warm or decaying dark matter, where free-streaming likewise breaks the clean split between space and time.
- The paper's own figures suggest the Lyman-$alpha$ forest slope constraint excludes more of the parameter space than cluster counts or damped systems do, an early hint that small-scale intergalactic medium data rather than cluster abundance would end up driving neutrino mass limits.
Original abstract
The bulk of recent cosmological research has focused on the adiabatic cold dark matter model and its simple extensions. Here we present an accurate fitting formula that describes the matter transfer functions of all common variants, including mixed dark matter models. The result is a function of wavenumber, time, and six cosmological parameters: the massive neutrino density, number of neutrino species degenerate in mass, baryon density, Hubble constant, cosmological constant, and spatial curvature. We show how observational constraints---e.g. the shape of the power spectrum, the abundance of clusters and damped Lyman-alpha systems, and the properties of the Lyman-alpha forest--- can be extended to a wide range of cosmologies, including variations in the neutrino and baryon fractions in both high-density and low-density universes.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
The bulk of recent cosmological research has focused on the adiabatic cold dark matter model and its simple extensions. Here we present an accurate fitting formula that describes the matter transfer functions of all common variants, including mixed dark matter models. The result is a function of wavenumber, time, and six cosmological parameters: the massive neutrino density, number of neutrino species degenerate in mass, baryon density, Hubble constant, cosmological constant, and spatial curvature. We show how observational constraints---e.g. the shape of the power spectrum, the abundance of clusters and damped Lyman-alpha systems, and the properties of the Lyman-alpha forest--- can be extended to a wide range of cosmologies, including variations in the neutrino and baryon fractions in both high-density and low-density universes.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | The bulk of recent cosmological research has focused on the adiabatic cold dark matter model and its simple extensions. | background 1.00 | 0.03 | 0.03 | asserted 1.00 |
| s1 | Here we present an accurate fitting formula that describes the matter transfer functions of all common variants, including mixed dark matter models. | contribution 0.99 | 0.96 | 0.52 | asserted 1.00 |
| s2 | The result is a function of wavenumber, time, and six cosmological parameters: the massive neutrino density, number of neutrino species degenerate in mass, baryon density, Hubble constant, cosmological constant, and spatial curvature. | result 0.55 | 0.91 | 0.23 | asserted 1.00 |
| s3 | We show how observational constraints---e.g. the shape of the power spectrum, the abundance of clusters and damped Lyman-alpha systems, and the properties of the Lyman-alpha forest--- can be extended to a wide range of cosmologies, including variations in the neutrino and baryon fractions in both high-density and low-density universes. | contribution 0.44 | 0.82 | 0.23 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:astro-ph/9710252 by the authors.
References
How later papers used it
176 later papers cite it
8 later papers by other authors are listed below: 8 used its method or tool.
Uses its method or tool 8 of 85 papers with this label
-
Used this paper's cold dark matter power spectrum parametrization for its warm dark matter models and applied a warm dark matter thermal cutoff from another work.
The median density of the Universe (2017) · 0.98 · paraphrased: the licence does not allow quoting
-
Computed the linear matter power spectrum using this paper's fitting formula, with the non-linear part from another work's Halofit code.
Biased Dark Energy Constraints from Neglecting Reduced Shear in Weak Lensing Surveys (2008) · 1.00 · paraphrased: the licence does not allow quoting
-
Used another work's halo mass function form with this paper's transfer function to generate the halo mass function, the SMF and other quantities.
The Dual Role of Outflows in Quenching Satellites of Low-Mass Hosts: NGC 3109 (2022) · 0.87 · paraphrased: the licence does not allow quoting
-
Obtained the power spectrum through this paper's analytic prescription, with initial conditions generated by another work.
Painting halos from cosmic density fields of dark matter with physically motivated neural networks (2019) · 1.00 · paraphrased: the licence does not allow quoting
-
Simulated scaled log-normal fields with this paper's linear matter power spectrum as its P_LN, implemented in another work's module.
Lossless, Scalable Implicit Likelihood Inference for Cosmological Fields (2021) · 0.98 · paraphrased: the licence does not allow quoting
-
The linear matter power spectrum used in PyCosmo is calculated using the Eisenstein & Hu fitting formulas [ref] , describing the matter transfer functions, including massive neutrinos for a degenerate mass hierarchy.
Combined $13\times2$-point analysis of the Cosmic Microwave Background and Large-Scale Structure: implications for the $S_8$-tension and neutrino mass constraints (2021) · Analytical Predictions · 0.99
-
In practice, the growth module selects among these four solutions according to the cosmological content of the model: the [ref] approximation for standard flat $Λ$CDM, the Heath-Peebles integral for non-flat geometries, a radiation-aware hybrid for cases where relativistic species are included, and a Linder-Jenkins ODE for evolving dark energy.
PHANTOM: A MATLAB and Octave Toolbox Connecting Linear Field Statistics to Dark Matter Halo Observables (2026) · Linear growth factor · 0.97
-
Adopted this paper's form of the linear matter power spectrum, along with another work's Sheth-Tormen halo mass function and another work's halo bias form.
The power spectrum of extended [C II] halos around high redshift galaxies (2023) · 0.98 · paraphrased: the licence does not allow quoting
How the literature uses it
95 arXiv papers cite this in 121 sentences
Papers citing it · 55 sentences without a confident label
Builds on it 1 paper
-
Extensions to the fitting formula of Ref. [ref] to include a non-zero $f_ν$ and $w ≠-1$ have been investigated in Ref. [ref] and [ref] respectively.
Cosmological parameters from large scale structure - geometric versus shape information · Halo power spectrum data in extended models · 0.84
Uses its method or tool 43 papers
-
In this approximation, if $𝒯_NCDM$ and $𝒯_CDM$ are the transfer functions of NCDM and CDM, respectively, the non-linear NCDM power spectrum is [ref] [equation] Both the transfer functions, relating the primordial and the present day power spectra [ref] , and the non-linear CDM power spectrum $P_δ, CDM$ are obtained using CLASS [ref] , with parameters $m_a$, $Ω_a$ from Eq. [ref] and $T_a$ from Eq. [ref], computed using the value of $g_*,s$ satisfying the structure formation bound.
Probing the Blue Axion with Cosmic Optical Background Anisotropies · Modified Freeze-out · 1.00
-
This analysis is performed using [<https://github.com/santiagocasas/fishermathica>] [ref] , a multipurpose code designed for Fisher analysis. works in cooperation with [<https://github.com/santiagocasas/cosmomathica>] [ref] , a wrapper for Einstein-Boltzmann codes [ref] and [ref] , fitting function codes [ref] , and [ref] . has been extensively tested and validated within the Euclid IST:Forecasts code comparison effort.
Quintessential $\alpha$-attractor inflation: forecasts for Stage IV galaxy surveys · Fiducial parameters and Fisher matrix analysis · 1.00
-
I compute the linear matter power spectrum using the fitting formula of [ref] , and the non-linear part is computed from the “Halofit” code of [ref] which has also been calibrated by numerical simulations.
Biased Dark Energy Constraints from Neglecting Reduced Shear in Weak Lensing Surveys · Perturbative Correction · 1.00
Uses its data or values 3 papers
-
We employ the transfer function for CDM plus baryons from Ref. [ref] .
The Halo Model of Large Scale Structure for Warm Dark Matter · Cold Dark Matter · 0.94
-
For the matter distribution, we used the best-fit cosmological parameters from WMAP5 [ref] and the linear power spectrum of Eisenstein and Hu [ref] with neutrino free-streaming and gravitational wave effects neglected.
Effects of P-wave Annihilation on the Angular Power Spectrum of Extragalactic Gamma-rays from Dark Matter Annihilation · The point distribution of halos · 0.95
-
Cosmological parameters used were from the WMAP7 data [ref] , and we use the linear power spectrum of [ref] to describe large-scale fluctuations in the matter distribution.
Extragalactic and galactic gamma-rays and neutrinos from annihilating dark matter · Gamma-ray Results · 1.00
Cites it as background 8 papers
-
In Ref. [ref] , Eisenstein and Hu provided a fitting formula when baryons, CDM and massive neutrinos are considered.
Using machine learning to compress the matter transfer function $T(k)$ · Previous Fitting Formulae · 0.83
-
It is well known that the CMB and matter power spectra constrain the sum of the neutrino masses $∑m_ν$ because while neutrinos freestream, they suppress power on scales smaller than the horizon when neutrinos become nonrelativistic [ref] .
Massive neutrinos and interacting dark matter look alike through the lens of lensing · Introduction · 0.90
-
Readers can refer to Ref. [ref] for more complex nonlinear evolution of matter trispectrum in the context of primordial parity violation.], the matter power spectrum could be written as [ref] , [equation] Here $𝒯(k)$ is the transfer function, $D(χ)$ is the growth factor, $𝒫_ℛ(k)$ denotes the primordial power spectrum of the comoving curvature perturbations and $δ^(3) _D$ is the Dirac delta function.
Probing Parity Violation with Weak Lensing Trispectrum · Weak Lensing Convergence Power Spectrum · 0.86
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.