# CHEX-MATE: exploring the kinematical properties of Planck galaxy clusters

Lorenzo Pizzuti, Rafael Barrena, Mauro Sereno, Alina Streblyanska, Antonio Ferragamo, Sophie Maurogordato, Alberto Cappi, Stefano Ettori, Gabriel W Pratt, Gianluca Castignani, et al.

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

To cite this version:

Lorenzo Pizzuti, Rafael Barrena, Mauro Sereno, Alina Streblyanska, Antonio Ferragamo, et al. CHEX-MATE: exploring the kinematical properties of Planck galaxy clusters. Astronomy & Astrophysics - A&A, 2025, 699, pp.A88. 10.1051/0004-6361/202555417. hal-05077881

HAL Id: hal-05077881

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

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# CHEX-MATE: Exploring the kinematical properties of Planck galaxy clusters*

Lorenzo Pizzuti<fnref n="1" />, ** ID, Rafael Barrena<fnref n="2" />,<fnref n="3" /> ID, Mauro Sereno<fnref n="4" />,<fnref n="5" /> ID, Alina Streblyanska<fnref n="2" />,<fnref n="3" /> ID, Antonio Ferragamo<fnref n="6" /> ID, Sophie Maurogordato<fnref n="7" /> ID, Alberto Cappi<fnref n="4" />,<fnref n="7" /> ID, Stefano Ettori<fnref n="4" />,<fnref n="5" /> ID, Gabriel W. Pratt<fnref n="8" />, Gianluca Castignani<fnref n="4" /> ID, Megan Donahue<fnref n="9" />, Dominique Eckert<fnref n="4" />,<fnref n="5" /> ID, Fabio Gastaldello<fnref n="10" />, Raphael Gavazzi<fnref n="11" />,<fnref n="12" /> ID, Christopher P. Haines<fnref n="13" />,<fnref n="14" /> ID, Scott T. Kay<fnref n="15" /> ID, Lorenzo Lovisari<fnref n="10" />,<fnref n="16" /> ID, Ben J. Maughan<fnref n="17" /> ID, Etienne Pointecouteau<fnref n="18" /> ID, Elena Rasia<fnref n="19" />,<fnref n="20" />,<fnref n="21" /> ID, Mario Radovich<fnref n="22" /> ID, and Jack Sayers<fnref n="23" />

Received 7 May 2025 / Accepted 23 May 2025

## ABSTRACT

We analysed the kinematical properties of the CHEX-MATE galaxy cluster sample. Our study is based on the radial velocities retrieved from the SDSS DR18, DESI, and NED spectroscopic databases and new data obtained with the 10.4 m GTC and ESO-NTT telescopes. We derived cluster mass profiles for 75 clusters using the MG-MAMPOSST procedure, which recovers the gravitational potential and the anisotropy profiles from line-of-sight velocities and projected positions of galaxy members. The standard Navarro–Frenk–White (NFW) model and the Burkert model, with flatter cores than the NFW, both adequately fit the kinematic data, with only a marginal statistical preference for one model over the other. An estimation of the mass bias ($1 - B_1$) = $M_{500}^{SZ}/M_{500}^M$ was performed via a comparison with Sunyaev-Zel'dovich–X-ray-calibrated mass estimates, resulting in a value of 0.54 ± 0.11 when four evidently disturbed clusters are removed from the sample. We assessed the dynamical state of the clusters by inferring the Anderson-Darling coefficient ($A^2$) and the fraction of galaxies in substructures ($f_{sub}$). Except for a few cases, we find relatively low values for $A^2$, which suggests that CHEX-MATE clusters are not too far from relaxation. Moreover, no significant trends emerge between $A^2$ and $f_{sub}$, nor between the log-masses estimated by MG-MAMPOSST and those based on the Sunyaev–Zel'dovich effect calibrated through X-rays measurements. We studied the concentration–mass relation for the sample; despite the large scatter, we observe signs of an increasing trend for high-mass clusters, in agreement with recent theoretical expectations. Finally, our analysis of the radial anisotropy profiles of member galaxies – stacked in five bins of mass and redshift – reveals that orbits tend to be isotropic at the centre and more radial towards the edge, as found in previous studies. A slight trend of increasing radial orbits at $r_{200}$ is observed in clusters with larger velocity dispersions.

**Key words.** galaxies: clusters: general – galaxies: kinematics and dynamics

<footnote n="1">¹ Dipartimento di Fisica G. Occhialini, Università degli Studi di Milano Bicocca, Piazza della Scienza 3, I-20126 Milano, Italy</footnote>

<footnote n="2">² Instituto de Astrofísica de Canarias, C/ Vía Láctea s/n, E-38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="3">³ Universidad de La Laguna, Departamento de Astrofísica, E-38206 La Laguna, Tenerife, Spain</footnote>

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

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

<footnote n="6">⁶ Dipartimento di Fisica 'E. Pancini', Università degli Studi di Napoli Federico II, Via Cinthia, 21, I-80126 Napoli, Italy</footnote>

<footnote n="7">⁷ Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire Lagrange, Bd de l'Observatoire, CS 34229, 06304 Nice Cedex 4, France</footnote>

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

<footnote n="9">⁹ Michigan State University, Physics and Astronomy Department, 567 Wilson Road, East Lansing, Michigan 48824, USA</footnote>

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

<footnote n="11">¹¹ Laboratoire d'Astrophysique de Marseille, Aix-Marseille Univ., CNRS, CNES, 13013 Marseille, France</footnote>

<footnote n="12">¹² Institut d'Astrophysique de Paris, UMR7095 CNRS & Sorbonne Université, 98bis Bd Arago, 75014 Paris, France</footnote>

<footnote n="13">¹³ Instituto de Astronomía y Ciencias Planetarias, Universidad de Atacama, Copayapu 485, Copiapó, Chile</footnote>

<footnote n="14">¹⁴ Millennium Nucleus for Galaxies (MINGAL), Valparaíso, Chile</footnote>

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

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

<footnote n="17">¹⁷ HH Wills Physics Laboratory, University of Bristol, Tyndall Ave, Bristol BS8 1TL, UK</footnote>

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

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

<footnote n="20">²⁰ IFPU, Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy</footnote>

<footnote n="21">²¹ Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA</footnote>

<footnote n="22">²² INAF – Osservatorio Astronomico di Padova, Vicolo Osservatorio 5, I-35122 Padova, Italy</footnote>

<footnote n="23">²³ California Institute of Technology, Pasadena, CA 91125, USA</footnote>

<footnote>* Based in part on observations collected at the European Southern Observatory under ESO programmes 0110.A-4192 and 0111.1-0186.</footnote>

<footnote>** Corresponding author.</footnote>
# 1. Introduction

Galaxy clusters are the most massive virialised structures in the Universe and, as such, are excellent natural laboratories to study the formation, structure assembly, and evolution of mass halos and sub-halos across cosmic time. The combination of optical data, X-ray observations, and the signal from the Sunyaev–Zel'dovich (SZ) effect (Sunyaev & Zeldovich 1972) is ideal for studying galaxy clusters. A multi-wavelength analysis allows us to acquire a complete picture of the physics behind different components – dark matter (Zwicky 1933), gas (Cavaliere & Fusco-Femiano 1981), and galaxies (e.g. Abell 1958) – and see how they interact to configure these massive objects.

In this work we studied dynamical cluster mass profiles, the presence of substructures, and velocity anisotropies by analysing the projected-phase-space (PPS) distribution of member galaxies (e.g. Carlberg et al. 1997) under the assumption of dynamical equilibrium. As galaxies enter into clusters through hierarchical accretion, their orbits contain important information on the processes that lead to cluster mass assembly and internal dynamical relaxation (e.g. Biviano & Katgert 2004). Moreover, knowledge of the orbits of cluster galaxies is crucial to understanding environmental effects that differentiate galaxy evolution in clusters from those in the field (Lotz et al. 2019; Tonnesen 2019).

Our studies were carried out using the MG-MAMPOSST (Modelling Anisotropy and Mass Profile of Spherical Observed Systems) code from Pizzuti et al. (2023). This technique allows us to reconstruct cluster mass profiles and velocity anisotropy profiles for several parametrisations of the total gravitational potential, in general relativity and in general dark energy or modified gravity frameworks. In other words, MG-MAMPOSST jointly obtains the best-fit parameters for the cluster mass distribution and for the velocity anisotropy profile of cluster members (Munari et al. 2014; Zarattini et al. 2021); this analysis allows us to retrieve information about the orbits of galaxies by distinguishing, statistically, between populations with radial and circular configurations.

We developed our analysis using the Cluster HERitage project with XMM-Newton – Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) cluster sample (CHEX-MATE Collaboration 2021). This project is based on the three-mega-second Multi-Year Heritage Programme, whose goal is to obtain X-ray observations of a minimally biased, signal-to-noise-limited (S/N > 6.5) sample of 118 galaxy clusters detected by Planck through the SZ effect. The CHEX-MATE census contains two subsamples: Tier 1 (T1), which includes a census of the population at the most recent time ($0.05 < z < 0.2$), and Tier 2 (T2), the most massive objects to have formed at $z < 0.6$ (see Fig. 1 in CHEX-MATE Collaboration 2021). Therefore, the CHEX-MATE clusters offer us the opportunity to study the kinematical and dynamical properties of the most massive clusters of the Universe detected through their SZ signal.

The masses of galaxy clusters can be obtained directly under specific assumptions (see e.g. Pratt et al. 2019). Usually, X-ray mass estimates can be recovered either by assuming hydrostatic equilibrium (see e.g. Ettori et al. 2013) or through scaling relations among integrated quantities, such as gas mass and temperature (e.g. Kravtsov et al. 2006). The SZ effect can also be used to estimate masses through the Comptonisation parameter, $Y_{SZ}$.

We adopted a dynamical method that uses the line-of-sight (LOS) velocities of member galaxies and their projected positions with respect to the cluster centre to directly measure

masses. Different mass estimation techniques can be biased due to various effects, such as the violation of the assumption of hydrostatic or dynamical equilibrium, or asymmetric, complicated true geometries. To account for these deviations for all these methods, we studied the mass bias parameter $(1 - B_1) = M_{500}^{SZ}/M_{500}^{M}$ (Ruel et al. 2014; Penna-Lima et al. 2017; Amodeo et al. 2017) to characterise the cluster mass estimations when these are performed using the SZ–X-ray signature or other dynamical techniques. We refer to the SZ–X-ray-calibrated mass estimates as $M^{SZ}$, and the mass derived by the MG-MAMPOSST analysis as $M^{M}$.

Our main aims can be summarised in three points: (1) obtain the full kinematic mass profiles for the CHEX-MATE clusters, which allows future comparisons with those obtained from X-ray, SZ, and weak-lensing (WL) techniques; (2) study the scaling relations and the presence of possible biases in the mass estimates derived from different techniques; and (3) statistically study the main kinematical properties of the cluster members by disentangling their velocity anisotropy and orbits, which depend on the mass of the clusters.

Throughout this paper, we adopt a Λ cold dark matter cosmology with $\Omega_m = 0.3$, $\Omega_\Lambda = 0.7$, and $H_0 = 70 \text{ km s}^{-1} \text{ Mpc}^{-1}$. We define $r_\Delta = r_{\Delta c}$ as the radius of a spherical overdensity within which the mean density is $\Delta$ times the critical density of the Universe. The corresponding enclosed mass is $M_\Delta$ (e.g. $r_{200}$, $M_{200}$ for $\Delta = 200$, and so on). Logarithms in base 10 are defined as $\text{Log } X = \log_{10} X$ for any variable $X$, while natural logarithms are indicated by $\ln X$.

# 2. Dataset

Our analysis is based on the radial velocities of cluster members. However, not all the 118 clusters of the CHEX-MATE sample have spectroscopic data of sufficient quality to enable the analysis. So, we restricted our study to 75 clusters, which correspond to a subset of galaxy systems with at least 50 confirmed members within an aperture radius $R = r_{\text{ap,200}}$ – as estimated by Sereno et al. (2025) – with spectroscopic redshifts, which is a reasonable minimum number of members to obtain an accurate enough analysis of the mass and dynamical properties (see e.g. Ferragamo et al. 2020). As stated in Sereno et al. (2025), the cluster centre has been identified as the X-ray luminosity peak. As shown in Fig. 1, this subsample represents the 66% of the full CHEX-MATE clusters, with no bias in redshift or mass. In this figure the sample is divided into redshift bins (with width equal to 0.05) and mass bins (with width equal to $10^{14} M_\odot$), our analysis covers more than 50% of the CHEX-MATE clusters in each bin, with the exception of those with $z > 0.55$. In other words, the full CHEX-MATE sample is well represented in the present work. Most of the clusters present a large number of members with spectroscopic redshifts, with more than 100 confirmed members within $4.0 \text{ Mpc}$ in the majority of the cases, and even more than 300 for some of them. PSZ2 G006.49+50.56 (Abell 2029) and PSZ2 G044.20+48.66 (Abell 2142) represent outstanding cases with more than 1000 total spectroscopic redshifts available.

The cluster member selection for the sample used in this work was developed by Sereno et al. (2025) using the CLEAN method (fully described in Mamon et al. 2013; see Appendix B of that paper), which is based on galaxy location in the PPS. So, the radial velocities of cluster members are the same as used in Sereno et al. (2025), who based their analysis on the Sloan Digital Sky Survey (SDSS) DR18, the Dark Energy Spectroscopic Instrument (DESI) databases, and the NASA/IPAC
<figure><img src="image_3.png" /><figcaption>Fig. 1. Cluster redshift and $M_{500}^{SZ}$ mass distributions of the full CHEX-MATE sample (red) and the subsample (blue). $M_{500}^{SZ}$ has been derived iteratively from the $Y_{SZ}-M_{500}^{SZ}$ relation calibrated using hydrostatic masses obtained from XMM-Newton data (see CHEX-MATE Collaboration 2021).</figcaption></figure>

Extragalactic Database (NED), in addition to some private catalogues and proprietary data. Moreover, we added new spectroscopic observations<fnref n="1" /> for three CHEX-MATE clusters, PSZ2 G113.91–37.01, PSZ2G046.10+27.18, and PSZ2G083.29–31.03, all of which are in the redshift range $0.37 < z < 0.41$. We observed galaxies in these fields at the 10.4 m Gran Telescopio Canarias (GTC) telescope with the Optical Spectrograph and InfraRed Imager System (OSIRIS) in multi-object observing mode acquiring 2400 s exposures and two mask configurations for each cluster. For more details on the instrumental setup, data reduction process, and redshift estimates, see Ferragamo et al. (2020) and Aguado-Barahona et al. (2022). The procedures applied in this work are exactly the same as those used there. So, we obtained radial velocities for 80, 61, and 66 cluster members, respectively, with a precision of about 70 km s⁻<fnref n="1" />. These spectroscopic observations are limited, showing magnitude limits of $r' \sim 21.5$, and low sample completeness, below 30% within $r_{500}$.

For these three clusters, we selected members by applying a $2.7\sigma$ clipping in the phase-space to minimise the fraction of interlopers<fnref n="2" /> (Mamon et al. 2010). However, we find that only PSZ2G046.10+27.18 has more than $N = 50$ members within the estimated $r_{200}$. From here on, we refer to $r_{ap,200}$ as the radius provided by the CLEAN procedure, to distinguish it from the values obtained by the MG-MAMPOSST fit.

In summary, including the GTC observations, we created 77 catalogues of cluster member galaxies (one for each cluster), with a total of 17 521 radial velocities, i.e. a mean of about 220 members per cluster. Of these, 75 clusters have been found to have $N \ge 50$ confirmed members within the aperture radius $R = r_{ap,200}$, after the application of the CLEAN algorithm (Sereno et al. 2025). In Fig. E.1 we plot the number of cluster members identified within $R = r_{ap,200}$ after the cleaning as a function of the aperture velocity dispersion $\sigma_{ap,200}$ estimated by Sereno et al. (2025). Considering this spectroscopic dataset, we present in the following sections our procedure to analyse the main kinematical and dynamical properties of the CHEX-MATE cluster sample.

### 3. Kinematics analysis setup

We inferred the cluster mass profiles of the CHEX-MATE sample by employing the MG-MAMPOSST method (Pizzuti et al. 2023), which is based on the MAMPOSST (Mamon et al. 2013) algorithm to simultaneously recover the (total) gravitational

potential and the orbit anisotropy profile of clusters from kinematic analyses of their member galaxies.

MAMPOSST is based on the assumption that the system is spherically symmetric and that galaxies are collisionless tracers of the gravitational potential. In this case, the dynamics is determined by the so-called (spherical) Jeans equation:

$$ \frac{d(v\sigma_r^2)}{dr} + 2\beta(r) \frac{v\sigma_r^2}{r} = -v(r) \frac{d\Phi}{dr}, \quad (1) $$

