Pith. Sign in

REVIEW

Cosmological Perturbation Theory in the Synchronous and Conformal Newtonian Gauges

Astrophysics · Submitted 11 Jun 1995 · 40 references

Pith's verdict: The standard reference for linear perturbation theory in two gauges, and it holds up — the only thing worth pressing is that the 10^-3 accuracy is an estimate, not a bound.

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 · 9 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/9506072 (version 1)
Pith IDPPFIXYDP, Pith's permanent record of this paper and its review
CitationChung-Pei Ma, Edmund Bertschinger. "Cosmological Perturbation Theory in the Synchronous and Conformal Newtonian Gauges." arXiv:astro-ph/9506072 (1995). https://pith.science/paper/PPFIXYDP
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 paper derives the linear theory of cosmic structure formation twice over — once in the synchronous gauge and once in the conformal Newtonian gauge — for every particle species that matters, and uses the agreement between the two to make

The reading

The paper asks a bookkeeping-shaped question with real physical stakes: if you write the linear theory of gravitational instability in two different coordinate choices — the synchronous gauge used since Lifshitz, and the conformal Newtonian gauge — do you get the same universe? Ma and Bertschinger answer by deriving side by side the complete set of Einstein, Boltzmann, and fluid equations for every species that matters in a flat cold dark matter, hot dark matter, or mixed model: CDM, baryons, photons, massless neutrinos, and massive neutrinos, with a cosmological constant allowed. The payoff is a numerical pipeline that can be run twice, once per gauge, and checked against itself; the two independent integrations agree to better than 0.1 percent. On that foundation they produce the first angular power spectra for mixed dark matter models that retain photon polarization, the full neutrino multipole hierarchy, and helium recombination. For a reader who wants to know why cosmic microwave background predictions can be trusted at the percent level, the paper's real contribution is the demonstration that a gauge choice is a computational convenience rather than a physical assumption.

What carries the argument

Three pieces carry the argument. The first is the gauge transformation of equation (18), which converts between $(h, eta)$ and $(psi, varphi)$ and is what makes a two-code cross-check possible at all. The second is the parallel set of Boltzmann hierarchies: the photon and massless-neutrino distribution functions are expanded in Legendre polynomials and reduced to coupled ordinary differential equations in equations (49), (50), (63), and (64), while the massive neutrinos keep an explicit momentum dependence and are tracked on a grid of 128 $q$-points per wavenumber with their own hierarchy in equations (56) and (57). The third is the set of truncation rules — equations (51), (58), and (65) — which close the infinite hierarchies by extrapolating the spherical-Bessel-like behavior of the high multipoles rather than simply setting them to zero, plus the tight-coupling approximation of §5.7 that replaces the stiff Thomson-scattering terms with a sound-speed equation for the combined photon-baryon fluid. Together these turn an intractable integro-differential problem into a finite system that a Runge-Kutta integrator can carry from redshift $z sim 10^9$ to today.

Extended reading notes

Core claim

The central claim is that a single, complete linear perturbation theory can be written down for all the species relevant to flat CDM, HDM, and CDM+HDM cosmologies, and that solving it in two different gauges gives the same answer. The bridge between the gauges is an explicit coordinate transformation, equation (18), which expresses the conformal Newtonian potentials $psi$ and $varphi$ in terms of the synchronous fields $h$ and $eta$ and their time derivatives. With that bridge in hand, the authors derive super-horizon isentropic initial conditions in both gauges — equations (96) and (98) — showing that the same physical mode appears as a growing density perturbation in the synchronous gauge and as a constant one in the conformal Newtonian gauge, with the difference traceable entirely to the coordinate choice. They then integrate the coupled system numerically for a standard CDM model and two mixed models with neutrino mass densities $Omega_nu = 0.2$ and $0.3$ (neutrino masses 4.7 eV and 7 eV), assuming a scale-invariant adiabatic spectrum. The two gauge integrations agree to better than 0.1 percent, and the resulting spectra are, on the authors' account, the first for mixed dark matter models that include photon polarization, higher neutrino multipole moments, and helium recombination.

Load-bearing premise

The quoted accuracy of the spectra rests on an estimate — grids, multipole ceilings, and timesteps doubled until the answer stopped moving — rather than on a proven error bound, so everything hangs on the assumption that the closure rules of equations (51), (58), and (65) and the tight-coupling switch temperatures of §5.7 converge for every wavenumber and multipole the paper uses.

Editorial extensions

