# Linking a universal gas density profile to the core-excised X-ray luminosity in galaxy clusters up to z 1.1

G.W. Pratt, M. Arnaud, B.J. Maughan, J.-B. Melin

► To cite this version:

G.W. Pratt, M. Arnaud, B.J. Maughan, J.-B. Melin. Linking a universal gas density profile to the core-excised X-ray luminosity in galaxy clusters up to z 1.1. Astronomy and Astrophysics - A&A, 2022, 665, pp.A24. 10.1051/0004-6361/202243074. hal-03706338

HAL Id: hal-03706338

https://hal.science/hal-03706338

Submitted on 11 Apr 2023

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Linking a universal gas density profile to the core-excised X-ray luminosity in galaxy clusters up to $z \sim 1.1$

G. W. Pratt<fnref n="1" />, M. Arnaud<fnref n="1" />, B. J. Maughan<fnref n="2" />, and J.-B. Melin<fnref n="3" />

Received 10 January 2022 / Accepted 26 May 2022

## ABSTRACT

We investigate the regularity of galaxy cluster gas density profiles and the link to the relation between core-excised luminosity, $L_{Xc}$, and mass from the $Y_X$ proxy, $M_{Y_X}$, for 93 objects selected through their Sunyaev-Zeldovich effect (SZE) signal. The sample spans a mass range of $M_{500} = [0.5-20] \times 10^{14} M_\odot$, and lies at redshifts $0.05 < z < 1.13$. To investigate differences in X-ray and SZE selection, we compare to the local X-ray-selected REXCESS sample. Using XMM-Newton observations, we derive an average intra-cluster medium (ICM) density profile for the SZE-selected systems and determine its scaling with mass and redshift. This average profile exhibits an evolution that is slightly stronger than self-similar ($\alpha_Z = 2.09 \pm 0.02$), and a significant dependence on mass ($\alpha_M = 0.22 \pm 0.01$). Deviations from this average scaling with radius, which we quantify, indicate different evolution for the core regions as compared to the bulk. We measure the radial variation of the intrinsic scatter in scaled density profiles, finding a minimum of ~20% at $R \sim [0.5-0.7] R_{500}$ and a value of ~40% at $R_{500}$; moreover, the scatter evolves slightly with redshift. The average profile of the SZE-selected systems adequately describes the X-ray-selected systems and their intrinsic scatter at low redshift, except in the very central regions. We examine the evolution of the scaled core properties over time, which are positively skewed at later times, suggesting an increased incidence of centrally peaked objects at lower redshifts. The relation between core-excised luminosity, $L_{Xc}$, and mass is extremely tight, with a measured logarithmic intrinsic scatter of $\sigma_{\ln L_{Xc}|M_{Y_X}} \sim 0.13$. Using extensive simulations, we investigate the impact of selection effects, intrinsic scatter, and covariance between quantities on this relation. The slope is insensitive to selection and intrinsic scatter between quantities; however, the scatter is very dependent on the covariance between $L_{Xc}$ and $Y_X$. Accounting for our use of the $Y_X$ proxy to determine the mass, for observationally motivated values of covariance we estimate an upper limit to the logarithmic intrinsic scatter with respect to the true mass of $\sigma_{\ln L_{Xc}|M} \sim 0.22$. We explicitly illustrate the connection between the scatter in density profiles and that in the $L_{Xc} - M$ relation. Our results are consistent with the overall conclusion that the ICM bulk evolves approximately self-similarly, with the core regions evolving separately. They indicate a systematic variation of the gas content with mass. They also suggest that the core-excised X-ray luminosity, $L_{Xc}$, has a tight and well-understood relation to the underlying mass.

**Key words.** X-rays: galaxies: clusters – galaxies: clusters: intracluster medium – large-scale structure of Universe

## 1. Introduction

In a Λ cold dark matter Universe, halo assembly is driven by the hierarchical gravitational collapse of the dominant dark matter component. To first order, this process is self-similar and scale-free because gravity has no characteristic scale. If the baryon content remains constant, power-law relations link the baryonic observable properties – such as X-ray luminosity $L_X$, or Sunyaev-Zeldovich effect (SZE) signal $Y_{SZ}$, or total optical richness $\Lambda$ – to the cluster mass and redshift. Each of these observables is sensitive to a different underlying intrinsic physical characteristic (e.g., the distribution of gas or the number of red sequence galaxies above a given threshold). However, the detection of a baryonic observable also depends on the intrinsic properties of the signal itself. For example, while the SZE signal is proportional to the gas density, the X-ray emission is proportional to the square of the gas density. This means that X-ray measurements are very sensitive to the physical conditions in the core regions, while SZE measurements are much less so.

The X-ray luminosity, $L_X$, is an attractive quantity because it can be measured from very few source counts once the

redshift is known. For a virialised galaxy cluster where the intra-cluster medium (ICM) is in hydrostatic equilibrium, $L_X$ depends only on the halo mass, $M$, redshift, $z$, and the distribution of the gas in the dark matter potential. Power-law $L_X - M$ relations have indeed been observed (e.g., Maughan 2007; Rykoff et al. 2008; Zhang et al. 2008; Connor et al. 2014; Pratt et al. 2009; Vikhlinin et al. 2009; Schellenberger & Reiprich 2017; Lovisari et al. 2020). However, these relations exhibit a large intrinsic scatter (~40 per cent), linked to the presence of cool cores and merging activity (e.g., Pratt et al. 2009). Exclusion of the core regions, by measuring $L_X$ in an annulus excluding the cluster centre, significantly reduces the intrinsic scatter (Fabian et al. 1994; Maughan 2007; Pratt et al. 2009; Mantz et al. 2018).

Radial gas density profiles have been a key means of obtaining information about the ICM since the advent of X-ray imaging. The high spatial resolution observations afforded by XMM-Newton and Chandra observations have revealed the complexity of both the core regions, which are strongly affected by non-gravitational processes (Croston et al. 2008; Pratt et al. 2010), and the outskirts, where the gas distribution becomes

<footnote n="1">¹ Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM de Paris-Saclay, 91191 Gif-sur-Yvette, France
e-mail: gabriel.pratt@cea.fr</footnote>

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

<footnote n="3">³ IRFU, CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Sample properties. Left: redshift-mass distribution of the clusters used in this paper. The SZE-selected clusters comprise a subset of 44 systems from the Planck Early SZ sample (Planck Collaboration VIII 2011) at $z < 0.5$ and a further 49 clusters at $z > 0.5$. REXCESS (Böhringer et al. 2007) is an X-ray-selected sample of 31 objects at $z < 0.25$. Right: stacked histogram of the mass distribution. The REXCESS sample has a lower median mass than the SZE-selected samples.</figcaption></figure>

progressively more inhomogeneous (e.g., Eckert et al. 2015). Key open questions are how the ICM evolves over time in the dark matter potential, and how this connects to the formation and evolution of cool cores.

The advent of SZE surveys has resulted in the detection of large numbers of clusters at $z > 0.5$ (Hasselfield et al. 2013; Bleem et al. 2015; Planck Collaboration XXVII 2016; Hilton et al. 2021), extending the redshift leverage for studies of how the population changes over time. The suggestion that X-ray-selected and SZE-selected samples may not have the same distribution of dynamical states (e.g., Planck Collaboration IX 2011; Rossetti et al. 2016; Andrade-Santos et al. 2017; Lovisari et al. 2017) has prompted examination of the relationship between the baryon signatures and the true underlying cluster population. Indeed, the dynamical state may well be as fundamental a characteristic as the mass or the redshift (Bartalucci et al. 2019). At the same time, McDonald et al. (2017) found that the new SZE-selected samples suggest that cool cores are in place very early in the history of a cluster, and have not changed in size, density, or total mass up to the present. These authors further suggest that much of what was thought to be cool core growth over time is in fact due to the self-similar evolution of the cluster bulk around this static core.

Here we investigate the properties of the gas density profiles and the core-excised X-ray luminosity, $L_{Xc}$, of 31 X-ray-selected clusters at $z < 0.2$ and 93 SZE-selected clusters at $z < 1.13$. We describe a universal gas density profile for the SZE-selected objects and quantify its variation with redshift and mass. We quantify the radial variation in scaled profiles with respect to the best-fitting evolution and the best-fitting mass dependence. Outside the cores, the median scaled gas density profiles are remarkably similar, showing no dependence on selection. We obtain the radial variation of the intrinsic scatter in scaled profiles and investigate the evolution of this scatter with redshift. We find that $L_{Xc}$, measured in the [0.15–1] $R_{500}$ region<fnref n="1" />, is an extremely well-behaved mass proxy that does not depend on cluster selection, and shows little evolution beyond self-similar in the broad redshift range of the sample.

Throughout the paper we assume a flat $\Lambda$CDM cosmology with $\Omega_m = 0.3$, $\Omega_\Lambda = 0.7$, and $H_0 = 70 \text{ km s}^{-1} \text{ Mpc}^{-1}$. The Sunyaev-Zeldovich flux in units of square arcminutes is denoted $Y_{SZ}$; the quantity $D_A^2 Y_{SZ}$, in units of square megaparsecs, is the spherically integrated Compton parameter within $R_{500}$, where $D_A$ is the angular diameter distance of the cluster. Unless stated otherwise, logarithmic quantities, including scatter, are given to base e, and uncertainties are quoted at the 68% confidence level.

## 2. Data and analysis

### 2.1. Dataset

The dataset consists of 31 X-ray-selected clusters at $0.05 < z < 0.2$ and 93 SZE-selected objects at $0.08 < z < 1.13$, with six systems in common. The local X-ray-selected dataset is REXCESS (Böhringer et al. 2007; Croston et al. 2008; Pratt et al. 2009). The SZE-selected systems consist of a local sample comprising a subset of 44 objects at $z < 0.5$ from the Planck Early SZ sample (ESZ; Planck Collaboration VIII 2011), which Bartalucci et al. (2019) show is representative of the full ESZ; these were complemented by a further 49 distant objects observed by XMM-Newton in a series of three large programmes obtained as part of the M2C project. These cover the redshift ranges $0.5 < z < 0.7$ (LP1, ID 069366, 072378), $0.7 < z < 0.9$ (LP2, ID 078388), and $z > 0.9$ (LP3, ID 074440). The LP1 sample was selected from objects detected at a signal-to-noise $S/R > 4$ in the Planck SZ catalogue, and confirmed by Autumn 2011 to be at $z > 0.5$. The LP2 sample consists of clusters with estimated masses $M_{500} > 5 \times 10^{14} \text{ M}_\odot$ at $0.7 < z < 0.9$ in the PSZ2 catalogue. The LP3 sample is derived from the five highest-mass objects at $z > 0.9$ from the combined Planck and SPT catalogues (Bartalucci et al. 2017, 2018). Full observation details for the sample can be found in Table A.1.

The left-hand panel of Fig. 1 shows the distribution of the 118 clusters in the redshift-mass plane. Here and in the following, we group the SZE-selected systems into three sub-samples in the redshift ranges $z < 0.3$ (blue), $0.3 \le z < 0.6$ (light green), and $z \ge 0.6$ (dark green), containing 37, 27, and 29 objects, respectively. The right-hand panel of Fig. 1 shows the stacked mass histogram. This plot makes clear that the REXCESS sample

<footnote n="1">¹ $R_\Delta$ is the radius within which the mean density is $\Delta$ times the critical density at the redshift of the object.</footnote>
(light blue) has a lower median mass ($M_{500} = 2.7 \times 10^{14} M_{\odot}$) than any of the SZE-selected sub-samples ($M_{500} = 6.4, 7.9$, and $5.0 \times 10^{14} M_{\odot}$, respectively). The LP2 sub-sample is subject to significant Eddington bias in the Planck signal, leading in the most extreme case to an estimated mass of only $6.5 \times 10^{13} M_{\odot}$ for PSZ2 G208.57–44.31. We show below that this has a negligible effect on our results.

## 2.2. Analysis

As our aim was to compare the SZE-selected clusters to the low-redshift X-ray-selected REXCESS systems, we followed the X-ray data reduction and analysis procedures described in Croston et al. (2008), Pratt et al. (2009, 2010). Event files were reprocessed with the XMM-Newton Science analysis System v15 and associated calibration files. Standard filtering for clean events (PATTERN < 4 and < 13 for MOS1/2 and pn detectors, respectively, and FLAG = 0) and soft proton flares was applied. The instrumental and particle background was obtained from custom stacked, recast data files derived from observations obtained with the filter wheel in the CLOSED position (FWC), renormalised using the count rate in a high energy band free of cluster emission.

Vignetting-corrected, background-subtracted [0.3–2] keV surface brightness profiles were extracted in annular bins centred on the X-ray peak. Temperature profiles were produced using the procedures described in Pratt et al. (2010). These were extracted in logarithmically spaced annular bins centred on the X-ray peak, with a binning of $R_{out}/R_{in} = 1.33–1.5$ depending on data quality. After subtraction of the FWC spectra, all spectra were grouped to a minimum of 25 counts per bin. The FWC-subtracted spectrum of the region external to the cluster was fitted with a model consisting of two MEKAL components plus an absorbed power law with a fixed slope of $\Gamma = 1.4$. The spectra were fitted in the [0.5–10] keV range using $\chi^2$ statistics, excluding the [1.4–1.6] keV band (due to the Al line in all three detectors), and, in the pn, the [7.45–9.0] keV band (due to the strong Cu line complex). In these fits the MEKAL models were unabsorbed and have solar abundances, and the temperature and normalisations are free parameters; the powerlaw component is absorbed by the Galactic absorption. Since it has a fixed slope, only its normalisation is an additional free parameter in the fit. This best-fitting model was added as an extra component to the annular spectral fits, with its normalisation rescaled to the ratio of the areas of the extraction regions (corrected for bad pixels, chip gaps, etc). In the annular spectral fits, the temperature and metallicity of the cluster component were left free, and the absorption was fixed to the HI value (Kalberla et al. 2005). The metallicity was fixed to a value of $Z = 0.3 Z_{\odot}$ when its relative uncertainty exceeded 30%.

### 2.2.1. Luminosity

The core-excised X-ray luminosity $L_{Xc}$ was measured in the [0.15–1] $R_{500}$ region for all objects. Here the only change with respect to the analysis in Pratt et al. (2009) was the use of an updated $M_{500} - Y_X$ relation from Arnaud et al. (2010) to estimate the relevant masses $M_{500}$ and scaled apertures $R_{500}$. We show below that this change has a negligible impact on the results. The core-excised X-ray luminosity $L_{Xc}$ was calculated both in the bolometric ([0.01–100] keV) and soft ([0.5–2] keV) bands for comparison to previous work. As in Pratt et al. (2009), the luminosities were calculated from the [0.3–2] keV band surface brightness profile count rates, using the best-fitting spectral

model estimated in the [0.15–1] $R_{500}$ aperture to convert from count rates to luminosity. In cases where the surface brightness profile did not extend to $R_{500}$ (seven systems), we extrapolated using a power law with a slope measured from the data at large radius. Errors on $L_{Xc}$ take into account the uncertainties in the spectral model, the count rates, and the value of $R_{500}$, and were estimated from Monte Carlo realisations in which the luminosity calculation was derived for 100 surface brightness profiles, the profiles and $R_{500}$ values each being randomised according to the observed uncertainties. Once obtained, the luminosities were further corrected for point spread function (PSF) effects by calculating the ratio of the observed to PSF-corrected count rates in each aperture (see below). Table B.1 lists the resulting masses and luminosities for the present sample, in addition to the name, redshift, RA, and Dec.

### 2.2.2. Density profiles

The vignetting-corrected background-subtracted [0.3–2] keV surface brightness profiles were used to obtain the deprojected, PSF-corrected density profiles using the regularised, non-parametric technique described in Croston et al. (2006), and applied to the REXCESS sample in Croston et al. (2008). The surface brightness profiles were converted to gas density by calculating an emissivity profile $\Lambda(\theta)$ in XSPEC, taking into account the absorption and instrumental response, and using a parameterised model of the projected temperature and abundance profiles (see e.g., Pratt & Arnaud 2003). The Croston et al. (2006) method uses the parametric PSF model of Ghizzardi (2001) as a function of the energy and angular offsets, the parameters of which can be found in EPIC-MCT-TN-011² and EPIC-MCT-TN-012³. In Bartalucci et al. (2017), the deprojected density profiles from XMM-Newton observations of a number of clusters obtained using this method were compared to Chandra observations, for which the PSF can be neglected. It was shown that the results obtained with the Croston et al. (2006) method reproduced the deprojected Chandra density profiles accurately down to an effective resolution limit of ~5 arcseconds (Fig. 6 of Bartalucci et al. 2017). The gas density $\rho_{gas} = n_e \times (\mu_e m_p)$, where $n_e$ is the electron density measured in X-rays, $m_p$ is the proton mass, and $\mu_e = 1.148$ is the mean molecular weight per free electron:

$$ \rho_{\text{gas}} = 1.92 \times 10^{-24} \left( \frac{n_e}{\text{cm}^{-3}} \right) \text{g cm}^{-3}. \quad (1) $$

## 3. Gas density profiles

### 3.1. Model

In the self-similar model, a cluster can be completely defined by only two parameters: its mass, $M_V$, and its redshift, $z$. A fundamental property of this model is the cluster overdensity, $\Delta$, with respect to the reference density of the Universe $\rho_{\text{Uni}}(z)$, from which the virial mass $M_V = \Delta\rho_{\text{Uni}}(z)(4\pi/3)R_V^3$, and radius $R_V$, can thereafter be defined. Cluster profiles then exhibit a universal form when the radii are scaled to $R_V$.

The original, and simplest, self-similar model concerns top-hat spherical collapse in the Standard CDM ($\Omega = 1$) cosmology. Here a cluster at redshift $z$ is represented by a spherical

<footnote n="2">² http://www.iasf-milano.inaf.it/~simona/pub/EPIC-MCT/EPIC-MCT-TN-011.pdf</footnote>

<footnote n="3">³ http://www.iasf-milano.inaf.it/~simona/pub/EPIC-MCT/EPIC-MCT-TN-012.pdf</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 2. Deprojected, PSF-corrected density profiles for 118 galaxy clusters, normalised by the critical density $\rho_{crit}$ and $R_{500}$. Fully self-similar clusters would trace the same locus in this plot. The profiles are colour-coded by mass $M_{500}$ in the left-hand panel, and by redshift $z$ in the right-hand panel. There are clear trends with respect to both quantities.</figcaption></figure>

perturbation that has just collapsed, with $\rho_{Univ}(z)$ being the critical density of the Universe, $\rho_c = 3H^2(z)/(8\pi G)$, and $\Delta = 178$. Of course, the hierarchical formation of structure in a $\Lambda$CDM Universe is a very complex dynamical process: objects continuously accrete matter along large-scale filaments, and there is no strict boundary that would separate a virialised region from the infall zone. The definitions of the mass and the corresponding overdensity are therefore ambiguous, as is the choice of $\rho_{Univ}(z)$, as one can use either the critical density or the mean density (see Voit 2005, for a review).

Using numerical simulations, Lau et al. (2015) showed that the structure of the inner part of clusters that is typically covered by X-ray observations is more self-similar when scaling by fixed overdensities with respect to the critical density $\rho_c(z)$. The zone in question corresponds to overdensities of $\Delta \gtrsim 200$. As a scaling radius, we therefore chose an $R_{500}$ corresponding to $\Delta = 500$, the radius within which the mean matter density is $\rho_{500} = 500\rho_c(z)$. The corresponding total mass within this radius, $M_{500}$, is

$$M_{500} = (4\pi/3) R_{500}^3 \rho_{500}, \qquad (2)$$

with

$$\rho_{500} = 9.205 \times 10^{-30} E(z)^2 \text{ g cm}^{-3}, \qquad (3)$$

and where $E(z) = [\Omega_m(1+z)^3 + \Omega_A]^{1/2}$ is the evolution of the Hubble parameter with redshift in a flat cosmology. The scaled gas density profile expressed as a function of scaled radius $x$ is then

$$\rho(x) = \frac{\rho_{\text{gas}}(R)}{\rho_{500}} \quad \text{where} \quad x = \frac{R}{R_{500}}. \qquad (4)$$

