We present the weak lensing mass calibration of the stellar mass based $\mu_{\star}$ mass proxy for redMaPPer galaxy clusters in the Dark Energy Survey Year 1. For the first time we are able to perform a calibration of $\mu_{\star}$ at high redshifts, $z>0.33$. In a blinded analysis, we use $\sim 6,000$ clusters split into 12 subsets spanning the ranges $0.1 \leqslant z<0.65$ and $\mu_{\star}$ up to $\sim 5.5 \times 10^{13} M_{\odot}$, and infer the average masses of these subsets through modelling of their stacked weak lensing signal. In our model we account for the following sources of systematic uncertainty: shear measurement and photometric redshift errors, miscentring, cluster-member contamination of the source sample, deviations from the NFW halo profile, halo triaxiality and projection effects. We use the inferred masses to estimate the joint mass--$\mu_{\star}$--$z$ scaling relation given by $\langle M_{200c} | \mu_{\star},z \rangle = M_0 (\mu_{\star}/5.16\times 10^{12} \mathrm{M_{\odot}})^{F_{\mu_{\star}}} ((1+z)/1.35)^{G_z}$. We find $M_0= (1.14 \pm 0.07) \times 10^{14} \mathrm{M_{\odot}}$ with $F_{\mu_{\star}}= 0.76 \pm 0.06$ and $G_z= -1.14 \pm 0.37$. We discuss the use of $\mu_{\star}$ as a complementary mass proxy to the well-studied richness $\lambda$ for: $i)$ exploring the regimes of low $z$, $\lambda<20$ and high $\lambda$, $z \sim 1$; $ii)$ testing systematics such as projection effects for applications in cluster cosmology.

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

We present here the weak gravitational lensing detection of four nearby galaxy clusters in the southern sky: Abell 2029, Abell 85, Abell 1606 and Abell 2457. The weak lensing detections of Abell 1606 and Abell 2457 are the first in the literature. This work capitalizes on the wide field of view of the Dark Energy Camera at the Cerro Tololo Inter-American Observatory, which we use to obtain deep, multi-wavelength imaging of all targets. We publish maps of the clusters' projected mass distributions, and obtain the $M_{200}$ of their clusters through NFW profile fits to the two-dimensional tangential ellipticity signal.

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

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

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$\sigma$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters.

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.

This paper presents multiwavelength photometric catalogues of the last two Hubble Frontier Fields (HFF), the massive galaxy clusters Abell 370 and RXC J2248.7-4431. The photometry ranges from imaging performed on the Hubble Space Telescope (HST) to ground based Very Large Telescope (VLT) and Spitzer/IRAC, in collaboration with the ASTRODEEP team, and using the ASTRODEEP pipeline. While the main purpose of this paper is to release the catalogues, we also perform, as a proof of concept, a brief analysis of z > 6 objects selected using drop-out method, as well as spectroscopically confirmed sources and multiple images in both clusters. While dropout methods yield a sample of high-z galaxies, the addition of longer wavelength data reveals that as expected the samples have substantial contamination at the ~30-45% level by dusty galaxies at lower redshifts. Furthermore, we show that spectroscopic redshifts are still required to unambiguously determine redshifts of multiply imaged systems. Finally, the now publicly available ASTRODEEP catalogues were combined for all HFFs and used to explore stellar properties of a large sample of 20,000 galaxies across a large photometric redshift range. The powerful magnification provided by the HFF clusters allows us an exploration of the properties of galaxies with intrinsic stellar masses as low as $M_* \gtrsim 10^7M_{\odot}$ and intrinsic star formation rates $\mbox{SFRs}\sim 0.1\mbox{-}1M_\odot/\mbox yr$ at z > 6.

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

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

Community access to deep (i ~ 25), highly-multiplexed optical and near-infrared multi-object spectroscopy (MOS) on 8-40m telescopes would greatly improve measurements of cosmological parameters from LSST. The largest gain would come from improvements to LSST photometric redshifts, which are employed directly or indirectly for every major LSST cosmological probe; deep spectroscopic datasets will enable reduced uncertainties in the redshifts of individual objects via optimized training. Such spectroscopy will also determine the relationship of galaxy SEDs to their environments, key observables for studies of galaxy evolution. The resulting data will also constrain the impact of blending on photo-z's. Focused spectroscopic campaigns can also improve weak lensing cosmology by constraining the intrinsic alignments between the orientations of galaxies. Galaxy cluster studies can be enhanced by measuring motions of galaxies in and around clusters and by testing photo-z performance in regions of high density. Photometric redshift and intrinsic alignment studies are best-suited to instruments on large-aperture telescopes with wider fields of view (e.g., Subaru/PFS, MSE, or GMT/MANIFEST) but cluster investigations can be pursued with smaller-field instruments (e.g., Gemini/GMOS, Keck/DEIMOS, or TMT/WFOS), so deep MOS work can be distributed amongst a variety of telescopes. However, community access to large amounts of nights for surveys will still be needed to accomplish this work. In two companion white papers we present gains from shallower, wide-area MOS and from single-target imaging and spectroscopy.

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.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

Weak lensing source galaxy catalogs used in estimating the masses of galaxy clusters can be heavily contaminated by cluster members, prohibiting accurate mass calibration. In this study we test the performance of an estimator for the extent of cluster member contamination based on decomposing the photometric redshift $P(z)$ of source galaxies into contaminating and background components. We perform a full scale mock analysis on a simulated sky survey approximately mirroring the observational properties of the Dark Energy Survey Year One observations (DES Y1), and find excellent agreement between the true number profile of contaminating cluster member galaxies in the simulation and the estimated one. We further apply the method to estimate the cluster member contamination for the DES Y1 redMaPPer cluster mass calibration analysis, and compare the results to an alternative approach based on the angular correlation of weak lensing source galaxies. We find indications that the correlation based estimates are biased by the selection of the weak lensing sources in the cluster vicinity, which does not strongly impact the $P(z)$ decomposition method. Collectively, these benchmarks demonstrate the strength of the $P(z)$ decomposition method in alleviating membership contamination and enabling highly accurate cluster weak lensing studies without broad exclusion of source galaxies, thereby improving the total constraining power of cluster mass calibration via weak lensing.

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

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

The galaxy phase-space distribution in galaxy clusters provides insights into the formation and evolution of cluster galaxies, and it can also be used to measure cluster mass profiles. We present a dynamical study based on $\sim$3000 passive, non-emission line cluster galaxies drawn from 110 galaxy clusters. The galaxy clusters were selected using the Sunyaev-Zel'dovich effect (SZE) in the 2500~deg$^2$ SPT-SZ survey and cover the redshift range $0.2 < z < 1.3$. We model the clusters using the Jeans equation, while adopting NFW mass profiles and a broad range of velocity dispersion anisotropy profiles. The data prefer velocity dispersion anisotropy profiles that are approximately isotropic near the center and increasingly radial toward the cluster virial radius, and this is true for all redshifts and masses we study. The pseudo-phase-space density profile of the passive galaxies is consistent with expectations for dark matter particles and subhalos from cosmological $N$-body simulations. The dynamical mass constraints are in good agreement with external mass estimates of the SPT cluster sample from either weak lensing, velocity dispersions, or X-ray $Y_X$ measurements. However, the dynamical masses are lower (at the 2.2$\sigma$ level) when compared to the mass calibration favored when fitting the SPT cluster data to a $\Lambda$CDM model with external cosmological priors, including CMB anisotropy data from Planck. The discrepancy grows with redshift, where in the highest redshift bin the ratio of dynamical to SPT+Planck masses is $\eta=0.63^{+0.13}_{-0.08}\pm0.06$ (statistical and systematic), corresponding to a $2.6\sigma$ discrepancy.

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

