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

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

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne telescope mission to search for inflationary gravitational waves from the early universe. PIPER employs two 32x40 arrays of superconducting transition-edge sensors, which operate at 100 mK. An open bucket dewar of liquid helium maintains the receiver and telescope optics at 1.7 K. We describe the thermal design of the receiver and sub-kelvin cooling with a continuous adiabatic demagnetization refrigerator (CADR). The CADR operates between 70-130 mK and provides ~10 uW cooling power at 100 mK, nearly five times the loading of the two detector assemblies. We describe electronics and software to robustly control the CADR, overall CADR performance in flight-like integrated receiver testing, and practical considerations for implementation in the balloon float environment.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

Voyage 2050 White Paper highlighting the unique science opportunities using spectral distortions of the cosmic microwave background (CMB). CMB spectral distortions probe many processes throughout the history of the Universe. Precision spectroscopy, possible with existing technology, would provide key tests for processes expected within the cosmological standard model and open an enormous discovery space to new physics. This offers unique scientific opportunities for furthering our understanding of inflation, recombination, reionization and structure formation as well as dark matter and particle physics. A dedicated experimental approach could open this new window to the early Universe in the decades to come, allowing us to turn the long-standing upper distortion limits obtained with COBE/FIRAS some 25 years ago into clear detections of the expected standard distortion signals.

Experience suggests that structural issues in how institutional Astrophysics approaches data-driven science and the development of discovery technology may be hampering the community's ability to respond effectively to a rapidly changing environment in which increasingly complex, heterogeneous datasets are challenging our existing information infrastructure and traditional approaches to analysis. We stand at the confluence of a new epoch of multimessenger science, remote co-location of data and processing power and new observing strategies based on miniaturized spacecraft. Significant effort will be required by the community to adapt to this rapidly evolving range of possible discovery moduses. In the suggested creation of a new Astrophysics element, Advanced Astrophysics Discovery Technology, we offer an affirmative solution that places the visibility of discovery technologies at a level that we suggest is fully commensurate with their importance to the future of the field.

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.

New determinations are presented of the cosmic infrared background monopole brightness in the Planck HFI bands from 100 GHz to 857 GHz. Planck was not designed to measure the monopole component of sky brightness, so cross-correlation of the 2015 HFI maps with COBE/FIRAS data is used to recalibrate the zero level of the HFI maps. For the HFI 545 and 857 GHz maps, the brightness scale is also recalibrated. Correlation of the recalibrated HFI maps with a linear combination of Galactic H I and H alpha data is used to separate the Galactic foreground emission and determine the cosmic infrared background brightness in each of the HFI bands. We obtain CIB values of 0.007 +- 0.014, 0.010 +- 0.019, 0.060 +- 0.023, 0.149 +- 0.017, 0.371 +- 0.018, and 0.576 +- 0.034 MJy/sr at 100, 143, 217, 353, 545, and 857 GHz, respectively. The estimated uncertainties for the 353 to 857 GHz bands are about 3 to 6 times smaller than those of previous direct CIB determinations at these frequencies. Our results are compared with integrated source brightness results from selected recent submillimeter and millimeter wavelength imaging surveys.

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

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

CII is one of the brightest emission lines from star-forming galaxies and is an excellent tracer for star formation. Recent work measured the CII emission line amplitude for redshifts $2<z<3.2$ by cross-correlating Planck High Frequency Instrument emission maps with tracers of overdensity from the Baryon Oscillation Spectroscopic Sky Survey, finding $I_{\mathrm{CII}}=6.6^{+5.0}_{-4.8}\times 10^4$ Jy/sr at $95\%$ confidence. In this paper, we present a refinement of this earlier work by improving the mask weighting in each of the Planck bands and the precision in the covariance matrix. We report a detection of excess emission in the 545 GHz Planck band separate from the cosmic infrared background (CIB) present in the 353-857 GHz Planck bands. This excess is consistent with redshifted CII emission, in which case we report $b_{\mathrm{CII}}I_{\mathrm{CII}} =2.0^{+1.2}_{-1.1}\times 10^5$ Jy/sr at $95\%$ confidence, which strongly favors many collisional excitation models of CII emission. Our detection shows strong evidence for a model with a non-zero CII parameter, though line intensity mapping observations at high spectral resolution will be needed to confirm this result.

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.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

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.

Line-Intensity Mapping is an emerging technique which promises new insights into the evolution of the Universe, from star formation at low redshifts to the epoch of reionization and cosmic dawn. It measures the integrated emission of atomic and molecular spectral lines from galaxies and the intergalactic medium over a broad range of frequencies, using instruments with aperture requirements that are greatly relaxed relative to surveys for single objects. A coordinated, comprehensive, multi-line intensity-mapping experimental effort can efficiently probe over 80% of the volume of the observable Universe - a feat beyond the reach of other methods. Line-intensity mapping will uniquely address a wide array of pressing mysteries in galaxy evolution, cosmology, and fundamental physics. Among them are the cosmic history of star formation and galaxy evolution, the compositions of the interstellar and intergalactic media, the physical processes that take place during the epoch of reionization, cosmological inflation, the validity of Einstein's gravity theory on the largest scales, the nature of dark energy and the origin of dark matter.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

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.

