# CHEX-MATE: The Impact of Triaxiality and Orientation on Planck SZ Cluster Selection and Weak Lensing Mass Measurements

H Saxena, J Sayers, A Gavidia, J.-B Melin, E.T Lau, J Kim, L Chappuis, D Eckert, S Ettori, M Gaspari, et al.

► To cite this version:

H Saxena, J Sayers, A Gavidia, J.-B Melin, E.T Lau, et al., CHEX-MATE: The Impact of Triaxiality and Orientation on Planck SZ Cluster Selection and Weak Lensing Mass Measurements. Astronomy & Astrophysics - A&A, 2025, 700, pp.A128. 10.1051/0004-6361/202555719. hal-05108916

HAL Id: hal-05108916

https://hal.science/hal-05108916v1

Submitted on 16 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# CHEX-MATE: The impact of triaxiality and orientation on Planck SZ cluster selection and weak lensing mass measurements

H. Saxena<fnref n="1" />,*, J. Sayers<fnref n="1" />, A. Gavidia<fnref n="1" />, J.-B. Melin<fnref n="2" />, E. T. Lau<fnref n="3" />, J. Kim<fnref n="4" />, L. Chappuis<fnref n="2" />,<fnref n="17" />, D. Eckert<fnref n="17" />, S. Ettori<fnref n="6" />,<fnref n="7" />, M. Gaspari<fnref n="12" />, F. Gastaldello<fnref n="8" />, C. Giocoli<fnref n="6" />,<fnref n="7" />, S. Kay<fnref n="18" />, L. Lovisari<fnref n="8" />,<fnref n="9" />, B. Maughan<fnref n="19" />, F. Oppizzi<fnref n="11" />, M. De Petris<fnref n="13" />, G. W. Pratt<fnref n="5" />, E. Pointecouteau<fnref n="10" />, E. Rasia<fnref n="14" />,<fnref n="15" />,<fnref n="16" />, M. Rossetti<fnref n="8" />, and M. Sereno<fnref n="6" />,<fnref n="7" />

Received 29 May 2025 / Accepted 10 June 2025

## ABSTRACT

**Context.** Galaxy cluster abundance measurements are a valuable tool for constraining cosmological parameters, such as the mass density ($\Omega_m$) and the density fluctuation amplitude ($\sigma_8$). Wide-area surveys detect clusters based on observables, such as the total integrated Sunyaev-Zel'dovich effect signal ($Y_{SZ}$) in the case of *Planck*. Quantifying the survey selection function is necessary for cosmological analyses, with completeness representing the probability of detecting a cluster as a function of its intrinsic properties, such as $Y_{SZ}$ and an angular scale $\theta_{500}$.

**Aims.** We determine the completeness of the *Planck*-selected CHEX-MATE cluster catalog using mock observations of clusters with triaxial shapes and random orientations, with physically-motivated distributions of axial ratios. From these mocks, we derive the distribution of shapes and orientations of the detected clusters, along with any associated bias in weak-lensing-derived mass ($M_{WL}$) due to this orientation-dependent selection (denoted as $1 - b_\chi$).

**Methods.** Employing a Monte Carlo method, we injected triaxial cluster profiles into random positions within the *Planck* all-sky maps and subsequently determined the completeness as a function of both geometry and SZ brightness. This result was then used to generate 1000 mock CHEX-MATE cluster catalogs. We computed $M_{WL}$ for these mock CHEX-MATE clusters and for equal-sized samples of randomly selected clusters with similar mass and redshift distributions.

**Results.** Cluster orientation impacts completeness, with a higher probability of detecting clusters elongated along the line of sight (LOS). This leads to $1 - b_\chi$ values of 0–4% for CHEX-MATE clusters relative to a random population. The largest increase in $M_{WL}$ is observed in the lowest-mass objects, which are most impacted by orientation-related selection bias.

**Conclusions.** Clusters in *Planck* SZ-selected catalogs are preferentially elongated along the LOS and have an average bias in $M_{WL}$ relative to randomly selected cluster samples. This bias is relevant for upcoming SZ surveys such as CMB-S4, and should be considered for surveys utilizing other probes for cluster detection, such as *Euclid*.

**Key words.** gravitational lensing: weak – galaxies: clusters: general – galaxies: clusters: intracluster medium – cosmology: observations

## 1. Introduction

Galaxy cluster abundance is a powerful cosmological probe, offering incisive tests of the ΛCDM model and precise measurements of fundamental cosmological parameters (e.g., Voit 2005; Allen et al. 2011; Lesci et al. 2022a). As the most massive gravitationally bound structures in the Universe, galaxy clusters trace the late stages of structure formation, representing the evolution of the highest peaks in the primordial matter density field

<footnote n="1">¹ California Institute of Technology, 1200 East California Boulevard, Pasadena, California, USA</footnote>

<footnote n="2">² IRFU, CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="3">³ Smithsonian Astrophysical Observatory, Cambridge, Massachusetts, USA</footnote>

<footnote n="4">⁴ Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea</footnote>

<footnote n="5">⁵ Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="6">⁶ INAF – Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, Via Piero Gobetti 93/3, I-40129 Bologna, Italy</footnote>

<footnote n="7">⁷ INFN, Sezione di Bologna, Viale Berti Pichat 6/2, I-40127 Bologna, Italy</footnote>

<footnote n="8">⁸ INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica di Milano, Via A. Corti 12, 20133 Milano, Italy</footnote>

<footnote n="9">⁹ Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA</footnote>

<footnote n="10">¹⁰ IRAP, CNRS, Université de Toulouse, CNES, Toulouse, France</footnote>

<footnote n="11">¹¹ INFN-Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy</footnote>

<footnote n="12">¹² Department of Physics, Informatics and Mathematics, University of Modena and Reggio Emilia, 41125 Modena, Italy</footnote>

<footnote n="13">¹³ Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Rome, Italy</footnote>

<footnote n="14">¹⁴ INAF, Osservatorio Astronomico di Trieste, Via Tiepolo 11, I-34131 Trieste, Italy</footnote>

<footnote n="15">¹⁵ IFPU, Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy</footnote>

<footnote n="16">¹⁶ Department of Physics; University of Michigan, Ann Arbor, MI 48109, USA</footnote>

<footnote n="17">¹⁷ Department of Astronomy, University of Geneva, Ch. d'Ecogia 16, CH-1290 Versoix, Switzerland</footnote>

<footnote n="18">¹⁸ Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, UK</footnote>

<footnote n="19">¹⁹ HH Wills Physics Laboratory, University of Bristol, Bristol, UK</footnote>

<footnote>* Corresponding author.</footnote>
(Despali et al. 2016). Their abundance is particularly useful for constraining both the mean matter density ($\Omega_m$) and the amplitude of matter fluctuations ($\sigma_8$), defined as the root mean square (rms) of linear density perturbations on scales of $8 \hbar^{-1}$ Mpc. Recent studies using cluster catalogs selected from X-ray, optical, and millimeter-wavelength surveys have yielded increasingly precise constraints on these parameters (Ghirardini et al. 2024; DES Collaboration 2025; Bocquet et al. 2024), although challenges related to halo mass calibration and selection biases still remain.

