Search

Upcoming advances in galaxy surveys and cosmic microwave background data will enable measurements of the anisotropic distribution of diffuse gas in filaments and superclusters at redshift $z=1$ and beyond, observed through the thermal Sunyaev-Zel'dovich (tSZ) effect. These measurements will help distinguish between different astrophysical feedback models, account for baryons that appear to be 'missing' from the cosmic census, and present opportunities for using locally-anisotropic tSZ statistics as cosmological probes. This study seeks to guide such future measurements by analysing whether diffuse intergalactic gas is a major contributor to anisotropic tSZ signal in The Three Hundred Gizmo-Simba hydrodynamic simulations. We apply multiple different halo boundary and temperature criteria to divide concentrated from diffuse gas at $z=1$, then create mock Compton-$y$ maps for the separated components. The maps from 98 simulation snapshots are centred on massive galaxy clusters, oriented by the most prominent galaxy filament axis, and stacked. Results vary significantly depending on the definition used for diffuse gas, indicating that assumptions should be clearly defined when claiming observations of the warm-hot intergalactic medium. In all cases, the diffuse gas is important, contributing 25-60% of the tSZ signal in the far field ($>4 h^{-1}$ comoving Mpc) from the stacked clusters. The gas 1-2 virial radii from halo centres is especially key. Oriented stacking and environmental selections help to amplify the signal from the warm-hot intergalactic medium, which is aligned but less concentrated along the filament axis than the hot halo gas., Published in MNRAS. 18 pages, 13 figures, 2 tables. Minor revisions for clarity; results unchanged

We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,\mu$K$^2$ greatly improves on the indirect upper limit of $420\,\mu$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,\mu$K$^2$, and the cosmic molecular gas density, $\rho_\text{H2}, Comment: Paper 5 of 7 in series. 17 pages + appendix and bibliography (30 pages total); 15 figures, 6 tables; accepted for publication in ApJ; v3 reflects the accepted version with minor changes and additions to text

The $\rho$ Oph molecular cloud is one of the best examples of spinning dust emission, first detected by the Cosmic Background Imager (CBI). Here we present 4.5 arcmin observations with CBI 2 that confirm 31 GHz emission from $\rho$ Oph W, the PDR exposed to B-type star HD 147889, and highlight the absence of signal from S1, the brightest IR nebula in the complex. In order to quantify an association with dust-related emission mechanisms, we calculated correlations at different angular resolutions between the 31 GHz map and proxies for the column density of IR emitters, dust radiance and optical depth templates. We found that the 31 GHz emission correlates best with the PAH column density tracers, while the correlation with the dust radiance improves when considering emission that is more extended (from the shorter baselines), suggesting that the angular resolution of the observations affects the correlation results. A proxy for the spinning dust emissivity reveals large variations within the complex, with a dynamic range of 25 at 3$\sigma$ and a variation by a factor of at least 23, at 3$\sigma$, between the peak in $\rho$ Oph W and the location of S1, which means that environmental factors are responsible for boosting spinning dust emissivities locally., Comment: 12 pages, 7 figures, 8 tables. Accepted for publication in MNRAS

We introduce COMAP-EoR, the next generation of the Carbon Monoxide Mapping Array Project aimed at extending CO intensity mapping to the Epoch of Reionization. COMAP-EoR supplements the existing 30 GHz COMAP Pathfinder with two additional 30 GHz instruments and a new 16 GHz receiver. This combination of frequencies will be able to simultaneously map CO(1--0) and CO(2--1) at reionization redshifts ($z\sim5-8$) in addition to providing a significant boost to the $z\sim3$ sensitivity of the Pathfinder. We examine a set of existing models of the EoR CO signal, and find power spectra spanning several orders of magnitude, highlighting our extreme ignorance about this period of cosmic history and the value of the COMAP-EoR measurement. We carry out the most detailed forecast to date of an intensity mapping cross-correlation, and find that five out of the six models we consider yield signal to noise ratios (S/N) $\gtrsim20$ for COMAP-EoR, with the brightest reaching a S/N above 400. We show that, for these models, COMAP-EoR can make a detailed measurement of the cosmic molecular gas history from $z\sim2-8$, as well as probe the population of faint, star-forming galaxies predicted by these models to be undetectable by traditional surveys. We show that, for the single model that does not predict numerous faint emitters, a COMAP-EoR-type measurement is required to rule out their existence. We briefly explore prospects for a third-generation Expanded Reionization Array (COMAP-ERA) capable of detecting the faintest models and characterizing the brightest signals in extreme detail., Comment: Paper 7 of 7 in series. 19 pages, 10 figures, to be submitted to ApJ