where $r$ indicates the 3D radial distance from the cluster centre, $v(r)$ is the number density profile of galaxies, $\sigma_r^2$ is the velocity dispersion along the radial direction, and $\Phi$ is the total gravitational potential. The quantity $\beta \equiv 1 - (\sigma_\theta^2 + \sigma_\phi^2)/2\sigma_r^2$ is the velocity anisotropy, with $\sigma_\theta^2, \sigma_\phi^2$ the velocity dispersion components along the tangential and azimuthal directions, respectively. In spherical symmetry, we have $\sigma_\theta^2 = \sigma_\phi^2$.

MAMPOSST and MG-MAMPOSST work in the so-called PPS, i.e. $R, v_{los}$, where $R$ is the projected radial distance from the cluster centre – defined to be at the peak of the projected X-ray emission (Sereno et al. 2025) – and $v_{los}$ is the LOS velocity of each galaxy, computed in the rest frame of the cluster. Given as input the parametric expressions of the gravitational potential, the number density profile, and the velocity anisotropy profile, the code implements a maximum likelihood fit to determine the best set of parameters describing the PPS distribution of cluster members.

MG-MAMPOSST extends the gravitational potential models available in the original MAMPOSST, further including generalised models of the velocity anisotropy profile, which we explore here. MG-MAMPOSST also includes popular modified gravity and alternative dark energy scenarios (see Pizzuti et al. 2021 for a description of the physics), viable at the cosmological level, which we do not analyse in this paper. Moreover, a recent update of the code implemented a module to include data of the velocity dispersion of the brightest cluster galaxy (BCG) in the fit, when available. Inclusion of the BCG velocity dispersion would allow the method to provide tight constraints on the slope of the dark matter profiles in clusters (see Biviano et al. 2023). However, in this work, we did not use the BCG module.

Therefore, we only considered galaxies outside of a minimum projected radius $R_{min} = 0.05$ Mpc. The central region was excluded to avoid contamination by the BCG, which dominates cluster dynamics at very small radii. Furthermore, in its current implementation, MG-MAMPOSST assumes negligible streaming motions and a Gaussian 3D velocity distribution of the tracers.

We ran MG-MAMPOSST over our CHEX-MATE subsample of 75 clusters, which includes members up to a maximum projected radius $R_{max} = 1.1 r_{200}^D$. Here $r_{200}^D$ is the estimate of the (dynamical) virial radius, either from the analysis of Sereno et al. (2025, $r_{ap,200}$) or from the velocity dispersion (see Sect. 5 in Barrena et al. 2024) for those clusters for which the former is not available. This definition of maximum radius for cluster membership is a common choice to ensure the validity of the Jeans equation, as $r_{200}$ is close to the virial radius which determines the region where the assumptions of Eq. (1) are satisfied.

We sampled the parameter space with a Markov chain Monte Carlo (MCMC) Metropolis-Hastings algorithm for 110 000 points, considering the first 10 000 as the burn-in phase. Our final chain contains 100 000 samples for each run; this number has been chosen has a reasonable value to save computational time while ensuring convergence. In more detail, we checked convergence by considering the four clusters with the fewest members

<footnote n="1">¹ These new dataset will be made publicly available soon, but it is already available upon request. An extract can be found in Appendix D.</footnote>

<footnote n="2">² Galaxies not belonging to the clusters falling inside the PPS.</footnote>
$(N < 60)$ within $r_{ap,200}$; and running $n = 5$ test chains for each of them. We then computed the corresponding Gelman-Rubin diagnostic coefficients $\hat{R}$ Gelman & Rubin (1992)<fnref n="3" /> checking that the requirement $\hat{R} \le 1.1$ is always satisfied, as done in Pizzuti et al. (2024).

For the main analysis, a total of eight MG-MAMPOSST chains were generated per cluster in the sample. We explored combinations of two models for the number density distribution and four models for the total mass profile. The first, the Navarro-Frenk-White (NFW), model can be written as

$$M_{\text{NFW}}(r) = M_{200} \frac{\ln(1 + r/r_s) - r/r_s/(1 + r/r_s)}{\ln(1 + c_{200}) - c_{200}/(1 + c_{200})}, \quad (2)$$

where $M_{200} = 200(H^2(z)/2G)r_{200}^3$ is the enclosed mass at $r_{200}$ and $c_{200} = r_{200}/r_{-2}$ is the concentration, with $r_{-2}$ being the radius at which the logarithmic slope of the density profile is -2. For the NFW model, this radius corresponds to $r_s$.

The Hernquist profile has a steeper decline of the density at large radii, and the mass enclosed in a radius $r$ reads (Hernquist 1990)

$$M_{\text{Her}}(r) = M_{200} \frac{r_s^2 + r_{200}^2}{r_{200}^2} \frac{r^2}{r + r_{200}^2}, \quad (3)$$

where $r_s = 2r_{-2}$.

The Burkert model (Burkert 1995) has a central core in the density distribution, and the mass is given by

$$M_{\text{Bur}}(r) = M_{200} \frac{\ln((r/r_s)^2 + 1) + 2\ln(r/r_s + 1) - 2\tan^{-1}(r/r_s)}{\ln((r_{200}/r_s)^2 + 1) + 2\ln(r_{200}/r_s + 1) - 2\tan^{-1}(r_{200}/r_s)}, \quad (4)$$

with $r_s \approx r_{-2}/1.5$.

Finally, the Einasto model Einasto (1965) can be written as

$$M_{\text{Eis}}(r) = M_{200} \frac{\gamma \left[ \frac{3}{m}, \frac{2}{m} \left( \frac{r}{r_v} \right)^m \right]}{\gamma \left[ \frac{3}{m}, \frac{2}{m} (c_{200})^m \right]}, \quad (5)$$

where we set the shape parameter<fnref n="4" /> $m = 3$ and $r_s = r_{-2}$ as for the NFW profile. Therefore, all models are described by two free parameters, the mass $M_{200}$ (or equivalently, the radius $r_{200}$) and the scale radius $r_s$, which will be fitted within the MG-MAMPOSST algorithm.

As mentioned above, we further considered two models for the number density profile: the projected NFW (pNFW) and the projected Hernquist (pHer). These models are defined as the integral of Eqs. (2) and (3) along the LOS. While the profiles are characterised by two free parameters, only the scale radius $r_v$<fnref n="5" /> is relevant in the analysis, as the other factors out in the solution of the Jeans equation implemented in MG-MAMPOSST. Note that for $r_v$, we first performed a preliminary fit to the member galaxy distribution in the phase space, which does not require the binning of data (Sarazin 1980), assuming a constant completeness of the sample in the projected radial distance explored.

This last assumption is rather strong, since our dataset exhibits incomplete and inhomogeneous spatial distributions; however, we carried out a detailed analysis on a few clusters to check that reasonable variation of the completeness does not provide relevant changes in the reconstructed mass profile. We then used the 95% confidence level (CL) limits of the inferred $r_v$ as the upper and lower bounds for the flat prior range in the MG-MAMPOSST fit (see below). We further verified that varying this interval produces negligible effects in the reconstructed mass profile.

As for the velocity anisotropy profile, we adopted a generalised Tiret (gT) model (see Mamon et al. 2019), which is an extension of the Tiret profile (Tiret et al. 2007) and is able to capture a broad ranges of possibilities for the orbits of member galaxies in clusters:

$$\beta_{gT}(r) = \beta_0 + \beta_{\infty}r/(r + r_{\beta}), \quad (6)$$

where $\beta_0$ and $\beta_{\infty}$ are the inner and outer values of the anisotropy, respectively, and $r_{\beta}$ is a characteristic scale radius. In all cases, we set $r_{\beta} = r_{-2}$ (i.e. the scale radius of the mass profile) as suggested by numerical simulations (Mamon et al. 2010). In MG-MAMPOSST, we worked with the rescaled parameters $\mathcal{A}_{0,\infty} = (1 - \beta_{0,\infty})^{-1/2}$, which are equal to one for completely isotropic orbits.

We adopted flat priors on all the fitted parameters, $r_{200}/\text{Mpc} \in [0.5, 5.0]$, $r_s/\text{Mpc} \in [0.04, 4.0]$, $\mathcal{A}_{0,\infty} \in [0.4, 7.0]$, $r_v \in [r_v^{\text{low}}, r_v^{\text{up}}]$, where the last range is between the 95% lower and upper limits from the external fit of the projected number density distribution. Except for $r_v$, the priors are non-informative.

For each posterior, we applied the Bayesian information criterion (BIC; see e.g. Mamon et al. 2019) to quickly select the best model among the eight combinations analysed:

$$\text{BIC} = k \ln N - 2 \ln \mathcal{L}, \quad (7)$$

where $k$ is the number of free parameters, $N$ the number of data points, and $\mathcal{L}$ the likelihood at the best fit. Generally, a model is said to be strongly preferred when its BIC is lower by a factor of ~6 with respect to the others. For the models with the same number of free parameters, the BIC comparison is equivalent to the chi-square $\chi^2$ comparison. It is important to stress that MG-MAMPOSST only indicates which model is favoured relative to the others, but not the best model in an absolute sense. We further quantified the systematics induced by the modelling as half of the maximum variation between all the best-fit values of each parameter.

## 4. Results

Here we present the main outcomes from the MG-MAMPOSST analysis for the CHEX-MATE sample of 75 clusters. Among the mass profile models assumed, all of them adequately fit the PPS distribution, with the NFW profile performing better for 45% of the clusters, followed by Burkert (32.5%), Einasto (12.5%), and Hernquist (8%). This distribution may indicate a slight preference towards steeper central profiles with respect to those with cores. While it is true that the presence of interlopers that survived the CLEAN procedure can produce a steeper profile, thus biasing in favour of NFW-like models, the statistical significance of this is generally negligible, with a very small $\Delta$BIC ($\le 1.0$ on average). In contrast, the fits to the number density distribution of galaxies notably prefer the projected Hernquist model compared to the NFW model for more than 70% of the analysed clusters, with several cases where $\Delta$BIC $\ge 6.0$. This preference

<footnote n="3">³ The Gelman-Rubin coefficient is defined as $\hat{R} = [W(L-1)/L + B/L]/W$, where $L$ is the length of the chains, $B$ the variance of the means of the chains, and $W$ the averaged variances of the individual chains across all chains.</footnote>

<footnote n="4">⁴ We checked that varying the exponent between $m = 2$ and $m = 5$ does not lead to statistically relevant changes in the final posterior.</footnote>

<footnote n="5">⁵ We defined the scale radius of the number density profile as $r_v$ to avoid confusion with that of the mass profile.</footnote>
seems to be in contrast to other studies (e.g. van der Burg et al. 2016), which indicate a preference of the pNFW model. However, one should note here that the assumption of a constant completeness may play a significant role. Indeed, while it is true that the constraints on a mass profile distribution are not strongly affected by the choice of the model of the number density profile in the MG-MAMPOSST fit, the slope of $v(r)$ and the value or the scale radius $r_v$ depends on the completeness. In order to deal with this issue, one should compute the completeness, cluster by cluster, by comparing photometric and spectroscopic catalogues. However, this analysis is beyond the scope of the paper. As an example, we considered the case of the cluster PSZ2 G049.22+30.8, for which the completeness as a function of the projected radius has been estimated according to the model

$$C(R) = A \exp(BR), \quad (8)$$

with $A = 0.817$, $B = -0.413$. While in this case the best-fit projected number density is always a NFW profile, the value of the scale radius changes from $r_v = 0.31^{+0.18}_{-0.10}$ when accounting for the correct completeness, to $r_v = 0.59^{+0.10}_{-0.13}$, when a constant completeness is assumed (uncertainties are given at the 95% CL).

In Table A.1 we report the regression results. In the second and third columns, we list the best-fit models (i.e. the models that give the highest posterior) for mass and density, respectively. For each entry, the first uncertainties refer to the (statistical) 95% confidence intervals, while the second values are the systematic uncertainties computed as discussed above. The last four columns show the value of the log posterior $\ln P$ of the best fit, the cluster redshift $z$, the value of the Anderson-Darling coefficient $A^2$ (see Sect. 4.3), and the estimated LOS velocity dispersion at $r_{200}$, $\sigma_{ap,200}$. Note that the MG-MAMPOSST procedure is robust in estimating $r_{200}$ against the choice of the mass, anisotropy and number density models, as shown by the small systematic uncertainties compared with the statistical errors.

In Fig. 2 we show an example of the mass and anisotropy profiles obtained from MG-MAMPOSST for one rich cluster ($N_{gal} = 291$ spectroscopic members within $r_{ap,200}$) in the sample, PSZ2 G067.17+67.46. The mass is parametrised as an NFW model, which gives the best-fit profile, with a pNFW model for the number density of the member galaxies.

We then computed the mass inside the virial radius as $M_{200} = (100 H^2(z)/G) \times r_{200}^3$, and we performed a first comparison with the results of Sereno et al. (2025). Figure 3 displays the values of Log $M_{200}$ as obtained by MG-MAMPOSST (vertical axis), and Log $M_{\sigma_{ap,200}}$ estimated from the aperture velocity dispersion $\sigma_{ap,200}$ with the method of Sereno et al. (2025, horizontal axis). The superscript (M) refers to the MG-MAMPOSST estimates.

The results show an overall good agreement among the two estimates, confirming the robustness of the simpler analytical approach of Sereno et al. (2025). However, for seven clusters MG-MAMPOSST tends to overestimate $M_{200}$ with respect to the findings of Sereno et al. (2025). For those clusters we find an average logarithmic difference Log $M_{200}^M$ - Log $M_{\sigma_{ap,200}}$ > 0.3 dex. If we look at the relation between $M_{200}^M$ and the aperture velocity dispersion $\sigma_{ap,200}$ (Fig. 4) obtained by Sereno et al. (2025) using the CLEAN method (Mamon et al. 2013), it is easy to see that six of the seven clusters are characterised by large values of $\sigma_{ap,200} > 10^3$ km/s. This pattern suggests that clusters in the CHEX-MATE sample with very large velocity dispersions suffer a stronger contamination of interlopers. This fact was first found by Mamon et al. (2010) who, working with simulated halos, found that the fraction of interlopers should

<figure><img src="image_4.png" /><figcaption>Fig. 2. Mass profile (top) and radial velocity anisotropy profile (bottom) of the cluster PSZ2 G067.17+67.46, as obtained via the MG-MAMPOSST analysis with $N_{gal} = 299$ galaxies and adopting a NFW model and a gT model, respectively. In both plots, darker and lighter shaded regions refer to 68% and 95% CLs, respectively. The vertical dashed line in the top panel identifies the best fit for $r_{200} = 2.19$ Mpc.</figcaption></figure>

be (slightly) more important in the more massive halos since they have (slightly) lower concentrations (Navarro et al. 1997; Macciò et al. 2008) and they typically live in richer environments. That is, the larger fraction of interlopers that survived the membership selection performed by Sereno et al. (2025) in high-mass clusters could translate into a higher mass estimation by MG-MAMPOSST – which is a more complex model with respect to that of Sereno et al. (2025) – for the most massive systems.

Moreover, in MG-MAMPOSST the concentration $c_{200} = r_{200}/r_{-2}$ is a free parameter (or better, it is a quantity derived from two parameters optimised within the fit), while in Sereno et al. (2025) it is fixed as a function of $r_{200}$ using the concentration-mass relation proposed by Diemer & Joyce (2019). As shown by the blue points in Fig. 8, the off-trend clusters are characterised by a concentration estimated by MG-MAMPOSST that is lower than the prediction of Diemer & Joyce (2019, solid brown line), consequently boosting the reconstructed $M_{200}$.

The blue line in Fig. 4 shows the best orthogonal distance regression (ODR) fit for the linear model

$$\mathrm{Log}\left(F_z \frac{M_{200}^M}{10^{14} M_{\odot}}\right) = \alpha_{\sigma} \mathrm{Log}\left(\frac{\sigma_{ap,200}}{10^3 \mathrm{km s}^{-1}}\right) + \beta_{\sigma} + \gamma_{\sigma} \mathrm{Log}F_z, \quad (9)$$

when fixing the redshift to the pivot value of the sample $z_p = 0.2$; in Eq. (9), $F_z = H(z)/H(z_p)$. We find

$$\alpha_{\sigma} = 3.61 \pm 0.17, \quad \beta_{\sigma} = 1.07 \pm 0.01, \quad \gamma_{\sigma} = -0.83 \pm 0.70, \quad (10)$$
<figure><img src="image_5.png" /><figcaption>Fig. 3. Scatter plot of Log $M_{200}$, inferred from the kinematic analysis with MG-MAMPOSST (vertical axis), and Log $M_{\sigma,200c}$ obtained by Sereno et al. (2025, horizontal axis). The error bars represent the 68% CL, and the dashed black line the bisector. Darker points highlight those clusters for which Log ($M_{200c}^M/M_{\sigma,200c}$) > 0.3.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 4. Logarithm of the aperture velocity dispersion within R = $r_{200}$, $\sigma_{ap,200}$ compared to the logarithm of $M_{200}^M$ estimated from MG-MAMPOSST. The bars refer to 68% uncertainties. The blue line indicates the best-fit linear model of Eq. (9) when z = 0.2. The solid black line identifies the $\sigma \propto M^3$ relation.</figcaption></figure>

