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

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

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

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

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

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

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

The cosmic microwave background (CMB) encodes information about the content and evolution of the universe. The presence of light, weakly interacting particles impacts the expansion history of the early universe, which alters the temperature and polarization anisotropies of the CMB. In this way, current measurements of the CMB place interesting constraints on the neutrino energy density and mass, as well as on the abundance of other possible light relativistic particle species. We present the status of an on-going 1500 sq. deg. survey with the SPT-3G receiver, a new mm-wavelength camera on the 10-m diameter South Pole Telescope (SPT). The SPT-3G camera consists of 16,000 superconducting transition edge sensors, a 10x increase over the previous generation camera, which allows it to map the CMB with an unprecedented combination of sensitivity and angular resolution. We highlight projected constraints on the abundance of sterile neutrinos and the sum of the neutrino masses for the SPT-3G survey, which could help determine the neutrino mass hierarchy.

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

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

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

Measurements of the Hubble constant, and more generally measurements of the expansion rate and distances over the interval $0 < z < 1$, appear to be inconsistent with the predictions of the standard cosmological model ($\Lambda$CDM) given observations of cosmic microwave background temperature and polarization anisotropies. Here we consider a variety of types of departures from $\Lambda$CDM that could, in principle, restore concordance among these datasets, and we explain why we find almost all of them unlikely to be successful. We single out the set of solutions that increase the expansion rate in the decade of scale factor expansion just prior to recombination as the least unlikely. These solutions are themselves tightly constrained by their impact on photon diffusion and on the gravitational driving of acoustic oscillations of the modes that begin oscillating during this epoch -- modes that project on to angular scales that are very well measured. We point out that a general feature of such solutions is a residual to fits to $\Lambda$CDM, like the one observed in Planck power spectra. This residual drives the modestly significant inferences of angular-scale dependence to the matter density and anomalously high lensing power, puzzling aspects of a data set that is otherwise extremely well fit by $\Lambda$CDM.

Separating galactic foreground emission from maps of the cosmic microwave background (CMB), and quantifying the uncertainty in the CMB maps due to errors in foreground separation are important for avoiding biases in scientific conclusions. Our ability to quantify such uncertainty is limited by our lack of a model for the statistical distribution of the foreground emission. Here we use a Deep Convolutional Generative Adversarial Network (DCGAN) to create an effective non-Gaussian statistical model for intensity of emission by interstellar dust. For training data we use a set of dust maps inferred from observations by the Planck satellite. A DCGAN is uniquely suited for such unsupervised learning tasks as it can learn to model a complex non-Gaussian distribution directly from examples. We then use these simulations to train a second neural network to estimate the underlying CMB signal from dust-contaminated maps. We discuss other potential uses for the trained DCGAN, and the generalization to polarized emission from both dust and synchrotron.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

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

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 the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

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

Type Ia Supernovae, calibrated by classical distance ladder methods, can be used, in conjunction with galaxy survey two-point correlation functions, to empirically determine the size of the sound horizon $r_{\rm s}$. Assumption of the $\Lambda$CDM model, together with data to constrain its parameters, can also be used to determine the size of the sound horizon. Using a variety of cosmic microwave background (CMB) datasets to constrain $\Lambda$CDM parameters, we find the model-based sound horizon to be larger than the empirically-determined one with a statistical significance of between 2 and 3$\sigma$, depending on the dataset. If reconciliation requires a change to the cosmological model, we argue that change is likely to be important in the two decades of scale factor evolution prior to recombination. Future CMB observations will therefore likely be able to test any such adjustments; e.g., a third generation CMB survey like SPT-3G can achieve a three-fold improvement in the constraints on $r_{\rm s}$ in the $\Lambda$CDM model extended to allow additional light degrees of freedom.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

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

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

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.

With low redshift probes reaching unprecedented precision, uncertainty of the CMB optical depth is expected to be the limiting factor for future cosmological neutrino mass constraints. In this paper, we discuss to what extent combinations of CMB lensing and galaxy surveys measurements at low redshifts $z\sim 0.5-5$ will be able to make competitive neutrino mass measurements without relying on any optical depth constraints. We find that the combination of LSST galaxies and CMB-S4 lensing should be able to achieve constraints on the neutrino mass sum of 25meV without optical depth information, an independent measurement that is competitive with or slightly better than the constraint of 30meV possible with CMB-S4 and present-day optical depth measurements. These constraints originate both in structure growth probed by cross-correlation tomography over a wide redshift range as well as, most importantly, the shape of the galaxy power spectrum measured over a large volume. We caution that possible complications such as higher-order biasing and systematic errors in the analysis of high redshift galaxy clustering are only briefly discussed and may be non-negligible. Nevertheless, our results show that new kinds of high-precision neutrino mass measurements at and beyond the present-day optical depth limit may be possible.

