We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions.

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA's Euclid, and NASA's Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.

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

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

CCAT-prime is a new 6 m crossed Dragone telescope designed to characterize the Cosmic Microwave Background (CMB) polarization and foregrounds, measure the Sunyaev-Zel'dovich effects of galaxy clusters, map the [CII] emission intensity from the Epoch of Reionization (EoR), and monitor accretion luminosity over multi-year timescales of hundreds of protostars in the Milky Way. CCAT-prime will make observations from a 5,600 m altitude site on Cerro Chajnantor in the Atacama Desert of northern Chile. The novel optical design of the telescope combined with high surface accuracy ($<$10 $\mu$m) mirrors and the exceptional atmospheric conditions of the site will enable sensitive broadband, polarimetric, and spectroscopic surveys at sub-mm to mm wavelengths. Prime-Cam, the first light instrument for CCAT-prime, consists of a 1.8 m diameter cryostat that can house seven individual instrument modules. Each instrument module, optimized for a specific science goal, will use state-of-the-art kinetic inductance detector (KID) arrays operated at $\sim$100 mK, and Fabry-Perot interferometers (FPI) for the EoR science. Prime-Cam will be commissioned with staged deployments to populate the seven instrument modules. The full instrument will consist of 60,000 polarimetric KIDs at a combination of 220/280/350/410 GHz, 31,000 KIDS at 250/360 GHz coupled with FPIs, and 21,000 polarimetric KIDs at 850 GHz. Prime-Cam is currently being built, and the CCAT-prime telescope is designed and under construction by Vertex Antennentechnik GmbH to achieve first light in 2021. CCAT-prime is also a potential telescope platform for the future CMB Stage-IV observations.

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

We present a new pipeline for the efficient generation of synthetic observations of the extragalactic microwave sky, tailored to large ground-based CMB experiments such as the Simons Observatory, Advanced ACTPol, SPT-3G, and CMB-S4. Such simulated observations are a key technical challenge in cosmology because of the dynamic range and accuracy required. The first part of the pipeline generates a random cosmological realization in the form of a dark matter halo catalog and matter displacement field, as seen from a given position. The halo catalog and displacement field are modeled with ellipsoidal collapse dynamics and Lagrangian perturbation theory, respectively. In the second part, the cosmological realization is converted into a set of intensity maps over the range 10 - 10^3 GHz using models based on existing observations and hydrodynamical simulations. These maps include infrared emission from dusty star forming galaxies (CIB), Comptonization of CMB photons by hot gas in groups and clusters through the thermal Sunyaev-Zel'dovich effect (tSZ), Doppler boosting by Thomson scattering of the CMB by bulk flows through the kinetic Sunyaev-Zel'dovich effect (kSZ), and weak gravitational lensing of primary CMB anisotropies by the large-scale distribution of matter in the universe. After describing the pipeline and its implementation, we present the Websky maps, created from a realization of the cosmic web on our past light cone in the redshift interval 0<z<4.6 over the full-sky and a volume of ~(600 Gpc/h)^3 resolved with ~10^12 resolution elements. The Websky maps and halo catalog are publicly available at mocks.cita.utoronto.ca/websky.

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

The small-scale cosmic microwave background (CMB) is dominated by anisotropies from the kinematic Sunyaev-Zeldovich (kSZ) effect, and upcoming experiments will measure it very precisely, but the optical depth degeneracy limits the cosmological information that can be extracted. At the same time, fast radio bursts (FRBs) are an exciting new frontier for astrophysics, but their usefulness as cosmological probes is currently unclear. We show that FRBs are uniquely suited for breaking the kSZ optical depth degeneracy. This opens up new possibilities for constraining cosmology with the kSZ effect, and new cosmological applications for FRBs.

One emerging application of machine learning methods is the inference of galaxy cluster masses. In this note, machine learning is used to directly combine five simulated multiwavelength measurements in order to find cluster masses. This is in contrast to finding mass estimates for each observable, normally by using a scaling relation, and then combining these scaling law based mass estimates using a likelihood. We also illustrate how the contributions of each observable to the accuracy of the resulting mass measurement can be compared via model-agnostic Importance Permutation values. Thirdly, as machine learning relies upon the accuracy of the training set in capturing observables, their correlations, and the observational selection function, and as the machine learning training set originates from simulations, two tests of whether a simulation's correlations are consistent with observations are suggested and explored as well.

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

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

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.

The galaxy cluster Zwicky 3146 is a sloshing cool core cluster at $z=0.291$ that in X-ray imaging does not appear to exhibit significant pressure substructure in the intracluster medium (ICM). The published $M_{500}$ values range between $3.88^{+0.62}_{-0.58}$ to $22.50 \pm 7.58 \times 10^{14}$ M$_{\odot}$, where ICM-based estimates with reported errors $<20$\% suggest that we should expect to find a mass between $6.53^{+0.44}_{-0.44} \times 10^{14}$ M$_{\odot}$ (from Planck, with an $8.4\sigma$ detection) and $8.52^{+1.77}_{-1.47} \times 10^{14}$ M$_{\odot}$ (from ACT, with a $14\sigma$ detection). This broad range of masses is suggestive that there is ample room for improvement for all methods. Here, we investigate the ability to estimate the mass of Zwicky 3146 via the Sunyaev-Zel'dovich (SZ) effect with data taken at 90 GHz by MUSTANG-2 to a noise level better than $15\ \mu$K at the center, and a cluster detection of $104\sigma$. We derive a pressure profile from our SZ data which is in excellent agreement with that derived from X-ray data. From our SZ-derived pressure profiles, we infer $M_{500}$ and $M_{2500}$ via three methods -- $Y$-$M$ scaling relations, the virial theorem, and hydrostatic equilibrium -- where we employ X-ray constraints from \emph{XMM-Newton} on the electron density profile when assuming hydrostatic equilibrium. Depending on the model and estimation method, our $M_{500}$ estimates range from $6.23 \pm 0.59$ to $10.6 \pm 0.95 \times 10^{14}$ M$_{\odot}$, where our estimate from hydrostatic equilibrium, is $8.29^{+1.93}_{-1.24}$ ($\pm 19.1$\% stat) ${}^{+0.74}_{-0.68}$ ($\pm 8.6$\% sys, calibration) $\times 10^{14}$ M$_{\odot}$. Our fiducial mass, derived from a $Y$-$M$ relation is $8.16^{+0.44}_{-0.54}$ ($\pm 5.5$\% stat) ${}^{+0.46}_{-0.43}$ ($\pm 5.5$\% sys, $Y$-$M$) ${}^{+0.59}_{-0.55}$ ($\pm 7.0$\% sys, cal.) $\times 10^{14}$ M$_{\odot}$.

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

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

Upcoming millimeter experiments that probe the cosmic microwave background (CMB) will observe tens of thousands of galaxy clusters through the thermal Sunyaev-Zeldovich (tSZ) effect. tSZ selected clusters are powerful probes of cosmological models, as they trace the late-time growth of structure. Late-time structure growth is highly sensitive to extensions to the standard cosmological model ($\Lambda$CDM), such as the sum of the neutrino masses, the dark energy equation of state, and modifications to general relativity. The nominal statistic used for cluster observations is their abundances as a function of redshift. We investigate what additional cosmological information is gained after including the clustering signal of clusters, the cluster power spectrum. We forecast the cluster power spectra for the upcoming Simons Observatory and a CMB Stage-4-like experiment and find that the cluster power spectrum reduces marginalized constraints on the dark energy equation of state by $3$--$12\%$ and the growth index by $5$--$8\%$, for example. We present the constraints using a generalized figure of merit and find improvements ranging from $12$--$18\%$ for extensions, $12$--$19\%$ for the astrophysical nuisance parameters, and $14$--$25\%$ for $\Lambda$CDM parameters. We also find that if the bias of clusters as a tracer of the matter density can be measured to within $3\%$ or better, these improvements can be increased by up to a factor of $6$. We discuss the possibility of utilizing the clustering signal to address specific systematic uncertainties present in cluster abundance measurements.

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

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

