Clusters of galaxies are outstanding laboratories for understanding the physics of supermassive black hole feedback. Here, we present the first \textit{Chandra}, Karl G. Janksy Very Large Array and \textit{Hubble Space Telescope} analysis of MACS J1447.4+0827 ($z = 0.3755$), one of the strongest cool core clusters known, in which extreme feedback from its central supermassive black hole is needed to prevent the hot intracluster gas from cooling. Using this multiwavelength approach, including 70 ks of \textit{Chandra} X-ray observations, we detect the presence of collimated jetted-outflows that coincides with a southern and a northern X-ray cavity. The total mechanical power associated with these outflows ($P_{\mathrm{cav}} \approx 6 \times 10^{44}$ erg s$^{-1}$) is roughly consistent with the energy required to prevent catastrophic cooling of the hot intracluster gas ($L_{\mathrm{cool}} = 1.71 \pm 0.01 \times 10^{45}$ erg s$^{-1}$ for t$_\mathrm{cool}$ = 7.7 Gyrs); implying that powerful supermassive black hole feedback has been in place several Giga-years ago in MACS J1447.7+0827. In addition, we detect the presence of a radio mini-halo that extends over 300 kpc in diameter ($P_{1.4 \mathrm{GHz}} = 3.0 \pm 0.3 \times 10^{24}$ W Hz$^{-1}$). The X-ray observations also reveal a $\sim20$ kpc plume-like structure that coincides with optical dusty filaments that surround the central galaxy. Overall, this study demonstrates that the various physical phenomena occurring in the most nearby clusters of galaxies are also occurring in their more distant analogues.

We present results from a 577 ks XMM-Newton observation of SPT-CL J0459-4947, the most distant cluster detected in the South Pole Telescope 2500 square degree (SPT-SZ) survey, and currently the most distant cluster discovered through its Sunyaev-Zel'dovich effect. The data confirm the cluster's high redshift, $z=1.71 \pm 0.02$, in agreement with earlier, less precise optical/IR photometric estimates. From the gas density profile, we estimate a characteristic mass of $M_{500}=(1.8 \pm 0.2) \times 10^{14}M_{Sun}$; cluster emission is detected above the background to a radius of $\sim 2.2 r_{500}$, or approximately the virial radius. The intracluster gas is characterized by an emission-weighted average temperature of $7.2 \pm 0.3$ keV and metallicity with respect to Solar of $0.37 \pm 0.08$. For the first time at such high redshift, this deep data set provides a measurement of metallicity outside the cluster center; at radii $r > 0.3 r_{500}$, we find it to be $0.33 \pm 0.17$, in good agreement with precise measurements at similar radii in the most nearby clusters, supporting an early enrichment scenario in which the bulk of the cluster gas is enriched to a universal metallicity prior to cluster formation, with little to no evolution thereafter. The leverage provided by the high redshift of this cluster tightens by a factor of 2 constraints on evolving metallicity models, when combined with previous measurements at lower redshifts.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

Galaxy clusters are the most massive virialized structures in the Universe and are formed through the gravitational accretion of matter over cosmic time. The discovery of an evolved galaxy cluster at redshift z=2, corresponding to a look-back time of 10.4 billion years, provides an opportunity to study its properties. The galaxy cluster XLSSC 122 was originally detected as a faint, extended X-ray source in the XMM Large Scale Structure survey and was revealed to be coincident with a compact over-density of galaxies with photometric redshifts of 1.9 +/- 0.2. Subsequent observations at millimetre wavelengths detected a Sunyaev-Zel'dovich decrement along the line of sight to XLSSC 122, thus confirming the existence of hot intracluster gas, while deep imaging spectroscopy from the European Space Agency's X-ray Multi-Mirror Mission (XMM-Newton) revealed an extended, X-ray bright gaseous atmosphere with a virial temperature of 60 million Kelvin, enriched with metals to the same extent as are local clusters. Here we report rest frame optical spectroscopic observations of XLSSC 122 and identify 37 member galaxies at a mean redshift of 1.98, corresponding to a look-back time of 10.4 billion years. We use photometry to determine a mean, dust-free stellar age of 2.98 billion years, indicating that star formation commenced in these galaxies at a mean redshift of 12, when the Universe was only 370 million years old. The full range of inferred formation redshifts, including the effects of dust, covers the interval from 7 to 13. These observations confirm that XLSSC 122 is a remarkably mature galaxy cluster with both evolved stellar populations in the member galaxies and a hot, metal-rich gas composing the intracluster medium.

We improve upon the cosmological constraints derived from the abundance and weak-lensing data of redMaPPer clusters detected in the Sloan Digital Sky Survey (SDSS). Specifically, we derive gas mass data using Chandra X-ray follow-up of a complete sample of the 30 richest SDSS redMaPPer clusters with $z\in[0.1,0.3]$, and use these additional data to improve upon the original analysis by Costanzi et al. (2019b). We simultaneously fit for the parameters of the richness-mass relation, the cluster gas mass-mass relation, and cosmology. By including our X-ray cluster sample in the SDSS cluster cosmology analysis, we measure $\Omega_{\rm m} = 0.25 \pm 0.04$ and $\sigma_8 = 0.85^{+0.06}_{-0.08}$. These constraints represent a 25.5% and 29.8% reduction in the size of the 68% confidence intervals of $\Omega_{\rm m}$ and $\sigma_8$ respectively, relative to the constraints published in Costanzi et al. (2019b). Our cosmological constraints are in agreement with early universe results from Planck. As a byproduct of our analysis, we also perform an independent calibration of the amplitude of the $\langle M_{\rm gas}^{\rm true}|M_{\rm 500c}\rangle$ scaling relation. Our calibration is consistent with and of comparable precision to that of Mantz et al. (2016b).

Weak-lensing measurements of the masses of galaxy clusters are commonly based on the assumption of spherically symmetric density profiles. Yet, the cold dark matter model predicts the shapes of dark matter halos to be triaxial. Halo triaxiality, and the orientation of the major axis with respect to the line of sight, are expected to be the leading cause of intrinsic scatter in weak-lensing mass measurements. The shape of central cluster galaxies (Brightest Cluster Galaxies; BCGs) is expected to follow the shape of the dark matter halo. Here we investigate the use of BCG ellipticity as predictor of the weak-lensing mass bias in individual clusters compared to the mean. Using weak lensing masses $M^{\rm WL}_{500}$ from the Weighing the Giants project, and $M_{500}$ derived from gas masses as low-scatter mass proxy, we find that, on average, the lensing masses of clusters with the roundest / most elliptical 25% of BCGs are biased $\sim 20$% high / low compared to the average, as qualitatively predicted by the cold dark matter model. For cluster cosmology projects utilizing weak-lensing mass estimates, the shape of the BCG can thus contribute useful information on the effect of orientation bias in weak lensing mass estimates as well as on cluster selection bias.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

The center determination of a galaxy cluster from an optical cluster finding algorithm can be offset from theoretical prescriptions or $N$-body definitions of its host halo center. These offsets impact the recovered cluster statistics, affecting both richness measurements and the weak lensing shear profile around the clusters. This paper models the centering performance of the \RM~cluster finding algorithm using archival X-ray observations of \RM-selected clusters. Assuming the X-ray emission peaks as the fiducial halo centers, and through analyzing their offsets to the \RM~centers, we find that $\sim 75\pm 8 \%$ of the \RM~clusters are well centered and the mis-centered offset follows a Gamma distribution in normalized, projected distance. These mis-centering offsets cause a systematic underestimation of cluster richness relative to the well-centered clusters, for which we propose a descriptive model. Our results enable the DES Y1 cluster cosmology analysis by characterizing the necessary corrections to both the weak lensing and richness abundance functions of the DES Y1 redMaPPer cluster catalog.

We derive cosmological constraints using a galaxy cluster sample selected from the 2500~deg$^2$ SPT-SZ survey. The sample spans the redshift range $0.25< z<1.75$ and contains 343 clusters with SZ detection significance $\xi>5$. The sample is supplemented with optical weak gravitational lensing measurements of 32 clusters with $0.29<z<1.13$ (from Magellan and HST) and X-ray measurements of 89 clusters with $0.25<z<1.75$ (from Chandra). We rely on minimal modeling assumptions: i) weak lensing provides an accurate means of measuring halo masses, ii) the mean SZ and X-ray observables are related to the true halo mass through power-law relations in mass and dimensionless Hubble parameter $E(z)$ with a-priori unknown parameters, iii) there is (correlated, lognormal) intrinsic scatter and measurement noise relating these observables to their mean relations. We simultaneously fit for these astrophysical modeling parameters and for cosmology. Assuming a flat $\nu\Lambda$CDM model, in which the sum of neutrino masses is a free parameter, we measure $\Omega_\mathrm{m}=0.276\pm0.047$, $\sigma_8=0.781\pm0.037$, and $\sigma_8(\Omega_\mathrm{m}/0.3)^{0.2}=0.766\pm0.025$. The redshift evolution of the X-ray $Y_\mathrm{X}$-mass and $M_\mathrm{gas}$-mass relations are both consistent with self-similar evolution to within $1\sigma$. The mass-slope of the $Y_\mathrm{X}$-mass relation shows a $2.3\sigma$ deviation from self-similarity. Similarly, the mass-slope of the $M_\mathrm{gas}$-mass relation is steeper than self-similarity at the $2.5\sigma$ level. In a $\nu w$CDM cosmology, we measure the dark energy equation of state parameter $w=-1.55\pm0.41$ from the cluster data. We perform a measurement of the growth of structure since redshift $z\sim1.7$ and find no evidence for tension with the prediction from General Relativity. We provide updated redshift and mass estimates for the SPT sample. (abridged)

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

We present measurements of the stellar mass fractions ($f_\star$) for a sample of high-redshift ($0.93 \le z \le 1.32$) infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS) and compare them to the stellar mass fractions of Sunyaev-Zel'dovich (SZ) effect-selected clusters in a similar mass and redshift range from the South Pole Telescope (SPT)-SZ Survey. We do not find a significant difference in mean $f_\star$ between the two selection methods, though we do find an unexpectedly large range in $f_\star$ for the SZ-selected clusters. In addition, we measure the luminosity function of the MaDCoWS clusters and find $m^*= 19.41\pm0.07$, similar to other studies of clusters at or near our redshift range. Finally, we present SZ detections and masses for seven MaDCoWS clusters and new spectroscopic redshifts for five MaDCoWS clusters. One of these new clusters, MOO J1521+0452 at $z=1.31$, is the most distant MaDCoWS cluster confirmed to date.

Using archival X-ray observations and a log-normal population model, we estimate constraints on the intrinsic scatter in halo mass at fixed optical richness for a galaxy cluster sample identified in Dark Energy Survey Year-One (DES-Y1) data with the redMaPPer algorithm. We examine the scaling behavior of X-ray temperatures, $T_X$, with optical richness, $\lambda_{RM}$, for clusters in the redshift range $0.2<z<0.7$. X-ray temperatures are obtained from Chandra and XMM observations for 58 and 110 redMaPPer systems, respectively. Despite non-uniform sky coverage, the $T_X$ measurements are $> 50\%$ complete for clusters with $\lambda_{RM} > 130$. Regression analysis on the two samples produces consistent posterior scaling parameters, from which we derive a combined constraint on the residual scatter, $\sigma_{\ln Tx | \lambda} = 0.275 \pm 0.019$. Joined with constraints for $T_X$ scaling with halo mass from the Weighing the Giants program and richness--temperature covariance estimates from the LoCuSS sample, we derive the richness-conditioned scatter in mass, $\sigma_{\ln M | \lambda} = 0.30 \pm 0.04\, _{({\rm stat})} \pm 0.09\, _{({\rm sys})}$, at an optical richness of approximately 70. Uncertainties in external parameters, particularly the slope and variance of the $T_X$--mass relation and the covariance of $T_X$ and $\lambda_{RM}$ at fixed mass, dominate the systematic error. The $95\%$ confidence region from joint sample analysis is relatively broad, $\sigma_{\ln M | \lambda} \in [0.14, \, 0.55]$, or a factor ten in variance.

Modern galaxy cluster science is a multi-wavelength endeavor with cornerstones provided by X-ray, optical/IR, mm, and radio measurements. In combination, these observations enable the construction of large, clean, complete cluster catalogs, and provide precise redshifts and robust mass calibration. The complementary nature of these multi-wavelength data dramatically reduces the impact of systematic effects that limit the utility of measurements made in any single waveband. The future of multi-wavelength cluster science is compelling, with cluster catalogs set to expand by orders of magnitude in size, and extend, for the first time, into the high-redshift regime where massive, virialized structures first formed. Unlocking astrophysical and cosmological insight from the coming catalogs will require new observing facilities that combine high spatial and spectral resolution with large collecting areas, as well as concurrent advances in simulation modeling campaigns. Together, future multi-wavelength observations will resolve the thermodynamic structure in and around the first groups and clusters, distinguishing the signals from active and star-forming galaxies, and unveiling the interrelated stories of galaxy evolution and structure formation during the epoch of peak cosmic activity.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

Outbursts from active galactic nuclei (AGN) can inflate cavities in the intracluster medium (ICM) of galaxy clusters and are believed to play the primary role in offsetting radiative cooling in the ICM. However, the details of how the energy from AGN feedback thermalizes to heat the ICM is not well understood, partly due to the unknown composition and energetics of the cavities. The Sunyaev-Zel'dovich (SZ) effect, a measure of the integrated pressure along the line of sight, provides a means of measuring the thermal contents of the cavities, to discriminate between thermal, nonthermal, and other sources of pressure support. Here we report measurements of the SZ effect at 30 GHz toward the galaxy cluster MS 0735.6+7421 (MS0735), using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). MS0735 hosts the most energetic AGN outburst known and lobes of radio synchrotron emission coincident with a pair of giant X-ray cavities $\sim 200$ across. Our CARMA maps show a clear deficit in the SZ signal coincident with the X-ray identified cavities, when compared to a smooth X-ray derived pressure model. We find that the cavities have very little SZ-contributing material, suggesting that they are either supported by very diffuse thermal plasma with temperature in excess of hundreds of keV, or are not supported thermally. Our results represent the first detection (with $4.4 \sigma$ significance) of this phenomenon with the SZ effect.

