Towards a characterization of X-ray galaxy clusters for cosmology

TL;DR

This paper investigates galaxy cluster X-ray surface brightness profiles, comparing traditional models with a novel wavelet-based approach, to improve understanding of cluster morphologies and selection effects for cosmological studies.

Contribution

It introduces a wavelet decomposition method for modeling cluster SB profiles, highlighting biases in traditional $eta$-models and analyzing correlations with cluster properties.

Findings

Single $eta$-models cause a ~3% bias in outer SB measurements.

Fixing $eta$ to 2/3 increases residuals to over 40%.

No systematic core-type differences in outskirts profiles.

Abstract

In the framework of the hierarchical model the intra-cluster medium properties of galaxy clusters are tightly linked to structure formation, which makes X-ray surveys well suited for cosmological studies. To constrain cosmological parameters accurately by use of galaxy clusters in X-ray surveys, a better understanding of selection effects related to the detection method is needed. We aim at a better understanding of galaxy cluster morphologies to include corrections between the different core types and covariances with X-ray luminosities in selection functions. We stress the morphological deviations between a newly described surface brightness (SB) profile characterization and a commonly used single $\beta$-model. We investigate a novel approach to describe SB profiles, where the excess cool-core emission in the galaxy cluster centres is modelled using wavelet decomposition. Morphological parameters and the residuals are compared to classical single $\beta$-models. Using single $\beta$-models to describe the ensemble of overall SB profiles leads on average to a non-zero bias ($0.032 \pm 0.003$) in the outer part of the clusters, i.e. a $\sim 3\%$ systematic difference in the SB at large radii. In addition $\beta$-models show a general trend towards underestimating the flux in the outskirts for smaller core radii. Fixing the $\beta$ parameter to $2/3$ doubles the bias and increases the residuals from a single $\beta$-model up to more than $40\%$. Modelling the core region in the fitting procedure reduces the impact of these two effects significantly. We find a positive scaling between shape parameters and temperature, as well as a negative correlation ($\sim-0.4$) between extent and luminosity. Our non-parametric analysis of the self-similar scaled emission measure profiles indicates no systematic core-type differences of median profiles in the galaxy clusters outskirts.

Figures

Tables

Table 1. Table 1 : Slopes, m 𝑚 m , normalizations, b 𝑏 b , and intrinsic scatter values for the different y 𝑦 y -temperature scaling relations for single β 𝛽 \beta -model and core-modelled single β 𝛽 \beta -model parameters.

	Single $\beta$-model
				$\beta =2/3$
$𝒚$	$𝒎$	$𝒃$	${𝝀}_{𝒚}$	$𝒎$	$𝒃$	${𝝀}_{𝒚}$
$\beta$	${0.19}_{-0.06}^{+0.06}$	$-{0.82}_{-0.09}^{+0.09}$	${0.13}_{-0.01}^{+0.02}$	-	-	-
$r_c$	${1.04}_{-0.37}^{+0.37}$	${2.94}_{-0.60}^{+0.61}$	${0.88}_{-0.09}^{+0.10}$	${0.77}_{-0.14}^{+0.14}$	${3.73}_{-0.25}^{+0.25}$	${0.62}_{-0.05}^{+0.06}$
$E{\left(z\right)}^{-1}{L}_{0.1-2.4\mathrm{keV}}$	${1.90}_{-0.23}^{+0.23}$	${99.16}_{-0.38}^{+0.37}$	${0.49}_{-0.06}^{+0.07}$	${2.26}_{-0.14}^{+0.14}$	${98.58}_{-0.24}^{+0.24}$	${0.52}_{-0.05}^{+0.06}$
	Core-modelled single $\beta$-model
				$\beta =2/3$
$𝒚$	$𝒎$	$𝒃$	${𝝀}_{𝒚}$	$𝒎$	$𝒃$	${𝝀}_{𝒚}$
$\beta$	${0.16}_{-0.07}^{+0.07}$	$-{0.63}_{-0.11}^{+0.11}$	${0.16}_{-0.02}^{+0.02}$	-	-	-
$r_c$	${0.75}_{-0.22}^{+0.22}$	${4.10}_{-0.36}^{+0.36}$	${0.52}_{-0.05}^{+0.06}$	${0.50}_{-0.16}^{+0.16}$	${4.48}_{-0.26}^{+0.26}$	${0.40}_{-0.04}^{+0.05}$
$E{\left(z\right)}^{-1}{L}_{0.1-2.4\mathrm{keV}}$	${1.86}_{-0.23}^{+0.22}$	${99.22}_{-0.38}^{+0.37}$	${0.50}_{-0.06}^{+0.07}$	${1.90}_{-0.23}^{+0.22}$	${99.17}_{-0.37}^{+0.37}$	${0.49}_{-0.06}^{+0.07}$

Table 2. Table 2 : Correlation coefficients, λ x ​ y subscript 𝜆 𝑥 𝑦 \lambda_{xy} , between the different galaxy cluster parameters x 𝑥 x and y 𝑦 y for single β 𝛽 \beta -model and core-modelled single β 𝛽 \beta -model parameters.

		Single $\beta$-model		Core-modelled single $\beta$-model
			$\beta =2/3$		$\beta =2/3$
$𝒙$	$𝒚$	${𝝀}_{𝒙𝒚}$		${𝝀}_{𝒙𝒚}$
$\beta$	$r_c$	${0.60}_{-0.11}^{+0.09}$	-	${0.67}_{-0.10}^{+0.08}$	-
$\beta$	$E{\left(z\right)}^{-1}{L}_{0.1-2.4\mathrm{keV}}$	$-{0.02}_{-0.17}^{+0.17}$	-	${0.02}_{-0.16}^{+0.16}$	-
$r_c$	$E{\left(z\right)}^{-1}{L}_{0.1-2.4\mathrm{keV}}$	$-{0.43}_{-0.13}^{+0.16}$	$-{0.43}_{-0.13}^{+0.15}$	$-{0.32}_{-0.15}^{+0.16}$	$-{0.50}_{-0.13}^{+0.15}$

Table 3. Table 3 : List of parameters and their priors. We note that the prior of the scaling relation slope is assumed to be uniform in sin ⁡ ( Θ ) Θ \sin(\Theta) (VanderPlas 2016 ) , with Θ Θ \Theta being the angle between the best-fit line and the x 𝑥 x -axis. The term ”pos-normal” refers to a probability distribution that follows an ordinary normal distribution but is set to zero for negative parameter values, i.e. the parameter is restricted to be positive.

Parameter	Description	Prior
$\beta$	$\beta$-model slope	pos-normal(0.67,10)
$r_{\text{c}}$	$\beta$-model core radius in arcminutes	pos-normal(0,100)
$m$	Scaling relation slope	$\propto {\left(1+m^2\right)}^{-3/2}$
$b$	Scaling relation intercept	non-informative uniform
${\lambda }_y$	Intrinsic scatter in $y$	pos-normal(0,10)
${\lambda }_{xy}$	Correlation coefficient between $x$ and $y$	uniform[-1,1]

Table 4. Table 4 : Galaxy cluster parameters. Column (1) gives the cluster name. Columns (2) and (3) list the equatorial coordinates of the cluster center in decimal degrees based on the large scale wavelet image. Column (4) gives the offset to an iteratively determined two-dimensional ’center of mass’ using an aperture radius of 3   ′ times 3 arcminute 3\text{\,}\mathrm{\SIUnitSymbolArcminute} (Reiprich & Böhringer 2002 ) . Column (5) and (6) list the cluster redshift (Reiprich & Böhringer 2002 ) and core-excised temperature (Hudson et al. 2010 ) , respectively. Column (7) and (8) give the β 𝛽 \beta -model slope and core radius for a core-modelled fit. Column (9) lists the luminosity in the 0.1 – 2.4   keV times range 0.1 2.4 kiloelectronvolt 0.12.4\text{\,}\mathrm{keV} energy range. The cool-core classification according to Hudson et al. ( 2010 ) is given in column (10). Column (11) lists the characteristic radius where the density corresponds to 500 times the critical density. Galaxy clusters, whose r 500 subscript 𝑟 500 r_{500} value is marked with a † don’t have SZ mass estimates (Planck Collaboration et al. 2016 ) and the Schellenberger & Reiprich ( 2017 ) mass estimate is used to determine the wavelet small scales. These clusters are excluded from further analysis steps which comprise characteristic radii. Column (12) gives the measured large scale ellipticity. The physical to angular scale conversion at the cluster redshift is given in column (13).