which is slightly steeper than $\alpha \sim 3$, the value reported in, for example, Munari et al. (2013) and Ferragamo et al. (2020). If we consider the T1 and T2 clusters separately, the ODR fit of Eq. (9) provides a slightly lower $\alpha^{T1} = 3.31 \pm 0.17$ for the T1 sample – close to the ~ $\sigma^3$ relation and a higher value $\alpha^{T2} = 4.81 \pm 0.64$ for T2. This reflects the fact that T2 contains massive clusters by construction, and that among the seven off-trend objects, five belong to T2. The best ODR fits for $z_p$ are shown as the dashed red and green lines for T1 and T2, respectively.

## 4.1. Comparison with SZ mass estimates

To compare our kinematic masses with the ones obtained from the SZ effect, we further derived the distribution of $M_{500} = 500 H^2(z)/(2G) \times r_{500}^3$, as follows: for each value of $r_{200}$, $r_s$ in the MCMC chain we obtained the corresponding parametric mass

profile $M(r|r_{200}, r_s)$. The radius $r_{500}$ is then given by the numerical solution of the equation

$$ \frac{M(x)}{x^3} = \frac{500 H^2(z)}{2G} \qquad (11) $$

From that, the computation of $M_{500}$ is straightforward.

In Fig. 5 we plot the values of Log $M_{500}$ as obtained by MG-MAMPOSST (vertical axis) and inferred by SZ (horizontal axis), where the superscript (M) and (SZ) refers to the MG-MAMPOSST and SZ estimates, respectively. We find a strong correlation between the two estimates (Pearson coefficient $r = 0.68$, p-value < 10⁻⁵); we fitted the linear relation

$$ \log\left(\frac{M_{500}^{\text{SZ}}}{10^{14} \text{ M}_{\odot}}\right) = \alpha_{SZ} \log\left(\frac{M_{500}^{M}}{10^{14} \text{ M}_{\odot}}\right) + \log(1 - B), \qquad (12) $$

where the mass-dependent bias (1 − B) is defined similarly to that done, for example, in Ferragamo et al. (2021) and Aguado-Barahona et al. (2022). This way, the usual bias (1 − $B_1$) = $M_{500}^{\text{SZ}}/M_{500}^M$ is recovered when $\alpha_{SZ} = 1$.

We obtain $\alpha_{SZ} = 0.84 \pm 0.05$, and an overall mass-dependent bias Log(1 − B) = −0.18 ± 0.06 dex when all the galaxy clusters in the sample are considered, corresponding to (1 − B) = 0.67 ± 0.09; this means that the kinematic results from MG-MAMPOSST produce larger masses with respect to those derived from the Planck SZ signal. This value seems to be in agreement within errors with those obtained by Planck Collaboration VI (2020) and Lesci et al. (2023) from Baryon acoustic oscillation and 3D clustering of the PSZ2 clusters. In addition, our estimate is similar to that found by the fit of Sereno et al. (2025), although slightly larger than that of Ferragamo et al. (2021) and Aguado-Barahona et al. (2022), who obtained a bias of (1 − B) ~ 0.80 comparing Planck SZ-X-ray calibrated masses of 270 PSZ1 and 388 PSZ2 clusters, respectively, with the kinematic mass estimated through the total velocity dispersion. When forcing $\alpha_{SZ}$ to unity, however, the difference between the two mass estimates is even larger, almost doubled: (1 − $B_1$) = 0.43 ± 0.08, which is not consistent with the abovementioned works.

It is interesting to note that the sample includes a few clusters for which the kinematic mass is about one order of magnitude larger than the one given by the SZ proxy. We mark the four clusters for which

$$ \delta\eta = \log M_{500}^{M} - \log M_{500}^{\text{SZ}} > 0.8, \qquad (13) $$

as red points in Fig. 5. However, the PPS distributions of those clusters show evident sign of a disturbed dynamical state, as displayed in the left plots of Fig. B.1, with very large velocity dispersions and presence of possible important structures that are not part of the main cluster, as in the case of PSZ2G056.77+36.32, and PSZ2G040.03+74.95.

The disturbed dynamical state of these four clusters is further indicated by the distributions of the velocities along the LOS (central plots of Fig. B.1); those for PSZ2G056.77+36.32 and PSZ2G040.03+74.95 exhibit noticeable bimodal distributions. This is further indicated by relatively large values of the Anderson Darling coefficient, $A^2$ – a measure of the degree of unrelaxation that is discussed in Sect. 4.3 – for three of the four outliers. Note that here we considered the members selected by running the CLEAN algorithm in its default settings, which do not employ an identification of substructures. However, as already done in Sereno et al. (2025), a refined analysis can be performed to study these systems. In more detail, PSZ2G040.03+74.95 (also known as Abell 1831) has been already identified as a double system: it is made of two clusters, A and B at z = 0.063 and
<figure><img src="image_7.png" /></figure>

**Fig. 5.** Scatter plot of $M_{500}$ inferred from the SZ effect by Planck (on the vertical axis) and from the kinematic analysis with MG-MAMPOSST (horizontal axis). The error bars represent the 68% CL, and the red points show those clusters for which the logarithmic difference ($\delta\eta$) is larger than 0.8 dex. Dark green points highlight clusters belonging to T2. The purple line is the ODR fit that considers all the clusters, while the blue one is the linear fit done by discarding the four outliers. The best-fit (with the outliers removed) linear models for T1 and T2 are shown separately with the dashed red and dark green lines, respectively.

0.076, respectively, nearly aligned along the LOS, as discussed in Sereno et al. (2025). A similar situation holds for the object PSZ2G218.81+35.51 (see Sereno et al. (2025) and also the histogram of the LOS velocities in Fig. 6). We find a logarithmic difference $\delta\eta = 0.68 \pm 0.42$ for PSZ2G218.81+35.51, where the large uncertainties reflect the inefficiency of MG-MAMPOSST in fitting the double system with a single mass model.

By considering the two sub-groups in PSZ2G040.03+74.95 and PSZ2G218.81+35.51 as independent clusters<fnref n="6" />, the kinematical analysis gives $M_{200} < 7 \times 10^{14} M_\odot$ and $M_{200} = 53_{-26}^{+9} \times 10^{14} M_\odot$ for PSZ2G040.03+74.95 A, B, and $M_{200} = 2.4_{-2.2}^{+1.8} M_\odot$, $M_{200} = 6.1_{-2.7}^{+3.6} M_\odot$ for PSZ2G218.81+35.51 A, B. Note that for PSZ2G040.03+74.95 A we are only able to set an upper limit with MG-MAMPOSST, given the small number of members considered in the fit (only 25 members lie within $R = r_{200}$, well below the limit of 50 we imposed on the full sample) and the relatively small $\sigma_{ap,200} = 461 \text{ km s}^{-1}$, compared to the error in the LOS velocities. It is important to mention that the results for these bi-modal systems are in agreement (albeit with larger uncertainties) with the values obtained by the simple halo model of Sereno et al. (2025).

For the most massive sub-components in the bimodal clusters (which are identified as the CHEX-MATE clusters: see Sereno et al. (2025) and references therein), the logarithmic difference is now reduced to $\delta\eta = 0.8 \pm 0.3$ (PSZ2G040.03+74.95) and $\delta\eta = 0.15 \pm 0.11$ (PSZ2G218.81+35.51), instead of 1.4 and 0.68. These interesting cases, as well the other three clusters mentioned above, deserve a separate detailed analysis of their properties. Such an analysis is beyond the scope of this paper and it will be performed in a separate work. We consider only the main sub-components of the double systems in the forthcoming analyses.

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

**Fig. 6.** LOS velocity distribution for the system PSZ2G218.81+35.51, made of two separate clusters, as mentioned in Sereno et al. (2025).

If we remove the four clusters with large $\delta\eta$ and fit Eq. (12) again (dashed blue line in Fig. 5), we remarkably find a negligible mass-dependent bias within the uncertainties $\text{Log}(1-B) = 0.02 \pm 0.04$, which corresponds to $1-B = 1.04 \pm 0.09$, with a slightly flatter slope $\alpha_{SZ} = 0.71 \pm 0.04$. Note that, however, the bias is still large when fixing the slope to 1, $1-B_1 = 0.54 \pm 0.11$, from which we can conclude that kinematic mass estimates are a factor of two larger, on average, with respect to SZ probes. Note that this value is now in agreement with the findings of, for example, Lesci et al. (2023). Here we have ignored the possible redshift dependence of the sample and neglected the Malmquist bias, both were already considered in Sereno et al. (2025). SZ masses are biased high, which can produce an underestimation of B. When considering T1 and T2 separately (after removing the four biased clusters), the fit is still consistent with a negligible mass-dependent bias: $(1-B) = 1.12 \pm 0.13$ for the former and $(1-B) = 1.32 \pm 0.44$ for the latter. The larger uncertainties in T2 reflect a wide scatter in $\delta\eta$, as shown by the green points in Fig. 5. For $B_1$, we find $(1-B_1) = 0.64 \pm 0.04$ for T1 and $(1-B_1) = 0.40 \pm 0.04$ for T2.

Note that $\delta\eta$ does not exhibit a strong correlation with redshift (see Fig. 7), with a scattered distribution. However, if we compute the fraction of galaxies with $\delta\eta > 0.5$ in two redshift bins, $z < 0.2$ and $z > 0.2$, we find 0.14 and 0.28, respectively. This difference should not come as a surprise, since clusters with $z > 0.2$ belong to the T2 selection, which contains massive clusters by construction (thus characterised by larger velocity dispersions; see e.g. Sereno et al. 2025). Those clusters exhibit overall larger discrepancies between $M_{500}^{SZ}$ and $M_{500}^M$, attributed to a potentially stronger contribution of interlopers in the kinematics analysis, as already mentioned above, in clusters with high velocity dispersions.

## 4.2. Concentration-mass relation

In Fig. 8 we plot the relation between $c_{200}$ and $M_{200}$ of the sample, as derived by the MG-MAMPOSST analysis. For the two double systems, PSZ2G056.77+36.32 and PSZ2G218.81+35.51, we considered the sub-components separately and ended up discarding PSZ2G056.77+36.32 A, given its limited statistics. Although the points appear quite scattered and have large uncertainties, a mild, positive correlation can be observed (Spearman rank coefficient 0.33, corresponding to a p-value of $2 \times 10^{-3}$).

<footnote n="6">⁶ The selection of members for the subsystems is performed as discussed in Sereno et al. (2025).</footnote>
<figure><img src="image_9.png" /><figcaption>Fig. 7. $\delta\eta \equiv \text{Log}(M_{200}^M/M_{500}^{SZ})$ as a function of the redshift. The error bars indicate the 68% statistical uncertainties. Points are colour-coded according to the values of the $A^2$ coefficient. The red vertical line indicates the pivot redshift $z = 0.2$. Note that we have included the four outliers ($\delta\eta > 0.8$) as well as the full double systems PSZ2G040.03+74.95 and PSZ2G218.81+35.51 to highlight the values of $A^2$.</figcaption></figure>

For a quantitative analysis, we performed a Bayesian regression. We followed the CoMaLit (Comparing Masses from Literature) scheme described in Sereno & Ettori (2015b,a, 2017) and Sereno et al. (2015a), which we refer to for details, and implemented in the R-package LIRA (Sereno 2016)<fnref n="7" />.

The mass-concentration relation is modelled with a linear relation in log-space, similarly to Eq. (9),

$$ \text{Log}(c_{200}) = \alpha_C + \beta_C \text{Log}(M_{200}^M/M_{\text{pivot}}) + \gamma_C \text{Log} F_z, \quad (14) $$

where $F_z = H(z)/H(z_{\text{ref}})$ is the redshift-dependent Hubble parameter normalised to the reference redshift of the sample, $z_{\text{ref}} = 0.2$, and $M_{\text{pivot}} = 10^{15} M_{\odot}$ is the pivot mass. The concentration at the reference mass and redshift is given by $c_{200,\text{ref}} = 10^{\alpha_C}$.

Statistical biases and correlation effects have to be carefully dealt with (Sereno et al. 2015b). For our analysis, we assumed that both $X = \text{Log}(M_{200c}^M/M_{\text{pivot}})$, that is, the covariate in Eq. (14), and $Y = \text{Log}(c_{200})$, that is, the response, are intrinsically scattered with respect to the true mass, $Z$, whose distribution we modelled as a Gaussian with a redshift-dependent mean and variance (Sereno & Ettori 2015a). We assumed that the covariate is scattered but unbiased with respect to the true mass. We considered heteroscedastic and correlated uncertainties for mass and concentration. For central location and scale, we considered the bi-weight estimators (Beers et al. 1990). We applied non-informative priors (Sereno & Ettori 2015a, 2017): the uniform distribution for $\alpha_C$ and the mean of the distribution of $Z$; the Student's $t_1$ distribution with one degree of freedom for the slopes $\beta_C$ and $\gamma_C$; and the Gamma distribution for the inverse of either variances or intrinsic scatters.

The regression favours a flattening and upturn of the relation at large masses, in agreement with theoretical predictions (Diemer & Joyce 2019). The best-fit values of the linear relation are

$$ \alpha_C = 0.49 \pm 0.07, \quad \beta_C = 0.20 \pm 0.31, \quad \gamma_C = -0.23 \pm 1.6. \quad (15) $$

To test the robustness of the results on the regression scheme, we also fitted the relation with a simple linear model, similarly to

<figure><img src="image_10.png" /><figcaption>Fig. 8. Concentration–mass relation for the CHEX-MATE sample. The error bars indicate 68% uncertainties in the parameters; we highlight in light blue the clusters with $\text{Log}(M_{200}^M/M_{\text{cr},200c}) > 0.3$. The dashed blue line is the best-fit Bayesian linear regression derived using the Sereno & Ettori (2015b) method; the shaded region represents the corresponding $1\sigma$ contour. Purple and red crosses indicate the values obtained by stacking clusters in bin of SZ mass and $\sigma_{\text{ap},200}$, as explained in Sect. 5. The solid brown line shows the theoretical prediction of Diemer & Joyce (2019).</figcaption></figure>

that done, for instance, in Biviano et al. (2017), using the ODR method and neglecting the redshift dependence,

$$ \text{Log} c_{200} = A_1 \text{Log} \left( \frac{M_{200}^M}{10^{14} M_{\odot}} \right) + A_2. \quad (16) $$

We find $A_1 = 0.27 \pm 0.08$, $A_2 = 0.20 \pm 0.11$, which confirms the trend of rising concentration for high mass clusters. The same is valid if we fit the T1 and T2 clusters separately, in particular $A_1^{T1} = 0.35 \pm 0.12$ and $A_1^{T2} = 0.21 \pm 0.14$. Our findings might appear to be in tension with other theoretical expectations (e.g. Dutton & Macciò 2014; Ragagnin et al. 2020) and observational determinations of the $c-M$ relation (see e.g. Merten et al. 2015 and Biviano et al. 2017; the latter study is also based on the MG-MAMPOSST technique). However, this statement should be considered cautiously. The CHEX-MATE sample covers the very massive end of the halo mass function at intermediate redshifts, where the upturn is visible. Selection effects as well as the uncertainties in the kinematic analysis can play a significant role; in particular, X-ray- and SZ-selected samples tend to favour more centrally concentrated and relaxed clusters, biasing the inferred $c-M$ relation upwards (e.g. Meneghetti et al. 2014). Moreover, as mentioned before, the presence of residual interlopers may affect the high-mass end of the sample, giving rise to the observed trend.

### 4.3. Substructures and un-relaxation

To further investigate the overall kinematical state of the cluster sample, we computed the Anderson-Darling coefficient $A^2$ (Anderson & Darling 1952), which measures deviation from Gaussianity of the LOS velocity distribution, which is a good indicator of the dynamical state of a cluster (e.g. Roberts et al. 2018). Large values of $A^2$ have been shown to be associated with deviations from relaxation (e.g. Pizzuti et al. 2020; Barrena et al.

<footnote n="7">⁷ The package LIRA (Linear Regression in Astronomy) is publicly available from the Comprehensive R Archive Network at https://cran.r-project.org/web/packages/lira/index.html.</footnote>
<figure><img src="image_11.png" /><figcaption>**Fig. 9.** $\delta\eta = \text{Log } M_{500}^M - \text{Log } M_{500}^{SZ}$ as a function of the fraction of galaxies in substructures, $f_{\text{sub}}$, computed using the DS+ code from Benavides et al. (2023). Colours and error bars are as in Fig. 7.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Fraction of galaxies in substructures, $f_{\text{sub}}$, as a function of the logarithm of the Anderson-Darling coefficient, $A^2$. The darker stars indicate the clusters for which $\text{Log}(M_{200}^M / M_{\sigma,200c}) > 0.3$.</figcaption></figure>

2024). The colours of each point in Fig. 7 indicates the corresponding value of the (logarithm of the) Anderson-Darling coefficient, computed using all galaxies within $R = r_{\text{ap,200c}}$ as given by CLEAN (Sereno et al. 2025). Note that three of the four outliers with $\delta\eta > 0.8$ are characterised by large values of $A^2 = 0.64, 0.87$, and 4.63 for PSZ2G056.77+36.32, PSZ2G008.31-64.74, and PSZ2G040.03+74.95. Note that for PSZ2G056.77+36.32 the group below $v_z = -2000$ km/s is at large radii from cluster centre (see the top left plot of Fig. B.1), and thus is not considered in the computation of the Anderson-Darling coefficient. If all galaxies in the phase space are included, we find $A^2 = 2.49$ While PSZ2G075.71+13.51 exhibits a low value $A^2 = 0.21$, the velocity distribution is negatively skewed ($s \sim -0.11$) with a further negative excess of kurtosis ($k-3 = -0.43$).

