Testing Sunyaev-Zel'dovich measurements of the hot gas content of dark matter haloes using synthetic skies

The thermal Sunyaev-Zel'dovich (tSZ) effect offers a means of probing the hot gas in and around massive galaxies and galaxy groups and clusters, which is thought to constitute a large fraction of the baryon content of the Universe. The Planck collaboration recently performed a stacking analysis of a large sample of `locally brightest galaxies' (LBGs) and, surprisingly, inferred an approximately self-similar relation between the tSZ flux and halo mass. At face value, this implies that the hot gas mass fraction is independent of halo mass, a result which is in apparent conflict with resolved X-ray observations. We test the robustness of the inferred trend using synthetic tSZ maps generated from cosmological hydrodynamical simulations and using the same tools and assumptions applied in the Planck study. We show that, while the detection and the estimate of the `total' flux (within $5 r_{500}$) is reasonably robust, the inferred flux originating from within $r_{500}$ (i.e. the limiting radius to which X-ray observations typically probe) is highly sensitive to the assumed pressure distribution of the gas. Using our most realistic simulations with AGN feedback, that reproduce a wide variety of X-ray and optical properties of groups and clusters, we estimate that the derived tSZ flux within $r_{500}$ is biased high by up to to an order of magnitude for haloes with masses $M_{500} \sim 10^{13}$ M$_{\odot}$. Moreover, we show that the AGN simulations are consistent with the total tSZ flux-mass relation observed with Planck, whereas a self-similar model is ruled out.Comment: 16 pages, 10 figures, MNRAS, accepted after minor revisio

Non-thermal X-ray Emission: An Alternative to Cluster Cooling Flows?

We report the results of experiments aimed at reducing the major problem with cooling flow models of rich cluster X-ray sources: the fact that most of the cooled gas or its products have not been found. Here we show that much of the X-ray emission usually attributed to cooling flows can, in fact, be modeled by a power-law component which is indicative of a source(s) other than thermal bremsstrahlung from the intracluster medium. We find that adequate simultaneous fits to ROSAT PSPCB and ASCA GIS/SIS spectra of the central regions of ten clusters are obtained for two-component models that include a thermal plasma component that is attributable to hot intracluster gas and a power-law component that is likely generated by compact sources and/or extended non-thermal emission. For five of the clusters that purportedly have massive cooling flows, the best-fit models have power-law components that contribute $\sim$ 30 % of the total flux (0.14 - 10.0 keV) within the central 3 arcminutes. Because cooling flow mass deposition rates are inferred from X-ray fluxes, our finding opens the possibility of significantly reducing cooling rates.Comment: 11 pages, 3 figures, emulateapj style. Accepted for publication in Ap

The scatter and evolution of the global hot gas properties of simulated galaxy cluster populations

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations to investigate the scatter and evolution of the global hot gas properties of large simulated populations of galaxy groups and clusters. Our aim is to compare the predictions of different physical models and to explore the extent to which commonly-adopted assumptions in observational analyses (e.g. self-similar evolution) are violated. We examine the relations between (true) halo mass and the X-ray temperature, X-ray luminosity, gas mass, Sunyaev-Zel'dovich (SZ) flux, the X-ray analogue of the SZ flux ($Y_X$) and the hydrostatic mass. For the most realistic models, which include AGN feedback, the slopes of the various mass-observable relations deviate substantially from the self-similar ones, particularly at late times and for low-mass clusters. The amplitude of the mass-temperature relation shows negative evolution with respect to the self-similar prediction (i.e. slower than the prediction) for all models, driven by an increase in non-thermal pressure support at higher redshifts. The AGN models predict strong positive evolution of the gas mass fractions at low halo masses. The SZ flux and $Y_X$ show positive evolution with respect to self-similarity at low mass but negative evolution at high mass. The scatter about the relations is well approximated by log-normal distributions, with widths that depend mildly on halo mass. The scatter decreases significantly with increasing redshift. The exception is the hydrostatic mass-halo mass relation, for which the scatter increases with redshift. Finally, we discuss the relative merits of various hot gas-based mass proxies.Comment: 31 pages (21 before appendices), 19 figures, 12 tables, accepted by MNRAS after minor revisio

Models of the ICM with Heating and Cooling: Explaining the Global and Structural X-ray Properties of Clusters

(Abridged) Theoretical models that include only gravitationally-driven processes fail to match the observed mean X-ray properties of clusters. As a result, there has recently been increased interest in models in which either radiative cooling or entropy injection play a central role in mediating the properties of the intracluster medium. Both sets of models give reasonable fits to the mean properties of clusters, but cooling only models result in fractions of cold baryons in excess of observationally established limits and the simplest entropy injection models do not treat the "cooling core" structure present in many clusters and cannot account for entropy profiles revealed by recent X-ray observations. We consider models that marry radiative cooling with entropy injection, and confront model predictions for the global and structural properties of massive clusters with the latest X-ray data. The models successfully and simultaneously reproduce the observed L-T and L-M relations, yield detailed entropy, surface brightness, and temperature profiles in excellent agreement with observations, and predict a cooled gas fraction that is consistent with observational constraints. The model also provides a possible explanation for the significant intrinsic scatter present in the L-T and L-M relations and provides a natural way of distinguishing between clusters classically identified as "cooling flow" clusters and dynamically relaxed "non-cooling flow" clusters. The former correspond to systems that had only mild levels (< 300 keV cm^2) of entropy injection, while the latter are identified as systems that had much higher entropy injection. This is borne out by the entropy profiles derived from Chandra and XMM-Newton.Comment: 20 pages, 15 figures, accepted for publication in the Astrophysical Journa

The impact of baryonic processes on the two-point correlation functions of galaxies, subhaloes and matter

The observed clustering of galaxies and the cross-correlation of galaxies and mass provide important constraints on both cosmology and models of galaxy formation. Even though the dissipation and feedback processes associated with galaxy formation are thought to affect the distribution of matter, essentially all models used to predict clustering data are based on collisionless simulations. Here, we use large hydrodynamical simulations to investigate how galaxy formation affects the autocorrelation functions of galaxies and subhaloes, as well as their cross-correlation with matter. We show that the changes due to the inclusion of baryons are not limited to small scales and are even present in samples selected by subhalo mass. Samples selected by subhalo mass cluster ~10% more strongly in a baryonic run on scales r > 1Mpc/h, and this difference increases for smaller separations. While the inclusion of baryons boosts the clustering at fixed subhalo mass on all scales, the sign of the effect on the cross-correlation of subhaloes with matter can vary with radius. We show that the large-scale effects are due to the change in subhalo mass caused by the strong feedback associated with galaxy formation and may therefore not affect samples selected by number density. However, on scales r < r_vir significant differences remain after accounting for the change in subhalo mass. We conclude that predictions for galaxy-galaxy and galaxy-mass clustering from models based on collisionless simulations will have errors greater than 10% on sub-Mpc scales, unless the simulation results are modified to correctly account for the effects of baryons on the distributions of mass and satellites.Comment: 15 pages, 9 figures. Replaced to match the version accepted by MNRA