The hot intracluster medium (ICM) emits X-rays via thermal bremsstrahlung (Sarazin 1986), producing a bright, diffuse signal that makes X-ray surveys a valuable tool for detecting galaxy clusters. Early X-ray halo abundance studies were based on catalogs of approximately 100 clusters detected in the ROSAT All-Sky Survey (RASS), using deep X-ray follow-up from *Chandra* to determine halo masses via hydrostatic equilibrium-based estimates (Vikhlinin et al. 2009; Mantz et al. 2010). Because the bias associated with such mass estimates is difficult to quantify, later work improved on these results by utilizing weak-lensing (WL) observations of a subsample of clusters to calibrate the halo mass scale (Mantz et al. 2015). The more recent XXL survey from *XMM-Newton*, which is much narrower and deeper than the RASS, delivered a comparably sized cluster sample with WL again calibrating the mass scale for a halo abundance measurement (Pacaud et al. 2018). The current state of the art in X-ray cluster cosmology was obtained from the eROSITA (extended ROentgen Survey with an Imaging Telescope Array) All-Sky Survey (eRASS), based on over 5000 X-ray-selected clusters (Ghirardini et al. 2024). Cluster masses were calibrated using WL, in this case using WL mass ($M_{WL}$) estimates from three different optical surveys (Chiu et al. 2022; Kleinebreil et al. 2025; Grandis et al. 2024), enabling a robust scaling between X-ray observables and total halo mass (see also Okabe et al. 2025).

In the optical regime, cluster detection is most often based on searching for overdensities of galaxies, for example, using algorithms that attempt to identify a red sequence of cluster-member galaxies (e.g., Rykoff et al. 2014; Bellagamba et al. 2017). For instance, commensurate with the RASS-based studies noted above, the maxBCG cluster sample from the Sloan Digital Sky Survey (SDSS) was used by Rozo et al. (2010) to derive cosmological constraints from over $10^4$ clusters and groups at low redshifts. To connect the detection observable, optical richness, with the underlying halo mass, stacked $M_{WL}$ values within bins of fixed richness were used. Subsequent work by Costanzi et al. (2019) combines abundance and WL measurements from the SDSS DR8 to jointly constrain cosmology and the richness–mass relation for a photometrically selected sample of red-sequence Matched-filter Probabilistic Percolation (redMaPPer) clusters. Illustrating the challenges related to selection and mass calibration, the Dark Energy Survey (DES) released its initial halo abundance measurement soon after (Abbott et al. 2020), with an anomalously low value of the matter density that was later found to be due to biases in $M_{WL}$ values for the lowest-richness clusters (To et al. 2021; Costanzi et al. 2021). The current state of the art in optical cluster cosmology is presented in the DES Y3 results and the third release of the Kilo-Degree Survey (DES Collaboration 2025; Lesci et al. 2022b), who find cosmological parameter values that are fully consistent with those obtained from primary CMB anisotropy measurements from Planck.

Scattering of CMB photons with electrons in the ICM produces an average boost in photon energy known as the

Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zeldovich 1972), which enables the detection of clusters in millimeter-wave surveys. Recent large-area surveys from the Atacama Cosmology Telescope (Hilton et al. 2021), the South Pole Telescope (Bleem et al. 2015, 2024; Kornoelje et al. 2025), and Planck (Planck Collaboration XXVII 2016, hereafter P16) have now provided catalogs with thousands of clusters, which have in turn been used to constrain cosmological parameters (Sehgal et al. 2011; Bocquet et al. 2024; Planck Collaboration XXIV 2016). Notably, the SZ effect surface brightness is independent of redshift, making these surveys approximately mass-limited across most redshifts, and thus providing catalogs that are particularly attractive for halo abundance measurements. However, results from these surveys are generally limited by systematic uncertainties in halo mass calibration. For example, Planck detected fewer clusters than predicted by the base $\Lambda$CDM model, as constrained by its primary CMB anisotropy measurement, when hydrostatic mass estimates are used for the cluster measurement (Planck Collaboration XXIV 2016; Planck Collaboration VI 2020). This discrepancy motivated numerous WL studies to improve the mass calibration (e.g., von der Linden et al. 2014; Hoekstra et al. 2015), although most of these calibrations still suggest at least mild tension between the cluster and CMB measurements.

As noted above, a central challenge in cluster cosmology, both for halo abundance and clustering, is accurately calibrating the scaling relations between mass and observables—such as X-ray brightness, optical richness, and SZ flux ($Y_{SZ}$)—which are needed to connect cluster detections based on these observables to the Halo Mass Function (HMF). In general, systematic uncertainties in this calibration limit the precision of the derived cosmological constraints (e.g., Pratt et al. 2019). To improve mass calibration, nearly all recent efforts focus on $M_{WL}$, which is sensitive to the total projected mass distribution and is independent of the dynamical state of the ICM. However, $M_{WL}$ values exhibit significant cluster-to-cluster scatter, along with an average bias, primarily from projection effects due to halo asphericity, orientation, and miscentering (e.g., Meneghetti et al. 2010; Giocoli et al. 2025; Euclid Collaboration: Ragagnin et al. 2025). To quantify this scatter and bias, researchers commonly perform mock WL analyses of simulated clusters with random orientations (e.g., Euclid Collaboration: Giocoli et al. 2024). A limitation of such studies arises if the observable used to detect the cluster sample is sensitive to halo orientation, resulting in a selected population with a preferred geometry. In this case, the scatter and bias in $M_{WL}$ could be misestimated.

Understanding the selected population of galaxy clusters requires quantifying survey completeness, which measures the probability of detecting a cluster based on its observables. For the Planck SZ survey, completeness, $\chi(Y_{5R500}, \theta_{500}, l, b)$, is defined as a function of four cluster parameters, where $Y_{5R500}$ is the integrated SZ signal within $5R_{500}$, $\theta_{500}$ is the corresponding angular scale, and $(l, b)$ are the cluster's galactic coordinates (P16)<fnref n="1" />. Detectability is assessed using the observed signal-to-noise ratio (S/N) obtained from, e.g., the Multi-Matched Filtering (MMF3) algorithm (P16, see also Melin et al. 2006, 2012), which must exceed a threshold value. A robust Monte Carlo (MC) method for characterizing completeness involves injecting a large number of model clusters with a range of values for $Y_{5R500}$ and $\theta_{500}$ into random locations within the Planck maps,

<footnote n="1">¹ $R_{500}$ denotes the spherical radius enclosing an average overdensity 500 times larger than the critical density of the Universe at the cluster redshift.</footnote>
then computing the S/N from the MMF3 for each injected cluster. The Planck team did not consider asphericity in this calculation and injected clusters based only on spherically symmetric models. Subsequent work by Gallo et al. (2024) estimated the Planck completeness using simulated clusters with complicated aspherical shapes more closely resembling reality, finding that steeper SZ signal profiles yield a higher probability of detection. The impact of halo ellipticity was minimal. Overall, the change in completeness obtained by Gallo et al. (2024) results in approximately $1\sigma$ shifts in the derived values of $\Omega_m$ and $\sigma_8$ relative to the nominal Planck result.