No strong correlations can be highlighted among $A^2$ and $\delta\eta$, but dividing again the sample into two subsamples according to the median value $A^2 = 0.44$, the average logarithmic differences in the two bins, computed by 1000 bootstrap resampling, are $\langle|\delta\eta|\rangle = 0.23 \pm 0.08$ and $0.13 \pm 0.07$ for $A^2 > 0.44$ and $A^2 < 0.44$, respectively. This small difference in mean $|\delta\eta|$ shows a weak preference of lower biases for clusters with a more regular LOS velocity dispersion.

The MG-MAMPOSST method may be sensitive to the presence of relevant substructures in clusters, which are connected to recent merging activity and thus indicate a possible disturbed dynamical state (e.g. Wen & Han 2013; Parekh et al. 2015; Kimmig et al. 2023). Thus, we searched for possible correlations between the number of substructures in the cluster and the observed logarithmic difference, $\delta\eta$. We estimated the fraction of galaxies in substructures, $f_{\text{sub}}$, over the total number of galaxies within $R = r_{200}$ by means of the DS+ method<fnref n="8" /> from Benavides et al. (2023). This DS+ method extends the original algorithm of Dressler & Shectman (1988) for the identification of sub-groups in clusters by looking for galaxy clumps whose mean velocities and dispersion deviate from the global cluster values. DS+ employs a combination of projected position and LOS velocities to identify substructures and estimate the probability of each galaxy belonging to each of those groups.

For each cluster, we ran 1000 Monte Carlo simulations and we a *p*-value of 0.01 as the threshold for the identification of a substructure (see also Barrena et al. 2024). We worked in the so-called non-overlapping mode, i.e. each galaxy is attributed to at most one group. We then defined $f_{\text{sub}} = N_{\text{sub}}/N_{200}$, where $N_{\text{sub}}$ is the total number of galaxies in all the substructures found by DS+, and $N_{200}$ is the number of members in the projected cylinder $R < r_{200}$. In Fig. 9 the distribution of $\delta\eta$ as a function of $f_{\text{sub}}$ is shown, colour-coded according to the values of $A^2$. Again, no relevant trends can be drawn from the sample, with a more scattered distribution across the range of $f_{\text{sub}}$. It is worth pointing out that a $\delta\eta$ close to zero does not indicate that the cluster mass estimations are unbiased, but it may suggest that the SZ and kinematic masses are biased towards the same direction. Figure 10 further shows the relation between $f_{\text{sub}}$ and $\text{Log } A^2$. It is evident that, despite the correlation not being strong (Pearson coefficient $r = 0.43$), a trend of having large $A^2$ for large values of $f_{\text{sub}}$ is present.

A more accurate analysis would involve masses estimated from WL, as well as the combination of different morphological parameters (e.g. those derived from X-ray analyses; see Campitiello et al. 2022) to have a clearer indication of the equilibrium state of the clusters. This will be the subject of follow-up work.

## 5. Orbits of galaxies as a function of mass and redshift

The orbital anisotropy of member galaxies in clusters is one of the key quantities in kinematic analyses. It has been shown that the nature of anisotropy can be linked to the dynamical state of the cluster and its formation history (e.g. Mamon et al. 2019); in general, studies of observed clusters and cluster-size halos in cosmological simulations have identified an overall trend of isotropic orbits in the cluster core, and an increasing radial anisotropy towards the outskirts. Moreover, galaxies of different morphological types in clusters exhibit different orbital behaviours. In particular, early-type galaxies that are believed to have fallen into the cluster at early times may prefer more isotropic orbits, while spiral galaxies are thought to follow fairly radial paths (e.g. Munari et al. 2014; Mamon et al. 2019; Biviano et al. 2024). The reconstruction of anisotropy profiles in clusters is generally

<footnote n="8">⁸ The code is freely available at https://github.com/josegit88/MilaDS.</footnote>
<figure><img src="image_13.png" /><figcaption>**Fig. 11.** Radial velocity anisotropy profiles for stacked clusters, where the stacking is in velocity dispersion. The bands show the 68% and 95% regions, and the black curve is the median profile. Top: $z < 0.2$. Bottom: $z > 0.2$. The bottom-right plot shows the value of $\beta(r = r_{200})$ with 68% error bars (blue and red for $z < 0.2$ and $z > 0.2$, respectively). In all plots, the vertical dashed lines indicate the average $r_{ap,200}$ in the corresponding velocity dispersion bin.</figcaption></figure>

a complicated issue; this is because $\beta(r)$ is not a direct observable and it is degenerate with the total mass (Binney & Mamon 1982), constituting a possible source of systematic error in the reconstruction of dynamical mass profiles. To address this, both parametric approaches (e.g. Mamon et al. 2013; Read et al. 2021) and non-parametric methods (see e.g. Mamon et al. 2019 and references therein) have been developed. The MG-MAMPOSST technique (as does the original MAMPOSST) relies on a joint parametric modelling of both mass and anisotropy profiles, which has been shown to work robustly and adequately well in a wide range of cases (Mamon et al. 2013; Pizzuti et al. 2020; Biviano et al. 2023). However, the uncertainties in the reconstructed radial profile $\beta(r)$ are large for a single cluster, even for a considerably large number of members in the MG-MAMPOSST fit (see the example of PSZ2 G067.17+67.46 in the bottom plot of Fig. 2, with $N_{gal} = 299$ galaxies used in the analysis).

This issue can be solved in two ways. One possibility is to assume prior constraints on the mass profile, as provided, for example, by hydrostatic equilibrium or WL analyses (e.g. Annunziatella et al. 2016; Aguerri et al. 2017). On the other hand, one can infer information on the anisotropy profile as a function of cluster mass and redshift by populating the PPS with a stacking procedure. We adopted the latter approach and performed two different stacking procedures. First, we identified clusters with similar velocity dispersions in two redshift bins, namely $z \le 0.2$ and $0.2 < z < 0.6$ (which mostly splits clusters into the T1 or the T2 subsample). In the first bin, we further divided the clusters into three sub-bins of velocity dispersion, centred on $\sigma_{ap,200}/[\text{km s}^{-1}] = 750, 980, \text{ and } 1360$, respectively. We then generated the three stacked clusters (one for each bin of $\sigma_{ap,200}$) by considering all the member galaxies of CHEX-MATE clusters belonging to the same bin. The choices of the bins and the construction ensure that each stack contains $N \sim 3000$ galaxies.

For the high-redshift bin, we end up with only two stacked clusters, with $\sigma_{ap,200}/[\text{km s}^{-1}] = 1040, 1600$.

The second stacking has been performed considering five bins for clusters with similar SZ–X-ray calibrated masses – three for $z < 0.2$, $M_{500}^{SZ}/10^{14} \text{ M}_\odot = 3.23, 5.48, 7.74$, and two for $z > 0.2$, $M_{500}^{SZ}/10^{14} \text{ M}_\odot = 7.91, 12.10$. This second set is to consider a stacking criterion ($M_{500}^{SZ}$) that is independent of kinematic variables. As an example, the stacked PPS for this case are shown in Appendix C.

For each of the ten stacked clusters, we considered galaxies in the projected range $[0.05 \text{ Mpc}, 1.1 \langle r_{ap,200} \rangle]$, where $\langle r_{ap,200} \rangle$ is the average of $r_{ap,200}$ over all the clusters in the bin. We applied the MG-MAMPOSST method assuming a NFW model for the mass distribution and a pNFW for the projected number density distribution; as before, we performed a MCMC sampling of the parameter space $\{r_{200}, r_s, r_\nu, A_0, A_\infty\}$, adopting the same flat prior listed in Sect. 3. Figures 11 and 12 show the radial profiles of the velocity anisotropy obtained from the MG-MAMPOSST analysis of the stacked clusters ($\sigma_{ap,200}$ and SZ, respectively) at $z < 0.2$ (top panels in each plot) and $z > 0.2$ (bottom panels in each plot). For the stacking in velocity dispersion, we see that all the profiles (except the one for the most massive cluster, $\sigma_{ap,200} \sim 1600 \text{ km s}^{-1}$) are consistent with central isotropic orbits within one $\sigma$, and with an increase in the anisotropy towards the cluster outskirts, as found by other studies in the literature (Mamon et al. 2019).

In the bottom-right panel, we can observe a slight, although not very significant, trend to have more radial orbits at the virial radius for clusters with a larger velocity dispersion. The most massive stacked cluster exhibits a more ‘extreme’ behaviour, passing from tangential orbits ($\beta < 0$) at the centre, to largely radial ones in the outskirts. However, this bin contains clusters with $\delta\eta > 0.8$, which, as we have discussed, are characterised by disturbed dynamical states. The same trend of increasing radial
<figure><img src="image_14.png" /><figcaption>**Fig. 12.** Radial velocity anisotropy profiles for the stacked clusters, where the stacking is in bins of $M_{500}^{SZ}$. The structure is exactly as in Fig. 11. In all plots, the vertical dashed lines indicate $r_{500}$ as obtained from the average $M_{500}^{SZ}$ of the corresponding mass bin.</figcaption></figure>

orbits at virial radius is found in terms of the SZ mass, as highlighted in Fig. 12. In this case, however, only the first mass sub-bins in each redshift cut ($z < 0.2$ and $z > 0.2$) exhibit a central anisotropy consistent with zero at the 68% CL, while the others are in agreement with no anisotropy at the 95% CL, with a slight preference for central tangential orbits that is not very statistically significant.

## 6. Conclusions

We have studied the kinematical and dynamical properties of 75 massive galaxy clusters of the CHEX-MATE sample for which high-precision spectroscopic measurements of LOS velocities and projected positions are available. By means of the MG-MAMOOSST code, we jointly reconstructed mass profiles and velocity anisotropy profiles for eight combinations of mass and number density models. We find that the NFW prescription is preferred by the kinematic analysis for about half of the sample; the second most preferred is the Burkert model. This confirms the assumptions of the halo model used in the dynamical mass determinations of Sereno et al. (2025). Moreover, the constraints on the virial radius $r_{200}$ (and on the mass $M_{200}$) are robust against variations in the mass and number density models, as shown by the relatively small systematic uncertainties (Table A.1).

We compared the masses obtained from the SZ–X-ray-calibrated signal, finding $(1 - B_1) = 0.54 \pm 0.11$ when clusters with evident disturbed dynamical states are removed. A negligible bias is instead obtained when fitting with a linear relation that allows for mass dependence, $(1 - B) = 1.04 \pm 0.09$. When the clusters with a logarithmic mass difference $\delta\eta > 0.8$ dex (which are high-velocity dispersions $\sigma_{ap,200} > 1.4 \times 10^3$ km/s) are also considered, an overall positive mass-dependent bias, $(1 - B) = 0.67$, is found. Note that SZ masses may be overestimated due to the Malmquist bias, as discussed in Sereno et al. (2025); this might affect our findings, leading to a larger value of $(1 - B)$. A refined

analysis is beyond the scope of the current work, but one is presented in Sereno et al. (2025). We investigated the concentration–mass relation and compared it to theoretical expectations. By performing a Bayesian linear regression by means of the LIRA package, we found weak evidence of a rising trend at high masses, as also predicted by the Diemer & Joyce (2019) model and recently found by Butt et al. (2025) via a kinematic analysis of the Hectospec Cluster Survey (HeCS) cluster sample. Interestingly, the results obtained for the sample of CHEX-MATE clusters are fully consistent with the prediction of the Λ cold dark matter scenario, given an average value of $\langle c_{200} \rangle = 3.5$ with a scatter of 1.8. The results of the various scaling relations investigated in this work are summarised in Table 1.

We studied the dynamical state of the clusters by comparing $\delta\eta$ with two quantities that can be related with dynamical un-relaxation and disturbed morphology, namely the Anderson-Darling coefficient, $A^2$, and the fraction of galaxies in substructures, $f_{\text{sub}}$. No evident relations are found for $A^2$ and $f_{\text{sub}}$. Nevertheless, a lower averaged $\delta\eta$ is obtained for clusters with $A^2 < 0.44$ (0.44 being the median value of the sample). Note, however, that a more complete analysis should investigate the differences with the WL mass determinations and account for the X-ray morphology to better identify substructures and departures from equilibrium. This will be the subject of an upcoming work. Finally, we studied the shape of the radial velocity anisotropy profile by dividing the sample into two bins of redshift and collecting clusters within three and two sub-bins of $\sigma_{\text{ap,200}}$. We stacked clusters in the same bin and applied MG-MAMPOSST to those stacked systems to reconstruct the orbits from the centre to the edges. We performed a similar stacking procedure with bins of $M_{500}^{SZ}$, which are not directly related to the kinematical properties of the clusters. Consistent with the literature, orbits tend to be isotropic or tangential at the centre and more radial towards the edge of the clusters, where galaxies are thought to be falling into the cluster’s potential for the first time (Mamon et al. 2019).
<table><caption>**Table 1.** Results for the parameters of the linear fits of the relations discussed in this paper.</caption><thead><tr><th>Relation</th><th>Reference</th><th>Constraints</th></tr></thead><tbody><tr><td>Log M<sup>M</sup><sub>200</sub> - σ<sub>ap,200</sub></td><td>Eq. (9)</td><td>α<sub>σ</sub> = 3.61 ± 0.17, β<sub>σ</sub> = 1.07 ± 0.01, γ<sub>σ</sub> = -0.83 ± 0.70</td></tr><tr><td>Log M<sup>M</sup><sub>500</sub> - Log M<sup>SZ</sup><sub>500</sub> (outliers removed)</td><td>Eq. (12)<br>α<sub>SZ</sub> fixed to 1</td><td>α<sub>SZ</sub> = 0.71 ± 0.04, 1 - B = 1.04 ± 0.09<br>(1 - B<sub>1</sub>) = 0.54 ± 0.11</td></tr><tr><td>Log c<sub>200</sub> - Log M<sup>M</sup><sub>200</sub> (outliers removed)</td><td>Eq. (14)</td><td>α<sub>C</sub> = 0.49 ± 0.07, β<sub>C</sub> = 0.20 ± 0.31, γ<sub>C</sub> = -0.23 ± 1.6</td></tr></tbody></table>

There is a slight trend towards increasingly radial orbits at $r_{200}$ with larger $\sigma_{ap,200}$ (which, in turn, is connected to higher dynamical masses); we find the same trend with increasing $M_{500}^{SZ}$.

More broadly, our results indicate that the most massive galaxy clusters at $z < 0.6$ are generally close to dynamical equilibrium. Nonetheless, we identify some cases where even the most massive systems appear to deviate from full virialisation. As previously mentioned, a comprehensive assessment of the dynamical state of the CHEX-MATE sample will be carried out in a forthcoming study, which will incorporate a more complete set of morphological and kinematic indicators, as well as a comparison with mass estimates from WL and X-ray data.

**Acknowledgements.** This project has been partially funded by the Spanish Ministerio de Ciencia, Innovación y Universidades through the programme Generación del Concierto 2021, with code PID2021-122665OB-I00. R. Barrena acknowledges support by the Severo Ochoa 2020 research programme of the Instituto de Astrofísica de Canarias. S.E., M.S., acknowledge 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). M.S. acknowledges financial contributions from INAF Theory Grant 2023: Gravitational lensing detection of matter distribution at galaxy cluster boundaries and beyond (1.05.23.06.17). S.E. acknowledges the financial contribution from the European Union's Horizon 2020 Programme under the AHEAD2020 project (grant agreement n. 871158). DE acknowledges support from the Swiss National Science Foundation (SNSF) through grant agreement 200021_212576. C.P.H. acknowledges support from ANID through Fondecyt Regular project number 1252233, and from MILENIO - NCN2024_112. JS was supported by NASA Astrophysics Data Analysis Program (ADAP) Grant 80NSSC21K1571. BJM acknowledges support from Science and Technology Facilities Council grants ST/V000454/1 and ST/Y002008/1. GWP acknowledges long-term support from CNES, the French space agency. GC acknowledges the support from the Next Generation EU funds within the National Recovery and Resilience Plan (PNRR), Mission 4 – Education and Research, Component 2 – From Research to Business (M4C2), Investment Line 3.1 – Strengthening and creation of Research Infrastructures, Project IR0000012 – “CTA+ – Cherenkov Telescope Array Plus”. MD acknowledges the support of two NASA programs: NASA award 80NSSC19K0116/ SAO sub-award SV9-89010 and NASA award 80NSSC 22K0476.

