Extragalactic foregrounds in temperature maps of the Cosmic Microwave Background (CMB) severely limit the ability of standard estimators to reconstruct the weak lensing potential. These foregrounds are not fully removable by multi-frequency cleaning or masking and can lead to large biases if not properly accounted for. For foregrounds made of a number of unclustered point sources, an estimator for the source amplitude can be derived and deprojected, removing any bias to the lensing reconstruction. We show with simulations that all of the extragalactic foregrounds in temperature can be approximated by a collection of sources with identical profiles, and that a simple bias hardening technique is effective at reducing any bias to lensing, at a minimal noise cost. We compare the performance and bias to other methods such as "shear-only" reconstruction, and discuss how to jointly deproject any arbitrary number of foregrounds, each with an arbitrary profile. In particular, for a Simons Observatory-like experiment foreground-hardened estimators allow us to extend the maximum multipole used in the reconstruction, increasing the overall statistical power by $\sim 50\%$ over the standard quadratic estimator, both in auto and cross-correlation. We conclude that source hardening outperforms the standard lensing quadratic estimator both in auto and cross-correlation, and in terms of lensing signal-to-noise and foreground bias.

The epoch of reionization is one of the major phase transitions in the history of the universe, and is a focus of ongoing and upcoming cosmic microwave background (CMB) experiments with improved sensitivity to small-scale fluctuations. Reionization also represents a significant contaminant to CMB-derived cosmological parameter constraints, due to the degeneracy between the Thomson-scattering optical depth, $\tau$, and the amplitude of scalar perturbations, $A_s$. This degeneracy subsequently hinders the ability of large-scale structure data to constrain the sum of the neutrino masses, a major target for cosmology in the 2020s. In this work, we explore the kinematic Sunyaev-Zel'dovich (kSZ) effect as a probe of reionization, and show that it can be used to mitigate the optical depth degeneracy with high-sensitivity, high-resolution data from the upcoming CMB-S4 experiment. We discuss the dependence of the kSZ power spectrum on physical reionization model parameters, as well as on empirical reionization parameters, namely $\tau$ and the duration of reionization, $\Delta z$. We show that by combining the kSZ two-point function and the reconstructed kSZ four-point function, degeneracies between $\tau$ and $\Delta z$ can be strongly broken, yielding tight constraints on both parameters. We forecast $\sigma(\tau) = 0.003$ and $\sigma(\Delta z) = 0.25$ for a combination of CMB-S4 and Planck data, including detailed treatment of foregrounds and atmospheric noise. The constraint on $\tau$ is nearly identical to the cosmic-variance limit that can be achieved from large-angle CMB polarization data. The kSZ effect thus promises to yield not only detailed information about the reionization epoch, but also to enable high-precision cosmological constraints on the neutrino mass.

MB lensing tomography, or the cross-correlation between CMB lensing maps and large-scale structure tracers over a well-defined redshift range, has the potential to map the amplitude and growth of structure over cosmic time, provide some of the most stringent tests of gravity, and break important degeneracies between cosmological parameters. In this work, we use the unWISE galaxy catalog to provide three samples at median redshifts $z \sim 0.6, 1.1$ and 1.5, fully spanning the Dark Energy dominated era, together with the most recent Planck CMB lensing maps. We obtain a combined cross-correlation significance $S/N = 79.3$ over the range of scales $100 < \ell < 1000$. We measure the redshift distribution of unWISE sources by a combination of cross-matching with the COSMOS photometric catalog and cross-correlation with BOSS galaxies and quasars and eBOSS quasars. We also show that magnification bias must be included in our analysis and perform a number of null tests. In a companion paper, we explore the derived cosmological parameters by modeling the non-linearities and propagating the redshift distribution uncertainties.

