{"id":"c9a3f91e-077a-4003-818a-91decd75112c","arxiv_id":"1908.09200","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"MUSTANG-2 SZ observations of Zwicky 3146 recover a pressure profile in agreement with X-rays and yield M500 = 8.16 x 10^14 solar masses from a Y-M scaling relation.","lead":"A 90 GHz Sunyaev-Zel'dovich map of the galaxy cluster Zwicky 3146, made with the MUSTANG-2 camera on the Green Bank Telescope, yields a cluster mass of about 8 x 10^14 solar masses. The measurement demonstrates a new data-processing pipeline that recovers cluster pressure profiles beyond the telescope field of view, sharpening comparison of cluster mass measurement methods.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The outermost non-parametric pressure bin is plausibly biased low by ~2.1σ (Appendix B); correcting it shifts M500 estimates by up to ~20%, yet Section 7 leaves this systematic unquantified, so the fiducial mass error bar is understated.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the outermost non-parametric pressure bin is assumed unbiased, but the paper's own Appendix B shows it is plausibly biased low and that this biases the mass estimates. I agree this is the central issue. The concern is load-bearing because both the headline pressure-profile agreement and the fiducial mass derive from the NP profile, and the outer bin contributes to the integrated Y and the hydrostatic mass profile. The paper is transparent: Section 7 explicitly lists an unquantified extra systematic, and Appendix B quantifies the sensitivity with Table 5. This is not an internal inconsistency or a disagreement with consensus; it is a missing uncertainty term in an otherwise careful measurement. The conditional verdict is appropriate: the measurement is accept-shaped, but the archived precision on M500 cannot be taken at face value until the outer-bin bias is either refuted by an injection test or folded into the error budget. If the injection test clears the bin, the concern is resolved; if it confirms the bias, the fiducial mass error bar must be enlarged and the central claim revised. No change to the reader's verdict is needed.","tokens_in":32997,"tokens_out":4449,"duration_ms":46052,"concrete_test":"Run an end-to-end injection test: add a known gNFW pressure profile with outer-bin pressure anchored to the XMM-Newton profile into the MUSTANG-2 TODs, reduce with Minkasi using the same binning and noise model, and check whether the recovered outer-bin value is biased low by ~2.1σ. If the bias reproduces, apply the empirically measured transfer correction to the real data and recompute Table 5 masses with full covariance; if any M500 shift exceeds the originally quoted total uncertainty, the fiducial error bar must be expanded accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central mass claim rests on the non-parametric pressure profile, whose outermost bin (207″ to infinity, Section 4.2) is weakly constrained and appears biased low. Appendix B shows this explicitly: requiring the entropy profile to continue the fitted r^1.34 power law raises the outer-bin pressure by a factor of 2.3, a ~2.1σ shift. Table 5 then shows that this single-bin correction moves M_HE from 8.29 to 9.95×10^14 M☉ and the fiducial (Y-M)_A10 estimate from 8.16 to 8.52×10^14 M☉, with even larger excursions for the M12 relation. Section 7 concedes that this 'additional systematic error' cannot yet be quantified. Because Y_sph(R500) and the hydrostatic mass profile (Eq. 10) both integrate over this bin, the incompleteness of the reported error budget is not a peripheral issue: the 5.5% statistical and 7% calibration figures do not capture a systematic that can move M500 by up to ~20% in an unfavourable but physically motivated direction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents MUSTANG-2 90 GHz Sunyaev-Zel'dovich observations of the relaxed cool-core cluster Zwicky 3146. Using a new maximum-likelihood processing pipeline (Minkasi) that fits surface-brightness annuli directly to time-ordered data, the authors derive non-parametric and gNFW/A10 pressure profiles over 5\" to 300\", compare them with an XMM-Newton pressure profile, and estimate M500 and M2500 through three routes: Y-M scaling relations (A10, M12, P17, C11), the virial theorem, and hydrostatic equilibrium with an external X-ray electron density profile. The fiducial Y-M mass is 8.16(+0.44,-0.54)(+0.46,-0.43)(+0.59,-0.55)x10^14 M_sun, and the hydrostatic mass is 8.29(+1.93,-1.24)(+0.74,-0.68)x10^14 M_sun. The paper also investigates residuals, temperature, entropy, and gas fraction profiles, and it explicitly analyzes the possibility that the outermost pressure bin is biased low.","tokens_in":33296,"tokens_out":4136,"duration_ms":45165,"significance":"If the central mass estimates held at their quoted precision, this would be a valuable demonstration that a ground-based SZ instrument can constrain a cluster pressure profile non-parametrically beyond its radial field of view, and it would provide a useful multi-method mass comparison for a relaxed cluster that has a wide spread of literature masses. The paper's strengths are its detailed treatment of the Minkasi covariance matrix, the MCMC fitting, the dual-pipeline (MIDAS/Minkasi) comparison, the careful point-source modeling, and the unusually honest Appendix B, which tests the outermost pressure bin against a physically motivated entropy-continuation prior. However, the same appendix shows that this single bin shifts mass estimates by up to about 20%, and Section 7 admits that the associated systematic error cannot yet be quantified. Because the fiducial mass and the hydrostatic mass both integrate over this bin, the quoted error budget is incomplete and the central mass claims are not yet established at the stated precision.","major_comments":[{"comment":"The outermost non-parametric pressure bin, spanning 207 arcseconds to infinity, is load-bearing for the integrated quantities Y_sph(R500) and the hydrostatic mass profile (Eqs. 6, 7, and 10), yet its value is weakly constrained. Appendix B shows that requiring the entropy profile to continue the fitted r^1.34 power law raises this bin's pressure by a factor of 2.3, a 2.1-sigma shift, and Table 5 reports that this shift moves M_HE from 8.29 to 9.95 x 10^14 M_sun and the fiducial (Y-M)_A10 estimate from 8.16 to 8.52 x 10^14 M_sun, with larger excursions for some other relations. Section 7 concedes that the associated systematic error is currently unquantifiable. As a result, the reported 5.5% statistical and 7% calibration error bars on the fiducial mass do not include a systematic that can change M500 by roughly 20% in a physically motivated direction. The authors should either quantify this systematic with a prior-based analysis or restate the central mass as a range that encompasses the entropy-continuation adjustment.","section":"Sections 4.2, Appendix B, Table 5"},{"comment":"The paper finds that the hydrostatic mass is larger than the Y-M masses, implying a negative hydrostatic mass bias, in direct tension with the expected positive bias of 0.1-0.3. This is not a peripheral discrepancy: it is a central consistency check for the three mass estimators. Appendix B shows that the adjustment that resolves the entropy turnover makes the bias more negative (b = -0.17 to -0.48 for the A10 relation), so the tension is not cured by the apparently plausible outer-bin correction. The manuscript should directly examine whether the negative bias points to a bias in the MUSTANG-2 pressure profile, a bias in the XMM electron density profile, or a problem in the adopted Y-M relations. As written, Section 7's statement that there is an additional problem to be solved leaves the origin unresolved and weakens the paper's claim that the mass estimates are self-consistent.","section":"Sections 5.2, 6.3, and 7"},{"comment":"The hydrostatic mass estimate uses an external X-ray electron density profile, but the quoted errors for M_HE include only statistical and flux-calibration terms; uncertainties in the XMM density profile (background subtraction, absolute calibration, and possible clumping) are not propagated. Given that the hydrostatic mass is one of the headline results and is compared against the Y-M masses, the absence of a ne-related systematic term makes the error budget incomplete. The authors should either quote M_HE with a density-profile systematic term or explicitly present the hydrostatic estimate as a consistency check rather than a primary mass measurement.","section":"Section 4.3.2, Eq. (10)"}],"minor_comments":[{"comment":"The abstract states that the SZ-derived pressure profile is in excellent agreement with the X-ray pressure profile, but Section 6.2 and Appendix B describe conflicts in derived products such as entropy and hydrostatic mass. I recommend softening the abstract or quantifying the agreement with a goodness-of-fit statistic.","section":"Abstract and Section 6.2"},{"comment":"The relativistic SZ correction assumes kBTe = 7 keV, which is reasonable for this cluster, but the sensitivity of the fitted pressure profile and masses to this assumed temperature is not stated. A short statement that the effect is negligible would be sufficient.","section":"Section 4.2"},{"comment":"The gNFW entropy slope is quoted with an uncertainty of 4e-5, which is artificially small because the pressure shape parameters are fixed to the A10 values. The text should note that this uncertainty does not include shape-parameter or calibration systematics.","section":"Section 5.4.2"},{"comment":"The MIDAS pressure profile shows its two outermost points about 2 sigma below the Minkasi and X-ray profiles, which reinforces the outer-bin concern in Appendix B. Cross-referencing this behavior in Section 4.2 or Appendix B would make the systematic issue easier for the reader to track.","section":"Appendix C"},{"comment":"There are several typographical and wording issues, including 'quadropole' (Section 5.3), 'The conversion is comes from' (Section 4.2), and 'completeleness' (Appendix C). I recommend a careful proofreading pass.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"This is a technically careful paper whose central issue is the unquantified systematic associated with the outermost pressure bin. The authors are commendably transparent about it, and Appendix B provides a concrete test, but the mass estimates and the paper's central claims are sufficiently sensitive to this bin that the error budget must be revised before publication. I do not see a novelty or citation-practice concern; the MUSTANG-2/Minkasi pressure-profile demonstration is valuable. I recommend major revision rather than rejection because the issue appears addressable within the manuscript's scope, for example by adopting the entropy-continuation profile as a systematic envelope or by re-framing the fiducial mass as a range."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a careful single-cluster SZ measurement with a real methodological payoff. The Minkasi pipeline fits annuli directly to the timestreams, avoids the high-pass filtering that limits the older MIDAS pipeline, and recovers the pressure profile to about twice the field of view. That beyond-FOV recovery is the genuinely new capability, and the internal checks — agreement with XMM, comparison with MIDAS, residual analysis — make the claim credible.\n\nThe measurement itself is also useful. Zw3146 has literature masses scattered from roughly 4 to 22 x 10^14 M_sun; the MUSTANG-2 data give a 104-sigma detection and a set of ICM-based masses that cluster around 8 x 10^14, which narrows the picture. The treatment of point sources, covariance, and calibration is above average for this literature.\n\nThe soft spot is real and it is in the paper. The outermost non-parametric bin, 207 arcsec to infinity, is weakly constrained. Appendix B shows that if the entropy profile is forced to continue its fitted power law, the bin rises by a factor of 2.3, a roughly 2.1-sigma shift; M_HE goes from 8.29 to 9.95 x 10^14, and the fiducial (Y-M)_A10 estimate goes from 8.16 to about 8.5, with bigger moves for other relations. Section 7 then says there is an additional systematic that cannot yet be quantified. That sentence is honest, but it also means the paper's headline 5.5% statistical and 7% calibration errors do not cover the dominant systematic, which can move M500 by about 20%. The negative hydrostatic mass bias is the same symptom: something about the outer pressure or the Y-M calibration is off. The paper speculates the problem may lie in the scaling relations, which is possible, but it does not resolve which side is wrong.\n\nI do not think this kills the paper. It is one cluster, and the central mass remains inside the broad literature range. But the error budget should be revised, or the caveat made much more prominent, before the precision claim is taken at face value.\n\nThis is a methods-plus-measurement paper for cluster cosmologists and SZ pipeline builders. It deserves a serious referee. I would send it out, and I would ask for the outer-bin systematic to be quantified or the quoted uncertainties widened accordingly.","headline":"Solid single-cluster SZ measurement with a genuinely useful pipeline advance; the mass is plausible, but the quoted error bar omits a known outer-bin systematic that can move M500 by roughly 20%.","tokens_in":33872,"tokens_out":3131,"would_cite":false,"duration_ms":34172,"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":"This paper shows that deep, high-resolution 90 GHz Sunyaev-Zel'dovich observations of the relaxed cluster Zwicky 3146 recover the intracluster pressure profile non-parametrically and pin its $M_{500}$ at $8.16^{+0.44}_{-0.54}$ (stat.)…","keywords":["galaxy clusters","Sunyaev-Zel'dovich effect","intracluster medium","pressure profiles","mass estimation","hydrostatic equilibrium","Zwicky 3146","Compton y parameter"],"falsifier":"Re-observe Zwicky 3146 with a larger scan radius (or use an independent X-ray pressure measurement beyond 200 arcseconds) and compare the recovered outermost pressure bin with the value assumed here. A rise of roughly a factor of 2.3, as the paper's entropy-power-law test predicts, would confirm the suspected bias and shift the hydrostatic mass estimate upward by about $1.7 \\times 10^{14}\\,M_\\odot$.","tokens_in":32834,"feed_emoji":"🔭","tokens_out":9193,"duration_ms":80461,"temperature":0.7,"pith_summary":"High-resolution Sunyaev-Zel'dovich (SZ) observations of the galaxy cluster Zwicky 3146, taken at 90 GHz to a central noise below 15 $\\mu$K, are used to recover the thermal pressure profile of the intracluster gas without assuming a fixed analytic shape. The paper's central claim is that a model-fitting data pipeline can constrain that profile beyond the telescope's field of view, and that the resulting pressure profile agrees with the one obtained from XMM-Newton X-ray data. From that profile the authors derive mass estimates through three independent routes -- $Y$-$M$ scaling relations, hydrostatic equilibrium, and the virial theorem -- obtaining $M_{500}$ values from $6.23 \\pm 0.59$ to $10.6 \\pm 0.95 \\times 10^{14}\\,M_\\odot$. Their fiducial mass, from the $Y$-$M$ relation, is $8.16^{+0.44}_{-0.54}$ (5.5% statistical), with 5.5% systematic uncertainty from the scaling relation and 7.0% from calibration, in units of $10^{14}\\,M_\\odot$. The result matters because published masses for this cluster scattered from $3.88$ to $22.5 \\times 10^{14}\\,M_\\odot$, and cluster-cosmology constraints hinge on exactly this kind of mass calibration.","feed_headline":"Deep SZ data pin Zwicky 3146 at 8.16 x 10^14 Msun","feed_subtitle":"High-resolution 90 GHz observations recover the gas pressure profile, matching X-rays and tightening a mass range that varied fivefold.","key_machinery":"The central object is the Compton $y$ parameter, the line-of-sight integral of thermal electron pressure, which is what the 90 GHz observations measure through the SZ temperature decrement. The argument is carried by three pieces: (1) a non-parametric model that bins pressure into twelve logarithmically spaced radii and fits those bins directly to the detector timestreams with a maximum-likelihood pipeline, recovering scales beyond the 4.25-arcminute field of view; (2) a generalized NFW pressure profile with A10 shape parameters as a parametric cross-check; and (3) three mass estimators -- integrated $Y_{\\rm sph}$ compared with published $Y$-$M$ relations, the hydrostatic equation $M_{\\rm HE} = -(d\\ln P_e/d\\ln r)\\,P_e\\,r/(n_e\\,\\mu m_p G)$ using X-ray electron densities, and a virial-theorem expression relating thermal energy to an NFW gravitational potential. All mass estimates are made self-consistent by finding where the derived mass curve crosses the reference $M_{500}(<r)$ curve.","core_discovery":"For a relaxed, cool-core cluster, the paper demonstrates that SZ data alone can deliver a deconvolved, non-parametric pressure profile spanning radii from about 5 arcseconds to beyond 300 arcseconds, with the outermost bin extending to infinity. That profile is statistically consistent with the pressure profile derived from XMM-Newton, and it anchors the cluster's total mass: the fiducial estimate $M_{500} = 8.16^{+0.44}_{-0.54}$ (stat.) $^{+0.46}_{-0.43}$ (sys., $Y$-$M$) $^{+0.59}_{-0.55}$ (sys., calibration) $\\times 10^{14}\\,M_\\odot$. The paper also finds that the hydrostatic mass estimate, $8.29^{+1.93}_{-1.24}$ (stat.) $^{+0.74}_{-0.68}$ (sys., calibration) $\\times 10^{14}\\,M_\\odot$, sits above the $Y$-$M$ masses, implying a negative hydrostatic mass bias rather than the usual positive one, and it investigates whether residual SZ substructure or an underestimated outer pressure bin can explain the inconsistencies.","pith_inferences":["A direct extension, not made in the paper, is to apply the same timestream-fitting pressure recovery to a sample of relaxed clusters spanning a range of masses and redshifts; if the outer-bin low bias is generic, current SZ mass estimates from ground-based dishes could be systematically low and $Y$-$M$ calibrations would need revisiting.","The paper's suspected link between the residual SZ decrement and the radio minihalo could be tested with high-resolution spectral-index mapping: a thermal SZ component tracing the minihalo would indicate that sloshing redistributes pressure as well as gas density on small scales.","If the negative hydrostatic bias persists across a larger sample, the standard practice of applying a 10-30% positive correction for hydrostatic bias in cluster cosmology would have to be reexamined, directly affecting mass-function estimates."],"forward_implications":["If the non-parametric recovery beyond the field of view is unbiased, single-dish SZ instruments can measure cluster masses to a few percent without relying on X-ray or weak-lensing calibration for the pressure shape.","The fiducial mass of $\\sim 8.16 \\times 10^{14}\\,M_\\odot$ agrees with ACT, Planck, and X-ray based estimates near $8 \\times 10^{14}\\,M_\\odot$, sharpening the picture against weak-lensing estimates that prefer lower masses.","The negative hydrostatic mass bias, if real, suggests that non-thermal pressure support alone cannot explain the offset between hydrostatic and $Y$-$M$ masses for this cluster, and that calibration of the $Y$-$M$ relations themselves may be implicated.","The agreement between SZ and X-ray pressure profiles implies that high-resolution SZ data can supply the pressure side of thermodynamic products such as temperature, entropy, and gas fraction when X-ray or other data supply the density."],"supporting_citations":[{"why":"Supplies the A10 universal pressure profile shape and the fiducial $Y$-$M$ relation used to derive $M_{500}$.","marker":"Arnaud et al. 2010"},{"why":"Provides a weak-lensing-calibrated $Y$-$M$ relation used as one of the three $R_{500}$ mass estimators.","marker":"Marrone et al. 2012"},{"why":"Provides a simulation-calibrated $Y$-$M$ relation used for a second $R_{500}$ mass estimate.","marker":"Planelles et al. 2017"},{"why":"Provides the $Y$-$M$ relation used for the $R_{2500}$ mass estimates.","marker":"Comis et al. 2011"},{"why":"Derives the virial-theorem mass estimator with an external surface-pressure term that the paper applies to its pressure profile.","marker":"Mroczkowski 2011"},{"why":"Supplies the XMM-Newton electron density profile used in the hydrostatic equilibrium mass calculation.","marker":"Ghirardini et al. 2018"},{"why":"Provides the weak-lensing mass estimate that the paper compares against and finds in tension with ICM-based masses.","marker":"Okabe & Smith 2016"}],"fun_headline_variants":["MUSTANG-2 SZ data pin Zwicky 3146 at 8.16e14 Msun","High-res SZ pressure profile matches X-rays, trims mass range","104σ SZ detection yields tight cluster mass estimate","Negative hydrostatic bias hinted in Zwicky 3146 masses","SZ maps narrow Zwicky 3146 mass spread from fivefold"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the outermost non-parametric pressure bin, spanning 207 arcseconds to infinity, is unbiased; if that bin is low, the paper's own entropy-continuity test raises it by about 2.1 $\\sigma$ and pushes $M_{500}$ estimates, such as the hydrostatic one from $8.29$ to $9.95 \\times 10^{14}\\,M_\\odot$, substantially higher.","fun_headline_variants_meta":{"raw":{"variants":["MUSTANG-2 SZ data pin Zwicky 3146 at 8.16e14 Msun","High-res SZ pressure profile matches X-rays, trims mass range","104σ SZ detection yields tight cluster mass estimate","Negative hydrostatic bias hinted in Zwicky 3146 masses","SZ maps narrow Zwicky 3146 mass spread from fivefold"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000773,"raw_usage":{"total_tokens":3648,"prompt_tokens":1396,"completion_tokens":2252,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":1012,"completion_tokens_details":{"reasoning_tokens":2151}},"tokens_in":1012,"tokens_out":2252,"duration_ms":20066,"temperature":1.0,"reasoning_tokens":2151,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:18:28.337761+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observe Zwicky 3146 with a larger scan radius (or use an independent X-ray pressure measurement beyond 200 arcseconds) and compare the recovered outermost pressure bin with the value assumed here. A rise of roughly a factor of 2.3, as the paper's entropy-power-law test predicts, would confirm the suspected bias and shift the hydrostatic mass estimate upward by about $1.7 \\times 10^{14}\\,M_\\odot$.","supporting_citations":[{"cited_title":"2011, , 728, L35, 10.1088/2041-8205/728/2/L35","cited_arxiv_id":null,"evidence_quote":"Derives the virial-theorem mass estimator with an external surface-pressure term that the paper applies to its pressure profile."}],"review_version":1}