## References

Abell, G. O. 1958, ApJS, 3, 211

Aguado-Barahona, A., Rubiño-Martín, J. A., Ferragamo, A., et al. 2022, A&A, 659, A126

Aguerri, J. A. L., Agulli, I., Diaferio, A., & Dalla Vecchia, C. 2017, MNRAS, 468, 364

Amodeo, S., Mei, S., Stanford, S. A., et al. 2017, ApJ, 844, 101

Anderson, T. W., & Darling, D. A. 1952, Ann. Math. Stat., 23, 193

Annunziatella, M., Mercurio, A., Biviano, A., et al. 2016, A&A, 585, A160

Barrena, R., Pizzuti, L., Chon, G., & Böhringer, H. 2024, A&A, 691, A135

Beers, T. C., Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32

Benavides, J. A., Biviano, A., & Abadi, M. G. 2023, A&A, 669, A147

Binney, J., & Mamon, G. A. 1982, MNRAS, 200, 361

Biviano, A., & Katgert, P. 2004, A&A, 424, 779

Biviano, A., Moretti, A., Paccagnella, A., et al. 2017, A&A, 607, A81

Biviano, A., Pizzuti, L., Mercurio, A., et al. 2023, ApJ, 958, 148

Biviano, A., Poggianti, B. M., Jaffé, Y., et al. 2024, ApJ, 965, 117

Burkert, A. 1995, ApJ, 447, L25

Butt, M. A., Haridasu, S., Diaferio, A., et al. 2025, arXiv e-prints [arXiv:2504.16685]

Campitiello, M. G., Ettori, S., Lovisari, L., et al. 2022, A&A, 665, A117

Carlb erg, R. G., Yee, H. K. C., & Ellingson, E. 1997, ApJ, 478, 462

Cavaliere, A., & Fusco-Femiano, R. 1981, A&A, 100, 194

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

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

Dressler, A., & Shectman, S. A. 1988, AJ, 95, 985

Dutton, A. A., & Macciò, A. V. 2014, MNRAS, 441, 3359

Einasto, J. 1965, Trudy Astrofizicheskogo Instituta Alma-Ata, 5, 87

Ettori, S., Donnarumma, A., Pointecouteau, E., et al. 2013, Space Sci. Rev., 177, 119

Ferragamo, A., Rubiño-Martín, J. A., Betancort-Rijo, J., et al. 2020, A&A, 641, A11

Ferragamo, A., Barrena, R., Rubiño-Martín, J. A., et al. 2021, A&A, 655, A115

Gelman, A., & Rubin, D. B. 1992, Stat. Sci., 7, 457

Hernquist, L. 1990, ApJ, 356, 359

Kimmig, L. C., Remus, R.-S., Dolag, K., & Biffi, V. 2023, ApJ, 949, 92

Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128

Lesci, G. F., Veropalumbo, A., Sereno, M., et al. 2023, A&A, 674, A80

Lotz, M., Remus, R.-S., Dolag, K., Biviano, A., & Burkert, A. 2019, MNRAS, 488, 5370

Macciò, A. V., Dutton, A. A., & van den Bosch, F. C. 2008, MNRAS, 391, 1940

Mamon, G. A., Biviano, A., & Murante, G. 2010, A&A, 520, A30

Mamon, G. A., Biviano, A., & Boué, G. 2013, MNRAS, 429, 3079

Mamon, G. A., Cava, A., Biviano, A., et al. 2019, A&A, 631, A131

Meneghetti, M., Rasia, E., Vega, J., Merten, J., & Postman, M. E. A. 2014, ApJ, 797, 34

Merten, J., Meneghetti, M., Postman, M., et al. 2015, ApJ, 806, 4

Munari, E., Biviano, A., Borgani, S., Murante, G., & Fabjan, D. 2013, MNRAS, 430, 2638

Munari, E., Biviano, A., & Mamon, G. A. 2014, A&A, 566, A68

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493

Parekh, V., van der Heyden, K., Ferrari, C., Angus, G., & Holwerda, B. 2015, A&A, 575, A127

Penna-Lima, M., Bartlett, J. G., Rozo, E., et al. 2017, A&A, 604, A89

Pizzuti, L., Sartoris, B., Borgani, S., & Biviano, A. 2020, JCAP, 04, 024

Pizzuti, L., Saltas, I. D., & Amendola, L. 2021, MNRAS, 506, 595

Pizzuti, L., Saltas, I. D., Biviano, A., Mamon, G., & Amendola, L. 2023, JOSS, 8, 4800

Pizzuti, L., Boumechta, Y., Haridasu, S., et al. 2024, J. Cosmol. Astropart. Phys., 2024, 014

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

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

Ragagnin, A., Saro, A., Singh, P., & Dolag, K. 2020, MNRAS, 500, 5056

Read, J. I., Mamon, G. A., Vasiliev, E., et al. 2021, MNRAS, 501, 978

Roberts, I. D., Parker, L. C., & Hlavacek-Larrondo, J. 2018, MNRAS, 475, 4704

Ruel, J., Bazin, G., Bayliss, M., et al. 2014, ApJ, 792, 45

Sarazin, C. L. 1980, ApJ, 236, 75

Sereno, M. 2016, MNRAS, 455, 2149

Sereno, M., & Ettori, S. 2015a, MNRAS, 450, 3675 (CoMaLit-IV)

Sereno, M., & Ettori, S. 2015b, MNRAS, 450, 3633 (CoMaLit-I)

Sereno, M., & Ettori, S. 2017, MNRAS, 468, 3322 (CoMaLit-V)

Sereno, M., Ettori, S., & Moscardini, L. 2015a, MNRAS, 450, 3649 (CoMaLit-II)

Sereno, M., Giocoli, C., Ettori, S., & Moscardini, L. 2015b, MNRAS, 449, 2024

Sereno, M., Maurogordato, S., Cappi, A., et al. 2025, A&A, 693, A2

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

Tiret, O., Combes, F., Angus, G. W., Famaey, B., & Zhao, H. S. 2007, A&A, 476, L1

Tonnesen, S. 2019, ApJ, 874, 161

van der Burg, R. F. J., Muzzin, A., & Hoekstra, H. 2016, A&A, 590, A20

Wen, Z. L., & Han, J. L. 2013, MNRAS, 436, 275

Zarattini, S., Biviano, A., Aguerri, J. A. L., Girardi, M., & D’Onghia, E. 2021, A&A, 655, A103

Zwicky, F. 1933, Helvetica Physica Acta, 6, 110
## Appendix A: Results from MG-MAMPOSST

The constraints on the mass profile, number density and velocity anisotropy parameters are summarised in Table A.1, along with cluster redshift, Anderson-Darling coefficient $A^2$, fraction of galaxies in substructures $f_{sub}$, number of galaxies $N_{200}$ within a projected cylinder $R < r_{ap,200}$ aperture velocity dispersion. Here we also include the MG-MAMPOSST analysis of clusters PSZ2G113.91-37.01 and PSZ2G083.29-31.03, which have fewer than $N = 50$ members within $r_{ap,200}$.

