The impact of strong feedback on galaxy group scaling relations
Pith reviewed 2026-05-21 17:52 UTC · model grok-4.3
The pith
Highly ejective feedback models underpredict the X-ray luminosity of galaxy groups at fixed mass.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Highly ejective feedback models under-predict the X-ray luminosity of galaxy groups in the mass range 10^13 to 10^14 solar masses at 5.7 sigma significance. This mismatch is found in a sample of 44 groups with high-quality XMM-Newton data and remains after selection effects are considered, using quantities that are directly measurable and minimally correlated.
What carries the argument
The X-ray luminosity at fixed halo mass for galaxy groups, used to test feedback models calibrated on baryon fractions.
If this is right
- Feedback models in cosmological simulations should be calibrated on observable scaling relations rather than on derived baryon fractions.
- The efficiency of gas ejection by active galactic nuclei in groups may need to be lowered to reproduce local X-ray properties.
- Stacking heterogeneous galaxy group samples to estimate average gas fractions can introduce large systematic errors.
- Tension between different baryon fraction constraints, such as those from kinetic Sunyaev-Zel'dovich measurements and optical group catalogs, requires resolution.
Where Pith is reading between the lines
- If the claim holds, large-volume simulations may require mass-dependent adjustments to feedback prescriptions to avoid underpredicting hot gas in groups.
- Upcoming wide-field X-ray surveys could enlarge the sample size and test whether the same luminosity deficit appears at higher redshifts or in different environments.
- The result underscores the value of comparing simulations directly to raw observables before converting them into physical quantities like gas mass.
Load-bearing premise
The mismatch between observed X-ray luminosities and simulation predictions arises chiefly from the ejective strength of feedback rather than from unmodeled astrophysics, measurement biases, or how representative the 44-group sample is.
What would settle it
Precise X-ray luminosity and mass measurements for a larger, volume-limited sample of galaxy groups that instead match the luminosities predicted by highly ejective feedback simulations.
Figures
read the original abstract
Feedback from active supermassive black holes alters the distribution of matter in the Universe by injecting energy in the neighbouring hot gaseous medium, which leads to ejection of gas from the halos of galaxy groups and massive galaxies. Recent cosmological simulations such as FLAMINGO calibrate their feedback model on the baryon fractions of galaxy groups to tune the efficiency of gas ejection. However, recent observational constraints from optically selected groups and the kinetic Sunyaev-Zel'dovich effect yield lower baryon fractions than previous studies, which indicates that feedback may be more ejective than previously thought. Here we show that models involving highly ejective feedback are inconsistent with the scaling relations of local galaxy groups in the mass range $10^{13}-10^{14}M_\odot$. We study the X-ray luminosity-temperature relation in a sample of 44 galaxy groups with high-quality XMM-Newton observations. We show that highly ejective models under-predict the luminosity of galaxy groups at fixed mass at high significance ($5.7\sigma$). This conclusion is robust against selection effects and is obtained from directly measurable and minimally correlated quantities. We point out that turning observable quantities into gas fraction estimates is challenging, especially in the context of stacking large samples of heterogeneous systems. We argue that calibrating feedback models on baryon fractions is prone to systematic uncertainties and that observable scaling relations are better suited for this task.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the X-ray luminosity-temperature scaling relation for a sample of 44 local galaxy groups with high-quality XMM-Newton data in the halo mass range 10^13-10^14 M_⊙. It concludes that cosmological simulations with highly ejective AGN feedback (e.g., FLAMINGO variants calibrated on baryon fractions) under-predict observed group X-ray luminosities at fixed mass at 5.7σ significance. The authors argue this discrepancy is robust to selection effects, relies on directly measurable and minimally correlated quantities, and implies that calibrating feedback models on baryon fractions is prone to systematics while observable scaling relations are preferable.
Significance. If the central result holds, the work would meaningfully constrain feedback prescriptions in large-volume cosmological simulations by showing that highly ejective models fail to reproduce observed X-ray scaling relations in the group regime. The emphasis on directly observable quantities rather than derived gas fractions is a constructive contribution, as it sidesteps some stacking and conversion uncertainties highlighted in recent kSZ and optical-group studies. This could shift calibration practices away from baryon-fraction tuning toward scaling-relation matching.
major comments (2)
- [Abstract and §4] Abstract and §4 (results on Lx at fixed mass): The 5.7σ discrepancy is presented as arising from directly measurable, minimally correlated quantities, yet the manuscript does not quantify how residual covariance between the X-ray luminosity and the mass proxy (e.g., via hydrostatic equilibrium assumptions or gas clumping in the 10^13-10^14 M_⊙ regime) propagates into the significance. A bias in M that correlates with Lx at even the 10-20% level could shift the comparison horizontally and reduce the tension below 5σ; an explicit error-budget table or Monte-Carlo test of this effect is needed to support the quoted significance.
- [§3] §3 (sample and mass estimation): The claim of robustness to selection effects is stated, but the text does not demonstrate that the 44-group XMM-Newton sample is representative of the underlying population at fixed true mass rather than at fixed observed Lx or T. If the selection function correlates with the same gas properties that set Lx, the under-prediction by ejective models could be partly selection-driven; a completeness simulation or comparison to a volume-limited catalog would strengthen the result.
minor comments (2)
- [Figure 2] Figure 2 or equivalent (Lx-T or Lx-M plot): Axis labels and error bars should explicitly indicate whether masses are derived from X-ray hydrostatic equilibrium or from an independent proxy; the current presentation leaves the degree of correlation ambiguous.
- [§2] §2 (simulation comparison): The specific FLAMINGO feedback variants used (e.g., the exact AGN heating temperature or ejection efficiency parameter) should be tabulated with their baryon-fraction calibration values for direct reproducibility.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments. These have prompted us to strengthen the presentation of our error analysis and selection robustness. We respond point-by-point below and indicate the revisions we will make.
read point-by-point responses
-
Referee: [Abstract and §4] Abstract and §4 (results on Lx at fixed mass): The 5.7σ discrepancy is presented as arising from directly measurable, minimally correlated quantities, yet the manuscript does not quantify how residual covariance between the X-ray luminosity and the mass proxy (e.g., via hydrostatic equilibrium assumptions or gas clumping in the 10^13-10^14 M_⊙ regime) propagates into the significance. A bias in M that correlates with Lx at even the 10-20% level could shift the comparison horizontally and reduce the tension below 5σ; an explicit error-budget table or Monte-Carlo test of this effect is needed to support the quoted significance.
Authors: We agree that an explicit quantification of residual covariance between the mass proxy and L_X would further support the quoted significance. While L_X and T are directly measured with low correlation, the mass estimate (derived from the X-ray data) can in principle introduce covariance through hydrostatic assumptions or clumping. In the revised manuscript we will add a dedicated Monte-Carlo error-budget section that injects 10–20 % correlated biases in the mass proxy and recomputes the offset significance. Our preliminary tests show the tension remains above 5σ, but the full table will be included so readers can evaluate the effect directly. revision: yes
-
Referee: [§3] §3 (sample and mass estimation): The claim of robustness to selection effects is stated, but the text does not demonstrate that the 44-group XMM-Newton sample is representative of the underlying population at fixed true mass rather than at fixed observed Lx or T. If the selection function correlates with the same gas properties that set Lx, the under-prediction by ejective models could be partly selection-driven; a completeness simulation or comparison to a volume-limited catalog would strengthen the result.
Authors: We acknowledge that the manuscript states robustness to selection without a dedicated demonstration. The 44 groups comprise all systems in the local volume with high-quality XMM-Newton data meeting our quality cuts; the selection is therefore driven by data quality rather than by L_X or T thresholds. To address the referee’s concern, the revised version will include a direct comparison of the sample’s mass and temperature distributions against a larger, volume-limited group catalog from the literature, together with a brief discussion of how any residual selection bias would affect the L_X–M comparison. revision: yes
Circularity Check
No circularity: independent X-ray observables tested against external simulation outputs
full rationale
The paper derives its central claim by comparing directly measured X-ray luminosity and temperature from a sample of 44 galaxy groups observed with XMM-Newton against predictions from cosmological simulations (e.g., FLAMINGO) that adopt varying feedback efficiencies. The reported 5.7σ under-prediction of luminosity at fixed mass uses quantities the paper explicitly describes as minimally correlated and robust to selection. No equation or step reduces the discrepancy to a fitted parameter, self-defined quantity, or load-bearing self-citation; simulation results function as external inputs calibrated on baryon fractions, while the observational scaling relations serve as an independent test. The derivation chain therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The X-ray luminosity-temperature relation in the 44-group sample is minimally affected by selection effects and can be directly compared to simulation predictions.
Forward citations
Cited by 2 Pith papers
-
The limits of feedback from active galactic nuclei
AGN feedback creates a mass-independent entropy ceiling that allows outflows to escape halos only below M_200m = 10^13.7 M_sun, explaining depleted gas in groups versus near-cosmic fractions in clusters.
-
Signatures of Suppressed Matter Clustering revealed by Fast Radio Bursts
FRB dispersion measures directly constrain suppression of the matter power spectrum due to feedback at k ~ 0.1-3 h/Mpc, reduce posterior variance by a factor of ~8 at k~1 h/Mpc, and exclude extreme large-scale feedbac...
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
- [3]
-
[4]
E., Bulbul , E., Ghirardini , V., et al
Bahar , Y. E., Bulbul , E., Ghirardini , V., et al. 2024, arXiv e-prints, arXiv:2401.17276
- [5]
-
[6]
A., Ir s i c , V., Amon , A., & Sijacki , D
Bigwood , L., Bourne , M. A., Ir s i c , V., Amon , A., & Sijacki , D. 2025, [ [arXiv] 2501.16983 ]
-
[7]
Booth , C. M. & Schaye , J. 2009, , 398, 53
work page 2009
-
[8]
Braspenning , J., Schaye , J., Schaller , M., et al. 2024, , 533, 2656
work page 2024
-
[9]
Chisari , N. E., Mead , A. J., Joudaki , S., et al. 2019, The Open Journal of Astrophysics, 2, 4
work page 2019
- [10]
-
[11]
Damsted , S., Finoguenov , A., Lietzen , H., et al. 2024, , 690, A52
work page 2024
- [12]
-
[13]
Eckert , D., Ettori , S., Pointecouteau , E., van der Burg , R. F. J., & Loubser , S. I. 2022, , 662, A123
work page 2022
-
[14]
2020, The Open Journal of Astrophysics, 3, 12
Eckert , D., Finoguenov , A., Ghirardini , V., et al. 2020, The Open Journal of Astrophysics, 3, 12
work page 2020
-
[15]
Eckert , D., Gaspari , M., Gastaldello , F., Le Brun , A. M. C., & O'Sullivan , E. 2021, Universe, 7, 142
work page 2021
-
[16]
Eckert , D., Gastaldello , F., Lovisari , L., et al. 2025, , 701, A127
work page 2025
-
[17]
Eckert , D., Gastaldello , F., O'Sullivan , E., Finoguenov , A., & Brienza , M. 2024, Galaxies, 12, 24
work page 2024
-
[18]
Gianfagna , G., De Petris , M., Yepes , G., et al. 2021, , 502, 5115
work page 2021
- [19]
-
[20]
Hadzhiyska , B., Ferraro , S., Ried Guachalla , B., et al. 2025, , 112, 083509
work page 2025
-
[21]
A., Puchwein , E., Shen , S., & Sijacki , D
Henden , N. A., Puchwein , E., Shen , S., & Sijacki , D. 2018, , 479, 5385
work page 2018
-
[22]
Lovisari , L., Ettori , S., Gaspari , M., & Giles , P. A. 2021, Universe, 7, 139
work page 2021
-
[23]
Lovisari , L., Reiprich , T. H., & Schellenberger , G. 2015, , 573, A118
work page 2015
-
[24]
G., Amon , A., Schaye , J., et al
McCarthy , I. G., Amon , A., Schaye , J., et al. 2025, , 540, 143
work page 2025
-
[25]
G., Schaye , J., Bird , S., & Le Brun , A
McCarthy , I. G., Schaye , J., Bird , S., & Le Brun , A. M. C. 2017, , 465, 2936
work page 2017
-
[26]
McCarthy , I. G., Schaye , J., Ponman , T. J., et al. 2010, , 406, 822
work page 2010
-
[27]
Munari , E., Biviano , A., Borgani , S., Murante , G., & Fabjan , D. 2013, , 430, 2638
work page 2013
-
[28]
Planck Collaboration , Ade , P. A. R., Aghanim , N., et al. 2016, , 594, A13
work page 2016
-
[29]
2024 a , arXiv e-prints, arXiv:2411.16555
Popesso , P., Biviano , A., Marini , I., et al. 2024 a , arXiv e-prints, arXiv:2411.16555
-
[30]
2024 b , arXiv e-prints, arXiv:2411.17120
Popesso , P., Marini , I., Dolag , K., et al. 2024 b , arXiv e-prints, arXiv:2411.17120
-
[31]
2025, arXiv e-prints, arXiv:2503.19870
Ried Guachalla , B., Schaan , E., Hadzhiyska , B., et al. 2025, arXiv e-prints, arXiv:2503.19870
-
[32]
Robotham , A. S. G., Norberg , P., Driver , S. P., et al. 2011, , 416, 2640
work page 2011
-
[33]
Rossetti , M., Eckert , D., Gastaldello , F., et al. 2024, , 686, A68
work page 2024
-
[34]
Schaan , E., Ferraro , S., Amodeo , S., et al. 2021, , 103, 063513
work page 2021
- [35]
- [36]
- [37]
- [38]
-
[39]
Spinelli , C., Veronica , A., Pacaud , F., et al. 2025, , 700, A220
work page 2025
-
[40]
2012, New Journal of Physics, 14, 045004
Sun , M. 2012, New Journal of Physics, 14, 045004
work page 2012
-
[41]
Sun , M., Voit , G. M., Donahue , M., et al. 2009, , 693, 1142
work page 2009
-
[42]
Tempel , E., Tuvikene , T., Kipper , R., & Libeskind , N. I. 2017, , 602, A100
work page 2017
- [43]
- [44]
- [45]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.