We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2\sigma$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.

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

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

Modern galaxy cluster science is a multi-wavelength endeavor with cornerstones provided by X-ray, optical/IR, mm, and radio measurements. In combination, these observations enable the construction of large, clean, complete cluster catalogs, and provide precise redshifts and robust mass calibration. The complementary nature of these multi-wavelength data dramatically reduces the impact of systematic effects that limit the utility of measurements made in any single waveband. The future of multi-wavelength cluster science is compelling, with cluster catalogs set to expand by orders of magnitude in size, and extend, for the first time, into the high-redshift regime where massive, virialized structures first formed. Unlocking astrophysical and cosmological insight from the coming catalogs will require new observing facilities that combine high spatial and spectral resolution with large collecting areas, as well as concurrent advances in simulation modeling campaigns. Together, future multi-wavelength observations will resolve the thermodynamic structure in and around the first groups and clusters, distinguishing the signals from active and star-forming galaxies, and unveiling the interrelated stories of galaxy evolution and structure formation during the epoch of peak cosmic activity.

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.

Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$\alpha$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$\alpha$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario.

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

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.

Achieving a precise understanding of galaxy formation in a cosmological context is one of the great challenges in theoretical astrophysics, due to the vast range of spatial scales involved in the relevant physical processes. Observations in the millimeter bands, particularly those using the cosmic microwave background (CMB) radiation as a "backlight", provide a unique probe of the thermodynamics of these processes, with the capability to directly measure the density, pressure, and temperature of ionized gas. Moreover, these observations have uniquely high sensitivity into the outskirts of the halos of galaxies and clusters, including systems at high redshift. In the next decade, the combination of large spectroscopic and photometric optical galaxy surveys and wide-field, low-noise CMB surveys will transform our understanding of galaxy formation via these probes.

Recent advances in cosmic observations have brought us to the verge of discovery of the absolute scale of neutrino masses. Nonzero neutrino masses are known evidence of new physics beyond the Standard Model. Our understanding of the clustering of matter in the presence of massive neutrinos has significantly improved over the past decade, yielding cosmological constraints that are tighter than any laboratory experiment, and which will improve significantly over the next decade, resulting in a guaranteed detection of the absolute neutrino mass scale.

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

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

