The morphology-density relationship states that dense cosmic environments such as galaxy clusters have an overabundance of quiescent elliptical galaxies, but it is unclear at which redshift this relationship is first established. We study the morphology of 4 clusters with $1.2<z<1.8$ using HST imaging and the morphology computation code statmorph. By comparing median morphology of cluster galaxies to CANDELS field galaxies using Monte Carlo analysis, we find that 2 out of 4 clusters (at z=1.19 and z=1.75) have an established morphology-density relationship with more than $3\sigma$ significance. $\sim$50% of galaxies in these clusters are bulge-dominated compared to $\sim$30% in the field, and they are significantly more compact. This result is more significant for low-mass galaxies with $\log M/M_\odot \lessapprox 10.5$, showing that low-mass galaxies are affected the most in clusters. We also find an intriguing system of two z $\approx$ 1.45 clusters at a unusually small separation 2D separation of $3'$ and 3D separation of $\approx73$ Mpc that exhibit no morphology-density relationship but have enhanced merger signatures. We conclude that the environmental mechanism responsible for the morphology-density relationship is 1) already active as early as z=1.75, 2) forms compact, bulge-dominated galaxies and 3) affects primarily low-mass galaxies. However, there is a significant degree of intracluster variance that may depend on the larger cosmological environment in which the cluster is embedded.

Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present ``total light" stacking techniques spanning 3.4-500{\mu}m aimed at revealing the total cluster IR emission, including low mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our stacking techniques to 232 well-studied massive (log M200/Msun~13.8) clusters across multiple z bins, recovering extended cluster emission at all wavelengths, typically at >5sigma. We measure the averaged near- to far-IR radial profiles and SEDs, quantifying the total stellar and dust content. The near-IR radial profiles are well described by an NFW model with a high (c~7) concentration parameter. Dust emission is similarly concentrated, albeit suppressed at small radii (r<0.2Mpc). The measured SEDs lack warm dust, consistent with the colder SEDs expected for low mass galaxies. We derive total stellar masses consistent with the theoretical Mhalo-M_star relation and specific-star formation rates that evolve strongly with redshift, echoing that of massive (log Mstar/Msun>10) cluster galaxies. Separating out the massive galaxy population reveals that the majority of cluster far-IR emission (~70-80%) is provided by the low mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low mass cluster galaxies.

We present a multi-wavelength investigation of the radio galaxy population in the galaxy cluster MOO J1506+5137 at $z$=1.09$\pm$0.03, which in previous work we identified as having multiple complex radio sources. The combined dataset used in this work includes data from the Low-Frequency Array Two-metre Sky Survey (LoTSS), NSF's Karl G. Jansky Very Large Array (VLA), the Robert C. Byrd Green Bank Telescope (GBT), the Spitzer Space Telescope, and the Dark Energy Camera Legacy Survey (DECaLS). We find that there are five radio sources which are all located within 500 kpc ($\sim$1$^{\prime}$) of the cluster center and have radio luminosities $P_{\mathrm{1.4GHz}}$ > 1.6$\times$10$^{24}$ W Hz$^{-1}$. The typical host galaxies are among the highest stellar mass galaxies in the cluster. The exceptional radio activity among the massive galaxy population appears to be linked to the dynamical state of the cluster. The galaxy distribution suggests an ongoing merger, with a subgroup found to the northwest of the main cluster. Further, two of the five sources are classified as bent-tail sources with one being a potential wide-angle tail (WAT)/hybrid morphology radio source (HyMoRS) indicating a dynamic environment. The cluster also lies in a region of the mass-richness plane occupied by other merging clusters in the Massive and Distant Clusters of WISE Survey (MaDCoWS). The data suggest that during the merger phase radio activity can be dramatically enhanced, which would contribute to the observed trend of increased radio activity in clusters with increasing redshift.

The properties of galaxy clusters as a function of redshift can be utilized as an important cosmological tool. We present initial results from a program of follow-up observations of the Sunyaev-Zeldovich effect (SZE) in high redshift galaxy clusters detected at infrared wavelengths in the Massive and Distant Clusters of \textit{WISE} Survey (MaDCoWS). Using typical on-source integration times of 3--4 hours per cluster, MUSTANG2 on the Green Bank Telescope was able to measure strong detections of SZE decrements and statistically significant masses on 14 out of 16 targets. On the remaining two, weaker (3.7$\sigma$) detections of the SZE signal and strong upper limits on the masses were obtained. In this paper we present masses and pressure profiles of each target and outline the data analysis used to recover these quantities. Of the clusters with strong detections, three show significantly flatter pressure profiles while, from the MUSTANG2 data, five others show signs of disruption at their cores. However, outside of the cores of the clusters, we were unable to detect significant amounts of asymmetry. Finally, there are indications that the relationship between optical richness used by MaDCoWS and SZE-inferred mass may be significantly flatter than indicated in previous studies.

The Massive and Distant Clusters of WISE Survey (MaDCoWS) provides a catalog of high-redshift ($0.7\lesssim z\lesssim 1.5$) infrared-selected galaxy clusters. However, the verification of the ionized intracluster medium, indicative of a collapsed and nearly virialized system, is made challenging by the high redshifts of the sample members. The main goal of this work is to test the capabilities of the Atacama Compact Array (ACA; also known as the Morita Array) Band 3 observations, centered at about 97.5 GHz, to provide robust validation of cluster detections via the thermal Sunyaev-Zeldovich (SZ) effect. Using a pilot sample that comprises ten MaDCoWS galaxy clusters, accessible to ACA and representative of the median sample richness, we infer the masses of the selected galaxy clusters and respective detection significance by means of a Bayesian analysis of the interferometric data. Our test of the "Verification with the ACA - Localization and Cluster Analysis" (VACA LoCA) program demonstrates that the ACA can robustly confirm the presence of the virialized intracluster medium in galaxy clusters previously identified in full-sky surveys. In particular, we obtain a significant detection of the SZ effect for seven out of the ten VACA LoCA clusters. We note that this result is independent of the assumed pressure profile. However, the limited angular dynamic range of the ACA in Band 3 alone, short observational integration times, and possible contamination from unresolved sources limit the detailed characterization of the cluster properties and the inference of the cluster masses within scales appropriate for the robust calibration of mass-richness scaling relations.

We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel'dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confirmed 81 candidates as galaxy clusters. We use these follow-up images and archival images to estimate photometric redshifts for 66 galaxy clusters and spectroscopic observations to obtain redshifts for 13 systems. An additional 2 galaxy clusters are confirmed using the overdensity of near-infrared galaxies only, and are presented without redshifts. We find that 15 candidates (18% of the total sample) are at redshift of $z \geq 1.0$, with a maximum confirmed redshift of $z_{\rm{max}} = 1.38 \pm 0.10$. We expect this catalog to contain every galaxy cluster with $M_{500c} > 2.6 \times 10^{14} M_\odot h^{-1}_{70}$ and $z > 0.25$ in the survey area. The mass threshold is approximately constant above $z = 0.25$, and the complete catalog has a median mass of approximately $ M_{500c} = 2.7 \times 10^{14} M_\odot h^{-1}_{70}$. Compared to previous SPT works, the increased depth of the millimeter-wave data (11.2 and 6.5 $\mu$K-arcmin at 95 and 150 GHz, respectively) makes it possible to find more galaxy clusters at high redshift and lower mass.

We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $\xi \ge 5$ and an additional 204 systems at $4<\xi<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $\xi>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($\lambda-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $\sigma$ level---with the relations intersecting at $\lambda=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses.

We present ~103 ks of Chandra observations of the galaxy cluster SPT-CLJ0528-5300 (SPT0528, z=0.768). This cluster harbors the most radio-loud (L_1.4GHz = 1.01 x 10^33 erg/s/Hz) central AGN of any cluster in the South Pole Telescope (SPT) SZ survey with available X-ray data. We find evidence of AGN-inflated cavities in the X-ray emission, which are consistent with the orientation of the jet direction revealed by ATCA radio data. The combined probability that two such depressions -- each at ~1.4-1.8sigma significance, oriented ~180 degrees apart and aligned with the jet axis -- would occur by chance is 0.1%. At >10^61 erg, the outburst in SPT0528 is among the most energetic known in the universe, and certainly the most powerful known at z>0.25. This work demonstrates that such powerful outbursts can be detected even in shallow X-ray exposures out to relatively high redshifts (z~0.8), providing an avenue for studying the evolution of extreme AGN feedback. The ratio of the cavity power (P_cav = 9.4+/-5.8 x 10^45 erg/s) to the cooling luminosity (L_cool = 1.5+/-0.5 x 10^44 erg/s) for SPT0528 is among the highest measured to date. If, in the future, additional systems are discovered at similar redshifts with equally high P_cav/L_cool ratios, it would imply that the feedback/cooling cycle was not as gentle at high redshifts as in the low-redshift universe.

We present the results from a study with NSF's Karl G. Jansky Very Large Array (VLA) to determine the radio morphologies of extended radio sources and the properties of their host galaxies in 50 massive galaxy clusters at z~1. We find a majority of the radio morphologies to be Fanaroff-Riley (FR) type IIs. By analyzing the infrared counterparts of the radio sources, we find that ~40% of the host galaxies are the candidate brightest cluster galaxy (BCG) and ~83% are consistent with being one of the top six most massive galaxies in the cluster. We investigate the role of environmental factors on the radio-loud AGN population by examining correlations between environmental and radio-galaxy properties. We find that the highest stellar mass hosts ($M_{*} \gtrsim$ 4$\times 10^{11} M_{\odot}$) are confined to the cluster center and host compact jets. There is evidence for an increase in the size of the jets with cluster-centric radius, which may be attributed to the decreased ICM pressure confinement with increasing radius. Besides this correlation, there are no other significant correlations between the properties of the radio-AGN (luminosity, morphology, or size) and environmental properties (cluster richness and location within the cluster). The fact that there are more AGN in the cluster environment than the field at this epoch, combined with the lack of strong correlation between galaxy and environmental properties, argues that the cluster environment fosters radio activity but does not solely drive the evolution of these sources at this redshift.

We constrain the evolution of the brightest cluster galaxy plus intracluster light (BCG+ICL) using an ensemble of 42 galaxy groups and clusters that span redshifts of z = 0.05-1.75 and masses of $M_{500,c}=2\times10^{13}-10^{15}$ M$_\odot$ Specifically, we measure the relationship between the BCG+ICL stellar mass $M_\star$ and $M_{500,c}$ at projected radii 10 < r < 100 kpc for three different epochs. At intermediate redshift (z = 0.40), where we have the best data, we find $M_\star\propto M_{500,c}^{0.48\pm0.06}$. Fixing the exponent of this power law for all redshifts, we constrain the normalization of this relation to be $2.08\pm0.21$ times higher at z = 0.40 than at high redshift (z = 1.55). We find no change in the relation from intermediate to low redshift (z = 0.10). In other words, for fixed $M_{500,c}$, $M_\star$ at 10 < r < 100 kpc increases from z = 1.55 to z = 0.40 and not significantly thereafter. Theoretical models predict that the physical mass growth of the cluster from z = 1.5 to z = 0 within $r_{500,c}$ is a factor of 1.4, excluding evolution due to definition of $r_{500,c}$. We find that $M_\star$ within the central 100 kpc increases by a factor of 3.8 over the same period. Thus, the growth of $M_\star$ in this central region is more than a factor of two greater than the physical mass growth of the cluster as a whole. Furthermore, the concentration of the BCG+ICL stellar mass, defined by the ratio of stellar mass within 10 kpc to the total stellar mass within 100 kpc, decreases with increasing $M_{500,c}$ at all redshift. We interpret this result as evidence for inside-out growth of the BCG+ICL over the past ten Gyrs, with stellar mass assembly occuring at larger radii at later times.

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

Here we present new red sequence overdensity measurements for 77 fields in the high-$z$ Clusters Occupied by Bent Radio AGN (COBRA) survey, based on $r$- and $i$-band imaging taken with Lowell Observatory's Discovery Channel Telescope. We observe 38 COBRA fields in $r$-band and 90 COBRA fields in $i$-band. By combining the $r$- and $i$-band photometry with our 3.6$\mu$m and 4.5$\mu$m $Spitzer$ IRAC observations, we identify 39 red sequence cluster candidates that host a strong overdensity of galaxies when measuring the excess of red sequence galaxies relative to a background field. We initially treat the radio host as the cluster center and then determine a new cluster center based on the surface density of red sequence sources. Using our color selection, we identify which COBRA cluster candidates have strong red sequence populations. By removing foreground and background contaminants, we more securely determine which fields include cluster candidates with a higher significance than our single-band observations. Additionally, of the 77 fields we analyze with a redshift estimate, 26 include newly estimated photometric redshifts.

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

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2\sigma$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.

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)

We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets.

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

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

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

We present the results from a pilot study with the Karl G. Jansky Very Large Array (JVLA) to determine the radio morphologies of extended radio sources and the properties of their host-galaxies in 10 massive galaxy clusters at z~1, an epoch in which clusters are assembling rapidly. These clusters are drawn from a parent sample of WISE-selected galaxy clusters that were cross-correlated with the VLA Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST) to identify extended radio sources within 1$^{\prime}$ of the cluster centers. Out of the ten targeted sources, six are FR II sources, one is an FR I source, and three sources have undetermined morphologies. Eight radio sources have associated Spitzer data, 75% presenting infrared counterparts. A majority of these counterparts are consistent with being massive galaxies. The angular extent of the FR sources exhibits a strong correlation with the cluster-centric radius, which warrants further investigation with a larger sample.

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

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 present an analysis of the radial distribution of Active Galactic Nuclei (AGN) in $2300$ galaxy clusters from the Massive and Distant Clusters of {\it WISE} Survey (MaDCoWS). MaDCoWS provides the largest coverage of the extragalactic sky for a cluster sample at $z\sim1$. We use literature catalogs of AGN selected via optical, mid-infrared (MIR), and radio data, and by optical-to-MIR (OIR) color. Stacking the radial distribution of AGN within the $6\arcmin$ of the centers of MaDCoWS galaxy clusters, we find a distinct overdensity of AGN within $1\arcmin$ of the galaxy cluster center for AGN of all selection methods. The fraction of red galaxies that host AGN as a function of clustercentric distance is, however, dependent on the AGN selection. The fraction of red galaxies in cluster environments that host AGN selected by optical signatures or blue OIR color is at a deficit compared to the field, while MIR-selected and red OIR color AGN are enhanced in the centers of clusters when compared to field levels. The radio-selected AGN fraction is more than $2.5$ times that of the field, implying that the centers of clusters are conducive to the triggering of radio emission in AGN. We do not find a statistically significant change in the AGN fraction as a function of cluster richness. We also investigate the correlation of central radio activity with other AGN in galaxy clusters. Clusters with radio activity have more central AGN than radio-inactive clusters, implying that central cluster radio activity and AGN triggering may be linked.

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

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

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

The incidence of Active Galactic Nuclei (AGN) with local environment is a potentially valuable probe of the mechanisms that trigger and provide fuel for accretion onto supermassive black holes. While the correlation between AGN fraction and environment has been well-studied in the local universe, AGN fractions have been measured for relatively few dense environments at high redshift. In this paper we present a measurement of the X-ray AGN fraction in the USS 1558-003 protocluster associated with the z=2.53 radio galaxy 4C-00.62. Our measurement is based on a 100ks Chandra observation, follow-up spectroscopy from the Multi-Object Double Spectrograph on the Large Binocular Telescope, and broad and narrow band photometry. These data are sensitive to AGN more luminous than $L_{X}>2\times10^{43}$ erg s$^{-1}$ in the rest-frame hard X-ray band (2-10 keV). We have identified two X-ray AGN at the redshift of USS 1558-003, one of which is the radio galaxy. We have determined that $2.0^{+2.6}_{-1.3}$% of the H$\alpha$ emitters in the protocluster are X-ray AGN. Unlike most other high-redshift cluster progenitors studied with similar techniques, USS 1558-003 does not have a significantly higher fraction of AGN than field galaxies at similar redshifts. This lower AGN fraction is inconsistent with the expectation that the higher gas fractions at high redshift, combined with the high galaxy densities and modest relative velocities in protoclusters, should produce higher AGN fractions.

We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of 2.8+0.6-0.5 (1.10 +- 0.23 mag)---compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from LambdaCDM simulations---making it the most amplified SN Ia discovered behind a galaxy cluster.

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today.

We report spectroscopic results from our 40-orbit $Hubble~Space~Telescope$ slitless grism spectroscopy program observing the 20 densest CARLA candidate galaxy clusters at $1.4 < z < 2.8$. These candidate rich structures, among the richest and most distant known, were identified on the basis of $[3.6]-[4.5]$ color from a $408~$hr multi-cycle $Spitzer$ program targeting $420$ distant radio-loud AGN. We report the spectroscopic confirmation of $16$ distant structures at $1.4 < z < 2.8$ associated with the targeted powerful high-redshift radio-loud AGN. We also report the serendipitous discovery and spectroscopic confirmation of seven additional structures at $0.87 < z < 2.12$ not associated with the targeted radio-loud AGN. We find that $10^{10} - 10^{11}\,M_{\odot}$ member galaxies of our confirmed CARLA structures form significantly fewer stars than their field counterparts at all redshifts within $1.4\leq z\leq 2$. We also observe higher star-forming activity in the structure cores up to $z = 2$, finding similar trends as cluster surveys at slightly lower redshifts ($1.0 < z < 1.5$). By design, our efficient strategy of obtaining just two grism orbits per field only obtains spectroscopic confirmation of emission-line galaxies. Deeper spectroscopy will be required to study the population of evolved, massive galaxies in these (forming) clusters. Lacking multi-band coverage of the fields, we adopt a very conservative approach of calling all confirmations "structures", although we note that a number of features are consistent with some of them being bona fide galaxy clusters. Together, this survey represents a unique and large homogenous sample of spectroscopically confirmed structures at high redshifts, potentially more than doubling the census of confirmed, massive clusters at $z > 1.4$.

According to the current understanding of cosmic structure formation, the precursors of the most massive structures in the Universe began to form shortly after the Big Bang, in regions corresponding to the largest fluctuations in the cosmic density field. Observing these structures during their period of active growth and assembly - the first few hundred million years of the Universe - is challenging because it requires surveys that are sensitive enough to detect the distant galaxies that act as signposts for these structures and wide enough to capture the rarest objects. As a result, very few such objects have been detected so far. Here we report observations of a far-infrared-luminous object at redshift 6.900 (less than 800 Myr after the Big Bang) that was discovered in a wide-field survey. High-resolution imaging reveals this source to be a pair of extremely massive star-forming galaxies. The larger of these galaxies is forming stars at a rate of 2900 solar masses per year, contains 270 billion solar masses of gas and 2.5 billion solar masses of dust, and is more massive than any other known object at a redshift of more than 6. Its rapid star formation is probably triggered by its companion galaxy at a projected separation of just 8 kiloparsecs. This merging companion hosts 35 billion solar masses of stars and has a star-formation rate of 540 solar masses per year, but has an order of magnitude less gas and dust than its neighbor and physical conditions akin to those observed in lower-metallicity galaxies in the nearby Universe. These objects suggest the presence of a dark-matter halo with a mass of more than 400 billion solar masses, making it among the rarest dark-matter haloes that should exist in the Universe at this epoch.