In the self-similar model, $\rho(x)$ follows a universal shape and its normalisation is independent of mass and redshift. In such a case $\langle\rho_{gas}\rangle \propto \rho_{500}$, where the angle brackets denote the average within $R_{500}$, as expected for gas evolution purely driven by gravitation. Figure 2 shows the scaled density profiles of all 118 systems. If the clusters were perfectly self-similar, they would trace the same locus in this plot. This is clearly not the

case. The colour-coding by mass and by redshift highlights that at large radius, the scaled profiles of the higher-mass, higher-redshift systems lie systematically above those of lower-mass, lower redshift-objects. These trends suggest a dependence of the scaling on $M_{500}$ and/or redshift.

To better understand this dependence, we fitted the observed scaled profiles with a model consisting of a median analytical profile, the normalisation of which is allowed to vary with $z$ and $M_{500}$, with a radially varying intrinsic scatter. The median profile was expressed as

$$\rho_m(x, z, M_{500}) = A(z, M_{500})f(x), \qquad (5)$$

where $f(x)$ is the function describing the profile shape. Here we adopted a generalised Navarro-Frenk-White (GNFW) model (Nagai et al. 2007):

$$f(x) = \frac{f_0}{(x/x_s)^\alpha [1 + (x/x_s)^\gamma]^{(3\beta-\alpha)/\gamma}}, \qquad (6)$$

where $x_s$ is the scaling radius, and the parameters $(\alpha, \gamma, 3\beta)$ are the central ($x \ll x_s$), intermediate ($x \sim x_s$), and outer ($x \gg x_s$) slopes, respectively. The case $\alpha = 0$ and $\gamma = 2$ corresponds to the standard $\beta$ model (Cavaliere & Fusco-Femiano 1976), while the case $\gamma = 2$ corresponds to the AB model introduced by Pratt & Arnaud (2002). The latter was used to model the median density profile of the REXCESS sample (Piffaretti et al. 2011).

The normalisation is given by the product $f_0 A(z, M_{500})$, where $A(z, M_{500})$ describes the departure from standard self-similarity in terms of a possible mass and/or redshift dependence of the scaled gas density. For this we assumed a power-law dependence on $M_{500}$ and $E(z)$:

$$A(z, M_{500}) = E(z)^{\alpha z} \left[ \frac{M_{500}}{5 \times 10^{14} M_{\odot}} \right]^{\alpha M}. \qquad (7)$$

The standard self-similar model corresponds to $\alpha_M = 0, \alpha_z = 0$. We expect $\alpha_M > 0$, as it is well established that the gas mass fraction of local clusters decreases with decreasing mass due to non-gravitational effects (e.g., Pratt et al. 2010). The model
<figure><img src="image_4.png" /><figcaption>Fig. 3. Marginalised posterior likelihood for the parameters of the best-fitting density profile model detailed in Sect. 3.1.</figcaption></figure>

above allows us to disentangle mass dependence and possible evolution.

Equations (6) and (7) translate into a gas mass fraction within $R_{500}$, $f_{g,500}$, which varies with mass and redshift as a function of $\alpha_M$ and $\alpha_z$:

$$f_{g,500} = \frac{M_{\text{gas}}(<R_{500})}{M_{500}} = \frac{1}{3} \int_{0}^{1} \rho(x)x^2 dx \quad (8)$$

$$= f_0 A(z, M_{500}) I(x_s, \alpha, \gamma, 3\beta), \quad (9)$$

where $M_{\text{gas}}(<R_{500})$ is the gas mass within $R_{500}$ and we have used Eqs. (2) and (4). The quantity $I(x_s, \alpha, \gamma, 3\beta)$ is the three dimensional integral value for $f_0 = 1$, which depends solely on the shape parameters.

We introduced a radially varying intrinsic scatter term around the model profile, assuming a log-normal distribution at

each radius. Taking into account measurement errors, the probability of measuring a given scaled gas density $\rho$ at given scaled radius $x$ for a cluster of mass $M_{500}$ at redshift $z$ is then

$$p(\rho|x, z, M_{500}) = N[\log \rho_m(x, z, M_{500}), \sigma^2(x)] \quad (10)$$

$$\sigma^2(x) = \sigma_{\text{int}}^2(x) + \sigma_{\text{stat}}^2, \quad (11)$$

where $N$ is the log-normal distribution. The variance term, $\sigma^2(x)$, is the quadratic sum of the statistical error, $\sigma_{\text{stat}}$ on the measured log $\rho$ and of the intrinsic scatter on log $\rho_m$ at radius $x$, $\sigma_{\text{int}}(x)$.

We expect the intrinsic scatter to increase towards the centre, as observed in Fig. 2, due to the increasing effect of non-gravitational physics on the density profiles. Ghirardini et al. (2019) studied the intrinsic scatter of massive local clusters in the XCOP sample, modelling the radially varying scatter $\sigma_{\text{int}}(x)$
<figure><img src="image_5.png" /><figcaption>Fig. 4. The universal cluster ICM density profile. Left: scaled density profiles of the SZE-selected clusters (grey points), overplotted with the best-fitting GNFW model with free evolution and mass dependence: $\rho_{\text{gas}}/\rho_{500}(R/R_{500}) \propto E(z)^{\alpha z} M^{\alpha M}$ with $\alpha_z = 2.09 \pm 0.02$ and $\alpha_M = 0.22 \pm 0.01$ (orange line). The model includes a radially varying intrinsic scatter term (orange envelope). Right: comparison of the best-fitting model, defined on the SZE-selected sample, to the best-fitting model for the X-ray-selected REXCESS sample (light blue). Here, the points with error bars are the REXCESS sample.</figcaption></figure>

with a log-parabola function. However, we found that such an analytical form significantly overestimates the scatter in the inner core, which was not covered by their data. To allow for more freedom we used a non-analytical form for the intrinsic scatter, where $\sigma_{\text{int}}(x)$ is defined at $n$ equally spaced points in $\log(x)$ in the typical observed radial range, $[x_{\text{min}}-x_{\text{max}}]$. The scatter, $\sigma_{\text{int}}(x)$ at other radii is computed by spline interpolation. We used $n = 7$ between $x_{\text{min}} = 0.01$ and $x_{\text{max}} = 1$.

The likelihood of a set of scaled density profiles measured for a sample of $i = 1, N_c$ clusters of mass $M_{500,i}$ and redshift $z_i$ is:

$$ \mathcal{L} = \prod_{i=0}^{N_c} \prod_{j=0}^{N_R[i]} p(\rho_{i,j} | x_{i,j}, z_i, M_{500,i}), \qquad (12) $$

where $N_R[i]$ is the number of points of the profile of cluster $i$, and the quantity $\rho_{i,j} = \rho_{\text{gas}[i,j]} / \rho_{500}(z_i)$ is the scaled density measured at each scaled radius $x_{i,j} = r[i,j] / R_{500}(z_i, M_{500,i})$, with $\rho_{\text{gas}[i,j]}$ and $r[i,j]$ being the physical gas density and radius. The statistical error on $\log \rho_{i,j}$ is $\sigma_{\text{stat},i,j}$.

We fitted the data (i.e. the set of observed $\rho_{i,j}$) using Bayesian maximum likelihood estimation with Markov chain Monte Carlo (MCMC) sampling. Using the emcee package developed by Foreman-Mackey et al. (2013), we maximised the log of the likelihood, which reads (up to an additive constant)

$$ \ln \mathcal{L} = -0.5 \sum_{i,j} \left[ \ln \sigma_{i,j}^2 + \frac{(\log \rho_{i,j} - \log \rho_{\text{m},i,j})^2}{\sigma_{i,j}^2} \right] \qquad (13) $$

$$ \rho_{\text{m},i,j} = \rho_{\text{m}}(x_{i,j}, z_i, M_{500,i}) \qquad (14) $$

$$ \sigma_{i,j}^2 = \sigma_{\text{int}}^2(x_{i,j}) + \sigma_{\text{stat},i,j}^2. \qquad (15) $$

The fit marginalises over a total of fourteen parameters: four describing the shape of the median profile $(x_s, \alpha, \gamma, 3\beta)$, a global

normalisation, $f_0$, the slopes $\alpha_M$ and $\alpha_z$ that describe the non-standard mass and evolution dependences, and seven additional parameters describing the intrinsic scatter profile. We used flat priors on all parameters.

## 3.2. Results

To establish a baseline, we fitted the model described above to the 93 SZE-selected systems. The resulting best-fitting model is

$$ \rho_{\text{m}}(x, z, M_{500}) = A(z, M_{500}) f(x), \qquad (16) $$

with

$$ A(z, M_{500}) = E(z)^{2.09 \pm 0.02} \times \left[ \frac{M_{500}}{5 \times 10^{14} h_{70}^{-1} M_{\odot}} \right]^{0.22 \pm 0.01} \qquad (17) $$

and

$$ f(x) = \frac{f_0}{(x/x_s)^\alpha [1 + (x/x_s)^\gamma]^{(3\beta - \alpha)/\gamma}}, \qquad (18) $$

where

$$ f_0 = 0.13 \pm 0.01, $$

$$ x_s = 0.28 \pm 0.01, $$

$$ \alpha = 0.42 \pm 0.06, $$

$$ \beta = 0.78 \pm 0.03, \text{ and} $$

$$ \gamma = 1.52 \pm 0.16. $$

Figure 3 shows the marginalised posterior likelihood for the parameters of the best-fitting density profile model detailed in Sect. 3.1. All parameters are well constrained: in particular, we note that the mass and evolution parameters $\alpha_M$ and $\alpha_z$ do not show any degeneracies, implying that we clearly separate the mass and redshift effects.
<table><caption>**Table 1.** Numerical values for the best-fitting intrinsic scatter term, measured at seven equally spaced points in $\log(R/R_{500})$, in the range $[0.01-1]R_{500}$.</caption><thead><tr><td>Radius</td><td>σ<sub>int</sub></td></tr></thead><tbody><tr><td>0.010</td><td>0.98 ± 0.06</td></tr><tr><td>0.021</td><td>0.75 ± 0.03</td></tr><tr><td>0.046</td><td>0.60 ± 0.02</td></tr><tr><td>0.100</td><td>0.38 ± 0.01</td></tr><tr><td>0.215</td><td>0.19 ± 0.01</td></tr><tr><td>0.464</td><td>0.12 ± 0.01</td></tr><tr><td>1.000</td><td>0.19 ± 0.01</td></tr></tbody></table>

The left-hand panel of Fig. 4 shows the density profiles of the SZE-selected clusters together with the best-fitting model. The intrinsic scatter term is represented by the orange envelope; the numerical values for this term are given in Table 1. The right-hand panel of Fig. 4 shows the best-fitting model for the SZE-selected systems compared to the profiles from the X-ray-selected REXCESS sample. The agreement is excellent beyond the core; in the inner regions, there is a hint that the X-ray-selected systems may show more dispersion. We will return to this point below in Sect. 5.1.5.

# 4. Luminosity scaling relations

We now turn to the scaling relation between $L_{Xc}$ and the mass $M_{500}$. The bolometric X-ray luminosity of a cluster can be written (Arnaud & Evrard 1999)

$$L(T) = f_{\text{gas}}^2(T)[M(T)\Lambda(T)]\hat{Q}(T), \quad (19)$$

where $f_{\text{gas}} = M_{\text{gas}}/M$ is the gas mass fraction, and $\Lambda(T)$ is the cooling function. The quantity $\hat{Q}(T) = \langle\rho_{\text{gas}}^2\rangle/\langle\rho_{\text{gas}}\rangle^2$ is a dimensionless structure factor that depends on the spatial distribution of the gas density (e.g., clumpiness at small scale, shape at large scale, etc.). Further assuming (i) virial equilibrium of the gas in the dark matter potential [$M \propto T^{3/2}$]; (ii) simple Bremsstrahlung emission [$\Lambda(T) \propto T^{1/2}$]; (iii) similar internal structure [$\hat{Q}(T) = \text{const.}$]; (iv) a constant gas mass fraction [$f_{\text{gas}}^2(T) = \text{const.}$], the standard self-similar relation between bolometric X-ray luminosity and mass, $L_X \propto M^{4/3}$, can be obtained. Similar arguments can be used to obtain the soft-X-ray luminosity-mass relation of $L_X \propto M$.

## 4.1. Fitting method

We fitted the data with a power-law relation of the form

$$E(z)^n L/L_0 = A_L (X/X_0)^{B_L}, \quad (20)$$

where $L_0 = 1 \times 10^{44} \text{ erg s}^{-1}$ and $5 \times 10^{44} \text{ erg s}^{-1}$ for the soft and bolometric bands, respectively, and $X_0 = 5 \times 10^{14} M_{\odot} \text{ keV}$ and $4 \times 10^{14} M_{\odot} \text{ keV}$ for $Y_X$ and $M$, respectively. Fitting was undertaken using linear regression in the log-log plane, taking uncertainties in both variables into account, and including the intrinsic scatter. We fitted the data using a Bayesian maximum likelihood estimation approach with Markov chain Monte Carlo (MCMC) sampling. We write the likelihood as defined by

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Relation between $L_{Xc}$ and $Y_X$. The blue envelope is the best-fitting relation given in Eq. (23), and the results from Maughan (2007) are also shown for comparison.</figcaption></figure>

Robotham & Obreschkow (2015)

$$\ln \mathcal{L} = \sum_{i=1}^{N} \left[ \ln \frac{B_L^2 + 1}{\sigma_i^2} - \frac{(\ln(L_i/L_0) - \ln A_L - B_L \ln(X_i/X_0))^2}{\sigma_i^2} \right], \quad (21)$$

with

$$\sigma_i^2 = \sigma^2 + B_L^2 \sigma_{X,i}^2 + \sigma_{L,i}^2 \quad (22)$$

and the intrinsic scatter, $\sigma^2$, as a free parameter. MCMC sampling was undertaken using the `emcee` package developed by Foreman-Mackey et al. (2013), with flat priors in the ranges $[-2.0, 3.0]$ and $[0.0, 1.0]$ for $B_L$ and $\sigma$, respectively. The results, reported in Table 2, were compared to those obtained with the `LINMIX` (Kelly 2007) Bayesian regression package: these were indistinguishable and so are not reported here.

## 4.2. Results

### 4.2.1. $L_{Xc} - Y_X$

We first fitted the relation between the bolometric $L_{Xc}$ and the mass proxy $Y_X$ for the SZE-selected systems only. With the evolution term left free, the best-fitting relation, shown in Fig. 5, is

$$E(z)^{-1.79 \pm 0.08} (L_{Xc}/L_0) = (1.49 \pm 0.03) (Y_X/Y_0)^{0.97 \pm 0.02}, \quad (23)$$

with an intrinsic scatter of $\sigma_{L_{Xc}|Y_X} = 0.09 \pm 0.01$. This result yields evolution and mass dependences that are in excellent agreement with the self-similar predictions of -9/5 and 1.0, respectively. It is also in excellent agreement with the REXCESS only results of Pratt et al. (2009) and that of Maughan (2007), the latter of which was estimated from *Chandra* data and is overplotted on the figure.
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Relation between the core-excised X-ray luminosity $L_{Xc}$ and mass estimated from the $Y_X$ proxy, for the bolometric and soft-band luminosities of 118 systems. **Left:** data points with the best fitting relation to the X-ray-selected REXCESS sample with the evolution factor fixed to the self-similar value of $n = -2$ (grey envelope). The dark grey envelope shows the best fitting relation to the 37 SZE-selected systems at $z < 0.3$ with $n = -2$. **Right:** histogram of the log space residuals from the best fitting relation to the SZE-selected objects at $z < 0.3$.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Relation between the core-excised X-ray luminosity $L_{Xc}$ and mass estimated from the $Y_X$ proxy, for the bolometric and soft-band luminosities of 118 systems. **Left:** best fitting relation (grey envelope) to the full sample (data points) with the evolution factor $n$ left free to vary. The best-fitting values of $n$ are given in Table 2. **Right:** histogram of the log space residuals from the best fitting relation. Solid lines show the best-fitting Gaussian distributions with $\sigma$ corresponding to the best-fitting intrinsic scatter in log space (Table 2).</figcaption></figure>

### 4.2.2. $L_{Xc}$ – $M$: Low redshift with fixed evolution

We initially fixed the evolution factor, $n$, to the self-similar values of $-2$ and $-7/3$ for the soft and bolometric bands, respectively. We first fitted the bolometric $L_{Xc}-M$ for the REXCESS data only, using a mass pivot of $2 \times 10^{14} M_\odot$, as used by Pratt et al. (2009). The resulting normalisation, $A_L = 1.06 \pm 0.04$, slope $B_L = 1.73 \pm 0.05$, and intrinsic scatter $\sigma_{\ln L_{Xc}} = 0.16 \pm 0.03$, are in excellent agreement with those found by Pratt et al. (2009)

using orthogonal Bivariate Correlated Errors and intrinsic Scatter (BCES; Akritas & Bershady 1996) fitting. Similarly good agreement was found for the soft-band $L_{Xc}-M$ relation, showing that the scaling relation parameters that are obtained for the X-ray-selected are robust to the change in underlying $M-Y_X$ relation used to estimate the mass, and also to differences in fitting method.

We then fitted the SZE-selected clusters at $z < 0.3$ (37 systems). The results are given in Table 2 and show that the
<table><caption>**Table 2.** Fits to the core-excised X-ray luminosity $L_{Xc} - M$ relation.</caption><thead><tr><th>Relation</th><th>Selection</th><th>Redshift</th><th>n</th><th>L<sub>0</sub><br>(10<sup>44</sup> erg s<sup>-1</sup>)</th><th>A<sub>L</sub></th><th>B<sub>L</sub></th><th>σ<sub>ln L</sub></th></tr></thead><tbody><tr><td rowspan="3">bolometric</td><td>REXCESS</td><td>0.05 &lt; z &lt; 0.20</td><td>-7/3</td><td>5</td><td>0.70 ± 0.03</td><td>1.73 ± 0.05</td><td>0.16 ± 0.03</td></tr><tr><td>SZ</td><td>z &lt; 0.3</td><td>-7/3</td><td>5</td><td>0.73 ± 0.02</td><td>1.69 ± 0.05</td><td>0.09 ± 0.02</td></tr><tr><td>All SZ</td><td>z &gt; 0.05</td><td>-2.50 ± 0.09</td><td>5</td><td>0.70 ± 0.02</td><td>1.74 ± 0.02</td><td>0.12 ± 0.01</td></tr><tr><td rowspan="3">0.5–2 keV</td><td>REXCESS</td><td>0.05 &lt; z &lt; 0.20</td><td>-2</td><td>1</td><td>1.06 ± 0.04</td><td>1.48 ± 0.05</td><td>0.17 ± 0.03</td></tr><tr><td>SZ</td><td>z &lt; 0.30</td><td>-2</td><td>1</td><td>1.12 ± 0.04</td><td>1.37 ± 0.06</td><td>0.13 ± 0.02</td></tr><tr><td>All SZ</td><td>z &gt; 0.05</td><td>-2.23 ± 0.10</td><td>1</td><td>1.05 ± 0.03</td><td>1.46 ± 0.03</td><td>0.15 ± 0.02</td></tr></tbody></table>

**Notes.** $L_{Xc}$ is measured in the [0.15–1] $R_{500}$ region. The data were fitted in log-log space with a relation of the form $E(z)^n (L/L_0) = A_L (M/4 \times 10^{14} M_\odot)^{B_L}$.

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

**Fig. 8.** Scaled density profiles. Left: scaled density profiles (points) and best-fitting model (orange envelope) for SZE-selected systems in our sample compared to the median and 68% dispersion from the X-COP sample (Ghirardini et al. 2019, magenta envelope). Right: comparison of best fitting intrinsic scatter model (blue) with that found by Ghirardini et al. (2019, magenta). The best-fitting intrinsic scatter obtained from our sample when the evolution factor is forced to the self-similar value of $E(z)^2$ is also shown in grey.

normalisation and slope for this sub-sample are in agreement within $1\sigma$ with those found for REXCESS. This indicates that there is no difference in the scaling relation between the local X-ray and SZE-selected samples, once the core region has been excised. There is a slight hint that the intrinsic scatter of the SZE-selected sample about the best-fitting relation is lower than that for REXCESS, although this is only a ~2$\sigma$ effect for the bolometric luminosity, and is less significant for the soft-band. The data and best-fitting relations, including the $1\sigma$ scatter envelopes, are shown in the left-hand panel of Fig. 6.

