We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. In this work, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant ($z > 1.2$) SPT-selected clusters, and compare them with those of the nearby clusters ($z<0.1$) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

We simulate stacked observations of nearby hot X-ray coronae associated with galaxies in the EAGLE and Illustris-TNG hydrodynamic simulations. A forward modeling pipeline is developed to predict 4-year eROSITA observations and stacked image analysis, including the effects of instrumental and astrophysical backgrounds. We propose an experiment to stack z~0.01 galaxies separated by specific star-formation rate (sSFR) to examine how the hot (T>=10^6 K) circumgalactic medium (CGM) differs for high- and low-sSFR galaxies. The simulations indicate that the hot CGM of low-mass (M_*~10^{10.5} Msol), high-sSFR (defined as the top one-third ranked by sSFR) central galaxies will be detectable to a galactocentric radius r~30-50 kpc. Both simulations predict lower luminosities at fixed stellar mass for the low-sSFR galaxies (the lower third of sSFR) with Illustris-TNG predicting 3x brighter coronae around high-sSFR galaxies than EAGLE. Both simulations predict detectable emission out to r~150-200 kpc for stacks centered on high-mass (M_*~10^{11.0} Msol) galaxies, with EAGLE predicting brighter X-ray halos. The extended soft X-ray luminosity correlates strongly and positively with the mass of circumgalactic gas within the virial radius (f_{CGM}). Prior analyses of both simulations have established that f_{CGM} is reduced by expulsive feedback driven mainly by black hole growth, which quenches galaxy growth by inhibiting replenishment of the ISM. Both simulations predict that eROSITA stacks should not only conclusively detect and resolve the hot CGM around L^* galaxies for the first time, but provide a powerful probe of how the baryon cycle operates, for which there remains an absence of consensus between state-of-the-art simulations.

The particle-induced background of X-ray observatories is produced by Galactic Cosmic Ray (GCR) primary protons, electrons, and He ions. Events due to direct interaction with the detector are usually removed by on board processing. The interactions of these primary particles with the detector environment produce secondary particles that mimic X-ray events from celestial sources and are much more difficult to identify. The filter wheel closed data from the XMM-Newton EPIC-pn camera in small window mode (SWM) contains both the X-ray-like background events and the events due to direct interactions with the primary particles. From this data we demonstrate that X-ray-like background events are spatially correlated with the primary particle interaction. This result can be used to further characterise and reduce the non-X-ray background in silicon-based X-ray detectors in current and future missions. We also show that spectrum and pattern fractions of secondary particle events are different from those produced by cosmic X-rays.

One key ingredient in using galaxy clusters (GCs) as a precision cosmological probe in large X-ray surveys is to understand selection effects. The dependence of the X-ray emission on the square of the gas density leads to a predominant role of cool cores in the detection of GCs. The contribution of cool cores to the X-ray luminosity does not scale with GC mass and cosmology and therefore affects the use of X-ray GCs in producing cosmological constraints. One of the main science goals of the eROSITA mission is to constrain cosmology with a wide X-ray survey. We propose an eROSITA GC detection scheme that avoids the use of X-ray GC centers in detection. We calculate theoretical expectations and characterize the performance of this scheme by simulations. Performing realistic simulations of point sources (PSs) in survey mode we search for spatial scales where the extended signal is uncontaminated by the PS flux. We derive a combination of scales and thresholds, which result in a clean extended source catalog. We design the output of the GC detection which enables calibrating the core-excised luminosity using external mass measurements. We provide a way to incorporate the results of this calibration in the production of final core-excised luminosity. Similarly to other GC detection pipelines, we sample the flux - core radius detection space of our method and find many similarities with the pipeline used in the 400d survey. Both detection methods require large statistics on compact GCs, in order to reduce the contamination from PSs. The benefit of our pipeline consists in the sensitivity to the outer GC shapes, which are characterized by large core sizes with little GC to GC variation at a fixed total mass. GC detection through cluster outskirts improves the GC characterization using eROSITA survey data and is expected to yield well characterized GC catalogs having simple selection functions.

