Search

The X-ray morphologies of clusters of galaxies display significant variations, reflecting their dynamical histories and the nonlinear dependence of X-ray emissivity on the density of the intracluster gas. Qualitative and quantitative assessments of X-ray morphology have long been considered a proxy for determining whether clusters are dynamically active or "relaxed." Conversely, the use of circularly or elliptically symmetric models for cluster emission can be complicated by the variety of complex features realized in nature, spanning scales from Mpc down to the resolution limit of current X-ray observatories. In this work, we use mock X-ray images from simulated clusters from THE THREE HUNDRED project to define a basis set of cluster image features. We take advantage of clusters' approximate self similarity to minimize the differences between images before encoding the remaining diversity through a distribution of high order polynomial coefficients. Principal component analysis then provides an orthogonal basis for this distribution, corresponding to natural perturbations from an average model. This representation allows novel, realistically complex X-ray cluster images to be easily generated, and we provide code to do so. The approach provides a simple way to generate training data for cluster image analysis algorithms, and could be straightforwardly adapted to generate clusters displaying specific types of features, or selected by physical characteristics available in the original simulations., Comment: 12 pages, 9 figures. To appear in ApJ

The most dynamically relaxed clusters of galaxies play a special role in cosmological studies as well as astrophysical studies of the intracluster medium (ICM) and active galactic nucleus feedback. While high spatial resolution imaging of the morphology of the ICM has long been the gold standard for establishing a cluster's dynamical state, such data are not available from current or planned surveys, and thus require separate, pointed follow-up observations. With optical and/or near-IR photometric imaging, and red-sequence cluster finding results from those data, expected to be ubiquitously available for clusters discovered in upcoming optical and mm-wavelength surveys, it is worth asking how effectively photometric data alone can identify relaxed cluster candidates, before investing in, e.g., high-resolution X-ray observations. Here we assess the ability of several simple photometric measurements, based on the redMaPPer cluster finder run on Sloan Digital Sky Survey data, to reproduce X-ray classifications of dynamical state for an X-ray selected sample of massive clusters. We find that two simple metrics contrasting the Bright Central Galaxy (BCG) to other cluster members can identify a complete sample of relaxed clusters with a purity of ~40 per cent in our data set. Including minimal ICM information in the form of a center position increases the purity to ~60 per cent. However, all three metrics depend critically on correctly identifying the BCG, which is presently a challenge for optical red-sequence cluster finders., Comment: 10 pages, 6 figures. Accepted version; to appear in ApJ

We present the discovery of the most distant, dynamically relaxed cool core cluster, SPT-CL J2215-3537 (SPT2215) and its central brightest cluster galaxy (BCG) at z=1.16. Using new X-ray observations, we demonstrate that SPT2215 harbors a strong cool core, with a central cooling time of 200 Myr (at 10 kpc) and a maximal intracluster medium cooling rate of 1900+/-400 Msun/yr. This prodigious cooling may be responsible for fueling extended, star-forming filaments observed in Hubble Space Telescope imaging. Based on new spectrophotometric data, we detect bright [O II] emission in the BCG, implying an unobscured star formation rate (SFR) of 320^{+230}_{-140} Msun/yr. The detection of a weak radio source (2.0+/-0.8 mJy at 0.8 GHz) suggests ongoing feedback from an active galactic nucleus (AGN), though the implied jet power is less than half the cooling luminosity of the hot gas, consistent with cooling overpowering heating. The extreme cooling and SFR of SPT2215 is rare among known cool core clusters, and it is even more remarkable that we observe these at such high redshift, when most clusters are still dynamically disturbed. The high mass of this cluster, coupled with the fact that it is dynamically relaxed with a highly-isolated BCG, suggests that it is an exceptionally rare system that must have formed very rapidly in the early Universe. Combined with the high SFR, SPT2215 may be a high-z analog of the Phoenix cluster, potentially providing insight into the limits of AGN feedback and star formation in the most massive galaxies., Comment: 14 pages, 5 figures, 1 table. Accepted in ApJ

The relationship linking a galaxy cluster's total mass with the concentration of its mass profile and its redshift is a fundamental prediction of the Cold Dark Matter (CDM) paradigm of cosmic structure formation. However, confronting those predictions with observations is complicated by the fact that simulated clusters are not representative of observed samples where detailed mass profile constraints are possible. In this work, we calculate the Symmetry-Peakiness-Alignment (SPA) morphology metrics for maps of X-ray emissivity from THE THREE HUNDRED project hydrodynamical simulations of galaxy clusters at four redshifts, and thereby select a sample of morphologically relaxed, simulated clusters, using observational criteria. These clusters have on average earlier formation times than the full sample, confirming that they are both morphologically and dynamically more relaxed than typical. We constrain the concentration-mass-redshift relation of both the relaxed and complete sample of simulated clusters, assuming power-law dependences on mass ($\kappa_m$) and $1+z$ ($\kappa_\zeta$), finding $\kappa_m = -0.12 \pm 0.07$ and $\kappa_\zeta = -0.27 \pm 0.19$ for the relaxed subsample. From an equivalently selected sample of massive, relaxed clusters observed with ${\it Chandra}$, we find $\kappa_m = -0.12 \pm 0.08$ and $\kappa_\zeta = -0.48 \pm 0.19$, in good agreement with the simulation predictions. The simulated and observed samples also agree well on the average concentration at a pivot mass and redshift providing further validation of the $\Lambda$CDM paradigm in the properties of the largest gravitationally collapsed structures observed. This also represents the first clear detection of decreasing concentration with redshift, a longstanding prediction of simulations, in data., Comment: 10 pages, 10 figures. Accepted by MNRAS

