{"id":"e0147609-a335-4c62-a6a2-58c53f8e04f3","arxiv_id":"2411.19911","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Cluster abundance with weak-lensing masses sets the dark matter-dark radiation temperature ratio to ξ_DR < 0.098 at 95% credibility when combined with CMB and BAO data.","lead":"This paper uses the number of galaxy clusters seen by the South Pole Telescope, with masses measured from light bending by the Dark Energy Survey and Hubble, to test a dark matter model where some dark matter particles interact with a dark radiation fluid. The key measurement is an upper limit on how much interacting dark radiation can exist, roughly three times tighter than previous limits when combined with cosmic microwave background and galaxy data.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cluster-abundance limit on xi_DR rests on an unvalidated mapping from the linear IDM-DR power spectrum to the halo mass function, so the quoted 95% upper bounds may be biased without N-body tests.","rationale":"The paper is careful in several respects: it uses the same analysis framework as Bocquet et al. (2024), checks robustness to fIDM and a_dark in Appendix A, and consistently translates xi_DR into Delta_N_eff. The load-bearing assumption is not the linear-theory calculation itself, which is standard, but the translation from linear power to halo abundance and to weak-lensing masses. The 10% IDM fraction is small but not negligible; if that component is pressure-supported at cluster scales, the effective collapse density differs from what the total-matter sigma(M) encodes, and the CDM-calibrated Tinker/emulator HMF may be biased in a way that directly changes the predicted cluster counts. The same physics can alter halo concentrations and thereby the MWL-Mhalo relation used to set the absolute mass scale. Because the headline result is a tight upper limit on xi_DR, a systematic error in this step could plausibly move the limit by an amount comparable to the improvement from adding cluster data. This does not make the paper wrong, but it makes the central claim conditional on a simulation calibration that is not currently provided. The reader's verdict of CONDITIONAL is therefore appropriate, and no change in verdict is recommended beyond requiring or clearly caveating this validation.","tokens_in":27131,"tokens_out":8533,"duration_ms":91058,"concrete_test":"Run N-body simulations of the ETHOS n=0 model (e.g., a modified GADGET/Arepo implementation of the coupled IDM-DR fluid equations) for fIDM=0.1, log10[a_dark/Mpc^-1]=8, and xi_DR=0.1 and 0.2, with the same cosmology as the SPT cluster analysis. Measure the M200c halo mass function over 0.25<z<0.8 and the concentration-mass relation, and compare with the Tinker/emulator prediction computed from the CLASS linear P(k) used in Eq. (24). If the ratio deviates by more than the statistical precision of the 1,005-cluster sample, re-run the likelihood with the simulation-calibrated HMF and MWL-Mhalo relation; a shift in the xi_DR upper limit beyond ~0.01 would weaken the 'three times tighter' claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result—xi_DR<0.166 from SPT-clusters x WL and xi_DR<0.098 when combined with CMB+BAO—depends on Eq. (24), where the IDM-DR halo mass function is computed by feeding the CLASS linear matter power spectrum into the standard gravity-only Tinker/emulator HMF. This is an unvalidated universality assumption. In the tight-coupling limit (log10[a_dark/Mpc^-1]=8, fIDM=0.1), the 10% IDM component is pressure-supported by the DR fluid and may not collapse into halos on cluster scales; the total-matter sigma(M) used by the HMF is then not the correct collapse variable, and the mapping from P(k) to dN/dM can be biased even though the linear spectrum is computed correctly. The same issue affects the weak-lensing mass calibration: the MWL-Mhalo relation in Eq. (19) is calibrated to standard CDM/hydro simulations, but IDM-DR suppression changes the concentration-mass relation, so the mass scale entering the abundance could shift. No N-body or calibrated emulator test of this step is presented; Appendix A varies fIDM and a_dark within the same unvalidated HMF framework, so it does not test the assumption. If the true HMF differs from the Tinker+CLASS prediction at the ~10-20% level at M~3e14 h^-1 Mpc, the 95% upper limits on xi_DR would shift by an amount comparable to the claimed improvement from adding clusters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives cosmological constraints on the ETHOS n=0 interacting dark matter–dark radiation (IDM–DR) model using the abundance of 1,005 SPT clusters with weak-lensing mass calibration from DES Y3 and HST. The key result is an upper limit on the dark radiation temperature ratio, ξ_DR < 0.166 at 95% credibility from SPT clusters × WL alone, and ξ_DR < 0.098 at 95% credibility when combined with Planck 2018 CMB and BOSS BAO data. In an SU(3) realization this corresponds to ΔN_eff < 0.003 for tightly coupled, fluid-like dark radiation. The paper also examines the impact of massive neutrinos, verifies robustness to freeing f_IDM and a_dark (Appendix A), and reports a statistically insignificant preference for non-zero ξ_DR. The analysis pipeline is the same validated cluster-abundance and mass-calibration framework used in recent SPT/DES work, with the IDM–DR matter power spectrum computed by CLASS.","tokens_in":27473,"tokens_out":5228,"duration_ms":54942,"significance":"If the underlying modeling assumption holds, this is a significant result: it would be the tightest cluster-abundance constraint on interacting dark radiation and demonstrates that tSZE cluster counts with weak-lensing mass calibration provide a probe complementary to CMB and BAO for ETHOS-type dark sectors. The paper is careful and honest in reporting: the Poisson abundance likelihood is well established, the weak-lensing mass calibration is anchored to external hydrodynamical simulations, Appendix A explicitly tests freedom in adark and fIDM, and the neutrino-mass check in Sec. V.B is a useful robustness test. The main caveat is that the central bound is obtained by feeding the CLASS linear power spectrum into a gravity-only Tinker/emulator halo mass function without a dedicated validation for the IDM–DR model; this is a load-bearing assumption that needs to be addressed before the quoted limits can be taken at face value.","major_comments":[{"comment":"The central ξ_DR constraints assume that the halo mass function for the IDM–DR model is obtained by inserting the CLASS linear matter power spectrum into the standard gravity-only Tinker/emulator HMF. This is not obviously valid when 10% of the DM is tightly coupled to a dark radiation fluid: in that regime the IDM component is pressure-supported and may not participate in halo collapse on cluster scales, so the total-matter σ(M) is not an established collapse variable for this model. A 10–20% bias in dN/dM at M ~ 3×10^14 h^-1 Mpc could shift the 95% upper limits by an amount comparable to the improvement claimed from adding clusters to CMB+BAO. Appendix A varies adark and fIDM within the same unvalidated HMF framework, so it does not test this assumption. I request either a simulation-based validation of the HMF mapping (e.g., N-body or calibrated emulator tests) or an explicit conservative systematic term folded into the reported limits.","section":"Sec. IV.B–IV.C, Eq. (24)"},{"comment":"The MWL–Mhalo relation used for weak-lensing mass calibration is calibrated to standard CDM hydrodynamical simulations and assumed unchanged in the IDM–DR model. However, the IDM–DR power suppression alters the concentration–mass relation of collapsed halos, and the NFW-derived lensing mass bias MWL/M200c may differ from the CDM-calibrated value. Since the abundance likelihood is very sensitive to the absolute mass scale, a few-percent shift in this bias could bias ξ_DR and also contribute to the reported slight preference for non-zero ξ_DR. The paper should quantify the expected change in the mass calibration for IDM–DR halos, or at least provide an estimate of the resulting systematic uncertainty on the ξ_DR limits.","section":"Sec. IV.A, Eq. (19)"}],"minor_comments":[{"comment":"The abstract states that the combined limit is 'around three times tighter' while Sec. VI says the ξ_DR bound itself is '30% tighter'; both statements are correct because the DR energy density scales as ξ_DR^4, but the wording is easy to misread. Please state explicitly which quantity (ξ_DR versus ΔN_eff or ρ_DR) is improved by which factor.","section":"Abstract and Sec. VI"},{"comment":"The non-standard combination S8^opt = σ8(Ωm/0.3)^0.2 is defined in the figure caption and Table II, but it would help to define it in the main text at first use, since most readers will expect the standard S8 exponent of 0.5.","section":"Sec. V.A and Fig. 3"},{"comment":"The 'slight preference for a non-zero value of ξ_DR' is described as robust to several variations, but no Bayes factor or other model-comparison statistic is reported. A posterior mode away from zero is not by itself evidence for the model; please add a quantitative model comparison or explicitly label the statement as a posterior-only hint.","section":"Sec. V.A"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the observational analysis is careful, but the central claim rests on an unvalidated mapping from the IDM–DR linear power spectrum to the halo mass function and on a CDM-calibrated weak-lensing mass relation. Both are in principle addressable with simulations or with a conservative systematic term. I would not reject the paper on this basis, but the limits as currently quoted are not yet robust enough for publication without either validation or a clearly quantified caveat."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First real cluster-abundance constraint on the ETHOS n=0 IDM-DR model, carefully executed, but the headline xi_DR bound rests on an unvalidated mapping from linear power spectrum to halo mass function; treat the number as provisional.\n\nWhat is new is exactly what the title says: the first application of actual SPT cluster abundance data, with DES Y3 and HST weak-lensing mass calibration, to the ETHOS n=0 interacting dark matter-dark radiation model. The main result is real: xi_DR < 0.166 at 95% from clusters alone, and < 0.098 combined with Planck 2018 CMB and BOSS BAO, about three times tighter than CMB+BAO alone. The analysis is careful. The Poisson abundance likelihood is standard, the weak-lensing mass calibration is empirical and anchored to hydro simulations, and the robustness checks in Appendix A (freeing fIDM and a_dark) plus the neutrino check in Sec. V.B are honest and useful. The data handling follows the companion papers, and the citations to the ETHOS model and previous constraints look appropriate. They also label the slight preference for non-zero xi_DR as not statistically significant, which is the right call.\n\nThe soft spot is the one the stress-test flags. Equation (24) computes the IDM-DR halo mass function by feeding the CLASS linear total-matter power spectrum into the standard gravity-only Tinker/emulator HMF. That presupposes universality. With fIDM = 0.1 tightly coupled to the DR fluid, the interacting fraction is pressure-supported and may not collapse into halos on cluster scales; the total-matter sigma(M) may not be the right collapse variable. The MWL-Mhalo relation in Eq. (19) is calibrated on standard CDM/hydro simulations, and if the concentration-mass relation shifts, the mass scale entering the abundance shifts too. No N-body test of this step is presented, and Appendix A does not test it because it varies fIDM and a_dark within the same HMF framework. If the true HMF differs at the 10-20% level near the pivot mass, the quoted upper limits could move by an amount comparable to the improvement from adding clusters.\n\nThis is a caveat, not a fatal flaw. The direction and size of any bias are unknown without simulations. The qualitative conclusion that cluster abundance with weak-lensing mass calibration is a competitive probe for this class of models would likely survive even if the mapping needs correction. But the quantitative limit should carry that uncertainty.\n\nThe paper is for cluster cosmologists, dark sector model builders, and anyone tracking the S8 tension. It deserves a serious referee. I would send it to review, with the HMF validation as the explicit central request: either run N-body tests or clearly delimit where the universality assumption is reliable.","headline":"First real cluster-abundance limit on ETHOS n=0, carefully done, but the quoted xi_DR bound depends on an unvalidated halo mass function mapping that could shift the result.","tokens_in":28371,"tokens_out":6376,"would_cite":true,"duration_ms":55851,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using 1,005 South Pole Telescope clusters with DES and HST weak-lensing masses, the paper bounds the interacting dark-radiation temperature ratio to below 10% when combined with CMB and BAO data.","keywords":["interacting dark matter","dark radiation","ETHOS","galaxy cluster abundance","weak lensing mass calibration","Sunyaev-Zeldovich effect","dark sector cosmology","S8 tension"],"falsifier":"A cosmological N-body simulation suite evolving the ETHOS $n=0$ model with $f_{\\rm IDM}=10\\%$ and tight coupling, comparing its halo mass function to the Tinker/emulator prediction used in Eq. (24), would settle the central claim: if cluster counts differ by more than the analysis's systematic error budget, the $\\xi_{\\rm DR}$ upper limits are biased.","tokens_in":26959,"feed_emoji":"🔭","tokens_out":10895,"duration_ms":91430,"temperature":0.7,"pith_summary":"This paper asks how much dark matter can interact with a bath of dark radiation without changing the number of galaxy clusters that form. Using 1,005 Sunyaev-Zeldovich-selected clusters from the South Pole Telescope with masses calibrated by DES Year 3 and HST weak-lensing data, the authors compare the observed abundance to the predictions of an interacting dark-matter–dark-radiation model (ETHOS $n=0$). They obtain an upper limit on the dark-radiation-to-CMB temperature ratio of $\\xi_{\\rm DR}<0.166$ at 95% credibility from clusters alone, tightening to $\\xi_{\\rm DR}<0.098$ when Planck 2018 CMB and BOSS BAO data are added, roughly three times tighter than CMB+BAO alone. If correct, this would be the strongest cluster-based constraint on interacting dark radiation and would leave room for such interactions to ease the $S_8$ tension, with a weak preference for a small nonzero $\\xi_{\\rm DR}$.","feed_headline":"Cluster counts cap dark radiation temperature at 10 percent","feed_subtitle":"Galaxy clusters from SPT with DES and HST weak-lensing masses set the tightest limit yet on interacting dark radiation.","key_machinery":"The load-bearing machinery is the ETHOS perturbation system for a tightly coupled dark-matter–dark-radiation fluid, in which the Euler equation for the interacting dark matter carries a drag term $\\Gamma_{\\rm IDM-DR}(\\theta_{\\rm IDM}-\\theta_{\\rm DR})$ with $\\Gamma_{\\rm IDM-DR}(z)\\propto(1+z)$ for $n=0$, while the dark radiation is treated as a self-interacting fluid with zero shear stress. This drag damps density perturbations on a wide range of scales, with $\\xi_{\\rm DR}$ setting where the suppression begins. The cluster abundance likelihood of Eq. (24) then maps the halo mass function, computed from the CLASS linear power spectrum for this model, into predicted counts in tSZE significance, richness, weak-lensing shear, and redshift, with the observable–mass relations calibrated empirically by DES Y3 and HST weak-lensing data.","core_discovery":"The paper's central claim is that weak-lensing-calibrated galaxy cluster abundance is sharply sensitive to the temperature ratio $\\xi_{\\rm DR}=T_{\\rm DR}/T_{\\rm CMB}$ in ETHOS $n=0$ interacting dark-matter–dark-radiation models, and that this sensitivity produces the tightest existing upper limit on that ratio. Within the microscopically motivated $SU(N)$ dark-sector realization, a fraction $f_{\\rm IDM}=10\\%$ of the dark matter is tightly coupled to a self-interacting dark-radiation fluid, and the interaction suppresses the matter power spectrum over a broad range of scales whose onset is set by $\\xi_{\\rm DR}$. The paper reports $\\xi_{\\rm DR}<0.166$ at 95% credibility from SPT clusters with DES Y3 and HST masses alone, and $\\xi_{\\rm DR}<0.098$ when combined with Planck 2018 and BOSS BAO data, corresponding in an $SU(3)$ realization to $\\Delta N_{\\rm eff}<0.003$. It also finds $S_8=0.793\\pm0.032$ from clusters alone and a mild, statistically insignificant preference for nonzero $\\xi_{\\rm DR}$ that persists across analysis choices, which the authors interpret as leaving a physical route toward resolving the $S_8$ tension.","pith_inferences":["A direct test of the central assumption would be N-body simulations of the ETHOS $n=0$ model; if the simulated halo mass function deviates from the Tinker/emulator prediction by more than the analysis's systematic budget, the reported $\\xi_{\\rm DR}$ limits would be biased.","Combining cluster abundance with cosmic-shear measurements should break the degeneracy between $\\xi_{\\rm DR}$ and the interaction parameters $f_{\\rm IDM}$ and $a_{\\rm dark}$, which cluster counts alone cannot constrain.","An independent cluster sample selected at other frequencies or wavelengths could confirm or refute the mild preference for nonzero $\\xi_{\\rm DR}$, since that preference is currently based on one dataset combination.","If the preference is physical, the same model predicts a characteristic scale-dependent shape change in the halo mass function that future lensing-calibrated cluster surveys should detect as a deviation from $\\Lambda$CDM at the few-percent level."],"forward_implications":["SPT cluster counts alone constrain $\\xi_{\\rm DR}<0.166$ at 95% credibility, demonstrating that cluster abundance is a competitive standalone probe of interacting dark radiation.","Adding Planck 2018 and BOSS BAO data tightens this to $\\xi_{\\rm DR}<0.098$, which for a minimal $SU(3)$ dark sector corresponds to $\\Delta N_{\\rm eff}<0.003$, about three times stronger than the CMB+BAO-only bound.","Cluster abundance is most sensitive to scales $k\\in[0.02,0.08]\\,h/\\mathrm{Mpc}$, making it complementary to small-scale probes for testing the full ETHOS parameter space.","The measured $S_8=0.793\\pm0.032$ is consistent with both Planck and cosmic-shear measurements, and the posterior retains a tail toward lower values, leaving room to address the $S_8$ tension.","Marginalizing over neutrino mass shows no degeneracy with $\\xi_{\\rm DR}$ and yields $\\sum m_\\nu<0.096$ eV from clusters+CMB+BAO, tighter than CMB+BAO alone."],"supporting_citations":[{"why":"Supplies the SPT cluster abundance and DES/HST weak-lensing analysis framework and the Lambda CDM comparison.","marker":"[5]"},{"why":"Provides the forecast and sensitivity argument motivating cluster abundance as a probe of xi_DR and the tight-coupling setup.","marker":"[46]"},{"why":"Derives the SU(N) microscopic model mapping and the interaction-rate formula used for the ETHOS n=0 setup.","marker":"[44]"},{"why":"Previous Planck+BAO analysis of IDM-DR that defines the model parameter space and prior bounds.","marker":"[42]"},{"why":"Documents the cluster cosmology methods, likelihood, and mass-calibration pipeline used here.","marker":"[60]"},{"why":"Calibrates the weak-lensing mass-halo mass relation that turns WL profiles into halo masses.","marker":"[83]"},{"why":"Provides the HST high-redshift weak-lensing data for 39 clusters.","marker":"[54]"},{"why":"Supplies the universal gravity-only halo mass function that converts the IDM-DR power spectrum into predicted cluster counts.","marker":"[84]"},{"why":"Provides the Planck 2018 CMB likelihood used in the combined analysis.","marker":"[7]"},{"why":"Provides the BOSS DR12 BAO measurements used in the combined analysis.","marker":"[8]"}],"fun_headline_variants":["Cluster counts cap dark radiation temperature at 10%","Dark radiation temperature limited to 10% by clusters","SPT clusters set tightest limit on interacting dark radiation","Cluster data bound dark radiation to 10% of CMB temperature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on assuming that the dark-matter–dark-radiation power suppression maps directly onto the standard gravity-only halo mass function, a mapping that this paper does not check with simulations.","fun_headline_variants_meta":{"raw":{"variants":["Cluster counts cap dark radiation temperature at 10%","Dark radiation temperature limited to 10% by clusters","SPT clusters set tightest limit on interacting dark radiation","Cluster data bound dark radiation to 10% of CMB temperature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000335,"raw_usage":{"total_tokens":1981,"prompt_tokens":1193,"completion_tokens":788,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":809,"completion_tokens_details":{"reasoning_tokens":721}},"tokens_in":809,"tokens_out":788,"duration_ms":7827,"temperature":1.0,"reasoning_tokens":721,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:41:58.429715+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A cosmological N-body simulation suite evolving the ETHOS $n=0$ model with $f_{\\rm IDM}=10\\%$ and tight coupling, comparing its halo mass function to the Tinker/emulator prediction used in Eq. (24), would settle the central claim: if cluster counts differ by more than the analysis's systematic error budget, the $\\xi_{\\rm DR}$ upper limits are biased.","supporting_citations":[{"cited_title":"Sunyaev-Zel'dovich Effect and X-ray Scaling Relations from Weak-Lensing Mass Calibration of 32 SPT Selected Galaxy Clusters","cited_arxiv_id":"1711.05344","evidence_quote":"Supplies the universal gravity-only halo mass function that converts the IDM-DR power spectrum into predicted cluster counts."}],"review_version":1}