Pith. sign in

REVIEW 2 major objections 6 minor 24 references

Galaxy velocities on cluster outskirts can pin cluster masses to sub-percent precision with DESI spectra, matching Stage-IV weak lensing.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-14 20:13 UTC pith:DEY7XVT2

load-bearing objection Solid reduced-parameter pure-infall velocity model and transparent DESI Fisher forecast; the sub-percent claim is an ideal-calibration lower bound, which the authors themselves flag. the 2 major comments →

arxiv 2603.22670 v2 pith:DEY7XVT2 submitted 2026-03-24 astro-ph.CO

Cluster Infall for Mass Calibration in the Stage-IV Era

classification astro-ph.CO
keywords galaxy clustersmass calibrationinfall regionsline-of-sight velocitiesDESIweak lensing comparisonpairwise velocitiesStage-IV cosmology
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Cluster masses are hard to measure cleanly because baryonic physics muddies the dense cores. This paper argues that the far outskirts, where galaxies are still falling in, are cleaner: they still feel the cluster's gravity and the expansion of the Universe, but they are largely free of feedback. The authors build a simulation-calibrated model for how radial and tangential speeds of infalling galaxies depend on distance and cluster mass, then project that model along the line of sight so it can be compared with real spectroscopic surveys. With that model in hand they forecast that DESI spectra alone can recover median cluster masses to better than one percent, competitive with or better than stacked weak-lensing forecasts from Stage-IV surveys. The claim matters because a mass calibrator that is both precise and robust to baryons would strengthen the use of cluster abundances for cosmology.

Core claim

A carefully parameterized model of the joint radial-and-tangential velocity distribution of infalling galaxies outside 5 h^{-1} Mpc, once projected along the line of sight, yields P(v_LOS|R,M) accurate enough that DESI-like spectroscopic samples can constrain cluster masses at the sub-percent level, matching or exceeding Stage-IV weak-lensing forecasts.

What carries the argument

The joint velocity distribution P(v_r, v_tan|r,M) written as the product of a Johnson-SU radial distribution and a Student-t conditional tangential distribution, both with mass- and radius-dependent parameters calibrated on simulations, then LOS-projected with the infalling halo-galaxy density profile to produce the observable P(v_LOS|R,M).

Load-bearing premise

The forecast freezes every model parameter at the simulation cosmology and assumes galaxies fall in exactly like dark matter, with no residual velocity bias or imperfect calibration.

What would settle it

Apply the same P(v_LOS|R,M) model to DESI cluster-galaxy pairs in a mass bin already measured by weak lensing or caustic methods; a statistically significant offset between the spectroscopically recovered mass and the independent mass would falsify the claimed sub-percent accuracy.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 6 minor

Summary. The paper constructs a smooth, mass-dependent model for the joint radial–tangential peculiar-velocity distribution of infalling galaxies around clusters on scales r ≥ 5 h⁻¹ Mpc, using MDPL2 + UniverseMachine. The joint PDF is written as P(v_r|r,M) imes P(v_t|v_r,r,M), with a two-parameter Johnson SU for the radial velocities (after exploiting linear relations among the four JSU parameters) and a Student-t (dof = 5) whose scale σ_vt is a cubic in v_r. Both pieces, together with an infalling halo–galaxy correlation function taken from Salazar et al., are projected to obtain P(v_LOS|R,M). The model recovers the simulated joint and projected distributions to roughly 5 %. A Fisher forecast that freezes every model coefficient at the simulation MAP values then claims that DESI BGS/LRG spectra can deliver sub-percent mass precision, competitive with (and in some bins better than) Stage-IV stacked weak-lensing forecasts.

Significance. If the modeling accuracy and the idealized Fisher numbers survive a more realistic error budget, the work supplies a practical, baryon-robust mass-calibration channel that is complementary to weak lensing and that can be applied directly to DESI spectra. The reduction of the joint-velocity model to a compact set of power-law and linear mass trends (Tables I–II), the explicit projection machinery, and the side-by-side comparison with Stage-IV WL forecasts are concrete technical advances over earlier GIK-style analyses. The authors themselves flag the principal caveats (perfect calibration, frozen cosmology, possible galaxy–matter velocity bias), so the paper already points to the next necessary steps.