Intensity mapping has attracted significant interest as an approach to measure the properties of the interstellar medium in typical galaxies at high redshift. Intensity mapping measures the statistics of surface brightness as a function of frequency, making it sensitive not only to all line emission of interest but also radiation from all other sources. Significant effort has gone into developing approaches that reject foreground contamination. Additionally, the target galaxies have multiple sources of emission that can complicate the interpretation of the line brightness. We describe the problem of jointly estimating correlated continuum emission and cleaning uncorrelated continuum emission, such as from the Milky Way. We apply these considerations to a cross-correlation of Planck data with BOSS quasars for a determination of CII for 2 < z < 3.2. Intensity mapping surveys with few bands have unique challenges for treating foregrounds and avoiding bias from correlated continuum emission. We show how a future intensity mapping survey with many bands can separate line from continuum emission in cross-correlation.

We report on measurements of the polarization of extragalactic sources at 148 GHz made during the first two seasons of the Atacama Cosmology Telescope Polarization (ACTPol) survey. The survey covered 680 deg$^{2}$ of the sky on the celestial equator. Polarization measurements of 169 intensity-selected sources brighter than 30 mJy, that are predominantly Active Galactic Nuclei, are presented. Above a total flux of 215 mJy where the noise bias removal in the polarization measurement is reliable, we detect 26 sources, 14 of which have a detection of linear polarization at greater than 3$\sigma_{p}$ significance. The distribution of the fractional polarization as a function of total source intensity is analyzed. Our result is consistent with the scenario that the fractional polarization of our measured radio source population is independent of total intensity down to the limits of our measurements and well described by a Gaussian distribution with a mean fractional polarization $p=0.028\pm$0.005 and standard deviation $\sigma_{\mathrm{p}}=0.054$, truncated at $p=0$. Extrapolating this model for the distribution of source polarization below the ACTPol detection threshold, we predict that one could get a clean measure of the E-mode polarization power spectrum of the microwave background out to $\ell=6000$ with 1 $\mu$K-arcminute maps over 10$\%$ of the sky from a future survey. We also study the spectral energy distribution of the total and polarized source flux densities by cross-matching with low radio frequency catalogs. We do not find any correlation between the spectral indices for total flux and polarized flux.

The sky-averaged (global) highly redshifted 21-cm spectrum from neutral hydrogen is expected to appear in the VHF range of $\sim20-200$ MHz and its spectral shape and strength are determined by the heating properties of the first stars and black holes, by the nature and duration of reionization, and by the presence or absence of exotic physics. Measurements of the global signal would therefore provide us with a wealth of astrophysical and cosmological knowledge. However, the signal has not yet been detected because it must be seen through strong foregrounds weighted by a large beam, instrumental calibration errors, and ionospheric, ground and radio-frequency-interference effects, which we collectively refer to as "systematics". Here, we present a signal extraction method for global signal experiments which uses Singular Value Decomposition (SVD) of "training sets" to produce systematics basis functions specifically suited to each observation. Instead of requiring precise absolute knowledge of the systematics, our method effectively requires precise knowledge of how the systematics can vary. After calculating eigenmodes for the signal and systematics, we perform a weighted least square fit of the corresponding coefficients and select the number of modes to include by minimizing an information criterion. We compare the performance of the signal extraction when minimizing various information criteria and find that minimizing the Deviance Information Criterion (DIC) most consistently yields unbiased fits. The methods used here are built into our widely applicable, publicly available Python package, $\texttt{pylinex}$, which analytically calculates constraints on signals and systematics from given data, errors, and training sets.

We report results from 21-cm intensity maps acquired from the Parkes radio telescope and cross-correlated with galaxy maps from the 2dF galaxy survey. The data span the redshift range $0.057<z<0.098$ and cover approximately 1,300 square degrees over two long fields. Cross correlation is detected at a significance of $5.18\sigma$. The amplitude of the cross-power spectrum is low relative to the expected dark matter power spectrum, assuming a neutral hydrogen (HI) bias and mass density equal to measurements from the ALFALFA survey. The decrement is pronounced and statistically significant at small scales. At $k\sim1.5$ $ h \mathrm{Mpc^{-1}}$, the cross power spectrum is more than a factor of 6 lower than expected, with a significance of $14.8\,\sigma$. This decrement indicates either a lack of clustering of neutral hydrogen (HI), a small correlation coefficient between optical galaxies and HI, or some combination of the two. Separating 2dF into red and blue galaxies, we find that red galaxies are much more weakly correlated with HI on $k\sim1.5$ $h \mathrm{Mpc^{-1}}$ scales, suggesting that HI is more associated with blue star-forming galaxies and tends to avoid red galaxies.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