### 4.2.3. $L_{Xc} - M$: Free evolution

The right-hand panel of Fig. 6 shows the histogram of the residuals of the full SZE-selected sample (93 systems, $0.08 < z < 1.13$) with respect to the best-fitting relation to the systems at $z < 0.3$. The peak is offset by $\Delta \ln L \le 0.1$, indicating that some evolution beyond self-similar is in fact needed.

We then fitted the full SZE-selected sample with a power-law relation, including a free evolution factor, $n$. The results are given in Table 2 and the best-fitting relations are shown in the left-hand panel of Fig. 7; the right-hand panel of Fig. 7 shows

the residual histograms. The latter are well-centred on zero. The best-fitting evolution terms, $n = -2.23 \pm 0.09$ for the soft-band and $n = -2.50 \pm 0.10$ for the bolometric luminosity, suggest that stronger than self-similar evolution is significant at the ~2$\sigma$ level.

## 5. Discussion

### 5.1. Gas density profiles

#### 5.1.1. Comparison with previous work

Pioneering work on parametric models of scaled density profiles (Neumann & Arnaud 1999) obtained from ROSAT allowed the dispersion in radial slopes to be constrained. Croston et al. (2008) studied the scaled density profiles of the REXCESS sample, obtaining for the first time constraints on the radial dependence of the intrinsic scatter. The scaled density profile of the X-COP sample, obtained assuming a self-similar evolution factor $E(z)^2$, was presented in Ghirardini et al. (2019). Figure 8 compares the scaled density profiles and the best-fitting model from our SZE-selected sample to their median scaled density
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Deviations from the average scaling with radius. **Left:** redshift-mass distribution of the SZE-selected sample used in this work. The shaded regions indicate cuts for two sub-samples: a large redshift range at nearly fixed mass, and a large mass range at nearly fixed redshift. **Middle:** degree to which the radial ICM density profile evolves as a function of redshift at nearly fixed mass. The dotted line shows the self-similar expectation ($\alpha_z = 2$). The dashed line shows the best-fitting evolution, which varies from slower than self-similar in the centre ($\alpha_z \sim 0.3$) to faster than self-similar around $R_{500}$ ($\alpha_z \sim 2.4$). Envelopes show the 1 and $2\sigma$ uncertainties. **Right:** degree to which the radial ICM density profile scales with mass at nearly fixed redshift. The dotted line shows the self-similar expectation ($\alpha_z = 0$). The dashed line shows the best-fitting mass dependence of $\alpha_z = 0.22$. The scaled density at nearly fixed redshift does not depend on radius.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Scaled profiles and scatter. **Left:** median scaled profiles (solid lines) and 68% scatter (envelopes) for REXCESS and the SZE-selected sample split into three redshift bins. Beyond ~0.2 $R_{500}$ the scaled profiles are almost indistinguishable. **Right:** radial profile of the intrinsic scatter for the various sub-samples. The best-fitting intrinsic scatter model obtained from the SZE-selected sample is also shown. Intrinsic scatter is less than 20% between 0.2 ≲ $R_{500}$ ≲ 1.0. The gold line shows the model intrinsic scatter profile corrected for the covariance between $M_{\text{gas}}$ and $R_{500}$ (see Sect. 5.1.4), which results in a suppression of the scatter by a factor of about two at $R_{500}$.</figcaption></figure>

profile and 68% dispersion. The agreement is good out to the maximum X-COP radius of ~2 $R_{500}$, although with subtle differences in the inner regions ($R < 0.1 R_{500}$). Their profile is less peaked, likely due to their not having corrected the density profiles for PSF effects, and has a smaller dispersion than our sample, which may be linked to their more limited mass coverage. The right-hand panel compares the intrinsic scatter measurements, which also agree quite well, although the scatter of the present sample is better constrained. The best-fitting intrinsic scatter obtained from our sample when the evolution factor is forced to the self-similar value of $E(z)^2$ is also shown in grey.

We can also compare to the results obtained by Mantz et al. (2016), who modelled the evolution and mass dependence of the scaled density profiles of a morphologically relaxed cluster sample of 40 systems at $0.08 < z < 1.06$. While their evolution dependence of $2.0 \pm 0.2$ is in agreement with our results, they find a mass dependence that is consistent with zero ($0.03 \pm 0.06$). The difference with respect to our results may be due simply to cluster selection. They studied dynamically relaxed, hot systems, for which the mass leverage is more limited. Once scaled, they found a scatter in scaled density of ≲20% at $R_{2500}$. For the typical mass of the present sample, $R_{2500} \approx 0.45 R_{500}$, where our intrinsic scatter measurements are in good agreement with theirs.
### 5.1.2. Radial dependence of the scaling

Once the best fitting model was obtained, we quantified how well the model represents the data by calculating the variation of scaled density at different scaled radii. To better disentangle redshift and mass evolution, we extracted two sub-samples: one covering a large redshift range at nearly constant mass, and another covering a large mass range at nearly constant redshift. These sub-samples are illustrated in the $z-M$ plane by the green and blue regions in the left-hand panel of Fig. 9. We defined ten radial bins in terms of $x = r/R_{500}$ and measured the gas density in each bin. We then scaled the density by the best-fitting model and fitted a function of the form $\rho_{gas}(x)/\rho_{crit,0} \propto E(z)^{\alpha_z} M^{\alpha_M}$. For a self-similarly evolving population, the gas density scales with the critical density $\rho_{crit}$ and there is no mass dependence, and so $\alpha_z = 2$ and $\alpha_M = 0$.

The middle- and right-hand panels of Fig. 9 show the degree to which the two sub-samples vary from self-similar scaling and from the best-fitting model scaling, as a function of scaled radius. Uncertainties are large because of the reduced number of data points and the scatter in the data. At nearly constant mass the overall variation with redshift is slightly greater than self-similar. However, the density evolves differently in the core and in the outer regions. The density evolution is consistent with zero in the core: $\alpha_z = 0.28 \pm 1.10$ at $x = 0.01$, a value that is in good agreement the result found by McDonald et al. (2017) although with large uncertainties. However, the density in the outer regions appears to evolve more strongly than self-similar: $\alpha_z = 2.42 \pm 0.22$ at $x = 1.20$, a result that is significant at slightly more than $1\sigma$. At nearly fixed redshift, the density varies with mass in agreement with the $\alpha_M = 0.22 \pm 0.01$ scaling established above. At $R_{500}$, this mass dependence is significant at $\sim 2\sigma$, but does not depend on radius.

### 5.1.3. Median and scatter of scaled profiles

We now turn to the ensemble properties of the scaled density profiles. The left-hand panel of Fig. 10 shows the median scaled profile, obtained in the log-log plane, and 68% dispersion for REXCESS X-ray-selected sample and the SZE-selected sample split into three redshift bins. It is clear that once scaled, the four sub-samples are remarkably similar beyond $0.2 R_{500}$. In the core region, the median central density decreases progressively with redshift. The median scaled central densities of the REXCESS and of the SZE-selected sample at $z < 0.3$ are virtually indistinguishable. We will return to the central regions in Sect. 5.1.5 below.

The radial variation of the intrinsic scatter about these median profiles is quantified in the right-hand panel of Fig. 10, together with that of the best-fitting intrinsic scatter model obtained above (from the SZE-selected clusters only). The intrinsic scatter of this model falls below 20% in the radial range $0.2 \le R_{500} \le 1.0$. There is excellent agreement between the best-fitting intrinsic scatter model and the observed profiles, which all follow broadly the same trend with scaled radius: a steep decrease with a minimum at $\sim 0.5-0.7 R_{500}$, followed by an increase towards larger radii.

Within $R \le 0.2 R_{500}$ the intrinsic scatter in the scaled gas density profiles climbs steeply towards the centre. This increase is intimately linked to the complex physics of the core regions, dynamical activity, and to the presence or absence of cool core systems in the various samples. In this connection, the sample with the largest intrinsic scatter in the central regions is REXCESS, reflecting the presence of cool core systems in this dataset.

Beyond $\sim 0.2 R_{500}$ the relative dispersion of all samples dips below 20%, and at $R_{500}$ the dispersion in profiles is $\sim 15\%$. In the SZE-selected sub-samples, there is a clear evolution, in the sense that the low-redshift systems exhibit the lowest intrinsic scatter values while the high-redshift systems show higher values. The scatter will be related to intrinsic cluster-to-cluster variations linked to inhomogeneities that will depend on the mass accretion rate and associated dynamical state, together with a component due to uncertainties in the total mass. It is possible that both of these effects conspire to produce higher intrinsic scatter values for higher-redshift systems: one expects an increase in dynamical activity with redshift, while uncertainties in the cluster mass measurement will also increase in the same sense.

### 5.1.4. Suppression of scatter due to covariance between $M_{gas}$ and $R_{500}$

The observed scatter in the density profiles may be suppressed by the use of $M_{gas}$ in the computation of $R_{500}$ when scaling the radial coordinate. All other things being equal, a cluster with a higher than average $\rho_{gas}$ (for its mass) at some radius, will have a higher than average $M_{gas}$ and hence $Y_X$ relative to its mass. Since $Y_X$ is then used to estimate $R_{500}$ assuming a mean scaling relation, one would then overestimate $R_{500}$ for this cluster. The radial scaling for this cluster would then be too large, which would move its density profile back towards the mean profile, reducing the apparent scatter. The reduction in scatter will depend on the slope of the density profile (i.e. the reduction is larger where the profile is steeper), so will be radially dependent.

The possible magnitude of this effect was estimated by generating synthetic cluster density profiles with a known amount of scatter, and then scaling them in radius following the method used for the observed clusters in order to test how much the scatter was changed. In more detail, the best-fitting median profile presented in Sect. 3.2 was normalised to match a 6 keV cluster at $z = 0.15$ (i.e. when integrated to $R_{500}$, the gas mass and $Y_X$ were consistent with the scaling relations used for the observed clusters). For these reference values, the 'true' $R_{500}$ is 1210 kpc.

A large number of realisations of this median profile were then generated by resampling the normalisation from a lognormal distribution with a standard deviation of 0.2 (approximating a constant 20% scatter in $\rho_{gas}$ at all radii). For each realisation, $R_{500}$ was then computed in the same manner as for the observed clusters; the profile was integrated to compute $Y_X$ (assuming a fixed temperature of 6 keV) and hence $R_{500}$, with the process performed iteratively until $R_{500}$ converged. The profile was then scaled in radius by $R_{500}$ and the process was repeated for each realisation of the density profile.

When the distribution of densities in the realisations was measured at the 'true' value of $R_{500} = 1210$ kpc, the input scatter of 20% was recovered. However, when the profiles were each scaled in radius by the value of $R_{500}$ estimated for each realisation from the $M-Y_X$ relation, the scatter at a scaled radius of unity was found to be 10%; that is, the scatter is suppressed by a factor of two at around $R_{500}$ due to the dependence of $R_{500}$ on $M_{gas}$ in our analysis. The same factor of two suppression was found for different values of the input scatter. We calculated the reduction in scatter at different scaled radii, obtaining a radial profile of the suppression factor. The intrinsic scatter profile corrected for this suppression factor is plotted in gold in Fig. 10. This method makes a number of simplifying assumptions (e.g., there is no scatter in $T$, the profile is fixed at the median form), but based on this analysis, we estimate that the scatter measured
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Central density. Left: histogram of central densities for the SZE-selected systems, scaled according to the best-fitting model (Eq. (17)) derived in Sect. 3.2, measured at $0.05 R_{500}$. The solid line is a kernel density plot with a smoothing width of 0.15. Middle: histogram of scaled central densities for the $z < 0.3$ SZE-selected systems compared to the X-ray-selected sample, measured at $0.015 R_{500}$. The solid line is a kernel density plot with a smoothing width of 0.15. Right: histogram of central densities for the SZE-selected sample at 40 kpc.</figcaption></figure>

for the observed clusters is likely to be underestimated by a factor of approximately two at $R_{500}$.

### 5.1.5. Change of central regions over time

McDonald et al. (2017) showed that while the ICM outside the core regions of their SZE-selected sample, covering the redshift range $0.25 < z < 1.2$, evolved self-similarly with redshift, the central absolute median density (i.e. expressed in units of $cm^{-3}$) did not. They interpreted this result as being due to an un-evolving core component embedded in a self-similarly evolving bulk.

Our sample is of comparable size to that of McDonald et al. (2017), but the evolution with mass and redshift has been decoupled and quantified (Sect. 3.1). The effective 5'' resolution of XMM-Newton after PSF correction (Bartalucci et al. 2017) allows us to measure the density of the SZE-selected sample across all redshifts down to a scaled radius of $R \sim 0.05 R_{500}$. The left-hand panel of Fig. 11 shows a histogram of the resulting values<fnref n="4" /> scaled according to the best-fitting model (Eq. (17)) established in Sect. 3.2. The SZE-selected sample was further divided into three redshift bins to better visualise how the sample changes over time. This histogram is characterised by a strong peak centred on a scaled central density of –0.2 (in log space), which is clearly visible, and coincident, in all three SZE-selected sub-samples. While the histogram of the $z > 0.6$ sub-sample has no detectable skewness, the histogram of the $z < 0.3$ sub-sample exhibits a distinct tail to higher scaled central densities, which is characterised by a moderately large positive skewness of $G_1 = 0.86$ that is significant at $> 90\%$ (Doane & Seward 2011). This may indicate the gradual appearance of objects with more peaked scaled central densities towards lower redshifts.

At $z < 0.3$, the effective 5'' resolution of XMM-Newton after PSF correction allows us to measure the density down to a scaled radius of $R \sim 0.015 R_{500}$. The middle panel of Fig. 11 shows a histogram of the scaled central density of the SZE-selected clusters at $z < 0.3$ compared to that of REXCESS. The positive skewness of the SZE-selected systems is confirmed at greater significance ($G_1 = 1.06$), while the histogram of the REXCESS sample exhibits two peaks in scaled central density: a main peak that is coincident with the peak of the SZE-selected sub-samples,

and a secondary peak at a scaled density of 0.7 (in log space). The latter peak is due to cool core systems and may indicate that centrally peaked systems are over-represented in X-ray-selected samples, as has been argued by Rossetti et al. (2017) from their comparison of the image concentration parameter in Planck clusters to those for X-ray-selected systems, and also by Andrade-Santos et al. (2017).

The right-hand panel of Fig. 11 shows the histogram of the central density of the SZE-selected sample at 40 kpc, measured in physical units ($cm^{-3}$). There is a broad maximum at $n_{e,40 kpc} \sim 0.01 cm^{-3}$, and the histograms of the three sub-samples coincide. A Kolmogorov-Smirnov test indicates that all three sub-samples come from the same parent distribution. This result suggests, in agreement with McDonald et al. (2017), that the absolute central density remains constant over the redshift range probed by the current sample.

### 5.2. The $L_{Xc} - M$ relations

#### 5.2.1. Comparison with other work

Figure 12 shows the best-fitting bolometric $L_{Xc} - M$ relation for the SZE-selected clusters in the present sample compared to a number of results from the literature (Maughan 2007; Mantz et al. 2010; Bulbul et al. 2019; Lovisari et al. 2020). With the exception of those obtained by Mantz et al. (2010), these studies generally find slopes that are steeper than the self-similar expectation of 4/3, ranging from 1.63 to 1.92. Studies that put constraints on the evolution with redshift (Mantz et al. 2010; Bulbul et al. 2019; Lovisari et al. 2020) generally find good agreement with self-similar expectations (although with large uncertainties). However, any measurement of the dependence of a quantity on the mass will be affected strongly by the sample selection and on how the mass itself has been measured. Data fidelity and sample sizes are now such that systematic effects are starting to become dominant over measurement uncertainties.

Concerning the sample selection, the results in Sect. 4.2 show that the $L_{Xc} - M$ relations of X-ray- and SZE-selected systems are in good agreement, suggesting that once the core regions are excluded, effects due to detection methods relying on the ICM do not have any impact. Similarly, Lovisari et al. (2020) showed that their relaxed and disturbed samples had similar $L_{Xc} - M$ relation slopes and normalisations. This suggests

<footnote n="4">⁴ The measurement requires a small extrapolation, by <1'' in the log-log plane, for seven systems.</footnote>
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Comparison of the bolometric $L_{Xc}-M$ relation to previous work (Maughan 2007; Mantz et al. 2010; Bulbul et al. 2019; Lovisari et al. 2020).</figcaption></figure>

that the $L_{Xc} - M$ relation may also be relatively robust to selection effects linked to cluster dynamical state, likely due to the small intrinsic scatter.

A more fundamental issue is the mass measurement itself (Pratt et al. 2019). In the present work we have used $Y_X$ as a mass proxy; the works listed above use variously $Y_X$ (Maughan 2007), the gas mass (Mantz et al. 2010), the SZE signal-to-noise (Bulbul et al. 2019), and the hydrostatic mass (Lovisari et al. 2020), as proxies. In this context, the shallower slope of the Mantz et al. (2010) relation compared to the others can be fully explained by their assumption of a constant gas mass fraction in the mass calculation (e.g., Rozo et al. 2014b).

All of the above mass estimates are derived from ICM observables, and all except Lovisari et al. (2020) use scaling laws that have been calibrated on X-ray hydrostatic mass estimates. Independent mass measurements, such as those available from lensing, galaxy velocity dispersions, or caustic measurements (e.g., Maughan et al. 2016), are critical to making progress on this issue. In this connection, weak-lensing mass measurements for individual clusters have been carried out by several projects, such as the Local Cluster Substructure Survey (LoCuSS; Okabe et al. 2010, 2013; Okabe & Smith 2016), the Canadian Cluster Comparison Project (CCCP; Hoekstra et al. 2012, 2015), the Cluster Lensing And Supernova survey with Hubble (CLASH; Merten et al. 2015; Umetsu et al. 2014, 2016), Weighing the Giants (WtG; von der Linden et al. 2014; Kelly et al. 2014; Applegate et al. 2014), and CHEX-MATE (CHEX-MATE Collaboration 2021). The cluster community is undertaking a major ongoing effort to critically compare various mass estimates, obtained from mass proxies and from direct X-ray, lensing, or velocity dispersion analyses (e.g., Rozo et al. 2014a; Sereno & Ettori 2015, 2017; Sereno 2015; Groener et al. 2016). Ultimately, this effort will help to better constrain the parameters of the scaling relations.

### 5.2.2. Link between density profile and $L_{Xc}$

As noted in for example Maughan et al. (2008), the similarity of the ICM density profiles outside the core implies a low scatter in $L_{Xc}$, as is indeed observed here. In order to explore how much of the scatter in $L_{Xc}$ is due to the variation in density profiles, we computed a 'pseudo luminosity' $L'_{Xc} = \int \rho_{\text{gas}}^2 T_X^{1/2} dV$ for each cluster. For this calculation, we used the measured density profile for each cluster and assumed isothermality at the measured core-excised temperature for the cluster. The integral was performed over a cylindrical volume from projected radii of 0.15 $R_{500}$ to $R_{500}$. For 16 of 118 clusters, the density profiles did not reach $R_{500}$, so the integrals were truncated at the maximum observed radius. In all cases the profiles reached to ≈ 90% of $R_{500}$, and the contribution to the luminosity of the outer parts of the profile is very small, so the effect of this truncation is negligible. The scatter in $L'_{Xc}$ then provides an estimate of the scatter in the bolometric $L_{Xc}$ due only to the scatter in density profiles.

The intrinsic scatter about the best fitting relation to $L'_{Xc} E(z)^{-2.5}$ versus $M_{500}$ was then measured (assuming that the fractional statistical error on $L'_{Xc}$ is the same as that in $L_{Xc}$ for each cluster), giving a value of 11%. This implies that most or all of the intrinsic scatter in the bolometric $L_{Xc}$ at fixed mass can be explained by the variation in the ICM density profiles.