The CO Mapping Array Project (COMAP) aims to use line intensity mapping of carbon monoxide (CO) to trace the distribution and global properties of galaxies over cosmic time, back to the Epoch of Reionization (EoR). To validate the technologies and techniques needed for this goal, a Pathfinder instrument has been constructed and fielded. Sensitive to CO(1-0) emission from $z=2.4$-$3.4$ and a fainter contribution from CO(2-1) at $z=6$-8, the Pathfinder is surveying $12$ deg$^2$ in a 5-year observing campaign to detect the CO signal from $z\sim3$. Using data from the first 13 months of observing, we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\mathrm{K}^2 \mathrm{Mpc}^3$ on scales $k=0.051-0.62 \mathrm{Mpc}^{-1}$ - the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Based on these observations alone, we obtain a constraint on the amplitude of the clustering component (the squared mean CO line temperature-bias product) of $\langle Tb\rangle^2, Comment: Paper 1 of 7 in series. 18 pages, 16 figures, submitted to ApJ

Diverse astrophysical observations suggest the existence of cold dark matter that interacts only gravitationally with radiation and ordinary baryonic matter. Any nonzero coupling between dark matter and baryons would provide a significant step towards understanding the particle nature of dark matter. Measurements of the cosmic microwave background (CMB) provide constraints on such a coupling that complement laboratory searches. In this work we place upper limits on a variety of models for dark matter elastic scattering with protons and electrons by combining large-scale CMB data from the Planck satellite with small-scale information from Atacama Cosmology Telescope (ACT) DR4 data. In the case of velocity-independent scattering, we obtain bounds on the interaction cross section for protons that are 40\% tighter than previous constraints from the CMB anisotropy. For some models with velocity-dependent scattering we find best-fitting cross sections with a 2$\sigma$ deviation from zero, but these scattering models are not statistically preferred over $\Lambda$CDM in terms of model selection., Comment: 8 pages, 4 figures