The primary difficulty in measuring dynamical masses of galaxy clusters from galaxy data lies in the separation between true cluster members from interloping galaxies along the line of sight. We study the impact of membership contamination and incompleteness on cluster mass estimates obtained with 25 commonly used techniques applied to nearly 1000 mock clusters. We show that all methods overestimate or underestimate cluster masses when applied to contaminated or incomplete galaxy samples respectively. This appears to be the main source of the intrinsic scatter in the mass scaling relation. Applying corrections based on a prior knowledge of contamination and incompleteness can reduce the scatter to the level of shot noise expected for poorly sampled clusters. We establish an empirical model quantifying the effect of imperfect membership on cluster mass estimation and discuss its universal and method-dependent features. We find that both imperfect membership and the response of the mass estimators depend on cluster mass, effectively causing a flattening of the estimated - true mass relation. Imperfect membership thus alters cluster counts determined from spectroscopic surveys, hence the cosmological parameters that depend on such counts.

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

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1\sigma$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal.

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.

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

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

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

With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.

Exploiting the data of the Grism Lens-Amplified Survey from Space (GLASS), we characterize the spatial distribution of star formation in 76 high star forming galaxies in 10 clusters at 0.3< z <0.7. All these galaxies are likely restricted to first infall. In a companion paper we contrast the properties of field and cluster galaxies, whereas here we correlate the properties of H{\alpha} emitters to a number of tracers of the cluster environment to investigate its role in driving galaxy transformations. H{\alpha} emitters are found in the clusters out to 0.5 virial radii, the maximum radius covered by GLASS. The peak of the H{\alpha} emission is offset with respect to the peak of the UV-continuum. We decompose this offsets into a radial and tangential component. The radial compo- nent points away from the cluster center in 60% of the cases, with 95% confidence. The decompositions agree with cosmological simulations, i.e. the H{\alpha} emission offset correlates with galaxy velocity and ram-pressure stripping signatures. Trends between H{\alpha} emitter properties and surface mass density distributions and X-ray emissions emerge only for unrelaxed clusters. The lack of strong correlations with the global environment does not allow us to identify a unique environmental effect originating from the cluster center. In contrast, correla- tions between H{\alpha} morphology and local number density emerge. We conclude that local effects, uncorrelated to the cluster-centric radius, play a more important role in shaping galaxy properties.

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

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

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.

In Hubble Space Telescope (HST) imaging taken on 10 November 2014, four images of supernova (SN) "Refsdal" (redshift z=1.49) appeared in an Einstein-cross--like configuration (images S1-S4) around an early-type galaxy in the cluster MACS J1149.5+2223 (z=0.54). Almost all lens models of the cluster have predicted that the SN should reappear within a year in a second host-galaxy image created by the cluster's potential. In HST observations taken on 11 December 2015, we find a new source at the predicted position of the new image of SN Refsdal approximately 8" from the previous images S1-S4. This marks the first time the appearance of a SN at a particular time and location in the sky was successfully predicted in advance! We use these data and the light curve from the first four observed images of SN Refsdal to place constraints on the relative time delay and magnification of the new image (SX), compared to images S1-S4. This enables us, for the first time, to test "blind" lens model predictions of both magnifications and time delays for a lensed SN. We find that the timing and brightness of the new image are consistent with the blind predictions of a fraction of the models. The reappearance illustrates the discriminatory power of this blind test and its utility to uncover sources of systematic uncertainty. From planned HST photometry, we expect to reach a precision of 1-2% on the time delay between S1-S4 and SX.

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

We study the stellar population properties of the IRAC-detected $6 \lesssim z \lesssim 10$ galaxy candidates from the Spitzer UltRa Faint SUrvey Program (SURFS UP). Using the Lyman Break selection technique, we find a total of 16 new galaxy candidates at $6 \lesssim z \lesssim 10$ with $S/N \geq 3$ in at least one of the IRAC $3.6\mu$m and $4.5\mu$m bands. According to the best mass models available for the surveyed galaxy clusters, these IRAC-detected galaxy candidates are magnified by factors of $\sim 1.2$--$5.5$. We find that the IRAC-detected $6 \lesssim z \lesssim 10$ sample is likely not a homogeneous galaxy population: some are relatively massive (stellar mass as high as $4 \times 10^9\,M_{\odot}$) and evolved (age $\lesssim 500$ Myr) galaxies, while others are less massive ($M_{\text{stellar}}\sim 10^8\,M_{\odot}$) and very young ($\sim 10$ Myr) galaxies with strong nebular emission lines that boost their rest-frame optical fluxes. We identify two Ly$\alpha$ emitters in our sample from the Keck DEIMOS spectra, one at $z_{\text{Ly}\alpha}=6.76$ (in RXJ1347) and one at $z_{\text{Ly}\alpha}=6.32$ (in MACS0454). We show that IRAC $[3.6]-[4.5]$ color, when combined with photometric redshift, can be used to identify galaxies likely with strong nebular emission lines within certain redshift windows.

We present the first study of the spatial distribution of star formation in z~0.5 cluster galaxies. The analysis is based on data taken with the Wide Field Camera 3 as part of the Grism Lens-Amplified Survey from Space (GLASS). We illustrate the methodology by focusing on two clusters (MACS0717.5+3745 and MACS1423.8+2404) with different morphologies (one relaxed and one merging) and use foreground and background galaxies as field control sample. The cluster+field sample consists of 42 galaxies with stellar masses in the range 10^8-10^11 M_sun, and star formation rates in the range 1-20 M_sun/yr. Both in clusters and in the field, H{\alpha} is more extended than the rest-frame UV continuum in 60% of the cases, consistent with diffuse star formation and inside out growth. In ~20% of the cases, the H{\alpha} emission appears more extended in cluster galaxies than in the field, pointing perhaps to ionized gas being stripped and/or star formation being enhanced at large radii. The peak of the H{\alpha} emission and that of the continuum are offset by less than 1 kpc. We investigate trends with the hot gas density as traced by the X-ray emission, and with the surface mass density as inferred from gravitational lens models and find no conclusive results. The diversity of morphologies and sizes observed in H_alpha illustrates the complexity of the environmental process that regulate star formation. Upcoming analysis of the full GLASS dataset will increase our sample size by almost an order of magnitude, verifying and strengthening the inference from this initial dataset.

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

We give an overview of the Grism Lens Amplified Survey from Space (GLASS), a large Hubble Space Telescope program aimed at obtaining grism spectroscopy of the fields of ten massive clusters of galaxies at redshift z=0.308-0.686, including the Hubble Frontier Fields (HFF). The Wide Field Camera 3 yields near infrared spectra of the cluster cores, covering the wavelength range 0.81-1.69mum through grisms G102 and G141, while the Advanced Camera for Surveys in parallel mode provides G800L spectra of the infall regions of the clusters. The WFC3 spectra are taken at two almost orthogonal position angles in order to minimize the effects of confusion. After summarizing the scientific drivers of GLASS, we describe the sample selection as well as the observing strategy and data processing pipeline. We then utilize MACSJ0717.5+3745, a HFF cluster and the first one observed by GLASS, to illustrate the data quality and the high-level data products. Each spectrum brighter than H_AB=23 is visually inspected by at least two co-authors and a redshift is measured when sufficient information is present in the spectra. Furthermore, we conducted a thorough search for emission lines through all the GLASS WFC3 spectra with the aim of measuring redshifts for sources with continuum fainter than H_AB=23. We provide a catalog of 139 emission-line based spectroscopic redshifts for extragalactic sources, including three new redshifts of multiple image systems (one probable, two tentative). In addition to the data itself we also release software tools that are helpful to navigate the data.

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.

We cross-match galaxy cluster candidates selected via their Sunyaev-Zel'dovich effect (SZE) signatures in 129.1 deg$^2$ of the South Pole Telescope 2500d SPT-SZ survey with optically identified clusters selected from the Dark Energy Survey (DES) science verification data. We identify 25 clusters between $0.1\lesssim z\lesssim 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 $\lambda$-mass relation with the following function $\langle\ln\lambda|M_{500}\rangle\propto B_\lambda\ln M_{500}+C_\lambda\ln E(z)$ and use SPT-SZ cluster masses and RM richnesses $\lambda$ to constrain the parameters. We find $B_\lambda= 1.14^{+0.21}_{-0.18}$ and $C_\lambda=0.73^{+0.77}_{-0.75}$. The associated scatter in mass at fixed richness is $\sigma_{\ln M|\lambda} = 0.18^{+0.08}_{-0.05}$ at a characteristic richness $\lambda=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 sub-dominant population characterized by larger offsets. We also cross-match the RM catalog with SPT-SZ candidates below the official catalog threshold significance $\xi=4.5$, using the RM catalog to provide optical confirmation and redshifts for additional low-$\xi$ SPT-SZ candidates.In this way, we identify 15 additional clusters with $\xi\in [4,4.5]$ over the redshift regime explored by RM in the overlapping region between DES science verification data and the SPT-SZ survey.