CMB-S4 is a proposed experiment to map the polarization of the Cosmic Microwave Background (CMB) to nearly the cosmic variance limit for angular scales that are accessible from the ground. The science goals and capabilities of CMB-S4 in illuminating cosmic inflation, measuring the sum of neutrino masses, searching for relativistic relics in the early universe, characterizing dark energy and dark matter, and mapping the matter distribution in the universe have been described in the CMB-S4 Science Book. This Technology Book is a companion volume to the Science Book. The ambitious science goals of CMB-S4, a "Stage-4" experiment, require a step forward in experimental capability from the current Stage=II experiments. To guide this process, we summarize the current state of CMB instrumentation technology, and identify R&D efforts necessary to advance it for use in CMB-S4. The book focuses on technical challenges in four broad areas: Telescope Design; Receiver Optics; Focal-Plane Optical Coupling; and Focal-Plane Sensor and Readout.

The redshifted 21-cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies ($10<z<35$). The global 21-cm signal is sensitive to the thermal and ionization state of hydrogen gas and thus provides a tracer of sources of energetic photons -- primarily hot stars and accreting black holes -- which ionize and heat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the global spectrum with a realizable instrument placed in a low altitude lunar orbit, performing night-time 40-120 MHz spectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption, and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of the redshifted 21-cm spectrum are distinct from the spectrally featureless, spatially-varying, and polarized emission from the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. For signal extraction, we model the foreground, instrument, and 21-cm spectrum with eigenmodes calculated via Singular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameter space defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured and how astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence of foregrounds using the Dark Ages Radio Explorer (DARE).

A full account of galaxy evolution in the context of LCDM cosmology requires measurements of the average star-formation rate (SFR) and cold gas abundance across cosmic time. Emission from the CO ladder traces cold gas, and [CII] fine structure emission at 158 um traces the SFR. Intensity mapping surveys the cumulative surface brightness of emitting lines as a function of redshift, rather than individual galaxies. CMB spectral distortion instruments are sensitive to both the mean and anisotropy of the intensity of redshifted CO and [CII] emission. Large-scale anisotropy is proportional to the product of the mean surface brightness and the line luminosity-weighted bias. The bias provides a connection between galaxy evolution and its cosmological context, and is a unique asset of intensity mapping. Cross-correlation with galaxy redshift surveys allows unambiguous measurements of redshifted line brightness despite residual continuum contamination and interlopers. Measurement of line brightness through cross-correlation also evades cosmic variance and suggests new observation strategies. Galactic foreground emission is $\sim 10^3$ times larger than the expected signals, and this places stringent requirements on instrument calibration and stability. Under a range of assumptions, a linear combination of bands cleans continuum contamination sufficiently that residuals produce a modest penalty over the instrumental noise. For PIXIE, the $2 \sigma$ sensitivity to CO and [CII] emission scales from $\sim 5 \times 10^{-2}$ kJy/sr at low redshift to ~2 kJy/sr by reionization.

We present the first application of a new foreground removal pipeline to the current leading HI intensity mapping dataset, obtained by the Green Bank Telescope (GBT). We study the 15hr and 1hr field data of the GBT observations previously presented in Masui et al (2013) and Switzer et al (2013), covering about 41 square degrees at 0.6<z<1.0, for which cross-correlations may be measured with the galaxy distribution of the WiggleZ Dark Energy Survey. In the presented pipeline, we subtract the Galactic foreground continuum and the point source contamination using an independent component analysis technique (fastica), and develop a Fourier-based optimal estimator to compute the temperature power spectrum of the intensity maps and cross-correlation with the galaxy survey data. We show that fastica is a reliable tool to subtract diffuse and point-source emission through the non-Gaussian nature of their probability distributions. The temperature power spectra of the intensity maps is dominated by instrumental noise on small scales which fastica, as a conservative subtraction technique of non-Gaussian signals, can not mitigate. However, we determine similar GBT-WiggleZ cross-correlation measurements to those obtained by the Singular Value Decomposition (SVD) method, and confirm that foreground subtraction with fastica is robust against 21cm signal loss, as seen by the converged amplitude of these cross-correlation measurements. We conclude that SVD and fastica are complementary methods to investigate the foregrounds and noise systematics present in intensity mapping datasets.

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

The B-mode polarization of the cosmic microwave background (CMB) provides a unique window into tensor perturbations from inflationary gravitational waves. Survey effects complicate the estimation and description of the power spectrum on the largest angular scales. The pixel-space likelihood yields parameter distributions without the power spectrum as an intermediate step, but it does not have the large suite of tests available to power spectral methods. Searches for primordial B-modes must rigorously reject and rule out contamination. Many forms of contamination vary or are uncorrelated across epochs, frequencies, surveys, or other data treatment subsets. The cross power and the power spectrum of the difference of subset maps provide approaches to reject and isolate excess variance. We develop an analogous joint pixel-space likelihood. Contamination not modeled in the likelihood produces parameter-dependent bias and complicates the interpretation of the difference map. We describe a null test that consistently weights the difference map. Excess variance should either be explicitly modeled in the covariance or be removed through reprocessing the data.

