Search

We present the splashback radius measurements around the SRG/eROSITA eFEDS X-ray selected galaxy clusters by cross-correlating them with HSC S19A photometric galaxies. The X-ray selection is expected to be less affected by systematics related to projection that affects optical cluster finder algorithms. We use a nearly volume-limited sample of 109 galaxy clusters selected in 0.5-2.0 keV band having luminosity $L_X > 10^{43.5}\,{\rm erg s^{-1} h^{-2}}$ within the redshift $z, Comment: 15 pages, 10 figures

We present our determination of the baryon budget for an X-ray-selected XXL sample of 136 galaxy groups and clusters spanning nearly two orders of magnitude in mass ($M_{500}\sim 10^{13}-10^{15}M_\odot$) and the redshift range $0< z < 1$. Our joint analysis is based on the combination of HSC-SSP weak-lensing mass measurements, XXL X-ray gas mass measurements, and HSC and SDSS multiband photometry. We carry out a Bayesian analysis of multivariate mass-scaling relations of gas mass, galaxy stellar mass, stellar mass of brightest cluster galaxies (BCGs), and soft-band X-ray luminosity, by taking into account the intrinsic covariance between cluster properties, selection effect, weak-lensing mass calibration, and observational error covariance matrix. The mass-dependent slope of the gas mass--total mass ($M_{500}$) relation is found to be $1.29_{-0.10}^{+0.16}$, which is steeper than the self-similar prediction of unity, whereas the slope of the stellar mass--total mass relation is shallower than unity, $0.85_{-0.09}^{+0.12}$. The BCG stellar mass weakly depends on cluster mass with a slope of $0.49_{-0.10}^{+0.11}$. The baryon, gas mass, and stellar mass fractions as a function of $M_{500}$ agree with the results from numerical simulations and previous observations. We successfully constrain the full intrinsic covariance of the baryonic contents. The BCG stellar mass shows the larger intrinsic scatter at a given halo total mass, followed in order by stellar mass and gas mass. We find a significant positive intrinsic correlation coefficient between total (and satellite) stellar mass and BCG stellar mass and no evidence for intrinsic correlation between gas mass and stellar mass. All the baryonic components show no redshift evolution., 34 pages, 16 figures, 8 tables; accepted for the publication in PASJ

We study the angular correlation function of star-forming galaxies and properties of their host dark matter halos at z>1 using the Hyper-Suprime Cam (HSC) SSP survey. We use [OII] emitters identified using two narrow-band (NB) filters, NB816 and NB921, in the Deep/UltraDeep layers, which respectively cover large angular areas of 16.3 deg^2 and 16.9 deg^2. Our sample contains 8302 and 9578 [OII] emitters at z=1.19 (NB816) and z=1.47 (NB921), respectively. We detect a strong clustering signal over a wide angular range, 0.001 < \theta < 1 [deg], with the bias $b=1.61^{+0.13}_{-0.11}$ (z=1.19) and $b=2.09^{+0.17}_{-0.15}$ (z=1.47). We also find a clear deviation of the correlation from a simple power-law form. To interpret the measured clustering signal, we adopt a halo occupation distribution (HOD) model that is constructed to explain the spatial distribution of galaxies selected by a star formation rate. The observed correlation function and number density are simultaneously explained by the best-fitting HOD model. From the constrained HOD model, the average mass of halos hosting the [OII] emitters is derived to be $\log{M_{eff}/(h^{-1}M_\odot)}=12.70^{+0.09}_{-0.07}$ and $12.61^{+0.09}_{-0.05}$ at z=1.19 and 1.47, respectively, which will become halos with the present-day mass, $M\sim 1.5 \times 10^{13}h^{-1}M_\odot$. The satellite fraction of the [OII] emitter sample is found to be $f_{sat}\sim 0.15$. All these values are consistent with the previous studies of similar samples, but we obtain tighter constraints even in a larger parameter space due to the larger sample size from the HSC. The results obtained for host halos of [OII] emitters in this paper enable the construction of mock galaxy catalogs and the systematic forecast study of cosmological constraints from upcoming emission line galaxy surveys such as the Subaru PFS survey., Comment: 21 pages, 13 figures, 5 tables, accepted to PASJ