The South Pole Telescope (SPT) is a millimeter-wavelength telescope designed for high-precision measurements of the cosmic microwave background (CMB). The SPT measures both the temperature and polarization of the CMB with a large aperture, resulting in high resolution maps sensitive to signals across a wide range of angular scales on the sky. With these data, the SPT has the potential to make a broad range of cosmological measurements. These include constraining the effect of massive neutrinos on large-scale structure formation as well as cleaning galactic and cosmological foregrounds from CMB polarization data in future searches for inflationary gravitational waves. The SPT began observing in January 2017 with a new receiver (SPT-3G) containing $\sim$16,000 polarization-sensitive transition-edge sensor bolometers. Several key technology developments have enabled this large-format focal plane, including advances in detectors, readout electronics, and large millimeter-wavelength optics. We discuss the implementation of these technologies in the SPT-3G receiver as well as the challenges they presented. In late 2017 the implementations of all three of these technologies were modified to optimize total performance. Here, we present the current instrument status of the SPT-3G receiver.

The Probe of Inflation and Cosmic Origins (PICO) is a probe-class mission concept currently under study by NASA. PICO will probe the physics of the Big Bang and the energy scale of inflation, constrain the sum of neutrino masses, measure the growth of structures in the universe, and constrain its reionization history by making full sky maps of the cosmic microwave background with sensitivity 80 times higher than the Planck space mission. With bands at 21-799 GHz and arcmin resolution at the highest frequencies, PICO will make polarization maps of Galactic synchrotron and dust emission to observe the role of magnetic fields in Milky Way's evolution and star formation. We discuss PICO's optical system, focal plane, and give current best case noise estimates. The optical design is a two-reflector optimized open-Dragone design with a cold aperture stop. It gives a diffraction limited field of view (DLFOV) with throughput of 910 square cm sr at 21 GHz. The large 82 square degree DLFOV hosts 12,996 transition edge sensor bolometers distributed in 21 frequency bands and maintained at 0.1 K. We use focal plane technologies that are currently implemented on operating CMB instruments including three-color multi-chroic pixels and multiplexed readouts. To our knowledge, this is the first use of an open-Dragone design for mm-wave astrophysical observations, and the only monolithic CMB instrument to have such a broad frequency coverage. With current best case estimate polarization depth of 0.65 microK(CMB}-arcmin over the entire sky, PICO is the most sensitive CMB instrument designed to date.

The Probe of Inflation and Cosmic Origins (PICO) is a NASA-funded study of a Probe-class mission concept. The top-level science objectives are to probe the physics of the Big Bang by measuring or constraining the energy scale of inflation, probe fundamental physics by measuring the number of light particles in the Universe and the sum of neutrino masses, to measure the reionization history of the Universe, and to understand the mechanisms driving the cosmic star formation history, and the physics of the galactic magnetic field. PICO would have multiple frequency bands between 21 and 799 GHz, and would survey the entire sky, producing maps of the polarization of the cosmic microwave background radiation, of galactic dust, of synchrotron radiation, and of various populations of point sources. Several instrument configurations, optical systems, cooling architectures, and detector and readout technologies have been and continue to be considered in the development of the mission concept. We will present a snapshot of the baseline mission concept currently under development.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

Temperature and polarization variations across the microwave sky include the fingerprints of quantum fluctuations in the early universe. They may soon reveal physics at unprecedented energy scales.

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

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

Clusters of galaxies gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint in the CMB on arcminute scales. Measurement of this effect offers a promising way to constrain the masses of galaxy clusters, particularly those at high redshift. We use CMB maps from the South Pole Telescope Sunyaev-Zel'dovich (SZ) survey to measure the CMB lensing signal around galaxy clusters identified in optical imaging from first year observations of the Dark Energy Survey. The cluster catalog used in this analysis contains 3697 members with mean redshift of $\bar{z} = 0.45$. We detect lensing of the CMB by the galaxy clusters at $8.1\sigma$ significance. Using the measured lensing signal, we constrain the amplitude of the relation between cluster mass and optical richness to roughly $17\%$ precision, finding good agreement with recent constraints obtained with galaxy lensing. The error budget is dominated by statistical noise but includes significant contributions from systematic biases due to the thermal SZ effect and cluster miscentering.

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.

