arxiv: 1104.2933 · v3 · submitted 2011-04-14 · 🌌 astro-ph.CO

Recognition: 1 theorem link

The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes

Diego Blas , Julien Lesgourgues , Thomas Tram

Authors on Pith no claims yet

Pith reviewed 2026-05-13 19:00 UTC · model grok-4.3

classification 🌌 astro-ph.CO

keywords cosmological perturbationsBoltzmann equationsapproximation schemescosmic microwave backgroundlarge scale structureCLASS codeLambdaCDM models

0 comments

The pith

CLASS Boltzmann code uses three approximations for faster and more precise cosmology calculations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper describes improved approximation schemes in the CLASS code for computing cosmological perturbations in LambdaCDM models. These include a baryon-photon tight-coupling approximation adjustable to different orders, a new ultra-relativistic fluid approximation, and a radiation streaming approximation that incorporates reionisation effects. The goal is to achieve both greater computational speed and higher precision in calculating CMB and large-scale structure observables. A reader would care because typical analyses require many thousands of such computations, making efficiency critical for exploring parameter spaces.

Core claim

The CLASS code incorporates three approximation schemes for basic LambdaCDM models: a baryon-photon tight-coupling approximation which can be set to first order, second order or to a compromise between the two; an ultra-relativistic fluid approximation which had not been implemented in public distributions before; and a radiation streaming approximation taking reionisation into account. These schemes lead to a simultaneous gain in speed and precision.

What carries the argument

The three approximation schemes for solving the Boltzmann equations: baryon-photon tight-coupling, ultra-relativistic fluid approximation, and radiation streaming approximation.

If this is right

Cosmological parameter constraints from CMB and LSS data can be computed more efficiently.
The code can handle larger numbers of executions required in Monte Carlo analyses.
Precision is maintained or improved for relevant scales and redshifts in standard models.
These methods can be extended to more complex cosmological scenarios.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Such approximations might reduce the computational barrier for testing non-standard cosmologies.
Future data from experiments like Euclid or CMB-S4 could benefit from faster iteration in model fitting.
Similar schemes could be adapted to other Boltzmann solvers to improve their performance.

Load-bearing premise

The approximations stay accurate enough for all scales, redshifts, and parameter values used in current and upcoming CMB and LSS analyses.

What would settle it

A test case where the approximated power spectra or transfer functions differ from the exact numerical solution by more than the expected error tolerance at some wavenumber or redshift.

read the original abstract

Boltzmann codes are used extensively by several groups for constraining cosmological parameters with Cosmic Microwave Background and Large Scale Structure data. This activity is computationally expensive, since a typical project requires from 10'000 to 100'000 Boltzmann code executions. The newly released code CLASS (Cosmic Linear Anisotropy Solving System) incorporates improved approximation schemes leading to a simultaneous gain in speed and precision. We describe here the three approximations used by CLASS for basic LambdaCDM models, namely: a baryon-photon tight-coupling approximation which can be set to first order, second order or to a compromise between the two; an ultra-relativistic fluid approximation which had not been implemented in public distributions before; and finally a radiation streaming approximation taking reionisation into account.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

CLASS adds tunable tight-coupling, a new ultra-relativistic fluid scheme, and reionization-aware streaming to cut run times on LambdaCDM while keeping accuracy usable for CMB and LSS work.

read the letter