We present the constraints on the helium abundance in 12 X-ray luminous galaxy clusters that have been mapped in their X-ray and Sunyaev-Zeldovich (SZ) signals out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). The unprecedented precision available for the estimate of $H_0$ allows us to investigate how much the reconstructed X-ray and SZ signals are consistent with the expected ratio $x$ between helium and proton densities of 0.08-0.1. We find that a $H_0$ around 70 km/s/Mpc is preferred from our measurements, with lower values of $H_0$ as requested from the Planck collaboration (67 km/s/Mpc) requiring a 34% higher value of $x$. On the other hand, higher values of $H_0$, as obtained by measurements in the local universe, impose $x$, from the primordial nucleosynthesis calculations and current solar abundances, reduced by 37--44\%.

The galaxy cluster Zwicky 3146 is a sloshing cool core cluster at $z=0.291$ that in X-ray imaging does not appear to exhibit significant pressure substructure in the intracluster medium (ICM). The published $M_{500}$ values range between $3.88^{+0.62}_{-0.58}$ to $22.50 \pm 7.58 \times 10^{14}$ M$_{\odot}$, where ICM-based estimates with reported errors $<20$\% suggest that we should expect to find a mass between $6.53^{+0.44}_{-0.44} \times 10^{14}$ M$_{\odot}$ (from Planck, with an $8.4\sigma$ detection) and $8.52^{+1.77}_{-1.47} \times 10^{14}$ M$_{\odot}$ (from ACT, with a $14\sigma$ detection). This broad range of masses is suggestive that there is ample room for improvement for all methods. Here, we investigate the ability to estimate the mass of Zwicky 3146 via the Sunyaev-Zel'dovich (SZ) effect with data taken at 90 GHz by MUSTANG-2 to a noise level better than $15\ \mu$K at the center, and a cluster detection of $104\sigma$. We derive a pressure profile from our SZ data which is in excellent agreement with that derived from X-ray data. From our SZ-derived pressure profiles, we infer $M_{500}$ and $M_{2500}$ via three methods -- $Y$-$M$ scaling relations, the virial theorem, and hydrostatic equilibrium -- where we employ X-ray constraints from \emph{XMM-Newton} on the electron density profile when assuming hydrostatic equilibrium. Depending on the model and estimation method, our $M_{500}$ estimates range from $6.23 \pm 0.59$ to $10.6 \pm 0.95 \times 10^{14}$ M$_{\odot}$, where our estimate from hydrostatic equilibrium, is $8.29^{+1.93}_{-1.24}$ ($\pm 19.1$\% stat) ${}^{+0.74}_{-0.68}$ ($\pm 8.6$\% sys, calibration) $\times 10^{14}$ M$_{\odot}$. Our fiducial mass, derived from a $Y$-$M$ relation is $8.16^{+0.44}_{-0.54}$ ($\pm 5.5$\% stat) ${}^{+0.46}_{-0.43}$ ($\pm 5.5$\% sys, $Y$-$M$) ${}^{+0.59}_{-0.55}$ ($\pm 7.0$\% sys, cal.) $\times 10^{14}$ M$_{\odot}$.

In this work, we investigate the relation between the radially-resolved thermodynamic quantities of the intracluster medium in the X-COP cluster sample, aiming to assess the stratification properties of the ICM. We model the relations between radius, gas temperature, density and pressure using a combination of power-laws, also evaluating the intrinsic scatter in these relations. We show that the gas pressure is remarkably well correlated to the density, with very small scatter. Also, the temperature correlates with gas density with similar scatter. The slopes of these relations have values that show a clear transition from the inner cluster regions to the outskirts. This transition occurs at the radius $r_t = 0.19(\pm0.04)R_{500}$ and electron density $n_t = (1.91\pm0.21)\cdot10^{-3} cm^{-3} E^2 (z)$. We find that above 0.2 $R_{500}$ the radial thermodynamic profiles are accurately reproduced by a well defined and physically motivated framework, where the dark matter follows the NFW potential and the gas is represented by a polytropic equation of state. By modeling the gas temperature dependence upon both the gas density and radius, we propose a new method to reconstruct the hydrostatic mass profile based only on the quite inexpensive measurement of the gas density profile.

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

As the largest virialized structures in the universe, galaxy clusters continue to grow and accrete matter from the cosmic web. Due to the low gas density in the outskirts of clusters, measurements are very challenging, requiring extremely sensitive telescopes across the entire electromagnetic spectrum. Observations using X-rays, the Sunyaev-Zeldovich effect, and weak lensing and galaxy distributions from the optical band, have over the last decade helped to unravel this exciting new frontier of cluster astrophysics, where the infall and virialization of matter takes place. Here, we review the current state of the art in our observational and theoretical understanding of cluster outskirts, and discuss future prospects for exploration using newly planned and proposed observatories.