We present first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg$^2$ South Pole Telescope Sunyaev Zel'dovich effect (SPT-SZ) survey. The cluster selection is to first order independent of galaxy properties, making the SPT-SZ sample particularly well suited for cluster galaxy population studies. We carry out a 4-band imaging campaign with the {\it Hubble} and {\it Spitzer} Space Telescopes of the five $z\gtrsim 1.4$, S/N$_{SZE}>$5 clusters, that are among the rarest most massive clusters known at this redshift. All five show clear overdensities of red galaxies whose colors agree with the initial cluster redshift estimates. The highest redshift cluster in this sample, SPT-CLJ0459-4947 at $z\sim1.72$, is the most distant $M_{500}>10^{14}~M_{\odot}$ ICM-selected cluster discovered thus far, and is one of only three known clusters in this mass range at $z\gtrsim 1.7$, regardless of selection. Based on UVJ-like photometric classification of quiescent and star-forming galaxies, the passive fraction in the cluster central regions ($r/r_{500}<0.7$) is higher than in the field at the same redshift, with corresponding environmental quenching efficiencies typically in the range $\sim0.5-0.8$ for stellar masses $\log(M/M_{\odot})>10.85$. We have explored the impact of emission from star formation on the selection of this sample, concluding that all five clusters studied here would still have been detected with S/N$_{SZE}>$5, even if they had the same passive fraction as measured in the field. Our results thus point towards an efficient suppression of star formation in massive galaxies in the central regions of the most massive clusters, occurring already earlier than $z\sim1.5$. [Abridged]

The finite sensitivity of instruments or detection methods means that data sets in many areas of astronomy, for example cosmological or exoplanet surveys, are necessarily systematically incomplete. Such data sets, where the population being investigated is of unknown size and only partially represented in the data, are called "truncated" in the statistical literature. Truncation can be accounted for through a relatively straightforward modification to the model being fitted in many circumstances, provided that the model can be extended to describe the population of undetected sources. Here I examine the problem of regression using truncated data in general terms, and use a simple example to show the impact of selecting a subset of potential data on the dependent variable, on the independent variable, and on a second dependent variable that is correlated with the variable of interest. Special circumstances in which selection effects are ignorable are noted. I also comment on computational strategies for performing regression with truncated data, as an extension of methods that have become popular for the non-truncated case, and provide some general recommendations.

We present the Massive and Distant Clusters of WISE Survey (MaDCoWS), a search for galaxy clusters at 0.7<z<1.5 based upon data from the Wide-field Infrared Survey Explorer (WISE) mission. MaDCoWS is the first cluster survey capable of discovering massive clusters at these redshifts over the full extragalactic sky. The search is divided into two regions -- the region of the extragalactic sky covered by Pan-STARRS (Dec>-30 degrees) and the remainder of the southern extragalactic sky at Dec<-30 degrees for which shallower optical data from SuperCOSMOS Sky Survey are available. In this paper we describe the search algorithm, characterize the sample, and present the first MaDCoWS data release -- catalogs of the 2433 highest amplitude detections in the WISE--Pan-STARRS region and the 250 highest amplitude detections in the WISE--SuperCOSMOS region. A total of 1723 of the detections from the WISE--Pan-STARRS sample have also been observed with the Spitzer Space Telescope, providing photometric redshifts and richnesses, and an additional 64 detections within the WISE--SuperCOSMOS region also have photometric redshifts and richnesses. Spectroscopic redshifts for 38 MaDCoWS clusters with IRAC photometry demonstrate that the photometric redshifts have an uncertainty of $\sigma_z/(1+z)\sim0.036$. Combining the richness measurements with Sunyaev-Zel'dovich observations of MaDCoWS clusters, we also present a preliminary mass-richness relation that can be used to infer the approximate mass distribution of the full sample. The estimated median mass for the WISE--Pan-STARRS catalog is $M_{500}=1.6^{+0.7}_{-0.8}\times10^{14} \mathrm{M}_\odot$, with the Sunyaev-Zel'dovich data confirming that we detect clusters with masses up to $M_{500}\sim5\times10^{14} \mathrm{M}_\odot$ $(M_{200}\sim10^{15} \mathrm{M}_\odot)$.

Uncertainty in the mass-observable scaling relations is currently the limiting factor for galaxy cluster based cosmology. Weak gravitational lensing can provide a direct mass calibration and reduce the mass uncertainty. We present new ground-based weak lensing observations of 19 South Pole Telescope (SPT) selected clusters at redshifts $0.29 \leq z \leq 0.61$ and combine them with previously reported space-based observations of 13 galaxy clusters at redshifts $0.576 \leq z \leq 1.132$ to constrain the cluster mass scaling relations with the Sunyaev-Zel'dovich effect (SZE), the cluster gas mass \mgas, and \yx, the product of \mgas\ and X-ray temperature. We extend a previously used framework for the analysis of scaling relations and cosmological constraints obtained from SPT-selected clusters to make use of weak lensing information. We introduce a new approach to estimate the effective average redshift distribution of background galaxies and quantify a number of systematic errors affecting the weak lensing modelling. These errors include a calibration of the bias incurred by fitting a Navarro-Frenk-White profile to the reduced shear using $N$-body simulations. We blind the analysis to avoid confirmation bias. We are able to limit the systematic uncertainties to 5.6% in cluster mass (68% confidence). Our constraints on the mass--X-ray observable scaling relations parameters are consistent with those obtained by earlier studies, and our constraints for the mass--SZE scaling relation are consistent with the simulation-based prior used in the most recent SPT-SZ cosmology analysis. We can now replace the external mass calibration priors used in previous SPT-SZ cosmology studies with a direct, internal calibration obtained on the same clusters.

We perform the first blind analysis of cluster abundance data. Specifically, we derive cosmological constraints from the abundance and weak-lensing signal of \redmapper\ clusters of richness $\lambda\geq 20$ in the redshift range $z\in[0.1,0.3]$ as measured in the Sloan Digital Sky Survey (SDSS). We simultaneously fit for cosmological parameters and the richness--mass relation of the clusters. For a flat $\Lambda$CDM cosmological model with massive neutrinos, we find $S_8 \equiv \sigma_{8}(\Omega_m/0.3)^{0.5}=0.79^{+0.05}_{-0.04}$. This value is both consistent and competitive with that derived from cluster catalogues selected in different wavelengths. Our result is also consistent with the combined probes analyses by the Dark Energy Survey (DES) and the Kilo-Degree Survey (KiDS), and with the Cosmic Microwave Background (CMB) anisotropies as measured by \planck. We demonstrate that the cosmological posteriors are robust against variation of the richness--mass relation model and to systematics associated with the calibration of the selection function. In combination with Baryon Acoustic Oscillation (BAO) data and Big-Bang Nucleosynthesis (BBN) data, we constrain the Hubble rate to be $h=0.66\pm 0.02$, independent of the CMB. Future work aimed at improving our understanding of the scatter of the richness--mass relation has the potential to significantly improve the precision of our cosmological posteriors. The methods described in this work were developed for use in the forthcoming analysis of cluster abundances in the DES. Our SDSS analysis constitutes the first part of a staged-unblinding analysis of the full DES data set.

We present a multi-wavelength analysis of the four most relaxed clusters in the South Pole Telescope 2500 deg^2 survey, which lie at 0.55 < z < 0.75. This study, which utilizes new, deep data from Chandra and Hubble, along with ground-based spectroscopy from Gemini and Magellan, improves significantly on previous studies in both depth and angular resolution, allowing us to directly compare to clusters at z~0. We find that the temperature, density, and entropy profiles of the intracluster medium (ICM) are very similar among the four clusters, and share similar shapes to clusters at z~0. Specifically, we find no evidence for deviations from self similarity in the temperature profile over the radial range 10kpc < r < 1Mpc, implying that the processes responsible for preventing runaway cooling over the past >6 Gyr are, at least roughly, preserving self similarity. We find typical metallicities of ~0.3 Zsun in the bulk of the ICM, rising to ~0.5 Zsun in the inner ~100 kpc, and reaching ~1 Zsun at r < 10kpc. This central excess is similar in magnitude to what is observed in the most relaxed clusters at z~0, suggesting that both the global metallicity and the central excess that we see in cool core clusters at z~0 were in place very early in the cluster lifetime and, specifically, that the central excess is not due to late-time enrichment by the central galaxy. Consistent with observations at z~0, we measure a diversity of stellar populations in the central brightest cluster galaxies of these four clusters, with star formation rates spanning a factor of ~500, despite the similarity in cooling time, cooling rate, and central entropy. These data suggest that, while the details vary dramatically from system to system, runaway cooling has been broadly regulated in relaxed clusters over the past 6 Gyr.

We constrain the mass--richness scaling relation of redMaPPer galaxy clusters identified in the Dark Energy Survey Year 1 data using weak gravitational lensing. We split clusters into $4\times3$ bins of richness $\lambda$ and redshift $z$ for $\lambda\geq20$ and $0.2 \leq z \leq 0.65$ and measure the mean masses of these bins using their stacked weak lensing signal. By modeling the scaling relation as $\langle M_{\rm 200m}|\lambda,z\rangle = M_0 (\lambda/40)^F ((1+z)/1.35)^G$, we constrain the normalization of the scaling relation at the 5.0 per cent level as $M_0 = [3.081 \pm 0.075 ({\rm stat}) \pm 0.133 ({\rm sys})] \cdot 10^{14}\ {\rm M}_\odot$ at $\lambda=40$ and $z=0.35$. The richness scaling index is constrained to be $F=1.356 \pm 0.051\ ({\rm stat})\pm 0.008\ ({\rm sys})$ and the redshift scaling index $G=-0.30\pm 0.30\ ({\rm stat})\pm 0.06\ ({\rm sys})$. These are the tightest measurements of the normalization and richness scaling index made to date. We use a semi-analytic covariance matrix to characterize the statistical errors in the recovered weak lensing profiles. Our analysis accounts for the following sources of systematic error: shear and photometric redshift errors, cluster miscentering, cluster member dilution of the source sample, systematic uncertainties in the modeling of the halo--mass correlation function, halo triaxiality, and projection effects. We discuss prospects for reducing this systematic error budget, which dominates the uncertainty on $M_0$. Our result is in excellent agreement with, but has significantly smaller uncertainties than, previous measurements in the literature, and augurs well for the power of the DES cluster survey as a tool for precision cosmology and upcoming galaxy surveys such as LSST, Euclid and WFIRST.

The Science Products Module (SPM), a US contribution to the Athena Wide Field Imager, is a highly capable secondary CPU that performs special processing on the science data stream. The SPM will have access to both accepted X-ray events and those that were rejected by the on-board event recognition processing. It will include two software modules. The Transient Analysis Module will perform on-board processing of the science images to identify and characterize variability of the prime target and/or detection of serendipitous transient X-ray sources in the field of view. The Background Analysis Module will perform more sophisticated flagging of potential background events as well as improved background characterization, making use of data that are not telemetered to the ground, to provide improved background maps and spectra. We present the preliminary design of the SPM hardware as well as a brief overview of the software algorithms under development.

The cosmological utility of galaxy cluster catalogues is primarily limited by our ability to calibrate the relation between halo mass and observable mass proxies such as cluster richness, X-ray luminosity or the Sunyaev-Zeldovich signal. Projection effects are a particularly pernicious systematic effect that can impact observable mass proxies; structure along the line of sight can both bias and increase the scatter of the observable mass proxies used in cluster abundance studies. In this work, we develop an empirical method to characterize the impact of projection effects on redMaPPer cluster catalogues. We use numerical simulations to validate our method and illustrate its robustness. We demonstrate that modeling of projection effects is a necessary component for cluster abundance studies capable of reaching $\approx 5\%$ mass calibration uncertainties (e.g. the Dark Energy Survey Year 1 sample). Specifically, ignoring the impact of projection effects in the observable--mass relation --- i.e. marginalizing over a log-normal model only --- biases the posterior of the cluster normalization condition $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{1/2}$ by $\Delta S_8 =0.05$, more than twice the uncertainty in the posterior for such an analysis.

We present spectroscopic confirmation of five galaxy clusters at $1.25 < \textit{z} < 1.5$, discovered in the $2500$ deg$^{2}$ South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. These clusters, taken from a mass-limited sample with a nearly redshift independent selection function, have multi-wavelength follow-up imaging data from the X-ray to near-infrared, and currently form the most homogeneous massive high-redshift cluster sample known. We identify $44$ member galaxies, along with $25$ field galaxies, among the five clusters, and describe the full set of observations and data products from Magellan/LDSS3 multi-object spectroscopy of these cluster fields. We briefly describe the analysis pipeline, and present ensemble analyses of cluster member galaxies that demonstrate the reliability of the measured redshifts. We report $\textit{z} = 1.259, 1.288, 1.316, 1.401$ and $1.474$ for the five clusters from a combination of absorption-line (Ca II H$\&$K doublet - $3968,3934$ {\AA}) and emission-line ([OII] $3727,3729$ {\AA}) spectral features. Moreover, the calculated velocity dispersions yield dynamical cluster masses in good agreement with SZ masses for these clusters. We discuss the velocity and spatial distributions of passive and [OII]-emitting galaxies in these clusters, showing that they are consistent with velocity segregation and biases observed in lower redshift SPT clusters. We identify modest [OII] emission and pronounced CN and H$\delta$ absorption in a stacked spectrum of $28$ passive galaxies with Ca II H$\&$K-derived redshifts. This work increases the number of spectroscopically-confirmed SZ-selected galaxy clusters at $\textit{z} > 1.25$ from three to eight, further demonstrating the efficacy of SZ selection for the highest redshift massive clusters, and enabling detailed study of these systems.

We jointly analyze Bolocam Sunyaev-Zeldovich (SZ) effect and Chandra X-ray data for a set of 45 clusters to derive gas density and temperature profiles without using spectroscopic information. The sample spans the mass and redshift range $3 \times 10^{14} M_{\odot} \le M_{500} \le 25 \times 10^{14} M_{\odot}$ and $0.15\le z \le 0.89$. We define cool-core (CC) and non-cool core (NCC) subsamples based on the central X-ray luminosity, and 17/45 clusters are classified as CC. In general, the profiles derived from our analysis are found to be in good agreement with previous analyses, and profile constraints beyond $r_{500}$ are obtained for 34/45 clusters. In approximately 30% of the CC clusters our analysis shows a central temperature drop with a statistical significance of $>3\sigma$; this modest detection fraction is due mainly to a combination of coarse angular resolution and modest S/N in the SZ data. Most clusters are consistent with an isothermal profile at the largest radii near $r_{500}$, although 9/45 show a significant temperature decrease with increasing radius. The sample mean density profile is in good agreement with previous studies, and shows a minimum intrinsic scatter of approximately 10% near $0.5 \times r_{500}$. The sample mean temperature profile is consistent with isothermal, and has an intrinsic scatter of approximately 50% independent of radius. This scatter is significantly higher compared to earlier X-ray-only studies, which find intrinsic scatters near 10%, likely due to a combination of unaccounted for non-idealities in the SZ noise, projection effects, and sample selection.

