Systematics in the X-ray cluster mass estimators

We examine the systematics affecting the X-ray mass estimators applied to a set of five galaxy clusters resolved at high resolution in hydrodynamic simulations, including cooling, star formation and feedback processes. These simulated objects are processed through the X-ray Map Simulator, X-MAS, to provide Chandra-like long exposures that are analysed to reconstruct the gas temperature, density and mass profiles used as input. These clusters have different dynamic state: we consider a hot cluster with temperature T = 11.4 keV, a perturbed cluster with T = 3.9 keV, a merging object with T = 3.6 keV, and two relaxed systems with T = 3.3 keV and T = 2.7 keV, respectively. These systems are located at z = 0.175 so that their emission fits within the Chandra ACIS-S3 chip between 0.6 and 1.2 R-500. We find that the mass profile obtained via a direct application of the hydrostatic equilibrium ( HE) equation is dependent upon the measured temperature profile. An irregular radial distribution of the temperature values, with associated large errors, induces a significant scatter on the reconstructed mass measurements. At R-2500, the actual mass is recovered within 1 sigma, although we notice this estimator shows high statistical errors due to high level of Chandra background. Instead, the poorness of the beta-model in describing the gas density profile makes the evaluated masses to be underestimated by similar to 40 per cent with respect to the true mass, both with an isothermal and a polytropic temperature profile. We also test ways to recover the mass by adopting an analytic mass model, such as those proposed by Nvarro, Frenk & White and Rasia, Tormen & Moscardini, and fitting the temperature profile expected from the HE equation to the observed one. We conclude that the methods of the HE equation and those of the analytic fits provide a more robust mass estimation than the ones based on the beta-model. In the present work, the main limitation for a precise mass reconstruction is to ascribe to the relatively high level of the background chosen to reproduce the Chandra one. After artificially reducing the total background by a factor of 100, we find that the estimated mass significantly underestimates the true mass profiles. This is manly due (i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to the bias towards lower values of the X-ray temperature measurements because of the complex thermal structure of the emitting plasma.