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

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

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

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

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

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

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

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

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

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

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

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

We perform a detailed study of the bispectrum of the Sunyaev-Zel'dovich effect. Using an analytical model for the pressure profiles of the intracluster medium, we demonstrate the SZ bispectrum to be a sensitive probe of the amplitude of the matter power spectrum parameter sigma_8. We find that the bispectrum amplitude scales as B_SZ ~ sigma_8^{11-12}, compared to that of the power spectrum, which scales as A_tSZ ~ sigma_8^{7-9}. We show that the SZ bispectrum is principally sourced by massive clusters at redshifts around z~0.4, which have been well-studied observationally. This is in contrast to the SZ power spectrum, which receives a significant contribution from less-well understood low-mass and high-redshift groups and clusters. Therefore, the amplitude of the bispectrum at l~3000 is less sensitive to astrophysical uncertainties than the SZ power spectrum. We show that current high resolution CMB experiments should be able to detect the SZ bispectrum amplitude with high significance, in part due to the low contamination from extra-galactic foregrounds. A combination of the SZ bispectrum and the power spectrum can sharpen the measurements of thermal and kinetic SZ components and help distinguish cosmological and astrophysical information from high-resolution CMB maps.

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

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

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

We present a detailed investigation of the impact of astrophysical processes on the shape and amplitude of the kinetic Sunyaev-Zel'dovich (kSZ) power spectrum from the post-reionization epoch. This is achieved by constructing a new model of the kSZ power spectrum which we calibrate to the results of hydrodynamic simulations. By construction, our method accounts for all relevant density and velocity modes and so is unaffected by the limited box size of our simulations. We find that radiative cooling and star-formation can reduce the amplitude of the kSZ power spectrum by up to 33%, or 1 uK^2 at ell = 3000. This is driven by a decrease in the mean gas density in groups and clusters due to the conversion of gas into stars. Variations in the redshifts at which helium reionization occurs can effect the amplitude by a similar fraction, while current constraints on cosmological parameters (namely sigma_8) translate to a further +-15% uncertainty on the kSZ power spectrum. We demonstrate how the models presented in this work can be constrained -- reducing the astrophysical uncertainty on the kSZ signal -- by measuring the redshift dependence of the signal via kSZ tomography. Finally, we discuss how the results of this work can help constrain the duration of reionization via measurements of the kinetic SZ signal sourced by inhomogeneous (or patchy) reionization.

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

One of the central goals of multi-wavelength galaxy cluster cosmology is to unite all cluster observables to form a consistent understanding of cluster mass. Here, we study the impact of systematic effects from optical cluster catalogs on stacked SZ signals. We show that the optically predicted Y-decrement can vary by as much as 50% based on the current 2 sigma systematic uncertainties in the observed mass-richness relationship. Mis-centering and impurities will suppress the SZ signal compared to expectations for a clean and perfectly centered optical sample, but to a lesser degree. We show that the level of these variations and suppression is dependent on the amount of systematics in the optical cluster catalogs. We also study X-ray luminosity-dependent sub-sampling of the optical catalog and find that it creates Malmquist bias increasing the observed Y-decrement of the stacked signal. We show that the current Planck measurements of the Y-decrement around SDSS optical clusters and their X-ray counterparts are consistent with expectations after accounting for the 1 sigma optical systematic uncertainties using the Johnston mass richness relation.

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

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

In this work, we examine the effects of mergers on the hydrostatic mass estimate of galaxy clusters using high-resolution Eulerian cosmological simulations. We utilize merger trees to isolate the last merger for each cluster in our sample and follow the time evolution of the hydrostatic mass bias as the systems relax. We find that during a merger, a shock propagates outward from the parent cluster, resulting in an overestimate in the hydrostatic mass bias. After the merger, as a cluster relaxes, the bias in hydrostatic mass estimate decreases but remains at a level of -5-10% with 15-20% scatter within r500. We also investigate the post-merger evolution of the pressure support from bulk motions, a dominant cause of this residual mass bias. At r500, the contribution from random motions peaks at 30% of the total pressure during the merger and quickly decays to \sim 10-15% as a cluster relaxes. Additionally, we use a measure of the random motion pressure to correct the hydrostatic mass estimate. We discover that 4 Gyr after mergers, the direct effects of the merger event on the hydrostatic mass bias have become negligible. Thereafter, the mass bias is primarily due to residual bulk motions in the gas which are not accounted for in the hydrostatic equilibrium equation. We present a hydrostatic mass bias correction method that can recover the unbiased cluster mass for relaxed clusters with 9% scatter at r500 and 11% scatter in the outskirts, within r200.

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

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

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

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

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

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

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

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

