Parametrization effects in the analysis of AMI Sunyaev-Zel'dovich observations

• Published in: Monthly notices of the Royal Astronomical Society 2012-04, Vol.421 (2), p.1136-1154
• Main Authors: Olamaie, Malak, Rodríguez-Gonzálvez, Carmen, Davies, Matthew L., Feroz, Farhan, Franzen, Thomas M. O., Grainge, Keith J. B., Hobson, Michael P., Hurley-Walker, Natasha, Lasenby, Anthony N., Pooley, Guy G., Saunders, Richard D. E., Scaife, Anna M. M., Schammel, Michel, Scott, Paul F., Shimwell, Timothy W., Titterington, David J., Waldram, Elizabeth M., Zwart, Jonathan T. L.
• Format: Article
• Language: English
• Subjects: Astrophysics; cosmic background radiation; Cosmology; cosmology: observations; galaxies: clusters: general; methods: data analysis; Simulation; Stars & galaxies
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1111/j.1365-2966.2011.20374.x

Most Sunyaev-Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the cluster total mass (M T(r)) and/or gas mass (M g(r)) using parametrized models derived from both simulations and imaging observations, and assumptions of spherical symmetry and hydrostatic equilibrium. By numerically exploring the probability distributions of the cluster parameters given the simulated interferometric SZ data in the context of Bayesian methods, and assuming a β-model for the electron number density n e(r) described by two shape parameters β and r c, we investigate the capability of this model and analysis to return the simulated cluster input quantities via three parametrizations. In parametrization I we assume that the gas temperature is an independent free parameter and assume hydrostatic equilibrium, spherical geometry and an ideal gas equation of state. We find that parametrization I can hardly constrain the cluster parameters and fails to recover the true values of the simulated cluster. In particular it overestimates M T(r 200) and T g(r 200) (M T(r 200) = (6.43 ± 5.43) × 1015 M⊙ and T g(r 200) = (10.61 ± 5.28) keV) compared to the corresponding values of the simulated cluster (M T(r 200) = 5.83 × 1014 M⊙ and T g(r 200) = 5 keV). We then investigate parametrizations II and III in which f g(r 200) replaces temperature as a main variable; we do this because f g may vary significantly less from cluster to cluster than temperature. In parametrization II we relate M T(r 200) and T g assuming hydrostatic equilibrium. We find that parametrization II can constrain the cluster physical parameters but the temperature estimate is biased low (M T(r 200) = (6.8 ± 2.1) × 1014 M⊙ and T g(r 200) = (3.0 ± 1.2) keV). In parametrization III, the virial theorem (plus the assumption that all the kinetic energy of the cluster is the internal energy of the gas) replaces the hydrostatic equilibrium assumption because we consider it more robust both in theory and in practice. We find that parametrization III results in unbiased estimates of the cluster properties (M T(r 200) = (4.68 ± 1.56) × 1014 M⊙ and T g(r 200) = (4.3 ± 0.9) keV). We generate a second simulated cluster using a generalized Navarro-Frenk-White pressure profile and analyse it with an entropy-based model to take into account the temperature gradient in our analysis and improve the cluster gas density distribution. This model also constrains the cluster physical parameters and the results show a radial decl