In principle, the use of $M_{YX}$ to determine $R_{500}$ and hence the aperture within which $L_{Xc}$ and $L'_{Xc}$ are measured, could introduce additional scatter in the $L_{Xc} - M_{500}$ relation. If $M_{YX}$ were scattered high relative to the true mass of a cluster, then $R_{500}$ would be overestimated and the 0.15–1 $R_{500}$ aperture would be shifted to larger radii. This shift would reduce $L_{Xc}$ since more of the luminosity comes from the inner edge of the aperture than the outer edge. Hence, if $M_{YX}$ were scattered high, then $L_{Xc}$ and $L'_{Xc}$ would be scattered low, adding to the observed scatter in the $L'_{Xc} - M_{500}$ relation. We examined the impact of this effect by adding 10% scatter to $M_{YX}$ when computing $L'_{Xc}$. This increased the measured scatter in $L'_{Xc}$ by less than 1%. The dependence of the luminosity aperture on $M_{YX}$ leads to a negligible contribution to the scatter in the $L_{Xc} - M_{YX}$ relation.

We therefore conclude that the measured intrinsic scatter in the bolometric $L_{Xc}$ at fixed mass is dominated almost entirely by the variation in the ICM density profiles. The results do not depend on the aperture for reasonable assumptions on the scatter between $M_{YX}$ and the true mass $M_{500}$. Any residual scatter will come from inhomogeneity and/or substructure in the density distribution, or from the effects of structure in the temperature or metallicity distribution.

### 5.2.3. Impact of selection bias and covariance

We found above that the bolometric $L_{Xc} - M$ relation has a slope $B_L \sim 1.7$, which is steeper than the self-similar value of 1.3, and has a small scatter of $\sigma_{\ln L_{Xc}} \sim 0.15$. A similar difference in slope of $\Delta B_L \sim +0.4$ is seen relative to the self-similar soft-band $L_{Xc} - M$ relation. One might wonder whether and to what extent selection effects and covariance may impact these results. Firstly, in survey data the observable that is used for cluster detection, $O_{\text{det}}$, can be affected by so-called Malmquist bias. This happens because clusters that are scattered to higher values of $O_{\text{det}}$, whether by noise or by the intrinsic scatter between the observable and the mass, will be preferentially detected, leading to a positive bias in the average value of $O_{\text{det}}$. This is a particular concern near the detection threshold, and will affect the apparent slope of any scaling relation with $O_{\text{det}}$. Secondly, the intrinsic scatter of each quantity around the mean relation may well be
correlated, leading to a reduction in the measured scatter with respect to the true underlying value.

Since we showed in Sect. 4.2 that the results for the X-ray-selected sample are compatible within $1\sigma$ with those for the local SZE-selected sample, we focus here on the SZE selection. As detailed in Sect. 2.1, the SZE-selected sample is composed of four different sub-samples and with different selections: the ESZ, LP1, LP2 and LP3. Although the ESZ sub-sample was generated from high signal-to-noise ratio ($S/N > 6$) detections, Andrade-Santos et al. (2021) showed that its selection is dominated by the instrumental and astrophysical scatter, and that the intrinsic scatter in the $Y_{SZ}-M$ relation has a negligible effect on the recovered scaling parameters. The three other SZE sub-samples (LP1, LP2, LP3) are selected at an $S/N$ lower than that of the ESZ sample, and are thus even more affected by the instrumental and astrophysical scatter. We restricted in consequence our selection study to the ESZ. This will provide upper limits to the effects of selection on the results for the $L_{Xc}-M$ of the SZE-selected systems in the current study.

Appendix C describes in detail the simulations we used, which were similar to those produced for the Andrade-Santos et al. (2021) study. Simulated clusters, modelled with the Arnaud et al. (2010) pressure profile and drawn from a Tinker mass function (Tinker et al. 2008), were injected into the Planck Early SZ maps in the ‘cosmological’ mask region. The Multi-Matched Filter extraction algorithm (Melin et al. 2006) was then applied, to obtain SZ detections at $S/N > 6$, corresponding to the threshold for the ESZ sample. We then matched the injected and recovered clusters to produce a mock ESZ catalogue, doing this twenty times to generate 3188 detections in total. To investigate selection and covariance effects, we assumed a Gaussian lognormal correlated distribution for $Y_{SZ}$, $Y_X$, and $L_{Xc}$ at fixed true mass $M$, with a covariance matrix to describe correlations between parameters.

A first-order estimate of the expected scatter in the $L_{Xc}-M_{YX}$ relation can be obtained by assuming that selection effects are negligible. In this case,

$$ \sigma^2_{\ln L_{Xc}|M_{YX}} = \sigma^2_{\ln L_{Xc}|M} + (B_L^2/B_Y^2) \sigma^2_{\ln Y_{X}|M} - 2t \sigma_{\ln L_{Xc}|M} \sigma_{\ln Y_{X}|M}. \quad (24) $$

For a measured intrinsic scatter $\sigma_{\ln L_{Xc}|M_{YX}} \sim 0.13$, a covariance between $L_{Xc}$ and $Y_X$ of $t = 0.85$ (Farahi et al. 2019), and assuming a scatter between $Y_X$ and true mass $M$ of $\sigma_{\ln Y_{X}|M} \sim 0.16$ motivated by recent numerical simulations (Planelles et al. 2014; Le Brun et al. 2017; Truong et al. 2018), the resulting dispersion between core-excised luminosity and true mass is $\sigma_{\ln L_{Xc}|M} \sim 0.22$. This suggests that the measured scatter is underestimated by a factor of $0.22/0.13 \sim 1.7$ due to the use of $Y_X$ as a mass proxy.

We generated a series of simulations using the above values for true scatter and covariance between quantities, and used them to investigate the effects of selection, intrinsic scatter, and covariance on the results. We found:

- The Malmquist bias induced on the $L_{Xc}-Y_{SZ}$ relation due to intrinsic scatter in the $D_A^2 Y_{SZ}-M$ relation is completely negligible, with the slope changing by less than 0.5%. The resulting impact on the slope of the $L_{Xc}-M$ relation is $\Delta B_L = -0.01$. This is important, because it implies that selection effects due to intrinsic scatter in the $Y_{SZ}$ observable cannot account for the observed steeper slope of the $L_{Xc}-M$ relations seen here.

- Use of $Y_X$ as a mass proxy introduces additional intrinsic scatter with respect to the underlying mass. The net impact

on the recovered slope of the $L_{Xc}-M$ relation is minimal, with $\Delta B_L = -0.04$.

- Intrinsic scatter in the $L_{Xc}-Y_X$ relation has the effect of redressing the slope towards its original value, with $\Delta B_L = +0.02$.

- Covariance between $Y_X$ and $L_{Xc}$ again changes the slope. For $t = 0.85$ (Farahi et al. 2019), the impact on the slope of the $L_{Xc}-M$ relation is $\Delta B_L = +0.04$.

We thus conclude that the slope of the $L_{Xc}-M_{YX}$ relation found here is robust to selection effects due to intrinsic scatter in the $Y_{SZ}$ and $Y_X$ proxies, and to covariance between quantities. The dispersion, however, is very sensitive to the covariance. We estimate that the measured dispersion in the relation between the core-excised luminosity and true mass, $L_{Xc}-M$, is underestimated by a factor of $\sim 1.7$ due to the use of $Y_X$ as a mass proxy.

### 5.3. Link to $f_{gas}-M$

The $L_{Xc}-M$ relations derived above have a steeper dependence than self-similar, which cannot be explained by selection effects, intrinsic scatter, or covariance. There is evidence for a dependence of the gas content with mass, which has the effect of suppressing the luminosity preferentially in lower-mass systems, leading to the observed steepening of the relation (see discussion in Pratt et al. 2009). Interestingly, with the assumption of a standard dependence of temperature on mass ($T \propto M^{2/3}$, e.g., Arnaud et al. 2005; Mantz et al. 2016), use of Eq. (19) with the observed bolometric $L_{Xc} \propto M^{1.7}$ relation yields $f_{gas} \propto M^{0.21}$.

The gas density profile model detailed in Sect. 3.1 yields an alternative method to obtain the dependence of the gas mass fraction with mass, as the model includes a dependence of $f_{gas} = M_{gas}/M$ on the total mass $M$, via Eq. (8). The best-fitting gas density model (Eq. (17)) suggests $f_{gas} \propto M^{0.22 \pm 0.01}$. This dependence of the gas content on total mass is therefore in good agreement with the expected relation given the observed $L_{Xc}-M$ dependence, and with previous findings (e.g., Pratt et al. 2009; Lovisari et al. 2015; Ettori 2015) based primarily on X-ray-selected clusters.

## 6. Conclusions

We have examined the gas density profiles and the relation between the core excised X-ray luminosity $L_{Xc}$ and the total mass derived from the $Y_X$ mass proxy for 118 X-ray and SZE-selected objects covering a mass range of $M_{500} = [0.5-20] \times 10^{14} M_{\odot}$ and extending in redshift up to $z \sim 1.13$. We first examined the scaled density profiles:

- The gas density profiles do not scale perfectly self-similar, exhibiting subtle trends in mass and redshift.

- Motivated by this finding, we fitted an analytic gas density model to the 93 SZE-selected systems. The analytic model is based on a generalised NFW profile, and correctly reproduces the scaled gas density profile and the radial variation of its intrinsic dispersion. Combined with the empirical mass scaling of the profiles, this analytic model defines the gas density profile of SZE-selected clusters as a function of mass and redshift. This model is given in Eqs. (16)–(18).

- The intrinsic dispersion in scaled profiles is greatest in the central regions, declining to a minimum at $\sim 0.5-0.7 R_{500}$, and increasing thereafter. The dispersion is similar for X-ray-selected clusters and for local SZE-selected clusters, except in the centre, where the X-ray-selected systems have a higher
dispersion. There is a hint for an evolution of the disper-
sion with redshift, which may be linked to an increase in
perturbed clusters at higher redshifts.

- We investigated the effect of covariance between $M_{gas}$ and $R_{500}$ due to the use of $M_{YX}$ as a mass proxy, obtaining a radial profile of the scatter suppression factor. Taking into account this suppression factor, we estimated a scatter in scaled density profiles of approximately 40% at $R_{500}$.

- We quantified deviations from the average scaling with radius. These show no variation with mass, but which show a significant variation with redshift, in the sense that the core regions clearly evolve differently as compared to the bulk.

- We examined the scaled central density measured at $R = 0.05 R_{500}$ for the SZE-selected systems, finding that only the $z < 0.3$ sample is skewed. This skewness is positive, and may indicate the increased presence of centrally peaked systems at later times.

- We measured the scaled central density at $R = 0.015 R_{500}$ for the X-ray and SZE-selected systems at $z < 0.3$. The scaled central density of the local X-ray-selected sample exhibits two peaks. The main peak, corresponding to non-cool core systems in the X-ray-selected sample, is slightly offset to higher scaled central density from that of the local SZE-selected sample. The secondary peak in the X-ray-selected sample, corresponding to the cool core systems, is not seen in the SZE-selected sample, although the latter does exhibit a clear tail to higher scaled central density as confirmed by the strongly positively skewed distribution.

- The absolute value of the central density in the SZE-selected sample measured at 40 kpc does not appear to evolve with redshift, consistent with the findings of McDonald et al. (2017).

We then examined the relation between the core excised X-ray luminosity $L_{Xc}$ and the total mass derived from the $Y_{X,500}$ mass proxy, $M_{YX}$.

- This relation is extremely tight, with a logarithmic intrinsic scatter of $\sigma_{\ln L_{Xc}/M_{YX}} \le 0.15$ depending on sub-sample and band in which the luminosity is measured. Importantly, at low redshift, the best-fitting parameters of this relation do not depend on whether the sample was selected in X-rays or through the SZE, suggesting that $L_{Xc}$ is a selection-independent quantity.

- The slope of the bolometric relation fitted to the SZE-selected clusters, $B \sim 1.74 \pm 0.02$, is significantly steeper than self-similar. When left free to vary, the evolution of $n = -2.5 \pm 0.09$ is in agreement with the self-similar value of $-7/3$ within $<2\sigma$.

- We thoroughly examined the impact of selection bias and covariance on the relation. We found that the slope of the $L_{Xc} - M_{YX}$ relation is robust to selection effects due to intrinsic scatter in the $Y_{SZ}$ and $Y_X$ proxies, and to covariance between quantities. The dispersion, however, is very sensitive to the covariance. For reasonable values of covariance, we estimate that the measured dispersion in the $L_{Xc} - M$ relation is underestimated by a factor of at most ~1.7 due to the use of $Y_X$ as a mass proxy, implying a true scatter of $\sigma_{\ln L_{Xc}/M} \sim 0.22$.

- We show explicitly that the scatter in the $L_{Xc} - M$ relation can be accounted for almost entirely by object-to-object variations in gas density profiles.

With our study we have examined the mass and redshift depen-
dence of the ICM gas density profile, and made quantitative com-
parisons between X-ray- and SZE-selected samples. Our overall
conclusion is consistent with the view that the ICM bulk evolves

approximately self-similarly, with the core regions evolving sepa-
rately due to cooling and feedback from the central active galac-
tic nucleus. Indeed, it suggests potentially subtle differences in
the core regions between X-ray- and SZE-selected systems. It
also supports a view where the ICM gas mass fraction depends
on mass up to high redshift, with a dependence $f_{gas} \propto M_{500}^{0.22 \pm 0.02}$
for the present sample. Further progress can be undoubtedly be
made by bringing to bear fully independent mass estimates, such
as those that can be obtained from weak lensing and/or galaxy
velocity dispersions. Such studies are one of the goals of the
CHEX-MATE project (CHEX-MATE Collaboration 2021).

*Acknowledgements.* G.W.P., M.A., and J.-B.M. acknowledge funding from the European Research Council under the European Union's Seventh Framework Programme (FP72007-2013) ERC grant agreement no. 340519, and from the French space agency, CNES. B.J.M. acknowledges support from the Science and Technology Facilities Council (grant number ST/V000454/1). The results reported in this article are based on data obtained from the *XMM-Newton* observatory, an ESA science mission with instruments and contributions directly funded by ESA Member States and NASA.

References

Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706

Andrade-Santos, F., Jones, C., Forman, W. R., et al. 2017, ApJ, 843, 76

Andrade-Santos, F., Pratt, G. W., Melin, J.-B., et al. 2021, ApJ, 914, 58

Applegate, D. E., von der Linden, A., Kelly, P. L., et al. 2014, MNRAS, 439, 48

Arnaud, M., & Evrard, A. E. 1999, MNRAS, 305, 631

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893

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

Bartalucci, I., Arnaud, M., Pratt, G. W., et al. 2017, A&A, 598, A61

Bartalucci, I., Arnaud, M., Pratt, G. W., & Le Brun, A. M. C. 2018, A&A, 617, A64

Bartalucci, I., Arnaud, M., Pratt, G. W., Démoclès, J., & Lovisari, L. 2019, A&A, 628, A86

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

Böhringer, H., Schuecker, P., Pratt, G. W., et al. 2007, A&A, 469, 363

Bulbul, E., Chiu, I. N., Mohr, J. J., et al. 2019, ApJ, 871, 50

Cavaliere, A., & Fusco-Femiano, R. 1976, A&A, 500, 95

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

Connor, T., Donahue, M., Sun, M., et al. 2014, ApJ, 794, 48

Croston, J. H., Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2006, A&A, 459, 1007

Croston, J. H., Pratt, G. W., Böhringer, H., et al. 2008, A&A, 487, 431

Doane, D. P., & Seward, L. E. 2011, J. Stat. Educ., 19

Eckert, D., Roncarelli, M., Ettori, S., et al. 2015, MNRAS, 447, 2198

Ettori, S. 2015, MNRAS, 446, 2629

Fabian, A. C., Crawford, C. S., Edge, A. C., & Mushotzky, R. F. 1994, MNRAS, 267, 779

Farahi, A., Mulroy, S. L., Evrard, A. E., et al. 2019, Nat. Commun., 10, 2504

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Ghirardini, V., Eckert, D., Ettori, S., et al. 2019, A&A, 621, A41

Ghizzardi, S. 2001, in Flight calibration of the PSF for the MOS1 and MOS2 cameras, XMM-SOC-CAL-TN-0022

Groener, A. M., Goldberg, D. M., & Sereno, M. 2016, MNRAS, 455, 892

Hasselfield, M., Hilton, M., Marriage, T. A., et al. 2013, JCAP, 7, 008

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

Hoekstra, H., Mahdavi, A., Babul, A., & Bildfell, C. 2012, MNRAS, 427, 1298

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

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775

Kay, S. T., Peel, M. W., Short, C. J., et al. 2012, MNRAS, 422, 1999

Kelly, B. C. 2007, ApJ, 665, 1489

Kelly, P. L., von der Linden, A., Applegate, D. E., et al. 2014, MNRAS, 439, 28

Lau, E. T., Nagai, D., Avestruz, C., Nelson, K., & Vikhlinin, A. 2015, ApJ, 806, 68

Le Brun, A. M. C., McCarthy, I. G., Schaye, J., & Ponman, T. J. 2017, MNRAS, 466, 4442

Lovisari, L., Reiprich, T. H., & Schellenberger, G. 2015, A&A, 573, A118

Lovisari, L., Forman, W. R., Jones, C., et al. 2017, ApJ, 846, 51

Lovisari, L., Schellenberger, G., Sereno, M., et al. 2020, ApJ, 892, 102

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

Mantz, A. B., Allen, S. W., & Morris, R. G. 2016, MNRAS, 462, 681
Mantz, A. B., Allen, S. W., Morris, R. G., & von der Linden, A. 2018, MNRAS, 473, 3072

Maughan, B. J. 2007, ApJ, 668, 772

Maughan, B. J., Jones, C., Forman, W., & Van Speybroeck, L. 2008, ApJS, 174, 117

Maughan, B. J., Giles, P. A., Rines, K. J., et al. 2016, MNRAS, 461, 4182

McDonald, M., Allen, S. W., Bayliss, M., et al. 2017, ApJ, 843, 28

Melin, J. B., Bartlett, J. G., & Delabrouille, J. 2006, A&A, 459, 341

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

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

Nagarajan, A., Pacaud, F., Sommer, M., et al. 2019, MNRAS, 488, 1728

Neumann, D. M., & Arnaud, M. 1999, A&A, 348, 711

Okabe, N., & Smith, G. P. 2016, MNRAS, 461, 3794

Okabe, N., Takada, M., Umetsu, K., Futamase, T., & Smith, G. P. 2010, PASJ, 62, 811

Okabe, N., Smith, G. P., Umetsu, K., Takada, M., & Futamase, T. 2013, ApJ, 769, L35

Piffaretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & Melin, J.-B. 2011, A&A, 534, A109

Planck Collaboration VIII. 2011, A&A, 536, A8

Planck Collaboration IX. 2011, A&A, 536, A9

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

Planelles, S., Borgani, S., Fabjan, D., et al. 2014, MNRAS, 438, 195

Pratt, G. W., & Arnaud, M. 2002, A&A, 394, 375

Pratt, G. W., & Arnaud, M. 2003, A&A, 408, 1

Pratt, G. W., Croston, J. H., Arnaud, M., & Böhringer, H. 2009, A&A, 498, 361

Pratt, G. W., Arnaud, M., Piffaretti, R., et al. 2010, A&A, 511, A85

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

Robotham, A. S. G., & Obreschkow, D. 2015, PASA, 32, e033

Rossetti, M., Gastaldello, F., Ferioli, G., et al. 2016, MNRAS, 457, 4515

Rossetti, M., Gastaldello, F., Eckert, D., et al. 2017, MNRAS, 468, 1917

Rozo, E., Rykoff, E. S., Bartlett, J. G., & Evrard, A. 2014a, MNRAS, 438, 49

Rozo, E., Bartlett, J. G., Evrard, A. E., & Rykoff, E. S. 2014b, MNRAS, 438, 78

Rykoff, E. S., Evrard, A. E., McKay, T. A., et al. 2008, MNRAS, 387, L28

Schellenberger, G., & Reiprich, T. H. 2017, MNRAS, 469, 3738

Sereno, M. 2015, MNRAS, 450, 3665

Sereno, M., & Ettori, S. 2015, MNRAS, 450, 3633

Sereno, M., & Ettori, S. 2017, MNRAS, 468, 3322

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

Truong, N., Rasia, E., Mazzotta, P., et al. 2018, MNRAS, 474, 4089

Umetsu, K., Medezinski, E., Nonino, M., et al. 2014, ApJ, 795, 163