Using the latest physical modeling and constrained by the most recent data, we develop a phenomenological parameterized model of the contributions to intensity and polarization maps at millimeter wavelengths from external galaxies and Sunyaev-Zeldovich effects. We find such modeling to be necessary for estimation of cosmological parameters from Planck data. For example, ignoring the clustering of the infrared background would result in a bias in n_s of 7 sigma. We show that the simultaneous marginalization over a full foreground model can eliminate such biases, while increasing the statistical uncertainty in cosmological parameters by less than 20%. The small increases in uncertainty can be significantly reduced with the inclusion of higher-resolution ground-based data. The multi-frequency analysis we employ involves modeling 46 total power spectra and marginalization over 17 foreground parameters. We show that we can also reduce the data to a best estimate of the CMB power spectra, and just two principal components (with constrained amplitudes) describing residual foreground contamination.

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

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

We use an analytic model to investigate the theoretical uncertainty on the thermal Sunyaev-Zel'dovich (SZ) power spectrum due to astrophysical uncertainties in the thermal structure of the intracluster medium. Our model accounts for star formation and energy feedback (from supernovae and active galactic nuclei) as well as radially dependent non-thermal pressure support due to random gas motions, the latter calibrated by recent hydrodynamical simulations. We compare the model against X-ray observations of low redshift clusters, finding excellent agreement with observed pressure profiles. Varying the levels of feedback and non-thermal pressure support can significantly change both the amplitude and shape of the thermal SZ power spectrum. Increasing the feedback suppresses power at small angular scales, shifting the peak of the power spectrum to lower ell. On the other hand, increasing the non-thermal pressure support has the opposite effect, significantly reducing power at large angular scales. In general, including non-thermal pressure at the level measured in simulations has a large effect on the power spectrum, reducing the amplitude by 50% at angular scales of a few arcminutes compared to a model without a non-thermal component. Our results demonstrate that measurements of the shape of the power spectrum can reveal useful information on important physical processes in groups and clusters, especially at high-redshift where there exists little observational data. Comparing with the recent South Pole Telescope measurements of the small-scale cosmic microwave background power spectrum, we find our model reduces the tension between the values of sigma_8 measured from the SZ power spectrum and from cluster abundances.

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

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

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

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

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

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

We investigate the impact of non-Gaussian scatter in the cluster mass-observable scaling relation on the mass and redshift distribution of clusters detected by wide area surveys. We parameterize non-Gaussian scatter by incorporating the third and forth moments (skewness and kurtosis) into the distribution P(Mobs|M). We demonstrate that for low scatter mass proxies the higher order moments do not significantly affect the observed cluster mass and redshift distributions. However, for high scatter mass indicators it is necessary for the survey limiting mass threshold to be less than 10^14 h^-1 Msol to prevent the skewness from having a significant impact on the observed number counts, particularly at high redshift. We also show that an unknown level of non-Gaussianity in the scatter is equivalent to an additional uncertainty on the variance in P(Mobs|M) and thus may limit the constraints that can be placed on the dark energy equation of state parameter w. Furthermore, positive skewness flattens the mass function at the high mass end, and so one must also account for skewness in P(Mobs|M) when using the shape of the mass function to constrain cluster scaling-relations.

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

Using both halo model calculations and a large sample of simulated SZ maps, we demonstrate that high-mass clusters add significant non-Gaussian variance to measurements of the SZ power spectrum amplitude. The difficulty in correctly accounting theoretically for the contribution of these objects to the uncertainty in C_l leads to a reduced sensitivity to \sigma_8. We show that a simple solution is to mask out the brightest clusters in the map before measuring the power spectrum. We demonstrate that fairly conservative masking can reduce the variance and Gaussianize the statistics significantly, thus increasing the sensitivity to cosmological parameters. Choosing which objects to mask is non-trivial; we found that using a fixed sky density produced a well-defined and well-behaved estimate that can easily be applied to real maps. For example, masking the 10 (90) brightest clusters in a 100 deg^2 SZ map will improve the sensitivity to C_l by a factor of two at l = 1000 (2000) and 1.5 at l = 2000 (4000). We show that even in the presence of astrophysical foregrounds (primary CMB and point sources) and instrument noise, one can increase the precision on measurements of \sigma_8 by masking up to 0.9 clusters/deg^2.

We present a detailed analysis of the intrinsic scatter in the integrated SZ effect - cluster mass (Y-M) relation, using semi-analytic and simulated cluster samples. Specifically, we investigate the impact on the Y-M relation of energy feedback, variations in the host halo concentration and substructure populations, and projection effects due to unresolved clusters along the line of sight (the SZ background). Furthermore, we investigate at what radius (or overdensity) one should measure the integrated SZE and define cluster mass so as to achieve the tightest possible scaling. We find that the measure of Y with the least scatter is always obtained within a smaller radius than that at which the mass is defined; e.g. for M_{200} (M_{500}) the scatter is least for Y_{500} (Y_{1100}). The inclusion of energy feedback in the gas model significantly increases the intrinsic scatter in the Y-M relation due to larger variations in the gas mass fraction compared to models without feedback. We also find that variations in halo concentration for clusters of a given mass may partly explain why the integrated SZE provides a better mass proxy than the central decrement. Substructure is found to account for approximately 20% of the observed scatter in the Y-M relation. Above M_{200} = 2x10^{14} h^{-1} msun, the SZ background does not significantly effect cluster mass measurements; below this mass, variations in the background signal reduce the optimal angular radius within which one should measure Y to achieve the tightest scaling with M_{200}.