The Primordial Inflation Polarization ExploreR (PIPER) is a balloon-borne telescope designed to measure the polarization of the Cosmic Microwave Background on large angular scales. PIPER will map 85% of the sky at 200, 270, 350, and 600 GHz over a series of 8 conventional balloon flights from the northern and southern hemispheres. The first science flight will use two 32x40 arrays of backshort-under-grid transition edge sensors, multiplexed in the time domain, and maintained at 100 mK by a Continuous Adiabatic Demagnetization Refrigerator. Front-end cryogenic Variable-delay Polarization Modulators provide systematic control by rotating linear to circular polarization at 3 Hz. Twin telescopes allow PIPER to measure Stokes I, Q, U, and V simultaneously. The telescope is maintained at 1.5 K in an LHe bucket dewar. Cold optics and the lack of a warm window permit sensitivity at the sky-background limit. The ultimate science target is a limit on the tensor-to-scalar ratio of r ~ 0.007, from the reionization bump to l ~ 300. PIPER's first flight will be from the Northern hemisphere, and overlap with the CLASS survey at lower frequencies. We describe the current status of the PIPER instrument.

Intensity mapping experiments survey the spectrum of diffuse line radiation rather than detect individual objects at high signal-to-noise. Spectral maps of unresolved atomic and molecular line radiation contain three-dimensional information about the density and environments of emitting gas, and efficiently probe cosmological volumes out to high redshift. Intensity mapping survey volumes also contain all other sources of radiation at the frequencies of interest. Continuum foregrounds are typically ~10^2-10^3 times brighter than the cosmological signal. The instrumental response to bright foregrounds will produce new spectral degrees of freedom that are not known in advance, nor necessarily spectrally smooth. The intrinsic spectra of foregrounds may also not be well-known in advance. We describe a general class of quadratic estimators to analyze data from single-dish intensity mapping experiments, and determine contaminated spectral modes from the data itself. The key attribute of foregrounds is not that they are spectrally smooth, but instead that they have fewer bright spectral degrees of freedom than the cosmological signal. Spurious correlations between the signal and foregrounds produce additional bias. Compensation for signal attenuation must estimate and correct this bias. A successful intensity mapping experiment will control instrumental systematics that spread variance into new modes, and it must observe a large enough volume that contaminant modes can be determined independently from the signal on scales of interest.

Variable-delay Polarization Modulators (VPMs) are currently being implemented in experiments designed to measure the polarization of the cosmic microwave background on large angular scales because of their capability for providing rapid, front-end polarization modulation and control over systematic errors. Despite the advantages provided by the VPM, it is important to identify and mitigate any time-varying effects that leak into the synchronously modulated component of the signal. In this paper, the effect of emission from a $300$ K VPM on the system performance is considered and addressed. Though instrument design can greatly reduce the influence of modulated VPM emission, some residual modulated signal is expected. VPM emission is treated in the presence of rotational misalignments and temperature variation. Simulations of time-ordered data are used to evaluate the effect of these residual errors on the power spectrum. The analysis and modeling in this paper guides experimentalists on the critical aspects of observations using VPMs as front-end modulators. By implementing the characterizations and controls as described, front-end VPM modulation can be very powerful for mitigating $1/f$ noise in large angular scale polarimetric surveys. None of the systematic errors studied fundamentally limit the detection and characterization of B-modes on large scales for a tensor-to-scalar ratio of $r=0.01$. Indeed, $r<0.01$ is achievable with commensurately improved characterizations and controls.

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

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

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

