Recognition: unknown
Testing template-fitting models for the multipoles of the two-point clustering of galaxy clusters
Pith reviewed 2026-05-07 14:57 UTC · model grok-4.3
The pith
The dispersion model for galaxy cluster clustering multipoles yields unbiased fσ8 inferences down to 10 h^{-1} Mpc in Euclid DR1-like conditions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors apply the dispersion, Scoccimarro, and Taruya-Nishimichi-Saito models to the multipoles of the redshift-space two-point correlation function extracted from PINOCCHIO simulations of clusters with virial mass above 10^{14} h^{-1} M_odot over 0 < z < 2 in a DR1-like footprint. After constructing covariance matrices from 1000 realizations and imposing permissive and conservative scale cuts, they infer fσ8 and find that the dispersion model returns unbiased values down to 10 h^{-1} Mpc under both realistic and optimistic photometric redshift uncertainty assumptions, while all three models produce statistically indistinguishable goodness-of-fit metrics.
What carries the argument
Template-fitting models of redshift-space distortions applied to the even multipoles of the two-point correlation function of galaxy clusters, including modeling of nonlinear evolution, halo bias, and photometric redshift errors.
If this is right
- The dispersion model can be used for unbiased fσ8 inference from Euclid cluster clustering multipoles down to 10 h^{-1} Mpc.
- All three models deliver comparable goodness-of-fit metrics once DR1-like photometric redshift uncertainties are included.
- The validated multipole cut-off criteria allow template fitting to be applied reliably in DR1-like settings.
- Template fitting of cluster multipoles provides a viable complement to full-shape methods for cosmological inference.
Where Pith is reading between the lines
- If the dispersion model performs equally well on real data, analyses for early Euclid releases could rely on the simplest implementation without loss of accuracy.
- The result suggests that model choice is unlikely to dominate systematics in DR1 cluster cosmology, freeing effort for other observational uncertainties.
- The framework could be tested on lower-mass thresholds or combined with galaxy clustering to tighten constraints on modified gravity models.
Load-bearing premise
The PINOCCHIO simulations with third-order Lagrangian perturbation theory accurately represent the clustering properties and photometric redshift uncertainties of real optically-selected galaxy clusters in the Euclid footprint.
What would settle it
Recovering a statistically significant bias in fσ8 when the dispersion model is fitted to the measured multipoles from the actual Euclid DR1 cluster catalog at scales below 10 h^{-1} Mpc.
Figures
read the original abstract
The \textit{Euclid} satellite will deliver a catalogue of optically-selected galaxy clusters spanning from around $2000$ deg$^2$ in Data Release (DR) 1 to around $14\,000$ deg$^2$ in DR3. We assess the validity of cluster clustering (CC) models for template-fitting, which complements the full-shape methodology in providing cosmological information from the anisotropy of the redshift-space two-point correlation function (2PCF). Both will be used to analyse the cluster 2PCF multipoles in \textit{Euclid}. We analyse the multipoles of the two-point redshift-space clustering of galaxy clusters simulated with the semi-analytic \code{PINOCCHIO} code using third-order Lagrangian perturbation theory, assuming a \textit{Euclid} DR1-like footprint of 500 deg$^2$ in the Northern Hemisphere and 1400 deg$^2$ in the Southern Hemisphere. We estimate the first three even multipoles of the 2PCF and associated covariance matrix from 1000 DR1-like synthetic catalogues. We study the impact of the modelling of nonlinearities, halo bias, and photometric redshift uncertainties on the 2PCF. We apply three clustering models to the mock catalogues at $0<z<2$ and virial mass $M_{\rm vir}>10^{14}\;h^{-1}\,M_\odot$ under realistic and optimistic photometric redshift uncertainty scenarios. We formulate a set of permissive and conservative criteria that should be fulfilled by the multipole cut-off scales and validate against 100 mock catalogues via inference of the growth rate times the matter power spectrum normalisation parameter, $f\sigma_8$. We test the dispersion, Scoccimarro, and Taruya--Nishimichi--Saito models. We find that the simplest of the three -- the dispersion model -- yields unbiased inferences on $f\sigma_8$ from CC down to $10$ $h^{-1}$ Mpc in a DR1-like setting. All clustering models provide very similar goodness-of-fit metrics in the presence of DR1-like cluster redshift uncertainties.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper tests three template-fitting models (dispersion, Scoccimarro, and Taruya-Nishimichi-Saito) for the even multipoles of the redshift-space two-point correlation function of galaxy clusters using 1000 PINOCCHIO mocks with 3LPT in a Euclid DR1-like footprint. It estimates multipoles and covariances, incorporates photometric redshift uncertainties, defines permissive and conservative cut-off criteria validated via fσ8 recovery tests, and concludes that the dispersion model recovers unbiased fσ8 down to 10 h^{-1} Mpc while all models yield similar goodness-of-fit under DR1-like photo-z errors.
Significance. If the PINOCCHIO mocks adequately capture the relevant clustering and error properties, the result supports adopting the simple dispersion model for template fitting of cluster multipoles in Euclid analyses, providing a computationally efficient complement to full-shape methods for growth-rate constraints. The large mock ensemble and explicit recovery tests on fσ8 are strengths that make the validation concrete and falsifiable within the simulated setting.
major comments (2)
- [Mock generation and simulation setup (implicit in § on PINOCCHIO catalogues)] The central claim that the dispersion model yields unbiased fσ8 inferences down to 10 h^{-1} Mpc rests entirely on the fidelity of the PINOCCHIO 3LPT mocks to real optically-selected cluster clustering (including higher-order bias, velocity dispersion, and the detailed form of photo-z errors). No comparison to full N-body or hydrodynamical simulations is presented to quantify residual systematics in the multipoles at these scales; this is load-bearing for extrapolating the result beyond the specific mock suite.
- [Validation criteria and fσ8 recovery tests] The permissive and conservative criteria for multipole cut-off scales are validated only on the same 1000-mock ensemble used to build the covariance; it is unclear whether the covariance estimation (from 1000 realisations) remains stable and invertible at the smallest scales tested, or whether the criteria would hold under a different mock realisation or covariance regularisation scheme.
minor comments (2)
- [Abstract] The abstract and introduction should explicitly state the mass and redshift cuts (M_vir > 10^{14} h^{-1} M_⊙, 0 < z < 2) when summarising the main result on the dispersion model.
- [Clustering models section] Notation for the three models is introduced without a compact table comparing their functional forms for the multipoles; adding such a table would improve readability.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive feedback on our manuscript. We address each major comment in detail below, proposing revisions where appropriate to enhance the clarity and robustness of our results.
read point-by-point responses
-
Referee: The central claim that the dispersion model yields unbiased fσ8 inferences down to 10 h^{-1} Mpc rests entirely on the fidelity of the PINOCCHIO 3LPT mocks to real optically-selected cluster clustering (including higher-order bias, velocity dispersion, and the detailed form of photo-z errors). No comparison to full N-body or hydrodynamical simulations is presented to quantify residual systematics in the multipoles at these scales; this is load-bearing for extrapolating the result beyond the specific mock suite.
Authors: We agree that the fidelity of the simulation suite is crucial for the generalizability of our findings. Our study is designed to test the performance of template-fitting models in a Euclid DR1-like setting using mocks that incorporate the relevant effects, including 3LPT for nonlinear evolution and photometric redshift uncertainties. While we do not perform a direct comparison to full N-body simulations in this work, we note that PINOCCHIO has been extensively validated against N-body results for halo clustering in the literature. In the revised manuscript, we will expand the discussion section to explicitly address potential limitations of the 3LPT mocks, such as missing higher-order effects or baryonic physics, and caution readers about extrapolating beyond the mock setup. This constitutes a partial revision as we cannot undertake new large-scale simulations at this time. revision: partial
-
Referee: The permissive and conservative criteria for multipole cut-off scales are validated only on the same 1000-mock ensemble used to build the covariance; it is unclear whether the covariance estimation (from 1000 realisations) remains stable and invertible at the smallest scales tested, or whether the criteria would hold under a different mock realisation or covariance regularisation scheme.
Authors: The covariance matrix is estimated from 1000 independent mock catalogues, which is a standard and sufficient number for the dimensionality of our data vector (three even multipoles). We have confirmed the numerical stability by ensuring the matrix is positive definite and invertible, with condition numbers that do not indicate instability at the smallest scales. The cut-off criteria were validated using an independent set of 100 mocks. To further alleviate concerns, in the revision we will add a subsection discussing the robustness of the covariance estimation, including checks with varying numbers of mocks and alternative regularisation methods if applicable. This will show that our conclusions are not sensitive to these choices. revision: partial
- Direct quantitative comparison of PINOCCHIO 3LPT multipoles to those from full N-body or hydrodynamical simulations at scales down to 10 h^{-1} Mpc.
Circularity Check
No circularity: forward simulation validates models against independent mocks
full rationale
The paper generates 1000 DR1-like mock catalogs using the independent PINOCCHIO code with 3LPT, computes the first three even multipoles of the redshift-space 2PCF and covariance, then fits three template models (dispersion, Scoccimarro, TNS) and checks recovery of input fσ8 via inference on scales down to 10 h^{-1} Mpc. No equation or step reduces the reported unbiased recovery to a fitted parameter by construction, nor does any load-bearing premise rest on self-citation. The central claim is an empirical validation result conditional on the mocks, which is self-contained and externally falsifiable.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Third-order Lagrangian perturbation theory in PINOCCHIO produces sufficiently accurate mock catalogues for testing cluster clustering models at the scales and redshifts considered
Reference graph
Works this paper leans on
-
[1]
Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2020, Phys. Rev. D, 102, 023509 Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2025a, Phys. Rev. D, 112, 083535 Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2025b, Phys. Rev. D, 112, 083535 Abell, P. A., Allison, J., Anderson, S. F., et al. 2009, arXiv:0912.0201 Adame, A. G., Aguilar, J., Ahlen, S....
-
[2]
Contribution of the hexadecapole10 For the investigation of the contribution ofℓ=4 we compare thefσ 8 posteriors with and without including it in the inference
A.1. Contribution of the hexadecapole10 For the investigation of the contribution ofℓ=4 we compare thefσ 8 posteriors with and without including it in the inference. We work in the highest-S/N redshift bin as above, and consider rmin,4 =r min,0 =10h −1Mpc, following the assumptions of Eu- clid Collaboration: Kärcher et al. (2026). In Fig. A.1 we present o...
2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.