Most secondary sources of cosmic microwave background anisotropy (radio sources, dusty galaxies -DSFGs-, thermal Sunyaev Zel'dovich -tSZ- distortions from hot gas, and gravitational lensing) are highly non-Gaussian. Statistics beyond the power spectrum are therefore potentially important sources of information about the physics of these processes. We combine data from the Atacama Cosmology Telescope and with data from the Planck satellite (only using Planck data in the overlapping region) to constrain the amplitudes of a set of theoretical bispectrum templates from the tSZ effect, DSFGs, gravitational lensing, and radio galaxies. We make a strong detection of radio galaxies (>5$\sigma$) and have hints of non-Gaussianity arising from the tSZ effect, DSFGs, from cross-correlations between the tSZ effect and DSFGs and from cross-correlations among the tSZ effect, DSFGs and radio galaxies. These results suggest that the same halos host radio sources, DSFGs, and have tSZ signal. We present a new method to calculate the non-Gaussian contributions to the template covariances. Using this method we find significant non-Gaussian contributions to the variance and covariance of our templates, with templates involving the tSZ effect most effected. Strong degeneracies exist between the various sources at the current noise levels. In light of these degeneracies, combined with theoretical uncertainty in the templates, these results are a demonstration of this technique. With these caveats, we demonstrate the utility of future bispectrum measurements by using the tSZ bispectrum measurement to constrain a combination of the amplitude of matter fluctuations and the matter density to be $\sigma_8 \Omega_m^{0.17}=0.65^{+0.05}_{-0.06}$. Improvements in signal to noise from upcoming observations will enable the separation of bispectrum components and robust constraints on cosmological parameters.

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

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

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

We present the detailed science case, and brief descriptions of the telescope design, site, and first light instrument plans for a new ultra-wide field submillimeter observatory, CCAT-prime, that we are constructing at a 5600 m elevation site on Cerro Chajnantor in northern Chile. Our science goals are to study star and galaxy formation from the epoch of reionization to the present, investigate the growth of structure in the Universe, improve the precision of B-mode CMB measurements, and investigate the interstellar medium and star formation in the Galaxy and nearby galaxies through spectroscopic, polarimetric, and broadband surveys at wavelengths from 200 um to 2 mm. These goals are realized with our two first light instruments, a large field-of-view (FoV) bolometer-based imager called Prime-Cam (that has both camera and an imaging spectrometer modules), and a multi-beam submillimeter heterodyne spectrometer, CHAI. CCAT-prime will have very high surface accuracy and very low system emissivity, so that combined with its wide FoV at the unsurpassed CCAT site our telescope/instrumentation combination is ideally suited to pursue this science. The CCAT-prime telescope is being designed and built by Vertex Antennentechnik GmbH. We expect to achieve first light in the spring of 2021.

CCAT-prime will be a 6-meter aperture telescope operating from sub-mm to mm wavelengths, located at 5600 meters elevation on Cerro Chajnantor in the Atacama Desert in Chile. Its novel crossed-Dragone optical design will deliver a high throughput, wide field of view capable of illuminating much larger arrays of sub-mm and mm detectors than can existing telescopes. We present an overview of the motivation and design of Prime-Cam, a first-light instrument for CCAT-prime. Prime-Cam will house seven instrument modules in a 1.8 meter diameter cryostat, cooled by a dilution refrigerator. The optical elements will consist of silicon lenses, and the instrument modules can be individually optimized for particular science goals. The current design enables both broadband, dual-polarization measurements and narrow-band, Fabry-Perot spectroscopic imaging using multichroic transition-edge sensor (TES) bolometers operating between 190 and 450 GHz. It also includes broadband kinetic induction detectors (KIDs) operating at 860 GHz. This wide range of frequencies will allow excellent characterization and removal of galactic foregrounds, which will enable precision measurements of the sub-mm and mm sky. Prime-Cam will be used to constrain cosmology via the Sunyaev-Zeldovich effects, map the intensity of [CII] 158 $\mu$m emission from the Epoch of Reionization, measure Cosmic Microwave Background polarization and foregrounds, and characterize the star formation history over a wide range of redshifts. More information about CCAT-prime can be found at www.ccatobservatory.org.

The initial conditions of cosmological simulations are commonly drawn from a Gaussian ensemble. The limited number of modes inside a simulation volume gives rise to statistical fluctuations known as \textit{sample variance}, limiting the accuracy of simulation predictions. Fixed fields offer an alternative initialization strategy; they have the same power spectrum as standard Gaussian fields but without intrinsic amplitude scatter at linear order. Paired fixed fields consists of two fixed fields with opposite phases that cancel phase correlations which otherwise induce second-order scatter in the non-linear power spectrum. We study the statistical properties of those fields for 19 different quantities at different redshifts through a large set of 600 N-body and 506 state-of-the-art magneto-hydrodynamic simulations covering a wide range of scales, mass and spatial resolutions. We find that paired fixed simulations do not introduce a bias on any of the examined quantities. We quantify the statistical improvement brought by these simulations, over standard ones, on different power spectra such as matter, halos, CDM, gas, stars, black-holes and magnetic fields, finding that they can reduce their variance by factors as large as $10^6$. We quantify the improvement achieved by fixing and by pairing, showing that sample variance in some quantities can be highly suppressed by pairing after fixing. Paired fixed simulations do not change the scatter in quantities such as the probability distribution function of matter density, or the halo, void or stellar mass functions. We argue that procedures aiming at reducing the sample variance of those quantities are unlikely to work. Our results show that paired fixed simulations do not affect either mean relations or scatter of galaxy properties, and suggest that the information embedded in 1-pt statistics is highly complementary to that in clustering.

Upcoming cosmic microwave background (CMB) surveys will soon make the first detection of the polarized Sunyaev-Zel'dovich effect, the linear polarization generated by the scattering of CMB photons on the free electrons present in collapsed objects. Measurement of this polarization along with knowledge of the electron density of the objects allows a determination of the quadrupolar temperature anisotropy of the CMB as viewed from the space-time location of the objects. Maps of these remote temperature quadrupoles have several cosmological applications. Here we propose a new application: reconstruction of the cosmological reionization history. We show that with quadrupole measurements out to redshift 3, constraints on the mean optical depth can be improved by an order of magnitude beyond the CMB cosmic variance limit.

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

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present a catalog of 182 galaxy clusters detected through the Sunyaev-Zel'dovich effect by the Atacama Cosmology Telescope in a contiguous 987.5 deg$^{2}$ field. The clusters were detected as SZ decrements by applying a matched filter to 148 GHz maps that combine the original ACT equatorial survey with data from the first two observing seasons using the ACTPol receiver. Optical/IR confirmation and redshift measurements come from a combination of large public surveys and our own follow-up observations. Where necessary, we measured photometric redshifts for clusters using a pipeline that achieves accuracy $\Delta z/(1 + z)=0.015$ when tested on SDSS data. Under the assumption that clusters can be described by the so-called Universal Pressure Profile and its associated mass-scaling law, the full signal-to-noise > 4 sample spans the mass range $1.6 < M^{\rm UPP}_{\rm 500c}/10^{14}{\rm M}_{\odot}<9.1$, with median $M^{\rm UPP}_{\rm 500c}=3.1 \times 10^{14}$ M$_{\odot}$. The sample covers the redshift range $0.1 < z < 1.4$ (median $z = 0.49$) and 28 clusters are new discoveries (median $z = 0.80$). We compare our catalog with other overlapping cluster samples selected using the SZ, optical,and X-ray wavelengths. We find the ratio of the UPP-based SZ mass to richness-based weak-lensing mass is $\langle M^{\rm UPP}_{\rm 500c} \rangle / \langle M^{\rm \lambda WL}_{\rm 500c} \rangle = 0.68 \pm 0.11$. After applying this calibration, the mass distribution for clusters with $M_{\rm 500c} > 4 \times 10^{14}$ M$_{\odot}$ is consistent with the number of such clusters found in the South Pole Telescope SZ survey.

Future high-resolution measurements of the cosmic microwave background (CMB) will produce catalogs of tens of thousands of galaxy clusters through the thermal Sunyaev-Zel'dovich (tSZ) effect. We forecast how well different configurations of a CMB Stage-4 experiment can constrain cosmological parameters, in particular the amplitude of structure as a function of redshift $\sigma_8(z)$, the sum of neutrino masses $\Sigma m_{\nu}$, and the dark energy equation of state $w(z)$. A key element of this effort is calibrating the tSZ scaling relation by measuring the lensing signal around clusters. We examine how the mass calibration from future optical surveys like the Large Synoptic Survey (LSST) compares with a purely internal calibration using lensing of the CMB itself. We find that, due to its high-redshift leverage, internal calibration gives constraints on cosmological parameters comparable to the optical calibration, and can be used as a cross-check of systematics in the optical measurement. We also show that in contrast to the constraints using the CMB lensing power spectrum, lensing-calibrated tSZ cluster counts can detect a minimal $\Sigma m_{\nu}$ at the 3-5$\sigma$ level even when the dark energy equation of state is freed up.

The improving sensitivity of measurements of the kinetic Sunyaev-Zel'dovich (SZ) effect opens a new window into the thermodynamic properties of the baryons in halos. We propose a methodology to constrain these thermodynamic properties by combining the kinetic SZ, which is an unbiased probe of the free electron density, and the thermal SZ, which probes their thermal pressure. We forecast that our method constrains the average thermodynamic processes that govern the energetics of galaxy evolution like energetic feedback across all redshift ranges where viable halos sample are available. Current Stage-3 cosmic microwave background (CMB) experiments like AdvACT and SPT-3G can measure the kSZ and tSZ to greater than 100$\sigma$ if combined with a DESI-like spectroscopic survey. Such measurements translate into percent-level constraints on the baryonic density and pressure profiles and on the feedback and non-thermal pressure support parameters for a given ICM model. This in turn will provide critical thermodynamic tests for sub-grid models of feedback in cosmological simulations of galaxy formation. The high fidelity measurements promised by the next generation CMB experiment, CMB-S4, allow one to further sub-divide these constraints beyond redshift into other classifications, like stellar mass or galaxy type.

The use of large, X-ray selected galaxy cluster catalogues for cosmological analyses requires a thorough understanding of the X-ray mass estimates. Weak gravitational lensing is an ideal method to shed light on such issues, due to its insensitivity to the cluster dynamical state. We perform a weak lensing calibration of 166 galaxy clusters from the REFLEX and BCS cluster catalogue and compare our results to the X-ray masses based on scaled luminosities from that catalogue. To interpret the weak lensing signal in terms of cluster masses, we compare the lensing signal to simple theoretical Navarro-Frenk-White models and to simulated cluster lensing profiles, including complications such as cluster substructure, projected large-scale structure, and Eddington bias. We find evidence of underestimation in the X-ray masses, as expected, with $\langle M_{\mathrm{X}}/M_{\mathrm{WL}}\rangle = 0.75 \pm 0.07$ stat. $\pm 0.05$ sys. for our best-fit model. The biases in cosmological parameters in a typical cluster abundance measurement that ignores this mass bias will typically exceed the statistical errors.

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

There are well-observed differences between cool-core (CC) and non-cool-core (NCC) clusters, but the origin of this distinction is still largely unknown. Competing theories can be divided into internal (inside-out), in which internal physical processes transform or maintain the NCC phase, and external (outside-in), in which the cluster type is determined by its initial conditions, which in turn lead to different formation histories (i.e., assembly bias). We propose a new method that uses the relative assembly bias of CC to NCC clusters, as determined via the two-point cluster-galaxy cross-correlation function (CCF), to test whether formation history plays a role in determining their nature. We apply our method to 48 ACCEPT clusters, which have well resolved central entropies, and cross-correlate with the SDSS-III/BOSS LOWZ galaxy catalog. We find that the relative bias of NCC over CC clusters is $b = 1.42 \pm 0.35$ ($1.6\sigma$ different from unity). Our measurement is limited by the small number of clusters with core entropy information within the BOSS footprint, 14 CC and 34 NCC. Future compilations of X-ray cluster samples, combined with deep all-sky redshift surveys, will be able to better constrain the relative assembly bias of CC and NCC clusters and determine the origin of the bimodality.

The recent emergence of detections of the kinetic Sunyaev-Zel'dovich (kSZ) effect through cross-correlation techniques is encouraging for the prospects of future cosmic microwave background (CMB) experiments. Extracting information on the large-scale velocity fields and constraining cosmological parameters from such kSZ measurements requires an understanding of the optical depth to CMB photons through halos. Using cosmological hydrodynamic simulations we find that there exists a low-scatter relation between the optical depth and thermal Sunyaev-Zel'dovich (tSZ) signal of halos within a physical aperture. We propose that such a relation can be used to break the degeneracy between optical depth and line-of-sight velocity in kSZ measurements. The limiting factors in our proposal are systematic uncertainties associated with the sub-grid physics models in the simulations, which we calculate to be less than 10 percent. We discuss future observational measurements that could potentially be used to mitigate the systematic uncertainties in this scaling relation.

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

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

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

The kinematic Sunyaev-Zel'dovich (kSZ) effect --- the Doppler boosting of cosmic microwave background (CMB) photons due to Compton-scattering off free electrons with non-zero bulk velocity --- probes the abundance and distribution of baryons in the Universe. All kSZ measurements to date have explicitly required spectroscopic redshifts. Here, we implement a novel estimator for the kSZ -- large-scale structure cross-correlation based on projected fields: it does not require redshift estimates for individual objects, allowing kSZ measurements from large-scale imaging surveys. We apply this estimator to cleaned CMB temperature maps constructed from Planck and Wilkinson Microwave Anisotropy Probe data and a galaxy sample from the Wide-field Infrared Survey Explorer (WISE). We measure the kSZ effect at 3.8-4.5$\sigma$ significance, depending on the use of additional WISE galaxy bias constraints. We verify that our measurements are robust to possible dust emission from the WISE galaxies. Assuming the standard $\Lambda$CDM cosmology, we directly constrain $( {f_{b}}/{0.158} ) ( {f_{\rm free}}/{1.0} ) = 1.48 \pm 0.19$ (statistical error only) at redshift $z \approx 0.4$, where $f_{b}$ is the fraction of matter in baryonic form and $f_{\rm free}$ is the free electron fraction. This is the tightest kSZ-derived constraint reported to date on these parameters. The consistency between the $f_{b}$ value found here and the values inferred from analyses of the primordial CMB and Big Bang nucleosynthesis verifies that baryons approximately trace the dark matter distribution down to $\sim$Mpc scales. While our projected-field estimator is already competitive with other kSZ approaches when applied to current datasets (because we are able to use the full-sky WISE photometric survey), it will yield enormous signal-to-noise when applied to upcoming high-resolution, multi-frequency CMB surveys.

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

We show that the gravitational wave source counts distribution can test how gravitational radiation propagates on cosmological scales. This test does not require obtaining redshifts for the sources. If the signal-to-noise ratio (SNR, $\rho$) from a gravitational wave source is proportional to the strain then it falls as $R^{-1}$, thus we expect the source counts to follow $dN/d\rho \propto \rho^{-4}$. However, if gravitational waves decay as they propagate or propagate into other dimensions, then there can be deviations from this generic prediction. We consider the possibility that the strain falls as $R^{-\gamma}$, where $\gamma=1$ recovers the expected predictions in a Euclidean uniformly-filled universe, and forecast the sensitivity of future observations to deviations from standard General Relativity. We first consider the case of few objects, 7 sources, with a signal-to-noise from 8 to 24, and impose a lower limit on $\gamma$, finding $\gamma>0.33$ at $95\%$ confidence level. The distribution of our simulated sample is very consistent with the distribution of the trigger events reported by Advanced LIGO. Future measurements will improve these constraints: with 100 events, we estimate that $\gamma$ can be measured with an uncertainty of $15\%$. We generalize the formalism to account for a range of chirp masses and the possibility that the signal falls as $\exp(-R/R_0)/R^\gamma$.

SPHEREx is a proposed SMEX mission selected for Phase A. SPHEREx will carry out the first all-sky spectral survey and provide for every 6.2" pixel a spectra between 0.75 and 4.18 $\mu$m [with R$\sim$41.4] and 4.18 and 5.00 $\mu$m [with R$\sim$135]. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. It is readily apparent, however, that many other questions in astrophysics and planetary sciences could be addressed with the SPHEREx data. The SPHEREx team convened a community workshop in February 2016, with the intent of enlisting the aid of a larger group of scientists in defining these questions. This paper summarizes the rich and varied menu of investigations that was laid out. It includes studies of the composition of main belt and Trojan/Greek asteroids; mapping the zodiacal light with unprecedented spatial and spectral resolution; identifying and studying very low-metallicity stars; improving stellar parameters in order to better characterize transiting exoplanets; studying aliphatic and aromatic carbon-bearing molecules in the interstellar medium; mapping star formation rates in nearby galaxies; determining the redshift of clusters of galaxies; identifying high redshift quasars over the full sky; and providing a NIR spectrum for most eROSITA X-ray sources. All of these investigations, and others not listed here, can be carried out with the nominal all-sky spectra to be produced by SPHEREx. In addition, the workshop defined enhanced data products and user tools which would facilitate some of these scientific studies. Finally, the workshop noted the high degrees of synergy between SPHEREx and a number of other current or forthcoming programs, including JWST, WFIRST, Euclid, GAIA, K2/Kepler, TESS, eROSITA and LSST.

We present galaxy velocity dispersions and dynamical mass estimates for 44 galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect by the Atacama Cosmology Telescope. Dynamical masses for 18 clusters are reported here for the first time. Using \Nbody\ simulations, we model the different observing strategies used to measure the velocity dispersions and account for systematic effects resulting from these strategies. We find that the galaxy velocity distributions may be treated as isotropic, and that an aperture correction of up to 7 per cent in the velocity dispersion is required if the spectroscopic galaxy sample is sufficiently concentrated towards the cluster centre. Accounting for the radial profile of the velocity dispersion in simulations enables consistent dynamical mass estimates regardless of the observing strategy. Cluster masses $M_{200}$ are in the range $(1-15)\times10^{14}M_\odot$. Comparing with masses estimated from the SZ distortion assuming a gas pressure profile derived from X-ray observations gives a mean SZ-to-dynamical mass ratio of $1.10\pm0.13$, but there is an additional 0.14 systematic uncertainty due to the unknown velocity bias; the statistical uncertainty is dominated by the scatter in the mass-velocity dispersion scaling relation. This ratio is consistent with previous determinations at these mass scales.

The kinematic Sunyaev-Zel'dovich (kSZ) signal is a powerful probe of the cosmic baryon distribution. The kSZ signal is proportional to the integrated free electron momentum rather than the electron pressure (which sources the thermal SZ signal). Since velocities should be unbiased on large scales, the kSZ signal is an unbiased tracer of the large-scale electron distribution, and thus can be used to detect the "missing baryon" that evade most observational techniques. While most current methods for kSZ extraction rely on the availability of very accurate redshifts, we revisit a method that allows measurements even in the absence of redshift information for individual objects. It involves cross-correlating the square of an appropriately filtered cosmic microwave background (CMB) temperature map with a projected density map constructed from a sample of large-scale structure tracers. We show that this method will achieve high signal-to-noise when applied to the next generation of high-resolution CMB experiments, provided that component separation is sufficiently effective at removing foreground contamination. Considering statistical errors only, we forecast that this estimator can yield $S/N \approx $ 3, 120 and over 150 for Planck, Advanced ACTPol, and hypothetical Stage-IV CMB experiments, respectively, in combination with a galaxy catalog from WISE, and about 20% larger $S/N$ for a galaxy catalog from the proposed SPHEREx experiment. This work serves as a companion paper to the first kSZ measurement with this method, where we used CMB temperature maps constructed from Planck and WMAP data, together with galaxies from the WISE survey, to obtain a 3.8 - 4.5$\sigma$ detection of the kSZ$^2$ amplitude.

Using a radio-quiet subsample of the Sloan Digital Sky Survey spectroscopic quasar catalogue, spanning redshifts 0.5-3.5, we derive the mean millimetre and far-infrared quasar spectral energy distributions (SEDs) via a stacking analysis of Atacama Cosmology Telescope and Herschel-Spectral and Photometric Imaging REceiver data. We constrain the form of the far-infrared emission and find 3$\sigma$-4$\sigma$ evidence for the thermal Sunyaev-Zel'dovich (SZ) effect, characteristic of a hot ionized gas component with thermal energy $(6.2 \pm 1.7)\times 10^{60}$ erg. This amount of thermal energy is greater than expected assuming only hot gas in virial equilibrium with the dark matter haloes of $(1-5)\times 10^{12}h^{-1}$M$_\odot$ that these systems are expected to occupy, though the highest quasar mass estimates found in the literature could explain a large fraction of this energy. Our measurements are consistent with quasars depositing up to $(14.5 \pm 3.3)~\tau_8^{-1}$ per cent of their radiative energy into their circumgalactic environment if their typical period of quasar activity is $\tau_8\times~10^8$ yr. For high quasar host masses, $\sim10^{13}h^{-1}$M$_\odot$, this percentage will be reduced. Furthermore, the uncertainty on this percentage is only statistical and additional systematic uncertainties enter at the 40 per cent level. The SEDs are dust dominated in all bands and we consider various models for dust emission. While sufficiently complex dust models can obviate the SZ effect, the SZ interpretation remains favoured at the 3$\sigma$-4$\sigma$ level for most models.

The cosmic microwave background (CMB) energy spectrum is a near-perfect blackbody. The standard model of cosmology predicts small spectral distortions to this form, but no such distortion of the sky-averaged CMB spectrum has yet been measured. We calculate the largest expected distortion, which arises from the inverse Compton scattering of CMB photons off hot, free electrons, known as the thermal Sunyaev-Zel'dovich (tSZ) effect. We show that the predicted signal is roughly one order of magnitude below the current bound from the COBE-FIRAS experiment, but can be detected at enormous significance ($\gtrsim 1000\sigma$) by the proposed Primordial Inflation Explorer (PIXIE). Although cosmic variance reduces the effective signal-to-noise to $230\sigma$, this measurement will still yield a sub-percent constraint on the total thermal energy of electrons in the observable universe. Furthermore, we show that PIXIE can detect subtle relativistic effects in the sky-averaged tSZ signal at $30\sigma$, which directly probe moments of the optical depth-weighted intracluster medium electron temperature distribution. These effects break the degeneracy between the electron density and temperature in the mean tSZ signal, allowing a direct inference of the mean baryon density at low redshift. Future spectral distortion probes will thus determine the global thermodynamic properties of ionized gas in the universe with unprecedented precision. These measurements will impose a fundamental "integral constraint" on models of galaxy formation and the injection of feedback energy over cosmic time.

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

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

We use the Planck full mission temperature maps to examine the stacked thermal Sunyaev-Zel'dovich (tSZ) signal of 188042 "locally brightest galaxies'" (LBGs) selected from the Sloan Digital Sky Survey Data Release 7. Our LBG sample closely matches that of Planck Collaboration XI (2013, PCXI) but our analysis differs in several ways. We work directly in terms of physically observable quantities, requiring minimal assumptions about the gas pressure profile. We explicitly model the dust emission from each LBG and simultaneously measure both the stacked tSZ and dust signals as a function of stellar mass $M_*$. There is a small residual bias in stacked tSZ measurements; we measure this bias and subtract it from our results, finding that the effects are non-negligible at the lowest masses in the LBG sample. Finally, we compare our measurements with two pressure profile models, finding that the profile from Battaglia et al. (2012b) provides a better fit to the results than the Arnaud et al. (2010) "universal pressure profile." However, within the uncertainties, we find that the data are consistent with a self-similar scaling with mass --- more precise measurements are needed to detect the relatively small deviations from self-similarity predicted by these models. Consistent with PCXI, we measure the stacked tSZ signal from LBGs with stellar masses down to $\log_{10}(M_*/M_{\odot}) \sim 11.1-11.3$. For lower stellar masses, however, we do not see evidence for a stacked tSZ signal. We note that the stacked dust emission is comparable to, or larger than, the stacked tSZ signal for $\log_{10}(M_*/M_{\odot}) \lesssim 11.3$. Future tSZ analyses with larger samples and lower noise levels should be able to probe deviations from self-similarity and thus provide constraints on models of feedback and the evolution of hot halo gas over cosmic time.

Whether among the myriad tiny proto-galaxies there exists a population with similarities to present day galaxies is an open question. We show, using BlueTides, the first hydrodynamic simulation large enough to resolve the relevant scales, that the first massive galaxies to form are %in fact predicted to have extensive rotationally-supported disks. Although their morphology resembles in some ways Milky-way types seen at much lower redshifts, these high-redshift galaxies are smaller, denser, and richer in gas than their low redshift counterparts. From a kinematic analysis of a statistical sample of 216 galaxies at redshift $z=8-10$ we have found that disk galaxies make up 70\% of the population of galaxies with stellar mass $10^{10} M_\odot$ or greater. Cold Dark Matter cosmology therefore makes specific predictions for the population of large galaxies 500 million years after the Big Bang. We argue that wide-field satellite telescopes (e.g. WFIRST) will in the near future discover these first massive disk galaxies. The simplicity of their structure and formation history should make possible new tests of cosmology.

We present a modern machine learning approach for cluster dynamical mass measurements that is a factor of two improvement over using a conventional scaling relation. Different methods are tested against a mock cluster catalog constructed using halos with mass >= 10^14 Msolar/h from Multidark's publicly-available N-body MDPL halo catalog. In the conventional method, we use a standard M(sigma_v) power law scaling relation to infer cluster mass, M, from line-of-sight (LOS) galaxy velocity dispersion, sigma_v. The resulting fractional mass error distribution is broad, with width=0.87 (68% scatter), and has extended high-error tails. The standard scaling relation can be simply enhanced by including higher-order moments of the LOS velocity distribution. Applying the kurtosis as a correction term to log(sigma_v) reduces the width of the error distribution to 0.74 (16% improvement). Machine learning can be used to take full advantage of all the information in the velocity distribution. We employ the Support Distribution Machines (SDMs) algorithm that learns from distributions of data to predict single values. SDMs trained and tested on the distribution of LOS velocities yield width=0.46 (47% improvement). Furthermore, the problematic tails of the mass error distribution are effectively eliminated. Decreasing cluster mass errors will improve measurements of the growth of structure and lead to tighter constraints on cosmological parameters.

We introduce the BlueTides simulation and report initial results for the luminosity functions of the first galaxies and AGN, and their contribution to reionization. BlueTides was run on the BlueWaters cluster at NCSA from $z=99$ to $z=8.0$ and includes 2$\times$7040$^3$ particles in a $400$Mpc/h per side box, making it the largest hydrodynamic simulation ever performed at high redshift. BlueTides includes a pressure-entropy formulation of smoothed particle hydrodynamics, gas cooling, star formation (including molecular hydrogen), black hole growth and models for stellar and AGN feedback processes. The star formation rate density in the simulation is a good match to current observational data at $z\sim 8-10$. We find good agreement between observations and the predicted galaxy luminosity function in the currently observable range $-18\le M_{\mathrm UV} \le -22.5$ with some dust extinction required to match the abundance of brighter objects. BlueTides implements a patchy reionization model that produces a fluctuating UV background. BlueTides predicts number counts for galaxies fainter than current observational limits which are consistent with extrapolating the faint end slope of the luminosity function with a power law index $\alpha\sim -1.8$ at $z\sim 8$ and redshift dependence of $\alpha\sim (1+z)^{-0.4}$. The AGN population has a luminosity function well fit by a power law with a slope $\alpha\sim -2.4$ that compares favourably with the deepest CANDELS-Goods fields. We investigate how these luminosity functions affect the progress of reionization, and find that a high Lyman-$\alpha$ escape fraction ($f_\mathrm{esc} \sim 0.5$) is required if galaxies dominate the ionising photon budget during reionization. Smaller galaxy escape fractions imply a large contribution from faint AGN (down to $M_\mathrm{UV}=-12$) which results in a rapid reionization, disfavoured by current observations.

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

We present Southern African Large Telescope (SALT) follow-up observations of seven massive clusters detected by the Atacama Cosmology Telescope (ACT) on the celestial equator using the Sunyaev-Zel'dovich (SZ) effect. We conducted multi-object spectroscopic observations with the Robert Stobie Spectrograph in order to measure galaxy redshifts in each cluster field, determine the cluster line-of-sight velocity dispersions, and infer the cluster dynamical masses. We find that the clusters, which span the redshift range 0.3 < z < 0.55, range in mass from (5 -- 20) x 10$^{14}$ solar masses (M200c). Their masses, given their SZ signals, are similar to those of southern hemisphere ACT clusters previously observed using Gemini and the VLT. We note that the brightest cluster galaxy in one of the systems studied, ACT-CL J0320.4+0032 at z = 0.38, hosts a Type II quasar. Only a handful of such systems are currently known, and therefore ACT-CL J0320.4+0032 may be a rare example of a very massive halo in which quasar-mode feedback is actively taking place.

Recent first detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) signal in Planck cosmic microwave background (CMB) temperature maps with gravitational lensing maps inferred from the Planck CMB data and the CFHTLenS galaxy survey provide new probes of the relationship between baryons and dark matter. Using cosmological hydrodynamics simulations, we show that these cross-correlation signals are dominated by contributions from hot gas in the intracluster medium (ICM), rather than diffuse, unbound gas located beyond the virial radius (the "missing baryons"). Thus, these cross-correlations offer a tool with which to study the ICM over a wide range of halo masses and redshifts. In particular, we show that the tSZ - CMB lensing cross-correlation is more sensitive to gas in lower-mass, higher-redshift halos and gas at larger cluster-centric radii than the tSZ - galaxy lensing cross-correlation. Combining these measurements with primary CMB data will constrain feedback models through their signatures in the ICM pressure profile. We forecast the ability of ongoing and future experiments to constrain such ICM parameters, including the mean amplitude of the pressure - mass relation, the redshift evolution of this amplitude, and the mean outer logarithmic slope of the pressure profile. The results are promising, with $\approx 5-20$% precision constraints achievable with upcoming experiments, even after marginalizing over cosmological parameters.

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

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

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

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

