Detecting shocked intergalactic gas with X-ray and radio observations

Detecting the thermal and non-thermal emission from the shocked cosmic gas surrounding large-scale structures represents a challenge for observations, as well as a unique window into the physics of the warm-hot intergalactic medium. In this work, we present synthetic radio and X-ray surveys of large cosmological simulations in order to assess the chances of jointly detecting the cosmic web in both frequency ranges. We then propose best observing strategies tailored for existing (LOFAR, MWA and XMM) or future instruments (SKA-LOW and SKA-MID, ATHENA and eROSITA). We find that the most promising targets are the extreme peripheries of galaxy clusters in an early merging stage, where the merger causes the fast compression of warm-hot gas onto the virial region. By taking advantage of a detection in the radio band, future deep X-ray observations will probe this gas in emission, and help us to study plasma conditions in the dynamic warm-hot intergalactic medium with unprecedented detail.Comment: 22 pages, 25 Figures. A\&A accepted, in press. Moderate revision compared to version 1, with a few new figure

Large-scale inhomogeneities of the intracluster medium: improving mass estimates using the observed azimuthal scatter

Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study the ICM of inhomogeneities, focusing on the ones on the large scale that, unlike clumps, are the most difficult to identify. To this purpose we introduce the concept of residual clumpiness, C_R, that quantifies the large-scale inhomogeneity of the ICM. After showing that this quantity can be robustly defined for relaxed systems, we characterize how it varies with radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces an overestimate in the determination of the density profile from the X-ray emission, which translates into a systematic overestimate of 6 (12)% in the measurement of M_gas at R_200 for our relaxed (perturbed) cluster sample. At the same time, the increase of C_R with radius introduces also a ~2% systematic underestimate in the measurement of the hydrostatic-equilibrium mass (M_he), which adds to the previous one generating a systematic ~8.5% overestimate in f_gas in our relaxed sample. Since the residual clumpiness of the ICM is not directly observable, we study its correlation with the azimuthal scatter in the X-ray surface brightness of the halo and in the y-parameter profiles. We find that their correlation is highly significant (r_S = 0.6-0.7), allowing to define the azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as robust proxies of C_R. After providing a function that connects the two quantities, we obtain that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias in the measurement of M_gas for relaxed objects, which becomes (0+/-2)% up to 2R_200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of M_he and f_gas, although a significant halo to halo scatter remains. (abridged)Comment: 18 pages, 17 figures, 3 tables. Submitted to MNRAS, revised after referee's comment

Simulated X-ray galaxy clusters at the virial radius: slopes of the gas density, temperature and surface brightness profiles

Using a set of hydrodynamical simulations of 9 galaxy clusters with masses in the range 1.5 10^{14} M_sun < M_vir < 3.4 10^{15} M_sun, we have studied the density, temperature and X-ray surface brightness profiles of the intracluster medium in the regions around the virial radius. We have analyzed the profiles in the radial range well above the cluster core, the physics of which are still unclear and matter of tension between simulated and observed properties, and up to the virial radius and beyond, where present observations are unable to provide any constraints. We have modeled the radial profiles between 0.3 R_200 and 3 R_200 with power laws with one index, two indexes and a rolling index. The simulated temperature and [0.5-2] keV surface brightness profiles well reproduce the observed behaviours outside the core. The shape of all these profiles in the radial range considered depends mainly on the activity of the gravitational collapse, with no significant difference among models including extraphysics. The profiles steepen in the outskirts, with the slope of the power-law fit that changes from -2.5 to -3.4 in the gas density, from -0.5 to -1.8 in the gas temperature, and from -3.5 to -5.0 in the X-ray soft surface brightness. We predict that the gas density, temperature and [0.5-2] keV surface brightness values at R_200 are, on average, 0.05, 0.60, 0.008 times the measured values at 0.3 R_200. At 2 R_200, these values decrease by an order of magnitude in the gas density and surface brightness, by a factor of 2 in the temperature, putting stringent limits on the detectable properties of the intracluster-medium (ICM) in the virial regions.Comment: 13 pages, 6 figures; added reference and other minor change

Gas clumping in galaxy clusters

The reconstruction of galaxy cluster's gas density profiles is usually performed by assuming spherical symmetry and averaging the observed X-ray emission in circular annuli. In the case of a very inhomogeneous and asymmetric gas distribution, this method has been shown to return biased results in numerical simulations because of the $n^2$ dependence of the X-ray emissivity. We propose a method to recover the true density profiles in the presence of inhomogeneities, based on the derivation of the azimuthal median of the surface brightness in concentric annuli. We demonstrate the performance of this method with numerical simulations, and apply it to a sample of 31 galaxy clusters in the redshift range 0.04-0.2 observed with ROSAT/PSPC. The clumping factors recovered by comparing the mean and the median are mild and show a slight trend of increasing bias with radius. For $R<R_{500}$, we measure a clumping factor $\sqrt{C}<1.1$, which indicates that the thermodynamic properties and hydrostatic masses measured in this radial range are only mildly affected by this effect. Comparing our results with three sets of hydrodynamical numerical simulations, we found that non-radiative simulations significantly overestimate the level of inhomogeneities in the ICM, while the runs including cooling, star formation, and AGN feedback reproduce the observed trends closely. Our results indicate that most of the accretion of X-ray emitting gas is taking place in the diffuse, large-scale accretion patterns rather than in compact structures.Comment: 12 pages, 11 figures, accepted for publication in MNRAS. Largely-improved version compared to v1, method and comparison with simulations update