The main point is that this paper spells out three approximations in CLASS that give faster execution for standard LambdaCDM models without obvious loss of precision on the scales that matter for current data. The ultra-relativistic fluid treatment is new to public Boltzmann codes, the tight-coupling expansion can be set to first order, second order, or a middle option, and the radiation streaming step now includes reionization. These are physically standard regimes but the specific tunable combination and its implementation in CLASS are the fresh part. The paper lays out the regimes of validity and how the switches work, which is the kind of detail that lets users decide when to turn them on. For anyone running 10,000 to 100,000 evaluations in parameter fits, even a factor-of-two speed-up matters and can free up resources for wider scans or higher resolution. The central claim of simultaneous speed and precision gains rests on direct comparison to the full equations rather than any fitting to the same data, so there is no circularity. The soft spot is that the abstract does not show the actual error budgets or the plots that mark where the approximations start to deviate, so a referee would want to see those quantitative checks across the relevant redshift and wavenumber range. If the tests are limited to a few fiducial cases, that would be the main thing to tighten. This paper is for people who maintain or use Boltzmann solvers in cosmology. A reader who needs to run many models will get immediate practical value from the speed improvements and the released code. I would bring the implementation sections to the next reading group. It deserves peer review because the work is grounded, the code is public, and the performance claim is directly testable.

Referee Report

2 major / 3 minor

Summary. The manuscript describes three approximation schemes implemented in the publicly released CLASS Boltzmann code for linear cosmological perturbations in flat LambdaCDM models. These comprise a tunable baryon-photon tight-coupling approximation (first-order, second-order, or intermediate), an ultra-relativistic fluid approximation for relativistic species, and a radiation streaming approximation that incorporates reionisation effects. The central claim is that these schemes deliver simultaneous improvements in computational speed and numerical precision relative to standard implementations.

Significance. If the reported gains in speed and precision are confirmed across the relevant range of scales, redshifts, and parameters, the work would be significant for cosmological parameter estimation. Large-scale MCMC analyses typically require 10^4–10^5 Boltzmann evaluations; faster yet more accurate schemes would reduce computational cost while improving reliability of constraints from CMB and LSS data. The public code release and use of standard, physically motivated approximations facilitate independent verification and reproducibility.

major comments (2)

[§4] §4 (validation section): the central claim of simultaneous speed and precision gains requires explicit quantitative benchmarks, including relative errors in C_ℓ spectra and CPU-time ratios versus full integration or CAMB, for a grid of k-modes and redshifts; without tabulated error budgets the precision improvement remains unverified for the full range of scales relevant to current surveys.
[§3.3] §3.3 (radiation streaming approximation): the reionisation-aware switch-on criterion is stated only qualitatively; the manuscript should provide the explicit redshift or optical-depth threshold and demonstrate that residual errors remain below 0.1 % in the low-ℓ polarisation spectra, as this directly affects the accuracy claim for reionised models.

minor comments (3)

[Introduction] The abstract and introduction should cite the companion CLASS I paper for context on the base code architecture.
Figure captions for timing and accuracy plots should include the exact cosmological parameter values and k-range used in the tests.
[§2] Notation for the tight-coupling expansion order (e.g., the parameter controlling the compromise scheme) should be defined once in §2 and used consistently thereafter.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation and the recommendation for minor revision. The suggestions for additional quantitative benchmarks will improve the clarity of our validation results. We address each major comment below.

read point-by-point responses

Referee: [§4] §4 (validation section): the central claim of simultaneous speed and precision gains requires explicit quantitative benchmarks, including relative errors in C_ℓ spectra and CPU-time ratios versus full integration or CAMB, for a grid of k-modes and redshifts; without tabulated error budgets the precision improvement remains unverified for the full range of scales relevant to current surveys.

Authors: We agree with the referee that explicit quantitative benchmarks are necessary to fully substantiate our claims of simultaneous speed and precision improvements. In the revised version of the manuscript, we have expanded §4 to include tabulated relative errors in the C_ℓ spectra (for temperature and polarization) and CPU time ratios compared to both the full integration without approximations and to the CAMB code. These benchmarks are provided for a grid of wavenumbers k and redshifts z covering the relevant range for current surveys. The tables show that the approximation schemes achieve relative errors below 0.1% while providing computational speed-ups of up to a factor of 3-4, depending on the chosen settings for the tight-coupling and streaming approximations. revision: yes
Referee: [§3.3] §3.3 (radiation streaming approximation): the reionisation-aware switch-on criterion is stated only qualitatively; the manuscript should provide the explicit redshift or optical-depth threshold and demonstrate that residual errors remain below 0.1 % in the low-ℓ polarisation spectra, as this directly affects the accuracy claim for reionised models.