We examine the importance of baryonic feedback effects on the matter power spectrum on small scales, and the implications for the precise measurement of neutrino masses through gravitational weak lensing. Planned large galaxy surveys such as the Large Synoptic Sky Telescope (LSST) and Euclid are expected to measure the sum of neutrino masses to extremely high precision, sufficient to detect non-zero neutrino masses even in the minimal mass normal hierarchy. We show that weak lensing of galaxies while being a very good probe of neutrino masses, is extremely sensitive to baryonic feedback processes. We use publicly available results from the Overwhelmingly Large Simulations (OWLS) project to investigate the effects of active galactic nuclei feedback, the nature of the stellar initial mass function, and gas cooling rates, on the measured weak lensing shear power spectrum. Using the Fisher matrix formalism and priors from CMB+BAO data, we show that when one does not account for feedback, the measured neutrino mass may be substantially larger or smaller than the true mass, depending on the dominant feedback mechanism, with the mass error |\Delta m_nu| often exceeding the mass m_nu itself. We also consider gravitational lensing of the cosmic microwave background (CMB) and show that it is not sensitive to baryonic feedback on scales l < 2000, although CMB experiments that aim for sensitivities sigma(m_nu) < 0.02 eV will need to include baryonic effects in modeling the CMB lensing potential. A combination of CMB lensing and galaxy lensing can help break the degeneracy between neutrino masses and baryonic feedback processes. We conclude that future large galaxy lensing surveys such as LSST and Euclid can only measure neutrino masses accurately if the matter power spectrum can be measured to similar accuracy.

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.