Finding a sample of the most massive clusters with redshifts $z>0.6$ can provide an interesting consistency check of the $\Lambda$ cold dark matter ($\Lambda$CDM) model. Here we present results from our search for clusters with $0.6\lesssim z\lesssim1.0$ where the initial candidates were selected by cross-correlating the RASS faint and bright source catalogues with red galaxies from the Sloan Digital Sky Survey DR8. Our survey thus covers $\approx10,000\,\rm{deg^2}$, much larger than previous studies of this kind. Deeper follow-up observations in three bands using the William Herschel Telescope and the Large Binocular Telescope were performed to confirm the candidates, resulting in a sample of 44 clusters for which we present richnesses and red sequence redshifts, as well as spectroscopic redshifts for a subset. At least two of the clusters in our sample are comparable in richness to RCS2-$J$232727.7$-$020437, one of the richest systems discovered to date. We also obtained new observations with the Combined Array for Research in Millimeter Astronomy for a subsample of 21 clusters. For 11 of those we detect the Sunyaev-Zel'dovich effect signature. The Sunyaev-Zel'dovich signal allows us to estimate $M_{200}$ and check for tension with the cosmological standard model. We find no tension between our cluster masses and the $\Lambda$CDM model.

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

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and ICL, its high redshift ($z=0.6986$) makes it ideal for studying background galaxies. Using new ACS and WFC3/IR HST data, we identify 16 multiple images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple image, likely confirming the redshift of that system to $z=2.083$. With a highly magnified ($\mu\gtrsim2$) source plane area of $\sim0.7$ arcmin$^2$ at $z=7$, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified $z\sim7$ candidate Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for magnification, we find that all four galaxies are fainter than $0.5 L_{\star}$. One candidate is detected at ${>10\sigma}$ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyman-alpha emission line from any of the four candidates. From the MOSFIRE spectra we place median upper limits on the Lyman-alpha flux of $5-14 \times 10^{-19}\, \mathrm{erg \,\, s^{-1} cm^{-2}}$ ($5\sigma$).

In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging, we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The cluster's gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.

We analyse the spatially-resolved stellar populations of 9 local ($z<0.1$) Brightest Cluster Galaxies (BCGs) observed with VIMOS in IFU mode. Our sample is composed of 7 slow-rotating and 2 fast-rotating BCGs. We do not find a connection between stellar kinematics and stellar populations in this small sample. The BCGs have shallow metallicity gradients (median $\Delta$[Fe/H] $= -0.11\pm0.1$), high central metallicities (median $[$Fe/H]$_{[\alpha/Fe]=0} = 0.13\pm0.07$), and a wide range of central ages (from 5 to 15 Gyr). We propose that the reason for this is diverse evolutionary paths in BCGs. 67 per cent of the sample (6/9) show $\sim 7$ Gyr old central ages, which reflects an active accretion history, and 33 per cent of the sample (3/9) have central ages older than 11 Gyr, which suggest no star formation since $z=2$. The BCGs show similar central stellar populations and stellar population gradients to early-type galaxies of similar mass (M$_{dyn}> 10^{11.3}$M$_{\odot}$) from the ATLAS$^{3D}$ survey (median [Z/H] $= 0.04\pm0.07$, $\Delta$[Z/H] $= -0.19\pm0.1$). However, massive early-type galaxies from ATLAS$^{3D}$ have consistently old ages (median Age $=12.0\pm3.8$Gyr). We also analyse the close massive companion galaxies of two of the BCGs. These galaxies have similar stellar populations to their respective BCGs.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

This is the first in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here we present a new, automated method for identifying relaxed clusters based on their morphologies in X-ray imaging data. While broadly similar to others in the literature, the morphological quantities that we measure are specifically designed to provide a fair basis for comparison across a range of data quality and cluster redshifts, to be robust against missing data due to point-source masks and gaps between detectors, and to avoid strong assumptions about the cosmological background and cluster masses. Based on three morphological indicators - Symmetry, Peakiness and Alignment - we develop the SPA criterion for relaxation. This analysis was applied to a large sample of cluster observations from the Chandra and ROSAT archives. Of the 361 clusters which received the SPA treatment, 57 (16 per cent) were subsequently found to be relaxed according to our criterion. We compare our measurements to similar estimators in the literature, as well as projected ellipticity and other image measures, and comment on trends in the relaxed cluster fraction with redshift, temperature, and survey selection method. Code implementing our morphological analysis will be made available on the web.

We employ robust weak gravitational lensing measurements to improve cosmological constraints from measurements of the galaxy cluster mass function and its evolution, using X-ray selected clusters detected in the ROSAT All-Sky Survey. Our lensing analysis constrains the absolute mass scale of such clusters at the 8 per cent level, including both statistical and systematic uncertainties. Combining it with the survey data and X-ray follow-up observations, we find a tight constraint on a combination of the mean matter density and late-time normalization of the matter power spectrum, $\sigma_8(\Omega_m/0.3)^{0.17}=0.81\pm0.03$, with marginalized, one-dimensional constraints of $\Omega_m=0.26\pm0.03$ and $\sigma_8=0.83\pm0.04$. For these two parameters, this represents a factor of two improvement in precision with respect to previous work, primarily due to the reduced systematic uncertainty in the absolute mass calibration provided by the lensing analysis. Our new results are in good agreement with constraints from cosmic microwave background (CMB) data, both WMAP and Planck (plus WMAP polarization), under the assumption of a flat $\Lambda$CDM cosmology with minimal neutrino mass. Consequently, we find no evidence for non-minimal neutrino mass from the combination of cluster data with CMB, supernova and baryon acoustic oscillation measurements, regardless of which all-sky CMB data set is used (and independent of the recent claimed detection of B-modes on degree scales). We also present improved constraints on models of dark energy (both constant and evolving), modifications of gravity, and primordial non-Gaussianity. Assuming flatness, the constraints for a constant dark energy equation of state from the cluster data alone are at the 15 per cent level, improving to $\sim 6$ per cent when the cluster data are combined with other leading probes.

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

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

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

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

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

We present Spitzer/IRAC observations of nine $z'$-band dropouts highly magnified (2<mu<12) by the Bullet Cluster. We combine archival imaging with our Exploratory program (SURFS UP), which results in a total integration time of ~30 hr per IRAC band. We detect (>3sigma) in both IRAC bands the brightest of these high-redshift galaxies, with [3.6]=23.80+-0.28 mag, [4.5]=23.78+-0.25 mag, and (H-[3.6])=1.17+-0.32 mag. The remaining eight galaxies are undetected to [3.6]~26.4 mag and [4.5]~26.0 mag with stellar masses of ~5x10^7 M_sol. The detected galaxy has an estimated magnification of mu=12+-4, which implies this galaxy has an ultraviolet luminosity of L_1500~0.3 L*_{z=7} --- the lowest luminosity individual source detected in IRAC at z>7. By modeling the broadband photometry, we estimate the galaxy has an intrinsic star-formation rate of SFR~1.3 M_sol/yr and stellar mass of M~2x10^9 M_sol, which gives a specific star-formation rate of sSFR~0.7 Gyr^-1. If this galaxy had sustained this star-formation rate since z~20, it could have formed the observed stellar mass (to within a factor of ~2), we also discuss alternate star-formation histories and argue the exponentially-increasing model is unlikely. Finally, based on the intrinsic star-formation rate, we estimate this galaxy has a likely [C II] flux of <f_[C II]> = 10^{-17} erg/s/cm2.

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)