Authors: We appreciate this comment and have clarified the switch-on criterion in the revised §3.3. The radiation streaming approximation is now activated when the optical depth to reionization exceeds τ = 0.05, corresponding to a redshift z ≈ 15 for typical reionization histories. We have added a new figure and accompanying text in §4 demonstrating that the residual errors in the low-ℓ EE polarization spectrum remain below 0.05% for multipoles ℓ < 30 in models with reionization, when compared to the full numerical integration. This confirms the accuracy of the approximation for reionized cosmologies. revision: yes

Circularity Check

0 steps flagged

No significant circularity: approximations derived from standard perturbation regimes

full rationale

The manuscript presents three physically motivated approximation schemes (baryon-photon tight-coupling to variable order, ultra-relativistic fluid, and reionisation-aware radiation streaming) for the CLASS Boltzmann solver. These are obtained from standard fluid and perturbation theory expansions with stated regimes of validity; none reduce by construction to fitted parameters drawn from the same CMB/LSS data the code is intended to analyze. The central performance claim is therefore independent of self-citation chains or internal redefinitions and remains externally falsifiable by direct numerical comparison.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard linearized cosmological perturbation theory and the validity of fluid approximations in well-defined epochs; no new free parameters are introduced beyond user choices of approximation order.

axioms (2)

domain assumption Linearized Einstein-Boltzmann equations remain valid for the scales and redshifts of interest in LambdaCDM
All approximations are applied inside the standard linear perturbation framework for basic LambdaCDM models.
standard math Fluid descriptions are accurate when mean free paths are short or particles are ultra-relativistic
The tight-coupling and ultra-relativistic fluid approximations are standard limiting cases of the Boltzmann hierarchy.

pith-pipeline@v0.9.0 · 5422 in / 1415 out tokens · 46402 ms · 2026-05-13T19:00:14.272950+00:00 · methodology

discussion (0)

Forward citations

Cited by 27 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Nonlinear Relativistic Effects on Cosmological Redshift Drift
astro-ph.CO 2026-04 unverdicted novelty 8.0

Second-order relativistic effects on redshift drift are computed, showing distortions appear only at this order with enhanced nonlinear bispectrum contributions at low redshift and large momenta.
Effective Field Theory of Large Scale Structure and Newtonian Motion Gauges
astro-ph.CO 2026-05 unverdicted novelty 7.0

Newtonian motion gauges extend the validity of Newtonian EFTofLSS to scale-dependent growth and GR effects by transforming linear equations to Newtonian form, computing nonlinear clustering there, and transforming res...
Cobaya: Code for Bayesian Analysis of hierarchical physical models
astro-ph.IM 2020-05 accept novelty 7.0

Cobaya is a modular Bayesian analysis code that exploits model interdependencies via automatic caching and a novel parameter-blocking algorithm to minimize sampling cost.
New Isocurvature Constraints from JWST UV Luminosity Function
astro-ph.CO 2026-05 unverdicted novelty 6.0

First UVLF-based constraints on model-agnostic isocurvature power spectra for CDM, baryon, neutrino, and dark radiation modes yield consistent 95% credible envelopes over k ~ 0.5-10 Mpc^{-1}.
CMB Limits on the Absorption of Light Vector and Axial-Vector Dark Matter
astro-ph.CO 2026-05 conditional novelty 6.0

Planck CMB data set upper limits on vector and axial-vector dark matter-electron couplings for masses 100 eV to 100 keV via energy injection from inelastic scattering and hydrogen absorption.
Post-Recombination Fluctuations from a Sequestered Dark Sector
astro-ph.CO 2026-05 unverdicted novelty 6.0