We present predictions for the 21 cm brightness temperature power spectrum during the Epoch of Reionization (EoR). We discuss the implications of the "light cone" effect, which incorporates evolution of the neutral hydrogen fraction and 21 cm brightness temperature along the line of sight. Using a novel method calibrated against radiation-hydrodynamic simulations, we model the neutral hydrogen density field and 21 cm signal in large volumes ($L = 2$ Gpc/$h$). The inclusion of the light cone effect leads to a relative decrease of about 50% in the 21 cm power spectrum on all scales. We also find that the effect is more prominent at the midpoint of reionization and later. The light cone effect also can introduce an anisotropy along the line of sight. By decomposing the 3D power spectrum into components perpendicular to and along the line of sight, we find that in our fiducial reionization model, there is no significant anisotropy. However, parallel modes can contribute up to 40% more power for shorter reionization scenarios. The scales on which the light cone effect is relevant are comparable to scales where one measures the baryon acoustic oscillation. We argue that due to its large comoving scale and introduction of anisotropy, the light cone effect is important when considering redshift space distortions and future application to the Alcock-Paczynski test for the determination of cosmological parameters.

Understanding the outskirts of galaxy clusters at the virial radius (R200) and beyond is critical for an accurate determination of cluster masses and to ensure unbiased cosmological parameter estimates from cluster surveys. This problem has drawn renewed interest due to recent determinations of gas mass fractions beyond R200, which appear to be considerably larger than the cosmic mean, and because the clusters' total Sunyaev-Zel'dovich flux receives a significant contribution from these regions. Here, we use a large suite of cosmological hydrodynamical simulations to study the clumpiness of density and pressure and employ different variants of simulated physics, including radiative gas physics and thermal feedback by active galactic nuclei. We find that density and pressure clumping closely trace each other as a function of radius, but the bias on density remains on average < 20% within the virial radius R200. At larger radius, clumping increases steeply due to the continuous infall of coherent structures that have not yet passed the accretion shock. Density and pressure clumping increase with cluster mass and redshift, which probes on average dynamically younger objects that are still in the process of assembling. The angular power spectra of gas density and pressure show that the clumping signal is dominated by comparably large substructures with scales >R200/5, signaling the presence of gravitationally-driven "super-clumping". In contrast, the angular power spectrum of the dark matter (DM) shows an almost uniform size distribution due to unimpeded subhalos. The quadrupolar anisotropy dominates the signal and correlates well across different radii as a result of the prolateness of the DM potential. We provide a synopsis of the radial dependence of the clusters' non-equilibrium measures (kinetic pressure support, ellipticity, and clumping) that all increase sharply beyond R200.