Cluster	R.A.	Dec.	Offset	$𝒛$	$𝑻$	$𝜷$	${𝒓}_{\text{c}}$	$𝑳$	CCC	${𝒓}_{\mathrm{𝟓𝟎𝟎}}$	Ellipticity	Scale
name	J2000		$″$		$\mathrm{keV}$		$\mathrm{kpc}$	${10}^{44}\mathrm{erg}{s}^{-1}$		$\mathrm{Mpc}$		$\mathrm{kpc}/″$
A0085	10.4600	-9.3110	23.5	0.0556	${5.50}_{-0.10}^{+0.10}$	${0.731}_{-0.009}^{+0.010}$	${283.33}_{-6.43}^{+6.72}$	$5.034\pm 0.030$	SCC	1.18	0.102	1.08
A0119	14.0610	-1.2490	14.1	0.0440	${5.26}_{-0.26}^{+0.30}$	${0.681}_{-0.017}^{+0.017}$	${371.94}_{-13.22}^{+13.49}$	$1.725\pm 0.016$	NCC	1.05	0.177	0.87
A0133	15.6740	-21.8800	1.9	0.0569	${3.72}_{-0.09}^{+0.07}$	${0.634}_{-0.007}^{+0.008}$	${115.90}_{-3.59}^{+3.79}$	$1.514\pm 0.012$	SCC	1.00	0.240	1.10
NGC507	20.9080	33.2560	9.1	0.0165	${1.44}_{-0.10}^{+0.08}$	${0.420}_{-0.005}^{+0.005}$	${34.35}_{-2.10}^{+2.17}$	$0.127\pm 0.002$	SCC	${0.47}^{†}$	0.047	0.34
A0262	28.1960	36.1600	26.2	0.0161	${2.39}_{-0.04}^{+0.03}$	${0.483}_{-0.006}^{+0.007}$	${57.07}_{-2.45}^{+2.55}$	$0.534\pm 0.020$	SCC	0.73	0.303	0.33
A0400	44.4110	6.0050	44.8	0.0240	${2.23}_{-0.11}^{+0.09}$	${0.544}_{-0.010}^{+0.010}$	${126.97}_{-4.91}^{+5.06}$	$0.352\pm 0.004$	NCC	${0.54}^{†}$	0.371	0.48
A0399	44.4680	13.0440	9.3	0.0715	${6.11}_{-0.12}^{+0.12}$	${0.551}_{-0.013}^{+0.014}$	${120.71}_{-12.93}^{+13.25}$	$3.634\pm 0.196$	NCC	1.20	0.051	1.36
A0401	44.7380	13.5790	3.7	0.0748	${7.70}_{-0.19}^{+0.30}$	${0.527}_{-0.004}^{+0.004}$	${69.55}_{-3.68}^{+3.83}$	$6.452\pm 0.071$	NCC	1.30	0.346	1.42
A3112	49.4910	-44.2370	0.5	0.0750	${4.39}_{-0.11}^{+0.10}$	${0.646}_{-0.008}^{+0.009}$	${107.81}_{-4.03}^{+4.26}$	$3.832\pm 0.042$	SCC	0.99	0.246	1.42
FORNAX	54.6490	-35.3640	202.9	0.0046	${1.34}_{-0.00}^{+0.00}$	${0.613}_{-0.023}^{+0.025}$	${101.54}_{-4.02}^{+4.12}$	$0.042\pm 0.002$	SCC	${0.38}^{†}$	0.061	0.09
2A0335	54.6660	9.9710	9.5	0.0349	${3.34}_{-0.11}^{+0.09}$	${0.731}_{-0.011}^{+0.012}$	${141.70}_{-4.13}^{+4.23}$	$2.462\pm 0.020$	SCC	0.92	0.177	0.69
IIIZw54	55.3150	15.3970	46.4	0.0311	${2.44}_{-0.05}^{+0.04}$	${0.736}_{-0.078}^{+0.120}$	${166.98}_{-34.00}^{+44.75}$	$0.427\pm 0.033$	WCC	${0.56}^{†}$	0.104	0.62
A3158	55.7260	-53.6300	4.0	0.0590	${4.62}_{-0.06}^{+0.06}$	${0.697}_{-0.018}^{+0.020}$	${234.54}_{-11.41}^{+11.72}$	$2.899\pm 0.043$	NCC	1.12	0.207	1.14
A0478	63.3550	10.4660	1.3	0.0900	${6.67}_{-0.17}^{+0.16}$	${0.731}_{-0.006}^{+0.007}$	${186.43}_{-3.58}^{+3.65}$	$9.083\pm 0.054$	SCC	1.31	0.240	1.68
NGC1550	64.9130	2.4070	37.9	0.0123	${1.34}_{-0.00}^{+0.00}$	${0.615}_{-0.038}^{+0.047}$	${55.14}_{-9.43}^{+11.39}$	$0.155\pm 0.008$	SCC	${0.44}^{†}$	0.054	0.25
EXO0422	66.4670	-8.5620	16.4	0.0390	${2.81}_{-0.11}^{+0.11}$	${0.688}_{-0.026}^{+0.030}$	${70.53}_{-7.91}^{+8.55}$	$1.036\pm 0.064$	SCC	${0.69}^{†}$	0.205	0.77
A3266	67.8810	-61.4260	86.3	0.0594	${8.52}_{-0.32}^{+0.31}$	${0.933}_{-0.018}^{+0.018}$	${560.18}_{-12.26}^{+12.57}$	$4.483\pm 0.031$	WCC	1.30	0.397	1.15
A0496	68.4080	-13.2590	8.2	0.0328	${4.50}_{-0.05}^{+0.04}$	${0.651}_{-0.010}^{+0.010}$	${146.49}_{-4.67}^{+4.75}$	$1.972\pm 0.014$	SCC	0.97	0.217	0.65
A3376	90.5000	-39.9670	52.0	0.0455	${3.58}_{-0.09}^{+0.10}$	${0.675}_{-0.026}^{+0.029}$	${375.21}_{-21.94}^{+22.79}$	$1.118\pm 0.016$	NCC	0.93	0.495	0.89
A3391	96.5800	-53.6980	30.2	0.0531	${5.30}_{-0.32}^{+0.27}$	${0.762}_{-0.034}^{+0.038}$	${292.76}_{-19.82}^{+21.42}$	$1.379\pm 0.026$	NCC	0.98	0.398	1.03
A3395s	96.6970	-54.5460	11.2	0.0498	${4.47}_{-0.23}^{+0.23}$	${0.699}_{-0.086}^{+0.127}$	${286.63}_{-57.77}^{+77.07}$	$1.096\pm 0.042$	NCC	1.03	0.330	0.97
A0576	110.3530	55.7930	103.9	0.0381	${3.83}_{-0.09}^{+0.07}$	${0.716}_{-0.072}^{+0.107}$	${210.87}_{-39.80}^{+52.62}$	$0.962\pm 0.065$	WCC	0.90	0.255	0.76
A0754	137.3110	-9.6830	82.6	0.0528	${9.99}_{-0.38}^{+0.34}$	${0.767}_{-0.015}^{+0.016}$	${405.13}_{-11.60}^{+11.92}$	$2.052\pm 0.033$	NCC	1.32	0.455	1.03
HYDRA-A	139.5230	-12.0940	2.7	0.0538	${3.27}_{-0.08}^{+0.07}$	${0.719}_{-0.007}^{+0.007}$	${126.94}_{-2.63}^{+2.68}$	$3.049\pm 0.018$	SCC	${0.84}^{†}$	0.237	1.05
A1060	159.1680	-27.5280	41.9	0.0114	${3.02}_{-0.04}^{+0.04}$	${0.602}_{-0.007}^{+0.007}$	${98.62}_{-1.97}^{+1.92}$	$0.284\pm 0.009$	WCC	${0.68}^{†}$	0.062	0.23
A1367	176.2100	19.7030	68.4	0.0216	${3.38}_{-0.05}^{+0.05}$	${0.670}_{-0.021}^{+0.022}$	${346.02}_{-12.76}^{+13.68}$	$0.619\pm 0.005$	NCC	0.83	0.089	0.44
MKW4	181.1020	1.8990	38.1	0.0200	${2.01}_{-0.04}^{+0.04}$	${0.566}_{-0.011}^{+0.011}$	${79.67}_{-3.89}^{+4.05}$	$0.200\pm 0.003$	SCC	${0.51}^{†}$	0.275	0.41
ZwCl1215	184.4200	3.6590	7.2	0.0750	${5.74}_{-0.25}^{+0.28}$	${0.748}_{-0.017}^{+0.018}$	${248.82}_{-10.45}^{+10.78}$	$2.693\pm 0.035$	NCC	${1.01}^{†}$	0.268	1.42
NGC4636	190.7310	2.7600	272.0	0.0037	${0.90}_{-0.02}^{+0.02}$	${2.043}_{-0.682}^{+1.127}$	${149.01}_{-35.75}^{+45.65}$	$0.012\pm 0.001$	SCC	${0.34}^{†}$	0.176	0.08
A3526	192.2300	-41.3000	89.3	0.0103	${3.68}_{-0.02}^{+0.02}$	${0.551}_{-0.003}^{+0.003}$	${110.89}_{-1.29}^{+1.30}$	$0.637\pm 0.014$	SCC	0.77	0.280	0.21
A1644	194.2900	-17.3990	14.0	0.0474	${4.70}_{-0.08}^{+0.08}$	${0.688}_{-0.062}^{+0.084}$	${296.25}_{-54.79}^{+68.92}$	$1.993\pm 0.102$	SCC	1.07	0.050	0.93
A1650	194.6700	-1.7590	6.4	0.0845	${5.33}_{-0.06}^{+0.05}$	${0.642}_{-0.038}^{+0.048}$	${171.84}_{-26.11}^{+31.73}$	$3.754\pm 0.248$	WCC	1.13	0.160	1.59
A1651	194.8390	-4.1950	9.0	0.0860	${5.80}_{-0.24}^{+0.24}$	${0.696}_{-0.012}^{+0.012}$	${180.93}_{-6.44}^{+6.75}$	$4.109\pm 0.049$	WCC	1.18	0.152	1.61
COMA	194.9480	27.9300	30.6	0.0232	${8.26}_{-0.14}^{+0.15}$	${0.874}_{-0.004}^{+0.004}$	${449.21}_{-2.39}^{+2.37}$	$4.067\pm 0.057$	NCC	1.35	0.295	0.47
NGC5044	198.8640	-16.3860	38.8	0.0090	${1.22}_{-0.04}^{+0.03}$	${0.664}_{-0.013}^{+0.014}$	${65.97}_{-2.13}^{+2.19}$	$0.099\pm 0.000$	SCC	${0.41}^{†}$	0.063	0.18
A1736	201.7500	-27.1580	108.3	0.0461	${2.98}_{-0.11}^{+0.10}$	${0.591}_{-0.048}^{+0.062}$	${337.18}_{-56.53}^{+70.00}$	$1.657\pm 0.104$	NCC	0.99	0.065	0.91
A3558	202.0050	-31.5130	57.3	0.0480	${4.58}_{-0.13}^{+0.11}$	${0.656}_{-0.006}^{+0.007}$	${265.23}_{-4.54}^{+4.67}$	$3.401\pm 0.017$	WCC	1.17	0.289	0.94
A3562	203.3960	-31.6650	13.1	0.0499	${4.13}_{-0.14}^{+0.18}$	${0.585}_{-0.010}^{+0.011}$	${192.27}_{-7.30}^{+7.73}$	$1.603\pm 0.014$	WCC	0.94	0.230	0.98
A3571	206.8670	-32.8520	12.3	0.0397	${6.38}_{-0.11}^{+0.11}$	${0.692}_{-0.010}^{+0.010}$	${199.73}_{-5.15}^{+5.33}$	$4.181\pm 0.029$	WCC	1.16	0.336	0.79
A1795	207.2190	26.5980	14.7	0.0616	${5.57}_{-0.06}^{+0.06}$	${0.732}_{-0.003}^{+0.003}$	${166.95}_{-1.49}^{+1.52}$	$5.205\pm 0.016$	SCC	1.14	0.262	1.19
A3581	211.8790	-27.0180	22.3	0.0214	${1.97}_{-0.07}^{+0.07}$	${0.509}_{-0.017}^{+0.018}$	${25.47}_{-4.78}^{+5.36}$	$0.337\pm 0.011$	SCC	${0.63}^{†}$	0.210	0.43
MKW8	220.1690	3.4710	35.7	0.0270	${2.88}_{-0.11}^{+0.11}$	${0.543}_{-0.044}^{+0.061}$	${100.58}_{-26.45}^{+34.75}$	$0.405\pm 0.034$	NCC	0.72	0.151	0.54
A2029	227.7310	5.7440	7.8	0.0767	${7.48}_{-0.08}^{+0.08}$	${0.693}_{-0.006}^{+0.006}$	${160.89}_{-3.15}^{+3.15}$	$8.897\pm 0.053$	SCC	1.32	0.282	1.45
A2052	229.1880	7.0260	19.5	0.0348	${3.18}_{-0.02}^{+0.02}$	${0.791}_{-0.020}^{+0.022}$	${163.91}_{-6.39}^{+6.83}$	$1.259\pm 0.013$	SCC	${0.67}^{†}$	0.207	0.69
MKW3S	230.4660	7.7060	4.5	0.0450	${3.66}_{-0.08}^{+0.08}$	${0.707}_{-0.008}^{+0.009}$	${102.22}_{-2.47}^{+2.57}$	$1.473\pm 0.015$	SCC	0.86	0.277	0.88
A2065	230.6080	27.7170	18.8	0.0721	${4.98}_{-0.10}^{+0.18}$	${1.147}_{-0.209}^{+0.467}$	${510.85}_{-103.71}^{+182.02}$	$2.858\pm 0.174$	WCC	1.10	0.275	1.37
A2063	230.7790	8.6110	18.5	0.0354	${3.55}_{-0.05}^{+0.05}$	${0.570}_{-0.010}^{+0.011}$	${60.75}_{-4.99}^{+5.21}$	$1.168\pm 0.015$	WCC	0.86	0.160	0.70
A2142	239.6000	27.2220	70.9	0.0899	${7.60}_{-0.67}^{+0.88}$	${0.923}_{-0.013}^{+0.013}$	${556.15}_{-10.08}^{+10.25}$	$10.960\pm 0.099$	WCC	1.41	0.388	1.68
A2147	240.5640	15.9550	14.4	0.0351	${3.81}_{-0.11}^{+0.10}$	${0.359}_{-0.017}^{+0.021}$	${85.48}_{-19.42}^{+24.73}$	$1.500\pm 0.048$	NCC	1.06	0.088	0.70
A2163	243.9410	-6.1470	13.6	0.2010	${14.17}_{-0.71}^{+0.71}$	${0.738}_{-0.014}^{+0.015}$	${368.75}_{-13.65}^{+14.01}$	$17.179\pm 0.258$	NCC	1.67	0.151	3.31
A2199	247.1580	39.5480	2.0	0.0302	${4.07}_{-0.06}^{+0.06}$	${0.751}_{-0.005}^{+0.005}$	${182.12}_{-1.93}^{+1.96}$	$2.140\pm 0.039$	SCC	0.99	0.201	0.60
A2204	248.1960	5.5760	10.2	0.1523	${8.06}_{-0.53}^{+0.63}$	${0.707}_{-0.018}^{+0.020}$	${173.91}_{-10.33}^{+10.89}$	$13.715\pm 0.219$	SCC	1.33	0.046	2.65
A2244	255.6740	34.0590	4.9	0.0970	${5.31}_{-0.10}^{+0.09}$	${0.642}_{-0.014}^{+0.015}$	${126.79}_{-7.43}^{+7.91}$	$4.345\pm 0.091$	WCC	1.12	0.145	1.80
A2256	255.9750	78.6400	31.0	0.0601	${6.91}_{-0.55}^{+0.57}$	${1.041}_{-0.013}^{+0.014}$	${547.10}_{-7.35}^{+7.58}$	$4.793\pm 0.067$	NCC	1.27	0.341	1.16
A2255	258.2080	64.0650	26.5	0.0800	${5.33}_{-0.18}^{+0.17}$	${0.851}_{-0.027}^{+0.029}$	${502.79}_{-20.19}^{+21.65}$	$2.829\pm 0.034$	NCC	1.21	0.209	1.51
A3667	303.1450	-56.8430	17.3	0.0560	${5.84}_{-0.04}^{+0.04}$	${0.567}_{-0.006}^{+0.007}$	${260.96}_{-6.82}^{+7.05}$	$4.949\pm 0.035$	WCC	1.33	0.341	1.09
S1101	348.4940	-42.7270	1.3	0.0580	${2.50}_{-0.11}^{+0.11}$	${0.783}_{-0.011}^{+0.011}$	${115.72}_{-3.15}^{+3.27}$	$1.850\pm 0.017$	SCC	0.79	0.178	1.12
A2589	350.9880	16.7720	10.7	0.0416	${3.65}_{-0.04}^{+0.04}$	${0.671}_{-0.016}^{+0.017}$	${148.38}_{-6.86}^{+7.35}$	$0.989\pm 0.013$	WCC	0.84	0.285	0.82
A2597	351.3320	-12.1250	1.9	0.0852	${3.79}_{-0.06}^{+0.06}$	${0.702}_{-0.013}^{+0.014}$	${106.77}_{-5.02}^{+5.42}$	$3.535\pm 0.042$	SCC	0.93	0.192	1.60
A2634	354.6090	27.0130	61.9	0.0312	${3.04}_{-0.10}^{+0.10}$	${0.680}_{-0.027}^{+0.030}$	${297.03}_{-16.72}^{+17.64}$	$0.518\pm 0.008$	WCC	0.80	0.554	0.62
A2657	356.2290	9.1940	16.4	0.0404	${3.33}_{-0.10}^{+0.11}$	${0.555}_{-0.006}^{+0.006}$	${96.03}_{-3.00}^{+3.05}$	$0.911\pm 0.008$	WCC	0.79	0.125	0.80
A4038	356.9320	-28.1450	13.3	0.0283	${3.00}_{-0.04}^{+0.03}$	${0.620}_{-0.013}^{+0.014}$	${98.38}_{-4.88}^{+5.08}$	$1.005\pm 0.013$	WCC	0.80	0.313	0.57
A4059	359.2520	-34.7580	7.1	0.0460	${3.94}_{-0.03}^{+0.03}$	${0.705}_{-0.018}^{+0.020}$	${168.53}_{-8.62}^{+9.33}$	$1.477\pm 0.019$	SCC	0.94	0.255	0.90

