We present a catalog of emissive point sources detected in the SPT-SZ survey, a contiguous 2530-square-degree area surveyed with the South Pole Telescope (SPT) from 2008 - 2011 in three bands centered at 95, 150, and 220 GHz. The catalog contains 4845 sources measured at a significance of 4.5 sigma or greater in at least one band, corresponding to detections above approximately 9.8, 5.8, and 20.4 mJy in 95, 150, and 220 GHz, respectively. Spectral behavior in the SPT bands is used for source classification into two populations based on the underlying physical mechanisms of compact, emissive sources that are bright at millimeter wavelengths: synchrotron radiation from active galactic nuclei and thermal emission from dust. The latter population includes a component of high-redshift sources often referred to as submillimeter galaxies (SMGs). In the relatively bright flux ranges probed by the survey, these sources are expected to be magnified by strong gravitational lensing. The survey also contains sources consistent with protoclusters, groups of dusty galaxies at high redshift undergoing collapse. We cross-match the SPT-SZ catalog with external catalogs at radio, infrared, and X-ray wavelengths and identify available redshift information. The catalog splits into 3980 synchrotron-dominated and 865 dust-dominated sources and we determine a list of 506 SMGs. Ten sources in the catalog are identified as stars. We calculate number counts for the full catalog, and synchrotron and dusty components, using a bootstrap method and compare our measured counts with models. This paper represents the third and final catalog of point sources in the SPT-SZ survey.

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

Precision measurements of cosmic microwave background (CMB) polarization require extreme control of instrumental systematics. In a companion paper we have presented cosmological constraints from observations with the BICEP2 and Keck Array experiments up to and including the 2015 observing season (BK15), resulting in the deepest CMB polarization maps to date and a statistical sensitivity to the tensor-to-scalar ratio of $\sigma(r) = 0.020$. In this work we characterize the beams and constrain potential systematic contamination from main beam shape mismatch at the three BK15 frequencies (95, 150, and 220 GHz). Far-field maps of 7,360 distinct beam patterns taken from 2010-2015 are used to measure differential beam parameters and predict the contribution of temperature-to-polarization leakage to the BK15 B-mode maps. In the multifrequency, multicomponent likelihood analysis that uses BK15, Planck, and WMAP maps to separate sky components, we find that adding this predicted leakage to simulations induces a bias of $\Delta r = 0.0027 \pm 0.0019$. Future results using higher-quality beam maps and improved techniques to detect such leakage in CMB data will substantially reduce this uncertainty, enabling the levels of systematics control needed for BICEP Array and other experiments that plan to definitively probe large-field inflation.

We present an internal consistency test of South Pole Telescope (SPT) measurements of the cosmic microwave background (CMB) temperature anisotropy using three-band data from the SPT-SZ survey. These measurements are made from observations of ~2500 deg^2 of sky in three frequency bands centered at 95, 150, and 220 GHz. We combine the information from these three bands into six semi-independent estimates of the CMB power spectrum (three single-frequency power spectra and three cross-frequency spectra) over the multipole range 650 < l < 3000. We subtract an estimate of foreground power from each power spectrum and evaluate the consistency among the resulting CMB-only spectra. We determine that the six foreground-cleaned power spectra are consistent with the null hypothesis, in which the six cleaned spectra contain only CMB power and noise. A fit of the data to this model results in a chi-squared value of 236.3 for 235 degrees of freedom, and the probability to exceed this chi-squared value is 46%.

Correlations between tracers of the matter density field and gravitational lensing are sensitive to the evolution of the matter power spectrum and the expansion rate across cosmic time. Appropriately defined ratios of such correlation functions, on the other hand, depend only on the angular diameter distances to the tracer objects and to the gravitational lensing source planes. Because of their simple cosmological dependence, such ratios can exploit available signal-to-noise down to small angular scales, even where directly modeling the correlation functions is difficult. We present a measurement of lensing ratios using galaxy position and lensing data from the Dark Energy Survey, and CMB lensing data from the South Pole Telescope and Planck, obtaining the highest precision lensing ratio measurements to date. Relative to the concordance $\Lambda$CDM model, we find a best fit lensing ratio amplitude of $A = 1.1 \pm 0.1$. We use the ratio measurements to generate cosmological constraints, focusing on the curvature parameter. We demonstrate that photometrically selected galaxies can be used to measure lensing ratios, and argue that future lensing ratio measurements with data from a combination of LSST and Stage-4 CMB experiments can be used to place interesting cosmological constraints, even after considering the systematic uncertainties associated with photometric redshift and galaxy shear estimation.

We present measurements of the power spectra of cosmic infrared background (CIB) and cosmic microwave background (CMB) fluctuations in six frequency bands. Maps at the lower three frequency bands, 95, 150, and 220 GHz (3330, 2000, 1360 $\mu$m) are from the South Pole Telescope, while the upper three frequency bands, 600, 857, and 1200 GHz (500, 350, 250 $\mu$m) are observed with Herschel/SPIRE. From these data, we produce 21 angular power spectra (six auto- and fifteen cross-frequency) spanning the multipole range $600 \le \ell \le 11,000$. Our measurements are the first to cross-correlate measurements near the peak of the CIB spectrum with maps at 95 GHz, complementing and extending the measurements from Planck Collaboration et al. (2014) at 218, 550, and 857 GHz. The observed fluctuations originate largely from clustered, infrared-emitting, dusty star-forming galaxies, the CMB, and to a lesser extent radio galaxies, active galactic nuclei, and the Sunyaev-Zel'dovich effect.

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

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

We present three maps of the millimeter-wave sky created by combining data from the South Pole Telescope (SPT) and the Planck satellite. We use data from the SPT-SZ survey, a survey of 2540 deg$^2$ of the the sky with arcminute resolution in three bands centered at 95, 150, and 220 GHz, and the full-mission Planck temperature data in the 100, 143, and 217 GHz bands. A linear combination of the SPT-SZ and Planck data is computed in spherical harmonic space, with weights derived from the noise of both instruments. This weighting scheme results in Planck data providing most of the large-angular-scale information in the combined maps, with the smaller-scale information coming from SPT-SZ data. A number of tests have been done on the maps. We find their angular power spectra to agree very well with theoretically predicted spectra and previously published results.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

We measure the cross-correlation between redMaGiC galaxies selected from the Dark Energy Survey (DES) Year-1 data and gravitational lensing of the cosmic microwave background (CMB) reconstructed from South Pole Telescope (SPT) and Planck data over 1289 sq. deg. When combining measurements across multiple galaxy redshift bins spanning the redshift range of $0.15<z<0.90$, we reject the hypothesis of no correlation at 19.9$\sigma$ significance. When removing small-scale data points where thermal Sunyaev-Zel'dovich signal and nonlinear galaxy bias could potentially bias our results, the detection significance is reduced to 9.9$\sigma$. We perform a joint analysis of galaxy-CMB lensing cross-correlations and galaxy clustering to constrain cosmology, finding $\Omega_{\rm m} = 0.276^{+0.029}_{-0.030}$ and $S_{8}=\sigma_{8}\sqrt{\mathstrut \Omega_{\rm m}/0.3} = 0.800^{+0.090}_{-0.094}$. We also perform two alternate analyses aimed at constraining only the growth rate of cosmic structure as a function of redshift, finding consistency with predictions from the concordance $\Lambda$CDM model. The measurements presented here are part of a joint cosmological analysis that combines galaxy clustering, galaxy lensing and CMB lensing using data from DES, SPT and Planck.

