We present a search for anisotropic cosmic birefringence in 500 deg$^2$ of southern sky observed at 150 GHz with the SPTpol camera on the South Pole Telescope. We reconstruct a map of cosmic polarization rotation anisotropies using higher-order correlations between the observed cosmic microwave background (CMB) $E$ and $B$ fields. We then measure the angular power spectrum of this map, which is found to be consistent with zero. The non-detection is translated into an upper limit on the amplitude of the scale-invariant cosmic rotation power spectrum, $L(L+1)C_L^{\alpha\alpha}/2\pi < 0.10 \times 10^{-4}$ rad$^2$ (0.033 deg$^2$, 95% C.L.). This upper limit can be used to place constraints on the strength of primordial magnetic fields, $B_{1 \rm Mpc} < 17 {\rm nG} $ (95% C.L.), and on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma} < 4.0 \times 10^{-2}/H_I $ (95% C.L.), where $H_I$ is the inflationary Hubble scale. For the first time, we also cross-correlate the CMB temperature fluctuations with the reconstructed rotation angle map, a signal expected to be non-vanishing in certain theoretical scenarios, and find no detectable signal. We perform a suite of systematics and consistency checks and find no evidence for contamination.

We report a B-mode power spectrum measurement from the cosmic microwave background (CMB) polarization anisotropy observations made using the SPTpol instrument on the South Pole Telescope. This work uses 500 deg$^2$ of SPTpol data, a five-fold increase over the last SPTpol B-mode release. As a result, the bandpower uncertainties have been reduced by more than a factor of two, and the measurement extends to lower multipoles: $52 < \ell < 2301$. Data from both 95 and 150 GHz are used, allowing for three cross-spectra: 95 GHz x 95 GHz, 95 GHz x 150 GHz, and 150 GHz x 150 GHz. B-mode power is detected at very high significance; we find $P(BB < 0) = 5.8 \times 10^{-71}$, corresponding to a $18.1 \sigma$ detection of power. An upper limit is set on the tensor-to-scalar ratio, $r < 0.44$ at 95% confidence (the expected $1 \sigma$ constraint on $r$ given the measurement uncertainties is 0.22). We find the measured B-mode power is consistent with the Planck best-fit $\Lambda$CDM model predictions. Scaling the predicted lensing B-mode power in this model by a factor Alens, the data prefer Alens = $1.17 \pm 0.13$. These data are currently the most precise measurements of B-mode power at $\ell > 320$.

We report new measurements of millimeter-wave power spectra in the angular multipole range $2000 \le \ell \le 11,000$ (angular scales $5^\prime \gtrsim \theta \gtrsim 1^\prime$). By adding 95 and 150\,GHz data from the low-noise 500 deg$^2$ SPTpol survey to the SPT-SZ three-frequency 2540 deg$^2$ survey, we substantially reduce the uncertainties in these bands. These power spectra include contributions from the primary cosmic microwave background, cosmic infrared background, radio galaxies, and thermal and kinematic Sunyaev-Zel'dovich (SZ) effects. The data favor a thermal SZ (tSZ) power at 143\,GHz of $D^{\rm tSZ}_{3000} = 3.42 \pm 0.54~ \mu {\rm K}^2$ and a kinematic SZ (kSZ) power of $D^{\rm kSZ}_{3000} = 3.0 \pm 1.0~ \mu {\rm K}^2$. This is the first measurement of kSZ power at $\ge 3\,\sigma$. We study the implications of the measured kSZ power for the epoch of reionization, finding the duration of reionization to be $\Delta z_{re} = 1.0^{+1.6}_{-0.7}$ ($\Delta z_{re}< 4.1$ at 95% confidence), when combined with our previously published tSZ bispectrum measurement.