Equations

$$T_{\mathrm{4.1.1}}=0.875\cdot T_{\mathrm{3.2.1}}+0.251,$$

$$r_{500}={\left(\frac{3M_{500}}{4\pi \cdot 500{\rho }_{\text{crit},z}}\right)}^{1/3}.$$

$$s_{\textnormal{X}}(R)=\sum_{j=1}^{N}s_{0,j}\left[1+\left(\frac{R}{r_{\textnormal{c},j}}\right)^{2}\right]^{-3\beta_{j}+$0.5$}.$$

$$\ln \mathcal{L}\propto \sum _i{O}_i\ln \left(M_i\right)-M_i,$$

$$\ln \mathcal{L}\propto \sum _i{O}_i\ln \left(M_i+B_{\text{wv},i}+B_{\text{p},i}\right)-\left(M_i+B_{\text{wv},i}+B_{\text{p},i}\right).$$

$$s_{i}(R_{i})=s_{0}\left[1+\left(\frac{R_{i}}{r_{\textnormal{c}}}\right)^{2}\right]^{-3\beta+$0.5$}+b_{\textnormal{c}}.$$

$$M_i=s_i(R_i)\cdot {\alpha }_i\cdot {ϵ}_{\text{mean},i}.$$

$$\frac{10^{-14}}{4\pi [D_{\text{A}}(1+z){]}^2}\int n_{\text{e}}{n}_{\text{H}}\text{d}V.$$

$$EM(R)=\int n_{\text{e}}{n}_{\text{H}}\text{d}l,$$

$$T\propto {\left(E(z)M_{500}\right)}^{2/3}$$

$$\Delta_{\textnormal{z}}^{3/2}(1+z)^{9/2}\left(\frac{E(z)M_{500}}{$2\text{\cdot}{10}^{15}\text{\,}\mathrm{\textnormal{M}_{\odot}}$}\right)^{1/3},$$

$${\Delta }_{\text{z}}={\Delta }_{\text{c}}{\Omega }_{\text{m}}/\left(18{\pi }^2{\Omega }_{\text{z}}\right).$$

$${\Delta }_{\text{c}}=18{\pi }^2+82w-39w^2,$$

$$\log y/n_y=m\cdot \log x/n_x+b.$$

$$\mathcal{L}=\prod _{n=1}^N\frac{1}{2\pi \sqrt{\det ({\Sigma }_n+\Lambda )}}\exp \left(-\frac{1}{2}{\tilde{r}}_n^{\mathrm{T}}({\Sigma }_n+\Lambda {)}^{-1}{\tilde{r}}_n\right),$$

$${\tilde{r}}_n$$