We constrain cold dark energy of negligible sound speed using galaxy cluster abundance observations. In contrast to standard quasi-homogeneous dark energy, negligible sound speed implies clustering of the dark energy fluid at all scales, allowing us to measure the effects of dark energy perturbations at cluster scales. We compare those models and set the stage for using non-linear information from semi-analytical modelling in cluster growth data analyses. For this, we recalibrate the halo mass function with non-linear characteristic quantities, the spherical collapse threshold and virial overdensity, that account for model and redshift dependent behaviours, as well as an additional mass contribution for cold dark energy. We present the first constraints from this cold dark matter plus cold dark energy mass function using our cluster abundance likelihood, which self-consistently accounts for selection effects, covariances and systematic uncertainties. We combine cluster growth data with CMB, SNe Ia and BAO data, and find a shift between cold versus quasi-homogeneous dark energy of up to $1\sigma$. We make a Fisher matrix forecast of constraints attainable with cluster growth data from the on-going Dark Energy Survey (DES). For DES, we predict $\sim$50$\%$ tighter constraints on $\left(\Omega_\mathrm{m},w \right)$ for cold dark energy versus $w$CDM models, with the same free parameters. Overall, we show that cluster abundance analyses are sensitive to cold dark energy, an alternative, viable model that should be routinely investigated alongside the standard dark energy scenario.

We present an HST/ACS weak gravitational lensing analysis of 13 massive high-redshift (z_median=0.88) galaxy clusters discovered in the South Pole Telescope (SPT) Sunyaev-Zel'dovich Survey. This study is part of a larger campaign that aims to robustly calibrate mass-observable scaling relations over a wide range in redshift to enable improved cosmological constraints from the SPT cluster sample. We introduce new strategies to ensure that systematics in the lensing analysis do not degrade constraints on cluster scaling relations significantly. First, we efficiently remove cluster members from the source sample by selecting very blue galaxies in V-I colour. Our estimate of the source redshift distribution is based on CANDELS data, where we carefully mimic the source selection criteria of the cluster fields. We apply a statistical correction for systematic photometric redshift errors as derived from Hubble Ultra Deep Field data and verified through spatial cross-correlations. We account for the impact of lensing magnification on the source redshift distribution, finding that this is particularly relevant for shallower surveys. Finally, we account for biases in the mass modelling caused by miscentring and uncertainties in the concentration-mass relation using simulations. In combination with temperature estimates from Chandra we constrain the normalisation of the mass-temperature scaling relation ln(E(z) M_500c/10^14 M_sun)=A+1.5 ln(kT/7.2keV) to A=1.81^{+0.24}_{-0.14}(stat.) +/- 0.09(sys.), consistent with self-similar redshift evolution when compared to lower redshift samples. Additionally, the lensing data constrain the average concentration of the clusters to c_200c=5.6^{+3.7}_{-1.8}.

The cosmological constraining power of modern galaxy cluster catalogs can be improved by obtaining low-scatter mass proxy measurements for even a small fraction of sources. In the context of large upcoming surveys that will reveal the cluster population down to the group scale and out to high redshifts, efficient strategies for obtaining such mass proxies will be valuable. In this work, we use high-quality weak lensing and X-ray mass estimates for massive clusters in current X-ray selected catalogs to revisit the scaling relations of the projected, center-excised X-ray luminosity ($L_{ce}$), which previous work suggests correlates tightly with total mass. Our data confirm that this is the case, with $L_{ce}$ having an intrinsic scatter at fixed mass comparable to that of gas mass, temperature or $Y_X$. Compared to these other proxies, however, $L_{ce}$ is less susceptible to systematic uncertainties due to background modeling, and can be measured precisely with shorter exposures. This opens up the possibility of using $L_{ce}$ to estimate masses for large numbers of clusters discovered by new X-ray surveys (e.g. eROSITA) directly from the survey data, as well as for clusters discovered at other wavelengths, with relatively short follow-up observations. We describe a simple procedure for making such estimates from X-ray surface brightness data, and comment on the spatial resolution required to apply this method as a function of cluster mass and redshift. We also explore the potential impact of Chandra and XMM-Newton follow-up observations over the next decade on dark energy constraints from new cluster surveys.

We use Chandra X-ray data to measure the metallicity of the intracluster medium (ICM) in 245 massive galaxy clusters selected from X-ray and Sunyaev-Zel'dovich (SZ) effect surveys, spanning redshifts $0<z<1.2$. Metallicities were measured in three different radial ranges, spanning cluster cores through their outskirts. We explore trends in these measurements as a function of cluster redshift, temperature, and surface brightness "peakiness" (a proxy for gas cooling efficiency in cluster centers). The data at large radii (0.5--1 $r_{500}$) are consistent with a constant metallicity, while at intermediate radii (0.1-0.5 $r_{500}$) we see a late-time increase in enrichment, consistent with the expected production and mixing of metals in cluster cores. In cluster centers, there are strong trends of metallicity with temperature and peakiness, reflecting enhanced metal production in the lowest-entropy gas. Within the cool-core/sharply peaked cluster population, there is a large intrinsic scatter in central metallicity and no overall evolution, indicating significant astrophysical variations in the efficiency of enrichment. The central metallicity in clusters with flat surface brightness profiles is lower, with a smaller intrinsic scatter, but increases towards lower redshifts. Our results are consistent with other recent measurements of ICM metallicity as a function of redshift. They reinforce the picture implied by observations of uniform metal distributions in the outskirts of nearby clusters, in which most of the enrichment of the ICM takes place before cluster formation, with significant later enrichment taking place only in cluster centers, as the stellar populations of the central galaxies evolve.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We present the statistical analysis of X-ray surface brightness and gas density fluctuations in cool cores of ten, nearby and bright galaxy clusters that have deep Chandra observations and show observational indications of radio-mechanical AGN feedback. Within the central parts of cool cores the total variance of fluctuations is dominated by isobaric and/or isothermal fluctuations on spatial scales ~ 10-60 kpc, which are likely associated with slow gas motions and bubbles of relativistic plasma. Adiabatic fluctuations associated with weak shocks constitute less than 10 per cent of the total variance in all clusters. The typical amplitude of density fluctuations is small, ~ 10 per cent or less on scales of ~ 10-15 kpc. Subdominant contribution of adiabatic fluctuations and small amplitude of density fluctuations support a model of gentle AGN feedback as opposed to periodically explosive scenarios which are implemented in some numerical simulations. Measured one-component velocities of gas motions are typically below 100-150 km/s on scales < 50 kpc, and can be up to ~ 300 km/s on ~ 100 kpc scales. The non-thermal energy is < 12 per cent of the thermal energy. Regardless of the source that drives these motions the dissipation of the energy in such motions provides heat that is sufficient to balance radiative cooling on average, albeit the uncertainties are large. Presented results here support previous conclusions based on the analysis of the Virgo and Perseus Clusters, and agree with the Hitomi measurements. With next generation observatories like Athena and Lynx, these techniques will be yet more powerful.

Suzaku measurements of a homogeneous metal distribution of $Z\sim0.3$ Solar in the outskirts of the nearby Perseus cluster suggest that chemical elements were deposited and mixed into the intergalactic medium before clusters formed, likely over 10 billion years ago. A key prediction of this early enrichment scenario is that the intracluster medium in all massive clusters should be uniformly enriched to a similar level. Here, we confirm this prediction by determining the iron abundances in the outskirts ($r>0.25r_{200}$) of a sample of ten other nearby galaxy clusters observed with Suzaku for which robust measurements based on the Fe-K lines can be made. Across our sample the iron abundances are consistent with a constant value, $Z_{\rm Fe}=0.316\pm0.012$ Solar ($\chi^2=28.85$ for 25 degrees of freedom). This is remarkably similar to the measurements for the Perseus cluster of $Z_{\rm Fe}=0.314\pm0.012$ Solar, using the Solar abundance scale of Asplund et al. (2009).

