In order to extract cosmological information from observations of the millimeter and submillimeter sky, foreground components must first be removed to produce an estimate of the cosmic microwave background (CMB). We developed a machine learning approach for doing so for full-sky temperature maps of the millimeter and submillimeter sky. We constructed a Bayesian spherical convolutional neural network architecture to produce a model that captures both spectral and morphological aspects of the foregrounds. Additionally, the model outputs a per-pixel error estimate that incorporates both statistical and model uncertainties. The model was then trained using simulations that incorporated knowledge of these foreground components that was available at the time of the launch of the Planck satellite. On simulated maps, the CMB is recovered with a mean absolute difference of $<4\mu$K over the full sky after masking map pixels with a predicted standard error of $>50\mu$K; the angular power spectrum is also accurately recovered. Once validated with the simulations, this model was applied to Planck temperature observations from its 70GHz through 857GHz channels to produce a foreground-cleaned CMB map at a Healpix map resolution of NSIDE=512. Furthermore, we demonstrate the utility of the technique for evaluating how well different simulations match observations, particularly in regard to the modeling of thermal dust.

We present a catalog of 510 radio-loud active galactic nuclei (AGN, primarily blazars) and 287 dusty star-forming galaxies (DSFGs) detected by the Atacama Cosmology Telescope at $>5\sigma$ significance in bands centered on 148 GHz (2 mm), 218 GHz (1.4 mm) and 277 GHz (1.1 mm), from a 480 square degrees strip on the celestial equator with additional (360 square degrees) shallower fields. Combining the deepest available 218 GHz wide-field imaging, 277 GHz data, and multi-band filtering yields the most sensitive wide-field millimeter-wave DSFG selection to date with rms noise referenced to 218 GHz reaching $<2$ mJy. We developed techniques to remove Galactic contamination from the extragalactic catalog, yielding 321 additional Galactic sources. We employ a new flux debiasing method that handles the heterogeneous sample selection due to Galactic cuts. We present spectral properties and source counts of the AGN and DSFGs. The DSFG spectra depart from an optically thin modified blackbody between 218 GHz and 277 GHz, consistent with optically thick emission or an additional cold dust component. For bright AGN, the inter-year RMS fractional deviation in flux density from source variability is $\sim40\%$. We report 8$-$2870 mJy source counts for AGN and 8$-$90 mJy source counts for DSFGs, the latter probing both the brighter, lensed population and the fainter, unlensed population. At 277 GHz we report the first source counts measurements at these flux densities, finding an excess above most model count predictions. Finally, we select thirty of the brightest DSFGs for multi-frequency study as candidate high-$z$ lensed systems.