<table><caption>**Table A.1.** Constraints on the free parameters from the kinematic analysis with MG-MAMPOSST, shown for the combination of mass and number density model that provides the best BIC.</caption><thead><tr><th>PSZ2 name</th><th>Mass</th><th>N(R)</th><th>r<sub>200</sub><br>[Mpc]</th><th>r<sub>v</sub><br>[Mpc]</th><th>r<sub>s</sub><br>[Mpc]</th><th>A<sub>0</sub></th><th>-ln P</th><th>z</th><th>A<sup>2</sup></th><th>f<sub>sub</sub></th><th>N<sub>200</sub></th><th>σ<sub>ap,200</sub><br>[Km/s]</th></tr></thead><tbody><tr><td>G049.22+30.87</td><td>NFW</td><td>pNFW</td><td>2.13<sup>+0.25</sup><sub>-0.27</sub> ± 0.06</td><td>0.59<sup>+0.10</sup><sub>-0.13</sub> ± 0.28</td><td>0.81<sup>+1.22</sup><sub>-0.72</sub> ± 0.73</td><td>1.78<sup>+1.76</sup><sub>-1.33</sub> ± 0.10</td><td>3403.919</td><td>0.1606</td><td>0.27</td><td>0.29</td><td>368</td><td>1040 ± 40</td></tr><tr><td>G226.18+76.79</td><td>Bur</td><td>pHer</td><td>2.29<sup>+0.53</sup><sub>-0.44</sub> ± 0.04</td><td>0.67<sup>+0.17</sup><sub>-0.13</sub> ± 0.24</td><td>0.30<sup>+0.41</sup><sub>-0.28</sub> ± 0.38</td><td>1.64<sup>+1.94</sup><sub>-1.18</sub> ± 0.04</td><td>1171.266</td><td>0.1415</td><td>0.44</td><td>0.09</td><td>131</td><td>1027 ± 75</td></tr><tr><td>G053.53+59.52</td><td>NFW</td><td>pHer</td><td>2.40<sup>+0.36</sup><sub>-0.32</sub> ± 0.08</td><td>0.87<sup>+0.14</sup><sub>-0.15</sub> ± 0.30</td><td>0.85<sup>+1.21</sup><sub>-0.73</sub> ± 0.66</td><td>1.58<sup>+1.90</sup><sub>-1.10</sub> ± 0.12</td><td>2100.914</td><td>0.1136</td><td>0.78</td><td>0.12</td><td>222</td><td>1041 ± 46</td></tr><tr><td>G056.77+36.32</td><td>NFW</td><td>pNFW</td><td>3.41<sup>+0.81</sup><sub>-0.76</sub> ± 0.15</td><td>0.85<sup>+0.29</sup><sub>-0.16</sub> ± 0.74</td><td>2.09<sup>+1.64</sup><sub>-1.54</sub> ± 0.46</td><td>1.38<sup>+2.14</sup><sub>-0.91</sub> ± 0.14</td><td>981.519</td><td>0.09696</td><td>0.64</td><td>0.19</td><td>96</td><td>1364 ± 76</td></tr><tr><td>G042.81+56.61</td><td>NFW</td><td>pNFW</td><td>2.48<sup>+0.31</sup><sub>-0.29</sub> ± 0.08</td><td>0.53<sup>+0.17</sup><sub>-0.15</sub> ± 0.49</td><td>0.63<sup>+0.76</sup><sub>-0.52</sub> ± 0.77</td><td>1.63<sup>+1.78</sup><sub>-1.06</sub> ± 0.08</td><td>2014.210</td><td>0.07225</td><td>0.54</td><td>0.30</td><td>210</td><td>1148 ± 47</td></tr><tr><td>G033.81+77.18</td><td>Bur</td><td>pNFW</td><td>1.79<sup>+0.31</sup><sub>-0.25</sub> ± 0.05</td><td>0.34<sup>+0.12</sup><sub>-0.12</sub> ± 0.34</td><td>0.40<sup>+0.56</sup><sub>-0.37</sub> ± 0.61</td><td>1.48<sup>+1.82</sup><sub>-1.01</sub> ± 0.08</td><td>1293.172</td><td>0.06272</td><td>0.56</td><td>0.23</td><td>154</td><td>786 ± 49</td></tr><tr><td>G041.45+29.10</td><td>Bur</td><td>pHer</td><td>2.40<sup>+0.52</sup><sub>-0.44</sub> ± 0.15</td><td>0.38<sup>+0.10</sup><sub>-0.08</sub> ± 0.14</td><td>0.12<sup>+0.15</sup><sub>-0.08</sub> ± 0.15</td><td>1.74<sup>+1.78</sup><sub>-1.22</sub> ± 0.09</td><td>961.440</td><td>0.1771</td><td>0.28</td><td>0.37</td><td>110</td><td>1223 ± 85</td></tr><tr><td>G238.69+63.26</td><td>Her</td><td>pHer</td><td>1.83<sup>+0.28</sup><sub>-0.26</sub> ± 0.07</td><td>1.95<sup>+0.20</sup><sub>-0.29</sub> ± 0.52</td><td>1.56<sup>+1.74</sup><sub>-1.30</sub> ± 0.46</td><td>1.51<sup>+2.00</sup><sub>-1.05</sub> ± 0.08</td><td>1770.555</td><td>0.1672</td><td>0.45</td><td>0.21</td><td>229</td><td>871 ± 46</td></tr><tr><td>G080.41-33.24</td><td>NFW</td><td>pHer</td><td>2.13<sup>+0.62</sup><sub>-0.55</sub> ± 0.18</td><td>0.36<sup>+0.08</sup><sub>-0.08</sub> ± 0.14</td><td>0.99<sup>+1.81</sup><sub>-0.97</sub> ± 0.69</td><td>1.34<sup>+1.74</sup><sub>-1.62</sub> ± 0.11</td><td>1124.055</td><td>0.1101</td><td>0.54</td><td>0.10</td><td>137</td><td>797 ± 50</td></tr><tr><td>G048.10+57.16</td><td>Bur</td><td>pHer</td><td>1.87<sup>+0.36</sup><sub>-0.33</sub> ± 0.06</td><td>1.11<sup>+0.24</sup><sub>-0.15</sub> ± 0.50</td><td>0.64<sup>+1.10</sup><sub>-0.63</sub> ± 0.73</td><td>1.49<sup>+1.84</sup><sub>-1.03</sub> ± 0.09</td><td>1546.255</td><td>0.07782</td><td>0.48</td><td>0.22</td><td>165</td><td>811 ± 52</td></tr><tr><td>G287.46+81.12</td><td>Eis</td><td>pHer</td><td>1.69<sup>+0.31</sup><sub>-0.29</sub> ± 0.06</td><td>1.72<sup>+0.36</sup><sub>-0.46</sub> ± 0.50</td><td>0.64<sup>+1.47</sup><sub>-0.64</sub> ± 0.73</td><td>1.85<sup>+1.83</sup><sub>-1.32</sub> ± 0.09</td><td>964.668</td><td>0.07194</td><td>0.25</td><td>0.28</td><td>109</td><td>813 ± 57</td></tr><tr><td>G031.93+78.71</td><td>Bur</td><td>pNFW</td><td>1.28<sup>+0.56</sup><sub>-0.73</sub> ± 0.06</td><td>0.35<sup>+0.31</sup><sub>-0.18</sub> ± 0.50</td><td>0.82<sup>+1.33</sup><sub>-0.82</sub> ± 0.73</td><td>1.76<sup>+1.85</sup><sub>-1.27</sub> ± 0.09</td><td>400.254</td><td>0.07533</td><td>0.28</td><td>0.20</td><td>57</td><td>560 ± 57</td></tr><tr><td>G067.17+67.46</td><td>NFW</td><td>pNFW</td><td>2.19<sup>+0.28</sup><sub>-0.29</sub> ± 0.06</td><td>0.74<sup>+0.20</sup><sub>-0.19</sub> ± 0.50</td><td>0.85<sup>+1.34</sup><sub>-0.75</sub> ± 0.73</td><td>1.55<sup>+1.72</sup><sub>-1.01</sub> ± 0.09</td><td>2880.857</td><td>0.1667</td><td>1.57</td><td>0.28</td><td>291</td><td>1077 ± 43</td></tr><tr><td>G006.49+50.56</td><td>NFW</td><td>pNFW</td><td>2.13<sup>+0.16</sup><sub>-0.15</sub> ± 0.06</td><td>1.08<sup>+0.19</sup><sub>-0.14</sub> ± 0.50</td><td>0.40<sup>+0.35</sup><sub>-0.27</sub> ± 0.73</td><td>2.02<sup>+1.58</sup><sub>-0.18</sub> ± 0.09</td><td>6711.942</td><td>0.07854</td><td>0.89</td><td>0.26</td><td>750</td><td>1083 ± 28</td></tr><tr><td>G124.20-36.48</td><td>Her</td><td>pHer</td><td>2.61<sup>+0.39</sup><sub>-0.37</sub> ± 0.06</td><td>1.22<sup>+0.11</sup><sub>-0.15</sub> ± 0.50</td><td>1.23<sup>+1.65</sup><sub>-1.10</sub> ± 0.73</td><td>1.55<sup>+1.86</sup><sub>-1.08</sub> ± 0.09</td><td>2031.290</td><td>0.1511</td><td>0.95</td><td>0.08</td><td>74</td><td>1240 ± 130</td></tr><tr><td>G021.10+33.24</td><td>Eis</td><td>pNFW</td><td>2.34<sup>+0.52</sup><sub>-0.75</sub> ± 0.06</td><td>1.79<sup>+1.38</sup><sub>-0.98</sub> ± 0.50</td><td>0.70<sup>+2.01</sup><sub>-0.68</sub> ± 0.73</td><td>1.59<sup>+2.05</sup><sub>-1.12</sub> ± 0.09</td><td>643.999</td><td>0.1511</td><td>0.45</td><td>0.35</td><td>232</td><td>1248 ± 63</td></tr><tr><td>G044.20+48.66</td><td>NFW</td><td>pHer</td><td>2.38<sup>+0.26</sup><sub>-0.28</sub> ± 0.06</td><td>2.33<sup>+0.14</sup><sub>-0.21</sub> ± 0.50</td><td>1.60<sup>+1.59</sup><sub>-1.20</sub> ± 0.73</td><td>1.42<sup>+1.74</sup><sub>-0.88</sub> ± 0.09</td><td>7619.565</td><td>0.08987</td><td>0.44</td><td>0.33</td><td>778</td><td>1102 ± 29</td></tr><tr><td>G075.71+13.51</td><td>Her</td><td>pHer</td><td>3.97<sup>+0.88</sup><sub>-0.85</sub> ± 0.06</td><td>0.50<sup>+0.13</sup><sub>-0.11</sub> ± 0.50</td><td>1.64<sup>+1.47</sup><sub>-1.18</sub> ± 0.73</td><td>1.26<sup>+1.95</sup><sub>-0.79</sub> ± 0.09</td><td>1249.919</td><td>0.05455</td><td>0.21</td><td>0.12</td><td>134</td><td>1586 ± 92</td></tr><tr><td>G324.04+48.79</td><td>Eis</td><td>pHer</td><td>1.79<sup>+0.48</sup><sub>-0.42</sub> ± 0.06</td><td>0.92<sup>+0.42</sup><sub>-0.29</sub> ± 0.50</td><td>0.97<sup>+2.04</sup><sub>-0.98</sub> ± 0.73</td><td>1.61<sup>+1.95</sup><sub>-1.14</sub> ± 0.09</td><td>634.364</td><td>0.4505</td><td>0.12</td><td>0.00</td><td>78</td><td>937 ± 80</td></tr><tr><td>G284.41+52.45</td><td>NFW</td><td>pHer</td><td>2.07<sup>+0.23</sup><sub>-0.26</sub> ± 0.04</td><td>2.19<sup>+0.40</sup><sub>-0.45</sub> ± 0.24</td><td>0.90<sup>+1.56</sup><sub>-0.85</sub> ± 0.38</td><td>1.49<sup>+1.52</sup><sub>-0.85</sub> ± 0.04</td><td>2983.187</td><td>0.4399</td><td>0.40</td><td>0.33</td><td>333</td><td>1217 ± 50</td></tr><tr><td>G349.46-59.95</td><td>NFW</td><td>pHer</td><td>3.46<sup>+1.16</sup><sub>-1.09</sub> ± 0.04</td><td>0.62<sup>+0.27</sup><sub>-0.22</sub> ± 0.24</td><td>1.35<sup>+1.92</sup><sub>-1.34</sub> ± 0.38</td><td>1.31<sup>+1.87</sup><sub>-0.85</sub> ± 0.04</td><td>719.670</td><td>0.3465</td><td>0.54</td><td>0.19</td><td>81</td><td>1440 ± 110</td></tr><tr><td>G008.94-81.22</td><td>Bur</td><td>pHer</td><td>3.05<sup>+0.32</sup><sub>-0.30</sub> ± 0.04</td><td>0.90<sup>+0.09</sup><sub>-0.06</sub> ± 0.24</td><td>0.26<sup>+0.18</sup><sub>-0.17</sub> ± 0.38</td><td>2.24<sup>+1.62</sup><sub>-1.41</sub> ± 0.04</td><td>4631.269</td><td>0.3064</td><td>0.87</td><td>0.39</td><td>491</td><td>1657 ± 51</td></tr><tr><td>G266.04-21.25</td><td>NFW</td><td>pHer</td><td>2.98<sup>+0.83</sup><sub>-0.80</sub> ± 0.04</td><td>0.46<sup>+0.10</sup><sub>-0.15</sub> ± 0.24</td><td>1.25<sup>+1.90</sup><sub>-1.21</sub> ± 0.38</td><td>1.43<sup>+1.85</sup><sub>-0.98</sub> ± 0.04</td><td>838.472</td><td>0.2965</td><td>0.25</td><td>0.15</td><td>93</td><td>1200 ± 94</td></tr><tr><td>G107.10+65.32</td><td>Bur</td><td>pHer</td><td>3.00<sup>+0.34</sup><sub>-0.32</sub> ± 0.04</td><td>1.74<sup>+0.10</sup><sub>-0.15</sub> ± 0.24</td><td>0.69<sup>+0.67</sup><sub>-0.51</sub> ± 0.38</td><td>1.62<sup>+1.92</sup><sub>-1.14</sub> ± 0.04</td><td>3463.596</td><td>0.2772</td><td>1.05</td><td>0.41</td><td>348</td><td>1475 ± 56</td></tr></tbody></table>
<table><caption>Table A.1. continued.</caption><thead><tr><th>PSZ2 name</th><th>Mass</th><th>N(R)</th><th>r<sup>200</sup> [Mpc]</th><th>r<sub>v</sub> [Mpc]</th><th>r<sub>s</sub> [Mpc]</th><th>A<sub>∞</sub></th><th>A<sub>0</sub></th><th>-ln P</th><th>z</th><th>A<sup>2</sup></th><th>f<sub>sub</sub></th><th>N<sub>200</sub></th><th>σ<sub>ap,200</sub> [km/s]</th></tr></thead><tbody><tr><td>G340.36+60.58</td><td>Eis</td><td>pHer</td><td>2.84<sup>+0.42</sup><sub>-0.49</sub> ± 0.04</td><td>3.56<sup>+0.92</sup><sub>-0.75</sub> ± 0.24</td><td>1.53<sup>+1.77</sup><sub>-1.17</sub> ± 0.38</td><td>1.16<sup>+1.73</sup><sub>-0.69</sub> ± 0.04</td><td>1.72<sup>+2.47</sup><sub>-1.26</sub> ± 0.14</td><td>3070.106</td><td>0.251</td><td>2.38</td><td>0.38</td><td>368</td><td>1452 ± 62</td></tr><tr><td>G072.62+41.46</td><td>NFW</td><td>pHer</td><td>2.52<sup>+0.33</sup><sub>-0.33</sub> ± 0.04</td><td>2.13<sup>+0.42</sup><sub>-0.34</sub> ± 0.24</td><td>1.54<sup>+1.81</sup><sub>-1.39</sub> ± 0.38</td><td>1.23<sup>+1.55</sup><sub>-0.72</sub> ± 0.04</td><td>2.11<sup>+2.30</sup><sub>-1.44</sub> ± 0.14</td><td>3316.267</td><td>0.2255</td><td>0.34</td><td>0.19</td><td>368</td><td>1271 ± 51</td></tr><tr><td>G055.59+31.85</td><td>NFW</td><td>pHer</td><td>2.03<sup>+0.51</sup><sub>-0.50</sub> ± 0.04</td><td>3.74<sup>+1.12</sup><sub>-1.28</sub> ± 0.24</td><td>1.44<sup>+1.90</sup><sub>-1.36</sub> ± 0.38</td><td>1.55<sup>+1.94</sup><sub>-1.09</sub> ± 0.04</td><td>1.38<sup>+2.68</sup><sub>-0.93</sub> ± 0.14</td><td>1095.847</td><td>0.2255</td><td>0.37</td><td>0.21</td><td>126</td><td>1060 ± 63</td></tr><tr><td>G159.91-73.50</td><td>Eis</td><td>pHer</td><td>2.70<sup>+0.31</sup><sub>-0.29</sub> ± 0.08</td><td>1.59<sup>+0.32</sup><sub>-0.27</sub> ± 0.30</td><td>0.74<sup>+0.78</sup><sub>-0.57</sub> ± 0.66</td><td>1.74<sup>+1.91</sup><sub>-1.28</sub> ± 0.12</td><td>1.04<sup>+2.12</sup><sub>-0.56</sub> ± 0.16</td><td>2603.277</td><td>0.2095</td><td>1.19</td><td>0.36</td><td>272</td><td>1331 ± 65</td></tr><tr><td>G195.75-24.32</td><td>NFW</td><td>pHer</td><td>2.15<sup>+0.33</sup><sub>-0.30</sub> ± 0.08</td><td>1.32<sup>+0.24</sup><sub>-0.29</sub> ± 0.30</td><td>0.98<sup>+1.74</sup><sub>-0.99</sub> ± 0.66</td><td>1.49<sup>+1.84</sup><sub>-1.02</sub> ± 0.12</td><td>1.79<sup>+2.37</sup><sub>-1.34</sub> ± 0.16</td><td>1430.192</td><td>0.2007</td><td>0.60</td><td>0.32</td><td>153</td><td>1051 ± 60</td></tr><tr><td>G313.33+61.13</td><td>Eis</td><td>pHer</td><td>2.87<sup>+0.36</sup><sub>-0.32</sub> ± 0.08</td><td>0.93<sup>+0.12</sup><sub>-0.09</sub> ± 0.30</td><td>0.28<sup>+0.20</sup><sub>-0.18</sub> ± 0.66</td><td>2.24<sup>+1.62</sup><sub>-1.40</sub> ± 0.12</td><td>1.16<sup>+1.74</sup><sub>-0.71</sub> ± 0.16</td><td>4277.167</td><td>0.1848</td><td>1.00</td><td>0.40</td><td>492</td><td>1620 ± 63</td></tr><tr><td>G208.80-30.67</td><td>Bur</td><td>pHer</td><td>2.78<sup>+1.09</sup><sub>-0.76</sub> ± 0.08</td><td>0.83<sup>+0.31</sup><sub>-0.26</sub> ± 0.30</td><td>0.68<sup>+1.46</sup><sub>-0.65</sub> ± 0.66</td><td>1.61<sup>+1.93</sup><sub>-1.13</sub> ± 0.12</td><td>1.44<sup>+2.52</sup><sub>-0.98</sub> ± 0.16</td><td>1141.373</td><td>0.2468</td><td>0.89</td><td>0.27</td><td>130</td><td>1177 ± 91</td></tr><tr><td>G172.98-53.55</td><td>NFW</td><td>pHer</td><td>3.08<sup>+0.77</sup><sub>-0.67</sub> ± 0.08</td><td>0.24<sup>+0.03</sup><sub>-0.03</sub> ± 0.30</td><td>0.39<sup>+0.45</sup><sub>-0.32</sub> ± 0.66</td><td>1.41<sup>+1.96</sup><sub>-0.95</sub> ± 0.12</td><td>1.48<sup>+2.47</sup><sub>-1.03</sub> ± 0.16</td><td>2412.888</td><td>0.374</td><td>0.44</td><td>0.36</td><td>272</td><td>1526 ± 63</td></tr><tr><td>G205.93-39.46</td><td>NFW</td><td>pNFW</td><td>2.57<sup>+0.69</sup><sub>-0.62</sub> ± 0.08</td><td>0.57<sup>+0.46</sup><sub>-0.28</sub> ± 0.30</td><td>1.48<sup>+1.91</sup><sub>-1.43</sub> ± 0.66</td><td>1.33<sup>+1.85</sup><sub>-0.87</sub> ± 0.12</td><td>1.99<sup>+2.42</sup><sub>-1.52</sub> ± 0.16</td><td>544.822</td><td>0.4409</td><td>0.36</td><td>0.07</td><td>65</td><td>1450 ± 150</td></tr><tr><td>G273.59+63.27</td><td>Bur</td><td>pHer</td><td>2.63<sup>+0.35</sup><sub>-0.33</sub> ± 0.15</td><td>1.74<sup>+0.48</sup><sub>-0.34</sub> ± 0.74</td><td>0.49<sup>+0.43</sup><sub>-0.33</sub> ± 0.46</td><td>2.08<sup>+1.71</sup><sub>-1.42</sub> ± 0.14</td><td>0.96<sup>+1.21</sup><sub>-0.49</sub> ± 0.39</td><td>1685.399</td><td>0.1335</td><td>0.44</td><td>0.19</td><td>180</td><td>1275 ± 56</td></tr><tr><td>G057.61+34.93</td><td>NFW</td><td>pHer</td><td>2.19<sup>+0.35</sup><sub>-0.32</sub> ± 0.15</td><td>1.13<sup>+0.28</sup><sub>-0.28</sub> ± 0.74</td><td>0.62<sup>+1.30</sup><sub>-0.62</sub> ± 0.46</td><td>1.52<sup>+1.73</sup><sub>-0.98</sub> ± 0.14</td><td>1.93<sup>+2.41</sup><sub>-1.47</sub> ± 0.39</td><td>1232.084</td><td>0.08492</td><td>0.50</td><td>0.22</td><td>141</td><td>993 ± 49</td></tr><tr><td>G046.88+56.48</td><td>NFW</td><td>pHer</td><td>2.34<sup>+0.43</sup><sub>-0.43</sub> ± 0.15</td><td>2.69<sup>+0.80</sup><sub>-0.54</sub> ± 0.74</td><td>2.29<sup>+1.46</sup><sub>-1.36</sub> ± 0.46</td><td>1.60<sup>+1.96</sup><sub>-1.13</sub> ± 0.14</td><td>0.84<sup>+0.54</sup><sub>-0.36</sub> ± 0.39</td><td>1834.999</td><td>0.1137</td><td>0.61</td><td>0.27</td><td>204</td><td>1069 ± 46</td></tr><tr><td>G243.64+67.74</td><td>Bur</td><td>pNFW</td><td>1.41<sup>+0.56</sup><sub>-0.86</sub> ± 0.15</td><td>0.61<sup>+0.43</sup><sub>-0.25</sub> ± 0.74</td><td>0.84<sup>+1.52</sup><sub>-0.84</sub> ± 0.46</td><td>1.61<sup>+1.96</sup><sub>-1.15</sub> ± 0.14</td><td>2.30<sup>+2.33</sup><sub>-1.81</sub> ± 0.39</td><td>488.051</td><td>0.08108</td><td>0.55</td><td>0.30</td><td>54</td><td>781 ± 86</td></tr><tr><td>G008.31-64.74</td><td>NFW</td><td>pHer</td><td>4.43<sup>+0.57</sup><sub>-0.61</sub> ± 0.15</td><td>0.68<sup>+0.11</sup><sub>-0.12</sub> ± 0.74</td><td>1.39<sup>+1.14</sup><sub>-0.89</sub> ± 0.46</td><td>1.04<sup>+1.80</sup><sub>-0.57</sub> ± 0.14</td><td>1.85<sup>+2.43</sup><sub>-1.33</sub> ± 0.39</td><td>2369.679</td><td>0.3153</td><td>0.87</td><td>0.26</td><td>249</td><td>1884 ± 79</td></tr><tr><td>G285.63+72.75</td><td>Eis</td><td>pHer</td><td>2.12<sup>+0.36</sup><sub>-0.33</sub> ± 0.15</td><td>1.36<sup>+0.33</sup><sub>-0.36</sub> ± 0.74</td><td>0.95<sup>+1.78</sup><sub>-0.89</sub> ± 0.46</td><td>1.58<sup>+1.95</sup><sub>-1.12</sub> ± 0.14</td><td>1.74<sup>+2.48</sup><sub>-1.29</sub> ± 0.39</td><td>1243.294</td><td>0.1667</td><td>0.34</td><td>0.14</td><td>140</td><td>1000 ± 70</td></tr><tr><td>G077.90-26.63</td><td>Bur</td><td>pNFW</td><td>1.86<sup>+0.61</sup><sub>-0.88</sub> ± 0.08</td><td>0.45<sup>+0.24</sup><sub>-0.15</sub> ± 0.49</td><td>0.73<sup>+1.60</sup><sub>-0.73</sub> ± 0.77</td><td>1.71<sup>+1.83</sup><sub>-1.24</sub> ± 0.08</td><td>2.29<sup>+2.28</sup><sub>-1.78</sub> ± 0.10</td><td>725.684</td><td>0.1449</td><td>0.29</td><td>0.14</td><td>81</td><td>930 ± 89</td></tr><tr><td>G114.79-33.71</td><td>NFW</td><td>pNFW</td><td>2.06<sup>+0.37</sup><sub>-0.34</sub> ± 0.08</td><td>0.44<sup>+0.19</sup><sub>-0.14</sub> ± 0.49</td><td>0.87<sup>+1.51</sup><sub>-0.80</sub> ± 0.77</td><td>1.49<sup>+1.87</sup><sub>-1.02</sub> ± 0.08</td><td>1.64<sup>+2.49</sup><sub>-1.19</sub> ± 0.10</td><td>1316.109</td><td>0.09447</td><td>0.20</td><td>0.21</td><td>155</td><td>927 ± 50</td></tr><tr><td>G224.00+69.33</td><td>NFW</td><td>pHer</td><td>2.28<sup>+0.50</sup><sub>-0.44</sub> ± 0.05</td><td>0.96<sup>+0.31</sup><sub>-0.24</sub> ± 0.34</td><td>0.88<sup>+1.96</sup><sub>-0.88</sub> ± 0.61</td><td>1.45<sup>+1.88</sup><sub>-0.99</sub> ± 0.08</td><td>2.06<sup>+2.36</sup><sub>-1.58</sub> ± 0.17</td><td>983.873</td><td>0.1917</td><td>0.35</td><td>0.20</td><td>114</td><td>1038 ± 69</td></tr><tr><td>G113.91-37.01</td><td>Eis</td><td>pHer</td><td>2.87<sup>+1.07</sup><sub>-0.96</sub> ± 0.05</td><td>0.59<sup>+0.13</sup><sub>-0.17</sub> ± 0.34</td><td>1.15<sup>+1.97</sup><sub>-1.15</sub> ± 0.61</td><td>1.60<sup>+1.91</sup><sub>-1.14</sub> ± 0.08</td><td>2.12<sup>+2.38</sup><sub>-1.65</sub> ± 0.17</td><td>553.347</td><td>0.3677</td><td>0.57</td><td>0.00</td><td>19</td><td>1370 ± 200</td></tr><tr><td>G067.52+34.75</td><td>NFW</td><td>pHer</td><td>2.03<sup>+0.44</sup><sub>-0.39</sub> ± 0.05</td><td>1.96<sup>+0.68</sup><sub>-0.62</sub> ± 0.34</td><td>1.27<sup>+1.87</sup><sub>-1.24</sub> ± 0.61</td><td>1.54<sup>+1.90</sup><sub>-1.06</sub> ± 0.08</td><td>1.63<sup>+2.34</sup><sub>-1.19</sub> ± 0.17</td><td>940.716</td><td>0.1822</td><td>0.41</td><td>0.17</td><td>116</td><td>1041 ± 61</td></tr><tr><td>G113.29-29.69</td><td>NFW</td><td>pNFW</td><td>1.86<sup>+0.30</sup><sub>-0.28</sub> ± 0.05</td><td>0.41<sup>+0.15</sup><sub>-0.16</sub> ± 0.34</td><td>0.42<sup>+0.83</sup><sub>-0.41</sub> ± 0.61</td><td>1.57<sup>+1.81</sup><sub>-1.06</sub> ± 0.08</td><td>1.84<sup>+2.42</sup><sub>-1.39</sub> ± 0.17</td><td>958.728</td><td>0.103</td><td>0.45</td><td>0.14</td><td>114</td><td>892 ± 51</td></tr><tr><td>G192.18+56.12</td><td>Bur</td><td>pHer</td><td>1.66<sup>+0.79</sup><sub>-1.06</sub> ± 0.05</td><td>2.20<sup>+0.63</sup><sub>-0.73</sub> ± 0.34</td><td>1.72<sup>+1.35</sup><sub>-0.88</sub> ± 0.61</td><td>1.35<sup>+2.00</sup><sub>-0.88</sub> ± 0.08</td><td>1.42<sup>+2.25</sup><sub>-0.98</sub> ± 0.17</td><td>834.688</td><td>0.1268</td><td>0.70</td><td>0.25</td><td>104</td><td>824 ± 59</td></tr><tr><td>G179.09+60.12</td><td>Bur</td><td>pNFW</td><td>2.53<sup>+0.93</sup><sub>-1.26</sub> ± 0.05</td><td>0.67<sup>+0.40</sup><sub>-0.30</sub> ± 0.34</td><td>1.00<sup>+0.99</sup><sub>-0.96</sub> ± 0.61</td><td>1.57<sup>+1.92</sup><sub>-1.10</sub> ± 0.08</td><td>2.11<sup>+2.30</sup><sub>-1.63</sub> ± 0.17</td><td>808.742</td><td>0.1397</td><td>0.14</td><td>0.12</td><td>76</td><td>1130 ± 100</td></tr><tr><td>G283.91+73.87</td><td>Bur</td><td>pHer</td><td>1.79<sup>+0.58</sup><sub>-0.83</sub> ± 0.15</td><td>1.46<sup>+0.65</sup><sub>-0.42</sub> ± 0.14</td><td>0.77<sup>+1.66</sup><sub>-0.77</sub> ± 0.15</td><td>1.56<sup>+1.89</sup><sub>-1.09</sub> ± 0.09</td><td>1.80<sup>+2.37</sup><sub>-1.35</sub> ± 0.10</td><td>671.203</td><td>0.08635</td><td>0.30</td><td>0.28</td><td>88</td><td>793 ± 53</td></tr><tr><td>G149.39-36.84</td><td>Eis</td><td>pHer</td><td>1.65<sup>+0.53</sup><sub>-0.66</sub> ± 0.15</td><td>3.15<sup>+1.49</sup><sub>-1.32</sub> ± 0.14</td><td>1.42<sup>+1.99</sup><sub>-1.40</sub> ± 0.15</td><td>1.46<sup>+1.94</sup><sub>-1.00</sub> ± 0.09</td><td>2.12<sup>+2.35</sup><sub>-1.63</sub> ± 0.10</td><td>657.580</td><td>0.1699</td><td>0.44</td><td>0.19</td><td>100</td><td>912 ± 85</td></tr><tr><td>G040.03+74.95</td><td>NFW</td><td>pHer</td><td>4.26<sup>+0.65</sup><sub>-0.64</sub> ± 0.15</td><td>1.64<sup>+0.38</sup><sub>-0.25</sub> ± 0.14</td><td>0.96<sup>+1.28</sup><sub>-0.86</sub> ± 0.15</td><td>1.49<sup>+1.76</sup><sub>-0.95</sub> ± 0.09</td><td>1.67<sup>+2.36</sup><sub>-1.23</sub> ± 0.10</td><td>1645.232</td><td>0.07067</td><td>4.63</td><td>0.39</td><td>241</td><td>2041 ± 66</td></tr><tr><td>G040.03+74.95a</td><td>NFW</td><td>pHer</td><td>2.76<sup>+0.49</sup><sub>-0.45</sub> ± 0.12</td><td>1.42<sup>+0.41</sup><sub>-0.37</sub> ± 0.29</td><td>1.31<sup>+1.21</sup><sub>-1.15</sub> ± 0.60</td><td>1.61<sup>+2.65</sup><sub>-1.15</sub> ± 0.11</td><td>1.50<sup>+2.58</sup><sub>-1.05</sub> ± 0.08</td><td>1186.989</td><td>0.07558</td><td>0.54</td><td>0.24</td><td>125</td><td>1180 ± 80</td></tr></tbody></table>
<table><caption>Table A.1. continued.</caption><thead><tr><th>PSZ2 name</th><th>Mass</th><th>N(R)</th><th>r<sub>200</sub><br>[Mpc]</th><th>r<sub>γ</sub><br>[Mpc]</th><th>r<sub>s</sub><br>[Mpc]</th><th>A<sub>∞</sub></th><th>A<sub>0</sub></th><th>-ln P</th><th>z</th><th>A<sup>2</sup></th><th>f<sub>sub</sub></th><th>N<sub>200</sub></th><th>σ<sub>20,200</sub><br>[km/s]</th></tr></thead><tbody><tr><td>G217.09+40.15</td><td>NFW</td><td>pHer</td><td>1.83<sup>+0.35</sup><sub>-0.32</sub> ± 0.15</td><td>1.47<sup>+0.09</sup><sub>-0.09</sub> ± 0.14</td><td>1.09<sup>+1.98</sup><sub>-1.09</sub> ± 0.15</td><td>1.40<sup>+1.73</sup><sub>-0.94</sub> ± 0.09</td><td>2.11<sup>+2.35</sup><sub>-1.61</sub> ± 0.10</td><td>1050.259</td><td>0.1382</td><td>1.42</td><td>0.33</td><td>117</td><td>893 ± 66</td></tr><tr><td>G172.74+65.30</td><td>NFW</td><td>pNFW</td><td>1.57<sup>-0.47</sup><sub>-0.42</sub> ± 0.15</td><td>0.52<sup>+0.32</sup><sub>-0.19</sub> ± 0.14</td><td>1.25<sup>+1.96</sup><sub>-1.24</sub> ± 0.15</td><td>1.28<sup>+1.83</sup><sub>-0.83</sub> ± 0.09</td><td>2.05<sup>+2.34</sup><sub>-1.57</sub> ± 0.10</td><td>577.073</td><td>0.07529</td><td>0.30</td><td>0.20</td><td>74</td><td>702 ± 50</td></tr><tr><td>G105.55+77.21</td><td>NFW</td><td>pHer</td><td>1.72<sup>+0.41</sup><sub>-0.36</sub> ± 0.15</td><td>1.34<sup>+0.58</sup><sub>-0.42</sub> ± 0.14</td><td>1.33<sup>+1.97</sup><sub>-1.32</sub> ± 0.15</td><td>1.31<sup>+1.87</sup><sub>-0.85</sub> ± 0.09</td><td>1.98<sup>+2.45</sup><sub>-1.51</sub> ± 0.10</td><td>700.214</td><td>0.07242</td><td>0.36</td><td>0.17</td><td>85</td><td>799 ± 62</td></tr><tr><td>G066.68+68.44</td><td>Bur</td><td>pHer</td><td>1.64<sup>+0.35</sup><sub>-0.29</sub> ± 0.15</td><td>1.63<sup>+0.63</sup><sub>-0.53</sub> ± 0.14</td><td>0.36<sup>+0.73</sup><sub>-0.34</sub> ± 0.15</td><td>1.79<sup>+1.82</sup><sub>-1.24</sub> ± 0.09</td><td>1.77<sup>+2.41</sup><sub>-1.52</sub> ± 0.10</td><td>845.978</td><td>0.1625</td><td>0.15</td><td>0.11</td><td>97</td><td>830 ± 60</td></tr><tr><td>G050.40+31.17</td><td>NFW</td><td>pNFW</td><td>1.89<sup>-0.43</sup><sub>-0.43</sub> ± 0.07</td><td>0.63<sup>+0.37</sup><sub>-0.24</sub> ± 0.52</td><td>1.18<sup>+1.83</sup><sub>-1.14</sub> ± 0.46</td><td>1.33<sup>+1.92</sup><sub>-0.87</sub> ± 0.08</td><td>1.69<sup>+2.27</sup><sub>-1.25</sub> ± 0.16</td><td>837.047</td><td>0.1605</td><td>0.36</td><td>0.07</td><td>95</td><td>932 ± 59</td></tr><tr><td>G057.78+52.32</td><td>Eis</td><td>pHer</td><td>1.65<sup>+0.29</sup><sub>-0.29</sub> ± 0.07</td><td>0.94<sup>+0.29</sup><sub>-0.25</sub> ± 0.52</td><td>0.32<sup>+0.55</sup><sub>-0.30</sub> ± 0.46</td><td>1.79<sup>+1.83</sup><sub>-1.25</sub> ± 0.08</td><td>2.01<sup>+2.39</sup><sub>-1.55</sub> ± 0.16</td><td>848.210</td><td>0.06628</td><td>0.23</td><td>0.14</td><td>104</td><td>792 ± 57</td></tr><tr><td>G040.58+77.12</td><td>NFW</td><td>pHer</td><td>1.43<sup>+0.48</sup><sub>-0.51</sub> ± 0.07</td><td>2.52<sup>+1.64</sup><sub>-1.52</sub> ± 0.52</td><td>1.58<sup>+1.92</sup><sub>-1.54</sub> ± 0.46</td><td>1.21<sup>+1.77</sup><sub>-0.75</sub> ± 0.08</td><td>2.20<sup>+2.32</sup><sub>-1.64</sub> ± 0.16</td><td>442.840</td><td>0.0754</td><td>0.34</td><td>0.00</td><td>60</td><td>697 ± 57</td></tr><tr><td>G028.63+50.15</td><td>Her</td><td>pHer</td><td>1.84<sup>+0.49</sup><sub>-0.44</sub> ± 0.07</td><td>0.73<sup>+0.36</sup><sub>-0.31</sub> ± 0.52</td><td>1.55<sup>+1.90</sup><sub>-1.49</sub> ± 0.46</td><td>1.43<sup>+1.94</sup><sub>-0.97</sub> ± 0.08</td><td>1.95<sup>+2.45</sup><sub>-1.49</sub> ± 0.16</td><td>416.306</td><td>0.09083</td><td>0.32</td><td>0.06</td><td>50</td><td>784 ± 76</td></tr><tr><td>G186.37+37.26</td><td>Bur</td><td>pNFW</td><td>2.25<sup>+0.28</sup><sub>-0.11</sub> ± 0.18</td><td>0.87<sup>+0.24</sup><sub>-0.14</sub> ± 0.14</td><td>0.36<sup>+0.38</sup><sub>-0.28</sub> ± 0.69</td><td>1.79<sup>+1.78</sup><sub>-1.17</sub> ± 0.11</td><td>1.32<sup>+2.12</sup><sub>-0.87</sub> ± 0.09</td><td>2213.302</td><td>0.2818</td><td>0.48</td><td>0.21</td><td>243</td><td>1232 ± 55</td></tr><tr><td>G228.16+75.20</td><td>Bur</td><td>pHer</td><td>2.54<sup>+0.34</sup><sub>-0.31</sub> ± 0.18</td><td>0.75<sup>+0.09</sup><sub>-0.19</sub> ± 0.14</td><td>0.29<sup>+0.24</sup><sub>-0.19</sub> ± 0.69</td><td>1.93<sup>+1.79</sup><sub>-1.42</sub> ± 0.11</td><td>0.95<sup>+1.17</sup><sub>-0.48</sub> ± 0.09</td><td>2705.849</td><td>0.5423</td><td>0.72</td><td>0.31</td><td>293</td><td>1467 ± 68</td></tr><tr><td>G286.98+32.90</td><td>Her</td><td>pNFW</td><td>3.21<sup>+0.77</sup><sub>-0.70</sub> ± 0.18</td><td>0.31<sup>+0.23</sup><sub>-0.18</sub> ± 0.14</td><td>1.24<sup>+1.88</sup><sub>-1.22</sub> ± 0.69</td><td>1.59<sup>+1.92</sup><sub>-1.12</sub> ± 0.11</td><td>1.95<sup>+2.41</sup><sub>-1.50</sub> ± 0.09</td><td>619.176</td><td>0.386</td><td>0.16</td><td>0.00</td><td>73</td><td>1790 ± 160</td></tr><tr><td>G028.89+60.13</td><td>Her</td><td>pNFW</td><td>1.70<sup>+0.32</sup><sub>-0.30</sub> ± 0.18</td><td>0.60<sup>+0.30</sup><sub>-0.23</sub> ± 0.14</td><td>1.35<sup>+1.85</sup><sub>-1.22</sub> ± 0.69</td><td>1.46<sup>+1.87</sup><sub>-1.00</sub> ± 0.11</td><td>1.82<sup>+2.49</sup><sub>-1.37</sub> ± 0.09</td><td>1001.433</td><td>0.1519</td><td>0.19</td><td>0.09</td><td>120</td><td>812 ± 59</td></tr><tr><td>G087.03-57.37</td><td>Bur</td><td>pNFW</td><td>1.86<sup>+0.69</sup><sub>-0.30</sub> ± 0.06</td><td>0.72<sup>+0.39</sup><sub>-0.30</sub> ± 0.50</td><td>1.28<sup>+1.53</sup><sub>-1.12</sub> ± 0.73</td><td>1.43<sup>+1.87</sup><sub>-0.97</sub> ± 0.09</td><td>2.08<sup>+2.29</sup><sub>-1.58</sub> ± 0.10</td><td>1085.206</td><td>0.2761</td><td>0.33</td><td>0.32</td><td>97</td><td>1055 ± 84</td></tr><tr><td>G071.63+29.78</td><td>Bur</td><td>pHer</td><td>1.82<sup>+0.89</sup><sub>-1.12</sub> ± 0.06</td><td>1.50<sup>+0.49</sup><sub>-0.49</sub> ± 0.50</td><td>1.97<sup>+1.63</sup><sub>-1.46</sub> ± 0.73</td><td>1.22<sup>+1.95</sup><sub>-0.76</sub> ± 0.09</td><td>1.88<sup>+2.35</sup><sub>-1.34</sub> ± 0.10</td><td>807.446</td><td>0.1575</td><td>0.41</td><td>0.33</td><td>91</td><td>851 ± 58</td></tr><tr><td>G073.97-27.82</td><td>NFW</td><td>pNFW</td><td>2.48<sup>+0.33</sup><sub>-0.22</sub> ± 0.06</td><td>0.72<sup>+0.20</sup><sub>-0.17</sub> ± 0.50</td><td>0.46<sup>+0.46</sup><sub>-0.34</sub> ± 0.73</td><td>1.95<sup>+1.68</sup><sub>-1.19</sub> ± 0.09</td><td>1.04<sup>+1.37</sup><sub>-0.58</sub> ± 0.10</td><td>3368.643</td><td>0.2282</td><td>0.74</td><td>0.34</td><td>353</td><td>1339 ± 52</td></tr><tr><td>G092.71+73.46</td><td>NFW</td><td>pHer</td><td>2.74<sup>+0.39</sup><sub>-0.36</sub> ± 0.06</td><td>1.95<sup>+0.37</sup><sub>-0.40</sub> ± 0.50</td><td>1.16<sup>+1.64</sup><sub>-1.04</sub> ± 0.73</td><td>1.50<sup>+1.85</sup><sub>-0.97</sub> ± 0.09</td><td>1.69<sup>+2.39</sup><sub>-1.23</sub> ± 0.10</td><td>2246.441</td><td>0.2317</td><td>0.97</td><td>0.18</td><td>261</td><td>1342 ± 51</td></tr><tr><td>G187.53+21.92</td><td>Bur</td><td>pNFW</td><td>1.63<sup>+0.33</sup><sub>-0.37</sub> ± 0.06</td><td>0.71<sup>+0.31</sup><sub>-0.23</sub> ± 0.50</td><td>0.37<sup>+0.58</sup><sub>-0.34</sub> ± 0.73</td><td>1.59<sup>+1.86</sup><sub>-1.11</sub> ± 0.09</td><td>1.65<sup>+2.55</sup><sub>-1.19</sub> ± 0.10</td><td>1407.205</td><td>0.1707</td><td>0.53</td><td>0.24</td><td>168</td><td>851 ± 48</td></tr><tr><td>G278.58+39.16</td><td>NFW</td><td>pHer</td><td>2.18<sup>+0.34</sup><sub>-0.52</sub> ± 0.06</td><td>1.75<sup>+0.35</sup><sub>-0.81</sub> ± 0.50</td><td>0.81<sup>+1.74</sup><sub>-0.81</sub> ± 0.73</td><td>1.62<sup>+1.78</sup><sub>-1.10</sub> ± 0.09</td><td>1.95<sup>+2.39</sup><sub>-1.50</sub> ± 0.10</td><td>1170.971</td><td>0.3063</td><td>0.24</td><td>0.26</td><td>130</td><td>1140 ± 70</td></tr><tr><td>G218.81+35.51</td><td>Bur</td><td>pHer</td><td>3.70<sup>+0.57</sup><sub>-0.61</sub> ± 0.06</td><td>2.24<sup>+0.38</sup><sub>-0.37</sub> ± 0.50</td><td>1.85<sup>+1.86</sup><sub>-1.88</sub> ± 0.73</td><td>2.50<sup>+1.54</sup><sub>-2.03</sub> ± 0.09</td><td>3.63<sup>+1.50</sup><sub>-2.20</sub> ± 0.10</td><td>4820.696</td><td>0.1663</td><td>22.14</td><td>0.32</td><td>486</td><td>1821 ± 84</td></tr><tr><td>G218.81+35.51 a</td><td>Bur</td><td>pNFW</td><td>1.03<sup>+0.42</sup><sub>-0.48</sub> ± 0.06</td><td>3.42<sup>+0.36</sup><sub>-0.55</sub> ± 0.40</td><td>1.46<sup>+1.92</sup><sub>-1.37</sub> ± 0.65</td><td>1.87<sup>+2.62</sup><sub>-1.39</sub> ± 0.11</td><td>2.09<sup>+2.55</sup><sub>-1.59</sub> ± 0.10</td><td>1357.077</td><td>0.1643</td><td>0.9179</td><td>0.35</td><td>151</td><td>671 ± 40</td></tr><tr><td>G218.81+35.51 b</td><td>Bur</td><td>pNFW</td><td>1.52<sup>+0.46</sup><sub>-0.43</sub> ± 0.07</td><td>0.62<sup>+0.23</sup><sub>-0.25</sub> ± 0.22</td><td>0.49<sup>+0.80</sup><sub>-0.44</sub> ± 0.72</td><td>2.08<sup>+2.61</sup><sub>-1.61</sub> ± 0.09</td><td>1.36<sup>+2.71</sup><sub>-0.90</sub> ± 0.09</td><td>843.035</td><td>0.1776</td><td>0.70</td><td>0.32</td><td>89</td><td>769 ± 45</td></tr><tr><td>G049.32+44.37</td><td>NFW</td><td>pHer</td><td>1.74<sup>+0.38</sup><sub>-0.37</sub> ± 0.06</td><td>1.57<sup>+0.56</sup><sub>-0.40</sub> ± 0.50</td><td>1.71<sup>+1.83</sup><sub>-1.51</sub> ± 0.73</td><td>1.27<sup>+1.84</sup><sub>-0.81</sub> ± 0.09</td><td>1.94<sup>+2.46</sup><sub>-1.44</sub> ± 0.10</td><td>1124.578</td><td>0.09606</td><td>0.20</td><td>0.24</td><td>134</td><td>797 ± 51</td></tr><tr><td>G046.10+27.18</td><td>Her</td><td>pHer</td><td>3.29<sup>+0.96</sup><sub>-0.92</sub> ± 0.06</td><td>0.57<sup>+0.29</sup><sub>-0.23</sub> ± 0.50</td><td>1.96<sup>+1.72</sup><sub>-1.48</sub> ± 0.73</td><td>1.28<sup>+1.93</sup><sub>-0.81</sub> ± 0.09</td><td>1.69<sup>+2.54</sup><sub>-1.24</sub> ± 0.10</td><td>560.821</td><td>0.3878</td><td>0.57</td><td>0.05</td><td>63</td><td>1600 ± 250</td></tr><tr><td>G262.73+40.92</td><td>NFW</td><td>pHer</td><td>2.79<sup>+0.96</sup><sub>-0.86</sub> ± 0.06</td><td>0.59<sup>+0.22</sup><sub>-0.23</sub> ± 0.50</td><td>1.39<sup>+1.92</sup><sub>-1.34</sub> ± 0.73</td><td>1.24<sup>+2.06</sup><sub>-0.78</sub> ± 0.09</td><td>1.98<sup>+2.39</sup><sub>-1.48</sub> ± 0.10</td><td>530.355</td><td>0.4218</td><td>0.44</td><td>0.18</td><td>67</td><td>1191 ± 92</td></tr><tr><td>G206.45+13.89</td><td>Bur</td><td>pHer</td><td>3.54<sup>+1.10</sup><sub>-0.91</sub> ± 0.06</td><td>0.98<sup>+0.47</sup><sub>-0.39</sub> ± 0.50</td><td>1.12<sup>+1.57</sup><sub>-0.90</sub> ± 0.73</td><td>1.51<sup>+1.96</sup><sub>-0.90</sub> ± 0.09</td><td>1.35<sup>+2.37</sup><sub>-0.90</sub> ± 0.10</td><td>553.161</td><td>0.4054</td><td>0.37</td><td>0.07</td><td>59</td><td>1520 ± 140</td></tr><tr><td>G083.29-31.03</td><td>NFW</td><td>pHer</td><td>2.44<sup>+0.64</sup><sub>-0.57</sub> ± 0.20</td><td>0.74<sup>+0.14</sup><sub>-0.18</sub> ± 0.27</td><td>1.56<sup>+1.90</sup><sub>-1.48</sub> ± 0.51</td><td>1.32<sup>+1.90</sup><sub>-0.86</sub> ± 0.14</td><td>1.86<sup>+2.42</sup><sub>-1.40</sub> ± 0.19</td><td>546.380</td><td>0.4102</td><td>0.40</td><td>0.00</td><td>21</td><td>1170 ± 170</td></tr><tr><td>G201.50-27.31</td><td>NFW</td><td>pHer</td><td>2.35<sup>+0.31</sup><sub>-0.29</sub> ± 0.06</td><td>1.31<sup>+0.31</sup><sub>-0.25</sub> ± 0.50</td><td>1.13<sup>+1.55</sup><sub>-0.98</sub> ± 0.73</td><td>1.44<sup>+1.72</sup><sub>-0.90</sub> ± 0.09</td><td>1.65<sup>+2.37</sup><sub>-1.19</sub> ± 0.10</td><td>2281.875</td><td>0.539</td><td>0.93</td><td>0.31</td><td>241</td><td>1330 ± 50</td></tr><tr><td>G057.25-45.34</td><td>Bur</td><td>pHer</td><td>2.52<sup>+1.04</sup><sub>-1.77</sub> ± 0.14</td><td>4.72<sup>+3.45</sup><sub>-2.61</sub> ± 1.22</td><td>1.15<sup>+2.19</sup><sub>-1.14</sub> ± 0.20</td><td>1.65<sup>+1.86</sup><sub>-1.19</sub> ± 0.08</td><td>1.96<sup>+2.45</sup><sub>-1.50</sub> ± 0.06</td><td>479.187</td><td>0.3977</td><td>0.41</td><td>0.21</td><td>54</td><td>1540 ± 160</td></tr><tr><td>G111.61-45.71</td><td>Her</td><td>pHer</td><td>2.66<sup>+0.45</sup><sub>-0.42</sub> ± 0.06</td><td>0.98<sup>+0.21</sup><sub>-0.16</sub> ± 0.50</td><td>1.40<sup>+1.66</sup><sub>-1.16</sub> ± 0.73</td><td>1.37<sup>+1.97</sup><sub>-0.91</sub> ± 0.09</td><td>1.54<sup>+2.45</sup><sub>-1.10</sub> ± 0.10</td><td>2157.449</td><td>0.5477</td><td>0.76</td><td>0.23</td><td>255</td><td>1765 ± 84</td></tr></tbody></table>