We cross-correlate galaxy weak lensing measurements from the Dark Energy Survey (DES) year-one (Y1) data with a cosmic microwave background (CMB) weak lensing map derived from South Pole Telescope (SPT) and Planck data, with an effective overlapping area of 1289 deg$^{2}$. With the combined measurements from four source galaxy redshift bins, we reject the hypothesis of no lensing with a significance of $10.8\sigma$. When employing angular scale cuts, this significance is reduced to $6.8\sigma$, which remains the highest signal-to-noise measurement of its kind to date. We fit the amplitude of the correlation functions while fixing the cosmological parameters to a fiducial $\Lambda$CDM model, finding $A = 0.99 \pm 0.17$. We additionally use the correlation function measurements to constrain shear calibration bias, obtaining constraints that are consistent with previous DES analyses. Finally, when performing a cosmological analysis under the $\Lambda$CDM model, we obtain the marginalized constraints of $\Omega_{\rm m}=0.261^{+0.070}_{-0.051}$ and $S_{8}\equiv \sigma_{8}\sqrt{\Omega_{\rm m}/0.3} = 0.660^{+0.085}_{-0.100}$. These measurements are used in a companion work that presents cosmological constraints from the joint analysis of two-point functions among galaxies, galaxy shears, and CMB lensing using DES, SPT and Planck data.

We perform a joint analysis of the auto and cross-correlations between three cosmic fields: the galaxy density field, the galaxy weak lensing shear field, and the cosmic microwave background (CMB) weak lensing convergence field. These three fields are measured using roughly 1300 sq. deg. of overlapping optical imaging data from first year observations of the Dark Energy Survey and millimeter-wave observations of the CMB from both the South Pole Telescope Sunyaev-Zel'dovich survey and Planck. We present cosmological constraints from the joint analysis of the two-point correlation functions between galaxy density and galaxy shear with CMB lensing. We test for consistency between these measurements and the DES-only two-point function measurements, finding no evidence for inconsistency in the context of flat $\Lambda$CDM cosmological models. Performing a joint analysis of five of the possible correlation functions between these fields (excluding only the CMB lensing autospectrum) yields $S_{8}\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.782^{+0.019}_{-0.025}$ and $\Omega_{\rm m}=0.260^{+0.029}_{-0.019}$. We test for consistency between these five correlation function measurements and the Planck-only measurement of the CMB lensing autospectrum, again finding no evidence for inconsistency in the context of flat $\Lambda$CDM models. Combining constraints from all six two-point functions yields $S_{8}=0.776^{+0.014}_{-0.021}$ and $\Omega_{\rm m}= 0.271^{+0.022}_{-0.016}$. These results provide a powerful test and confirmation of the results from the first year DES joint-probes analysis.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

We study the consistency of 150 GHz data from the South Pole Telescope (SPT) and 143 GHz data from the Planck satellite over the patch of sky covered by the SPT-SZ survey. We first visually compare the maps and find that the residuals appear consistent with noise after accounting for differences in angular resolution and filtering. We then calculate (1) the cross-spectrum between two independent halves of SPT data, (2) the cross-spectrum between two independent halves of Planck data, and (3) the cross-spectrum between SPT and Planck data. We find the three cross-spectra are well-fit (PTE = 0.30) by the null hypothesis in which both experiments have measured the same sky map up to a single free calibration parameter---i.e., we find no evidence for systematic errors in either data set. As a by-product, we improve the precision of the SPT calibration by nearly an order of magnitude, from 2.6% to 0.3% in power. Finally, we compare all three cross-spectra to the full-sky Planck power spectrum and find marginal evidence for differences between the power spectra from the SPT-SZ footprint and the full sky. We model these differences as a power law in spherical harmonic multipole number. The best-fit value of this tilt is consistent among the three cross-spectra in the SPT-SZ footprint, implying that the source of this tilt is a sample variance fluctuation in the SPT-SZ region relative to the full sky. The consistency of cosmological parameters derived from these datasets is discussed in a companion paper.

We report constraints on cosmological parameters from the angular power spectrum of a cosmic microwave background (CMB) gravitational lensing potential map created using temperature data from 2500 deg$^2$ of South Pole Telescope (SPT) data supplemented with data from Planck in the same sky region, with the statistical power in the combined map primarily from the SPT data. We fit the corresponding lensing angular power spectrum to a model including cold dark matter and a cosmological constant ($\Lambda$CDM), and to models with single-parameter extensions to $\Lambda$CDM. We find constraints that are comparable to and consistent with constraints found using the full-sky Planck CMB lensing data. Specifically, we find $\sigma_8 \Omega_{\rm m}^{0.25}=0.598 \pm 0.024$ from the lensing data alone with relatively weak priors placed on the other $\Lambda$CDM parameters. In combination with primary CMB data from Planck, we explore single-parameter extensions to the $\Lambda$CDM model. We find $\Omega_k = -0.012^{+0.021}_{-0.023}$ or $M_{\nu}< 0.70$eV both at 95% confidence, all in good agreement with results that include the lensing potential as measured by Planck over the full sky. We include two independent free parameters that scale the effect of lensing on the CMB: $A_{L}$, which scales the lensing power spectrum in both the lens reconstruction power and in the smearing of the acoustic peaks, and $A^{\phi \phi}$, which scales only the amplitude of the CMB lensing reconstruction power spectrum. We find $A^{\phi \phi} \times A_{L} =1.01 \pm 0.08$ for the lensing map made from combined SPT and Planck temperature data, indicating that the amount of lensing is in excellent agreement with what is expected from the observed CMB angular power spectrum when not including the information from smearing of the acoustic peaks.

The Planck cosmic microwave background (CMB) temperature data are best fit with a LCDM model that is in mild tension with constraints from other cosmological probes. The South Pole Telescope (SPT) 2540 $\text{deg}^2$ SPT-SZ survey offers measurements on sub-degree angular scales (multipoles $650 \leq \ell \leq 2500$) with sufficient precision to use as an independent check of the Planck data. Here we build on the recent joint analysis of the SPT-SZ and Planck data in \citet{hou17} by comparing LCDM parameter estimates using the temperature power spectrum from both data sets in the SPT-SZ survey region. We also restrict the multipole range used in parameter fitting to focus on modes measured well by both SPT and Planck, thereby greatly reducing sample variance as a driver of parameter differences and creating a stringent test for systematic errors. We find no evidence of systematic errors from such tests. When we expand the maximum multipole of SPT data used, we see low-significance shifts in the angular scale of the sound horizon and the physical baryon and cold dark matter densities, with a resulting trend to higher Hubble constant. When we compare SPT and Planck data on the SPT-SZ sky patch to Planck full-sky data but keep the multipole range restricted, we find differences in the parameters $n_s$ and $A_se^{-2\tau}$. We perform further checks, investigating instrumental effects and modeling assumptions, and we find no evidence that the effects investigated are responsible for any of the parameter shifts. Taken together, these tests reveal no evidence for systematic errors in SPT or Planck data in the overlapping sky coverage and multipole range and, at most, weak evidence for a breakdown of LCDM or systematic errors influencing either the Planck data outside the SPT-SZ survey area or the SPT data at $\ell >2000$.

We present a cosmic microwave background (CMB) lensing map produced from a linear combination of South Pole Telescope (SPT) and \emph{Planck} temperature data. The 150 GHz temperature data from the $2500\ {\rm deg}^{2}$ SPT-SZ survey is combined with the \emph{Planck} 143 GHz data in harmonic space, to obtain a temperature map that has a broader $\ell$ coverage and less noise than either individual map. Using a quadratic estimator technique on this combined temperature map, we produce a map of the gravitational lensing potential projected along the line of sight. We measure the auto-spectrum of the lensing potential $C_{L}^{\phi\phi}$, and compare it to the theoretical prediction for a $\Lambda$CDM cosmology consistent with the \emph{Planck} 2015 data set, finding a best-fit amplitude of $0.95_{-0.06}^{+0.06}({\rm Stat.})\! _{-0.01}^{+0.01}({\rm Sys.})$. The null hypothesis of no lensing is rejected at a significance of $24\,\sigma$. One important use of such a lensing potential map is in cross-correlations with other dark matter tracers. We demonstrate this cross-correlation in practice by calculating the cross-spectrum, $C_{L}^{\phi G}$, between the SPT+\emph{Planck} lensing map and Wide-field Infrared Survey Explorer (\emph{WISE}) galaxies. We fit $C_{L}^{\phi G}$ to a power law of the form $p_{L}=a(L/L_{0})^{-b}$ with $a=2.15 \times 10^{-8}$, $b=1.35$, $L_{0}=490$, and find $\eta^{\phi G}=0.94^{+0.04}_{-0.04}$, which is marginally lower, but in good agreement with $\eta^{\phi G}=1.00^{+0.02}_{-0.01}$, the best-fit amplitude for the cross-correlation of \emph{Planck}-2015 CMB lensing and \emph{WISE} galaxies over $\sim67\%$ of the sky. The lensing potential map presented here will be used for cross-correlation studies with the Dark Energy Survey (DES), whose footprint nearly completely covers the SPT $2500\ {\rm deg}^2$ field.