If this is right

  • Adding a massive neutrino at $Omega_nu = 0.2$ or $0.3$ raises the anisotropy at high multipoles by roughly 10 percent at the second and third acoustic peaks, so a percent-level measurement of peak heights separates mixed dark matter from pure CDM.
  • Free-streaming of the massive neutrino suppresses its own perturbation growth and slows the growth of CDM in turn, with the suppression scale set by $k_{fs}(a) propto a^{1/2} m_ u$, so the shape of the matter power spectrum becomes a handle on neutrino mass.
  • Because the two gauges agree to 0.1 percent, the gauge-invariant content of the calculation is confirmed, and a user may choose whichever gauge simplifies an extension — for example, integrating neutrino geodesics to build N-body initial conditions.
  • The released code makes the entire pipeline reproducible and extendable to a cosmological constant and to additional species without rederiving the perturbation equations.

Reading between the lines

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

  • If two independently written Boltzmann codes ever disagree by more than about 0.1 percent on a mixed dark matter spectrum, that disagreement is more likely a truncation artifact or a coding error than new physics — the paper effectively sets a benchmark for what 'agreeing codes' should mean.
  • The weakest link in the numerical chain is the hierarchy closure, which rests on an asymptotic analogy to spherical Bessel functions; a genuinely independent check would come from a line-of-sight integration that never truncates the multipole expansion at all.
  • The sensitivity of the acoustic peak heights to $Omega_ nu$ runs in both directions: the same machinery, inverted, turns a measured peak ratio into a constraint on the sum of neutrino masses.
  • The gauge transformation itself is a reusable tool beyond this paper — any calculation performed in the synchronous gauge can be mapped into the conformal Newtonian gauge where the potentials have a Newtonian interpretation, which is exactly what is needed when feeding linear results into Newtonian N-body simulations.
Original abstract

This paper presents a systematic treatment of the linear theory of scalar gravitational perturbations in the synchronous gauge and the conformal Newtonian (or longitudinal) gauge. It differs from others in the literature in that we give, in both gauges, a complete discussion of all particle species that are relevant to any flat cold dark matter (CDM), hot dark matter (HDM), or CDM+HDM models (including a possible cosmological constant). The particles considered include CDM, baryons, photons, massless neutrinos, and massive neutrinos (an HDM candidate), where the CDM and baryons are treated as fluids while a detailed phase-space description is given to the photons and neutrinos. Particular care is applied to the massive neutrino component, which has been either ignored or approximated crudely in previous works. Isentropic initial conditions on super-horizon scales are derived. The coupled, linearized Boltzmann, Einstein and fluid equations that govern the evolution of the metric and density perturbations are then solved numerically in both gauges for the standard CDM model and two CDM+HDM models with neutrino mass densities $\onu=0.2$ and 0.3, assuming a scale-invariant, adiabatic spectrum of primordial fluctuations. We also give the full details of the cosmic microwave background anisotropy, and present the first accurate calculations of the angular power spectra in the two CDM+HDM models including photon polarization, higher neutrino multipole moments, and helium recombination. The numerical programs for both gauges are available at http://arcturus.mit.edu/cosmics/ .

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, 6,594 input tokens, $0.0003 reported.

What the abstract itself states, as probabilities: claims first 0.44 · states limitation 0.08 · human participants 0.01 · controlled comparison 0.05 · randomized assignment 0.01 · states quantitative result 0.65 · compares with prior results 0.84 · mentions code or data availability 0.98.

The abstract, sentence by sentence

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