We explore the fundamental limits to which reionization histories can be constrained using only large-scale cosmic microwave background (CMB) anisotropy measurements. The redshift distribution of the fractional ionization $x_e(z)$ affects the angular distribution of CMB polarization. We project constraints on the reionization history of the universe using low-noise full-sky temperature and E-mode measurements of the CMB. We show that the measured TE power spectrum, $\hat C_\ell^\mathrm{TE}$, has roughly one quarter of the constraining power of $\hat C_\ell^\mathrm{EE}$ on the reionization optical depth $\tau$, and its addition improves the precision on $\tau$ by 20% over using $\hat C_\ell^\mathrm{EE}$ only. We also use a two-step reionization model with an additional high redshift step, parametrized by an early ionization fraction $x_e^\mathrm{min}$, and a late reionization step at $z_\mathrm{re}$. We find that future high signal-to-noise measurements of the multipoles $10\leqslant\ell<20$ are especially important for breaking the degeneracy between $x_e^\mathrm{min}$ and $z_\mathrm{re}$. In addition, we show that the uncertainties on these parameters determined from a map with sensitivity $10\,\mathrm{\mu K\,arcmin}$ are less than 5% larger than the uncertainties in the noiseless case, making this noise level a natural target for future large sky area E-mode measurements.

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (covering approximately 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel'dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. The tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper will be made available on LAMBDA.

We report circular polarization measurements from the first two years of observation with the 40 GHz polarimeter of the Cosmology Large Angular Scale Surveyor (CLASS). CLASS is conducting a multi-frequency survey covering 75% of the sky from the Atacama Desert designed to measure the cosmic microwave background (CMB) linear E and B polarization on angular scales $1^\circ \lesssim \theta \leq 90^\circ$, corresponding to a multipole range of $2 \leq \ell \lesssim 200$. The modulation technology enabling measurements of linear polarization at the largest angular scales from the ground, the Variable-delay Polarization Modulator, is uniquely designed to provide explicit sensitivity to circular polarization (Stokes $V$). We present a first detection of circularly polarized atmospheric emission at 40 GHz that is well described by a dipole with an amplitude of $124\pm4\,\mathrm{\mu K}$ when observed at an elevation of $45^\circ$, and discuss its potential impact as a foreground to CMB experiments. Filtering the atmospheric component, CLASS places a 95% C.L. upper limit of $0.4\,\mathrm{\mu K}^2$ to $13.5\,\mathrm{\mu K}^2$ on $\ell(\ell+1)C_\ell^{VV}/(2\pi)$ between $1 \leq \ell \leq 120$, representing a two-orders-of-magnitude improvement over previous limits.

In this paper we probe the hot, post-shock gas component of quasar-driven winds through the thermal Sunyaev-Zel'dovich (tSZ) effect. Combining datasets from the Atacama Cosmology Telescope, the $\textit{Herschel}$ Space Observatory, and the Very Large Array, we measure average spectral energy distributions (SEDs) of 109,829 optically-selected, radio quiet quasars from 1.4~GHz to 3000~GHz in six redshift bins between $0.3<z<3.5$. We model the emission components in the radio and far-infrared, plus a spectral distortion from the tSZ effect. At $z>1.91$, we measure the tSZ effect at $3.8\sigma$ significance with an amplitude corresponding to a total thermal energy of $3.1\times10^{60}$ ergs. If this energy is due to virialized gas, then our measurement implies quasar host halo masses are $\sim6\times10^{12}~h^{-1}$M$_\odot$. Alternatively, if the host dark matter halo masses are $\sim2\times10^{12}~h^{-1}$M$_\odot$ as some measurements suggest, then we measure a $>$90 per cent excess in the thermal energy over that expected due to virialization. If the measured SZ effect is primarily due to hot bubbles from quasar-driven winds, we find that $(5^{+1.2}_{-1.3}$) per cent of the quasar bolometric luminosity couples to the intergalactic medium over a fiducial quasar lifetime of 100 Myr. An additional source of tSZ may be correlated structure, and further work is required to separate the contributions. At $z\leq1.91$, we detect emission at 95 and 148~GHz that is in excess of thermal dust and optically thin synchrotron emission. We investigate potential sources of this excess emission, finding that CO line emission and an additional optically thick synchrotron component are the most viable candidates.

Consistency checks of cosmological data sets are an important tool because they may suggest systematic errors or the type of modifications to \lcdm\ necessary to resolve current tensions. In this work, we derive an analytic method for calculating the level of correlations between model parameters from two correlated cosmological data sets, which complements more computationally expensive simulations. This method is an extension of the Fisher analysis that assumes a Gaussian likelihood and a known data covariance matrix. We apply this method to the SPTpol temperature and polarization CMB spectra (TE and EE). We find weak correlations between \lcdm\ parameters with a 9$\%$ correlation between the TE-only and EE-only constraints on $H_0$ and a 25$\%$ and 32$\%$ correlation for log($A_s$) and $n_s$ respectively. Despite the negative correlations between the TE and EE power spectra, the correlations in the parameters are positive. The TE-EE parameter differences are consistent with zero, with a PTE of 0.53, in contrast to the PTE of 0.017 reported by SPTpol for the consistency of the TE and EE power spectra with \lcdm. Using simulations we find that the results of these two tests are independent and that this difference can arise simply from statistical fluctuations. Ignoring correlations in the TT-TE and TE-EE comparisons biases the $\chi^2$ low, artificially making parameters look more consistent. Therefore, we conclude that these correlations need to be accounted for when performing internal consistency checks of the TT vs TE vs EE power spectra for future CMB analyses.

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

Angular power spectra computed from Planck HFI 353 GHz intensity and polarization maps produce a TB correlation that can be approximated by a power-law (Planck Collaboration XI 2018). Whether the observed TB correlation is an induced systematic feature or a physical property of Galactic dust emission is of interest both for cosmological and Galactic studies. We investigate the large angular scale E and B-mode morphology of microwave polarized thermal dust emission, and relate it to physical quantities of polarization angle and polarization fraction. We use empirical models of polarized dust to show that dust polarization angle is a key factor in producing the TB correlation. A small sample of both simulated and observed polarization angle maps are combined with 353 GHz intensity and dust polarization fraction to produce a suite of maps from which we compute TB and EB. Model realizations that produce a positive TB correlation are common and can result from large-scale (>5 degree) structure in the polarization angle. The TB correlation appears robust to introduction of individual intensity, polarization angle and polarization fraction model components that are independent of the 353 GHz observations. We conclude that the observed TB correlation is likely the result of large-scale Galactic dust polarization properties.

This is an overview of the data products and other resources available through NASA's LAMBDA site https://lambda.gsfc.nasa.gov/. An up-to-date version of this document, along with code tools actively maintained and developed by LAMBDA staff, can be found on the LAMBDA GitHub page at https://github.com/nasa-lambda/lambda_overview. New data products and other updates are announced on LAMBDA's twitter account at https://twitter.com/NASA_LAMBDA. If you have questions or suggestions relating to LAMBDA, or are interested in joining a LAMBDA advisory group, please contact us using the form here: https://lambda.gsfc.nasa.gov/contact/contact.cfm.

We perform forecasts for how baryon acoustic oscillation (BAO) scale and redshift-space distortion (RSD) measurements from future spectroscopic emission line galaxy (ELG) surveys such as Euclid are degraded in the presence of spectral line misidentification. Using analytic calculations verified with mock galaxy catalogs from log-normal simulations we find that constraints are degraded in two ways, even when the interloper power spectrum is modeled correctly in the likelihood. Firstly, there is a loss of signal-to-noise ratio for the power spectrum of the target galaxies, which propagates to all cosmological constraints and increases with contamination fraction, $f_c$. Secondly, degeneracies can open up between $f_c$ and cosmological parameters. In our calculations this typically increases BAO scale uncertainties at the 10-20% level when marginalizing over parameters determining the broadband power spectrum shape. External constraints on $f_c$, or parameters determining the shape of the power spectrum, for example from cosmic microwave background (CMB) measurements, can remove this effect. There is a near-perfect degeneracy between $f_c$ and the power spectrum amplitude for low $f_c$ values, where $f_c$ is not well determined from the contaminated sample alone. This has the potential to strongly degrade RSD constraints. The degeneracy can be broken with an external constraint on $f_c$, for example from cross-correlation with a separate galaxy sample containing the misidentified line, or deeper sub-surveys.

Current cosmological tensions motivate investigating extensions to the standard $\Lambda$CDM model. Additional model parameters are typically varied one or two at a time, in a series of separate tests. The purpose of this paper is to highlight that information is lost by not also examining the correlations between these additional parameters, which arise when their effects on model predictions are similar, even if the parameters are not varied simultaneously. We show how these correlations can be quantified with simulations and Markov Chain Monte Carlo (MCMC) methods. As an example, we assume that $\Lambda$CDM is the true underlying model, and calculate the correlations expected between the phenomenological lensing amplitude parameter, $A_L$, the running of the spectral index, $n_{\rm run}$, and the primordial helium mass fraction, $Y_P$, when these parameters are varied one at a time along with the $\Lambda$CDM parameters in fits to the $\textit{Planck}$ 2015 temperature power spectrum. These correlations are not small, ranging from 0.31 ($A_L-n_{\rm run}$) to $-0.93$ ($n_{\rm run}-Y_P$). We find that the values of these three parameters from the $\textit{Planck}$ data are consistent with $\Lambda$CDM expectations within $0.9\sigma$ when the correlations are accounted for. This does not explain the 1.8-2.7$\sigma$ $\textit{Planck}$ preference for $A_L>1$, but provides an additional $\Lambda$CDM consistency test. For example, if $A_L>1$ was a symptom of an underlying systematic error or some real but unknown physical effect that also produced spurious correlations with $n_{\rm run}$ or $Y_P$ our test might have revealed this. We recommend that future cosmological analyses examine correlations between additional model parameters in addition to investigating them separately, one a time.

The reionization optical depth is the most poorly determined of the six $\Lambda$CDM parameters fit to CMB anisotropy data. Instrumental noise and systematics have prevented uncertainties from reaching their cosmic variance limit. At present, the datasets providing the most statistical constraining power are the WMAP, Planck LFI, and Planck HFI full-sky polarization maps. As the reprocessed HFI data with reduced systematics are not yet publicly unavailable, we examine determinations of $\tau$ using 9-year WMAP and 2015 Planck LFI data, with an emphasis on characterizing potential systematic bias resulting from foreground template and masking choices. We find evidence for a low-level systematic in the LFI polarization data with a roughly common-mode morphology across the LFI frequencies and a spectrum consistent with leakage of intensity signal into the polarization channels. We demonstrate significant bias in the optical depth derived when using the LFI 30 GHz map as a template to clean synchrotron from WMAP data, and recommend against use of the 2015 LFI 30 GHz polarization data as a foreground template for non-LFI datasets. We find an inconsistency between versions of the 2015 polarized 353 GHz dust templates reconstructed from the Planck likelihood and those from delivered maps, which can affect $\tau$ at the 1$\sigma$ level. The spread in $\tau$ values over the ensemble of data combinations we study suggests that systematic uncertainties still contribute significantly to the current uncertainty in $\tau$, but all values are consistent with the range of $\tau$ = 0.07 +/- 0.02.

We perform a comparison of WMAP 9-year (WMAP9) and Planck 2015 cosmic microwave background (CMB) temperature power spectra across multipoles $30\leq\ell\leq1200$. We generate simulations to estimate the correlation between the two datasets due to cosmic variance from observing the same sky. We find that their spectra are consistent within $1\sigma$. While we do not implement the optimal "$C^{-1}$" estimator on WMAP maps as in the WMAP9 analysis, we demonstrate that the change of pixel weighting only shifts our results at most at the $0.66\sigma$ level. We also show that changing the fiducial power spectrum for simulations only impacts the comparison at around $0.1\sigma$ level. We exclude $\ell<30$ both because WMAP9 data were included in the Planck 2015 $\ell<30$ analysis, and because the cosmic variance uncertainty on these scales is large enough that any remaining systematic difference between the experiments is extremely unlikely to affect cosmological constraints. The consistency shown in our analysis provides high confidence in both the WMAP9 temperature power spectrum and the overlapping multipole region of Planck 2015's, virtually independent of any assumed cosmological model. Our results indicate that cosmological model differences between Planck and WMAP do not arise from measurement differences, but from the high multipoles not measured by WMAP.

We revisit the degeneracy between the Hubble constant, $H_0$, and matter density, $\Omega_m$, for current cosmic microwave background (CMB) observations within the standard $\Lambda CDM$ model. We show that Planck, Wilkinson Microwave Anisotropy Probe (WMAP), South Pole Telescope (SPT), and Atacama Cosmology Telescope Polarimeter (ACTPol) temperature power spectra produce different values of the exponent $x$ from minimizing the variance of the product $\Omega_mH_0^x$. The distribution of $x$ from the different data sets does not follow the Markov Chain Monte Carlo (MCMC) best-fit values for $H_0$ or $\Omega_m$. Particularly striking is the difference between Planck multipoles $\ell\leq800$ ($x=2.81$), and WMAP ($x = 2.94$), despite very similar best-fit cosmologies. We use a Fisher matrix analysis to show that, in fact, this range in exponent values is exactly as expected in $\Lambda CDM$ given the multipole coverage and power spectrum uncertainties for each experiment. We show that the difference in $x$ from the Planck $\ell \leq 800$ and WMAP data is explained by a turning point in the relationship between $x$ and the maximum effective multipole, at around $\ell=700$. The value of $x$ is determined by several physical effects, and we highlight the significant impact of gravitational lensing for the high-multipole measurements. Despite the spread of $H_0$ values from different CMB experiments, the experiments are consistent with their sampling of the $\Omega_m-H_0$ degeneracy and do not show evidence for the need for new physics or for the presence of significant underestimated systematics according to these tests. The Fisher calculations can be used to predict the $\Omega_m-H_0$ degeneracy of future experiments.

We examine the impact of baryon acoustic oscillation (BAO) scale measurements on the discrepancy between the value of the Hubble constant ($H_0$) inferred from the local distance ladder and from Planck cosmic microwave background (CMB) data. While the BAO data alone cannot constrain $H_0$, we show that combining the latest BAO results with WMAP, Atacama Cosmology Telescope (ACT), or South Pole Telescope (SPT) CMB data produces values of $H_0$ that are $2.4-3.1\sigma$ lower than the distance ladder, independent of Planck, and that this downward pull was less apparent in some earlier analyses that used only angle-averaged BAO scale constraints rather than full anisotropic information. At the same time, the combination of BAO and CMB data also disfavors the lower values of $H_0$ preferred by the Planck high-multipole temperature power spectrum. Combining galaxy and Lyman-$\alpha$ forest (Ly$\alpha$) BAO with a precise estimate of the primordial deuterium abundance produces $H_0=66.98\pm1.18$ km s$^{-1}$ Mpc$^{-1}$ for the flat $\Lambda$CDM model. This value is completely independent of CMB anisotropy constraints and is $3.0\sigma$ lower than the latest distance ladder constraint, although $2.4\sigma$ tension also exists between the galaxy BAO and Ly$\alpha$ BAO. These results show that it is not possible to explain the $H_0$ disagreement solely with a systematic error specific to the Planck data. The fact that tensions remain even after the removal of any single data set makes this intriguing puzzle all the more challenging to resolve.

We report a measurement of the power spectrum of cosmic microwave background (CMB) lensing from two seasons of Atacama Cosmology Telescope Polarimeter (ACTPol) CMB data. The CMB lensing power spectrum is extracted from both temperature and polarization data using quadratic estimators. We obtain results that are consistent with the expectation from the best-fit Planck LCDM model over a range of multipoles L=80-2100, with an amplitude of lensing A_lens = 1.06 +/- 0.15 (stat.) +/- 0.06 (sys.) relative to Planck. Our measurement of the CMB lensing power spectrum gives sigma_8 Omega_m^0.25 = 0.643 +/- 0.054; including baryon acoustic oscillation scale data, we constrain the amplitude of density fluctuations to be sigma_8 = 0.831 +/- 0.053. We also update constraints on the neutrino mass sum. We verify our lensing measurement with a number of null tests and systematic checks, finding no evidence of significant systematic errors. This measurement relies on a small fraction of the ACTPol data already taken; more precise lensing results can therefore be expected from the full ACTPol dataset.

We present the temperature and polarization angular power spectra measured by the Atacama Cosmology Telescope Polarimeter (ACTPol). We analyze night-time data collected during 2013-14 using two detector arrays at 149 GHz, from 548 deg$^2$ of sky on the celestial equator. We use these spectra, and the spectra measured with the MBAC camera on ACT from 2008-10, in combination with Planck and WMAP data to estimate cosmological parameters from the temperature, polarization, and temperature-polarization cross-correlations. We find the new ACTPol data to be consistent with the LCDM model. The ACTPol temperature-polarization cross-spectrum now provides stronger constraints on multiple parameters than the ACTPol temperature spectrum, including the baryon density, the acoustic peak angular scale, and the derived Hubble constant. Adding the new data to planck temperature data tightens the limits on damping tail parameters, for example reducing the joint uncertainty on the number of neutrino species and the primordial helium fraction by 20%.

We examine the internal consistency of the Planck 2015 cosmic microwave background (CMB) temperature anisotropy power spectrum. We show that tension exists between cosmological constant cold dark matter (LCDM) model parameters inferred from multipoles l<1000 (roughly those accessible to Wilkinson Microwave Anisotropy Probe), and from l>=1000, particularly the CDM density, Omega_ch^2, which is discrepant at 2.5 sigma for a Planck-motivated prior on the optical depth, tau=0.07+/-0.02. We find some parameter tensions to be larger than previously reported because of inaccuracy in the code used by the Planck Collaboration to generate model spectra. The Planck l>=1000 constraints are also in tension with low-redshift data sets, including Planck's own measurement of the CMB lensing power spectrum (2.4 sigma), and the most precise baryon acoustic oscillation (BAO) scale determination (2.5 sigma). The Hubble constant predicted by Planck from l>=1000, H_0=64.1+/-1.7 km/s/Mpc, disagrees with the most precise local distance ladder measurement of 73.0+/-2.4 km/s/Mpc at the 3.0 sigma level, while the Planck value from l<1000, 69.7+/-1.7 km/s/Mpc, is consistent within 1 sigma. A discrepancy between the Planck and South Pole Telescope (SPT) high-multipole CMB spectra disfavors interpreting these tensions as evidence for new physics. We conclude that the parameters from the Planck high-multipole spectrum probably differ from the underlying values due to either an unlikely statistical fluctuation or unaccounted-for systematics persisting in the Planck data.

We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2{\sigma}, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements.

We present a measurement of the gravitational lensing of the Cosmic Microwave Background (CMB) temperature and polarization fields obtained by cross-correlating the reconstructed convergence signal from the first season of ACTPol data at 146 GHz with Cosmic Infrared Background (CIB) fluctuations measured using the Planck satellite. Using an overlap area of 206 square degrees, we detect gravitational lensing of the CMB polarization by large-scale structure at a statistical significance of 4.5 sigma. Combining both CMB temperature and polarization data gives a lensing detection at 9.1 sigma significance. A B-mode polarization lensing signal is present with a significance of 3.2 sigma. We also present the first measurement of CMB lensing--CIB correlation at small scales corresponding to l > 2000. Null tests and systematic checks show that our results are not significantly biased by astrophysical or instrumental systematic effects, including Galactic dust. Fitting our measurements to the best-fit lensing-CIB cross power spectrum measured in Planck data, scaled by an amplitude A, gives A=1.02 +0.12/-0.18 (stat.) +/-0.06(syst.), consistent with the Planck results.

We present a measurement of the one-point probability distribution function (PDF) of the thermal Sunyaev-Zel'dovich (tSZ) decrement in the pixel temperature histogram of filtered 148 GHz sky maps from the Atacama Cosmology Telescope (ACT). The PDF includes the signal from all galaxy clusters in the map, including objects below the signal-to-noise threshold for individual detection, making it a particularly sensitive probe of the amplitude of matter density perturbations, $\sigma_8$. We use a combination of analytic halo model calculations and numerical simulations to compute the theoretical tSZ PDF and its covariance matrix, accounting for all noise sources and including relativistic corrections. From the measured ACT 148 GHz PDF alone, we find $\sigma_8 = 0.793 \pm 0.018$, with additional systematic errors of $\pm 0.017$ due to uncertainty in intracluster medium gas physics and $\pm 0.006$ due to uncertainty in infrared point source contamination. Using effectively the same data set, the statistical error here is a factor of two lower than that found in ACT's previous $\sigma_8$ determination based solely on the skewness of the tSZ signal. In future temperature maps with higher sensitivity, the tSZ PDF will break the degeneracy between intracluster medium gas physics and cosmological parameters.

We present a statistical analysis of the millimeter-wavelength properties of 1.4 GHz-selected sources and a detection of the Sunyaev-Zel'dovich (SZ) effect associated with the halos that host them. The Atacama Cosmology Telescope (ACT) has conducted a survey at 148 GHz, 218 GHz and 277 GHz along the celestial equator. Using samples of radio sources selected at 1.4 GHz from FIRST and NVSS, we measure the stacked 148, 218 and 277 GHz flux densities for sources with 1.4 GHz flux densities ranging from 5 to 200 mJy. At these flux densities, the radio source population is dominated by active galactic nuclei (AGN), with both steep and flat spectrum populations, which have combined radio-to-millimeter spectral indices ranging from 0.5 to 0.95, reflecting the prevalence of steep spectrum sources at high flux densities and the presence of flat spectrum sources at lower flux densities. The thermal SZ effect associated with the halos that host the AGN is detected at the 5$\sigma$ level through its spectral signature. When we compare the SZ effect with weak lensing measurements of radio galaxies, we find that the relation between the two is consistent with that measured by Planck for local bright galaxies. We present a detection of the SZ effect in some of the lowest mass halos (average $M_{200}\approx10^{13}$M$_{\odot}h_{70}^{-1}$) studied to date. This detection is particularly important in the context of galaxy evolution models, as it confirms that galaxies with radio AGN also typically support hot gaseous halos. With Herschel observations, we show that the SZ detection is not significantly contaminated by dust. We show that 5 mJy$<S_{1.4}<$200 mJy radio sources contribute $\ell(\ell+1)C_{\ell}/(2\pi)=0.37\pm0.03\mu$K$^2$ to the angular power spectrum at $\ell=3000$ at 148 GHz, after accounting for the SZ effect associated with their host halos.

We report on measurements of the cosmic microwave background (CMB) and celestial polarization at 146 GHz made with the Atacama Cosmology Telescope Polarimeter (ACTPol) in its first three months of observing. Four regions of sky covering a total of 270 square degrees were mapped with an angular resolution of $1.3'$. The map noise levels in the four regions are between 11 and 17 $\mu$K-arcmin. We present TT, TE, EE, TB, EB, and BB power spectra from three of these regions. The observed E-mode polarization power spectrum, displaying six acoustic peaks in the range $200<\ell<3000$, is an excellent fit to the prediction of the best-fit cosmological models from WMAP9+ACT and Planck data. The polarization power spectrum, which mainly reflects primordial plasma velocity perturbations, provides an independent determination of cosmological parameters consistent with those based on the temperature power spectrum, which results mostly from primordial density perturbations. We find that without masking any point sources in the EE data at $\ell<9000$, the Poisson tail of the EE power spectrum due to polarized point sources has an amplitude less than $2.4$ $\mu$K$^2$ at $\ell = 3000$ at 95\% confidence. Finally, we report that the Crab Nebula, an important polarization calibration source at microwave frequencies, has 8.7\% polarization with an angle of $150.7^\circ \pm 0.6^\circ$ when smoothed with a $5'$ Gaussian beam.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

We present a catalog of 191 extragalactic sources detected by the Atacama Cosmology Telescope (ACT) at 148 GHz and/or 218 GHz in the 2008 Southern survey. Flux densities span 14-1700 mJy, and we use source spectral indices derived using ACT-only data to divide our sources into two sub-populations: 167 radio galaxies powered by central active galactic nuclei (AGN), and 24 dusty star-forming galaxies (DSFGs). We cross-identify 97% of our sources (166 of the AGN and 19 of the DSFGs) with those in currently available catalogs. When combined with flux densities from the Australian Telescope 20 GHz survey and follow-up observations with the Australia Telescope Compact Array, the synchrotron-dominated population is seen to exhibit a steepening of the slope of the spectral energy distribution from 20 to 148 GHz, with the trend continuing to 218 GHz. The ACT dust-dominated source population has a median spectral index of 3.7+0.62-0.86, and includes both local galaxies and sources with redshifts as great as 5.6. Dusty sources with no counterpart in existing catalogs likely belong to a recently discovered subpopulation of DSFGs lensed by foreground galaxies or galaxy groups.

We present the temperature power spectrum of the Cosmic Microwave Background obtained by cross-correlating maps from the Atacama Cosmology Telescope (ACT) at 148 and 218 GHz with maps from the Planck satellite at 143 and 217 GHz, in two overlapping regions covering 592 square degrees. We find excellent agreement between the two datasets at both frequencies, quantified using the variance of the residuals between the ACT power spectra and the ACTxPlanck cross-spectra. We use these cross-correlations to calibrate the ACT data at 148 and 218 GHz, to 0.7% and 2% precision respectively. We find no evidence for anisotropy in the calibration parameter. We compare the Planck 353 GHz power spectrum with the measured amplitudes of dust and cosmic infrared background (CIB) of ACT data at 148 and 218 GHz. We also compare planet and point source measurements from the two experiments.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

We present constraints on cosmological and astrophysical parameters from high-resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and secondary foreground parameters is fit to the map power spectra and lensing deflection power spectrum, including contributions from both the thermal Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, radio sources, and the correlation between the tSZ effect and infrared sources. The power ell^2 C_ell/2pi of the thermal SZ power spectrum at 148 GHz is measured to be 3.4 +\- 1.4 muK^2 at ell=3000, while the corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6 muK^2. Combining ACT power spectra with the WMAP 7-year temperature and polarization power spectra, we find excellent consistency with the LCDM model. We constrain the number of effective relativistic degrees of freedom in the early universe to be Neff=2.79 +\- 0.56, in agreement with the canonical value of Neff=3.046 for three massless neutrinos. We constrain the sum of the neutrino masses to be Sigma m_nu < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to be Yp = 0.225 +\- 0.034, and measure no variation in the fine structure constant alpha since recombination, with alpha/alpha0 = 1.004 +/- 0.005. We also find no evidence for any running of the scalar spectral index, dns/dlnk = -0.004 +\- 0.012.

The Atacama Cosmology Telescope has measured the angular power spectra of microwave fluctuations to arcminute scales at frequencies of 148 and 218 GHz, from three seasons of data. At small scales the fluctuations in the primordial Cosmic Microwave Background (CMB) become increasingly obscured by extragalactic foregounds and secondary CMB signals. We present results from a nine-parameter model describing these secondary effects, including the thermal and kinematic Sunyaev-Zel'dovich (tSZ and kSZ) power; the clustered and Poisson-like power from Cosmic Infrared Background (CIB) sources, and their frequency scaling; the tSZ-CIB correlation coefficient; the extragalactic radio source power; and thermal dust emission from Galactic cirrus in two different regions of the sky. In order to extract cosmological parameters, we describe a likelihood function for the ACT data, fitting this model to the multi-frequency spectra in the multipole range 500<ell<10000. We extend the likelihood to include spectra from the South Pole Telescope at frequencies of 95, 150, and 220 GHz. Accounting for different radio source levels and Galactic cirrus emission, the same model provides an excellent fit to both datasets simultaneously, with chi2/dof= 675/697 for ACT, and 96/107 for SPT. We then use the multi-frequency likelihood to estimate the CMB power spectrum from ACT in bandpowers, marginalizing over the secondary parameters. This provides a simplified `CMB-only' likelihood in the range 500<ell<3500 for use in cosmological parameter estimation.

We constrain cosmological parameters using combined measurements of the baryon acoustic oscillation (BAO) feature in the correlation function of galaxies and Ly-\alpha absorbers that together cover 0.1 < z < 2.4. The BAO position measurements alone -- without fixing the absolute sound horizon `standard ruler' length with cosmic microwave background (CMB) data -- constrain \Omega_m = 0.303 +/- 0.040 (68 per cent confidence) for a flat \Lambda CDM model, and w = -1.06^{+0.33}_{-0.32}, \Omega_m = 0.292^{+0.045}_{-0.040} for a flat wCDM model. Adding other large-scale structure (LSS) clustering constraints -- correlation function shape, the Alcock-Paczynski test and growth rate information -- to the BAO considerably tightens constraints (\Omega_m = 0.290 +/- 0.019, H_0 = 67.5 +/- 2.8 km s^{-1} Mpc^{-1}, \sigma_8 = 0.80 +/- 0.05 for \Lambda CDM, and w = -1.14 +/- 0.19 for wCDM). The LSS data mildly prefer a lower value of H_0, and a higher value of \Omega_m, than local distance ladder and type IA supernovae (SNe) measurements, respectively. While tension in the combined CMB, SNe and distance ladder data appears to be relieved by allowing w < -1, this freedom introduces tension with the LSS \sigma_8 constraint from the growth rate of matter fluctuations. The combined constraint on w from CMB, BAO and LSS clustering for a flat wCDM model is w = -1.03 +/- 0.06.

We perform a joint fit to differential number counts from Spitzer's MIPS and Herschel's SPIRE instruments, and angular power spectra of cosmic infrared background (CIB) anisotropies from SPIRE, Planck, the Atacama Cosmology Telescope, and the South Pole Telescope, which together span 220 < \nu / GHz < 4300 (70 < \lambda / \mu m < 1400). We simultaneously constrain the dust luminosity function, thermal dust spectral energy distribution (SED) and clustering properties of CIB sources, and the evolution of these quantities over cosmic time. We find that the data strongly require redshift evolution in the thermal dust SED. In our adopted parametrization, this evolution takes the form of an increase in graybody dust temperature at high redshift, but it may also be related to a temperature - dust luminosity correlation or evolution in dust opacity. The counts and spectra together constrain the evolution of the thermal dust luminosity function up to z ~ 2.5-3, complementing approaches relying on rest-frame mid-infrared observations of the rarest bright objects. We are able to fit the power spectra without requiring a complex halo model approach, and show that neglecting scale-dependent halo bias may be impairing analyses that do use this framework.

[Abridged] We present a catalog of 68 galaxy clusters, of which 19 are new discoveries, detected via the Sunyaev-Zel'dovich effect (SZ) at 148 GHz in the Atacama Cosmology Telescope (ACT) survey of 504 square degrees on the celestial equator. A subsample of 48 clusters within the 270 square degree region overlapping SDSS Stripe 82 is estimated to be 90% complete for M_500c > 4.5e14 Msun and 0.15 < z < 0.8. While matched filters are used to detect the clusters, the sample is studied further through a "Profile Based Amplitude Analysis" using a single filter at a fixed \theta_500 = 5.9' angular scale. This new approach takes advantage of the "Universal Pressure Profile" (UPP) to fix the relationship between the cluster characteristic size (R_500) and the integrated Compton parameter (Y_500). The UPP scalings are found to be nearly identical to an adiabatic model, while a model incorporating non-thermal pressure better matches dynamical mass measurements and masses from the South Pole Telescope. A high signal to noise ratio subsample of 15 ACT clusters is used to obtain cosmological constraints. We first confirm that constraints from SZ data are limited by uncertainty in the scaling relation parameters rather than sample size or measurement uncertainty. We next add in seven clusters from the ACT Southern survey, including their dynamical mass measurements based on galaxy velocity dispersions. In combination with WMAP7 these data simultaneously constrain the scaling relation and cosmological parameters, yielding \sigma_8 = 0.829 \pm 0.024 and \Omega_m = 0.292 \pm 0.025. The results include marginalization over a 15% bias in dynamical mass relative to the true halo mass. In an extension to LCDM that incorporates non-zero neutrino mass density, we combine our data with WMAP7+BAO+Hubble constant measurements to constrain \Sigma m_\nu < 0.29 eV (95% C. L.).

We present the first dynamical mass estimates and scaling relations for a sample of Sunyaev-Zel'dovich effect (SZE) selected galaxy clusters. The sample consists of 16 massive clusters detected with the Atacama Cosmology Telescope (ACT) over a 455 sq. deg. area of the southern sky. Deep multi-object spectroscopic observations were taken to secure intermediate-resolution (R~700-800) spectra and redshifts for ~60 member galaxies on average per cluster. The dynamical masses M_200c of the clusters have been calculated using simulation-based scaling relations between velocity dispersion and mass. The sample has a median redshift z=0.50 and a median mass M_200c~12e14 Msun/h70 with a lower limit M_200c~6e14 Msun/h70, consistent with the expectations for the ACT southern sky survey. These masses are compared to the ACT SZE properties of the sample, specifically, the match-filtered central SZE amplitude y, the central Compton parameter y0, and the integrated Compton signal Y_200c, which we use to derive SZE-Mass scaling relations. All SZE estimators correlate with dynamical mass with low intrinsic scatter (<~20%), in agreement with numerical simulations. We explore the effects of various systematic effects on these scaling relations, including the correlation between observables and the influence of dynamically disturbed clusters. Using the 3-dimensional information available, we divide the sample into relaxed and disturbed clusters and find that ~50% of the clusters are disturbed. There are hints that disturbed systems might bias the scaling relations but given the current sample sizes these differences are not significant; further studies including more clusters are required to assess the impact of these clusters on the scaling relations.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

We cross-correlate the 100 um Improved Reprocessing of the IRAS Survey (IRIS) map and galaxy clusters at 0.1 < z < 0.3 in the maxBCG catalogue taken from the Sloan Digital Sky Survey, measuring an angular cross-power spectrum over multipole moments 150 < l < 3000 at a total significance of over 40 sigma. The cross-spectrum, which arises from the spatial correlation between unresolved dusty galaxies that make up the cosmic infrared background (CIB) in the IRIS map and the galaxy clusters, is well-fit by a single power law with an index of -1.28 +/- 0.12, similar to the clustering of unresolved galaxies from cross-correlating far-infrared and submillimetre maps at longer wavelengths. Using a recent, phenomenological model for the spectral and clustering properties of the IRIS galaxies, we constrain the large-scale bias of the maxBCG clusters to be 2.6 +/- 1.4, consistent with existing analyses of the real-space cluster correlation function. The success of our method suggests that future CIB-optical cross-correlations using Planck and Herschel data will significantly improve our understanding of the clustering and redshift distribution of the faint CIB sources.

We present the measured Sunyaev-Zel'dovich (SZ) flux from 474 optically-selected MaxBCG clusters that fall within the Atacama Cosmology Telescope (ACT) Equatorial survey region. The ACT Equatorial region used in this analysis covers 510 square degrees and overlaps Stripe 82 of the Sloan Digital Sky Survey. We also present the measured SZ flux stacked on 52 X-ray-selected MCXC clusters that fall within the ACT Equatorial region and an ACT Southern survey region covering 455 square degrees. We find that the measured SZ flux from the X-ray-selected clusters is consistent with expectations. However, we find that the measured SZ flux from the optically-selected clusters is both significantly lower than expectations and lower than the recovered SZ flux measured by the Planck satellite. Since we find a lower recovered SZ signal than Planck, we investigate the possibility that there is a significant offset between the optically-selected brightest cluster galaxies (BCGs) and the SZ centers, to which ACT is more sensitive due to its finer resolution. Such offsets can arise due to either an intrinsic physical separation between the BCG and the center of the gas concentration or from misidentification of the cluster BCG. We find that the entire discrepancy for both ACT and Planck can be explained by assuming that the BCGs are offset from the SZ maxima with a uniform random distribution between 0 and 1.5 Mpc. Such large offsets between gas peaks and BCGs for optically-selected cluster samples seem unlikely given that we find the physical separation between BCGs and X-ray peaks for an X-ray-selected subsample of MaxBCG clusters to have a much narrower distribution that peaks within 0.2 Mpc. It is possible that other effects are lowering the ACT and Planck signals by the same amount, with offsets between BCGs and SZ peaks explaining the remaining difference between measurements. (Abridged)

We present the temperature power spectra of the cosmic microwave background (CMB) derived from the three seasons of data from the Atacama Cosmology Telescope (ACT) at 148 GHz and 218 GHz, as well as the cross-frequency spectrum between the two channels. We detect and correct for contamination due to the Galactic cirrus in our equatorial maps. We present the results of a number of tests for possible systematic error and conclude that any effects are not significant compared to the statistical errors we quote. Where they overlap, we cross-correlate the ACT and the South Pole Telescope (SPT) maps and show they are consistent. The measurements of higher-order peaks in the CMB power spectrum provide an additional test of the Lambda CDM cosmological model, and help constrain extensions beyond the standard model. The small angular scale power spectrum also provides constraining power on the Sunyaev-Zel'dovich effects and extragalactic foregrounds. We also present a measurement of the CMB gravitational lensing convergence power spectrum at 4.6-sigma detection significance.

We present a detection of the unnormalized skewness <T^3> induced by the thermal Sunyaev-Zel'dovich (tSZ) effect in filtered Atacama Cosmology Telescope (ACT) 148 GHz cosmic microwave background temperature maps. Contamination due to infrared and radio sources is minimized by template subtraction of resolved sources and by constructing a mask using outlying values in the 218 GHz (tSZ-null) ACT maps. We measure <T^3>= -31 +- 6 \mu K^3 (measurement error only) or +- 14 \mu K^3 (including cosmic variance error) in the filtered ACT data, a 5-sigma detection. We show that the skewness is a sensitive probe of sigma_8, and use analytic calculations and tSZ simulations to obtain cosmological constraints from this measurement. From this signal alone we infer a value of sigma_8= 0.79 +0.03 -0.03 (68 % C.L.) +0.06 -0.06 (95 % C.L.). Our results demonstrate that measurements of non-Gaussianity can be a useful method for characterizing the tSZ effect and extracting the underlying cosmological information.

We measure the cross-correlation of Atacama Cosmology Telescope CMB lensing convergence maps with quasar maps made from the Sloan Digital Sky Survey DR8 SDSS-XDQSO photometric catalog. The CMB lensing-quasar cross-power spectrum is detected for the first time at a significance of 3.8 sigma, which directly confirms that the quasar distribution traces the mass distribution at high redshifts z>1. Our detection passes a number of null tests and systematic checks. Using this cross-power spectrum, we measure the amplitude of the linear quasar bias assuming a template for its redshift dependence, and find the amplitude to be consistent with an earlier measurement from clustering; at redshift z ~ 1.4, the peak of the distribution of quasars in our maps, our measurement corresponds to a bias of b = 2.5 +/- 0.6. With the signal-to-noise ratio on CMB lensing measurements likely to improve by an order of magnitude over the next few years, our results demonstrate the potential of CMB lensing cross-correlations to probe astrophysics at high redshifts.

We show how the correlation between the thermal Sunyaev Zel'dovich effect (tSZ) from galaxy clusters and dust emission from cosmic infrared background (CIB) sources can be calculated in a halo model framework. Using recent tSZ and CIB models, we find that the size of the tSZ x CIB cross-correlation is approximately 20 per cent at 150 GHz. The contribution to the total angular power spectrum is of order -2 \mu K^2 at ell=3000, however, this value is uncertain by a factor of two to three, primarily because of CIB source modelling uncertainties. We expect the large uncertainty in this component to degrade upper limits on the kinematic Sunyaev Zel'dovich effect (kSZ), due to similarity in the frequency dependence of the tSZ x CIB and kSZ across the frequency range probed by current Cosmic Microwave Background missions. We also find that the degree of tSZ x CIB correlation is higher for mm x sub-mm spectra than mm x mm, because more of the sub-mm CIB originates at lower redshifts (z<2), where most tSZ clusters are found.

We present a description of the data reduction and mapmaking pipeline used for the 2008 observing season of the Atacama Cosmology Telescope (ACT). The data presented here at 148 GHz represent 12% of the 90 TB collected by ACT from 2007 to 2010. In 2008 we observed for 136 days, producing a total of 1423 hours of data (11 TB for the 148 GHz band only), with a daily average of 10.5 hours of observation. From these, 1085 hours were devoted to a 850 deg^2 stripe (11.2 hours by 9.1 deg) centered on a declination of -52.7 deg, while 175 hours were devoted to a 280 deg^2 stripe (4.5 hours by 4.8 deg) centered at the celestial equator. We discuss sources of statistical and systematic noise, calibration, telescope pointing, and data selection. Out of 1260 survey hours and 1024 detectors per array, 816 hours and 593 effective detectors remain after data selection for this frequency band, yielding a 38% survey efficiency. The total sensitivity in 2008, determined from the noise level between 5 Hz and 20 Hz in the time-ordered data stream (TOD), is 32 micro-Kelvin sqrt{s} in CMB units. Atmospheric brightness fluctuations constitute the main contaminant in the data and dominate the detector noise covariance at low frequencies in the TOD. The maps were made by solving the least-squares problem using the Preconditioned Conjugate Gradient method, incorporating the details of the detector and noise correlations. Cross-correlation with WMAP sky maps, as well as analysis from simulations, reveal that our maps are unbiased at multipoles ell > 300. This paper accompanies the public release of the 148 GHz southern stripe maps from 2008. The techniques described here will be applied to future maps and data releases.

Using high-resolution microwave sky maps made by the Atacama Cosmology Telescope, we for the first time present strong evidence for motions of galaxy clusters and groups via microwave background temperature distortions due to the kinematic Sunyaev-Zel'dovich effect. Galaxy clusters are identified by their constituent luminous galaxies observed by the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. We measure the mean pairwise momentum of clusters, with a probability of the signal being due to random errors of 0.002, and the signal is consistent with the growth of cosmic structure in the standard model of cosmology.

We perform a combined fit to angular power spectra of unresolved infrared (IR) point sources from the Planck satellite (at 217, 353, 545 and 857 GHz, over angular scales 100 < l < 2200), the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST; 250, 350 and 500 um; 1000 < l < 9000), and from correlating BLAST and Atacama Cosmology Telescope (ACT; 148 and 218 GHz) maps. We find that the clustered power over the range of angular scales and frequencies considered is well fit by a simple power law of the form C_l \propto l^-n with n = 1.25 +/- 0.06. While the IR sources are understood to lie at a range of redshifts, with a variety of dust properties, we find that the frequency dependence of the clustering power can be described by the square of a modified blackbody, nu^beta B(nu,T_eff), with a single emissivity index beta = 2.20 +/- 0.07 and effective temperature T_eff = 9.7 K. Our predictions for the clustering amplitude are consistent with existing ACT and South Pole Telescope results at around 150 and 220 GHz, as is our prediction for the effective dust spectral index, which we find to be alpha_150-220 = 3.68 +/- 0.07 between 150 and 220 GHz. Our constraints on the clustering shape and frequency dependence can be used to model the IR clustering as a contaminant in Cosmic Microwave Background anisotropy measurements. The combined Planck and BLAST data also rule out a linear bias clustering model.

We present constraints on the primordial power spectrum of adiabatic fluctuations using data from the 2008 Southern Survey of the Atacama Cosmology Telescope (ACT). The angular resolution of ACT provides sensitivity to scales beyond \ell = 1000 for resolution of multiple peaks in the primordial temperature power spectrum, which enables us to probe the primordial power spectrum of adiabatic scalar perturbations with wavenumbers up to k \simeq 0.2 Mpc^{-1}. We find no evidence for deviation from power-law fluctuations over two decades in scale. Matter fluctuations inferred from the primordial temperature power spectrum evolve over cosmic time and can be used to predict the matter power spectrum at late times; we illustrate the overlap of the matter power inferred from CMB measurements (which probe the power spectrum in the linear regime) with existing probes of galaxy clustering, cluster abundances and weak lensing constraints on the primordial power. This highlights the range of scales probed by current measurements of the matter power spectrum.

We present measurements of the auto- and cross-frequency correlation power spectra of the cosmic (sub)millimeter background at: 250, 350, and 500 um (1200, 860, and 600 GHz) from observations made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST; and at 1380 and 2030 um (218 and 148 GHz) from observations made with the Atacama Cosmology Telescope, ACT. The overlapping observations cover 8.6 deg^2 in an area relatively free of Galactic dust near the south ecliptic pole (SEP). The ACT bands are sensitive to radiation from the CMB, the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio and dusty star-forming galaxies (DSFGs), while the dominant contribution to the BLAST bands is from DSFGs. We confirm and extend the BLAST analysis of clustering with an independent pipeline, and also detect correlations between the ACT and BLAST maps at over 25sigma significance, which we interpret as a detection of the DSFGs in the ACT maps. In addition to a Poisson component in the cross-frequency power spectra, we detect a clustered signal at >4sigma, and using a model for the DSFG evolution and number counts, we successfully fit all our spectra with a linear clustering model and a bias that depends only on redshift and not on scale. Finally, the data are compared to, and generally agree with, phenomenological models for the DSFG population. This study represents a first of its kind, and demonstrates the constraining power of the cross-frequency correlation technique to constrain models for the DSFGs. Similar analyses with more data will impose tight constraints on future models.