{"id":"06c37edc-fe24-4555-bd12-a6b723f38c17","arxiv_id":"2504.20159","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Simulated and observed cluster baryon fraction scaling slopes agree within 1 to 2 sigma, and the baryons missing from the observational count are mostly warm gas and intracluster light.","lead":"This paper compares how the fraction of ordinary matter in galaxy clusters changes with cluster mass in the IllustrisTNG simulations versus SPT-SZ survey observations, and finds the slopes agree within 1 to 2 sigma. It then tallies which baryon components are absent from the observational measurement (mainly warm gas and intracluster light) and links them to the missing baryon problem.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1–2σ slope agreement with C18 is likely a mass-range artifact: the paper's own BPL fits yield a high-mass slope of ~0.05–0.06 across the C18 mass range (log10 M500c ≈ 14.4–15.1), versus C18's 0.35, so the SPL agreement rests on fitting lower masses C18 never probed.","rationale":"The paper has genuine strengths: it uses public IllustrisTNG data, its component decomposition (Table 3, Appendix A) is transparent and reproducible, the BPL/GDPL formalism is borrowed from P22, and the conclusion that warm gas plus ICL, with cold gas negligible, constitute the baryons missed by C18-type measurements is concrete and partly independent of the slope comparison. The main comparison, however, does not establish what the abstract claims. The reader identified the proxy equivalence (hot gas + rop stars versus Chandra β-model plus SED stellar masses) as the weakest assumption; that is a fair concern, but it is not the most load-bearing. The most decisive issue is that the two slopes being compared are fitted over different mass ranges, and the paper's own Table 5 already provides the numbers showing the mismatch: above the ~10^14 pivot, which sits inside the C18 range, the high-mass slope is 0.05–0.06, in strong tension with C18's 0.35. The full-range SPL slopes that appear to agree are steepened by the 63–69% of haloes below 10^14 that C18 never observes. Because this affects the primary claim rather than the written record, it is more serious than the internal inconsistencies the reader catalogued (1605 versus 1604 sample sizes, swapped 17.67/17.80 percentages, Fig. 2 slope labels disagreeing with the text), which matter for reliability but not for the physics. A single re-fit restricted to the C18 mass range settles the question. If the restricted slopes confirm ~0.05–0.15, the headline agreement should be withdrawn, while the component analysis, breakpoint location, and FoF discussion would remain as contributions. I therefore keep the verdict CONDITIONAL, with the matched-mass-range refit as the explicit condition; confirmation of the discrepancy would move the verdict to REJECT.","tokens_in":28055,"tokens_out":17900,"duration_ms":170351,"concrete_test":"Re-fit the SPL to the TNG100 and TNG300 haloes restricted to the C18 mass range, log10(M500c/Msun) ∈ [14.42, 15.10], using the same non-linear least-squares procedure as Section 4.1, and compare the restricted slopes with αChiu = 0.350 ± 0.076. If the restricted slopes lie within 2σ of 0.35, the concern is refuted; if they land near the Table 5 high-mass slopes (~0.05–0.15), the reported 1–2σ agreement is a fitting-range artifact and the central claim fails. As a cross-check, recompute the full-range SPL slope with each simulation halo weighted by the C18 sample's mass distribution; a drop from ~0.26 toward ~0.15 or below would confirm the mass-range dependence.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline claim — that IllustrisTNG baryon-fraction scaling slopes agree with C18 within 1–2σ — appears to be an artifact of fitting over different mass ranges. C18's slope, αChiu = 0.350 ± 0.076, is measured over their 91 clusters spanning log10(M500c/Msun) ∈ [14.42, 15.10] (Section 2.2). The simulation SPL slopes (αSPL = 0.248 ± 0.018 for TNG100, 0.261 ± 0.006 for TNG300, Section 4.1) are fitted over M200c ≥ 7×10^13 Msun, reaching down to log10 M500c ≈ 13.5, roughly 0.9 dex below the C18 range. The samples are dominated by low-mass haloes (137/218 and 1115/1605 below log10 M500c = 14), where the local slope is steep (Table 5: αBPL,1 = 0.79 for TNG100, 0.35 for TNG300). Above the pivots (log10 Mpivot = 13.80 and 14.07), which lie below the C18 range, the paper's own BPL/GDPL fits give αBPL,2 = 0.06 ± 0.03 and 0.05 ± 0.01; the GDPL transition (δ ≈ 0.25) is sharp, so the simulation slope across 14.4–15.1 is ~0.05–0.10, a ~3–4σ discrepancy with C18's 0.35. The apparent agreement is produced by extrapolating the SPL into a mass regime C18 does not observe and weighting it with many low-mass haloes. The instability of the simulation slope itself confirms the problem: 0.27/0.16 when fitted to mass-binned medians versus 0.248/0.261 when fitted to individual haloes. Section 5.1's unpropagated 14–20% mass-bias systematic is a related but secondary issue.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper compares baryonic-mass-fraction scaling relations of galaxy clusters in the IllustrisTNG simulations (TNG100 and TNG300) with the SPT-SZ/Chandra sample of Chiu et al. (2018, C18). It fits simple power law (SPL), broken power law (BPL), and general double power law (GDPL) models to fbar,500c versus M500c and versus redshift, using two measurement proxies: hot gas (T≥10^7 K) plus galaxy stellar mass within a 2× half-mass-radius aperture, and all gas plus all stars within R500c. The paper reports approximately null redshift slopes, a 1–2σ agreement of SPL slopes with C18, a breakpoint at log10(M500c/Msun)≈14 where self-similarity is lost, and a 'missing' baryon fraction of about 13.8–14.1% dominated by warm gas and intracluster light. It also examines FoF-based scaling relations and the role of mass accretion in the baryon deficit.","tokens_in":28488,"tokens_out":10783,"duration_ms":100949,"significance":"If the claimed slope agreement were robust, it would strengthen confidence in IllustrisTNG's feedback model as a tool for interpreting cluster baryon fractions and the missing-baryon problem. The paper provides a useful component-by-component decomposition of baryons in simulated clusters and a clear measurement of the BPL/GDPL breakpoint, which are of interest to the cluster scaling-relation community. However, the central comparison to C18 is compromised by a mass-range mismatch: the simulation SPL is fitted over a much wider range that includes many low-mass haloes, while C18's slope is measured in a narrow high-mass range. The paper itself reports high-mass BPL slopes of 0.05–0.06, which differ from C18's 0.35 by several sigma. In addition, the paper identifies but does not propagate a 14–20% mass-bias systematic. These issues must be resolved before the main claim can be accepted.","major_comments":[{"comment":"","section":"Sec. 4.1, Fig. 3, and Table 5"},{"comment":"","section":"Sec. 5.1"},{"comment":"","section":"Abstract and Sec. 4.2, Table 3"},{"comment":"","section":"Sec. 4.1, Fig. 2, and Abstract"}],"minor_comments":[{"comment":"","section":"Abstract, Sec. 2.1, and Sec. 6"},{"comment":"","section":"Sec. 6, summary bullet"},{"comment":"","section":"Sec. 4.1 and Fig. 2"},{"comment":"","section":"Sec. 4.4 and Table 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript contains useful simulation-side measurements, particularly the component decomposition and the BPL/GDPL breakpoint, but the central comparison to C18 is not currently supported because of the mass-range mismatch and the unpropagated mass-bias systematic. I recommend major revision with a request for a mass-restricted analysis and a clear statement of which fitted slope is used for the headline claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read this one. First, the headline claim — that IllustrisTNG baryon fraction slopes agree with Chiu et al. (2018) within 1–2σ — is a mass-range artifact. Second, the baryon component budget underneath it is real and useful, and worth keeping after the headline gets fixed.\n\nThe SPL slopes in Section 4.1 are fitted over M200c ≥ 7×10^13 Msun, reaching down to log M500c ≈ 13.5, while C18's sample spans log M500c = 14.42–15.10. The paper's own Table 5 gives BPL high-mass slopes of 0.06 and 0.05, with pivots at 13.80 and 14.07 — both below the C18 range, and the GDPL transition is sharp (δ = 0.25). Across the mass range C18 actually observed, TNG therefore predicts a near-flat baryon fraction, ~0.06, against C18's 0.35 ± 0.076. That is a 3–4σ discrepancy. The apparent agreement comes from the many low-mass haloes dominating the TNG SPL fit. The reader's report was too generous on this point: the paper's own fits refute its abstract, so the issue is more than conditional.\n\nWhat is genuinely good: the decomposition into cold, warm, and hot gas, galaxy stars (both rop and Subfind apertures), and ICL is the most complete I have seen for this comparison. The missing-baryon budget — mostly warm gas plus ICL, larger in low-mass haloes — is a concrete statement people will want to quote, and the breakpoint at log M ~ 14 lines up with Pop et al. (2022a) and Ayromlou et al. (2023). The fits are straightforward, the engagement with the literature is honest, and the analysis is reproducible: TNG is public and the C18 numbers are in their tables.\n\nThe soft spots beyond the mass-range issue: Section 5.1 computes that SZ/X-ray mass biases would overestimate observed baryon fractions by 14–20 per cent, then drops it without propagating into the main comparison — a significant unquantified systematic, not a footnote. There are also sloppy internal inconsistencies: TNG300 is 1605 haloes in the abstract and Section 2.1, 1604 in the summary; the summary swaps the low-mass missing-baryon percentages relative to Table 3; a couple of redshift slopes differ between text and figure labels. All fixable, but they undermine confidence in the written record. And the abstract's WHIM framing overstates what was measured: the warm gas sits inside R500c, which the text itself later correctly calls WCGM.\n\nWho is this for? Anyone working on cluster scaling relations or the missing baryon problem. The component budget is worth keeping; the slope comparison needs to be redone honestly, fit over the C18 mass range and with the bias correction folded in. I would send it to a serious referee — the material is real and the fixes are achievable — but the mass-range issue should come back as major.","headline":"The headline 1–2σ slope agreement is a mass-range artifact — over the mass range C18 actually observed, the paper's own BPL fits give a slope near 0.06 versus C18's 0.35 — but the baryonic component budget is a genuinely useful new result.","tokens_in":29158,"tokens_out":9334,"would_cite":false,"duration_ms":80707,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The baryon fraction scaling slopes of simulated clusters match the observed SPT-SZ/Chandra slopes within 1-2 sigma, and the missing baryons are predominantly warm gas and intracluster light.","keywords":["galaxy clusters","baryon fraction","scaling relations","missing baryons","intracluster light","warm-hot intergalactic medium","IllustrisTNG","self-similarity"],"falsifier":"Measure the warm gas ($10^5\\lesssim T\\lesssim10^7$ K) mass in low-mass clusters ($\\log_{10}(M_{500c}/M_\\odot)\\lesssim14$) using O VII absorption stacking or thermal Sunyaev-Zel'dovich observations and compare with the roughly 10 per cent baryonic fraction predicted here; a detection below about 5 per cent would rule out warm gas as the dominant missing component and contradict the paper's missing-baryon budget.","tokens_in":27753,"feed_emoji":"🌌","tokens_out":8779,"duration_ms":76110,"temperature":0.7,"pith_summary":"This paper tries to establish that the IllustrisTNG cosmological simulations reproduce the observed scaling of cluster baryon fraction with mass, so the well-known missing baryon deficit can be attributed to specific undetected components. The authors compare 218 TNG100 and 1605 TNG300 haloes against 91 SPT-SZ/Chandra clusters from Chiu et al. (2018) using simple, broken, and general double power-law fits. They find baryon fraction slopes agree within 1-2 $\\sigma$, identify a self-similarity break at $\\log_{10}(M_{500c}/M_\\odot)\\sim 14$, and quantify that about 13.8-14.1 per cent of baryons are missing in the observational measurement. The missing budget is mostly warm gas at $10^5$-$10^7$ K in low-mass haloes and intracluster light at all masses, which would resolve the tension between cluster baryon fractions and the hierarchical formation picture.","feed_headline":"Simulations pin cluster missing baryons to warm gas and starlight","feed_subtitle":"Simulation matches observed baryon scaling; missing baryons are warm gas and diffuse starlight.","key_machinery":"The load-bearing object is the baryon mass fraction $f_{\\rm bar,500c}=M_{\\rm bar,500c}/M_{500c}$ measured in two ways: an observational proxy (hot gas with $T\\ge10^7$ K within $R_{500c}$ plus galaxy stellar mass within $2\\times R_{M_\\star/2}$, excluding warm gas, cold gas, and intracluster light) and a full baryon budget (all gas and all stars within $R_{500c}$). The argument runs through fits of the self-similar scaling relation $f_{\\rm bar}\\propto M^{\\alpha}$ using a simple power law, a broken power law, and a general double power law with a smoothness parameter $\\delta$; the breakpoint $M_{\\rm pivot}$ identifies the mass where self-similarity is lost. The mechanism that explains the break is AGN feedback pushing baryons into the cluster outskirts, and the friends-of-friends mass definition is used to test whether the universal baryon fraction is recovered when all gravitationally bound material is counted.","core_discovery":"The central discovery is that the observed baryon fraction scaling is reproduced by IllustrisTNG once the same measurement cuts are applied: hot gas with $T\\geq 10^7$ K inside $R_{500c}$ plus galaxy stellar mass within an aperture of $2\\times R_{M_\\star/2}$. With this proxy, TNG100 and TNG300 give SPL slopes $\\alpha_{\\rm SPL}=0.248\\pm0.018$ and $0.261\\pm0.006$, respectively, consistent at 1-2 $\\sigma$ with the slope $\\alpha_{\\rm Chiu}=0.350\\pm0.076$ reported by Chiu et al. (2018) for the SPT-SZ sample. The baryon fraction is consistent with zero redshift evolution (slopes $\\sim0.00$-$0.07$), and the missing baryons, defined as cold gas, warm gas, and intracluster light excluded by the observational proxy, amount to approximately 13.8-14.1 per cent of the baryon budget, with warm gas dominating in low-mass haloes (~10.5 per cent) and ICL contributing roughly 5 per cent at all masses. Fitting a broken or general double power law places the self-similarity break at $\\log_{10}(M_{\\rm pivot}/M_\\odot)\\sim13.80$-$14.07$, the mass below which AGN feedback redistributes baryons to the outskirts; when halo mass is instead defined by the friends-of-friends algorithm, the slope flattens toward zero, consistent with hierarchical assembly.","pith_inferences":["An extension would be to stack thermal Sunyaev-Zel'dovich or soft X-ray data on low-mass clusters to measure the warm gas component; if the warm gas mass is below roughly 5 per cent of baryons in that mass range, the identification of warm gas as the dominant missing reservoir would be falsified.","The ICL fraction of about 5 per cent could be tested with deep surface-photometry surveys that separate intracluster light from galaxy light; the paper's proxy choice alone shifts the inferred missing fraction by about 5 per cent, so the definition of a galaxy in future surveys may matter as much as the sensitivity.","If baryons are truly redistributed to 1.5-2.5 $R_{200c}$ in low-mass haloes, then cluster outskirts become the natural place to search for missing baryons, and cluster mass calibration via SZ or X-ray may require a radius-dependent correction.","The authors' result that about 21 per cent of baryons are missing even when all bound material is counted suggests the issue is not just where the baryons are, but how efficiently haloes accrete gas versus dark matter; comparing with particle-based gas simulations could settle whether the deficit is physical or numerical."],"forward_implications":["If the slope agreement holds, the IllustrisTNG feedback model can be used to predict the baryon budget of clusters and the exact components an observational survey is missing.","The missing-baryon deficit in observed low-mass clusters should be recoverable as warm gas ($10^5$-$10^7$ K) and diffuse intracluster light rather than cold gas.","The breakpoint at $M_{500c}\\sim10^{14}M_\\odot$ marks a physically distinct regime: above it clusters are self-similar, below it AGN feedback redistributes baryons beyond $R_{500c}$.","With friends-of-friends masses the baryon fraction slope approaches zero, so high-mass clusters can form from low-mass ones and the apparent contradiction with hierarchical growth is resolved.","The SZ and X-ray mass biases estimated from mock observations imply real observed baryon fractions could be lower by 14-20 per cent than currently estimated, strengthening the need for the missing warm component."],"supporting_citations":[{"why":"Supplies the observational sample of 91 SPT-SZ/Chandra clusters and the measured baryon fraction slopes that the simulations are compared against.","marker":"C18"},{"why":"Provides the general double power-law model and the prior demonstration that a single breakpoint captures AGN feedback effects in IllustrisTNG.","marker":"Pop et al. (2022a)"},{"why":"Supplies the closure radius idea and the 1.5-2.5 $R_{200c}$ redistribution range used to explain missing baryons outside $R_{500c}$.","marker":"Ayromlou et al. (2023)"},{"why":"Defines the self-similar scaling framework that predicts zero slope for the baryon fraction as a function of mass and redshift.","marker":"Kaiser (1986)"},{"why":"Provides the universal baryon fraction used as the reference value and the 3-sigma band in the scaling relation plots.","marker":"Planck Collaboration et al. (2016a)"},{"why":"Provides the SZ and X-ray mass bias values used to estimate how much observational mass errors would inflate the baryon fraction.","marker":"Pop et al. (2022b)"},{"why":"Supplies the gas phase decomposition in TNG300 used to identify warm circumcluster gas as the dominant missing phase.","marker":"Gouin et al. (2022)"},{"why":"Provides caustic-based mass accretion rate measurements in TNG300 used to discuss why dark matter may accrete more efficiently than gas.","marker":"Pizzardo et al. (2023b)"}],"fun_headline_variants":["Simulations pinpoint missing cluster baryons as warm gas and starlight","Simulations match cluster baryon scaling; missing baryons are warm gas and starlight","Cluster simulations link missing baryons to warm gas and starlight","Simulations resolve cluster missing baryons as warm gas and starlight"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulated proxy for the observed baryon fraction (hot gas above $10^7$ K within $R_{500c}$ plus galaxy stars in a $2R_{M_\\star/2}$ aperture) is taken to be directly comparable to the Chandra/SED measurement, so any systematic mismatch in how stars are assigned to galaxies or gas to the hot phase would shift the fitted slopes and the inferred missing budget.","fun_headline_variants_meta":{"raw":{"variants":["Simulations pinpoint missing cluster baryons as warm gas and starlight","Simulations match cluster baryon scaling; missing baryons are warm gas and starlight","Cluster simulations link missing baryons to warm gas and starlight","Simulations resolve cluster missing baryons as warm gas and starlight"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00102,"raw_usage":{"total_tokens":4436,"prompt_tokens":1210,"completion_tokens":3226,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":826,"completion_tokens_details":{"reasoning_tokens":3144}},"tokens_in":826,"tokens_out":3226,"duration_ms":19588,"temperature":1.0,"reasoning_tokens":3144,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:35:59.744639+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the warm gas ($10^5\\lesssim T\\lesssim10^7$ K) mass in low-mass clusters ($\\log_{10}(M_{500c}/M_\\odot)\\lesssim14$) using O VII absorption stacking or thermal Sunyaev-Zel'dovich observations and compare with the roughly 10 per cent baryonic fraction predicted here; a detection below about 5 per cent would rule out warm gas as the dominant missing component and contradict the paper's missing-baryon budget.","supporting_citations":[],"review_version":1}