We present 610 MHz and 2.1 GHz imaging of the massive SZE-selected z=0.870 cluster merger ACT-CL J0102-4915 (El Gordo), obtained with the GMRT and the ATCA, respectively. We detect two complexes of radio relics separated by 3.4' (1.6 Mpc) along the system's NW-to-SE collision axis that have high integrated polarizations (33%) and steep spectral indices, consistent with creation via Fermi acceleration by shocks in the ICM. From the spectral index of the relics, we compute a Mach number of 2.5^{+0.7}_{-0.3} and shock speed of 2500^{+400}_{-300} km/s. With our ATCA data, we compute the Faraday depth across the NW relic and find a mean value of 11 rad/m^2 and standard deviation of 6 rad/m^2. With the integrated line-of-sight gas density derived from new Chandra observations, our Faraday depth measurement implies B_parallel~0.01 \mu G in the cluster outskirts. The extremely narrow shock widths in the relics (<23 kpc) prevent us from placing a meaningful constraint on |B| using cooling time arguments. In addition to the relics, we detect a large (1.1 Mpc radius), powerful (log L_1.4[W/Hz]= 25.66+-0.12) radio halo with a Bullet-like morphology. The spectral-index map of the halo shows the synchrotron spectrum is flattest near the relics, along the collision axis, and in regions of high T_gas, all locations associated with recent energy injection. The spatial and spectral correlation between the halo emission and cluster X-ray properties supports primary-electron processes like turbulent reacceleration as the halo production mechanism. The halo's integrated 610 MHz to 2.1 GHz spectral index is 1.2+-0.1, consistent with the cluster's high T_gas in view of previously established global scaling relations. El Gordo is the highest-redshift cluster known to host a radio halo and/or radio relics, and provides new constraints on the non-thermal physics in clusters at z>0.6. [abridged]

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

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

