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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 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/9506072 by the authors.
References
How later papers used it
545 later papers cite it
32 later papers by other authors are listed below: 32 used its method or tool.
Uses its method or tool 32 of 190 papers with this label
-
Modified the public code Class to implement the coupled Boltzmann equations of this paper.
Residual Non-Abelian Dark Matter and Dark Radiation (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
Computed the Newtonian potential from code written in synchronous gauge using a transformation involving the synchronous gauge potentials of this paper.
The effects of the small-scale behaviour of dark matter power spectrum on CMB spectral distortion (2017) · 1.00 · paraphrased: the licence does not allow quoting
-
Used the closing relations defined in this paper for both the first and the second order.
The cosmic microwave background bispectrum from the non-linear evolution of the cosmological perturbations (2010) · 1.00 · paraphrased: the licence does not allow quoting
-
Followed the truncation method introduced in this paper as implemented in its code.
Dark matter self-interactions in the matter power spectrum (2022) · 1.00 · paraphrased: the licence does not allow quoting
-
Employed this paper's general formalism for scalar perturbations to study how the growth of perturbations is reduced by scattering.
Partially Acoustic Dark Matter, Interacting Dark Radiation, and Large Scale Structure (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
Truncated the Boltzmann hierarchy at order l_max=5 using the approximation outlined in this paper.
Cosmological Signatures of a Mirror Twin Higgs (2018) · 1.00 · paraphrased: the licence does not allow quoting
-
Used the truncation prescription of this paper.
The Cosmological Evolution of Self-interacting Dark Matter (2021) · 1.00 · paraphrased: the licence does not allow quoting
-
Used a re-implementation of both the Newtonian and the comoving synchronous equations of this paper for its full linear theory calculations.
General relativistic 'screening' in cosmological simulations (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
Employed this paper's general formalism for scalar perturbations in the conformal Newtonian gauge.
Reconciling Large And Small-Scale Structure In Twin Higgs Models (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
We truncate the Boltzmann hierarchy with the traditional method in Ref. [ref] .
Cosmic microwave background spectral distortions from Rayleigh scattering at second order (2022) · Boltzmann Hierarchy equations · 1.00
-
Followed this paper to compute the fluid limit to the perturbed equations in k-space in the synchronous gauge for WDM.
Analytic Fluid Approximation for Warm Dark Matter (2019) · 1.00 · paraphrased: the licence does not allow quoting
-
Since we will keep perturbation in the these fluids to linear order, we will consider these perturbations in momentum space following the method derived in Ref. [ref] .
A novel integrated Sachs-Wolfe effect from Early Dark Energy (2023) · Linear and Nonlinear Numerical Methods · 1.00
-
Chose l_max=8 and obtained closed differential equations using the truncation relation of this paper.
Isocurvature Induced Gravitational Waves at Pulsar Timing Arrays (2025) · 1.00 · paraphrased: the licence does not allow quoting
-
Derived the initial conditions for the gravitational potentials including anisotropic stress contributions using the method described in this paper.
The Phase of the BAO on Observable Scales (2020) · 1.00 · paraphrased: the licence does not allow quoting
-
Effectively truncated the Boltzmann hierarchy by adopting the scheme of this paper.
Cosmological constraints on Lorentz invariance violation in the neutrino sector (2012) · 0.99 · paraphrased: the licence does not allow quoting
-
Followed the procedure given by this paper, in which the stress-energy tensor to linear order is given by an equation with density, pressure, anisotropic shear and fluid velocity perturbations.
Dark Matter Decaying into a Fermi Sea of Neutrinos (2010) · 0.99 · paraphrased: the licence does not allow quoting
-
Employed this paper's formalism to study the more subtle effects of dark gluons on the evolution of DM overdensities.
Non-Abelian dark matter and dark radiation (2015) · 0.99 · paraphrased: the licence does not allow quoting
-
Gave the linear perturbation theory equations in conformal Newtonian gauge for DE-DM interaction fluids according to this paper.
Observational constraints on dark matter-dark energy scattering cross section (2017) · 0.81 · paraphrased: the licence does not allow quoting
-
Numerically solved the coupled linearized Einstein, Boltzmann and fluid equations of this paper to describe the spatial evolution of gas around the black hole from BBN until decoupling.
Nucleosynthetic signatures of primordial origin around supermassive black holes (2021) · 0.98 · paraphrased: the licence does not allow quoting
-
Followed the evolution of cosmological perturbations using equations with the gravitational potential in the synchronous gauge of this paper.
Effects of electrically charged dark matter on cosmic microwave background anisotropies (2016) · 0.85 · paraphrased: the licence does not allow quoting
-
The conformal Newtonian gauge is parametrized as, [equation] To go from the synchronous gauge to the conformal Newtonian gauge, we use the transformation [ref] [equation] Once again, to get the leading behavior in the superhorizon limit, we approximate $η≫h$, [equation] implying [equation] Since we showed above that the sign of $σ(<0 or >0)$ depends upon the nature of DR (CDR or FDR), and we also showed $σ≈-3η$ on superhorizon scales, we see that $ϕ+ψ>0 (<0)$ for CDR (FDR).
Free-streaming and Coupled Dark Radiation Isocurvature Perturbations: Constraints and Application to the Hubble Tension (2021) · Effects of shear · 0.90
-
Solved numerically its equations together with the standard Boltzmann equations describing the evolution of baryons, photons and neutrinos from this paper.
The Cosmology of Atomic Dark Matter (2012) · 0.98 · paraphrased: the licence does not allow quoting
-
Analyzed the evolution equations of the linear perturbations in Fourier space and in the synchronous gauge of this paper.
The Early-Time Evolution of the Cosmological Perturbations in f(R) Gravity (2011) · 0.85 · paraphrased: the licence does not allow quoting
-
Worked in synchronous gauge following the notation and conventions of this paper.
Constraining Dark Matter-Baryon Scattering with Linear Cosmology (2013) · 0.98 · paraphrased: the licence does not allow quoting
-
Used equations from this paper governing the evolution of the density contrast and fluid velocity of unparticles in a spatially flat universe.
Emergent Unparticles Dark Energy can restore cosmological concordance (2023) · 0.86 · paraphrased: the licence does not allow quoting
-
We use the definitions for stress-energy components according to [ref] , which account for the curvature of FLRW space since the negative pressure of dark energy is a large-scale effect, and we quantize the field below with this in mind as well.
Deriving the Dark Matter-Dark Energy Interaction Term in the Continuity Equation from the Boltzmann Equation (2019) · Introduction · 0.94
-
Transformed an equation into k space and otherwise followed the standard procedure of linear perturbation theory of this paper to arrive at a hierarchy.
SU(2)$_{\rm CMB}$ and the cosmological model: angular power spectra (2018) · 0.97 · paraphrased: the licence does not allow quoting
-
Put the synchronous gauge metric into an equation of this paper to obtain four equations.
Minimal parameterizations for modified gravity (2013) · 0.88 · paraphrased: the licence does not allow quoting
-
Used the recursion relation between higher photon multipole moments of this paper to obtain expressions for F_gamma2 accurate to second order in tau_c.
Photons and Baryons before Atoms: Improving the Tight-Coupling Approximation (2010) · 0.98 · paraphrased: the licence does not allow quoting
-
Closely followed this paper to describe the perturbation equations for massive neutrinos.
Bounds on Neutrino Mass in Viscous Cosmology (2017) · 0.90 · paraphrased: the licence does not allow quoting
-
Gave the modified Boltzmann hierarchy for massless neutrinos following the same notation of this paper, also for the definition of the harmonic expansion coefficients.
Hints of Neutrino Dark Matter scattering in the CMB? Constraints from the Marginalized and Profile Distributions (2023) · 0.84 · paraphrased: the licence does not allow quoting
-
We have adopted the convention $r_≡a^4 ρ̅_/ρ_c,0$, as in Ref. [ref] , which in the 2-body decay scenario under study leads to: [equation] The first three multipoles are given by: [equation] We choose the maximum multipole $ℓ_max$ to truncate the hierarchies of equations according to the scheme proposed in Ref. [ref] for both massless and massive neutrinos, i.e.
Linear cosmological constraints on 2-body decaying dark matter scenarios and the $S_8$ tension (2021) · Linear perturbation equations · 0.93
How the literature uses it
555 arXiv papers cite this in 1,115 sentences
Papers citing it · 730 sentences without a confident label
Contests it 1 paper
-
The time dependence of the terms proportional to $c_3$ and $c_4$ would not satisfy the second order differential equation given in [ref] , which therefore does not admit, strangely enough, the standard inflationary solution for $π_ν$ given for example in [ref] .
CMB temperature anisotropy at large scales induced by a causal primordial magnetic field · Resolution for the neutrino anisotropic stress at leading order · contest 0.84
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
-
Within this local framework, the scattering amplitude relies solely on Lorentz invariants and is hence equivalent to the flat-space outcome [ref] .
Boltzmann Dynamics in K-essence Cosmology: Photon Propagation in an Emergent Spacetime · Local nature of collision terms and validity of flat-space cross sections · 0.96
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.