Develops a short-timescale injection formalism for post-recombination fluctuations, derives CMB impacts from scalar/vector/tensor perturbations, and constrains sequestered dark sector phase transitions to permille-lev...
Geometric Constraints on the Pre-Recombination Expansion History from the Hubble Tension
astro-ph.CO 2026-04 unverdicted novelty 6.0

Model-independent reconstruction shows that early-universe modifications resolving the Hubble tension exist at the background level, requiring a smooth ~15% pre-recombination expansion rate enhancement.
Double the axions, half the tension: multi-field early dark energy eases the Hubble tension
astro-ph.CO 2026-04 unverdicted novelty 6.0

Two-field axion-like early dark energy reduces Hubble tension to 1.5 sigma residual and improves high-ell CMB fits over single-field models.
Into the Gompverse: A robust Gompertzian reionization model for CMB analyses
astro-ph.CO 2026-04 unverdicted novelty 6.0

A Gompertzian reionization model with three nuisance parameters demotes optical depth to a derived quantity, reducing its uncertainty by a factor of three and revealing potential neutrino mass tension in CMB analyses.
A unified harmonic framework for dark siren cosmology
astro-ph.CO 2026-03 unverdicted novelty 6.0

The GW-galaxy cross-correlation method, unified with spectral sirens in a harmonic framework, can measure H0 to 1% and Omega_m to 5% precision with 2 years of data from next-generation detectors like Einstein Telescop...
Primordial Magnetic Fields at Cosmic Dawn: 21-cm Forecasts with HERA and SKA
astro-ph.CO 2026-05 unverdicted novelty 5.0

Primordial magnetic fields enhance low-mass halo abundance at Cosmic Dawn, advancing Lyman-alpha coupling, X-ray heating and reionization and imprinting detectable signatures on global and fluctuating 21-cm signals ob...
Alleviating the Hubble Tension Using $\Lambda$sCDM Model: A Coupled Dark Energy - Dark Matter Interaction
astro-ph.CO 2026-05 unverdicted novelty 5.0

The ΛsCDM model with coupled dark sectors reduces the Hubble tension to 1.2σ via late-time expansion changes while keeping the early-universe sound horizon nearly unchanged.
Cosmological constraints on TeV-scale dark matter subcomponents decaying between recombination and reionisation
astro-ph.CO 2026-04 conditional novelty 5.0

Future 21-cm observations may constrain TeV-scale decaying dark matter subcomponents more tightly than CMB data for lifetimes above 10^15 s, with strongest sensitivity for neutrino decay channels due to differences in...
A barotropic alternative to Early Dark Energy for alleviating the $H_0$ tension
astro-ph.CO 2026-04 unverdicted novelty 5.0

A barotropic fluid with ω_s ≈ 0.29 and Ω_s ≈ 1.5×10^{-5} raises the inferred H0 to match SH0ES while remaining consistent with Planck CMB, DESI BAO, and Pantheon data.
Cosmology of Inelastic Self-Interacting Dark Matter: Linear Evolution and Observational Constraints
astro-ph.CO 2026-04 unverdicted novelty 5.0

Inelastic self-interacting dark matter with small mass splitting produces a cutoff in the matter power spectrum at k > 1 h Mpc^{-1} whose location depends on cross-section normalization, velocity dependence, dark matt...
Reconstructing inflationary features on large scales using genetic algorithm
astro-ph.CO 2026-04 unverdicted novelty 5.0

Genetic algorithm reconstructs single-field inflationary models with features in the scalar power spectrum that fit Planck 2018 CMB data better by Δχ² ≲ -10 and suggest alternative background parameters.
Is the $w_0w_a$CDM cosmological parameterization evidence for dark energy dynamics partially caused by the excess smoothing of Planck PR4 CMB anisotropy data?
astro-ph.CO 2026-04 conditional novelty 5.0

Planck PR4 CMB data mildly favors dynamical dark energy, but this preference weakens when accounting for possible excess smoothing, indicating the signal may partly arise from data processing issues.
Exploring the interplay of late-time dynamical dark energy and new physics before recombination
astro-ph.CO 2026-03 unverdicted novelty 5.0