Umetsu, K., Zitrin, A., Gruen, D., et al. 2016, ApJ, 821, 116

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

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

von der Linden, A., Allen, M. T., Applegate, D. E., et al. 2014, MNRAS, 439, 2

Zhang, Y. Y., Finoguenov, A., Böhringer, H., et al. 2008, A&A, 482, 451
## Appendix A: Sample observation details

Table A.1 contains the sample observation details, including: cluster name(s), redshift, coordinates, column density, exposure time, and XMM-Newton OBSID used for the analysis.

<table><caption>**Table A.1.** Sample observation details. Columns: (1) Planck name; (2) SPT name; (3) ACT name (4) alternative name (non-exhaustive), REX-CESS clusters in bold; (5) redshift; (6) RA; (7) Dec.; (8) neutral hydrogen column density integrated along the line of sight determined from the LAB survey (Kalberla et al. 2005); (9) observation exposure time for MOS1 and pn detectors, in ks; (10) XMM-Newton OBSID.</caption><thead><tr><th>Planck name</th><th>SPT name</th><th>ACT name</th><th>Alt. name</th><th>z</th><th>RA</th><th>Dec.</th><th>N<sub>H</sub></th><th>Exp. MOS1-2, PN [ks]</th><th>Obs. Id</th></tr><tr><td></td><td></td><td></td><td></td><td></td><td>[J2000]</td><td>[J2000]</td><td>[10<sup>20</sup>cm<sup>-3</sup>]</td><td></td><td></td></tr></thead><tbody><tr><td>PSZ1 G246.01-51.76</td><td></td><td></td><td>RXC J2023.0-2056, AS0868</td><td>0.056</td><td>305.7450</td><td>-20.9485</td><td>5.60</td><td>17, 9</td><td>201902301</td></tr><tr><td></td><td></td><td></td><td>RXC J2157.4-0747, A2399</td><td>0.058</td><td>329.3673</td><td>-7.8046</td><td>3.50</td><td>10, 7</td><td>201902801</td></tr><tr><td></td><td></td><td></td><td>RXC J0345.7-4112, AS0384</td><td>0.060</td><td>56.4428</td><td>-41.2042</td><td>1.90</td><td>17, 8</td><td>201900801</td></tr><tr><td></td><td></td><td></td><td>RXC J0225.1-2928</td><td>0.060</td><td>36.2877</td><td>-29.4773</td><td>1.70</td><td>20, 17</td><td>302610601</td></tr><tr><td></td><td></td><td></td><td>RXC J1236.7-3354, AS0700</td><td>0.080</td><td>189.1712</td><td>-33.9260</td><td>5.60</td><td>24, 18</td><td>302610701</td></tr><tr><td></td><td></td><td></td><td>RXC J2129.8-5048, A3771</td><td>0.080</td><td>322.4271</td><td>-50.8167</td><td>2.20</td><td>23, 13</td><td>201902501</td></tr><tr><td></td><td></td><td></td><td>A2255</td><td>0.081</td><td>258.2029</td><td>64.0636</td><td>2.50</td><td>6, 4</td><td>011226081</td></tr><tr><td></td><td></td><td></td><td>RXC J0821.8 + 0112, A653</td><td>0.082</td><td>125.4614</td><td>1.1967</td><td>4.20</td><td>11, 7</td><td>201903601</td></tr><tr><td></td><td></td><td></td><td>A1651</td><td>0.084</td><td>194.8438</td><td>-4.1983</td><td>1.81</td><td>7, 4</td><td>020302011</td></tr><tr><td></td><td></td><td></td><td>A1650</td><td>0.084</td><td>194.6734</td><td>-1.7622</td><td>0.72</td><td>34, 28</td><td>009320011</td></tr><tr><td></td><td></td><td></td><td>RXC J1302.8-0230, A1663</td><td>0.085</td><td>195.7217</td><td>-2.51695</td><td>1.70</td><td>26, 16</td><td>201901801</td></tr><tr><td></td><td></td><td></td><td>RXC J0003.8 +0203, A2700</td><td>0.092</td><td>0.9572</td><td>2.0657</td><td>3.00</td><td>26, 19</td><td>201900101</td></tr><tr><td></td><td></td><td></td><td>A3921</td><td>0.094</td><td>342.4917</td><td>-64.4294</td><td>1.61</td><td>30, 23</td><td>011224011</td></tr><tr><td></td><td></td><td></td><td>A3911</td><td>0.097</td><td>341.5711</td><td>-52.7261</td><td>1.50</td><td>23, 10</td><td>014967031</td></tr><tr><td></td><td></td><td></td><td>RXC J2319.6 - 7313, A3992</td><td>0.098</td><td>349.9168</td><td>-73.2277</td><td>1.90</td><td>10, 6</td><td>201903301</td></tr><tr><td></td><td></td><td></td><td>A3827</td><td>0.099</td><td>330.4720</td><td>-59.9454</td><td>2.09</td><td>21, 11</td><td>014967011</td></tr><tr><td></td><td></td><td></td><td>RXC J0211.4 - 4017, A2984</td><td>0.101</td><td>32.8528</td><td>-40.2915</td><td>1.40</td><td>29, 22</td><td>201900601</td></tr><tr><td></td><td></td><td></td><td>RXC J0049.4 - 2931, AS0084</td><td>0.108</td><td>12.3459</td><td>-29.5206</td><td>1.80</td><td>20, 13</td><td>201900401</td></tr><tr><td></td><td></td><td></td><td>A2034</td><td>0.113</td><td>227.5528</td><td>33.5104</td><td>1.54</td><td>10, 6</td><td>030393011</td></tr><tr><td></td><td></td><td></td><td>RXC J0006.0 - 3443, A2721</td><td>0.115</td><td>1.5015</td><td>-34.7224</td><td>1.20</td><td>12, 6</td><td>201903801</td></tr><tr><td></td><td></td><td></td><td>RXC J0616.8 - 4748</td><td>0.116</td><td>94.2158</td><td>-47.7950</td><td>4.80</td><td>23, 19</td><td>302610401</td></tr><tr><td>PSZ2 G053.53+59.52</td><td></td><td></td><td>RXC J0145.0-5300, A2941</td><td>0.117</td><td>26.2433</td><td>-53.02085</td><td>2.30</td><td>1, 0</td><td>201900501</td></tr><tr><td>PSZ2 G352.28+77.66</td><td></td><td></td><td>RXC J1516.3 +0005, A2050</td><td>0.118</td><td>229.0747</td><td>0.0893</td><td>4.60</td><td>27, 21</td><td>201902001</td></tr><tr><td>PSZ2 G255.64-25.30</td><td></td><td></td><td>RXC J2149.1-3041, A3814</td><td>0.118</td><td>327.2817</td><td>-30.7013</td><td>2.30</td><td>25, 18</td><td>201902601</td></tr><tr><td>PSZ2 G285.52-62.23</td><td>SPT-CLJ0145-5301</td><td>ACT-CL J0145-5301</td><td>RXC J1141.4-1216, A1348</td><td>0.119</td><td>175.3513</td><td>-12.2776</td><td>3.30</td><td>28, 22</td><td>201901601</td></tr><tr><td></td><td></td><td></td><td>RXC J1516.5-0056, A2051</td><td>0.120</td><td>229.1842</td><td>-0.9696</td><td>5.50</td><td>29, 22</td><td>201902101</td></tr><tr><td></td><td></td><td></td><td>RXC J1044.5-0704, A1084</td><td>0.134</td><td>161.1370</td><td>-7.0687</td><td>3.40</td><td>26, 18</td><td>201901501</td></tr><tr><td></td><td></td><td></td><td>RXC J0605.8-3518, A3378</td><td>0.139</td><td>91.4752</td><td>-35.3023</td><td>4.02</td><td>14, 11</td><td>020190101</td></tr><tr><td></td><td></td><td></td><td>RXC J0020.7-2542, A0022</td><td>0.141</td><td>5.1755</td><td>-25.7080</td><td>2.3</td><td>15, 11</td><td>201900301</td></tr><tr><td></td><td></td><td></td><td>RXC J2218.6-3853, A3856</td><td>0.141</td><td>334.6677</td><td>-38.9018</td><td>1.13</td><td>14, 4</td><td>020190301</td></tr><tr><td></td><td></td><td></td><td>A1413</td><td>0.143</td><td>178.8250</td><td>23.4049</td><td>1.84</td><td>62, 42</td><td>050269021</td></tr><tr><td>PSZ2 G277.38+47.07</td><td></td><td></td><td>RXC J2048.1-1750, A2328</td><td>0.147</td><td>312.0419</td><td>-17.8413</td><td>4.80</td><td>25, 19</td><td>201902401</td></tr><tr><td>PSZ2 G000.04+45.13</td><td></td><td></td><td>RXC J0547.6-3152, A3364</td><td>0.148</td><td>86.9081</td><td>-31.8727</td><td>2.07</td><td>22, 14</td><td>201900901</td></tr><tr><td>PSZ2 G256.38+44.04</td><td></td><td></td><td>RXC J2217.7-3543, A3854</td><td>0.149</td><td>334.4400</td><td>-35.7260</td><td>1.20</td><td>22, 14</td><td>201902901</td></tr><tr><td>PSZ2 G241.79-24.01</td><td></td><td></td><td>RXC J2234.5-3744, A3888</td><td>0.151</td><td>338.6125</td><td>-37.7360</td><td>1.32</td><td>20, 13</td><td>040491081</td></tr><tr><td>PSZ2 G042.81-82.97</td><td></td><td></td><td>A2204</td><td>0.152</td><td>248.1959</td><td>5.5754</td><td>6.97</td><td>14, 8</td><td>030649021</td></tr><tr><td>PSZ2 G002.77-56.16</td><td></td><td></td><td>RXC J0945.4-0839, A0868</td><td>0.153</td><td>146.3575</td><td>-8.6557</td><td>3.59</td><td>8, 5</td><td>001754011</td></tr><tr><td>PSZ2 G226.18+76.79</td><td></td><td></td><td>RXC J2014.8-2430</td><td>0.154</td><td>303.7154</td><td>-24.5059</td><td>7.40</td><td>25, 16</td><td>201902201</td></tr><tr><td>PSZ2 G028.77-33.56</td><td></td><td></td><td>RX J1720.1+2638</td><td>0.164</td><td>260.0417</td><td>26.6250</td><td>5.65</td><td>14, 8</td><td>050067041</td></tr><tr><td>PSZ2 G236.92-26.65</td><td></td><td></td><td>RXC J0645.4-5413, A3404</td><td>0.164</td><td>101.3712</td><td>-54.2273</td><td>5.60</td><td>11, 7</td><td>040491041</td></tr><tr><td>PSZ2 G008.48-56.34</td><td></td><td></td><td>RXC J0958.3-1103, A0907</td><td>0.167</td><td>149.5929</td><td>-11.0644</td><td>5.10</td><td>9, 5</td><td>201903501</td></tr><tr><td>PSZ2 G009.39-59.41</td><td></td><td></td><td>A2218</td><td>0.171</td><td>248.9597</td><td>66.2125</td><td>2.60</td><td>17, 11</td><td>011298011</td></tr><tr><td>PSZ2 G021.10+33.24</td><td></td><td></td><td>A1914</td><td>0.171</td><td>216.5105</td><td>37.8243</td><td>1.06</td><td>15, 7</td><td>011223021</td></tr><tr><td>PSZ2 G244.71+32.50</td><td></td><td></td><td>A0665</td><td>0.182</td><td>127.7462</td><td>65.8398</td><td>4.24</td><td>5, 2</td><td>010989041</td></tr><tr><td>PSZ2 G018.32-28.50</td><td></td><td></td><td>RXC J1311.4-0120, A1689</td><td>0.183</td><td>197.8726</td><td>-1.3417</td><td>1.52</td><td>36, 27</td><td>009303011</td></tr><tr><td>PSZ2 G049.22+30.87</td><td></td><td></td><td>A520</td><td>0.203</td><td>73.5402</td><td>2.9212</td><td>3.74</td><td>20, 9</td><td>020151011</td></tr><tr><td>PSZ2 G263.68-22.55</td><td>SPT-CLJ0645-5413</td><td>ACT-CL J0645-5413</td><td>A2163</td><td>0.203</td><td>243.9399</td><td>-6.1491</td><td>16.50</td><td>10, 6</td><td>011223061</td></tr><tr><td>PSZ2 G249.38+33.26</td><td></td><td></td><td>A963</td><td>0.206</td><td>154.2653</td><td>39.0482</td><td>1.25</td><td>18, 12</td><td>008423071</td></tr><tr><td>PSZ2 G097.72+38.12</td><td></td><td></td><td>A773</td><td>0.217</td><td>139.4765</td><td>51.7315</td><td>1.28</td><td>13, 14</td><td>008423061</td></tr><tr><td>PSZ2 G067.17+67.46</td><td></td><td></td><td>A1763</td><td>0.223</td><td>203.8298</td><td>41.0001</td><td>0.94</td><td>12, 9</td><td>008423091</td></tr><tr><td>PSZ2 G149.75+34.68</td><td></td><td></td><td>A2261</td><td>0.224</td><td>260.6127</td><td>32.1324</td><td>3.19</td><td>2, 1</td><td>0093031001</td></tr><tr><td>PSZ2 G313.33+61.13</td><td></td><td></td><td>A2219</td><td>0.228</td><td>250.0837</td><td>46.7107</td><td>1.76</td><td>12, 6</td><td>060500051</td></tr><tr><td>PSZ2 G195.75-24.32</td><td></td><td></td><td>A2390</td><td>0.231</td><td>328.4031</td><td>17.6949</td><td>8.66</td><td>10, 8</td><td>011127011</td></tr><tr><td>PSZ2 G006.76+30.45</td><td></td><td></td><td>RXCJ0303.8-7752</td><td>0.274</td><td>45.9366</td><td>-77.8784</td><td>8.73</td><td>11, 8</td><td>020533011</td></tr><tr><td>PSZ2 G182.59+55.83</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PSZ2 G166.09+43.38</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PSZ2 G092.71+73.46</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PSZ2 G055.59+31.85</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PSZ2 G072.62+41.46</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PSZ2 G073.97+27.82</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PSZ2 G294.68-37.01</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>
<table><caption>Table A.1. continued.</caption><thead><tr><th>Planck name</th><th>SPT name</th><th>ACT name</th><th>Alt. name</th><th>z</th><th>RA</th><th>Dec.</th><th>N<sub>H</sub></th><th>Exp. MOS1-2, PN</th><th>Obs. Id</th></tr><tr><td></td><td></td><td></td><td></td><td></td><td>[J2000]</td><td>[J2000]</td><td>[10<sup>20</sup>cm<sup>-3</sup>]</td><td>[ks]</td><td></td></tr></thead><tbody><tr><td>PSZ2 G241.76-30.88</td><td></td><td></td><td>RXCJ0532.9-3701</td><td>0.275</td><td>83.2323</td><td>-37.0270</td><td>2.90</td><td>11, 6</td><td>004234181</td></tr><tr><td>PSZ2 G259.98-63.43</td><td>SPT-CLJ0232-4421</td><td></td><td>RXCJ0232.2-4420</td><td>0.284</td><td>38.0772</td><td>-44.3464</td><td>2.49</td><td>12, 7</td><td>004234031</td></tr><tr><td>PSZ2 G244.37-32.15</td><td></td><td></td><td>RXC J0528.9-3927</td><td>0.284</td><td>82.2211</td><td>-39.4714</td><td>2.13</td><td>7, 4</td><td>004234081</td></tr><tr><td>PSZ2 G106.87-83.23</td><td></td><td></td><td>RXC J0043.4-2037, A2813</td><td>0.292</td><td>10.8519</td><td>-20.6229</td><td>1.54</td><td>11, 5</td><td>004234021</td></tr><tr><td>PSZ2 G262.27-35.38</td><td>SPT-CL J0516-5430</td><td>ACT-CL J0516-5430</td><td>ACO S520</td><td>0.295</td><td>79.1536</td><td>-54.5089</td><td>2.56</td><td>3, 1</td><td>0042340701</td></tr><tr><td>PSZ2 G266.04-21.25</td><td>SPT-CLJ0658-5556</td><td>ACT-CL J0658-5557</td><td>1ES0657-558</td><td>0.296</td><td>104.6277</td><td>-55.9434</td><td>4.17</td><td>21, 14</td><td>011298021</td></tr><tr><td>PSZ2 G195.60+44.06</td><td></td><td></td><td>A781</td><td>0.298</td><td>140.1018</td><td>30.5028</td><td>1.94</td><td>57, 47</td><td>040117011</td></tr><tr><td>PSZ2 G125.71+53.86</td><td></td><td></td><td>A1576</td><td>0.302</td><td>189.2441</td><td>63.1871</td><td>1.68</td><td>6, 1</td><td>040225011</td></tr><tr><td>PSZ2 G008.94-81.22</td><td></td><td></td><td>A2744</td><td>0.307</td><td>3.5775</td><td>-30.3863</td><td>1.60</td><td>14, 10</td><td>004234011</td></tr><tr><td>PSZ2 G278.58+39.16</td><td></td><td></td><td>A1300</td><td>0.308</td><td>172.9775</td><td>-19.9285</td><td>4.50</td><td>11, 9</td><td>004234101</td></tr><tr><td>PSZ1 G103.58+24.78</td><td></td><td></td><td></td><td>0.334</td><td>286.1007</td><td>72.4622</td><td>7.48</td><td>28, 15</td><td>0693660501<br>0693663101</td></tr><tr><td>PSZ2 G349.46-59.95</td><td>SPT-CLJ2248-4431</td><td></td><td>AS 1063</td><td>0.347</td><td>342.1824</td><td>-44.5305</td><td>1.84</td><td>25, 15</td><td>050463011</td></tr><tr><td>PSZ2 G083.29-31.03</td><td></td><td></td><td>RXCJ2228+2037</td><td>0.412</td><td>337.1405</td><td>20.6204</td><td>4.26</td><td>24, 15</td><td>014789011</td></tr><tr><td>PSZ2 G284.41+52.45</td><td></td><td></td><td>MACSJ1206.2-0848</td><td>0.441</td><td>181.5521</td><td>-8.8002</td><td>4.35</td><td>30, 21</td><td>050243041</td></tr><tr><td>PSZ2 G056.93-55.08</td><td></td><td></td><td>MACSJ2243.3-0935</td><td>0.444</td><td>340.8387</td><td>-9.5947</td><td>3.11</td><td>103, 78</td><td>050349021</td></tr><tr><td>PSZ2 G254.08-58.45</td><td>SPT-CL J0304-4401</td><td></td><td>SMACS J0304.3-4402</td><td>0.458</td><td>46.0705</td><td>-44.0257</td><td>1.29</td><td>17, 12</td><td>0700182201</td></tr><tr><td>PSZ2 G265.10-59.50</td><td>SPT-CLJ0243-4833</td><td></td><td>RXCJ0243.6-4834</td><td>0.500</td><td>40.9129</td><td>-48.5611</td><td>2.15</td><td>15, 6</td><td>0672090501<br>0723780801</td></tr><tr><td>PSZ2 G044.77-51.30</td><td></td><td></td><td>MACSJ2214.9-1359</td><td>0.503</td><td>333.7383</td><td>-14.0045</td><td>2.88</td><td>17, 7</td><td>0693661901</td></tr><tr><td>PSZ2 G211.21+38.66</td><td></td><td></td><td>MACSJ0911.2+1746</td><td>0.505</td><td>137.7970</td><td>17.7760</td><td>3.28</td><td>33, 25</td><td>0693662501</td></tr><tr><td>PSZ2 G004.45-19.55</td><td></td><td></td><td></td><td>0.516</td><td>289.2692</td><td>-33.5228</td><td>5.90</td><td>14, 5</td><td>0656201001</td></tr><tr><td>PSZ2 G110.28-87.48</td><td></td><td></td><td></td><td>0.520</td><td>12.2939</td><td>-24.6792</td><td>1.50</td><td>26, 4</td><td>0693662101<br>0723780201</td></tr><tr><td>PSZ2 G212.44+63.19</td><td></td><td></td><td>RMJ105252.4+241530.0</td><td>0.529</td><td>163.2159</td><td>24.2584</td><td>1.89</td><td>34, 20</td><td>0693660701<br>0723780701</td></tr><tr><td>PSZ2 G201.50-27.31</td><td></td><td></td><td>MACSJ0454.1-0300</td><td>0.538</td><td>73.5471</td><td>-3.0162</td><td>3.92</td><td>22, 17</td><td>0205670101</td></tr><tr><td>PSZ2 G094.56+51.03</td><td></td><td></td><td>WHL J227.050+57.90</td><td>0.539</td><td>227.0821</td><td>57.9164</td><td>1.50</td><td>26, 20</td><td>0693660101<br>0723780501</td></tr><tr><td>PSZ2 G228.16+75.20</td><td></td><td></td><td>MACSJ1149.5+2223</td><td>0.544</td><td>177.3976</td><td>22.4011</td><td>1.92</td><td>13, 2</td><td>0693661701</td></tr><tr><td>PSZ2 G111.61-45.71</td><td></td><td></td><td>CL0016+16</td><td>0.546</td><td>4.6399</td><td>16.4362</td><td>3.99</td><td>33, 24</td><td>0111000101<br>0111000201</td></tr><tr><td>PSZ2 G180.25+21.03</td><td></td><td></td><td>MACSJ0717.5+3745</td><td>0.546</td><td>109.3800</td><td>37.7587</td><td>6.63</td><td>156, 116</td><td>0672420101<br>0672420201<br>0672420301</td></tr><tr><td>PSZ2 G183.90+42.99</td><td></td><td></td><td>WHL J137.713+38.83</td><td>0.559</td><td>137.7032</td><td>38.8357</td><td>1.63</td><td>14, 8</td><td>0723780101</td></tr><tr><td>PSZ2 G155.27-68.42</td><td></td><td></td><td>WHL J24.3324-8.477</td><td>0.567</td><td>24.3536</td><td>-8.4557</td><td>3.52</td><td>28, 18</td><td>0693662801<br>0700180201</td></tr><tr><td>PSZ2 G046.13+30.72</td><td></td><td></td><td>WHL J171705.5+240424</td><td>0.569</td><td>259.2742</td><td>24.0737</td><td>5.18</td><td>26, 2</td><td>0693661401</td></tr><tr><td>PSZ2 G239.93-39.97</td><td></td><td></td><td></td><td>0.580</td><td>71.6966</td><td>-37.0625</td><td>1.44</td><td>34, 23</td><td>0679181001<br>0693661201</td></tr><tr><td>PSZ2 G254.64-45.20</td><td>SPT-CLJ0417-4748</td><td></td><td></td><td>0.581</td><td>64.3464</td><td>-47.8134</td><td>1.34</td><td>20, 9</td><td>0700182401</td></tr><tr><td>PSZ2 G144.83+25.11</td><td></td><td></td><td>MACSJ20647.7+7015</td><td>0.584</td><td>101.9590</td><td>70.2481</td><td>5.40</td><td>72, 46</td><td>0551850401<br>0551851301</td></tr><tr><td>PSZ2 G045.32-38.46</td><td></td><td></td><td>MACSJ2129.4-0741</td><td>0.589</td><td>322.3591</td><td>-7.6913</td><td>4.32</td><td>9, 3</td><td>0700182001</td></tr><tr><td>PSZ2 G070.89+49.26</td><td></td><td></td><td>WHL J155625.2+444042</td><td>0.602</td><td>239.1098</td><td>44.6772</td><td>1.23</td><td>47, 25</td><td>0693661301</td></tr><tr><td>PSZ2 G045.87+57.70</td><td></td><td></td><td>WHL J151820.6+292740</td><td>0.609</td><td>229.5866</td><td>29.4603</td><td>2.12</td><td>23, 15</td><td>0693661101</td></tr><tr><td>PSZ2 G073.31+67.52</td><td></td><td></td><td>WHL J215.168+39.91</td><td>0.609</td><td>215.1709</td><td>39.9187</td><td>0.82</td><td>26, 19</td><td>0693661001</td></tr><tr><td>PSZ2 G099.86+58.45</td><td></td><td></td><td>WHL J213.697+54.78</td><td>0.615</td><td>213.6952</td><td>54.7840</td><td>1.50</td><td>21, 9</td><td>0693660601</td></tr></tbody></table>
<table><caption>Table A.1. continued.</caption><thead><tr><th>Planck name</th><th>SPT name</th><th>ACT name</th><th>Alt. name</th><th>z</th><th>RA</th><th>Dec.</th><th>N<sub>H</sub></th><th>Exp. MOS1-2, PN</th><th>Obs. Id</th></tr><tr><td></td><td></td><td></td><td></td><td></td><td>[J2000]</td><td>[J2000]</td><td>[10<sup>20</sup>cm<sup>-3</sup>]</td><td>[ks]</td><td></td></tr></thead><tbody><tr><td>PSZ2 G193.31-46.13</td><td></td><td></td><td></td><td>0.634</td><td>53.9644</td><td>-6.9758</td><td>4.15</td><td>67.48</td><td>0693662701<br>0723780301</td></tr><tr><td>PLCK G147.32-16.59</td><td></td><td></td><td></td><td>0.645</td><td>44.1056</td><td>40.2885</td><td>8.29</td><td>59.39</td><td>0658200401<br>0693661501</td></tr><tr><td>PLCK G260.7-26.3</td><td>SPT-CLJ0616-5227</td><td>ACT-CL J0616-5227</td><td></td><td>0.680</td><td>94.1429</td><td>-52.4518</td><td>4.25</td><td>35.20</td><td>0693662301</td></tr><tr><td>PSZ2 G097.52+51.70</td><td></td><td></td><td></td><td>0.700</td><td>223.8379</td><td>58.8718</td><td>1.06</td><td>24.17</td><td>0783881301</td></tr><tr><td>PSZ2 G219.89-34.39</td><td></td><td></td><td></td><td>0.700</td><td>73.6894</td><td>-20.2851</td><td>3.25</td><td>51.34</td><td>0679180501<br>0693660301</td></tr><tr><td>PSZ1 G080.66-57.87</td><td></td><td>ACT-CL J2327.4-0204</td><td>RCS2 J2327-0204</td><td>0.699</td><td>351.8654</td><td>-2.0771</td><td>4.62</td><td>25.75, 71</td><td>7355<br>14025<br>14361</td></tr><tr><td>PSZ2 G208.61-74.39</td><td></td><td></td><td></td><td>0.711</td><td>30.0695</td><td>-24.9132</td><td>1.37</td><td>47.28</td><td>0693662901<br>0723780601</td></tr><tr><td>PSZ1 G226.65+28.43</td><td></td><td></td><td></td><td>0.724</td><td>134.0858</td><td>1.7803</td><td>3.38</td><td>44.27</td><td>0783881001</td></tr><tr><td>PSZ2 G084.10+58.72</td><td></td><td></td><td></td><td>0.731</td><td>222.2535</td><td>48.5569</td><td>2.29</td><td>72.58</td><td>0783880901</td></tr><tr><td>PSZ2 G087.39+50.92</td><td></td><td></td><td></td><td>0.748</td><td>231.6379</td><td>54.1524</td><td>1.21</td><td>20.11</td><td>0783881201</td></tr><tr><td>PSZ2 G088.98+55.07</td><td></td><td></td><td></td><td>0.754</td><td>224.7448</td><td>52.8167</td><td>1.83</td><td>44.30</td><td>0783880801</td></tr><tr><td>PSZ2 G086.93+53.18</td><td></td><td></td><td></td><td>0.771</td><td>228.5022</td><td>52.8035</td><td>1.62</td><td>57.37</td><td>0783880701</td></tr><tr><td>PLCK G079.95+46.96</td><td></td><td></td><td></td><td>0.790</td><td>240.5461</td><td>51.0615</td><td>1.56</td><td>81.51</td><td>0783880601</td></tr><tr><td>PSZ2 G352.05-24.01</td><td></td><td></td><td></td><td>0.798</td><td>290.2490</td><td>-45.8500</td><td>6.03</td><td>52.28</td><td>0679180201<br>0693660401</td></tr><tr><td>PLCK G227.99+38.11</td><td></td><td></td><td></td><td>0.810</td><td>143.0891</td><td>5.6899</td><td>3.40</td><td>26.10</td><td>0783880501</td></tr><tr><td>PSZ2 G091.83+26.11</td><td></td><td></td><td></td><td>0.822</td><td>277.7925</td><td>62.2429</td><td>4.30</td><td>23.13</td><td>0762801001</td></tr><tr><td>PSZ2G 208.57-44.31</td><td></td><td></td><td></td><td>0.830</td><td>60.6477</td><td>-15.6783</td><td>2.13</td><td>60.47</td><td>0783880301</td></tr><tr><td>PSZ2 G071.82-56.55</td><td></td><td></td><td></td><td>0.870</td><td>347.3914</td><td>-4.1707</td><td>3.59</td><td>46.27</td><td>0783880201</td></tr><tr><td>PSZ2 G160.83+81.66</td><td></td><td>WARPS J1226.9+3332</td><td></td><td>0.892</td><td>186.7420</td><td>33.5463</td><td>1.83</td><td>63.47</td><td>0200340101</td></tr><tr><td>PSZ1 G254.58-32.16</td><td>SPT-CL J0535-4801<br>SPT-CLJ2146-4633</td><td></td><td></td><td>0.900</td><td>83.9566</td><td>-48.0247</td><td>3.11</td><td>56.30</td><td>0783880101</td></tr><tr><td></td><td></td><td></td><td></td><td>0.933</td><td>326.6447</td><td>-46.5475</td><td>1.64</td><td>153.102</td><td>0744400501<br>0744401301</td></tr><tr><td>PSZ2 G266.54-27.31</td><td>SPT-CLJ0615-5746<br>SPT-CLJ2341-5119</td><td></td><td></td><td>0.972</td><td>93.9660</td><td>-57.7796</td><td>4.32</td><td>12.3</td><td>0658200101</td></tr><tr><td></td><td></td><td></td><td></td><td>1.003</td><td>355.3010</td><td>-51.3286</td><td>1.21</td><td>91.45</td><td>0744400401<br>0763670201</td></tr><tr><td></td><td>SPT-CLJ0546-5345</td><td>ACT-CL J0546-5345</td><td></td><td>1.066</td><td>86.6551</td><td>-53.7596</td><td>6.79</td><td>127.113</td><td>0744400201<br>0744400301</td></tr><tr><td></td><td>SPT-CLJ2106-5844</td><td></td><td></td><td>1.132</td><td>316.5221</td><td>-58.7421</td><td>4.33</td><td>26.16</td><td>0763670301</td></tr></tbody></table>
## Appendix B: Sample data

