Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

A Fast and Accurate Implementation of the Effective Fluid Approximation for Ultralight Axions

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A fast CAMB implementation of the Passaglia–Hu effective fluid approximation claims sub-percent CMB accuracy for ultralight axions and reports tighter dark-matter limits.

desk verdict Useful CAMB implementation of the PH EFA with honest new constraints, but the headline accuracy claim rests on a self-referential validation. read the letter →

arxiv 2501.13662 v1 pith:TCAFIH4O submitted 2025-01-23 astro-ph.CO hep-ph

Adam Moss , Lauren Gaughan , Anne M. Green This is my paper · ORCID
classification astro-ph.COhep-ph
keywords ultralightaxionseffectivefluidapproximationPassaglia-HuCMBpowerspectrumKlein-GordonequationcosmologicalparameterconstraintsaxiondarkmatterAxiCAMB
topics Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to establish that ultralight axions with masses around $10^{-28}$ to $10^{-24}$ eV can be evolved accurately and quickly by replacing the rapidly oscillating scalar field with a calibrated effective fluid. It presents an implementation of the Passaglia–Hu effective fluid approximation inside the Boltzmann code CAMB, in which two slowly varying envelope fields carry the axion dynamics instead of the full Klein–Gordon oscillation. The central claim is that this method matches the CMB temperature power spectrum of the exact Klein–Gordon equation to sub-percent accuracy across that mass range, while allowing an early and cheap switch from field to fluid evolution. Using Planck PR4 CMB data and DESI BAO data, the paper reports $2\sigma$ upper limits $f_{\rm ax} < 0.0082$ and $\Omega_{\rm ax} h^2 < 0.0010$ at $m = 10^{-28}$ eV, about 50% tighter than earlier standard-EFA analyses. If correct, this removes a known systematic from axion constraints at precisely the masses where current and near-future surveys are most sensitive.

What carries the argument

The load-bearing object is the Passaglia–Hu decomposition of the axion field into two slowly varying auxiliary fields, $\phi(\tau) = \phi_c(\tau)\cos(\tau-\tau_*) + \phi_s(\tau)\sin(\tau-\tau_*)$ with $\tau = mt$, which filters the fast oscillations out of the equations of motion. From these envelope fields the paper constructs an effective fluid energy density, pressure, and heat flux, and closes the system with a calibrated equation of state $w_{\rm ax} \simeq \frac{3}{2}(m/H)^{-2}$ and a sound speed that carries a $(m/H)^{-2}$ correction. Matching conditions at the switch time suppress spurious oscillatory modes in both the background and perturbations. This machinery is what lets the switch from Klein–Gordon to fluid evolution happen at $m/H_* \approx 50$ with far less sensitivity to the switch parameter than the standard effective fluid approximation.

What would settle it

Compute the CMB TT spectrum with the production settings ($m/H_* = 50$, accuracy=1) and compare it with a direct numerical integration of the perturbed Klein–Gordon equation at $m = 5 \times 10^{-27}$ eV; if the relative difference at $\ell \sim 2000$ exceeds roughly one percent, the sub-percent claim is refuted.

Watch

Extended reading notes

Core claim

The paper's contribution is a concrete numerical claim: the Passaglia–Hu effective fluid approximation, coded into CAMB as AxiCAMB, reproduces the observable cosmological effects of ultralight axions to sub-percent accuracy in the CMB temperature power spectrum for $m \in [10^{-28}, 10^{-24}]$ eV. The method decomposes the axion field into amplitude and phase envelope fields whose equations of motion are free of the fast oscillations, then maps them onto an effective fluid with an equation of state and sound speed calibrated to exact Klein–Gordon solutions. The validation in the body compares against a late-switch, high-accuracy version of the same scheme as a ground-truth reference, and finds that the production setting $m/H_* \approx 50$ with default integrator accuracy stays close to that reference across the mass range, whereas the standard effective fluid approximation requires much later switches. The resulting MCMC analysis yields $2\sigma$ upper limits $f_{\rm ax} < 0.0082$ and $\Omega_{\rm ax} h^2 < 0.0010$ at $m = 10^{-28}$ eV, with comparable limits at higher masses improving on earlier standard-EFA results.

Load-bearing premise

The sub-percent accuracy claim rests on treating a late-switch run of the same method at $m/H_* = 300$ with accuracy=3 as ground truth, because the paper does not compare directly against exact Klein–Gordon solutions in the body.