This paper presents a systematic treatment of the linear theory of scalar gravitational perturbations in the synchronous gauge and the conformal Newtonian (or longitudinal) gauge. It differs from others in the literature in that we give, in both gauges, a complete discussion of all particle species that are relevant to any flat cold dark matter (CDM), hot dark matter (HDM), or CDM+HDM models (including a possible cosmological constant). The particles considered include CDM, baryons, photons, massless neutrinos, and massive neutrinos (an HDM candidate), where the CDM and baryons are treated as fluids while a detailed phase-space description is given to the photons and neutrinos. Particular care is applied to the massive neutrino component, which has been either ignored or approximated crudely in previous works. Isentropic initial conditions on super-horizon scales are derived. The coupled, linearized Boltzmann, Einstein and fluid equations that govern the evolution of the metric and density perturbations are then solved numerically in both gauges for the standard CDM model and two CDM+HDM models with neutrino mass densities $\onu=0.2$ and 0.3, assuming a scale-invariant, adiabatic spectrum of primordial fluctuations. We also give the full details of the cosmic microwave background anisotropy, and present the first accurate calculations of the angular power spectra in the two CDM+HDM models including photon polarization, higher neutrino multipole moments, and helium recombination. The numerical programs for both gauges are available at http://arcturus.mit.edu/cosmics/ .

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 This paper presents a systematic treatment of the linear theory of scalar gravitational perturbations in the synchronous gauge and the conformal Newtonian (or longitudinal) gauge. contribution 0.99 0.93 0.07 asserted 1.00
s1 It differs from others in the literature in that we give, in both gauges, a complete discussion of all particle species that are relevant to any flat cold dark matter (CDM), hot dark matter (HDM), or CDM+HDM models (including a possible cosmological constant). contribution 0.98 0.85 0.07 asserted 1.00
s2 The particles considered include CDM, baryons, photons, massless neutrinos, and massive neutrinos (an HDM candidate), where the CDM and baryons are treated as fluids while a detailed phase-space description is given to the photons and neutrinos. contribution 0.56 0.86 0.04 asserted 1.00
s3 Particular care is applied to the massive neutrino component, which has been either ignored or approximated crudely in previous works. gap 0.89 0.23 0.05 asserted 0.99
s4 Isentropic initial conditions on super-horizon scales are derived. contribution 0.59 0.94 0.06 asserted 0.99
s5 The coupled, linearized Boltzmann, Einstein and fluid equations that govern the evolution of the metric and density perturbations are then solved numerically in both gauges for the standard CDM model and two CDM+HDM models with neutrino mass densities $\onu=0.2$ and 0.3, assuming a scale-invariant, adiabatic spectrum of primordial fluctuations. method 0.94 0.97 0.09 asserted 1.00
s6 We also give the full details of the cosmic microwave background anisotropy, and present the first accurate calculations of the angular power spectra in the two CDM+HDM models including photon polarization, higher neutrino multipole moments, and helium recombination. contribution 0.58 0.96 0.17 asserted 1.00
s7 The numerical programs for both gauges are available at http://arcturus.mit.edu/cosmics/ . contribution 0.69 0.72 0.02 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 0.53 · evidence from simulation or computation 0.98.

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. — Top: The angular power spectrum Cl of the photon anisotropy (including polarization) in the CDM (solid), the Ων = 0.2 CDM+HDM (dashed), and the Ων = 0.3 CDM+HDM (dotted) models. Middle: The ratio of Cl for the CDM+HDM models relative to the CDM (dashed for Ων = 0.2; dotted for Ων = 0.3). Bottom: The fractional difference in Cl with and without polarization for the three models
Figure 2. — The scalar metric perturbations φ(k,τ ) and ψ(k,τ ) in the conformal Newtonian gauge in the Ων = 0.2 CDM+HDM model. The 31 curves from left to right correspond to 31 values of k between 10.0 Mpc−1 and 0.01 Mpc−1 . The labels τnr, τeq and τrec indicate, respectively, the time the 4.7 eV neutrinos become non-relativistic, the matter-radiation equality time, and the recombination time
Figure 3. — Evolution of the density fields in the synchronous ga
Figure 4. — The z = 0 power spectra for the pure CDM (solid), the Ων = 0.2 CDM+HDM (dashed), and the Ων = 0.3 CDM+HDM (dotted) models. In the mixed models, the upper curve is for the CDM component and the lower one is for the HDM. The COBE-compatible quadrupole moment of Qrms−PS = 17.6 µ K is used

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

References
How later papers used it

545 later papers cite it

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

Pith read 40 of the 1,100 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 32 of 190 papers with this label
How the literature uses it

555 arXiv papers cite this in 1,115 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 · 730 sentences without a confident label

Contests it 1 paper
Builds on it 6 papers
  • The treatment of paper II generalises and extends the standard spatial Boltzmann hierarchy for early-universe perturbations [ref] to describe the full spectro-spatial evolution of the photon field.

    Spectro-spatial evolution of the CMB III: transfer functions, power spectra and Fisher forecasts · Brief recap of the important equations from paper I+II · 0.91

  • This is a straightforward generalization, suitable to the Newtonian evolution (i.e. let $ϵ= √(q^2+a^2m^2) →am$), of the truncation scheme found in [ref] .

    Neutrinos in N-body simulations · Anisotropies of Newtonian and non-relativistic NCDM components · 0.98

  • The first-order equation of motion ([ref]) for the perturbed phase space density $F_i$ can also be recast as an equation of motion for the phase space density contrast $Ψ_i (𝐤,𝐪_i,τ)= F_i(𝐤,𝐪_i,τ)/f̅_i(|𝐪_i|,τ)$: [equation] This is a generalisation of equation (40) of reference [ref] to include a time-dependent background distribution.

    Invisible neutrino decay in precision cosmology · First-order perturbation equations · 0.89