In this work, we provide a formalism to investigate how halo triaxiality and orientation impact cluster detection and SZ selection, focusing on the Planck-selected Cluster HERITAGE project with Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) sample. By utilizing smooth cluster models, rather than the simulated clusters of Gallo et al. (2024), we isolate the impact of triaxiality and orientation from other characteristics such as substructure and dynamical state. In addition, because our study is not limited to a relatively small number of simulated objects, we are also able to obtain sufficient statistics to characterize sub-percent level biases. The CHEX-MATE project (CHEX-MATE Collaboration 2021) is an XMM-Newton Heritage Project designed to study ICM physics and improve cluster-based cosmological constraints by refining links between observables and halo mass<fnref n="2" />. It includes 118 clusters, divided into a volume-limited sample (Tier 1, 0.05 < z < 0.2) and a mass-limited sample (Tier 2, z < 0.6, M > 7.25×10¹⁴ M☉), each containing 61 clusters, with four clusters in common. Furthermore, all clusters in the catalog have a Planck S/N > 6.5, and Tier 1 is restricted to the northern sky, Dec > 0°. This deep, uniform, and high-quality multi-probe dataset enables precise measurements of triaxial 3D cluster shapes and orientations (Kim et al. 2024). Our goal is to model the SZ selection function of the CHEX-MATE sample while incorporating cluster triaxiality through informative priors. The end result is a CHEX-MATE completeness estimate that depends on the triaxial geometry, along with a determination of the expected distribution of the geometries for the CHEX-MATE sample. We further explore the impact of this selection on the values of $M_{WL}$ derived for the CHEX-MATE clusters.

This paper is organized as follows. In Section 2, we describe the details of the MMF3 algorithm used to compute cluster detection and selection. In Section 3, we describe the triaxial model assumed for the ICM and the total mass distribution of the clusters. In Section 4, we describe the MC injection and detection procedure used to calculate the completeness for a set of both spherical and triaxial clusters, and we quantify how the completeness depends on the triaxial shape and orientation parameters. In Section 5, we discuss the algorithm for generating a mock CHEX-MATE catalog, beginning with assembling a mock cluster sample from an HMF, injecting triaxial clusters from this sample based on informative priors on their shape and orientation, and subsequently detecting these mock clusters using the MMF3 algorithm. In Section 6, we estimate $M_{WL}$ for these mock CHEX-MATE samples to quantify any WL mass bias, defined as $1 - b_{WL} = M_{WL}/M_{true}$, due to preferentially selected shapes and orientations. Finally, in Section 7, we summarize the main results of this work.

All of the code developed for this work, and some of the data products generated are publicly available at this link.

## 2. Multi-Matched filtering (MMF3) algorithm

The CHEX-MATE sample was selected from the Planck PSZ2 catalog, and we thus follow P16 in their construction of the MMF3. We independently recreated the MMF3 algorithm, with an implementation as close as possible to that used by Planck, with our verification of this implementation detailed below.

The total signal in the Planck maps at each of the High Frequency Instrument (HFI) frequencies at a given position in the sky is modeled as

$$ \mathbf{m}(\mathbf{x}) = y_0 \mathbf{t}_{\theta_s}(\mathbf{x}) + \mathbf{n}(\mathbf{x}), \quad (1) $$

where $\mathbf{x}$ is the 2D pixel coordinate within the map, $\mathbf{n}(\mathbf{x})$ is the unwanted astrophysical signals and instrumental noise, and $\mathbf{t}_{\theta_s}(\mathbf{x})$ is the SZ signal from a cluster with angular size $\theta_s$ and Comptonization parameter $y_0 = k\sigma_T \int n_e T_e dl/(m_e c^2)$. To obtain this SZ signal profile, we constructed the normalized cluster profile, $\tau_{\theta_s}$, which is described by the Generalized Navarro-Frank-White (GNFW) pressure profile (Nagai et al. 2007; Arnaud et al. 2010):

$$ \tau(X) = \frac{P_0}{(c_{500} * X)^\gamma [1 + (c_{500} * X)^\alpha]^{(\beta - \gamma)/\alpha}}, \quad (2) $$

with $P_0$ giving the amplitude, $c_{500} = 1.177$, $\alpha = 1.0510$, $\beta = 5.490$, and $\gamma = 0.3081$. Here, $X = \theta/\theta_{500}$, where $\theta$ is the angular radius and $\theta_s = \theta_{500}/c_{500}$. This 3D profile was then integrated along the LOS to give the projected 2D signal. This 2D map was then convolved with the Planck beam in each HFI channel, $b_i$, which is assumed to have Gaussian profiles as given in P16, with the final observed SZ signal given by

$$ \mathbf{t}_{\theta_s}(\mathbf{x})_i = j_v(v_i)[b_i * \tau_{\theta_s}](\mathbf{x}). \quad (3) $$

The spectral shape of the SZ signal is

$$ j_v(x_v) = x_v T_{CMB} \left[ \frac{e^{x_v} + 1}{e^{x_v} - 1} - 4 \right], \quad (4) $$

where $x_v = \frac{h\nu}{k_B T_{CMB}}$ is the dimensionless frequency. Following P16, we do not consider here the effect of relativistic corrections to the SZ effect. For MMF3, we applied the detection algorithm to square cutouts of the all-sky maps centered on the cluster position, with sides of 2.5°. The noise matrix was estimated directly from these local maps, by mean-subtracting the maps, applying a Hanning window to ensure periodicity at the map edges, and then computing the cross-power spectra for the 6×6 Planck HFI frequencies. An azimuthal averaging procedure was applied, and then this cross-channel matrix was inverted for each Fourier-space pixel, $\mathbf{k}$, defining the inverse noise covariance matrix, $\mathbf{P}^{-1}(\mathbf{k})$. The uncertainty on the SZ signal detection was then estimated from

$$ \sigma_{\theta_s} = \left[ \sum_k t_{\theta_s}(\mathbf{k}) \mathbf{P}^{-1}(\mathbf{k}) t_{\theta_s}(\mathbf{k}) \right]^{-1/2}. \quad (5) $$

The multi-matched filter assumes the spatial and spectral form for the SZ signal defined above. The filter ($\Psi_{\theta_s}$) is uniquely specified by demanding a minimum variance estimate (Melin et al. 2006), and constructed as

$$ \Psi_{\theta_s}(\mathbf{k}) = \sigma_{\theta_s}^2 \mathbf{P}^{-1}(\mathbf{k}) t_{\theta_s}(\mathbf{k}). \quad (6) $$

These are circularly symmetric filters and show typical Fourier ringing, as demonstrated in Melin et al. (2006), as a result of

<footnote n="2">² http://xmm-heritage.oas.inaf.it/</footnote>
down-weighting the large angular scale modes contaminated by primary CMB anisotropy and dust. The filter uses both spatial and frequency weighting to optimally extract the cluster signal from the maps. This allowed us to recover an estimate of the signal, quantified by the central Comptonization parameter as