Reconstruction is becoming a crucial procedure of galaxy clustering analysis for future spectroscopic redshift surveys to obtain subpercent level measurement of the baryon acoustic oscillation scale. Most reconstruction algorithms rely on an estimation of the displacement field from the observed galaxy distribution. However, the displacement reconstruction degrades near the survey boundary due to incomplete data and the boundary effects extend to $\sim100\,\mathrm{Mpc}/h$ within the interior of the survey volume. We study the possibility of using radial velocities measured from the cosmic microwave background observation through the kinematic Sunyaev-Zel'dovich effect to improve performance near the boundary. We find that the boundary effect can be reduced to $\sim30-40\,\mathrm{Mpc}/h$ with the velocity information from Simons Observatory. This is especially helpful for dense low redshift surveys where the volume is relatively small and a large fraction of total volume is affected by the boundary.

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

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

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.

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.

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.

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

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

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

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.

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

CMB lensing from current and upcoming wide-field CMB experiments such as AdvACT, SPT-3G and Simons Observatory relies heavily on temperature (vs. polarization). In this regime, foreground contamination to the temperature map produces significant lensing biases, which cannot be fully controlled by multi-frequency component separation, masking or bias hardening. In this letter, we split the standard CMB lensing quadratic estimator into a new set of optimal "multipole" estimators. On large scales, these multipole estimators reduce to the known magnification and shear estimators, and a new shear B-mode estimator. We leverage the different symmetries of the lensed CMB and extragalactic foregrounds to argue that the shear-only estimator should be approximately immune to extragalactic foregrounds. We build a new method to compute separately and without noise the primary, secondary and trispectrum biases to CMB lensing from foreground simulations. Using this method, we demonstrate that the shear estimator is indeed insensitive to extragalactic foregrounds, even when applied to a single-frequency temperature map contaminated with CIB, tSZ, kSZ and radio point sources. This dramatic reduction in foreground biases allows us to include higher temperature multipoles than with the standard quadratic estimator, thus increasing the total lensing signal-to-noise beyond the quadratic estimator. In addition, magnification-only and shear B-mode estimators provide useful diagnostics for potential residuals. Our Python code LensQuEst to forecast the signal-to-noise of the various estimators, generate mock maps, lense them, and apply the various lensing estimators to them is publicly available at https://github.com/EmmanuelSchaan/LensQuEst .

Over the past decade, the gravitational lensing of the Cosmic Microwave Background (CMB) has become a powerful tool for probing the matter distribution in the Universe. The standard technique used to reconstruct the CMB lensing signal employs the quadratic estimator (QE) method, which has recently been shown to be suboptimal for lensing measurements on very small scales in temperature and polarization data. We implement a simple, more optimal method for the small-scale regime, which involves taking the direct inverse of the background gradient. We derive new techniques to make continuous maps of lensing using this "Gradient-Inversion" (GI) method and validate our method with simulated data, finding good agreement with predictions. For idealized simulations of lensing cross- and autospectra that neglect foregrounds, we demonstrate that our method performs significantly better than previous quadratic estimator methods in temperature; at $L=5000-9000$, it reduces errors on the lensing auto-power spectrum by a factor of $\sim 4$ for both idealized CMB-S4 and Simons Observatory-like experiments and by a factor of $\sim 2.6$ for cross-correlations of CMB-S4-like lensing reconstruction and the true lensing field. We caution that the level of the neglected small-scale foreground power, while low in polarization, is very high in temperature; though we briefly outline foreground mitigation methods, further work on this topic is required. Nevertheless, our results show the future potential for improved small-scale CMB lensing measurements, which could provide stronger constraints on cosmological parameters and astrophysics at high redshifts.