This paper is the first in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our primary aim is to test the performance of these cluster mass estimation techniques on a diverse set of models that will increase in complexity. We begin by providing participating methods with data from a simple model that delivers idealised clusters, enabling us to quantify the underlying scatter intrinsic to these mass estimation techniques. The mock catalogue is based on a Halo Occupation Distribution (HOD) model that assumes spherical Navarro, Frenk and White (NFW) haloes truncated at R_200, with no substructure nor colour segregation, and with isotropic, isothermal Maxwellian velocities. We find that, above 10^14 M_solar, recovered cluster masses are correlated with the true underlying cluster mass with an intrinsic scatter of typically a factor of two. Below 10^14 M_solar, the scatter rises as the number of member galaxies drops and rapidly approaches an order of magnitude. We find that richness-based methods deliver the lowest scatter, but it is not clear whether such accuracy may simply be the result of using an over-simplistic model to populate the galaxies in their haloes. Even when given the true cluster membership, large scatter is observed for the majority non-richness-based approaches, suggesting that mass reconstruction with a low number of dynamical tracers is inherently problematic.

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey.

The Grism Lens-Amplified Survey from Space (GLASS) is a Hubble Space Telescope (HST) Large Program, which will obtain 140 orbits of grism spectroscopy of the core and infall regions of 10 galaxy clusters, selected to be among the very best cosmic telescopes. Extensive HST imaging is available from many sources including the CLASH and Frontier Fields programs. We introduce the survey by analyzing spectra of faint multiply-imaged galaxies and $z\gtrsim6$ galaxy candidates obtained from the first seven orbits out of fourteen targeting the core of the Frontier Fields cluster MACS0717.5+3745. Using the G102 and G141 grisms to cover the wavelength range 0.8-1.7$\mu$m, we confirm 4 strongly lensed systems by detecting emission lines in each of the images. For the 9 $z\gtrsim6$ galaxy candidates clear from contamination, we do not detect any emission lines down to a seven-orbit 1$\sigma$ noise level of $\sim$5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. Taking lensing magnification into account, our flux sensitivity reaches $\sim$0.2-5$\times$10$^{-18}$erg s$^{-1}$cm$^{-2}$. These limits over an uninterrupted wavelength range rule out the possibility that the high-$z$ galaxy candidates are instead strong line emitters at lower redshift. These results show that by means of careful modeling of the background - and with the assistance of lensing magnification - interesting flux limits can be reached for large numbers of objects, avoiding pre-selection and the wavelength restrictions inherent to ground-based multi-slit spectroscopy. These observations confirm the power of slitless HST spectroscopy even in fields as crowded as a cluster core.

We present a multi-wavelength study of the interstellar medium in eight nearby, X-ray and optically bright, giant elliptical galaxies. Using Herschel PACS, we map the cold gas in the lines of [CII], [OI], and [OIb]. Additionally, we present Ha+[NII] imaging of warm ionized gas with the SOAR telescope, and a study of the hot X-ray emitting plasma with Chandra. All systems with extended Ha emission in our sample (6/8 galaxies) display significant [CII] line emission indicating the presence of cold gas. This emission is co-spatial with the Ha+[NII] emitting nebulae and the lowest entropy X-ray emitting plasma. The entropy profiles of the hot galactic atmospheres show a clear dichotomy, with the systems displaying extended emission line nebulae having lower entropies beyond r~1 kpc than the cold-gas-poor systems. We show that while the hot atmospheres of the cold-gas-poor galaxies are thermally stable outside of their innermost cores, the atmospheres of the cold-gas-rich systems are prone to cooling instabilities. This result indicates that the cold gas is produced chiefly by thermally unstable cooling from the hot phase. We show that cooling instabilities may develop more easily in rotating systems and discuss an alternative condition for thermal instability for this case. The hot atmospheres of cold-gas-rich galaxies display disturbed morphologies indicating that the accretion of clumpy multiphase gas in these systems may result in variable power output of the AGN jets, potentially triggering sporadic, larger outbursts. In the two cold-gas-poor, X-ray morphologically relaxed galaxies of our sample, NGC 1399 and NGC 4472, powerful AGN outbursts may have destroyed or removed most of the cold gas from the cores, allowing the jets to propagate and deposit most of their energy further out, increasing the entropy of the hot galactic atmospheres and leaving their cores relatively undisturbed.

Using the VIMOS Integral Field Unit (IFU) spectrograph on the Very Large Telescope (VLT), we have spatially mapped the kinematic properties of 10 nearby Brightest Cluster Galaxies (BCGs) and 4 BCG companion galaxies located within a redshift of $z=0.1$. In the hierarchical formation model, these massive galaxies $(10^{10.5} M_{\odot} < M_{dyn} < 10^{11.9} M_{\odot})$ are expected to undergo more mergers than lower mass galaxies, and simulations show that dry minor mergers can remove angular momentum. We test whether BCGs have low angular momenta by using the $\lambda_{Re}$ parameter developed by the SAURON and ATLAS\textsuperscript{3D} teams and combine our kinematics with Sloan Digital Sky Survey (SDSS) photometry to analyze the BCGs' merger status. We find that 30% (3/10) of the BCGs and 100% of the companion galaxies (4/4) are fast rotators as defined by the ATLAS\textsuperscript{3D} criteria. Our fastest rotating BCG has a $\lambda_{Re}=0.35\pm0.05$. We increase the number of BCGs analyzed from 1 in the combined SAURON and ATLAS\textsuperscript{3D} surveys to 11 BCGs total and find that above $M_{dyn}\sim11.5 M_{\odot}$, virtually all galaxies regardless of environment are slow rotators. To search for signs of recent merging, we analyze the photometry of each system and use the $G-M_{20}$ selection criteria to identify mergers. We find that $40\pm20$% of our BCGs are currently undergoing or have recently undergone a merger (within 0.2 Gyrs). Surprisingly, we find no correlation between galaxies with high angular momentum and morphological signatures of merging.

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

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We present a detailed Chandra, XMM-Newton, VLA and HST analysis of one of the strongest cool core clusters known, RX J1532.9+3021 (z=0.3613). Using new, deep 90 ks Chandra observations, we confirm the presence of a western X-ray cavity or bubble, and report on a newly discovered eastern X-ray cavity. The total mechanical power associated with these AGN-driven outflows is (22+/-9)*10^44 erg/s, and is sufficient to offset the cooling, indicating that AGN feedback still provides a viable solution to the cooling flow problem even in the strongest cool core clusters. Based on the distribution of the optical filaments, as well as a jet-like structure seen in the 325 MHz VLA radio map, we suggest that the cluster harbours older outflows along the north to south direction. The jet of the central AGN is therefore either precessing or sloshing-induced motions have caused the outflows to change directions. There are also hints of an X-ray depression to the north aligned with the 325 MHz jet-like structure, which might represent the highest redshift ghost cavity discovered to date. We further find evidence of a cold front (r=65kpc) that coincides with the outermost edge of the western X-ray cavity and the edge of the radio mini-halo. The common location of the cold front with the edge of the radio mini-halo supports the idea that the latter originates from electrons being reaccelerated due to sloshing induced turbulence. Alternatively, its coexistence with the edge of the X-ray cavity may be due to cool gas being dragged out by the outburst. We confirm that the central AGN is highly sub-Eddington and conclude that a >10^10M_Sun or a rapidly spinning black hole is favoured to explain both the radiative-inefficiency of the AGN and the powerful X-ray cavities.

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

Major astrophysical questions related to the formation and evolution of structures, and more specifically of galaxy groups and clusters, will still be open in the coming decade and beyond: what is the interplay of galaxy, supermassive black hole, and intergalactic gas evolution in the most massive objects in the Universe - galaxy groups and clusters? What are the processes driving the evolution of chemical enrichment of the hot diffuse gas in large-scale structures? How and when did the first galaxy groups in the Universe, massive enough to bind more than 10^7 K gas, form? Focussing on the period when groups and clusters assembled (0.5<z<2.5), we show that, due to the continuum and line emission of this hot intergalactic gas at X-ray wavelengths, Athena+, combining high sensitivity with excellent spectral and spatial resolution, will deliver breakthrough observations in view of the aforementioned issues. Indeed, the physical and chemical properties of the hot intra-cluster gas, and their evolution across time, are a key to understand the co-evolution of galaxy and supermassive black hole within their environments.