$$ \hat{y}_0 = \sum_x \Psi_{\theta_s}(x)m(x), \qquad (7) $$

and thus, finally, measuring the S/N of the cluster as $\hat{y}_0/\sigma_{\theta_s}$.

To verify our implementation of the MMF3 algorithm relative to the one utilized by P16, we compared our derived S/N values and integrated SZ signals to their published results for the CHEX-MATE sample. To do this, we first obtained 2.5° × 2.5° cutouts centered on the MMF3 position reported in P16 for each CHEX-MATE cluster from the six all-sky HFI channel Planck maps. We then applied $\Psi_{\theta_s}$ to each of these maps, with the angular scale of the filter spanning $\theta_s$ values between 0.8' and 32'. We defined the cluster size as the filter scale that maximizes the S/N at the location of the cluster, and the best-fit SZ signal was given by the corresponding $\hat{y}_0$ parameter. The integrated signal within a radial extent of $5R_{500}$ was then defined as

$$ Y_{5R_{500}} = \hat{y}_0 \int_{r<5R_{500}} \tau_{\theta_s}(r) r^2 dr. \qquad (8) $$

Figure 1 shows a comparison of our derived values of S/N for the CHEX-MATE sample to those obtained from P16. On average, we find very good agreement, with slight differences due to inherent instabilities in the MMF3 algorithm. At large angular scales, corresponding to small values of $|k|$, the primary CMB anisotropy dominates the lower frequency HFI maps. Since this signal is coherent across frequencies, it introduces strong inter-channel correlations in the noise. Similarly, in the high frequency HFI maps, Galactic dust emission becomes the dominant signal, and this emission is also spatially coherent and highly correlated across these frequencies. Due to these large correlations, even minor differences in how the noise is modeled can propagate nonlinearly when inverting the noise covariance matrix, resulting in measurable discrepancies in the final S/N. These discrepancies can be quantified by the difference between S/N values obtained from our MMF3 implementation and the one utilized by P16. For the CHEX-MATE sample, the median difference is 0.21, with a standard deviation of 0.72, see Fig. 1. Compared to a Gaussian distribution with this mean and scale factor, a Kolmogorov-Smirnov (KS) test yields a p-value of 0.56, indicating that the observed difference in S/N is consistent with this distribution. We note that this level of agreement is approximately a factor of two better than that found between the detection algorithms described in P16, suggesting that our implementation of the MMF3 is consistent with the one utilized by P16. Furthermore, applying our MMF3 algorithm results in only three out of 118 CHEX-MATE clusters falling below the detection threshold of S/N > 6.5. As a result, we expect the overall impact on sample selection from our MMF3 implementation relative to that of P16, to be minimal.

### 3. Triaxial cluster shapes

The matter distribution within galaxy clusters is better approximated by a triaxial ellipsoid than by a sphere (Despali et al. 2014; Bonamigo et al. 2015; Vega-Ferrero et al. 2017), primarily due to the anisotropic accretion of matter along cosmic filaments and the complex dynamics of structure formation.

<figure><img src="image_2.png" /><figcaption>Fig. 1. Histogram of the difference in S/N obtained from our implementation of the MMF3 detection algorithm and that used by P16 for the CHEX-MATE clusters.</figcaption></figure>

Simulations suggest that modeling clusters as 3D triaxial ellipsoids results in lower bias and less intrinsic scatter in estimating intrinsic quantities (Becker & Kravtsov 2011). Observations also strongly support the triaxial nature of galaxy clusters. Ellipticity has been detected in a variety of observational probes: X-ray surface brightness maps reveal elongated morphologies (Lau et al. 2012), SZ maps show asymmetric pressure distributions (Sayers et al. 2011), and weak lensing reconstructions uncover elliptical mass distributions (Oguri et al. 2012), among others. These complementary data sets highlight the motivation for 3D triaxial modeling in both theoretical and observational contexts. As examples of such modeling, Limousin et al. (2013) present a general parametric framework within a triaxial basis to simultaneously fit complementary data sets from X-ray, SZ and lensing, and Sereno et al. (2018) demonstrate cluster shape constraints from a multi-probe 3D analysis of the X-ray regular Cluster Lensing and Supernova survey with Hubble (CLASH) clusters, using a combined analysis of strong and weak lensing, X-ray photometry and spectroscopy, and SZ.

Following these works, we assumed a geometry described by a triaxial ellipsoid with the ICM pressure following a distribution given by a radial profile in this triaxial basis, consistent with the formalism described in Kim et al. (2024). The main parameters involved in the construction of this profile are related to the intrinsic versus observed coordinate system of a triaxial ellipsoid. The projection of the ellipsoid on the sky plane is defined by an ellipse with a semi-major axis $l_{p1}$ and eccentricity $q_p$, while the projection of the ellipsoid along the LOS is defined by a semi-axis $l_{los}$. The parameters $l_{p1}$, $q_p$, and $l_{los}$ have complicated descriptions in terms of the viewing angle and geometry of the ellipsoid ($q_1, q_2, \cos\theta, \phi, \psi$), where $q_1$ and $q_2$ are the minor-to-major and intermediate-to-major axial ratios of the ellipsoid, respectively, with the angles defined as the Euler angles between the intrinsic and observed coordinate system. These are shown in Fig. 2 for clarity.

We modeled the pressure as a radial profile in the elliptical basis parameterized by the ellipsoidal radius $\zeta$, where $\zeta^2 = (x_1/q_1)^2 + (x_2/q_2)^2 + x_3^2$, which is used in place of the spherical radius in the GNFW model given in Eq. (2). Within this framework, the 2D maps were calculated from the 3D ellipsoidal model projected onto the sky plane. The model Compton-y
<figure><img src="image_3.png" /><figcaption>Fig. 2. Triaxial ellipsoid model used in the analysis. The intrinsic coordinate system of the ellipsoid is indicated by the dotted arrows, and the black arrows correspond to the observer's coordinate system. The red ellipse shows the projection of the ellipsoid on the sky plane, with $l_{p1}$ as its semi-major axis. The green ellipse shows the projection of the ellipsoid onto the plane perpendicular to the sky plane, and $l_{los}$ is the half size of the ellipse along the observer line of sight (Kim et al. 2024).</figcaption></figure>

parameter is

$$y_{\text{model}}(x_{\xi}; l_{p1}) = (2l_{p1}e_{\parallel}) \left(\frac{\sigma_{\text{T}}}{m_e c^2}\right) \int_{x_{\xi}}^{\infty} P_e(x_{\xi}) \frac{x_{\xi}}{\sqrt{x_{\xi}^2 - x_{\xi}^2}} dx_{\xi}, \quad (9)$$

where $y$ has been described in terms of the elliptical radius on the sky $x_{\xi}^2 = (x_1^2 + (x_2/q_p)^2)(l_s/l_{p1})^2$, $x_{\xi} = \zeta/l_s$, $q_p$ is the minor-to-major axial ratio of the observed projected isophote, $l_s$ is the semi-major axis of the ellipsoid, $e_{\parallel} = l_{los}/l_{p1}$ is the elongation parameter of the ellipsoid and $P_e$ is obtained from the GNFW profile