We report a significant detection of the hot intergalactic medium in the filamentary bridge connecting the galaxy clusters Abell 399 and Abell 401. This result is enabled by a low-noise, high-resolution map of the thermal Sunyaev-Zeldovich signal from the Atacama Cosmology Telescope (ACT) and Planck satellite. The ACT data provide the $1.65'$ resolution that allows us to clearly separate the profiles of the clusters, whose centres are separated by $37'$, from the gas associated with the filament. A model that fits for only the two clusters is ruled out compared to one that includes a bridge component at $>5\sigma$. Using a gas temperature determined from Suzaku X-ray data, we infer a total mass of $(3.3\pm0.7)\times10^{14}\,\mathrm{M}_{\odot}$ associated with the filament, comprising about $8\%$ of the entire Abell 399-Abell 401 system. We fit two phenomenological models to the filamentary structure; the favoured model has a width transverse to the axis joining the clusters of ${\sim}1.9\,\mathrm{Mpc}$. When combined with the Suzaku data, we find a gas density of $(0.88\pm0.24)\times10^{-4}\,\mathrm{cm}^{-3}$, considerably lower than previously reported. We show that this can be fully explained by a geometry in which the axis joining Abell 399 and Abell 401 has a large component along the line of sight, such that the distance between the clusters is significantly greater than the $3.2\,\mathrm{Mpc}$ projected separation on the plane of the sky. Finally, we present initial results from higher resolution ($12.7"$ effective) imaging of the bridge with the MUSTANG-2 receiver on the Green Bank Telescope., Comment: 23 pages, 10 figures, 3 tables. This is a pre-copyedited, author-produced PDF of an article accepted for publication in the Monthly Notices of the Royal Astronomical Society following peer review. The version of record is available online at: https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/stab3391/6442294

We present the power spectrum methodology used for the first-season COMAP analysis, and assess the quality of the current data set. The main results are derived through the Feed-feed Pseudo-Cross-Spectrum (FPXS) method, which is a robust estimator with respect to both noise modeling errors and experimental systematics. We use effective transfer functions to take into account the effects of instrumental beam smoothing and various filter operations applied during the low-level data processing. The power spectra estimated in this way have allowed us to identify a systematic error associated with one of our two scanning strategies, believed to be due to residual ground or atmospheric contamination. We omit these data from our analysis and no longer use this scanning technique for observations. We present the power spectra from our first season of observing and demonstrate that the uncertainties are integrating as expected for uncorrelated noise, with any residual systematics suppressed to a level below the noise. Using the FPXS method, and combining data on scales $k=0.051-0.62 \,\mathrm{Mpc}^{-1}$ we estimate $P_\mathrm{CO}(k) = -2.7 \pm 1.7 \times 10^4\mu\textrm{K}^2\mathrm{Mpc}^3$, the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum in the literature., Comment: Paper 4 of 7 in series. 18 pages, 11 figures, as accepted in ApJ

Line-intensity mapping observations will find fluctuations of integrated line emission are attenuated by varying degrees at small scales due to the width of the line emission profiles. This attenuation may significantly impact estimates of astrophysical or cosmological quantities derived from measurements. We consider a theoretical treatment of the effect of line broadening on both the clustering and shot-noise components of the power spectrum of a generic line-intensity power spectrum using a halo model. We then consider possible simplifications to allow easier application in analysis, particularly in the context of inferences that require numerous, repeated, fast computations of model line-intensity signals across a large parameter space. For the CO Mapping Array Project (COMAP) and the CO(1-0) line-intensity field at $z\sim3$ serving as our primary case study, we expect a $\sim10\%$ attenuation of the spherically averaged power spectrum on average at relevant scales of $k\approx0.2$-$0.3$ Mpc$^{-1}$, compared to $\sim25\%$ for the interferometric Millimetre-wave Intensity Mapping Experiment (mmIME) targeting shot noise from CO lines at $z\sim1$-$5$ at scales of $k\gtrsim1$ Mpc$^{-1}$. We also consider the nature and amplitude of errors introduced by simplified treatments of line broadening, and find that while an approximation using a single effective velocity scale is sufficient for spherically-averaged power spectra, a more careful treatment is necessary when considering other statistics such as higher multipoles of the anisotropic power spectrum or the voxel intensity distribution., 24 pages + appendix and bibliography (33 pages total), 16 figures, 2 tables; accepted for publication in ApJ

Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg^2 of sky; derived power spectra and cosmological likelihoods; a catalog of over 4,000 galaxy clusters; and supporting ancillary products including beam functions and masks. The data and products are described in a suite of ACT papers; here we provide a summary. In order to facilitate ease of access to these data we present a set of Jupyter IPython notebooks developed to introduce users to DR4, DR5, and the tools needed to analyze these data. The data products (excluding simulations) and the set of notebooks are publicly available on the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA); simulation products are available on the National Energy Research Scientific Computing Center (NERSC)., Accepted to ApJS. 21 pages, 8 figures. Data and notebooks available on LAMBDA https://lambda.gsfc.nasa.gov/product/act/

Line intensity mapping (LIM) is a new technique for tracing the global properties of galaxies over cosmic time. Detection of the very faint signals from redshifted carbon monoxide (CO), a tracer of star formation, pushes the limits of what is feasible with a total-power instrument. The CO Mapping Project (COMAP) Pathfinder is a first-generation instrument aiming to prove the concept and develop the technology for future experiments, as well as delivering early science products. With 19 receiver channels in a hexagonal focal plane arrangement on a 10.4 m antenna, and an instantaneous 26-34 GHz frequency range with 2 MHz resolution, it is ideally suited to measuring CO($J$=1-0) from $z\sim3$. In this paper we discuss strategies for designing and building the Pathfinder and the challenges that were encountered. The design of the instrument prioritized LIM requirements over those of ancillary science. After a couple of years of operation, the instrument is well understood, and the first year of data is already yielding useful science results. Experience with this Pathfinder will drive the design of the next generations of experiments., Comment: Paper 2 of 7 in series. 27 pages, 28 figures, submitted to ApJ

We report on the serendipitous discovery of three transient millimeter-wave sources using data from the Atacama Cosmology Telescope. The first, detected at R.A. = 273.8138, decl. = −49.4628 at ~50σ total, brightened from less than 5 mJy to at least 1100 mJy at 150 GHz with an unknown rise time shorter than 13 days, during which the increase from 250 mJy to 1100 mJy took only 8 minutes. Maximum flux was observed on 2019 November 8. The source's spectral index in flux between 90–150 GHz was positive, α = 1.5 ± 0.2. The second, detected at R.A. = 105.1584, decl. = −11.2434 at ~20σ total, brightened from less than 20 mJy to at least 300 mJy at 150 GHz with an unknown rise time shorter than 8 days. Maximum flux was observed on 2019 December 15. Its spectral index was also positive, α = 1.8 ± 0.2. The third, detected at R.A. = 301.9952, decl. = 16.1652 at ~40σ total, brightened from less than 8 mJy to at least 300 mJy at 150 GHz over a day or less but decayed over a few days. Maximum flux was observed on 2018 September 11. Its spectrum was approximately flat, with a spectral index of α = −0.2 ± 0.1. None of the sources were polarized to the limits of these measurements. The two rising-spectrum sources are coincident in position with M and K stars, while the third is coincident with a G star.