Cosmic Microwave Background (CMB) lensing is a powerful probe of the matter distribution in the Universe. The standard quadratic estimator, which is typically used to measure the lensing signal, is known to be suboptimal for low-noise polarization data from next-generation experiments. In this paper we explain why the quadratic estimator will also be suboptimal for measuring lensing on very small scales, even for measurements in temperature where this estimator typically performs well. Though maximum likelihood methods could be implemented to improve performance, we explore a much simpler solution, revisiting a previously proposed method to measure lensing which involves a direct inversion of the background gradient. An important application of this simple formalism is the measurement of cluster masses with CMB lensing. We find that directly applying a gradient inversion matched filter to simulated lensed images of the CMB can tighten constraints on cluster masses compared to the quadratic estimator. While the difference is not relevant for existing surveys, for future surveys it can translate to significant improvements in mass calibration for distant clusters, where galaxy lensing calibration is ineffective due to the lack of enough resolved background galaxies. Improvements can be as large as $\sim 50\%$ for a cluster at $z = 2$ and a next-generation CMB experiment with 1$\mu$K-arcmin noise, and over an order of magnitude for lower noise levels. For future surveys, this simple matched-filter or gradient inversion method approaches the performance of maximum likelihood methods, at a fraction of the computational cost.

Galaxy redshift surveys, such as 2dF, SDSS, 6df, GAMA and VIPERS, have shown that the spatial distribution of matter forms a rich web, known as the cosmic web. The majority of galaxy survey analyses measure the amplitude of galaxy clustering as a function of scale, ignoring information beyond a small number of summary statistics. Since the matter density field becomes highly non-Gaussian as structure evolves under gravity, we expect other statistical descriptions of the field to provide us with additional information. One way to study the non-Gaussianity is to study filaments, which evolve non-linearly from the initial density fluctuations produced in the primordial Universe. In our study, we report the first detection of CMB (Cosmic Microwave Background) lensing by filaments and we apply a null test to confirm our detection. Furthermore, we propose a phenomenological model to interpret the detected signal and we measure how filaments trace the matter distribution on large scales through filament bias, which we measure to be around 1.5. Our study provides a new scope to understand the environmental dependence of galaxy formation. In the future, the joint analysis of lensing and Sunyaev-Zel'dovich observations might reveal the properties of `missing baryons', the vast majority of the gas which resides in the intergalactic medium and has so far evaded most observations.

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.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

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.

Direct measurement of acceleration is a key scientific goal for both cosmology and exoplanets. For cosmology, the concept of redshift drift (more than 60 years old by the 2020s) could directly establish the Friedmann-Lema{\^\i}tre-Robertson-Walker model. It would increase the dark energy figure of merit by a factor of 3 beyond Stage 4 experiments, in combination with cosmic microwave background measurements. For exoplanets, the same technology required provides unprecedented radial velocity accuracy, enabling detection of Earth mass planets in the habitable zone. Other science cases include mapping the Milky Way gravitational potential and testing its dark matter distribution.

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.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

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

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

We study the cross-correlation between the Planck CMB lensing convergence map and the eBOSS quasar overdensity obtained from the Sloan Digital Sky Survey (SDSS) IV, in the redshift range $0.9 < z < 2.2$. We detect the CMB lensing convergence-quasar cross power spectrum at $5.4 \sigma$ significance. The cross power spectrum provides a quasar clustering bias measurement that is expected to be particularly robust against systematic effects. The redshift distribution of the quasar sample has a median redshift $z \approx 1.55$, and an effective redshift about $1.51$. The best fit bias of the quasar sample is $b_q = 2.43 \pm 0.45$, corresponding to a host halo mass of $\log_{10}\left( \frac{M}{h^{-1} M_\odot} \right) = 12.54^{+0.25}_{-0.36}$. This is broadly consistent with the previous literature on quasars with a similar redshift range and selection. Since our constraint on the bias comes from the cross-correlation between quasars and CMB lensing, we expect it to be robust to a wide range of possible systematic effects that may contaminate the auto correlation of quasars. We checked for a number of systematic effects from both CMB lensing and quasar overdensity, and found that all systematics are consistent with null within $2 \sigma$. The data is not sensitive to a possible scale dependence of the bias at present, but we expect that as the number of quasars increases (in future surveys such as DESI), it is likely that strong constraints on the scale dependence of the bias can be obtained.