Pizzuti, L., et al.: A&A, 699, A88 (2025)
**Appendix B: Disturbed clusters**

In Fig. B.1 we show the PPSs, the histograms of the LOS velocities, and the 2D spatial distributions for the four clusters characterised by a value of $\delta\eta > 0.8$, to highlight the evident disturbed dynamical state. From top to bottom: PSZ2G056.77+36.32 (skewness $s = -0.01$, kurtosis excess $k-3 = -0.54$), PSZ2G075.71+13.51 ($s = -0.11$, $k-3 = -0.43$), PSZ2G008.31-64.74 ($s = -0.15$, $k-3 = -0.51$) and PSZ2G040.03+74.95 ($s = -0.13$, $k-3 = -1.01$).

<figure><img src="image_15.png" /><figcaption>**Fig. B.1.** Each of the four clusters in the sample with a logarithmic scatter $\delta\eta > 0.8$ dex. Left: PPS. Central: LOS velocity distribution. Right: Projected positions with respect to the cluster centre.</figcaption></figure>

**Appendix C: Stacked projected phase spaces in bins of $M_{500}^{SZ}$**

In Fig. C.1 we plot the PPS of the five stacked clusters obtained by collecting member galaxies in bins of SZ masses, for T1 (top plots, blue points) and T2 (bottom plots, green points).
<figure><img src="image_16.png" /><figcaption>**Fig. C.1.** PPS parameters (radius, *R*, and radial velocity, *v*) for the SZ-stacked clusters. Top, blue points: T1. Bottom, green points: T2. The vertical dashed lines indicate the average *r*<sub>500</sub> of all the clusters in each stacking.</figcaption></figure>