Using the full, three-dimensional potential of galaxy cluster halos (drawn from an N-body simulation of the current, most favored cosmology), the distribution of the X-ray emitting gas is found by assuming a polytropic equation of state and hydrostatic equilibrium, with constraints from conservation of energy and pressure balance at the cluster boundary. The resulting properties of the gas for these simulated redshift zero clusters (the temperature distribution, mass-temperature and luminosity-temperature relations, and the gas fraction) are compared with observations in the X-ray of nearby clusters. The observed properties are reproduced only under the assumption that substantial energy injection from non-gravitational sources has occurred. Our model does not specify the source, but star formation and AGN may be capable of providing this energy, which amounts to 3 to 5 x10^{-5} of the rest mass in stars (assuming ten percent of the gas initially in the cluster forms stars). With the method described here it is possible to generate realistic X-ray and Sunyaev-Zel'dovich cluster maps and catalogs from N-body simulations, with the distributions of internal halo properties (and their trends with mass, location, and time) taken into account.

We present a new definition of subhalos in dissipationless dark matter N-body simulations, based on the coherent identification of their dynamically bound constituents. Whereas previous methods of determining the energetically bound components of a subhalo ignored the contribution of all the remaining particles in the halo (those not geometrically or dynamically associated with the subhalo), our method allows for all the forces, both internal and external, exerted on the subhalo. We demonstrate, using the output of a simulation at different timesteps, that our new method is more accurate at identifying the bound mass of a subhalo. We then compare our new method to previously adopted means of identifying subhalos by applying each to a sample of 1838 virialized halos extracted from a high resolution cosmological simulation. We find that the subhalo distributions are similar in each case, and that the increase in the binding energy of a subhalo from including all the particles located within it is almost entirely balanced by the losses due to the external forces; the net increase in the mass fraction of subhalos is roughly 10%, and the extra substructures tending to reside in the inner parts of the system. Finally, we compare the subhalo populations of halos to the sub-subhalo populations of subhalos, finding the two distributions to be similar. This is a new and interesting result, suggesting a self-similarity in the hierarchy substructures within cluster mass halos.

We have identified over 2000 well resolved cluster halos, and also their associated bound subhalos, from the output of a 1024^3 particle cosmological N-body simulation (of box size 320 h^-1 Mpc and softening length 3.2 h^-1 kpc). This has allowed us to measure halo quantities in a statistically meaningful way, and for the first time analyse their distribution for a large and well resolved sample. We characterize each halo in terms of its morphology, concentration, spin, circular velocity and the fraction of their mass in substructure. We also identify those halos that have not yet reached a state of dynamical equilibrium using the virial theorem with an additional correction to account for the surface pressure at the boundary. These amount to 3.4% of our initial sample. For the virialized halos, we find a median of 5.6% of halo mass is contained within substructure, with the distribution ranging between no identified subhalos to 65%. The fraction of mass in substructure increases with halo mass with logarithmic slope of 0.44 +- 0.06. Halos tend to have a prolate morphology, becoming more so with increasing mass. Subhalos have a greater orbital angular momentum per unit mass than their host halo. Furthermore, their orbital angular momentum is typically well aligned with that of their host. Overall, we find that dimensionless properties of dark matter halos do depend on their mass, thereby demonstrating a lack of self-similarity.

We present an algorithm which is designed to allow the efficient identification and preliminary dynamical analysis of thousands of structures and substructures in large N-body simulations. First we utilise a refined density gradient system (based on DENMAX) to identify the structures, and then apply an iterative approximate method to identify unbound particles, allowing fast calculation of bound substructures. After producing a catalog of separate energetically bound substructures we check to see which of these are energetically bound to adjacent substructures. For such bound complex subhalos, we combine components and check if additional free particles are also bound to the union, repeating the process iteratively until no further changes are found. Thus our subhalos can contain more than one density maximum, but the scheme is stable: starting with a small smoothing length initially produces small structures which must be combined later, and starting with a large smoothing length produces large structures within which sub-substructure is found. We apply this algorithm to three simulations. Two which are using the TPM algorithm by Bode et al. (2000) and one on a simulated halo by Diemand (2004). For all these halos we find about 5-8% of the mass in substructures.