$${\Sigma }_n$$

$$\left(\begin{array}[]{c}\hat{x}_{n}\\ \hat{y}_{n}\end{array}\right)\sim\mathcal{N}\left(\left(\begin{array}[]{c}\tilde{x}_{n}\\ m\,\tilde{x}_{n}-b\end{array}\right),\Lambda\right).$$

$$\mathrm{cov}(x,y):=\frac{1}{N-1}\sum _{i=0}^N(x_i-{\bar{x}}_i)(y_i-{\bar{y}}_i).$$

$${\lambda }_{xy}:=\frac{\mathrm{cov}(x,y)}{{\lambda }_x{\lambda }_y}.$$

$${\Lambda }_{l,m}:=\mathrm{cov}(X_l,X_m).$$

$$\Lambda=\pmatrix{}{cov}(x,x)&\mathrm{cov}(y,x)\\ \mathrm{cov}(x,y)\mathrm{cov}(y,y)=\pmatrix{\lambda}_{x}^{2}&\lambda_{xy}\lambda_{x}\lambda_{y}\\ \lambda_{xy}\lambda_{x}\lambda_{y}\lambda_{y}^{2}.$$

$$e=1-e_2/e_1$$

$$M_{500,\text{GS17}}^{\ast }\left(<r_{500,\text{Planck}}\right)=M_{500,\text{GS17}}\frac{Y(\frac{r_{500,\text{Planck}}}{r_{500,\text{GS17}}}{c}_{500,\text{GS17}})}{Y(c_{500,\text{GS17}})},$$

Full text

*Aims. *We aim at a better understanding of the morphology of galaxy clusters to include corrections between the different core types and covariances with X-ray luminosities in selection functions. In particular we stress the morphological deviations between a newly described surface brightness profile characterization and a commonly used single $\beta$-model.

*Methods. *We investigate a novel approach to describe surface brightness profiles, where the excess cool-core emission in the centres of the galaxy clusters is modelled using wavelet decomposition. Morphological parameters and the residuals are compared to classical single $\beta$-models, fitted to the overall surface brightness profiles.

*Results. *Using single $\beta$-models to describe the ensemble of overall surface brightness profiles leads on average to a non-zero bias ($0.032\pm 0.003$) in the outer part of the clusters, i.e. a $\sim 3\%$ systematic difference in the surface brightness at large radii. In addition $\beta$-models show a general trend towards underestimating the flux in the outskirts for smaller core radii. Fixing the $\beta$ parameter to $2/3$ doubles the bias and increases the residuals from a single $\beta$-model up to more than $40\%$. Modelling the core region in the fitting procedure reduces the impact of these two effects significantly. In addition, we find a positive scaling between shape parameters and temperature, as well as a negative correlation ($\sim-0.4$) between extent and luminosity.

*Conclusions. *We demonstrate the caveats in modelling galaxy clusters with single $\beta$-models and recommend the usage of them with caution, especially when not taking the systematics into account that arise using them. Our non-parametric analysis of the self-similar scaled emission measure profiles indicates no systematic core-type differences of median profiles in the galaxy clusters outskirts.

††offprints: F. Käfer,

1 Introduction

Clusters of galaxies are formed from the collapse of initial density fluctuations in the early Universe and grow hierarchically to the densest regions of the large-scale structure. This makes them the most massive (M_{\text{tot}}\sim$${10}^{14}{10}^{15}\text{\,}\mathrm{\textnormal{M}_{\odot}}) gravitationally bound structures in our universe and their virialisation timescales are less than the Hubble time. The gas between the galaxies, the intra-cluster medium (ICM), has been heated to temperatures111X-ray gas temperatures are often expressed in kiloelectronvolt ($\text{\,}\mathrm{keV}$). Using the Boltzmann constant $k_{\textnormal{B}}$, 1\text{\,}\mathrm{keV}$/k_{\textnormal{B}}\approx$1.16\text{\cdot}{10}^{7}\text{\,}\mathrm{K}. of several ${10}^{7}\text{\,}\mathrm{K}$ by gravitational collapse. The primary emission mechanism of this hot, fully-ionised thermal plasma is thermal bremsstrahlung and line emission of heavy elements, e.g. iron. The majority ($\sim$85 %) of the baryonic component is in the form of the hot ICM. Therefore the most massive visible component can be traced by X-ray emission, which makes X-ray astronomy a great and important tool to study galaxy clusters. However, flux-limited galaxy cluster samples compiled from X-ray surveys suffer from selection effects like Malmquist bias, i.e. the preferential detection of intrinsically brighter sources (a more detailed discussion of different selection effect biases is compiled in e.g. Hudson et al. 2010; Giodini et al. 2013). Another form of selection effect arises from the different core-types of galaxy clusters. In the central regions of galaxy clusters, gas is able to cool more efficiently compared to the outskirts. Several diagnostics were proposed to identify and categorize galaxy clusters according to their different core-types, e.g. a central temperature decrease (Sanderson et al., 2006), mass-deposition rates (Chen et al., 2007), cuspiness (Vikhlinin et al., 2007) or surface brightness concentration (Santos et al., 2008). Galaxy clusters exhibiting cool-cores show centrally-peaked surface brightness profiles, whereas non-cool-core clusters have flat profiles. In surveys differently shaped profiles are detected with different efficiencies. Even for the same brightness, cool-core clusters may be more easily detected, since their surface brightness profiles are more peaked and thus the central emission sticks out more above the background. The preferential detection of cool-core objects close to the detection threshold of flux-limited samples leads to the so-called cool-core bias (e.g. Eckert et al., 2011; Rossetti et al., 2017). It is crucial to take such selection effects into account in cosmological studies to obtain unbiased results. One possibility to quantify these biases is running the source detection chain on well-defined simulations. The quantification of the completeness, i.e. the fraction of detected clusters as function of mass and redshift, requires an accurate galaxy cluster model as input for such simulations.

Outside the core regions, scaled radial profiles (e.g. temperature, pressure or entropy profiles) of galaxy clusters show a so-called ”self-similar” behaviour (e.g. Zhang et al., 2007; Ghirardini et al., 2018). It is believed that this is the result of a similar formation process of galaxy clusters, namely that tiny density perturbations in the early universe are amplified by gravitational instabilities and grow hierarchically, yielding the large scale structure observed today. Galaxy clusters are then believed to correspond to the densest regions of the large scale structure. This formation history motivated the theoretical consideration of the self-similar model (e.g. Kaiser, 1986), where all galaxy clusters share the same average density and evolve with redshift and mass according to prescriptions given by spherical gravitational collapse. Therefore galaxy cluster observables such as X-ray luminosity, spectral temperature or gas mass are correlated to the total cluster mass. Assuming gravity is the dominant process, the self-similar model predicts simple power-law relations between those cluster observables and the total mass, so-called scaling relations (e.g. Maughan, 2007; Pratt et al., 2009; Mantz et al., 2010, 2016; Maughan et al., 2012).

In this work, we aim towards a proper characterization of galaxy cluster shapes using different surface brightness parameterisations. We investigate scaling relations between surface brightness parameters and temperature by use of the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS), a statistically complete, X-ray-selected and X-ray flux-limited sample of $64$ galaxy clusters compiled from the ROSAT all-sky survey (RASS, Voges et al. 1999). In addition we study the covariances between shape and other galaxy cluster parameters. The impact of the different core-types on the obtained scaling relations and covariances are quantified. The goal of this paper is to improve our understanding of galaxy cluster shapes. This serves as a basis for simulations quantifying selection effects, inter alia for the future X-ray all-sky survey performed by the extended Roentgen Survey with an Imaging Telescope Array (eROSITA, Merloni et al. 2012; Predehl et al. 2018). In addition, the obtained covariance matrices can be implemented in current cosmological studies using e.g. the COnstrain Dark Energy with X-ray galaxy clusters (CODEX) sample.

Throughout this paper a flat $\Lambda$CDM cosmology is assumed. The matter density, vacuum energy density and Hubble constant are assumed to be $\Omega_{\textnormal{m}}=$0.3, $\Omega_{\Lambda}=$0.7 and $H_{0}=$ $70\text{\,}\mathrm{km}\text{\,}{\mathrm{s}}^{-1}\text{\,}{\mathrm{Mpc}}^{-1}$ with $h_{70}\coloneqq H_{0}/$70\text{,}\mathrm{km}\text{,}{\mathrm{s}}^{-1}\text{,}{\mathrm{Mpc}}^{-1}$=$1$$, respectively. ’ln’ refers to the natural logarithm and ’log’ is the logarithm to base 10. All errors are $1\sigma$ unless otherwise stated.

2 Data

This section describes the underlying galaxy cluster sample used for this study, as well as a description of the data analysis.

2.1 The sample

The HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS, Reiprich & Böhringer, 2002; Hudson et al., 2010) comprises $64$ galaxy clusters, constructed from highly complete cluster catalogs based on the ROSAT all-sky survey (RASS). The final flux limit of f_{\text{X}}\text{(0.12.4\text{,}\mathrm{keV})}=\text{20\text{\cdot}{10}^{-12}\text{,}\mathrm{erg}\text{,}{\mathrm{s}}^{-1}\text{,}{\mathrm{cm}}^{-2}} defines the X-ray-selected and X-ray flux-limited sample of the brightest galaxy clusters away from the Galactic plane. Although statically complete, HIFLUGCS is not necessarily representative or unbiased with respect to the cluster morphology (Hudson et al., 2010; Mittal et al., 2011). Eckert et al. (2011) calculated a significant bias in the selection of X-ray clusters of about 29 percent in favour of centrally peaked, cool-core objects compared to non-cool-core (NCC) clusters. We minimize these kind of selection effects by restricting our study to objects above a temperature of $3\text{\,}\mathrm{keV}$. This temperature threshold excludes the low mass galaxy groups, which are closer to the HIFLUGCS flux threshold and therefore have a high cool-core fraction as discussed in Sect. 1. Originally, the cluster RX J1504.1-0248 was not included in HIFLUGCS. This object is a strong cool-core (SCC) cluster that appears only marginally extended in the RASS, i.e. its extent is comparable to the ROSAT survey point-spread function (PSF). To avoid biasing our results because of the small extent of this system compared to the ROSAT PSF, this cluster is excluded in our full analysis. In total, we consider $49$ galaxy clusters above our selected temperature threshold of $3\text{\,}\mathrm{keV}$.