Uses its method or tool 194 papers
  • Considering the drag force [ref] as the only non-gravitational force between DE and DM, the linear perturbation theory equations (in conformal newtonian gauge) for DE-DM interaction fluids are given by [ref]

    Observational constraints on dark matter-dark energy scattering cross section · Dark matter-dark energy scattering models · 0.81

  • we have modified the equations for the time evolution of the density, velocity, shear and higher anisotropic moments of the phase space distribution in the synchronous gauge [ref] and we have included them in CLASS [ref] , using MontePython [ref] to deduce the cosmological implications of this scenario.

    Cosmology of an Axion-Like Majoron · Cosmological evolution: $T<10$ eV · 1.00

  • To visualize DM-DR's effects on matter power spectrum, we have modified the public code Class [ref] to implement the coupled Boltzmann equations [ref] , [equation] where dot means derivative over conformal time $dτ≡dt/a$ ($a$ is the scale factor), $θ_r$ and $θ_m$ are velocity divergences of DR $A^a$ and DM $A^m$, $k$ is the co-moving wavenumber, $Ψ$ is the gravitational potential, $δ_r$ and $σ_r$ ($σ_r=0$ for $A^a$ as perfect fluid) are the density perturbation and the anisotropic stress potential of $A^a$, and $ℋ ≡ȧ/a$ is the conformal Hubble parameter.

    Residual Non-Abelian Dark Matter and Dark Radiation · Numeric Results · 1.00

Compares with it 12 papers
  • With $π_Ai = ∂_i π_A$, the spatial parts are recast into [equation] where we write $Σ$ and $ξ_Σ$ as [equation] Then, the second-order action yields the first-order continuity and Euler equations in the presence of the stochastic noise: [equation] where a hat implies [equation] Setting the noise to zero, Eqs.[ref] and [ref] are consistent with the linearized continuity and Euler equations presented in Ref. [ref] .

    Fluctuation-dissipation relation in cosmic microwave background · Equation of motion with noise · 0.95

  • In the pure CDM limit, where $γ_0 = 0$, we can set $f_1 = 1$ and $f_2 = γ_E - 1/2$, which satisfies Eq. ([ref]) and the initial condition of $δ_χ = -1/2 C (k τ)^2$ for $x →0$ and $θ_χ = 0$, coinciding with that for the adiabatic mode given in Eq. (96) of Ref. [ref] .

    Dark matter kinetic decoupling with a light particle · Evolution of the cosmological perturbations · 0.86

  • The difference here with the formulas from the Ref. [ref] is that we have fixed the Newtonian gauge from the very beginning, neglected the vector and tensor perturbations, and changed the names of the Newtonian potentials, $ϕ↔ψ$, in order to be in accordance with the Ref. [ref] .

    Toward a concordance teleparallel Cosmology II: Linear perturbation · Linear perturbations in $f(T)$ gravity · 0.89

Cites it as support 1 paper
Cites it as background 93 papers
  • The formalism describing the evolution of matter and dark energy density perturbations without [ref] and with dark matter and dark energy interaction [ref] is well established.

    Testing the Interaction between Dark Energy and Dark Matter with Planck Data · The phenomenological model on the interaction between dark sectors · 0.90

  • In linear cosmological perturbation theory, the Boltzmann equation of NCDM particles depends on the amplitude of the particle momentum and the angle between the momentum and the Fourier wavevector $k⃗$ [ref] .[The case for massless particles is simpler, since the dependence on the norm of the momentum can be factored out.] It is typically rewritten as an infinite set of coupled differential equations for the momentum-dependent multipole moments of the perturbed phase-space-density distribution function.

    Rapid and accurate numerical evolution of linear cosmological perturbations with non-cold relics · Motivation – Limitations of Boltzmann hierarchy approach · 0.93

  • Among other features, offers a fast solver for Einstein-Boltzmann equations, a large system of ODEs which describes the evolution of linear perturbations in the Universe [ref] .

    sympy2c: from symbolic expressions to fast C/C++ functions and ODE solvers in Python · Introduction · 0.91

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.