We present the first cluster catalog extracted from combined space-based (
Planck
) and ground-based (South Pole Telescope; SPT-SZ) millimeter data. We developed and applied a matched multi-filter ...(MMF) capable of dealing with the different transfer functions and resolutions of the two datasets. We verified that it produces results consistent with publications from
Planck
and SPT collaborations when applied on the datasets individually. We also verified that
Planck
and SPT-SZ cluster fluxes are consistent with each other. When applied blindly to the combined dataset, the MMF generated a catalog of 419 detections (
S
/
N
> 5), of which 323 are already part of the SPT-SZ or PSZ2 catalogs; 54 are new SZ detections, which have been identified in other catalogs or surveys; and 42 are new unidentified candidates. The MMF takes advantage of the complementarity of the two datasets,
Planck
being particularly useful for detecting clusters at a low redshift (
z
< 0.3), while SPT is efficient at finding higher redshift (
z
> 0.3) sources. This work represents a proof of concept that blind cluster extraction can be performed on combined, inhomogeneous millimeter datasets acquired from space and ground. This result is of prime importance for planned ground-based cosmic microwave background (CMB) experiments (e.g., Simons Observatory, CMB-S4) and envisaged CMB space missions (e.g., PICO, Backlight) that will detect hundreds of thousands of clusters in the low mass regime (
M
500
≤ 10
14
M
⊙
), for which the various sources of intra-cluster emission (gas, dust, synchrotron) will be of the same order of magnitude and hence require broad ground and space frequency coverage with a comparable spatial resolution for adequate separation.
We study the dynamical state and the integrated total mass profiles of 75 massive (M500 > 5 × 1014 M⊙) Sunyaev–Zeldovich(SZ)-selected clusters at 0.08 < z < 1.1. The sample is built from the Planck ...catalogue, with the addition of four SPT clusters at z > 0.9. Using XMM-Newton imaging observations, we characterise the dynamical state with the centroid shift ⟨w⟩, the concentration CSB, and their combination, M, which simultaneously probes the core and the large-scale gas morphology. Using spatially resolved spectroscopy and assuming hydrostatic equilibrium, we derive the total integrated mass profiles. The mass profile shape is quantified by the sparsity, that is the ratio of M500 to M2500, the masses at density contrasts of 500 and 2500, respectively. We study the correlations between the various parameters and their dependence on redshift. We confirm that SZ-selected samples, thought to most accurately reflect the underlying cluster population, are dominated by disturbed and non-cool core objects at all redshifts. There is no significant evolution or mass dependence of either the cool core fraction or the centroid shift parameter. The M parameter evolves slightly with z, having a correlation coefficient of ρ = −0.2 ± 0.1 and a null hypothesis p-value of 0.01. In the high-mass regime considered here, the sparsity evolves minimally with redshift, increasing by 10% between z < 0.2 and z > 0.55, an effect that is significant at less than 2σ. In contrast, the dependence of the sparsity on dynamical state is much stronger, increasing by a factor of ∼60% from the one third most relaxed to the one third most disturbed objects, an effect that is significant at more than 3σ. This is the first observational evidence that the shape of the integrated total mass profile in massive clusters is principally governed by the dynamical state and is only mildly dependent on redshift. We discuss the consequences for the comparison between observations and theoretical predictions.
Abstract We investigate the evolution of the dark matter density profiles of the most massive galaxy clusters in the Universe. Using a ‘zoom-in’ procedure on a large suite of cosmological simulations ...of total comoving volume of 3 (h − 1 Gpc)3, we study the 25 most massive clusters in four redshift slices from z ∼ 1 to the present. The minimum mass is M500 > 5.5 × 1014 M⊙ at z = 1. Each system has more than two million particles within r500. Once scaled to the critical density at each redshift, the dark matter profiles within r500 are strikingly similar from z ∼ 1 to the present day, exhibiting a low dispersion of 0.15 dex, and showing little evolution with redshift in the radial logarithmic slope and scatter. They have the running power-law shape typical of the Navarro–Frenk–White type profiles, and their inner structure, resolved to 3.8 h−1 comoving kpc at z = 1, shows no signs of converging to an asymptotic slope. Our results suggest that this type of profile is already in place at z > 1 in the highest-mass haloes in the Universe, and that it remains exceptionally robust to merging activity.
In high-energy astronomy, spectro-imaging instruments such as X-ray detectors allow investigation of the spatial and spectral properties of extended sources including galaxy clusters, galaxies, ...diffuse interstellar medium, supernova remnants, and pulsar wind nebulae. In these sources, each physical component possesses a different spatial and spectral signature, but the components are entangled. Extracting the intrinsic spatial and spectral information of the individual components from this data is a challenging task. Current analysis methods do not fully exploit the 2D-1D (x, y, E) nature of the data, as spatial information is considered separately from spectral information. Here we investigate the application of a blind source separation (BSS) algorithm that jointly exploits the spectral and spatial signatures of each component in order to disentangle them. We explore the capabilities of a new BSS method (the general morphological component analysis; GMCA), initially developed to extract an image of the cosmic microwave background from Planck data, in an X-ray context. The performance of the GMCA on X-ray data is tested using Monte-Carlo simulations of supernova remnant toy models designed to represent typical science cases. We find that the GMCA is able to separate highly entangled components in X-ray data even in high-contrast scenarios, and can extract the spectrum and map of each physical component with high accuracy. A modification of the algorithm is proposed in order to improve the spectral fidelity in the case of strongly overlapping spatial components, and we investigate a resampling method to derive realistic uncertainties associated to the results of the algorithm. Applying the modified algorithm to the deep Chandra observations of Cassiopeia A, we are able to produce detailed maps of the synchrotron emission at low energies (0.6–2.2 keV), and of the red- and blueshifted distributions of a number of elements including Si and Fe K.
Context.
Galaxy clusters grow through the accretion of mass over cosmic time. Their observed properties are then shaped by how baryons distribute and energy is diffused. Thus, a better understanding ...of spatially resolved, projected thermodynamic properties of the intra-cluster medium (ICM) may provide a more consistent picture of how mass and energy act locally in shaping the X-ray observed quantities of these massive virialized or still collapsing structures.
Aims.
We study the perturbations in the temperature (and density) distribution to evaluate and characterize the level of inhomogeneities and the related dynamical state of the ICM.
Methods.
We obtain and analyze the temperature and density distribution for 28 clusters (2.4 × 10
14
M
⊙
<
M
500
< 1.2 × 10
15
M
⊙
; 0.07 <
z
< 0.45) selected from the CHEX-MATE sample. We use these spatially resolved two-dimensional distributions to measure the global and radial scatter and identify the regions that deviate the most from the average distribution. During this process, we introduce three dynamical state estimators and produce “clean” temperature profiles after removing the most deviant regions.
Results.
We find that the temperature distribution of most of the clusters is skewed towards high temperatures and is well described by a log-normal function. There is no indication that the number of regions deviating more than 1
σ
from the azimuthal value is correlated with the dynamical state inferred from morphological estimators. The removal of these regions leads to local temperature variations up to 10–20% and an average increase of ∼5% in the overall cluster temperatures. The measured relative intrinsic scatter within
R
500
,
σ
T, int
/
T
, has values of 0.17
−0.05
+0.08
, and is almost independent of the cluster mass and dynamical state. Comparing the scatter of temperature and density profiles to hydrodynamic simulations, we constrain the average Mach number regime of the sample to
Ṁ
3D
= 0.36
−0.09
+0.16
. We infer the ratio between the energy in turbulence and the thermal energy, and translate this ratio in terms of a predicted hydrostatic mass bias
b
, estimating an average value of
b
∼ 0.11 (covering a range between 0 and 0.37) within
R
500
.
Conclusions.
This study provides detailed temperature fluctuation measurements for 28 CHEX-MATE clusters which can be used to study turbulence, derive the mass bias, and make predictions on the scaling relation properties.
We present a study of the structural and scaling properties of the gas distributions in the intracluster medium (ICM) of 31 nearby ($z < 0.2$) clusters observed with XMM-Newton, which together ...comprise the Representative XMM-Newton Cluster Structure Survey (REXCESS). In contrast to previous studies, this sample is unbiased with respect to X-ray surface brightness and cluster dynamical state, and it fully samples the cluster X-ray luminosity function. The clusters cover a temperature range of 2.0-8.5 keV and possess a variety of morphologies. The sampling strategy allows us to compare clusters with a wide range of central cooling times on an equal footing. We applied a recently developed technique for the deprojection and PSF-deconvolution of X-ray surface brightness profiles to obtain non-parametric gas-density profiles out to distances ranging between $0.8~R_{500}$ and $1.5~R_{500}$. We scaled the gas density distributions to allow for the systems' differing masses and redshifts. The central gas densities differ greatly from system to system, with no clear correlation with system temperature. At intermediate radii (~$0.3~R_{500}$), the scaled density profiles show much less scatter, with a clear dependence on system temperature. We find that the density at this radius scales proportionally to the square root of temperature, consistent with the presence of an entropy excess as suggested in previous literature. However, at larger scaled radii this dependence becomes weaker: clusters with $kT > 3$ keV scale self-similarly, with no temperature dependence of gas-density normalisation. The REXCESS sample allows us to investigate the correlations between cluster properties and dynamical state. We find no evidence of correlations between cluster dynamical state and either the gas density slope in the inner regions or temperature, but do find some evidence of a correlation between dynamical state and outer gas density slope. We also find a weak correlation between dynamical state and both central gas normalisation and inner cooling times, but this is only significant at the 10% level. We conclude that, for the X-ray cluster population as a whole, both the central gas properties and the angle-averaged, large-scale gas properties are linked to the cluster dynamical state. We also investigate the central cooling times of the clusters. While the cooling times span a wide range, we find no evidence of a significant bimodality in the distributions of central density, density gradient, or cooling time. Finally, we present the gas mass-temperature relation for the REXCESS sample, finding that $h(z)M_{\rm gas} \propto T^{1.99\pm0.11}$, which is consistent with the expectation of self-similar scaling modified by the presence of an entropy excess in the inner regions of the cluster and consistent with earlier work on relaxed cluster samples. We measure a logarithmic intrinsic scatter in this relation of ~$10\%$, which should be a good measure of the intrinsic scatter in the Mgas-T relation for the cluster population as a whole.
The “Cluster HEritage project with XMM-Newton : Mass Assembly and Thermodynamics at the End point of structure formation” (CHEX-MATE) is a multi-year heritage program to obtain homogeneous XMM-Newton ...observations of a representative sample of 118 galaxy clusters. The observations are tuned to reconstruct the distribution of the main thermodynamic quantities of the intra-cluster medium up to R 500 and to obtain individual mass measurements, via the hydrostatic-equilibrium equation, with a precision of 15−20%. Temperature profiles are a necessary ingredient for the scientific goals of the project and it is thus crucial to derive the best possible temperature measurements from our data. This is why we have built a new pipeline for spectral extraction and analysis of XMM-Newton data, based on a new physically motivated background model and on a Bayesian approach with Markov chain Monte Carlo methods, which we present in this paper for the first time. We applied this new method to a subset of 30 galaxy clusters representative of the CHEX-MATE sample and show that we can obtain reliable temperature measurements up to regions where the source intensity is as low as 20% of the background, keeping systematic errors below 10%. We compare the median profile of our sample and the best-fit slope at large radii with literature results and we find a good agreement with other measurements based on XMM-Newton data. Conversely, when we exclude the most contaminated regions, where the source intensity is below 20% of the background, we find significantly flatter profiles, in agreement with predictions from numerical simulations and independent measurements with a combination of Sunyaev–Zeldovich and X-ray imaging data.
Aims. Our goal is to provide a robust estimate of the metal content of the intracluster medium (ICM) in massive clusters. Methods. We made use of published abundance profiles for a sample of ~60 ...nearby systems. We included in our estimate uncertainties associated with the measurement process and with the almost total lack of information in cluster outskirts. Results. We performed a first, albeit rough, census of metals and find that the mean abundance of the ICM within r180 is very poorly constrained, 0.06 Z⊙ ≲ Z ≲ 0.26 Z⊙ and presents no disagreement with expectations. Similarly, whether and how the bulk of the metal content in clusters varies with cosmic time are very much open questions. Conclusions. A solid estimate of abundances in cluster outskirts could be achieved by combining observations of the two experiments that will operate on board Athena, the XIFU and the WFI, provided they do not fall victim to the de-scoping process that has afflicted several space observatories over the last decade.
Using XMM-Newton observations, we investigate the scaling and structural properties of the ICM entropy in a sample of 10 nearby ($z < 0.2$) morphologically relaxed galaxy clusters in the temperature ...range 2-9 keV. We derive the local entropy-temperature ($S{-}T$) relation at $R = 0.1, 0.2, 0.3$ and $0.5R_{200}$. The logarithmic slope of the relation is the same within the $1\sigma$ error at all scaled radii. However, the intrinsic dispersion about the best fitting relation is significantly higher at $0.1R_{200}$. The slope is $0.64\pm0.11$ at $0.3\,R_{200}$, in excellent agreement with previous work. We also investigate the entropy-mass relation at density contrasts $\delta=5000, 2500$ and 1000. We find a shallower slope than that expected in simple self-similar models, which is in agreement with the observed empirically-determined entropy-temperature and mass-temperature scaling. The dispersion is smaller than for the $S{-}T$ relation. Once scaled appropriately, the entropy profiles appear similar beyond ~$0.1R_{200}$, with an intrinsic dispersion of ~15 per cent and a shape consistent with gravitational heating ($S(r) \lower 3pt \hbox{$\, \buildrel {\textstyle \propto}\over {\textstyle \sim}\,$}r^{1.1}$). However, the scatter in scaled entropy profiles increases with smaller scaled radius, to more than 60 per cent at $R \la 0.05 R_{200}$. Our results are in qualitative agreement with models which boost entropy production at the accretion shock. However, localised entropy modification may be needed to explain the dispersion in the inner regions.