2.2 Data analysis

Pointed ROSAT observations were used whenever available. The pointed data is reduced with the ROSAT Extended Source Analysis Software package (Snowden et al., 1994) as described in Eckert et al. (2012). Otherwise RASS data from the public archive are used. There are $4$ objects without pointed observations above our representative temperature threshold of $3\text{\,}\mathrm{keV}$ and excluding those from the analysis does not change our results significantly. Therefore, we do not expect our results to be affected by the use of heterogeneous data. All images for count rate measurements are restricted to the ROSAT hard energy band (channels $42201$ $\approx$ $0.42.0\text{\,}\mathrm{keV}$) due to higher background levels in the soft band. Luminosities are taken from Reiprich & Böhringer 2002, i.e. are based on ROSAT data and are not corrected for cooling flows. The central temperature drop of the ICM in cool-core clusters biases the estimation of the cluster virial temperature, i.e. the temperature of the hot gas which is in hydrostatic equilibrium with the potential well of the cluster. This bias is a source of scatter in scaling relations related to temperature and can be minimized by excluding the central region for the temperature fitting. Since we are interested in how the galaxy cluster shape parameters scale with temperature, we adapt core-excised HIFLUGCS temperatures measured by Hudson et al. 2010 using Chandra’s Advanced CCD Imaging Spectrometer (ACIS) data. We re-scale all temperatures greater than $2\text{\,}\mathrm{keV}$ due to the Chandra calibration package updates according to the Mittal et al. (2011) best-fit relation

[TABLE]

which links temperature measurements (in $\text{\,}\mathrm{keV}$) between the Calibration Database (CALDB) 3.2.1 and CALDB 4.1.1. The reader is referred to the aforementioned papers in this subsection for a detailed description of individual data analysis steps.

2.3 Masses

We adapt $M_{500}$ values of the ”Union catalogue” (Planck Collaboration et al., 2016) calculated by Planck Sunyaev-Zel’dovich (SZ) observations, which contains detections with a minimum signal-to-noise of $4.5$. We note that SZ mass uncertainties are small due to being purely statistical and are not propagated when rescaling radii. The advantage of using masses calculated by the use of the $Y_{\textnormal{SZ}}-M$ relation in this study is that they are statistically not covariant with X-ray parameters and less affected by galaxy cluster core states (Lin et al., 2015). Above our selected temperature threshold of $3\text{\,}\mathrm{keV}$ there are 4 galaxy clusters without counterpart in the Planck catalogue (Hydra-A, A1060, ZwCl1215, A2052). These clusters are rejected for studies that require characteristic masses or radii. Assuming spherical symmetry, the galaxy cluster masses can be transformed into $r_{500}$ values according to

[TABLE]

3 Analysis

3.1 Surface brightness profiles

Using King’s (1962) analytical approximation of an isothermal sphere, measured X-ray surface brightness profiles of galaxy clusters are well described222Note that the assumption of a single $\beta$-models is that the hot gas and the galaxies are in hydrostatic equilibrium and isothermal. by a $\beta$-model (Cavaliere & Fusco-Femiano, 1976)

[TABLE]

For each component $j$, $s_{0,j}$ is the central surface brightness, i.e. at projected radial distance $R=0$, $r_{\textnormal{c},j}$ is the core radius of the gas distribution and the slope $\beta_{j}$ is motivated by the ratio of the specific energy in galaxies to the specific energy in the hot gas. For galaxy clusters exhibiting a central excess emission due to the presence of cool cores, a double ($N=2$) $\beta$-model can improve the agreement between model and data as one component accounts for the central excess emission while the other accounts for the overall cluster emission. However, the two components are highly degenerate and except for very nearby galaxy clusters, the ROSAT point-spread function is insufficient to resolve the core regions since the apparent size of the objects is smaller. Therefore, a single ($N=1$) beta model is used to describe the galaxy cluster emission and the central excess emission is included in the background model. Simulations (Navarro et al., 1995; Bartelmann & Steinmetz, 1996) indicate that the measured $\beta$ values are biased systematically low if the range of radii used for fitting is less than the virial radius of the cluster. The advantage of using ROSAT PSPC data to determine the surface brightness profiles is in the large field-of-view and the low background, allowing to trace the galaxy cluster emission to relatively large radii.

3.1.1 Wavelet decomposition

We use a wavelet decomposition technique as described in Vikhlinin et al. (1998). The technique is implemented as wvdecomp task of the publicly available ZHTOOLS333Please contact A. Vikhlinin for the latest version of ZHTOOLS (avikhlinin@cfa.harvard.edu). package. The basic idea is to convolve the input image with a kernel which allows isolation of structures of given angular size. Particular angular sizes are isolated by varying the scale of the kernel. The wavelet kernel on scale $i$ used in wvdecomp is approximately the difference of two Gaussians, isolating structures in the convolved image of a characteristic scale of $\approx 2^{i-1}$. The input image is convolved with a series of kernels with varying scales, starting with the smallest scale. In each step, significantly detected features of the particular scale are subtracted from the input image before going to the next scale. This allows, among other things, to decompose structures of different sizes into their components, e.g. in the case of point-like sources in the vicinity of an extended object. Wavelet kernels have the advantage of a simple linear back transformation, i.e. the original image is the sum of the different scales. We define a scale around $0.2\text{\,}\mathrm{\textnormal{r}_{500}}$ up to which all emission from smaller scales is classified as contamination and is included in our background modelling for the core-modelled single $\beta$-model approach. The galaxy clusters $0.2\text{\,}\mathrm{\textnormal{r}_{500}}$ wavelet scales are around $35$ ($23$) for pointed (survey) observations. This corresponds to $416\text{\,}\mathrm{pixel}$ ($24\text{\,}\mathrm{pixel}$), with a pixel size of $15\text{\,}\arcsec$ ($45\text{\,}\arcsec$). The detection threshold of a wavelet kernel convolved image is the level above which all maxima are statistically significant. Vikhlinin et al. (1998) performed Monte-Carlo simulations of flat Poisson background to define detection thresholds such that one expects on average $1/3$ false detections per scale in a $512\text{\times}512\text{\,}\mathrm{pixel}$ image. We adapt a slightly more stringent threshold of $5\sigma$.

3.1.2 Likelihood function

Under the assumption that the observed counts are Poisson distributed, the maximum-likelihood estimation statistic to estimate the surface brightness profile parameters is chosen to be the Poisson likelihood. The so-called Cash statistic (Cash, 1979) is derived by taking the logarithm of the Poisson likelihood function and neglecting the constant factorial term of the observed counts

[TABLE]

where $M_{i}$ and $O_{i}$ are the model and observed counts in bin $i$, respectively. The model counts of the background sources using wavelet decomposition, $B_{\textnormal{wv},i}$, are not Poissonian. We assume this background component without error, i.e. just the total amount of counts show dispersion. Thus, we can add this background component to the model counts (Greiner et al., 2016). In the same way, we add an additional particle background component, $B_{\textnormal{p},i}$, to equation (4) for pointed observations, i.e. the likelihood function becomes

[TABLE]

A single $\beta$- plus constant background model is used to describe the surface brightness of each cluster (see equation (3), using $N=1$ and dropping the index $j$)

[TABLE]

The projected radii $R_{i}$ are placed at the center of the bins. By use of the exposure map, we calculate the proper area, $\alpha_{i}$, and the vignetting corrected mean exposure time, $\epsilon_{\textnormal{mean},i}$. The model counts in each bin are then calculated by multiplying equation (6) with the corresponding area and exposure time

[TABLE]

3.2 Point-spread function

The ability of an X-ray telescope to focus photons, i.e. its response to a point source, is characterized by its point-spread function. More peaked cool-core objects are affected more by PSF effects compared to non-cool-core objects. The ROSAT PSF depends amongst others on photon energy, off-axis angle and observation mode. A detailed description of the ROSAT PSF functions is presented in Boese (2000). We use the Python package pyproffit444https://github.com/domeckert/pyproffit to calculate PSF mixing matrices based on equations (7) and (30) of Boese (2000) for pointed and survey observations, respectively. These matrices are folded in our surface brightness profile fitting method to obtain PSF unconvolved parameters.

3.3 Emission measure profiles

This subsection describes our approach to obtain background subtracted self-similar scaled emission measure profiles. First, the outer significance radius and background level of each galaxy cluster are iteratively determined using the growth curve analysis method (Böhringer et al., 2000; Reiprich & Böhringer, 2002). The outer significance radius determines the maximum radius out to which galaxy cluster emission is detected and thus to which radius each profile is extracted. Background-subtracted and logarithmically binned surface-brightness profiles are converted into emission measure profiles using the normalisation of a partially absorbed Astrophysical Plasma Emission Code (APEC) model

[TABLE]

The total weighted hydrogen column density (calculated with the method of Willingale et al. (2013))555http://www.swift.ac.uk/analysis/nhtot/index.php is used to describe the absorption by the atomic and molecular Galactic column density of hydrogen. Metallicities are fixed to $0.35\text{\,}\mathrm{\textnormal{Z}_{\odot}}$ and the abundance table compiled by Anders & Grevesse (1989) is used. The emission measure along the line-of-sight,

[TABLE]

is self-similar scaled according to Arnaud et al. (2002) and

[TABLE]

by

[TABLE]

where $\Delta_{\textnormal{z}}$ is calculated using the density contrast, $\Delta_{\textnormal{c}}$, and matter density parameter at redshift $z$, $\Omega_{\textnormal{z}}=\Omega_{\textnormal{m}}(1+z)^{3}/E(z)^{2}$, according to