Editorial extensions

If this is right

  • Axion constraints in the mass range $10^{-28}$ to $10^{-24}$ eV can be computed without the switch-time systematic that affects the standard effective fluid approximation.
  • Production MCMC runs can use an early switch with default integrator settings, making percent-level axion limits computationally affordable.
  • For Planck PR4 plus DESI BAO data, the method gives $f_{\rm ax} < 0.0082$ and $\Omega_{\rm ax} h^2 < 0.0010$ at $m = 10^{-28}$ eV.
  • The implementation is compatible with the current Boltzmann-code base, so it can be combined with new CMB and large-scale structure likelihoods as they appear.

Reading between the lines

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

  • The same auxiliary-field construction is the natural starting point for the full cosine potential and self-interacting axion-like fields; an extension there would let 'extreme axion' scenarios be sampled with the same speed.
  • Because the central accuracy claim is checked against a late-switch version of the same approximation rather than an independent exact solver, a direct comparison with a true Klein–Gordon integrator would independently test the sub-percent statement.
  • If the switch-time robustness extends above $10^{-24}$ eV, the method could sharpen forecasts for CMB-S4 and intensity mapping experiments, whose axion sensitivity moves to smaller scales.
  • The tighter limits combine a more accurate fluid model with newer datasets, so a fixed-data comparison between standard and PH EFA would separate the method's contribution from the data's contribution.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper describes AxiCAMB, a public implementation of the Passaglia–Hu effective fluid approximation (PH EFA) for ultralight axions in the CAMB Boltzmann code. The authors derive the background and perturbation equations, detail the numerical switching between Klein–Gordon evolution and the fluid description, and show reduced sensitivity to the switch time compared to the standard EFA. They then perform MCMC analyses with Planck PR4 and DESI BAO data to obtain upper limits on the axion fraction and physical density for fixed masses from 10^-28 to 10^-24 eV, reporting improved constraints relative to previous standard-EFA analyses. The abstract claims sub-percent accuracy relative to exact Klein–Gordon solutions, and that claim is the focus of my concerns.

Significance. If the sub-percent accuracy claim is fully supported, AxiCAMB would be a valuable public tool for fast and accurate ULA parameter estimation, and the demonstration of switch-time robustness is a useful contribution. The code is open source, and the MCMC pipeline is described in sufficient detail for reproduction. However, the headline accuracy claim is not directly validated against exact Klein–Gordon solutions in the body of the paper; the validation instead compares the PH EFA with an early switch to the same approximation with a late switch. This leaves the cosmological constraints vulnerable to systematic errors common to the approximation, so the significance of the results depends on closing that validation gap.

major comments (2)
  1. [III.B and III.C] The paper uses the phrase 'ground truth' inconsistently. In Section III.B, the m/H*=1000 case is called 'a ground-truth benchmark', while Section III.C adopts the PH EFA with m/H*=300 and accuracy=3 as 'a ground-truth reference'. Neither is an exact solution; both are the same approximation at different switch times. This is not merely a wording issue: it means the central accuracy claim is not tied to an external standard, and the choice of reference switch time can itself affect the inferred chi^2. The authors should state explicitly what the reference is, and justify why m/H*=300 is sufficient to serve as a proxy for the exact solution (e.g., by showing convergence of the CMB spectrum with increasing m/H* for the PH EFA).
  2. [III.A, Eq. (30)] The weighting function defined in Eq. (30), W(a)=exp(alpha(ln a* - ln a)) with alpha>0, is not bounded above; for a << a*, W(a) is larger than 1 and the combination rho_ax = W rho_KG + (1-W)rho_fluid involves a negative weight on the fluid density. The accompanying text states that 'For a << a*, the evolution is dominated by the exact field solution', which is not what this formula implies unless rho_fluid is negligible in that regime. Please clarify the domain of validity of Eq. (30) and either use a bounded weighting (e.g., a sigmoid or a normalized average) or justify the unbounded form and show that it does not introduce artifacts in the background and perturbation evolution.