Table B.1 contains the name, redshift, RA, Dec., mass from the $Y_X$ proxy, and soft-band and bolometric core-excised X-ray luminosity for the sample.

<table><caption>Table B.1. Sample observational data.</caption><thead><tr><th>Planck name</th><th>REXCESS name</th><th>z</th><th>RA [deg]</th><th>Dec. [deg]</th><th>M<sub>500</sub> [10<sup>14</sup> M<sub>☉</sub>]</th><th>L<sub>Xc</sub> [0.5-2] keV [10<sup>44</sup> erg s<sup>-1</sup>]</th><th>L<sub>Xc</sub> bol [10<sup>44</sup> erg s<sup>-1</sup>]</th></tr></thead><tbody><tr><td></td><td>RXC J2023.0 - 2056</td><td>0.056</td><td>305.7450</td><td>-20.9485</td><td>1.21<sup>+0.03</sup><sub>-0.03</sub></td><td>0.16<sup>+0.00</sup><sub>-0.00</sub></td><td>0.40<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td></td><td>RXC J2157.4 - 0747</td><td>0.058</td><td>329.3673</td><td>-7.8046</td><td>1.29<sup>+0.03</sup><sub>-0.03</sub></td><td>0.15<sup>+0.00</sup><sub>-0.00</sub></td><td>0.37<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td></td><td>RXC J0345.7 - 4112</td><td>0.060</td><td>56.4428</td><td>-41.2042</td><td>0.97<sup>+0.02</sup><sub>-0.02</sub></td><td>0.15<sup>+0.00</sup><sub>-0.00</sub></td><td>0.37<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td></td><td>RXC J0225.1 - 2928</td><td>0.060</td><td>36.2877</td><td>-29.4773</td><td>0.96<sup>+0.04</sup><sub>-0.04</sub></td><td>0.12<sup>+0.00</sup><sub>-0.00</sub></td><td>0.31<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td></td><td>RXC J1236.7 - 3354</td><td>0.080</td><td>189.1712</td><td>-33.9260</td><td>1.33<sup>+0.02</sup><sub>-0.02</sub></td><td>0.24<sup>+0.00</sup><sub>-0.00</sub></td><td>0.61<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td></td><td>RXC J2129.8 - 5048</td><td>0.080</td><td>322.4271</td><td>-50.8167</td><td>2.26<sup>+0.06</sup><sub>-0.06</sub></td><td>0.41<sup>+0.01</sup><sub>-0.01</sub></td><td>1.19<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>PSZ2 G093.92+34.92</td><td></td><td>0.081</td><td>258.2029</td><td>64.0636</td><td>5.00<sup>+0.21</sup><sub>-0.15</sub></td><td>1.24<sup>+0.01</sup><sub>-0.01</sub></td><td>4.78<sup>+0.10</sup><sub>-0.10</sub></td></tr><tr><td>PSZ2 G222.52+20.58</td><td>RXC J0821.8 + 0112</td><td>0.082</td><td>125.4614</td><td>1.1967</td><td>1.31<sup>+0.03</sup><sub>-0.04</sub></td><td>0.22<sup>+0.00</sup><sub>-0.00</sub></td><td>0.54<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td>PSZ2 G306.77+58.61</td><td></td><td>0.084</td><td>194.8438</td><td>-4.1983</td><td>4.51<sup>+0.12</sup><sub>-0.11</sub></td><td>1.48<sup>+0.01</sup><sub>-0.01</sub></td><td>5.02<sup>+0.07</sup><sub>-0.07</sub></td></tr><tr><td>PSZ2 G306.66+61.06</td><td></td><td>0.084</td><td>194.6734</td><td>-1.7622</td><td>4.21<sup>+0.05</sup><sub>-0.05</sub></td><td>1.27<sup>+0.00</sup><sub>-0.00</sub></td><td>4.25<sup>+0.03</sup><sub>-0.03</sub></td></tr><tr><td>PSZ2 G308.64+60.26</td><td>RXC J1302.8 - 0230</td><td>0.085</td><td>195.7216</td><td>-2.5170</td><td>1.89<sup>+0.03</sup><sub>-0.03</sub></td><td>0.31<sup>+0.00</sup><sub>-0.00</sub></td><td>0.83<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td>PSZ2 G099.57-58.64</td><td>RXC J0003.8 + 0203</td><td>0.092</td><td>0.9572</td><td>2.0657</td><td>2.11<sup>+0.04</sup><sub>-0.04</sub></td><td>0.40<sup>+0.00</sup><sub>-0.00</sub></td><td>1.15<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td>PSZ2 G321.98-47.96</td><td></td><td>0.094</td><td>342.4917</td><td>-64.4294</td><td>3.73<sup>+0.06</sup><sub>-0.06</sub></td><td>1.07<sup>+0.00</sup><sub>-0.00</sub></td><td>3.51<sup>+0.03</sup><sub>-0.03</sub></td></tr><tr><td>PSZ2 G336.60-55.43</td><td></td><td>0.097</td><td>341.5711</td><td>-52.7261</td><td>3.99<sup>+0.08</sup><sub>-0.08</sub></td><td>1.21<sup>+0.01</sup><sub>-0.01</sub></td><td>3.95<sup>+0.04</sup><sub>-0.04</sub></td></tr><tr><td>PSZ2 G311.62-42.31</td><td>RXC J2319.6 - 7313</td><td>0.098</td><td>349.9167</td><td>-73.2277</td><td>1.56<sup>+0.03</sup><sub>-0.03</sub></td><td>0.38<sup>+0.01</sup><sub>-0.01</sub></td><td>0.97<sup>+0.02</sup><sub>-0.01</sub></td></tr><tr><td>PSZ2 G332.23-46.37</td><td></td><td>0.099</td><td>330.4720</td><td>-59.9454</td><td>5.64<sup>+0.09</sup><sub>-0.09</sub></td><td>1.76<sup>+0.01</sup><sub>-0.01</sub></td><td>6.61<sup>+0.06</sup><sub>-0.06</sub></td></tr><tr><td></td><td>RXC J0211.4 - 4017</td><td>0.101</td><td>32.8528</td><td>-40.2915</td><td>1.00<sup>+0.02</sup><sub>-0.02</sub></td><td>0.20<sup>+0.00</sup><sub>-0.00</sub></td><td>0.48<sup>+0.01</sup><sub>-0.00</sub></td></tr><tr><td></td><td>RXC J0049.4 - 2931</td><td>0.108</td><td>12.3460</td><td>-29.5206</td><td>1.62<sup>+0.04</sup><sub>-0.04</sub></td><td>0.38<sup>+0.01</sup><sub>-0.01</sub></td><td>1.00<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>PSZ2 G053.53+59.52</td><td></td><td>0.113</td><td>227.5528</td><td>33.5104</td><td>5.64<sup>+0.14</sup><sub>-0.14</sub></td><td>1.71<sup>+0.01</sup><sub>-0.01</sub></td><td>6.34<sup>+0.10</sup><sub>-0.10</sub></td></tr><tr><td>PSZ2 G352.28-77.66</td><td>RXC J0006.0 - 3443</td><td>0.115</td><td>1.5015</td><td>-34.7224</td><td>3.95<sup>+0.12</sup><sub>-0.12</sub></td><td>0.97<sup>+0.01</sup><sub>-0.01</sub></td><td>3.17<sup>+0.05</sup><sub>-0.05</sub></td></tr><tr><td>PSZ2 G255.64-25.30</td><td>RXC J0616.8 - 4748</td><td>0.116</td><td>94.2158</td><td>-47.7950</td><td>2.70<sup>+0.05</sup><sub>-0.05</sub></td><td>0.60<sup>+0.01</sup><sub>-0.01</sub></td><td>1.88<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>PSZ2 G285.52-62.23</td><td>RXC J0145.0 - 5300</td><td>0.117</td><td>26.2433</td><td>-53.0208</td><td>4.37<sup>+0.08</sup><sub>-0.08</sub></td><td>1.08<sup>+0.01</sup><sub>-0.01</sub></td><td>3.87<sup>+0.03</sup><sub>-0.03</sub></td></tr><tr><td></td><td>RXC J1516.3 + 0005</td><td>0.118</td><td>229.0747</td><td>0.0893</td><td>3.28<sup>+0.07</sup><sub>-0.04</sub></td><td>0.89<sup>+0.00</sup><sub>-0.01</sub></td><td>2.76<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td></td><td>RXC J2149.1 - 3041</td><td>0.118</td><td>327.2817</td><td>-30.7013</td><td>2.25<sup>+0.03</sup><sub>-0.03</sub></td><td>0.55<sup>+0.00</sup><sub>-0.00</sub></td><td>1.58<sup>+0.02</sup><sub>-0.01</sub></td></tr><tr><td>PSZ2 G277.38+47.07</td><td>RXC J1141.4 - 1216</td><td>0.119</td><td>175.3513</td><td>-12.2776</td><td>2.27<sup>+0.02</sup><sub>-0.02</sub></td><td>0.60<sup>+0.00</sup><sub>-0.00</sub></td><td>1.69<sup>+0.01</sup><sub>-0.01</sub></td></tr><tr><td>PSZ2 G000.04+45.13</td><td>RXC J1516.5 - 0056</td><td>0.120</td><td>229.1842</td><td>-0.9696</td><td>2.59<sup>+0.05</sup><sub>-0.04</sub></td><td>0.62<sup>+0.01</sup><sub>-0.01</sub></td><td>1.77<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>PSZ2 G256.38+44.04</td><td>RXC J1044.5 - 0704</td><td>0.134</td><td>161.1370</td><td>-7.0687</td><td>2.69<sup>+0.02</sup><sub>-0.02</sub></td><td>1.04<sup>+0.01</sup><sub>-0.01</sub></td><td>2.99<sup>+0.02</sup><sub>-0.02</sub></td></tr><tr><td>PSZ2 G241.79-24.01</td><td>RXC J0605.8 - 3518</td><td>0.139</td><td>91.4752</td><td>-35.3023</td><td>3.89<sup>+0.09</sup><sub>-0.09</sub></td><td>1.27<sup>+0.01</sup><sub>-0.01</sub></td><td>4.23<sup>+0.16</sup><sub>-0.03</sub></td></tr><tr><td>PSZ2 G042.81-82.97</td><td>RXC J0020.7 - 2542</td><td>0.141</td><td>5.1755</td><td>-25.7080</td><td>3.84<sup>+0.06</sup><sub>-0.06</sub></td><td>1.17<sup>+0.01</sup><sub>-0.01</sub></td><td>4.05<sup>+0.04</sup><sub>-0.03</sub></td></tr><tr><td>PSZ2 G002.77-56.16</td><td>RXC J2218.6 - 3853</td><td>0.141</td><td>334.6677</td><td>-38.9018</td><td>4.92<sup>+0.11</sup><sub>-0.11</sub></td><td>1.50<sup>+0.02</sup><sub>-0.02</sub></td><td>5.56<sup>+0.06</sup><sub>-0.06</sub></td></tr><tr><td>PSZ2 G226.18+76.79</td><td></td><td>0.143</td><td>178.8250</td><td>23.4049</td><td>6.24<sup>+0.06</sup><sub>-0.06</sub></td><td>2.27<sup>+0.01</sup><sub>-0.01</sub></td><td>8.84<sup>+0.06</sup><sub>-0.06</sub></td></tr><tr><td>PSZ2 G028.77-33.56</td><td>RXC J2048.1 - 1750</td><td>0.147</td><td>312.0419</td><td>-17.8413</td><td>4.31<sup>+0.07</sup><sub>-0.07</sub></td><td>1.36<sup>+0.01</sup><sub>-0.01</sub></td><td>4.40<sup>+0.03</sup><sub>-0.03</sub></td></tr><tr><td>PSZ2 G236.92-26.65</td><td>RXC J0547.6 - 3152</td><td>0.148</td><td>86.9081</td><td>-31.8727</td><td>5.01<sup>+0.10</sup><sub>-0.10</sub></td><td>1.58<sup>+0.01</sup><sub>-0.01</sub></td><td>5.73<sup>+0.38</sup><sub>-0.04</sub></td></tr><tr><td>PSZ2 G008.47-56.34</td><td>RXC J2217.7 - 3543</td><td>0.149</td><td>334.4400</td><td>-35.7260</td><td>3.61<sup>+0.05</sup><sub>-0.05</sub></td><td>1.15<sup>+0.01</sup><sub>-0.01</sub></td><td>3.68<sup>+0.03</sup><sub>-0.03</sub></td></tr><tr><td>PSZ2 G003.93-59.41</td><td>RXC J2234.5 - 3744</td><td>0.151</td><td>338.6125</td><td>-37.7360</td><td>7.36<sup>+0.09</sup><sub>-0.09</sub></td><td>3.05<sup>+0.02</sup><sub>-0.02</sub></td><td>12.28<sup>+0.12</sup><sub>-0.10</sub></td></tr><tr><td>PSZ2 G021.10+33.24</td><td></td><td>0.152</td><td>248.1959</td><td>5.5754</td><td>8.06<sup>+0.23</sup><sub>-0.22</sub></td><td>2.96<sup>+0.02</sup><sub>-0.02</sub></td><td>12.51<sup>+0.17</sup><sub>-0.17</sub></td></tr><tr><td>PSZ2 G244.71+32.50</td><td>RXC J0945.4 - 0839</td><td>0.153</td><td>146.3575</td><td>-8.6557</td><td>3.91<sup>+0.15</sup><sub>-0.11</sub></td><td>1.43<sup>+0.02</sup><sub>-0.02</sub></td><td>4.66<sup>+0.08</sup><sub>-0.08</sub></td></tr><tr><td>PSZ2 G018.32-28.50</td><td>RXC J2014.8 - 2430</td><td>0.154</td><td>303.7154</td><td>-24.5059</td><td>5.38<sup>+0.07</sup><sub>-0.07</sub></td><td>2.08<sup>+0.02</sup><sub>-0.02</sub></td><td>7.47<sup>+0.07</sup><sub>-0.06</sub></td></tr><tr><td>PSZ2 G049.22+30.87</td><td></td><td>0.164</td><td>260.0417</td><td>26.6250</td><td>5.30<sup>+0.11</sup><sub>-0.11</sub></td><td>2.06<sup>+0.02</sup><sub>-0.02</sub></td><td>7.49<sup>+0.08</sup><sub>-0.08</sub></td></tr><tr><td>PSZ2 G263.68-22.55</td><td>RXC J0645.4 - 5413</td><td>0.164</td><td>101.3712</td><td>-54.2273</td><td>7.38<sup>+0.18</sup><sub>-0.18</sub></td><td>2.81<sup>+0.02</sup><sub>-0.02</sub></td><td>11.33<sup>+0.25</sup><sub>-0.08</sub></td></tr><tr><td>PSZ2 G249.38+33.26</td><td>RXC J0958.3 - 1103</td><td>0.167</td><td>149.5930</td><td>-11.0644</td><td>4.17<sup>+0.22</sup><sub>-0.15</sub></td><td>1.42<sup>+0.04</sup><sub>-0.04</sub></td><td>5.21<sup>+0.15</sup><sub>-0.14</sub></td></tr><tr><td>PSZ2 G097.72+38.12</td><td></td><td>0.171</td><td>248.9597</td><td>66.2125</td><td>5.37<sup>+0.10</sup><sub>-0.10</sub></td><td>2.20<sup>+0.01</sup><sub>-0.01</sub></td><td>7.76<sup>+0.09</sup><sub>-0.09</sub></td></tr></tbody></table>
<table><caption>Table B.1. continued.</caption><thead><tr><th>Planck name</th><th>REXCESS name</th><th>z</th><th>RA [deg]</th><th>Dec. [deg]</th><th>M<sub>500</sub> [10<sup>14</sup> M<sub>☉</sub>]</th><th>L<sub>Xc</sub> [0.5-2] keV [10<sup>44</sup> erg s<sup>-1</sup>]</th><th>L<sub>Xc</sub> bol [10<sup>44</sup> erg s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>PSZ2 G067.17+67.46</td><td></td><td>0.171</td><td>216.5105</td><td>37.8243</td><td>8.25<sup>+0.19</sup><sub>-0.19</sub></td><td>3.23<sup>+0.02</sup><sub>-0.02</sub></td><td>14.34<sup>+0.18</sup><sub>-0.18</sub></td></tr><tr><td>PSZ2 G149.75+34.68</td><td></td><td>0.182</td><td>127.7462</td><td>65.8398</td><td>8.07<sup>+0.53</sup><sub>-0.42</sub></td><td>3.04<sup>+0.06</sup><sub>-0.06</sub></td><td>13.15<sup>+0.41</sup><sub>-0.41</sub></td></tr><tr><td>PSZ2 G313.33+61.13</td><td>RXC J1311.4-0120</td><td>0.183</td><td>197.8727</td><td>-1.3415</td><td>8.41<sup>+0.08</sup><sub>-0.08</sub></td><td>3.36<sup>+0.01</sup><sub>-0.01</sub></td><td>14.93<sup>+0.08</sup><sub>-0.07</sub></td></tr><tr><td>PSZ2 G195.75-24.32</td><td></td><td>0.203</td><td>73.5402</td><td>2.9212</td><td>8.31<sup>+0.22</sup><sub>-0.21</sub></td><td>3.46<sup>+0.03</sup><sub>-0.03</sub></td><td>14.83<sup>+0.22</sup><sub>-0.22</sub></td></tr><tr><td>PSZ2 G006.76+30.45</td><td></td><td>0.203</td><td>243.9399</td><td>-6.1491</td><td>20.10<sup>+0.61</sup><sub>-0.62</sub></td><td>9.50<sup>+0.08</sup><sub>-0.08</sub></td><td>55.06<sup>+1.00</sup><sub>-1.00</sub></td></tr><tr><td>PSZ2 G182.59+55.83</td><td></td><td>0.206</td><td>154.2653</td><td>39.0482</td><td>5.02<sup>+0.10</sup><sub>-0.10</sub></td><td>2.15<sup>+0.02</sup><sub>-0.02</sub></td><td>7.74<sup>+0.10</sup><sub>-0.10</sub></td></tr><tr><td>PSZ2 G166.09+43.38</td><td></td><td>0.217</td><td>139.4765</td><td>51.7315</td><td>6.61<sup>+0.14</sup><sub>-0.14</sub></td><td>2.75<sup>+0.02</sup><sub>-0.02</sub></td><td>10.85<sup>+0.12</sup><sub>-0.12</sub></td></tr><tr><td>PSZ2 G092.71+73.46</td><td></td><td>0.223</td><td>203.8298</td><td>41.0001</td><td>7.38<sup>+0.18</sup><sub>-0.18</sub></td><td>3.58<sup>+0.02</sup><sub>-0.02</sub></td><td>13.63<sup>+0.16</sup><sub>-0.16</sub></td></tr><tr><td>PSZ2 G055.59+31.85</td><td></td><td>0.224</td><td>260.6127</td><td>32.1324</td><td>6.33<sup>+0.59</sup><sub>-0.51</sub></td><td>3.15<sup>+0.13</sup><sub>-0.13</sub></td><td>12.33<sup>+0.65</sup><sub>-0.65</sub></td></tr><tr><td>PSZ2 G072.62+41.46</td><td></td><td>0.228</td><td>250.0837</td><td>46.7107</td><td>11.62<sup>+0.33</sup><sub>-0.24</sub></td><td>6.31<sup>+0.05</sup><sub>-0.05</sub></td><td>30.96<sup>+0.54</sup><sub>-0.54</sub></td></tr><tr><td>PSZ2 G073.97-27.82</td><td></td><td>0.231</td><td>328.4031</td><td>17.6949</td><td>10.67<sup>+0.28</sup><sub>-0.28</sub></td><td>5.61<sup>+0.05</sup><sub>-0.05</sub></td><td>25.89<sup>+0.38</sup><sub>-0.38</sub></td></tr><tr><td>PSZ2 G294.68-37.01</td><td></td><td>0.274</td><td>45.9366</td><td>-77.8784</td><td>7.54<sup>+0.31</sup><sub>-0.32</sub></td><td>3.16<sup>+0.04</sup><sub>-0.04</sub></td><td>13.35<sup>+0.28</sup><sub>-0.28</sub></td></tr><tr><td>PSZ2 G241.76-30.88</td><td></td><td>0.275</td><td>83.2323</td><td>-37.0270</td><td>6.40<sup>+0.23</sup><sub>-0.23</sub></td><td>2.96<sup>+0.04</sup><sub>-0.04</sub></td><td>11.77<sup>+0.23</sup><sub>-0.23</sub></td></tr><tr><td>PSZ2 G259.98-63.43</td><td></td><td>0.284</td><td>38.0772</td><td>-44.3464</td><td>6.99<sup>+0.21</sup><sub>-0.20</sub></td><td>4.34<sup>+0.04</sup><sub>-0.04</sub></td><td>16.28<sup>+0.27</sup><sub>-0.27</sub></td></tr><tr><td>PSZ2 G244.37-32.15</td><td></td><td>0.284</td><td>82.2211</td><td>-39.4714</td><td>6.86<sup>+0.30</sup><sub>-0.33</sub></td><td>4.18<sup>+0.05</sup><sub>-0.05</sub></td><td>15.29<sup>+0.37</sup><sub>-0.37</sub></td></tr><tr><td>PSZ2 G106.87-83.23</td><td></td><td>0.292</td><td>10.8519</td><td>-20.6229</td><td>5.91<sup>+0.25</sup><sub>-0.18</sub></td><td>3.38<sup>+0.04</sup><sub>-0.04</sub></td><td>12.36<sup>+0.25</sup><sub>-0.25</sub></td></tr><tr><td>PSZ2 G262.27-35.38</td><td></td><td>0.295</td><td>79.1536</td><td>-54.5090</td><td>7.81<sup>+0.81</sup><sub>-0.73</sub></td><td>3.90<sup>+0.17</sup><sub>-0.17</sub></td><td>15.98<sup>+0.95</sup><sub>-0.95</sub></td></tr><tr><td>PSZ2 G266.04-21.25</td><td></td><td>0.296</td><td>104.6277</td><td>-55.9434</td><td>13.92<sup>+0.29</sup><sub>-0.28</sub></td><td>8.72<sup>+0.06</sup><sub>-0.06</sub></td><td>46.88<sup>+0.53</sup><sub>-0.53</sub></td></tr><tr><td>PSZ2 G195.60+44.06</td><td></td><td>0.298</td><td>140.1018</td><td>30.5028</td><td>5.36<sup>+0.09</sup><sub>-0.09</sub></td><td>2.63<sup>+0.02</sup><sub>-0.02</sub></td><td>9.46<sup>+0.11</sup><sub>-0.11</sub></td></tr><tr><td>PSZ2 G125.71+53.86</td><td></td><td>0.302</td><td>189.2441</td><td>63.1871</td><td>5.67<sup>+0.44</sup><sub>-0.39</sub></td><td>2.99<sup>+0.09</sup><sub>-0.09</sub></td><td>11.43<sup>+0.53</sup><sub>-0.53</sub></td></tr><tr><td>PSZ2 G008.94-81.22</td><td></td><td>0.307</td><td>3.5775</td><td>-30.3863</td><td>9.81<sup>+0.24</sup><sub>-0.24</sub></td><td>6.70<sup>+0.05</sup><sub>-0.05</sub></td><td>27.83<sup>+0.32</sup><sub>-0.32</sub></td></tr><tr><td>PSZ2 G278.58+39.16</td><td></td><td>0.308</td><td>172.9775</td><td>-19.9285</td><td>8.61<sup>+0.33</sup><sub>-0.30</sub></td><td>4.68<sup>+0.05</sup><sub>-0.05</sub></td><td>18.93<sup>+0.31</sup><sub>-0.31</sub></td></tr><tr><td>PSZ1 G103.58+24.78</td><td></td><td>0.330</td><td>286.1007</td><td>72.4622</td><td>4.89<sup>+0.24</sup><sub>-0.23</sub></td><td>1.45<sup>+0.03</sup><sub>-0.03</sub></td><td>6.09<sup>+0.15</sup><sub>-0.15</sub></td></tr><tr><td>PSZ2 G349.46-59.95</td><td></td><td>0.347</td><td>342.1824</td><td>-44.5305</td><td>12.66<sup>+0.27</sup><sub>-0.27</sub></td><td>8.49<sup>+0.06</sup><sub>-0.06</sub></td><td>44.19<sup>+0.63</sup><sub>-0.63</sub></td></tr><tr><td>PSZ2 G083.29-31.03</td><td></td><td>0.412</td><td>337.1405</td><td>20.6204</td><td>8.70<sup>+0.28</sup><sub>-0.28</sub></td><td>5.24<sup>+0.05</sup><sub>-0.05</sub></td><td>23.01<sup>+0.43</sup><sub>-0.43</sub></td></tr><tr><td>PSZ2 G284.41+52.45</td><td></td><td>0.441</td><td>181.5521</td><td>-8.8002</td><td>10.76<sup>+0.30</sup><sub>-0.30</sub></td><td>6.79<sup>+0.06</sup><sub>-0.06</sub></td><td>33.82<sup>+0.65</sup><sub>-0.65</sub></td></tr><tr><td>PSZ2 G056.93-55.08</td><td></td><td>0.444</td><td>340.8387</td><td>-9.5947</td><td>9.10<sup>+0.12</sup><sub>-0.12</sub></td><td>7.30<sup>+0.03</sup><sub>-0.03</sub></td><td>31.65<sup>+0.26</sup><sub>-0.26</sub></td></tr><tr><td>PSZ2 G254.08-58.45</td><td></td><td>0.458</td><td>46.0705</td><td>-44.0257</td><td>5.84<sup>+0.33</sup><sub>-0.33</sub></td><td>3.19<sup>+0.06</sup><sub>-0.06</sub></td><td>11.74<sup>+0.32</sup><sub>-0.32</sub></td></tr><tr><td>PSZ2 G265.10-59.50</td><td></td><td>0.500</td><td>40.9129</td><td>-48.5611</td><td>6.13<sup>+0.63</sup><sub>-0.45</sub></td><td>3.10<sup>+0.08</sup><sub>-0.08</sub></td><td>12.46<sup>+0.57</sup><sub>-0.57</sub></td></tr><tr><td>PSZ2 G044.77-51.30</td><td></td><td>0.503</td><td>333.7383</td><td>-14.0045</td><td>7.95<sup>+0.44</sup><sub>-0.43</sub></td><td>5.29<sup>+0.08</sup><sub>-0.08</sub></td><td>22.64<sup>+0.66</sup><sub>-0.66</sub></td></tr><tr><td>PSZ2 G211.21+38.66</td><td></td><td>0.505</td><td>137.7970</td><td>17.7760</td><td>5.48<sup>+0.22</sup><sub>-0.22</sub></td><td>3.03<sup>+0.04</sup><sub>-0.04</sub></td><td>11.51<sup>+0.26</sup><sub>-0.26</sub></td></tr><tr><td>PSZ2 G004.45-19.55</td><td></td><td>0.516</td><td>289.2692</td><td>-33.5228</td><td>8.73<sup>+0.48</sup><sub>-0.48</sub></td><td>6.28<sup>+0.15</sup><sub>-0.15</sub></td><td>28.93<sup>+1.05</sup><sub>-1.05</sub></td></tr><tr><td>PSZ2 G110.28-87.48</td><td></td><td>0.520</td><td>12.2939</td><td>-24.6792</td><td>4.83<sup>+0.46</sup><sub>-0.36</sub></td><td>2.30<sup>+0.08</sup><sub>-0.08</sub></td><td>8.25<sup>+0.38</sup><sub>-0.38</sub></td></tr><tr><td>PSZ2 G212.44+63.19</td><td></td><td>0.529</td><td>163.2159</td><td>24.2584</td><td>4.15<sup>+0.23</sup><sub>-0.23</sub></td><td>1.93<sup>+0.04</sup><sub>-0.04</sub></td><td>7.08<sup>+0.25</sup><sub>-0.25</sub></td></tr><tr><td>PSZ2 G201.50-27.31</td><td></td><td>0.538</td><td>73.5471</td><td>-3.0162</td><td>7.90<sup>+0.30</sup><sub>-0.29</sub></td><td>6.47<sup>+0.07</sup><sub>-0.07</sub></td><td>28.31<sup>+0.57</sup><sub>-0.57</sub></td></tr><tr><td>PSZ2 G094.56+51.03</td><td></td><td>0.539</td><td>227.0821</td><td>57.9164</td><td>6.15<sup>+0.25</sup><sub>-0.24</sub></td><td>3.87<sup>+0.05</sup><sub>-0.05</sub></td><td>14.86<sup>+0.34</sup><sub>-0.34</sub></td></tr><tr><td>PSZ2 G228.16+75.20</td><td></td><td>0.544</td><td>177.3976</td><td>22.4011</td><td>9.36<sup>+0.64</sup><sub>-0.62</sub></td><td>7.08<sup>+0.17</sup><sub>-0.17</sub></td><td>31.46<sup>+1.15</sup><sub>-1.15</sub></td></tr><tr><td>PSZ2 G180.25+21.03</td><td></td><td>0.546</td><td>109.3800</td><td>37.7587</td><td>12.83<sup>+0.17</sup><sub>-0.17</sub></td><td>11.41<sup>+0.04</sup><sub>-0.04</sub></td><td>56.07<sup>+0.41</sup><sub>-0.41</sub></td></tr><tr><td>PSZ2 G111.61-45.71</td><td></td><td>0.546</td><td>4.6399</td><td>16.4362</td><td>9.21<sup>+0.24</sup><sub>-0.24</sub></td><td>7.49<sup>+0.07</sup><sub>-0.07</sub></td><td>34.24<sup>+0.52</sup><sub>-0.52</sub></td></tr><tr><td>PSZ2 G183.90+42.99</td><td></td><td>0.559</td><td>137.7032</td><td>38.8357</td><td>8.44<sup>+0.60</sup><sub>-0.53</sub></td><td>4.79<sup>+0.10</sup><sub>-0.10</sub></td><td>22.54<sup>+0.78</sup><sub>-0.78</sub></td></tr><tr><td>PSZ2 G155.27-68.42</td><td></td><td>0.567</td><td>24.3536</td><td>-8.4557</td><td>8.01<sup>+0.46</sup><sub>-0.38</sub></td><td>4.20<sup>+0.06</sup><sub>-0.06</sub></td><td>18.41<sup>+0.51</sup><sub>-0.51</sub></td></tr><tr><td>PSZ2 G046.13+30.72</td><td></td><td>0.569</td><td>259.2742</td><td>24.0737</td><td>3.17<sup>+0.22</sup><sub>-0.22</sub></td><td>1.26<sup>+0.04</sup><sub>-0.04</sub></td><td>4.40<sup>+0.20</sup><sub>-0.20</sub></td></tr><tr><td>PSZ2 G239.93-39.97</td><td></td><td>0.580</td><td>71.6966</td><td>-37.0625</td><td>5.73<sup>+0.26</sup><sub>-0.23</sub></td><td>3.36<sup>+0.04</sup><sub>-0.04</sub></td><td>13.89<sup>+0.31</sup><sub>-0.31</sub></td></tr><tr><td>PSZ2 G254.64-45.20</td><td></td><td>0.581</td><td>64.3464</td><td>-47.8134</td><td>4.97<sup>+0.27</sup><sub>-0.27</sub></td><td>3.02<sup>+0.06</sup><sub>-0.06</sub></td><td>11.69<sup>+0.40</sup><sub>-0.40</sub></td></tr><tr><td>PSZ2 G144.83+25.11</td><td></td><td>0.584</td><td>101.9590</td><td>70.2481</td><td>7.78<sup>+0.21</sup><sub>-0.20</sub></td><td>4.80<sup>+0.04</sup><sub>-0.04</sub></td><td>22.72<sup>+0.34</sup><sub>-0.34</sub></td></tr><tr><td>PSZ2 G045.32-38.46</td><td></td><td>0.589</td><td>322.3591</td><td>-7.6913</td><td>7.39<sup>+0.66</sup><sub>-0.65</sub></td><td>4.95<sup>+0.16</sup><sub>-0.16</sub></td><td>22.33<sup>+1.18</sup><sub>-1.18</sub></td></tr><tr><td>PSZ2 G070.89+49.26</td><td></td><td>0.602</td><td>239.1098</td><td>44.6772</td><td>5.02<sup>+0.20</sup><sub>-0.21</sub></td><td>2.93<sup>+0.05</sup><sub>-0.05</sub></td><td>11.71<sup>+0.31</sup><sub>-0.31</sub></td></tr><tr><td>PSZ2 G045.87+57.70</td><td></td><td>0.609</td><td>229.5866</td><td>29.4603</td><td>5.82<sup>+0.22</sup><sub>-0.22</sub></td><td>5.07<sup>+0.06</sup><sub>-0.06</sub></td><td>19.56<sup>+0.35</sup><sub>-0.35</sub></td></tr><tr><td>PSZ2 G073.31+67.52</td><td></td><td>0.609</td><td>215.1709</td><td>39.9187</td><td>6.15<sup>+0.26</sup><sub>-0.25</sub></td><td>4.34<sup>+0.07</sup><sub>-0.07</sub></td><td>17.97<sup>+0.47</sup><sub>-0.47</sub></td></tr><tr><td>PSZ2 G099.86+58.45</td><td></td><td>0.615</td><td>213.6952</td><td>54.7840</td><td>7.09<sup>+0.42</sup><sub>-0.42</sub></td><td>5.01<sup>+0.09</sup><sub>-0.09</sub></td><td>22.03<sup>+0.79</sup><sub>-0.79</sub></td></tr><tr><td>PSZ2 G193.31-46.13</td><td></td><td>0.620</td><td>53.9644</td><td>-6.9758</td><td>5.49<sup>+0.30</sup><sub>-0.32</sub></td><td>3.25<sup>+0.04</sup><sub>-0.04</sub></td><td>12.88<sup>+0.26</sup><sub>-0.26</sub></td></tr></tbody></table>
<table><caption>Table B.1. continued.</caption><thead><tr><th>Planck name</th><th>REXCESS name</th><th>z</th><th>RA [deg]</th><th>Dec. [deg]</th><th>M<sub>500</sub> [10<sup>14</sup> M<sub>☉</sub>]</th><th>L<sub>Xc</sub> [0.5-2] keV [10<sup>44</sup> erg s<sup>-1</sup>]</th><th>L<sub>Xc</sub> bol [10<sup>44</sup> erg s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>PLCK G147.32-16.59</td><td></td><td>0.645</td><td>44.1056</td><td>40.2885</td><td>6.51<sup>+0.29</sup><sub>-0.28</sub></td><td>4.26<sup>+0.05</sup><sub>-0.05</sub></td><td>17.69<sup>+0.38</sup><sub>-0.38</sub></td></tr><tr><td>PLCK G260.7 - 26.3</td><td></td><td>0.680</td><td>94.1429</td><td>-52.4518</td><td>4.56<sup>+0.34</sup><sub>-0.32</sub></td><td>2.86<sup>+0.06</sup><sub>-0.06</sub></td><td>10.57<sup>+0.35</sup><sub>-0.35</sub></td></tr><tr><td>PSZ2 G097.52 + 51.70</td><td></td><td>0.700</td><td>223.8379</td><td>58.8718</td><td>4.08<sup>+0.23</sup><sub>-0.23</sub></td><td>3.17<sup>+0.07</sup><sub>-0.07</sub></td><td>10.63<sup>+0.39</sup><sub>-0.39</sub></td></tr><tr><td>PSZ2 G219.89-34.39</td><td></td><td>0.700</td><td>73.6894</td><td>-20.2851</td><td>6.77<sup>+0.33</sup><sub>-0.29</sub></td><td>4.04<sup>+0.06</sup><sub>-0.06</sub></td><td>19.31<sup>+0.53</sup><sub>-0.53</sub></td></tr><tr><td>PSZ1 G080.66 - 57.87</td><td></td><td>0.700</td><td>351.8654</td><td>-2.0771</td><td>8.49<sup>+0.70</sup><sub>-0.55</sub></td><td>10.56<sup>+0.22</sup><sub>-0.22</sub></td><td>43.59<sup>+1.53</sup><sub>-1.53</sub></td></tr><tr><td>PSZ2 G208.61-74.39</td><td></td><td>0.711</td><td>30.0695</td><td>-24.9132</td><td>5.23<sup>+0.23</sup><sub>-0.23</sub></td><td>3.83<sup>+0.06</sup><sub>-0.06</sub></td><td>15.88<sup>+0.37</sup><sub>-0.37</sub></td></tr><tr><td>PSZ1 G226.65 + 28.43</td><td></td><td>0.724</td><td>134.0858</td><td>1.7803</td><td>3.44<sup>+0.20</sup><sub>-0.20</sub></td><td>2.55<sup>+0.04</sup><sub>-0.04</sub></td><td>8.39<sup>+0.25</sup><sub>-0.25</sub></td></tr><tr><td>PSZ2 G084.10 + 58.72</td><td></td><td>0.731</td><td>222.2535</td><td>48.5569</td><td>2.93<sup>+0.17</sup><sub>-0.18</sub></td><td>1.84<sup>+0.03</sup><sub>-0.03</sub></td><td>5.65<sup>+0.16</sup><sub>-0.16</sub></td></tr><tr><td>PSZ2 G087.39 + 50.92</td><td></td><td>0.748</td><td>231.6379</td><td>54.1523</td><td>3.29<sup>+0.30</sup><sub>-0.29</sub></td><td>2.41<sup>+0.07</sup><sub>-0.07</sub></td><td>7.55<sup>+0.34</sup><sub>-0.34</sub></td></tr><tr><td>PSZ2 G088.98 + 55.07</td><td></td><td>0.754</td><td>224.7448</td><td>52.8167</td><td>1.46<sup>+0.23</sup><sub>-0.20</sub></td><td>0.54<sup>+0.03</sup><sub>-0.03</sub></td><td>1.57<sup>+0.12</sup><sub>-0.12</sub></td></tr><tr><td>PSZ2 G086.93 + 53.18</td><td></td><td>0.771</td><td>228.5022</td><td>52.8035</td><td>5.21<sup>+0.19</sup><sub>-0.19</sub></td><td>4.63<sup>+0.06</sup><sub>-0.06</sub></td><td>17.55<sup>+0.42</sup><sub>-0.42</sub></td></tr><tr><td>PLCK G079.95 + 46.96</td><td></td><td>0.790</td><td>240.5461</td><td>51.0615</td><td>1.99<sup>+0.23</sup><sub>-0.20</sub></td><td>0.98<sup>+0.03</sup><sub>-0.03</sub></td><td>2.64<sup>+0.11</sup><sub>-0.11</sub></td></tr><tr><td>PSZ2 G352.05-24.01</td><td></td><td>0.798</td><td>290.2490</td><td>-45.8500</td><td>5.50<sup>+0.32</sup><sub>-0.31</sub></td><td>4.40<sup>+0.07</sup><sub>-0.07</sub></td><td>20.12<sup>+0.51</sup><sub>-0.51</sub></td></tr><tr><td>PLCK G227.99 + 38.11</td><td></td><td>0.810</td><td>143.0891</td><td>5.6900</td><td>2.55<sup>+0.59</sup><sub>-0.41</sub></td><td>1.15<sup>+0.07</sup><sub>-0.07</sub></td><td>4.15<sup>+0.48</sup><sub>-0.48</sub></td></tr><tr><td>PSZ2 G091.83+26.11</td><td></td><td>0.816</td><td>277.7924</td><td>62.2429</td><td>9.83<sup>+0.52</sup><sub>-0.50</sub></td><td>10.62<sup>+0.20</sup><sub>-0.20</sub></td><td>53.33<sup>+1.87</sup><sub>-1.87</sub></td></tr><tr><td>PSZ2 G208.57 - 44.31</td><td></td><td>0.830</td><td>60.6477</td><td>-15.6784</td><td>0.65<sup>+0.22</sup><sub>-0.14</sub></td><td>0.35<sup>+0.04</sup><sub>-0.04</sub></td><td>0.66<sup>+0.07</sup><sub>-0.07</sub></td></tr><tr><td>PSZ2 G071.82 - 56.55</td><td></td><td>0.870</td><td>347.3914</td><td>-4.1707</td><td>4.42<sup>+0.27</sup><sub>-0.21</sub></td><td>4.96<sup>+0.09</sup><sub>-0.09</sub></td><td>17.64<sup>+0.50</sup><sub>-0.50</sub></td></tr><tr><td>PSZ2 G160.83 + 81.66</td><td></td><td>0.890</td><td>186.7420</td><td>33.5463</td><td>5.95<sup>+0.23</sup><sub>-0.22</sub></td><td>5.06<sup>+0.07</sup><sub>-0.07</sub></td><td>22.03<sup>+0.49</sup><sub>-0.49</sub></td></tr><tr><td>PSZ1 G254.58 - 32.16</td><td></td><td>0.900</td><td>83.9566</td><td>-48.0247</td><td>3.80<sup>+0.31</sup><sub>-0.31</sub></td><td>2.38<sup>+0.06</sup><sub>-0.06</sub></td><td>8.59<sup>+0.32</sup><sub>-0.32</sub></td></tr><tr><td>SPT-CL J2146 - 4633</td><td></td><td>0.933</td><td>326.6447</td><td>-46.5475</td><td>3.24<sup>+0.15</sup><sub>-0.16</sub></td><td>2.98<sup>+0.04</sup><sub>-0.04</sub></td><td>9.75<sup>+0.23</sup><sub>-0.23</sub></td></tr><tr><td>PLCK G266.6 - 27.3</td><td></td><td>0.972</td><td>93.9660</td><td>-57.7796</td><td>8.93<sup>+0.40</sup><sub>-0.40</sub></td><td>10.27<sup>+0.18</sup><sub>-0.18</sub></td><td>54.62<sup>+1.65</sup><sub>-1.65</sub></td></tr><tr><td>SPT-CL J2341 - 5119</td><td></td><td>1.003</td><td>355.3010</td><td>-51.3286</td><td>4.17<sup>+0.21</sup><sub>-0.20</sub></td><td>3.68<sup>+0.07</sup><sub>-0.07</sub></td><td>15.06<sup>+0.43</sup><sub>-0.43</sub></td></tr><tr><td>SPT-CL J0546 - 5345</td><td></td><td>1.066</td><td>86.6551</td><td>-53.7596</td><td>4.33<sup>+0.21</sup><sub>-0.19</sub></td><td>3.67<sup>+0.05</sup><sub>-0.05</sub></td><td>15.25<sup>+0.32</sup><sub>-0.32</sub></td></tr><tr><td>SPT-CL J2106 - 5844</td><td></td><td>1.132</td><td>316.5221</td><td>-58.7421</td><td>7.01<sup>+0.49</sup><sub>-0.48</sub></td><td>9.33<sup>+0.23</sup><sub>-0.23</sub></td><td>47.17<sup>+1.95</sup><sub>-1.95</sub></td></tr></tbody></table>
## Appendix C: Sample selection bias and covariance tests

We used similar simulations as those produced for the Andrade-Santos et al. (2021) study. Simulated clusters, modelled with the Arnaud et al. (2010) pressure profile and drawn from a Tinker mass function (Tinker et al. 2008), were injected into the Planck Early SZ maps in the ‘cosmological’ mask region. The $Y_{SZ}$ value for each object was drawn from the $Y_{SZ} - M$ relation of Arnaud et al. (2010), with a bias between the X-ray calibrated mass and the true mass of $(1-b) = 0.65$, a value used to obtain the observed cluster number counts in the Planck cosmology. As the slope of the $Y_{SZ} - Y_X$ relation is expected to be close to unity, and we are only interested in slope variations, we assume that the normalisation and slope of the $Y_X - M$ and $Y_{SZ} - M$ relations are the same.

We draw the $Y_{SZ}$, $Y_X$ and $L_X$ quantities associated to each simulated cluster following a correlated Gaussian distribution with covariance matrix. If $Q_M$ is the latent value of $Q$ at a given mass, obtained if there were no scatter, one can write

$$Y_M = A_Y E(z)^{2/3} M^{B_Y} \qquad (C.1)$$

$$L_M = A_L E(z)^2 M^{B_L} \qquad (C.2)$$

$$P(D_A^2 Y_{SZ}, Y_X, L_{Xc} | Y_M, Y_M, L_M) = N[(Y_M, Y_M, L_M), V_\sigma], \qquad (C.3)$$

where $Y_M$ and $L_M$ are the latent values of $Y$ and $L$ at a given mass (the values obtained if there were no scatter), and $D_A^2 Y_{SZ}$, $Y_X$, and $L_{Xc}$ are the true values. For these simulations we assume $B_Y = 1.79$ (Arnaud et al. 2010) and $B_L = 1.37$ (Table 2). $N$ is a Gaussian log-normal correlated distribution at fixed mass, where

$$V_{\sigma} = \begin{pmatrix} \sigma_{\ln Y_{SZ}}^2 & r \sigma_{\ln Y_{SZ}} \sigma_{\ln Y_X} & s \sigma_{\ln Y_{SZ}} \sigma_{\ln L_{Xc}} \\ r \sigma_{\ln Y_{SZ}} \sigma_{\ln Y_X} & \sigma_{\ln Y_X}^2 & t \sigma_{\ln Y_X} \sigma_{\ln L_{Xc}} \\ s \sigma_{\ln Y_{SZ}} \sigma_{\ln L_{Xc}} & t \sigma_{\ln Y_X} \sigma_{\ln L_{Xc}} & \sigma_{\ln L_{Xc}}^2 \end{pmatrix}. \qquad (C.4)$$

For the covariance matrix $V_\sigma$, we assume $\sigma_{\ln Y_{SZ}} = 0.12$ (Kay et al. 2012; Le Brun et al. 2017); $\sigma_{\ln Y_X} = 0.16$ (Planelles et al. 2014; Le Brun et al. 2017; Truong et al. 2018); $r = 0.4$ (Farahi et al. 2019; Nagarajan et al. 2019); $s = 0.4$ (Farahi et al. 2019); and $t = 0.85$ (Farahi et al. 2019). As detailed in Sect. 5.2.3, with these assumptions, Eqn. 24 gives a first-order estimate of $\sigma_{\ln L_{Xc}} = 0.22$ in the absence of selection effects.

The Multi-Matched Filter extraction algorithm (Melin et al. 2006) was then applied, to obtain SZ detections at $S/N > 6$, corresponding to the threshold for the ESZ sample. We then matched the injected and recovered clusters to produce a mock ESZ catalogue, doing this twenty times, resulting in a total of 3188 detections. Measurement errors were estimated from the maps as described in Melin et al. (2006). Fits to the simulated data were performed as described in Sect. 4.2. Our baseline simulation in the following assumes $r = s = t = 0$ (i.e. zero covariance between quantities).

We first fitted the $Y_M - L_M$ relation, which assumes that there is zero intrinsic scatter in either of the observables with respect to the mass. The scatter in the $Y$-axis seen in Fig. C.1 is then entirely due to the observational uncertainties in the SZE measurements. The resulting slope of $0.759 \pm 0.001$ implies a change

<figure><img src="image_14.png" /><figcaption>Fig. C.1. Simulated $Y_{SZ}$ reextracted from PSZ2 maps, plotted as a function of $L_{Xc}$. The quantity $Q_M$ is the latent variable with respect to the mass $M$, obtained if there were no scatter in the relation. The Figure shows the effect of progressively adding scatter in $Y_{SZ}$, $Y_X$, and $L_{Xc}$ with $M$. The blue points include covariance of $t = 0.85$ between $Y_X$ and $L_{Xc}$.</figcaption></figure>

in the $L_{Xc} - M$ relation of $\Delta B_L = -0.01$. Adding intrinsic scatter of $\sigma_{\ln Y_{SZ}} = 0.12$ (Kay et al. 2012; Le Brun et al. 2017), the slope of the $Y_{SZ} - L_M$ relation is $0.759 \pm 0.001$, again implying a negligible change of $\Delta B_L = -0.01$ on the slope of the $L_{Xc} - M$ relation. These results imply that Malmquist bias in the $Y_{SZ}$ observable are negligible, and cannot account for the significantly steeper than self-similar slope we find for the $L_{Xc} - M$ relation.

We next studied the robustness of the recovery of the $L_{Xc} - M$ relation slope to intrinsic scatter in the $Y_X$ proxy, due to $Y_X$ being a scattered estimates of $Y_{SZ}$. An intrinsic scatter of $\sigma_{\ln Y_X} = 0.16$ (Planelles et al. 2014; Le Brun et al. 2017; Truong et al. 2018) pushes the slope of the $Y_X - L_M$ relation to $0.741 \pm 0.002$, in turn changing the $L_{Xc} - M$ relation slope by $\Delta B_L = -0.04$.

We then added an intrinsic scatter of $\sigma_{\ln L_{Xc}} = 0.22$, finding that this redresses the slope to $0.775 \pm 0.005$, implying a change of $\Delta B_L = +0.02$ in the slope of the $L_{Xc} - M$ relation. Finally, we added a covariance of $t = 0.85$ between $Y_X$ and $L_{Xc}$ (Farahi et al. 2019). This pushes the slope to a slightly steeper value of $0.785 \pm 0.003$, changing the $L_{Xc} - M$ relation slope by $\Delta B_L = +0.04$. The covariance significantly reduces the dispersion in the $Y_X - L_{Xc}$ relation, by a factor of two. This can clearly be seen in the difference in dispersion between the red and blue points in Fig. C.1.

In conclusion, we find that selection effects and intrinsic scatter have a negligible effect on the slope of the $L_{Xc} - M$ relation, and cannot account for the significantly steeper than self-similar value that we find in this work. The dispersion of the $L_{Xc} - M$ relation that we derive is significantly underestimated, most likely by a factor of $\sim 1.7$, due to the covariance between $L_{Xc}$ and the mass proxy, $Y_X$.