[TABLE]

Under the assumption that the cluster has just virialized, Bryan & Norman (1998) derived an analytical approximation of $\Delta_{\textnormal{c}}$ for a flat universe from the solution to the collapse of a spherical top-hot perturbation

[TABLE]

with $w=\Omega_{\textnormal{z}}-1$.

3.4 Scaling relations

In this subsection we describe the basic principle of our linear regression routine to obtain scaling relations. A set of two variates, $x/y$, is fitted by a power-law relation according to

[TABLE]

The pivot elements, $n_{x/y}$, are set to the median along a given axis, such that the results of the slope and normalisation are approximately uncorrelated.

3.4.1 Likelihood function

Linear regression of the scaling relations is performed using a Markov chain Monte Carlo (MCMC) posterior sampling technique. We adapt an $N$ dimensional Gaussian likelihood function

[TABLE]

extended compared to Kelly (2007) to account for intrinsic scatter correlation. The intrinsic scatter tensor, $\Lambda$, is described in more detail in Sect. 3.4.2. The uncertainty tensors $\Sigma_{n}$ account for measurement errors in the independent and the dependent variables and $\tilde{r}_{n}$ denote the residual vectors. For illustration purposes, this is how these two objects would look like in a bivariate example:

[TABLE]

In this study, the correlation between different measurement errors in the uncertainty tensor is set to zero. The ”true” coordinate $\tilde{x}_{n}$ is normal-scattered according to the intrinsic scatter tensor via

[TABLE]

We integrate out, i.e. marginalize over, $\hat{x}_{n}$ and $\hat{y}_{n}$. The scatter along the independent axis, $\lambda_{x}$, of the intrinsic scatter tensor is fixed to avoid degeneracies. This means that for this study the intrinsic scatter in temperature is fixed to $20\text{\,}\mathrm{\char 37\relax}$, i.e. $\lambda_{\mathrm{T}}=0.11$ (Kravtsov et al., 2006). The correlation between the intrinsic scatter values of the two variates $x$ and $y$, $\lambda_{xy}$, is of particular interest for this study and will be described in more detail in Sect. 3.4.2. We use the emcee algorithm and implementation (Foreman-Mackey et al., 2013) for optimization. A chain is considered as converged when the integrated autocorrelation time is greater than one-hundredth of the chain length.

3.4.2 Covariance

The linear relationship and thus the joint variability between two or more sets of random variables can be quantified by the covariance between those variates. In the simple case of two variables $x$ and $y$, each with a sample size of $N$ and expected values $\bar{x}$ and $\bar{y}$, the covariance is given by

[TABLE]

The degree of correlation can be calculated by normalizing the covariance to the maximum possible dispersion of the single standard deviations $\lambda_{x}$ and $\lambda_{y}$, the so-called Pearson correlation coefficient:

[TABLE]

The Pearson correlation coefficient can take values between $-1$ and $+1$, where [math] means no linear correlation and $+1$ ($-1$) means total positive (negative) linear correlation. In the general case of $n$ sets of variables $\{X_{1}\},\dots,\{X_{n}\}$, the covariances can be displayed in a matrix, where the first-order covariance matrix is defined by

[TABLE]

In the previous example of 2 variables $x$ and $y$, the covariance matrix reads

[TABLE]

The latter equality makes use of equation (23), which implies that the covariance of a variate with itself, i.e. $\mathrm{cov}(x,x)$, reduces to the variance of $x$ or the square of the standard deviation of $x$. The off-axis elements are rewritten by solving equation (24) for the covariance and using the symmetry $\lambda_{xy}=\lambda_{yx}$.

Calculating the Pearson correlation coefficient between the ranked variables is a non-parametric measure of a monotonic relationship between the variables and is called the Spearman rank-order correlation coefficient.

3.4.3 Selection effects

As already discussed in Sect. 1, centrally peaked galaxy clusters are more likely to enter an X-ray selected sample due to their enhanced central emission. Mittal et al. (2011) investigated this effect by applying the HIFLUGCS flux limit to Monte Carlo simulated samples. Assuming HIFLUGCS being complete, one can vary the input fractions of different core-types in the simulations to match the observed ones. The intrinsic scatter increases the normalization of the luminosity-temperature relation because up-scattered clusters have a higher chance of lying above the flux threshold. In this study, we are not trying to determine the true luminosity-temperature relation but are interested in the residuals of the sample with respect to the mean to study the intrinsic scatter covariances. Therefore we are neglecting Malmquist bias in the parameter optimization, although it is present in HIFLUGCS. To investigate the effect of Malmquist bias on the best-fit shape-temperature relation parameters and the intrinsic scatter correlation coefficients, we artificially decrease the luminosity-temperature relation normalization and find that the differences are insignificant.

3.5 Cool-core classification

Hudson et al. (2010) used HIFLUGCS to compare $16$ different techniques to differentiate cool-core and non-cool-core clusters. The central cooling time, $t_{\textnormal{cool}}$, was found to be suited best and used to classify clusters into three categories. Clusters with central cooling times shorter than $1\text{\,}\mathrm{Gyr}$ are classified as strong-cool-core (SCC) clusters. They usually show characteristic temperature drops towards the center and low central entropies. Clusters exhibiting high central entropies and cooling times greater than $7.7\text{\,}\mathrm{Gyr}$ are classified as non-cool-core (NCC) clusters. In intermediate class with cooling times in between those of SCC and NCC clusters are classified as weak-cool-core (WCC) clusters. We adapt the Hudson et al. (2010) classification scheme and categorization of HIFLUGCS clusters for this study. There are $45$ galaxy clusters above our selected temperature threshold of $3\text{\,}\mathrm{keV}$ with mass estimates in the Planck ”Union catalogue”. The amount of each core-type category is $15$, $16$ and $14$ for SCC, WCC and NCC, respectively. For $1$ of the SCC and $3$ of the WCC objects, no ROSAT pointed observations are available and RASS data is used.

4 Emission measure profiles

This section describes a model-independent way to test the surface brightness (SB) behaviour between different core-type populations by comparing their background subtracted self-similar scaled emission measure (EM) profiles. Fig. 1 shows the SB and self-similar scaled EM profiles. Outside the cluster core the galaxy clusters show self-similar behaviour, as already discussed in several previous studies (e.g. Vikhlinin et al. 1999). The self-similar scaled EM profiles show a smaller intrinsic scatter compared to the surface brightness profiles in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ range. It is reduced by $28$% ($33$%) with respect to the median (weighted mean). To investigate possible differences between the core-type populations, the emission measure profiles are stacked according to their core types by calculating the weighted mean and the median, as shown in the upper panel of Fig. 2. The lower panel of Fig. 2 shows the ratio of the average EM profiles sorted according to their core types. The statistical errors on the weighted means or medians are very small. The error bars indicated on the plot correspond to uncertainties calculated by bootstrapping the data, i.e. by shuffling the profiles with repetition, repeating the operation $10\,000$ times and computing the median and percentiles of the output distribution. The bootstrap errors thus include information on the sample variance and non-Gaussianity of the underlying distribution.

4.1 Discussion

The weighted mean profiles reveal in a model-independent way the existence of subtle differences between the galaxy cluster populations. The amplitude of this effect between SCC and NCC clusters in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range is up to $30$% and confirms the finding of Eckert et al. (2012). Compared to the heterogeneous sample of Eckert et al. (2012) HIFLUGCS is statistically complete, which confirms the result in a more robust way. If true, this finding implies that the outskirts are affected by the core-type and a detection algorithm tailored to the galaxy cluster outskirts will be more sensitive to the more abundant NCC objects, which needs to be taken into account in selection functions. In addition one could determine the statistical likelihood for the core-state of a particular galaxy cluster. However, the asymmetric bootstrap errors indicate an underlying non-Gaussian distribution or that the sample is affected by outliers. The median profiles, which are more robust against outliers, do not reveal the same trend of the emission measure ratios. As a test, we exclude the strong cool-core profile that deviates most with respect to the median within $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ (Abell 3526). This cluster also shows the smallest statistical errors in this radial range. The weighted mean ratios resembles well the trend of the median when excluding this single cluster from the analysis. This reveals that the weighted mean is driven by an outlier with small statistical errors. Therefore, we conclude that there is no indication of a systematic core-type differences in the galaxy clusters $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range. Nevertheless, this comparison is useful because outliers like Abell 3526 will affect the selection function. The investigation of a possible redshift evolution of this analysis is left to a future study. Beyond $0.5\text{\,}\mathrm{\textnormal{r}_{500}}$ the difference between SCC and NCC clusters become larger. Eckert et al. (2012) discussed gas redistribution between the core region and the outskirts as possible explanation. In this scenario, the injected energy due to a merging event flattens the density profile of interacting objects. Assuming that NCC clusters are more likely to have experienced a recent merging event, their self-similar scaled EM profiles would be different compared to CC objects. Another explanation could be the current accretion of large scale blobs. Clusters with higher mass accretion rates show a larger fraction of non-thermal pressure in simulations (Nelson et al., 2014). Again assuming that NCC objects are merging clusters, the discrepancies can be explained by the different non-thermal energy content. However, this scenario seems unlikely since we expect to detect such structures in our wavelet images. An additional explanation can be that the dark matter halos of NCC and CC objects have different shapes and thus a different concentration at a given radius.

5 Large scale center and ellipticity

The wavelet decomposition allows us to study galaxy cluster parameters on large scales. The one scale around $0.2\text{\,}\mathrm{\textnormal{r}_{500}}$ is used to determine the center and ellipticities using the SExtractor program (Bertin & Arnouts, 1996). This minimizes the impact of the different core states on these parameters. The chosen scale contains most of the cluster counts outside the core region and is therefore a good tracer for large scale properties. Including scales up to around $0.5\text{\,}\mathrm{\textnormal{r}_{500}}$ (an additional $13$ scales) shows that the median and mean difference in ellipticity is just about 10%. The ellipticity, $e$, of an ellipse with major-axis, $e_{1}$, and minor-axis, $e_{2}$, is defined as