We investigate the effect of patchy reionization on the cosmic microwave background temperature. An anisotropic optical depth tau (theta) alters the TT power spectrum on small scales l > 2000. We make use of the correlation between the matter density and the reionization redshift fields to construct full sky maps of tau(theta). Patchy reionization transfers CMB power from large scales to small scales, resulting in a non-zero cross correlation between large and small angular scales. We show that the patchy tau correlator is sensitive to small root mean square values tau_rms ~ 0.003 seen in our maps. We include other secondary anisotropies such as CMB lensing, kinetic and thermal Sunyaev-Zel'dovich terms, as well as the infrared and point source background, and show that patchy reionization may be detected in the low frequency channels ~ 90 GHz, particularly for extended reionization histories. If frequency dependent secondaries can be minimized by a multi-frequency analysis, we show that even small degrees of patchiness can be detected at high significance. In this case, we show that models with different values of tau_rms are distinguishable, even for the same mean value <tau>. The patchy tau correlator thus provides information that is complementary to what may be obtained from the polarization and the kinetic Sunyaev-Zel'dovich power spectra.

We measure a significant correlation between the thermal Sunyaev-Zel'dovich effect in the Planck and WMAP maps and an X-ray cluster map based on ROSAT. We use the 100, 143 and 343 GHz Planck maps and the WMAP 94 GHz map to obtain this cluster cross spectrum. We check our measurements for contamination from dusty galaxies using the cross correlations with the 220, 545 and 843 GHz maps from Planck. Our measurement yields a direct characterization of the cluster power spectrum over a wide range of angular scales that is consistent with large cosmological simulations. The amplitude of this signal depends on cosmological parameters that determine the growth of structure (\sigma_8 and \Omega_M) and scales as \sigma_8^7.4 and \Omega_M^1.9 around the multipole (ell) ~ 1000. We constrain \sigma_8 and \Omega_M from the cross-power spectrum to be \sigma_8 (\Omega_M/0.30)^0.26 = 0.8 +/- 0.02. Since this cross spectrum produces a tight constraint in the \sigma_8 and \Omega_M plane the errors on a \sigma_8 constraint will be mostly limited by the uncertainties from external constraints. Future cluster catalogs, like those from eRosita and LSST, and pointed multi-wavelength observations of clusters will improve the constraining power of this cross spectrum measurement. In principle this analysis can be extended beyond \sigma_8 and \Omega_M to constrain dark energy or the sum of the neutrino masses.

Gas masses tightly correlate with the virial masses of galaxy clusters, allowing for a precise determination of cosmological parameters by means of large-scale X-ray surveys. However, according to recent Suzaku X-ray measurements, gas mass fractions, f_gas, appear to be considerably larger than the cosmic mean at the virial radius, R_200, questioning the accuracy of the cosmological parameter estimations. Here, we use a large suite of cosmological hydrodynamical simulations to study measurement biases of f_gas. We employ different variants of simulated physics, including radiative gas physics, star formation, and thermal feedback by active galactic nuclei. Computing the mass profiles in 48 angular cones, whose footprints partition the sphere, we find anisotropic gas and total mass distributions that imply an angular variance of f_gas at the level of 30%. This anisotropic distribution originates from the recent formation epoch of clusters and from the strong internal baryon-to-dark-matter density bias. In the most extreme cones, f_gas can be biased high by a factor of two at R_200 in massive clusters, thereby providing a potential explanation for high f_gas measurements by Suzaku. While projection lowers this factor, there are other measurement biases that may (partially) compensate. We find that at R_200, f_gas is biased high by 20% when assuming hydrostatic equilibrium masses, i.e., neglecting the kinetic pressure, and by another ~10-20% due to the presence of density clumping. At larger radii, both measurement biases increase dramatically. While the cluster sample variance of the true f_gas decreases to a level of 5% at R_200, the sample variance that includes both measurement biases remains fairly constant at the level of 10-20%. The constant redshift evolution of f_gas within R_500 for massive clusters is encouraging for using gas masses to derive cosmological parameters.

We present a first measurement of the stellar mass component of galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect, using 3.6 um and 4.5 um photometry from the Spitzer Space Telescope. Our sample consists of 14 clusters detected by the Atacama Cosmology Telescope (ACT), which span the redshift range 0.27 < z < 1.07 (median z = 0.50), and have dynamical mass measurements, accurate to about 30 per cent, with median M500 = 6.9 x 10^{14} MSun. We measure the 3.6 um and 4.5 um galaxy luminosity functions, finding the characteristic magnitude (m*) and faint-end slope (alpha) to be similar to those for IR-selected cluster samples. We perform the first measurements of the scaling of SZ-observables (Y500 and y0) with both brightest cluster galaxy (BCG) stellar mass and total cluster stellar mass (M500star). We find a significant correlation between BCG stellar mass and Y500 (E(z)^{-2/3} DA^2 Y500 ~ M*^{1.2 +/- 0.6}), although we are not able to obtain a strong constraint on the slope of the relation due to the small sample size. Additionally, we obtain E(z)^{-2/3} DA^2 Y500 ~ M500star^{1.0 +/- 0.6} for the scaling with total stellar mass. The mass fraction in stars spans the range 0.006-0.034, with the second ranked cluster in terms of dynamical mass (ACT-CL J0237-4939) having an unusually low total stellar mass and the lowest stellar mass fraction. For the five clusters with gas mass measurements available in the literature, we see no evidence for a shortfall of baryons relative to the cosmic mean value.

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

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

We study primordial magnetic field effects on the matter perturbations in the Universe. We assume magnetic field generation prior to the big bang nucleosynthesis (BBN), i.e. during the radiation dominated epoch of the Universe expansion, but do not limit analysis by considering a particular magnetogenesis scenario. Contrary to previous studies, we limit the total magnetic field energy density and not the smoothed amplitude of the magnetic field at large (order of 1 Mpc) scales. We review several cosmological signatures, such as halos abundance, thermal Sunyaev Zel'dovich (tSZ) effect, and Lyman-$\alpha$ data. For a cross check we compare our limits with that obtained through the CMB faraday rotation effect and BBN. The limits are ranging between 1.5 nG and 4.5 nG for $n_B \in (-3; -1.5)$.

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

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

