Properties of gas clumps and gas clumping factor in the intra cluster medium

The spatial distribution of gas matter inside galaxy clusters is not completely smooth, but may host gas clumps associated with substructures. These overdense gas substructures are generally a source of unresolved bias of X-ray observations towards high density gas, but their bright luminosity peaks may be resolved sources within the ICM, that deep X-ray exposures may be (already) capable to detect. In this paper we aim at investigating both features, using a set of high-resolution cosmological simulations with ENZO. First, we monitor how the bias by unresolved gas clumping may yield incorrect estimates of global cluster parameters and affects the measurements of baryon fractions by X-ray observations. We find that based on X-ray observations of narrow radial strips, it is difficult to recover the real baryon fraction to better than 10 - 20 percent uncertainty. Second, we investigated the possibility of observing bright X-ray clumps in the nearby Universe (z<=0.3). We produced simple mock X-ray observations for several instruments (XMM, Suzaku and ROSAT) and extracted the statistics of potentially detectable bright clumps. Some of the brightest clumps predicted by simulations may already have been already detected in X- ray images with a large field of view. However, their small projected size makes it difficult to prove their existence based on X-ray morphology only. Preheating, AGN feedback and cosmic rays are found to have little impact on the statistical properties of gas clumps.Comment: 17 pages, 11 figures. MNRAS accepte

The cluster gas mass fraction as a cosmological probe: a revised study

(Abriged) We present the analysis of the baryonic content of 52 X-ray luminous galaxy clusters observed with Chandra in the redshift range 0.3-1.273. We use the deprojected X-ray surface brightness profiles and the measured values of the gas temperature to recover the gas and total mass profiles. By assuming that galaxy clusters are representative of the cosmic baryon budget, the distribution of the cluster baryon fraction in the hottest (T> 4 keV) systems as a function of redshift is used to constrain the cosmological parameters. We discuss how our constraints are affected by several systematics, namely the isothermality, the assumed baryon fraction in stars, the depletion parameter and the sample selection. By using only the cluster baryon fraction as a proxy for the cosmological parameters, we obtain that Omega is very well constrained at the value of 0.35 with a relative statistical uncertainty of 11% (1 sigma level; w=-1) and a further systematic error of about (-6,+7)%. On the other hand, constraints on Lambda (without the prior of flat geometry) and w (using the prior of flat geometry) are definitely weaker due to the presence of larger statistical and systematic uncertainties (of the order of 40 per cent on Lambda and larger than 50 per cent on w). If the WMAP 5-year best-fit results are assumed to fix the cosmological parameters, we limit the contributions expected from non-thermal pressure support and ICM clumpiness to be lower than about 10 per cent, leaving also room to accommodate baryons not accounted for either in the X-ray emitting plasma or in stars of the order of 18 per cent of the total cluster baryon budget.Comment: A&A in press. Accepted on March 28, 2009. Revised to match version in prin

On the evolution of cooling cores in X-ray galaxy clusters

(Abridged) To define a framework for the formation and evolution of the cooling cores in X-ray galaxy clusters, we study how the physical properties change as function of the cosmic time in the inner regions of a 4 keV and 8 keV galaxy cluster under the action of radiative cooling and gravity only. The cooling radius, R_cool, defined as the radius at which the cooling time equals the Universe age at given redshift, evolves from ~0.01 R200 at z>2, where the structures begin their evolution, to ~0.05 R200 at z=0. The values measured at 0.01 R200 show an increase of about 15-20 per cent per Gyr in the gas density and surface brightness and a decrease with a mean rate of 10 per cent per Gyr in the gas temperature. The emission-weighted temperature diminishes by about 25 per cent and the bolometric X-ray luminosity rises by a factor ~2 after 10 Gyrs when all the cluster emission is considered in the computation. On the contrary, when the core region within 0.15 R500 is excluded, the gas temperature value does not change and the X-ray luminosity varies by 10-20 per cent only. The cooling time and gas entropy radial profiles are well represented by power-law functions. The behaviour of the inner slopes of the gas temperature and density profiles are the most sensitive and unambiguous tracers of an evolving cooling core. Their values after 10 Gyrs of radiative losses, T_gas ~ r^0.4 and n_gas ~ r^(-1.2) for the hot (cool) object, are remarkably in agreement with the observational constraints available for nearby X-ray luminous cooling core clusters. Because our simulations do not consider any AGN heating, they imply that the feedback process does not greatly alter the gas density and temperature profiles as generated by radiative cooling alone.Comment: 8 pages. MNRAS in pres

The evolution of the spatially-resolved metal abundance in galaxy clusters up to z=1.4

We present the combined analysis of the metal content of 83 objects in the redshift range 0.09-1.39, and spatially-resolved in the 3 bins (0-0.15, 0.15-0.4, >0.4) R500, as obtained with similar analysis using XMM-Newton data in Leccardi & Molendi (2008) and Baldi et al. (2012). We use the pseudo-entropy ratio to separate the Cool-Core (CC) cluster population, where the central gas density tends to be relatively higher, cooler and more metal rich, from the Non-Cool-Core systems. The average, redshift-independent, metal abundance measured in the 3 radial bins decrease moving outwards, with a mean metallicity in the core that is even 3 (two) times higher than the value of 0.16 times the solar abundance in Anders & Grevesse (1989) estimated at r>0.4 R500 in CC (NCC) objects. We find that the values of the emission-weighted metallicity are well-fitted by the relation $Z(z) = Z_0 (1+z)^{-\gamma}$ at given radius. A significant scatter, intrinsic to the observed distribution and of the order of 0.05-0.15, is observed below 0.4 R500. The nominal best-fit value of $\gamma$ is significantly different from zero in the inner cluster regions ($\gamma = 1.6 \pm 0.2$) and in CC clusters only. These results are confirmed also with a bootstrap analysis, which provides a still significant negative evolution in the core of CC systems (P>99.9 per cent). No redshift-evolution is observed when regions above the core (r > 0.15 R500) are considered. A reasonable good fit of both the radial and redshift dependence is provided from the functional form $Z(r,z)=Z_0 (1+(r/0.15 R500)^2)^{-\beta} (1+z)^{-\gamma}$, with $(Z_0, \beta, \gamma) = (0.83 \pm 0.13, 0.55 \pm 0.07, 1.7 \pm 0.6)$ in CC clusters and $(0.39 \pm 0.04, 0.37 \pm 0.15, 0.5 \pm 0.5)$ for NCC systems. Our results represent the most extensive study of the spatially-resolved metal distribution in the cluster plasma as function of redshift.Comment: 5 pages. Research Note accepted for publication in A&

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

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

