{"id":"483e9bd8-9595-44ab-a425-f8a0f48d8c18","arxiv_id":"2507.10857","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":15,"one_line_summary":"A triaxial model combining X-ray, SZ, and weak lensing data finds Abell 1689 is elongated along the line of sight and yields M200c = 13.7e14 Msun, lower than the spherical value of 17.8e14 Msun.","lead":"This paper adds weak lensing data to an existing triaxial fitting pipeline that uses X-ray and microwave observations, and applies it to the galaxy cluster Abell 1689. It finds the cluster is elongated mostly toward us, which lowers the inferred mass by about 30 percent relative to a spherical fit.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Quoted RLP uncertainty (0.02) excludes documented systematic shifts up to 0.20, so the claimed 30% triaxial mass reduction is not yet robust.","rationale":"The central claim is that the triaxial fit yields a WL mass roughly 30% lower than the spherical fit because Abell 1689 is elongated along the line of sight. The decisive quantity is RLP, because the mass ratio tracks it: 17.77/13.69 ≈ 1.30 versus RLP = 1.27. The paper quotes RLP = 1.27 ± 0.02, but Appendix A is an internal record of shifts up to 0.20 from plausible alternative treatments, notably the one-dimensional versus two-dimensional X-ray temperature weighting. The final RLP is stable only by cancellation of these shifts, and cancellation does not remove the uncertainty in which treatment is correct. This directly affects the headline mass comparison and the word 'significantly.' I did not select the qpot = qICM assumption as the primary concern even though it is real: the mock tests in Section 4 generate clusters from the same model and therefore cannot validate that assumption against realistic potential shapes. The spherical fit being WL-only is less decisive, because in this pipeline the gas observables mainly constrain geometry and gas profiles, while M200c is constrained by WL in both fits. The reader's CONDITIONAL verdict already captures the need for revision, so this stress-test does not change it.","tokens_in":34792,"tokens_out":10118,"duration_ms":129383,"concrete_test":"Re-run the Abell 1689 fit exactly as in Section 6 but with the K24 two-dimensional X-ray temperature map instead of the one-dimensional projected profile from Rossetti et al. (2024), keeping priors, masks, covariances, and all other data fixed. Record RLP and M200c,triax. If RLP changes by at least 0.1 or M200c shifts by more than the quoted 68% credible interval, the headline claim requires a systematic error term and the reader's CONDITIONAL verdict stands as the appropriate outcome.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix A documents that plausible alternative treatments move RLP by amounts far exceeding the quoted statistical uncertainty: switching from the one-dimensional projected X-ray temperature profile to the two-dimensional temperature map changes RLP by about 0.20; the ARF choice by 0.07; the X-ray SB profile method by 0.05; and integration tolerance by 0.03. Some listed updates (APEC normalization, vignetting) are corrections of earlier errors and may not represent remaining uncertainty, but the temperature-weighting and ARF choices are genuine modeling ambiguities. The final RLP is stable only because these shifts partially cancel; cancellation among modeling choices is not a systematic error budget. Because the mass comparison is driven by RLP, with M200c,sph/M200c,triax = 17.77/13.69 ≈ 1.30 close to RLP = 1.27, an RLP error of ±0.2 changes M200c,triax by roughly ±2×10^14 Msun, comparable to or larger than the quoted ±1.5×10^14 uncertainty. The mock validation in Section 4 does not cover this, because the mocks are generated from the same model choices and do not vary temperature weighting or instrument response. A secondary assumption, qpot = qICM in Section 2.3.1, is also not stress-tested by those mocks, but the RLP systematic is the more direct threat to the stated mass comparison.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an extension of the triaxial cluster modeling pipeline of Kim et al. (2024) to include weak-lensing shear data, and applies it to Abell 1689 using X-ray surface brightness and temperature, SZ maps from Planck/ACT, and Subaru weak-lensing maps. The main results are an inferred line-of-sight elongation R_LP = 1.27 +/- 0.02, a triaxial mass M_200c = (13.69^{+1.56}_{-1.41}) x 10^14 Msun, a concentration c_200c = 8.55^{+2.20}_{-1.61}, and a non-thermal pressure fraction of roughly 20-30% between 0.18 and 1.37 Mpc. The paper reports that the spherical-fit mass is (17.77^{+2.00}_{-1.75}) x 10^14 Msun and interprets the difference as a projection bias. A mock-data validation of the mass and concentration recovery is included.","tokens_in":35123,"tokens_out":5440,"duration_ms":62408,"significance":"If the central claims hold, the paper would provide a valuable demonstration that joint triaxial modeling of X-ray, SZ, and WL data can remove orientation-dependent bias in weak-lensing mass estimates of individual clusters. The work uses a rich, publicly available CHEX-MATE dataset, documents its methodological updates transparently in Appendix A, and includes mock validation with full instrumental response. However, several load-bearing aspects of the analysis currently prevent full confidence in the headline results: the spherical comparison is not apples-to-apples, the concentration measurement is prior-dominated, and the elongation parameter carries unquantified modeling systematics that are comparable to the quoted statistical uncertainties.","major_comments":[{"comment":"The abstract states that the spherical fit 'otherwise employs the same methodology' as the triaxial fit, but Section 6 says the spherical fit was performed 'to the WL data' only. The triaxial mass, in contrast, comes from a joint X-ray, SZ, and WL fit. The comparison of M_200c = 17.77 vs 13.69 therefore conflates a change in geometry with a change in the data used in the fit, and the claim that the mass difference is due to line-of-sight elongation is not established by this comparison. Please either rerun the spherical fit on the full multi-probe dataset or state clearly that the spherical reference is WL-only and discuss the implications.","section":"Abstract and Section 6"},{"comment":"The concentration measurement c_200c = 8.55 is dominated by the informative log-normal prior centered on the Diemer & Joyce (2019) mass-concentration relation. The paper itself states in Section 6 and Figure 5 that the upper edge of the posterior is set by the prior and that a flat prior returns a posterior that hits the prior boundary. The mock tests in Section 4 (Table 2) also show a measurable prior-induced bias for input values c_200c = 2 and 10. The conclusion that the high concentration is intrinsic and 'not due to triaxiality and orientation' is therefore not supported by the data. Please present the flat-prior result and limit the concentration claim to what the data actually constrain.","section":"Sections 2.5.1, 4, and 6"},{"comment":"The quoted elongation R_LP = 1.27 +/- 0.02 is a statistical uncertainty only. Appendix A documents that alternative modeling choices shift R_LP by up to 0.20 (1D vs 2D temperature profile), 0.18 (APEC normalization), 0.07 (ARF), 0.05 (X-ray SB profile method), and 0.03 (integration tolerance). Because the mass ratio M_200c,sph / M_200c,triax ~ 1.30 closely tracks R_LP, an unquantified R_LP systematic of +/-0.2 changes the triaxial mass by roughly +/-2 x 10^14 Msun, comparable to or larger than the quoted mass uncertainty. The mock validation in Section 4 does not address this, since the mocks are generated with the same modeling choices. Please provide a systematic error budget for R_LP or weaken the claim that the triaxial mass is 'significantly lower' than the spherical value.","section":"Appendix A and Section 6"},{"comment":"The assumption that the gravitational potential has exactly the same axial ratios and orientation as the ICM (q_pot = q_ICM, co-aligned) is load-bearing for the triaxial mass estimate, but it is not stress-tested. The mock observations in Section 4 were generated from the same model, so they cannot validate this assumption. If the potential is rounder than the gas or is misaligned, the projected WL signal and the inferred M_200c will change. A sensitivity test with a rounder or misaligned potential, or at least a quantitative discussion of the expected bias based on simulations, is needed before the mass comparison can be considered robust.","section":"Section 2.3.1"}],"minor_comments":[{"comment":"There is a typo 'Custer Lensing And Supernovae survey' and 'hearafter' should be 'hereafter'.","section":"Introduction"},{"comment":"The sentence 'A grid of values for Lambda(Te,Z) is precalculated using using the Python package pyproffit' contains a duplicated 'using'.","section":"Section 2.2.2"},{"comment":"The phrase 'an the initial spherically-symmetric rho(r)' is ungrammatical and should be revised.","section":"Section 2.3.1"},{"comment":"The word 'Arbritary' before Eq. (17) should be 'Arbitrary'.","section":"Section 2.4"},{"comment":"The sentence 'We see that for the Compton-y and X-ray SB profiles, the model follows the the data closely' contains a duplicated 'the'.","section":"Section 6"},{"comment":"The heading 'Logrithmically-Spaced LOS Projection Integral' should be 'Logarithmically-Spaced'.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a credible and potentially useful pipeline, and the authors are transparent about their modeling changes in Appendix A. The main issues are internal consistency (the spherical comparison is not what the abstract claims), a prior-dominated concentration result that the text itself acknowledges, and an unquantified systematic error on the elongation parameter. These are addressable with additional analysis and careful rewriting, so I recommend major revision rather than rejection. If the authors can provide an apples-to-apples spherical comparison and a systematic error budget for R_LP, the paper would be a valuable contribution to the cluster-lensing literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: the headline result—triaxial mass about 30% lower than spherical—is not yet robust. The quoted RLP = 1.27 ± 0.02 carries only statistical error, but Appendix A shows that switching between 1D and 2D temperature weighting moves RLP by ~0.20, the ARF choice by ~0.07, and several other modeling choices by similar amounts. The shifts happen to cancel in the final value, but cancellation is not a systematic error budget. Since the mass ratio 17.77/13.69 ≈ 1.30 is essentially RLP, a ±0.2 uncertainty in RLP translates into a ±2 × 10^14 M_sun change in M200c—comparable to or larger than the quoted uncertainty. The mock validation doesn't cover this, because the mocks use the same modeling choices.\n\nWhat the paper does well: it integrates WL into the K24 gas-only triaxial formalism, applies it to A1689, validates mass recovery on mocks, and produces a non-thermal pressure profile consistent with Chappuis et al. and simulations. The comparison with prior triaxial analyses is thorough. This is a real pipeline that will be used on the CHEX-MATE sample.\n\nThe fixable problems: (1) The abstract says the spherical fit 'otherwise employs the same methodology,' but Section 6 shows it uses WL only. That is misleading. (2) The concentration result is prior-dominated—the log-normal prior centered on the Diemer & Joyce relation sets the posterior's upper edge, and a flat prior hits the boundary. So the claim that the high concentration is intrinsic is not supported. (3) The RLP systematic must be quantified, either by marginalizing over the weighting/response choices or by reporting an error an order of magnitude larger than 0.02.\n\nThe qpot = qICM assumption is a separate dependency that the mocks don't stress-test; it's secondary to the RLP issue.\n\nWho this is for: cluster cosmologists and anyone doing triaxial multi-probe fitting. It's a methods demonstration on one cluster, not a cosmological result. I'd send it to peer review, but the revision must address the systematic error budget and the spherical-fit comparison honestly.","headline":"The 30% triaxial mass reduction is not yet robust because the RLP systematic error is an order of magnitude larger than the quoted statistical uncertainty, but the pipeline itself is a solid methods contribution.","tokens_in":35726,"tokens_out":3175,"would_cite":true,"duration_ms":36178,"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":"Modeling Abell 1689 as elongated along the line of sight lowers its weak-lensing mass by about 30 percent relative to a spherical fit.","keywords":["galaxy clusters","triaxial modeling","weak gravitational lensing","Sunyaev-Zeldovich effect","X-ray astronomy","Abell 1689","non-thermal pressure","cluster mass bias"],"falsifier":"Measure the cluster's three-dimensional shape with an independent tracer that does not assume $q_{pot} = q_{ICM}$, for example strong-lensing arc positions combined with central-galaxy stellar kinematics, and check whether the inferred potential axial ratio matches the gas axial ratio. If the potential is significantly rounder, the 30 percent mass correction is an overestimate; alternatively, repeat the fit using the two-dimensional temperature map instead of the one-dimensional profile, which the appendix says shifts $\\mathcal{R}_{LP}$ by about 0.20, and see whether the mass difference disappears.","tokens_in":34582,"feed_emoji":"🌌","tokens_out":8809,"duration_ms":85235,"temperature":0.7,"pith_summary":"Galaxy clusters are triaxial, not spherical, so fitting them as spheres can bias the masses measured from weak lensing. This paper builds a multiprobe triaxial analysis that combines X-ray surface brightness and temperature, Sunyaev-Zeldovich effect, and weak-lensing shear in a single geometric model, and applies it to the massive cluster Abell 1689. It finds the cluster is stretched along the line of sight, with elongation $\\mathcal{R}_{LP} = 1.27 \\pm 0.02$, which lowers the weak-lensing mass from about $17.8 \\times 10^{14}$ to $13.7 \\times 10^{14}$ solar masses relative to an otherwise identical spherical fit. The paper argues that the cluster's unusually high concentration is not a projection artifact, and it derives a non-thermal pressure fraction rising from about 20 percent at 600 kpc to near 30 percent at the largest radius. If the approach is correct, spherical models systematically overestimate the masses of line-of-sight elongated clusters by tens of percent.","feed_headline":"Triaxial fit cuts Abell 1689's mass by ~30 percent","feed_subtitle":"The cluster is elongated toward us; spherical weak-lensing models overestimate its mass.","key_machinery":"The central object is the triaxial ellipsoidal radius $\\zeta$ defined by $\\zeta^2 = x_1^2/q_1^2 + x_2^2/q_2^2 + x_3^2$, with two axial ratios $q_1$ and $q_2$ and three Euler angles relating the ellipsoid to the observer. All ICM and gravitational-potential profiles are written as functions of $\\zeta$, and the projection integral $F_{2D}(x_\\xi) = 2 l_p e_\\parallel \\int_{x_\\xi}^\\infty F_{3D}(x_\\zeta) x_\\zeta / \\sqrt{x_\\zeta^2 - x_\\xi^2} \\, dx_\\zeta$ converts them into predicted two-dimensional maps of X-ray surface brightness, Compton-$y$, and lensing potential. The mass model works through the gravitational potential $\\Phi(\\zeta)$, which is assumed to have the same ellipsoidal shape and orientation as the X-ray gas; Poisson's equation is inverted numerically to connect $\\Phi$ to the total density, characterized by $M_{200c}$ and $c_{200c}$. The elongation parameter $\\mathcal{R}_{LP}$ (ratio of line-of-sight extent to average projected extent) sets the projection correction, and the generalized hydrostatic-equilibrium equation $dP_{tot}/d\\zeta = -\\rho_{gas}\\,d\\Phi/d\\zeta$ is integrated outside the fit to obtain the non-thermal pressure fraction.","core_discovery":"The paper's central claim is that combining X-ray, SZ, and weak-lensing observations with a triaxial ellipsoidal model removes a major projection bias in cluster mass measurement. For Abell 1689, the fit yields $M_{200c} = (13.69_{-1.41}^{+1.56}) \\times 10^{14}\\,M_\\odot$ with concentration $c_{200c} = 8.55_{-1.61}^{+2.20}$, whereas a spherical fit to the same weak-lensing data gives $M_{200c} = (17.77_{-1.75}^{+2.00}) \\times 10^{14}\\,M_\\odot$ and $c_{200c} = 9.99_{-1.78}^{+2.26}$. The mass difference is attributed to the cluster's line-of-sight elongation $\\mathcal{R}_{LP} = 1.27 \\pm 0.02$. Because the triaxial fit retains the high concentration, the paper concludes that the high concentration is intrinsic rather than a projection effect. The analysis also uses the triaxial gravitational potential, assumed to share the gas axial ratios, to derive the non-thermal pressure support, which runs from roughly 20 percent at intermediate radii to near 30 percent at the edges of the fitted range, at about $\\pm5$ percent precision. The pipeline was first tested on mock observations with known input parameters, and it recovered masses consistent with the inputs within noise.","pith_inferences":["If the potential is rounder than the gas, as some simulations suggest, the quoted mass correction is an upper limit; a rounder potential would produce a smaller triaxial-versus-spherical mass difference.","The appendix's sensitivity tests show individual modeling choices shift $\\mathcal{R}_{LP}$ by up to about 0.20, so the quoted statistical uncertainty of 0.02 likely understates the real uncertainty; a systematic shift that large would move the mass by several $\\times 10^{13}\\,M_\\odot$.","Applying this pipeline to a full sample should reveal a population-level correlation between inferred $\\mathcal{R}_{LP}$ and the spherical-versus-triaxial mass offset, which can be checked against independent shape indicators such as strong-lensing morphology."],"forward_implications":["Spherical weak-lensing fits overestimate the mass of line-of-sight elongated clusters by roughly 30 percent for a cluster like Abell 1689, and the size of the bias tracks the elongation $\\mathcal{R}_{LP}$.","The unexpectedly high concentration of Abell 1689 is intrinsic to the cluster rather than a viewing-angle artifact.","The multiprobe triaxial pipeline recovers input masses in mock observations, so it can be applied to the rest of the CHEX-MATE sample to map orientation biases.","Abell 1689 carries roughly 20–30 percent of its pressure in non-thermal form over 0.18–1.37 Mpc, with about $\\pm5$ percent precision."],"supporting_citations":[{"why":"Supplies the gas-only X-ray and SZ triaxial fitting formalism and the geometric definitions that this paper extends with weak lensing.","marker":"Kim et al. (2024)"},{"why":"Supplies the weak-lensing mass modeling methodology, including co-aligned constant axial ratios for the total matter distribution.","marker":"Sereno et al. (2017)"},{"why":"Provides the comparison analysis of the same cluster with the same raw data but a different modeling formalism, used to check mass, concentration, and non-thermal pressure.","marker":"Chappuis et al. (2025)"},{"why":"Provides the two-dimensional X-ray surface brightness map and the azimuthal-median radial profile used in the fit.","marker":"Bartalucci et al. (2023)"},{"why":"Provides the projected X-ray temperature profile used in the fit.","marker":"Rossetti et al. (2024)"},{"why":"Supplies the mass-concentration relation and intrinsic scatter that set the log-normal prior on c200c.","marker":"Diemer & Joyce (2019)"},{"why":"Provides a previous triaxial WL/SL analysis of Abell 1689 with flat priors, against which the paper compares its mass and concentration.","marker":"Umetsu et al. (2015)"}],"fun_headline_variants":["Triaxial fit trims Abell 1689's mass by a quarter","Cluster elongation toward us lowers mass estimate by ~23%","Spherical weak-lensing models overestimate Abell 1689's mass by 23%","3D cluster fit reduces Abell 1689's mass to 13.7e14 Msun","New triaxial analysis cuts Abell 1689's mass by a quarter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The mass and projection correction assume the dark-matter potential has exactly the same axial ratios as the X-ray gas and is aligned with it; if the potential is rounder or tilted, the inferred line-of-sight elongation and the mass shift change.","fun_headline_variants_meta":{"raw":{"variants":["Triaxial fit trims Abell 1689's mass by a quarter","Cluster elongation toward us lowers mass estimate by ~23%","Spherical weak-lensing models overestimate Abell 1689's mass by 23%","3D cluster fit reduces Abell 1689's mass to 13.7e14 Msun","New triaxial analysis cuts Abell 1689's mass by a quarter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00187,"raw_usage":{"total_tokens":7516,"prompt_tokens":1300,"completion_tokens":6216,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":916,"completion_tokens_details":{"reasoning_tokens":6109}},"tokens_in":916,"tokens_out":6216,"duration_ms":54307,"temperature":1.0,"reasoning_tokens":6109,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:24:31.516223+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the cluster's three-dimensional shape with an independent tracer that does not assume $q_{pot} = q_{ICM}$, for example strong-lensing arc positions combined with central-galaxy stellar kinematics, and check whether the inferred potential axial ratio matches the gas axial ratio. If the potential is significantly rounder, the 30 percent mass correction is an overestimate; alternatively, repeat the fit using the two-dimensional temperature map instead of the one-dimensional profile, which the appendix says shifts $\\mathcal{R}_{LP}$ by about 0.20, and see whether the mass difference disappears.","supporting_citations":[{"cited_title":"2024, , 686, A97","cited_arxiv_id":null,"evidence_quote":"Supplies the gas-only X-ray and SZ triaxial fitting formalism and the geometric definitions that this paper extends with weak lensing."},{"cited_title":"2017, Monthly Notices of the Royal Astronomical Society, 467, 3801–3826","cited_arxiv_id":null,"evidence_quote":"Supplies the weak-lensing mass modeling methodology, including co-aligned constant axial ratios for the total matter distribution."},{"cited_title":"CHEX-MATE: Multi-probe analysis of Abell 1689","cited_arxiv_id":"2503.22316","evidence_quote":"Provides the comparison analysis of the same cluster with the same raw data but a different modeling formalism, used to check mass, concentration, and non-thermal pressure."},{"cited_title":"2015, , 806, 207","cited_arxiv_id":null,"evidence_quote":"Provides a previous triaxial WL/SL analysis of Abell 1689 with flat priors, against which the paper compares its mass and concentration."}],"review_version":1}