major comments (2)
  1. Sec. V.A, Eq. (35) and Fig. 7: the Fisher matrix treats only the mass-bin centers M_α as free parameters; every coefficient of P(v_r,v_t|r,M) (Table I) and of the infalling density profile (Table II) is held fixed at the MDPL2 MAP values, and the cosmology is likewise frozen. The abstract and Fig. 7 therefore present optimistic lower bounds rather than realistic DESI uncertainties. Sec. VI already notes this idealization; the forecast section should either (i) marginalize a representative subset of the dominant nuisance parameters (or add a systematic floor) or (ii) re-label the quoted numbers as ideal-calibration limits and show how they degrade under plausible residual systematics.
  2. Sec. VI and the comparison in Fig. 7: the paper acknowledges possible galaxy-versus-matter velocity bias but does not quantify its impact on the recovered mass. Because the entire calibration rests on the assumption that the simulated galaxy velocities faithfully trace the mass, even a few-percent coherent bias would shift the mass scale at a level comparable to the claimed statistical precision. A short test (e.g., rescaling the mean infall velocity or σ_vr by a few percent and re-running the Fisher) would make the robustness claim quantitative rather than qualitative.
minor comments (6)
  1. Abstract and Sec. I: the phrase “can be used to for cluster mass calibration” contains a duplicated preposition; correct to “can be used for”.
  2. Sec. II: the stellar-mass cut is written M_* = 10^10 h^{-1} M_⊙; confirm whether the h-scaling is intentional or should be M_* ≥ 10^{10} M_⊙/h.
  3. Fig. 3 caption: the text mentions 68 %, 95 %, and 97 % contours while the body text refers to 99.7 %; align the numbers.
  4. Eq. (3) and surrounding text: the ± convention for the Hubble term is clear, but a one-sentence reminder that the sign is chosen so that the Hubble flow always points away from the cluster would help readers less familiar with the distant-observer setup.
  5. Table I: several parameters (e.g., A, C_1,c) carry units that are easy to misread; a short column header or footnote listing units for every entry would improve usability.
  6. Sec. IV: the surface-density integral (Eq. 32) assumes a Gaussian LOS velocity distribution with a fixed σ_LOS = 532 km s^{-1}; a brief statement of how sensitive Σ_inf is to that choice would be useful.

Circularity Check

0 steps flagged

No load-bearing circularity: simulation-calibrated velocity/density model is projected and Fisher-forecasted under frozen parameters; the DESI sub-percent claim is an idealized forecast, not a tautology of the fit.

full rationale

The derivation chain is: (i) fit a reduced-parameter JSU + Student-t model for P(v_r, v_t | r, M) and the Salazar et al. infalling density profile to MDPL2+UniverseMachine (Tables I–II, Figs. 1–4); (ii) project via the standard LOS integral (Eqs. 29–33) to obtain P(v_LOS | R, M); (iii) form a Poisson Fisher matrix (Eq. 35) whose only free parameters are the mass-bin centers M_α, with all velocity and density coefficients held fixed at the simulation MAP values and N_pairs rescaled to DESI number densities. The projected distributions are then compared to the same simulation (Figs. 5–6) as a consistency check, not as an independent prediction. Prior works (Zu & Weinberg, Aung et al., Salazar et al.) supply functional forms that are re-calibrated here; none is invoked as a uniqueness theorem that forces the mass-precision claim. The authors themselves flag that the forecast “assumes perfect model calibration and neglects cosmology dependence” (Sec. VI). That is an optimistic assumption, not a circular reduction of the result to its inputs. Score 1 only for the minor, non-load-bearing reliance on the same research lineage for the starting ansatz.

Axiom & Free-Parameter Ledger

9 free parameters · 5 axioms · 0 invented entities

The central forecast rests on a large set of free parameters fitted to one simulation snapshot, on the assumption that UniverseMachine galaxies trace the relevant velocity field, on a fixed Planck cosmology, and on the idealizations of the Fisher matrix (perfect model, no cosmology dependence, no velocity bias). No new physical entities are postulated; the work is phenomenological calibration plus forecasting.