We present maps of the Large and Small Magellanic Clouds from combined South Pole Telescope (SPT) and Planck data. The Planck satellite observes in nine bands, while the SPT data used in this work were taken with the three-band SPT-SZ camera, The SPT-SZ bands correspond closely to three of the nine Planck bands, namely those centered at 1.4, 2.1, and 3.0 mm. The angular resolution of the Planck data ranges from 5 to 10 arcmin, while the SPT resolution ranges from 1.0 to 1.7 arcmin. The combined maps take advantage of the high resolution of the SPT data and the long-timescale stability of the space-based Planck observations to deliver robust brightness measurements on scales from the size of the maps down to ~1 arcmin. In each band, we first calibrate and color-correct the SPT data to match the Planck data, then we use noise estimates from each instrument and knowledge of each instrument's beam to make the inverse-variance-weighted combination of the two instruments' data as a function of angular scale. We create maps assuming a range of underlying emission spectra and at a range of final resolutions. We perform several consistency tests on the combined maps and estimate the expected noise in measurements of features in the maps. We compare maps from this work to maps from the Herschel HERITAGE survey, finding general consistency between the datasets. All data products described in this paper are available for download from the NASA Legacy Archive for Microwave Background Data Analysis server.

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.

BICEP3, the latest telescope in the BICEP/Keck program, started science observations in March 2016. It is a 550mm aperture refractive telescope observing the polarization of the cosmic microwave background at 95 GHz. We show the focal plane design and detector performance, including spectral response, optical efficiency and preliminary sensitivity of the upgraded BICEP3. We demonstrate 9.72$\mu$K$\sqrt{\textrm{s}}$ noise performance of the BICEP3 receiver.

BICEP3 is a small-aperture refracting cosmic microwave background (CMB) telescope designed to make sensitive polarization maps in pursuit of a potential B-mode signal from inflationary gravitational waves. It is the latest in the BICEP/Keck Array series of CMB experiments at the South Pole, which has provided the most stringent constraints on inflation to date. For the 2016 observing season, BICEP3 was outfitted with a full suite of 2400 optically coupled detectors operating at 95 GHz. In these proceedings we report on the far field beam performance using calibration data taken during the 2015-2016 summer deployment season in situ with a thermal chopped source. We generate high-fidelity per-detector beam maps, show the array-averaged beam profile, and characterize the differential beam response between co-located, orthogonally polarized detectors which contributes to the leading instrumental systematic in pair differencing experiments. We find that the levels of differential pointing, beamwidth, and ellipticity are similar to or lower than those measured for BICEP2 and Keck Array. The magnitude and distribution of BICEP3's differential beam mismatch - and the level to which temperature-to-polarization leakage may be marginalized over or subtracted in analysis - will inform the design of next-generation CMB experiments with many thousands of detectors.

BICEP3 is a 520 mm aperture, compact two-lens refractor designed to observe the polarization of the cosmic microwave background (CMB) at 95 GHz. Its focal plane consists of modularized tiles of antenna-coupled transition edge sensors (TESs), similar to those used in BICEP2 and the Keck Array. The increased per-receiver optical throughput compared to BICEP2/Keck Array, due to both its faster f/1.7 optics and the larger aperture, more than doubles the combined mapping speed of the BICEP/Keck program. The BICEP3 receiver was recently upgraded to a full complement of 20 tiles of detectors (2560 TESs) and is now beginning its second year of observation (and first science season) at the South Pole. We report on its current performance and observing plans. Given its high per-receiver throughput while maintaining the advantages of a compact design, BICEP3-class receivers are ideally suited as building blocks for a 3rd-generation CMB experiment, consisting of multiple receivers spanning 35 GHz to 270 GHz with total detector count in the tens of thousands. We present plans for such an array, the new "BICEP Array" that will replace the Keck Array at the South Pole, including design optimization, frequency coverage, and deployment/observing strategies.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

Mega-pixel charge-integrating detectors are common in near-IR imaging applications. Optimal signal-to-noise ratio estimates of the photocurrents, which are particularly important in the low-signal regime, are produced by fitting linear models to sequential reads of the charge on the detector. Algorithms that solve this problem have a long history, but can be computationally intensive. Furthermore, the cosmic ray background is appreciable for these detectors in Earth orbit, particularly above the Earth's magnetic poles and the South Atlantic Anomaly, and on-board reduction routines must be capable of flagging affected pixels. In this paper we present an algorithm that generates optimal photocurrent estimates and flags random transient charge generation from cosmic rays, and is specifically designed to fit on a computationally restricted platform. We take as a case study the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer (SPHEREx), a NASA Small Explorer astrophysics experiment concept, and show that the algorithm can easily fit in the resource-constrained environment of such a restricted platform. Detailed simulations of the input astrophysical signals and detector array performance are used to characterize the fitting routines in the presence of complex noise properties and charge transients. We use both Hubble Space Telescope Wide Field Camera-3 and Wide-field Infrared Survey Explorer to develop an empirical understanding of the susceptibility of near-IR detectors in low earth orbit and build a model for realistic cosmic ray energy spectra and rates. We show that our algorithm generates an unbiased estimate of the true photocurrent that is identical to that from a standard line fitting package, and characterize the rate, energy, and timing of both detected and undetected transient events.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

(abridged) We present cosmological constraints obtained from galaxy clusters identified by their Sunyaev-Zel'dovich effect signature in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich survey. We consider the 377 cluster candidates identified at z>0.25 with a detection significance greater than five, corresponding to the 95% purity threshold for the survey. We compute constraints on cosmological models using the measured cluster abundance as a function of mass and redshift. We include additional constraints from multi-wavelength observations, including Chandra X-ray data for 82 clusters and a weak lensing-based prior on the normalization of the mass-observable scaling relations. Assuming a LCDM cosmology, where the species-summed neutrino mass has the minimum allowed value (mnu = 0.06 eV) from neutrino oscillation experiments, we combine the cluster data with a prior on H0 and find sigma_8 = 0.797+-0.031 and Omega_m = 0.289+-0.042, with the parameter combination sigma_8(Omega_m/0.27)^0.3 = 0.784+-0.039. These results are in good agreement with constraints from the CMB from SPT, WMAP, and Planck, as well as with constraints from other cluster datasets. Adding mnu as a free parameter, we find mnu = 0.14+-0.08 eV when combining the SPT cluster data with Planck CMB data and BAO data, consistent with the minimum allowed value. Finally, we consider a cosmology where mnu and N_eff are fixed to the LCDM values, but the dark energy equation of state parameter w is free. Using the SPT cluster data in combination with an H0 prior, we measure w = -1.28+-0.31, a constraint consistent with the LCDM cosmological model and derived from the combination of growth of structure and geometry. When combined with primarily geometrical constraints from Planck CMB, H0, BAO and SNe, adding the SPT cluster data improves the w constraint from the geometrical data alone by 14%, to w = -1.023+-0.042.