We present the reconstruction of hydrostatic mass profiles in 13 X-ray luminous galaxy clusters that have been mapped in their X-ray and SZ signal out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). Using profiles of the gas temperature, density and pressure that have been spatially resolved out to (median value) 0.9 $R_{500}$, 1.8 $R_{500}$, and 2.3 $R_{500}$, respectively, we are able to recover the hydrostatic gravitating mass profile with several methods and using different mass models. The hydrostatic masses are recovered with a relative (statistical) median error of 3% at $R_{500}$ and 6% at $R_{200}$. By using several different methods to solve the equation of the hydrostatic equilibrium, we evaluate some of the systematic uncertainties to be of the order of 5% at both $R_{500}$ and $R_{200}$. A Navarro-Frenk-White profile provides the best-fit in nine cases out of 13, with the remaining four cases that do not show a statistically significant tension with it. The distribution of the mass concentration follows the correlations with the total mass predicted from numerical simulations with a scatter of 0.18 dex, with an intrinsic scatter on the hydrostatic masses of 0.15 dex. We compare them with the estimates of the total gravitational mass obtained through X-ray scaling relations applied to $Y_X$, gas fraction and $Y_{SZ}$, and from weak lensing and galaxy dynamics techniques, and measure a substantial agreement with the results from scaling laws, from WL at both $R_{500}$ and $R_{200}$ (with differences below 15%), from cluster velocity dispersions, but a significant tension with the caustic masses that tend to underestimate the hydrostatic masses by 40% at $R_{200}$. We also compare these measurements with predictions from alternative models to the Cold Dark Matter, like the "Emergent Gravity" and MOND scenarios.

The hot plasma in galaxy clusters is expected to be heated to high temperatures through shocks and adiabatic compression. The thermodynamical properties of the gas encode information on the processes leading to the thermalization of the gas in the cluster's potential well as well as non-gravitational processes such as gas cooling, AGN feedback and kinetic energy. In this work we present the radial profiles of the thermodynamic properties of the intracluster medium (ICM) out to the virial radius for a sample of 12 galaxy clusters selected from the Planck all-sky survey. We determine the universal profiles of gas density, temperature, pressure, and entropy over more than two decades in radius. We exploit jointly X-ray information from XMM and Sunyaev-Zel'dovich constraints from Planck to recover thermodynamic properties out to 2 R500. We provide average functional forms for the radial dependence of the main quantities and quantify the slope and intrinsic scatter of the population as a function of radius. We find that gas density and pressure profiles steepen steadily with radius, in excellent agreement with previous observational results. Entropy profiles beyond R500 closely follow the predictions for the gravitational collapse of structures. The scatter in all thermodynamical quantities reaches a minimum in the range [0.2-0.8] R500 and increases outwards. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region R > 0.3 R500) follow remarkably well the predictions of simple gravitational collapse and require little non-gravitational corrections.

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons.

We present a new cosmological probe for galaxy clusters, the halo sparsity. This characterises halos in terms of the ratio of halo masses measured at two different radii and carries cosmological information encoded in the halo mass profile. Building upon the work of Balmes et al. (2014) we test the properties of the sparsity using halo catalogs from a numerical N-body simulation of ($2.6$ Gpc/h)$^3$ volume with $4096^3$ particles. We show that at a given redshift the average sparsity can be predicted from prior knowledge of the halo mass function. This provides a quantitative framework to infer cosmological parameter constraints using measurements of the sparsity of galaxy clusters. We show this point by performing a likelihood analysis of synthetic datasets with no systematics, from which we recover the input fiducial cosmology. We also perform a preliminary analysis of potential systematic errors and provide an estimate of the impact of baryonic effects on sparsity measurements. We evaluate the sparsity for a sample of 104 clusters with hydrostatic masses from X-ray observations and derive constraints on the cosmic matter density $\Omega_m$ and the normalisation amplitude of density fluctuations at the $8$ Mpc h$^{-1}$ scale, $\sigma_8$. Assuming no systematics, we find $\Omega_m=0.42\pm 0.17$ and $\sigma_8=0.80\pm 0.31$ at $1\sigma$, corresponding to $S_8\equiv \sigma_8\sqrt{\Omega_m}=0.48\pm 0.11$. Future cluster surveys may provide opportunities for precise measurements of the sparsity. A sample of a few hundreds clusters with mass estimate errors at a few percent level can provide competitive cosmological parameter constraints complementary to those inferred from other cosmic probes.