Spectral measurements of the 21 cm monopole background have the promise of revealing the bulk energetic properties and ionization state of our universe from z ~ 6-30. Synchrotron foregrounds are orders of magnitude larger than the cosmological signal, and are the principal challenge faced by these experiments. While synchrotron radiation is thought to be spectrally smooth and described by relatively few degrees of freedom, the instrumental response to bright foregrounds may be much more complex. To deal with such complexities, we develop an approach that discovers contaminated spectral modes using spatial fluctuations of the measured data. This approach exploits the fact that foregrounds vary across the sky while the signal does not. The discovered modes are projected out of each line-of-sight of a data cube. An angular weighting then optimizes the cosmological signal amplitude estimate by giving preference to lower-noise regions. Using this method, we show that it is essential for the passband to be stable to at least ~10^{-4}. In contrast, the constraints on the spectral smoothness of the absolute calibration are mainly aesthetic if one is able to take advantage of spatial information. To the extent it is understood, controlling polarization to intensity leakage at the ~10^{-2} level will also be essential to rejecting Faraday rotation of the polarized synchrotron emission.

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.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne cosmic microwave background (CMB) polarimeter designed to search for evidence of inflation by measuring the large-angular scale CMB polarization signal. BICEP2 recently reported a detection of B-mode power corresponding to the tensor-to-scalar ratio r = 0.2 on ~2 degree scales. If the BICEP2 signal is caused by inflationary gravitational waves (IGWs), then there should be a corresponding increase in B-mode power on angular scales larger than 18 degrees. PIPER is currently the only suborbital instrument capable of fully testing and extending the BICEP2 results by measuring the B-mode power spectrum on angular scales $\theta$ = ~0.6 deg to 90 deg, covering both the reionization bump and recombination peak, with sensitivity to measure the tensor-to-scalar ratio down to r = 0.007, and four frequency bands to distinguish foregrounds. PIPER will accomplish this by mapping 85% of the sky in four frequency bands (200, 270, 350, 600 GHz) over a series of 8 conventional balloon flights from the northern and southern hemispheres. The instrument has background-limited sensitivity provided by fully cryogenic (1.5 K) optics focusing the sky signal onto four 32x40-pixel arrays of time-domain multiplexed Transition-Edge Sensor (TES) bolometers held at 140 mK. Polarization sensitivity and systematic control are provided by front-end Variable-delay Polarization Modulators (VPMs), which rapidly modulate only the polarized sky signal at 3 Hz and allow PIPER to instantaneously measure the full Stokes vector (I, Q, U, V) for each pointing. We describe the PIPER instrument and progress towards its first flight.

We describe the design, construction, and initial validation of the variable-delay polarization modulator (VPM) designed for the PIPER cosmic microwave background polarimeter. The VPM modulates between linear and circular polarization by introducing a variable phase delay between orthogonal linear polarizations. Each VPM has a diameter of 39 cm and is engineered to operate in a cryogenic environment (1.5 K). We describe the mechanical design and performance of the kinematic double-blade flexure and drive mechanism along with the construction of the high precision wire grid polarizers.

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

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

The large-scale distribution of neutral hydrogen in the Universe will be luminous through its 21 cm emission. Here, for the first time, we use the auto-power spectrum of 21 cm intensity fluctuations to constrain neutral hydrogen fluctuations at z~0.8. Our data were acquired with the Green Bank Telescope and span the redshift range 0.6 < z < 1 over two fields totalling ~41 deg. sq. and 190 h of radio integration time. The dominant synchrotron foregrounds exceed the signal by ~10^3, but have fewer degrees of freedom and can be removed efficiently. Even in the presence of residual foregrounds, the auto-power can still be interpreted as an upper bound on the 21 cm signal. Our previous measurements of the cross-correlation of 21 cm intensity and the WiggleZ galaxy survey provide a lower bound. Through a Bayesian treatment of signal and foregrounds, we can combine both fields in auto- and cross-power into a measurement of Omega_HI b_HI = [0.62^{+0.23}_{-0.15}] * 10^{-3} at 68% confidence with 9% systematic calibration uncertainty, where Omega_HI is the neutral hydrogen (HI) fraction and b_HI is the HI bias parameter. We describe observational challenges with the present data set and plans to overcome them.

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

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

Future high resolution, high sensitivity Sunyaev-Zeldovich (SZ) observations of individual clusters will provide an exciting opportunity to answer specific questions about the dynamical state of the intra-cluster medium (ICM). In this paper we develop a new method that clearly shows the connection of the SZ signal with the underlying cluster model. We include relativistic temperature and kinematic corrections in the single-scattering approximation, allowing studies of hot clusters. In our approach, particular moments of the temperature and velocity field along the line-of-sight determine the precise spectral shape and morphology of the SZ signal. We illustrate how to apply our method to different cluster models, highlighting parameter degeneracies and instrumental effects that are important for interpreting future high-resolution SZ data. Our analysis shows that line-of-sight temperature variations can introduce significant biases in the derived SZ temperature and peculiar velocity. We furthermore discuss how the position of the SZ null is affected by the cluster's temperature and velocity structure. Our computations indicate that the SZ signal around the null alone is rather insensitive to different cluster models and that high frequency channels add a large leverage in this respect. We also apply our method to recent high sensitivity SZ data of the Bullet cluster, showing how the results can be linked to line-of-sight variations in the electron temperature. The tools developed here as part of SZpack should be useful for analyzing high-resolution SZ data and computing SZ maps from simulated clusters.

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.