[TABLE]

and is shown as function of the cooling time in Fig. 3. The two parameters do not show a significant correlation (Spearman rank-order correlation coefficient of $0.17$). The medians and sample standard deviations in the 3 bins are ($0.24$, $0.28$, $0.25$) and ($0.08$, $0.12$, $0.14$), respectively. There are no highly elliptical clusters with short cooling times, i.e. selecting clusters above an ellipticity of \sim$$0.3 creates a sample without SCC objects. The universality of this result needs to be confirmed with galaxy cluster samples of larger sizes.

We perform a similar stacking analysis as in Sect. 4, where the sample is divided into 3 sub-classes according to ellipticity, rather than core-state. The weighted mean and median profiles are shown in Fig. 4. The median profiles do not show a difference between the sub-classes except in the core and the very outskirts.

We quantify the covariance between ellipticity and core radius in kpc. For a core-modelled single $\beta$-model and fixed $\beta$-parameter, the Spearman rank-order correlation coefficient is $0.20$, i.e. no strong indication for a correlation. Allowing the $\beta$-parameter to vary introduces a small positive correlation coefficient of $0.33$. A similar behaviour is seen for best-fit core radii using a single $\beta$-model.

5.1 Discussion

Clusters with large cooling times show a large range of 2D ellipticities. Assuming that they originate from the same population in 3D, the large scatter of ellipticity might be explained by projection effects, i.e. triaxial halos with random orientations will yield a wide range of observed 2D ellipticities. The verification of this effect using simulations is left to a future study. Considering ellipticity adds supplementary information since it is not significantly correlated with the core radius. The ellipticity in the galaxy cluster outskirts traces the amount of baryon dissipation (Lau et al., 2012). Thus, ellipticity is linked and may help to constrain cosmological parameters like the amplitude of the matter power spectrum ($\sigma_{8}$). Halos form later for lower $\sigma_{8}$ values and are therefore on average more elliptical (Allgood et al., 2006; Macciò et al., 2008). In addition, measuring ellipticity on several scales allows to indirectly study large-scale gas rotation in galaxy clusters (Bianconi et al., 2013). This allows us to test if the ICM is in hydrostatic equilibrium, an often made assumption in X-ray mass measurements. This makes the ellipticity an interesting survey measure for the future eROSITA X-ray all-sky survey, where we expect a range of ellipticities. The ratio of the weighted mean profiles between objects with low and high ellipticity shows a slight offset. Since this offset is not visible in the median profile, it can’t solely be explained by the larger non-Gaussianity per bin of elliptical surface brightness distributions extracted in circular annuli.

Hashimoto et al. (2007) studied the relationships between the X-ray morphology and several other cluster properties using a heterogeneous sample of $101$ clusters taken from the Chandra archive, out of which $18$ objects are represented in HIFLUGCS. The ellipticity measurements are in agreement with our study within a factor of $\sim$2, except for A0399. For this cluster, the ellipticity calculated by Hashimoto et al. (2007) is a factor of $5.6$ larger ($0.284$ compared to $0.051$). We note that the semi-major and semi-minor axis are calculated in the same way, i.e. by the 2nd-order moments but we subtracted the central excess emission of the cluster before. Including the central excess emission brings our results in better agreement with Hashimoto et al. (2007), especially in the case of A0399. For a single $\beta$-model description of the surface brightness profile, they find a slightly larger Spearman rank-order correlation coefficient of $0.37$ (compared to $0.26$) between ellipticity and core radius in kpc.

6 Analysis of the residuals

In this section we quantify how well the different model parameterisations represent the underlying HIFLUGCS surface brightness profiles. All profiles are fitted over the full extracted radial range (see Fig. D.1) using an MCMC posterior sampling technique, taking uncertainties in the measured surface brightness into account. In case of the core-modelled single $\beta$-model all contaminating sources, including the excess emission in the cluster core, are modelled according to equation (5). For the single $\beta$-model fitting procedure we perform a classical approach of masking contaminating sources that are detected on the wavelet images. This includes point-like sources or extended emission like substructure but not emission from the cluster core. Thus the likelihood functions for the two cases are different because for a single $\beta$-model, $B_{\textnormal{wv},i}$ in equation (5) is equal to zero for all $i$. The choice of priors is discussed in appendix B. The background is modelled with an additive constant. We note that not subtracting the background before fitting introduces a positive degeneracy between the best-fit slope of the $\beta$-model and the background level. In this framework, residuals are defined as (data-model)/model, i.e. positive (negative) residuals indicate that the model under- (over-) predicts the observed surface brightnesses. We focus particularly on the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range because most of the galaxy cluster counts outside the non-scalable core regions are expected there.

6.1 Discussion: Single $\beta$-models

A single $\beta$-model is a widely used description to fit surface brightness profiles, especially in the low statistics regime. Therefore it is also commonly used to detect extended sources in X-ray surveys. Fig. 5 (upper left panel) shows the fractional median of the residuals from a single $\beta$-model in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range as function of core radius. There is a clear trend of positive residuals towards smaller core radii. This means that the flux in the outskirts is systematically under-predicted, especially for SCC objects (see Sect. 6.3). In addition, the core radii of objects that exhibit cool-cores are systematically biased low since a single $\beta$-model lacks degrees of freedom to model the central excess emission. Due to higher photon statistics in the core region, a single $\beta$-model fit tends to be driven by the inner radial bins. As an additional test we tried to reduce the weight of the core region in the fitting procedure by assuming that the variance in the Gaussian likelihood function is underestimated by a given fractional amount $f=0.1$. The qualitative behaviour of the residuals remains the same. In cases where the cluster outskirts are poorly fitted, the background level is not determined properly as well because its level compensates for the poor fit in the outskirts. This influences the $\beta$ value determination as already discussed. The median of the $\beta$ parameters is $0.59$ and thus smaller than the often assumed generic value of $2/3$. The best-fit values of the surface brightness parameters are in good agreement with Reiprich & Böhringer (2002).

The single $\beta$-model residuals for individual clusters have a wavy form. As additional test we investigated different functional forms to describe the common deviations from a single $\beta$-model, but no significant common form of the residuals could be found. The scatter around [math] and bias in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial bin are $0.092$ and $0.032\pm 0.003$, respectively. The non-zero bias reflects the same finding as discussed above.

The residuals in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range from a single $\beta$-model with fixed $\beta$ parameter ($\beta=2/3$) is shown in the lower left panel of Fig. 5. The amplitude of the residuals for SCC objects increases up to over 40%. Fixing the $\beta$ parameter increases the scatter and bias to $0.161$ and $0.066\pm 0.003$, respectively.

6.2 Discussion: Core-modelled single $\beta$-models

The negative effects of a single $\beta$-model as described in Sect. 6.1 are reduced when modelling the excess core emission by adding the counts on the small scales of the wavelet decomposition to the model counts in the Poisson likelihood function (see Sect. 3.1.2). The residuals in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range are shown in Fig. 5 (upper right panel). The measured residuals scatter around zero and there is no bias between individual core types visible. The scatter around [math] in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial bin is slightly reduced compared to the single $\beta$-model ($0.074$). Most importantly, the bias gets more consistent with zero ($-0.004\pm 0.003$). In case of a fixed $\beta$ parameter (lower right panel of Fig. 5), the bias of a core-modelled single $\beta$-model is consistent with zero ($-0.001\pm 0.003$) and the scatter is slightly increased ($0.087$). The median of the $\beta$ parameters is larger compared to the single $\beta$-model case and with $0.696$ close to generic value of $2/3$.

To verify our core-modelling method and to confirm its results we excise the core region ($<$0.2\text{,}\mathrm{\textnormal{r}_{500}}$$) in the single $\beta$-model fitting procedure. This is an independent way to avoid the single $\beta$-model fit to be driven by the non-scalable core emission and to reduce residuals between model and data in the cluster outskirts. In more than 90% of the cases, we find that this approach delivers comparable best-fit parameters for $\beta$ and core-radii as when modelling the core emission. Due to excluding data, the constraining power is reduced and the degeneracy of the shape parameters with the $\beta$-model normalization is larger. In several cases, this results in a larger mismatch of the $\beta$-model fluxes compared to the real flux of the objects.

6.3 Flux comparison

This section compares the overall model flux, as well as the model flux in the cluster outskirts for the single $\beta$-model and core-modelled single $\beta$-model to the measured flux in the corresponding radial range. For each model, we calculate the cluster count rate by integrating a single $\beta$-model with the corresponding best-fit parameters in a given radial range, i.e. [math]-$$r_{\textnormal{x}} (the outer significance radius) for the overall model flux and $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ for the flux in the outskirts. We are interested in flux ratios, in which the count rate to flux conversion factors of individual clusters cancel each other out. We calculate the total flux of the single $\beta$-model including the central excess emission and subtract the wavelet-detected central excess emission in the core-modelled case. The median of the measured fluxes in the outskirts and the overall fluxes of the NCC and WCC objects agree very well with each other. The measured overall single $\beta$-model fluxes of the SCC objects are on average \sim{\color[rgb]{0,0,0}23}% larger compared to the measured overall core-modelled fluxes and the flux ratio has an intrinsic scatter of \sim{\color[rgb]{0,0,0}14}% around the median. The median accuracy between model and measured total flux is within \sim{\color[rgb]{0,0,0}4}% (\sim{\color[rgb]{0,0,0}2}%) for the (core-modelled) single $\beta$-model, regardless of the core-type or if the $\beta$ value is free to vary or fixed to $2/3$. The flux in the outskirts of the core-modelled single $\beta$-model has the same median accuracy of \sim{\color[rgb]{0,0,0}2}%. For the single $\beta$-model, the accuracy stays at the \sim{\color[rgb]{0,0,0}4}% level for WCC and NCC objects. For SCC objects, for the single beta-model case, the flux in the outskirts is biased low by \sim{\color[rgb]{0,0,0}6}% and increases to \sim{\color[rgb]{0,0,0}10}% when $\beta$ is fixed to $2/3$; i.e. biased $35$ times larger than the for the core-modelled case. In all cases, the intrinsic scatter values of the ratios are below \sim{\color[rgb]{0,0,0}6}%. The Spearman rank-order correlation coefficients do not reveal a significant correlation between the flux underestimation and the cluster temperature.