The effect of feedback on the emission properties of the Warm-Hot Intergalactic Medium

At present, 30-40 per cent of the baryons in the local Universe is still undetected. According to theoretical predictions, this gas should reside in filaments filling the large-scale structure (LSS) in the form of a Warm-Hot Intergalactic Medium (WHIM), at temperatures of 10^5 - 10^7 K, thus emitting in the soft X-ray energies via free-free interaction and line emission from heavy elements. In this work we characterize the properties of the X-ray emission of the WHIM, and the LSS in general, focusing on the influence of different physical mechanisms, namely galactic winds (GWs), black-hole feedback and star-formation, and providing estimates of possible observational constraints. To this purpose we use a set of cosmological hydrodynamical simulations that include a self-consistent treatment of star-formation and chemical enrichment of the intergalactic medium, that allows us to follow the evolution of different metal species. We construct a set of simulated light-cones to make predictions of the emission in the 0.3-10 keV energy range. We obtain that GWs increase by a factor of 2 the emission of both galaxy clusters and WHIM. The amount of oxygen at average temperature and, consequently, the amount of expected bright Ovii and Oviii lines is increased by a factor of 3 due to GWs and by 20 per cent when assuming a top-heavy IMF. We compare our results with current observational constraints and find that the emission from faint groups and WHIM should account from half to all of the unresolved X-ray background in the 1-2 keV band.Comment: 15 pages, 8 figures, 4 tables. Accepted for publication in the MNRAS. Minor changes after referee repor

The XMM Cluster Outskirts Project (X-COP): Physical conditions to the virial radius of Abell 2142

Context. Galaxy clusters are continuously growing through the accretion of matter in their outskirts. This process induces inhomogeneities in the gas density distribution (clumping) which need to be taken into account to recover the physical properties of the intracluster medium (ICM) at large radii. Aims. We studied the thermodynamic properties in the outskirts (R > R500) of the massive galaxy cluster Abell 2142 by combining the Sunyaev Zel'dovich (SZ) effect with the X-ray signal. Methods. We combined the SZ pressure profile measured by Planck with the XMM-Newton gas density profile to recover radial profiles of temperature, entropy and hydrostatic mass out to 2R500. We used a method that is insensitive to clumping to recover the gas density, and we compared the results with traditional X-ray measurement techniques. Results. When taking clumping into account, our joint SZ/X-ray entropy profile is consistent with the predictions from pure gravitational collapse, whereas a significant entropy flattening is found when the effect of clumping is neglected. The hydrostatic mass profile recovered using joint X-ray/SZ data agrees with that obtained from spectroscopic X-ray measurements and with mass reconstructions obtained through weak lensing and galaxy kinematics. Conclusions. We found that clumping can explain the entropy flattening observed by Suzaku in the outskirts of several clusters. When using a method insensitive to clumping for the reconstruction of the gas density, the thermodynamic properties of Abell 2142 are compatible with the assumption that the thermal gas pressure sustains gravity and that the entropy is injected at accretion shocks, with no need to evoke more exotic physics. Our results highlight the need for X-ray observations with sufficient spatial resolution, and large collecting area, to understand the processes at work in cluster outer regions.Comment: 22 pages, 32 figures, accepted in the journal A&

The Sunyaev-Zel'dovich effects from a cosmological hydrodynamical simulation: large-scale properties and correlation with the soft X-ray signal

Using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model, we study the global properties of the Sunyaev-Zel'dovich (SZ) effects, both considering the thermal (tSZ) and the kinetic (kSZ) component. The simulation follows gravitation and gas dynamics and includes also several physical processes that affect the baryonic component, like a simple reionization scenario, radiative cooling, star formation and supernova feedback. Starting from the outputs of the simulation we create mock maps of the SZ signals due to the large structures of the Universe integrated in the range 0 < z < 6. We predict that the Compton y-parameter has an average value of (1.19 +/- 0.32) 10^-6 and is lognormally distributed in the sky; half of the whole signal comes from z < 1 and about 10 per cent from z > 2. The Doppler b-parameter shows approximately a normal distribution with vanishing mean value and a standard deviation of 1.6 10^-6, with a significant contribution from high-redshift (z > 3) gas. We find that the tSZ effect is expected to dominate the primary CMB anisotropies for l >~ 3000 in the Rayleigh-Jeans limit, while interestingly the kSZ effect dominates at all frequencies at very high multipoles (l >~ 7 10^4). We also analyse the cross-correlation between the two SZ effects and the soft (0.5-2 keV) X-ray emission from the intergalactic medium and we obtain a strong correlation between the three signals, especially between X-ray emission and tSZ effect (r_l ~ 0.8-0.9) at all angular scales.Comment: 12 pages, 15 figures. Accepted for publication in MNRAS. Minor changes, added reference