We present a multi-wavelength study of the emission-line nebulae located southeast of the nucleus of M87, the central dominant galaxy of the Virgo Cluster. We report the detection of far-infrared (FIR) [CII] line emission from the nebulae using observations made with Herschel PACS. The infrared line emission is extended and cospatial with optical H{\alpha}+[NII], far-ultraviolet CIV lines, and soft X-ray emission. The filamentary nebulae evidently contain multi-phase material spanning a temperature range of at least 5 orders of magnitude, from ~100 K to ~10^7 K. This material has most likely been uplifted by the AGN from the center of M87. The thermal pressure of the 10^4 K phase appears to be significantly lower than that of the surrounding hot intra-cluster medium (ICM) indicating the presence of additional turbulent and magnetic pressure in the filaments. If the turbulence in the filaments is subsonic then the magnetic field strength required to balance the pressure of the surrounding ICM is B~30-70 {\mu}G. The spectral properties of the soft X-ray emission from the filaments indicate that it is due to thermal plasma with kT~0.5-1 keV, which is cooling by mixing with the cold gas and/or radiatively. Charge exchange can be ruled out as a significant source of soft X-rays. Both cooling and mixing scenarios predict gas with a range of temperatures. This is at first glance inconsistent with the apparent lack of X-ray emitting gas with kT<0.5 keV. However, we show that the missing very soft X-ray emission could be absorbed by the cold gas in the filaments with an integrated hydrogen column density of ~1.6x10^21 cm^-2, providing a natural explanation for the apparent temperature floor to the X-ray emission at kT~0.5 keV. The FIR through ultra-violet line emission is most likely primarily powered by the ICM particles penetrating the cold gas following a shearing induced mixing process.

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

We present new spectroscopic observations in a field containing the highest redshift cluster of the ESO Distant Cluster Survey (EDisCS). We measure galaxy redshifts and determine the velocity dispersions of the galaxy structures located in this field. Together with the main cluster Cl1103.7$-$1245 (z=0.9580; sigma_{clus} = 522 +/- 111 km/s) we find a secondary structure at z=0.9830, Cl1103.7-1245c. We then characterize the galaxy properties in both systems, and find that they contain very different galaxy populations. The cluster Cl1103.7-1245 hosts a mixture of passive elliptical galaxies and star-forming spirals and irregulars. In the secondary structure Cl1103.7-1245c all galaxies are lower-mass star-forming irregulars and peculiars. In addition, we compare the galaxy populations in the Cl1103.7-1245 z=0.9580 cluster with those in lower redshift EDisCS clusters with similar velocity dispersions. We find that the properties of the galaxies in Cl1103.7-1245 follow the evolutionary trends found at lower redshifts: the number of cluster members increases with time in line with the expected growth in cluster mass, and the fraction of passive early-type galaxies increases with time while star-forming late types become less dominant. Finally, we find that the mean stellar masses are similar in all clusters, suggesting that massive cluster galaxies were already present at z~1.

We present Integral Field Unit (IFU) spectroscopy of four brightest cluster galaxies (BCGs) at z~0.1. Three of the BCGs have close companions within a projected radius of 20 kpc and one has no companion within that radius. We calculate the dynamical masses of the BCGs and their companions to be 1.4x10^11<M_dyn (M_solar)<1.5x10^12. We estimate the probability that the companions of the BCGs are bound using the observed masses and velocity offsets. We show that the lowest mass companion (1:4) is not bound while the two nearly equal mass (1:1.45 and 1:1.25) companions are likely to merge with their host BCGs in 0.35 Gyr in major, dry mergers. We conclude that some BCGs continue to grow from major merging even at z~0. We analyse the stellar kinematics of these systems using the \lambda_R parameter developed by the SAURON team. This offers a new and unique means to measure the stellar angular momentum of BCGs and make a direct comparison to other early-type galaxies. The BCGs and their companions have similar ellipticities to those of other early-type galaxies but are more massive. We find that not all these massive galaxies have low \lambda_R_e as one might expect. One of the four BCGs and the two massive companions are found to be fast-rotating galaxies with high angular momentum, thereby providing a new test for models of galaxy evolution and the formation of Intra-Cluster Light.

We report on a deep, multiwavelength study of the galaxy cluster MACS J1931.8-2634 using Chandra X-ray, Subaru optical, and VLA 1.4 GHz radio data. This cluster (z=0.352) harbors one of the most X-ray luminous cool cores yet discovered, with an equivalent mass cooling rate within the central 50 kpc is approximately 700 solar masses/yr. Unique features observed in the central core of MACSJ1931.8-2634 hint to a wealth of past activity that has greatly disrupted the original cool core. We observe a spiral of relatively cool, dense, X-ray emitting gas connected to the cool core, as well as highly elongated intracluster light (ICL) surrounding the cD galaxy. Extended radio emission is observed surrounding the central AGN, elongated in the east-west direction, spatially coincident with X-ray cavities. The power input required to inflate these `bubbles' is estimated from both the X-ray and radio emission to reside between 4 and 14e45 erg/s, putting it among the most powerful jets ever observed. This combination of a powerful AGN outburst and bulk motion of the cool core have resulted in two X-ray bright ridges to form to the north and south of the central AGN at a distance of approximately 25 kpc. The northern ridge has spectral characteristics typical of cool cores and is consistent with being a remnant of the cool core after it was disrupted by the AGN and bulk motions. It is also the site of H-alpha filaments and young stars. The X-ray spectroscopic cooling rate associated with this ridge is approximately 165 solar masses/yr, which agrees with the estimate of the star formation rate from broad-band optical imaging (170 solar masses/yr). MACS J1931.8-2634 appears to harbor one of most profoundly disrupted low entropy cores observed in a cluster, and offers new insights into the survivability of cool cores in the context of hierarchical structure formation.

We investigate how the specific star formation rates of galaxies of different masses depend on cluster-centric radius and on the central/satellite dichotomy in both field and cluster environments. Recent data from a variety of sources, including the cluster catalogue of von der Linden et al. are compared to the semi-analytic models of De Lucia & Blaizot. We find that these models predict too many passive satellite galaxies in clusters, too few passive central galaxies with low stellar masses, and too many passive central galaxies with high masses. We then outline a series of modifications to the model necessary to solve these problems: a) Instead of instantaneous stripping of the external gas reservoir after a galaxy becomes a satellite, the gas supply is assumed to decrease at the same rate that the surrounding halo loses mass due to tidal stripping, b) The AGN feedback efficiency is lowered to bring the fraction of massive passive centrals in better agreement with the data. We also allow for radio mode AGN feedback in satellite galaxies. c) We assume that satellite galaxies residing in host haloes with masses below 10^12 M_sun do not undergo any stripping. We highlight the fact that in low mass galaxies, the external reservoir is composed primarily of gas that has been expelled from the galactic disk by supernovae driven winds. This gas must remain available as a future reservoir for star formation, even in satellite galaxies. Finally, we present a simple recipe for the stripping of gas and dark matter in satellites that can be used in models where subhalo evolution is not followed in detail.

Using a combination of deep 574ks Chandra data, XMM-Newton high-resolution spectra, and optical Halpha+NII images, we study the nature and spatial distribution of the multiphase plasma in M87. Our results provide direct observational evidence of `radio mode' AGN feedback in action, stripping the central galaxy of its lowest entropy gas and preventing star-formation. This low entropy gas was entrained with and uplifted by the buoyantly rising relativistic plasma, forming long "arms". These arms are likely oriented within 15-30 degrees of our line-of-sight. The mass of the uplifted gas in the arms is comparable to the gas mass in the approximately spherically symmetric 3.8 kpc core, demonstrating that the AGN has a profound effect on its immediate surroundings. The coolest X-ray emitting gas in M87 has a temperature of ~0.5 keV and is spatially coincident with Halpha+NII nebulae, forming a multiphase medium where the cooler gas phases are arranged in magnetized filaments. We place strong upper limits of 0.06 Msun/yr on the amount of plasma cooling radiatively from 0.5 keV and show that a uniform, volume-averaged heating mechanism could not be preventing the cool gas from further cooling. All of the bright Halpha filaments appear in the downstream region of the <3 Myr old shock front, at smaller radii than ~0.6'. We suggest that shocks induce shearing around the filaments, thereby promoting mixing of the cold gas with the ambient hot ICM via instabilities. By bringing hot thermal particles into contact with the cool, line-emitting gas, mixing can supply the power and ionizing particles needed to explain the observed optical spectra. Mixing of the coolest X-ray emitting plasma with the cold optical line emitting filamentary gas promotes efficient conduction between the two phases, allowing non-radiative cooling which could explain the lack of X-ray gas with temperatures under 0.5 keV.