XMM-Newton observation of the interacting cluster Abell 3528

We analyze the XMM dataset of the interacting cluster of galaxies Abell 3528 located westward in the core of the Shapley Supercluster, the largest concentration of mass in the nearby Universe. A3528 is formed by two interacting clumps (A3528-N at North and A3528-S at South) separated by 0.9 $h_{70}^{-1}$ Mpc at redshift 0.053. XMM data describe these clumps as relaxed structure with an overall temperature of $4.14 \pm 0.09$ and $4.29 \pm 0.07$ keV in A3528-N and A3528-S, respectively, and a core cooler by a factor 1.4--1.5 and super-solar metal abundance in the inner 30 arcsec. These clumps are connected by a X-ray soft, bridge-like emission and present asymmetric surface brightness with significant excess in the North--West region of A3528-N and in the North--East area of A3528-S. However, we do not observe any evidence of shock heated gas, both in the surface brightness and in the temperature map. Considering also that the optical light distribution is more concentrated around A3528-N and makes A3528-S barely detectable, we do not find support to the originally suggested head-on pre-merging scenario and conclude that A3528 is in a off-axis post-merging phase, where the closest cores encounter happened about 1--2 Gyrs ago.Comment: 15 pages, 16 figures (10 in colors), accepted for publication on A&A. For a Poscript version with high resolution figures see http://www.mi.iasf.cnr.it/~gasta/www/a3528.htm

Mass profiles and concentration-dark matter relation in X-ray luminous galaxy clusters

(Abriged) Assuming that the hydrostatic equilibrium holds between the intracluster medium and the gravitational potential, we constrain the NFW profiles in a sample of 44 X-ray luminous galaxy clusters observed with XMM-Newton in the redshift range 0.1-0.3. We evaluate several systematic uncertainties that affect our reconstruction of the X-ray masses. We measure the concentration c200, the dark mass M200 and the gas mass fraction within R500 in all the objects of our sample, providing the largest dataset of mass parameters for galaxy clusters in this redshift range. We confirm that a tight correlation between c200 and M200 is present and in good agreement with the predictions from numerical simulations and previous observations. When we consider a subsample of relaxed clusters that host a Low-Entropy-Core (LEC), we measure a flatter c-M relation with a total scatter that is lower by 40 per cent. From the distribution of the estimates of c200 and M200, with associated statistical (15-25%) and systematic (5-15%) errors, we use the predicted values from semi-analytic prescriptions calibrated through N-body numerical runs and measure sigma_8*Omega_m^(0.60+-0.03)= 0.45+-0.01 (at 2 sigma level, statistical only) for the subsample of the clusters where the mass reconstruction has been obtained more robustly, and sigma_8*Omega_m^(0.56+-0.04) = 0.39+-0.02 for the subsample of the 11 more relaxed LEC objects. With the further constraint from the fgas distribution in our sample, we break the degeneracy in the sigma_8-Omega_m plane and obtain the best-fit values sigma_8~1.0+-0.2 (0.75+-0.18 when the subsample of the more relaxed objects is considered) and Omega_m = 0.26+-0.01.Comment: 21 pages. A&A in press. Minor revisions to match accepted version. Corrected 2nd and 3rd column in Table 3, and equation (A.4

On the connection between turbulent motions and particle acceleration in galaxy clusters

Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected throughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{\rm radio}\propto\sigma_v^{3.3\pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.Comment: Submitted to ApJ Letter

A Chandra archival study of the temperature and metal abundance profiles in hot Galaxy Clusters at 0.1 < z < 0.3

We present the analysis of the temperature and metallicity profiles of 12 galaxy clusters in the redshift range 0.1--0.3 selected from the Chandra archive with at least ~20,000 net ACIS counts and kT>6 keV. We divide the sample between 7 Cooling-Core (CC) and 5 Non-Cooling-Core (NCC) clusters according to their central cooling time. We find that single power-laws can describe properly both the temperature and metallicity profiles at radii larger than 0.1 r_180 in both CC and NCC systems, showing the NCC objects steeper profiles outwards. A significant deviation is only present in the inner 0.1 r_180. We perform a comparison of our sample with the De Grandi & Molendi BeppoSAX sample of local CC and NCC clusters, finding a complete agreement in the CC cluster profile and a marginally higher value (at ~1sigma) in the inner regions of the NCC clusters. The slope of the power-law describing kT(r) within 0.1 r_180 correlates strongly with the ratio between the cooling time and the age of the Universe at the cluster redshift, being the slope >0 and tau_c/tau_age<=0.6 in CC systems.Comment: 12 pages, 6 figures, Accepted for publication by the Astrophysical Journa