We present a weak-lensing study of SPT-CLJ2040-4451 and IDCSJ1426+3508 at z=1.48 and 1.75, respectively. The two clusters were observed in our "See Change" program, a HST survey of 12 massive high-redshift clusters aimed at high-z supernova measurements and weak-lensing estimation of accurate cluster masses. We detect weak but significant galaxy shape distortions using IR images from the WFC3, which has not yet been used for weak-lensing studies. Both clusters appear to possess relaxed morphology in projected mass distribution, and their mass centroids agree nicely with those defined by both the galaxy luminosity and X-ray emission. Using an NFW profile, for which we assume that the mass is tightly correlated with the concentration parameter, we determine the masses of SPT-CL J2040-4451 and IDCS J1426+3508 to be M_{200}=8.6_{-1.4}^{+1.7}x10^14 M_sun and 2.2_{-0.7}^{+1.1}x10^14 M_sun, respectively. The weak-lensing mass of SPT-CLJ2040-4451 shows that the cluster is clearly a rare object. Adopting the central value, the expected abundance of such a massive cluster at z>1.48 is only ~0.07 in the parent 2500 sq. deg. survey. However, it is yet premature to claim that the presence of this cluster creates a serious tension with the current LCDM paradigm unless that tension will remain in future studies after marginalizing over many sources of uncertainties such as the accuracy of the mass function and the mass-concentration relation at the high mass end. The mass of IDCSJ1426+3508 is in excellent agreement with our previous ACS-based weak-lensing result while the much higher source density from our WFC3 imaging data makes the current statistical uncertainty ~40% smaller.

Several recent studies have reported different intrinsic correlations between the AGN mid-IR luminosity ($L_{MIR}$) and the rest-frame 2-10 keV luminosity ($L_{X}$) for luminous quasars. To understand the origin of the difference in the observed $L_{X}-L_{MIR}$ relations, we study a sample of 3,247 spectroscopically confirmed type 1 AGNs collected from Bo\"{o}tes, XMM-COSMOS, XMM-XXL-North, and the SDSS quasars in the Swift/XRT footprint spanning over four orders of magnitude in luminosity. We carefully examine how different observational constraints impact the observed $L_{X}-L_{MIR}$ relations, including the inclusion of X-ray non-detected objects, possible X-ray absorption in type 1 AGNs, X-ray flux limits, and star formation contamination. We find that the primary factor driving the different $L_{X}-L_{MIR}$ relations reported in the literature is the X-ray flux limits for different studies. When taking these effects into account, we find that the X-ray luminosity and mid-IR luminosity (measured at rest-frame $6\mu m$, or $L_{6\mu m}$) of our sample of type 1 AGNs follow a bilinear relation in the log-log plane: $\log L_X =(0.84\pm0.03)\times\log L_{6\mu m}/10^{45}{\rm erg\;s^{-1}} + (44.60\pm0.01)$ for $L_{6\mu m} < 10^{44.79}{\rm erg\;s^{-1}} $, and $\log L_X = (0.40\pm0.03)\times\log L_{6\mu m}/10^{45}{\rm erg\;s^{-1}} +(44.51\pm0.01)$ for $L_{6\mu m} \geq 10^{44.79}{\rm erg\;s^{-1}} $. This suggests that the luminous type 1 quasars have a shallower $L_{X}-L_{MIR}$ correlation than the approximately linear relations found in local Seyfert galaxies. This result is consistent with previous studies reporting a luminosity-dependent $L_{X}-L_{MIR}$ relation, and implies that assuming a linear $L_{X}-L_{MIR}$ relation to infer the neutral gas column density for X-ray absorption might overestimate the column densities in luminous quasars.

We present 190 galaxy cluster candidates (most at high redshift) based on galaxy overdensity measurements in the \spitzer/IRAC imaging of the fields surrounding 646 bent, double-lobed radio sources drawn from the Clusters Occupied by Bent Radio AGN (COBRA) Survey. The COBRA sources were chosen as objects in the VLA FIRST survey that lack optical counterparts in the Sloan Digital Sky Survey (SDSS) to a limit of $m_r=22$, making them likely to lie at high redshift. This is confirmed by our observations: the redshift distribution of COBRA sources with estimated redshifts peaks near $z=1$, and extends out to $z\approx3$. Cluster candidates were identified by comparing our target fields to a background field and searching for statistically significant ($\ge2\sigma$) excesses in the galaxy number counts surrounding the radio sources; 190 fields satisfy the $\ge2\sigma$ limit. We find that 530 fields (82.0\%) have a net positive excess of galaxies surrounding the radio source. Many of the fields with positive excesses but below the $2\sigma$ cutoff are likely to be galaxy groups. Forty-one COBRA sources are quasars with known spectroscopic redshifts, which may be tracers of some of the most distant clusters known.

We present a study of the relation between galaxy stellar age and mass for 14 members of the $z=1.62$ protocluster IRC 0218, using multiband imaging and HST G102 and G141 grism spectroscopy. Using $UVJ$ colors to separate galaxies into star forming and quiescent populations, we find that at stellar masses $M_* \geq 10^{10.85} M_{\odot}$, the quiescent fraction in the protocluster is $f_Q=1.0^{+0.00}_{-0.37}$, consistent with a $\sim 2\times $ enhancement relative to the field value, $f_Q=0.45^{+0.03}_{-0.03}$. At masses $10^{10.2} M_{\odot} \leq M_* \leq 10^{10.85} M_{\odot}$, $f_Q$ in the cluster is $f_Q=0.40^{+0.20}_{-0.18}$, consistent with the field value of $f_Q=0.28^{+0.02}_{-0.02}$. Using galaxy $D_{n}(4000)$ values derived from the G102 spectroscopy, we find no relation between galaxy stellar age and mass. These results may reflect the impact of merger-driven mass redistribution, which is plausible as this cluster is known to host many dry mergers. Alternately, they may imply that the trend in $f_Q$ in IRC 0218 was imprinted over a short timescale in the protocluster's assembly history. Comparing our results with those of other high-redshift studies and studies of clusters at $z\sim 1$, we determine that our observed relation between $f_Q$ and stellar mass only mildly evolves between $z\sim 1.6$ and $z \sim 1$, and only at stellar masses $M_* \leq 10^{10.85} M_{\odot}$. Both the $z\sim 1$ and $z\sim 1.6$ results are in agreement that the red sequence in dense environments was already populated at high redshift, $z \ge 3$, placing constraints on the mechanism(s) responsible for quenching in dense environments at $z\ge 1.5$

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

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

We report new ALMA observations of the CO(3-2) line emission from the $2.1\pm0.3\times10^{10}\rm\thinspace M_{\odot}$ molecular gas reservoir in the central galaxy of the Phoenix cluster. The cold molecular gas is fuelling a vigorous starburst at a rate of $500-800\rm\thinspace M_{\odot}\rm\; yr^{-1}$ and powerful black hole activity in the form of both intense quasar radiation and radio jets. The radio jets have inflated huge bubbles filled with relativistic plasma into the hot, X-ray atmospheres surrounding the host galaxy. The ALMA observations show that extended filaments of molecular gas, each $10-20\rm\; kpc$ long with a mass of several billion solar masses, are located along the peripheries of the radio bubbles. The smooth velocity gradients and narrow line widths along each filament reveal massive, ordered molecular gas flows around each bubble, which are inconsistent with gravitational free-fall. The molecular clouds have been lifted directly by the radio bubbles, or formed via thermal instabilities induced in low entropy gas lifted in the updraft of the bubbles. These new data provide compelling evidence for close coupling between the radio bubbles and the cold gas, which is essential to explain the self-regulation of feedback. The very feedback mechanism that heats hot atmospheres and suppresses star formation may also paradoxically stimulate production of the cold gas required to sustain feedback in massive galaxies.

We explore 7.5 billion years of evolution in the star formation activity of massive ($M_{\star}>10^{10.1}\,M_{\odot}$) cluster galaxies using a sample of 25 clusters over $0.15<z<1$ from the Cluster Lensing And Supernova survey with Hubble and 11 clusters over $1<z<1.5$ from the IRAC Shallow Cluster Survey. Galaxy morphologies are determined visually using high-resolution Hubble Space Telescope images. Using the spectral energy distribution fitting code CIGALE, we measure star formation rates, stellar masses, and 4000 \AA\ break strengths. The latter are used to separate quiescent and star-forming galaxies (SFGs). From $z\sim1.3$ to $z\sim0.2$, the specific star formation rate (sSFR) of cluster SFGs and quiescent galaxies decreases by factors of three and four, respectively. Over the same redshift range, the sSFR of the entire cluster population declines by a factor of 11, from $0.48\pm0.06\;\mathrm{Gyr}^{-1}$ to $0.043\pm0.009\;\mathrm{Gyr}^{-1}$. This strong overall sSFR evolution is driven by the growth of the quiescent population over time; the fraction of quiescent cluster galaxies increases from $28^{+8}_{-19}\%$ to $88^{+5}_{-4}\%$ over $z\sim1.3\rightarrow0.2$. The majority of the growth occurs at $z\gtrsim0.9$, where the quiescent fraction increases by 0.41. While the sSFR of the majority of star-forming cluster galaxies is at the level of the field, a small subset of cluster SFGs have low field-relative star formation activity, suggestive of long-timescale quenching. The large increase in the fraction of quiescent galaxies above $z\sim0.9$, coupled with the field-level sSFRs of cluster SFGs, suggests that higher redshift cluster galaxies are likely being quenched quickly. Assessing those timescales will require more accurate stellar population ages and star formation histories.

We present the stellar mass - size relation for 49 galaxies within the $z$ = 1.067 cluster SPT-CL J0546$-$5345, with FWHM $\sim$80-120 mas $K_{\mathrm s}$-band data from the Gemini multi-conjugate adaptive optics system (GeMS/GSAOI). This is the first such measurement in a cluster environment, performed at sub-kpc resolution at rest-frame wavelengths dominated by the light of the underlying old stellar populations. The observed stellar mass - size relation is offset from the local relation by 0.21 dex, corresponding to a size evolution proportional to $(1+z)^{-1.25}$, consistent with the literature. The slope of the stellar mass - size relation $\beta$ = 0.74 $\pm$ 0.06, consistent with the local relation. The absence of slope evolution indicates that the amount of size growth is constant with stellar mass. This suggests that galaxies in massive clusters such as SPT-CL J0546$-$5345 grow via processes that increase the size without significant morphological interference, such as minor mergers and/or adiabatic expansion. The slope of the cluster stellar mass - size relation is significantly shallower if measured in $HST$/ACS imaging at wavelengths blueward of the Balmer break, similar to rest-frame UV relations at $z$ = 1 in the literature. The stellar mass - size relation must be measured at redder wavelengths, which are more sensitive to the old stellar population that dominates the stellar mass of the galaxies. The slope is unchanged when GeMS $K_s$-band imaging is degraded to the resolution of $K$-band HST/NICMOS resolution but dramatically affected when degraded to $K_s$-band Magellan/FourStar resolution. Such measurements must be made with AO in order to accurately characterise the sizes of compact, $z$ = 1 galaxies.

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.

