{"id":"e6693f68-3393-4a0b-9201-31cfc1614586","arxiv_id":"2507.20938","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A revised escape velocity mass estimator yields weak lensing and escape masses for 46 clusters that agree within 0.02 dex, with a correlation of 0.68.","lead":"The authors compare two independent ways of weighing galaxy clusters, weak lensing and escape velocity, and find they now agree once the escape velocity measurement is corrected for sparse sampling. The result suggests previous disagreements came from an outdated caustic technique, not from missing cluster physics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The suppression model Z_v, the primary escape-mass systematic (Eq. 5-10), is calibrated on analytic spherical phase-spaces and validated only on gravity-only Millennium halos; a few-percent baryonic bias would shift escape masses by ~0.05-0.1 dex and could produce or erase the claimed 0.02 dex…","rationale":"The reader's weakest-assumption analysis points to the hydrodynamical validation gap for Z_v, and I agree that this is the most load-bearing concern. The paper's central claim is not simply that escape and lensing masses correlate, but that they agree with zero mean bias and scatter fully explained by individual errors. That quantitative statement inherits every systematic in the suppression model. The authors themselves call Z_v the primary systematic of the technique (Section 2) and explicitly note in Section 4.3.9 that they have not tested it in hydrodynamical simulations; the order-of-magnitude energy estimate they offer addresses wholesale galaxy displacement, not the subtle re-population of the phase-space edge by tides, harassment, or feedback-modulated orbits. Because M200 enters through v_esc^2, even a few percent error in the suppressed edge is a 0.05-0.1 dex mass error, comparable to the claimed 0.02 dex bias. The Millennium validation is genuine independent support—it shows the model is robust to asphericity, substructure, and interlopers in a gravity-only universe—but it does not cover baryonic physics. A hydrodynamical mock test would settle the issue. Since the concern is a missing calibration rather than a demonstrated contradiction, and since the paper's internal tests and the data comparison are otherwise coherent, the existing CONDITIONAL verdict is appropriate; I would not escalate to rejection or lower confidence without running the proposed test.","tokens_in":38933,"tokens_out":8151,"duration_ms":100845,"concrete_test":"Re-run the full pipeline (AGAMA Z_v table, MCMC likelihood, shifting-gapper interloper rejection, 5-bin monotonic edge extraction) on mock galaxy catalogs from a hydrodynamical simulation such as IllustrisTNG or The Three Hundred, matching the observed sample in N (median ~100; range 50-200), spectroscopic velocity errors (30 km/s), and cluster mass/redshift range. For 100+ independent lines of sight, compare inferred escape M200 to the true M200 and to mock weak-lensing masses. Settle the concern by checking (a) mean bias <0.03 dex with z-scores consistent with a unit Gaussian, and (b) the empirical Z_v distribution at fixed N matches the AGAMA-scaled skew-normal within a few percent. If either fails, the Section 4.2 concordance cannot be cleanly attributed to unbiased escape masses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—zero mean bias and no extra intrinsic scatter between escape and weak-lensing masses (Section 4.2, Figure 8)—rests on the accuracy of the suppression function Z_v in Equation 5. Z_v is the primary systematic of the technique (the authors' language in Section 2), yet it is calibrated purely from AGAMA phase-spaces of isolated, spherical Dehnen potentials (Section 3.1) and validated only against gravity-only Millennium halos with semi-analytic galaxies (Section 3.2). The paper explicitly does not test Z_v against hydrodynamical simulations containing baryonic feedback, MHD, star formation, and the galaxy-scale processes that set the high-velocity tail of the observed phase-space (Section 4.3.9). The substitution of an order-of-magnitude energy estimate (~10^60 erg to displace one galaxy) does not address the actual sensitivity of the measured edge: lower-energy processes such as tidal disruption, harassment, and orbital redistribution can change which galaxies populate the edge at fixed N without requiring a single galaxy to be displaced by AGN output. Because M200 scales roughly as v_esc^2 r, a 2-3% error in the suppressed edge translates into ~0.05-0.1 dex in mass, comparable to the quoted 0.02±0.02 dex bias. The Millennium test demonstrates insensitivity to gravity-only complexities (asphericity, substructure, interlopers), which is genuine but insufficient; baryonic effects remain uncalibrated. If Z_v is biased at the few-percent level, the observed concordance with weak lensing could be partly spurious or partly masked by a compensating lensing bias.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a revised comparison between escape-velocity masses and weak-lensing masses for 46 galaxy clusters at 0.05<z<0.3. The authors model the observed down-sampled escape edge with a skew-normal suppression function Z_v calibrated on AGAMA phase-spaces of spherical potentials, and they validate the mass-inference pipeline on Millennium N-body halos, finding unbiased recovery. Applied to clusters drawn from Herbonnet et al. (2020), Okabe & Smith (2016), and AS1063, the method yields a mean bias of 0.02±0.02 dex, a scatter of 0.17 dex, and a correlation of 0.679^{+0.046}_{-0.049} with weak-lensing masses, with no additional intrinsic scatter required. The authors contrast this with the prior caustic-based comparison, which showed negligible correlation and a 0.25 dex bias, and they provide a detailed discussion of systematic uncertainties.","tokens_in":39275,"tokens_out":7418,"duration_ms":90443,"significance":"If the result is correct, it is significant: it would demonstrate that a carefully modeled escape-edge technique yields cluster masses statistically consistent with independent weak-lensing measurements, overturning the previously reported lack of correlation and establishing escape masses as a viable independent mass probe with potential cosmological applications. The paper's strengths include the explicit treatment of the skewness of the suppression distribution, the end-to-end validation on Millennium halos with a range of dynamical states, the outlier analysis, and a clear systematics budget. However, the central claim is not yet fully secured: the suppression function Z_v is calibrated on gravity-only simulations and the paper explicitly does not test it against hydrodynamical simulations, so a few-percent model bias could shift the escape masses by an amount comparable to the quoted zero bias. The result is therefore promising but conditional on additional validation.","major_comments":[{"comment":"The suppression function Z_v, which the paper identifies as the primary systematic of the escape technique, is calibrated on AGAMA phase-spaces of isolated spherical potentials and validated only against gravity-only Millennium halos. The paper explicitly does not test Z_v against hydrodynamical simulations, and the order-of-magnitude energy argument (~10^60 erg to displace a single galaxy) does not bound the effect of lower-energy processes such as tidal stripping, harassment, or orbital redistribution, which can change which galaxies populate the phase-space edge at fixed N without requiring AGN displacement of a single galaxy. Because M200 scales roughly as v_esc^2 r, a 2–3% bias in the suppressed edge translates to ~0.05–0.1 dex in mass, comparable to the quoted bias of 0.02±0.02 dex. I recommend testing Z_v on a hydrodynamical simulation (e.g., IllustrisTNG or the Three Hundred project) or otherwise providing a quantitative demonstration that baryonic processes do not bias the edge population.","section":"§4.3.9 and Eqs. (5)–(10)"},{"comment":"The escape masses are not fully independent of the weak-lensing masses because the weak-lensing M200 is used to define the initial r200 that sets the radial binning and the phase-space count N entering Z_v. The AGAMA test in §3.5.3 shows the induced correlation is weak, and the Millennium test with 0.6 dex mass errors gives zero bias, which is reassuring. However, in the real sample, the weak-lensing masses are themselves the comparison quantity, and if those masses carry a correlated bias (e.g., from orientation or photometric-redshift systematics), part of that bias could propagate into the escape masses through the N-dependent suppression. The quoted 0.01 dex binning systematic in §4.3.2 does not obviously include this correlated component. I ask the authors to either propagate the full weak-lensing uncertainty through the escape-mass pipeline in a joint analysis or demonstrate with a realistic mock that the covariance is negligible at the 0.01 dex level.","section":"§4.2 and §3.5.3"},{"comment":"The claim that the 0.17 dex observed scatter requires no additional intrinsic component rests on z-score histograms in which the escape-mass errors appear consistent with σ=1 but the weak-lensing errors have σ≈1.4. Underestimated lensing errors directly weaken the inference about intrinsic scatter. The conclusion should be tested with a joint likelihood that includes both error sets and a free intrinsic-scatter parameter, with a posterior on that parameter, rather than by comparing z-score widths informally.","section":"§4.2 and Figure 8"}],"minor_comments":[{"comment":"The abstract states the sample spans 0.05≤z≤0.3, but Table 3 includes AS1063 at z=0.345; the redshift range should be updated or the sample definition made consistent.","section":"Abstract and Table 3"},{"comment":"The 'chance of observing this correlation due to random chance' (1.25% and 32.09%) is not defined; please specify the null model and whether the p-value comes from a permutation test or from the Monte Carlo realizations under a null hypothesis.","section":"§4.2"},{"comment":"The statement that non-uniform sampling variations up to 30% have no effect on the measured velocity dispersion is cited from A. Rodriguez et al. (2024); since the present work concerns the phase-space edge rather than the velocity dispersion, please clarify whether this test also applies to the edge measurement and Z_v.","section":"§4.3.9"},{"comment":"There are several typos: 'quantity' should be 'quantify' in §3.1, 'Herbonett' should be 'Herbonnet' in §5, and 'affect' should be 'effect' in §4.3.1; the manuscript also references 'Figure A' in §3.1 instead of 'Figure A1'.","section":"§3.1 and §5"},{"comment":"In the sentence discussing the Hubble constant, '±3 km−1s' should be '±3 km s^−1'.","section":"§4.3.7"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and presents an interesting result, but the core claim of zero bias and no intrinsic scatter between escape and weak-lensing masses currently rests on the suppression model Z_v, which has not been tested against hydrodynamical simulations. The authors' energy-based argument is not sufficient to rule out a few-percent bias in the edge population. I encourage a hydrodynamical simulation test or a more direct demonstration of baryonic insensitivity, along with a small additional analysis of the intrinsic-scatter question. The redshift-range inconsistency in the abstract should also be fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a serious paper with a plausible but not yet proven central claim. The new thing is the application of the revised escape-velocity technique to 46 clusters, finding a bias of 0.02±0.02 dex and scatter 0.17 dex against weak-lensing masses. Previous caustic work showed no correlation; here the correlation is 0.679, which is a real improvement. If it holds, it resolves a long-standing tension and opens a new dynamical probe of cosmic expansion.\n\nWhat the paper does well: the simulation tests are the strongest part. The authors show unbiased mass recovery in the Millennium simulation, including with scattered input masses, and they quantify the impact of the shifting-gapper parameters, redshift errors, and binning. The systematics section (4.3) is thorough, and they explicitly test the effect of choosing different cosmologies. The Monte Carlo error propagation for the correlation is appropriate. They also honestly flag the lack of hydrodynamical validation in Section 4.3.9.\n\nThe soft spots, in proportion. The biggest is that the suppression function Z_v, which is the primary systematic of the escape-mass estimate, is calibrated on analytic spherical AGAMA phase-spaces and validated only against gravity-only Millennium halos. A few-percent baryonic bias in Z_v would shift masses by 0.05–0.1 dex, comparable to the quoted bias of 0.02 dex. The stress-test concern lands: the order-of-magnitude energy estimate in 4.3.9 doesn't address how baryonic processes change which galaxies populate the high-velocity edge. This doesn't kill the paper, but it means the concordance is tentative until Z_v is tested with hydro simulations. The circularity from using lensing masses to set the initial r200 is real but I think they've largely addressed it; the 200% mass-error test shows the effect is weak. Calling a 0.679 correlation 'excellent' is an overstatement; it's moderate and improved, but the language should be toned down. The lack of public code/data is a genuine reproducibility gap, though the underlying data are public.\n\nBottom line: this deserves a serious referee. I'd recommend acceptance after major revision, pushing for a hydro test of Z_v or at least a clear statement that the concordance is conditional on the suppression model. The paper is honest about its limitations, and the central result is interesting enough to warrant referee time.","headline":"A plausible resolution of the caustic–WL discrepancy, but the escape-mass suppression model needs hydro validation before I'd treat the concordance as established.","tokens_in":39883,"tokens_out":3189,"would_cite":false,"duration_ms":36408,"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":"Galaxy cluster masses inferred from escape velocities now agree with weak-lensing masses to within 0.02 dex, ending a decade of apparent disagreement.","keywords":["galaxy clusters","escape velocity","weak lensing","phase-space","mass estimation","suppression function","cluster dynamics","cosmology"],"falsifier":"Run the same escape-mass pipeline on mock cluster phase-spaces drawn from a cosmological hydrodynamical simulation with baryonic feedback and realistic survey selection: if the recovered masses differ from the true halo masses by more than about 0.05 dex in the mean, or if re-applying a hydro-calibrated $Z_v$ to the 46 real clusters pushes the lensing bias outside $0.02 \\pm 0.02$ dex, the concordance would be exposed as a calibration artifact rather than a physical agreement.","tokens_in":38662,"feed_emoji":"🌌","tokens_out":5285,"duration_ms":50012,"temperature":0.7,"pith_summary":"The paper argues that the long-standing disagreement between dynamical and weak-lensing cluster masses was a measurement problem, not a physics problem. By re-deriving the escape-velocity edge from the effective potential in an accelerating universe and modeling how sparse galaxy sampling suppresses the observed edge, the authors infer masses for 46 clusters that agree with weak lensing masses to within $0.02 \\pm 0.02$ dex, with a correlation coefficient of $0.679^{+0.046}_{-0.049}$ and a scatter of $0.17$ dex that matches the reported individual uncertainties. If this holds, two independent mass measurement techniques — light deflection and orbital speeds — now give consistent cluster masses, opening a cleaner route to cluster cosmology and tests of gravity on megaparsec scales.","feed_headline":"Escape-velocity masses now match weak lensing for 46 clusters","feed_subtitle":"A skew-normal correction for sparse phase-space sampling erases the old caustic bias and scatter.","key_machinery":"The load-bearing object is the suppression function $Z_v$, the factor by which the true three-dimensional escape velocity profile is pulled inward when the phase-space is sparsely sampled. The paper models the distribution of $Z_v$ as a skew-normal whose location, scale, and skewness depend on the projected tracer count $N$, mass, and redshift, calibrated on analytic spherical phase-spaces and validated against a gravity-only N-body simulation; the observed edge is taken as the maximum absolute line-of-sight velocity in five radial bins, and the cluster mass is recovered by an MCMC comparison of the bin edges to the suppressed theoretical escape profile built from a Dehnen potential in a flat $\\Lambda$CDM cosmology.","core_discovery":"The central discovery is that galaxy cluster masses measured from projected radius–velocity phase-space data, via the down-sampled escape velocity profile corrected by a skew-normal suppression function $Z_v$, are statistically indistinguishable from weak lensing masses. For 46 clusters at $0.05 \\le z \\le 0.3$ spanning $14.4 \\le \\log_{10} M/M_\\odot \\le 15.4$, the mean logarithmic bias is $0.02 \\pm 0.02$ dex, the scatter is $0.17$ dex, and the correlation is $0.679^{+0.046}_{-0.049}$; the scatter is fully consistent with the sum of individual measurement errors, so no additional intrinsic scatter is required. This contrasts with the same comparison made with caustic-inferred masses, which shows a correlation consistent with zero and a $0.25 \\pm 0.05$ dex bias.","pith_inferences":["A hydrodynamical test of $Z_v$ would settle the main residual worry: if baryonic feedback redistributes galaxies in phase-space, the suppression calibration could shift and the reported bias could be hiding a cancellation between errors.","The same framework could be extended to lower-mass groups or higher-redshift clusters, where the strong $N$ dependence of $Z_v$ predicts that sampling sparsity, not internal dynamics, will dominate the error budget.","If the escape technique constrains $qH^2$ as the paper suggests, combining it with independent $H_0$ probes could provide a purely dynamical cross-check of the current Hubble tension.","A stacked or ensemble application of the method, rather than per-cluster mass estimates, might sharpen the cosmological constraint from the same 46 clusters."],"forward_implications":["If the concordance is real, phase-space escape masses can be used as an independent check on weak lensing systematics at the ~0.1 dex level without relying on hydrodynamical simulations.","The skew-normal treatment of the suppression function is essential: ignoring its skewness biases inferred masses low by about 0.1 dex, so future escape-velocity work must include it.","Because the escape edge depends on cosmology through $qH^2$, the same data can in principle constrain the late-universe expansion rate and acceleration, with $H_0$ uncertainty currently the dominant systematic in escape masses.","Concordance holds for clusters with evidence of disturbed dynamical states as well as relaxed ones, implying the method is not limited to equilibrium systems.","The bias between lensing and escape masses shifts with the assumed cosmology, meaning the technique is sensitive enough to distinguish competing $H_0$ values in a larger sample."],"supporting_citations":[{"why":"Introduces the suppression function $Z_v$ and its dependence on phase-space sampling, which this paper re-derives and models as skew-normal.","marker":"V. Halenka et al. (2022)"},{"why":"Supplies the edge-measurement technique, the binning scheme, and the single-cluster validation on Abell S1063 that this paper scales to 46 clusters.","marker":"A. Rodriguez et al. (2024)"},{"why":"Establishes the escape-velocity potential profile for halos in an accelerating universe and tests unbound-particle behavior in simulations.","marker":"C. J. Miller et al. (2016)"},{"why":"Provides the theoretical interpretation of the escape profile for cluster mass inference on which this work builds.","marker":"A. Stark et al. (2016a)"},{"why":"Supplies the weak lensing masses for the majority of the sample and the caustic comparison showing zero correlation.","marker":"R. Herbonnet et al. (2020)"},{"why":"Supplies the independent weak lensing masses used for the second subsample.","marker":"N. Okabe & G. P. Smith (2016)"},{"why":"Provides the Millennium simulation halo and galaxy sample used to validate the suppression model against N-body data.","marker":"D. Gifford et al. (2013)"},{"why":"Provides the potential-density pair used to parametrize the escape velocity profile.","marker":"W. Dehnen (1993)"}],"fun_headline_variants":["Escape-velocity masses now match weak lensing for 46 clusters","Old caustic bias gone: escape masses agree with lensing","Corrected escape-velocity masses indistinguishable from lensing","No extra scatter: escape and lensing masses agree on 46 clusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The suppression function $Z_v$, calibrated on analytic spherical model clusters and validated only against gravity-only N-body halos, is assumed to correct the observed escape edge in real clusters that contain baryonic feedback, substructure, and interlopers without additional bias; a few-percent error in $Z_v$ would shift escape masses enough to make the reported lensing concordance spurious.","fun_headline_variants_meta":{"raw":{"variants":["Escape-velocity masses now match weak lensing for 46 clusters","Old caustic bias gone: escape masses agree with lensing","Corrected escape-velocity masses indistinguishable from lensing","No extra scatter: escape and lensing masses agree on 46 clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000692,"raw_usage":{"total_tokens":3145,"prompt_tokens":974,"completion_tokens":2171,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":2097}},"tokens_in":590,"tokens_out":2171,"duration_ms":19750,"temperature":1.0,"reasoning_tokens":2097,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T13:05:40.472362+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same escape-mass pipeline on mock cluster phase-spaces drawn from a cosmological hydrodynamical simulation with baryonic feedback and realistic survey selection: if the recovered masses differ from the true halo masses by more than about 0.05 dex in the mean, or if re-applying a hydro-calibrated $Z_v$ to the 46 real clusters pushes the lensing bias outside $0.02 \\pm 0.02$ dex, the concordance would be exposed as a calibration artifact rather than a physical agreement.","supporting_citations":[{"cited_title":"J., & Vansickle , P","cited_arxiv_id":null,"evidence_quote":"Introduces the suppression function $Z_v$ and its dependence on phase-space sampling, which this paper re-derives and models as skew-normal."},{"cited_title":"J., Halenka , V., & Kremin , A","cited_arxiv_id":null,"evidence_quote":"Supplies the edge-measurement technique, the binning scheme, and the single-cluster validation on Abell S1063 that this paper scales to 46 clusters."}],"review_version":1}