We present results from an analysis of all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season. This includes the first Keck Array observations at 95 GHz. The maps reach a depth of 50 nK deg in Stokes $Q$ and $U$ in the 150 GHz band and 127 nK deg in the 95 GHz band. We take auto- and cross-spectra between these maps and publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353 GHz. An excess over lensed-LCDM is detected at modest significance in the 95x150 $BB$ spectrum, and is consistent with the dust contribution expected from our previous work. No significant evidence for synchrotron emission is found in spectra such as 23x95, or for correlation between the dust and synchrotron sky patterns in spectra such as 23x353. We take the likelihood of all the spectra for a multi-component model including lensed-LCDM, dust, synchrotron and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio $r$), using priors on the frequency spectral behaviors of dust and synchrotron emission from previous analyses of WMAP and Planck data in other regions of the sky. This analysis yields an upper limit $r_{0.05}<0.09$ at 95% confidence, which is robust to variations explored in analysis and priors. Combining these $B$-mode results with the (more model-dependent) constraints from Planck analysis of CMB temperature plus BAO and other data, yields a combined limit $r_{0.05}<0.07$ at 95% confidence. These are the strongest constraints to date on inflationary gravitational waves.

BICEP3 is a $550~mm$ aperture telescope with cold, on-axis, refractive optics designed to observe at the $95~GHz$ band from the South Pole. It is the newest member of the BICEP/Keck family of inflationary probes specifically designed to measure the polarization of the cosmic microwave background (CMB) at degree-angular scales. BICEP3 is designed to house 1280 dual-polarization pixels, which, when fully-populated, totals to $\sim$9$\times$ the number of pixels in a single Keck $95~GHz$ receiver, thus further advancing the BICEP/Keck program's $95~GHz$ mapping speed. BICEP3 was deployed during the austral summer of 2014-2015 with 9 detector tiles, to be increased to its full capacity of 20 in the second season. After instrument characterization measurements were taken, CMB observation commenced in April 2015. Together with multi-frequency observation data from Planck, BICEP2, and the Keck Array, BICEP3 is projected to set upper limits on the tensor-to-scalar ratio to $r$ $\lesssim 0.03$ at $95\%$ C.L..

In a companion paper we have reported a $>5\sigma$ detection of degree scale $B $-mode polarization at 150 GHz by the BICEP2 experiment. Here we provide a detailed study of potential instrumental systematic contamination to that measurement. We focus extensively on spurious polarization that can potentially arise from beam imperfections. We present a heuristic classification of beam imperfections according to their symmetries and uniformities, and discuss how resulting contamination adds or cancels in maps that combine observations made at multiple orientations of the telescope about its boresight axis. We introduce a technique, which we call "deprojection", for filtering the leading order beam-induced contamination from time ordered data, and show that it removes power from BICEP2's $BB$ spectrum consistent with predictions using high signal-to-noise beam shape measurements. We detail the simulation pipeline that we use to directly simulate instrumental systematics and the calibration data used as input to that pipeline. Finally, we present the constraints on $BB$ contamination from individual sources of potential systematics. We find that systematics contribute $BB$ power that is a factor $\sim10\times$ below BICEP2's 3-year statistical uncertainty, and negligible compared to the observed $BB$ signal. The contribution to the best-fit tensor/scalar ratio is at a level equivalent to $r=(3-6)\times10^{-3}$.

The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.

BICEP2 and the Keck Array are polarization-sensitive microwave telescopes that observe the cosmic microwave background (CMB) from the South Pole at degree angular scales in search of a signature of inflation imprinted as B-mode polarization in the CMB. BICEP2 was deployed in late 2009, observed for three years until the end of 2012 at 150 GHz with 512 antenna-coupled transition edge sensor bolometers, and has reported a detection of B-mode polarization on degree angular scales. The Keck Array was first deployed in late 2010 and will observe through 2016 with five receivers at several frequencies (95, 150, and 220 GHz). BICEP2 and the Keck Array share a common optical design and employ the field-proven BICEP1 strategy of using small-aperture, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. Here we describe the optical design of both instruments and report a full characterization of the optical performance and beams of BICEP2 and the Keck Array at 150 GHz.

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 report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.

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 have developed antenna-coupled transition-edge sensor (TES) bolometers for a wide range of cosmic microwave background (CMB) polarimetry experiments, including BICEP2, Keck Array, and the balloon borne SPIDER. These detectors have reached maturity and this paper reports on their design principles, overall performance, and key challenges associated with design and production. Our detector arrays repeatedly produce spectral bands with 20%-30% bandwidth at 95, 150, or 220~GHz. The integrated antenna arrays synthesize symmetric co-aligned beams with controlled side-lobe levels. Cross-polarized response on boresight is typically ~0.5%, consistent with cross-talk in our multiplexed readout system. End-to-end optical efficiencies in our cameras are routinely 35% or higher, with per detector sensitivities of NET~300 uKrts. Thanks to the scalability of this design, we have deployed 2560 detectors as 1280 matched pairs in Keck Array with a combined instantaneous sensitivity of ~9 uKrts, as measured directly from CMB maps in the 2013 season. Similar arrays have recently flown in the SPIDER instrument, and development of this technology is ongoing.

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 secondary cosmic microwave background (CMB) anisotropies and cosmic infrared background (CIB) fluctuations using data from the South Pole Telescope (SPT) covering the complete 2540 sq.deg. SPT-SZ survey area. Data in the three SPT-SZ frequency bands centered at 95, 150, and 220 GHz, are used to produce six angular power spectra (three single-frequency auto-spectra and three cross-spectra) covering the multipole range 2000 < ell < 11000 (angular scales 5' > \theta > 1'). These are the most precise measurements of the angular power spectra at ell > 2500 at these frequencies. The main contributors to the power spectra at these angular scales and frequencies are the primary CMB, CIB, thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ), and radio galaxies. We include a constraint on the tSZ power from a measurement of the tSZ bispectrum from 800 sq.deg. of the SPT-SZ survey. We measure the tSZ power at 143 GHz to be DtSZ = 4.08 +0.58 -0.67 \mu K^2 and the kSZ power to be DkSZ = 2.9 +- 1.3 \mu K^2. The data prefer positive kSZ power at 98.1% CL. We measure a correlation coefficient of \xi = 0.113 +0.057 -0.054 between sources of tSZ and CIB power, with \xi < 0 disfavored at a confidence level of 99.0%. The constraint on kSZ power can be interpreted as an upper limit on the duration of reionization. When the post-reionization homogeneous kSZ signal is accounted for, we find an upper limit on the duration \Delta z < 5.4 at 95% CL.

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 report on the design and performance of the BICEP2 instrument and on its three-year data set. BICEP2 was designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 1 to 5 degrees ($\ell$=40-200), near the expected peak of the B-mode polarization signature of primordial gravitational waves from cosmic inflation. Measuring B-modes requires dramatic improvements in sensitivity combined with exquisite control of systematics. The BICEP2 telescope observed from the South Pole with a 26~cm aperture and cold, on-axis, refractive optics. BICEP2 also adopted a new detector design in which beam-defining slot antenna arrays couple to transition-edge sensor (TES) bolometers, all fabricated on a common substrate. The antenna-coupled TES detectors supported scalable fabrication and multiplexed readout that allowed BICEP2 to achieve a high detector count of 500 bolometers at 150 GHz, giving unprecedented sensitivity to B-modes at degree angular scales. After optimization of detector and readout parameters, BICEP2 achieved an instrument noise-equivalent temperature of 15.8 $\mu$K sqrt(s). The full data set reached Stokes Q and U map depths of 87.2 nK in square-degree pixels (5.2 $\mu$K arcmin) over an effective area of 384 square degrees within a 1000 square degree field. These are the deepest CMB polarization maps at degree angular scales to date. The power spectrum analysis presented in a companion paper has resulted in a significant detection of B-mode polarization at degree scales.

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

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

