A New Robust Low-Scatter X-ray Mass Indicator for Clusters of Galaxies

We present comparison of X-ray proxies for the total cluster mass, including the spectral temperature (Tx), gas mass measured within r500 (Mg), and the new proxy, Yx, which is a simple product of Tx and Mg and is related to the total thermal energy of the ICM. We use mock Chandra images constructed for a sample of clusters simulated with the eulerian N-body+gasdynamics adaptive mesh refinement ART code in the concordance LCDM cosmology. The simulations achieve high spatial and mass resolution and include radiative cooling, star formation, and other processes accompanying galaxy formation. Our analysis shows that simulated clusters exhibit a high degree of regularity and tight correlations between the considered observables and total mass. The normalizations of the M-Tx, Mg-Tx, and M-Yx relations agree to better than 10-15% with the current observational measurements of these relations. Our results show that Yx is the best mass proxy with a remarkably low scatter of only ~5-7% in M500 for a fixed Yx, at both low and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx-M500 relation is close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the cluster mass function from X-ray selected samples.Comment: submitted to ApJ; 9 pages, 6 figures, uses emulateap

The Radial Distribution of Galaxies in LCDM clusters

We study the radial distribution of subhalos and galaxies using high-resolution cosmological simulations of galaxy clusters formed in the concordance LCDM cosmology. In agreement with previous studies, we find that the radial distribution of subhalos is significantly less concentrated than that of the dark matter, when subhalos are selected using their present-day gravitationally bound mass. We show that the difference in the radial distribution is not a numerical artifact and is due to tidal stripping. The subhalos in the cluster core lose more than 70% of their initial mass since accretion, while the average tidal mass loss for halos near the virial radius is ~30%. This introduces a radial bias in the spatial distribution of subhalos when they are selected using their tidally truncated mass. We demonstrate that the radial bias disappears almost entirely if subhalos are selected using their mass or circular velocity at the accretion epoch. The comparisons of the results of dissipationless simulations to the observed distribution of galaxies in clusters are therefore sensitive to the selection criteria used to select subhalo samples. Using the simulations that include cooling and starformation, we show that the radial distribution of subhalos is in reasonable agreement with the observed radial distribution of galaxies in clusters for 0.1<R/R200<2.0, if subhalos are selected using the stellar mass of galaxies. The radial bias is minimized in this case because the stars are located in the centers of dark matter subhalos and are tightly bound. The stellar mass of an object is therefore approximately conserved as the dark matter is stripped from the outer regions. Nevertheless, the concentration of the radial distribution of galaxies is systematically lower than that of the dark matter.Comment: submitted to ApJ, 12 pages, 12 figure

Simulating the Formation of Galaxy Clusters

We study the effects of radiative cooling, star formation and stellar feedback on the properties and evolution of galaxy clusters using high-resolution Adaptive Mesh Refinement N-body+gasdynamics simulations of clusters forming in the LCDM universe. Cooling leads to the condensation of gas in the inner regions of clusters, which in turn leads to steepening of the dark matter profile. The cooling gas is replaced by the higher-entropy gas from the outer regions, which raises the entropy and temperature of gas in the cluster core. The magnitude of these effects is likely overestimated in the current simulations because they suffer from the overcooling problem: a much larger fraction of baryons is in the form of cold gas and stars than is observed. We find that the thermal stellar feedback alone does not remedy this problem. Additional ad-hoc preheating can lower the amount of cold gas but a simple uniform preheating results in incorrect star formation history, as it delays the bulk of star formation until z<1. Our analysis shows that the overcooling in a cluster as a whole is really the overcooling in the central galaxy and its progenitors at high redshifts. This indicates that an additional heating mechanism that can continuously heat the gas in the cluster core is required to reproduce the observed cluster properties. Energy injection by the Active Galactic Nuclei, which may provide such heating, may thus be an important missing ingredient in the current theoretical models of cluster formation.Comment: 3 pages, 2 figures. Use iaus style. To appear in the proceedings of IAU Colloquium 195 - "Outskirts of galaxy clusters: intense life in the suburbs", Torino, Italy, March 12-16, 200

Testing X-ray Measurements of Galaxy Clusters with Cosmological Simulations

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this paper, we present mock Chandra analyses of cosmological cluster simulations and assess X-ray measurements of galaxy cluster properties using a model and procedure essentially identical to that used in real data analysis. We show that reconstruction of three-dimensional ICM density and temperature profiles is excellent for relaxed clusters, but still reasonably accurate for unrelaxed systems. The total ICM mass is measured quite accurately (<6%) in all clusters, while the hydrostatic estimate of the gravitationally bound mass is biased low by about 5%-20% through the virial region, primarily due to additional pressure support provided by subsonic bulk motions in the ICM, ubiquitous in our simulations even in relaxed systems. Gas fraction determinations are therefore biased high; the bias increases toward cluster outskirts and depends sensitively on its dynamical state, but we do not observe significant trends of the bias with cluster mass or redshift. We also find that different average ICM temperatures, such as the X-ray spectroscopic Tspec and gas-mass-weighted Tmg, are related to each other by a constant factor with a relatively small object-to-object scatter and no systematic trend with mass, redshift or the dynamical state of clusters. We briefly discuss direct applications of our results for different cluster-based cosmological tests.Comment: 11 pages, 6 figures, submitted to Ap

Shapes of Gas, Gravitational Potential and Dark Matter in Lambda-CDM Clusters

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological simulations of the Lambda-CDM cosmology and compare it to the shape of dark matter distribution and the shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation and stellar feedback (CSF), intracluster gas outside the cluster core is more spherical compared to non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in NR simulations the difference between gas and potential shape at small radii is due to random gas motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters against observations can therefore constrain the amount of cooling of the intracluster medium and the presence of random gas motions in cluster cores.Comment: 11 pages, 8 figures, 3 tables, updated to match the version accepted for publication in the Astrophysical Journa

Effects of Galaxy Formation on Thermodynamics of the Intracluster Medium

We present detailed comparisons of the intracluster medium (ICM) in cosmological Eulerian cluster simulations with deep Chandra observations of nearby relaxed clusters. To assess the impact of galaxy formation, we compare two sets of simulations, one performed in the non-radiative regime and another with radiative cooling and several physical processes critical to various aspects of galaxy formation: star formation, metal enrichment and stellar feedback. We show that the observed ICM properties outside cluster cores are well-reproduced in the simulations that include cooling and star formation, while the non-radiative simulations predict an overall shape of the ICM profiles inconsistent with observations. In particular, we find that the ICM entropy in our runs with cooling is enhanced to the observed levels at radii as large as half of the virial radius. We also find that outside cluster cores entropy scaling with the mean ICM temperature in both simulations and Chandra observations is consistent with being self-similar within current error bars. We find that the pressure profiles of simulated clusters are also close to self-similar and exhibit little cluster-to-cluster scatter. The X-ray observable-total mass relations for our simulated sample agree with the Chandra measurements to \~10%-20% in normalization. We show that this systematic difference could be caused by the subsonic gas motions, unaccounted for in X-ray hydrostatic mass estimates. The much improved agreement of simulations and observations in the ICM profiles and scaling relations is encouraging and the existence of tight relations of X-ray observables, such as Yx, and total cluster mass and the simple redshift evolution of these relations hold promise for the use of clusters as cosmological probes.Comment: 14 pages, 6 figures. Matches version accepted to Ap