free parameters (9)
  • JSU shape intercepts and slopes (γ-bar, δ-bar linear relations)
    Two linear relations that collapse four JSU parameters to two; coefficients fitted to MDPL2 velocity histograms.
  • v_r,peak and σ²_vr power-law amplitudes and slopes (v_p,p, v_p,s, v_s,p, v_s,s, σ_p,p, σ_p,s, σ_s,p, σ_s,s)
    Eight parameters controlling radial dependence of the peak and variance of P(v_r); MAP values in Table I.
  • Δm, Δb (mass-dependent δ correction)
    Extra linear mass term added post-hoc to improve the smooth-model shape at high mass.
  • σ²_r,c (large-scale radial velocity variance floor)
    Mass-independent constant fixed at 2.13e5 (km/s)².
  • A, B_p, B_s, μ0,c/p, μ1,c/p, C1,c/p (cubic σ_vt model)
    Parameters of the cubic conditional tangential-velocity width and its mass/radius dependence; Table I.
  • Student-t degrees of freedom = 5
    Fixed by hand, independent of mass and radius, following earlier literature.
  • Infalling density-profile parameters (r_h,p/s, b_p/s, γ_p/s, η0, ηm, ησ, r_inf, μ, Δ)
    Eleven parameters of the Salazar et al. infalling ξ model, re-fitted to the same simulation (Table II).
  • σ_LOS = 532 km/s (for surface-density noise integral)
    Measured once from the simulation and held fixed in the projection.
  • DESI number densities n_g,BGS and n_g,LRG
    Taken from an external figure and held constant across broad redshift ranges for the Fisher forecast.
axioms (5)
  • domain assumption UniverseMachine galaxies with M* > 10^10 h^-1 M⊙ faithfully sample the dark-matter velocity field on r ≥ 5 h^-1 Mpc for the purpose of mass calibration.
    All velocity histograms and density profiles are measured on UM galaxies; residual galaxy–matter velocity bias is acknowledged only as a future systematic.
  • domain assumption A single MDPL2 snapshot at a=0.8376 (z≈0.194) is representative for DESI BGS/LRG forecasts across 0.1 < z < 0.8.
    Model calibrated at one redshift; no redshift evolution of the velocity PDFs is modeled.
  • ad hoc to paper Fisher matrix with Poisson covariance and all nuisance parameters fixed yields a realistic mass uncertainty.
    Sec. V.A freezes the model at the simulation cosmology; Sec. VI notes this neglects cosmology dependence and model uncertainty.
  • domain assumption Orbiting galaxies contribute negligibly for r ≥ 5 h^-1 Mpc (or R ≥ 4–5 h^-1 Mpc).
    Justified by the Salazar et al. decomposition and Fig. 4; sets the analysis floor.
  • standard math Distant-observer approximation and periodic-box geometry introduce negligible bias for the projected LOS distributions.
    Standard simulation practice; not re-validated here.

pith-pipeline@v1.1.0-grok45 · 19336 in / 3787 out tokens · 28720 ms · 2026-07-14T20:13:29.934887+00:00 · methodology

0 comments
read the original abstract

The outskirts of galaxy clusters present a promising avenue for constraining cluster masses in a way that is robust to the impact of baryonic physics. We assess the accuracy to which the cluster infall regions can be used for cluster mass calibration. Building on previous work, we parameterize the velocity distribution $P(v_{\rm r},v_{\rm tan}|r,M)$ of dark matter halos on scales $r \geq 5\ h^{-1}\ \rm{Mpc}$ as the product of the marginalized distribution $P(v_{\rm r}|r,M)$ and the conditional distribution $P(v_{\rm tan}|v_{\rm r},r,M)$, calibrating the radial and mass dependence of these distributions in numerical simulations. We then project our model along the line-of-sight to obtain accurate predictions for the distributions of line-of-sight velocities at a given projected radius and cluster mass $P(v_{\rm LOS}|R,M)$, which we can observe with spectroscopic survey data. With our model, we forecast that spectra from the Dark Energy Spectroscopic Instrument (DESI) can constrain cluster masses with sub-percent level precision, comparable to that of Stage IV weak lensing surveys.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