We present cosmological constraints based on the cosmic microwave background (CMB) lensing potential power spectrum measurement from the recent 500 deg$^2$ SPTpol survey, the most precise CMB lensing measurement from the ground to date. We fit a flat $\Lambda$CDM model to the reconstructed lensing power spectrum alone and in addition with other data sets: baryon acoustic oscillations (BAO) as well as primary CMB spectra from Planck and SPTpol. The cosmological constraints based on SPTpol and Planck lensing band powers are in good agreement when analysed alone and in combination with Planck full-sky primary CMB data. With weak priors on the baryon density and other parameters, the CMB lensing data alone provide a 4\% constraint on $\sigma_8\Omega_m^{0.25} = 0.0593 \pm 0.025$.. Jointly fitting with BAO data, we find $\sigma_8=0.779 \pm 0.023$, $\Omega_m = 0.368^{+0.032}_{-0.037}$, and $H_0 = 72.0^{+2.1}_{-2.5}\,\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1} $, up to $2\,\sigma$ away from the central values preferred by Planck lensing + BAO. However, we recover good agreement between SPTpol and Planck when restricting the analysis to similar scales. We also consider single-parameter extensions to the flat $\Lambda$CDM model. The SPTpol lensing spectrum constrains the spatial curvature to be $\Omega_K = -0.0007 \pm 0.0025$ and the sum of the neutrino masses to be $\sum m_{\nu} < 0.23$ eV at 95\% C.L. (with Planck primary CMB and BAO data), in good agreement with the Planck lensing results. With the differences in the $S/N$ of the lensing modes and the angular scales covered in the lensing spectra, this analysis represents an important independent check on the full-sky Planck lensing measurement.

We study the polarisation properties of extragalactic sources at 95 and 150 GHz in the SPTpol 500 deg$^2$ survey. We estimate the polarised power by stacking maps at known source positions, and correct for noise bias by subtracting the mean polarised power at random positions in the maps. We show that the method is unbiased using a set of simulated maps with similar noise properties to the real SPTpol maps. We find a flux-weighted mean-squared polarisation fraction $\langle p^2 \rangle= [8.9\pm1.1] \times 10^{-4}$ at 95 GHz and $[6.9\pm1.1] \times 10^{-4}$ at 150~GHz for the full sample. This is consistent with the values obtained for a sub-sample of active galactic nuclei. For dusty sources, we find 95 per cent upper limits of $\langle p^2 \rangle_{\rm 95}<16.9 \times 10^{-3}$ and $\langle p^2 \rangle_{\rm 150}<2.6 \times 10^{-3}$. We find no evidence that the polarisation fraction depends on the source flux or observing frequency. The 1-$\sigma$ upper limit on measured mean squared polarisation fraction at 150 GHz implies that extragalactic foregrounds will be subdominant to the CMB E and B mode polarisation power spectra out to at least $\ell\lesssim5700$ ($\ell\lesssim4700$) and $\ell\lesssim5300$ ($\ell\lesssim3600$), respectively at 95 (150) GHz.

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 a measurement of the cosmic microwave background (CMB) lensing potential using 500 deg$^2$ of 150 GHz data from the SPTpol receiver on the South Pole Telescope. The lensing potential is reconstructed with signal-to-noise per mode greater than unity at lensing multipoles $L \lesssim 250$, using a quadratic estimator on a combination of CMB temperature and polarization maps. We report measurements of the lensing potential power spectrum in the multipole range of $100< L < 2000$ from sets of temperature-only, polarization-only, and minimum-variance estimators. We measure the lensing amplitude by taking the ratio of the measured spectrum to the expected spectrum from the best-fit $\Lambda$CDM model to the $\textit{Planck}$ 2015 TT+lowP+lensing dataset. For the minimum-variance estimator, we find $A_{\rm{MV}} = 0.944 \pm 0.058{\rm (Stat.)}\pm0.025{\rm (Sys.)}$; restricting to only polarization data, we find $A_{\rm{POL}} = 0.906 \pm 0.090 {\rm (Stat.)} \pm 0.040 {\rm (Sys.)}$. Considering statistical uncertainties alone, this is the most precise polarization-only lensing amplitude constraint to date (10.1 $\sigma$), and is more precise than our temperature-only constraint. We perform null tests and consistency checks and find no evidence for significant contamination.

We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $\lambda \ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8\sigma$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements.

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

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

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

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

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