The primordial CMB at small angular scales is sensitive to the ionization and expansion history of the universe around the time of recombination. This dependence has been exploited to constrain the helium abundance and the effective number of relativistic species. Here we focus on allowed ionization fraction trajectories, $\Xe (z)$, by constraining low-order principal components of perturbations to the standard recombination scenario ($\Xe$-eigenmodes) in the circa 2011 SPT, ACT and WMAP7 data. Although the trajectories are statistically consistent with the standard recombination, we find that there is a tension similar to that found by varying the helium fraction. As this paper was in press, final SPT and ACT datasets were released and we applied our framework to them: we find the tension continues, with slightly higher significance, in the new 2012 SPT data, but find no tension with the standard model of recombination in the new 2012 ACT data. We find that the prior probabilities on the eigenamplitudes are substantially influenced by the requirement that $\Xe$ trajectories conserve electron number. We propose requiring a sufficient entropy decrease between posterior and prior marginalized distributions be used as an $\Xe$-mode selection criterion. We find that in the case of the 2011 SPT/ACT+WMAP7 data only two modes are constrainable, but upcoming ACTPol, Planck and SPTPol data will be able to test more modes and more precisely address the current tension.

In this letter, 21 cm intensity maps acquired at the Green Bank Telescope are cross-correlated with large-scale structure traced by galaxies in the WiggleZ Dark Energy Survey. The data span the redshift range 0.6 < z < 1 over two fields totaling ~41 deg. sq. and 190 hours of radio integration time. The cross-correlation constrains Omega_HI b_HI r = [0.43 \pm 0.07 (stat.) \pm 0.04(sys.)] x 10^-3, where Omega_HI is the neutral hydrogen HI fraction, r is the galaxy-hydrogen correlation coefficient, and b_HI is the HI bias parameter. This is the most precise constraint on neutral hydrogen density fluctuations in a challenging redshift range. Our measurement improves the previous 21 cm cross-correlation at z ~ 0.8 both in its precision and in the range of scales probed.

We present the optical and X-ray properties of 68 galaxy clusters selected via the Sunyaev-Zel'dovich Effect at 148 GHz by the Atacama Cosmology Telescope (ACT). Our sample, from an area of 504 square degrees centered on the celestial equator, is divided into two regions. The main region uses 270 square degrees of the ACT survey that overlaps with the co-added ugriz imaging from the Sloan Digital Sky Survey (SDSS) over Stripe 82 plus additional near-infrared pointed observations with the Apache Point Observatory 3.5-meter telescope. We confirm a total of 49 clusters to z~1.3, of which 22 (all at z>0.55) are new discoveries. For the second region the regular-depth SDSS imaging allows us to confirm 19 more clusters up to z~0.7, of which 10 systems are new. We present the optical richness, photometric redshifts, and separation between the SZ position and the brightest cluster galaxy (BCG). We find no significant offset between the cluster SZ centroid and BCG location and a weak correlation between optical richness and SZ-derived mass. We also present X-ray fluxes and luminosities from the ROSAT All Sky Survey which confirm that this is a massive sample. One of the newly discovered clusters, ACT-CL J0044.4+0113 at z=1.1 (photometric), has an integrated XMM-Newton X-ray temperature of kT_x=7.9+/-1.0 keV and combined mass of M_200a=8.2(-2.5,+3.3)x10^14 M_sun/h70 placing it among the most massive and X-ray-hot clusters known at redshifts beyond z=1. We also highlight the optically-rich cluster ACT-CL J2327.4-0204 (RCS2 2327) at z=0.705 (spectroscopic) as the most significant detection of the whole equatorial sample with a Chandra-derived mass of M_200a=1.9(-0.4,+0.6)x10^15 M_sun/h70, comparable to some of the most massive known clusters like "El Gordo" and the Bullet Cluster.

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.

In this paper we refine the calculation of primordial helium recombination, accounting for several additional effects that were neglected or treated more approximately in previous studies. These include consideration of (i) time-dependent radiative transfer interaction between the 2^1 P_1 - 1^1 S_0 and 2^3 P_1 - 1^1 S_0 resonances; (ii) time-dependent radiative transfer for the partially overlapping n^1 P_1 - 1^1 S_0, n^1 D_2 - 1^1 S_0 and n^3 P_1 - 1^1 S_0 series with 3 <= n <= 10; (iii) electron scattering within a kernel approach. We also briefly discuss the effect of electron scattering and HI quadrupole lines on the recombination of hydrogen. Although the physics of all the considered processes is interesting and subtle, for the standard cosmology, with Y_p \sim 0.24, the overall correction to the ionization history during helium recombination with respect to the previous implementation of CosmoRec remains smaller than |DeltaNe/Ne| \sim 0.05%. The dominant improvement is caused by consistent inclusion of resonance scattering for the 2^1 P_1 - 1^1 S_0 resonance. For cosmologies with a large helium fraction, Y_p \sim 0.4, the difference reaches |D N_e/N_e|\sim 0.22% at z\sim1800, however, the overall correction to the CMB power spectra is small, exceeding |D C_l/C_l| \sim 0.05 % only at l>\sim 3000. In comparison to the current version of Recfast the difference reaches |D C_l/C_l|\sim0.4% at l\sim 3000 for Y_p \sim 0.4, and also for the standard value Y_p \sim 0.24 we find differences |DC_l/C_l|>\sim 0.1% at l\sim 2500. The new processes are now included by the cosmological recombination code CosmoRec and can be activated as needed for most settings without affecting its runtime significantly.

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.

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