We describe the first season COMAP analysis pipeline that converts raw detector readouts to calibrated sky maps. This pipeline implements four main steps: gain calibration, filtering, data selection, and map-making. Absolute gain calibration relies on a combination of instrumental and astrophysical sources, while relative gain calibration exploits real-time total-power variations. High efficiency filtering is achieved through spectroscopic common-mode rejection within and across receivers, resulting in nearly uncorrelated white noise within single-frequency channels. Consequently, near-optimal but biased maps are produced by binning the filtered time stream into pixelized maps; the corresponding signal bias transfer function is estimated through simulations. Data selection is performed automatically through a series of goodness-of-fit statistics, including $\chi^2$ and multi-scale correlation tests. Applying this pipeline to the first-season COMAP data, we produce a dataset with very low levels of correlated noise. We find that one of our two scanning strategies (the Lissajous type) is sensitive to residual instrumental systematics. As a result, we no longer use this type of scan and exclude data taken this way from our Season 1 power spectrum estimates. We perform a careful analysis of our data processing and observing efficiencies and take account of planned improvements to estimate our future performance. Power spectrum results derived from the first-season COMAP maps are presented and discussed in companion papers., Comment: Paper 3 of 7 in series. 26 pages, 23 figures, submitted to ApJ

We present a detailed overview of the science goals and predictions for the Prime-Cam direct detection camera/spectrometer being constructed by the CCAT-prime collaboration for dedicated use on the Fred Young Submillimeter Telescope (FYST). The FYST is a wide-field, 6-m aperture submillimeter telescope being built (first light in mid-2024) by an international consortium of institutions led by Cornell University and sited at more than 5600 meters on Cerro Chajnantor in northern Chile. Prime-Cam is one of two instruments planned for FYST and will provide unprecedented spectroscopic and broadband measurement capabilities to address important astrophysical questions ranging from Big Bang cosmology through reionization and the formation of the first galaxies to star formation within our own Milky Way galaxy. Prime-Cam on the FYST will have a mapping speed that is over ten times greater than existing and near-term facilities for high-redshift science and broadband polarimetric imaging at frequencies above 300 GHz. We describe details of the science program enabled by this system and our preliminary survey strategies., Comment: 61 pages, 16 figures. Resubmitted to ApJSS July 11, 2022

The anisotropy or triaxiality of massive dark matter haloes largely defines the structure of the cosmic web, in particular the filaments that join the haloes together. Here we investigate such oriented correlations in mass-Peak Patch halo catalogues by using the initial strain tensor of spherical proto-halo regions to orient the haloes. To go beyond the spherically averaged two-point correlation function of haloes we use oriented stacks to compute oriented two-point correlations: we explicitly break isotropy by imposing a local frame set by the strain tensor of the reference halo before stacking neighbouring haloes. Beyond the exclusion zone of the reference halo, clustering is found to be strongly enhanced along the major direction of the strain tensor as expected. This anisotropic clustering of haloes along filaments is further quantified by using a spherical harmonics decomposition. Furthermore, we compute the evolution of cluster-scale halo principal directions relative to those of their neighbours and show that there are strong correlations extending up to very large scales. In order to provide calculations more suitable to observational confrontations, we also utilize 2D projected versions of some equivalent correlation functions. Finally, we show that the multipole structure of the mass-peak patch halo's anisotropic clustering can be qualitatively captured in an analytic treatment based on peak theory. Though highly informative, analytic evaluation involves extensive use of Monte Carlo methods, which is also what the simulated catalogue uses, taking into account as they do the adaptive nature of the mass-peak patch mass hierarchy and all non-local complexities associated with the exclusion of smaller haloes overlapping with larger ones: there is no substitute for the mass-Peak Patch simulation-based determination of oriented and anisotropic correlations., Comment: 22 pages, 15 figures, accepted by MNRAS

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 are available on LAMBDA., Comment: 24 pages, 11 figures, matches version accepted by PRD. Data products are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/act_dr4_derived_maps_get.cfm

