This is the first in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here we present a new, automated method for identifying relaxed clusters based on their morphologies in X-ray imaging dat
a. While broadly similar to others in the literature, the morphological quantities that we measure are specifically designed to provide a fair basis for comparison across a range of data quality and cluster redshifts, to be robust against missing data due to point-source masks and gaps between detectors, and to avoid strong assumptions about the cosmological background and cluster masses. Based on three morphological indicators - Symmetry, Peakiness and Alignment - we develop the SPA criterion for relaxation. This analysis was applied to a large sample of cluster observations from the Chandra and ROSAT archives. Of the 361 clusters which received the SPA treatment, 57 (16 per cent) were subsequently found to be relaxed according to our criterion. We compare our measurements to similar estimators in the literature, as well as projected ellipticity and other image measures, and comment on trends in the relaxed cluster fraction with redshift, temperature, and survey selection method. Code implementing our morphological analysis will be made available on the web.
We employ robust weak gravitational lensing measurements to improve cosmological constraints from measurements of the galaxy cluster mass function and its evolution, using X-ray selected clusters detected in the ROSAT All-Sky Survey. Our lensing anal
ysis constrains the absolute mass scale of such clusters at the 8 per cent level, including both statistical and systematic uncertainties. Combining it with the survey data and X-ray follow-up observations, we find a tight constraint on a combination of the mean matter density and late-time normalization of the matter power spectrum, $sigma_8(Omega_m/0.3)^{0.17}=0.81pm0.03$, with marginalized, one-dimensional constraints of $Omega_m=0.26pm0.03$ and $sigma_8=0.83pm0.04$. For these two parameters, this represents a factor of two improvement in precision with respect to previous work, primarily due to the reduced systematic uncertainty in the absolute mass calibration provided by the lensing analysis. Our new results are in good agreement with constraints from cosmic microwave background (CMB) data, both WMAP and Planck (plus WMAP polarization), under the assumption of a flat $Lambda$CDM cosmology with minimal neutrino mass. Consequently, we find no evidence for non-minimal neutrino mass from the combination of cluster data with CMB, supernova and baryon acoustic oscillation measurements, regardless of which all-sky CMB data set is used (and independent of the recent claimed detection of B-modes on degree scales). We also present improved constraints on models of dark energy (both constant and evolving), modifications of gravity, and primordial non-Gaussianity. Assuming flatness, the constraints for a constant dark energy equation of state from the cluster data alone are at the 15 per cent level, improving to $sim 6$ per cent when the cluster data are combined with other leading probes.
We present cosmological constraints from measurements of the gas mass fraction, $f_{gas}$, for massive, dynamically relaxed galaxy clusters. Our data set consists of Chandra observations of 40 such clusters, identified in a comprehensive search of th
e Chandra archive, as well as high-quality weak gravitational lensing data for a subset of these clusters. Incorporating a robust gravitational lensing calibration of the X-ray mass estimates, and restricting our measurements to the most self-similar and accurately measured regions of clusters, significantly reduces systematic uncertainties compared to previous work. Our data for the first time constrain the intrinsic scatter in $f_{gas}$, $(7.4pm2.3)$% in a spherical shell at radii 0.8-1.2 $r_{2500}$, consistent with the expected variation in gas depletion and non-thermal pressure for relaxed clusters. From the lowest-redshift data in our sample we obtain a constraint on a combination of the Hubble parameter and cosmic baryon fraction, $h^{3/2}Omega_b/Omega_m=0.089pm0.012$, that is insensitive to the nature of dark energy. Combined with standard priors on $h$ and $Omega_b h^2$, this provides a tight constraint on the cosmic matter density, $Omega_m=0.27pm0.04$, which is similarly insensitive to dark energy. Using the entire cluster sample, extending to $z>1$, we obtain consistent results for $Omega_m$ and interesting constraints on dark energy: $Omega_Lambda=0.65^{+0.17}_{-0.22}$ for non-flat $Lambda$CDM models, and $w=-0.98pm0.26$ for flat constant-$w$ models. Our results are both competitive and consistent with those from recent CMB, SNIa and BAO data. We present constraints on models of evolving dark energy from the combination of $f_{gas}$ data with these external data sets, and comment on the possibilities for improved $f_{gas}$ constraints using current and next-generation X-ray observatories and lensing data. (Abridged)
(Abridged) This is the second in a series of papers in which we derive simultaneous constraints on cosmology and X-ray scaling relations using observations of massive, X-ray flux-selected galaxy clusters. The data set consists of 238 clusters drawn f
rom the ROSAT All-Sky Survey with 0.1-2.4 keV luminosities >2.5e44 erg/second, and incorporates extensive follow-up observations using the Chandra X-ray Observatory. Our analysis accounts self-consistently for all selection effects, covariances and systematic uncertainties. Here we describe the reduction of the follow-up X-ray observations, present results on the cluster scaling relations, and discuss their implications. Our constraints on the luminosity-mass and temperature-mass relations, measured within r_500, lead to three important results. First, the data support the conclusion that excess heating of the intracluster medium has altered its thermodynamic state from that expected in a simple, gravitationally dominated system; however, this excess heating is primarily limited to the central regions of clusters (r<0.15r_500). Second, the intrinsic scatter in the center-excised luminosity-mass relation is remarkably small, being undetected at the <10% level in current data; for the hot, massive clusters under investigation, this scatter is smaller than in either the temperature-mass or Y_X-mass relations (10-15%). Third, the evolution with redshift of the scaling relations is consistent with the predictions of simple, self-similar models of gravitational collapse, indicating that the mechanism responsible for heating the central regions of clusters was in operation before redshift 0.5 (the limit of our data) and that its effects on global cluster properties have not evolved strongly since then.
(Abridged) This is the first of a series of papers in which we derive simultaneous constraints on cosmological parameters and X-ray scaling relations using observations of the growth of massive, X-ray flux-selected galaxy clusters. Our data set consi
sts of 238 clusters drawn from the ROSAT All-Sky Survey, and incorporates extensive follow-up observations using the Chandra X-ray Observatory. Here we describe and implement a new statistical framework required to self-consistently produce simultaneous constraints on cosmology and scaling relations from such data, and present results on models of dark energy. In spatially flat models with a constant dark energy equation of state, w, the cluster data yield Omega_m=0.23 +- 0.04, sigma_8=0.82 +- 0.05, and w=-1.01 +- 0.20, marginalizing over conservative allowances for systematic uncertainties. These constraints agree well and are competitive with independent data in the form of cosmic microwave background (CMB) anisotropies, type Ia supernovae (SNIa), cluster gas mass fractions (fgas), baryon acoustic oscillations (BAO), galaxy redshift surveys, and cosmic shear. The combination of our data with current CMB, SNIa, fgas, and BAO data yields Omega_m=0.27 +- 0.02, sigma_8=0.79 +- 0.03, and w=-0.96 +- 0.06 for flat, constant w models. For evolving w models, marginalizing over transition redshifts in the range 0.05-1, we constrain the equation of state at late and early times to be respectively w_0=-0.88 +- 0.21 and w_et=-1.05 +0.20 -0.36. The combined data provide constraints equivalent to a DETF FoM of 15.5. Our results highlight the power of X-ray studies to constrain cosmology. However, the new statistical framework we apply to this task is equally applicable to cluster studies at other wavelengths.
We present constraints on the mean matter density, Omega_m, the normalization of the density fluctuation power spectrum, sigma_8, and the dark-energy equation-of-state parameter, w, obtained from measurements of the X-ray luminosity function of the l
argest known galaxy clusters at redshifts z<0.7, as compiled in the Massive Cluster Survey (MACS) and the local BCS and REFLEX galaxy cluster samples. Our analysis employs an observed mass-luminosity relation, calibrated by hydrodynamical simulations, including corrections for non-thermal pressure support and accounting for the presence of intrinsic scatter. Conservative allowances for all known systematic uncertainties are included, as are standard priors on the Hubble constant and mean baryon density. We find Omega_m=0.28 +0.11 -0.07 and sigma_8=0.78 +0.11 -0.13 for a spatially flat, cosmological-constant model, and Omega_m=0.24 +0.15 -0.07, sigma_8=0.85 +0.13 -0.20 and w=-1.4 +0.4 -0.7 for a flat, constant-w model. Future work improving our understanding of redshift evolution and observational biases affecting the mass--X-ray luminosity relation have the potential to significantly tighten these constraints. Our results are consistent with those from recent analyses of type Ia supernovae, cosmic microwave background anisotropies, the X-ray gas mass fraction of relaxed galaxy clusters, baryon acoustic oscillations and cosmic shear. Combining the new X-ray luminosity function data with current supernova, cosmic microwave background and cluster gas fraction data yields the improved constraints Omega_m=0.269 +- 0.016, sigma_8=0.82 +- 0.03 and w=-1.02 +- 0.06. (Abridged)