## Appendix D: GTC cluster data

Here we report an extract of the member galaxies with measured spectroscopic redshifts for the three clusters observed with the OSIRIS spectrograph at the GTC. These three catalogues are organised as follow: first column corresponds to an ID label; the second and third columns indicate the right ascension and declination coordinates in J2000 epoch (RA and Dec); the fourth is the measured redshift and the last column is the errors on the redshift. The full catalogues are available upon request to the authors.

<table><caption>**Table D.1.** Extract of the galaxies with measured spectroscopic redshifts in the field of PSZ2G083.29-31.03.</caption><thead><tr><th>ID</th><th>RA<br>(deg; J2000)</th><th>Dec</th><th>z</th><th>δz</th></tr></thead><tbody><tr><td>1</td><td>337.14050</td><td>20.62120</td><td>0.41158</td><td>0.00009</td></tr><tr><td>2</td><td>337.13602</td><td>20.62673</td><td>0.42092</td><td>0.00013</td></tr><tr><td>3</td><td>337.15583</td><td>20.61639</td><td>0.40267</td><td>0.00013</td></tr><tr><td>4</td><td>337.13124</td><td>20.59847</td><td>0.41477</td><td>0.00013</td></tr><tr><td>5</td><td>337.15079</td><td>20.59610</td><td>0.42112</td><td>0.00009</td></tr><tr><td>6</td><td>337.11054</td><td>20.61091</td><td>0.41220</td><td>0.00011</td></tr><tr><td>7</td><td>337.10786</td><td>20.61638</td><td>0.41187</td><td>0.00008</td></tr><tr><td>8</td><td>337.13583</td><td>20.59031</td><td>0.41003</td><td>0.00008</td></tr><tr><td>9</td><td>337.17094</td><td>20.63067</td><td>0.41885</td><td>0.00013</td></tr><tr><td>10</td><td>337.17612</td><td>20.60638</td><td>0.41726</td><td>0.00010</td></tr></tbody></table>
<table><caption>**Table D.2.** Extract of the galaxies with measured spectroscopic redshifts in the field of PSZ2G113.91-37.01.</caption><thead><tr><th>ID</th><th>RA</th><th>Dec</th><th>z</th><th>δz</th></tr><tr><td></td><td colspan="4">(deg; J2000)</td></tr></thead><tbody><tr><td>1</td><td>4.92222</td><td>25.30204</td><td>0.36999</td><td>0.00007</td></tr><tr><td>2</td><td>4.93030</td><td>25.30678</td><td>0.35787</td><td>0.00006</td></tr><tr><td>3</td><td>4.93419</td><td>25.28389</td><td>0.36852</td><td>0.00009</td></tr><tr><td>4</td><td>4.93971</td><td>25.31805</td><td>0.36313</td><td>0.00012</td></tr><tr><td>5</td><td>4.90339</td><td>25.31256</td><td>0.36536</td><td>0.00008</td></tr><tr><td>6</td><td>4.95634</td><td>25.30796</td><td>0.37480</td><td>0.00008</td></tr><tr><td>7</td><td>4.90230</td><td>25.27815</td><td>0.36741</td><td>0.00009</td></tr><tr><td>8</td><td>4.88488</td><td>25.28588</td><td>0.35860</td><td>0.00012</td></tr><tr><td>9</td><td>4.90188</td><td>25.34320</td><td>0.37568</td><td>0.00008</td></tr><tr><td>10</td><td>4.91762</td><td>25.25141</td><td>0.36569</td><td>0.00009</td></tr></tbody></table>

<table><caption>**Table D.3.** Extract of the galaxies with measured spectroscopic redshifts in the field of PSZ2G046.10+27.18.</caption><thead><tr><th>ID</th><th>RA</th><th>Dec</th><th>z</th><th>δz</th></tr><tr><td></td><td colspan="4">(deg; J2000)</td></tr></thead><tbody><tr><td>1</td><td>262.93623</td><td>22.80826</td><td>0.38764</td><td>0.00037</td></tr><tr><td>2</td><td>262.90941</td><td>22.81799</td><td>0.38311</td><td>0.00014</td></tr><tr><td>3</td><td>262.91685</td><td>22.81963</td><td>0.37632</td><td>0.00010</td></tr><tr><td>4</td><td>262.92559</td><td>22.82227</td><td>0.39349</td><td>0.00021</td></tr><tr><td>5</td><td>262.90121</td><td>22.82538</td><td>0.37398</td><td>0.00027</td></tr><tr><td>6</td><td>262.90976</td><td>22.82926</td><td>0.37562</td><td>0.00017</td></tr><tr><td>7</td><td>262.90436</td><td>22.83141</td><td>0.38326</td><td>0.00014</td></tr><tr><td>8</td><td>262.89507</td><td>22.83543</td><td>0.38148</td><td>0.00020</td></tr><tr><td>9</td><td>262.93955</td><td>22.83818</td><td>0.37660</td><td>0.00033</td></tr><tr><td>10</td><td>262.89412</td><td>22.83914</td><td>0.38374</td><td>0.00033</td></tr></tbody></table>

## Appendix E: Cluster members and velocity dispersion

In Fig. E.1 we show for each cluster the number of cluster members $N_{200}$ that remains after the CLEAN procedure within $R = r_{ap,200}$ (see also Fig. 4 of Sereno et al. 2025), plotted against the aperture velocity dispersion $\sigma_{ap,200}$. It is worth noticing that the majority of clusters have more than 100 members, making the CHEX-MATE catalogue a very good sample for kinematical analyses of galaxy cluster.

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

**Fig. E.1.** Number of members within $R = r_{ap,200}$ after the application of the CLEAN algorithm, in relation to the aperture velocity dispersion $\sigma_{ap,200}$.