In this paper, we conduct a search in the latest large-scale structure measurements for signatures of the dark matter-dark radiation interaction proposed by Buen-Abad et al. (2015). We show that prior claims of a detection of this interaction rely on a use of the SZ cluster mass function that ignores uncertainty in the mass-observable relationship. Including this uncertainty we find that the inferred level of interaction remains consistent with the data, but so does zero interaction; i.e., there is no longer a detection. We also point out that inference of the shape and amplitude of the matter power spectrum from Ly$\alpha$ forest measurements is highly inconsistent with the predictions of the $\Lambda$CDM model conditioned on Planck CMB temperature, polarization, and lensing power spectra, and that the dark matter-dark radiation model can restore that consistency. We also phenomenologically generalize the model of Buen-Abad et al. (2015) to allow for interaction rates with different scalings with temperature, and find that the original scaling is preferred by the data.

We present updated constraints on the free-streaming nature of cosmological neutrinos from cosmic microwave background (CMB) power spectra, baryonic acoustic oscillation data, and local measurements of the Hubble constant. Specifically, we consider a Fermi-like four-fermion interaction between massless neutrinos, characterized by an effective coupling constant $ G_{\rm eff}$, and resulting in a neutrino opacity $\dot{\tau}_\nu\propto G_{\rm eff}^2 T_\nu^5$. Using a conservative prior on the parameter $\log_{10}\left(G_{\rm eff} {\rm MeV}^2\right)$, we find a bimodal posterior distribution. The first of these modes is consistent with the standard $\Lambda$CDM cosmology and corresponds to neutrinos decoupling at redshift $z_{\nu,{\rm dec}} > 1.3\times10^5$. The other mode of the posterior, dubbed the "interacting neutrino mode", corresponds to neutrino decoupling occurring within a narrow redshift window centered around $z_{\nu,{\rm dec}}\sim8300$. This mode is characterized by a high value of the effective neutrino coupling constant, together with a lower value of the scalar spectral index and amplitude of fluctuations, and a higher value of the Hubble parameter. Using both a maximum likelihood analysis and the ratio of the two mode's Bayesian evidence, we find the interacting neutrino mode to be statistically disfavored compared to the standard $\Lambda$CDM cosmology. Interestingly, the addition of CMB polarization and direct Hubble constant measurements significantly raises the statistical significance of this secondary mode, indicating that new physics in the neutrino sector could help explain the difference between local measurements of $H_0$, and those inferred from CMB data. A robust consequence of our results is that neutrinos must be free streaming long before the epoch of matter-radiation equality.

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

One of the major targets for next-generation cosmic microwave background (CMB) experiments is the detection of the primordial B-mode signal. Planning is under way for Stage-IV experiments that are projected to have instrumental noise small enough to make lensing and foregrounds the dominant source of uncertainty for estimating the tensor-to-scalar ratio $r$ from polarization maps. This makes delensing a crucial part of future CMB polarization science. In this paper we present a likelihood method for estimating the tensor-to-scalar ratio $r$ from CMB polarization observations, which combines the benefits of a full-scale likelihood approach with the tractability of the quadratic delensing technique. This method is a pixel space, all order likelihood analysis of the quadratic delensed B modes, and it essentially builds upon the quadratic delenser by taking into account all order lensing and pixel space anomalies. Its tractability relies on a crucial factorization of the pixel space covariance matrix of the polarization observations which allows one to compute the full Gaussian approximate likelihood profile, as a function of $r$, at the same computational cost of a single likelihood evaluation.

