We present new constraints on anisotropic birefringence of the cosmic microwave background polarization using two seasons of data from the Atacama Cosmology Telescope covering $456$ square degrees of sky. The birefringence power spectrum, measured using a curved-sky quadratic estimator, is consistent with zero. Our results provide the tightest current constraint on birefringence over a range of angular scales between $5$ arcminutes and $9$ degrees. We improve previous upper limits on the amplitude of a scale-invariant birefringence power spectrum by a factor of between $2$ and $3$. Assuming a nearly-massless axion field during inflation, our result is equivalent to a $2\,\sigma$ upper limit on the Chern-Simons coupling constant between axions and photons of $g_{\alpha\gamma}<4.0\times 10^{-2}/H_I$ where $H_I$ is the inflationary Hubble scale.

We construct cosmic microwave background lensing mass maps using data from the 2014 and 2015 seasons of observations with the Atacama Cosmology Telescope (ACT). These maps cover 2100 square degrees of sky and overlap with a wide variety of optical surveys. The maps are signal dominated on large scales and have fidelity such that their correlation with the cosmic infrared background is clearly visible by eye. We also create lensing maps with thermal Sunyaev-Zel'dovich contamination removed using a novel cleaning procedure that only slightly degrades the lensing signal-to-noise ratio. The cross-spectrum between the cleaned lensing map and the BOSS CMASS galaxy sample is detected at $10$-$\sigma$ significance, with an amplitude of $A=1.02 \pm 0.10$ relative to the Planck best-fit LCDM cosmological model with fiducial linear galaxy bias. Our measurement lays the foundation for lensing cross-correlation science with current ACT data and beyond.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