One of the most important checks for systematic errors in CMB studies is the cross correlation of maps made by independent experiments. In this paper we report on the cross correlation between maps from the Atacama B-mode Search (ABS) and Atacama Cosmology Telescope (ACT) experiments in both temperature and polarization. These completely different measurements have a clear correlation with each other and with the Planck satellite in both the EE and TE spectra at $\ell, Comment: 7 pages, 5 figures. For an interactive demonstration of the methods, see https://colab.research.google.com/drive/1CnMGLn-J3pySv8A9ffNWDSlXlJWMrp2W

We present the temperature and polarization angular power spectra of the CMB measured by the Atacama Cosmology Telescope (ACT) from 5400 deg$^2$ of the 2013-2016 survey, which covers $>$15000 deg$^2$ at 98 and 150 GHz. For this analysis we adopt a blinding strategy to help avoid confirmation bias and, related to this, show numerous checks for systematic error done before unblinding. Using the likelihood for the cosmological analysis we constrain secondary sources of anisotropy and foreground emission, and derive a "CMB-only" spectrum that extends to $\ell=4000$. At large angular scales, foreground emission at 150 GHz is $\sim$1% of TT and EE within our selected regions and consistent with that found by Planck. Using the same likelihood, we obtain the cosmological parameters for $\Lambda$CDM for the ACT data alone with a prior on the optical depth of $\tau=0.065\pm0.015$. $\Lambda$CDM is a good fit. The best-fit model has a reduced $\chi^2$ of 1.07 (PTE=0.07) with $H_0=67.9\pm1.5$ km/s/Mpc. We show that the lensing BB signal is consistent with $\Lambda$CDM and limit the celestial EB polarization angle to $\psi_P =-0.07^{\circ}\pm0.09^{\circ}$. We directly cross correlate ACT with Planck and observe generally good agreement but with some discrepancies in TE. All data on which this analysis is based will be publicly released., Comment: 44 pages, 27 figures, products available on the NASA LAMBDA website, version accepted for publication in JCAP

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}, Comment: 18 pages, 3 figures, Accepted for publication in PRD

Many models of high energy physics suggest that the cosmological dark sector consists of not just one, but a spectrum of ultralight scalar particles with logarithmically distributed masses. To study the potential signatures of low concentrations of ultralight axion (also known as fuzzy) dark matter, we modify Lagrangian perturbation theory (LPT) by distinguishing between trajectories of different dark matter species. We further adapt LPT to include the effects of a quantum pressure, which is necessary to generate correct initial conditions for ultralight axion simulations. Based on LPT, our modified scheme is extremely efficient on large scales and it can be extended to an arbitrary number of particle species at very little computational cost. This allows for computation of self-consistent initial conditions in mixed dark matter models. Additionally, we find that shell-crossing is delayed for ultralight particles and that the deformation tensor extracted from LPT can be used to identify the range of redshifts and scales for which the Madelung formalism of fuzzy dark matter is a reliable approximation., Comment: Submitted to MNRAS, comments welcome

The Sunyaev–Zel’dovich (SZ) effect introduces a specific distortion of the blackbody spectrum of the cosmic microwave background (CMB) radiation when it scatters off hot gas in clusters of galaxies. The frequency dependence of the distortion is only independent of the cluster redshift when the evolution of the CMB radiation is adiabatic. Using 370 clusters within the redshift range 0.07 ≲ z ≲ 1.4 from the largest SZ-selected cluster sample to date from the Atacama Cosmology Telescope, we provide new constraints on the deviation of CMB temperature evolution from the standard model α = 0.017 − 0.032 + 0.029 , where T ( z ) = T 0 1 + z 1 − α . This result is consistent with no deviation from the standard adiabatic model. Combining it with previous, independent data sets we obtain a joint constraint of α = −0.001 ± 0.012. Attributing deviation from adiabaticity to the decay of dark energy, this result constrains its effective equation of state w eff = − 0.998 − 0.010 + 0.008 .

We correlate the positions of radio galaxies in the FIRST survey with the CMB lensing convergence estimated from the Atacama Cosmology Telescope over 470 square degrees to determine the bias of these galaxies. We remove optically cross-matched sources below redshift $z=0.2$ to preferentially select Active Galactic Nuclei (AGN). We measure the angular cross-power spectrum $C_l^{\kappa g}$ at $4.4\sigma$ significance in the multipole range $100, Comment: 11 pages, 9 figures. Submitted to MNRAS

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.

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 $3001500$, 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., Comment: 10 pages, 10 figures