Recent determination of the Hubble constant via Cepheid-calibrated supernovae by \citet{riess_2.4_2016} (R16) find $\sim 3\sigma$ tension with inferences based on cosmic microwave background temperature and polarization measurements from $Planck$. This tension could be an indication of inadequacies in the concordance $\Lambda$CDM model. Here we investigate the possibility that the discrepancy could instead be due to systematic bias or uncertainty in the Cepheid calibration step of the distance ladder measurement by R16. We consider variations in total-to-selective extinction of Cepheid flux as a function of line-of-sight, hidden structure in the period-luminosity relationship, and potentially different intrinsic color distributions of Cepheids as a function of host galaxy. Considering all potential sources of error, our final determination of $H_0 = 73.3 \pm 1.7~{\rm km/s/Mpc}$ (not including systematic errors from the treatment of geometric distances or Type Ia Supernovae) shows remarkable robustness and agreement with R16. We conclude systematics from the modeling of Cepheid photometry, including Cepheid selection criteria, cannot explain the observed tension between Cepheid-variable and CMB-based inferences of the Hubble constant. Considering a `model-independent' approach to relating Cepheids in galaxies with known distances to Cepheids in galaxies hosting a Type Ia supernova and finding agreement with the R16 result, we conclude no generalization of the model relating anchor and host Cepheid magnitude measurements can introduce significant bias in the $H_0$ inference.

As modern precision cosmological measurements continue to show agreement with the broad features of the standard $\Lambda$-Cold Dark Matter ($\Lambda$CDM) cosmological model, we are increasingly motivated to look for small departures from the standard model's predictions which might not be detected with standard approaches. While searches for extensions and modifications of $\Lambda$CDM have to date turned up no convincing evidence of beyond-the-standard-model cosmology, the list of models compared against $\Lambda$CDM is by no means complete and is often governed by readily-coded modifications to standard Boltzmann codes. Also, standard goodness-of-fit methods such as a naive $\chi^2$ test fail to put strong pressure on the null $\Lambda$CDM hypothesis, since modern datasets have orders of magnitudes more degrees of freedom than $\Lambda$CDM. Here we present a method of tuning goodness-of-fit tests to detect potential sub-dominant extra-$\Lambda$CDM signals present in the data through compressing observations in a way that maximizes extra-$\Lambda$CDM signal variation over noise and $\Lambda$CDM variation. This method, based on a Karhunen-Lo\`{e}ve transformation of the data, is tuned to be maximally sensitive to particular types of variations characteristic of the tuning model; but, unlike direct model comparison, the test is also sensitive to features that only partially mimic the tuning model. As an example of its use, we apply this method in the context of a nonstandard primordial power spectrum compared against the $2015$ $Planck$ CMB temperature and polarization power spectrum. We find weak evidence of extra-$\Lambda$CDM physics, conceivably due to known systematics in the 2015 Planck polarization release.

The angular two-point correlation function of the temperature of the cosmic microwave background (CMB), as inferred from nearly all-sky maps, is very close to zero on large angular scales. A statistic invented to quantify this feature, $S_{1/2}$, has a value sufficiently low that only about 7 in 1000 simulations generated assuming the standard cosmological model have lower values; i.e., it has a $p$-value of 0.007. As such, it is one of several unusual features of the CMB sky on large scales, including the low value of the observed CMB quadrupole, whose importance is unclear: are they multiple and independent clues about physics beyond the cosmological standard model, or an expected consequence of our ability to find signals in Gaussian noise? We find they are not independent: using only simulations with quadrupole values near the observed one, the $S_{1/2}$ $p$-value increases from 0.007 to 0.08. We also find strong evidence that corrections for a "look-elsewhere effect" are large. To do so, we use a one-dimensional generalization of the $S_{1/2}$ statistic, and select along the one dimension for the statistic that is most extreme. Subjecting our simulations to this process increases the $p$-value from 0.007 to 0.03; a result similar to that found in Planck XVI (2016). We argue that this optimization process along the one dimension provides an $underestimate$ of the look-elsewhere effect correction for the historical human process of selecting the $S_{1/2}$ statistic from a very high-dimensional space of alternative statistics $after$ having examined the data.

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.

The six parameters of the standard $\Lambda$CDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We investigate these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium $\tau$, the baryon density $\omega_{\rm b}$, the matter density $\omega_{\rm m}$, the angular size of the sound horizon $\theta_*$, the spectral index of the primordial power spectrum, $n_{\rm s}$, and $A_{\rm s}e^{-2\tau}$ (where $A_{\rm s}$ is the amplitude of the primordial power spectrum), we examine the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment $\ell<800$ in the Planck temperature power spectrum) and an all angular-scale data set ($\ell<2500$ Planck temperature power spectrum), each with a prior on $\tau$ of $0.07\pm0.02$. We find that the shifts, in units of the 1$\sigma$ expected dispersion for each parameter, are $\{\Delta \tau, \Delta A_{\rm s} e^{-2\tau}, \Delta n_{\rm s}, \Delta \omega_{\rm m}, \Delta \omega_{\rm b}, \Delta \theta_*\} = \{-1.7, -2.2, 1.2, -2.0, 1.1, 0.9\}$, with a $\chi^2$ value of 8.0. We find that this $\chi^2$ value is exceeded in 15% of our simulated data sets, and that a parameter deviates by more than 2.2$\sigma$ in 9% of simulated data sets, meaning that the shifts are not unusually large. Comparing $\ell<800$ instead to $\ell>800$, or splitting at a different multipole, yields similar results. We examine the $\ell<800$ model residuals in the $\ell>800$ power spectrum data and find that the features there... [abridged]

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.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.

We present the results of SPT-GMOS, a spectroscopic survey with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South. The targets of SPT-GMOS are galaxy clusters identified in the SPT-SZ survey, a millimeter-wave survey of 2500 squ. deg. of the southern sky using the South Pole Telescope (SPT). Multi-object spectroscopic observations of 62 SPT-selected galaxy clusters were performed between January 2011 and December 2015, yielding spectra with radial velocity measurements for 2595 sources. We identify 2243 of these sources as galaxies, and 352 as stars. Of the galaxies, we identify 1579 as members of SPT-SZ galaxy clusters. The primary goal of these observations was to obtain spectra of cluster member galaxies to estimate cluster redshifts and velocity dispersions. We describe the full spectroscopic dataset and resulting data products, including galaxy redshifts, cluster redshifts and velocity dispersions, and measurements of several well-known spectral indices for each galaxy: the equivalent width, W, of [O II] 3727,3729 and H-delta, and the 4000A break strength, D4000. We use the spectral indices to classify galaxies by spectral type (i.e., passive, post-starburst, star-forming), and we match the spectra against photometric catalogs to characterize spectroscopically-observed cluster members as a function of brightness (relative to m*). Finally, we report several new measurements of redshifts for ten bright, strongly-lensed background galaxies in the cores of eight galaxy clusters. Combining the SPT-GMOS dataset with previous spectroscopic follow-up of SPT-SZ galaxy clusters results in spectroscopic measurements for >100 clusters, or ~20% of the full SPT-SZ sample.

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

Cosmic Microwave Background (CMB) full-sky temperature data show a hemispherical asymmetry in power nearly aligned with the Ecliptic. In real space, this anomaly can be quantified by the temperature variance in the northern and southern Ecliptic hemispheres. In this context, the northern hemisphere displays an anomalously low variance while the southern hemisphere appears unremarkable (consistent with expectations from the best-fitting theory, $\Lambda$CDM). While this is a well established result in temperature, the low signal-to-noise ratio in current polarization data prevents a similar comparison. This will change with a proposed ground-based CMB experiment, CMB-S4. With that in mind, we generate realizations of polarization maps constrained by the temperature data and predict the distribution of the hemispherical variance in polarization considering two different sky coverage scenarios possible in CMB-S4: full Ecliptic north coverage and just the portion of the North that can be observed from a ground based telescope at the high Chilean Atacama plateau. We find that even in the set of realizations constrained by the temperature data, the low northern hemisphere variance observed in temperature is not expected in polarization. Therefore, an anomalously low variance detection in polarization would provide strong evidence against the hypothesis that the temperature anomaly is simply a statistical fluke. We show, within $\Lambda$CDM, how variance measurements in both sky coverage scenarios are related. We find that the variance makes for a good statistic in cases where the sky coverage is limited, however a full northern coverage is still preferable.

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

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $\Lambda$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5{\sigma} from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($\mu$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $\Lambda$CDM parameters consistent with those from temperature alone.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $\tau$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $\tau$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $\tau$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $\Lambda$CDM cosmological model and recent models of reionization.

The Planck collaboration has measured the temperature and polarization of the cosmic microwave background well enough to determine the locations of eight peaks in the temperature (TT) power spectrum, five peaks in the polarization (EE) power spectrum and twelve extrema in the cross (TE) power spectrum. The relative locations of these extrema give a striking, and beautiful, demonstration of what we expect from acoustic oscillations in the plasma; e.g., that EE peaks fall half way between TT peaks. We expect this because the temperature map is predominantly sourced by temperature variations in the last scattering surface, while the polarization map is predominantly sourced by gradients in the velocity field, and the harmonic oscillations have temperature and velocity 90 degrees out of phase. However, there are large differences in expectations for extrema locations from simple analytic models vs. numerical calculations. Here we quantitatively explore the origin of these differences in gravitational potential transients, neutrino free-streaming, the breakdown of tight coupling, the shape of the primordial power spectrum, details of the geometric projection from three to two dimensions, and the thickness of the last scattering surface. We also compare the peak locations determined from Planck measurements to expectations under the $\Lambda$CDM model. Taking into account how the peak locations were determined, we find them to be in agreement.

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

Our tightest upper limit on the sum of neutrino mass eigenvalues $M_\nu$ comes from cosmological observations that will improve substantially in the near future, enabling a detection. The combination of the Baryon Acoustic Oscillation feature measured from the Dark Energy Spectroscopic Instrument and a Stage-IV Cosmic Microwave Background experiment has been forecasted to achieve $\sigma(M_\nu) < 1/3$ of the lower limit on $M_\nu$ from atmospheric and solar neutrino oscillations \citep{2013arXiv1309.5383A,2012PhRvD..86a3012F}. Here we examine in detail the physical effects of neutrino mass on cosmological observables that make these constraints possible. We also consider how these constraints would be improved to ensure at least a $5\sigma$ detection.

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

The unimpeded relativistic propagation of cosmological neutrinos prior to recombination of the baryon-photon plasma alters gravitational potentials and therefore the details of the time-dependent gravitational driving of acoustic oscillations. We report here a first detection of the resulting shifts in the temporal phase of the oscillations, which we infer from their signature in the Cosmic Microwave Background (CMB) temperature power spectrum.

Our velocity relative to the rest frame of the cosmic microwave background (CMB) generates a dipole temperature anisotropy on the sky which has been well measured for more than 30 years, and has an accepted amplitude of v/c = 0.00123, or v = 369km/s. In addition to this signal generated by Doppler boosting of the CMB monopole, our motion also modulates and aberrates the CMB temperature fluctuations (as well as every other source of radiation at cosmological distances). This is an order 0.1% effect applied to fluctuations which are already one part in roughly one hundred thousand, so it is quite small. Nevertheless, it becomes detectable with the all-sky coverage, high angular resolution, and low noise levels of the Planck satellite. Here we report a first measurement of this velocity signature using the aberration and modulation effects on the CMB temperature anisotropies, finding a component in the known dipole direction, (l,b)=(264, 48) [deg], of 384km/s +- 78km/s (stat.) +- 115km/s (syst.). This is a significant confirmation of the expected velocity.

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

The anisotropies in the cosmic microwave background (CMB) provide our best laboratory for testing models of the formation and evolution of large-scale structure. The rich features in the cosmic microwave background anisotropy spectrum, in combination with highly precise observations and theoretical predictions, also allow us to simultaneously constrain a number of cosmological parameters. As observations have progressed, measurements at smaller angular scales have provided increasing leverage. These smaller angular scales provide sensitive measures of the matter density through the effect of gravitational lensing. In this work we provide an analytic explanation of the manner in which the lensing of CMB anisotropies depends on the matter density, finding that the dominant effect comes from the shape of the matter power spectrum set by the decay of small-scale potentials between horizon crossing and matter-radiation equality.

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

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

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $\sigma$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $\chi^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $\chi^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

We give updated constraints on hypothetical light bosons with a two-photon coupling such as axions or axion-like particles (ALPs). We focus on masses and lifetimes where decays happen near big bang nucleosynthesis (BBN), thus altering the baryon-to-photon ratio and number of relativistic degrees of freedom between the BBN epoch and the cosmic microwave background (CMB) last scattering epoch, in particular such that $N_{\rm eff}^{\rm CMB} < N_{\rm eff}^{\rm BBN}$ and $\eta^{\rm CMB} < \eta^{\rm BBN}$. New constraints presented here come from Planck measurements of the CMB power spectrum combined with the latest inferences of primordial $^4$He and D/H abundances. We find that a previously allowed region in parameter space near $m=1\,\rm MeV$ and $\tau=100\,\rm ms$, consistent with a QCD axion arising from a symmetry breaking near the electroweak scale, is now ruled out at $>3\sigma$ by the combination of CMB+D/H measurements if only ALPs and three thermalized neutrino species contribute to $N_{\rm eff}$. The bound relaxes if there are additional light degrees of freedom present which, in this scenario, have their contribution limited to $\Delta N_{\rm eff}=1.1\pm0.3$. We give forecasts showing that a number of experiments are expected to reach the sensitivity needed to further test this region, such as Stage-IV CMB and SUPER-KEKB, the latter a direct test insensitive to any extra degrees of freedom.