Rogue planets have been ejected from their planetary system. We investigate the possibility that a rogue planet could maintain a liquid ocean under layers of thermally-insulating water ice and frozen atmosphere as a result of geothermal heat flux. We find that a rogue planet of Earth-like composition and age could maintain a subglacial liquid ocean if it were ~3.5 times more massive than Earth. If a rogue planet had about ten times higher water mass fraction or a thick cryo-atmospheric layer, it would need to be only ~0.3 times the mass of Earth to maintain a liquid ocean. Such a planet could be detected from reflected solar radiation and its thermal emission could be characterized in the far-IR if it passed within O(1000) AU of Earth.

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

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

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 a method of reliably extracting the flux of individual sources from sky maps in the presence of noise and a source population in which number counts are a steeply falling function of flux. The method is an extension of a standard Bayesian procedure in the millimeter/submillimeter literature. As in the standard method, the prior applied to source flux measurements is derived from an estimate of the source counts as a function of flux, dN/dS. The key feature of the new method is that it enables reliable extraction of properties of individual sources, which previous methods in the literature do not. We first present the method for extracting individual source fluxes from data in a single observing band, then we extend the method to multiple bands, including prior information about the spectral behavior of the source population(s). The multi-band estimation technique is particularly relevant for classifying individual sources into populations according to their spectral behavior. We find that proper treatment of the correlated prior information between observing bands is key to avoiding significant biases in estimations of multi-band fluxes and spectral behavior, biases which lead to significant numbers of misclassified sources. We test the single- and multi-band versions of the method using simulated observations with observing parameters similar to that of the South Pole Telescope data used in Vieira, et al. (2010).

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 discuss the potential of using the HeI 584 A forest to detect and study HeII reionization. Significant 584 A absorption is expected from intergalactic HeII regions, whereas there should be no detectable absorption from low density gas in HeIII regions. Unlike HeII Ly-alpha absorption (the subject of much recent study), the difficulty with using this transition to study HeII reionization is not saturation but rather that the absorption is weak. The Gunn-Peterson optical depth for this transition is tau ~ 0.1 x_{HeII} Delta^2 [(1+z)/5]^{9/2}, where x_{HeII} is the fraction of helium in HeII and Delta is the density in units of the cosmic mean. In addition, HeI 584 A absorption is contaminated by lower redshift HI Ly-alpha absorption with a comparable flux decrement. We estimate the requirements for a definitive detection of redshifted HeI absorption from low density gas (Delta ~ 1), which would indicate that HeII reionization was occurring. We find that this objective can be accomplished (using coeval HI Ly-alpha absorption to mask dense regions and in cross correlation) with a spectral resolution of 10^4 and a signal-to-noise ratio per resolution element of ~ 10. Such specifications may be achievable on a few known z ~ 3.5 quasar sightlines with the Cosmic Origins Spectrograph on the Hubble Space Telescope. We also discuss how HeI absorption can be used to measure the hardness of the ionizing background above 13.6 eV.

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

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

The Atacama Cosmology Telescope was designed to measure small-scale anisotropies in the Cosmic Microwave Background and detect galaxy clusters through the Sunyaev-Zel'dovich effect. The instrument is located on Cerro Toco in the Atacama Desert, at an altitude of 5190 meters. A six-meter off-axis Gregorian telescope feeds a new type of cryogenic receiver, the Millimeter Bolometer Array Camera. The receiver features three 1000-element arrays of transition-edge sensor bolometers for observations at 148 GHz, 218 GHz, and 277 GHz. Each detector array is fed by free space mm-wave optics. Each frequency band has a field of view of approximately 22' x 26'. The telescope was commissioned in 2007 and has completed its third year of operations. We discuss the major components of the telescope, camera, and related systems, and summarize the instrument performance.

Motivated by recent interest in redshifted 21 cm emission of intergalactic hydrogen, we investigate the 8.7 GHz ^2S_{1/2} F=0-1 hyperfine transition of ^3He^+. While the primordial abundance of 3He relative to hydrogen is 10^-5, the hyperfine spontaneous decay rate is 680 times larger. Furthermore, the antenna temperature is much lower at the frequencies relevant for the ^3He^+ transition compared to that of z>6 21 cm emission. We find that the spin temperature of this 8.7 GHz line in the intergalactic medium is approximately the cosmic microwave background temperature, such that this transition is best observed in absorption against high-redshift, radio-bright quasars. We show that intergalactic 8.7 GHz absorption is a promising, unsaturated observable of the ionization history of intergalactic helium (for which HeII->HeIII reionization is believed to be complete at z~3) and of the primordial 3He abundance. Instruments must reach ~1 \mu Jy RMS noise in bands of 1 MHz on a 1 Jy source to directly resolve this absorption. However, in combination with HI Ly\alpha forest measurements, an instrument can statistically detect this absorption from z > 3 with 30 \mu Jy RMS noise in 0.1 MHz spectral bands over 100 MHz, which may be within the reach of present instruments.