We present results from Suzaku Key Project observations of the Virgo Cluster, the nearest galaxy cluster to us, mapping its X-ray properties along four long `arms' extending beyond the virial radius. The entropy profiles along all four azimuths increase with radius, then level out beyond $0.5r_{200}$, while the average pressure at large radii exceeds Planck Sunyaev-Zel'dovich measurements. These results can be explained by enhanced gas density fluctuations (clumping) in the cluster's outskirts. Using a standard Navarro, Frenk and White (1997) model, we estimate a virial mass, radius, and concentration parameter of $M_{200}=1.05\pm0.02\times10^{14}$ M$_\odot$, $r_{200}=974.1\pm5.7$ kpc, and $c = 8.8 \pm0.2$, respectively. The inferred cumulative baryon fraction exceeds the cosmic mean at $r\sim r_{200}$ along the major axis, suggesting enhanced gas clumping possibly sourced by a candidate large-scale structure filament along the north-south direction. The Suzaku data reveal a large-scale sloshing pattern, with two new cold fronts detected at radii of 233 kpc and 280 kpc along the western and southern arms, respectively. Two high-temperature regions are also identified 1 Mpc towards the south and 605 kpc towards the west of M87, likely representing shocks associated with the ongoing cluster growth. Although systematic uncertainties in measuring the metallicity for low temperature plasma remain, the data at large radii appear consistent with a uniform metal distribution on scales of $\sim 90\times180$ kpc and larger, providing additional support for the early chemical enrichment scenario driven by galactic winds at redshifts of 2-3.

We present results from a 100 ks XMM-Newton observation of galaxy cluster XLSSC 122, the first massive cluster discovered through its X-ray emission at $z\approx2$. The data provide the first precise constraints on the bulk thermodynamic properties of such a distant cluster, as well as an X-ray spectroscopic confirmation of its redshift. We measure an average temperature of $kT=5.0\pm0.7$ keV; a metallicity with respect to solar of $Z/Z_{\odot}=0.33^{+0.19}_{-0.17}$, consistent with lower-redshift clusters; and a redshift of $z=1.99^{+0.07}_{-0.06}$, consistent with the earlier photo-z estimate. The measured gas density profile leads to a mass estimate at $r_{500}$ of $M_{500}=(6.3\pm1.5)\times10^{13}M_{\odot}$. From CARMA 30 GHz data, we measure the spherically integrated Compton parameter within $r_{500}$ to be $Y_{500}=(3.6\pm0.4)\times10^{-12}$. We compare the measured properties of XLSSC 122 to lower-redshift cluster samples, and find good agreement when assuming the simplest (self-similar) form for the evolution of cluster scaling relations. While a single cluster provides limited information, this result suggests that the evolution of the intracluster medium in the most massive, well developed clusters is remarkably simple, even out to the highest redshifts where they have been found. At the same time, our data reaffirm the previously reported spatial offset between the centers of the X-ray and SZ signals for XLSSC 122, suggesting a disturbed configuration. Higher spatial resolution data could thus provide greater insights into the internal dynamics of this system.

We present the results of a Chandra X-ray survey of the 8 most massive galaxy clusters at z>1.2 in the South Pole Telescope 2500 deg^2 survey. We combine this sample with previously-published Chandra observations of 49 massive X-ray-selected clusters at 0<z<0.1 and 90 SZ-selected clusters at 0.25<z<1.2 to constrain the evolution of the intracluster medium (ICM) over the past ~10 Gyr. We find that the bulk of the ICM has evolved self similarly over the full redshift range probed here, with the ICM density at r>0.2R500 scaling like E(z)^2. In the centers of clusters (r<0.1R500), we find significant deviations from self similarity (n_e ~ E(z)^{0.1+/-0.5}), consistent with no redshift dependence. When we isolate clusters with over-dense cores (i.e., cool cores), we find that the average over-density profile has not evolved with redshift -- that is, cool cores have not changed in size, density, or total mass over the past ~9-10 Gyr. We show that the evolving "cuspiness" of clusters in the X-ray, reported by several previous studies, can be understood in the context of a cool core with fixed properties embedded in a self similarly-evolving cluster. We find no measurable evolution in the X-ray morphology of massive clusters, seemingly in tension with the rapidly-rising (with redshift) rate of major mergers predicted by cosmological simulations. We show that these two results can be brought into agreement if we assume that the relaxation time after a merger is proportional to the crossing time, since the latter is proportional to H(z)^(-1).

We present significantly improved measurements of turbulent velocities in the hot gaseous halos of nearby giant elliptical galaxies. Using deep XMM-Newton Reflection Grating Spectrometer (RGS ) observations and a combination of resonance scattering and direct line broadening methods, we obtain well bounded constraints for 13 galaxies. Assuming that the turbulence is isotropic, we obtain a best fit mean 1D turbulent velocity of ~110 km/s. This implies a typical 3D Mach number ~0.45 and a typical non-thermal pressure contribution of ~6 per cent in the cores of nearby massive galaxies. The intrinsic scatter around these values is modest - consistent with zero, albeit with large statistical uncertainty - hinting at a common and quasi-continuous mechanism sourcing the velocity structure in these objects. Using conservative estimates of the spatial scales associated with the observed turbulent motions, we find that turbulent heating can be sufficient to offset radiative cooling in the inner regions of these galaxies (<10 kpc, typically 2-3 kpc). The full potential of our analysis methods will be enabled by future X-ray microcalorimeter observations.

We present constraints on the scaling relations of galaxy cluster X-ray luminosity, temperature and gas mass (and derived quantities) with mass and redshift, employing masses from robust weak gravitational lensing measurements. These are the first such results obtained from an analysis that simultaneously accounts for selection effects and the underlying mass function, and directly incorporates lensing data to constrain total masses. Our constraints on the scaling relations and their intrinsic scatters are in good agreement with previous studies, and reinforce a picture in which departures from self-similar scaling laws are primarily limited to cluster cores. However, the data are beginning to reveal new features that have implications for cluster astrophysics and provide new tests for hydrodynamical simulations. We find a positive correlation in the intrinsic scatters of luminosity and temperature at fixed mass, which is related to the dynamical state of the clusters. While the evolution of the nominal scaling relations over the redshift range $0.0<z<0.5$ is consistent with self similarity, we find tentative evidence that the luminosity and temperature scatters respectively decrease and increase with redshift. Physically, this likely related to the development of cool cores and the rate of major mergers. We also examine the scaling relations of redMaPPer richness and Compton $Y$ from Planck. While the richness--mass relation is in excellent agreement with recent work, the measured $Y$--mass relation departs strongly from that assumed in the Planck cluster cosmology analysis. The latter result is consistent with earlier comparisons of lensing and Planck scaling-relation-derived masses.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We present the results of SPT-GMOS, a spectroscopic survey with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South. The targets of SPT-GMOS are galaxy clusters identified in the SPT-SZ survey, a millimeter-wave survey of 2500 squ. deg. of the southern sky using the South Pole Telescope (SPT). Multi-object spectroscopic observations of 62 SPT-selected galaxy clusters were performed between January 2011 and December 2015, yielding spectra with radial velocity measurements for 2595 sources. We identify 2243 of these sources as galaxies, and 352 as stars. Of the galaxies, we identify 1579 as members of SPT-SZ galaxy clusters. The primary goal of these observations was to obtain spectra of cluster member galaxies to estimate cluster redshifts and velocity dispersions. We describe the full spectroscopic dataset and resulting data products, including galaxy redshifts, cluster redshifts and velocity dispersions, and measurements of several well-known spectral indices for each galaxy: the equivalent width, W, of [O II] 3727,3729 and H-delta, and the 4000A break strength, D4000. We use the spectral indices to classify galaxies by spectral type (i.e., passive, post-starburst, star-forming), and we match the spectra against photometric catalogs to characterize spectroscopically-observed cluster members as a function of brightness (relative to m*). Finally, we report several new measurements of redshifts for ten bright, strongly-lensed background galaxies in the cores of eight galaxy clusters. Combining the SPT-GMOS dataset with previous spectroscopic follow-up of SPT-SZ galaxy clusters results in spectroscopic measurements for >100 clusters, or ~20% of the full SPT-SZ sample.

This is the fifth in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Our sample comprises 40 clusters identified as being dynamically relaxed and hot in Papers I and II of this series. Here we use constraints on cluster mass profiles from X-ray data to test some of the basic predictions of cosmological structure formation in the Cold Dark Matter (CDM) paradigm. We present constraints on the concentration--mass relation for massive clusters, finding a power-law mass dependence with a slope of $\kappa_m=-0.16\pm0.07$, in agreement with CDM predictions. For this relaxed sample, the relation is consistent with a constant as a function of redshift (power-law slope with $1+z$ of $\kappa_\zeta=-0.17\pm0.26$), with an intrinsic scatter of $\sigma_{\ln c}=0.16\pm0.03$. We investigate the shape of cluster mass profiles over the radial range probed by the data (typically $\sim50$kpc--1Mpc), and test for departures from the simple Navarro, Frenk & White (NFW) form, for which the logarithmic slope of the density profile tends to $-1$ at small radii. Specifically, we consider as alternatives the generalized NFW (GNFW) and Einasto parametrizations. For the GNFW model, we find an average value of (minus) the logarithmic inner slope of $\beta=1.02\pm0.08$, with an intrinsic scatter of $\sigma_\beta=0.22\pm0.07$, while in the Einasto case we constrain the average shape parameter to be $\alpha=0.29\pm0.04$ with an intrinsic scatter of $\sigma_\alpha=0.12\pm0.04$. Our results are thus consistent with the simple NFW model on average, but we clearly detect the presence of intrinsic, cluster-to-cluster scatter about the average.

We measure the SZ signal toward a set of 47 clusters with a median mass of $9.5 \times 10^{14}$ M$_{\odot}$ and a median redshift of 0.40 using data from Planck and the ground-based Bolocam receiver. When Planck XMM-like masses are used to set the scale radius $\theta_{\textrm{s}}$, we find consistency between the integrated SZ signal, $Y_{\textrm{5R500}}$, derived from Bolocam and Planck based on gNFW model fits using A10 shape parameters, with an average ratio of $1.069 \pm 0.030$ (allowing for the $\simeq 5$% Bolocam flux calibration uncertainty). We also perform a joint fit to the Bolocam and Planck data using a modified A10 model with the outer logarithmic slope $\beta$ allowed to vary, finding $\beta = 6.13 \pm 0.16 \pm 0.76$ (measurement error followed by intrinsic scatter). In addition, we find that the value of $\beta$ scales with mass and redshift according to $\beta \propto M^{0.077 \pm 0.026} \times (1+z)^{-0.06 \pm 0.09}$. This mass scaling is in good agreement with recent simulations. We do not observe the strong trend of $\beta$ with redshift seen in simulations, though we conclude that this is most likely due to our sample selection. Finally, we use Bolocam measurements of $Y_{500}$ to test the accuracy of the Planck completeness estimate. We find consistency, with the actual number of Planck detections falling approximately $1 \sigma$ below the expectation from Bolocam. We translate this small difference into a constraint on the the effective mass bias for the Planck cluster cosmology results, with $(1-b) = 0.93 \pm 0.06$.

(abridged) We present cosmological constraints obtained from galaxy clusters identified by their Sunyaev-Zel'dovich effect signature in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich survey. We consider the 377 cluster candidates identified at z>0.25 with a detection significance greater than five, corresponding to the 95% purity threshold for the survey. We compute constraints on cosmological models using the measured cluster abundance as a function of mass and redshift. We include additional constraints from multi-wavelength observations, including Chandra X-ray data for 82 clusters and a weak lensing-based prior on the normalization of the mass-observable scaling relations. Assuming a LCDM cosmology, where the species-summed neutrino mass has the minimum allowed value (mnu = 0.06 eV) from neutrino oscillation experiments, we combine the cluster data with a prior on H0 and find sigma_8 = 0.797+-0.031 and Omega_m = 0.289+-0.042, with the parameter combination sigma_8(Omega_m/0.27)^0.3 = 0.784+-0.039. These results are in good agreement with constraints from the CMB from SPT, WMAP, and Planck, as well as with constraints from other cluster datasets. Adding mnu as a free parameter, we find mnu = 0.14+-0.08 eV when combining the SPT cluster data with Planck CMB data and BAO data, consistent with the minimum allowed value. Finally, we consider a cosmology where mnu and N_eff are fixed to the LCDM values, but the dark energy equation of state parameter w is free. Using the SPT cluster data in combination with an H0 prior, we measure w = -1.28+-0.31, a constraint consistent with the LCDM cosmological model and derived from the combination of growth of structure and geometry. When combined with primarily geometrical constraints from Planck CMB, H0, BAO and SNe, adding the SPT cluster data improves the w constraint from the geometrical data alone by 14%, to w = -1.023+-0.042.

Kelly (2007, hereafter K07) described an efficient algorithm, using Gibbs sampling, for performing linear regression in the fairly general case where non-zero measurement errors exist for both the covariates and response variables, where these measurements may be correlated (for the same data point), where the response variable is affected by intrinsic scatter in addition to measurement error, and where the prior distribution of covariates is modeled by a flexible mixture of Gaussians rather than assumed to be uniform. Here I extend the K07 algorithm in two ways. First, the procedure is generalized to the case of multiple response variables. Second, I describe how to model the prior distribution of covariates using a Dirichlet process, which can be thought of as a Gaussian mixture where the number of mixture components is learned from the data. I present an example of multivariate regression using the extended algorithm, namely fitting scaling relations of the gas mass, temperature, and luminosity of dynamically relaxed galaxy clusters as a function of their mass and redshift. An implementation of the Gibbs sampler in the R language, called LRGS, is provided.

This is the third in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Our sample comprises 40 clusters identified as being dynamically relaxed and hot (i.e., massive) in Papers I and II of this series. Here we consider the thermodynamics of the intracluster medium, in particular the profiles of density, temperature and related quantities, as well as integrated measurements of gas mass, average temperature, total luminosity and center-excluded luminosity. We fit power-law scaling relations of each of these quantities as a function of redshift and cluster mass, which can be measured precisely and with minimal bias for these relaxed clusters using hydrostatic arguments. For the thermodynamic profiles, we jointly model the density and temperature and their intrinsic scatter as a function of radius, thus also capturing the behavior of the gas pressure and entropy. For the integrated quantities, we also jointly fit a multidimensional intrinsic covariance. Our results reinforce the view that simple hydrodynamical models provide a good description of relaxed clusters outside their centers, but that additional heating and cooling processes are important in the inner regions (radii $r < 0.5r_{2500} \approx 0.15r_{500}$). The thermodynamic profiles remain regular, with small intrinsic scatter, down to the smallest radii where deprojection is straightforward ($\sim 20$ kpc); within this radius, even the most relaxed systems show clear departures from spherical symmetry. Our results suggest that heating and cooling are continuously regulated in a tight feedback loop, allowing the cluster atmosphere to remain stratified on these scales.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

Traditionally, galaxy clusters have been expected to retain all the material accreted since their formation epoch. For this reason, their matter content should be representative of the Universe as a whole, and thus their baryon fraction should be close to the Universal baryon fraction. We make use of the sample of the 100 brightest galaxy clusters discovered in the XXL Survey to investigate the fraction of baryons in the form of hot gas and stars in the cluster population. We measure the gas masses of the detected halos and use a mass--temperature relation directly calibrated using weak-lensing measurements for a subset of XXL clusters to estimate the halo mass. We find that the weak-lensing calibrated gas fraction of XXL-100-GC clusters is substantially lower than was found in previous studies using hydrostatic masses. Our best-fit relation between gas fraction and mass reads $f_{\rm gas,500}=0.055_{-0.006}^{+0.007}\left(M_{\rm 500}/10^{14}M_\odot\right)^{0.21_{-0.10}^{+0.11}}$. The baryon budget of galaxy clusters therefore falls short of the Universal baryon fraction by about a factor of two at $r_{\rm 500}$. Our measurements require a hydrostatic bias $1-b=M_X/M_{\rm WL}=0.72_{-0.07}^{+0.08}$ to match the gas fraction obtained using lensing and hydrostatic equilibrium. Comparing our gas fraction measurements with the expectations from numerical simulations, our results favour an extreme feedback scheme in which a significant fraction of the baryons are expelled from the cores of halos. This model is, however, in contrast with the thermodynamical properties of observed halos, which might suggest that weak-lensing masses are overestimated. We note that a mass bias $1-b=0.58$ as required to reconcile Planck CMB and cluster counts should translate into an even lower baryon fraction, which poses a major challenge to our current understanding of galaxy clusters. [Abridged]

We present confirmation of the cluster MOO J1142+1527, a massive galaxy cluster discovered as part of the Massive and Distant Clusters of WISE Survey. The cluster is confirmed to lie at $z=1.19$, and using the Combined Array for Research in Millimeter-wave Astronomy we robustly detect the Sunyaev-Zel'dovich (SZ) decrement at 13.2$\sigma$. The SZ data imply a mass of $\mathrm{M}_{200m}=(1.1\pm0.2)\times10^{15}$ $\mathrm{M}_\odot$, making MOO J1142+1527 the most massive galaxy cluster known at $z>1.15$ and the second most massive cluster known at $z>1$. For a standard $\Lambda$CDM cosmology it is further expected to be one of the $\sim 5$ most massive clusters expected to exist at $z\ge1.19$ over the entire sky. Our ongoing Spitzer program targeting $\sim1750$ additional candidate clusters will identify comparably rich galaxy clusters over the full extragalactic sky.

This is the fourth in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here, we use measurements of weak gravitational lensing from the Weighing the Giants project to calibrate Chandra X-ray measurements of total mass that rely on the assumption of hydrostatic equilibrium. This comparison of X-ray and lensing masses provides a measurement of the combined bias of X-ray hydrostatic masses due to both astrophysical and instrumental sources. Assuming a fixed cosmology, and within a characteristic radius (r_2500) determined from the X-ray data, we measure a lensing to X-ray mass ratio of 0.96 +/- 9% (stat) +/- 9% (sys). We find no significant trends of this ratio with mass, redshift or the morphological indicators used to select the sample. In accordance with predictions from hydro simulations for the most massive, relaxed clusters, our results disfavor strong, tens-of-percent departures from hydrostatic equilibrium at these radii. In addition, we find a mean concentration of the sample measured from lensing data of c_200 = $3.0_{-1.8}^{+4.4}$. Anticipated short-term improvements in lensing systematics, and a modest expansion of the relaxed lensing sample, can easily increase the measurement precision by 30--50%, leading to similar improvements in cosmological constraints that employ X-ray hydrostatic mass estimates, such as on Omega_m from the cluster gas mass fraction.

We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344-4243). Deep optical imaging reveals previously-undetected filaments of star formation, extending to radii of ~50-100 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy reveals a massive (2x10^9 Msun)), young (~4.5 Myr) population of stars, consistent with a time-averaged star formation rate of 610 +/- 50 Msun/yr. We report a strong detection of OVI(1032,1038) which appears to originate primarily in shock-heated gas, but may contain a substantial contribution (>1000 Msun/yr) from the cooling intracluster medium. We confirm the presence of deep X-ray cavities in the inner ~10 kpc, which are amongst the most extreme examples of radio-mode feedback detected to date, implying jet powers of 2-7 x10^45 erg/s. We provide evidence that the AGN inflating these cavities may have only recently transitioned from "quasar-mode" to "radio-mode", and may currently be insufficient to completely offset cooling. A model-subtracted residual X-ray image reveals evidence for prior episodes of strong radio-mode feedback at radii of ~100 kpc, with extended "ghost" cavities indicating a prior epoch of feedback roughly 100 Myr ago. This residual image also exhibits significant asymmetry in the inner ~200 kpc (0.15R500), reminiscent of infalling cool clouds, either due to minor mergers or fragmentation of the cooling ICM. Taken together, these data reveal a rapidly evolving cool core which is rich with structure (both spatially and in temperature), is subject to a variety of highly energetic processes, and yet is cooling rapidly and forming stars along thin, narrow filaments.

The abundance of massive galaxy clusters is a powerful probe of departures from General Relativity (GR) on cosmic scales. Despite current stringent constraints placed by stellar and galactic tests, on larger scales alternative theories of gravity such as $f(R)$ can still work as effective theories. Here we present constraints on two popular models of $f(R)$, Hu-Sawicki and "designer", derived from a fully self-consistent analysis of current samples of X-ray selected clusters and accounting for all the covariances between cosmological and astrophysical parameters. Using cluster number counts in combination with recent data from the cosmic microwave background (CMB) and the CMB lensing potential generated by large scale structures, as well as with other cosmological constraints on the background expansion history and its mean matter density, we obtain the upper bounds $\log_{10}|f_{R0}| < 4.79$ and $\log_{10}B_0 < 3.75$ at the 95.4 per cent confidence level, for the Hu-Sawicki (with $n=1$) and designer models, respectively. The robustness of our results derives from high quality cluster growth data for the most massive clusters known out to redshifts $z \sim 0.5$, a tight control of systematic uncertainties including an accurate and precise mass calibration from weak gravitational lensing data, and the use of the full shape of the halo mass function over the mass range of our data.

Clusters of galaxies are expected to gravitationally lens the cosmic microwave background (CMB) and thereby generate a distinct signal in the CMB on arcminute scales. Measurements of this effect can be used to constrain the masses of galaxy clusters with CMB data alone. Here we present a measurement of lensing of the CMB by galaxy clusters using data from the South Pole Telescope (SPT). We develop a maximum likelihood approach to extract the CMB cluster lensing signal and validate the method on mock data. We quantify the effects on our analysis of several potential sources of systematic error and find that they generally act to reduce the best-fit cluster mass. It is estimated that this bias to lower cluster mass is roughly $0.85\sigma$ in units of the statistical error bar, although this estimate should be viewed as an upper limit. We apply our maximum likelihood technique to 513 clusters selected via their SZ signatures in SPT data, and rule out the null hypothesis of no lensing at $3.1\sigma$. The lensing-derived mass estimate for the full cluster sample is consistent with that inferred from the SZ flux: $M_{200,\mathrm{lens}} = 0.83_{-0.37}^{+0.38}\, M_{200,\mathrm{SZ}}$ (68% C.L., statistical error only).

We present deep, 170 ks, Chandra X-ray observations of Abell 2219 (z=0.23) one of the hottest and most X-ray luminous clusters known, and which is experiencing a major merger event. We discover a 'horseshoe' of high temperature gas surrounding the ram-pressure-stripped, bright, hot, X-ray cores. We confirm an X-ray shock front located north-west of the X-ray centroid and along the projected merger axis. We also find a second shock front to the south-east of the X-ray centroid making this only the second cluster where both the shock and reverse shock are confirmed with X-ray temperature measurements. We also present evidence for a sloshing cold front in the 'remnant tail' of one of the sub-cluster cores. The cold front and north-west shock front geometrically bound the radio halo and appear to be directly influencing the radio properties of the cluster.

We present scaling relations between the integrated Sunyaev-Zel'dovich Effect (SZE) signal, $Y_{\rm SZ}$, its X-ray analogue, $Y_{\rm X}\equiv M_{\rm gas}T_{\rm X}$, and total mass, $M_{\rm tot}$, for the 45 galaxy clusters in the Bolocam X-ray-SZ (BOXSZ) sample. All parameters are integrated within $r_{2500}$. $Y_{2500}$ values are measured using SZE data collected with Bolocam, operating at 140 GHz at the Caltech Submillimeter Observatory (CSO). The temperature, $T_{\rm X}$, and mass, $M_{\rm gas,2500}$, of the intracluster medium are determined using X-ray data collected with Chandra, and $M_{\rm tot}$ is derived from $M_{\rm gas}$ assuming a constant gas mass fraction. Our analysis accounts for several potential sources of bias, including: selection effects, contamination from radio point sources, and the loss of SZE signal due to noise filtering and beam-smoothing effects. We measure the $Y_{2500}$--$Y_{\rm X}$ scaling to have a power-law index of $0.84\pm0.07$, and a fractional intrinsic scatter in $Y_{2500}$ of $(21\pm7)\%$ at fixed $Y_{\rm X}$, both of which are consistent with previous analyses. We also measure the scaling between $Y_{2500}$ and $M_{2500}$, finding a power-law index of $1.06\pm0.12$ and a fractional intrinsic scatter in $Y_{2500}$ at fixed mass of $(25\pm9)\%$. While recent SZE scaling relations using X-ray mass proxies have found power-law indices consistent with the self-similar prediction of 5/3, our measurement stands apart by differing from the self-similar prediction by approximately 5$\sigma$. Given the good agreement between the measured $Y_{2500}$--$Y_{\rm X}$ scalings, much of this discrepancy appears to be caused by differences in the calibration of the X-ray mass proxies adopted for each particular analysis.

We present CARMA 30 GHz Sunyaev-Zel'dovich (SZ) observations of five high-redshift ($z \gtrsim 1$), infrared-selected galaxy clusters discovered as part of the all-sky Massive and Distant Clusters of WISE Survey (MaDCoWS). The SZ decrements measured toward these clusters demonstrate that the MaDCoWS selection is discovering evolved, massive galaxy clusters with hot intracluster gas. Using the SZ scaling relation calibrated with South Pole Telescope clusters at similar masses and redshifts, we find these MaDCoWS clusters have masses in the range $M_{200} \approx 2-6 \times 10^{14}$ $M_\odot$. Three of these are among the most massive clusters found to date at $z\gtrsim 1$, demonstrating that MaDCoWS is sensitive to the most massive clusters to at least $z = 1.3$. The added depth of the AllWISE data release will allow all-sky infrared cluster detection to $z \approx 1.5$ and beyond.

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 data. 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 present a catalog of galaxy clusters selected via their Sunyaev-Zel'dovich (SZ) effect signature from 2500 deg$^2$ of South Pole Telescope (SPT) data. This work represents the complete sample of clusters detected at high significance in the 2500-square-degree SPT-SZ survey, which was completed in 2011. A total of 677 (409) cluster candidates are identified above a signal-to-noise threshold of $\xi$ =4.5 (5.0). Ground- and space-based optical and near-infrared (NIR) imaging confirms overdensities of similarly colored galaxies in the direction of 516 (or 76%) of the $\xi$>4.5 candidates and 387 (or 95%) of the $\xi$>5 candidates; the measured purity is consistent with expectations from simulations. Of these confirmed clusters, 415 were first identified in SPT data, including 251 new discoveries reported in this work. We estimate photometric redshifts for all candidates with identified optical and/or NIR counterparts; we additionally report redshifts derived from spectroscopic observations for 141 of these systems. The mass threshold of the catalog is roughly independent of redshift above $z$~0.25 leading to a sample of massive clusters that extends to high redshift. The median mass of the sample is $M_{\scriptsize 500c}(\rho_\mathrm{crit})$ ~ 3.5 x 10$^{14} M_\odot h^{-1}$, the median redshift is $z_{med}$ =0.55, and the highest-redshift systems are at $z$>1.4. The combination of large redshift extent, clean selection, and high typical mass makes this cluster sample of particular interest for cosmological analyses and studies of cluster formation and evolution.

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 analysis 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.81\pm0.03$, with marginalized, one-dimensional constraints of $\Omega_m=0.26\pm0.03$ and $\sigma_8=0.83\pm0.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.

In light of the tension in cosmological constraints reported by the Planck team between their SZ-selected cluster counts and Cosmic Microwave Background (CMB) temperature anisotropies, we compare the Planck cluster mass estimates with robust, weak-lensing mass measurements from the Weighing the Giants (WtG) project. For the 22 clusters in common between the Planck cosmology sample and WtG, we find an overall mass ratio of $\left< M_{Planck}/M_{\rm WtG} \right> = 0.688 \pm 0.072$. Extending the sample to clusters not used in the Planck cosmology analysis yields a consistent value of $\left< M_{Planck}/M_{\rm WtG} \right> = 0.698 \pm 0.062$ from 38 clusters in common. Identifying the weak-lensing masses as proxies for the true cluster mass (on average), these ratios are $\sim 1.6\sigma$ lower than the default mass bias of 0.8 assumed in the Planck cluster analysis. Adopting the WtG weak-lensing-based mass calibration would substantially reduce the tension found between the Planck cluster count cosmology results and those from CMB temperature anisotropies, thereby dispensing of the need for "new physics" such as uncomfortably large neutrino masses (in the context of the measured Planck temperature anisotropies and other data). We also find modest evidence (at 95 per cent confidence) for a mass dependence of the calibration ratio and discuss its potential origin in light of systematic uncertainties in the temperature calibration of the X-ray measurements used to calibrate the Planck cluster masses. Our results exemplify the critical role that robust absolute mass calibration plays in cluster cosmology, and the invaluable role of accurate weak-lensing mass measurements in this regard.

We report the detection of the Sunyaev-Zel'dovich (SZ) effect of galaxy cluster XLSSU J021744.1-034536, using 30 GHz CARMA data. This cluster was discovered via its extended X-ray emission in the XMM-Newton Large Scale Structure survey, the precursor to the XXL survey. It has a photometrically determined redshift $z=1.91^{+0.19}_{-0.21}$, making it among the most distant clusters known, and nominally the most distant for which the SZ effect has been measured. The spherically integrated Comptonization is $Y_{500}=(3.0\pm0.4)\times 10^{-12}$, a measurement which is relatively insensitive to assumptions regarding the size and redshift of the cluster, as well as the background cosmology. Using a variety of locally calibrated cluster scaling relations extrapolated to z~2, we estimate a mass $M_{500} \sim (1$-$2)\times 10^{14}M_{sun}$ from the X-ray flux and SZ signal. The measured properties of this cluster are in good agreement with the extrapolation of an X-ray luminosity-SZ effect scaling relation calibrated from clusters discovered by the South Pole Telescope at higher masses and lower redshifts. The full XXL-CARMA sample will provide a more complete, multi-wavelength census of distant clusters in order to robustly extend the calibration of cluster scaling relations to these high redshifts.

SPT-CLJ2040-4451 -- spectroscopically confirmed at z = 1.478 -- is the highest redshift galaxy cluster yet discovered via the Sunyaev-Zel'dovich effect. SPT-CLJ2040-4451 was a candidate galaxy cluster identified in the first 720 deg^2 of the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey, and confirmed in follow-up imaging and spectroscopy. From multi-object spectroscopy with Magellan-I/Baade+IMACS we measure spectroscopic redshifts for 15 cluster member galaxies, all of which have strong [O II] 3727 emission. SPT-CLJ2040-4451 has an SZ-measured mass of M_500,SZ = 3.2 +/- 0.8 X 10^14 M_Sun/h_70, corresponding to M_200,SZ = 5.8 +/- 1.4 X 10^14 M_Sun/h_70. The velocity dispersion measured entirely from blue star forming members is sigma_v = 1500 +/- 520 km/s. The prevalence of star forming cluster members (galaxies with > 1.5 M_Sun/yr) implies that this massive, high-redshift cluster is experiencing a phase of active star formation, and supports recent results showing a marked increase in star formation occurring in galaxy clusters at z >1.4. We also compute the probability of finding a cluster as rare as this in the SPT-SZ survey to be >99%, indicating that its discovery is not in tension with the concordance Lambda-CDM cosmological model.

We present the results of a new analysis of the X-ray selected Active Galactic Nuclei (AGN) population in the vicinity of 135 of the most massive galaxy clusters in the redshift range of 0.2 < z < 0.9 observed with Chandra. With a sample of more than 11,000 X-ray point sources, we are able to measure, for the first time, evidence for evolution in the cluster AGN population beyond the expected evolution of field AGN. Our analysis shows that overall number density of cluster AGN scales with the cluster mass as $\sim M_{500}^{-1.2}$. There is no evidence for the overall number density of cluster member X-ray AGN depending on the cluster redshift in a manner different than field AGN, nor there is any evidence that the spatial distribution of cluster AGN (given in units of the cluster overdensity radius r_500) strongly depends on the cluster mass or redshift. The $M^{-1.2 \pm 0.7}$ scaling relation we measure is consistent with theoretical predictions of the galaxy merger rate in clusters, which is expected to scale with the cluster velocity dispersion, $\sigma$, as $ \sim \sigma^{-3}$ or $\sim M^{-1}$. This consistency suggests that AGN in clusters may be predominantly triggered by galaxy mergers, a result that is further corroborated by visual inspection of Hubble images for 23 spectroscopically confirmed cluster member AGN in our sample. A merger-driven scenario for the triggering of X-ray AGN is not strongly favored by studies of field galaxies, however, suggesting that different mechanisms may be primarily responsible for the triggering of cluster and field X-ray AGN.

We present optical spectroscopy of galaxies in clusters detected through the Sunyaev-Zel'dovich (SZ) effect with the South Pole Telescope (SPT). We report our own measurements of $61$ spectroscopic cluster redshifts, and $48$ velocity dispersions each calculated with more than $15$ member galaxies. This catalog also includes $19$ dispersions of SPT-observed clusters previously reported in the literature. The majority of the clusters in this paper are SPT-discovered; of these, most have been previously reported in other SPT cluster catalogs, and five are reported here as SPT discoveries for the first time. By performing a resampling analysis of galaxy velocities, we find that unbiased velocity dispersions can be obtained from a relatively small number of member galaxies ($\lesssim 30$), but with increased systematic scatter. We use this analysis to determine statistical confidence intervals that include the effect of membership selection. We fit scaling relations between the observed cluster velocity dispersions and mass estimates from SZ and X-ray observables. In both cases, the results are consistent with the scaling relation between velocity dispersion and mass expected from dark-matter simulations. We measure a $\sim$30% log-normal scatter in dispersion at fixed mass, and a $\sim$10% offset in the normalization of the dispersion-mass relation when compared to the expectation from simulations, which is within the expected level of systematic uncertainty.

(Abridged) We present the results of an X-ray analysis of 80 galaxy clusters selected in the 2500 deg^2 South Pole Telescope survey and observed with the Chandra X-ray Observatory. We divide the full sample into subsamples of ~20 clusters based on redshift and central density, performing an X-ray fit to all clusters in a subsample simultaneously, assuming self-similarity of the temperature profile. This approach allows us to constrain the shape of the temperature profile over 0<r<1.5R500, which would be impossible on a per-cluster basis, since the observations of individual clusters have, on average, 2000 X-ray counts. The results presented here represent the first constraints on the evolution of the average temperature profile from z=0 to z=1.2. We find that high-z (0.6<z<1.2) clusters are slightly (~40%) cooler both in the inner (r<0.1R500) and outer (r>R500) regions than their low-z (0.3<z<0.6) counterparts. Combining the average temperature profile with measured gas density profiles from our earlier work, we infer the average pressure and entropy profiles for each subsample. Overall, our observed pressure profiles agree well with earlier lower-redshift measurements, suggesting minimal redshift evolution in the pressure profile outside of the core. We find no measurable redshift evolution in the entropy profile at r<0.7R500. We observe a slight flattening of the entropy profile at r>R500 in our high-z subsample. This flattening is consistent with a temperature bias due to the enhanced (~3x) rate at which group-mass (~2 keV) halos, which would go undetected at our survey depth, are accreting onto the cluster at z~1. This work demonstrates a powerful method for inferring spatially-resolved cluster properties in the case where individual cluster signal-to-noise is low, but the number of observed clusters is high.

We report weak-lensing masses for 51 of the most X-ray luminous galaxy clusters known. This cluster sample, introduced earlier in this series of papers, spans redshifts 0.15 < z_cl < 0.7, and is well suited to calibrate mass proxies for current cluster cosmology experiments. Cluster masses are measured with a standard `color-cut' lensing method from three-filter photometry of each field. Additionally, for 27 cluster fields with at least five-filter photometry, we measure high-accuracy masses using a new method that exploits all information available in the photometric redshift posterior probability distributions of individual galaxies. Using simulations based on the COSMOS-30 catalog, we demonstrate control of systematic biases in the mean mass of the sample with this method, from photometric redshift biases and associated uncertainties, to better than 3%. In contrast, we show that the use of single-point estimators in place of the full photometric redshift posterior distributions can lead to significant redshift-dependent biases on cluster masses. The performance of our new photometric redshift-based method allows us to calibrate `color-cut` masses for all 51 clusters in the present sample to a total systematic uncertainty of ~7% on the mean mass, a level sufficient to significantly improve current cosmology constraints from galaxy clusters. Our results bode well for future cosmological studies of clusters, potentially reducing the need for exhaustive spectroscopic calibration surveys as compared to other techniques, when deep, multi-filter optical and near-IR imaging surveys are coupled with robust photometric redshift methods.