(abridged for arXiv) We report results from the BICEP2 experiment, a cosmic microwave background (CMB) polarimeter specifically designed to search for the signal of inflationary gravitational waves in the B-mode power spectrum around $\ell\sim80$. The telescope comprised a 26 cm aperture all-cold refracting optical system equipped with a focal plane of 512 antenna coupled transition edge sensor 150 GHz bolometers each with temperature sensitivity of $\approx300\mu\mathrm{K}_\mathrm{CMB}\sqrt{s}$. BICEP2 observed from the South Pole for three seasons from 2010 to 2012. A low-foreground region of sky with an effective area of 380 square deg was observed to a depth of 87 nK deg in Stokes $Q$ and $U$. We find an excess of $B$-mode power over the base lensed-LCDM expectation in the range $30< \ell< 150$, inconsistent with the null hypothesis at a significance of $> 5\sigma$. Through jackknife tests and simulations we show that systematic contamination is much smaller than the observed excess. We also examine a number of available models of polarized dust emission and find that at their default parameter values they predict power $\sim(5-10)\times$ smaller than the observed excess signal. However, these models are not sufficiently constrained to exclude the possibility of dust emission bright enough to explain the entire excess signal. Cross correlating BICEP2 against 100 GHz maps from the BICEP1 experiment, the excess signal is confirmed and its spectral index is found to be consistent with that of the CMB, disfavoring dust at $1.7\sigma$. The observed $B$-mode power spectrum is well fit by a lensed-LCDM + tensor theoretical model with tensor-to-scalar ratio $r=0.20^{+0.07}_{-0.05}$, with $r=0$ disfavored at $7.0\sigma$. Accounting for the contribution of foreground dust will shift this value downward by an amount which will be better constrained with upcoming data sets.

We present a measurement of the angular bispectrum of the millimeter-wave sky in observing bands centered at roughly 95, 150, and 220 GHz, on angular scales of $1^\prime \lesssim \theta \lesssim 10^\prime$ (multipole number $1000 \lesssim l \lesssim 10000$). At these frequencies and angular scales, the main contributions to the bispectrum are expected to be the thermal Sunyaev-Zel'dovich (tSZ) effect and emission from extragalactic sources, predominantly dusty, star-forming galaxies (DSFGs) and active galactic nuclei. We measure the bispectrum in 800 $\mathrm{deg}^2$ of three-band South Pole Telescope data, and we use a multi-frequency fitting procedure to separate the bispectrum of the tSZ effect from the extragalactic source contribution. We simultaneously detect the bispectrum of the tSZ effect at $>$10$\sigma$, the unclustered component of the extragalactic source bispectrum at $>$5$\sigma$ in each frequency band, and the bispectrum due to the clustering of DSFGs---i.e., the clustered cosmic infrared background (CIB) bispectrum---at $>$5$\sigma$. This is the first reported detection of the clustered CIB bispectrum. We use the measured tSZ bispectrum amplitude, compared to model predictions, to constrain the normalization of the matter power spectrum to be $\sigma_8 = 0.787 \pm 0.031$ and to predict the amplitude of the tSZ power spectrum at $l = 3000$. This prediction improves our ability to separate the thermal and kinematic contributions to the total SZ power spectrum. The addition of bispectrum data improves our constraint on the tSZ power spectrum amplitude by a factor of two compared to power spectrum measurements alone and demonstrates a preference for a nonzero kinematic SZ (kSZ) power spectrum, with a derived constraint on the kSZ amplitude at $l=3000$ of A_kSZ $ = 2.9 \pm 1.6 \ \mu$K$^2$, or A_kSZ $ = 2.6 \pm 1.8 \ \mu$K$^2$ if the default A_kSZ > 0 prior is removed.

We explore extensions to the $\Lambda$CDM cosmology using measurements of the cosmic microwave background (CMB) from the recent SPT-SZ survey, along with data from WMAP7 and measurements of $H_0$ and BAO. We check for consistency within $\Lambda$CDM between these datasets, and find some tension. The CMB alone gives weak support to physics beyond $\Lambda$CDM, due to a slight trend relative to $\Lambda$CDM of decreasing power towards smaller angular scales. While it may be due to statistical fluctuation, this trend could also be explained by several extensions. We consider running index (nrun), as well as two extensions that modify the damping tail power (the primordial helium abundance $Y_p$ and the effective number of neutrino species $N_{\rm eff}$) and one that modifies the large-scale power due to the ISW effect (the sum of neutrino masses $\sum m_\nu$). These extensions have similar observational consequences and are partially degenerate when considered simultaneously. Of the 6 one-parameter extensions considered, we find CMB to have the largest preference for nrun with -0.046<nrun<-0.003 at 95% confidence, which strengthens to a 2.7$\sigma$ indication of nrun<0 from CMB+BAO+$H_0$. Detectable non-zero nrun is difficult to explain in the context of single-field, slow-roll inflation models. We find $N_{\rm eff}=3.62\pm0.48$ for the CMB, which tightens to $N_{\rm eff}=3.71\pm0.35$ from CMB+BAO+$H_0$. Larger values of $N_{\rm eff}$ relieve the mild tension between CMB, BAO and $H_0$. When the SZ selected galaxy cluster abundances ($\rm{SPT_{CL}}$) data are also included, we obtain $N_{\rm eff}=3.29\pm0.31$. Allowing for $\sum m_\nu$ gives a 3$\sigma$ detection of $\sum m_\nu$>0 from CMB+BAO+$H_0$+$\rm{SPT_{CL}}$. The median value is $(0.32\pm0.11)$ eV, a factor of six above the lower bound set by neutrino oscillation observations. ... [abridged]

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.

We present a measurement of the cosmic microwave background (CMB) temperature power spectrum using data from the recently completed South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. This measurement is made from observations of 2540 deg$^2$ of sky with arcminute resolution at $150\,$GHz, and improves upon previous measurements using the SPT by tripling the sky area. We report CMB temperature anisotropy power over the multipole range $650<\ell<3000$. We fit the SPT bandpowers, combined with the seven-year Wilkinson Microwave Anisotropy Probe (WMAP7) data, with a six-parameter LCDM cosmological model and find that the two datasets are consistent and well fit by the model. Adding SPT measurements significantly improves LCDM parameter constraints; in particular, the constraint on $\theta_s$ tightens by a factor of 2.7. The impact of gravitational lensing is detected at $8.1\, \sigma$, the most significant detection to date. This sensitivity of the SPT+WMAP7 data to lensing by large-scale structure at low redshifts allows us to constrain the mean curvature of the observable universe with CMB data alone to be $\Omega_k=-0.003^{+0.014}_{-0.018}$. Using the SPT+WMAP7 data, we measure the spectral index of scalar fluctuations to be $n_s=0.9623 \pm 0.0097$ in the LCDM model, a $3.9\,\sigma$ preference for a scale-dependent spectrum with $n_s<1$. The SPT measurement of the CMB damping tail helps break the degeneracy that exists between the tensor-to-scalar ratio $r$ and $n_s$ in large-scale CMB measurements, leading to an upper limit of $r<0.18$ (95%,C.L.) in the LCDM+$r$ model. Adding low-redshift measurements of the Hubble constant ($H_0$) and the baryon acoustic oscillation (BAO) feature to the SPT+WMAP7 data leads to further improvements. The combination of SPT+WMAP7+$H_0$+BAO constrains $n_s=0.9538 \pm 0.0081$ in the LCDM model, a $5.7\,\sigma$ detection of $n_s < 1$, ... [abridged]

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

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

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

We present a point source catalog from 771 square degrees of the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey at 95, 150, and 220 GHz. We detect 1545 sources above 4.5 sigma significance in at least one band. Based on their relative brightness between survey bands, we classify the sources into two populations, one dominated by synchrotron emission from active galactic nuclei, and one dominated by thermal emission from dust-enshrouded star-forming galaxies. We find 1238 synchrotron and 307 dusty sources. We cross-match all sources against external catalogs and find 189 unidentified synchrotron sources and 189 unidentified dusty sources. The dusty sources without counterparts are good candidates for high-redshift, strongly lensed submillimeter galaxies. We derive number counts for each population from 1 Jy down to roughly 9, 5, and 11 mJy at 95, 150, and 220 GHz. We compare these counts with galaxy population models and find that none of the models we consider for either population provide a good fit to the measured counts in all three bands. The disparities imply that these measurements will be an important input to the next generation of millimeter-wave extragalactic source population models.

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

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

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

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