Searching for galaxy clusters in the Kilo-Degree Survey

In this paper, we present the tools used to search for galaxy clusters in the Kilo Degree Survey (KiDS), and our first results. The cluster detection is based on an implementation of the optimal filtering technique that enables us to identify clusters as over-densities in the distribution of galaxies using their positions on the sky, magnitudes, and photometric redshifts. The contamination and completeness of the cluster catalog are derived using mock catalogs based on the data themselves. The optimal signal to noise threshold for the cluster detection is obtained by randomizing the galaxy positions and selecting the value that produces a contamination of less than 20%. Starting from a subset of clusters detected with high significance at low redshifts, we shift them to higher redshifts to estimate the completeness as a function of redshift: the average completeness is ~ 85%. An estimate of the mass of the clusters is derived using the richness as a proxy. We obtained 1858 candidate clusters with redshift 0 < z_c < 0.7 and mass 13.5 < log(M500/Msun) < 15 in an area of 114 sq. degrees (KiDS ESO-DR2). A comparison with publicly available Sloan Digital Sky Survey (SDSS)-based cluster catalogs shows that we match more than 50% of the clusters (77% in the case of the redMaPPer catalog). We also cross-matched our cluster catalog with the Abell clusters, and clusters found by XMM and in the Planck-SZ survey; however, only a small number of them lie inside the KiDS area currently available.Comment: 13 pages, 15 figures. Accepted for publication on Astronomy & Astrophysic

Towards a Holistic View of the Heating and Cooling of the Intracluster Medium

(Abridged) X-ray clusters are conventionally divided into two classes: "cool core" (CC) clusters and "non-cool core" (NCC) clusters. Yet relatively little attention has been given to the origins of this dichotomy and, in particular, to the energetics and thermal histories of the two classes. We develop a model for the entropy profiles of clusters starting from the configuration established by gravitational shock heating and radiative cooling. At large radii, gravitational heating accounts for the observed profiles and their scalings well. However, at small and intermediate radii, radiative cooling and gravitational heating cannot be combined to explain the observed profiles of either type of cluster. The inferred entropy profiles of NCC clusters require that material is preheated prior to cluster collapse in order to explain the absence of low entropy (cool) material in these systems. We show that a similar modification is also required in CC clusters in order to match their properties at intermediate radii. In CC clusters, this modification is unstable, and an additional process is required to prevent cooling below a temperature of a few keV. We show that this can be achieved by adding a self-consistent AGN feedback loop in which the lowest-entropy, most rapidly cooling material is heated so that it rises buoyantly to mix with material at larger radii. The resulting model does not require fine tuning and is in excellent agreement with a wide variety of observational data. Some of the other implications of this model are briefly discussed.Comment: 27 pages, 13 figures, MNRAS accepted. Discussion of cluster heating energetics extended, results unchange

The Santa Fe Light Cone Simulation Project: I. Confusion and the WHIM in Upcoming Sunyaev-Zel'dovich Effect Surveys

We present the first results from a new generation of simulated large sky coverage (~100 square degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh refinement N-body/hydro code Enzo. We have simulated a very large (512^3h^{-3}Mpc^3) volume with unprecedented dynamic range. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. Unlike many previous studies of this type, our simulation includes unbound gas, where an appreciable fraction of the baryons in the universe reside. We have found that cluster line-of-sight overlap may be a significant issue in upcoming single-dish SZE surveys. Smaller beam surveys (~1 arcmin) have more than one massive cluster within a beam diameter 5-10% of the time, and a larger beam experiment like Planck has multiple clusters per beam 60% of the time. We explore the contribution of unresolved halos and unbound gas to the SZE signature at the maximum decrement. We find that there is a contribution from gas outside clusters of ~16% per object on average for upcoming surveys. This adds both bias and scatter to the deduced value of the integrated SZE, adding difficulty in accurately calibrating a cluster Y-M relationship. Finally, we find that in images where objects with M > 5x10^{13} M_{\odot} have had their SZE signatures removed, roughly a third of the total SZE flux still remains. This gas exists at least partially in the Warm Hot Intergalactic Medium (WHIM), and will possibly be detectable with the upcoming generation of SZE surveys.Comment: 14 pages, 13 figures, version accepted to ApJ. Major revisions mad