We present a Chandra study of surface brightness fluctuations in the diffuse intracluster medium of the Perseus Cluster. Our study utilizes deep, archival imaging of the cluster core as well as a new mosaic of 29 short 5 ks observations extending in 8 different directions out to radii of r_500 ~ 2.2r_2500. Under the assumption that the distribution of densities at a given radius is log-normally distributed, two important quantities can be derived from the width of the log-normal density distribution on a given spatial scale: the density bias, which is equal to the square root of the clumping factor C; and the one-component turbulent velocity, v_(k, 1D). We forward-model all contributions to the measured surface brightness, including astrophysical and particle background components, and account for the Poisson nature of the measured signal. Measuring the distribution of surface brightness fluctuations in 1 arcmin^2 regions, spanning the radial range 0.3-2.2 r_2500 (7.8-57.3 arcmin), we find a small to moderate average density bias of around 3% at radii below 1.6r_2500. We also infer an average turbulent velocity at these radii of v_1D <400 km s^-1. Direct confirmation of our results on turbulent velocities inferred from surface brightness fluctuations should be possible using the X-ray calorimeter spectrometers to be flown aboard the XRISM and Athena. observatories., Comment: 17 pages, 11 figures. to be published in MNRAS

We have combined X-ray observations from Chandra with Sunyaev-Zel'dovich (SZ) effect data from Planck and Bolocam to measure intra-cluster medium pressure profiles from 0.03R$_{500}$ $\le$ R $\le$ 5R$_{500}$ for a sample of 21 low-$z$ galaxy clusters with a median redshift $\langle z \rangle = 0.08$ and a median mass $\langle \textrm{M}_{500} \rangle = 6.1 \times 10^{14}$ M$_{\odot}$ and a sample of 19 mid-$z$ galaxy clusters with $\langle z \rangle = 0.50$ and $\langle \textrm{M}_{500} \rangle = 10.6 \times 10^{14}$ M$_{\odot}$. The mean scaled pressure in the low-$z$ sample is lower at small radii and higher at large radii, a trend that is accurately reproduced in similarly selected samples from The300 simulations. This difference appears to be primarily due to dynamical state at small radii, evolution at intermediate radii, and a combination of evolution and mass dependence at large radii. Furthermore, the overall flattening of the mean scaled pressure profile in the low-$z$ sample compared to the mid-$z$ sample is consistent with expectations due to differences in mass accretion rate and the fractional impact of feedback mechanisms. In agreement with previous studies, the fractional scatter about the mean scaled pressure profile reaches a minimum of $\simeq 20$ per cent near 0.5R$_{500}$. This scatter is consistent between the low-$z$ and mid-$z$ samples at all radii, suggesting it is not strongly impacted by sample selection, and this general behavior is reproduced in The300 simulations. Finally, analytic functions that approximately describe the mass and redshift trends in mean pressure profile shape are provided., Comment: Accepted for publication in ApJ. Updated to closely match the published version

We present updated cosmological constraints from measurements of the gas mass fractions ($f_{gas}$) of massive, dynamically relaxed galaxy clusters. Our new data set has greater leverage on models of dark energy, thanks to the addition of the Perseus Cluster at low redshifts, two new clusters at redshifts $z>0.97$, and significantly longer observations of four clusters at $0.6

We present a measurement of the relativistic corrections to the thermal Sunyaev-Zel'dovich (SZ) effect spectrum, the rSZ effect, toward the massive galaxy cluster RX J1347.5-1145 by combining sub-mm images from Herschel-SPIRE with mm-wave Bolocam maps. Our analysis simultaneously models the SZ effect signal, the population of cosmic infrared background (CIB) galaxies, and galactic cirrus dust emission in a manner that fully accounts for their spatial and frequency-dependent correlations. Gravitational lensing of background galaxies by RX J1347.5-1145 is included in our methodology based on a mass model derived from HST observations. Utilizing a set of realistic mock observations, we employ a forward modelling approach that accounts for the non-Gaussian covariances between observed astrophysical components to determine the posterior distribution of SZ effect brightness values consistent with the observed data. We determine a maximum a posteriori (MAP) value of the average Comptonization parameter of the intra-cluster medium (ICM) within R$_{2500}$ to be $\langle y \rangle_{2500} = 1.56 \times 10^{-4}$, with corresponding 68~per cent credible interval $[1.42,1.63] \times 10^{-4}$, and a MAP ICM electron temperature of $\langle \textrm{T}_{\textrm{sz}} \rangle_{2500} = 22.4$~keV with 68~per cent credible interval spanning $[10.4,33.0]$~keV. This is in good agreement with the pressure-weighted temperature obtained from {\it Chandra} X-ray observations, $\langle \textrm{T}_{\textrm{x,pw}}\rangle_{2500} = 17.4 \pm 2.3$~keV. We aim to apply this methodology to comparable existing data for a sample of 39 galaxy clusters, with an estimated uncertainty on the ensemble mean $\langle \textrm{T}_{\textrm{sz}} \rangle_{2500}$ at the $\simeq 1$~keV level, sufficiently precise to probe ICM physics and to inform X-ray temperature calibration., Comment: 20 pages, 9 figures

