Pith. Sign in

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

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/9710252 (version 1)
Pith IDNDT2PYZ2, Pith's permanent record of this paper and its review
CitationDaniel J. Eisenstein, Princeton), Wayne Hu (IAS. "Power Spectra for Cold Dark Matter and its Variants." arXiv:astro-ph/9710252 (1997). https://pith.science/paper/NDT2PYZ2
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

Working out how gravity amplified density ripples in the early universe normally means solving coupled Boltzmann equations mode by mode for every choice of cosmology. This paper sets out to replace that integration with a formula: an expression in wavenumber, redshift, and six cosmological parameters (massive neutrino density, number of degenerate neutrino species, baryon density, Hubble constant, cosmological constant, and spatial curvature) that reproduces the linear matter transfer function for adiabatic cold dark matter and its common variants, mixed dark matter included. If it holds, a theorist can sweep a large parameter space in seconds and read off the shape of the power spectrum, the abundance of clusters, the count of damped Lyman-$alpha$ halos, and the amplitude of Lyman-$alpha$ forest fluctuations without touching a Boltzmann code. The authors claim residuals below 5 percent in the transfer function (10 percent in power) for $0.06 < Omega_0 h^2 < 0.4$, baryon and neutrino fractions each below 0.3, and $z < 30$, improving to about 3 percent once the baryon fraction drops below 10 percent. They also demonstrate the payoff directly, mapping four observational constraints across two-dimensional slices of baryon and neutrino fraction in four fiducial cosmologies.

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 3,745 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.11 · human participants 0.01 · controlled comparison 0.05 · randomized assignment 0.01 · states quantitative result 0.21 · compares with prior results 0.22 · mentions code or data availability 0.04.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: method or tool 0.82 · evidence from simulation or computation 0.87.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. — Comparison of the fitting formula to the numerical re
Figure 2. — Same as Figure 1, but for different choices of cosmolo
Figure 3. — (a–d) Ratio of fluctuations within 50h −1 Mpc spheres to that within 8h −1 Mpc spheres as a function of neutrino and baryon fractions for various cosmologies. The shaded region shows the preferred range of σ50/σ8 = 0.151 ± 0.016. (a) Ω0 = 1, h = 0.5, n = 1, Nν = 1. (b) As (a), but with n = 0.95 and tensors as per power-law inflation. (c) As (a), but with two degenerate neutrino species (Nν = 2). (d) Ω0 = 0.35, ΩΛ = 0.65, h = 0.7, n = 1, Nν = 1. All are COBE-normalized. (e–h) Amplitude of fluctuations within 8h −1 Mpc spheres for the cosmologies given in (a–d), respectively. The…
Figure 4. — (a–d) Amplitude of fluctuations at z = 4 within a Gaussian window of mass corresponding to halos of 50 km s−1 Mpc−1 circular velocity. Cosmologies are as per

Figures reproduced from arXiv:astro-ph/9710252 by the authors.

References
How later papers used it

176 later papers cite it

Evidence: substantial. 171 of the citing papers share no author with it, and 92 of those carry a confident label from Pith's citation classifier.

Pith read 23 of the 215 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 8 of 85 papers with this label
How the literature uses it

95 arXiv papers cite this in 121 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 55 sentences without a confident label

Builds on it 1 paper
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
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.