24 extracted references · 13 linked inside Pith

  1. [1]

    di- mensionless

    [and therefore 3]. We then use the resulting velocity probability distributions to determine the accuracy with which pairwise-velocity data in the cluster-infall regions arXiv:2603.22670v1 [astro-ph.CO] 24 Mar 2026 2 can constrain clusters masses with DESI, updating the forecasts from [5]. We structure the remainder of this paper as follows: in Section II...

  2. [2]

    A. N. Salcedo, E. Rozo, H.-Y. Wu, D. H. Wein- berg, P. Chiploonkar, C.-H. To, S. Cao, E. S. Rykoff, N. Marcelina Gountanis, and C. Zhou, Cosmological Con- straints from Dark Energy Survey Year 1 Cluster Lens- 11 ing and Abundances with Simulation-based Forward- Modeling, arXiv e-prints , arXiv:2510.25706 (2025), arXiv:2510.25706 [astro-ph.CO]

  3. [3]

    Miyatake, Cosmology with galaxy clusters (2025), arXiv:2505.07697 [astro-ph.CO]

    H. Miyatake, Cosmology with galaxy clusters (2025), arXiv:2505.07697 [astro-ph.CO]

  4. [4]

    Zu and D

    Y. Zu and D. H. Weinberg, The redshift-space cluster- galaxy cross-correlation function - I. Modelling galaxy infall on to Millennium simulation clusters and SDSS groups, MNRAS431, 3319 (2013), arXiv:1211.1379 [astro-ph.CO]

  5. [5]

    Cuesta-Lazaro, B

    C. Cuesta-Lazaro, B. Li, A. Eggemeier, P. Zarrouk, C. M. Baugh, T. Nishimichi, and M. Takada, To- wards a non-gaussian model of redshift space distortions, Monthly Notices of the Royal Astronomical Society498, 1175 (2020), https://academic.oup.com/mnras/article- pdf/498/1/1175/33731014/staa2249.pdf

  6. [6]

    Hamabata, M

    A. Hamabata, M. Oguri, and T. Nishimichi, Con- straining cluster masses from the stacked phase space distribution at large radii, Monthly No- tices of the Royal Astronomical Society489, 1344 (2019), https://academic.oup.com/mnras/article- pdf/489/1/1344/29229772/stz2227.pdf

  7. [7]

    H. Aung, D. Nagai, E. Rozo, B. Wolfe, and S. Ad- hikari, Accurate model of the projected velocity dis- tribution of galaxies in dark matter haloes, Monthly Notices of the Royal Astronomical Society521, 3981 (2023), https://academic.oup.com/mnras/article- pdf/521/3/3981/49672134/stad601.pdf

  8. [8]

    T. Y. Lam, T. Nishimichi, F. Schmidt, and M. Takada, Testing Gravity with the Stacked Phase Space around Galaxy Clusters, Phys. Rev. Lett.109, 051301 (2012), arXiv:1202.4501 [astro-ph.CO]

  9. [9]

    Y. Zu, D. H. Weinberg, E. Jennings, B. Li, and M. Wyman, Galaxy infall kinematics as a test of modi- fied gravity, MNRAS445, 1885 (2014), arXiv:1310.6768 [astro-ph.CO]

  10. [10]

    Tomooka, E

    P. Tomooka, E. Rozo, E. L. Wagoner, H. Aung, D. Na- gai, and S. Safonova, Clusters have edges: the projected phase-space structure of SDSS redMaPPer clusters, MN- RAS499, 1291 (2020), arXiv:2003.11555 [astro-ph.CO]

  11. [11]

    Robertson, E

    A. Robertson, E. Huff, K. Markoviˇ c, and B. Li, Modelling the redshift-space cluster-galaxy correlation function on Mpc scales with emulation of the pairwise velocity dis- tribution, MNRAS533, 4081 (2024), arXiv:2406.01527 [astro-ph.CO]

  12. [12]

    Klypin, G

    A. Klypin, G. Yepes, S. Gottl¨ ober, F. Prada, and S. Heß, MultiDark simulations: the story of dark matter halo concentrations and density profiles, MNRAS457, 4340 (2016), arXiv:1411.4001 [astro-ph.CO]

  13. [13]

    E. S. Rykoff, E. Rozo, M. T. Busha, C. E. Cunha, A. Finoguenov, A. Evrard, J. Hao, B. P. Koester, A. Leauthaud, B. Nord, M. Pierre, R. Reddick, T. Sadibekova, E. S. Sheldon, and R. H. Wechsler, redMaPPer. I. Algorithm and SDSS DR8 Catalog, ApJ 785, 104 (2014), arXiv:1303.3562 [astro-ph.CO]

  14. [14]

    E. S. Rykoff, E. Rozo, D. Hollowood, A. Bermeo- Hernandez, T. Jeltema, J. Mayers, A. K. Romer, P. Rooney, A. Saro, C. Vergara Cervantes, R. H. Wechsler, H. Wilcox, T. M. C. Abbott, F. B. Ab- dalla, S. Allam, J. Annis, A. Benoit-L´ evy, G. M. Bern- stein, E. Bertin, D. Brooks, D. L. Burke, D. Capozzi, A. Carnero Rosell, M. Carrasco Kind, F. J. Castander, M...

  15. [15]

    P. S. Behroozi, R. H. Wechsler, and H.-Y. Wu, The ROCKSTAR Phase-space Temporal Halo Finder and the Velocity Offsets of Cluster Cores, ApJ762, 109 (2013), arXiv:1110.4372 [astro-ph.CO]

  16. [16]

    P. S. Behroozi, R. H. Wechsler, H.-Y. Wu, M. T. Busha, A. A. Klypin, and J. R. Primack, Gravitationally Con- sistent Halo Catalogs and Merger Trees for Precision Cosmology, ApJ763, 18 (2013), arXiv:1110.4370 [astro- ph.CO]

  17. [17]

    Behroozi, R

    P. Behroozi, R. H. Wechsler, A. P. Hearin, and C. Conroy, UNIVERSEMACHINE: The correlation between galaxy growth and dark matter halo assembly from z = 0- 10, MNRAS488, 3143 (2019), arXiv:1806.07893 [astro- ph.GA]

  18. [18]

    G. L. Bryan and M. L. Norman, Simulating X-ray Clus- ters with Adaptive Mesh Refinement, inAmerican As- tronomical Society Meeting Abstracts, American Astro- nomical Society Meeting Abstracts, Vol. 187 (1995) p. 95.04

  19. [19]

    H. Aung, D. Nagai, E. Rozo, and R. Garc´ ıa, The phase-space structure of dark matter haloes, Monthly Notices of the Royal Astronomical Society502, 1041 (2020), https://academic.oup.com/mnras/article- pdf/502/1/1041/36169658/staa3994.pdf

  20. [20]

    Foreman-Mackey, A

    D. Foreman-Mackey, A. Conley, W. Meierjurgen Farr, D. W. Hogg, D. Lang, P. Marshall, A. Price-Whelan, J. Sanders, and J. Zuntz, emcee: The MCMC Hammer, Astrophysics Source Code Library, record ascl:1303.002 (2013), ascl:1303.002

  21. [21]

    Scoccimarro, Redshift-space distortions, pairwise ve- locities, and nonlinearities, Phys

    R. Scoccimarro, Redshift-space distortions, pairwise ve- locities, and nonlinearities, Phys. Rev. D70, 083007 (2004)

  22. [22]

    E. M. Salazar, E. Rozo, R. Garc´ ıa, N. Kokron, S. Ad- hikari, B. Diemer, and C. Osinga, Dynamics-based halo model for large scale structure, Phys. Rev. D111, 043527 (2025)

  23. [23]

    D. H. Weinberg, M. J. Mortonson, D. J. Eisenstein, C. Hi- rata, A. G. Riess, and E. Rozo, Observational probes of cosmic acceleration, Physics Reports530, 87 (2013), ob- servational Probes of Cosmic Acceleration

  24. [24]

    Abdul Karim, J

    M. Abdul Karim, J. Aguilar, S. Ahlen, S. Alam, L. Allen, C. A. Prieto, O. Alves, A. Anand, U. Andrade, E. Ar- mengaud, A. Aviles, S. Bailey, C. Baltay, P. Bansal, A. Bault, J. Behera, S. BenZvi, D. Bianchi, C. Blake, S. Brieden, A. Brodzeller, D. Brooks, E. Buckley-Geer, E. Burtin, R. Calderon, R. Canning, A. C. Rosell, P. Carrilho, L. Casas, F. J. Castan...