We present the joint analysis of the X-ray and SZ signals in A2319, the galaxy cluster with the highest signal-to-noise ratio in Planck maps and that has been surveyed within our XMM Cluster Outskirts Project (X-COP). We recover the thermodynamical profiles by the geometrical deprojection of the X-ray surface brightness, of the SZ comptonization parameter, and an accurate and robust spectroscopic measurements of the temperature. We resolve the clumpiness of the density to be below 20 per cent demonstrating that most of this clumpiness originates from the ongoing merger and can be associated to large-scale inhomogeneities. This analysis is done in azimuthally averaged radial bins and in eight independent angular sectors, enabling us to study in details the azimuthal variance of the recovered properties. Given the exquisite quality of the X-ray and SZ datasets, we constrain at $R_{200}$ the total hydrostatic mass, modelled with a NFW profile, with very high precision ($M_{200} = 9.76 \pm 0.16^{stat.} \pm 0.31^{syst.} \times 10^{14} M_\odot$). We identify the ongoing merger and how it is affecting differently the gas properties in the resolved azimuthal sectors. We have several indications that the merger has injected a high level of non-thermal pressure in this system: the clumping free density profile is above the average profile obtained by stacking Rosat observations; the gas mass fraction exceeds the expected cosmic gas fraction beyond $R_{500}$; the pressure profile is flatter than the fit obtained by the Planck collaboration; the entropy profile is flatter than the mean one predicted from non-radiative simulations; the analysis in azimuthal sectors has revealed that these deviations occur in a preferred region of the cluster. All these tensions are resolved by requiring a relative support of about 40 per cent from non-thermal to the total pressure at $R_{200}$.

We present the radial distribution of the dark matter in two massive, X-ray luminous galaxy clusters, Abell~2142 and Abell~2319, and compare it with the quantity predicted as apparent manifestation of the baryonic mass in the context of the "Emergent Gravity" scenario, recently suggested from Verlinde (2016). Thanks to the observational strategy of the \xmm\ Cluster Outskirt Programme (X-COP), using the X-ray emission mapped with \xmm\ and the SZ signal in the Planck survey, we recover the gas density, temperature and thermal pressure profiles up to $\sim R_{200}$, allowing to constrain at unprecedented level the total mass through the hydrostatic equilibrium equation. We show that, also including systematic uncertainties related to the X-ray based mass modelling, the apparent "dark" matter shows a radial profile that has a shape different from the traditional dark matter distribution, with larger discrepancies (by a factor 2--3) in the inner ($r<200$ kpc) cluster's regions and a remarkable agreement only across $R_{500}$.

Galaxy clusters are the most recent products of hierarchical accretion over cosmological scales. The gas accreted from the cosmic field is thermalized inside the cluster halo. Gas entropy and pressure are expected to have a self-similar behaviour with their radial distribution following a power law and a generalized Navarro-Frenk-White profile, respectively. This has been shown also in many different hydrodynamical simulations. We derive the spatially-resolved thermodynamical properties of 47 X-ray galaxy clusters observed with Chandra in the redshift range 0.4 < z < 1.2, the largest sample investigated so far in this redshift range with X-rays spectroscopy, with a particular care in reconstructing the gas entropy and pressure radial profiles. We search for deviation from the self-similar behaviour and look for possible evolution with redshift. The entropy and pressure profiles lie very close to the baseline prediction from gravitational structure formation. We show that these profiles deviate from the baseline prediction as function of redshift, in particular at z > 0.75, where, in the central regions, we observe higher values of the entropy (by a factor of 2.2) and systematically lower estimates (by a factor of 2.5) of the pressure. The effective polytropic index, which retains informations about the thermal distribution of the gas, shows a slight linear positive evolution with the redshift and the concentration of the dark matter distribution. A prevalence of non-cool-core, disturbed systems, as we observe at higher redshifts, can explain such behaviours.