Model-independent reconstruction finds 96.7-98.5% probability of phantom crossing if recombination is standard, but early new physics to ease Hubble tension weakens this preference while requiring unrealistically high...
Planck 2018 results. VI. Cosmological parameters
astro-ph.CO 2018-07 accept novelty 5.0

Final Planck CMB data confirms the flat 6-parameter ΛCDM model with Ω_c h² = 0.120 ± 0.001, Ω_b h² = 0.0224 ± 0.0001, n_s = 0.965 ± 0.004, τ = 0.054 ± 0.007, H_0 = 67.4 ± 0.5 km/s/Mpc, and no strong evidence for extensions.
KiDS+VIKING-450 cosmology with Bayesian hierarchical model redshift distributions
astro-ph.CO 2026-05 conditional novelty 4.0

Bayesian hierarchical modeling of photometric redshifts in KiDS+VIKING-450 raises S8 to 0.756 ± 0.039 and reduces Planck tension to 1.9σ.
The Status of Gravitational Vector Perturbations with Recent CMB Data
astro-ph.CO 2026-05 unverdicted novelty 4.0

Recent CMB datasets tighten 95% CL upper bounds on vector-mode amplitude r_v to 1.3e-4 (neutrino isocurvature), 6.8 (octupole), and 4.2 (sourced) at k=0.05 Mpc^-1, with no significant detection.
Alleviating the Hubble Tension Using $\Lambda$sCDM Model: A Coupled Dark Energy - Dark Matter Interaction
astro-ph.CO 2026-05 unverdicted novelty 4.0

The ΛsCDM interacting dark sector model reduces the Hubble tension to 1.2σ via late-time energy transfer from dark matter to dark energy while leaving the sound horizon nearly unchanged.
In-depth analysis of the clustering of dark matter particles around primordial black holes. Part III: CMB constraints
astro-ph.CO 2026-04 unverdicted novelty 4.0

CMB data limits the s-wave annihilation cross section of thermal dark matter particles to ≲ 10^{-30} cm³/s scaled by PBH fraction and mass for PBHs heavier than ~10^{-10} solar masses.
Constraining Cosmological and Astrophysical Parameters with the Cosmic Star Formation History
astro-ph.CO 2026-04 unverdicted novelty 4.0

Cosmic star formation history provides complementary constraints on cosmological parameters, breaking degeneracies when combined with standard probes and yielding H0 = 68.28 ± 0.18 km s^{-1} Mpc^{-1} with DESI data.
Primordial magnetic fields in the light of upcoming post-EoR Lyman-$\alpha$ and 21-cm observations
astro-ph.CO 2026-04 unverdicted novelty 4.0

Forecasts show DESI-like and SKA1-Mid observations could constrain primordial magnetic field amplitude and spectral index to roughly 10% precision via 21cm auto-spectrum and Lyα-21cm cross-spectrum.
Mapping the redshift drift at various redshifts through cosmography
astro-ph.CO 2026-04 unverdicted novelty 4.0

Cosmographic Taylor and Padé models fitted to Pantheon+SH0ES+GRB+DESI BAO data yield redshift drift predictions compatible with ΛCDM and ω0ω1CDM at 1-2σ, with mock drift data tightening q0 and j0 bounds.
Extended Dark Energy analysis using DESI DR2 BAO measurements
astro-ph.CO 2025-03 conditional novelty 4.0

Extended analysis of DESI DR2 data confirms robust evidence for dynamical dark energy with phantom crossing preference, stable under parametric and non-parametric modeling.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages · cited by 26 Pith papers

[1]

Bertschinger, arXiv:astro-ph/9506070

E. Bertschinger, arXiv:astro-ph/9506070

work page arXiv
[2]

A Line of Sight Approach to Cosmic Microwave Background Anisotropies

U. Seljak, M. Zaldarriaga, Astrophys. J. 469 (1996) 437-444. [astro-ph/9603033]