minor comments (4)
  1. [III.B] There is a typo in the sentence 'the the full KG evolution' (should be 'the full KG evolution'), and 'reyonization' should be 'reionization' in the same section.
  2. [Fig. 6 caption] The caption reads 'As for Fig. 6 but for a switch m/H*=10.' It should refer to Fig. 5, not Fig. 6.
  3. [IV] The convergence criterion states 'we require R1 <= 0.05', which is inconsistent with the standard Gelman-Rubin statistic; this should likely read 'R-1 <= 0.05' (or equivalently R <= 1.05). Please clarify to ensure reproducibility.
  4. [Throughout] There are several small textual issues: 'max' is used instead of 'm_ax' in the Introduction and elsewhere, and the name 'Ureña-López' is rendered incorrectly in the bibliography. These should be corrected in a final polish.

Circularity Check

1 steps flagged · score 6.0 of 10

Accuracy validation is self-referential: the 'ground-truth' reference is the same PH EFA at a later switch, so the abstract's sub-percent claim versus exact Klein-Gordon is not actually tested in the body.

  1. self definitional [Section III.C, Eq. (31), Figure 7; also Section III.B, Figures 3-4]
    "As a ground-truth reference, we adopt the PH EFA with a switch at m/H∗ = 300 and accuracy=3. For earlier switch times and axion masses in the range 10−28 eV to 10−25 eV, we compute the resulting CMB temperature power spectra, C (model) ℓ , using accuracy=1. These are compared to the ground-truth spectra, C (ref) ℓ , via the mismatch measure: χ2 = ..."

    The reference spectrum C(ref)_ell is produced by the same PH EFA code at a later switch time, not by integrating the exact Klein-Gordon equation. Therefore Eq. (31) measures the switch-time convergence of PH EFA against itself; any systematic error common to all PH EFA switch choices (e.g., the approximate equation of state in Eq. 16, the sound speed in Eq. 26, or the matching and weighting implementation) is invisible to this test. The abstract's claim 'Compared to exact solutions of the Klein-Gordon equation, our method achieves sub-percent accuracy' is not established by this comparison, and no exact-KG CMB spectra are shown in the body.

full rationale

The circularity is partial and localized to the validation of the accuracy claim. The fluid equations themselves (the Passaglia-Hu EFA, including Eqs. 16 and 26 and the matching conditions) are imported from Ref. [23], which the paper states was calibrated against exact solutions of the Klein-Gordon equation; that prior work provides independent, parameter-free support and is not authored by the present authors, so the self-citation and uniqueness-import patterns do not apply. The MCMC constraints are data-driven fits to Planck PR4 and DESI BAO data, not derived from the approximation, and are therefore not circular. However, the paper's own headline claim of 'sub-percent accuracy' compared to exact Klein-Gordon solutions is not tested in the body: the designated ground truth is the PH EFA itself at a later switch time, so the central accuracy validation reduces by construction to a self-consistency check against the same approximation family. Errors common to the PH EFA at all switch times would not be detected. This warrants a partial-circularity score of 6 rather than a higher score, because the underlying method has external grounding in the cited Passaglia-Hu calibration and the implementation is publicly available for independent checking.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities. The auxiliary fields phi_c and phi_s are mathematical variables, not new particles. The main free parameters are numerical choices (switch time, weighting, accuracy). The central claim rests on the prior PH EFA calibration and on the unverified assumption that the late-switch PH EFA approximates exact KG solutions.

free parameters (3)
  • Switch threshold m/H* = 50 (used in MCMC)
    Chosen as a practical compromise between accuracy and speed; not fitted to data but is a user-defined numerical parameter.
  • Weighting sharpness alpha = 100
    Controls the smoothness of the KG-to-fluid transition in Eq. (30); chosen by hand.
  • CAMB accuracy setting = 1 (default) for MCMC, 3 for ground truth
    Integrator accuracy settings are numerical choices that affect run time and precision.
assumptions (5)
  • domain assumption The axion potential is approximated as quadratic, V(phi) ≈ 1/2 m^2 phi^2.
    Stated in Section II.A; assumes the field is near the minimum of the cosine potential, excluding 'extreme axion' scenarios with large initial field values.
  • domain assumption The Passaglia-Hu effective fluid approximation, including the calibrated sound speed and equation of state, is correct.
    The paper uses Eqs. (16) and (26) from Ref [23] without re-deriving them; the accuracy of this paper depends on that prior calibration.
  • ad hoc to paper PH EFA with switch m/H* = 300 and accuracy=3 serves as a valid ground-truth reference.
    In Section III.C this is adopted as the reference for computing chi^2; the body does not independently verify convergence to exact Klein-Gordon solutions.
  • standard math Standard synchronous-gauge perturbation equations and FLRW background are valid.
    Standard cosmological perturbation theory, used throughout Sections II and III.
  • domain assumption Neglecting the effect of modified dH/dz on recombination is acceptable.
    Section III.B states that empirically the effect on CMB anisotropies is negligible for the masses and abundances considered, verified by modifying recfast.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Fast and Accurate Implementation of the Effective Fluid Approximation for Ultralight Axions." pith.science (2026). https://pith.science/paper/TCAFIH4O