One of the most powerful cosmological datasets when it comes to constraining neutrino masses is represented by galaxy power spectrum measurements, $P_{gg}(k)$. The constraining power of $P_{gg}(k)$ is however severely limited by uncertainties in the modeling of the scale-dependent galaxy bias $b(k)$. In this Letter we present a new method to constrain $b(k)$ by using the cross-correlation between the Cosmic Microwave Background (CMB) lensing signal and galaxy maps ($C_\ell^{\rm \kappa g}$) using a simple but theoretically well-motivated parametrization for $b(k)$. We apply the method using $C_\ell^{\rm \kappa g}$ measured by cross-correlating Planck lensing maps and the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 (DR11) CMASS galaxy sample, and $P_{gg}(k)$ measured from the BOSS DR12 CMASS sample. We detect a non-zero scale-dependence at moderate significance, which suggests that a proper modeling of $b(k)$ is necessary in order to reduce the impact of non-linearities and minimize the corresponding systematics. The accomplished increase in constraining power of $P_{gg}(k)$ is demonstrated by determining a 95% C.L. upper bound on the sum of the three active neutrino masses $M_{\nu}$ of $M_{\nu}<0.19\, {\rm eV}$. This limit represents a significant improvement over previous bounds with comparable datasets. Our method will prove especially powerful and important as future large-scale structure surveys will overlap more significantly with the CMB lensing kernel providing a large cross-correlation signal.

Kinetic Sunyaev Zel'dovich (kSZ) tomography provides a powerful probe of the radial velocity field of matter in the Universe. By cross-correlating a high resolution CMB experiment like CMB S4 and a galaxy survey like DESI or LSST, one can measure the radial velocity field with very high signal to noise over a large volume of the universe. In this paper we show how this measurement can be used to improve constraints on primordial non-Gaussianities of the local type. The velocity field provides a measurement of the unbiased matter perturbations on large scales, which can be cross-correlated with the biased large-scale galaxy density field. This results in sample variance cancellation for a measurement of scale-dependent bias due to a non-zero $f_{NL}$. Using this method we forecast that CMB S4 and LSST combined reach a sensitivity $\sigma_{f_{NL}} \sim 0.5$, which is a factor of three improvement over the sensitivity using LSST alone (without internal sample variance cancellation). We take into account critical systematics like photometric redshifts, the kSZ optical depth degeneracy, and systematics affecting the shape of the galaxy auto-power spectrum and find that these have negligible impact, thus making kSZ tomography a robust probe for primordial non-Gaussianities. We also forecast the impact of mass binning on our constraints. The techniques proposed in this paper could be an important component of achieving the theoretically important threshold of $\sigma_{f_{NL}} \lesssim 1$ with future surveys.

Several statistics have been proposed for measuring the kSZ effect by combining the small-scale CMB with galaxy surveys. We review five such statistics, and show that they are all mathematically equivalent to the optimal bispectrum estimator of type $\langle ggT \rangle$. Reinterpreting these kSZ statistics as special cases of bispectrum estimation makes many aspects transparent, for example optimally weighting the estimator, or incorporating photometric redshift errors. We analyze the information content of the bispectrum and show that there are two observables: the small-scale galaxy-electron power spectrum $P_{ge}(k_S)$, and the large-scale galaxy-velocity power spectrum $P_{gv}(k)$. The cosmological constraining power of the kSZ arises from its sensitivity to fluctuations on large length scales, where its effective noise level can be much better than galaxy surveys.

With low redshift probes reaching unprecedented precision, uncertainty of the CMB optical depth is expected to be the limiting factor for future cosmological neutrino mass constraints. In this paper, we discuss to what extent combinations of CMB lensing and galaxy surveys measurements at low redshifts $z\sim 0.5-5$ will be able to make competitive neutrino mass measurements without relying on any optical depth constraints. We find that the combination of LSST galaxies and CMB-S4 lensing should be able to achieve constraints on the neutrino mass sum of 25meV without optical depth information, an independent measurement that is competitive with or slightly better than the constraint of 30meV possible with CMB-S4 and present-day optical depth measurements. These constraints originate both in structure growth probed by cross-correlation tomography over a wide redshift range as well as, most importantly, the shape of the galaxy power spectrum measured over a large volume. We caution that possible complications such as higher-order biasing and systematic errors in the analysis of high redshift galaxy clustering are only briefly discussed and may be non-negligible. Nevertheless, our results show that new kinds of high-precision neutrino mass measurements at and beyond the present-day optical depth limit may be possible.