work page Pith review arXiv 1996
[3]

Efficient Computation of CMB anisotropies in closed FRW models

A. Lewis, A. Challinor, A. Lasenby, Astrophys. J. 538 (2000) 473-476. [astro-ph/9911177]

work page Pith review arXiv 2000
[4]

Doran, JCAP 0510 (2005) 011

M. Doran, JCAP 0510 (2005) 011. [astro-ph/0302138]

work page arXiv 2005
[5]

Hamann, A

J. Hamann, A. Balbi, J. Lesgourgues, C. Quercellini, JCAP 0904 (2009) 011. [arXiv:0903.0382 [astro-ph.CO]]

work page arXiv 2009
[6]

The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview

J. Lesgourgues, “The Cosmic Linear Anisotropy Solving System ( CLASS) I: Overview,” [arXiv:1104.2932 [astro-ph.IM]]

work page Pith review arXiv
[7]

The Cosmic Linear Anisotropy Solving System (CLASS) IV: eﬃcient implementation of non-cold relics,

J. Lesgourgues, T. Tram, “The Cosmic Linear Anisotropy Solving System (CLASS) IV: eﬃcient implementation of non-cold relics,” [arXiv:1104.2935 [astro-ph.CO]]

work page arXiv
[8]

The Cosmic Linear Anisotropy Solving System ( CLASS) III: Comparison with CAMB for LambdaCDM,

J. Lesgourgues, “The Cosmic Linear Anisotropy Solving System ( CLASS) III: Comparison with CAMB for LambdaCDM,” [arXiv:1104.2934 [astro-ph.CO]]. 27

work page arXiv
[9]

P. J. E. Peebles and J. T. Yu, Astrophys. J. 162 (1970) 815

work page 1970
[10]

C. P. Ma and E. Bertschinger, Astrophys. J. 455 (1995) 7 [arXiv:astro-ph/9506072]

work page Pith review arXiv 1995
[11]

Doran, JCAP 0506 (2005) 011

M. Doran, JCAP 0506 (2005) 011. [astro-ph/0503277]

work page arXiv 2005
[12]

F. -Y. Cyr-Racine, K. Sigurdson, [arXiv:1012.0569 [astro-ph.CO ]]

work page arXiv
[13]

Hu, Astrophys

W. Hu, Astrophys. J. 506 (1998) 485-494. [astro-ph/9801234]

work page arXiv 1998
[14]

Pitrou, Phys

C. Pitrou, Phys. Lett. B 698 (2011) 1 [arXiv:1012.0546 [astro-ph.CO]]

work page arXiv 2011
[15]

S. Naoz, R. Barkana, Mon. Not. Roy. Astron. Soc. 362 (2005) 1047-1053. [astro-ph/0503196]

work page arXiv 2005
[16]

Lewis, Phys

A. Lewis, Phys. Rev. D76 (2007) 063001. [arXiv:0707.2727 [astro-ph]]

work page arXiv 2007
[17]

and Reichelt, Mark W

Shampine, Lawrence F. and Reichelt, Mark W. SIAM J. Sci. Compu t. vol. 18, 1 (1997)

work page 1997
[18]

Davis, Timothy A. SIAM J. Sci. Comput. vol. 18, 1 (1997) Direct M ethods for Sparse Linear Systems

work page 1997
[19]

Shoji, E

M. Shoji, E. Komatsu, Phys. Rev. D81 (2010) 123516. [arXiv:1003.0942 [astro-ph.CO]]

work page arXiv 2010
[20]

Bucher, K

M. Bucher, K. Moodley, N. Turok, Phys. Rev. D62 (2000) 083508. [astro-ph/9904231]. A. Stiﬀ integrator The standard numerical method for solving Ordinary Diﬀerent ial Equations (ODEs) is to use an adaptive step size Runge-Kutta solver. While this met hod is fast and accurate in simple cases, it may fail completely (or take a very large n umber of steps) wh...

work page arXiv 2000