We present improved methods for using stars found in astronomical exposures to calibrate both star and galaxy colors as well as to adjust the instrument flat field. By developing a spectroscopic model for the SDSS stellar locus in color-color space, synthesizing an expected stellar locus, and simultaneously solving for all unknown zeropoints when fitting to the instrumental locus, we increase the calibration accuracy of stellar locus matching. We also use a new combined technique to estimate improved flat-field models for the Subaru SuprimeCam camera, forming `star flats' based on the magnitudes of stars observed in multiple positions or through comparison with available SDSS magnitudes. These techniques yield galaxy magnitudes with reliable color calibration (< 0.01 - 0.02 mag accuracy) that enable us to estimate photometric redshift probability distributions without spectroscopic training samples. We test the accuracy of our photometric redshifts using spectroscopic redshifts z_s for ~5000 galaxies in 27 cluster fields with at least five bands of photometry, as well as galaxies in the COSMOS field, finding sigma((z_p - z_s)/(1 + z_s)) ~ 0.03 for the most probable redshift z_p. We show that the full posterior probability distributions for the redshifts of galaxies with five-band photometry exhibit good agreement with redshifts estimated from thirty-band photometry in the COSMOS field. The growth of shear with increasing distance behind each galaxy cluster shows the expected redshift-distance relation for a flat Lambda-CDM cosmology. Photometric redshifts and calibrated colors are used in subsequent papers to measure the masses of 51 galaxy clusters from their weak gravitational shear. We make our Python code for stellar locus matching available at http://big-macs-calibrate.googlecode.com; the code requires only a catalog and filter functions.

This is the first in a series of papers in which we measure accurate weak-lensing masses for 51 of the most X-ray luminous galaxy clusters known at redshifts 0.15<z<0.7, in order to calibrate X-ray and other mass proxies for cosmological cluster experiments. The primary aim is to improve the absolute mass calibration of cluster observables, currently the dominant systematic uncertainty for cluster count experiments. Key elements of this work are the rigorous quantification of systematic uncertainties, high-quality data reduction and photometric calibration, and the "blind" nature of the analysis to avoid confirmation bias. Our target clusters are drawn from RASS X-ray catalogs, and provide a versatile calibration sample for many aspects of cluster cosmology. We have acquired wide-field, high-quality imaging using the Subaru and CFHT telescopes for all 51 clusters, in at least three bands per cluster. For a subset of 27 clusters, we have data in at least five bands, allowing accurate photo-z estimates of lensed galaxies. In this paper, we describe the cluster sample and observations, and detail the processing of the SuprimeCam data to yield high-quality images suitable for robust weak-lensing shape measurements and precision photometry. For each cluster, we present wide-field color optical images and maps of the weak-lensing mass distribution, the optical light distribution, and the X-ray emission, providing insights into the large-scale structure in which the clusters are embedded. We measure the offsets between X-ray centroids and Brightest Cluster Galaxies in the clusters, finding these to be small in general, with a median of 20kpc. For offsets <100kpc, weak-lensing mass measurements centered on the BCGs agree well with values determined relative to the X-ray centroids; miscentering is therefore not a significant source of systematic uncertainty for our mass measurements. [abridged]

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 the 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.4\pm2.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.089\pm0.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.27\pm0.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.98\pm0.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)

We report our analysis of MACS J0717.5+3745 using 140 and 268 GHz Bolocam data collected at the Caltech Submillimeter Observatory. We detect extended Sunyaev-Zel'dovich (SZ) effect signal at high significance in both Bolocam bands, and we employ Herschel-SPIRE observations to subtract the signal from dusty background galaxies in the 268 GHz data. We constrain the two-band SZ surface brightness toward two of the sub-clusters of MACS J0717.5+3745: the main sub-cluster (named C), and a sub-cluster identified in spectroscopic optical data to have a line-of-sight velocity of +3200 km/s (named B). We determine the surface brightness in two separate ways: via fits of parametric models and via direct integration of the images. For both sub-clusters, we find consistent surface brightnesses from both analysis methods. We constrain spectral templates consisting of relativistically corrected thermal and kinetic SZ signals, using a jointly-derived electron temperature from Chandra and XMM-Newton under the assumption that each sub-cluster is isothermal. The data show no evidence for a kinetic SZ signal toward sub-cluster C, but they do indicate a significant kinetic SZ signal toward sub-cluster B. The model-derived surface brightnesses for sub-cluster B yield a best-fit line-of-sight velocity of v_z = +3450 +- 900 km/s, with (1 - Prob[v_z > 0]) = 1.3 x 10^-5 (4.2 sigma away from 0 for a Gaussian distribution). The directly integrated sub-cluster B SZ surface brightnesses provide a best-fit v_z = +2550 +- 1050 km/s, with (1 - Prob[v_z > 0]) = 2.2 x 10^-3 (2.9 sigma).

We describe a method for measuring the integrated Comptonization (YSZ) of clusters of galaxies from measurements of the Sunyaev-Zel'dovich (SZ) effect in multiple frequency bands and use this method to characterize a sample of galaxy clusters detected in South Pole Telescope (SPT) data. We test this method on simulated cluster observations and verify that it can accurately recover cluster parameters with negligible bias. In realistic simulations of an SPT-like survey, with realizations of cosmic microwave background anisotropy, point sources, and atmosphere and instrumental noise at typical SPT-SZ survey levels, we find that YSZ is most accurately determined in an aperture comparable to the SPT beam size. We demonstrate the utility of this method to measure YSZ and to constrain mass scaling relations using X-ray mass estimates for a sample of 18 galaxy clusters from the SPT-SZ survey. Measuring YSZ within a 0.75' radius aperture, we find an intrinsic log-normal scatter of 21+/-11% in YSZ at a fixed mass. Measuring YSZ within a 0.3 Mpc projected radius (equivalent to 0.75' at the survey median redshift z = 0.6), we find a scatter of 26+/-9%. Prior to this study, the SPT observable found to have the lowest scatter with mass was cluster detection significance. We demonstrate, from both simulations and SPT observed clusters, that YSZ measured within an aperture comparable to the SPT beam size is equivalent, in terms of scatter with cluster mass, to SPT cluster detection significance.

The adiabatic evolution of the temperature of the cosmic microwave background (CMB) is a key prediction of standard cosmology. We study deviations from the expected adiabatic evolution of the CMB temperature of the form $T(z) =T_0(1+z)^{1-\alpha}$ using measurements of the spectrum of the Sunyaev Zel'dovich Effect with the South Pole Telescope (SPT). We present a method for using the ratio of the Sunyaev Zel'dovich signal measured at 95 and 150 GHz in the SPT data to constrain the temperature of the CMB. We demonstrate that this approach provides unbiased results using mock observations of clusters from a new set of hydrodynamical simulations. We apply this method to a sample of 158 SPT-selected clusters, spanning the redshift range $0.05 < z < 1.35$, and measure $\alpha = 0.017^{+0.030}_{-0.028}$, consistent with the standard model prediction of $\alpha=0$. In combination with other published results, we constrain $\alpha = 0.011 \pm 0.016$, an improvement of $\sim 20\%$ over published constraints. This measurement also provides a strong constraint on the effective equation of state in models of decaying dark energy $w_\mathrm{eff} = -0.987^{+0.016}_{-0.017}$.

We present a measurement of the fraction of cluster galaxies hosting X-ray bright Active Galactic Nuclei (AGN) as a function of clustercentric distance scaled in units of $r_{500}$. Our analysis employs high quality Chandra X-ray and Subaru optical imaging for 42 massive X-ray selected galaxy cluster fields spanning the redshift range of $0.2 < z < 0.7$. In total, our study involves 176 AGN with bright ($R <23$) optical counterparts above a $0.5-8.0$ keV flux limit of $10^{-14} \rm{erg} \ \rm{cm}^{-2} \ \rm{s}^{-1}$. When excluding central dominant galaxies from the calculation, we measure a cluster-galaxy AGN fraction in the central regions of the clusters that is $\sim 3 $ times lower that the field value. This fraction increases with clustercentric distance before becoming consistent with the field at $\sim 2.5 r_{500}$. Our data exhibit similar radial trends to those observed for star formation and optically selected AGN in cluster member galaxies, both of which are also suppressed near cluster centers to a comparable extent. These results strongly support the idea that X-ray AGN activity and strong star formation are linked through their common dependence on available reservoirs of cold gas.

We present first results on the cooling properties derived from Chandra X-ray observations of 83 high-redshift (0.3 < z < 1.2) massive galaxy clusters selected by their Sunyaev-Zel'dovich signature in the South Pole Telescope data. We measure each cluster's central cooling time, central entropy, and mass deposition rate, and compare to local cluster samples. We find no significant evolution from z~0 to z~1 in the distribution of these properties, suggesting that cooling in cluster cores is stable over long periods of time. We also find that the average cool core entropy profile in the inner ~100 kpc has not changed dramatically since z ~ 1, implying that feedback must be providing nearly constant energy injection to maintain the observed "entropy floor" at ~10 keV cm^2. While the cooling properties appear roughly constant over long periods of time, we observe strong evolution in the gas density profile, with the normalized central density (rho_0/rho_crit) increasing by an order of magnitude from z ~ 1 to z ~ 0. When using metrics defined by the inner surface brightness profile of clusters, we find an apparent lack of classical, cuspy, cool-core clusters at z > 0.75, consistent with earlier reports for clusters at z > 0.5 using similar definitions. Our measurements indicate that cool cores have been steadily growing over the 8 Gyr spanned by our sample, consistent with a constant, ~150 Msun/yr cooling flow that is unable to cool below entropies of 10 keV cm^2 and, instead, accumulates in the cluster center. We estimate that cool cores began to assemble in these massive systems at z ~ 1, which represents the first constraints on the onset of cooling in galaxy cluster cores. We investigate several potential biases which could conspire to mimic this cool core evolution and are unable to find a bias that has a similar redshift dependence and a substantial amplitude.

In addition to cosmological tests based on the mass function and clustering of galaxy clusters, which probe the growth of cosmic structure, nature offers two independent ways of using clusters to measure cosmic distances. The first uses measurements of the X-ray emitting gas mass fraction, which is an approximately standard quantity, independent of mass and redshift, for the most massive clusters. The second uses combined millimeter (mm) and X-ray measurements of cluster pressure profiles. We review these methods, their current status and the prospects for improvements over the next decade. For the first technique, which currently provides comparable dark energy constraints to type Ia supernova studies, improvements of a factor of 6 or more should be readily achievable, together with tight constraints on the mean matter density that are largely independent of the cosmological model assumed. Realizing this potential will require a coordinated, multiwavelength approach, utilizing new cluster surveys, X-ray, optical and mm facilities, and a continued emphasis on improved hydrodynamical simulations.

We present the results from extensive, new observations of the Perseus Cluster of galaxies, obtained as a Suzaku Key Project. The 85 pointings analyzed span eight azimuthal directions out to 2 degrees = 2.6 Mpc, to and beyond the virial radius r_200 ~ 1.8 Mpc, offering the most detailed X-ray observation of the intracluster medium (ICM) at large radii in any cluster to date. The azimuthally averaged density profile for r>0.4r_200 is relatively flat, with a best-fit power-law index of 1.69+/-0.13 significantly smaller than expected from numerical simulations. The entropy profile in the outskirts lies systematically below the power-law behavior expected from large-scale structure formation models which include only the heating associated with gravitational collapse. The pressure profile beyond ~0.6r_200 shows an excess with respect to the best-fit model describing the SZ measurements for a sample of clusters observed with Planck. The inconsistency between the expected and measured density, entropy, and pressure profiles can be explained primarily by an overestimation of the density due to inhomogeneous gas distribution in the outskirts; there is no evidence for a bias in the temperature measurements within the virial radius. We find significant differences in thermodynamic properties of the ICM at large radii along the different arms. Along the cluster minor axis, we find a flattening of the entropy profiles outside ~0.6r_200, while along the major axis, the entropy rises all the way to the outskirts. Correspondingly, the inferred gas clumping factor is typically larger along the minor than along the major axis.

We present results from a large mosaic of Suzaku observations of the Coma Cluster, the nearest and X-ray brightest hot, dynamically active, non-cool core system, focusing on the thermodynamic properties of the ICM on large scales. For azimuths not aligned with an infalling subcluster towards the southwest, our measured temperature and X-ray brightness profiles exhibit broadly consistent radial trends, with the temperature decreasing from about 8.5 keV at the cluster center to about 2 keV at a radius of 2 Mpc, which is the edge of our detection limit. The SW merger significantly boosts the surface brightness, allowing us to detect X-ray emission out to ~2.2 Mpc along this direction. Apart from the southwestern infalling subcluster, the surface brightness profiles show multiple edges around radii of 30-40 arcmin. The azimuthally averaged temperature profile, as well as the deprojected density and pressure profiles, all show a sharp drop consistent with an outward propagating shock front located at 40 arcmin, corresponding to the outermost edge of the giant radio halo observed at 352 MHz with the WSRT. The shock front may be powering this radio emission. A clear entropy excess inside of r_500 reflects the violent merging events linked with these morphological features. Beyond r_500, the entropy profiles of the Coma Cluster along the relatively relaxed directions are consistent with the power-law behavior expected from simple models of gravitational large-scale structure formation. The pressure is also in agreement at these radii with the expected values measured from SZ data from the Planck satellite. However, due to the large uncertainties associated with the Coma Cluster measurements, we cannot yet exclude an entropy flattening in this system consistent with that seen in more relaxed cool core clusters.

We describe Sunyaev-Zel'dovich (SZ) effect measurements and analysis of the intracluster medium (ICM) pressure profiles of a set of 45 massive galaxy clusters imaged using Bolocam at the Caltech Submillimeter Observatory. We have used masses determined from Chandra X-ray observations to scale each cluster's profile by the overdensity radius R500 and the mass-and-redshift-dependent normalization factor P500. We deproject the average pressure profile of our sample into 13 logarithmically spaced radial bins between 0.07R500 and 3.5R500. We find that a generalized Navarro, Frenk, and White (gNFW) profile describes our data with sufficient goodness-of-fit and best-fit parameters (C500, alpha, beta, gamma, P0 = 1.18, 0.86, 3.67, 0.67, 4.29). We also use the X-ray data to define cool-core and disturbed subsamples of clusters, and we constrain the average pressure profiles of each of these subsamples. We find that given the precision of our data the average pressure profiles of disturbed and cool-core clusters are consistent with one another at R>~0.15R500, with cool-core systems showing indications of higher pressure at R<~0.15R500. In addition, for the first time, we place simultaneous constraints on the mass scaling of cluster pressure profiles, their ensemble mean profile, and their radius-dependent intrinsic scatter between 0.1R500 and 2.0R500. The scatter among profiles is minimized at radii between ~0.2R500 and ~0.5R500, with a value of ~20%. The best-fit mass scaling has a power-law slope of 0.49, which is shallower than the nominal prediction of 2/3 from self-similar hydrostatic equilibrium models. These results for the intrinsic scatter and mass scaling are largely consistent with previous analyses, most of which have relied heavily on X-ray derived pressures of clusters at significantly lower masses and redshifts compared to our sample.

We report constraints on primordial non-Gaussianity from the abundance of X-ray detected clusters. Our analytic prescription for adding non-Gaussianity to the cluster mass function takes into account moments beyond the skewness, and we demonstrate that those moments should not be ignored in most analyses of cluster data. We constrain the amplitude of the skewness for two scenarios that have different overall levels of non-Gaussianity, characterized by how amplitudes of higher cumulants scale with the skewness. We find that current data can constrain these one-parameter non-Gaussian models at a useful level, but are not sensitive to adding further details of the corresponding inflation scenarios. Combining cluster data with Cosmic Microwave Background constraints on the cosmology and power spectrum amplitude, we find the dimensionless skewness to be 1000*M3=-1+24-28 for one of our scaling scenarios, and 1000*M3=-4+/-7 for the other. These are the first constraints on non-Gaussianity from Large Scale Structure that can be usefully applied to any model of primordial non-Gaussianity. When applied to the standard local ansatz, the constraint corresponds to fNL=-3+78-91. When applied to a model with a local-shape bispectrum but higher cumulants that are relatively more important, the amplitude of the local-shape bispectrum is constrained to be fNL=-14+22-21. This sensitivity implies that cluster counts could be used to distinguish qualitatively different models for the primordial fluctuations that have identical bispectra. (Abstract Abridged.)

Combining galaxy cluster data from the ROSAT All-Sky Survey and the Chandra X-ray Observatory, cosmic microwave background data from the Wilkinson Microwave Anisotropy Probe, and galaxy clustering data from the WiggleZ Dark Energy Survey, the 6-degree Field Galaxy Survey and the Sloan Digital Sky Survey III, we test for consistency the cosmic growth of structure predicted by General Relativity (GR) and the cosmic expansion history predicted by the cosmological constant plus cold dark matter paradigm (LCDM). The combination of these three independent, well studied measurements of the evolution of the mean energy density and its fluctuations is able to break strong degeneracies between model parameters. We model the key properties of cosmic growth with the normalization of the matter power spectrum, sigma_8, and the cosmic growth index, gamma, and those of cosmic expansion with the mean matter density, Omega_m, the Hubble constant, H_0, and a kinematical parameter equivalent to that for the dark energy equation of state, w. For a spatially flat geometry, w=-1, and allowing for systematic uncertainties, we obtain sigma_8=0.785+-0.019 and gamma=0.570+0.064-0.063 (at the 68.3 per cent confidence level). Allowing both w and gamma to vary we find w=-0.950+0.069-0.070 and gamma=0.533+-0.080. To further tighten the constraints on the expansion parameters, we also include supernova, Cepheid variable and baryon acoustic oscillation data. For w=-1, we have gamma=0.616+-0.061. For our most general model with a free w, we measure Omega_m=0.278+0.012-0.011, H_0=70.0+-1.3 km s^-1 Mpc^-1 and w=-0.987+0.054-0.053 for the expansion parameters, and sigma_8=0.789+-0.019 and gamma=0.604+-0.078 for the growth parameters. These results are in excellent agreement with GR+LCDM (gamma~0.55; w=-1) and represent the tightest and most robust simultaneous constraint on cosmic growth and expansion to date.

We describe in detail our characterization of the compact radio source population in 140 GHz Bolocam observations of a set of 45 massive galaxy clusters. We use a combination of 1.4 and 30 GHz data to select a total of 28 probable cluster-member radio galaxies and also to predict their 140 GHz flux densities. All of these galaxies are steep-spectrum radio sources and they are found preferentially in the cool-core clusters within our sample. In particular, 11 of the 12 brightest cluster member radio sources are associated with cool-core systems. Although none of the individual galaxies are robustly detected in the Bolocam data, the ensemble-average flux density at 140 GHz is consistent with, but slightly lower than, the extrapolation from lower frequencies assuming a constant spectral index. In addition, our data indicate an intrinsic scatter of 30 percent around the power-law extrapolated flux densities at 140 GHz, although our data do not tightly constrain this scatter. For our cluster sample, which is composed of high-mass and moderate-redshift systems, we find that the maximum fractional change in the Sunyaev-Zel'dovich signal integrated over any single cluster due to the presence of these radio sources is 20 percent, and only 1/4 of the clusters show a fractional change of more than 1 percent. The amount of contamination is strongly dependent on cluster morphology, and nearly all of the clusters with more than 1 percent contamination are cool-core systems. This result indicates that radio contamination is not significant compared to current noise levels in 140 GHz images of massive clusters and is in good agreement with the level of radio contamination found in previous results based on lower frequency data or simulations.

We present an analysis of the X-ray bright point source population in 43 massive clusters of galaxies observed with the Chandra X-ray Observatory. We have constructed a catalog of 4210 rigorously selected X-ray point sources in these fields, which span a survey area of 4.2 square degrees. This catalog reveals a clear excess of sources when compared to deep blank-field surveys, which amounts to roughly 1 additional source per cluster, likely Active Galactic Nuclei (AGN) associated with the clusters. The excess sources are concentrated within the virial radii of the clusters, with the largest excess observed near the cluster centers. The average radial profile of the excess X-ray sources of the cluster are well described by a power law (N(r) ~ r^\beta) with an index of \beta ~ -0.5. An initial analysis using literature results on the mean profile of member galaxies in massive X-ray selected clusters indicates that the fraction of galaxies hosting X-ray AGN rises with increasing clustercentric radius, being approximately 5 to 10 times higher near the virial radius than in the central regions. This trend is qualitatively similar to that observed for star formation in cluster member galaxies.

We report the first investigation of cool-core properties of galaxy clusters selected via their Sunyaev--Zel'dovich (SZ) effect. We use 13 galaxy clusters uniformly selected from 178 deg^2 observed with the South Pole Telescope (SPT) and followed up by the Chandra X-ray Observatory. They form an approximately mass-limited sample (> 3 x 10^14 M_sun h^-1_70) spanning redshifts 0.3 < z < 1.1. Using previously published X-ray-selected cluster samples, we compare two proxies of cool-core strength: surface brightness concentration (cSB) and cuspiness ({\alpha}). We find that cSB is better constrained. We measure cSB for the SPT sample and find several new z > 0.5 cool-core clusters, including two strong cool cores. This rules out the hypothesis that there are no z > 0.5 clusters that qualify as strong cool cores at the 5.4{\sigma} level. The fraction of strong cool-core clusters in the SPT sample in this redshift regime is between 7% and 56% (95% confidence). Although the SPT selection function is significantly different from the X-ray samples, the high-z cSB distribution for the SPT sample is statistically consistent with that of X-ray-selected samples at both low and high redshifts. The cool-core strength is inversely correlated with the offset between the brightest cluster galaxy and the X-ray centroid, providing evidence that the dynamical state affects the cool-core strength of the cluster. Larger SZ-selected samples will be crucial in understanding the evolution of cluster cool cores over cosmic time.

The galaxy cluster SPT-CL J0205-5829 currently has the highest spectroscopically-confirmed redshift, z=1.322, in the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. XMM-Newton observations measure a core-excluded temperature of Tx=8.7keV producing a mass estimate that is consistent with the Sunyaev-Zel'dovich derived mass. The combined SZ and X-ray mass estimate of M500=(4.9+/-0.8)e14 h_{70}^{-1} Msun makes it the most massive known SZ-selected galaxy cluster at z>1.2 and the second most massive at z>1. Using optical and infrared observations, we find that the brightest galaxies in SPT-CL J0205-5829 are already well evolved by the time the universe was <5 Gyr old, with stellar population ages >3 Gyr, and low rates of star formation (<0.5Msun/yr). We find that, despite the high redshift and mass, the existence of SPT-CL J0205-5829 is not surprising given a flat LambdaCDM cosmology with Gaussian initial perturbations. The a priori chance of finding a cluster of similar rarity (or rarer) in a survey the size of the 2500 deg^2 SPT-SZ survey is 69%.

We use weak gravitational lensing to measure the masses of five galaxy clusters selected from the South Pole Telescope (SPT) survey, with the primary goal of comparing these with the SPT Sunyaev--Zel'dovich (SZ) and X-ray based mass estimates. The clusters span redshifts 0.28 < z < 0.43 and have masses M_500 > 2 x 10^14 h^-1 M_sun, and three of the five clusters were discovered by the SPT survey. We observed the clusters in the g'r'i' passbands with the Megacam imager on the Magellan Clay 6.5m telescope. We measure a mean ratio of weak lensing (WL) aperture masses to inferred aperture masses from the SZ data, both within an aperture of R_500,SZ derived from the SZ mass, of 1.04 +/- 0.18. We measure a mean ratio of spherical WL masses evaluated at R_500,SZ to spherical SZ masses of 1.07 +/- 0.18, and a mean ratio of spherical WL masses evaluated at R_500,WL to spherical SZ masses of 1.10 +/- 0.24. We explore potential sources of systematic error in the mass comparisons and conclude that all are subdominant to the statistical uncertainty, with dominant terms being cluster concentration uncertainty and N-body simulation calibration bias. Expanding the sample of SPT clusters with WL observations has the potential to significantly improve the SPT cluster mass calibration and the resulting cosmological constraints from the SPT cluster survey. These are the first WL detections using Megacam on the Magellan Clay telescope.

By combining large-scale mosaics of ROSAT PSPC, XMM-Newton, and Suzaku X-ray observations, we present evidence for large-scale motions in the intracluster medium of the nearby, X-ray bright Perseus Cluster. These motions are suggested by several alternating and interleaved X-ray bright, low-temperature, low-entropy arcs located along the east-west axis, at radii ranging from ~10 kpc to over a Mpc. Thermodynamic features qualitatively similar to these have previously been observed in the centers of cool core clusters, and were successfully modeled as a consequence of the gas sloshing/swirling motions induced by minor mergers. Our observations indicate that such sloshing/swirling can extend out to larger radii than previously thought, on scales approaching the virial radius.

In the cores of some galaxy clusters the hot intracluster plasma is dense enough that it should cool radiatively in the cluster's lifetime, leading to continuous "cooling flows" of gas sinking towards the cluster center, yet no such cooling flow has been observed. The low observed star formation rates and cool gas masses for these "cool core" clusters suggest that much of the cooling must be offset by astrophysical feedback to prevent the formation of a runaway cooling flow. Here we report X-ray, optical, and infrared observations of the galaxy cluster SPT-CLJ2344-4243 at z = 0.596. These observations reveal an exceptionally luminous (L_2-10 keV = 8.2 x 10^45 erg/s) galaxy cluster which hosts an extremely strong cooling flow (dM/dt = 3820 +/- 530 Msun/yr). Further, the central galaxy in this cluster appears to be experiencing a massive starburst (740 +/- 160 Msun/yr), which suggests that the feedback source responsible for preventing runaway cooling in nearby cool core clusters may not yet be fully established in SPT-CLJ2344-4243. This large star formation rate implies that a significant fraction of the stars in the central galaxy of this cluster may form via accretion of the intracluster medium, rather than the current picture of central galaxies assembling entirely via mergers.

We present a catalog of 224 galaxy cluster candidates, selected through their Sunyaev-Zel'dovich (SZ) effect signature in the first 720 deg2 of the South Pole Telescope (SPT) survey. This area was mapped with the SPT in the 2008 and 2009 austral winters to a depth of 18 uK-arcmin at 150 GHz; 550 deg2 of it was also mapped to 44 uK-arcmin at 95 GHz. Based on optical imaging of all candidates and near-infrared imaging of the majority of candidates, we have found optical and/or infrared counterparts for 158 clusters. Of these, 135 were first identified as clusters in SPT data, including 117 new discoveries reported in this work. This catalog triples the number of confirmed galaxy clusters discovered through the SZ effect. We report photometrically derived (and in some cases spectroscopic) redshifts for confirmed clusters and redshift lower limits for the remaining candidates. The catalog extends to high redshift with a median redshift of z = 0.55 and maximum redshift of z = 1.37. Based on simulations, we expect the catalog to be nearly 100% complete above M500 ~ 5e14 Msun h_{70}^-1 at z > 0.6. There are 121 candidates detected at signal-to-noise greater than five, at which the catalog purity is measured to be 95%. From this high-purity subsample, we exclude the z < 0.3 clusters and use the remaining 100 candidates to improve cosmological constraints following the method presented by Benson et al., 2011. Adding the cluster data to CMB+BAO+H0 data leads to a preference for non-zero neutrino masses while only slightly reducing the upper limit on the sum of neutrino masses to sum mnu < 0.38 eV (95% CL). For a spatially flat wCDM cosmological model, the addition of this catalog to the CMB+BAO+H0+SNe results yields sigma8=0.807+-0.027 and w = -1.010+-0.058, improving the constraints on these parameters by a factor of 1.4 and 1.3, respectively. [abbrev]

We demonstrate the Sunyaev-Zel'dovich (SZ) effect imaging capabilities of the Combined Array for Research in Millimeter-wave Astronomy (CARMA) by presenting an SZ map of the galaxy cluster RXJ1347.5-1145. By combining data from multiple CARMA bands and configurations, we are able to capture the structure of this cluster over a wide range of angular scales, from its bulk properties to its core morphology. We find that roughly 9% of this cluster's thermal energy is associated with sub-arcminute-scale structure imparted by a merger, illustrating the value of high-resolution SZ measurements for pursuing cluster astrophysics and for understanding the scatter in SZ scaling relations. We also find that the cluster's SZ signal is lower in amplitude than suggested by a spherically-symmetric model derived from X-ray data, consistent with compression along the line of sight relative to the plane of the sky. Finally, we discuss the impact of upgrades currently in progress that will further enhance CARMA's power as an SZ imaging instrument.

We perform a joint analysis of X-ray and Sunyaev Zel'dovich (SZ) effect data using an analytic model that describes the gas properties of galaxy clusters. The joint analysis allows the measurement of the cluster gas mass fraction profile and Hubble constant independent of cosmological parameters. Weak cosmological priors are used to calculate the overdensity radius within which the gas mass fractions are reported. Such an analysis can provide direct constraints on the evolution of the cluster gas mass fraction with redshift. We validate the model and the joint analysis on high signal-to-noise data from the Chandra X-ray Observatory and the Sunyaev-Zel'dovich Array for two clusters, Abell 2631 and Abell 2204.

We use measurements from the South Pole Telescope (SPT) Sunyaev Zel'dovich (SZ) cluster survey in combination with X-ray measurements to constrain cosmological parameters. We present a statistical method that fits for the scaling relations of the SZ and X-ray cluster observables with mass while jointly fitting for cosmology. The method is generalizable to multiple cluster observables, and self-consistently accounts for the effects of the cluster selection and uncertainties in cluster mass calibration on the derived cosmological constraints. We apply this method to a data set consisting of an SZ-selected catalog of 18 galaxy clusters at z > 0.3 from the first 178 deg2 of the 2500 deg2 SPT-SZ survey, with 14 clusters having X-ray observations from either Chandra or XMM. Assuming a spatially flat LCDM cosmological model, we find the SPT cluster sample constrain sigma_8 (Omega_m/0.25)^0.30 = 0.785 +- 0.037. In combination with measurements of the CMB power spectrum from the SPT and the seven-year WMAP data, the SPT cluster sample constrain sigma_8 = 0.795 +- 0.016 and Omega_m = 0.255 +- 0.016, a factor of 1.5 improvement on each parameter over the CMB data alone. We consider several extensions beyond the LCDM model by including the following as free parameters: the dark energy equation of state (w), the sum of the neutrino masses (sum mnu), the effective number of relativistic species (Neff), and a primordial non-Gaussianity (fNL). We find that adding the SPT cluster data significantly improves the constraints on w and sum mnu beyond those found when using measurements of the CMB, supernovae, baryon acoustic oscillations, and the Hubble constant. Considering each extension independently, we best constrain w=-0.973 +- 0.063 and the sum of neutrino masses sum mnu < 0.28 eV at 95% confidence, a factor of 1.25 and 1.4 improvement, respectively, over the constraints without clusters. [abbrev.]

Studies of the diffuse X-ray emitting gas in galaxy clusters have provided powerful constraints on cosmological parameters and insights into plasma astrophysics. However, measurements of the faint cluster outskirts have become possible only over the last few years. Here, we present results from Suzaku observations of the Perseus Cluster, which provide our best measurements of the thermodynamic properties of the ICM at large radii to date. In particular, we focus on the details of the data analysis procedure and discuss the evidence for a clumpy distribution of the gas in the outskirts, which is important for understanding the physics of the ongoing growth of clusters from the surrounding cosmic web.

Deriving the total masses of galaxy clusters from observations of the intracluster medium (ICM) generally requires some prior information, in addition to the assumptions of hydrostatic equilibrium and spherical symmetry. Often, this information takes the form of particular parametrized functions used to describe the cluster gas density and temperature profiles. In this paper, we investigate the implicit priors on hydrostatic masses that result from this fully parametric approach, and the implications of such priors for scaling relations formed from those masses. We show that the application of such fully parametric models of the ICM naturally imposes a prior on the slopes of the derived scaling relations, favoring the self-similar model, and argue that this prior may be influential in practice. In contrast, this bias does not exist for techniques which adopt an explicit prior on the form of the mass profile but describe the ICM non-parametrically. Constraints on the slope of the cluster mass--temperature relation in the literature show a separation based the approach employed, with the results from fully parametric ICM modeling clustering nearer the self-similar value. Given that a primary goal of scaling relation analyses is to test the self-similar model, the application of methods subject to strong, implicit priors should be avoided. Alternative methods and best practices are discussed.

Studies of galaxy clusters have proved crucial in helping to establish the standard model of cosmology, with a universe dominated by dark matter and dark energy. A theoretical basis that describes clusters as massive, multi-component, quasi-equilibrium systems is growing in its capability to interpret multi-wavelength observations of expanding scope and sensitivity. We review current cosmological results, including contributions to fundamental physics, obtained from observations of galaxy clusters. These results are consistent with and complementary to those from other methods. We highlight several areas of opportunity for the next few years, and emphasize the need for accurate modeling of survey selection and sources of systematic error. Capitalizing on these opportunities will require a multi-wavelength approach and the application of rigorous statistical frameworks, utilizing the combined strengths of observers, simulators and theorists.

Studies of the diffuse X-ray emitting gas in galaxy clusters have provided powerful constraints on cosmological parameters and insights into plasma astrophysics. However, measurements of the faint cluster outskirts have become possible only recently. Using data from the Suzaku X-ray telescope, we determined an accurate, spatially resolved census of the gas, metals, and dark matter out to the edge of the Perseus Cluster. Contrary to previous results, our measurements of the cluster baryon fraction are consistent with the expected universal value at half of the virial radius. The apparent baryon fraction exceeds the cosmic mean at larger radii, suggesting a clumpy distribution of the gas, which is important for understanding the ongoing growth of clusters from the surrounding cosmic web.

This is the fourth 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. Here we examine the constraints on neutrino properties that are enabled by the precise and robust constraint on the amplitude of the matter power spectrum at low redshift that is available from our data. In combination with cluster gas-mass fraction, cosmic microwave background, supernova and baryon acoustic oscillation data, and incorporating conservative allowances for systematic uncertainties, we limit the species-summed neutrino mass, M_nu, to <0.33 eV at 95.4 per cent confidence in a spatially flat, cosmological constant (LambdaCDM) model. In a flat LambdaCDM model where the effective number of neutrino species, N_eff, is allowed to vary, we find N_eff = 3.4 -0.5 +0.6 (68.3 per cent confidence, incorporating a direct constraint on the Hubble parameter from Cepheid and supernova data). We also obtain results with additional degrees of freedom in the cosmological model, in the form of global spatial curvature (Omega_k) and a primordial spectrum of tensor perturbations (r and n_t). The results are not immune to these generalizations; however, in the most general case we consider, in which M_nu, N_eff, curvature and tensors are all free, we still obtain M_nu < 0.70 eV and N_eff = 3.7 +- 0.7 (at respectively the same confidence levels as above). These results agree well with recent work using independent data, and highlight the importance of measuring cosmic structure and expansion at low as well as high redshifts. Although our cluster data extend to redshift z=0.5, the effect of neutrino mass on the growth of structure at late times is not yet detected at a significant level.

(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 from 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 consists 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 a statistically complete sample of very X-ray luminous galaxy clusters detected in the MAssive Cluster Survey (MACS). This second MACS release comprises all 34 MACS clusters with nominal X-ray fluxes in excess of 2x10^(-12) erg/s/cm^2 (0.1-2.4 keV) in the ROSAT Bright Source Catalogue; two thirds of them are new discoveries. Extending over the redshift range from 0.3 to 0.5, this subset complements the complete sample of the 12 most distant MACS clusters (z>0.5) published in 2007 and further exemplifies the efficacy of X-ray selection for the compilation of samples of intrinsically massive galaxy clusters. Extensive follow-up observations with Chandra/ACIS led to three additional MACS cluster candidates being eliminated as (predominantly) X-ray point sources. For another four clusters --- which, however, remain in our sample of 34 --- the point-source contamination was found to be about 50%. The median X-ray luminosity of 1.3x10^45 erg/s (0.1-2.4 keV, Chandra, within r_500) of the clusters in this subsample demonstrates the power of the MACS survey strategy to find the most extreme and rarest clusters out to significant redshift. A comparison of the optical and X-ray data for all clusters in this release finds a wide range of morphologies with no obvious bias in favour of either relaxed or merging systems.

This is the third 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 consists of 238 clusters drawn from the ROSAT All-Sky Survey, and incorporates extensive follow-up observations using the Chandra X-ray Observatory. Here we present improved constraints on departures from General Relativity (GR) on cosmological scales, using the growth index, gamma, to parameterize the linear growth rate of cosmic structure. Using the method of Mantz et al. (2009a), we simultaneously and self-consistently model the growth of X-ray luminous clusters and their observable-mass scaling relations, accounting for survey biases, parameter degeneracies and systematic uncertainties. We combine the cluster growth data with gas mass fraction, SNIa, BAO and CMB data. This combination leads to a tight correlation between gamma and sigma_8. Consistency with GR requires gamma~0.55. Under the assumption of self-similar evolution and constant scatter in the scaling relations, and for a flat LCDM model, we measure gamma(sigma_8/0.8)^6.8=0.55+0.13-0.10, with 0.79<sigma_8<0.89. Relaxing the assumptions on the scaling relations by introducing two additional parameters to model possible evolution in the normalization and scatter of the luminosity-mass relation, we obtain consistent constraints on gamma that are only ~20% weaker than those above. Allowing the dark energy equation of state, w, to take any constant value, we simultaneously constrain the growth and expansion histories, and find no evidence for departures from either GR or LCDM. Our results represent the most robust consistency test of GR on cosmological scales to date. (Abridged)

We use measurements of the growth of cosmic structure, as inferred from the observed evolution of the X-ray luminosity function (XLF) of galaxy clusters, to constrain departures from General Relativity (GR) on cosmological scales. We employ the popular growth rate parameterization, Omega_m(z)^gamma, for which GR predicts a growth index gamma~0.55. We use observations of the cosmic microwave background (CMB), type Ia supernovae (SNIa), and X-ray cluster gas-mass fractions (fgas), to simultaneously constrain the expansion history and energy content of the Universe, as described by the background model parameters: Omega_m, w, and Omega_k, i.e., the mean matter density, the dark energy equation of state parameter, and the mean curvature, respectively. Using conservative allowances for systematic uncertainties, in particular for the evolution of the mass-luminosity scaling relation in the XLF analysis, we find gamma=0.51+0.16-0.15 and Omega_m=0.27+-0.02 (68.3 per cent confidence limits), for a flat cosmological constant (LCDM) background model. Allowing w to be a free parameter, we find gamma=0.44+0.17-0.15. Relaxing the flatness prior in the LCDM model, we obtain gamma=0.51+0.19-0.16. When in addition to the XLF data we use the CMB data to constrain gamma through the ISW effect, we obtain a combined constraint of gamma=0.45+0.14-0.12 for the flat LCDM model. Our analysis provides the tightest constraints to date on the growth index. We find no evidence for departures from General Relativity on cosmological scales.

We examine the ability of a future X-ray observatory to constrain dark energy via measurements of the cluster X-ray gas mass fraction, fgas. We find that fgas measurements for a sample of ~500 hot, X-ray bright, dynamically relaxed clusters, to a precision of ~5 per cent, can be used to constrain dark energy with a Dark Energy Task Force (DETF) figure of merit of 15-40, with the possibility of boosting these values by 40 per cent or more by optimizing the redshift distribution of target clusters. Such constraints are comparable to those predicted by the DETF for other leading, planned dark energy experiments. A future fgas experiment will be preceded by a large X-ray or SZ survey that will find hot, X-ray luminous clusters out to high redshifts. Short `snapshot' observations with the new X-ray observatory should then be able to identify a sample of ~500 suitably relaxed systems. The redshift, temperature and X-ray luminosity range of interest has already been partially probed by existing X-ray cluster surveys which allow reasonable estimates of the fraction of clusters that will be suitably relaxed for fgas work. Our analysis uses a Markov Chain Monte Carlo method which fully captures the relevant degeneracies between parameters and facilitates the incorporation of priors and systematic uncertainties in the analysis. We explore the effects of such uncertainties for scenarios ranging from optimistic to pessimistic. We conclude that the fgas experiment will provide tight constraints on the mean matter and dark energy densities, with a peak sensitivity for dark energy work at redshifts midway between those of supernovae and baryon acoustic oscillation/weak lensing/cluster number counts experiments. In combination, these experiments should enable a precise measurement of the evolution of dark energy. (Abridged)

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 largest 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)

The galaxy cluster RX J1347-1145 is one of the most X-ray luminous and most massive clusters known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and cosmology. In this paper we present new high-resolution HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing analysis. Combining this information with shape measurements of weak lensing sources in the central regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map provides the best available quantification of the total mass of the central part of the cluster to date. We compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude that both mass estimates agree extremely well in the observed region, namely within 400 / h_70 kpc of the cluster center. In addition we study the major baryonic components (gas and stars) and hence derive the dark matter distribution in the center of the cluster. We find that the dark matter and baryons are both centered on the BCG within the uncertainties (alignment is better than <10 kpc). We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with both NFW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint the concentration parameter, and hence the inner slope of the dark matter profile.