Gravitational lensing deflects the paths of cosmic infrared background (CIB) photons, leaving a measurable imprint on CIB maps. The resulting statistical anisotropy can be used to reconstruct the matter distribution out to the redshifts of CIB sources. To this end, we generalize the CMB lensing quadratic estimator to any weakly non-Gaussian source field, by deriving the optimal lensing weights. We point out the additional noise and bias caused by the non-Gaussianity and the `self-lensing' of the source field. We propose methods to reduce, subtract or model these non-Gaussianities. We show that CIB lensing should be detectable with Planck data, and detectable at high significance for future CMB experiments like CCAT-Prime. The CIB thus constitutes a new source image for lensing studies, providing constraints on the amplitude of structure at intermediate redshifts between galaxies and the CMB. CIB lensing measurements will also give valuable information on the star formation history in the universe, constraining CIB halo models beyond the CIB power spectrum. By laying out a detailed treatment of lens reconstruction from a weakly non-Gaussian source field, this work constitutes a stepping stone towards lens reconstruction from continuum or line intensity mapping data, such as the Lyman-alpha emission, absorption, and the 21cm radiation.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. 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. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

Patchy reionization leaves a number of imprints on the small-scale cosmic microwave background (CMB) temperature fluctuations, the largest of which is the kinematic Sunyaev-Zel'dovich (kSZ), the Doppler shift of CMB photons scattering off moving electrons in ionized bubbles. It has long been known that in the CMB power spectrum, this imprint of reionization is largely degenerate with the kSZ signal produced by late-time galaxies and clusters, thus limiting our ability to constrain reionization. Following Smith & Ferraro (2017), it is possible to isolate the reionization contribution in a model independent way, by looking at the large scale modulation of the small scale CMB power spectrum. In this paper we extend the formalism to use the full shape information of the small scale power spectrum (rather than just its broadband average), and argue that this is necessary to break the degeneracy between the optical depth $\tau$ and parameters setting the duration of reionization. In particular, we show that the next generation of CMB experiments could achieve up to a factor of 3 improvement on the optical depth $\tau$ and at the same time, constrain the duration of reionization to $\sim$ 25 %. This can help tighten the constrains on neutrino masses, which will be limited by our knowledge of $\tau$, and shed light on the physical processes responsible for reionization.

Gravitational lensing of the cosmic microwave background (CMB) is expected to be amongst the most powerful cosmological tools for ongoing and upcoming CMB experiments. In this work, we investigate a bias to CMB lensing reconstruction from temperature anisotropies due to the kinematic Sunyaev-Zel'dovich (kSZ) effect, that is, the Doppler shift of CMB photons induced by Compton-scattering off moving electrons. The kSZ signal yields biases due to both its own intrinsic non-Gaussianity and its non-zero cross-correlation with the CMB lensing field (and other fields that trace the large-scale structure). This kSZ-induced bias affects both the CMB lensing auto-power spectrum and its cross-correlation with low-redshift tracers. Furthermore, it cannot be removed by multifrequency foreground separation techniques because the kSZ effect preserves the blackbody spectrum of the CMB. While statistically negligible for current datasets, we show that it will be important for upcoming surveys, and failure to account for it can lead to large biases in constraints on neutrino masses or the properties of dark energy. For a Stage 4 CMB experiment, the bias can be as large as $\approx$ 15% or 12% in cross-correlation with LSST galaxy lensing convergence or galaxy overdensity maps, respectively, when the maximum temperature multipole used in the reconstruction is $\ell_{\rm max} = 4000$, and about half of that when $\ell_{\rm max} = 3000$. Similarly, we find that the CMB lensing auto-power spectrum can be biased by up to several percent. These biases are many times larger than the expected statistical errors. Reducing $\ell_{\rm max}$ can significantly mitigate the bias at the cost of a decrease in the overall lensing reconstruction signal-to-noise. Polarization-only reconstruction may be the most robust mitigation strategy.

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.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

The CMB Stage 4 (CMB-S4) experiment is a next-generation, ground-based experiment that will measure the cosmic microwave background (CMB) polarization to unprecedented accuracy, probing the signature of inflation, the nature of cosmic neutrinos, relativistic thermal relics in the early universe, and the evolution of the universe. To advance the progress towards designing the instrument for CMB-S4, we have established a framework to optimize the instrumental configuration to maximize its scientific output. In this paper, we report our first results from this framework, using simplified instrumental and cost models. We have primarily studied two classes of instrumental configurations: arrays of large aperture telescopes with diameters ranging from 2-10 m, and hybrid arrays that combine small-aperture telescopes (0.5 m diameter) with large-aperture telescopes. We explore performance as a function of the telescope aperture size, the distribution of the detectors into different microwave frequencies, the survey strategy and survey area, the low-frequency noise performance, and the balance between small and large aperture telescopes for the hybrid configurations. We also examine the impact from the uncertainties of the instrumental model. There are several areas that deserve further improvement. In our forecasting framework, we adopt a simple two-component foregrounds model with spacially varying power-law spectral indices. We estimate delensing performance statistically and ignore possible non-idealities. Instrumental systematics, which is not accounted for in our study, may influence the design. Further study of the instrumental and cost models will be one of the main areas of study by the whole CMB-S4 community. We hope that our framework will be useful for estimating the influence of these improvement in future, and we will incorporate them in order to improve the optimization further.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

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

We present a Bayesian phase-space reconstruction of the cosmic large-scale matter density and velocity fields from the SDSS-III Baryon Oscillations Spectroscopic Survey Data Release 12 (BOSS DR12) CMASS galaxy clustering catalogue. We rely on a given $\Lambda$CDM cosmology, a mesh resolution in the range of 6-10 $h^{-1}$ Mpc, and a lognormal-Poisson model with a redshift dependent nonlinear bias. The bias parameters are derived from the data and a general renormalised perturbation theory approach. We use combined Gibbs and Hamiltonian sampling, implemented in the \textsc{argo} code, to iteratively reconstruct the dark matter density field and the coherent peculiar velocities of individual galaxies, correcting hereby for coherent redshift space distortions (RSD). Our tests relying on accurate $N$-body based mock galaxy catalogues, show unbiased real space power spectra of the nonlinear density field up to $k\sim0.2\, h$ Mpc$^{-1}$, and vanishing quadrupoles down to $r\sim20\,h^{-1}$ Mpc. We also demonstrate that the nonlinear cosmic web can be obtained from the tidal field tensor based on the Gaussian component of the reconstructed density field. We find that the reconstructed velocities have a statistical correlation coefficient compared to the true velocities of each individual lightcone mock galaxy of $r\sim0.68$ including about 10% of satellite galaxies with virial motions (about $r=0.75$ without satellites). The power spectra of the velocity divergence agree well with theoretical predictions up to $k\sim0.2\,h\,{\rm Mpc}^{-1}$. This work will be especially useful to improve, e.g. BAO reconstructions, kinematic Sunyaev-Zeldovich (kSZ), integrated Sachs-Wolfe (ISW) measurements, or environmental studies.

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.

Upcoming cosmic microwave background (CMB) experiments will measure temperature fluctuations on small angular scales with unprecedented precision. Small-scale CMB fluctuations are a mixture of late-time effects: gravitational lensing, Doppler shifting of CMB photons by moving electrons (the kSZ effect), and residual foregrounds. We propose a new statistic which separates the kSZ signal from the others, and also allows the kSZ signal to be decomposed in redshift bins. The decomposition extends to high redshift, and does not require external datasets such as galaxy surveys. In particular, the high-redshift signal from patchy reionization can be cleanly isolated, enabling future CMB experiments to make high-significance and qualitatively new measurements of the reionization era.

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.

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.

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.

Primordial non-Gaussianity of local type is known to produce a scale-dependent contribution to the galaxy bias. Several classes of multi-field inflationary models predict non-Gaussian bias which is stochastic, in the sense that dark matter and halos don't trace each other perfectly on large scales. In this work, we forecast the ability of next-generation Large Scale Structure surveys to constrain common types of primordial non-Gaussianity like $f_{NL}$, $g_{NL}$ and $\tau_{NL}$ using halo bias, including stochastic contributions. We provide fitting functions for statistical errors on these parameters which can be used for rapid forecasting or survey optimization. A next-generation survey with volume $V = 25 h^{-3}$Gpc$^3$, median redshift $z = 0.7$ and mean bias $b_g = 2.5$, can achieve $\sigma(f_{NL}) = 6$, $\sigma(g_{NL}) = 10^5$ and $\sigma(\tau_{NL}) = 10^3$ if no mass information is available. If halo masses are available, we show that optimally weighting the halo field in order to reduce sample variance can achieve $\sigma(f_{NL}) = 1.5$, $\sigma(g_{NL}) = 10^4$ and $\sigma(\tau_{NL}) = 100$ if halos with mass down to $M_{min} = 10^{11}$ $h^{-1} M_\odot $ are resolved, outperforming Planck by a factor of 4 on $f_{NL}$ and nearly an order of magnitude on $g_{NL}$ and $\tau_{NL}$. Finally, we study the effect of photometric redshift errors and discuss degeneracies between different non-Gaussian parameters, as well as the impact of marginalizing Gaussian bias and shot noise.

Photometric observations of galaxies at submillimetre to millimetre wavelengths (50 - 1000 GHz) are susceptible to spatial variations in both the background CMB temperature and CIB emission that can be comparable to the flux from the target galaxy. We quantify the residual uncertainty when background emission inside a circular aperture is estimated by the mean flux in a surrounding annular region, assumed to have no contribution from the source of interest. We present simple formulae to calculate this uncertainty as a function of wavelength and aperture size. Drawing on examples from the literature, we illustrate the use of our formalism in practice and highlight cases in which uncertainty in the background subtraction needs to be considered in the error analysis. We make the code used to calculate the uncertainties publicly available on the web.

The Integrated Sachs-Wolfe effect (ISW) measures the decay of the gravitational potential due to cosmic acceleration and is thus a direct probe of Dark Energy. In some of the earlier studies, the amplitude of the ISW effect was found to be in tension with the predictions of the standard $\Lambda$CDM model. We measure the cross-power of galaxies and AGN from the WISE mission with CMB temperature data from WMAP9 in order to provide an independent measurement of the ISW amplitude. Cross-correlations with the recently released Planck lensing potential maps are used to calibrate the bias and contamination fraction of the sources, thus avoiding systematic effects that could be present when using auto-spectra to measure bias. We find an amplitude of the cross-power of $\mathcal{A} = 1.24\pm 0.47$ from the galaxies and $\mathcal{A} = 0.88 \pm 0.74$ from the AGN, fully consistent with the $\Lambda$CDM prediction of $\mathcal{A} =1$. The ISW measurement signal-to-noise ratio is 2.7 and 1.2 respectively, giving a combined significance close to $3 \sigma$. Comparing the amplitudes of the galaxy and AGN cross-correlations, which arise from different redshifts, we find no evidence for redshift evolution in Dark Energy properties, consistent with a Cosmological Constant.

Several, apparently distinct, formalisms exist in the literature for predicting the clustering of dark matter halos. It has been noticed on a case-by-case basis that the predictions of these different methods agree in specific examples, but there is no general proof that they are equivalent. In this paper, we give a simple proof of the mathematical equivalence of barrier crossing, peak-background split, and local biasing.

In this article, we show that a stochastic form of scale-dependent halo bias arises in multi-source inflationary models, where multiple fields determine the initial curvature perturbation. We derive this effect for general non-Gaussian initial conditions and study various examples, such as curvaton models and quasi-single field inflation. We present a general formula for both the stochastic and the non-stochastic parts of the halo bias, in terms of the N-point cumulants of the curvature perturbation at the end of inflation. At lowest order, the stochasticity arises if the collapsed limit of the four-point function is boosted relative to the square of the three-point function in the squeezed limit. We derive all our results in two ways, using the barrier crossing formalism and the peak-background split method. In a companion paper, we prove that these two approaches are mathematically equivalent.

It has recently been shown that supersonic relative velocities between dark matter and baryonic matter can have a significant effect on formation of the first structures in the universe. If this effect is still non-negligible during the epoch of hydrogen reionization, it generates large-scale anisotropy in the free electron density, which gives rise to a CMB B-mode. We compute the B-mode power spectrum and find a characteristic shape with acoustic peaks at l ~ 200, 400, ... The amplitude of this signal is a free parameter which is related to the dependence of the ionization fraction on the relative baryon-CDM velocity during the epoch of reionization. However, we find that the B-mode signal is undetectably small for currently favored reionization models in which hydrogen is reionized promptly at z ~ 10, although constraints on this signal by future experiments may help constrain models in which partial reionization occurs at higher redshift, e.g. by accretion onto primordial black holes.

Cosmic inflation provides a mechanism for generating the early density perturbations that seeded the large-scale structures we see today. Primordial non-Gaussianity is among the most promising of few observational tests of physics at this epoch. At present non-Gaussianity is best constrained by the cosmic microwave background, but in the near term large-scale structure data may be competitive so long as the effects of primordial non-Gaussianity can be modeled through the non-linear process of structure formation. We discuss recent work modeling effects of a few types of primordial non-Gaussianity on the large-scale halo clustering and the halo mass function. More specifically, we compare analytic and N-body results for two variants of the curvaton model of inflation: (i) a "tauNL" scenario in which the curvaton and inflaton contribute equally to the primordial curvature perturbation and (ii) a "gNL" model where the usual quadratic fNL term in the potential cancels, but a large cubic term remains.

A wide range of multifield inflationary models generate non-Gaussian initial conditions in which the initial adiabatic fluctuation is of the form (zeta_G + g_{NL} zeta_G^3). We study halo clustering in these models using two different analytic methods: the peak-background split framework, and brute force calculation in a barrier crossing model, obtaining agreement between the two. We find a simple, theoretically motivated expression for halo bias which agrees with N-body simulations and can be used to constrain g_{NL} from observations. We discuss practical caveats to constraining g_{NL} using only observable properties of a tracer population, and argue that constraints obtained from populations whose observed bias is <~ 2.5 are generally not robust to uncertainties in modeling the halo occupation distribution of the population.

As one of the most powerful probes of cosmological structure formation, the abundance of massive galaxy clusters is a sensitive probe of modifications to gravity on cosmological scales. In this paper, we present results from N-body simulations of a general class of f(R) models, which self-consistently solve the non-linear field equation for the enhanced forces. Within this class we vary the amplitude of the field, which controls the range of the enhanced gravitational forces, both at the present epoch and as a function of redshift. Most models in the literature can be mapped onto the parameter space of this class. Focusing on the abundance of massive dark matter halos, we compare the simulation results to a simple spherical collapse model. Current constraints lie in the large-field regime, where the chameleon mechanism is not important. In this regime, the spherical collapse model works equally well for a wide range of models and can serve as a model-independent tool for placing constraints on f(R) gravity from cluster abundance. Using these results, we show how constraints from the observed local abundance of X-ray clusters on a specific f(R) model can be mapped onto other members of this general class of models.