We use cosmic microwave background (CMB) temperature maps from the 500 deg$^{2}$ SPTpol survey to measure the stacked lensing convergence of galaxy clusters from the Dark Energy Survey (DES) Year-3 redMaPPer (RM) cluster catalog. The lensing signal is extracted through a modified quadratic estimator designed to be unbiased by the thermal Sunyaev-Zel{'}dovich (tSZ) effect. The modified estimator uses a tSZ-free map, constructed from the SPTpol 95 and 150 GHz datasets, to estimate the background CMB gradient. For lensing reconstruction, we employ two versions of the RM catalog: a flux-limited sample containing 4003 clusters and a volume-limited sample with 1741 clusters. We detect lensing at a significance of 8.7$\sigma$(6.7$\sigma$) with the flux(volume)-limited sample. By modeling the reconstructed convergence using the Navarro-Frenk-White profile, we find the average lensing masses to be $M_{200m}$ = ($1.62^{+0.32}_{-0.25}$ [stat.] $\pm$ 0.04 [sys.]) and ($1.28^{+0.14}_{-0.18}$ [stat.] $\pm$ 0.03 [sys.]) $\times\ 10^{14}\ M_{\odot}$ for the volume- and flux-limited samples respectively. The systematic error budget is much smaller than the statistical uncertainty and is dominated by the uncertainties in the RM cluster centroids. We use the volume-limited sample to calibrate the normalization of the mass-richness scaling relation, and find a result consistent with the galaxy weak-lensing measurements from DES (Mcclintock et al. 2018).

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

We present measurements of the $E$-mode polarization angular auto-power spectrum ($EE$) and temperature-$E$-mode cross-power spectrum ($TE$) of the cosmic microwave background (CMB) using 150 GHz data from three seasons of SPTpol observations. We report the power spectra over the spherical harmonic multipole range $50 < \ell \leq 8000$, and detect nine acoustic peaks in the $EE$ spectrum with high signal-to-noise ratio. These measurements are the most sensitive to date of the $EE$ and $TE$ power spectra at $\ell > 1050$ and $\ell > 1475$, respectively. The observations cover 500 deg$^2$, a fivefold increase in area compared to previous SPTpol analyses, which increases our sensitivity to the photon diffusion damping tail of the CMB power spectra enabling tighter constraints on \LCDM model extensions. After masking all sources with unpolarized flux $>50$ mJy we place a 95% confidence upper limit on residual polarized point-source power of $D_\ell = \ell(\ell+1)C_\ell/2\pi <0.107\,\mu{\rm K}^2$ at $\ell=3000$, suggesting that the $EE$ damping tail dominates foregrounds to at least $\ell = 4050$ with modest source masking. We find that the SPTpol dataset is in mild tension with the $\Lambda CDM$ model ($2.1\,\sigma$), and different data splits prefer parameter values that differ at the $\sim 1\,\sigma$ level. When fitting SPTpol data at $\ell < 1000$ we find cosmological parameter constraints consistent with those for $Planck$ temperature. Including SPTpol data at $\ell > 1000$ results in a preference for a higher value of the expansion rate ($H_0 = 71.3 \pm 2.1\,\mbox{km}\,s^{-1}\mbox{Mpc}^{-1}$ ) and a lower value for present-day density fluctuations ($\sigma_8 = 0.77 \pm 0.02$).

We present a demonstration of delensing the observed cosmic microwave background (CMB) B-mode polarization anisotropy. This process of reducing the gravitational-lensing generated B-mode component will become increasingly important for improving searches for the B modes produced by primordial gravitational waves. In this work, we delens B-mode maps constructed from multi-frequency SPTpol observations of a 90 deg$^2$ patch of sky by subtracting a B-mode template constructed from two inputs: SPTpol E-mode maps and a lensing potential map estimated from the $\textit{Herschel}$ $500\,\mu m$ map of the CIB. We find that our delensing procedure reduces the measured B-mode power spectrum by 28% in the multipole range $300 < \ell < 2300$; this is shown to be consistent with expectations from theory and simulations and to be robust against systematics. The null hypothesis of no delensing is rejected at $6.9 \sigma$. Furthermore, we build and use a suite of realistic simulations to study the general properties of the delensing process and find that the delensing efficiency achieved in this work is limited primarily by the noise in the lensing potential map. We demonstrate the importance of including realistic experimental non-idealities in the delensing forecasts used to inform instrument and survey-strategy planning of upcoming lower-noise experiments, such as CMB-S4.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

The Hydrogen Epoch of Reionization Array (HERA) is a radio interferometer aiming to detect the power spectrum of 21 cm fluctuations from neutral hydrogen from the Epoch of Reionization (EOR). Drawing on lessons from the Murchison Widefield Array (MWA) and the Precision Array for Probing the Epoch of Reionization (PAPER), HERA is a hexagonal array of large (14 m diameter) dishes with suspended dipole feeds. Not only does the dish determine overall sensitivity, it affects the observed frequency structure of foregrounds in the interferometer. This is the first of a series of four papers characterizing the frequency and angular response of the dish with simulations and measurements. We focus in this paper on the angular response (i.e., power pattern), which sets the relative weighting between sky regions of high and low delay, and thus, apparent source frequency structure. We measure the angular response at 137 MHz using the ORBCOMM beam mapping system of Neben et al. We measure a collecting area of 93 m^2 in the optimal dish/feed configuration, implying HERA-320 should detect the EOR power spectrum at z~9 with a signal-to-noise ratio of 12.7 using a foreground avoidance approach with a single season of observations, and 74.3 using a foreground subtraction approach. Lastly we study the impact of these beam measurements on the distribution of foregrounds in Fourier space.

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 time-domain electromagnetic simulations to determine the spectral characteristics of the Hydrogen Epoch of Reionization Arrays (HERA) antenna. These simulations are part of a multi-faceted campaign to determine the effectiveness of the dish's design for obtaining a detection of redshifted 21 cm emission from the epoch of reionization. Our simulations show the existence of reflections between HERA's suspended feed and its parabolic dish reflector that fall below -40 dB at 150 ns and, for reasonable impedance matches, have a negligible impact on HERA's ability to constrain EoR parameters. It follows that despite the reflections they introduce, dishes are effective for increasing the sensitivity of EoR experiments at relatively low cost. We find that electromagnetic resonances in the HERA feed's cylindrical skirt, which is intended to reduce cross coupling and beam ellipticity, introduces significant power at large delays ($-40$ dB at 200 ns) which can lead to some loss of measurable Fourier modes and a modest reduction in sensitivity. Even in the presence of this structure, we find that the spectral response of the antenna is sufficiently smooth for delay filtering to contain foreground emission at line-of-sight wave numbers below $k_\parallel \lesssim 0.2$ $h$Mpc$^{-1}$, in the region where the current PAPER experiment operates. Incorporating these results into a Fisher Matrix analysis, we find that the spectral structure observed in our simulations has only a small effect on the tight constraints HERA can achieve on parameters associated with the astrophysics of reionization.

The millimeter transient sky is largely unexplored, with measurements limited to follow-up of objects detected at other wavelengths. High-angular-resolution telescopes designed for measurement of the cosmic microwave background offer the possibility to discover new, unknown transient sources in this band, particularly the afterglows of unobserved gamma-ray bursts. Here we use the 10-meter millimeter-wave South Pole Telescope, designed for the primary purpose of observing the cosmic microwave background at arcminute and larger angular scales, to conduct a search for such objects. During the 2012-2013 season, the telescope was used to continuously observe a 100 square degree patch of sky centered at RA 23h30m and declination -55 degrees using the polarization-sensitive SPTpol camera in two bands centered at 95 and 150 GHz. These 6000 hours of observations provided continuous monitoring for day- to month-scale millimeter-wave transient sources at the 10 mJy level. One candidate object was observed with properties broadly consistent with a gamma-ray burst afterglow, but at a statistical significance too low (p=0.01) to confirm detection.

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

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status.

(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 South Pole Telescope (SPT) is a high-resolution microwave-frequency telescope designed to observe the Cosmic Microwave Background (CMB). To date, two cameras have been installed on the SPT to conduct two surveys of the CMB, the first in intensity only (SPT-SZ) and the second in intensity and polarization (SPTpol). A third-generation polarization-sensitive camera is currently in development (SPT-3G). This thesis describes work spanning all three instruments on the SPT. I present my work in time-reversed order, to follow the canonical narrative of instrument development, deployment, and analysis. First, the development and testing of novel 3-band multichroic Transition Edge Sensor (TES) bolometers for the SPT-3G experiment is detailed, followed by the development and deployment of the frequency multiplexed cryogenic readout electronics for the SPTpol experiment, and concluding with the analysis of data taken by the SPT-SZ instrument. I describe the development of a Bayesian likelihood based method I developed for measuring the integrated Comptonization (Ysz) of galaxy clusters from the Sunyaev-Zel'dovich (SZ) effect, and constraining galaxy cluster Ysz-mass scaling relations.

We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 \sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 \sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7\sigma$ deviation.

We describe the cryogenic half-wave plate rotation mechanisms built for and used in Spider, a polarization-sensitive balloon-borne telescope array that observed the Cosmic Microwave Background at 95 GHz and 150 GHz during a stratospheric balloon flight from Antarctica in January 2015. The mechanisms operate at liquid helium temperature in flight. A three-point contact design keeps the mechanical bearings relatively small but allows for a large (305 mm) diameter clear aperture. A worm gear driven by a cryogenic stepper motor allows for precise positioning and prevents undesired rotation when the motors are depowered. A custom-built optical encoder system monitors the bearing angle to an absolute accuracy of +/- 0.1 degrees. The system performed well in Spider during its successful 16 day flight.

We present a measurement of the cosmic microwave background (CMB) gravitational lensing potential using data from the first two seasons of observations with SPTpol, the polarization-sensitive receiver currently installed on the South Pole Telescope (SPT). The observations used in this work cover 100 deg$^2$ of sky with arcminute resolution at 150 GHz. Using a quadratic estimator, we make maps of the CMB lensing potential from combinations of CMB temperature and polarization maps. We combine these lensing potential maps to form a minimum-variance (MV) map. The lensing potential is measured with a signal-to-noise ratio of greater than one for angular multipoles between $100< L <250$. This is the highest signal-to-noise mass map made from the CMB to date and will be powerful in cross-correlation with other tracers of large-scale structure. We calculate the power spectrum of the lensing potential for each estimator, and we report the value of the MV power spectrum between $100< L <2000$ as our primary result. We constrain the ratio of the spectrum to a fiducial $\Lambda$CDM model to be $A_{\rm MV}=0.92 \pm 0.14 {\rm\, (Stat.)} \pm 0.08 {\rm\, (Sys.)}$. Restricting ourselves to polarized data only, we find $A_{\rm POL}=0.92 \pm 0.24 {\rm\, (Stat.)} \pm 0.11 {\rm\, (Sys.)}$. This measurement rejects the hypothesis of no lensing at $5.9 \sigma$ using polarization data alone, and at $14 \sigma$ using both temperature and polarization data.

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

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 a measurement of the $B$-mode polarization power spectrum (the $BB$ spectrum) from 100 $\mathrm{deg}^2$ of sky observed with SPTpol, a polarization-sensitive receiver currently installed on the South Pole Telescope. The observations used in this work were taken during 2012 and early 2013 and include data in spectral bands centered at 95 and 150 GHz. We report the $BB$ spectrum in five bins in multipole space, spanning the range $300 \le \ell \le 2300$, and for three spectral combinations: 95 GHz $\times$ 95 GHz, 95 GHz $\times$ 150 GHz, and 150 GHz $\times$ 150 GHz. We subtract small ($< 0.5 \sigma$ in units of statistical uncertainty) biases from these spectra and account for the uncertainty in those biases. The resulting power spectra are inconsistent with zero power but consistent with predictions for the $BB$ spectrum arising from the gravitational lensing of $E$-mode polarization. If we assume no other source of $BB$ power besides lensed $B$ modes, we determine a preference for lensed $B$ modes of $4.9 \sigma$. After marginalizing over tensor power and foregrounds, namely polarized emission from galactic dust and extragalactic sources, this significance is $4.3 \sigma$. Fitting for a single parameter, $A_\mathrm{lens}$, that multiplies the predicted lensed $B$-mode spectrum, and marginalizing over tensor power and foregrounds, we find $A_\mathrm{lens} = 1.08 \pm 0.26$, indicating that our measured spectra are consistent with the signal expected from gravitational lensing. The data presented here provide the best measurement to date of the $B$-mode power spectrum on these angular scales.

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 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 measurements of $E$-mode polarization and temperature-$E$-mode correlation in the cosmic microwave background (CMB) using data from the first season of observations with SPTpol, the polarization-sensitive receiver currently installed on the South Pole Telescope (SPT). The observations used in this work cover 100~\sqdeg\ of sky with arcminute resolution at $150\,$GHz. We report the $E$-mode angular auto-power spectrum ($EE$) and the temperature-$E$-mode angular cross-power spectrum ($TE$) over the multipole range $500 < \ell \leq5000$. These power spectra improve on previous measurements in the high-$\ell$ (small-scale) regime. We fit the combination of the SPTpol power spectra, data from \planck\, and previous SPT measurements with a six-parameter \LCDM cosmological model. We find that the best-fit parameters are consistent with previous results. The improvement in high-$\ell$ sensitivity over previous measurements leads to a significant improvement in the limit on polarized point-source power: after masking sources brighter than 50\,mJy in unpolarized flux at 150\,GHz, we find a 95\% confidence upper limit on unclustered point-source power in the $EE$ spectrum of $D_\ell = \ell (\ell+1) C_\ell / 2 \pi < 0.40 \ \mu{\mbox{K}}^2$ at $\ell=3000$, indicating that future $EE$ measurements will not be limited by power from unclustered point sources in the multipole range $\ell < 3600$, and possibly much higher in $\ell.$

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.

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.

We describe the design of a new polarization sensitive receiver, SPT-3G, for the 10-meter South Pole Telescope (SPT). The SPT-3G receiver will deliver a factor of ~20 improvement in mapping speed over the current receiver, SPTpol. The sensitivity of the SPT-3G receiver will enable the advance from statistical detection of B-mode polarization anisotropy power to high signal-to-noise measurements of the individual modes, i.e., maps. This will lead to precise (~0.06 eV) constraints on the sum of neutrino masses with the potential to directly address the neutrino mass hierarchy. It will allow a separation of the lensing and inflationary B-mode power spectra, improving constraints on the amplitude and shape of the primordial signal, either through SPT-3G data alone or in combination with BICEP-2/KECK, which is observing the same area of sky. The measurement of small-scale temperature anisotropy will provide new constraints on the epoch of reionization. Additional science from the SPT-3G survey will be significantly enhanced by the synergy with the ongoing optical Dark Energy Survey (DES), including: a 1% constraint on the bias of optical tracers of large-scale structure, a measurement of the differential Doppler signal from pairs of galaxy clusters that will test General Relativity on ~200 Mpc scales, and improved cosmological constraints from the abundance of clusters of galaxies.

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

Gravitational lensing of the cosmic microwave background generates a curl pattern in the observed polarization. This "B-mode" signal provides a measure of the projected mass distribution over the entire observable Universe and also acts as a contaminant for the measurement of primordial gravity-wave signals. In this Letter we present the first detection of gravitational lensing B modes, using first-season data from the polarization-sensitive receiver on the South Pole Telescope (SPTpol). We construct a template for the lensing B-mode signal by combining E-mode polarization measured by SPTpol with estimates of the lensing potential from a Herschel-SPIRE map of the cosmic infrared background. We compare this template to the B modes measured directly by SPTpol, finding a non-zero correlation at 7.7 sigma significance. The correlation has an amplitude and scale-dependence consistent with theoretical expectations, is robust with respect to analysis choices, and constitutes the first measurement of a powerful cosmological observable.

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

The SPTpol camera is a two-color, polarization-sensitive bolometer receiver, and was installed on the 10 meter South Pole Telescope in January 2012. SPTpol is designed to study the faint polarization signals in the Cosmic Microwave Background, with two primary scientific goals. One is to constrain the tensor-to-scalar ratio of perturbations in the primordial plasma, and thus constrain the space of permissible inflationary models. The other is to measure the weak lensing effect of large-scale structure on CMB polarization, which can be used to constrain the sum of neutrino masses as well as other growth-related parameters. The SPTpol focal plane consists of seven 84-element monolithic arrays of 150 GHz pixels (588 total) and 180 individual 90 GHz single-pixel modules. In this paper we present the design and characterization of the 90 GHz modules.

We present the software system used to control and operate the South Pole Telescope. The South Pole Telescope is a 10-meter millimeter-wavelength telescope designed to measure anisotropies in the cosmic microwave background (CMB) at arcminute angular resolution. In the austral summer of 2011/12, the SPT was equipped with a new polarization-sensitive camera, which consists of 1536 transition-edge sensor bolometers. The bolometers are read out using 36 independent digital frequency multiplexing (\dfmux) readout boards, each with its own embedded processors. These autonomous boards control and read out data from the focal plane with on-board software and firmware. An overall control software system running on a separate control computer controls the \dfmux boards, the cryostat and all other aspects of telescope operation. This control software collects and monitors data in real-time, and stores the data to disk for transfer to the United States for analysis.

In January 2012, the 10m South Pole Telescope (SPT) was equipped with a polarization-sensitive camera, SPTpol, in order to measure the polarization anisotropy of the cosmic microwave background (CMB). Measurements of the polarization of the CMB at small angular scales (~several arcminutes) can detect the gravitational lensing of the CMB by large scale structure and constrain the sum of the neutrino masses. At large angular scales (~few degrees) CMB measurements can constrain the energy scale of Inflation. SPTpol is a two-color mm-wave camera that consists of 180 polarimeters at 90 GHz and 588 polarimeters at 150 GHz, with each polarimeter consisting of a dual transition edge sensor (TES) bolometers. The full complement of 150 GHz detectors consists of 7 arrays of 84 ortho-mode transducers (OMTs) that are stripline coupled to two TES detectors per OMT, developed by the TRUCE collaboration and fabricated at NIST. Each 90 GHz pixel consists of two antenna-coupled absorbers coupled to two TES detectors, developed with Argonne National Labs. The 1536 total detectors are read out with digital frequency-domain multiplexing (DfMUX). The SPTpol deployment represents the first on-sky tests of both of these detector technologies, and is one of the first deployed instruments using DfMUX readout technology. We present the details of the design, commissioning, deployment, on-sky optical characterization and detector performance of the complete SPTpol focal plane.

SPTpol is a dual-frequency polarization-sensitive camera that was deployed on the 10-meter South Pole Telescope in January 2012. SPTpol will measure the polarization anisotropy of the cosmic microwave background (CMB) on angular scales spanning an arcminute to several degrees. The polarization sensitivity of SPTpol will enable a detection of the CMB "B-mode" polarization from the detection of the gravitational lensing of the CMB by large scale structure, and a detection or improved upper limit on a primordial signal due to inflationary gravity waves. The two measurements can be used to constrain the sum of the neutrino masses and the energy scale of inflation. These science goals can be achieved through the polarization sensitivity of the SPTpol camera and careful control of systematics. The SPTpol camera consists of 768 pixels, each containing two transition-edge sensor (TES) bolometers coupled to orthogonal polarizations, and a total of 1536 bolometers. The pixels are sensitive to light in one of two frequency bands centered at 90 and 150 GHz, with 180 pixels at 90 GHz and 588 pixels at 150 GHz. The SPTpol design has several features designed to control polarization systematics, including: single-moded feedhorns with low cross-polarization, bolometer pairs well-matched to difference atmospheric signals, an improved ground shield design based on far-sidelobe measurements of the SPT, and a small beam to reduce temperature to polarization leakage. We present an overview of the SPTpol instrument design, project status, and science projections.

The SPTpol camera is a dichroic polarimetric receiver at 90 and 150 GHz. Deployed in January 2012 on the South Pole Telescope (SPT), SPTpol is looking for faint polarization signals in the Cosmic Microwave Background (CMB). The camera consists of 180 individual Transition Edge Sensor (TES) polarimeters at 90 GHz and seven 84-polarimeter camera modules (a total of 588 polarimeters) at 150 GHz. We present the design, dark characterization, and in-lab optical properties of the 150 GHz camera modules. The modules consist of photolithographed arrays of TES polarimeters coupled to silicon platelet arrays of corrugated feedhorns, both of which are fabricated at NIST-Boulder. In addition to mounting hardware and RF shielding, each module also contains a set of passive readout electronics for digital frequency-domain multiplexing. A single module, therefore, is fully functional as a miniature focal plane and can be tested independently. Across the modules tested before deployment, the detectors average a critical temperature of 478 mK, normal resistance R_N of 1.2 Ohm, unloaded saturation power of 22.5 pW, (detector-only) optical efficiency of ~ 90%, and have electrothermal time constants < 1 ms in transition.

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

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

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

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

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

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.