We show that, observationally, the projected local density distribution in high-z clusters is shifted towards higher values compared to clusters at lower redshift. To search for the origin of this evolution, we analyze a sample of haloes selected from the Millennium Simulation and populated using semi-analytic models, investigating the relation between observed projected density and physical 3D density, using densities computed from the 10 and 3 closest neighbours. Both observationally and in the simulations, we study the relation between number of cluster members and cluster mass, and number of members per unit of cluster mass. We find that the observed evolution of projected densities reflects a shift to higher values of the physical 3D density distribution. In turn, this must be related with the globally higher number of galaxies per unit of cluster volume N/V in the past. We show that the evolution of N/V is due to a combination of two effects: a) distant clusters were denser in dark matter (DM) simply because the DM density within R_{200} (~the cluster virial radius) is defined to be a fixed multiple of the critical density of the Universe, and b) the number of galaxies per unit of cluster DM mass is remarkably constant both with redshift and cluster mass if counting galaxies brighter than a passively evolving magnitude limit. Our results highlight that distant clusters were much denser environments than today's clusters, both in galaxy number and mass, and that the density conditions felt by galaxies in virialized systems do not depend on the system mass.

This paper reports the results obtained on the photometric redshifts measurement and accuracy, and cluster tomography in the ESO Distant Cluster Survey (EDisCS) fields. Photometric redshifts were computed using two independent codes (Hyperz and G. Rudnick's code). The accuracy of photometric redshifts was assessed by comparing our estimates with the spectroscopic redshifts of ~1400 galaxies in the 0.3<z<1.0 domain. The accuracy for galaxies fainter than the spectroscopic control sample was estimated using a degraded version of the photometric catalog for the spectroscopic sample. The accuracy of photometric redshifts is typically sigma(Delta z/(1+z)) ~ 0.05+/-0.01, depending on the field, the filter set, and the spectral type of the galaxies. The quality of the photometric redshifts degrades by a factor of two in sigma(Delta z/(1+z)) between the brightest (I~22) and the faintest (I~24-24.5) galaxies in the EDisCS sample. The photometric determination of cluster redshifts in the EDisCS fields using a simple algorithm based on photoz is in excellent agreement with the spectroscopic values (delta(z) ~0.03-0.04 in the high-z sample and ~0.05 in the low-z sample). We also developed a method that uses both photz codes jointly to reject interlopers at magnitudes fainter than the spectroscopic limit. When applied to the spectroscopic sample, this method rejects ~50-90% of all spectroscopically confirmed non-members, while retaining at least 90% of all confirmed members. Photometric redshifts are found to be particularly useful for the identification and study of clusters of galaxies in large surveys. They enable efficient and complete pre-selection of cluster members for spectroscopy, allow accurate determinations of the cluster redshifts, and provide a means of determining cluster membership. (Abridged)

We study the morphological content of a large sample of high-redshift clusters to determine its dependence on cluster mass and redshift. Quantitative morphologies are based on bulge+disk decompositions of cluster and field galaxies on deep VLT/FORS2 images of 18 optically-selected clusters at 0.45 < z < 0.80 from the ESO Distant Cluster Survey (EDisCS). Morphological content is given by the early-type galaxy fraction f_et, and early-type galaxies are selected based on their bulge fraction and image smoothness. A set of 158 SDSS clusters is analyzed exactly as the EDisCS sample to provide a robust local comparison. Our main results are: (1) f_et values for the SDSS and EDisCS clusters exhibit no clear trend as a function of sigma. (2) Mid-z EDisCS clusters around sigma = 500 km/s have f_et ~= 0.5 whereas high-z EDisCS clusters have f_et ~= 0.4 (~25% increase over 2 Gyrs). (3) There is a marked difference in the morphological content of EDisCS and SDSS clusters. None of the EDisCS clusters have f_et greater than 0.6 whereas half of the SDSS clusters lie above this value. This difference is seen in clusters of all velocity dispersions. (4) There is a strong correlation between morphology and star formation in SDSS and EDisCS clusters. This correlation holds independent of sigma and z even though the fraction of [OII] emitters decreases from z~0.8 to z~0.06 in all environments. Our results pose an interesting challenge to structural transformation and star formation quenching processes that strongly depend on the global cluster environment and suggest that cluster membership may be of lesser importance than other variables in determining galaxy properties. (ABRIDGED)

We investigate the recent and current star formation activity of galaxies as function of distance from the cluster center in a sample of 521 SDSS clusters at z<0.1. We show that when the BCGs are excluded from the galaxy sample, there is no evidence for mass segregation in the clusters, so that differences in cluster and field populations cannot simply be attributed to different mass functions. We find a marked star formation-radius relation in that almost all galaxies in the cluster core are quiescent, i.e. have terminated star formation a few Gyr ago. This star formation-radius relation is most pronounced for low-mass galaxies and is very weak or absent beyond the virial radius. The typical star formation rate of non-quiescent galaxies declines by approximately a factor of two towards the cluster center. However, the fraction of galaxies with young stellar populations indicating a recently completed starburst or a truncation of star formation does not vary significantly with radius. The fraction of star-forming galaxies which host a powerful optical AGN is also independent of clustercentric radius, indicating that the link between star formation and AGN in these galaxies operates independent of environment. The fraction of red galaxies which host a weak optical AGN decreases, however, towards the cluster center, with a similar timescale as the decline of star forming galaxies. Our results favor a scenario in which star formation is quenched slowly, on timescales similar to the cluster crossing time, and rule out significant contributions from more violent processes at least beyond cluster radii >0.1 R_200.

Many brightest cluster galaxies (BCGs) at the centers of X-ray selected clusters exhibit clear evidence for recent star formation. However, studies of BCGs in optically-selected clusters show that star formation is not enhanced when compared to control samples of non-BCGs of similar stellar mass. Here we analyze a sample of 113 BCGs in low redshift (z<0.1), optically-selected clusters, a matched control sample of non-BCGs, and a smaller sample of BCGs in X-ray selected clusters. We convolve the SDSS images of the BCGs to match the resolution of the GALEX data and we measure UV-optical colours in their inner and outer regions. We find that optically-selected BCGs exhibit smaller scatter in optical colours and redder inner NUV-r colours than the control galaxies, indicating that they are a homogenous population with very little ongoing star formation. The BCGs in the X-ray selected cluster sample span a similar range in optical colours, but have bluer NUV-r colours. Among X-ray selected BCGs, those located in clusters with central cooling times of less than 1 Gyr are significantly bluer than those located in clusters where the central gas cooling times are long. Our main conclusion is that the location of a galaxy at the centre of its halo is not sufficient to determine whether or not it is currently forming stars. One must also have information about the thermodynamic state of the gas in the core of the halo.

We present the rest-frame optical luminosity function (LF) of red sequence galaxies in 16 clusters at 0.4<z<0.8 drawn from the ESO Distant Cluster Survey (EDisCS). We compare our clusters to an analogous sample from the Sloan Digital Sky Survey (SDSS) and match the EDisCS clusters to their most likely descendants. We measure all LFs down to M M* + (2.5 - 3.5). At z<0.8, the bright end of the LF is consistent with passive evolution but there is a significant build-up of the faint end of the red sequence towards lower redshift. There is a weak dependence of the LF on cluster velocity dispersion for EDisCS but no such dependence for the SDSS clusters. We find tentative evidence that red sequence galaxies brighter than a threshold magnitude are already in place, and that this threshold evolves to fainter magnitudes toward lower redshifts. We compare the EDisCS LFs with the LF of co-eval red sequence galaxies in the field and find that the bright end of the LFs agree. However, relative to the number of bright red galaxies, the field has more faint red galaxies than clusters at 0.6<z<0.8 but fewer at 0.4<z<0.6, implying differential evolution. We compare the total light in the EDisCS cluster red sequences to the total red sequence light in our SDSS cluster sample. Clusters at 0.4<z<0.8 must increase their luminosity on the red sequence (and therefore stellar mass in red galaxies) by a factor of 1-3 by z=0. The necessary processes that add mass to the red sequence in clusters predict local clusters that are over-luminous as compared to those observed in the SDSS. The predicted cluster luminosities can be reconciled with observed local cluster luminosities by combining multiple previously known effects.

The ESO Distant Cluster Survey (EDisCS, P.I. Simon D.M. White, LP 166.A-0162) is an ESO large programme aimed at studying clusters and cluster galaxies at z=0.4-1. How different is the evolution of the star formation activity in clusters, in groups and in the field? Does it depend on cluster mass and/or the local galaxy density? How relevant are starburst and post-starburst galaxies in the different environments? Is there an evolution in the galaxies' structures, and if so, is this related to the changes in their star formation activity? These are some of the main questions that have been investigated using the EDisCS dataset.

Post-starburst (E+A or k+a) spectra, characterized by their exceptionally strong Balmer lines in absorption and the lack of emission lines, belong to galaxies in which the star formation activity ended abruptly sometime during the past Gyr. We perform a spectral analysis of galaxies in clusters, groups, poor groups and the field at z=0.4-0.8 based on the ESO Distant Cluster Survey. The incidence of k+a's at these redshifts depends strongly on environment. K+a's reside preferentially in clusters and, unexpectedly, in a subset of the sigma = 200-400 km/s groups, those that have a low fraction of [OII] emitters. In these environments, 20-30% of the recently star-forming galaxies have had their star formation activity recently truncated. In contrast, there are proportionally fewer k+a's in the field, the poor groups and groups with a high [OII] fraction. The incidence of k+a galaxies correlates with the cluster velocity dispersion: more massive clusters have higher proportions of k+a's. Spectra of dusty starburst candidates, with strong Balmer absorption and emission lines, present a very different environmental dependence from k+a's. They are numerous in all environments at z=0.4-0.8, but they are especially numerous in all types of groups, favoring the hypothesis of triggering by a merger. Our observations are consistent with previous suggestions that cluster k+a galaxies are observed in a transition phase as massive S0 and Sa galaxies, evolving from star-forming later types to passively evolving early-type galaxies. The correlation between k+a fraction and cluster sigma supports the hypothesis that k+a galaxies in clusters originate from processes related to the intracluster medium, while several possibilities are discussed for the origin of the k+a frequency in low-[OII] groups.(abr.)

We constrain the physical nature of dark matter using the newly identified massive merging galaxy cluster MACSJ0025.4-1222. As was previously shown by the example of the Bullet Cluster (1E0657-56), such systems are ideal laboratories for detecting isolated dark matter, and distinguishing between cold dark matter (CDM) and other scenarios (e.g. self-interacting dark matter, alternative gravity theories). MACSJ0025.4-1222 consists of two merging subclusters of similar richness at z=0.586. We measure the distribution of X-ray emitting gas from Chandra X-ray data and find it to be clearly displaced from the distribution of galaxies. A strong (information from highly distorted arcs) and weak (using weakly distorted background galaxies) gravitational lensing analysis based on Hubble Space Telescope observations and Keck arc spectroscopy confirms that the subclusters have near-equal mass. The total mass distribution in each of the subclusters is clearly offset (at >4sigma significance) from the peak of the hot X-ray emitting gas (the main baryonic component), but aligned with the distribution of galaxies. We measure the fractions of mass in hot gas (0.09^{+0.07}_{-0.03}) and stars (0.010^{+0.007}_{-0.004}), consistent with those of typical clusters, finding that dark matter is the dominant contributor to the gravitational field. Under the assumption that the subclusters experienced a head-on collision in the plane of the sky, we obtain an order-of-magnitude estimate of the dark matter self-interaction cross-section of sigma/m < 4cm^2/g, re-affirming the results from the Bullet Cluster on the collisionless nature of dark matter.

We use the ROSAT All-Sky Survey to study the X-ray properties of a sample of 625 groups and clusters of galaxies selected from the Sloan Digital Sky Survey. We stack clusters with similar velocity dispersions and investigate whether their average X-ray luminosities and surface brightness profiles vary with the radio activity level of their central galaxies. We find that at a given value of $\sigma$, clusters with a central radio AGN have more concentrated X-ray surface brightness profiles, larger central galaxy masses, and higher X-ray luminosities than clusters with radio-quiet central galaxies. The enhancement in X-ray luminosity is more than a factor of two, is detected with better than 6$\sigma$ significance, and cannot be explained by X-ray emission from the radio AGN itself. This difference is largely due to a subpopulation of radio-quiet, high velocity dispersion clusters with low mass central galaxies. These clusters are underluminous at X-ray wavelengths when compared to otherwise similar clusters where the central galaxy is radio-loud, more massive, or both.

We investigate how the [OII] properties and the morphologies of galaxies in clusters and groups at z=0.4-0.8 depend on projected local galaxy density, and compare with the field at similar redshifts and clusters at low-z. In both nearby and distant clusters, higher-density regions contain proportionally fewer star-forming galaxies, and the average [OII] equivalent width of star-forming galaxies is independent of local density. However, in distant clusters the average current star formation rate (SFR) in star-forming galaxies seems to peak at densities ~15-40 galaxies Mpc^{-2}. At odds with low-z results, at high-z the relation between star-forming fraction and local density varies from high- to low-mass clusters. Overall, our results suggest that at high-z the current star formation (SF) activity in star-forming galaxies does not depend strongly on global or local environment, though the possible SFR peak seems at odds with this conclusion. We find that the cluster SFR normalized by cluster mass anticorrelates with mass and correlates with the star-forming fraction. These trends can be understood given a) that the average star-forming galaxy forms about 1 Msun/yr in all clusters; b) that the total number of galaxies scales with cluster mass and c) the dependence of star-forming fraction on cluster mass. We present the morphology-density (MD) relation for our z=0.4-0.8 clusters, and uncover that the decline of the spiral fraction with density is entirely driven by galaxies of types Sc or later. For galaxies of a given Hubble type, we see no evidence that SF properties depend on local environment. In contrast with recent findings at low-z, in our distant clusters the SF-density relation and the MD-relation are equivalent, suggesting that neither of the two is more fundamental than the other.(abr.)

