REVIEW 2 major objections 5 minor 106 references
The youth of the intracluster medium. I. A non-parametric characterisation of the gas and electron number density profiles of $z \simeq 2$ protoclusters
T0 review · 2 major / 5 minor · reviewed 2026-07-12 · grok-4.5
Pith's one-line read At z~2, hot ionised gas in protocluster cores only dominates outside ~0.1–0.5 R500c and shows a strong double-β density shape that grows with mass and AGN activity.
desk verdict Solid large-N non-parametric maps of z≈2 proto-ICM density and ionisation from Magneticum; the double-β shape, ~70 ckpc dip, and secondary trends are real within the model, with the usual thermal-feedback and CIE caveats already flagged by the authors. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Non-parametric stacked radial profiles of gas density, temperature-binned mass fractions and electron number density, constructed for 3818 simulated protocluster cores and further stratified by M500c, the stellar mass ratio M12 of the two brightest galaxies, and the time-averaged Eddington ratio of the central black hole.
What would settle it
A statistically significant sample of z≃2 protoclusters with measured electron-density profiles (via deep X-ray or joint SZ+X-ray analysis) that lack the central double-β steepening or that show hot-ionised gas already dominant well inside 0.1 R500c would contradict the stacked simulation result.
Extended reading notes
Core claim
Protoclusters at z≃2 exhibit moderate departures from self-similar gas density and temperature structure, featuring a strong double-β profile that is most pronounced at high mass and high central AGN accretion; hot, ionised gas only dominates outside roughly 0.1–0.5 R500c, and its density at those radii correlates with both halo mass and dynamical disturbance.
Load-bearing premise
Electron densities rest on the assumption of collisional ionisation equilibrium with photo-ionisation neglected, so any substantial non-equilibrium or local AGN radiation field would reshape the reported ionisation fractions and the double-β electron-density profile.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a non-parametric characterisation of gas mass density, temperature structure, ionisation degree and electron number density profiles for 3818 protocluster cores at z ≃ 2 (M500c > 10^13 M⊙) drawn from Magneticum Box2b/hr. After re-centring with ASOHF, multiphase correction (Eq. 1), and CIE-based electron densities, the authors stack profiles versus mass and secondary indicators of assembly state (M12) and central AGN activity (integrated Eddington ratio f̄Edd). They report moderate departures from self-similarity, a ubiquitous double-β morphology that strengthens with mass and AGN accretion, and the dominance of hot, ionised gas only at intermediate radii (r ≳ [0.1–0.5] R500c). The work also supplies X-ray exposure-time estimates, a comparison to Spiderweb measurements, and an ANOVA-style variance decomposition showing that M500c, M12 and f̄Edd together explain up to ~30 % of the logarithmic scatter in ne(r).
Significance. If the reported trends hold, the paper supplies the first statistically large, simulation-based map of proto-ICM density and ionisation structure at z ≃ 2, directly usable for interpreting Faraday rotation, SZ and forthcoming X-ray constraints on high-redshift systems. Strengths include the carefully documented sample construction (ASOHF re-centring, multiphase treatment, CIE tables), robust stacking with bootstrap uncertainties, explicit Spearman correlations, and an open discussion of model limitations (thermal AGN feedback, Spiderweb tension). The non-parametric results are intended as the foundation for a forthcoming parametric calibration, which would be of clear practical value to observers.
major comments (2)
- §2.3 and §4.1: Electron densities rest on CIE and the neglect of local AGN photoionisation. The CIE timescale argument (4–6 orders of magnitude shorter than dynamical/cooling times) is persuasive for the bulk mass, and self-shielding is invoked inside ~0.1–0.2 R500c. However, the ~70 ckpc ionisation dip—central to the reported double-β ne shape and its correlation with f̄Edd—lies precisely where AGN radiation could raise the ionisation fraction. A quantitative upper bound (even a simple Strömgren-sphere or optically-thin estimate using the same SMBH accretion rates already measured) is needed to show that the dip and the secondary trends survive this uncertainty; without it the ionisation-structure claim remains only partially stress-tested.
- §4.3 and Fig. 7: The systematic under-prediction of Spiderweb central densities is acknowledged and partially mitigated by aperture and homogeneity corrections, yet none of the 3818 systems reach the observed mean ne inside 0.44 R500c. Given that the paper’s stated motivation is to provide density templates for interpreting high-z observables, the authors should either (i) quantify how much of the discrepancy is attributable to the purely thermal AGN feedback implementation (e.g., by citing or performing a controlled comparison with a kinetic-feedback run) or (ii) clearly demarcate the mass/radius regime in which the Magneticum profiles can be used as templates versus the regime in which they are known to be biased low.
minor comments (5)
- Fig. 1 (right) and Fig. 3 (bottom-right): Spearman coefficients are shown without indicating whether they are mass-corrected; the later panels of Fig. 4 do make this distinction. A uniform statement would improve clarity.
- §2.5: The integrated Eddington ratio averages over ~500 Myr between snapshots. A short sensitivity test (or at least a statement) using the instantaneous rate would reassure readers that the reported central-density correlations are not window-dependent.
- Appendix A / Fig. A.1: The example of Subfind–ASOHF centre mismatch is helpful; stating the median and 95th-percentile offsets in the main text (rather than only in the appendix) would strengthen the methods section.
- Fig. 6: The exposure-time forecast assumes Aeff = 500 cm2 and a fixed background; a brief note on how the contours scale with NewAthena’s larger effective area (even if angular resolution remains a limiting factor) would make the figure more forward-looking.
- Typographical: “ASOHF” is introduced with a footnote URL; a standard citation to Vallés-Pérez et al. (2022) in the main text would be cleaner. Occasional missing spaces after commas and inconsistent use of “ckpc” versus “kpc” appear in §§3.1–3.2.
Circularity Check
No circularity: direct non-parametric measurement of simulation profiles with independent secondary indicators; no fitted inputs renamed as predictions and no load-bearing self-citation chain.
full rationale
The paper extracts 3818 protocluster regions from the pre-existing Magneticum Box2b/hr run, recentres them with ASOHF, computes multiphase gas densities via the explicit rescaling of Eq. (1), obtains ne under stated CIE tables (Sect. 2.3), and stacks the resulting radial profiles versus M500c, M12 and fEdd (Sects. 2.4–2.5, 3). All reported trends (double-β shape, mass-dependent ionisation fractions, Spearman correlations, ANOVA variance fractions) are therefore direct statistics of the simulation particles; no free parameters are fitted to the target ne(r) or ρgas(r) profiles and then re-presented as predictions. M12 is defined from ASOHF stellar masses of BCG and first satellite; fEdd is the time-averaged BH accretion ratio between two snapshots—both independent of the density profiles they are later correlated against. Self-citations (Magneticum validation papers, prior ASOHF/Vallés-Pérez works) supply methodological background or low-z context and are not used to force the high-z characterisation itself. The Spiderweb comparison and X-ray exposure estimates are external benchmarks, not circular closures. The derivation chain is therefore self-contained against the simulation data and contains none of the six circular patterns.
Assumptions & free parameters
free parameters (2)
- Integrated Eddington-ratio averaging window (~500 Myr between snapshots) =
~500 Myr (snapshot spacing)
- Minimum particles per radial bin (Nbin_min_part=1000) and Δlog10 r=0.01 dex =
1000 particles, 0.01 dex
assumptions (4)
- domain assumption Collisional ionisation equilibrium (CIE) holds for radially-averaged ne profiles; photoionisation (UVB and local AGN) can be neglected inside the self-shielded core and is sub-dominant outside.
- domain assumption Cold sub-resolution ISM has negligible temperature and volume fraction, so diffuse gas properties are recovered by the simple rescaling mh=(1-fc)mt, Th=Tmw/(1-fc), ρh=(1-fc)ρt (Eq. 1).
- domain assumption Protocluster cores are adequately defined as spherical-overdensity peaks with M500c>10^13 M⊙ at z=1.98, without further environmental selection.
- domain assumption Magneticum Box2b/hr (WMAP7 cosmology, thermal AGN feedback, multiphase star formation) is a sufficiently faithful model of high-z baryonic physics for density and ionisation structure.
Cite this review
Pith. "Pith review of The youth of the intracluster medium. I. A non-parametric characterisation of the gas and electron number density profiles of $z \simeq 2$ protoclusters." pith.science (2026). https://pith.science/paper/GJBF7TWN
@misc{pith2026260702654,
author = {Pith},
title = {Pith review of: The youth of the intracluster medium. I. A non-parametric characterisation of the gas and electron number density profiles of $z \simeq 2$ protoclusters},
year = {2026},
howpublished = {\url{https://pith.science/paper/GJBF7TWN}},
note = {Machine review of arXiv:2607.02654}
}
abstract
Context. Protoclusters of galaxies are the earliest phase in the assembly of galaxy clusters and can provide invaluable information about plasma physics, cosmic magnetism, and cosmology. However, due to small angular sizes and cosmological dimming, observing the proto-intracluster medium (proto-ICM) associated with protocluster cores is far from trivial. Aims. We aim to provide a non-parametric description of the gas mass and electron number density profiles of the proto-ICM at $z = 2$, and to study their dependence on mass, dynamical state and central activity. Methods. We extract and analyse over $3800$ regions around protocluster cores with spherical-overdensity masses above $M_\mathrm{500c} > 10^{13} \, M_\odot$ out of a large simulated volume within the Magneticum suite. We study their density profiles, temperature structure, ionisation degree and electron number density as a function of mass and other secondary properties characterising dynamical state and central activity, extending from the central halo to the surrounding protocluster environment. Results. Protoclusters present moderate deviations from self-similarity in their density profiles and temperature structure, with a strong double-$\beta$ structure especially relevant at high masses and intense AGN accretion. Hot, ionised gas is only dominant at intermediate radii ($r \gtrsim [0.1-0.5] R_\mathrm{500c}$), where its density also correlates with mass and dynamical disturbance. Conclusions. These results constitute the basis for a forthcoming parametric calibration of proto-ICM density profiles, which could be useful for interpreting observables sensitive to the density and ionisation of the diffuse gas.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
The Spiderweb Protocluster is Being Magnetized by Its Central Radio Jet. , keywords =. doi:10.3847/1538-4357/ac7ec0 , archivePrefix =. 2207.03498 , primaryClass =
-
[2]
Scaling relations for galaxy clusters in the Millennium-XXL simulation. , keywords =. doi:10.1111/j.1365-2966.2012.21830.x , archivePrefix =. 1203.3216 , primaryClass =
-
[3]
Shock and splash: gas and dark matter halo boundaries around CDM galaxy clusters. , keywords =. doi:10.1093/mnras/stab2598 , archivePrefix =. 2012.00977 , primaryClass =
-
[4]
The SRG/eROSITA All-Sky Survey: Constraints on AGN feedback in galaxy groups. , keywords =. doi:10.1051/0004-6361/202449399 , archivePrefix =. 2401.17276 , primaryClass =
-
[5]
Kinetic AGN feedback effects on cluster cool cores simulated using SPH. , keywords =. doi:10.1093/mnras/stw1389 , archivePrefix =. 1605.08051 , primaryClass =
-
[6]
The DIANOGA simulations of galaxy clusters: characterizing star formation in proto-clusters
The DIANOGA simulations of galaxy clusters: characterising star formation in protoclusters. , keywords =. doi:10.1051/0004-6361/202038396 , archivePrefix =. 2006.13951 , primaryClass =
work page Pith review arXiv doi:10.1051/0004-6361/202038396 2006
-
[7]
An improved SPH scheme for cosmological simulations
An improved SPH scheme for cosmological simulations. , keywords =. doi:10.1093/mnras/stv2443 , archivePrefix =. 1502.07358 , primaryClass =
-
[8]
Measures of Location and Scale for Velocities in Clusters of Galaxies---A Robust Approach. , keywords =. doi:10.1086/115487 , adsurl =
Show all 106 references
- [9]
- [10]
-
[11]
An X-ray-emitting protocluster at z ≈5.7 reveals rapid structure growth , url =
Bogd. An X-ray-emitting protocluster at z ≈5.7 reveals rapid structure growth , url =. Nature , number =. 2026 , bdsk-url-1 =. doi:10.1038/s41586-025-09973-1 , id =
2026 doi
-
[12]
Exploring the properties of proto-clusters in the Simulated Euclid Wide Survey
Euclid preparation: LV. Exploring the properties of proto-clusters in the Simulated Euclid Wide Survey. , keywords =. doi:10.1051/0004-6361/202451683 , archivePrefix =. 2407.19919 , primaryClass =
- [13]
- [14]
-
[15]
, keywords =
X-rays from hot plasma in clusters of galaxies. , keywords =
- [16]
-
[17]
, keywords =
An Overabundance of Radio-AGN in the SPT2349-56 Protocluster: Preheating the Intra-Cluster Medium. , keywords =. doi:10.48550/arXiv.2511.17814 , archivePrefix =. 2511.17814 , primaryClass =
- [18]
- [19]
- [20]
- [21]
- [22]
- [23]
-
[24]
, keywords =
Signatures of dynamical activity in the hot gas profiles of groups and clusters in the FLAMINGO simulations. , keywords =. doi:10.1093/mnras/staf2073 , archivePrefix =. 2509.10230 , primaryClass =
-
[25]
, keywords =
Substructure in Clusters as a Cosmological Test. , keywords =. doi:10.1086/177626 , adsurl =
- [26]
-
[27]
, keywords =
Improving convergence in smoothed particle hydrodynamics simulations without pairing instability. , keywords =. doi:10.1111/j.1365-2966.2012.21439.x , archivePrefix =. 1204.2471 , primaryClass =
2012 doi
- [28]
- [29]
- [30]
- [31]
-
[32]
, keywords =
Substructures in hydrodynamical cluster simulations. , keywords =. doi:10.1111/j.1365-2966.2009.15034.x , archivePrefix =. 0808.3401 , primaryClass =
2009 doi
- [33]
- [34]
-
[35]
, keywords =
Encyclopedia Magneticum: Scaling Relations from Cosmic Dawn to Present Day. , keywords =. doi:10.48550/arXiv.2504.01061 , archivePrefix =. 2504.01061 , primaryClass =
- [36]
- [37]
- [38]
-
[39]
, keywords =
Simulating the effect of active galactic nuclei feedback on the metal enrichment of galaxy clusters. , keywords =. doi:10.1111/j.1365-2966.2009.15794.x , archivePrefix =. 0909.0664 , primaryClass =
2009 doi
-
[40]
Statistical methods for research workers
-
[41]
Atoms , year = 2020, month = aug, volume =
PyAtomDB: Extending the AtomDB Atomic Database to Model New Plasma Processes and Uncertainties. Atoms , year = 2020, month = aug, volume =. doi:10.3390/atoms8030049 , adsurl =
2020 doi
- [42]
- [43]
- [44]
- [45]
- [46]
- [47]
- [48]
- [49]
-
[50]
, keywords =
Chandra observation of the edge-on spiral NGC 5775: probing the hot galactic disc/halo connection. , keywords =. doi:10.1111/j.1365-2966.2008.13749.x , archivePrefix =. 0807.3587 , primaryClass =
2008 doi
-
[51]
, keywords =
The nature and space density of fossil groups of galaxies. , keywords =. doi:10.1046/j.1365-8711.2003.06702.x , archivePrefix =. astro-ph/0304257 , primaryClass =
2003 doi
- [52]
-
[53]
The structure of star clusters. III. Some simple dynamical models. , year = 1966, month = feb, volume =. doi:10.1086/109857 , adsurl =
1966 doi
- [54]
- [55]
- [56]
-
[57]
, keywords =
Merger rates in hierarchical models of galaxy formation. , keywords =. doi:10.1093/mnras/262.3.627 , adsurl =
- [58]
- [59]
-
[60]
, keywords =
Deep XMM-Newton observations of the most distant SPT-SZ galaxy cluster. , keywords =. doi:10.1093/mnras/staa1581 , archivePrefix =. 2006.02009 , primaryClass =
2006 doi
- [61]
- [62]
- [63]
- [64]
- [65]
- [66]
-
[67]
, keywords =
The statistics of CDM halo concentrations. , keywords =. doi:10.1111/j.1365-2966.2007.12381.x , archivePrefix =. 0706.2919 , primaryClass =
2007 doi
- [68]
- [69]
- [70]
- [71]
- [72]
-
[73]
Space Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray , year = 2018, editor =
The complicated evolution of the ACIS contamination layer over the mission life of the Chandra X-ray Observatory. Space Telescopes and Instrumentation 2018: Ultraviolet to Gamma Ray , year = 2018, editor =. doi:10.1117/12.2312748 , archivePrefix =. 1809.02225 , primaryClass =
- [74]
- [75]
- [76]
- [77]
-
[78]
Radiative processes in astrophysics
-
[79]
X-ray emission from clusters of galaxies
-
[80]
, keywords =
Simulating the formation of a protocluster at z -0.5ex 2. , keywords =. doi:10.1111/j.1365-2966.2008.14192.x , archivePrefix =. 0810.2791 , primaryClass =
2008 doi
-
[81]
, keywords =
A unified model for AGN feedback in cosmological simulations of structure formation. , keywords =. doi:10.1111/j.1365-2966.2007.12153.x , archivePrefix =. 0705.2238 , primaryClass =
2007 doi
- [82]
-
[83]
Results at z=0
Populating a cluster of galaxies - I. Results at z=0. , keywords =. doi:10.1046/j.1365-8711.2001.04912.x , archivePrefix =. astro-ph/0012055 , primaryClass =
2001 doi
-
[84]
, keywords =
Cosmological smoothed particle hydrodynamics simulations: a hybrid multiphase model for star formation. , keywords =. doi:10.1046/j.1365-8711.2003.06206.x , archivePrefix =. astro-ph/0206393 , primaryClass =
2003 doi
-
[85]
, keywords =
The cosmological simulation code GADGET-2. , keywords =. doi:10.1111/j.1365-2966.2005.09655.x , archivePrefix =. astro-ph/0505010 , primaryClass =
2005 doi
-
[86]
, keywords =
Modelling feedback from stars and black holes in galaxy mergers. , keywords =. doi:10.1111/j.1365-2966.2005.09238.x , archivePrefix =. astro-ph/0411108 , primaryClass =
2005 doi
- [87]
- [88]
-
[89]
, year = 1972, month = aug, volume =
Formation of Clusters of Galaxies; Protocluster Fragmentation and Intergalactic Gas Heating. , year = 1972, month = aug, volume =
1972
-
[90]
Comments on Astrophysics and Space Physics , keywords =
The Observations of Relic Radiation as a Test of the Nature of X-Ray Radiation from the Clusters of Galaxies. Comments on Astrophysics and Space Physics , keywords =
-
[91]
, keywords =
Chemical enrichment of galaxy clusters from hydrodynamical simulations. , keywords =. doi:10.1111/j.1365-2966.2007.12070.x , archivePrefix =. 0705.1921 , primaryClass =
2007 doi
- [92]
- [93]
- [94]
- [95]
- [96]
-
[97]
Impact of dark matter halo and intracluster medium properties on their full assembly histories
The eventful life journey of galaxy clusters: I. Impact of dark matter halo and intracluster medium properties on their full assembly histories. , keywords =. doi:10.1051/0004-6361/202453483 , archivePrefix =. 2505.02891 , primaryClass =
-
[98]
Impact of mass accretion on the thermodynamic structure of the ICM
The eventful life journey of galaxy clusters: II. Impact of mass accretion on the thermodynamic structure of the ICM. , keywords =. doi:10.1051/0004-6361/202557496 , adsurl =
- [99]
- [100]
-
[101]
Atomic Data and Nuclear Data Tables , year = 1997, month = jan, volume =
A Practical Fit Formula for Ionization Rate Coefficients of Atoms and Ions by Electron Impact: Z = 1-28. Atomic Data and Nuclear Data Tables , year = 1997, month = jan, volume =. doi:10.1006/adnd.1997.0732 , adsurl =
1997 doi
-
[102]
, keywords =
Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations. , keywords =. doi:10.1111/j.1365-2966.2009.15331.x , archivePrefix =. 0902.1535 , primaryClass =
2009 doi
- [103]
-
[104]
and Aravena, Manuel and Araya-Araya, Pablo and Archipley, Melanie and Cathey, Jared and Deane, Roger P
Zhou, Dazhi and Chapman, Scott C. and Aravena, Manuel and Araya-Araya, Pablo and Archipley, Melanie and Cathey, Jared and Deane, Roger P. and Di Mascolo, Luca and Gobat, Raphael and Greve, Thomas R. and Hill, Ryley and Kim, Seonwoo and Phadke, Kedar A. and Pillai, Vismaya R. a...
2026 doi
- [105]
- [106]
Reviewed July 12, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.