The epoch of reionization is a milestone of cosmological structure formation, marking the birth of the first objects massive enough to yield large numbers of ionizing photons. The mechanism and timescale of reionization remain largely unknown. Measurements of the CMB Doppler effect from ionizing bubbles embedded in large-scale velocity streams (the patchy kinetic Sunyaev-Zel'dovich effect) can constrain the duration of reionization. When combined with large-scale CMB polarization measurements, the evolution of the ionized fraction can be inferred. Using new multi-frequency data from the South Pole Telescope (SPT), we show that the ionized fraction evolved relatively rapidly. For our basic foreground model, we find the kinetic Sunyaev-Zel'dovich power sourced by reionization at l=3000 to be <= 2.1 micro K^2 at 95% CL. Using reionization simulations, we translate this to a limit on the duration of reionization of Delta z <= 4.4 (95% CL). We find that this constraint depends on assumptions about the angular correlation between the thermal Sunyaev-Zel'dovich power and the cosmic infrared background (CIB). Introducing the degree of correlation as a free parameter, we find that the limits on kSZ power weaken to <= 4.9 micro K^2, implying Delta z <= 7.9 (95% CL). We combine the SPT constraint on the duration of reionization with the WMAP7 measurement of the integrated optical depth to probe the cosmic ionization history. We find that reionization ended with 95% CL at z > 7.2 under the assumption of no tSZ-CIB correlation, and z>5.8 when correlations are allowed. Improved constraints from the full SPT data set in conjunction with upcoming Herschel and Planck data should detect extended reionization at >95% CL provided Delta z >= 4. (abbreviated)

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.

Between the BICEP2 and Keck Array experiments, we have deployed over 1500 dual polarized antenna coupled bolometers to map the Cosmic Microwave Background's polarization. We have been able to rapidly deploy these detectors because they are completely planar with an integrated phased-array antenna. Through our experience in these experiments, we have learned of several challenges with this technology- specifically the beam synthesis in the antenna- and in this paper we report on how we have modified our designs to mitigate these challenges. In particular, we discus differential steering errors between the polarization pairs' beam centroids due to microstrip cross talk and gradients of penetration depth in the niobium thin films of our millimeter wave circuits. We also discuss how we have suppressed side lobe response with a Gaussian taper of our antenna illumination pattern. These improvements will be used in Spider, Polar-1, and this season's retrofit of Keck Array.

The Keck Array (SPUD) began observing the cosmic microwave background's polarization in the winter of 2011 at the South Pole. The Keck Array follows the success of the predecessor experiments Bicep and Bicep2, using five on-axis refracting telescopes. These have a combined imaging array of 2500 antenna-coupled TES bolometers read with a SQUID-based time domain multiplexing system. We will discuss the detector noise and the optimization of the readout. The achieved sensitivity of the Keck Array is 11.5 {\mu}K_(CMB)*sqrt{s} in the 2012 configuration.

The Keck Array (SPUD) is a set of microwave polarimeters that observes from the South Pole at degree angular scales in search of a signature of Inflation imprinted as B-mode polarization in the Cosmic Microwave Background (CMB). The first three Keck Array receivers were deployed during the 2010-2011 Austral summer, followed by two new receivers in the 2011-2012 summer season, completing the full five-receiver array. All five receivers are currently observing at 150 GHz. The Keck Array employs the field-proven BICEP/BICEP2 strategy of using small, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. We describe our efforts to characterize the main beam shape and beam shape mismatch between co-located orthogonally-polarized detector pairs, and discuss the implications of measured differential beam parameters on temperature to polarization leakage in CMB analysis.

The BICEP2 and Keck Array experiments are designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 2-4 degrees (l=50-100). This is the region in which the B-mode signal, a signature prediction of cosmic inflation, is expected to peak. BICEP2 was deployed to the South Pole at the end of 2009 and is in the middle of its third year of observing with 500 polarization-sensitive detectors at 150 GHz. The Keck Array was deployed to the South Pole at the end of 2010, initially with three receivers--each similar to BICEP2. An additional two receivers have been added during the 2011-12 summer. We give an overview of the two experiments, report on substantial gains in the sensitivity of the two experiments after post-deployment optimization, and show preliminary maps of CMB polarization from BICEP2.

A technological milestone for experiments employing Transition Edge Sensor (TES) bolometers operating at sub-kelvin temperature is the deployment of detector arrays with 100s--1000s of bolometers. One key technology for such arrays is readout multiplexing: the ability to read out many sensors simultaneously on the same set of wires. This paper describes a frequency-domain multiplexed readout system which has been developed for and deployed on the APEX-SZ and South Pole Telescope millimeter wavelength receivers. In this system, the detector array is divided into modules of seven detectors, and each bolometer within the module is biased with a unique ~MHz sinusoidal carrier such that the individual bolometer signals are well separated in frequency space. The currents from all bolometers in a module are summed together and pre-amplified with Superconducting Quantum Interference Devices (SQUIDs) operating at 4 K. Room-temperature electronics demodulate the carriers to recover the bolometer signals, which are digitized separately and stored to disk. This readout system contributes little noise relative to the detectors themselves, is remarkably insensitive to unwanted microphonic excitations, and provides a technology pathway to multiplexing larger numbers of sensors.

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

We present the first three-frequency South Pole Telescope (SPT) cosmic microwave background (CMB) power spectra. The band powers presented here cover angular scales 2000 < ell < 9400 in frequency bands centered at 95, 150, and 220 GHz. At these frequencies and angular scales, a combination of the primary CMB anisotropy, thermal and kinetic Sunyaev-Zel'dovich (SZ) effects, radio galaxies, and cosmic infrared background (CIB) contributes to the signal. We combine Planck and SPT data at 220 GHz to constrain the amplitude and shape of the CIB power spectrum and find strong evidence for non-linear clustering. We explore the SZ results using a variety of cosmological models for the CMB and CIB anisotropies and find them to be robust with one exception: allowing for spatial correlations between the thermal SZ effect and CIB significantly degrades the SZ constraints. Neglecting this potential correlation, we find the thermal SZ power at 150 GHz and ell = 3000 to be 3.65 +/- 0.69 muK^2, and set an upper limit on the kinetic SZ power to be less than 2.8 muK^2 at 95% confidence. When a correlation between the thermal SZ and CIB is allowed, we constrain a linear combination of thermal and kinetic SZ power: D_{3000}^{tSZ} + 0.5 D_{3000}^{kSZ} = 4.60 +/- 0.63 muK^2, consistent with earlier measurements. We use the measured thermal SZ power and an analytic, thermal SZ model calibrated with simulations to determine sigma8 = 0.807 +/- 0.016. Modeling uncertainties involving the astrophysics of the intracluster medium rather than the statistical uncertainty in the measured band powers are the dominant source of uncertainty on sigma8 . We also place an upper limit on the kinetic SZ power produced by patchy reionization; a companion paper uses these limits to constrain the reionization history of the Universe.

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 South Pole Telescope data from 2008 and 2009 to detect the non-Gaussian signature in the cosmic microwave background (CMB) produced by gravitational lensing and to measure the power spectrum of the projected gravitational potential. We constrain the ratio of the measured amplitude of the lensing signal to that expected in a fiducial LCDM cosmological model to be 0.86 +/- 0.16, with no lensing disfavored at 6.3 sigma. Marginalizing over LCDM cosmological models allowed by the WMAP7 results in a measurement of A_lens=0.90+/-0.19, indicating that the amplitude of matter fluctuations over the redshift range 0.5 <~ z <~ 5 probed by CMB lensing is in good agreement with predictions. We present the results of several consistency checks. These include a clear detection of the lensing signature in CMB maps filtered to have no overlap in Fourier space, as well as a "curl" diagnostic that is consistent with the signal expected for LCDM. We perform a detailed study of bias in the measurement due to noise, foregrounds, and other effects and determine that these contributions are relatively small compared to the statistical uncertainty in the measurement. We combine this lensing measurement with results from WMAP7 to improve constraints on cosmological parameters when compared to those from WMAP7 alone: we find a factor of 3.9 improvement in the measurement of the spatial curvature of the Universe, Omega_k=-0.0014+/-0.0172; a 10% improvement in the amplitude of matter fluctuations within LCDM, sigma_8=0.810+/ 0.026; and a 5% improvement in the dark energy equation of state, w=-1.04+/-0.40. When compared with the measurement of w provided by the combination of WMAP7 and external constraints on the Hubble parameter, the addition of the lensing data improve the measurement of w by 15% to give w=-1.087+/-0.096.

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

The South Pole Telescope (SPT) has nearly completed a 2500-square-degree survey of the southern sky in three frequency bands. Here we present the first public release of SPT maps and associated data products. We present arcminute-resolution maps at 150 GHz and 220 GHz of an approximately 95-square-degree field centered at R.A. 82.7 degrees, decl. -55 degrees. The field was observed to a depth of approximately 17 micro-K arcmin at 150 GHz and 41 micro-K arcmin at 220 GHz during the 2008 austral winter season. Two variations on map filtering and map projection are presented, one tailored for producing catalogs of galaxy clusters detected through their Sunyaev-Zel'dovich effect signature and one tailored for producing catalogs of emissive sources. We describe the data processing pipeline, and we present instrument response functions, filter transfer functions, and map noise properties. All data products described in this paper are available for download at http://pole.uchicago.edu/public/data/maps/ra5h30dec-55 and from the NASA Legacy Archive for Microwave Background Data Analysis server. This is the first step in the eventual release of data from the full 2500-square-degree SPT survey.

We present a measurement of the angular power spectrum of the cosmic microwave background (CMB) using data from the South Pole Telescope (SPT). The data consist of 790 square degrees of sky observed at 150 GHz during 2008 and 2009. Here we present the power spectrum over the multipole range 650 < ell < 3000, where it is dominated by primary CMB anisotropy. We combine this power spectrum with the power spectra from the seven-year Wilkinson Microwave Anisotropy Probe (WMAP) data release to constrain cosmological models. We find that the SPT and WMAP data are consistent with each other and, when combined, are well fit by a spatially flat, LCDM cosmological model. The SPT+WMAP constraint on the spectral index of scalar fluctuations is ns = 0.9663 +/- 0.0112. We detect, at ~5-sigma significance, the effect of gravitational lensing on the CMB power spectrum, and find its amplitude to be consistent with the LCDM cosmological model. We explore a number of extensions beyond the LCDM model. Each extension is tested independently, although there are degeneracies between some of the extension parameters. We constrain the tensor-to-scalar ratio to be r < 0.21 (95% CL) and constrain the running of the scalar spectral index to be dns/dlnk = -0.024 +/- 0.013. We strongly detect the effects of primordial helium and neutrinos on the CMB; a model without helium is rejected at 7.7-sigma, while a model without neutrinos is rejected at 7.5-sigma. The primordial helium abundance is measured to be Yp = 0.296 +/- 0.030, and the effective number of relativistic species is measured to be Neff = 3.85 +/- 0.62. The constraints on these models are strengthened when the CMB data are combined with measurements of the Hubble constant and the baryon acoustic oscillation feature. Notable improvements include ns = 0.9668 +/- 0.0093, r < 0.17 (95% CL), and Neff = 3.86 +/- 0.42. The SPT+WMAP data show...

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.

The Keck Array is a cosmic microwave background (CMB) polarimeter that will begin observing from the South Pole in late 2010. The initial deployment will consist of three telescopes similar to BICEP2 housed in ultra-compact, pulse tube cooled cryostats. Two more receivers will be added the following year. In these proceedings we report on the design and performance of the Keck cryostat. We also report some initial results on the performance of antenna-coupled TES detectors operating in the presence of a pulse tube. We find that the performance of the detectors is not seriously impacted by the replacement of BICEP2's liquid helium cryostat with a pulse tube cooled cryostat.

We report measurements of the cosmic microwave background (CMB) power spectrum from the complete 2008 South Pole Telescope (SPT) data set. We analyze twice as much data as the first SPT power spectrum analysis, using an improved cosmological parameter estimator which fits multi-frequency models to the SPT 150 and $220\,$GHz bandpowers. We find an excellent fit to the measured bandpowers with a model that includes lensed primary CMB anisotropy, secondary thermal (tSZ) and kinetic (kSZ) Sunyaev-Zel'dovich anisotropies, unclustered synchrotron point sources, and clustered dusty point sources. In addition to measuring the power spectrum of dusty galaxies at high signal-to-noise, the data primarily constrain a linear combination of the kSZ and tSZ anisotropy contributions at $150\,$GHz and $\ell=3000$: $D^{tSZ}_{3000} + 0.5\,D^{kSZ}_{3000} = 4.5\pm 1.0 \,\mu{\rm K}^2$. The 95% confidence upper limits on secondary anisotropy power are $D^{tSZ}_{3000} < 5.3\,\mu{\rm K}^2$ and $D^{kSZ}_{3000} < 6.5\,\mu{\rm K}^2$. We also consider the potential correlation of dusty and tSZ sources, and find it incapable of relaxing the tSZ upper limit. These results increase the significance of the lower than expected tSZ amplitude previously determined from SPT power spectrum measurements. We find that models including non-thermal pressure support in groups and clusters predict tSZ power in better agreement with the SPT data. Combining the tSZ power measurement with primary CMB data halves the statistical uncertainty on $\sigma_8$. However, the preferred value of $\sigma_8$ varies significantly between tSZ models. Improved constraints on cosmological parameters from tSZ power spectrum measurements require continued progress in the modeling of the tSZ power.

The South Pole Telescope (SPT) is a 10 m diameter, wide-field, offset Gregorian telescope with a 966-pixel, multi-color, millimeter-wave, bolometer camera. It is located at the Amundsen-Scott South Pole station in Antarctica. The design of the SPT emphasizes careful control of spillover and scattering, to minimize noise and false signals due to ground pickup. The key initial project is a large-area survey at wavelengths of 3, 2 and 1.3 mm, to detect clusters of galaxies via the Sunyaev-Zeldovich effect and to measure the small-scale angular power spectrum of the cosmic microwave background (CMB). The data will be used to characterize the primordial matter power spectrum and to place constraints on the equation of state of dark energy. A second-generation camera will measure the polarization of the CMB, potentially leading to constraints on the neutrino mass and the energy scale of inflation.

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

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

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

BICEP2/Keck and SPIDER are cosmic microwave background (CMB) polarimeters targeting the B-mode polarization induced by primordial gravitational waves from inflation. They will be using planar arrays of polarization sensitive antenna-coupled TES bolometers, operating at frequencies between 90 GHz and 220 GHz. At 150 GHz each array consists of 64 polarimeters and four of these arrays are assembled together to make a focal plane, for a total of 256 dual-polarization elements (512 TES sensors). The detector arrays are integrated with a time-domain SQUID multiplexer developed at NIST and read out using the multi-channels electronics (MCE) developed at the University of British Columbia. Following our progress in improving detector parameters uniformity across the arrays and fabrication yield, our main effort has focused on improving detector arrays optical and noise performances, in order to produce science grade focal planes achieving target sensitivities. We report on changes in detector design implemented to optimize such performances and following focal plane arrays characterization. BICEP2 has deployed a first 150 GHz science grade focal plane to the South Pole in December 2009.

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

We use data from the first 100 square-degree field observed by the South Pole Telescope (SPT) in 2008 to measure the angular power spectrum of temperature anisotropies contributed by the background of dusty star-forming galaxies (DSFGs) at millimeter wavelengths. From the auto and cross-correlation of 150 and 220 GHz SPT maps, we significantly detect both Poisson distributed and, for the first time at millimeter wavelengths, clustered components of power from a background of DSFGs. The spectral indices between 150 and 220 GHz of the Poisson and clustered components are found to be 3.86 +- 0.23 and 3.8 +- 1.3 respectively, implying a steep scaling of the dust emissivity index beta ~ 2. The Poisson and clustered power detected in SPT, BLAST (at 600, 860, and 1200 GHz), and Spitzer (1900 GHz) data can be understood in the context of a simple model in which all galaxies have the same graybody spectrum with dust emissivity index of beta = 2 and dust temperature T_d = 34 K. In this model, half of the 150 GHz background light comes from redshifts greater than 3.2. We also use the SPT data to place an upper limit on the amplitude of the kinetic Sunyaev-Zel'dovich power spectrum at l = 3000 of 13 uK^2 at 95% confidence.

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

We present Sunyaev-Zel'dovich measurements of 15 massive X-ray selected galaxy clusters obtained with the South Pole Telescope. The Sunyaev-Zel'dovich (SZ) cluster signals are measured at 150 GHz, and concurrent 220 GHz data are used to reduce astrophysical contamination. Radial profiles are computed using a technique that takes into account the effects of the beams and filtering. In several clusters, significant SZ decrements are detected out to a substantial fraction of the virial radius. The profiles are fit to the beta model and to a generalized NFW pressure profile, and are scaled and stacked to probe their average behavior. We find model parameters that are consistent with previous studies: beta=0.86 and r_core/r_500 = 0.20 for the beta model, and (alpha, beta, gamma, c_500)=(1.0,5.5,0.5,1.0) for the generalized NFW model. Both models fit the SPT data comparably well, and both are consistent with the average SZ profile out to the virial radius. The integrated Compton-y parameter Y_SZ is computed for each cluster using both model-dependent and model-independent techniques, and the results are compared to X-ray estimates of cluster parameters. We find that Y_SZ scales with Y_X and gas mass with low scatter. Since these observables have been found to scale with total mass, our results point to a tight mass-observable relation for the SPT cluster survey.

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

We report cosmic microwave background (CMB) power spectrum measurements from the first 100 sq. deg. field observed by the South Pole Telescope (SPT) at 150 and 220 GHz. On angular scales where the primary CMB anisotropy is dominant, ell ~< 3000, the SPT power spectrum is consistent with the standard LambdaCDM cosmology. On smaller scales, we see strong evidence for a point source contribution, consistent with a population of dusty, star-forming galaxies. After we mask bright point sources, anisotropy power on angular scales of 3000 < ell < 9500 is detected with a signal-to-noise > 50 at both frequencies. We combine the 150 and 220 GHz data to remove the majority of the point source power, and use the point source subtracted spectrum to detect Sunyaev-Zel'dovich (SZ) power at 2.6 sigma. At ell=3000, the SZ power in the subtracted bandpowers is 4.2 +/- 1.5 uK^2, which is significantly lower than the power predicted by a fiducial model using WMAP5 cosmological parameters. This discrepancy may suggest that contemporary galaxy cluster models overestimate the thermal pressure of intracluster gas. Alternatively, this result can be interpreted as evidence for lower values of sigma8. When combined with an estimate of the kinetic SZ contribution, the measured SZ amplitude shifts sigma8 from the primary CMB anisotropy derived constraint of 0.794 +/- 0.028 down to 0.773 +/- 0.025. The uncertainty in the constraint on sigma8 from this analysis is dominated by uncertainties in the theoretical modeling required to predict the amplitude of the SZ power spectrum for a given set of cosmological parameters.

The South Pole Telescope (SPT) is conducting a Sunyaev-Zel'dovich (SZ) effect survey over large areas of the southern sky, searching for massive galaxy clusters to high redshift. In this preliminary study, we focus on a 40 square-degree area targeted by the Blanco Cosmology Survey (BCS), which is centered roughly at right ascension 5h30m, declination -53 degrees. Over two seasons of observations, this entire region has been mapped by the SPT at 95 GHz, 150 GHz, and 225 GHz. We report the four most significant SPT detections of SZ clusters in this field, three of which were previously unknown and, therefore, represent the first galaxy clusters discovered with an SZ survey. The SZ clusters are detected as decrements with greater than 5-sigma significance in the high-sensitivity 150 GHz SPT map. The SZ spectrum of these sources is confirmed by detections of decrements at the corresponding locations in the 95 GHz SPT map and non-detections at those locations in the 225 GHz SPT map. Multiband optical images from the BCS survey demonstrate significant concentrations of similarly colored galaxies at the positions of the SZ detections. Photometric redshift estimates from the BCS data indicate that two of the clusters lie at moderate redshift (z ~ 0.4) and two at high redshift (z >~ 0.8). One of the SZ detections was previously identified as a galaxy cluster using X-ray data from the ROSAT All-Sky Survey (RASS). Potential RASS counterparts (not previously identified as clusters) are also found for two of the new discoveries. These first four galaxy clusters are the most significant SZ detections from a subset of the ongoing SPT survey. As such, they serve as a demonstration that SZ surveys, and the SPT in particular, can be an effective means for finding galaxy clusters.

In this paper, we present results from the complete set of cosmic microwave background (CMB) radiation temperature anisotropy observations made with the Arcminute Cosmology Bolometer Array Receiver (ACBAR) operating at 150 GHz. We include new data from the final 2005 observing season, expanding the number of detector-hours by 210% and the sky coverage by 490% over that used for the previous ACBAR release. As a result, the band-power uncertainties have been reduced by more than a factor of two on angular scales encompassing the third to fifth acoustic peaks as well as the damping tail of the CMB power spectrum. The calibration uncertainty has been reduced from 6% to 2.1% in temperature through a direct comparison of the CMB anisotropy measured by ACBAR with that of the dipole-calibrated WMAP5 experiment. The measured power spectrum is consistent with a spatially flat, LambdaCDM cosmological model. We include the effects of weak lensing in the power spectrum model computations and find that this significantly improves the fits of the models to the combined ACBAR+WMAP5 power spectrum. The preferred strength of the lensing is consistent with theoretical expectations. On fine angular scales, there is weak evidence (1.1 sigma) for excess power above the level expected from primary anisotropies. We expect any excess power to be dominated by the combination of emission from dusty protogalaxies and the Sunyaev-Zel'dovich effect (SZE). However, the excess observed by ACBAR is significantly smaller than the excess power at ell > 2000 reported by the CBI experiment operating at 30 GHz. Therefore, while it is unlikely that the CBI excess has a primordial origin; the combined ACBAR and CBI results are consistent with the source of the CBI excess being either the SZE or radio source contamination.

A new 10 meter diameter telescope is being constructed for deployment at the NSF South Pole research station. The telescope is designed for conducting large-area millimeter and sub-millimeter wave surveys of faint, low contrast emission, as required to map primary and secondary anisotropies in the cosmic microwave background. To achieve the required sensitivity and resolution, the telescope design employs an off-axis primary with a 10m diameter clear aperture. The full aperture and the associated optics will have a combined surface accuracy of better than 20 microns rms to allow precision operation in the submillimeter atmospheric windows. The telescope will be surrounded with a large reflecting ground screen to reduce sensitivity to thermal emission from the ground and local interference. The optics of the telescope will support a square degree field of view at 2mm wavelength and will feed a new 1000-element micro-lithographed planar bolometric array with superconducting transition-edge sensors and frequency-multiplexed readouts. The first key project will be to conduct a survey over approximately 4000 degrees for galaxy clusters using the Sunyaev-Zel'dovich Effect. This survey should find many thousands of clusters with a mass selection criteria that is remarkably uniform with redshift. Armed with redshifts obtained from optical and infrared follow-up observations, it is expected that the survey will enable significant constraints to be placed on the equation of state of the dark energy.