7 Scaling relations

This section describes how the shape parameters scale with temperature and how they correlate with luminosity. Fig. 6 shows the core radius as a function of temperature, where the $\beta$ parameter is fixed to $2/3$ to avoid degeneracies between the surface brightness parameters. We note that the overall picture doesn’t change when fixing $\beta$ to the median of the full population or the medians of the individual core-type populations. To account for a possible temperature dependence, the core radii are divided by the square-root of the corresponding cluster temperatures (see Sect. 7.1 for more details). The figure emphasizes the systematic differences between the individual core types in the modelling with a single $\beta$-model. The discrepancies get less prominent when modelling the core region using the wavelet decomposition. In addition, the intrinsic scatter is reduced by $8$%/$11$%/$35$% for NCC/WCC/SCC objects, respectively. In both modelling cases, the Spearman rank-order correlation coefficients indicate a stronger negative ($\sim-0.7$), a mild positive ($\sim 0.3$) and no significant ($\sim 0$) correlation for the NCC, WCC and SCC populations, respectively.

We determine scaling relations as outlined in Sect. 3.4. The best-fit relations between shape parameters and temperature (Fig. 7), as well as luminosity and temperature are determined simultaneously. This allows studying the covariances between shape and luminosity from the joint fit. The best-fit values are shown in Table 1 and 2. Fig. 8 shows the correlation between the shape parameters and luminosity. The Spearman rank-order correlation coefficients between luminosity and $\beta$/core radius are 0.37/0.12 (0.32/0.22) for the (core-modelled) single $\beta$-model. When the $\beta$ value is fixed to $2/3$, the luminosity-core radius correlation coefficients become mildly negative with -0.23 (-0.22).

7.1 Discussion

The best-fit parameters of the single $\beta$-model and core-modelled single $\beta$-model agree on the $\sim$1$\sigma$ level, except that in the latter case the intrinsic scatter value of the core radius is reduced by almost a factor of $2$. In both modelling approaches the shape parameters show a positive correlation with temperature. The shape parameters of galaxy clusters are often fixed to generic values (e.g. $\beta=2/3$) or scaling relations (e.g. $r_{\textnormal{c}}\propto r_{500}$, Pacaud et al. (2018)). The assumption that the core-radius is proportional to $r_{500}$, together with equation (10) results in a self-similar scaling of $r_{\textnormal{c}}\propto T^{1/2}$. We find a core-radius-temperature relation with a marginally steeper slope of $1.04\pm 0.37$ ($0.75\pm 0.20$) for the (core-modelled) single $\beta$-model compared to this expectation. Fixing the value of $\beta$ to $2/3$ results in $0.77\pm 0.14$ ($0.50\pm 0.16$), consistent with the self-similar value when modelling the excess core emission. We study the correlation coefficients (Table 2) between the different galaxy cluster parameters by simultaneously fitting for the scaling relations and the intrinsic scatter tensor. We expect some degeneracy between the single $\beta$-model parameters in the fitting, amongst others that a larger core-radius is compensated by a steeper slope. This is reflected in the strong positive correlation between core-radius and $\beta$. We do not find a significant correlation between $\beta$ and luminosity. There is a strong negative correlation between core radius and luminosity, i.e. at a given temperature, more luminous objects tend to be more compact. For all modelling cases, the best-fit correlation coefficients are consistent on a 1$\sigma$ level. This implies, that the measured correlation between core radius and luminosity is not significantly affected by modelling the central excess emission. Neither this covariance nor the shape-temperature relations are taken into account in existing simulations for eROSITA but may play a crucial role in understanding selection effects related to the detection of clusters in X-ray surveys, which is a key ingredient for the use of X-ray galaxy clusters as a precision cosmological probe. The findings presented here can be used to perform more realistic simulations and a comparison between different sets of simulations allows to study the impact of these covariances on obtained cosmological parameters.

8 Summary

X-ray morphologies of galaxy clusters play a crucial role in the determination of the survey selection function. We compare self-similar scaled emission measure profiles of a well defined galaxy cluster sample (HIFLUGCS) above a representative temperature threshold of $3\text{\,}\mathrm{keV}$. One outlier (Abell 3526) with small statistical errors drives the weighted mean profiles of sub-populations according to different core properties towards a different behaviour of strong cool-core and non-cool-core objects in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range. Excluding this object from the analysis or calculating the median profiles reveals no systematic difference in the aforementioned radial range. We conclude that there is no indication for a correlation between the behaviour in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range and the core state, although the overall shapes of the SCC and NCC populations are different. The median SCC profile shows a larger normalization towards the center and is steeper compared to the median NCC profile. This leads to a turnover of the profile ratio at \sim$$0.3\text{\,}\mathrm{\textnormal{r}_{500}}. The difference in the center can be explained by the core state but the difference in the outskirts is still under debate. As discussed in Sect. 4.1 possible explanations are gas redistribution between the core region and the outskirts, the current accretion of large scale blobs or different shapes of the dark matter halos leading to different concentrations at a given radius. Characterizing galaxy cluster surface brightness profiles with single $\beta$-models is still state-of-the-art in the determination of selection functions. We investigate the residuals of a single $\beta$-model fit to the overall cluster profile, revealing that this description tends to underestimate the flux in the galaxy cluster outskirts for less extended clusters. Fixing the $\beta$ parameter to $2/3$ increases this effect dramatically (up to over $40\%$). In both cases, the core-radius measurement for SCC objects are biased low. In addition, the intrinsic scatter values with respect to the medians of the self-similar scaled extent parameters show a more than $1\sigma$ tension between strong and non-cool-core objects. These three effects can be minimized by adapting a wavelet decomposition based surface brightness modelling that is sensitive to the galaxy cluster outskirts and models the excess emission in the core region. Then, the fit is not driven by the local processes in the core. Compared to a single $\beta$-model approach the residuals in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ radial range are much smaller and the core radii depend much more mildly on the core state. Our method to model the excess core emission has very interesting applications for future galaxy cluster surveys, e.g. with eROSITA. The performance study for high redshift objects with small angular extent and establishing the most robust method for clusters where $r_{500}$ is unknown is left to a future work. Using wavelet decomposition allows to determine large scale ellipticities of the clusters. The ellipticity is an interesting new survey measure for eROSITA since it’s determination doesn’t require many photon counts and it adds additional information to the $\beta$-model shape parameters and core-excised luminosities. A detailed study of measuring galaxy cluster ellipticities with eROSITA and it’s implications is left to a future work.

We study how shape parameters and luminosity scale with temperature. There is no significant difference of the best-fit values between a single $\beta$-model and core-modelled single $\beta$-model, except that the intrinsic scatter of the core radius is almost twice as large for the single $\beta$-model case. The slope of the core radius-temperature relation is steeper than the self-similar prediction of $1/2$ but gets in agreement when fixing the $\beta$ parameter to $2/3$ in the surface brightness profile modelling. More interestingly, the shape parameters are covarient with luminosity, i.e. at a given temperature, more compact objects are more luminous. These covariances are usually neglected in simulations to determine the survey selection function (Pacaud et al., 2007; Clerc et al., 2018). In addition, these previous studies assumed a fixed $\beta$ value, while we find that $\beta$ is a function of temperature. Taking shape-temperature scaling relations and shape-luminosity covariances into account will lead to a more realistic set of simulated galaxy clusters and will provide a better understanding of the survey completeness.

Acknowledgements

This research was made possible by the International Max Planck Research School on Astrophysics at the Ludwig-Maximilians University Munich, as well as their funding received from the Max-Planck Society.

Appendix A Mass comparison

In this section we compare our previous emission measure ratio results of Sect. 4 to profiles that are re-scaled by a characteristic radii according to hydrostatic mass estimates by Schellenberger & Reiprich (2017). We adapt their preferred ”NFW Freeze” model, where a NFW profile (Navarro et al. 1996, 1997) is fit to the outermost measured mass profiles of Chandra observations and a concentration-mass relation is used to reduce the degrees of freedom. We are interested in the difference between individual core-types and not in the bias between the Planck and hydrostatic masses. Therefore we assume that the bias is constant for all clusters. To probe the masses at the same radii, we recalculate the hydrostatic masses at the Planck $r_{500}$ values according to Formula (31) of Schellenberger & Reiprich (2017)

[TABLE]

where $c_{500}$ denotes the NFW concentration parameter and $Y(u)=\ln(1+u)-u/(1+u)$. These recalculated masses are on average lower than the corresponding masses in the Planck catalogue. The median of the NCC cluster masses is $\sim$13% larger than that of the SCC objects (see Fig. 9), increasing the SCC to NCC weighted mean ratio effect in the $0.20.5\text{\,}\mathrm{\textnormal{r}_{500}}$ range to $40$%. This increase due to the dynamical states is also seen in simulations, e.g. in different fractions of non-thermal pressure for different mass accretion rates (Nelson et al. 2014). The differences can thus be explained by assuming that NCC clusters are merging objects and that they contain more non-thermal energy. An alternative explanation is that the differences come from the assumed concentration-mass relation. We expect different core-types to have different shapes of the dark matter halo and therefore different concentrations at a given radius. However, fitting for the concentration in Schellenberger & Reiprich (2017) results in some cases to unrealistic high or low masses, potentially because of limited radial coverage due to the relatively small Chandra field-of-view.

Appendix B Priors

Table 3 shows the priors used for this analysis. They are chosen to be weakly- or non-informative and varying them does not influence the results of this study significantly.

Appendix C HIFLUGCS parameters

666*†*

See pages - of ./anc/appendix.pdf