$$P_e(x_{\xi}) = \frac{P_0 P_{500}}{\left(c_{500} x_{\xi} \frac{l_s}{R_{500}}\right)^{\gamma_p} \left[1 + \left(c_{500} x_{\xi} \frac{l_s}{R_{500}}\right)^{\alpha_p}\right]^{(\beta_p - \gamma_p)/\alpha_p}}, \quad (10)$$

with the same GNFW profile parameters as used in Eq. (2). Thus, a 3D triaxial model of a cluster was generated and then projected into a 2D map based on its geometrical shape parameters ($q_1, q_2, \cos\theta, \phi, \psi$) and its pressure given by the GNFW model.

We defined the total SZ signal $Y_{5R500}$ and angular extent $\theta_{500}$ for a triaxial cluster by computing the spherically-averaged radial profile of the 3D triaxial model. The spherical GNFW model was then fit to this spherically-averaged profile model to determine the best-fit spherical $\theta_{500}$, which defined the angular scale of that cluster. Similarly, $Y_{5R500}$ ($Y_{500}$) was computed by integrating this spherically-averaged pressure profile within $5\theta_{500}$ ($\theta_{500}$). With this procedure, we ensured that the total integrated SZ signal $Y_{5R500}$ was the same regardless of triaxial geometry. Because we defined the elliptical radius relative to the major axis, the spherically-averaged $\theta_{500}$ is always smaller than its elliptical equivalent in Eq. (10).

## 4. Selection function

Estimating the Planck completeness requires introducing realistic noise into mock observations of cluster models with

<figure><img src="image_4.png" /><figcaption>Fig. 3. Comparison of our MC completeness, averaged over the unmasked sky, to the semi-analytical ERF completeness for $\theta_{500} = 1.19', 6.11', 17.51',$ and $31.41'$, assuming an S/N threshold of 8.5.</figcaption></figure>

known properties, such as $Y_{5R500}$ and $\theta_{500}$. To obtain representative noise, we directly sampled the Planck all-sky maps, and then applied an MC method to compute completeness. We masked regions according to the Planck 2015 cosmology mask (Planck Collaboration XI 2016), along with the point source mask (Planck Collaboration XXVI 2016), and an additional exclusion mask extending 1° from the center of all clusters in the PSZ2 catalog of P16. From the remaining unmasked sky regions, we randomly selected 100 locations, each separated by at least 5° in order to sample a representative range of noise environments. At each location, we injected clusters with varying $Y_{5R500}$ and $\theta_{500}$ values and determined the associated S/N from the MMF3 algorithm. Our assumed model for the cluster does not include any signal from correlated large-scale structure and, in general, neither did the random sky locations. Thus, the impact of such structure was not included in the derived S/N. The completeness function was then defined as the detection probability curve, which relates the likelihood of detection above a S/N threshold of 6.5 to $Y_{5R500}$ and $\theta_{500}$.

To validate our completeness pipeline, we first performed a direct comparison to the completeness function found by P16 for spherical clusters. We find good consistency, including percentile agreement with semi-analytical calculations obtained from the error function (ERF):

$$P(d|Y_{500}, \sigma_Y(\theta_{500}), q) = \frac{1}{2} \left[ 1 + \text{erf} \left( \frac{Y_{500} - q\sigma_Y(\theta_{500})}{\sqrt{2}\sigma_Y(\theta_{500})} \right) \right], \quad (11)$$

where $q$ is the detection threshold and $\sigma_Y(\theta_{500})$ is the estimate of the noise at that angular scale, see Fig. 3.

Following this validation, we extended our MC completeness calculation to triaxial clusters, constructing completeness curves to quantify how detection probability varies with triaxial axial ratios and orientation. In addition to the masks noted above, these samples included additional constraints: Tier 1 was limited to $\text{Dec} > 0^\circ$, while both Tier 1 and Tier 2 had a low XMM-Newton visibility mask to exclude regions with visibilities <55 ks per orbit (CHEX-MATE Collaboration 2021). From the unmasked sky regions, we again selected 100 widely separated map cutouts for both Tier 1 and Tier 2, ensuring good sampling of the available sky. We then inserted triaxial clusters with varying axial ratios, orientations, $Y_{5R500}$, and $\theta_{500}$ into these cutouts and computed completeness by averaging detection probabilities over all sky positions. To model triaxial shape and orientation
<figure><img src="image_5.png" /><figcaption>Fig. 4. Average MMF3 S/N as a function of $e'_\Delta$ for different triaxial axial ratios, for a cluster with an effective spherical $\theta_{500} = 5.63'$ and a fixed value of $Y_{5R500}$.</figcaption></figure>

effects on MMF3 S/N, we introduced a new parameter,

$$e'_{\Delta} = \frac{l_{\text{los}}}{(l_{\text{p1}} * l_{\text{p2}})^{0.5}}, \qquad (12)$$

where, as before, $l_{\text{los}}$ is the semi-major axis of the ellipse along the LOS and $(l_{\text{p1}} * l_{\text{p2}})^{0.5}$ is the effective spherical radius of the projected ellipse on the sky (see Fig. 2). We find that $e'_\Delta$ exhibits a complex dependence on $q_1, q_2, \cos\theta$, and $\phi$, but correlates monotonically with S/N–higher values of $e'_\Delta$ correspond to higher S/N for a fixed $(q_1, q_2)$, as shown in Fig. 4. We established this correlation for the full range of geometries considered in this paper, but we show only a subset in Fig. 4 for clarity. Thus, rather than using two separate orientation parameters ($\cos\theta, \phi$), we consolidate them into $e'_\Delta$ for each $q_1, q_2$. Using this framework, we constructed a comprehensive lookup table that quantifies detection probability as a function of SZ parameters and triaxial shape, thus describing $P(d|Y_{500}, \theta_{500}, q_1, q_2, e'_\Delta)$.

We present three slices of the five-parameter completeness function in Fig. 5, illustrating completeness as a function of triaxial axial ratios and orientations. We find a clear dependence on orientation, where clusters elongated along the LOS exhibit a higher probability of detection for the same $Y_{5R500}, \theta_{500}$, and shape parameters. Additionally, we observe that for the same effective spherical $\theta_{500}$ but for other orientations, there is no clear trend between detectability and cluster triaxiality. We thus find that the primary factor influencing detectability is the extent of the cluster along the LOS. For triaxial clusters whose major axes are not closely aligned with the LOS, varying triaxialities can still result in more or less centrally peaked SZ profiles relative to spherical clusters, due to projection effects. Consistent with the results reported by Gallo et al. (2024), we find that the more centrally peaked clusters from this projection effect are more easily detected.