We present the 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 present new predictions for temperature (on small angular scales) and polarization (on large angular scales) CMB anisotropies induced during the epoch of reionization (EoR). Using a novel method calibrated from Radiation-Hydrodynamic simulations we model the EoR in large volumes (L >~ 2 Gpc/h) in the context of galactic reionization. We find that the EoR contribution to the kinetic Sunyaev- Zel'dovich power spectrum (patchy kSZ) ranges between ~0.6 - 2.8 muK^2 at ell = 3000, for the parameter space we explored. These patchy kSZ power spectra are calculated from large 15 Deg x 15 Deg maps that are found to be necessary. Decreasing the size of these maps biases the overall patchy kSZ power to higher values. We find that the amplitude of the patchy kSZ power spectrum at ell = 3000 follows simple scalings of D_ell=3000^kSZ propto <z> and D_ell=3000^kSZ propto Delz^0.47 for the mean redshift (<z>) of reionization and the duration (dz). Using the constraints on <z> from WMAP 7-year results and the lower limit on dz from EDGES we find a lower limit of ~ 0.4 muK^2 on the kSZ at ell = 3000. Planck will constrain the mean redshift and the Thomson optical depth from the low-ell polarization power spectrum. Future measurements of the high-ell CMB power spectrum from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) should detect the patchy kSZ signal if the cross correlation between the cosmic infrared background and the thermal Sunyaev Zel'dovich effect is constrained. We show that the combination of temperature and polarization measurements constrains both <z> and dz. The patchy kSZ maps, power spectra templates and the polarization power spectra will be publicly available.

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048^3 dark matter particles, 2048^3 gas cells, and 17 billion adaptive rays in a L = 100 Mpc/h box, we show that the density and reionization-redshift fields are highly correlated on large scales (>~ 1 Mpc/h). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric function for the bias, which is then used to filter any large-scale density field to derive the corresponding spatially varying reionization-redshift field. The parametric model has three free parameters which can be reduced to one free parameter when we fit the two bias parameters to simulations results. We can differentiate degenerate combinations of the bias parameters by combining results for the global ionization histories and correlation length between ionized regions. Unlike previous semi-analytic models, the evolution of the reionization-redshift field in our model is directly compared cell by cell against simulations and preforms well in all tests. Our model maps the high-resolution, intermediate-volume radiation-hydrodynamic simulations onto lower-resolution, larger-volume N-body simulations (>~ 2 Gpc/h) in order to make mock observations and theoretical predictions.

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

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

Secondary anisotropies in the cosmic microwave background are a treasure-trove of cosmological information. Interpreting current experiments probing them are limited by theoretical uncertainties rather than by measurement errors. Here we focus on the secondary anisotropies resulting from the thermal Sunyaev-Zel'dovich (tSZ) effect; the amplitude of which depends critically on the average thermal pressure profile of galaxy groups and clusters. To this end, we use a suite of SPH simulations that include radiative cooling, star formation, supernova feedback, and energetic feedback from active galactic nuclei (AGN). We examine in detail how the pressure profile depends on cluster radius, mass, and redshift and provide an empirical fitting function. We employ three different approaches for calculating the tSZ power spectrum: an analytical approach that uses our pressure profile fit, a semi-analytical method of pasting our pressure fit onto simulated clusters, and a direct numerical integration of our simulated volumes. We demonstrate that the detailed structure of the intracluster medium and cosmic web affect the tSZ power spectrum. In particular, the substructure and asphericity of clusters increase the tSZ power spectrum by 10-20% at ell ~2000-8000, with most of the additional power being contributed by substructures. The contributions to the power spectrum from radii larger than R_500 is ~ 20% at ell = 3000, thus clusters interiors (r < R_500) dominate the power spectrum amplitude at these angular scales.

The utility of large Sunyaev Zel'dovich (SZ) surveys for determining cosmological parameters from cluster abundances is limited by the theoretical uncertainties in the integrated SZ-flux-to-mass relation, Y-M. We explore how non-thermal pressure and the anisotropic shape of the gas distribution of the intracluster medium (ICM) impacts Y-M scaling using a suite of SPH simulations of the cosmic web. We contrast results for models with different treatments of entropy injection and transport, varying radiative cooling, star formation and accompanying supernova feedback, cosmic rays, and energetic feedback from active galactic nuclei (AGN). We find that the gas kinetic-to-thermal pressure ratio from internal bulk motions depends on the cluster mass, and increases in the outer-cluster due to enhanced substructure, as does the asphericity of the ICM gas. With only a ~5-10% correction to projected (observable) ellipticities, we can infer the 3D ellipticities. Our simulated Y-M-slope roughly follows the self-similar prediction, except for a steepening due to a deficit of gas in lower mass clusters at low redshift in our AGN-feedback simulations. AGN feedback enhances slightly the overall Y-M-scatter, from ~11% to ~13%, a reflection of accretion history variations due to cluster merging. If we split the cluster system into lower, middle and upper bands of both P_kin/P_th and long-to-short axis ratio, we find a ~10% effect on Y-M. Identifying observable second parameters related to internal bulk flows and anisotropy for cluster-selection to minimize Y-M scatter in a "fundamental plane" would allow tighter cosmological parameter constraints.

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.

We explore how radiative cooling, supernova feedback, cosmic rays and a new model of the energetic feedback from active galactic nuclei (AGN) affect the thermal and kinetic Sunyaev-Zel'dovich (SZ) power spectra. To do this, we use a suite of hydrodynamical TreePM-SPH simulations of the cosmic web in large periodic boxes and tailored higher resolution simulations of individual galaxy clusters. Our AGN feedback simulations match the recent universal pressure profile and cluster mass scaling relations of the REXCESS X-ray cluster sample better than previous analytical or numerical approaches. For multipoles $\ell\lesssim 2000$, our power spectra with and without enhanced feedback are similar, suggesting theoretical uncertainties over that range are relatively small, although current analytic and semi-analytic approaches overestimate this SZ power. We find the power at high 2000-10000 multipoles which ACT and SPT probe is sensitive to the feedback prescription, hence can constrain the theory of intracluster gas, in particular for the highly uncertain redshifts $>0.8$. The apparent tension between $\sigma_8$ from primary cosmic microwave background power and from analytic SZ spectra inferred using ACT and SPT data is lessened with our AGN feedback spectra.

Recent improvements in the capabilities of low frequency radio telescopes provide a unique opportunity to study thermal and non-thermal properties of the cosmic web. We argue that the diffuse, polarized emission from giant radio relics traces structure formation shock waves and illuminates the large-scale magnetic field. To show this, we model the population of shock-accelerated relativistic electrons in high-resolution cosmological simulations of galaxy clusters and calculate the resulting radio synchrotron emission. We find that individual shock waves correspond to localized peaks in the radio surface brightness map which enables us to measure Mach numbers for these shocks. We show that the luminosities and number counts of the relics strongly depend on the magnetic field properties, the cluster mass and dynamical state. By suitably combining different cluster data, including Faraday rotation measures, we are able to constrain some macroscopic parameters of the plasma at the structure formation shocks, such as models of turbulence. We also predict upper limits for the properties of the warm-hot intergalactic medium, such as its temperature and density. We predict that the current generation of radio telescopes (LOFAR, GMRT, MWA, LWA) have the potential to discover a substantially larger sample of radio relics, with multiple relics expected for each violently merging cluster. Future experiments (SKA) should enable us to further probe the macroscopic parameters of plasma physics in clusters.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.