We examine the stacked thermal Sunyaev-Zel\text{'}dovich (SZ) signals for a sample of galaxy cluster candidates from the Spitzer-HETDEX Exploratory Large Area (SHELA) Survey, which are identified in combined optical and infrared SHELA data using the redMaPPer algorithm. We separate the clusters into three richness bins, with average photometric redshifts ranging from 0.70 to 0.80. The richest bin shows a clear temperature decrement at 148 GHz in the Atacama Cosmology Telescope data, which we attribute to the SZ effect. All richness bins show an increment at 220 GHz, which we attribute to dust emission from cluster galaxies. We correct for dust emission using stacked profiles from Herschel Stripe 82 data, and allow for synchrotron emission using stacked profiles created by binning source fluxes from NVSS data. We see dust emission in all three richness bins, but can only confidently detect the SZ decrement in the highest richness bin, finding $M_{500}$ = $8.7^{+1.7}_{-1.3} \times 10^{13} M_\odot$. Neglecting the correction for dust depresses the inferred mass by 26 percent, indicating a partial fill-in of the SZ decrement from thermal dust and synchrotron emission by the cluster member galaxies. We compare our corrected SZ masses to two redMaPPer mass--richness scaling relations and find that the SZ mass is lower than predicted by the richness. We discuss possible explanations for this discrepancy, and note that the SHELA richnesses may differ from previous richness measurements due to the inclusion of IR data in redMaPPer.

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 provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

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

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.

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

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

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

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

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

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

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

We use a complete sample of 38 richest maxBCG clusters to study the ICM-galaxy scaling relations and the halo mass selection properties of the maxBCG algorithm, based on X-ray and optical observations. The clusters are selected from the two largest bins of optical richness in the Planck stacking work with the maxBCG richness $N_{200} \geq 78$. We analyze their Chandra and XMM-Newton data to derive the X-ray properties of the ICM. We then use the distribution of $P(X|N)$, $X=T_X,\ L_X,\ Y_X$, to study the mass selection $P(M|N)$ of maxBCG. Compared with previous works based on the whole richness sample, a significant fraction of blended systems with boosted richness is skewed into this richest sample. Parts of the blended halos are picked apart by the redMaPPer, an updated red-sequence cluster finding algorithm with lower mass scatter. Moreover, all the optical blended halos are resolved as individual X-ray halos, following the established $L_X-T_X$ and $L_X-Y_X$ relations. We further discuss that the discrepancy between ICM-galaxy scaling relations, especially for future blind stacking, can come from several factors, including miscentering, projection, contamination of low mass systems, mass bias and covariance bias. We also evaluate the fractions of relaxed and cool core clusters in our sample. Both are smaller than those from SZ or X-ray selected samples. Moreover, disturbed clusters show a higher level of mass bias than relaxed clusters.

We present a method to measure the small-scale matter power spectrum using high-resolution measurements of the gravitational lensing of the Cosmic Microwave Background (CMB). To determine whether small-scale structure today is suppressed on scales below 10 kiloparsecs (corresponding to M < 10^9 M_sun), one needs to probe CMB-lensing modes out to L ~ 35,000, requiring a CMB experiment with about 20 arcsecond resolution or better. We show that a CMB survey covering 4,000 square degrees of sky, with an instrumental sensitivity of 0.5 uK-arcmin at 18 arcsecond resolution, could distinguish between cold dark matter and an alternative, such as 1 keV warm dark matter or 10^(-22) eV fuzzy dark matter with about 4-sigma significance. A survey of the same resolution with 0.1 uK-arcmin noise could distinguish between cold dark matter and these alternatives at better than 20-sigma significance; such high-significance measurements may also allow one to distinguish between a suppression of power due to either baryonic effects or the particle nature of dark matter, since each impacts the shape of the lensing power spectrum differently. CMB temperature maps yield higher signal-to-noise than polarization maps in this small-scale regime; thus, systematic effects, such as from extragalactic astrophysical foregrounds, need to be carefully considered. However, these systematic concerns can likely be mitigated with known techniques. Next-generation CMB lensing may thus provide a robust and powerful method of measuring the small-scale matter power spectrum.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

We present a method to delens the acoustic peaks of the CMB temperature and polarization power spectra internally, using lensing maps reconstructed from the CMB itself. We find that when delensing CMB acoustic peaks with a lensing potential map derived from the same CMB sky, a large bias arises in the delensed power spectrum. The cause of this bias is that the noise in the reconstructed potential map is derived from, and hence correlated with, the CMB map when delensing. This bias is more significant relative to the signal than an analogous bias found when delensing CMB B modes. We calculate the leading term of this bias, which is present even in the absence of lensing. We also demonstrate one method to remove this bias, using reconstructions from CMB angular scales within given ranges to delens CMB scales outside of those ranges. Some details relevant for a realistic analysis are also discussed, such as the importance of removing mask-induced effects for successful delensing, and a useful null test, obtained from randomizing the phases of the reconstructed potential. Our findings should help current and next-generation CMB experiments obtain tighter parameter constraints via the internal removal of lensing-induced smoothing from temperature and E-mode acoustic peaks.

We present a new measurement of the kinematic Sunyaev-Zeldovich effect using data from the Atacama Cosmology Telescope (ACT) and the Baryon Oscillation Spectroscopic Survey (BOSS). Using 600 square degrees of overlapping sky area, we evaluate the mean pairwise baryon momentum associated with the positions of 50,000 bright galaxies in the BOSS DR11 Large Scale Structure catalog. A non-zero signal arises from the large-scale motions of halos containing the sample galaxies. The data fits an analytical signal model well, with the optical depth to microwave photon scattering as a free parameter determining the overall signal amplitude. We estimate the covariance matrix of the mean pairwise momentum as a function of galaxy separation, using microwave sky simulations, jackknife evaluation, and bootstrap estimates. The most conservative simulation-based errors give signal-to-noise estimates between 3.6 and 4.1 for varying galaxy luminosity cuts. We discuss how the other error determinations can lead to higher signal-to-noise values, and consider the impact of several possible systematic errors. Estimates of the optical depth from the average thermal Sunyaev-Zeldovich signal at the sample galaxy positions are broadly consistent with those obtained from the mean pairwise momentum signal.

We measure the gravitational lensing shear signal around dark matter halos hosting CMASS galaxies using light sources at $z\sim 1$ (background galaxies) and at the surface of last scattering at $z\sim 1100$ (the cosmic microwave background). The galaxy shear measurement uses data from the CFHTLenS survey, and the microwave background shear measurement uses data from the {\it Planck} satellite. The ratio of shears from these cross-correlations provides a purely geometric distance measurement across the longest possible cosmological lever arm. This is because the matter distribution around the halos, including uncertainties in galaxy bias and systematic errors such as miscentering, cancels in the ratio for halos in thin redshift slices. We measure this distance ratio in three different redshift slices of the CMASS sample, and combine them to obtain a $17\%$ measurement of the distance ratio, $r=0.390^{+0.070}_{-0.062}$ at an effective redshift of $z=0.53$. This is consistent with the predicted ratio from the {\it Planck} best-fit $\Lambda$CDM cosmology of $r=0.419$.

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.

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

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

Mass calibration uncertainty is the largest systematic effect for using clusters of galaxies to constrain cosmological parameters. We present weak lensing mass measurements from the Canada-France-Hawaii Telescope Stripe 82 Survey for galaxy clusters selected through their high signal-to-noise thermal Sunyaev-Zeldovich (tSZ) signal measured with the Atacama Cosmology Telescope (ACT). For a sample of 9 ACT clusters with a tSZ signal-to-noise greater than five the average weak lensing mass is $\left(4.8\pm0.8\right)\,\times10^{14}\,\mathrm{M}_\odot$, consistent with the tSZ mass estimate of $\left(4.70\pm1.0\right)\,\times10^{14}\,\mathrm{M}_\odot$ which assumes a universal pressure profile for the cluster gas. Our results are consistent with previous weak-lensing measurements of tSZ-detected clusters from the Planck satellite. When comparing our results, we estimate the Eddington bias correction for the sample intersection of Planck and weak-lensing clusters which was previously excluded.

In recent years there have been significant improvements in the sensitivity and the angular resolution of the instruments dedicated to the observation of the Cosmic Microwave Background (CMB). ACTPol is the first polarization receiver for the Atacama Cosmology Telescope (ACT) and is observing the CMB sky with arcmin resolution over about 2000 sq. deg. Its upgrade, Advanced ACTPol (AdvACT), will observe the CMB in five frequency bands and over a larger area of the sky. We describe the optimization and implementation of the ACTPol and AdvACT surveys. The selection of the observed fields is driven mainly by the science goals, that is, small angular scale CMB measurements, B-mode measurements and cross-correlation studies. For the ACTPol survey we have observed patches of the southern galactic sky with low galactic foreground emissions which were also chosen to maximize the overlap with several galaxy surveys to allow unique cross-correlation studies. A wider field in the northern galactic cap ensured significant additional overlap with the BOSS spectroscopic survey. The exact shapes and footprints of the fields were optimized to achieve uniform coverage and to obtain cross-linked maps by observing the fields with different scan directions. We have maximized the efficiency of the survey by implementing a close to 24 hour observing strategy, switching between daytime and nighttime observing plans and minimizing the telescope idle time. We describe the challenges represented by the survey optimization for the significantly wider area observed by AdvACT, which will observe roughly half of the low-foreground sky. The survey strategies described here may prove useful for planning future ground-based CMB surveys, such as the Simons Observatory and CMB Stage IV surveys.

The Atacama Cosmology Telescope (ACT) is designed to make high angular resolution measurements of anisotropies in the Cosmic Microwave Background (CMB) at millimeter wavelengths. We describe ACTPol, an upgraded receiver for ACT, which uses feedhorn-coupled, polarization-sensitive detector arrays, a 3 degree field of view, 100 mK cryogenics with continuous cooling, and meta material anti-reflection coatings. ACTPol comprises three arrays with separate cryogenic optics: two arrays at a central frequency of 148 GHz and one array operating simultaneously at both 97 GHz and 148 GHz. The combined instrument sensitivity, angular resolution, and sky coverage are optimized for measuring angular power spectra, clusters via the thermal Sunyaev-Zel'dovich and kinetic Sunyaev-Zel'dovich signals, and CMB lensing due to large scale structure. The receiver was commissioned with its first 148 GHz array in 2013, observed with both 148 GHz arrays in 2014, and has recently completed its first full season of operations with the full suite of three arrays. This paper provides an overview of the design and initial performance of the receiver and related systems.

We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope (ACTPol) at 146 GHz, together with the Constant Mass CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Ze\v{l}dovich (kSZ) effect over the redshift range z = 0.4 - 0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the cosmic microwave background temperature at the location of each halo, weighted by the corresponding reconstructed velocity. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3sigma and 2.9sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Ze\v{l}dovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.

Advanced ACTPol is a polarization-sensitive upgrade for the 6 m aperture Atacama Cosmology Telescope (ACT), adding new frequencies and increasing sensitivity over the previous ACTPol receiver. In 2016, Advanced ACTPol will begin to map approximately half the sky in five frequency bands (28-230 GHz). Its maps of primary and secondary cosmic microwave background (CMB) anisotropies -- imaged in intensity and polarization at few arcminute-scale resolution -- will enable precision cosmological constraints and also a wide array of cross-correlation science that probes the expansion history of the universe and the growth of structure via gravitational collapse. To accomplish these scientific goals, the Advanced ACTPol receiver will be a significant upgrade to the ACTPol receiver, including four new multichroic arrays of cryogenic, feedhorn-coupled AlMn transition edge sensor (TES) polarimeters (fabricated on 150 mm diameter wafers); a system of continuously rotating meta-material silicon half-wave plates; and a new multiplexing readout architecture which uses superconducting quantum interference devices (SQUIDs) and time division to achieve a 64-row multiplexing factor. Here we present the status and scientific goals of the Advanced ACTPol instrument, emphasizing the design and implementation of the Advanced ACTPol cryogenic detector arrays.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

We present evidence of the gravitational lensing of the cosmic microwave background by $10^{13}$ solar mass dark matter halos. Lensing convergence maps from the Atacama Cosmology Telescope Polarimeter (ACTPol) are stacked at the positions of around 12,000 optically-selected CMASS galaxies from the SDSS-III/BOSS survey. The mean lensing signal is consistent with simulated dark matter halo profiles, and is favored over a null signal at 3.2 sigma significance. This result demonstrates the potential of microwave background lensing to probe the dark matter distribution in galaxy group and galaxy cluster halos.

We present a multi-wavelength analysis of eleven Sunyaev Zel'dovich effect (SZE)-selected galaxy clusters (ten with new data) from the Atacama Cosmology Telescope (ACT) southern survey. We have obtained new imaging from the Large APEX Bolometer Camera (345GHz; LABOCA) on the Atacama Pathfinder EXperiment (APEX) telescope, the Australia Telescope Compact Array (2.1GHz; ATCA), and the Spectral and Photometric Imaging Receiver (250, 350, and $500\,\rm\mu m$; SPIRE) on the Herschel Space Observatory. Spatially-resolved 345GHz SZE increments with integrated S/N > 5 are found in six clusters. We compute 2.1GHz number counts as a function of cluster-centric radius and find significant enhancements in the counts of bright sources at projected radii $\theta < \theta_{2500}$. By extrapolating in frequency, we predict that the combined signals from 2.1GHz-selected radio sources and 345GHz-selected SMGs contaminate the 148GHz SZE decrement signal by ~5% and the 345GHz SZE increment by ~18%. After removing radio source and SMG emission from the SZE signals, we use ACT, LABOCA, and (in some cases) new Herschel SPIRE imaging to place constraints on the clusters' peculiar velocities. The sample's average peculiar velocity relative to the cosmic microwave background is $153\pm 383\,\rm km\,s^{-1}$.

We correlate the positions of radio galaxies in the FIRST survey with the CMB lensing convergence estimated from the Atacama Cosmology Telescope over 470 square degrees to determine the bias of these galaxies. We remove optically cross-matched sources below redshift $z=0.2$ to preferentially select Active Galactic Nuclei (AGN). We measure the angular cross-power spectrum $C_l^{\kappa g}$ at $4.4\sigma$ significance in the multipole range $100<l<3000$, corresponding to physical scales between $\approx$ 2--60 Mpc at an effective redshift $z_{\rm eff}= 1.5$. Modelling the AGN population with a redshift-dependent bias, the cross-spectrum is well fit by the Planck best-fit $\Lambda$CDM cosmological model. Fixing the cosmology we fit for the overall bias model normalization, finding $b(z_{\rm eff}) = 3.5 \pm 0.8$ for the full galaxy sample, and $b(z_{\rm eff})=4.0\pm1.1 (3.0\pm1.1)$ for sources brighter (fainter) than 2.5 mJy. This measurement characterizes the typical halo mass of radio-loud AGN: we find $\log(M_{\rm halo} / M_\odot) = 13.6^{+0.3}_{-0.4}$.

We develop a general framework for generating estimators of a given quantity which are unbiased to a given order in the difference between the true value of the underlying quantity and the fiducial position in theory space around which we expand the likelihood. We apply this formalism to rederive the optimal quadratic estimator and show how the replacement of the second derivative matrix with the Fisher matrix is a generic way of creating an unbiased estimator (assuming choice of the fiducial model is independent of data). Next apply the approach to estimation of shear lensing, closely following the work of Bernstein and Armstrong (2014). Our first order estimator reduces to their estimator in the limit of zero shear, but it also naturally allows for the case of non-constant shear and the easy calculation of correlation functions or power spectra using standard methods. Both our first-order estimator and Bernstein and Armstrong's estimator exhibit a bias which is quadratic in true shear. Our third-order estimator is, at least in the realm of the toy problem of Bernstein and Armstrong, unbiased to 0.1% in relative shear errors $\Delta g/|g|$ for shears up to $|g|=0.2$.

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.

Updated constraints on dark matter cross section and mass are presented combining CMB power spectrum measurements from Planck, WMAP9, ACT, and SPT as well as several low-redshift datasets (BAO, HST, supernovae). For the CMB datasets, we combine WMAP9 temperature and polarization data for l <= 431 with Planck temperature data for 432 < l < 2500, ACT and SPT data for l > 2500, and Planck CMB four-point lensing measurements. We allow for redshift-dependent energy deposition from dark matter annihilation by using a `universal' energy absorption curve. We also include an updated treatment of the excitation, heating, and ionization energy fractions, and provide updated deposition efficiency factors (f_eff) for 41 different dark matter models. Assuming perfect energy deposition (f_eff = 1) and a thermal cross section, dark matter masses below 26 GeV are excluded at the 2-sigma level. Assuming a more generic efficiency of f_eff = 0.2, thermal dark matter masses below 5 GeV are disfavored at the 2-sigma level. These limits are a factor of ~2 improvement over those from WMAP9 data alone. These current constraints probe, but do not exclude, dark matter as an explanation for reported anomalous indirect detection observations from AMS-02/PAMELA and the Fermi Gamma-ray Inner Galaxy data. They also probe relevant models that would explain anomalous direct detection events from CDMS, CRESST, CoGeNT, and DAMA, as originating from a generic thermal WIMP. Projected constraints from the full Planck release should improve the current limits by another factor of ~2, but will not definitely probe these signals. The proposed CMB Stage IV experiment will more decisively explore the relevant regions and improve upon the Planck constraints by another factor of ~2.

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.

Future arcminute resolution polarization data from ground-based Cosmic Microwave Background (CMB) observations can be used to estimate the contribution to the temperature power spectrum from the primary anisotropies and to uncover the signature of reionization near $\ell=1500$ in the small angular-scale temperature measurements. Our projections are based on combining expected small-scale E-mode polarization measurements from Advanced ACTPol in the range $300<\ell<3000$ with simulated temperature data from the full Planck mission in the low and intermediate $\ell$ region, $2<\ell<2000$. We show that the six basic cosmological parameters determined from this combination of data will predict the underlying primordial temperature spectrum at high multipoles to better than $1\%$ accuracy. Assuming an efficient cleaning from multi-frequency channels of most foregrounds in the temperature data, we investigate the sensitivity to the only residual secondary component, the kinematic Sunyaev-Zel'dovich (kSZ) term. The CMB polarization is used to break degeneracies between primordial and secondary terms present in temperature and, in effect, to remove from the temperature data all but the residual kSZ term. We estimate a $15 \sigma$ detection of the diffuse homogeneous kSZ signal from expected AdvACT temperature data at $\ell>1500$, leading to a measurement of the amplitude of matter density fluctuations, $\sigma_8$, at $1\%$ precision. Alternatively, by exploring the reionization signal encoded in the patchy kSZ measurements, we bound the time and duration of the reionization with $\sigma(z_{\rm re})=1.1$ and $\sigma(\Delta z_{\rm re})=0.2$. We find that these constraints degrade rapidly with large beam sizes, which highlights the importance of arcminute-scale resolution for future CMB surveys.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

This is a report on the status and prospects of the quantification of neutrino properties through the cosmological neutrino background for the Cosmic Frontier of the Division of Particles and Fields Community Summer Study long-term planning exercise. Experiments planned and underway are prepared to study the cosmological neutrino background in detail via its influence on distance-redshift relations and the growth of structure. The program for the next decade described in this document, including upcoming spectroscopic galaxy surveys eBOSS and DESI and a new Stage-IV CMB polarization experiment CMB-S4, will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass measurement will produce a high significance detection of non-zero sum m_nu, whose lower bound derived from atmospheric and solar neutrino oscillation data is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this measurement will definitively rule out the inverted neutrino mass hierarchy, shedding light on one of the most puzzling aspects of the Standard Model of particle physics --- the origin of mass. This precise a measurement of N_eff will allow for high sensitivity to any light and dark degrees of freedom produced in the big bang and a precision test of the standard cosmological model prediction that N_eff = 3.046.

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.

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

The lensing power spectrum from cosmic microwave background (CMB) temperature maps will be measured with unprecedented precision with upcoming experiments, including upgrades to ACT and SPT. Achieving significant improvements in cosmological parameter constraints, such as percent level errors on sigma_8 and an uncertainty on the total neutrino mass of approximately 50 meV, requires percent level measurements of the CMB lensing power. This necessitates tight control of systematic biases. We study several types of biases to the temperature-based lensing reconstruction signal from foreground sources such as radio and infrared galaxies and the thermal Sunyaev-Zel'dovich effect from galaxy clusters. These foregrounds bias the CMB lensing signal due to their non-Gaussian nature. Using simulations as well as some analytical models we find that these sources can substantially impact the measured signal if left untreated. However, these biases can be brought to the percent level if one masks galaxies with fluxes at 150 GHz above 1 mJy and galaxy clusters with masses above M_vir = 10^14 M_sun. To achieve such percent level bias, we find that only modes up to a maximum multipole of l_max ~ 2500 should be included in the lensing reconstruction. We also discuss ways to minimize additional bias induced by such aggressive foreground masking by, for example, exploring a two-step masking and in-painting algorithm.

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.

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

[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 describe the measurement of the beam profiles and window functions for the Atacama Cosmology Telescope (ACT), which operated from 2007 to 2010 with kilo-pixel bolometer arrays centered at 148, 218, and 277 GHz. Maps of Saturn are used to measure the beam shape in each array and for each season of observations. Radial profiles are transformed to Fourier space in a way that preserves the spatial correlations in the beam uncertainty, to derive window functions relevant for angular power spectrum analysis. Several corrections are applied to the resulting beam transforms, including an empirical correction measured from the final CMB survey maps to account for the effects of mild pointing variation and alignment errors. Observations of Uranus made regularly throughout each observing season are used to measure the effects of atmospheric opacity and to monitor deviations in telescope focus over the season. Using the WMAP-based calibration of the ACT maps to the CMB blackbody, we obtain precise measurements of the brightness temperatures of the Uranus and Saturn disks at effective frequencies of 149 and 219 GHz. For Uranus we obtain thermodynamic brightness temperatures T_U^{149} = 106.7 \pm 2.2 K and T_U^{219} = 100.1 \pm 3.1 K. For Saturn, we model the effects of the ring opacity and emission using a simple model and obtain resulting (unobscured) disk temperatures of T_S^{149} = 137.3 \pm 3.2 K and T_S^{219} = 137.3 \pm 4.7 K.

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)

Recent data from the WMAP, ACT and SPT experiments provide precise measurements of the cosmic microwave background temperature power spectrum over a wide range of angular scales. The combination of these observations is well fit by the standard, spatially flat LCDM cosmological model, constraining six free parameters to within a few percent. The scalar spectral index, n_s = 0.9690 +/- 0.0089, is less than unity at the 3.6 sigma level, consistent with simple models of inflation. The damping tail of the power spectrum at high resolution, combined with the amplitude of gravitational lensing measured by ACT and SPT, constrains the effective number of relativistic species to be N_eff = 3.28 +/- 0.40, in agreement with the standard model's three species of light neutrinos.

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 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 present a detection of the Sunyaev-Zel'dovich (SZ) decrement associated with the Luminous Red Galaxy (LRG) sample of the Sloan Digital Sky Survey. The SZ data come from 148 GHz maps of the equatorial region made by the Atacama Cosmology Telescope (ACT). The LRG sample is divided by luminosity into four bins, and estimates for the central Sunyaev-Zel'dovich temperature decrement are calculated through a stacking process. We detect and account for a bias of the SZ signal due to weak radio sources. We use numerical simulations to relate the observed decrement to Y200 and clustering properties to relate the galaxy luminosity bins to mass. We also use a relation between brightest cluster galaxy luminosity and cluster mass based on stacked gravitational lensing measurements to estimate the characteristic halo masses. The masses are found to be around 1e14 M_sun.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

We present a detailed analysis from new multi-wavelength observations of the exceptional galaxy cluster ACT-CL J0102-4915 "El Gordo," likely the most massive, hottest, most X-ray luminous and brightest Sunyaev-Zeldovich (SZ) effect cluster known at z>0.6. The Atacama Cosmology Telescope collaboration discovered El Gordo as the most significant SZ decrement in a sky survey area of 755 deg^2. Our VLT/FORS2 spectra of 89 member galaxies yield a cluster redshift, z=0.870, and velocity dispersion, s=1321+/-106 km/s. Our Chandra observations reveal a hot and X-ray luminous system with an integrated temperature of Tx=14.5+/-1.0 keV and 0.5-2.0 keV band luminosity of Lx=(2.19+/-0.11)x10^45 h70^-2 erg/s. We obtain several statistically consistent cluster mass estimates; using mass scaling relations with velocity dispersion, X-ray Yx, and integrated SZ, we estimate a cluster mass of M200a=(2.16+/-0.32)x10^15 M_sun/h70. The Chandra and VLT/FORS2 optical data also reveal that El Gordo is undergoing a major merger between components with a mass ratio of approximately 2 to 1. The X-ray data show significant temperature variations from a low of 6.6+/-0.7 keV at the merging low-entropy, high-metallicity, cool core to a high of 22+/-6 keV. We also see a wake in the X-ray surface brightness caused by the passage of one cluster through the other. Archival radio data at 843 MHz reveal diffuse radio emission that, if associated with the cluster, indicates the presence of an intense double radio relic, hosted by the highest redshift cluster yet. El Gordo is possibly a high-redshift analog of the famous Bullet Cluster. Such a massive cluster at this redshift is rare, although consistent with the standard L-CDM cosmology in the lower part of its allowed mass range. Massive, high-redshift mergers like El Gordo are unlikely to be reproduced in the current generation of numerical N-body cosmological simulations.

We present follow-up observations with the Sunyaev-Zel'dovich Array (SZA) of optically-confirmed galaxy clusters found in the equatorial survey region of the Atacama Cosmology Telescope (ACT): ACT-CL J0022-0036, ACT-CL J2051+0057, and ACT-CL J2337+0016. ACT-CL J0022-0036 is a newly-discovered, massive (10^15 Msun), high-redshift (z=0.81) cluster revealed by ACT through the Sunyaev-Zel'dovich effect (SZE). Deep, targeted observations with the SZA allow us to probe a broader range of cluster spatial scales, better disentangle cluster decrements from radio point source emission, and derive more robust integrated SZE flux and mass estimates than we can with ACT data alone. For the two clusters we detect with the SZA we compute integrated SZE signal and derive masses from the SZA data only. ACT-CL J2337+0016, also known as Abell 2631, has archival Chandra data that allow an additional X-ray-based mass estimate. Optical richness is also used to estimate cluster masses and shows good agreement with the SZE and X-ray-based estimates. Based on the point sources detected by the SZA in these three cluster fields and an extrapolation to ACT's frequency, we estimate that point sources could be contaminating the SZE decrement at the <= 20% level for some fraction of clusters.

For the first time, measurements of the cosmic microwave background radiation (CMB) alone favor cosmologies with $w=-1$ dark energy over models without dark energy at a 3.2-sigma level. We demonstrate this by combining the CMB lensing deflection power spectrum from the Atacama Cosmology Telescope with temperature and polarization power spectra from the Wilkinson Microwave Anisotropy Probe. The lensing data break the geometric degeneracy of different cosmological models with similar CMB temperature power spectra. Our CMB-only measurement of the dark energy density $\Omega_\Lambda$ confirms other measurements from supernovae, galaxy clusters and baryon acoustic oscillations, and demonstrates the power of CMB lensing as a new cosmological tool.

We report the first detection of the gravitational lensing of the cosmic microwave background through a measurement of the four-point correlation function in the temperature maps made by the Atacama Cosmology Telescope. We verify our detection by calculating the levels of potential contaminants and performing a number of null tests. The resulting convergence power spectrum at 2-degree angular scales measures the amplitude of matter density fluctuations on comoving length scales of around 100 Mpc at redshifts around 0.5 to 3. The measured amplitude of the signal agrees with Lambda Cold Dark Matter cosmology predictions. Since the amplitude of the convergence power spectrum scales as the square of the amplitude of the density fluctuations, the 4-sigma detection of the lensing signal measures the amplitude of density fluctuations to 12%.

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.

We report on twenty-three clusters detected blindly as Sunyaev-Zel'dovich (SZ) decrements in a 148 GHz, 455 square-degree map of the southern sky made with data from the Atacama Cosmology Telescope 2008 observing season. All SZ detections announced in this work have confirmed optical counterparts. Ten of the clusters are new discoveries. One newly discovered cluster, ACT-CL J0102-4915, with a redshift of 0.75 (photometric), has an SZ decrement comparable to the most massive systems at lower redshifts. Simulations of the cluster recovery method reproduce the sample purity measured by optical follow-up. In particular, for clusters detected with a signal-to-noise ratio greater than six, simulations are consistent with optical follow-up that demonstrated this subsample is 100% pure. The simulations further imply that the total sample is 80% complete for clusters with mass in excess of 6x10^14 solar masses referenced to the cluster volume characterized by five hundred times the critical density. The Compton y -- X-ray luminosity mass comparison for the eleven best detected clusters visually agrees with both self-similar and non-adiabatic, simulation-derived scaling laws.

We report on extragalactic sources detected in a 455 square-degree map of the southern sky made with data at a frequency of 148 GHz from the Atacama Cosmology Telescope 2008 observing season. We provide a catalog of 157 sources with flux densities spanning two orders of magnitude: from 15 to 1500 mJy. Comparison to other catalogs shows that 98% of the ACT detections correspond to sources detected at lower radio frequencies. Three of the sources appear to be associated with the brightest cluster galaxies of low redshift X-ray selected galaxy clusters. Estimates of the radio to mm-wave spectral indices and differential counts of the sources further bolster the hypothesis that they are nearly all radio sources, and that their emission is not dominated by re-emission from warm dust. In a bright (>50 mJy) 148 GHz-selected sample with complete cross-identifications from the Australia Telescope 20 GHz survey, we observe an average steepening of the spectra between 5, 20, and 148 GHz with median spectral indices of $\alpha_{\rm 5-20} = -0.07 \pm 0.06$, $\alpha_{\rm 20-148} = -0.39 \pm0.04$, and $\alpha_{\rm 5-148} = -0.20 \pm 0.03$. When the measured spectral indices are taken into account, the 148 GHz differential source counts are consistent with previous measurements at 30 GHz in the context of a source count model dominated by radio sources. Extrapolating with an appropriately rescaled model for the radio source counts, the Poisson contribution to the spatial power spectrum from synchrotron-dominated sources with flux density less than 20 mJy is $C^{\rm Sync} = (2.8 \pm 0.3) \times 10^{-6} \micro\kelvin^2$.

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.

Recent measurements of the Sunyaev-Zel'dovich (SZ) angular power spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) demonstrate the importance of understanding baryon physics when using the SZ power spectrum to constrain cosmology. This is challenging since roughly half of the SZ power at l=3000 is from low-mass systems with 10^13 h^-1 M_sun < M_500 < 1.5x10^14 h^-1 M_sun, which are more difficult to study than systems of higher mass. We present a study of the thermal pressure content for a sample of local galaxy groups from Sun et al. (2009). The group Y_{sph, 500} - M_500 relation agrees with the one for clusters derived by Arnaud et al. (2010). The group median pressure profile also agrees with the universal pressure profile for clusters derived by Arnaud et al. (2010). With this in mind, we briefly discuss several ways to alleviate the tension between the measured low SZ power and the predictions from SZ templates.

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.

The Atacama Cosmology Telescope (ACT) is currently observing the cosmic microwave background with arcminute resolution at 148 GHz, 218 GHz, and 277 GHz. In this paper, we present ACT's first results. Data have been analyzed using a maximum-likelihood map-making method which uses B-splines to model and remove the atmospheric signal. It has been used to make high-precision beam maps from which we determine the experiment's window functions. This beam information directly impacts all subsequent analyses of the data. We also used the method to map a sample of galaxy clusters via the Sunyaev-Zel'dovich (SZ) effect, and show five clusters previously detected with X-ray or SZ observations. We provide integrated Compton-y measurements for each cluster. Of particular interest is our detection of the z = 0.44 component of A3128 and our current non-detection of the low-redshift part, providing strong evidence that the further cluster is more massive as suggested by X-ray measurements. This is a compelling example of the redshift-independent mass selection of the SZ effect.

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel'dovich-selected galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sample used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cosmology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives sigma_8 = 0.851 +/- 0.115 and w = -1.14 +/- 0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation obtained from numerical simulations and calibrated by X-ray observations, we find sigma_8 = 0.821 +/- 0.044 and w = -1.05 +/- 0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations plus type Ia supernoava which give sigma_8 = 0.802 +/- 0.038 and w = -0.98 +/- 0.053. A stacking analysis of the clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and the cosmological parameters derived from them are broadly consistent with current models.

We present optical and X-ray properties for the first confirmed galaxy cluster sample selected by the Sunyaev-Zel'dovich Effect from 148 GHz maps over 455 square degrees of sky made with the Atacama Cosmology Telescope. These maps, coupled with multi-band imaging on 4-meter-class optical telescopes, have yielded a sample of 23 galaxy clusters with redshifts between 0.118 and 1.066. Of these 23 clusters, 10 are newly discovered. The selection of this sample is approximately mass limited and essentially independent of redshift. We provide optical positions, images, redshifts and X-ray fluxes and luminosities for the full sample, and X-ray temperatures of an important subset. The mass limit of the full sample is around 8e14 Msun, with a number distribution that peaks around a redshift of 0.4. For the 10 highest significance SZE-selected cluster candidates, all of which are optically confirmed, the mass threshold is 1e15 Msun and the redshift range is 0.167 to 1.066. Archival observations from Chandra, XMM-Newton, and ROSAT provide X-ray luminosities and temperatures that are broadly consistent with this mass threshold. Our optical follow-up procedure also allowed us to assess the purity of the ACT cluster sample. Eighty (one hundred) percent of the 148 GHz candidates with signal-to-noise ratios greater than 5.1 (5.7) are confirmed as massive clusters. The reported sample represents one of the largest SZE-selected sample of massive clusters over all redshifts within a cosmologically-significant survey volume, which will enable cosmological studies as well as future studies on the evolution, morphology, and stellar populations in the most massive clusters in the Universe.

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg^2 with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500<l<10000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. The power from thermal and kinetic SZ at 148 GHz is estimated to be B_3000 = 6.8+-2.9 uK^2, where B_l=l(l+1)C_l/2pi. We estimate primary cosmological parameters from the 148 GHz spectrum, marginalizing over SZ and source power. The LCDM cosmological model is a good fit to the data, and LCDM parameters estimated from ACT+WMAP are consistent with the 7-year WMAP limits, with scale invariant n_s = 1 excluded at 99.7% CL (3sigma). A model with no CMB lensing is disfavored at 2.8sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6sigma detection of primordial helium, with Y_P = 0.313+-0.044, and a 4sigma detection of relativistic species, assumed to be neutrinos, with Neff = 5.3+-1.3 (4.6+-0.8 with BAO+H0 data). From the CMB alone the running of the spectral index is constrained to be dn/dlnk = -0.034 +- 0.018, the limit on the tensor-to-scalar ratio is r<0.25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension Gmu<1.6 \times 10^-7 (95% CL).

We present measurements of the cosmic microwave background (CMB) power spectrum made by the Atacama Cosmology Telescope at 148 GHz and 218 GHz, as well as the cross-frequency spectrum between the two channels. Our results clearly show the second through the seventh acoustic peaks in the CMB power spectrum. The measurements of these higher-order peaks provide an additional test of the {\Lambda}CDM cosmological model. At l > 3000, we detect power in excess of the primary anisotropy spectrum of the CMB. At lower multipoles 500 < l < 3000, we find evidence for gravitational lensing of the CMB in the power spectrum at the 2.8{\sigma} level. We also detect a low level of Galactic dust in our maps, which demonstrates that we can recover known faint, diffuse signals.

We present a new calibration method based on cross-correlations with WMAP and apply it to data from the Atacama Cosmology Telescope (ACT). ACT's observing strategy and map making procedure allows an unbiased reconstruction of the modes in the maps over a wide range of multipoles. By directly matching the ACT maps to WMAP observations in the multipole range of 400 < ell < 1000, we determine the absolute calibration with an uncertainty of 2% in temperature. The precise measurement of the calibration error directly impacts the uncertainties in the cosmological parameters estimated from the ACT power spectra. We also present a combined map based on ACT and WMAP data that has high signal-to-noise over a wide range of multipoles.

We present a measurement of the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 square degrees of the southern sky, in a 4.2-degree-wide strip centered on declination 53 degrees South. The CMB at arcminute angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio sources and dusty galaxies. After masking the 108 brightest point sources in our maps, we estimate the power spectrum between 600 < \ell < 8000 using the adaptive multi-taper method to minimize spectral leakage and maximize use of the full data set. Our absolute calibration is based on observations of Uranus. To verify the calibration and test the fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and recover the WMAP power spectrum from 250 < ell < 1150. The power beyond the Silk damping tail of the CMB is consistent with models of the emission from point sources. We quantify the contribution of SZ clusters to the power spectrum by fitting to a model normalized at sigma8 = 0.8. We constrain the model's amplitude ASZ < 1.63 (95% CL). If interpreted as a measurement of sigma8, this implies sigma8^SZ < 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to a 6-parameter LCDM model plus terms for point sources and the SZ effect is consistent with these results.

We use new kinematic data from the ultra-faint Milky Way satellite Segue 1 to model its dark matter distribution and derive upper limits on the dark matter annihilation cross-section. Using gamma-ray flux upper limits from the Fermi satellite and MAGIC, we determine cross-section exclusion regions for dark matter annihilation into a variety of different particles including charged leptons. We show that these exclusion regions are beginning to probe the regions of interest for a dark matter interpretation of the electron and positron fluxes from PAMELA, Fermi, and HESS, and that future observations of Segue 1 have strong prospects for testing such an interpretation. We additionally discuss prospects for detecting annihilation with neutrinos using the IceCube detector, finding that in an optimistic scenario a few neutrino events may be detected. Finally we use the kinematic data to model the Segue 1 dark matter velocity dispersion and constrain Sommerfeld enhanced models.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We create realistic, full-sky, half-arcminute resolution simulations of the microwave sky matched to the most recent astrophysical observations. The primary purpose of these simulations is to test the data reduction pipeline for the Atacama Cosmology Telescope (ACT) experiment; however, we have widened the frequency coverage beyond the ACT bands to make these simulations applicable to other microwave background experiments. Some of the novel features of these simulations are that the radio and infrared galaxy populations are correlated with the galaxy cluster populations, the CMB is lensed by the dark matter structure in the simulation via a ray-tracing code, the contribution to the thermal and kinetic Sunyaev-Zel'dovich (SZ) signals from galaxy clusters, groups, and the IGM has been included, and the gas prescription to model the SZ signals matches the most recent X-ray observations. Regarding the contamination of cluster SZ flux by radio galaxies, we find for 148 GHz (90 GHz) only 3% (4%) of halos have their SZ decrements contaminated at a level of 20% or more. We find the contamination levels higher for infrared galaxies. However, at 90 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<1.2 have their SZ decrements filled in at a level of 20% or more. At 148 GHz, less than 20% of clusters with M_{200} > 2.5 x 10^{14} Msun and z<0.8 have their SZ decrements filled in at a level of 50% or larger. Our models also suggest that a population of very high flux infrared galaxies, which are likely lensed sources, contribute most to the SZ contamination of very massive clusters at 90 and 148 GHz. These simulations are publicly available and should serve as a useful tool for microwave surveys to cross-check SZ cluster detection, power spectrum, and cross-correlation analyses.

We provide conservative bounds on the dark matter cross-section and lifetime from final state radiation produced by annihilation or decay into charged leptons, either directly or via an intermediate particle $\phi$. Our analysis utilizes the experimental gamma-ray flux upper limits from four Milky Way dwarf satellites: HESS observations of Sagittarius and VERITAS observations of Draco, Ursa Minor, and Willman 1. Using 90% confidence level lower limits on the integrals over the dark matter distributions, we find that these constraints are largely unable to rule out dark matter annihilations or decays as an explanation of the PAMELA and ATIC/PPB-BETS excesses. However, if there is an additional Sommerfeld enhancement in dwarfs, which have a velocity dispersion ~10 to 20 times lower than that of the local Galactic halo, then the cross-sections for dark matter annihilating through $\phi$'s required to explain the excesses are very close to the cross-section upper bounds from Willman 1. Dark matter annihilation directly into $\tau$'s is also marginally ruled out by Willman 1 as an explanation of the excesses, and the required cross-section is only a factor of a few below the upper bound from Draco. Finally, we make predictions for the gamma-ray flux expected from the dwarf galaxy Segue 1 for the Fermi Gamma-ray Space Telescope. We find that for a sizeable fraction of the parameter space in which dark matter annihilation into charged leptons explains the PAMELA excess, Fermi has good prospects for detecting a gamma-ray signal from Segue 1 after one year of observation.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

We compare X-ray and weak-lensing masses for four galaxy clusters that comprise the top-ranked shear-selected cluster system in the Deep Lens Survey. The weak-lensing observations of this system, which is associated with A781, are from the Kitt Peak Mayall 4-m telescope, and the X-ray observations are from both Chandra and XMM-Newton. For a faithful comparison of masses, we adopt the same matter density profile for each method, which we choose to be an NFW profile. Since neither the X-ray nor weak-lensing data are deep enough to well constrain both the NFW scale radius and central density, we estimate the scale radius using a fitting function for the concentration derived from cosmological hydrodynamic simulations and an X-ray estimate of the mass assuming isothermality. We keep this scale radius in common for both X-ray and weak-lensing profiles, and fit for the central density, which scales linearly with mass. We find that for three of these clusters, there is agreement between X-ray and weak-lensing NFW central densities, and thus masses. For the other cluster, the X-ray central density is higher than that from weak-lensing by 2 sigma. X-ray images suggest that this cluster may be undergoing a merger with a smaller cluster. This work serves as an additional step towards understanding the possible biases in X-ray and weak-lensing cluster mass estimation methods. Such understanding is vital to efforts to constrain cosmology using X-ray or weak-lensing cluster surveys to trace the growth of structure over cosmic time.

We study the ability of three-frequency, arcminute-resolution microwave measurements to detect galaxy clusters via their Sunyaev-Zel'dovich (SZ) distortion of the microwave background. For this purpose, we have constructed large-area simulations of the microwave sky, and we have made them publicly available to further investigations into optimal data reduction techniques for upcoming SZ cluster surveys. In these sky simulations, galaxy clusters are modeled using N-body simulated dark matter halos plus a gas prescription for the intracluster medium that allows the small scale cluster physics such as star formation and feedback to be realistically incorporated. We also model the primary microwave background, radio and infrared point sources, galactic dust emission, and the SZ flux including kinetic and relativistic contributions. We make use of these simulations to study the cluster scaling relation between integrated SZ flux and cluster mass and find our clusters fit a power-law well, with a power-law index that is steeper than that for self-similar cluster models. Some evolution of the power-law index and normalization with redshift is also observed. These simulations are also used to study cluster detection for the Atacama Cosmology Telescope (ACT). Using a multi-frequency Wiener filter to separate clusters from other microwave components, we find that ACT alone can recover a cluster sample that is ~90% complete above 3*10^{14} M_sun and ~85% pure.

Near-future SZ surveys such as ACT, SPT, APEX, and Planck will find thousands of galaxy clusters. Multi-frequency arcminute-resolution SZ observations can, in principle, determine each cluster's gas temperature (T_e), bulk velocity (v_pec), and optical depth (tau). However, the frequency bands and detector sensitivity employed by upcoming surveys will generally not be sufficient to disentangle the degeneracy between these three cluster parameters, even in the absence of SZ signal contamination from point sources and imperfect primary microwave background subtraction. Assuming contaminants can be removed, we find that near-future SZ surveys will be able to constrain well two cluster gas parameters that are linear combinations of tau*T_e, tau*v_z, and tau*T_e^2. Because the SZ intensity shift is nearly a linear function of tau*T_e, tau*v_z, and tau*T_e^2, a correspondence exists between the two effective gas parameters that SZ surveys can constrain and simple line-of-sight integrals through the three dimensional cluster. We illustrate the parameter constraints and correspondence to line-of-sight integrals using 3D Nbody + hydro cluster simulations and a Markov chain Monte Carlo method. We show that adding an independent T_e measurement to upcoming SZ data breaks the parameter degeneracy. The cluster effective velocity thus constrained is approximately the optical-depth-weighted velocity integrated along the cluster line of sight. A temperature prior with an error as large as 2 keV still gives bulk velocity errors < 100 km/sec, even for a more typical cluster with an electron temperature of 3 keV, for ACT-like SZ observations in the absence of signal contamination. The Markov chain constraints on v_pec and tau are more encouraging and likely more accurate than those obtained from Fisher matrices.