Upcoming precision measurements of the temperature anisotropy of the cosmic microwave background (CMB) at high multipoles will need to be complemented by a more complete understanding of recombination, which determines the damping of anisotropies on these scales. This is the third in a series of papers describing an accurate theory of HeI and HeII recombination. Here we describe the effect of Thomson scattering, the $^3$He isotope shift, the contribution of rare decays, collisional processes, and peculiar motion. These effects are found to be negligible: Thomson and $^3$He scattering modify the free electron fraction $x_e$ at the level of several $\times 10^{-4}$. The uncertainty in the $2^3P^o-1^1S$ rate is significant, and for conservative estimates gives uncertainties in $x_e$ of order $10^{-3}$. We describe several convergence tests for the atomic level code and its inputs, derive an overall $C_\ell$ error budget, and relate shifts in $x_e(z)$ to the changes in $C_\ell$, which are at the level of 0.5% at $\ell =3000$. Finally, we summarize the main corrections developed thus far. The remaining uncertainty from known effects is $\sim 0.3%$ in $x_e$.

Precision measurements of the cosmic microwave background temperature anisotropy on scales $\ell > 500$ will be available in the near future. Successful interpretation of these data is dependent on a detailed understanding of the damping tail and cosmological recombination of both hydrogen and helium. This paper and two companion papers are devoted to a precise calculation of helium recombination. We discuss several aspects of the standard recombination picture, and then include feedback, radiative transfer in HeI lines with partial redistribution, and continuum opacity from HI photoionization. In agreement with past calculations, we find that HeII recombination proceeds in Saha equilibrium, whereas HeI recombination is delayed relative to Saha due to the low rates connecting excited states of HeI to the ground state. However, we find that at $z<2200$ the continuum absorption by the rapidly increasing HI population becomes effective at destroying photons in the HeI $2^1P^o-1^1S$ line, causing HeI recombination to finish around $z\approx 1800$, much earlier than previously estimated.

Interpretation of precision measurements of the cosmic microwave background (CMB) will require a detailed understanding of the recombination era, which determines such quantities as the acoustic oscillation scale and the Silk damping scale. This paper is the second in a series devoted to the subject of helium recombination, with a focus on two-photon processes in He I. The standard treatment of these processes includes only the spontaneous two-photon decay from the 2^1S level. We extend this treatment by including five additional effects, some of which have been suggested in recent papers but whose impact on He I recombination has not been fully quantified. These are: (i) stimulated two-photon decays; (ii) two-photon absorption of redshifted HeI line radiation; (iii) two-photon decays from highly excited levels in HeI (n^1S and n^1D, with n>=3); (iv) Raman scattering; and (v) the finite width of the 2^1P^o resonance. We find that effect (iii) is highly suppressed when one takes into account destructive interference between different intermediate states contributing to the two-photon decay amplitude. Overall, these effects are found to be insignificant: they modify the recombination history at the level of several parts in 10^4.

The Atacama Cosmology Telescope is a 6-meter telescope designed to map the Cosmic Microwave Background simultaneously at 145 GHz, 215 GHz, and 280 GHz with arcminute resolution. Each frequency will have a 32 by 32 element focal plane array of TES bolometers. This paper describes the design of the telescope and the cold reimaging optics, which is optimized for millimeter-wave observations with these sensitive detectors.

It has been suggested that secondary CMB anisotropies generated by neutral lithium could open a new observational window into the universe around the redshift z ~ 400, and permit a determination of the primordial lithium abundance. The effect is due to resonant scattering in the allowed LiI doublet (2s ^2S_{1/2} - 2p ^2P_{1/2,3/2}), so its observability depends on the formation history of neutral lithium. Here we show that the ultraviolet photons produced during hydrogen recombination are sufficient to keep lithium in the LiII ionization stage in the relevant redshift range and suppress the neutral fraction by ~3 orders of magnitude from previous calculations, making the lithium signature unobservable.

Measurements of the linear power spectrum of galaxies have placed tight constraints on neutrino masses. We extend the framework of the halo model of cosmological nonlinear matter clustering to include the effect of massive neutrino infall into cold dark matter (CDM) halos. The magnitude of the effect of neutrino clustering for three degenerate mass neutrinos with m_nu=0.9 eV is of order ~1%, within the potential sensitivity of upcoming weak lensing surveys. In order to use these measurements to further constrain--or eventually detect--neutrino masses, accurate theoretical predictions of the nonlinear power spectrum in the presence of massive neutrinos will be needed, likely only possible through high-resolution multiple particle (neutrino, CDM and baryon) simulations.