We present the analysis of deep X-ray observations of 10 massive galaxy clusters at redshifts $1.05 < z < 1.71$, with the primary goal of measuring the metallicity of the intracluster medium (ICM) at intermediate radii, to better constrain models of the metal enrichment of the intergalactic medium. The targets were selected from X-ray and Sunyaev-Zel'dovich (SZ) effect surveys, and observed with both the \textit{XMM-Newton} and \textit{Chandra} satellites. For each cluster, a precise gas mass profile was extracted, from which the value of $r_{500}$ could be estimated. This allows us to define consistent radial ranges over which the metallicity measurements can be compared. In general, the data are of sufficient quality to extract meaningful metallicity measurements in two radial bins, $r<0.3r_{500}$ and $0.3

We use a sample of 14 massive, dynamically relaxed galaxy clusters to constrain the Hubble Constant, $H_0$, by combining X-ray and Sunyaev-Zel'dovich (SZ) effect signals measured with Chandra, Planck and Bolocam. This is the first such analysis to marginalize over an empirical, data-driven prior on the overall accuracy of X-ray temperature measurements, while our restriction to the most relaxed, massive clusters also minimizes astrophysical systematics. For a cosmological-constant model with $\Omega_m = 0.3$ and $\Omega_{\Lambda} = 0.7$, we find $H_0 = 67.3^{+21.3}_{-13.3}$ km/s/Mpc, limited by the temperature calibration uncertainty (compared to the statistically limited constraint of $H_0 = 72.3^{+7.6}_{-7.6}$ km/s/Mpc). The intrinsic scatter in the X-ray/SZ pressure ratio is found to be $13 \pm 4$ per cent ($10 \pm 3$ per cent when two clusters with significant galactic dust emission are removed from the sample), consistent with being primarily due to triaxiality and projection. We discuss the prospects for reducing the dominant systematic limitation to this analysis, with improved X-ray calibration and/or precise measurements of the relativistic SZ effect providing a plausible route to per cent level constraints on $H_0$., Comment: 16 pages, 9 figures, submitted to MNRAS

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., Comment: 12 pages, to be published in MNRAS

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., Comment: 9 pages, 3 figures, 3 tables

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., Comment: Science white paper submitted to the Astro2020 US Decadal Survey

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., Comment: Science white paper submitted for the Astro2020 decadal review, 5 pages + references

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., Comment: 11 pages, to appear in MNRAS

We present the Massive and Distant Clusters of WISE Survey (MaDCoWS), a search for galaxy clusters at 0.7-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)$., Comment: 26 pages, 17 figures, submitted to The Astrophysical Journal Supplement Series. Machine readable versions of tables 3-6 are included with the source files in the arXiv submission

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., Comment: 15 pages, 9 figures, submitted to ApJ Jun 2018, Accepted Dec 2018, Published Jan 2019. This is the version of the article before editing, as submitted by an author to ApJ. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at https://doi.org/10.3847/1538-4357/aaf888

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., Comment: 42 pages, 52 figures

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

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., Comment: 9 pages. Accepted version

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., Comment: 12 pages, A&A accepted

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., Comment: 14 pages, 6 figures, 3 tables. Accepted for publiction in MNRAS

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., Comment: 8 pages, to appear in MNRAS. Minor typos corrected in v2

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

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$., Comment: Updated to include one additional co-author. Also some minor changes to the text based on initial feedback

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., Comment: Submitted to ApJ Letters, 6 pages, 4 figures

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., Comment: Accepted version. 16 pages (+5 pages of appendix), 12 figures, 5 tables

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., Comment: 11 pages, 5 figures, 2 tables. Code is available on GitHub at https://github.com/abmantz/lrgs and from CRAN

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., Comment: MNRAS, in press. 43 pages in total, of which 17 are tables (please think twice before printing). 18 figures, 4 tables. Machine-readable tables will be available from the journal and at the url below; code will be posted at http://www.slac.stanford.edu/~amantz/work/morph14/

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., Comment: 11 pages, 1 Table, 3 figures, matches published version on PRD

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., Comment: 24 pages, 10 figures, 5 tables. v3: typo in table A1 corrected

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), Comment: 25 pages, 14 figures, 8 tables. Accepted by MNRAS. Code and data can be downloaded from http://www.slac.stanford.edu/~amantz/work/fgas14/ . v2: minor fix to table 1, updated bibliography

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., Comment: ApJ, in press. 9 pages, 4 figures, 4 tables

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., Comment: 53 pages, 21 figures. To appear in Annual Review of Astronomy & Astrophysics