[Abridged] We present K-band data for the brightest cluster galaxies (BCGs) from the ESO Distant Cluster Survey. These data are combined with photometry from Aragon-Salamanca et al. (1998) and a low-redshift comparison sample from von der Linden et al. (2007). The K-band Hubble diagram for BCGs exhibits very low scatter (~0.35mag) since z=1. The colour and $K$-band luminosity evolution of the BCGs are in good agreement with passively-evolving stellar populations formed at z>2. We do not detect any significant change in the stellar mass of the BCG since z~1. These results do not seem to depend on the velocity dispersion of the parent cluster. There is a correlation between the 1D velocity dispersion of the clusters and the K-band luminosity of the BCGs (after correcting for passive evolution). The clusters with large velocity dispersions tend to have brighter BCGs, i.e., BCGs with larger stellar masses. This dependency, although significant, is relatively weak: the stellar mass of the BCGs changes only by ~70% over a two-order-of-magnitude range in cluster mass. This dependency doesn't change significantly with redshift. The models of De Lucia & Blaizot (2007) predict colours which are in reasonable agreement with the observations because the growth in stellar mass is dominated by the accretion of old stars. However, the stellar mass in the model BCGs grows by a factor of 3-4 since z=1, a growth rate which seems to be ruled out by the observations. The models predict a dependency between the BCG's stellar mass and the velocity dispersion of the parent cluster in the same sense as the data, but the dependency is significantly stronger than observed. However, one major difficulty in this comparison is that we have measured fixed metric aperture magnitudes while the models compute total luminosities.