To ensure a well-sampled and smoothly varying lookup table connecting detection probability with the relevant parameters, we fitted the completeness function for each $(\theta_{500}, q_1, q_2, e'_\Delta)$ to an ERF of the form:

$$P(d) = \frac{1}{2} \left[ 1 + \text{erf} \left( \frac{Y_{500} - q * a(R_{500}, q_1, q_2, e'_\Delta)}{\sqrt{2} * b(R_{500}, q_1, q_2, e'_\Delta)} \right) \right]. \qquad (13)$$

Relative to the ERF in P16, this form incorporates modified mean (a) and variance (b) terms to capture non-Gaussian deviations introduced by cluster triaxiality. We then performed a grid interpolation over our lookup table to determine a and b as functions of $q_1, q_2$, and $e'_\Delta$.

## 5. Mock CHEXMATE catalogs

We used the triaxial completeness given by the ERF in Eq. (13) to construct mock CHEXMATE catalogs, allowing us to quantify the average distribution of orientations and geometries that arise from this sample selection. First, we generated mock universes with halo populations sampled from the Tinker et al. (2008) HMF, assuming the Planck Collaboration VI (2020) cosmology, considering redshifts $0 \le z \le 0.6$ and $M_{500}$ masses between $10^{14} M_\odot$ and $10^{16} M_\odot$, while also applying the sky fraction constraints for both Tier 1 and Tier 2 samples. We converted the generated halo masses ($M_{500}$) and redshifts ($z$) into $\theta_{500}$ and $Y_{500}$ using the established scaling relations from Planck Collaboration XX (2014):

$$E^{-\beta}(z) \left[ \frac{D_A^2(z) \bar{Y}_{500}}{10^{-4} \text{Mpc}^2} \right] = Y_* \left[ \frac{h}{0.7} \right]^{-2+\alpha} \left[ \frac{(1-b)M_{500}}{6 \times 10^{14} M_\odot} \right]^\alpha \qquad (14)$$

and

$$\bar{\theta}_{500} = \theta_* \left[ \frac{h}{0.7} \right]^{-2/3} \left[ \frac{(1-b)M_{500}}{3 \times 10^{14} M_{\odot}} \right]^{1/3} E^{-2/3}(z) \left[ \frac{D_A(z)}{500 \text{Mpc}} \right]^{-1}, \qquad (15)$$

where $E(z) = H(z)/H_0$, $D_A(z)$ is the angular diameter distance and $(1-b)$ is the mass bias factor, which is set to 0.62 for consistency with the cosmology obtained by Planck Collaboration VI (2020). The values of the scaling and normalization parameters were obtained from Planck Collaboration XXIV (2016). We then performed volume selections to form mock Tier 1 ($0.05 < z < 0.2$) and Tier 2 ($z < 0.6$) catalogs from halos generated using the HMF. To assign triaxial shapes and orientations to halos, we randomly draw axial ratio values for the gas from the empirical distributions obtained from the IllustrisTNG simulations presented in Machado Poletti Valle et al. (2021), along with random draws for $\cos\theta \in [0, 1]$ and $\phi \in [0^\circ, 90^\circ]$. Then, using the calibrated ERF, we computed the completeness for each halo and included it in the resulting catalog based on probabilistic sampling. Specifically, for each halo, we drew a uniform random value between 0 and 1, and included the halo in the catalog if the value of $P(d)$ was larger than this random value. We repeated this process for 1000 mock universes for both Tier 1 and Tier 2.

Since Tier 2 requires an additional selection based on MMF3-derived mass, we reprocessed the completeness-selected clusters through the MMF3 algorithm. This provided the extracted $Y_{500}$ for each cluster at every filter scale $\theta_{500}$ of MMF3. From this $Y-\theta$ relation we then determined the value of $M_{\text{MMF3}}$ that simultaneously satisfied Eq. (14) and (15). Since this process was repeated across 100 distinct map cutouts, we obtained 100 different values of $M_{\text{MMF3}}$ for each cluster, reflecting the variation in recovered mass due to noise fluctuations. To assign a single value of $M_{\text{MMF3}}$ to each cluster, we performed a probabilistic selection based on the cluster's detection completeness. Specifically, we randomly selected a value of $M_{\text{MMF3}}$ from among the $N$ highest values obtained from the 100 cutouts, where $N = P(d) \times 100$. This method effectively models the S/N-based mass boosting, as clusters with lower intrinsic S/N values are only included in the catalog when they fall within favorable
<figure><img src="image_6.png" /><figcaption>Fig. 6. Mass and redshift distributions for the observed Tier 1 and Tier 2 CHEX-MATE samples compared to the kernel density estimation (KDE) distribution of our 1000 mock CHEX-MATE selections. The distributions are consistent, indicating that the mock samples accurately reproduce the masses and redshifts of the observed samples.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 7. Distributions of orientation angles cos θ and φ for the mock Tier 1 and Tier 2 samples, along with randomly selected samples with similar distributions in mass and redshift. The CHEX-MATE samples show a clear bias toward larger values of cos θ and smaller values of φ compared to the random samples, indicating that SZ selection favors clusters elongated along the LOS.</figcaption></figure>

utilized for cosmological analyses, performed via a least-squares approach. Four different fits were performed using different values of c200, including fixing its value to either three or four, fixing its value based on the concentration-mass relation from Diemer & Joyce (2019), and allowing its value to float (within the range of one to ten). We performed these fits for the mock Tier 1 and Tier 2 samples, along with the randomly selected samples matching their redshift and mass distributions. We then computed the average ratio of the fitted M<sub>WL</sub> and the true halo mass, which we defined as 1-b<sub>WL</sub> = M<sub>WL</sub>/M<sub>True</sub>, within discrete mass bins for all cases (see Fig. 8). For both Tier 1 and Tier 2, the M<sub>WL</sub> of the CHEX-MATE-selected samples is higher than that obtained from the random samples at the low-mass end, trending toward agreement at the high-mass end. This is as expected, since the up-scatter in S/N due to clusters elongated along LOS is more pronounced in lower-mass (and thus lower S/N) clusters, with higher-mass clusters generally having a higher S/N, which allows for detection irrespective of orientation. Furthermore, for c = 3 and c = 4, we see a sharp increase in M<sub>WL</sub> for both the random and CHEX-MATE-selected samples at the high-mass end. This is a result of utilizing the relation of Diemer & Joyce (2019)

when generating the halos, as they find an increase in c200 at the high-mass end.

We note that the overall value of 1 − b<sub>WL</sub> depends strongly on our fitting methodology, i.e., the assumed prior on c200, as expected due to the strong anti-correlation between c and M. However, the ratio of 1 − b<sub>WL</sub> obtained from the CHEX-MATE-selected and random samples, which we define as 1 − b<sub>χ</sub>, appears to be independent of the fitting methodology, as illustrated in Fig. 9. In other words, the portion of 1 − b<sub>WL</sub> associated with the on-average orientation of an SZ selection is robust to the choice of fitting method. For both Tier 1 and Tier 2, we find that 1 − b<sub>χ</sub> is approximately 4% at the low-mass end, with an overall average value near 2%. Thus, for the CHEX-MATE samples– and likely any Planck SZ-selected sample–the correction for this bias in 1 − b<sub>WL</sub> due to orientation should be universal, regardless of the WL fitting methodology. To facilitate such a correction, we provide the parameters of a linear regression fit of the form

$$1 - b_{\chi} = \begin{cases} k \times (M_{\text{True}} - M_{\text{break}})/10^{14} M_{\odot} & \text{if } M_{\text{True}} < M_{\text{break}} \\ 0 & \text{if } M_{\text{True}} \ge M_{\text{break}} \end{cases}$$
<figure><img src="image_8.png" /><figcaption>**Fig. 8.** Comparison of $1 - b_{WL}$ as a function of true mass for the mock CHEX-MATE Tier 1 (top) and tier 2 (bottom) selected samples (purple lines and shaded 68% confidence level regions), compared to a random distribution of halos with matched mass and redshift distributions (orange lines and shaded 68% confidence level regions), using varying priors on the value of $c_{200}$ to obtain $M_{WL}$ from the mock shear data. Left: $c_{200} = 4$. Middle left: $c_{200} = 3$. Middle right: $c_{200}$ obtained from the $c-M$ relation of Diemer & Joyce (2019). Right: $c_{200}$ varied as a free parameter.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Ratio of $1 - b_{WL}$ for CHEX-MATE-selected samples (Left: Tier 1, Right: Tier 2), to $1 - b_{WL}$ for random samples with matching mass and redshift distributions, defined as $1 - b_{\chi}$. This plot isolates the bias due solely to orientation in an SZ-selected sample. This orientation bias is approximately independent of fitting methodology, as illustrated by the good agreement between all four choices used in this work. The black line shows a broken linear fit that reproduces the derived values of $1 - b_{\chi}$ to within an rms scatter of 0.003.</figcaption></figure>

to describe the selection-based orientation-induced bias in $M_{WL}$ for both the Tier 1 and Tier 2 samples, with the results summarized in Table 1. We find that these linear fits accurately reproduce the derived bias to within 0.3%, which is sufficient for percent-level mass calibration. Because of the different mass and redshift ranges for Tier 1 and Tier 2, along with the additional selection based on MMF3-derived mass for Tier 2, there is no expectation that the fitted slope and break point should be con-

sistent between the two samples, and we indeed find modest differences in their values.

## 7. Conclusions

In this paper, we investigated the impact of triaxiality and orientation on SZ selection and the associated $M_{WL}$ values, focusing on the specific example of CHEX-MATE. Quantifying these
<table><caption>**Table 1.** Parameters of the linear fit of $1 - b_\chi$ for the Tier 1 and Tier 2 CHEX-MATE samples, where $1 - b_\chi$ represents the bias in $M_{WL}$ due solely to orientation bias from the *Planck* SZ selection.</caption><thead><tr><th>CHEX-MATE</th><th><i>M</i><sub>break</sub></th><th><i>k</i></th></tr></thead><tbody><tr><td>Tier 1</td><td>15.03 × 10<sup>14</sup> M<sub>☉</sub></td><td>-0.0039</td></tr><tr><td>Tier 2</td><td>21.67 × 10<sup>14</sup> M<sub>☉</sub></td><td>-0.0026</td></tr></tbody></table>

selection effects is crucial for future cluster cosmology studies, as it enables corrections to the $M_{WL}$ values used for calibration of the halo masses.

To achieve this goal, we recreated the MMF3 algorithm used by P16 for the detection of clusters based on their SZ signal in the *Planck* all-sky maps. Using an MC method based on the injection of mock clusters into randomly sampled cutouts from *Planck* HFI maps within the unmasked region of the sky, we confirm that the values of S/N and $Y_{5R500}$ obtained from our implementation of the MMF3 are consistent with those obtained by P16 for spherically-symmetric clusters. By then applying this algorithm to mock observations of clusters with triaxial shapes, we determined the completeness as a function of SZ signal strength, angular size, and triaxial shape and orientation.

Using our measured triaxial completeness, we constructed mock CHEX-MATE catalogs to quantify selection biases introduced by triaxiality and orientation for SZ-detected clusters. We find a clear orientation bias, where clusters elongated along the LOS are preferentially detected due to their centrally boosted SZ signals. By fitting $M_{WL}$ to the mock Tier 1 and Tier 2 detected clusters, along with randomly selected clusters with comparable mass and redshift distributions, we identify a mass-dependent $M_{WL}$ bias due solely to the shape and orientation of the SZ-selected samples, $1 - b_\chi$. This bias reaches up to 4% in the lowest-mass clusters and averages approximately 2% across the full mass range. In addition, $1 - b_\chi$ appears to be independent of the exact WL fitting methodology. We also note that the value of $1 - b_\chi$ found here is expected to be a lower limit, since we do not account for correlated LSS in our selection or WL modeling. Such structures are expected to occur preferentially along the halo elongation direction, and would thus boost the SZ signal when viewed in that orientation, resulting in a stronger preference for LOS elongation in the CHEX-MATE-selected samples, as well as boosting the value of $M_{WL}$ at fixed LOS elongation.

We also consider the implications of our results in the context of cluster cosmology from halo abundance measurements. First, we find that $M_{WL}$ is systematically overestimated in SZ-selected samples from *Planck*, implying that true masses are, on average, lower than WL-derived masses for such samples. The discrepancy between halo abundance and primary CMB measurements from *Planck* requires higher halo masses for the cluster sample, so this correction acts in the opposite direction. In addition, we note that the average 2% systematic overestimate of $M_{WL}$ measured here is subdominant to current systematic uncertainties in $M_{WL}$ calibration. Thus it does not significantly impact previously published results. However, upcoming surveys like CMB-S4 and *Euclid* require percent-level mass calibration (Raghunathan et al. 2022; *Euclid Collaboration: Giocoli et al.* 2024; *Euclid Collaboration: Ragagnin et al.* 2025), and therefore such biases must be accounted for in their analyses.

*Acknowledgements.* This research was supported by the International Space Science Institute (ISSI) in Bern, through ISSI International Team project #565

(Multi-Wavelength Studies of the Culmination of Structure Formation in the Universe). H.S., J.S., A.G. and J.K. were supported by NASA Astrophysics Data Analysis Program (ADAP) Grant 80NSSC21K1571. SE, MS, MR acknowledges the financial contribution from the contracts Prin-MUR 2022 supported by Next Generation EU (M4.C2.1.1, n.20227RNLY3 The concordance cosmological model: stress-tests with galaxy clusters). MS acknowledges the financial contributions from contract INAF mainstream project 1.05.01.86.10 and INAF Theory Grant 2023: Gravitational lensing detection of matter distribution at galaxy cluster boundaries and beyond (1.05.23.06.17). SE, MR acknowledges the financial contributions from the European Union's Horizon 2020 Programme under the AHEAD2020 project (grant agreement n. 871158). GWP acknowledges long-term support from CNES, the French space agency. LL acknowledges the financial contribution from the INAF grant 1.05.12.04.01. EP acknowledges support from CNES, the French national space agency and from ANR the French Agence Nationale de la Recherche, under grant ANR-22-CE31-0010. J.K. acknowledges the support by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2019R1A6A1A10073887) and the 2025 KAIST-U.S. Joint Research Collaboration Open Track Project for Early-Career Researchers, supported by the International Office at the Korea Advanced Institute of Science and Technology (KAIST). M.G. acknowledges support from the ERC Consolidator Grant BlackHoleWeather (101086804). MDP acknowledges financial support from PRIN-MUR grant 20228B938N "Mass and selection biases of galaxy clusters: a multi-probe approach" funded by the European Union Next generation EU, Mission 4 Component 1 CUP C53D2300092 0006. DE acknowledges support from the Swiss National Science Foundation (SNSF) through grant agreement #200021_212576. BJM acknowledges support from Science and Technology Facilities Council grants ST/V000454/1 and ST/Y002008/1.

## References

*Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2020, Phys. Rev. D, 102, 023509*

*Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARA&A, 49, 409*

*Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92*

*Becker, M. R., & Kravtsov, A. V. 2011, ApJ, 740, 25*

*Bellagamba, F., Roncarelli, M., Maturi, M., & Moscardini, L. 2017, MNRAS, 473, 5221*

*Bleem, L. E., Stalder, B., de Haan, T., et al. 2015, ApJS, 216, 27*

*Bleem, L. E., Klein, M., Abbot, T. M. C., et al. 2024, Open J. Astrophys., 7, 13*

*Bocquet, S., Grandis, S., Bleem, L. E., et al. 2024, Phys. Rev. D, 110, 083510*

*Bonamigo, M., Despali, G., Limousin, M., et al. 2015, MNRAS, 449, 3171*

*CHEX-MATE Collaboration (Arnaud, M., et al.) 2021, A&A, 650, A104*

*Chiu, I. N., Ghirardini, V., Liu, A., et al. 2022, A&A, 661, A11*

*Costanzi, M., Rozo, E., Simet, M., et al. 2019, MNRAS, 488, 4779*

*Costanzi, M., Saro, A., Bocquet, S., et al. 2021, Phys. Rev. D, 103, 043522*

*DES Collaboration (Abbott, T. M. C.) 2025, ArXiv e-prints [arXiv:2503.13632]*

*Despali, G., Giocoli, C., & Tormen, G. 2014, MNRAS, 443, 3208*

*Despali, G., Giocoli, C., Angulo, R. E., et al. 2016, MNRAS, 456, 2486*

*Diemer, B., & Joyce, M. 2019, ApJ, 871, 168*

*Euclid Collaboration (Giocoli, C., et al.) 2024, A&A, 681, A67*

*Euclid Collaboration (Ragagnin, A., et al.) 2025, A&A, 695, A282*

*Gallo, S., Douspis, M., Soubrié, E., & Salvati, L. 2024, A&A, 686, A15*

*Ghirardini, V., Bulbul, E., Artis, E., et al. 2024, A&A, 689, A298*

*Giocoli, C., Despali, G., Meneghetti, M., et al. 2025, A&A, 697, A184*

*Grandis, S., Ghirardini, V., Bocquet, S., et al. 2024, A&A, 687, A178*

*Hilton, M., Sifón, C., Naess, S., et al. 2021, ApJS, 253, 3*

*Hoekstra, H., Herbonnet, R., Muzzin, A., et al. 2015, MNRAS, 449, 685*

*Kim, J., Sayers, J., Sereno, M., et al. 2024, A&A, 686, A97*

*Kleinebreil, F., Grandis, S., Schrabback, T., et al. 2025, A&A, 695, A216*

*Kornelje, K., Bleem, L. E., Rykoff, E. S., et al. 2025, ApJ, submitted [arXiv:2503.17271]*

*Lau, E. T., Nagai, D., Kravtsov, A. V., & Zentner, A. R. 2011, ApJ, 734, 93*

*Lau, E. T., Nagai, D., Kravtsov, A. V., Vikhlinin, A., & Zentner, A. R. 2012, ApJ, 755, 116*

*Lesci, G. F., Marulli, F., Moscardini, L., et al. 2022a, A&A, 659, A88*

*Lesci, G. F., Nanni, L., Marulli, F., et al. 2022b, A&A, 665, A100*

*Limousin, M., Morandi, A., Sereno, M., et al. 2013, Space Sci. Rev., 177, 155*

*Machado Poletti Valle, L. F., Avestruz, C., Barnes, D. J., et al. 2021, MNRAS, 507, 1468*

*Mantz, A., Allen, S. W., Rapetti, D., & Ebeling, H. 2010, MNRAS, 406, 1759*

*Mantz, A. B., von der Linden, A., Allen, S. W., et al. 2015, MNRAS, 446, 2205*

*Melin, J.-B., Bartlett, J. G., & Delabrouille, A., et al. 2006, A&A, 459, 341*

*Melin, J. B., Aghanim, N., Bartelmann, M., et al. 2012, A&A, 548, A51*

*Meneghetti, M., Rasia, E., Merten, J., et al. 2010, A&A, 514, A93*

*Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007, ApJ, 668, 1*

*Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493*
Oguri, M., Bayliss, M. B., Dahle, H., et al. 2012, MNRAS, 420, 3213

Okabe, N., Reiprich, T. H., Grandis, S., et al. 2025, A&A, 700, A46

Pacaud, F., Pierre, M., Melin, J.-B., et al. 2018, A&A, 620, A10

Planck Collaboration XX. 2014, A&A, 571, A20

Planck Collaboration XI. 2016, A&A, 594, A11

Planck Collaboration XXIV. 2016, A&A, 594, A24

Planck Collaboration XXVI. 2016, A&A, 594, A26

Planck Collaboration XXVII. 2016, A&A, 594, A27

Planck Collaboration VI. 2020, A&A, 641, A6

Pratt, G. W., Arnaud, M., Biviano, A., et al. 2019, Space Sci. Rev., 215, 25

Raghunathan, S., Whitehorn, N., Alvarez, M. A., et al. 2022, ApJ, 926, 172

Rozo, E., Wechsler, R. H., Rykoff, E. S., et al. 2010, ApJ, 708, 645

Rykoff, E. S., Rozo, E., Busha, M. T., et al. 2014, ApJ, 785, 104

Sarazin, C. L. 1986, Rev. Mod. Phys., 58, 1

Sayers, J., Golwala, S. R., Ameglio, S., & Pierpaoli, E. 2011, ApJ, 728, 39

Sehgal, N., Trac, H., Acquaviva, V., et al. 2011, ApJ, 732, 44

Sereno, M., Umetsu, K., Ettori, S., et al. 2018, ApJ, 860, L4

Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comments Astrophys. Space Phys., 4, 173

Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709

To, C., Krause, E., Rozo, E., et al. 2021, Phys. Rev. Lett., 126, 141301

Vega-Ferrero, J., Yepes, G., & Gottlöber, S. 2017, MNRAS, 467, 3226

Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 692, 1060

Voit, G. M. 2005, Rev. Mod. Phys., 77, 207

von der Linden, A., Mantz, A., Allen, S. W., et al. 2014, MNRAS, 443, 1973