The SRG/eROSITA all-sky survey: Constraints on Ultra-light Axion Dark Matter through Galaxy Cluster Number Counts
Pith reviewed 2026-05-23 03:47 UTC · model grok-4.3
The pith
Galaxy cluster counts from the eRASS1 survey exclude ultralight axion fractions above 0.0035 at masses near 10^{-27} eV.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using galaxy cluster number counts from the first eROSITA all-sky survey together with weak lensing data, the authors perform cosmological inference to place upper bounds on the ultralight axion relic density in logarithmic mass bins. They obtain the strongest constraints to date of Omega_a less than 0.0035 at m_a equals 10 to the minus 27 eV and less than 0.0079 at 10 to the minus 26 eV at 95 percent confidence level. These bounds tighten to 0.0030 and 0.0058 when CMB data are included. The work establishes cluster abundance as a probe for the ultralight axion fraction of dark matter for the first time.
What carries the argument
The response of the halo mass function to a nonzero ultralight axion fraction, used within a cosmological parameter inference pipeline on cluster counts and lensing data.
If this is right
- Improved theoretical modeling of halo abundances will allow future X-ray surveys to set even stronger limits on ultralight axions.
- The derived bounds restrict the viable parameter space for mixed dark matter models containing ultralight axions.
- Cluster number counts can now be added to the set of observables used to test fuzzy dark matter scenarios.
Where Pith is reading between the lines
- If these limits are confirmed, ultralight axions can constitute only a tiny fraction of dark matter at these masses, affecting models of structure formation on small scales.
- Extending the analysis to higher redshifts or combining with 21cm observations could test the same mass range through different physical effects.
Load-bearing premise
The halo mass function and its modification by ultralight axions are modeled accurately enough that errors do not shift the derived upper limits on the axion fraction.
What would settle it
Observing a higher number of galaxy clusters than allowed by the reported upper bounds on Omega_a in the relevant mass and redshift range would indicate that the exclusion is incorrect.
Figures
read the original abstract
Ultralight axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultralight axions constitute only a fraction of the dark matter has been limited to only a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam to constrain the fraction of ultralight axions in the mass range $10^{-32}$ eV to $10^{-24}$ eV. We put upper bounds on the ultralight axion relic density $\Omega_\mathrm{a}$ in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultralight axion mass regime with the tightest bounds reported so far in the mass bins around $m_\text{a} = 10^{-27}$ eV with $\Omega_\text{a} < 0.0035$ and $m_\text{a} = 10^{-26}$ eV with $\Omega_\text{a} < 0.0079$ (95% C.L.). When combined with cosmic microwave background probes, these bounds are tightened to $\Omega_\text{a} < 0.0030$ in the $m_\text{a} = 10^{-27}$ eV mass bin and $\Omega_\text{a} < 0.0058$ in the $m_\mathrm{a} = 10^{-26}$ eV mass bin (95% C.L.). This is the first time that constraints on ultralight axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultralight axions with better theoretical modeling of the abundance of halos.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses galaxy cluster number counts from the eRASS1 all-sky survey together with weak-lensing data from DES, KiDS and HSC to perform a full cosmological parameter inference and derive 95% CL upper limits on the ultralight axion relic density Ω_a in independent logarithmic mass bins spanning 10^{-32} eV to 10^{-24} eV. The tightest bounds are reported near m_a = 10^{-27} eV (Ω_a < 0.0035) and m_a = 10^{-26} eV (Ω_a < 0.0079), tightening further when combined with CMB data; the work claims to be the first application of cluster counts to fractional ULA constraints.
Significance. If the modeling assumptions hold, the result supplies the strongest existing limits on fractional ultralight axion dark matter in the intermediate-mass window and demonstrates that cluster abundance can serve as an independent growth-of-structure probe complementary to CMB and lensing. The explicit use of eRASS1 data for this purpose and the provision of binned limits constitute concrete advances.
major comments (3)
- [Methodology (HMF prescription)] Methodology section (halo-mass-function modeling): the predicted cluster abundance n(>M,z) that enters the likelihood is obtained from a mixed CDM+ULA halo mass function whose response to Ω_a ≲ 0.01 is taken from an unvalidated fitting function; because the reported 95% upper limits are driven by the difference between this prediction and the observed counts, any systematic bias in the suppression amplitude directly shifts the Ω_a bounds.
- [Abstract and methodology] Abstract and methodology: the manuscript itself states that the results 'pave the way … with better theoretical modeling of the abundance of halos,' confirming that the current HMF treatment is not demonstrated to be sub-dominant to statistical errors at the precision needed for the quoted limits.
- [Likelihood and parameter inference] Likelihood and covariance treatment: no information is supplied on the covariance matrix construction, mock validation, or marginalization over HMF nuisance parameters, all of which are required to establish that the Ω_a limits are robust rather than prior- or modeling-dominated.
minor comments (2)
- [Results] Notation for the mass bins and the precise definition of the logarithmic binning should be stated explicitly in the text and tables.
- [Figures] Figure captions should clarify whether the plotted limits include or exclude the CMB combination.
Simulated Author's Rebuttal
We thank the referee for their careful and constructive review of our manuscript. We address each major comment point by point below. Revisions have been made to improve the presentation of methodological details and to clarify limitations where appropriate.
read point-by-point responses
-
Referee: Methodology section (halo-mass-function modeling): the predicted cluster abundance n(>M,z) that enters the likelihood is obtained from a mixed CDM+ULA halo mass function whose response to Ω_a ≲ 0.01 is taken from an unvalidated fitting function; because the reported 95% upper limits are driven by the difference between this prediction and the observed counts, any systematic bias in the suppression amplitude directly shifts the Ω_a bounds.
Authors: We thank the referee for this important observation. The HMF response to ULAs is taken from the fitting function calibrated on N-body simulations in the referenced literature. We agree that independent validation at the precision needed for Ω_a ~ 0.01 would be desirable. In the revised manuscript we have expanded the methodology section with an explicit discussion of the fitting function's domain of applicability, added a quantitative estimate of the impact of plausible variations in the suppression amplitude on the final Ω_a limits, and included a dedicated paragraph on remaining theoretical uncertainties. We retain the original limits but now present them with this caveat more prominently. revision: partial
-
Referee: Abstract and methodology: the manuscript itself states that the results 'pave the way … with better theoretical modeling of the abundance of halos,' confirming that the current HMF treatment is not demonstrated to be sub-dominant to statistical errors at the precision needed for the quoted limits.
Authors: We acknowledge that the original abstract phrasing could be read as implying the HMF modeling uncertainties are not yet sub-dominant. We have revised the abstract to remove any ambiguity and added a new paragraph in Section 3 that quantifies the contribution of HMF modeling uncertainty relative to the statistical errors from the eRASS1 cluster counts. The revised text now states that, within the assumptions of the adopted fitting function, the modeling uncertainty is sub-dominant for the reported 95% CL bounds, while still noting that improved simulations will enable tighter future constraints. revision: yes
-
Referee: Likelihood and parameter inference: no information is supplied on the covariance matrix construction, mock validation, or marginalization over HMF nuisance parameters, all of which are required to establish that the Ω_a limits are robust rather than prior- or modeling-dominated.
Authors: We apologize for the omission. The revised manuscript now contains an expanded Section 4.2 that describes (i) the construction of the covariance matrix from a suite of mock catalogs, (ii) the mock-based validation of the likelihood pipeline, and (iii) the marginalization over the HMF nuisance parameters (including their priors). These additions demonstrate that the Ω_a upper limits remain stable under reasonable variations of the nuisance parameters and are not prior-dominated. revision: yes
Circularity Check
No significant circularity; limits derived from external data via standard inference
full rationale
The derivation uses observed eRASS1 cluster counts plus lensing data to place upper limits on Ω_a via full cosmological parameter inference in independent mass bins. The HMF response to fractional ULA is an input modeling choice whose accuracy is flagged by the authors as needing future improvement; it is not fitted to the present dataset and then re-used as a prediction. No self-definitional equations, fitted-input-as-prediction steps, or load-bearing self-citation chains appear in the abstract or described methodology. The central result remains an observational constraint rather than a tautology.
Axiom & Free-Parameter Ledger
free parameters (1)
- Ω_a per logarithmic mass bin
axioms (2)
- domain assumption The halo mass function and its dependence on ultralight axion fraction can be modeled accurately enough for cosmological inference
- standard math Standard cosmological parameter inference framework (flat priors, likelihood construction) applies without modification
invented entities (1)
-
ultralight axion
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We use axionCAMB … to compute a matter power spectrum … then compute the r.m.s density fluctuations … and use it as an input of the multiplicity function f(σ) from Tinker et al. (2008)
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the effect of ULAs on the HMF is illustrated … for different ULA masses ma and ULA relic densities Ωa
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 2 Pith papers
-
Artificial Precision Polarization Array: Sensitivity for the axion-like dark matter with clock satellites
Proposes satellite-based artificial pulsar polarization arrays (APPA) that simulations show can set tighter 95% C.L. upper limits on g_aγ than ground observations for axion masses 10^{-22} to 10^{-18} eV.
-
New constraints on cosmic anisotropy from galaxy clusters using an improved dipole fitting method
Galaxy cluster observations yield two preferred directions with cosmic anisotropy amplitude of about 5.3 times 10 to the minus 4 at roughly 1 sigma overall significance, though higher in the XMM-Newton subsample.
Reference graph
Works this paper leans on
-
[1]
Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2022, Phys. Rev. D, 105, 023520
work page 2022
-
[2]
Abbott, T. M. C., Allam, S., Andersen, P., et al. 2019, ApJ, 872, L30
work page 2019
- [3]
- [4]
- [5]
- [6]
-
[7]
2024a, arXiv e-prints, arXiv:2410.09499
Artis, E., Bulbul, E., Grandis, S., et al. 2024a, arXiv e-prints, arXiv:2410.09499
-
[8]
Asgari, M., Lin, C.-A., Joachimi, B., et al. 2021, A&A, 645, A104
work page 2021
- [9]
- [10]
- [11]
- [12]
-
[13]
E., Alonso, D., Krause, E., et al
Chisari, N. E., Alonso, D., Krause, E., et al. 2019, CCL: Core Cosmology Li- brary, Astrophysics Source Code Library, record ascl:1901.003
work page 2019
-
[14]
2024, arXiv e-prints, arXiv:2402.08457
Clerc, N., Comparat, J., Seppi, R., et al. 2024, arXiv e-prints, arXiv:2402.08457
-
[15]
2020, The Open Journal of As- trophysics, 3, 13
Comparat, J., Eckert, D., Finoguenov, A., et al. 2020, The Open Journal of As- trophysics, 3, 13
work page 2020
-
[16]
Conlon, J. P. 2006, Journal of High Energy Physics, 2006, 078
work page 2006
-
[17]
P., Day, F., Jennings, N., Krippendorf, S., & Rummel, M
Conlon, J. P., Day, F., Jennings, N., Krippendorf, S., & Rummel, M. 2017, J. Cosmology Astropart. Phys., 2017, 005
work page 2017
-
[18]
Costanzi, M., Saro, A., Bocquet, S., et al. 2021, Phys. Rev. D, 103, 043522
work page 2021
-
[19]
J., di Vecchia, P., Veneziano, G., & Witten, E
Crewther, R. J., di Vecchia, P., Veneziano, G., & Witten, E. 1979, Physics Letters B, 88, 123
work page 1979
- [20]
- [21]
-
[22]
Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10 de Jong, J. T. A., Verdoes Kleijn, G. A., Kuijken, K. H., & Valentijn, E. A. 2013, Experimental Astronomy, 35, 25 DES Collaboration, Abbott, T. M. C., Acevedo, M., et al. 2024, arXiv e-prints, arXiv:2401.02929 DESI Collaboration, Adame, A. G., Aguilar, J., et al. 2024, arXiv e-prints, ar...
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [23]
-
[24]
Du, X., Behrens, C., & Niemeyer, J. C. 2017, MNRAS, 465, 941
work page 2017
- [25]
-
[26]
Gatti, M., Sheldon, E., Amon, A., et al. 2021, MNRAS, 504, 4312
work page 2021
-
[27]
Gendler, N., Marsh, D. J. E., McAllister, L., & Moritz, J. 2024, J. Cosmology Astropart. Phys., 2024, 071
work page 2024
-
[28]
2024, arXiv e-prints, arXiv:2402.08458
Ghirardini, V ., Bulbul, E., Artis, E., et al. 2024, arXiv e-prints, arXiv:2402.08458
-
[29]
Giblin, B., Heymans, C., Asgari, M., et al. 2021, A&A, 645, A105
work page 2021
-
[30]
Grandis, S., Bocquet, S., Mohr, J. J., Klein, M., & Dolag, K. 2021, MNRAS, 507, 5671
work page 2021
-
[31]
2024, arXiv e-prints, arXiv:2402.08455
Grandis, S., Ghirardini, V ., Bocquet, S., et al. 2024, arXiv e-prints, arXiv:2402.08455
-
[32]
Grin, D., Marsh, D. J. E., & Hlozek, R. 2022, axionCAMB: Modification of the CAMB Boltzmann code, Astrophysics Source Code Library, record ascl:2203.026
work page 2022
-
[33]
Hildebrandt, H., van den Busch, J. L., Wright, A. H., et al. 2021, A&A, 647, A124
work page 2021
-
[34]
2013, ApJS, 208, 19 Hložek, R., Grin, D., Marsh, D
Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 Hložek, R., Grin, D., Marsh, D. J. E., & Ferreira, P. G. 2015, Phys. Rev. D, 91, 103512 Hložek, R., Marsh, D. J. E., & Grin, D. 2018, MNRAS, 476, 3063
work page 2013
- [35]
- [36]
-
[37]
Hui, L., Ostriker, J. P., Tremaine, S., & Witten, E. 2017, Phys. Rev. D, 95, 043541
work page 2017
-
[38]
Kim, J. E. & Carosi, G. 2010, Reviews of Modern Physics, 82, 557
work page 2010
-
[39]
2024, arXiv e-prints, arXiv:2402.08456
Kleinebreil, F., Grandis, S., Schrabback, T., et al. 2024, arXiv e-prints, arXiv:2402.08456
- [40]
-
[41]
Kobayashi, T., Murgia, R., De Simone, A., Irši ˇc, V ., & Viel, M. 2017, Phys. Rev. D, 96, 123514
work page 2017
- [42]
-
[43]
Kulkarni, M. & Ostriker, J. P. 2022, MNRAS, 510, 1425 Laguë, A., Bond, J. R., Hložek, R., et al. 2022, J. Cosmology Astropart. Phys., 2022, 049
work page 2022
-
[44]
F., Marulli, F., Moscardini, L., et al
Lesci, G. F., Marulli, F., Moscardini, L., et al. 2022, A&A, 659, A88
work page 2022
-
[45]
Lewis, A. & Challinor, A. 2011, CAMB: Code for Anisotropies in the Microwave
work page 2011
- [46]
- [47]
-
[48]
Marsh, D. J. E. 2016, Phys. Rep., 643, 1
work page 2016
-
[49]
Marsh, D. J. E., Macaulay, E., Trebitsch, M., & Ferreira, P. G. 2012, Phys. Rev. D, 85, 103514
work page 2012
- [50]
- [51]
-
[52]
Miyatake, H., Sugiyama, S., Takada, M., et al. 2023, Phys. Rev. D, 108, 123517
work page 2023
- [53]
- [54]
-
[55]
Passaglia, S. & Hu, W. 2022, Phys. Rev. D, 105, 123529
work page 2022
-
[56]
Peccei, R. D. & Quinn, H. R. 1977b, Phys. Rev. Lett., 38, 1440 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016a, A&A, 594, A13 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016b, A&A, 594, A24 Planck Collaboration, Aghanim, N., Akrami, Y ., et al. 2020, A&A, 641, A6
work page 2020
-
[57]
Predehl, P., Andritschke, R., Arefiev, V ., et al. 2021, A&A, 647, A1
work page 2021
-
[58]
Ragagnin, A., Saro, A., Singh, P., & Dolag, K. 2021, MNRAS, 500, 5056
work page 2021
-
[59]
W., Markevitch, M., Clowe, D., Gonzalez, A
Randall, S. W., Markevitch, M., Clowe, D., Gonzalez, A. H., & Bradaˇc, M. 2008, ApJ, 679, 1173
work page 2008
-
[60]
Reynolds, C. S., Marsh, M. C. D., Russell, H. R., et al. 2020, ApJ, 890, 59
work page 2020
-
[61]
K., Hložek, R., Laguë, A., et al
Rogers, K. K., Hložek, R., Laguë, A., et al. 2023, J. Cosmology Astropart. Phys., 2023, 023
work page 2023
-
[62]
Rogers, K. K. & Peiris, H. V . 2021, Phys. Rev. Lett., 126, 071302
work page 2021
-
[63]
Schive, H.-Y ., Chiueh, T., & Broadhurst, T. 2014, Nature Physics, 10, 496
work page 2014
-
[64]
Schwabe, B., Gosenca, M., Behrens, C., Niemeyer, J. C., & Easther, R. 2020, Phys. Rev. D, 102, 083518
work page 2020
- [65]
- [66]
-
[67]
Sevilla-Noarbe, I., Bechtol, K., Carrasco Kind, M., et al. 2021, ApJS, 254, 24
work page 2021
-
[68]
Stott, M. J. & Marsh, D. J. E. 2018, Phys. Rev. D, 98, 083006
work page 2018
-
[69]
Sunyaev, R., Arefiev, V ., Babyshkin, V ., et al. 2021, A&A, 656, A132
work page 2021
- [70]
-
[71]
Tinker, J., Kravtsov, A. V ., Klypin, A., et al. 2008, ApJ, 688, 709 V ogt, S. M. L., Marsh, D. J. E., & Laguë, A. 2023, Phys. Rev. D, 107, 063526
work page 2008
-
[72]
Winch, H., Rogers, K. K., Hložek, R., & Marsh, D. J. E. 2024, arXiv e-prints, arXiv:2404.11071
-
[73]
H., Hildebrandt, H., van den Busch, J
Wright, A. H., Hildebrandt, H., van den Busch, J. L., & Heymans, C. 2020, A&A, 637, A100
work page 2020
-
[74]
2022, MNRAS, 511, 5492 Article number, page 16 of 16
Zhao, C., Variu, A., He, M., et al. 2022, MNRAS, 511, 5492 Article number, page 16 of 16
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.