AIMS. We present spectroscopic observations of galaxies in 15 survey fields as part of the ESO Distant Cluster Survey (EDisCS). We determine the redshifts and velocity dispersions of the galaxy clusters located in these fields, and we test for possible substructure in the clusters. METHODS. We obtained multi-object mask spectroscopy using the FORS2 instrument at the VLT. We reduced the data with particular attention to the sky subtraction. We implemented the method of Kelson for performing sky subtraction prior to any rebinning/interpolation of the data. From the measured galaxy redshifts, we determine cluster velocity dispersions using the biweight estimator and test for possible substructure in the clusters using the Dressler-Shectman test. RESULTS. The method of subtracting the sky prior to any rebinning/interpolation of the data delivers photon-noise-limited results, whereas the traditional method of subtracting the sky after the data have been rebinned/interpolated results in substantially larger noise for spectra from tilted slits. Redshifts for individual galaxies are presented and redshifts and velocity dispersions are presented for 21 galaxy clusters. For the 9 clusters with at least 20 spectroscopically confirmed members, we present the statistical significance of the presence of substructure obtained from the Dressler-Shectman test, and substructure is detected in two of the clusters. CONCLUSIONS. Together with data from our previous paper, spectroscopy and spectroscopic velocity dispersions are now available for 26 EDisCS clusters with redshifts in the range 0.40-0.96 and velocity dispersions in the range 166-1080 km/s.

(Abridged) The prevalence of radio-loud AGN activity in present-day massive halos is determined using a sample of 625 nearby groups and clusters from the SDSS. Brightest group and cluster galaxies (BCGs) are more likely to host a radio-loud AGN than other galaxies of the same stellar mass (by less than a factor of two at a stellar mass of 5e11 M_sun, but over an order of magnitude below 1e11 M_sun). The distribution of radio luminosities for BCGs does not depend on stellar mass, however, and is similar to that of field galaxies of the same mass. Neither the radio-loud fraction nor the radio luminosity distribution of BCGs depends strongly on the cluster velocity dispersion. The radio-AGN fraction is also studied as a function of clustercentric distance: only within 0.2 r_200 do cluster galaxies exhibit an enhanced likelihood of radio-loud AGN activity, approaching that of the BCGs. The radio-loud AGN properties of both BCGs and non-BCGs can naturally be explained if this activity is fuelled by cooling from hot gas surrounding the galaxy. Using observational estimates of the mechanical output of the radio jets, the time-averaged energy output associated with recurrent radio source activity is estimated for all group/cluster galaxies. Within the cooling radius of the cluster, the radio-mode heating from the BCG dominates over that of all other galaxies combined. The scaling between total radio-AGN energy output and cluster velocity dispersion is considerably shallower than the sigma^4 scaling of the radiative cooling rate. Thus, unless either the mechanical-to-radio luminosity ratio or the efficiency of converting AGN mechanical energy into heating increases by 2-3 orders of magnitude between groups and rich clusters, radio-mode heating will not balance cooling in systems of all masses.

We use the Sloan Digital Sky Survey to construct a sample of 625 brightest group and cluster galaxies (BCGs) together with control samples of non-BCGs matched in stellar mass, redshift, and color. We investigate how the systematic properties of BCGs depend on stellar mass and on their privileged location near the cluster center. The groups and clusters that we study are drawn from the C4 catalogue of Miller et al. (2005) but we have developed improved algorithms for identifying the BCG and for measuring the cluster velocity dispersion. Since the SDSS photometric pipeline tends to underestimate the luminosities of large galaxies in dense environments, we have developed a correction for this effect which can be readily applied to the published catalog data. We find that BCGs are larger and have higher velocity dispersions than non-BCGs of the same stellar mass, which implies that BCGs contain a larger fraction of dark matter. In contrast to non-BCGs, the dynamical mass-to-light ratio of BCGs does not vary as a function of galaxy luminosity. Hence BCGs lie on a different fundamental plane than ordinary elliptical galaxies. BCGs also follow a steeper Faber-Jackson relation than non-BCGs, as suggested by models in which BCGs assemble via dissipationless mergers along preferentially radial orbits. We find tentative evidence that this steepening is stronger in more massive clusters. BCGs have similar mean stellar ages and metallicities to non-BCGs of the same mass, but they have somewhat higher alpha/Fe ratios, indicating that star formation may have occurred over a shorter timescale in the BCGs. Finally, we find that BCGs are more likely to host radio-loud active galactic nuclei than other galaxies of the same mass, but are less likely to host an optical AGN. The differences we find are more pronounced for the less massive BCGs.

Using galaxy clusters from the ESO Distant Cluster Survey, we study how the distribution of galaxies along the colour-magnitude relation has evolved since z~0.8. While red-sequence galaxies in all these clusters are well described by an old, passively evolving population, we confirm our previous finding of a significant evolution in their luminosity distribution as a function of redshift. When compared to galaxy clusters in the local Universe, the high redshift EDisCS clusters exhibit a significant "deficit" of faint red galaxies. Combining clusters in three different redshift bins, and defining as `faint' all galaxies in the range 0.4 > L/L* > 0.1, we find a clear decrease in the luminous-to-faint ratio of red galaxies from z~0.8 to z~0.4. The amount of such a decrease appears to be in qualitative agreement with predictions of a model where the blue bright galaxies that populate the colour-magnitude diagram of high redshift clusters, have their star formation suppressed by the hostile cluster environment. Although model results need to be interpreted with caution, our findings clearly indicate that the red-sequence population of high-redshift clusters does not contain all progenitors of nearby red-sequence cluster galaxies. A significant fraction of these must have moved onto the red-sequence below z~0.8.

We study how the proportion of star-forming galaxies evolves between z=0.8 and z=0 as a function of galaxy environment, using the [OII] line in emission as a signature of ongoing star formation. Our high-z dataset comprises 16 clusters, 10 groups and another 250 galaxies in poorer groups and the field at z=0.4-0.8 from the ESO Distant Cluster Survey, plus another 9 massive clusters at similar redshifts. As a local comparison, we use samples of galaxy systems selected from the Sloan Digital Sky Survey at 0.04< z < 0.08. At high-z most systems follow a broad anticorrelation between the fraction of star-forming galaxies and the system velocity dispersion. At face value, this suggests that at z=0.4-0.8 the mass of the system largely determines the proportion of galaxies with ongoing star formation. At these redshifts the strength of star formation (as measured by the [OII] equivalent width) in star-forming galaxies is also found to vary systematically with environment. Sloan clusters have much lower fractions of star-forming galaxies than clusters at z=0.4-0.8 and, in contrast with the distant clusters, show a plateau for velocity dispersions $ \ge 550 km s^-1$, where the fraction of galaxies with [OII] emission does not vary systematically with velocity dispersion. We quantify the evolution of the proportion of star-forming galaxies as a function of the system velocity dispersion and find it is strongest in intermediate-mass systems (sigma ~ 500-600 km s^-1 at z=0). To understand the origin of the observed trends, we use the Press-Schechter formalism and the Millennium Simulation and show that galaxy star formation histories may be closely related to the growth history of clusters and groups. We propose a scheme that is able to account for the observed relations between the star-forming fraction and \sigma [abridged].

Weak lensing surveys provide the possibility of identifying dark matter halos based on their total matter content rather than just the luminous matter content. On the basis of two sets of observations carried out with the CFHT, Erben et al. (2000) presented the first candidate dark clump, i.e. a dark matter concentration identified by its significant weak lensing signal without a corresponding galaxy overdensity or X-ray emission. We present a set of HST mosaic observations which confirms the presence of an alignment signal at the dark clump position. The signal strength, however, is weaker than in the ground-based data. It is therefore still unclear whether the signal is caused by a lensing mass or is just a chance alignment. We also present Chandra observations of the dark clump, which fail to reveal any significant extended emission. A comparison of the ellipticity measurements from the space-based HST data and the ground-based CFHT data shows a remarkable agreement on average, demonstrating that weak lensing studies from high-quality ground-based observations yield reliable results.