We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843~MHz radio AGN population in galaxy clusters from two large catalogs created using the Dark Energy Survey (DES): $\sim$11,800 optically selected RM-Y3 and $\sim$1,000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGN are highly concentrated around cluster centers to $z\sim1$. We measure the halo occupation number for cluster radio AGN above a threshold luminosity, finding that the number of radio AGN per cluster increases with cluster halo mass as $N\propto M^{1.2\pm0.1}$ ($N\propto M^{0.68\pm0.34}$) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad band colors and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGN luminosity function (LF). The LF exhibits a redshift trend of the form $(1+z)^\gamma$ in density and luminosity, respectively, of $\gamma_{\rm D}=3.0\pm0.4$ and $\gamma_{\rm P}=0.21\pm0.15$ in the RM-Y3 sample, and $\gamma_{\rm D}=2.6\pm0.7$ and $\gamma_{\rm P}=0.31\pm0.15$ in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGN, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core ($, Comment: 20 pages, 15 figures Replaced with published version

We present a detailed weak-lensing and X-ray study of the Frontier Fields galaxy cluster Abell 370, one of the most massive known lenses on the sky, using wide-field BRz Subaru/Sprime-Cam and Chandra X-ray observations. By combining two-dimensional (2D) shear and azimuthally averaged magnification constraints derived from Subaru data, we perform a lensing mass reconstruction in a free-form manner, which allows us to determine both radial structure and 2D morphology of the cluster mass distribution. In a triaxial framework assuming a Navarro-Frenk-White density profile, we constrain the intrinsic structure and geometry of the cluster halo by forward modeling the reconstructed mass map. We obtain a halo mass $M_{200}=(1.54 \pm 0.29)\times 10^{15}h^{-1}M_\odot$, a halo concentration $c_{200}=5.27 \pm 1.28$, and a minor-major axis ratio $q_a=0.62 \pm 0.23$ with uninformative priors. Using a prior on the line-of-sight alignment of the halo major axis derived from binary merger simulations constrained by multi-probe observations, we find that the data favor a more prolate geometry with lower mass and lower concentration. From triaxial lens modeling with the line-of-sight prior, we find a spherically enclosed gas mass fraction of $f_\mathrm{gas}=(8.4 \pm 1.0)\%$ at $0.7h^{-1}$ Mpc $\sim 0.7r_{500}$. When compared to the hydrostatic mass estimate from Chandra observations, our triaxial weak-lensing analysis yields spherically enclosed mass ratios of $1-b \equiv M_\mathrm{HE}/M_\mathrm{WL} = 0.56 \pm 0.09$ and $0.51 \pm 0.09$ at $0.7h^{-1}$ Mpc with and without using the line-of-sight prior, respectively. Since the cluster is in a highly disturbed dynamical state, this represents the likely maximum level of hydrostatic bias in galaxy clusters., Comment: Version accepted by ApJ; 26 pages (main body) + 8 pages (appendices and references), 20 figures, 9 tables

We present the first cosmological study of a sample of $eROSITA$ clusters, which were identified in the $eROSITA$ Final Equatorial Depth Survey (eFEDS). In a joint selection on X-ray and optical observables, the sample contains $455$ clusters within a redshift range of $0.1, Comment: Accepted for publication in MNRAS. Figures 18 and 19 contain the main results. Chains and cluster masses are at https://github.com/inonchiu/eFEDSproducts

We investigate the scaling relations between X-ray observables of the clusters detected in the eFEDS field using Spectrum-Roentgen-Gamma/eROSITA observations taking into account the selection effects and the distributions of observables with cosmic time. We extract X-ray observables (Lx, Lbol, T, Mgas, Yx) within R500 for the sample of 542 clusters in the eFEDS field. By applying detection and extent likelihoods, we construct a subsample of 265 clusters with a contamination level of 3sigma confidence level however, our results are in good agreement with the simulations including non-gravitational physics and the recent results that take into account selection effects. Strong deviations we find from the self-similar scenario indicate that the non-gravitational effects play an important role in shaping the observed physical state of clusters. This work extends the scaling relations to low mass, low luminosity galaxy cluster and group regime using eFEDS observations, demonstrating eROSITA's ability to measure ICM emission out to R500 with survey-depth exposures and constrain the scaling relations in a wide mass-luminosity-redshift range., Comment: 35 pages, 7 figures, 5 tables, Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) observed by Spectrum Roentgen Gamma/$eROSITA$ over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $\eta$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($\eta$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature $T_{\mathrm{X}}$, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged), Comment: Accepted for publication in A&A. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Using the extended halo-based group finder developed by Yang et al. (2021), which is able to deal with galaxies via spectroscopic and photometric redshifts simultaneously, we construct galaxy group and candidate protocluster catalogs in a wide redshift range ($0 < z < 6$) from the joint CFHT Large Area $U$-band Deep Survey (CLAUDS) and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) deep data set. Based on a selection of 5,607,052 galaxies with $i$-band magnitude $m_{i} < 26$ and a sky coverage of $34.41\ {\rm deg}^2$, we identify a total of 2,232,134 groups, within which 402,947 groups have at least three member galaxies. We have visually checked and discussed the general properties of those richest groups at redshift $z>2.0$. By checking the galaxy number distributions within a $5-7\ h^{-1}\mathrm{Mpc}$ projected separation and a redshift difference $Δz \le 0.1$ around those richest groups at redshift $z>2$, we identified a list of 761, 343 and 43 protocluster candidates in the redshift bins $2\leq z, 25 pages, 16 figures, 2 tables, accepted for publication in ApJ

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions: a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The main SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that the integrated Compton $y$ parameter and the temperature for the major merger are boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations., Comment: 36 pages, 25 figures, 5 tables; accepted for the publication in MNRAS

We present weak lensing (WL) mass constraints for a sample of massive galaxy clusters detected by the South Pole Telescope (SPT) via the Sunyaev-Zeldovich effect (SZE). We use $griz$ imaging data obtained from the Science Verification (SV) phase of the Dark Energy Survey (DES) to fit the WL shear signal of 33 clusters in the redshift range $0.25 \le z \le 0.8$ with NFW profiles and to constrain a four-parameter SPT mass-observable relation. To account for biases in WL masses, we introduce a WL mass to true mass scaling relation described by a mean bias and an intrinsic, log-normal scatter. We allow for correlated scatter within the WL and SZE mass-observable relations and use simulations to constrain priors on nuisance parameters related to bias and scatter from WL. We constrain the normalization of the $\zeta-M_{500}$ relation, $A_\mathrm{SZ}=12.0_{-6.7}^{+2.6}$ when using a prior on the mass slope $B_\mathrm{SZ}$ from the latest SPT cluster cosmology analysis. Without this prior, we recover $A_\mathrm{SZ}=10.8_{-5.2}^{+2.3}$ and $B_\mathrm{SZ}=1.30_{-0.44}^{+0.22}$. Results in both cases imply lower cluster masses than measured in previous work with and without WL, although the uncertainties are large. The WL derived value of $B_\mathrm{SZ}$ is $\approx 20\%$ lower than the value preferred by the most recent SPT cluster cosmology analysis. The method demonstrated in this work is designed to constrain cluster masses and cosmological parameters simultaneously and will form the basis for subsequent studies that employ the full SPT cluster sample together with the DES data., Comment: Revised version in response to referee report

Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($\beta_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($\beta_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($\beta_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $\rho_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $\rho_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$)., Comment: 19 pages; in press on MNRAS. One of two companion papers on HSC-XXL surveys (Keiichi Umetsu et al., arXiv:1909.10524)

The radial acceleration relation (RAR) in galaxies describes a tight empirical scaling law between the total acceleration $g_\mathrm{tot}(r)=GM_\mathrm{tot}(, Comment: Version accepted by ApJ. Textual changes to improve clarity. Semi-analytical model predictions improved (Figures 3, 4, and 5), better agreement with data. Results and conclusions unchanged

In the past decade, our understanding of how stars and galaxies formed during the first 5 billion years after the Big Bang has been revolutionized by observations that leverage gravitational lensing by intervening masses, which act as natural cosmic telescopes to magnify background sources. Previous studies have harnessed this effect to probe the distant universe at ultraviolet, optical, infrared and millimeter wavelengths. However, strong lensing studies of young, star-forming galaxies have never extended into X-ray wavelengths, which uniquely trace high-energy phenomena. Here we report an X-ray detection of star formation in a highly magnified, strongly lensed galaxy. This lensed galaxy, seen during the first third of the history of the Universe, is a low--mass, low--metallicity starburst with elevated X-ray emission, and is a likely analog to the first generation of galaxies. Our measurements yield insight into the role that X-ray emission from stellar populations in the first generation of galaxies may play in re-ionizing the Universe. This observation paves the way for future strong lensing-assisted X-ray studies of distant galaxies reaching orders of magnitude below the detection limits of current deep fields, and previews the depths that will be attainable with future X-ray observatories., Comment: Published in Nature Astronomy

We present the first results from a galaxy population study in the highest redshift galaxy clusters identified in the 2500 deg2 South Pole Telescope Sunyaev Zel’dovich effect (SPT-SZ) survey, which is sensitive to M500 ≳ 3 × 1014 M⊙ clusters from z ∼ 0.2 out to the highest redshifts where such massive structures exist. 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 carried out a four-band imaging campaign with the Hubble and Spitzer Space Telescopes of the five z ≳ 1.4, S/NSZE > 5 clusters, that are among the rarest most massive clusters known at this redshift. All five clusters show clear overdensities of red galaxies whose colors agree with the initial cluster redshift estimates, although one (SPT-CLJ0607–4448) shows a galaxy concentration much less prominent than the others. The highest redshift cluster in this sample, SPT-CLJ0459–4947 at z ∼ 1.72, is the most distant M500 > 1014 M⊙ cluster discovered thus far through its intracluster medium, and is one of only three known clusters in this mass range at z ≳ 1.7, regardless of selection. Based on UVJ-like photometric classification of quiescent and star-forming galaxies, we find that the quiescent fraction in the cluster central regions (r/r500 < 0.7) is higher than in the field at the same redshift, with corresponding environmental quenching efficiencies typically in the range ∼0.5 − 0.8 for stellar masses log(M/M⊙) > 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/NSZE> 5, even if they had the same quiescent fraction as measured in the field. Our results thus point towards an efficient suppression of star formation in the central regions of the most massive clusters, occurring already earlier than z ∼ 1.5.

We present a statistical weak-lensing magnification analysis on an optically selected sample of 3029 \texttt{CAMIRA} galaxy clusters with richness $N>15$ at redshift $0.2\leq z, 24 pages, 16 figures. Accepted for publication in the MNRAS

We present spectroscopic confirmation of five galaxy clusters at 1.25 < z < 1.5, discovered in the 2500 deg2 South Pole Telescope Sunyaev-Zel’dovich (SZ) survey. These clusters, taken from a mass-limited sample with a nearly redshift-independent selection function, have multiwavelength follow-up imaging data from the X-ray to near-IR 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 multiobject 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 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) and emission-line ([O II] λλ3727, 3729) spectral features. Moreover, the calculated velocity dispersions yield dynamical cluster masses in good agreement with the SZ masses for these clusters. We discuss the velocity and spatial distributions of passive and [O II]-emitting galaxies in these clusters, showing that they are consistent with velocity segregation and biases observed in lower redshift South Pole Telescope clusters. We identify modest [O II] emission and pronounced CN and Hδ 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 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.

© 2018 The Author(s). Uncertainty in mass-observable scaling relations is currently the limiting factor for galaxycluster-based cosmology. Weak gravitational lensing can provide 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 = z = 0.61 and combine them with previously reported space-based observations of 13 galaxy clusters at redshifts 0.576 = z = 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 relation 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 for the same clusters.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak gravitational lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster, assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle (PA) of the projected major axis. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the 5σ significance level within a scale of 2 Mpc h. The median projected axis ratio is 0.67 ± 0.07 at a virial mass of M = (15.2 ± 2.8) × 10 M, which is in agreement with theoretical predictions from recent numerical simulations of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of |ΔPA| = 21° ± 7°. We also conduct a stacked quadrupole shear analysis of the 20 clusters assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of 0.67 ± 0.10, suggesting again a tight alignment between the intracluster gas and dark matter. © 2018. The American Astronomical Society. All rights reserved., This work was made possible by the availability of high-quality weak-lensing data produced by the CLASH survey. We are grateful to the CLASH team who enabled us to carry out this work. We thank the anonymous referee for constructive suggestions and comments. We thank Masamune Oguri for making his simulated Subaru Suprime-Cam observations available to us. K.U. acknowledges support from the Ministry of Science and Technology of Taiwan (grants MOST 103-2112-M-001-030-MY3 and MOST 106-2628-M-001-003-MY3) and from the Academia Sinica Investigator Award. M.S. and S.E. acknowledge financial support from the contracts ASI-INAF I/009/10/0, NARO15 ASI-INAF I/037/12/0, ASI 2015-046-R.0 and ASI-INAF n.2017-14-H.0. Support for D.G. was provided by NASA through Einstein Postdoctoral Fellowship grant number PF5-160138 awarded by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for NASA under contract NAS8-03060. T.O. acknowledges support from the Ministry of Science and Technology of Taiwan under the grant MOST 106-2119-M-001-031-MY3. M.M., M.S., S.E., and J.S. acknowledge support from the Italian Ministry of Foreign Affairs and International Cooperation, Directorate General for Country Promotion (Project "Crack the lens"). J.S. was supported by NSF/AST-1617022.

We perform a three-dimensional triaxial analysis of 16 X-ray regular and 4 high-magnification galaxy clusters selected from the CLASH survey by combining two-dimensional weak-lensing and central strong-lensing constraints. In a Bayesian framework, we constrain the intrinsic structure and geometry of each individual cluster assuming a triaxial Navarro-Frenk-White halo with arbitrary orientations, characterized by the mass $M_{200\mathrm{c}}$, halo concentration $C_{200\mathrm{c}}$, and triaxial axis ratios ($q_{\mathrm{a}} \le q_{\mathrm{b}}$), and investigate scaling relations between these halo structural parameters. From triaxial modeling of the X-ray-selected subsample, we find that the halo concentration decreases with increasing cluster mass, with a mean concentration of $C_{200\mathrm{c}} = 4.82\pm0.30$ at the pivot mass $M_{200\mathrm{c}}=10^{15}M_{\odot}h^{-1}$. This is consistent with the result from spherical modeling, $C_{200\mathrm{c}}=4.51\pm 0.14$. Independently of the priors, the minor-to-major axis ratio $q_{\mathrm{a}}$ of our full sample exhibits a clear deviation from the spherical configuration ($q_{\mathrm{a}}=0.52 \pm 0.04$ at $10^{15}M_{\odot}h^{-1}$ with uniform priors), with a weak dependence on the cluster mass. Combining all 20 clusters, we obtain a joint ensemble constraint on the minor-to-major axis ratio of $q_{\mathrm{a}}=0.652^{+0.162}_{-0.078}$ and a lower bound on the intermediate-to-major axis ratio of $q_{\mathrm{b}}>0.63$ at the $2\sigma$ level from an analysis with uniform priors. Assuming priors on the axis ratios derived from numerical simulations, we constrain the degree of triaxiality for the full sample to be $\mathcal{T}=0.79 \pm 0.03$ at $10^{15}M_{\odot}h^{-1}$, indicating a preference for a prolate geometry of cluster halos. We find no statistical evidence for an orientation bias ($f_{\mathrm{geo}}=0.93 \pm 0.07$) (abridged), Comment: Accepted by the ApJ. This is one of the three companion papers, including Umetsu et al. 2018 (arXiv:1804.00664) and Sereno et al. 2018 (arXiv:1804.00667), in the CLUMP-3D project. Go go go, HoHo!

The ΛCDM model of structure formation makes strong predictions on the concentration and shape of dark matter (DM) halos, which are determined by mass accretion processes. Comparison between predicted shapes and observations provides a geometric test of the ΛCDM model. Accurate and precise measurements needs a full three-dimensional (3D) analysis of the cluster mass distribution. We accomplish this with a multi-probe 3D analysis of the X-ray regular Cluster Lensing and Supernova survey with Hubble (CLASH) clusters combining strong and weak lensing, X-ray photometry and spectroscopy, and the Sunyaev–Zel'dovich effect (SZe). The cluster shapes and concentrations are consistent with ΛCDM predictions. The CLASH clusters are randomly oriented, as expected given the sample selection criteria. Shapes agree with numerical results for DM-only halos, which hints at baryonic physics being less effective in making halos rounder.

We estimate total mass ($M_{500}$), intracluster medium (ICM) mass ($M_{\mathrm{ICM}}$) and stellar mass ($M_{\star}$) in a Sunyaev-Zel'dovich effect (SZE) selected sample of 91 galaxy clusters with masses $M_{500}\gtrsim2.5\times10^{14}M_{\odot}$ and redshift $0.2 < z < 1.25$ from the 2500 deg$^2$ South Pole Telescope SPT-SZ survey. The total masses $M_{500}$ are estimated from the SZE observable, the ICM masses $M_{\mathrm{ICM}}$ are obtained from the analysis of $Chandra$ X-ray observations, and the stellar masses $M_{\star}$ are derived by fitting spectral energy distribution templates to Dark Energy Survey (DES) $griz$ optical photometry and $WISE$ or $Spitzer$ near-infrared photometry. We study trends in the stellar mass, the ICM mass, the total baryonic mass and the cold baryonic fraction with cluster mass and redshift. We find significant departures from self-similarity in the mass scaling for all quantities, while the redshift trends are all statistically consistent with zero, indicating that the baryon content of clusters at fixed mass has changed remarkably little over the past $\approx9$ Gyr. We compare our results to the mean baryon fraction (and the stellar mass fraction) in the field, finding that these values lie above (below) those in cluster virial regions in all but the most massive clusters at low redshift. Using a simple model of the matter assembly of clusters from infalling groups with lower masses and from infalling material from the low density environment or field surrounding the parent halos, we show that the measured mass trends without strong redshift trends in the stellar mass scaling relation could be explained by a mass and redshift dependent fractional contribution from field material. Similar analyses of the ICM and baryon mass scaling relations provide evidence for the so-called "missing baryons" outside cluster virial regions., Comment: Accepted by MNRAS with further clarification and discussion of the mass calibration added

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology., Comment: 12 pages, 2 tables, 7 figures, accepted to PASJ special issue

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.115$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs., Comment: 18 pages, 18 figures, accepted for publication in PASJ; cluster catalogs are available at http://www.slac.stanford.edu/~oguri/cluster/

Multi-wavelength techniques can probe the distribution and the physical properties of baryons and dark matter in galaxy clusters from the inner regions out to the peripheries. We present a full three-dimensional analysis combining strong and weak lensing, X-ray surface brightness and temperature, and the Sunyaev-Zel'dovich effect. The method is applied to MACS J1206.2-0847, a remarkably regular, face-on, massive, M_{200}=(1.1+-0.2)*10^{15}M_Sun/h, cluster at z=0.44. The measured concentration, c_{200}=6.3+-1.2, and the triaxial shape are common to halos formed in a LCDM scenario. The gas has settled in and follows the shape of the gravitational potential, which is evidence of pressure equilibrium via the shape theorem. There is no evidence for significant non-thermal pressure and the equilibrium is hydrostatic., Comment: 26 pages, in press on MNRAS

The velocity distribution of galaxies in clusters is not universal; rather, galaxies are segregated according to their spectral type and relative luminosity. We examine the velocity distributions of different populations of galaxies within 89 Sunyaev Zel'dovich (SZ) selected galaxy clusters spanning $ 0.28 < z < 1.08$. Our sample is primarily draw from the SPT-GMOS spectroscopic survey, supplemented by additional published spectroscopy, resulting in a final spectroscopic sample of 4148 galaxy spectra---2868 cluster members. The velocity dispersion of star-forming cluster galaxies is $17\pm4$% greater than that of passive cluster galaxies, and the velocity dispersion of bright ($m < m^{*}-0.5$) cluster galaxies is $11\pm4$% lower than the velocity dispersion of our total member population. We find good agreement with simulations regarding the shape of the relationship between the measured velocity dispersion and the fraction of passive vs. star-forming galaxies used to measure it, but we find a small offset between this relationship as measured in data and simulations in which suggests that our dispersions are systematically low by as much as 3\% relative to simulations. We argue that this offset could be interpreted as a measurement of the effective velocity bias that describes the ratio of our observed velocity dispersions and the intrinsic velocity dispersion of dark matter particles in a published simulation result. Measuring velocity bias in this way suggests that large spectroscopic surveys can improve dispersion-based mass-observable scaling relations for cosmology even in the face of velocity biases, by quantifying and ultimately calibrating them out., Accepted to ApJ; 21 pages, 11 figures, 5 tables

The completeness of source detection is critical for analyzing the photometric and spatial properties of the population of interest observed by astronomical imaging. We present a software package $\mathtt{ComEst}$, which calculates the completeness of source detection on charge-coupled device (CCD) images of astronomical observations, especially for the optical and near-infrared (NIR) imaging of galaxies and point sources. The completeness estimator $\mathtt{ComEst}$ is designed for the source finder $\mathtt{SExtractor}$ used on the CCD images saved in the Flexible Image Transport System (FITS) format. Specifically, $\mathtt{ComEst}$ estimates the completeness of the source detection by deriving the detection rate of synthetic point sources and galaxies simulated on the observed CCD images. In order to capture any observational artifacts or noise properties while deriving the completeness, $\mathtt{ComEst}$ directly carries out the detection of simulated sources on the observed images. Given an observed CCD image saved in FITS format, $\mathtt{ComEst}$ derives the completeness of the source detection from end to end as a function of source flux (or magnitude) and CCD position. In addition, $\mathtt{ComEst}$ can also estimate the purity of the source detection by comparing the catalog of the detected sources to the input catalogs of the simulated sources. We run ComEst on the images from Blanco Cosmology Survey (BCS) and compare the derived completeness as a function of magnitude to the limiting magnitudes derived by using the Signal-to-Noise ratio (SNR) and number count histogram of the detected sources. $\mathtt{ComEst}$ is released as a Python package with an easy-to-use syntax and is publicly available at https://github.com/inonchiu/ComEst, This paper is dedicated to Chien-Ho Lin in Taiwan. Accepted for publication in Astronomy and Computing

We study the overdensity of point sources in the direction of X-ray-selected galaxy clusters from the meta-catalogue of X-ray-detected clusters of galaxies (MCXC; < z > = 0.14) at South Pole Telescope (SPT) and Sydney University Molonglo Sky Survey (SUMSS) frequencies. Flux densities at 95, 150 and 220 GHz are extracted from the 2500 deg(2) SPT-SZ survey maps at the locations of SUMSS sources, producing a multifrequency catalogue of radio galaxies. In the direction of massive galaxy clusters, the radio galaxy flux densities at 95 and 150 GHz are biased low by the cluster Sunyaev-Zel'dovich Effect (SZE) signal, which is negative at these frequencies. We employ a cluster SZE model to remove the expected flux bias and then study these corrected source catalogues. We find that the high-frequency radio galaxies are centrally concentrated within the clusters and that their luminosity functions (LFs) exhibit amplitudes that are characteristically an order of magnitude lower than the cluster LF at 843 MHz. We use the 150 GHz LF to estimate the impact of cluster radio galaxies on an SPT-SZ like survey. The radio galaxy flux typically produces a small bias on the SZE signal and has negligible impact on the observed scatter in the SZE mass-observable relation. If we assume there is no redshift evolution in the radio galaxy LF then 1.8 +/- 0.7 per cent of the clusters with detection significance xi >= 4.5 would be lost from the sample. Allowing for redshift evolution of the form ( 1 + z) (2.5) increases the incompleteness to 5.6 +/- 1.0 per cent. Improved constraints on the evolution of the cluster radio galaxy LF require a larger cluster sample extending to higher redshift., Gordon and Betty Moore Foundation (Grant GBMF-947)

We study the galaxy populations in 74 Sunyaev Zeldovich Effect (SZE) selected clusters from the South Pole Telescope (SPT) survey that have been imaged in the science verification phase of the Dark Energy Survey (DES). The sample extends up to $z\sim 1.1$ with $4 \times 10^{14} M_{\odot}\le M_{200}\le 3\times 10^{15} M_{\odot}$. Using the band containing the 4000~\AA\ break and its redward neighbor, we study the color-magnitude distributions of cluster galaxies to $\sim m_*+2$, finding: (1) the intrinsic rest frame $g-r$ color width of the red sequence (RS) population is $\sim$0.03 out to $z\sim0.85$ with a preference for an increase to $\sim0.07$ at $z=1$ and (2) the prominence of the RS declines beyond $z\sim0.6$. The spatial distribution of cluster galaxies is well described by the NFW profile out to $4R_{200}$ with a concentration of $c_{\mathrm{g}} = 3.59^{+0.20}_{-0.18}$, $5.37^{+0.27}_{-0.24}$ and $1.38^{+0.21}_{-0.19}$ for the full, the RS and the blue non-RS populations, respectively, but with $\sim40$\% to 55\% cluster to cluster variation and no statistically significant redshift or mass trends. The number of galaxies within the virial region $N_{200}$ exhibits a mass trend indicating that the number of galaxies per unit total mass is lower in the most massive clusters, and shows no significant redshift trend. The red sequence (RS) fraction within $R_{200}$ is $(68\pm3)$\% at $z=0.46$, varies from $\sim$55\% at $z=1$ to $\sim$80\% at $z=0.1$, and exhibits intrinsic variation among clusters of $\sim14$\%. We discuss a model that suggests the observed redshift trend in RS fraction favors a transformation timescale for infalling field galaxies to become RS galaxies of 2 to 3~Gyr., Comment: 22 pages, 14 figures, submitted to MNRAS

(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.

We present measurements of the X-ray observables of the intra-cluster medium (ICM), including luminosity $L_X$, ICM mass $M_{ICM}$, emission-weighted mean temperature $T_X$, and integrated pressure $Y_X$, that are derived from XMM-Newton X-ray observations of a Sunyaev-Zel'dovich Effect (SZE) selected sample of 59 galaxy clusters from the South Pole Telescope SPT-SZ survey that span the redshift range of $0.20 < z < 1.5$. We constrain the best-fit power law scaling relations between X-ray observables, redshift, and halo mass. The halo masses are estimated based on previously published SZE observable to mass scaling relations, calibrated using information that includes the halo mass function. Employing SZE-based masses in this sample enables us to constrain these scaling relations for massive galaxy clusters ($M_{500}\geq 3 \times10^{14}$ $M_\odot$) to the highest redshifts where these clusters exist without concern for X-ray selection biases. We find that the mass trends are steeper than self-similarity in all cases, and with $\geq 2.5{\sigma}$ significance in the case of $L_X$ and $M_{ICM}$. The redshift trends are consistent with the self-similar expectation, but the uncertainties remain large. Core-included scaling relations tend to have steeper mass trends for $L_X$. There is no convincing evidence for a redshift-dependent mass trend in any observable. The constraints on the amplitudes of the fitted scaling relations are currently limited by the systematic uncertainties on the SZE-based halo masses, however the redshift and mass trends are limited by the X-ray sample size and the measurement uncertainties of the X-ray observables., Comment: Accepted for publication in ApJ. 29 pages, 11 figures, 5 tables

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., 24 pages in eapj format. Accepted to the Astrophysical Journal Supplements

We present the results of an X-ray spectral analysis of 153 galaxy clusters observed with the Chandra, XMM-Newton, and Suzaku space telescopes. These clusters, which span 0 < z < 1.5, were drawn from a larger, mass-selected sample of galaxy clusters discovered in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. With a total combined exposure time of 9.1 Ms, these data yield the strongest constraints to date on the evolution of the metal content of the intracluster medium (ICM). We find no evidence for strong evolution in the global (r = 0.23 +/- 0.01 Zsun and a scatter of 0.08 +/- 0.01 Zsun. These results imply that >60% of the metals in the ICM were already in place at z=1 (at 95% confidence), consistent with the picture of an early (z>1) enrichment. We find, in agreement with previous works, a significantly higher mean value for the metallicity in the centers of cool core clusters versus non-cool core clusters. We find weak evidence for evolution in the central metallicity of cool core clusters (dZ/dz = -0.21 +/- 0.11 Zsun), which is sufficient to account for this enhanced central metallicity over the past ~10 Gyr. We find no evidence for metallicity evolution outside of the core (dZ/dz = -0.03 +/- 0.06 Zsun), and no significant difference in the core-excised metallicity between cool core and non-cool core clusters. This suggests that strong radio-mode AGN feedback does not significantly alter the distribution of metals at r>0.15R500. Given the limitations of current-generation X-ray telescopes in constraining the ICM metallicity at z>1, significant improvements on this work will likely require next-generation X-ray missions., Comment: 11 pages, 8 figures, 2 tables. Submitted to ApJ. Comments welcome!

We present a study of the optical properties of the 26 most massive galaxy clusters selected within the SPT-SZ 2500 deg$^2$ survey. This Sunyaev-Zel'dovich effect selected sample spans a redshift range of 0.10 < z < 1.13. We measure the galaxy radial profile, the luminosity function (LF), and the halo occupation number (HON) using optical data with a typical depth of $m^*$ + 2. The stacked radial profiles are consistent with a NFW profile with a concentration of $2.84^{+0.40}_{-0.37}$ for the red sequence (RS) and $2.36^{+0.38}_{-0.35}$ for the total population. Stacking the data in multiple redshift bins shows a hint of redshift evolution in the concentration when both the total population is used, and when only RS galaxies are used (at 2.1$\sigma$ and 2.8$\sigma$, respectively). The stacked LF shows a faint end slope $\alpha = -1.06^{+0.04}_{-0.03}$ for the total and $\alpha = -0.80^{+0.04}_{-0.03}$ for the RS population. The redshift evolution of $m^*$ is found to be consistent with a passively evolving Composite Stellar Population (CSP) model. By adopting the CSP model predictions, we explore the redshift evolution of the schechter parameters $\alpha$ and $\phi^*$. We find $\alpha$ for the total population to be consistent with no evolution (0.3$\sigma$), while evidence of evolution for the red galaxies is mildly significant (1.1-2.1$\sigma$). The data show that the density $\phi^*$/E$^2$(z) decreases with redshift, in tension with the self-similar expectation at a 2.4$\sigma$ level for the total population. The measured HON-mass relation has a lower normalization than previous studies at low redshift. Finally, our data support HON redshift evolution at a 2.1$\sigma$ level, with clusters at higher redshift containing fewer galaxies per unit mass to $m^*$ + 3 than their low-z counterparts [abridged]., Comment: 14 pages, 7 figures, 2 tables, submitted to MNRAS, abstract abridged for arXiv submission

We study the stellar, brightest cluster galaxy (BCG) and intracluster medium (ICM) masses of 14 South Pole Telescope (SPT) selected galaxy clusters with median redshift z = 0.9 and mass M-500 = 6 x 10(14) M-circle dot. We estimate stellar masses for each cluster and BCG using six photometric bands, the ICM mass using X-ray observations and the virial masses using the SPT Sunyaev-Zel'dovich effect signature. At z = 0.9, the BCG mass M-*(BCG) constitutes 0.12 +/- 0.01 per cent of the halo mass for a 6 x 10(14) M-circle dot cluster, and this fraction falls as M-500(-0.58 +/- 0.07). The cluster stellar mass function has a characteristic mass M-0 = 10(11.0 +/- 0.1) M-circle dot, and the number of galaxies per unit mass in clusters is larger than in the field by a factor of 1.65 +/- 0.20. We combine our SPT sample with previously published samples at low redshift and correct to a common initial mass function and for systematic virial mass differences. We then explore mass and redshift trends in the stellar fraction f(*), the ICM fraction f(ICM), the collapsed baryon fraction f(c) and the baryon fraction f(b). At a pivot mass of 6 x 10(14) M-circle dot and redshift z = 0.9, the characteristic values are f(*) = 1.1 +/- 0.1 per cent, f(ICM) = 9.6 +/- 0.5 per cent, f(c) = 10.7 +/- 1.1 per cent and f(b) = 10.7 +/- 0.6 per cent. These fractions all vary with cluster mass at high significance, with higher mass clusters having lower f(*) and f(c) and higher f(ICM) and f(b). When accounting for a 15 per cent systematic virial mass uncertainty, there is no statistically significant redshift trend at fixed mass. Our results support the scenario where clusters grow through accretion from subclusters (higher f(*), lower f(ICM)) and the field (lower f(*), higher f(ICM)), balancing to keep f(*) and f(ICM) approximately constant since z similar to 0.9.

We detect the kinematic Sunyaev-Zel'dovich (kSZ) effect with a statistical significance of $4.2 \sigma$ by combining a cluster catalogue derived from the first year data of the Dark Energy Survey (DES) with CMB temperature maps from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey. This measurement is performed with a differential statistic that isolates the pairwise kSZ signal, providing the first detection of the large-scale, pairwise motion of clusters using redshifts derived from photometric data. By fitting the pairwise kSZ signal to a theoretical template we measure the average central optical depth of the cluster sample, $\bar{\tau}_e = (3.75 \pm 0.89)\cdot 10^{-3}$. We compare the extracted signal to realistic simulations and find good agreement with respect to the signal-to-noise, the constraint on $\bar{\tau}_e$, and the corresponding gas fraction. High-precision measurements of the pairwise kSZ signal with future data will be able to place constraints on the baryonic physics of galaxy clusters, and could be used to probe gravity on scales $ \gtrsim 100$ Mpc., Comment: 23 pages, 14 figures; matches version published in MNRAS

We present a detection of the enhancement in the number densities of background galaxies induced from lensing magnification and use it to test the Sunyaev-Zel'dovich effect (SZE) inferred masses in a sample of 19 galaxy clusters with median redshift $z\simeq0.42$ selected from the South Pole Telescope SPT-SZ survey. Two background galaxy populations are selected for this study through their photometric colours; they have median redshifts ${z}_{\mathrm{median}}\simeq0.9$ (low-$z$ background) and ${z}_{\mathrm{median}}\simeq1.8$ (high-$z$ background). Stacking these populations, we detect the magnification bias effect at $3.3\sigma$ and $1.3\sigma$ for the low- and high-$z$ backgrounds, respectively. We fit NFW models simultaneously to all observed magnification bias profiles to estimate the multiplicative factor $\eta$ that describes the ratio of the weak lensing mass to the mass inferred from the SZE observable-mass relation. We further quantify systematic uncertainties in $\eta$ resulting from the photometric noise and bias, the cluster galaxy contamination and the estimations of the background properties. The resulting $\eta$ for the combined background populations with $1\sigma$ uncertainties is $0.83\pm0.24\mathrm{(stat)}\pm0.074\mathrm{(sys)}$, indicating good consistency between the lensing and the SZE-inferred masses. We use our best-fit $\eta$ to predict the weak lensing shear profiles and compare these predictions with observations, showing agreement between the magnification and shear mass constraints. This work demonstrates the promise of using the magnification as a complementary method to estimate cluster masses in large surveys., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

We present the stellar mass-halo mass scaling relation for 46 X-ray selected low-mass clusters or groups detected in the XMM-BCS survey with masses $2\times10^{13}M_{\odot}\lesssim M_{500}\lesssim2.5\times10^{14}M_{\odot}$ at redshift $0.1\le z \le1.02$. The cluster binding masses $M_{500}$ are inferred from the measured X-ray luminosities \Lx, while the stellar masses $M_{\star}$ of the galaxy populations are estimated using near-infrared imaging from the SSDF survey and optical imaging from the BCS survey. With the measured \Lx\ and stellar mass $M_{\star}$, we determine the best fit stellar mass-halo mass relation, accounting for selection effects, measurement uncertainties and the intrinsic scatter in the scaling relation. The resulting mass trend is $M_{\star}\propto M_{500}^{0.69\pm0.15}$, the intrinsic (log-normal) scatter is $\sigma_{\ln M_{\star}|M_{500}}=0.36^{+0.07}_{-0.06}$, and there is no significant redshift trend $M_{\star}\propto (1+z)^{-0.04\pm0.47}$, although the uncertainties are still large. We also examine $M_{\star}$ within a fixed projected radius of $0.5$~Mpc, showing that it provides a cluster binding mass proxy with intrinsic scatter of $\approx93\%$ (1$\sigma$ in $M_{500}$). We compare our $M_{\star}=M_{\star}(M_{500}, z)$ scaling relation from the XMM-BCS clusters with samples of massive, SZE-selected clusters ($M_{500}\approx6\times10^{14}M_{\odot}$) and low mass NIR-selected clusters ($M_{500}\approx10^{14}M_{\odot}$) at redshift $0.6\lesssim z \lesssim1.3$. After correcting for the known mass measurement systematics in the compared samples, we find that the scaling relation is in good agreement with the high redshift samples, suggesting that for both groups and clusters the stellar content of the galaxy populations within $R_{500}$ depends strongly on mass but only weakly on redshift out to $z\approx1$., Comment: 15 pages, 10 figures. Accepted for publication in MNRAS

We cross-match galaxy cluster candidates selected via their Sunyaev–Zel'dovich effect (SZE) signatures in 129.1 deg[superscript 2] of the South Pole Telescope 2500d SPT-SZ survey with optically identified clusters selected from the Dark Energy Survey science verification data. We identify 25 clusters between 0.1 ≲ z ≲ 0.8 in the union of the SPT-SZ and redMaPPer (RM) samples. RM is an optical cluster finding algorithm that also returns a richness estimate for each cluster. We model the richness λ-mass relation with the following function 〈ln λ|M500〉 ∝ Bλln M500 + Cλln E(z) and use SPT-SZ cluster masses and RM richnesses λ to constrain the parameters. We find Bλ=1.14[superscript +0.21][subscript −0.18 and Cλ=0.73[superscript +0.77][subscript −0.75]. The associated scatter in mass at fixed richness is σlnM|λ=0.18[superscript +0.08][subscript −0.05] at a characteristic richness λ = 70. We demonstrate that our model provides an adequate description of the matched sample, showing that the fraction of SPT-SZ-selected clusters with RM counterparts is consistent with expectations and that the fraction of RM-selected clusters with SPT-SZ counterparts is in mild tension with expectation. We model the optical-SZE cluster positional offset distribution with the sum of two Gaussians, showing that it is consistent with a dominant, centrally peaked population and a subdominant population characterized by larger offsets. We also cross-match the RM catalogue with SPT-SZ candidates below the official catalogue threshold significance ξ = 4.5, using the RM catalogue to provide optical confirmation and redshifts for 15 additional clusters with ξ ∈ [4, 4.5]., Kavli Foundation, Gordon and Betty Moore Foundation (Grant GBMF 947)

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)., Comment: 14 pages, 3 figures. Published in ApJ. Replaced to match published version

We present a velocity-dispersion-based mass calibration of the South Pole Telescope Sunyaev-Zel'dovich effect survey (SPT-SZ) galaxy cluster sample. Using a homogeneously selected sample of 100 cluster candidates from 720 deg^2 of the survey along with 63 velocity dispersion (σ_v) and 16 X-ray Y_X measurements of sample clusters, we simultaneously calibrate the mass-observable relation and constrain cosmological parameters. Our method accounts for cluster selection, cosmological sensitivity, and uncertainties in the mass calibrators. The calibrations using σ_v and Y_X are consistent at the 0.6σ level, with the σ_v calibration preferring ~16% higher masses. We use the full SP_(TCL) data set (SZ clusters+σ_v +Y_X) to measure σ_g(Ω_m/0.27)^(0.3) = 0.809 ± 0.036 within a flat ΛCDM model. The SPT cluster abundance is lower than preferred by either the WMAP9 or Planck+WMAP9 polarization (WP) data, but assuming that the sum of the neutrino masses is ∑m_ν = 0.06 eV, we find the data sets to be consistent at the 1.0σ level for WMAP9 and 1.5σ for Planck+WP. Allowing for larger ∑m_ν further reconciles the results. When we combine the SPT_(CL) and Planck+WP data sets with information from baryon acoustic oscillations and Type Ia supernovae, the preferred cluster masses are 1.9σ higher than the Y_X calibration and 0.8σ higher than the σ_v calibration. Given the scale of these shifts (~44% and ~23% in mass, respectively), we execute a goodness-of-fit test; it reveals no tension, indicating that the best-fit model provides an adequate description of the data. Using the multi-probe data set, we measure Ω_m = 0.299 ± 0.009 and σ_g = 0.829 ± 0.011. Within a νCDM model we find ∑m_ν = 0.148 ± 0.081 eV. We present a consistency test of the cosmic growth rate using SPT clusters. Allowing both the growth index γ and the dark energy equation-of-state parameter w to vary, we find γ = 0.73 ± 0.28 and w = –1.007 ± 0.065, demonstrating that the expansion and the growth histories are consistent with a ΛCDM universe (γ = 0.55; w = –1).

(Abridged) We use 95, 150, and 220GHz observations from the SPT to examine the SZE signatures of a sample of 46 X-ray selected groups and clusters drawn from ~6 deg^2 of the XMM-BCS. These systems extend to redshift z=1.02, have characteristic masses ~3x lower than clusters detected directly in the SPT data and probe the SZE signal to the lowest X-ray luminosities (>10^42 erg s^-1) yet. We develop an analysis tool that combines the SZE information for the full ensemble of X-ray-selected clusters. Using X-ray luminosity as a mass proxy, we extract selection-bias corrected constraints on the SZE significance- and Y_500-mass relations. The SZE significance- mass relation is in good agreement with an extrapolation of the relation obtained from high mass clusters. However, the fit to the Y_500-mass relation at low masses, while in good agreement with the extrapolation from high mass SPT clusters, is in tension at 2.8 sigma with the constraints from the Planck sample. We examine the tension with the Planck relation, discussing sample differences and biases that could contribute. We also present an analysis of the radio galaxy point source population in this ensemble of X-ray selected systems. We find 18 of our systems have 843 MHz SUMSS sources within 2 arcmin of the X-ray centre, and three of these are also detected at significance >4 by SPT. Of these three, two are associated with the group brightest cluster galaxies, and the third is likely an unassociated quasar candidate. We examine the impact of these point sources on our SZE scaling relation analyses and find no evidence of biases. We also examine the impact of dusty galaxies using constraints from the 220 GHz data. The stacked sample provides 2.8$��$ significant evidence of dusty galaxy flux, which would correspond to an average underestimate of the SPT Y_500 signal that is (17+-9) per cent in this sample of low mass systems., 15 pages, 7 figures

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}$., Comment: Submitted to MNRAS Letters

Since the discovery of the "bullet cluster" several similar cases have been uncovered suggesting relative velocities well beyond the tail of high speed collisions predicted by the concordance LCDM model. However, quantifying such post-merger events with hydrodynamical models requires a wide coverage of possible initial conditions. Here we show that it is simpler to interpret pre-merger cases, such as A1750, where the gas between the colliding clusters is modestly affected, so that the initial conditions are clear. We analyze publicly available Chandra data confirming a significant increase in the projected X-ray temperature between the two cluster centers in A1750 consistent with our expectations for a merging cluster. We model this system with a self-consistent hydrodynamical simulation of dark matter and gas using the FLASH code. Our simulations reproduce well the X-ray data, and the measured redshift difference between the two clusters in the phase before the first core passage viewed at an intermediate projection angle. The deprojected initial relative velocity derived using our model is 1460 km/sec which is considerably higher than the predicted mean impact velocity for simulated massive haloes derived by recent LCDM cosmological simulations, but it is within the allowed range. Our simulations demonstrate that such systems can be identified using a multi-wavelength approach and numerical simulations, for which the statistical distribution of relative impact velocities may provide a definitive examination of a broad range of dark matter scenarios., Comment: 11 pages, 10 Figures and 1 Table, accepted for publication in the Astrophysical Journal

We test the assumption of strict hydrostatic equilibrium in galaxy cluster MS2137.3-2353 (MS 2137) using the latest CHANDRA X-ray observations and results from a combined strong and weak lensing analysis based on optical observations. We deproject the two-dimensional X-ray surface brightness and mass surface density maps assuming spherical and spheroidal dark matter distributions. We find a significant, 40%-50%, contribution from non-thermal pressure in the core assuming a spherical model. This non-thermal pressure support is similar to what was found by Molnar et al. (2010) using a sample of massive relaxed clusters drawn from high resolution cosmological simulations. We have studied hydrostatic equilibrium in MS 2137 under the assumption of elliptical cluster geometry adopting prolate models for the dark matter density distribution with different axis ratios. Our results suggest that the main effect of ellipticity (compared to spherical models) is to decrease the non-thermal pressure support required for equilibrium at all radii without changing the distribution qualitatively. We find that a prolate model with an axis ratio of 1.25 (axis in the line of sight over perpendicular to it) provides a physically acceptable model implying that MS 2137 is close to hydrostatic equilibrium at about 0.04-0.15 Rvir and have an about 25% contribution from non-thermal pressure at the center. Our results provide further evidence that there is a significant contribution from non-thermal pressure in the core region of even relaxed clusters, i.e., the assumption of hydrostatic equilibrium is not valid in this region, independently of the assumed shape of the cluster., 11 pages, 4 figures, accepted for publication in ApJ

We present a multi-wavelength study of 90 brightest cluster galaxies (BCGs) in a sample of galaxy clusters selected via the Sunyaev Zel'dovich effect by the South Pole Telescope, utilizing data from various ground- and space-based facilities. We infer the star formation rate (SFR) for the BCG in each cluster, based on the UV and IR continuum luminosity, as well as the [O II] emission line luminosity in cases where spectroscopy is available, finding 7 systems with SFR > 100 Msun/yr. We find that the BCG SFR exceeds 10 Msun/yr in 31 of 90 (34%) cases at 0.25 < z < 1.25, compared to ~1-5% at z ~ 0 from the literature. At z > 1, this fraction increases to 92(+6)(-31)%, implying a steady decrease in the BCG SFR over the past ~9 Gyr. At low-z, we find that the specific star formation rate in BCGs is declining more slowly with time than for field or cluster galaxies, most likely due to the replenishing fuel from the cooling ICM in relaxed, cool core clusters. At z > 0.6, the correlation between cluster central entropy and BCG star formation - which is well established at z ~ 0 - is not present. Instead, we find that the most star-forming BCGs at high-z are found in the cores of dynamically unrelaxed clusters. We investigate the rest-frame near-UV morphology of a subsample of the most star-forming BCGs using data from the Hubble Space Telescope, finding complex, highly asymmetric UV morphologies on scales as large as ~50-60 kpc. The high fraction of star-forming BCGs hosted in unrelaxed, non-cool core clusters at early times suggests that the dominant mode of fueling star formation in BCGs may have recently transitioned from galaxy-galaxy interactions to ICM cooling., Comment: 20 pages, 10 figures. Submitted for publication in ApJ. Comments welcome