@misc{pith2026250113662,
  author       = {Pith},
  title        = {Pith review of: A Fast and Accurate Implementation of the Effective Fluid Approximation for Ultralight Axions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TCAFIH4O}},
  note         = {Machine review of arXiv:2501.13662}
}
abstract

We present a numerically efficient and accurate implementation of the Passaglia-Hu effective fluid approximation for ultralight axions (ULAs) within the Boltzmann code CAMB. This method is specifically designed to evolve the axion field accurately across cosmological timescales, mitigating the challenges associated with its rapid oscillations. Our implementation is based on the latest version of CAMB, ensuring compatibility with other cosmological codes., e.g. for calculating cosmological parameter constraints. Compared to exact solutions of the Klein-Gordon equation, our method achieves sub-percent accuracy in the CMB power spectrum across a broad range of axion masses, from $10^{-28}\,\mathrm{eV}$ to $10^{-24}\,\mathrm{eV}$. We perform Markov Chain Monte Carlo (MCMC) analyses incorporating our implementation, and find improved constraints on the axion mass and abundance compared to previous, simpler fluid-based approximations. For example, using \Planck\ PR4 and DESI BAO data, we find $2\sigma$ upper limits on the axion fraction $f_{\rm ax} < 0.0082$ and physical density $\Omega_{\rm ax}h^2 < 0.0010$ for $m=10^{-28}$ eV. The code is publicly available at \url{https://github.com/adammoss/AxiCAMB}.

Figures

Figures reproduced from arXiv: 2501.13662 by the authors.

Figure 1
Figure 1. FIG. 1. Violin plots of the posterior distributions for [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Evolution of the axion energy density around the time [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The axion density, [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. As for Fig. 3 but for a switch [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. The CMB temperature power spectra ( [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. As for Fig. 6 but for a switch [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Comparison of the mismatch measure [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Inflationary Axion Isocurvature in the CMB across All Ultralight Masses

    astro-ph.CO 2026-08 conditional novelty 6.0 of 10

    Accurate axion-isocurvature evolution across the full ultralight mass range converts Planck and BICEP bounds into inflation-scale limits and reveals a small discovery window near m_a ~ 10^-25 eV.

Reference graph

Works this paper leans on

30 extracted references · 2 canonical work pages · cited by 1 Pith paper

  1. [30]

    R. Liu, W. Hu, and D. Grin, (2024), arXiv:2412.15192 [astro-ph.CO]

    arXiv 2024

  2. [15]

    K. K. Rogers, R. Hloˇ zek, A. Lagu¨ e, M. M. Ivanov, O. H. E. Philcox, G. Cabass, K. Akitsu, and D. J. E. Marsh, JCAP 06, 023, arXiv:2301.08361 [astro-ph.CO]

    arXiv

  3. [1]

    Aghanim et al

    N. Aghanim et al. (Planck), Astron. Astrophys. 641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]

    arXiv 2020

  4. [2]

    Cirelli, A

    M. Cirelli, A. Strumia, and J. Zupan, (2024), arXiv:2406.01705 [hep-ph]

    arXiv 2024

  5. [3]

    E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006), arXiv:hep-th/0603057

    arXiv 2006

  6. [4]

    Svrcek and E

    P. Svrcek and E. Witten, JHEP 06, 051, arXiv:hep- th/0605206

    arXiv

  7. [5]

    Arvanitaki, S

    A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, and J. March-Russell, Phys. Rev. D 81, 123530 (2010), arXiv:0905.4720 [hep-th]

    arXiv 2010

  8. [6]

    D. J. E. Marsh, Phys. Rept. 643, 1 (2016), arXiv:1510.07633 [astro-ph.CO]

    arXiv 2016

Show all 30 references
  1. [7]

    W. Hu, R. Barkana, and A. Gruzinov, Phys. Rev. Lett. 85, 1158 (2000), arXiv:astro-ph/0003365

    2000 arXiv

  2. [8]

    J. A. Frieman, C. T. Hill, A. Stebbins, and I. Waga, Phys. Rev. Lett. 75, 2077 (1995), arXiv:astro-ph/9505060

    1995 arXiv

  3. [9]

    R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977)

    1977

  4. [10]

    Weinberg, Phys

    S. Weinberg, Phys. Rev. Lett. 40, 223 (1978)

    1978

  5. [11]

    Wilczek, Phys

    F. Wilczek, Phys. Rev. Lett. 40, 279 (1978)

    1978

  6. [12]

    Hlozek, D

    R. Hlozek, D. Grin, D. J. E. Marsh, and P. G. Ferreira, Phys. Rev. D 91, 103512 (2015), arXiv:1410.2896 [astro- ph.CO]

    2015 arXiv

  7. [13]

    Rosenberg, S

    E. Rosenberg, S. Gratton, and G. Efstathiou, Mon. Not. Roy. Astron. Soc. 517, 4620 (2022), arXiv:2205.10869 [astro-ph.CO]

    2022 arXiv

  8. [14]

    A. G. Adame et al. (DESI), (2024), arXiv:2404.03002 [astro-ph.CO]

    2024 arXiv

  9. [16]

    Abdalla et al., JHEAp 34, 49 (2022), arXiv:2203.06142 [astro-ph.CO]

    E. Abdalla et al., JHEAp 34, 49 (2022), arXiv:2203.06142 [astro-ph.CO]

    2022 arXiv

  10. [17]

    Cookmeyer, D

    T. Cookmeyer, D. Grin, and T. L. Smith, Phys. Rev. D 101, 023501 (2020), arXiv:1909.11094 [astro-ph.CO]

    2020 arXiv

  11. [18]

    L. A. Ure˜ na L´ opez and A. X. Gonzalez-Morales, JCAP 07, 048, arXiv:1511.08195 [astro-ph.CO]

    arXiv

  12. [19]

    Hlozek, D

    R. Hlozek, D. J. E. Marsh, and D. Grin, Mon. Not. Roy. Astron. Soc. 476, 3063 (2018), arXiv:1708.05681 [astro- ph.CO]

    2018 arXiv

  13. [20]

    Poulin, T

    V. Poulin, T. L. Smith, D. Grin, T. Karwal, and M. Kamionkowski, Phys. Rev. D 98, 083525 (2018), arXiv:1806.10608 [astro-ph.CO]

    2018 arXiv

  14. [21]

    Winch, R

    H. Winch, R. Hlozek, D. J. E. Marsh, D. Grin, and K. K. Rogers, Phys. Rev. D 110, 043517 (2024), arXiv:2311.02052 [astro-ph.CO]

    2024 arXiv

  15. [22]

    Abazajian et al., (2019), arXiv:1907.04473 [astro- ph.IM]

    K. Abazajian et al., (2019), arXiv:1907.04473 [astro- ph.IM]

    2019 arXiv

  16. [23]

    Passaglia and W

    S. Passaglia and W. Hu, Phys. Rev. D105, 123529 (2022), arXiv:2201.10238 [astro-ph.CO]

    2022 arXiv

  17. [24]

    Lewis, A

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

    2000 arXiv

  18. [25]

    Seager, D

    S. Seager, D. D. Sasselov, and D. Scott, Astrophys. J. Lett. 523, L1 (1999), arXiv:astro-ph/9909275. 12

    1999 arXiv

  19. [26]

    Scott and A

    D. Scott and A. Moss, Mon. Not. Roy. Astron. Soc. 397, 445 (2009), arXiv:0902.3438 [astro-ph.CO]

    2009 arXiv

  20. [27]

    Torrado and A

    J. Torrado and A. Lewis, JCAP05, 057, arXiv:2005.05290 [astro-ph.IM]

    2005 arXiv

  21. [28]

    Hloˇ zek, D

    R. Hloˇ zek, D. J. E. Marsh, D. Grin, R. Allison, J. Dunk- ley, and E. Calabrese, Phys. Rev. D 95, 123511 (2017), arXiv:1607.08208 [astro-ph.CO]

    2017 arXiv

  22. [29]

    J. B. Bauer, D. J. E. Marsh, R. Hloˇ zek, H. Padmanabhan, and A. Lagu¨ e, Mon. Not. Roy. Astron. Soc. 500, 3162 (2020), arXiv:2003.09655 [astro-ph.CO]

    2020 arXiv

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.