We use the Wide-field Infrared Survey Explorer (WISE) and the Sloan Digital Sky Survey (SDSS) to confirm a connection between dust-obscured active galactic nuclei (AGNs) and galaxy merging. Using a new, volume-limited (z <= 0.08) catalog of visually-selected major mergers and galaxy-galaxy interactions from the SDSS, with stellar masses above 2x10^10 Msun, we find that major mergers (interactions) are 5-17 (3-5) times more likely to have red [3.4]-[4.6] colors associated with dust-obscured or `dusty' AGNs, compared to non-merging galaxies with similar masses. Using published fiber spectral diagnostics, we map the [3.4]-[4.6] versus [4.6]-[12] colors of different emission-line galaxies and find one-quarter of Seyferts have colors indicative of a dusty AGN. We find that AGNs are five times more likely to be obscured when hosted by a merging galaxy, half of AGNs hosted by a merger are dusty, and we find no enhanced frequency of optical AGNs in merging over non-merging galaxies. We conclude that undetected AGNs missed at shorter wavelengths are at the heart of the ongoing AGN-merger connection debate. The vast majority of mergers hosting dusty AGNs are star-forming and located at the centers of Mhalo<10^13 Msun groups. Assuming plausibly short duration dusty-AGN phases, we speculate that a large fraction of gas-rich mergers experience a brief obscured AGN phase, in agreement with the strong connection between central star formation and black hole growth seen in merger simulations.

Using HST slitless grism data, we report the spectroscopic confirmation of two distant structures at $z \sim 2$ associated with powerful high-redshift radio-loud AGN. These rich structures, likely (forming) clusters, are among the most distant currently known and were identified on the basis of Spitzer/IRAC [3.6] - [4.5] color. We spectroscopically confirm 9 members in the field of MRC 2036-254, comprising eight star-forming galaxies and the targeted radio galaxy. The median redshift is $z = 2.000$. We spectroscopically confirm 10 members in the field of B3 0756+406, comprising eight star-forming galaxies and two AGN, including the targeted radio-loud quasar. The median redshift is $z = 1.986$. All confirmed members are within 500 kpc (1 arcmin) of the targeted AGN. We derive median (mean) star-formation rates of $\sim 35~M_{\odot}\rm ~ yr^{-1}$ ($\sim 50~M_{\odot}\rm ~ yr^{-1}$) for the confirmed star-forming members of both structures based on their [OIII]$\lambda5007$ luminosities, and estimate average galaxy stellar masses $\lesssim 1 \times 10^{11} ~M_{\odot}$ based on mid-infrared fluxes and SED modeling. Most of our confirmed members are located above the star-forming main-sequence towards starburst galaxies, consistent with clusters at these early epochs being the sites of significant levels of star formation. The structure around MRC 2036-254 shows an overdensity of IRAC-selected candidate galaxy cluster members consistent with being quiescent galaxies, while the structure around B3 0756+406 shows field values, albeit with many lower limits to colors that could allow an overdensity of faint red quiescent galaxies. The structure around MRC 2036-254 shows a red sequence of passive galaxy candidates.

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

We present a detailed, multi-wavelength study of star formation (SF) and AGN activity in 11 near-infrared (IR) selected, spectroscopically confirmed, massive ($\gtrsim10^{14}\,\rm{M_{\odot}}$) galaxy clusters at $1<z<1.75$. Using new, deep $Herschel$/PACS imaging, we characterize the optical to far-IR spectral energy distributions (SEDs) for IR-luminous cluster galaxies, finding that they can, on average, be well described by field galaxy templates. Identification and decomposition of AGN through SED fittings allows us to include the contribution to cluster SF from AGN host galaxies. We quantify the star-forming fraction, dust-obscured SF rates (SFRs), and specific-SFRs for cluster galaxies as a function of cluster-centric radius and redshift. In good agreement with previous studies, we find that SF in cluster galaxies at $z\gtrsim1.4$ is largely consistent with field galaxies at similar epochs, indicating an era before significant quenching in the cluster cores ($r<0.5\,$Mpc). This is followed by a transition to lower SF activity as environmental quenching dominates by $z\sim1$. Enhanced SFRs are found in lower mass ($10.1< \log \rm{M_{\star}}/\rm{M_{\odot}}<10.8$) cluster galaxies. We find significant variation in SF from cluster-to-cluster within our uniformly selected sample, indicating that caution should be taken when evaluating individual clusters. We examine AGN in clusters from $z=0.5-2$, finding an excess AGN fraction at $z\gtrsim1$, suggesting environmental triggering of AGN during this epoch. We argue that our results $-$ a transition from field-like to quenched SF, enhanced SF in lower mass galaxies in the cluster cores, and excess AGN $-$ are consistent with a co-evolution between SF and AGN in clusters and an increased merger rate in massive haloes at high redshift.

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

The Decadal IRAC Bootes Survey is a mid-IR variability survey of the ~9 sq. deg. of the NDWFS Bootes Field and extends the time baseline of its predecessor, the Spitzer Deep, Wide-Field Survey (SDWFS), from 4 to 10 years. The Spitzer Space Telescope visited the field five times between 2004 and 2014 at 3.6 and 4.5 microns. We provide the difference image analysis photometry for a half a million mostly extragalactic sources. In mid-IR color-color plane, sources with quasar colors constitute the largest variability class (75%), 16% of the variable objects have stellar colors and the remaining 9% have the colors of galaxies. Adding the fifth epoch doubles the number of variable active galactic nuclei (AGNs) for the same false positive rates as in SDWFS, or increases the number of sources by 20% while decreasing the false positive rates by factors of 2-3 for the same variability amplitude. We quantify the ensemble mid-IR variability of ~1500 spectroscopically confirmed AGNs using single power-law structure functions (SFs), which we find to be steeper (index $\gamma=0.45$) than in the optical ($\gamma=0.3$), leading to much lower amplitudes at short time-lags. This provides evidence for large emission regions, smoothing out any fast UV/optical variations, as the origin of infrared quasar variability. The mid-IR AGN SF slope $\gamma$ seems to be uncorrelated with both the luminosity and rest-frame wavelength, while the amplitude shows an anti-correlation with the luminosity and a correlation with the rest-frame wavelength.

We present a deep (100 ks) Chandra observation of IDCS J1426.5+3508, a spectroscopically confirmed, infrared-selected galaxy cluster at $z = 1.75$. This cluster is the most massive galaxy cluster currently known at $z > 1.5$, based on existing Sunyaev-Zel'dovich (SZ) and gravitational lensing detections. We confirm this high mass via a variety of X-ray scaling relations, including $T_X$-M, $f_g$-M, $Y_X$-M and $L_X$-M, finding a tight distribution of masses from these different methods, spanning M$_{500}$ = 2.3-3.3 $\times 10^{14}$ M$_{\odot}$, with the low-scatter $Y_X$-based mass $M_{500,Y_X} = 2.6^{+1.5}_{-0.5} \times 10^{14}$ M$_\odot$. IDCS J1426.5+3508 is currently the only cluster at $z > 1.5$ for which X-ray, SZ and gravitational lensing mass estimates exist, and these are in remarkably good agreement. We find a relatively tight distribution of the gas-to-total mass ratio, employing total masses from all of the aforementioned indicators, with values ranging from $f_{gas,500}$ = 0.087-0.12. We do not detect metals in the intracluster medium (ICM) of this system, placing a 2$\sigma$ upper limit of $Z(r < R_{500}) < 0.18 Z_{\odot}$. This upper limit on the metallicity suggests that this system may still be in the process of enriching its ICM. The cluster has a dense, low-entropy core, offset by $\sim$30 kpc from the X-ray centroid, which makes it one of the few "cool core" clusters discovered at $z > 1$, and the first known cool core cluster at $z > 1.2$. The offset of this core from the large-scale centroid suggests that this cluster has had a relatively recent ($\lesssim$500 Myr) merger/interaction with another massive system.

We report on the discovery of a z=1.58 mature cluster around the high-redshift radio galaxy 7C1753+6311, first identified in the Clusters Around Radio-Loud AGN survey. Two-thirds of the excess galaxies within the central 1Mpc lie on a red sequence with a colour that is consistent with an average formation redshift of zf~3. We show that 80+/-6% of the red sequence galaxies in the cluster core are quiescent, while the remaining 20% are red due to dusty star formation. We demonstrate that the cluster has an enhanced quiescent galaxy fraction that is three times that of the control field. We also show that this enhancement is mass dependent: 91+/-9% of the M* >10^{10.5}Msun cluster galaxies are quiescent, compared to only 36+/-2% of field galaxies, whereas the fraction of quiescent galaxies with lower masses is the same in the cluster and field environments. The presence of a dense core and a well-formed, quiescent red sequence suggest that this is a mature cluster. This means that distant radio galaxies do not solely reside in young, uncollapsed protoclusters, rather they can be found in clusters in a wide range of evolutionary states.

We probe star formation in the environments of massive $\sim10^{13}\,M_{\odot}$ dark matter halos at redshifts of $z$$\sim$$1$. This star formation is linked to a sub-millimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high-redshift ($z$$>$$2$) strongly-lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg$^2$ survey. The clustering signal has sub-millimetre colours which are consistent with the mean redshift of the foreground lensing halos ($z$$\sim$$1$). We report a mean excess of star formation rate (SFR) compared to the field, of $(2700\pm700)\,M_{\odot}\,{yr}^{-1}$ from all galaxies contributing to this clustering signal within a radius of 3.5' from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80$\%$ of the excess emission measured by Planck originates from galaxies lying in the neighbouring halos of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be ($370\pm40)$$\,M_{\odot}\,{yr}^{-1}$ per resolved, clustered source. Our findings suggest that the environments around these massive $z$$\sim$$1$ lensing halos host intense star formation out to about $2\,$Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

We study the stellar, Brightest Cluster Galaxy (BCG) and intracluster medium (ICM) masses of 14 South Pole Telescope (SPT) selected galaxy clusters with median redshift $z=0.9$ and median mass $M_{500}=6\times10^{14}M_{\odot}$. We estimate stellar masses for each cluster and BCG using six photometric bands spanning the range from the ultraviolet to the near-infrared observed with the VLT, HST and Spitzer. The ICM masses are derived from Chandra and XMM-Newton X-ray observations, and the virial masses are derived from the SPT Sunyaev-Zel'dovich Effect signature. At $z=0.9$ the BCG mass $M_{\star}^{\textrm{BCG}}$ constitutes $0.12\pm0.01$% of the halo mass for a $6\times10^{14}M_{\odot}$ cluster, and this fraction falls as $M_{500}^{-0.58\pm0.07}$. The cluster stellar mass function has a characteristic mass $M_{0}=10^{11.0\pm0.1}M_{\odot}$, and the number of galaxies per unit mass in clusters is larger than in the field by a factor $1.65\pm0.2$. Both results are consistent with measurements on group scales and at lower redshift. We combine our SPT sample with previously published samples at low redshift that we correct to a common initial mass function and for systematic differences in virial masses. We then explore mass and redshift trends in the stellar fraction (fstar), the ICM fraction (fICM), the cold baryon fraction (fc) and the baryon fraction (fb). At a pivot mass of $6\times10^{14}M_{\odot}$ and redshift $z=0.9$, the characteristic values are fstar=$1.1\pm0.1$%, fICM=$9.6\pm0.5$%, fc=$10.4\pm1.2$% and fb=$10.7\pm0.6$%. These fractions all vary with cluster mass at high significance, indicating that higher mass clusters have lower fstar and fc and higher fICM and fb. When accounting for a 15% systematic virial mass uncertainty, there is no statistically significant redshift trend at fixed mass in these baryon fractions. (abridged)

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

To understand cosmic mass assembly in the Universe at early epochs, we primarily rely on measurements of stellar mass and star formation rate of distant galaxies. In this paper, we present stellar masses and star formation rates of six high-redshift ($2.8\leq z \leq 5.7$) dusty, star-forming galaxies (DSFGs) that are strongly gravitationally lensed by foreground galaxies. These sources were first discovered by the South Pole Telescope (SPT) at millimeter wavelengths and all have spectroscopic redshifts and robust lens models derived from ALMA observations. We have conducted follow-up observations, obtaining multi-wavelength imaging data, using {\it HST}, {\it Spitzer}, {\it Herschel} and the Atacama Pathfinder EXperiment (APEX). We use the high-resolution {\it HST}/WFC3 images to disentangle the background source from the foreground lens in {\it Spitzer}/IRAC data. The detections and upper limits provide important constraints on the spectral energy distributions (SEDs) for these DSFGs, yielding stellar masses, IR luminosities, and star formation rates (SFRs). The SED fits of six SPT sources show that the intrinsic stellar masses span a range more than one order of magnitude with a median value $\sim$ 5 $\times 10^{10}M_{\Sun}$. The intrinsic IR luminosities range from 4$\times 10^{12}L_{\Sun}$ to 4$\times 10^{13}L_{\Sun}$. They all have prodigious intrinsic star formation rates of 510 to 4800 $M_{\Sun} {\rm yr}^{-1}$. Compared to the star-forming main sequence (MS), these six DSFGs have specific SFRs that all lie above the MS, including two galaxies that are a factor of 10 higher than the MS. Our results suggest that we are witnessing the ongoing strong starburst events which may be driven by major mergers.

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

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

We use a sample of 37 of the densest clusters and protoclusters across $1.3 \le z \le 3.2$ from the Clusters Around Radio-Loud AGN (CARLA) survey to study the formation of massive cluster galaxies. We use optical $i'$-band and infrared 3.6$\mu$m and 4.5$\mu$m images to statistically select sources within these protoclusters and measure their median observed colours; $\langle i'-[3.6] \rangle$. We find the abundance of massive galaxies within the protoclusters increases with decreasing redshift, suggesting these objects may form an evolutionary sequence, with the lower redshift clusters in the sample having similar properties to the descendants of the high redshift protoclusters. We find that the protocluster galaxies have an approximately unevolving observed-frame $i'-[3.6]$ colour across the examined redshift range. We compare the evolution of the $\langle i'-[3.6] \rangle$ colour of massive cluster galaxies with simplistic galaxy formation models. Taking the full cluster population into account, we show that the formation of stars within the majority of massive cluster galaxies occurs over at least 2Gyr, and peaks at $z \sim 2$-3. From the median $i'-[3.6]$ colours we cannot determine the star formation histories of individual galaxies, but their star formation must have been rapidly terminated to produce the observed red colours. Finally, we show that massive galaxies at $z>2$ must have assembled within 0.5Gyr of them forming a significant fraction of their stars. This means that few massive galaxies in $z>2$ protoclusters could have formed via dry mergers.

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.

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

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

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

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

We measure the star formation rates (SFRs) of massive ($M_{\star}>10^{10.1}M_{\odot}$) early-type galaxies (ETGs) in a sample of 11 high-redshift ($1.0 < z < 1.5$) galaxy clusters drawn from the IRAC Shallow Cluster Survey (ISCS). We identify ETGs visually from Hubble Space Telescope imaging and select likely cluster members as having either an appropriate spectroscopic redshift or red sequence color. Mid-infrared SFRs are measured using Spitzer 24 $\mu$m data for isolated cluster galaxies for which contamination by neighbors, and active galactic nuclei, can be ruled out. Cluster ETGs show enhanced specific star formation rates (sSFRs) compared to cluster galaxies in the local Universe, but have sSFRs more than four times lower than that of field ETGs at $1 < z < 1.5$. Relative to the late-type cluster population, isolated ETGs show substantially quenched mean SFRs, yet still contribute 12% of the overall star formation activity measured in $1 < z < 1.5$ clusters. We find that new ETGs are likely being formed in ISCS clusters; the fraction of cluster galaxies identified as ETGs increases from 34% to 56% from $z \sim 1.5 \rightarrow 1.25$. While the fraction of cluster ETGs that are highly star-forming ($\textrm{SFR}\geq26\ M_{\odot}$ yr$^{-1}$) drops from 27% to 10% over the same period, their sSFRs are roughly constant. All these factors taken together suggest that, particularly at $z\gtrsim1.25$, the events that created these distant cluster ETGs$-$likely mergers, at least among the most massive$-$were both recent and gas-rich.

We present a velocity dispersion-based mass calibration of the South Pole Telescope Sunyaev-Zel'dovich effect survey (SPT-SZ) galaxy cluster sample. Using a homogeneously selected sample of 100 cluster candidates from 720 deg2 of the survey along with 63 velocity dispersion ($\sigma_v$) and 16 X-ray Yx measurements of sample clusters, we simultaneously calibrate the mass-observable relation and constrain cosmological parameters. The calibrations using $\sigma_v$ and Yx are consistent at the $0.6\sigma$ level, with the $\sigma_v$ calibration preferring ~16% higher masses. We use the full cluster dataset to measure $\sigma_8(\Omega_ m/0.27)^{0.3}=0.809\pm0.036$. The SPT cluster abundance is lower than preferred by either the WMAP9 or Planck+WMAP9 polarization (WP) data, but assuming the sum of the neutrino masses is $\sum m_\nu=0.06$ eV, we find the datasets to be consistent at the 1.0$\sigma$ level for WMAP9 and 1.5$\sigma$ for Planck+WP. Allowing for larger $\sum m_\nu$ further reconciles the results. When we combine the cluster and Planck+WP datasets with BAO and SNIa, the preferred cluster masses are $1.9\sigma$ higher than the Yx calibration and $0.8\sigma$ higher than the $\sigma_v$ calibration. Given the scale of these shifts (~44% and ~23% in mass, respectively), we execute a goodness of fit test; it reveals no tension, indicating that the best-fit model provides an adequate description of the data. Using the multi-probe dataset, we measure $\Omega_ m=0.299\pm0.009$ and $\sigma_8=0.829\pm0.011$. Within a $\nu$CDM model we find $\sum m_\nu = 0.148\pm0.081$ eV. We present a consistency test of the cosmic growth rate. Allowing both the growth index $\gamma$ and the dark energy equation of state parameter $w$ to vary, we find $\gamma=0.73\pm0.28$ and $w=-1.007\pm0.065$, demonstrating that the expansion and the growth histories are consistent with a LCDM model ($\gamma=0.55; \,w=-1$).

We present an analysis of the clustering of high-redshift galaxies in the recently completed 94 deg$^2$ Spitzer-SPT Deep Field survey. Applying flux and color cuts to the mid-infrared photometry efficiently selects galaxies at $z\sim1.5$ in the stellar mass range $10^{10}-10^{11}M_\odot$, making this sample the largest used so far to study such a distant population. We measure the angular correlation function in different flux-limited samples at scales $>6^{\prime \prime}$ (corresponding to physical distances $>0.05$ Mpc) and thereby map the one- and two-halo contributions to the clustering. We fit halo occupation distributions and determine how the central galaxy's stellar mass and satellite occupation depend on the halo mass. We measure a prominent peak in the stellar-to-halo mass ratio at a halo mass of $\log(M_{\rm halo} / M_\odot) = 12.44\pm0.08$, 4.5 times higher than the $z=0$ value. This supports the idea of an evolving mass threshold above which star formation is quenched. We estimate the large-scale bias in the range $b_g=2-4$ and the satellite fraction to be $f_\mathrm{sat}\sim0.2$, showing a clear evolution compared to $z=0$. We also find that, above a given stellar mass limit, the fraction of galaxies that are in similar mass pairs is higher at $z=1.5$ than at $z=0$. In addition, we measure that this fraction mildly increases with the stellar mass limit at $z=1.5$, which is the opposite of the behavior seen at low-redshift.

We are conducting a large survey of distant clusters of galaxies using radio sources with bent jets and lobes as tracers. These radio sources are driven by AGN and achieve their bent morphologies through interaction with the surrounding gas found in clusters of galaxies. Based on low-redshift studies, these types of sources can be used to identify clusters very efficiently. We present initial results from our survey of 653 bent-double radio sources with optical hosts too faint to appear in the SDSS. The sample was observed in the infrared with Spitzer, and it has revealed $\sim$200 distant clusters or proto-clusters in the redshift range $z\sim0.7 - 3.0$. The sample of bent-doubles contains both quasars and radio galaxies enabling us to study both radiative and kinetic mode feedback in cluster and group environments at a wide range of redshifts.

We present clustering measurements and halo masses of star forming galaxies at 0.2 < z < 1.0. After excluding AGN, we construct a sample of 22553 24 {\mu}m sources selected from 8.42 deg^2 of the Spitzer MIPS AGN and Galaxy Evolution Survey of Bo\"otes. Mid-infrared imaging allows us to observe galaxies with the highest star formation rates (SFRs), less biased by dust obscuration afflicting the optical bands. We find that the galaxies with the highest SFRs have optical colors which are redder than typical blue cloud galaxies, with many residing within the green valley. At z > 0.4 our sample is dominated by luminous infrared galaxies (LIRGs, L_TIR > 10^11 Lsun) and is comprised entirely of LIRGs and ultra-luminous infrared galaxies (ULIRGs, L_TIR > 10^12 Lsun) at z > 0.6. We observe weak clustering of r_0 = 3-6 Mpc/h for almost all of our star forming samples. We find that the clustering and halo mass depend on L_TIR at all redshifts, where galaxies with higher L_TIR (hence higher SFRs) have stronger clustering. Galaxies with the highest SFRs at each redshift typically reside within dark matter halos of M_halo ~ 10^12.9 Msun/h. This is consistent with a transitional halo mass, above which star formation is largely truncated, although we cannot exclude that ULIRGs reside within higher mass halos. By modeling the clustering evolution of halos, we connect our star forming galaxy samples to their local descendants. Most star forming galaxies at z < 1.0 are the progenitors of L < 2.5L* blue galaxies in the local universe, but star forming galaxies with the highest SFRs (L_TIR >10^11.7 Lsun) at 0.6<z<1.0 are the progenitors of early-type galaxies in denser group environments.

X-ray cavities are key tracers of mechanical (or radio mode) heating arising from the active galactic nuclei (AGN) in brightest cluster galaxies. We report on a survey for X-ray cavities in 83 massive, high-redshift (0.4<z<1.2) clusters of galaxies selected by their Sunyaev-Zel'dovich signature in the South Pole Telescope data. Based on Chandra X-ray images, we find a total of 6 clusters having symmetric pairs of surface brightness depressions consistent with the picture of radio jets inflating X-ray cavities in the intracluster medium. The majority of these detections are of relatively low significance and require deeper follow-up data in order to be confirmed. Further, due to the limitations of Chandra at high redshift, this search misses small (<10 kpc), unresolved X-ray cavities at high (z>0.5) redshift. Despite these limitations, our results suggest that the power generated by AGN feedback in brightest cluster galaxies has remained unchanged for over half of the age of the Universe (>7 Gyrs at z=0.8). On average, the detected X-ray cavities have powers of 0.8-5*10^45 erg/s, enthalpies of 3-6*10^59 erg, and radii of 17 kpc. Integrating over 7 Gyrs, we find that the supermassive black holes in the brightest cluster galaxies may have accreted 10^8 to several 10^9M_sun of material to power these outflows. This level of accretion indicates that significant supermassive black hole growth may occur not only at early times, in the quasar era, but at late times as well. We also find that X-ray cavities at high-redshift may inject an excess heat of 0.1-1.0 keV per particle into the hot intracluster medium above and beyond the energy needed to offset cooling. This value is similar to the energy needed to preheat clusters, break self-similarity, and explain the excess entropy in hot atmospheres.

We present results from modeling the optical spectra of a large sample of quiescent galaxies between 0.1 < z < 0.7 from the Sloan Digital Sky Survey (SDSS) and the AGN and Galaxy Evolution Survey (AGES). We examine how the stellar ages and abundance patterns of galaxies evolve over time as a function of stellar mass from 10^{9.6}-10^{11.8} Msun. Galaxy spectra are stacked in bins of mass and redshift, and modeled over a wavelength range from 4000 A to 5500 A. Full spectrum stellar population synthesis modeling provides estimates of the age and the abundances of the elements Fe, Mg, C, N, and Ca. We find negligible evolution in elemental abundances at fixed stellar mass over roughly 7 Gyr of cosmic time. In addition, the increase in stellar ages with time for massive galaxies is consistent with passive evolution since z = 0.7. Taken together, these results favor a scenario in which the inner ~ 0.3-3 R_e of massive quiescent galaxies have been passively evolving over the last half of cosmic time. Interestingly, the derived stellar ages are considerably younger than the age of the universe at all epochs, consistent with an equivalent single-burst star formation epoch of z < 1.5. These young stellar population ages coupled with the existence of massive quiescent galaxies at z > 1 indicate the inhomogeneous nature of the z < 0.7 quiescent population. The data also permit the addition of newly-quenched galaxies at masses below ~10^{10.5} Msun at z < 0.7. Additionally, we analyze very deep Keck DEIMOS spectra of the two brightest quiescent galaxies in a cluster at z = 0.83. There is tentative evidence that these galaxies are older than their counterparts in low-density environments. In the Appendix, we demonstrate that our full spectrum modeling technique allows for accurate and reliable modeling of galaxy spectra to low S/N (~20 A^{-1}) and/or low spectral resolution (R ~ 500).

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

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

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

We explore the connection between absorption by neutral gas and extinction by dust in mid-infrared (IR) selected luminous quasars. We use a sample of 33 quasars at redshifts 0.7 < z < 3 in the 9 deg^2 Bo\"otes multiwavelength survey field that are selected using Spitzer Space Telescope Infrared Array Camera colors and are well-detected as luminous X-ray sources (with >150 counts) in Chandra observations. We divide the quasars into dust-obscured and unobscured samples based on their optical to mid-IR color, and measure the neutral hydrogen column density N_H through fitting of the X-ray spectra. We find that all subsets of quasars have consistent power law photon indices equal to 1.9 that are uncorrelated with N_H. We classify the quasars as gas-absorbed or gas-unabsorbed if N_H > 10^22 cm^-2 or N_H < 10^22 cm^-2, respectively. Of 24 dust-unobscured quasars in the sample, only one shows clear evidence for significant intrinsic N_H, while 22 have column densities consistent with N_H < 10^22 cm^-2. In contrast, of the nine dust-obscured quasars, six show evidence for intrinsic gas absorption, and three are consistent with N_H < 10^22 cm^-2. We conclude that dust extinction in IR-selected quasars is strongly correlated with significant gas absorption as determined through X-ray spectral fitting. These results suggest that obscuring gas and dust in quasars are generally co-spatial, and confirm the reliability of simple mid-IR and optical photometric techniques for separating quasars based on obscuration.

We present optical and infrared imaging and optical spectroscopy of galaxy clusters which were identified as part of an all-sky search for high-redshift galaxy clusters, the Massive and Distant Clusters of WISE Survey (MaDCoWS). The initial phase of MaDCoWS combined infrared data from the all-sky data release of the Wide-field Infrared Survey Explorer (WISE) with optical data from the Sloan Digital Sky Survey (SDSS) to select probable z ~ 1 clusters of galaxies over an area of 10,000 deg^2. Our spectroscopy confirms 19 new clusters at 0.7 < z < 1.3, half of which are at z > 1, demonstrating the viability of using WISE to identify high-redshift galaxy clusters. The next phase of MaDCoWS will use the greater depth of the AllWISE data release to identify even higher redshift cluster candidates.

We identify a strong lensing galaxy in the cluster IRC 0218 (also known as XMM-LSS J02182$-$05102) that is spectroscopically confirmed to be at $z=1.62$, making it the highest-redshift strong lens galaxy known. The lens is one of the two brightest cluster galaxies and lenses a background source galaxy into an arc and a counterimage. With Hubble Space Telescope (HST) grism and Keck/LRIS spectroscopy, we measure the source redshift to be $z_{\rm S}=2.26$. Using HST imaging in ACS/F475W, ACS/F814W, WFC3/F125W, and WFC3/F160W, we model the lens mass distribution with an elliptical power-law profile and account for the effects of the cluster halo and nearby galaxies. The Einstein radius is $\theta_{\rm E}=0.38^{+0.02}_{-0.01}$" ($3.2_{-0.1}^{+0.2}$ kpc) and the total enclosed mass is M$_{\rm tot} (< \theta_{\rm E})=1.8^{+0.2}_{-0.1}\times10^{11}~{\rm M}_{\odot}$. We estimate that the cluster environment contributes $\sim10$% of this total mass. Assuming a Chabrier IMF, the dark matter fraction within $\theta_{{\rm E}}$ is $f_{\rm DM}^{{\rm Chab}} = 0.3_{-0.3}^{+0.1}$, while a Salpeter IMF is marginally inconsistent with the enclosed mass ($f_{\rm DM}^{{\rm Salp}} = -0.3_{-0.5}^{+0.2}$). The total magnification of the source is $\mu_{\rm tot}=2.1_{-0.3}^{+0.4}$. The source has at least one bright compact region offset from the source center. Emission from Ly$\alpha$ and [O III] are likely to probe different regions in the source.

We present 279 galaxy cluster candidates at $z > 1.3$ selected from the 94 deg$^{2}$ Spitzer South Pole Telescope Deep Field (SSDF) survey. We use a simple algorithm to select candidate high-redshift clusters of galaxies based on Spitzer/IRAC mid-infrared data combined with shallow all-sky optical data. We identify distant cluster candidates in SSDF adopting an overdensity threshold that results in a high purity (80%) cluster sample based on tests in the Spitzer Deep, Wide-Field Survey of the Bo\"otes field. Our simple algorithm detects all three $1.4 < z \leq 1.75$ X-ray detected clusters in the Bo\"otes field. The uniqueness of the SSDF survey resides not just in its area, one of the largest contiguous extragalactic fields observed with Spitzer, but also in its deep, multi-wavelength coverage by the South Pole Telescope (SPT), Herschel/SPIRE and XMM-Newton. This rich dataset will allow direct or stacked measurements of Sunyaev-Zel'dovich effect decrements or X-ray masses for many of the SSDF clusters presented here, and enable systematic study of the most distant clusters on an unprecedented scale. We measure the angular correlation function of our sample and find that these candidates show strong clustering. Employing the COSMOS/UltraVista photometric catalog in order to infer the redshift distribution of our cluster selection, we find that these clusters have a comoving number density $n_c = (0.7^{+6.3}_{-0.6}) \times 10^{-7} h^{3} \mathrm{Mpc}^{-3}$ and a spatial clustering correlation scale length $r_0 = (32 \pm 7) h^{-1} \rm{Mpc}$. Assuming our sample is comprised of dark matter halos above a characteristic minimum mass, $M_{{\rm min}}$, we derive that at $z=1.5$ these clusters reside in halos larger than $M_{{\rm min}} = 1.5^{+0.9}_{-0.7} \times 10^{14} h^{-1} M_{\odot}$. (abridged)

The Blanco Cosmology Survey is 4-band (griz) optical-imaging survey that covers ~80 square degrees of the southern sky. The survey consists of two fields roughly centered at (RA,DEC) = (23h,-55d) and (5h30m,-53d) with imaging designed to reach depths sufficient for the detection of L* galaxies out to a redshift of one. In this paper we describe the reduction of the survey data, the creation of calibrated source catalogs and a new method for the separation of stars and galaxies. We search these catalogs for galaxy clusters at z< 0.75 by identifying spatial over-densities of red-sequence galaxies. We report the coordinates, redshift, and optical richness, Lambda, for 764 detected galaxy clusters at z < 0.75. This sample, >85% of which are new discoveries, has a median redshift of 0.52 and median richness Lambda(0.4L*) of 16.4. Accompanying this paper we also release data products including the reduced images and calibrated source catalogs. These products are available at http://data.rcc.uchicago.edu/dataset/blanco-cosmology-survey .

We present 4.5 {\mu}m luminosity functions for galaxies identified in 178 candidate galaxy clusters at 1.3 < z < 3.2. The clusters were identified as Spitzer/IRAC color-selected overdensities in the Clusters Around Radio-Loud AGN (CARLA) project, which imaged 421 powerful radio-loud AGN at z > 1.3. The luminosity functions are derived for different redshift and richness bins, and the IRAC imaging reaches depths of m*+2, allowing us to measure the faint end slopes of the luminosity functions. We find that {\alpha} = -1 describes the luminosity function very well in all redshifts bins and does not evolve significantly. This provides evidence that the rate at which the low mass galaxy population grows through star formation, gets quenched and is replenished by in-falling field galaxies does not have a major net effect on the shape of the luminosity function. Our measurements for m* are consistent with passive evolution models and high formation redshifts z_f ~ 3. We find a slight trend towards fainter m* for the richest clusters, implying that the most massive clusters in our sample could contain older stellar populations, yet another example of cosmic downsizing. Modelling shows that a contribution of a star-forming population of up to 40% cannot be ruled out. This value, found from our targeted survey, is significantly lower than the values found for slightly lower redshift, z ~ 1, clusters found in wide-field surveys. The results are consistent with cosmic downsizing, as the clusters studied here were all found in the vicinity of radio-loud AGNs -- which have proven to be preferentially located in massive dark matter halos in the richest environments at high redshift -- they may therefore be older and more evolved systems than the general protocluster population.

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.

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

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

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $\mu$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $z\sim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

We compare the star formation (SF) activity in cluster galaxies to the field from z=0.3-1.5 using $Herschel$ SPIRE 250$\mu$m imaging. We utilize 274 clusters from the IRAC Shallow Cluster Survey (ISCS) selected as rest-frame near-infrared overdensities over the 9 square degree Bootes field . This analysis allows us to quantify the evolution of SF in clusters over a long redshift baseline without bias against active cluster systems. Using a stacking analysis, we determine the average star formation rates (SFRs) and specific-SFRs (SSFR=SFR/M$_{\star}$) of stellar mass-limited (M>1.3x10$^{10}$ M$_{\odot}$), statistical samples of cluster and field galaxies, probing both the star forming and quiescent populations. We find a clear indication that the average SF in cluster galaxies is evolving more rapidly than in the field, with field SF levels at z>1.2 in the cluster cores (r<0.5 Mpc), in good agreement with previous ISCS studies. By quantifying the SF in cluster and field galaxies as an exponential function of cosmic time, we determine that cluster galaxies are evolving ~2 times faster than the field. Additionally, we see enhanced SF above the field level at z~1.4 in the cluster outskirts (r>0.5 Mpc). These general trends in the cluster cores and outskirts are driven by the lower mass galaxies in our sample. Blue cluster galaxies have systematically lower SSFRs than blue field galaxies, but otherwise show no strong differential evolution with respect to the field over our redshift range. This suggests that the cluster environment is both suppressing the star formation in blue galaxies on long time-scales and rapidly transitioning some fraction of blue galaxies to the quiescent galaxy population on short time-scales. We argue that our results are consistent with both strangulation and ram pressure stripping acting in these clusters, with merger activity occurring in the cluster outskirts.

We present Hubble Space Telescope near-IR spectroscopy for 18 galaxy clusters at 1.0 < $z$ < 1.5 in the IRAC Shallow Cluster Survey. We use Wide Field Camera 3 grism data to spectroscopically identify H$\alpha$ emitters in both the cores of galaxy clusters as well as in field galaxies. We find a large cluster-to-cluster scatter in the star formation rates within a projected radius of 500 kpc, and many of our clusters (~60%) have significant levels of star formation within a projected radius of 200 kpc. A stacking analysis reveals that dust reddening in these star-forming galaxies is positively correlated with stellar mass and may be higher in the field than the cluster at a fixed stellar mass. This may indicate a lower amount of gas in star-forming cluster galaxies than in the field population. Also, H$\alpha$ equivalent widths of star-forming galaxies in the cluster environment are still suppressed below the level of the field. This suppression is most significant for lower mass galaxies (log M$_{*}$ < 10.0 M$_{\odot}$). We therefore conclude that environmental effects are still important at 1.0 < $z$ < 1.5 for star-forming galaxies in galaxy clusters with log M$_{*}$ < 10.0 M$_{\odot}$.

The Spitzer-South Pole Telescope Deep Field (SSDF) is a wide-area survey using Spitzer's Infrared Array Camera (IRAC) to cover 94 square degrees of extragalactic sky, making it the largest IRAC survey completed to date outside the Milky Way midplane. The SSDF is centered at 23:30,-55:00, in a region that combines observations spanning a broad wavelength range from numerous facilities. These include millimeter imaging from the South Pole Telescope, far-infrared observations from Herschel/SPIRE, X-ray observations from the XMM XXL survey, near-infrared observations from the VISTA Hemisphere Survey, and radio-wavelength imaging from the Australia Telescope Compact Array, in a panchromatic project designed to address major outstanding questions surrounding galaxy clusters and the baryon budget. Here we describe the Spitzer/IRAC observations of the SSDF, including the survey design, observations, processing, source extraction, and publicly available data products. In particular, we present two band-merged catalogs, one for each of the two warm IRAC selection bands. They contain roughly 5.5 and 3.7 million distinct sources, the vast majority of which are galaxies, down to the SSDF 5-sigma sensitivity limits of 19.0 and 18.2 Vega mag (7.0 and 9.4 microJy) at 3.6 and 4.5 microns, respectively.

We measure the cross-power spectrum of the projected mass density as traced by the convergence of the cosmic microwave background lensing field from the South Pole Telescope (SPT) and a sample of Type 1 and 2 (unobscured and obscured) quasars at z~1 selected with the Wide-field Infrared Survey Explorer, over 2500 deg^2. The cross-power spectrum is detected at ~7-sigma, and we measure a linear bias b=1.67+/-0.24, consistent with clustering analyses. Using an independent lensing map, derived from Planck observations, to measure the cross-spectrum, we find excellent agreement with the SPT analysis. The bias of the combined sample of Type 1 and 2 quasars determined in this work is similar to that previously determined for Type 1 quasars alone; we conclude that that obscured and unobscured quasars must trace the matter field in a similar way. This result has implications for our understanding of quasar unification and evolution schemes.

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 observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

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

We present a measurement of the average supermassive black hole accretion rate (BHAR) as a function of star formation rate (SFR) for galaxies in the redshift range 0.25<z<0.8. We study a sample of 1,767 far-IR selected star-forming galaxies in the 9 deg^2 Bo\"otes multiwavelength survey field. The SFR is estimated using 250 micron observations from the Herschel Space Observatory, for which the contribution from the AGN is minimal. In this sample, 121 AGNs are directly identified using X-ray or mid-IR selection criteria. We combined these detected AGNs and an X-ray stacking analysis for undetected sources to study the average BHAR for all of the star-forming galaxies in our sample. We find an almost linear relation between the average BHAR (in M_sun/year) and the SFR (in M_sun/year) for galaxies across a wide SFR range 0.85<log SFR<2.56 : log BHAR=(-3.72\pm0.52)+(1.05\pm0.33) log SFR. This global correlation between SFR and average BHAR is consistent with a simple picture in which SFR and AGN activity are tightly linked over galaxy evolution timescales.

Stern et al.(2012) presented a study of WISE selection of AGN in the 2 deg^2 COSMOS field, finding that a simple criterion W1-W2>=0.8 provides a highly reliable and complete AGN sample for W2<15.05, where the W1 and W2 passbands are centered at 3.4 and 4.6 microns, respectively. Here we extend this study using the larger 9 deg^2 NOAO Deep Wide-Field Survey Bootes field which also has considerably deeper WISE observations than the COSMOS field, and find that this simple color-cut significantly loses reliability at fainter fluxes. We define a modified selection criterion combining the W1-W2 color and the W2 magnitude to provide highly reliable or highly complete AGN samples for fainter WISE sources. In particular, we define a color-magnitude cut that finds 130+/-4 deg^-2 AGN candidates for W2<17.11 with 90% reliability. Using the extensive UV through mid-IR broad-band photometry available in this field, we study the spectral energy distributions of WISE AGN candidates. As expected, the WISE AGN selection is biased towards objects where the AGN dominates the bolometric luminosity output, and that it can identify highly obscured AGN. We study the distribution of reddening in the AGN sample and discuss a formalism to account for sample incompleteness based on the step-wise maximum-likelihood method of Efstathiou et al.(1988). The resulting dust obscuration distributions depend strongly on AGN luminosity, consistent with the trend expected for a Simpson (2005) receding torus. At L_AGN~3x10^44 erg/s, 29+/-7% of AGN are observed as Type 1, while at ~4x10^45 erg/s the fraction is 64+/-13%. The distribution of obscuration values suggests that dust in the torus is present as both a diffuse medium and in optically thick clouds.

The details of the stellar mass assembly of brightest cluster galaxies (BCGs) remain an unresolved problem in galaxy formation. We have developed a novel approach that allows us to construct a sample of clusters that form an evolutionary sequence, and have applied it to the Spitzer IRAC Shallow Cluster Survey (ISCS) to examine the evolution of BCGs in progenitors of present-day clusters with mass of (2.5-4.5)x10^{14}Msun. We follow the cluster mass growth history extracted from a high resolution cosmological simulation, and then use an empirical method that infers the cluster mass based on the ranking of cluster luminosity to select high-z clusters of appropriate mass from ISCS to be progenitors of the given set of z=0 clusters. We find that, between z=1.5 and 0.5, the BCGs have grown in stellar mass by a factor of 2.3, which is well-matched by the predictions from a state-of-the-art semi-analytic model. Below z=0.5 we see hints of differences in behavior between the model and observation.

We report the detection of CO(1-0) line emission in the bright, lensed star-forming galaxies SPT-S 233227-5358.5 (z=2.73) and SPT-S 053816-5030.8 (z=2.78), using the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in a large-area millimeter survey with the South Pole Telescope (SPT) and found to be gravitationally lensed by intervening structures. The measured CO intensities imply galaxies with molecular gas masses of (3.2 \pm 0.5)x10^10 (mu/15)^{-1}(X_CO/0.8) M_sun and (1.7 \pm 0.3)x10^10 (mu/20)^{-1}(X_CO/0.8) M_sun, and gas depletion timescales of 4.9x10^7 (X_CO/0.8) yr and 2.6x10^7 (X_CO/0.8) yr, respectively, where mu corresponds to the lens magnification and X_CO is the CO luminosity to gas mass conversion factor. In the case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy distribution and an assumed synchrotron spectral index, we find that 42 \pm 10 % and 55 \pm 13 % of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. We find a radio-derived intrinsic star formation rate (SFR) of 470 \pm 170 M_sun/yr, consistent within the uncertainties with the infrared estimate. Based on the morphology of this object in the source plane, the derived gas mass and the possible flattening of the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy.

We report the first results from the Clusters Around Radio-Loud AGN (CARLA) program, a Cycle 7 and 8 Spitzer Space Telescope snapshot program to investigate the environments of a large sample of obscured and unobscured luminous radio-loud AGN at 1.2 < z < 3.2. These data, obtained for 387 fields, reach 3.6 and 4.5 um depths of [3.6] (AB) = 22.6 and [4.5] (AB) = 22.9 at the 95% completeness level, which is two to three times fainter than L* in this redshift range. By using the color cut [3.6]-[4.5] > -0.1 (AB), which efficiently selects high-redshift (z > 1.3) galaxies of all types, we identify galaxy cluster member candidates in the fields of the radio-loud AGN. The local density of these IRAC-selected sources is compared to the density of similarly selected sources in blank fields. We find that 92% of the radio-loud AGN reside in environments richer than average. The majority (55%) of the radio-loud AGN fields are found to be overdense at a > 2 {\sigma} level; 10% are overdense at a > 5 {\sigma} level. A clear rise in surface density of IRAC-selected sources towards the position of the radio-loud AGN strongly supports an association of the majority of the IRAC-selected sources with the radio-loud AGN. Our results provide solid statistical evidence that radio-loud AGN are likely beacons for finding high-redshift galaxy (proto-)clusters. We investigate how environment depends on AGN type (unobscured radio-loud quasars vs. obscured radio galaxies), radio luminosity and redshift, finding no correlation with either AGN type or radio luminosity. We find a decrease in density with redshift, consistent with galaxy evolution for this uniform, flux-limited survey. These results are consistent with expectations from the orientation-driven AGN unification model, at least for the high radio luminosity regimes considered in this sample.

We use a temperature map of the cosmic microwave background (CMB) obtained using the South Pole Telescope at 150 GHz to construct a map of the gravitational convergence to z ~ 1100, revealing the fluctuations in the projected mass density. This map shows individual features that are significant at the ~ 4 sigma level, providing the first image of CMB lensing convergence. We cross-correlate this map with Herschel/SPIRE maps covering 90 square degrees at wavelengths of 500, 350, and 250 microns. We show that these submillimeter-wavelength (submm) maps are strongly correlated with the lensing convergence map, with detection significances in each of the three submm bands ranging from 6.7 to 8.8 sigma. We fit the measurement of the cross power spectrum assuming a simple constant bias model and infer bias factors of b=1.3-1.8, with a statistical uncertainty of 15%, depending on the assumed model for the redshift distribution of the dusty galaxies that are contributing to the Herschel/SPIRE maps.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted a blind redshift survey in the 3 mm atmospheric transmission window for 26 strongly lensd dusty star-forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected to have S_1.4mm>20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio (S_843MHz<6mJy) or far-infrared counterparts (S_100um<1 Jy, S_60um<200mJy). We robustly detect 44 line features in our survey, which we identify as redshifted emission lines of 12CO, 13CO, [CI], H2O, and H2O+. We find one or more spectral features in 23 sources yielding a ~90% detection rate for this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single line. For the sources with only one detected line, we break the redshift degeneracy with additional spectroscopic observations if available, or infer the most likely line identification based on photometric data. This yields secure redshifts for ~70% of the sample. The three sources with no lines detected are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our sample is <z>=3.5. This finding is in contrast to the redshift distribution of radio-identified DSFGs, which have a significantly lower mean redshift of <z>=2.3 and for which only 10-15% of the population is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution and compare our measured redshift distribution to that of models in the literature.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

The fraction of cluster galaxies that host luminous AGN is an important probe of AGN fueling processes, the cold ISM at the centers of galaxies, and how tightly black holes and galaxies co-evolve. We present a new measurement of the AGN fraction in a sample of 13 clusters of galaxies (M >= 10^{14} Msun) at 1<z<1.5 selected from the Spitzer/IRAC Shallow Cluster Survey, as well as the field fraction in the immediate vicinity of these clusters, and combine these data with measurements from the literature to quantify the relative evolution of cluster and field AGN from the present to z~3. We estimate that the cluster AGN fraction at 1<z<1.5 is f_A = 3.0^{+2.4}_{-1.4}% for AGN with a rest-frame, hard X-ray luminosity greater than L_{X,H} >= 10^{44} erg/s. This fraction is measured relative to all cluster galaxies more luminous than M*_{3.6}(z)+1, where M*_{3.6}(z) is the absolute magnitude of the break in the galaxy luminosity function at the cluster redshift in the IRAC 3.6um bandpass. The cluster AGN fraction is 30 times greater than the 3sigma upper limit on the value for AGN of similar luminosity at z~0.25, as well as more than an order of magnitude greater than the AGN fraction at z~0.75. AGN with L_{X,H} >= 10^{43} erg/s exhibit similarly pronounced evolution with redshift. In contrast with the local universe, where the luminous AGN fraction is higher in the field than in clusters, the X-ray and MIR-selected AGN fractions in the field and clusters are consistent at 1<z<1.5. This is evidence that the cluster AGN population has evolved more rapidly than the field population from z~1.5 to the present. This environment-dependent AGN evolution mimics the more rapid evolution of star-forming galaxies in clusters relative to the field.

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

We measure the clustering of Extremely Red Objects (EROs) in ~8 deg^2 of the NOAO Deep Wide Field Survey Bo\"otes field in order to establish robust links between ERO z~1.2 and local galaxy z<0.1 populations. Three different color selection criteria from the literature are analyzed to assess the consequences of using different criteria for selecting EROs. Specifically, our samples are (R-K_s)>5.0 (28,724 galaxies), (I-K_s)>4.0 (22,451 galaxies) and (I-[3.6])>5.0 (64,370 galaxies). Magnitude-limited samples show the correlation length (r_0) to increase for more luminous EROs, implying a correlation with stellar mass. We can separate star-forming and passive ERO populations using the (K_s-[24]) and ([3.6]-[24]) colors to K_s=18.4 and [3.6]=17.5, respectively. Star-forming and passive EROs in magnitude limited samples have different clustering properties and host dark halo masses, and cannot be simply understood as a single population. Based on the clustering, we find that bright passive EROs are the likely progenitors of >4L^* elliptical galaxies. Bright EROs with ongoing star formation were found to occupy denser environments than star-forming galaxies in the local Universe, making these the likely progenitors of >L^* local ellipticals. This suggests that the progenitors of massive >4L^* local ellipticals had stopped forming stars by z>1.2, but that the progenitors of less massive ellipticals (down to L^*) can still show significant star formation at this epoch.

We present the results of the ground- and space-based optical and near-infrared (NIR) follow-up of 224 galaxy cluster candidates detected with the Sunyaev-Zel'dovich (SZ) effect in the 720 deg^2 of the South Pole Telescope (SPT) survey completed in the 2008 and 2009 observing seasons. We use the optical/NIR data to establish whether each candidate is associated with an overdensity of galaxies and to estimate the cluster redshift. Most photometric redshifts are derived through a combination of three different cluster redshift estimators using red-sequence galaxies, resulting in an accuracy of \Delta z/(1+z)=0.017, determined through comparison with a subsample of 57 clusters for which we have spectroscopic redshifts. We successfully measure redshifts for 158 systems and present redshift lower limits for the remaining candidates. The redshift distribution of the confirmed clusters extends to z=1.35 with a median of z_{med}=0.57. Approximately 18% of the sample with measured redshifts lies at z>0.8. We estimate a lower limit to the purity of this SPT SZ-selected sample by assuming that all unconfirmed clusters are noise fluctuations in the SPT data. We show that the cumulative purity at detection significance \xi>5 (\xi>4.5) is >= 95 (>= 70%). We present the red brightest cluster galaxy (rBCG) positions for the sample and examine the offsets between the SPT candidate position and the rBCG. The radial distribution of offsets is similar to that seen in X-ray-selected cluster samples, providing no evidence that SZ-selected cluster samples include a different fraction of recent mergers than X-ray-selected cluster samples.

We compare cosmic microwave background lensing convergence maps derived from South Pole Telescope (SPT) data with galaxy survey data from the Blanco Cosmology Survey, the Wide-field Infrared Survey Explorer, and a new large Spitzer/IRAC field designed to overlap with the SPT survey. Using optical and infrared catalogs covering between 17 and 68 square degrees of sky, we detect correlation between the SPT convergence maps and each of the galaxy density maps at >4 sigma, with zero cross-correlation robustly ruled out in all cases. The amplitude and shape of the cross-power spectra are in good agreement with theoretical expectations and the measured galaxy bias is consistent with previous work. The detections reported here utilize a small fraction of the full 2500 square degree SPT survey data and serve as both a proof of principle of the technique and an illustration of the potential of this emerging cosmological probe.

We measure the faint end slope of the galaxy luminosity function (LF) for cluster galaxies at 1<z<1.5 using Spitzer IRAC data. We investigate whether this slope, alpha, differs from that of the field LF at these redshifts, and with the cluster LF at low redshifts. The latter is of particular interest as low-luminosity galaxies are expected to undergo significant evolution. We use seven high-redshift spectroscopically confirmed galaxy clusters drawn from the IRAC Shallow Cluster Survey to measure the cluster galaxy LF down to depths of M* + 3 (3.6 microns) and M* + 2.5 (4.5 microns). The summed LF at our median cluster redshift (z=1.35) is well fit by a Schechter distribution with alpha[3.6] = -0.97 +/- 0.14 and alpha[4.5] = -0.91 +/- 0.28, consistent with a flat faint end slope and is in agreement with measurements of the field LF in similar bands at these redshifts. A comparison to alpha in low-redshift clusters finds no statistically significant evidence of evolution. Combined with past studies which show that M* is passively evolving out to z~1.3, this means that the shape of the cluster LF is largely in place by z~1.3. This suggests that the processes that govern the build up of the mass of low-mass cluster galaxies have no net effect on the faint end slope of the cluster LF at z<1.3.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

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

We present spectroscopic confirmation of a z=0.99 galaxy cluster discovered using data from the Wide-field Infrared Survey Explorer (WISE). This is the first z~1 cluster candidate from the Massive Distant Clusters of WISE Survey (MaDCoWS) to be confirmed. It was selected as an overdensity of probable z>~1 sources using a combination of WISE and SDSS-DR8 photometric catalogs. Deeper follow-up imaging data from Subaru and WIYN reveal the cluster to be a rich system of galaxies, and multi-object spectroscopic observations from Keck confirm five cluster members at z=0.99. The detection and confirmation of this cluster represents a first step towards constructing a uniformly-selected sample of distant, high-mass galaxy clusters over the full extragalactic sky using WISE data.

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

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

We present results for the assembly and star formation histories of massive (~L*) red sequence galaxies in 11 spectroscopically confirmed, infrared-selected galaxy clusters at 1.0 < z < 1.5, the precursors to present-day massive clusters with M ~ 10^15 M_sun. Using rest-frame optical photometry, we investigate evolution in the color and scatter of the red sequence galaxy population, comparing with models of possible star formation histories. In contrast to studies of central cluster galaxies at lower redshift (z < 1), these data are clearly inconsistent with the continued evolution of stars formed and assembled primarily at a single, much-earlier time. Specifically, we find that the colors of massive cluster galaxies at z = 1.5 imply that the bulk of star formation occurred at z ~ 3, whereas by z = 1 their colors imply formation at z ~ 2; therefore these galaxies exhibit approximately the same luminosity-weighted stellar age at 1 < z < 1.5. This likely reflects star formation that occurs over an extended period, the effects of significant progenitor bias, or both. Our results generally indicate that massive cluster galaxy populations began forming a significant mass of stars at z >~ 4, contained some red spheroids by z ~ 1.5, and were actively assembling much of their final mass during 1 < z < 2 in the form of younger stars. Qualitatively, the slopes of the cluster color-magnitude relations are consistent with no significant evolution relative to local clusters.

We report the discovery of an IR-selected galaxy cluster in the IRAC Distant Cluster Survey (IDCS). New data from the Hubble Space Telescope spectroscopically confirm IDCS J1433.2+3306 at z = 1.89 with robust spectroscopic redshifts for seven members, two of which are based on the 4000 Angstrom break. Detected emission lines such as [OII] and Hbeta indicate star formation rates of >20 solar masses per year for three galaxies within a 500 kpc projected radius of the cluster center. The cluster exhibits a red sequence with a scatter and color indicative of a formation redshift z > 3.5. The stellar age of the early-type galaxy population is approximately consistent with those of clusters at lower redshift (1 < z < 1.5) suggesting that clusters at these redshifts are experiencing ongoing or increasing star formation.

We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.

We report the discovery of an IR-selected massive galaxy cluster in the IRAC Distant Cluster Survey (IDCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm IDCS J1426+3508 at z=1.75. Moreover, the cluster is detected in archival Chandra data as an extended X-ray source, comprising 54 counts after the removal of point sources. We calculate an X-ray luminosity of L{0.5-2 keV} = (5.5 +/- 1.2) X 1e44 ergs/s within r = 60 arcsec (~1 Mpc diameter), which implies M_{200,L_x} = (5.6 +/- 1.6) X 1e14 Msun. IDCS J1426+3508 appears to be an exceptionally massive cluster for its redshift.

The galaxy cluster IDCS J1426.5+3508 at z = 1.75 is the most massive galaxy cluster yet discovered at z > 1.4 and the first cluster at this epoch for which the Sunyaev-Zel'Dovich effect has been observed. In this paper we report on the discovery with HST imaging of a giant arc associated with this cluster. The curvature of the arc suggests that the lensing mass is nearly coincident with the brightest cluster galaxy, and the color is consistent with the arc being a star-forming galaxy. We compare the constraint on M200 based upon strong lensing with Sunyaev-Zel'Dovich results, finding that the two are consistent if the redshift of the arc is z > 3. Finally, we explore the cosmological implications of this system, considering the likelihood of the existence of a strongly lensing galaxy cluster at this epoch in an LCDM universe. While the existence of the cluster itself can potentially be accomodated if one considers the entire volume covered at this redshift by all current high-redshift cluster surveys, the existence of this strongly lensed galaxy greatly exacerbates the long-standing giant arc problem. For standard LCDM structure formation and observed background field galaxy counts this lens system should not exist. Specifically, there should be no giant arcs in the entire sky as bright in F814W as the observed arc for clusters at z \geq 1.75, and only \sim 0.3 as bright in F160W as the observed arc. If we relax the redshift constraint to consider all clusters at z \geq 1.5, the expected number of giant arcs rises to \sim15 in F160W, but the number of giant arcs of this brightness in F814W remains zero. These arc statistic results are independent of the mass of IDCS J1426.5+3508. We consider possible explanations for this discrepancy.

We report 31 GHz CARMA observations of IDCS J1426.5+3508, an infrared-selected galaxy cluster at z = 1.75. A Sunyaev-Zel'dovich decrement is detected towards this cluster, indicating a total mass of M200 = (4.3 +/- 1.1) x 10^{14} Msun in agreement with the approximate X-ray mass of ~5 x 10^{14} Msun. IDCS J1426.5+3508 is by far the most distant cluster yet detected via the Sunyaev-Zel'dovich effect, and the most massive z >= 1.4 galaxy cluster found to date. Despite the mere ~1% probability of finding it in the 8.82 deg^2 IRAC Distant Cluster Survey, IDCS J1426.5+3508 is not completely unexpected in LCDM once the area of large, existing surveys is considered. IDCS J1426.5+3508 is, however, among the rarest, most extreme clusters ever discovered, and indeed is an evolutionary precursor to the most massive known clusters at all redshifts. We discuss how imminent, highly sensitive Sunyaev-Zel'dovich experiments will complement infrared techniques for statistical studies of the formation of the most massive galaxy clusters in the z > 1.5 Universe, including potential precursors to IDCS J1426.5+3508.

We present observations of CO J=2-1 line emission in infrared-luminous cluster galaxies at z~1 using the IRAM Plateau de Bure Interferometer. Our two primary targets are optically faint, dust-obscured galaxies (DOGs) found to lie within 2 Mpc of the centers of two massive (>10^14 Msun) galaxy clusters. CO line emission is not detected in either DOG. We calculate 3-sigma upper limits to the CO J=2-1 line luminosities, L'_CO < 6.08x10^9 and < 6.63x10^9 K km/s pc^2. Assuming a CO-to-H_2 conversion factor derived for ultraluminous infrared galaxies in the local Universe, this translates to limits on the cold molecular gas mass of M_H_2 < 4.86x10^9 Msun and M_H_2 < 5.30x10^9 Msun. Both DOGs exhibit mid-infrared continuum emission that follows a power-law, suggesting that an AGN contributes to the dust heating. As such, estimates of the star formation efficiencies in these DOGs are uncertain. A third cluster member with an infrared luminosity, L_IR < 7.4x10^11 Lsun, is serendipitously detected in CO J=2-1 line emission in the field of one of the DOGs located roughly two virial radii away from the cluster center. The optical spectrum of this object suggests that it is likely an obscured AGN, and the measured CO line luminosity is L'_CO = (1.94 +/- 0.35)x10^10 K km/s pc^2, which leads to an estimated cold molecular gas mass M_H_2 = (1.55+/-0.28)x10^10 Msun. A significant reservoir of molecular gas in a z~1 galaxy located away from the cluster center demonstrates that the fuel can exist to drive an increase in star-formation and AGN activity at the outskirts of high-redshift clusters.

Detailed analysis of the substructure of Lya nebulae can put important constraints on the physical mechanisms at work and the properties of galaxies forming within them. Using high resolution HST imaging of a Lya nebula at z~2.656, we have taken a census of the compact galaxies in the vicinity, used optical/near-infrared colors to select system members, and put constraints on the morphology of the spatially-extended emission. The system is characterized by (a) a population of compact, low luminosity (~0.1 L*) sources --- 17 primarily young, small (Re~1-2 kpc), disky galaxies including an obscured AGN --- that are all substantially offset (>20 kpc) from the line-emitting nebula; (b) the lack of a central galaxy at or near the peak of the Lya emission; and (c) several nearly coincident, spatially extended emission components --- Lya, HeII, and UV continuum --- that are extremely smooth. These morphological findings are difficult to reconcile with theoretical models that invoke outflows, cold flows, or resonant scattering, suggesting that while all of these physical phenomena may be occurring, they are not sufficient to explain the powering and large extent of Lya nebulae. In addition, although the compact galaxies within the system are irrelevant as power sources, the region is significantly overdense relative to the field galaxy population (by at least a factor of 4). These observations provide the first estimate of the luminosity function of galaxies within an individual Lya nebula system, and suggest that large Lya nebulae may be the seeds of galaxy groups or low-mass clusters.

We present a panoramic narrow-band imaging survey of [OII] emitters in and around the ClG J0218.3-0510 cluster at z=1.62 with Suprime-Cam on Subaru telescope. 352 [OII] emitters were identified on the basis of narrow-band excesses and photometric redshifts. We discovered a huge filamentary structure with some clumps traced by [OII] emitters and found that the ClG J0218.3-0510 cluster is embedded in an even larger super-structure than the one reported previously. 31 [OII] emitters were spectroscopically confirmed with the detection of H-alpha and/or [OIII] emission lines by FMOS observations. In the high density regions such as cluster core and clumps, star-forming [OII] emitters show a high overdensity by a factor of more than 10 compared to the field region. Although the star formation activity is very high even in the cluster core, some massive quiescent galaxies also exits at the same time. Furthermore, the properties of the individual [OII] emitters, such as star formation rates, stellar masses and specific star formation rates, do not show a significant dependence on the local density, either. Such lack of environmental dependence is consistent with our earlier result by Hayashi et al. (2011) on a z=1.5 cluster and its surrounding region. The fact that the star-forming activity of galaxies in the cluster core is as high as that in the field at z~1.6 may suggest that the star-forming galaxies are probably just in a transition phase from a starburst mode to a quiescent mode, and are thus showing comparable level of star formation rates to those in lower density environments. We may be witnessing the start of the reversal of the local SFR--density relation due to the "biased" galaxy formation and evolution in high density regions at high this redshift, beyond which massive galaxies would be forming vigorously in a more biased way in proto-cluster cores.

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

The South Pole Telescope (SPT) is a 10 meter telescope operating at mm wavelengths. It has recently completed a three-band survey covering 2500 sq. degrees. One of the survey's main goals is to detect galaxy clusters using Sunyaev-Zeldovich effect and use these clusters for a variety of cosmological and astrophysical studies such as the dark energy equation of state, the primordial non-gaussianity and the evolution of galaxy populations. Since 2005, we have been engaged in a comprehensive optical and near-infrared followup program (at wavelengths between 0.4 and 5 {\mu}m) to image high-significance SPT clusters, to measure their photometric redshifts, and to estimate the contamination rate of the candidate lists. These clusters are then used for various cosmological and astrophysical studies.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

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

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

We report on the serendipitous discovery in the Blanco Cosmology Survey (BCS) imaging data of a z = 0.9057 galaxy that is being strongly lensed by a massive galaxy cluster at a redshift of z=0.3838. The lens (BCS J2352-5452) was discovered while examining i- and z-band images being acquired in October 2006 during a BCS observing run. Follow-up spectroscopic observations with the GMOS instrument on the Gemini South 8m telescope confirmed the lensing nature of this system. Using weak plus strong lensing, velocity dispersion, cluster richness N_200, and fitting to an NFW cluster mass density profile, we have made three independent estimates of the mass M_200 which are all very consistent with each other. The combination of the results from the three methods gives M_200 = (5.1 +- 1.3) x 10^14 M_Sun, which is fully consistent with the individual measurements. The final NFW concentration c_200 from the combined fit is c_200 = 5.4(+1.4 -1.1). We have compared our measurements of M_200 and c_200 with predictions for (a) clusters from Lambda_CDM simulations, (b) lensing selected clusters from simulations, and (c) a real sample of cluster lenses. We find that we are most compatible with the predictions for Lambda_CDM simulations for lensing clusters, and we see no evidence based on this one system for an increased concentration compared to Lambda_CDM. Finally, using the flux measured from the [OII]3727 line we have determined the star formation rate (SFR) of the source galaxy and find it to be rather modest given the assumed lens magnification.

We present results from a systematic study of star formation in local galaxy clusters using 22 micron data from the Wide-field Infrared Survey Explorer (WISE). The 69 systems in our sample are drawn from the Cluster Infall Regions Survey (CIRS), and all have robust mass determinations. The all-sky WISE data enables us to quantify the amount of star formation, as traced by 22 micron, as a function of radius well beyond R200, and investigate the dependence of total star formation rate upon cluster mass. We find that the fraction of star-forming galaxies increases with cluster radius but remains below the field value even at 3 R200. We also find that there is no strong correlation between the mass-normalized total specific star formation rate and cluster mass, indicating that the mass of the host cluster does not strongly influence the total star formation rate of cluster members.

We present South Pole Telescope (SPT) observations of the five galaxy cluster candidates in the southern hemisphere which were reported as unconfirmed in the Planck Early Sunyaev-Zel'dovich (ESZ) sample. One cluster candidate, PLCKESZ G255.62-46.16, is located in the 2500-square-degree SPT SZ survey region and was reported previously as SPT-CL J0411-4819. For the remaining four candidates, which are located outside of the SPT SZ survey region, we performed short, dedicated SPT observations. Each of these four candidates was strongly detected in maps made from these observations, with signal-to-noise ratios ranging from 6.3 to 13.8. We have observed these four candidates on the Magellan-Baade telescope and used these data to estimate cluster redshifts from the red sequence. Resulting redshifts range from 0.24 to 0.46. We report measurements of Y_0.75', the integrated Comptonization within a 0.75' radius, for all five candidates. We also report X-ray luminosities calculated from ROSAT All-Sky Survey catalog counts, as well as optical and improved SZ coordinates for each candidate. The combination of SPT SZ measurements, optical red-sequence measurements, and X-ray luminosity estimates demonstrates that these five Planck ESZ cluster candidates do indeed correspond to real galaxy clusters with redshifts and observable properties consistent with the rest of the ESZ sample.

The South Pole Telescope (SPT) is currently surveying 2500 deg^2 of the southern sky to detect massive galaxy clusters out to the epoch of their formation using the Sunyaev-Zel'dovich (SZ) effect. This paper presents a catalog of the 26 most significant SZ cluster detections in the full survey region. The catalog includes 14 clusters which have been previously identified and 12 that are new discoveries. These clusters were identified in fields observed to two differing noise depths: 1500 deg^2 at the final SPT survey depth of 18 uK-arcmin at 150 GHz, and 1000 deg^2 at a depth of 54 uK-arcmin. Clusters were selected on the basis of their SZ signal-to-noise ratio (S/N) in SPT maps, a quantity which has been demonstrated to correlate tightly with cluster mass. The S/N thresholds were chosen to achieve a comparable mass selection across survey fields of both depths. Cluster redshifts were obtained with optical and infrared imaging and spectroscopy from a variety of ground- and space-based facilities. The redshifts range from 0.098 \leq z \leq 1.132 with a median of z_med = 0.40. The measured SZ S/N and redshifts lead to unbiased mass estimates ranging from 9.8 \times 10^14 M_sun/h_70 \leq M_200(rho_mean) \leq 3.1 \times 10^15 M_sun/h_70. Based on the SZ mass estimates, we find that none of the clusters are individually in significant tension with the LambdaCDM cosmological model. We also test for evidence of non-Gaussianity based on the cluster sample and find the data show no preference for non-Gaussian perturbations.

We present weak gravitational lensing analysis of 22 high-redshift (z >~1) clusters based on Hubble Space Telescope images. Most clusters in our sample provide significant lensing signals and are well detected in their reconstructed two-dimensional mass maps. Combining the current results and our previous weak-lensing studies of five other high-z clusters, we compare gravitational lensing masses of these clusters with other observables. We revisit the question whether the presence of the most massive clusters in our sample is in tension with the current LambdaCDM structure formation paradigm. We find that the lensing masses are tightly correlated with the gas temperatures and establish, for the first time, the lensing mass-temperature relation at z >~ 1. For the power law slope of the M-TX relation (M propto T^{\alpha}), we obtain \alpha=1.54 +/- 0.23. This is consistent with the theoretical self-similar prediction \alpha=3/2 and with the results previously reported in the literature for much lower redshift samples. However, our normalization is lower than the previous results by 20-30%, indicating that the normalization in the M-TX relation might evolve. After correcting for Eddington bias and updating the discovery area with a more conservative choice, we find that the existence of the most massive clusters in our sample still provides a tension with the current Lambda CDM model. The combined probability of finding the four most massive clusters in this sample after marginalization over current cosmological parameters is less than 1%.

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

[Abridged] We investigate the average physical properties and star formation histories of the most UV-luminous star-forming galaxies at z~3.7. Our results are derived from analyses of the average spectral energy distributions (SEDs), constructed from stacked optical to infrared photometry, of a sample of the 1,902 most UV-luminous star-forming galaxies found in 5.3 square degrees of the NOAO Deep Wide-Field Survey. We bin the sample according to UV luminosity, and find that the shape of the average SED in the rest-frame optical and infrared is fairly constant with UV luminosity: i.e., more UV luminous galaxies are, on average, also more luminous at longer wavelengths. In the rest-UV, however, the spectral slope (measured at 0.13-0.28 um) rises steeply with the median UV luminosity from -1.8 at L L* to -1.2 in the brightest bin (L~4-5L*). We use population synthesis analyses to derive the average physical properties of these galaxies and find that: (1) L_UV, and thus star formation rates (SFRs), scale closely with stellar mass such that more UV-luminous galaxies are also more massive; (2) The median ages indicate that the stellar populations are relatively young (200-400 Myr) and show little correlation with UV luminosity; and (3) More UV-luminous galaxies are dustier than their less-luminous counterparts, such that L~4-5L* galaxies are extincted up to A(1600)=2 mag while L L* galaxies have A(1600)=0.7-1.5 mag. Based on these observations, we argue that the average star formation histories of UV-luminous galaxies are better described by models in which SFR increases with time in order to simultaneously reproduce the tight correlation between the observed SFR and stellar mass, and the universally young ages of these galaxies. We demonstrate the potential of measurements of the SFR-M* relation at multiple redshifts to discriminate between simple models of star formation histories.

We present first results of an examination of the optical properties of the galaxy populations in SZE selected galaxy clusters. Using clusters selected by the South Pole Telescope survey and deep multiband optical data from the Blanco Cosmology Survey, we measure the radial profile, the luminosity function, the blue fraction and the halo occupation number of the galaxy populations of these four clusters with redshifts ranging from 0.3 to 1. Our goal is to understand whether there are differences among the galaxy populations of these SZE selected clusters and previously studied clusters selected in the optical and the X-ray. The radial distributions of galaxies in the four systems are consistent with NFW profiles with a galaxy concentration of 3 to 6. We show that the characteristic luminosities in $griz$ bands are consistent with passively evolving populations emerging from a single burst at redshift $z=3$. The faint end power law slope of the luminosity function is found to be on average $\alpha \approx -1.2$ in griz. Halo occupation numbers (to $m^*+2$) for these systems appear to be consistent with those based on X-ray selected clusters. The blue fraction estimated to $0.36L^*$, for the three lower redshift systems, suggests an increase with redshift, although with the current sample the uncertainties are still large. Overall, this pilot study of the first four clusters provides no evidence that the galaxy populations in these systems differ significantly from those in previously studied cluster populations selected in the X-ray or the optical.

We report the results of an 87 square-degree point-source survey centered at R.A. 5h30m, decl. -55 deg. taken with the South Pole Telescope (SPT) at 1.4 and 2.0 mm wavelengths with arc-minute resolution and milli-Jansky depth. Based on the ratio of flux in the two bands, we separate the detected sources into two populations, one consistent with synchrotron emission from active galactic nuclei (AGN) and one consistent with thermal emission from dust. We present source counts for each population from 11 to 640 mJy at 1.4 mm and from 4.4 to 800 mJy at 2.0 mm. The 2.0 mm counts are dominated by synchrotron-dominated sources across our reported flux range; the 1.4 mm counts are dominated by synchroton-dominated sources above ~15 mJy and by dust-dominated sources below that flux level. We detect 141 synchrotron-dominated sources and 47 dust-dominated sources at S/N > 4.5 in at least one band. All of the most significantly detected members of the synchrotron-dominated population are associated with sources in previously published radio catalogs. Some of the dust-dominated sources are associated with nearby (z << 1) galaxies whose dust emission is also detected by the Infrared Astronomy Satellite (IRAS). However, most of the bright, dust-dominated sources have no counterparts in any existing catalogs. We argue that these sources represent the rarest and brightest members of the population commonly referred to as sub-millimeter galaxies (SMGs). Because these sources are selected at longer wavelengths than in typical SMG surveys, they are expected to have a higher mean redshift distribution and may provide a new window on galaxy formation in the early universe.

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured quasars, using a sample in the redshift range 0.7 < z < 1.8 selected from the 9 deg^2 Bootes multiwavelength survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of obscured quasars that have been inferred from models of the X-ray background and supermassive black hole evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this study we use a sample of 806 mid-infrared selected quasars and ~250,000 galaxies to calculate the projected quasar-galaxy cross-correlation function w_p(R). The observed clustering yields characteristic dark matter halo masses of log (M_halo [h^-1 M_sun]) = 12.7^+0.4_-0.6 and 13.3^+0.3_-0.4 for unobscured quasars (QSO-1s) and obscured quasars (Obs-QSOs), respectively. The results for QSO-1s are in excellent agreement with previous measurements for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the QSO-1s. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and speculate on physical implications of stronger clustering for obscured quasars.

We present the J-band luminosity function of 1838 mid-infrared and X-ray selected AGNs in the redshift range 0<z<5.85. These luminosity functions are constructed by combining the deep multi-wavelength broad-band observations from the UV to the mid-IR of the NDWFS Bootes field with the X-ray observations of the XBootes survey and the spectroscopic observations of the same field by AGES. Our sample is primarily composed of IRAC-selected AGNs, targeted using modifications of the Stern et al.(2005) criteria, complemented by MIPS 24 microns and X-ray selected AGNs to alleviate the biases of IRAC mid-IR selection against z~4.5 quasars and AGNs faint with respect to their hosts. This sample provides an accurate link between low and high redshift AGN luminosity functions and does not suffer from the usual incompleteness of optical samples at z~3. We find that the space density of the brightest quasars strongly decreases from z=3 to z=0, while the space density of faint quasars is at least flat, and possibly increasing, over the same redshift range. At z>3 we observe a decrease in the space density of quasars of all brightnesses. We model the luminosity function by a double power-law and find that its evolution cannot be described by either pure luminosity or pure density evolution, but must be a combination of both. Our best-fit model has bright and faint power-law indices consistent with the low redshift measurements based on the 2QZ and 2SLAQ surveys and it generally agrees with the number of bright quasars predicted by other LFs at all redshifts. If we construct the QSO luminosity function using only the IRAC-selected AGNs, we find that the biases inherent to this selection method significantly modify the behavior of phi*(z) only for z<1 and have no significant impact upon the characteristic magnitude M*_J(z).

Using the South Pole Telescope (SPT), we have discovered the most massive known galaxy cluster at z > 1, SPT-CL J2106-5844. In addition to producing a strong Sunyaev-Zel'dovich effect signal, this system is a luminous X-ray source and its numerous constituent galaxies display spatial and color clustering, all indicating the presence of a massive galaxy cluster. VLT and Magellan spectroscopy of 18 member galaxies shows that the cluster is at z = 1.132^+0.002_-0.003. Chandra observations obtained through a combined HRC-ACIS GTO program reveal an X-ray spectrum with an Fe K line redshifted by z = 1.18 +/- 0.03. These redshifts are consistent with galaxy colors in extensive optical, near-infrared, and mid-infrared imaging. SPT-CL J2106-5844 displays extreme X-ray properties for a cluster, having a core-excluded temperature of kT = 11.0^+2.6_-1.9 keV and a luminosity (within r_500) of L_X (0.5 - 2.0 keV) = (13.9 +/- 1.0) x 10^44 erg/s. The combined mass estimate from measurements of the Sunyaev-Zel'dovich effect and X-ray data is M_200 = (1.27 +/- 0.21) x 10^15 M_sun. The discovery of such a massive gravitationally collapsed system at high redshift provides an interesting laboratory for galaxy formation and evolution, and is a powerful probe of extreme perturbations of the primordial matter density field. We discuss the latter, determining that, under the assumption of LambdaCDM cosmology with only Gaussian perturbations, there is only a 7% chance of finding a galaxy cluster similar to SPT-CL J2106-5844 in the 2500 deg^2 SPT survey region, and that only one such galaxy cluster is expected in the entire sky.

We report the X-ray detection of two z>1.4 infrared-selected galaxy clusters from the IRAC Shallow Cluster Survey (ISCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm cluster ISCS J1432.4+3250 at z=1.49, the most distant of 18 confirmed z>1 clusters in the ISCS to date. We also present new spectroscopy for ISCS J1438.1+3414, previously reported at z = 1.41, and measure its dynamical mass. Clusters ISCS J1432.4+3250 and ISCS J1438.1+3414 are detected in 36ks and 143ks Chandra exposures at significances of 5.2 sigma and 9.7 sigma, from which we measure total masses of log(M_{200,Lx}/Msun) = 14.4 +/- 0.2 and 14.35^{+0.14}_{-0.11}, respectively. The consistency of the X-ray and dynamical properties of these high redshift clusters further demonstrates that the ISCS is robustly detecting massive clusters to at least z = 1.5.

We present a detection-significance-limited catalog of 21 Sunyaev-Zel'dovich selected galaxy clusters. These clusters, along with 1 unconfirmed candidate, were identified in 178 deg^2 of sky surveyed in 2008 by the South Pole Telescope to a depth of 18 uK-arcmin at 150 GHz. Optical imaging from the Blanco Cosmology Survey (BCS) and Magellan telescopes provided photometric (and in some cases spectroscopic) redshift estimates, with catalog redshifts ranging from z=0.15 to z>1, with a median z = 0.74. Of the 21 confirmed galaxy clusters, three were previously identified as Abell clusters, three were presented as SPT discoveries in Staniszewski et al, 2009, and three were first identified in a recent analysis of BCS data by Menanteau et al, 2010; the remaining 12 clusters are presented for the first time in this work. Simulated observations of the SPT fields predict the sample to be nearly 100% complete above a mass threshold of M_200 ~ 5x10^14 M_sun/h at z = 0.6. This completeness threshold pushes to lower mass with increasing redshift, dropping to ~4x10^14 M_sun/h at z=1. The size and redshift distribution of this catalog are in good agreement with expectations based on our current understanding of galaxy clusters and cosmology. In combination with other cosmological probes, we use the cluster catalog to improve estimates of cosmological parameters. Assuming a standard spatially flat wCDM cosmological model, the addition of our catalog to the WMAP 7-year analysis yields sigma_8 = 0.81 +- 0.09 and w = -1.07 +- 0.29, a ~50% improvement in precision on both parameters over WMAP7 alone.

We present a study of the cosmic infrared background, which is a measure of the dust obscured activity in all galaxies in the Universe. We venture to isolate the galaxies responsible for the background at 1mm; with spectroscopic and photometric redshifts we constrain the redshift distribution of these galaxies. We create a deep 1.16mm map (sigma ~ 0.5mJy) by combining the AzTEC 1.1mm and MAMBO 1.2mm datasets in GOODS-N. This combined map contains 41 secure detections, 13 of which are new. By averaging the 1.16mm flux densities of individually undetected galaxies with 24um flux densities > 25uJy, we resolve 31--45 per cent of the 1.16mm background. Repeating our analysis on the SCUBA 850um map, we resolve a higher percentage (40--64 per cent) of the 850um background. A majority of the background resolved (attributed to individual galaxies) at both wavelengths comes from galaxies at z > 1.3. If the ratio of the resolved submillimeter to millimeter background is applied to a reasonable scenario for the origins of the unresolved submillimeter background, 60--88 per cent of the total 1.16mm background comes from galaxies at z > 1.3.

[NIRSS is one of three concepts that contributed to the Wide-Field Infrared Survey Telescope (WFIRST) mission advocated by the Decadal Survey.] Operating beyond the reaches of the Earth's atmosphere, free of its limiting absorption and thermal background, the Near-Infrared Sky Surveyor (NIRSS) will deeply map the entire sky at near-infrared wavelengths, thereby enabling new and fundamental discoveries ranging from the identification of extrasolar planets to probing the reionization epoch by identifying thousands of quasars at z>10. NIRSS will directly address the NASA scientific objective of studying cosmic origins by using a 1.5-meter telescope to reach full-sky 0.2 uJy (25.6 mag AB) sensitivities in four passbands from 1 to 4 microns in a 4-yr mission. At the three shorter passbands (1 - 2.5 microns), the proposed depth is comparable to the deepest pencil-beam surveys done to date and is 3000 times more sensitive than the only previous all-sky near-infrared survey, 2MASS. At the longest passband (3.5 micron), which is not feasible from the ground, NIRSS will be 500 times more sensitive than WISE. NIRSS fills a pivotal gap in our knowledge of the celestial sphere, is a natural complement to WISE, and is well matched to the next generation of deep (0.1 uJy), wide-area (>2 pi ster), ground-based optical surveys (LSST and Pan-Starrs). With the high thermal backgrounds of ground-based infrared observations, a near-infrared full sky survey at sub-uJy sensitivity is only feasible from space.

We measure the rest-frame colors (dust-corrected), infrared luminosities, star formation rates, and stellar masses of 92 galaxies in a Spitzer-selected cluster at z=1.62. By fitting spectral energy distributions (SEDs) to 10-band photometry (0.4 micron < lambda(obs) <8 micron) and measuring 24 micron fluxes for the 12 spectroscopically confirmed and 80 photometrically selected members, we discover an exceptionally high level of star formation in the cluster core of ~1700 Msun/yr per Mpc^2. The cluster galaxies define a strong blue sequence in (U-V) color and span a range in color. We identify 17 members with L(IR)>10^(11) Lsun, and these IR luminous members follow the same trend of increasing star formation with stellar mass that is observed in the field at z~2. Using rates derived from both the 24 micron imaging and SED fitting, we find that the relative fraction of star-forming members triples from the lowest to highest galaxy density regions, e.g. the IR luminous fraction increases from ~8% at Sigma~10 gal per Mpc^2 to ~25% at Sigma>100 gal per Mpc^2. The observed increase is a reversal of the well-documented trend at z<1 and signals that we have reached the epoch when massive cluster galaxies are still forming a substantial fraction of their stars.

We report the spectroscopic confirmation of SPT-CL J0546-5345 at <z> = 1.067. To date this is the most distant cluster to be spectroscopically confirmed from the 2008 South Pole Telescope (SPT) catalog, and indeed the first z > 1 cluster discovered by the Sunyaev-Zel'dovich Effect (SZE). We identify 21 secure spectroscopic members within 0.9 Mpc of the SPT cluster position, 18 of which are quiescent, early-type galaxies. From these quiescent galaxies we obtain a velocity dispersion of 1179^{+232}_{-167} km/s, ranking SPT-CL J0546-5345 as the most dynamically massive cluster yet discovered at z > 1. Assuming that SPT-CL J0546-5345 is virialized, this implies a dynamical mass of M_200 = 1.0^{+0.6}_{-0.4} x 10^{15} Msun, in agreement with the X-ray and SZE mass measurements. Combining masses from several independent measures leads to a best-estimate mass of M_200 = (7.95 +/- 0.92) x 10^{14} Msun. The spectroscopic confirmation of SPT-CL J0546-5345, discovered in the wide-angle, mass-selected SPT cluster survey, marks the onset of the high redshift SZE-selected galaxy cluster era.

We present composite 3.6 and 4.5 micron luminosity functions for cluster galaxies measured from the Spitzer Deep, Wide-Field Survey (SDWFS) for 0.3<z<2. We compare the evolution of m* for these luminosity functions to models for passively evolving stellar populations to constrain the primary epoch of star formation in massive cluster galaxies. At low redshifts (z < 1.3) our results agree well with models with no mass assembly and passively evolving stellar populations with a luminosity-weighted mean formation redshift zf=2.4 assuming a Kroupa initial mass function (IMF). We conduct a thorough investigation of systematic biases that might influence our results, and estimate systematic uncertainites of Delta zf=(+0.16-0.18) (model normalization), Delta zf=(+0.40-0.05) (alpha), and Delta zf=(+0.30-0.45) (choice of stellar population model). For a Salpeter type IMF, the typical formation epoch is thus strongly constrained to be z ~2-3. Higher formation redshifts can only be made consistent with the data if one permits an evolving IMF that is bottom-light at high redshift, as suggested by van Dokkum et al 2008. At high redshift (z > 1.3) we also witness a statistically significant (>5sigma) disagreement between the measured luminosity function and the continuation of the passive evolution model from lower redshifts. After considering potential systematic biases that might influence our highest redshift data points, we interpret the observed deviation as potential evidence for ongoing mass assembly at this epoch.

We present results of X-ray observations of a sample of 15 clusters selected via their imprint on the cosmic microwave background (CMB) from the thermal Sunyaev-Zel'dovich (SZ) effect. These clusters are a subset of the first SZ-selected cluster catalog, obtained from observations of 178 deg^2 of sky surveyed by the South Pole Telescope. Using X-ray observations with Chandra and XMM-Newton, we estimate the temperature, T_X, and mass, M_g, of the intracluster medium (ICM) within r_500 for each cluster. From these, we calculate Y_X=M_g T_X and estimate the total cluster mass using a M_500-Y_X scaling relation measured from previous X-ray studies. The integrated Comptonization, Y_SZ, is derived from the SZ measurements, using additional information from the X-ray measured gas density profiles and a universal temperature profile. We calculate scaling relations between the X-ray and SZ observables, and find results generally consistent with other measurements and the expectations from simple self-similar behavior. Specifically, we fit a Y_SZ-Y_X relation and find a normalization of 0.82 +- 0.07, marginally consistent with the predicted ratio of Y_SZ/Y_X=0.91+-0.01 that would be expected from the density and temperature models used in this work. Using the Y_X derived mass estimates, we fit a Y_SZ-M_500 relation and find a slope consistent with the self-similar expectation of Y_SZ ~ M^5/3 with a normalization consistent with predictions from other X-ray studies. We compare the X-ray mass estimates to previously published SZ mass estimates derived from cosmological simulations of the SPT survey. We find that the SZ mass estimates are lower by a factor of 0.89+-0.06, which is within the ~15% systematic uncertainty quoted for the simulation-based SZ masses.

We use the multi-epoch, mid-infrared Spitzer Deep, Wide-Field Survey to investigate the variability of 474,179 objects in 8.1 deg^2 of the NDWFS Bootes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 micron bands. We find that 1.1% of the studied sample meet our standard selection criteria for being classed as a variable source. We require that the 3.6 and 4.5 micron light-curves are strongly correlated (r>0.8) and that their joint variance exceeds that for all sources with the same magnitude by 2 sigma. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBootes survey, radio catalogs, 24 micron-selected AGN candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects. Most of the stellar, galaxy and unclassified sources are false positives. For our standard selection criteria, 11-12% of the mid-IR counterparts to X-ray sources, 24 micron-selected AGN candidates and spectroscopically identified AGNs show variability. Mid-IR AGN variability can be well described by a single power-law structure function with a power-law index of 0.5 at both 3.6 and 4.5 microns, and an amplitude of 0.1 mag on rest-frame time scales of 2 years. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities. (Abridged)

We present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria.

We report the discovery of a galaxy cluster at z=1.62 located in the Spitzer Wide-Area Infrared Extragalactic survey XMM-LSS field. This structure was selected solely as an overdensity of galaxies with red Spitzer/IRAC colors, satisfying [3.6]-[4.5] > -0.1 AB mag. Photometric redshifts derived from Subaru XMM Deep Survey (BViz-bands), UKIRT Infrared Deep Survey-Ultra-Deep Survey (UKIDSS-UDS, JK-bands), and from the Spitzer Public UDS survey (3.6-8.0 micron) show that this cluster corresponds to a surface density of galaxies at z ~ 1.6 that is more than 20 sigma above the mean at this redshift. We obtained optical spectroscopic observations of galaxies in the cluster region using IMACS on the Magellan telescope. We measured redshifts for seven galaxies in the range z=1.62-1.63 within 2.8 arcmin (<1.4 Mpc) of the astrometric center of the cluster. A posteriori analysis of the XMM data in this field reveal a weak (4 sigma) detection in the [0.5-2 keV] band compatible with the expected thermal emission from such a cluster. The color-magnitude diagram of the galaxies in this cluster shows a prominent red-sequence, dominated by a population of red galaxies with (z-J) > 1.7 mag. The photometric redshift probability distributions for the red galaxies are strongly peaked at z=1.62, coincident with the spectroscopically confirmed galaxies. The rest-frame (U-B) color and scatter of galaxies on the red-sequence are consistent with a mean luminosity-weighted age of 1.2 +/- 0.1 Gyr, yielding a formation redshift z_f = 2.35 +/- 0.10, and corresponding to the last significant star-formation period in these galaxies.

We have identified a sample of cool field brown dwarf candidates using IRAC data from the Spitzer Deep, Wide-Field Survey (SDWFS). The candidates were selected from 400,000 SDWFS sources with [4.5] <= 18.5 mag and required to have [3.6]-[4.5] >= 1.5 and [4.5] - [8.0] <= 2.0 on the Vega system. The first color requirement selects objects redder than all but a handful of presently known brown dwarfs with spectral classes later than T7, while the second eliminates 14 probable reddened AGN. Optical detection of 4 of the remaining 18 sources implies they are likely also AGN, leaving 14 brown dwarf candidates. For two of the brightest candidates (SDWFS J143524.44+335334.6 and SDWFS J143222.82+323746.5), the spectral energy distributions including near-infrared detections suggest a spectral class of ~ T8. The proper motion is < 0.25 "/yr, consistent with expectations for a luminosity inferred distance of >70 pc. The reddest brown dwarf candidate (SDWFS J143356.62+351849.2) has [3.6] - [4.5]=2.24 and H - [4.5] > 5.7, redder than any published brown dwarf in these colors, and may be the first example of the elusive Y-dwarf spectral class. Models from Burrows et al. (2003) predict larger numbers of cool brown dwarfs should be found for a Chabrier (2003) mass function. Suppressing the model [4.5] flux by a factor of two, as indicated by previous work, brings the Burrows models and observations into reasonable agreement. The recently launched Wide-field Infrared Survey Explorer (WISE) will probe a volume ~40x larger and should find hundreds of brown dwarfs cooler than T7.

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns based on the multi-wavelength photometric observations of the NOAO Deep-Wide Field Survey Bootes field and the spectroscopic observations of the AGN and Galaxy Evolution Survey. Our training sample is comprised of 14448 galaxies in the redshift range 0<~z<~1 and 5347 likely AGNs in the range 0<~z<~5.58. We use our templates to determine photometric redshifts for galaxies and AGNs. While they are relatively accurate for galaxies, their accuracies for AGNs are a strong function of the luminosity ratio between the AGN and galaxy components. Somewhat surprisingly, the relative luminosities of the AGN and its host are well determined even when the photometric redshift is significantly in error. We also use our templates to study the mid-IR AGN selection criteria developed by Stern et al.(2005) and Lacy et al.(2004). We find that the Stern et al.(2005) criteria suffers from significant incompleteness when there is a strong host galaxy component and at z =~ 4.5, when the broad Halpha emission line is redshifted into the [3.6] band, but that it is little contaminated by low and intermediate redshift galaxies. The Lacy et al.(2004) criterion is not affected by incompleteness at z =~ 4.5 and is somewhat less affected by strong galaxy host components, but is heavily contaminated by low redshift star forming galaxies. Finally, we use our templates to predict the color-color distribution of sources in the upcoming WISE mission and define a color criterion to select AGNs analogous to those developed for IRAC photometry. We estimate that in between 640,000 and 1,700,000 AGNs will be identified by these criteria, but will have serious completeness problems for z >~ 3.4.

We present redshifts and optical richness properties of 21 galaxy clusters uniformly selected by their Sunyaev-Zel'dovich signature. These clusters, plus an additional, unconfirmed candidate, were detected in a 178 square-degree area surveyed by the South Pole Telescope in 2008. Using griz imaging from the Blanco Cosmology Survey and from pointed Magellan telescope observations, as well as spectroscopy using Magellan facilities, we confirm the existence of clustered red-sequence galaxies, report red-sequence photometric redshifts, present spectroscopic redshifts for a subsample, and derive R_200 radii and M_200 masses from optical richness. The clusters span redshifts from 0.15 to greater than 1, with a median redshift of 0.74; three clusters are estimated to be at z > 1. Redshifts inferred from mean red-sequence colors exhibit 2% RMS scatter in sigma_z/(1+z) with respect to the spectroscopic subsample for z < 1. We show that M_200 cluster masses derived from optical richness correlate with masses derived from South Pole Telescope data and agree with previously derived scaling relations to within the uncertainties. Optical and infrared imaging is an efficient means of cluster identification and redshift estimation in large Sunyaev-Zel'dovich surveys, and exploiting the same data for richness measurements, as we have done, will be useful for constraining cluster masses and radii for large samples in cosmological analysis.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

In this work we present the first results of our imaging campaign at Keck Observatory to identify the host galaxies of "dark" gamma-ray bursts (GRBs), events with no detected optical afterglow or with detected optical flux significantly fainter than expected from the observed X-ray afterglow. We find that out of a uniform sample of 29 Swift bursts rapidly observed by the Palomar 60-inch telescope through March 2008 (14 of which we classify as dark), all events have either a detected optical afterglow, a probable optical host-galaxy detection, or both. Our results constrain the fraction of Swift GRBs coming from very high redshift (z > 7), such as the recent GRB 090423, to between 0.2-7 percent at 80% confidence. In contrast, a significant fraction of the sample requires large extinction columns (host-frame A_V > 1 mag, with several events showing A_V > 2-6 mag), identifying dust extinction as the dominant cause of the dark GRB phenomenon. We infer that a significant fraction of GRBs (and, by association, of high-mass star formation) occurs in highly obscured regions. However, the host galaxies of dark GRBs seem to have normal optical colors, suggesting that the source of obscuring dust is local to the vicinity of the GRB progenitor or highly unevenly distributed within the host galaxy.

Galaxy groups are key tracers of galaxy evolution, cluster evolution, and structure formation, yet they are difficult to study at even moderate redshift. We have undertaken a project to observe a flux-limited sample of intermediate-redshift (0.1 < z < 0.5) group candidates identified by the XBootes Chandra survey. When complete, this project will nearly triple the current number of groups with measured temperatures in this redshift range. Here we present deep Suzaku/XIS and Chandra/ACIS follow-up observations of the first 10 targets in this project; all are confirmed sources of diffuse, thermal emission with derived temperatures and luminosities indicative of rich groups/poor clusters. By exploiting the multi-wavelength coverage of the XBootes/NOAO Deep Wide Field Survey (NDWFS) field, we aim to (1) constrain non-gravitational effects that alter the energetics of the intragroup medium, and (2) understand the physical connection between the X-ray and optical properties of groups. We discuss the properties of the current group sample in the context of observed cluster scaling relations and group and cluster evolution and outline the future plans for this project.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

We are proposing to conduct a multicolor, synoptic infrared (IR) imaging survey of the Northern sky with a new, dedicated 6.5-meter telescope at San Pedro M\'artir (SPM) Observatory. This initiative is being developed in partnership with astronomy institutions in Mexico and the University of California. The 4-year, dedicated survey, planned to begin in 2017, will reach more than 100 times deeper than 2MASS. The Synoptic All-Sky Infrared (SASIR) Survey will reveal the missing sample of faint red dwarf stars in the local solar neighborhood, and the unprecedented sensitivity over such a wide field will result in the discovery of thousands of z ~ 7 quasars (and reaching to z > 10), allowing detailed study (in concert with JWST and Giant Segmented Mirror Telescopes) of the timing and the origin(s) of reionization. As a time-domain survey, SASIR will reveal the dynamic infrared universe, opening new phase space for discovery. Synoptic observations of over 10^6 supernovae and variable stars will provide better distance measures than optical studies alone. SASIR also provides significant synergy with other major Astro2010 facilities, improving the overall scientific return of community investments. Compared to optical-only measurements, IR colors vastly improve photometric redshifts to z ~ 4, enhancing dark energy and dark matter surveys based on weak lensing and baryon oscillations. The wide field and ToO capabilities will enable a connection of the gravitational wave and neutrino universe - with events otherwise poorly localized on the sky - to transient electromagnetic phenomena.

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we present low-resolution (64 < lambda / dlambda < 124), mid-infrared (20-38 micron) spectra of 23 high-redshift ULIRGs detected in the Bootes field of the NOAO Deep Wide-Field Survey. All of the sources were selected to have 1) fnu(24 micron) > 0.5 mJy; 2) R-[24] > 14 Vega mag; and 3) a prominent rest-frame 1.6 micron stellar photospheric feature redshifted into Spitzer's 3-8 micron IRAC bands. Of these, 20 show emission from polycyclic aromatic hydrocarbons (PAHs), usually interpreted as signatures of star formation. The PAH features indicate redshifts in the range 1.5 < z < 3.0, with a mean of <z>=1.96 and a dispersion of 0.30. Based on local templates, these sources have extremely large infrared luminosities, comparable to that of submillimeter galaxies. Our results confirm previous indications that the rest-frame 1.6 micron stellar bump can be efficiently used to select highly obscured starforming galaxies at z~2, and that the fraction of starburst-dominated ULIRGs increases to faint 24 micron flux densities. Using local templates, we find that the observed narrow redshift distribution is due to the fact that the 24 micron detectability of PAH-rich sources peaks sharply at z = 1.9. We can analogously explain the broader redshift distribution of Spitzer-detected AGN-dominated ULIRGs based on the shapes of their SEDs. Finally, we conclude that z~2 sources with a detectable 1.6 micron stellar opacity feature lack sufficient AGN emission to veil the 7.7 micron PAH band.

The Spitzer Deep, Wide-Field Survey (SDWFS) is a four-epoch infrared survey of ten square degrees in the Bootes field of the NOAO Deep Wide-Field Survey using the IRAC instrument on the Spitzer Space Telescope. SDWFS, a Cycle four Spitzer Legacy project, occupies a unique position in the area-depth survey space defined by other Spitzer surveys. The four epochs that make up SDWFS permit -- for the first time -- the selection of infrared-variable and high proper motion objects over a wide field on timescales of years. Because of its large survey volume, SDWFS is sensitive to galaxies out to z~3 with relatively little impact from cosmic variance for all but the richest systems. The SDWFS datasets will thus be especially useful for characterizing galaxy evolution beyond z~1.5. This paper explains the SDWFS observing strategy and data processing, presents the SDWFS mosaics and source catalogs, and discusses some early scientific findings. The publicly-released, full-depth catalogs contain 6.78, 5.23, 1.20, and 0.96 x 10e5 distinct sources detected to the average 5-sigma, 4" diameter, aperture-corrected limits of 19.77, 18.83, 16.50, and 15.82 Vega mag at 3.6, 4.5, 5.8, and 8.0 micron, respectively. The SDWFS number counts and color-color distribution are consistent with other, earlier Spitzer surveys. At the 6 min integration time of the SDWFS IRAC imaging, more than 50% of isolated FIRST radio sources and more than 80% of on-axis XBootes sources are detected out to 8.0 micron. Finally, we present the four highest proper motion IRAC-selected sources identified from the multi-epoch imaging, two of which are likely field brown dwarfs of mid-T spectral class.

We explore the connection between different classes of active galactic nuclei (AGNs) and the evolution of their host galaxies, by deriving host galaxy properties, clustering, and Eddington ratios of AGNs selected in the radio, X-ray, and infrared. We study a sample of 585 AGNs at 0.25 < z < 0.8 using redshifts from the AGN and Galaxy Evolution Survey (AGES) and data in the radio (WSRT 1.4 GHz), X-rays (Chandra XBootes), and mid-IR (IRAC Shallow Survey). The radio, X-ray, and IR AGN samples show modest overlap, indicating that to the flux limits of the survey, they represent largely distinct classes of AGNs. We derive host galaxy colors and luminosities, as well as Eddington ratios (lambda), for obscured or optically faint AGNs. We also measure the two-point cross-correlation between AGNs and galaxies on scales of 0.3-10 h^-1 Mpc, and derive typical dark matter halo masses. We find that: (1) radio AGNs are mainly found in luminous red galaxies, are strongly clustered (with M_halo ~ 3x10^13 h^-1 M_sun), and have very low lambda <~ 10^-3; (2) X-ray-selected AGNs are preferentially found in galaxies in the "green valley" of color-magnitude space and are clustered similarly to typical AGES galaxies (M_halo ~ 10^13 h^-1 M_sun), with 10^-3 <~ lambda <~ 1; (3) IR AGNs reside in slightly bluer, less luminous galaxies than X-ray AGNs, are weakly clustered (M_halo <~ 10^12 h^-1 M_sun), and have lambda > 10^-2. We interpret these results in terms of a simple model of AGN and galaxy evolution, whereby a "quasar" phase and the growth of the stellar bulge occurs when a galaxy's dark matter halo reaches a critical mass between ~10^12 and 10^13 M_sun. Subsequently, star formation ceases and AGN accretion shifts from radiatively efficient (optical- and IR- bright) to radiatively inefficient (optically-faint, radio-bright) modes.

We present galaxy luminosity functions at 3.6, 4.5, 5.8, and 8.0 micron measured by combining photometry from the IRAC Shallow Survey with redshifts from the AGN and Galaxy Evolution Survey of the NOAO Deep Wide-Field Survey Bootes field. The well-defined IRAC samples contain 3800-5800 galaxies for the 3.6-8.0 micron bands with spectroscopic redshifts and z < 0.6. We obtained relatively complete luminosity functions in the local redshift bin of z < 0.2 for all four IRAC channels that are well fit by Schechter functions. We found significant evolution in the luminosity functions for all four IRAC channels that can be fit as an evolution in M* with redshift, \Delta M* = Qz. While we measured Q=1.2\pm0.4 and 1.1\pm0.4 in the 3.6 and 4.5 micron bands consistent with the predictions from a passively evolving population, we obtained Q=1.8\pm1.1 in the 8.0 micron band consistent with other evolving star formation rate estimates. We compared our LFs with the predictions of semi-analytical galaxy formation and found the best agreement at 3.6 and 4.5 micron, rough agreement at 8.0 micron, and a large mismatch at 5.8 micron. These models also predicted a comparable Q value to our luminosity functions at 8.0 micron, but predicted smaller values at 3.6 and 4.5 micron. We also measured the luminosity functions separately for early and late-type galaxies. While the luminosity functions of late-type galaxies resemble those for the total population, the luminosity functions of early-type galaxies in the 3.6 and 4.5 micron bands indicate deviations from the passive evolution model, especially from the measured flat luminosity density evolution. Combining our estimates with other measurements in the literature, we found (53\pm18)% of the present stellar mass of early-type galaxies has been assembled at z=0.7.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

We present the surface density of luminous active galactic nuclei (AGN) associated with a uniformly selected galaxy cluster sample identified in the 8.5 square degree Bootes field of the NOAO Deep Wide-Field Survey. The clusters are distributed over a large range of redshift (0 < z < 1.5) and we identify AGN using three different selection criteria: mid-IR color, radio luminosity, and X-ray luminosity. Relative to the field, we note a clear overdensity of the number of AGN within 0.5 Mpc of the cluster centers at z > 0.5. The amplitude of this AGN overdensity increases with redshift. Although there are significant differences between the AGN populations probed by each selection technique, the rise in cluster AGN surface density generally increases more steeply than that of field quasars. In particular, X-ray selected AGN are at least three times more prevalent in clusters at 1 < z < 1.5 compared to clusters at 0.5 < z < 1. This effect is stronger than can be explained by the evolving median richness of our cluster sample. We thus confirm the existence of a Butcher-Oemler type effect for AGN in galaxy clusters, with the number of AGN in clusters increasing with redshift.

We present the angular autocorrelation function of 2603 dust-obscured galaxies (DOGs) in the Bootes field of the NOAO Deep Wide-Field Survey. DOGs are red, obscured galaxies, defined as having R-[24] \ge 14 (F_24/F_R \ga 1000). Spectroscopy indicates that they are located at 1.5 \la z \la 2.5. We find strong clustering, with r_0 = 7.40^{+1.27}_{-0.84} Mpc/h for the full F_24 > 0.3 mJy sample. The clustering and space density of the DOGs are consistent with those of submillimeter galaxies, suggestive of a connection between these populations. We find evidence for luminosity-dependent clustering, with the correlation length increasing to r_0 = 12.97^{+4.26}_{-2.64} Mpc/h for brighter (F_24 > 0.6 mJy) DOGs. Bright DOGs also reside in richer environments than fainter ones, suggesting these subsamples may not be drawn from the same parent population. The clustering amplitudes imply average halo masses of log M = 12.2^{+0.3}_{-0.2} Msun for the full DOG sample, rising to log M = 13.0^{+0.4}_{-0.3} Msun for brighter DOGs. In a biased structure formation scenario, the full DOG sample will, on average, evolve into ~ 3 L* present-day galaxies, whereas the most luminous DOGs may evolve into brightest cluster galaxies.

We use deep far-IR, submm, radio and X-ray imaging and mid-IR spectroscopy to explore the nature of a sample of Spitzer-selected dust-obscured galaxies (DOGs) in GOODS-N. A sample of 79 galaxies satisfy the criteria R-[24]>14 (Vega) down to S24>100 microJy. Twelve of these galaxies have IRS spectra available which we use to measure redshifts and classify these objects as being dominated by star formation or active galactic nuclei (AGN) activity in the mid-IR. The IRS spectra and Spitzer photometric redshifts confirm that the DOGs lie in a tight redshift distribution around z~2. Based on mid-IR colors, 80% of DOGs are likely dominated by star formation; the stacked X-ray emission from this sub-sample of DOGs is also consistent with star formation. Since only a small number of DOGs are individually detected at far-IR and submm wavelengths, we use a stacking analysis to determine the average flux from these objects and plot a composite IR (8-1000 microns) spectral energy distribution (SED). The average luminosity of these star forming DOGs is LIR~1e12 Lsun. We compare the average star forming DOG to the average bright (S850>5 mJy) submillimeter galaxy (SMG); the S24>100 microJy DOGs are 3 times more numerous but 8 times less luminous in the IR. The far-IR SED shape of DOGs is similar to that of SMGs (average dust temperature of around 30 K) but DOGs have a higher mid-IR to far-IR flux ratio. The average star formation-dominated DOG has a star formation rate of 200 Msun/yr which, given their space density, amounts to a contribution of 0.01 Msun/yr/Mpc3 (or 5-10%) to the star formation rate density at z~2. We use the composite SED to predict the average flux of DOGs in future Herschel/PACS 100 micron and SCUBA-2 450 micron surveys and show that the majority of them will be detected.

We have identified 335 galaxy cluster and group candidates spanning 0<z<2, using a 4.5um selected sample of galaxies in a 7.25 deg^2 region in the Spitzer/IRAC Shallow Survey. Using full redshift probability distributions for all galaxies, clusters were identified as 3-dimensional overdensities using a wavelet algorithm. To date 12 clusters at z>1, and over 60 at z < 0.5 have been spectroscopically confirmed. The mean I-[3.6] color for cluster galaxies up to z~1 is well matched by a z_f=3 passively evolving model. At z>1, a wider range of formation histories is needed, but higher formation redshifts (i.e. z_f \ge 4-5) are favored for most clusters. The cluster autocorrelation function, measured for the first time out to z=1.5, is found not to have evolved over the last 10 Gyr, in agreement with the prediction from LCDM. The average mass of the IRAC Shallow Cluster Survey sample, inferred from its clustering, is ~ 10^14 Msun.

We have identified 335 galaxy cluster and group candidates, 106 of which are at z > 1, using a 4.5 um selected sample of objects from a 7.25 deg^2 region in the Spitzer Infrared Array Camera (IRAC) Shallow Survey. Clusters were identified as 3-dimensional overdensities using a wavelet algorithm, based on photometric redshift probability distributions derived from IRAC and NOAO Deep Wide-Field Survey data. We estimate only ~10% of the detections are spurious. To date 12 of the z > 1 candidates have been confirmed spectroscopically, at redshifts from 1.06 to 1.41. Velocity dispersions of ~750 km/s for two of these argue for total cluster masses well above 10^14 M_sun, as does the mass estimated from the rest frame near infrared stellar luminosity. Although not selected to contain a red sequence, some evidence for red sequences is present in the spectroscopically confirmed clusters, and brighter galaxies are systematically redder than the mean galaxy color in clusters at all redshifts. The mean I - [3.6] color for cluster galaxies up to z ~ 1 is well matched by a passively evolving model in which stars are formed in a 0.1 Gyr burst starting at redshift z_f = 3. At z > 1, a wider range of formation histories is needed, but higher formation redshifts (i.e. z_f > 3) are favored for most clusters.

We have traced the past 7 Gyr of red galaxy stellar mass growth within dark matter halos. We have determined the halo occupation distribution, which describes how galaxies reside within dark matter halos, using the observed luminosity function and clustering of 40,696 0.2<z<1.0 red galaxies in Bootes. Half of 10^{11.9} Msun/h halos host a red central galaxy, and this fraction increases with increasing halo mass. We do not observe any evolution of the relationship between red galaxy stellar mass and host halo mass, although we expect both galaxy stellar masses and halo masses to evolve over cosmic time. We find that the stellar mass contained within the red population has doubled since z=1, with the stellar mass within red satellite galaxies tripling over this redshift range. In cluster mass halos most of the stellar mass resides within satellite galaxies and the intra-cluster light, with a minority of the stellar mass residing within central galaxies. The stellar masses of the most luminous red central galaxies are proportional to halo mass to the power of a third. We thus conclude that halo mergers do not always lead to rapid growth of central galaxies. While very massive halos often double in mass over the past 7 Gyr, the stellar masses of their central galaxies typically grow by only 30%.

We combine the results of the Spitzer IRAC Shallow Survey and the Chandra XBootes Survey of the 8.5 square degrees Bootes field of the NOAO Deep Wide- Field Survey to produce the largest comparison of mid-IR and X-ray sources to date. The comparison is limited to sources with X-ray fluxes >8x10-15 erg cm-2s-1 in the 0.5-7.0 keV range and mid-IR sources with 3.6 um fluxes brighter than 18.4 mag (12.3 uJy). In this most sensitive IRAC band, 85% of the 3086 X-ray sources have mid-IR counterparts at an 80% confidence level based on a Bayesian matching technique. Only 2.5% of the sample have no IRAC counterpart at all based on visual inspection. Even for a smaller but a significantly deeper Chandra survey in the same field, the IRAC Shallow Survey recovers most of the X-ray sources. A majority (65%) of the Chandra sources detected in all four IRAC bands occupy a well-defined region of IRAC [3.6] - [4.5] vs [5.8] - [8.0] color-color space. These X-ray sources are likely infrared luminous, unobscured type I AGN with little mid-infrared flux contributed by the AGN host galaxy. Of the remaining Chandra sources, most are lower luminosity type I and type II AGN whose mid-IR emission is dominated by the host galaxy, while approximately 5% are either Galactic stars or very local galaxies.

We present the redshift distribution of a complete, unbiased sample of 24 micron sources down to fnu(24 micron) = 300 uJy (5-sigma). The sample consists of 591 sources detected in the Bootes field of the NOAO Deep Wide-Field Survey. We have obtained optical spectroscopic redshifts for 421 sources (71%). These have a redshift distribution peaking at z~0.3, with a possible additional peak at z~0.9, and objects detected out to z=4.5. The spectra of the remaining 170 (29%) exhibit no strong emission lines from which to determine a redshift. We develop an algorithm to estimate the redshift distribution of these sources, based on the assumption that they have emission lines but that these lines are not observable due to the limited wavelength coverage of our spectroscopic observations. The redshift distribution derived from all 591 sources exhibits an additional peak of extremely luminous (L(8-1000 micron) > 3 x 10^{12} Lsun) objects at z~2, consisting primarily of sources without observable emission lines. We use optical line diagnostics and IRAC colors to estimate that 55% of the sources within this peak are AGN-dominated. We compare our results to published models of the evolution of infrared luminous galaxies. The models which best reproduce our observations predict a large population of star-formation dominated ULIRGs at z > 1.5 rather than the AGN-dominated sources we observe.

Observations with Spitzer Space Telescope have recently revealed a significant population of high-redshift z~2 dust-obscured galaxies (DOGs) with large mid-IR to UV luminosity ratios. These galaxies have been missed in traditional optical studies of the distant universe. We present a simple method for selecting this high-z population based solely on the ratio of the observed mid-IR 24um to optical R-band flux density. In the 8.6 sq.deg Bootes NDWFS Field, we uncover ~2,600 DOG candidates (= 0.089/sq.arcmin) with 24um flux densities F24>0.3mJy and (R-[24])>14 (i.e., F[24]/F[R] > 1000). These galaxies have no counterparts in the local universe, and become a larger fraction of the population at fainter F24, representing 13% of the sources at 0.3~mJy. DOGs exhibit evidence of both star-formation and AGN activity, with the brighter 24um sources being more AGN- dominated. We have measured spectroscopic redshifts for 86 DOGs, and find a broad z distribution centered at <z>~2.0. Their space density is 2.82E-5 per cubic Mpc, similar to that of bright sub-mm-selected galaxies at z~2. These redshifts imply very large luminosities LIR>~1E12-14 Lsun. DOGs contribute ~45-100% of the IR luminosity density contributed by all z~2 ULIRGs, suggesting that our simple selection criterion identifies the bulk of z~2 ULIRGs. DOGs may be the progenitors of ~4L* present-day galaxies seen undergoing a luminous,short- lived phase of bulge and black hole growth. They may represent a brief evolution phase between SMGs and less obscured quasars or galaxies. [Abridged]

We built an optimal basis of low resolution templates for galaxies over the wavelength range from 0.2 to 10 $\mu$m using a variant of the algorithm presented by Budavari et al. (2000). We derived them using eleven bands of photometry from the NDWFS, FLAMEX, zBo\"otes and IRAC Shallow surveys for 16033 galaxies in the NDWFS Bo\"otes field with spectroscopic redshifts measured by the AGN and Galaxy Evolution Survey. We also developed algorithms to accurately determine photometric redshifts, K corrections and bolometric luminosities using these templates. Our photometric redshifts have an accuracy of $\sigma_z/(1+z) = 0.04$ when clipped to the best 95%. We used these templates to study the spectral type distribution in the field and to estimate luminosity functions of galaxies as a function of redshift and spectral type. In particular, we note that the 5-8$\mu$m color distribution of galaxies is bimodal, much like the optical g--r colors.

We identify a population of 640 obscured and 839 unobscured AGNs at redshifts 0.7<z<~3 using multiwavelength observations of the 9 deg^2 NOAO Deep Wide-Field Survey (NDWFS) region in Bootes. We select AGNs on the basis of Spitzer IRAC colors obtained by the IRAC Shallow Survey. Redshifts are obtained from optical spectroscopy or photometric redshift estimators. We classify the IR-selected AGNs as IRAGN 1 (unobscured) and IRAGN 2 (obscured) using a simple criterion based on the observed optical to mid-IR color, with a selection boundary of R-[4.5]=6.1, where R and [4.5] are the Vega magnitudes in the R and IRAC 4.5 micron bands, respectively. We verify this selection using X-ray stacking analyses with data from the Chandra XBootes survey, as well as optical photometry from NDWFS and spectroscopy from MMT/AGES. We show that (1) these sources are indeed AGNs, and (2) the optical/IR color selection separates obscured sources (with average N_H~3x10^22 cm^-2 obtained from X-ray hardness ratios, and optical colors and morphologies typical of galaxies) and unobscured sources (with no X-ray absorption, and quasar colors and morphologies), with a reliability of >~80%. The observed numbers of IRAGNs are comparable to predictions from previous X-ray, optical, and IR luminosity functions, for the given redshifts and IRAC flux limits. We observe a bimodal distribution in R-[4.5] color, suggesting that luminous IR-selected AGNs have either low or significant dust extinction, which may have implications for models of AGN obscuration.

We present the galaxy cluster autocorrelation function of 277 galaxy cluster candidates with 0.25 \le z \le 1.5 in a 7 deg^2 area of the IRAC Shallow Cluster Survey. We find strong clustering throughout our galaxy cluster sample, as expected for these massive structures. Specifically, at <z> = 0.5 we find a correlation length of r_0 = 17.40^{+3.98}_{-3.10} h^-1 Mpc, in excellent agreement with the Las Campanas Distant Cluster Survey, the only other non-local measurement. At higher redshift, <z> = 1, we find that strong clustering persists, with a correlation length of r_0=19.14^{+5.65}_{-4.56} h^-1 Mpc. A comparison with high resolution cosmological simulations indicates these are clusters with halo masses of \sim 10^{14} Msun, a result supported by estimates of dynamical mass for a subset of the sample. In a stable clustering picture, these clusters will evolve into massive (10^{15} Msun) clusters by the present day.

A SST survey in the NOAO Deep-Wide Field in Bo\"otes provides a complete, 8-micron-selected sample of galaxies to a limiting (Vega) magnitude of 13.5. In the 6.88 deg$^2$ field sampled, 79% of the 4867 galaxies have spectroscopic redshifts, allowing an accurate determination of the local (z<0.3) galaxy luminosity function. Stellar and dust emission can be separated on the basis of observed galaxy colors. Dust emission (mostly PAH) accounts for 80% of the 8 micron luminosity, stellar photospheres account for 19%, and AGN emission accounts for roughly 1 %. A sub-sample of the 8 micron-selected galaxies have blue, early-type colors, but even most of these have significant PAH emission. The luminosity functions for the total 8 micron luminosity and for the dust emission alone are both well fit by Schechter functions. For the 8 micron luminosity function, the characteristic luminosity is \nu L_{\nu}^*(8.0 \micron) = 1.8 \times 10^{10}$ \Lsun while for the dust emission alone it is 1.6 x 10^{10}$ \Lsun\null. The average 8 \micron luminosity density at z<0.3 is 3.1 x 10^7 \Lsun Mpc^{-3}, and the average luminosity density from dust alone is 2.5 x 10^7 \Lsun Mpc^{-3}. This luminos ity arises predominantly from galaxies with 8 \micron luminosities ($\nu L_{\nu}$) between $2\times 10^9$ and $2 x 10^{10}$ \Lsun, i.e., normal galaxies, not LIRGs or ULIRGs.

We describe the angular power spectrum of resolved sources at 3.6 microns (L-band) in Spitzer imaging data of the GOODS HDF-N, the GOODS CDF-S, and the NDWFS Bootes field in several source magnitude bins. We also measure angular power spectra of resolved sources in the Bootes field at K_S and J-bands using ground-based IR imaging data. In the three bands, J, K_S, and L, we detect the clustering of galaxies on top of the shot-noise power spectrum at multipoles between ell ~ 10^2 and 10^5. The angular power spectra range from the large, linear scales to small, non-linear scales of galaxy clustering, and in some magnitude ranges, show departure from a power-law clustering spectrum. We consider a halo model to describe clustering measurements and to establish the halo occup ation number parameters of IR bright galaxies at redshifts around one. We also extend our clustering results and completeness-corrected faint source number counts in GOODS fields to understand the underlying nature of unresolved sources responsible for IR background (IRB) anisotropies that were detected in deep Spitzer images. While these unresolved fluctuations were measured at sub-arcminute angular scales, if a high-redshift diffuse component associated with first galaxies exists in the IRB, then it's clustering properties are best studied with shallow, wide-field images that allow a measurement of the clustering spectrum from a few degrees to arcminute angular scales.

We have measured the contribution of submillimeter and mid-infrared sources to the extragalactic background radiation at 70 and 160um. Specifically, we have stacked flux in 70 and 160um Spitzer Space Telescope (Spitzer) observations of the Canada-UK Deep Sub-millimeter Survey 14h field at the positions of 850um sources detected by SCUBA and also 8 and 24um sources detected by Spitzer. We find that per source, the SCUBA galaxies are the strongest and the 8um sources the weakest contributors to the background flux at both 70 and 160um. Our estimate of the contribution of the SCUBA sources is higher than previous estimates. However, expressed as a total contribution, the full 8um source catalogue accounts for twice the total 24um source contribution and ~10 times the total SCUBA source contribution. The 8um sources account for the majority of the background radiation at 160um with a flux of 0.87+/-0.16 MJy/sr and at least a third at 70um with a flux of 0.103+/-0.019 MJy/sr. These measurements are consistent with current lower limits on the background at 70 and 160um. Finally, we have investigated the 70 and 160um emission from the 8 and 24um sources as a function of redshift. We find that the average 70um flux per 24um source and the average 160um flux per 8 and 24um source is constant over all redshifts, up to z~4. In contrast, the low-redshift half (z<1) of the of 8um sample contributes approximately four times the total 70um flux of the high-redshift half. These trends can be explained by a single non-evolving SED.

We use observations from the 9 square degree multiwavelength survey in Bootes to identify hundreds of obscured active galactic nuclei (AGN) with high redshifts (z > 0.7), luminosities (L_bol > 10^45 ergs/s), and moderate obscuring columns (N_H > 10^22 cm^-2), and to measure the clustering properties of X-ray AGN at z > 1. In the Bootes region, shallow (5 ks) Chandra X-ray observations have detected ~4,000 X-ray sources, and the same region has been mapped with deep optical imaging and by Spitzer IRAC, which detects ~300,000 point sources, of which ~30,000 have detections in all four IRAC bands, for which we can select AGN on the basis of their mid-IR colors. With the MMT/Hectospec we have obtained modest resolution optical spectra for about half the X-ray sources (out to z > 3) and ~20,000 galaxies (out to z = 0.7). With this multiwavelength data we select >400 AGN per square degree (compared to 12 per square degree from SDSS). Among a sample of IRAC-selected AGN we identify 641 candidate obscured objects based on their R band and IRAC luminosities. We use X-ray stacking techniques to verify that they are obscured AGN and measure their absorbing column densities. We also measure the three-dimensional two-point correlation function for X-ray selected AGN.

We trace the assembly history of red galaxies since z=1, by measuring their evolving space density with the B-band luminosity function. Our sample of 39599 red galaxies, selected from 6.96 square degrees of imaging from the NOAO Deep Wide-Field and Spitzer IRAC Shallow surveys, is an order of magnitude larger, in size and volume, than comparable samples in the literature. We measure a higher space density of z=0.9 red galaxies than some of the recent literature, in part because we account for the faint yet significant galaxy flux which falls outside of our photometric aperture. The B-band luminosity density of red galaxies, which effectively measures the evolution of ~L* galaxies, increases by only 36 percent from z=0 to z=1. If red galaxy stellar populations have faded by 1.24 B-band magnitudes since z=1, the stellar mass contained within the red galaxy population has roughly doubled over the past 8 Gyr. This is consistent with star-forming galaxies being transformed into ~L* red galaxies after a decline in their star formation rates. In contrast, the evolution of 4L* red galaxies differs only slightly from a model with negligible star formation and no galaxy mergers since z=1. If this model approximates the luminosity evolution of red galaxy stellar populations, then 80 percent of the stellar mass contained within today's 4L* red galaxies was already in place at z=0.7. While red galaxy mergers have been observed, such mergers do not produce rapid growth of 4L* red galaxy stellar masses between z=1 and the present day.

We discuss color selection of rare objects in a wide-field, multiband survey spanning from the optical to the mid-infrared. Simple color criteria simultaneously identify and distinguish two of the most sought after astrophysical sources: the coolest brown dwarfs and the most distant quasars. We present spectroscopically-confirmed examples of each class identified in the IRAC Shallow Survey of the Bootes field of the NOAO Deep Wide-Field Survey. ISS J142950.9+333012 is a T4.5 brown dwarf at a distance of approximately 42 pc, and ISS J142738.5+331242 is a radio-loud quasar at redshift z=6.12. Our selection criteria identify a total of four candidates over 8 square degrees of the Bootes field. The other two candidates are both confirmed 5.5<z<6 quasars, previously reported by Cool et al. (2006). We discuss the implications of these discoveries and conclude that there are excellent prospects for extending such searches to cooler brown dwarfs and higher redshift quasars.

We compute accurate redshift distributions to I(AB) = 24 and R(AB) = 24.5 using photometric redshifts estimated from six-band UBVRIZ photometry in the Canada-France Deep Fields-Photometric Redshift Survey (CFDF-PRS). Our photometric redshift algorithm is calibrated using hundreds of CFRS spectroscopic redshifts in the same fields. The dispersion in redshift is \sigma/(1+z) \la 0.04 to the CFRS depth of I(AB) = 22.5, rising to \sigma/(1+z) \la 0.06 at our nominal magnitude and redshift limits of I(AB) = 24 and z \le 1.3, respectively. We describe a new method to compute N(z) that incorporates the full redshift likelihood functions in a Bayesian iterative analysis and we demonstrate in extensive Monte Carlo simulations that it is superior to distributions calculated using simple maximum likelihood redshifts. The field-to-field differences in the redshift distributions, while not unexpected theoretically, are substantial even on 30' scales. We provide I(AB) and R(AB) redshift distributions, median redshifts, and parametrized fits of our results in various magnitude ranges, accounting for both random and systematic errors in the analysis.