Search

The 21-cm brightness-temperature field of neutral hydrogen during the Epoch of Reionization and Cosmic Dawn is a rich source of cosmological and astrophysical information, primarily due to its significant non-Gaussian features. However, the complex, nonlinear nature of the underlying physical processes makes analytical modelling of this signal challenging. Consequently, studies often resort to semi-numerical simulations. Traditional analysis methods, which rely on a limited set of summary statistics, may not adequately capture the non-Gaussian content of the data, as the most informative statistics are not predetermined. This paper explores the application of machine learning (ML) to surpass the limitations of summary statistics by leveraging the inherent non-Gaussian characteristics of the 21-cm signal. We demonstrate that a well-trained neural network can independently reconstruct the hydrogen density, spin-temperature, and neutral-fraction fields with cross-coherence values exceeding 0.95 for $k$-modes below $0.5$ Mpc h$^{-1}$, based on a representative simulation at a redshift of $z \approx 15$. To achieve this, the neural network utilises the non-Gaussian information in brightness temperature images over many scales. We discuss how these reconstructed fields, which vary in their sensitivity to model parameters, can be employed for parameter inference, offering more direct insights into underlying cosmological and astrophysical processes only using limited summary statistics of the brightness temperature field, such as its power spectrum., Comment: 9 pages, 6 figures, prepared for submission to MNRAS

The redshifted 21-cm signal from the Cosmic Dawn and Epoch of Reionization carries invaluable information about the cosmology and astrophysics of the early Universe. Analyzing the data from a sky-averaged 21-cm signal experiment typically involves navigating through an intricate parameter space to accurately address various factors such as foregrounds, beam uncertainties, ionospheric distortions, and receiver noise for the search of the cosmological 21-cm signal. The traditional likelihood-based sampling methods for modeling these effects could become computationally demanding for such highly complex models, which makes it infeasible to include physically motivated 21-cm signal models in the analysis. Moreover, the inference with these traditional methods is driven by the assumed functional form of the likelihood function. This work demonstrates how Simulation-Based Inference through Truncated Marginal Neural Ratio Estimation (TMNRE) can naturally handle these issues at a significantly reduced computational cost than the likelihood-based methods. We estimate the posterior distribution on our model parameters with TMNRE for simulated mock observations, composed of beam-weighted foregrounds, physically motivated 21-cm signal, and radiometric noise. We find that maximizing the information content by simultaneously analyzing the data from multiple time slices and antennas significantly improves the parameter constraints and leads to a better exploration of the cosmological signal. We discuss the application of TMNRE for the current configuration of the REACH experiment and demonstrate how it can be utilized to explore potential avenues for REACH. The method presented here can be easily generalized for any sky-averaged 21-cm signal experiment., Comment: 13 pages, 10 figures

Inspired by the discussion in the community on possible hidden systematic errors in late universe cosmological probes and non-trivial physical models developed to reduce the Hubble tension, we investigate the Pantheon and Pantheon+ SNe samples for possible deviations from the original $\Lambda$CDM analysis. To simultaneously account for possible systematics or deviations from $\Lambda$CDM, we adopt Gaussian processes to model additional covariance while making no further assumptions on their origin. We explore both stationary and non-stationarity corrections to the covariance. While small changes in the inferred cosmological parameters $H_0$ and $\Omega_{m}$ can occur, we find no statistically significant evidence for missing covariance. We find an upper limit for the Gaussian processes amplitude $\sigma < 0.031$ mag with $95\%$ confidence, which corresponds to $20\%$ of the average statistical error in the Pantheon+ sample. The strongest effect we find on the inferred cosmological parameter posterior can reduce the statistical significance of the Hubble tension between Pantheon+ and Planck estimates from 5.3$\sigma$ to 4.5$\sigma$. Therefore, we conclude that the SN cosmological parameter inference is robust against the analysis modifications studied in this work., Comment: 16 pages, 6 figures, prepared for submission to JCAP

The 21-cm brightness temperature fluctuation from the Dark Ages ($z \simeq 30-100$) will allow us to probe the inflationary epoch on very small scales ($>0.1 \, \mbox{Mpc}^{-1}$), inaccessible to cosmic microwave background experiments. Combined with the possibility to collect information from different redshift slices, the 21-cm bispectrum has the potential to significantly improve constraints on primordial non-Gaussianity. However, recent work has shown secondary effects source off-diagonal terms in the covariance matrix which can significantly affect forecasted constraints, especially in signals that peak in the squeezed configuration, such as the local bispectrum. In this paper we propose the three-point $\langle 21-21-\rm CMB \rangle$ bispectrum cross-correlation as a new independent observational channel sensitive to local primordial non-Gaussianity. We find that, contrary to the 21-cm bispectrum, secondary contributions are subdominant to the primordial signal for values $f_{\rm NL}^{\rm loc} \sim 1$, resulting in negligible effects from off-diagonal terms in the covariance matrix. We forecast that CMB $T$ and $E$ modes cross-correlated with an ideal cosmic variance-limited 21-cm experiment with a $0.1$ MHz frequency and $0.1$ arc-minute angular resolution could reach $f_{\rm NL}^{\rm loc} \sim 6 \times 10^{-3}$. This forecast suggests cross-correlation between CMB and 21-cm experiments could provide a viable alternative to 21-cm auto-spectra in reaching unprecedented constraints on primordial local non-Gaussianities., Comment: 38 pages, 14 figures, code available at https://github.com/GiorgioOrlando/21cm_CMB_cross

We detail the the REACH radiometric system designed to enable measurements of the 21-cm neutral hydrogen line. Included is the radiometer architecture and end-to-end system simulations as well as a discussion of the challenges intrinsic to highly-calibratable system development. Following this, we share laboratory results based on the calculation of noise wave parameters utilising an over-constrained least squares approach demonstrating a calibration RMSE of 80 mK for five hours of integration on a custom-made source with comparable impedance to that of the antenna used in the field. This paper therefore documents the state of the calibrator and data analysis in December 2022 in Cambridge before shipping to South Africa., Comment: 30 pages, 19 figures

We study the use of U-Nets in reconstructing the linear dark matter density field and its consequences for constraining cosmological parameters, in particular primordial non-Gaussianity. Our network is able to reconstruct the initial conditions of redshift $z=0$ density fields from N-body simulations with $90\%$ accuracy out to $k \leq 0.4$ h/Mpc, competitive with state-of-the-art reconstruction algorithms at a fraction of the computational cost. We study the information content of the reconstructed $z=0$ density field with a Fisher analysis using the QUIJOTE simulation suite, including non-Gaussian initial conditions. Combining the pre- and post-reconstructed power spectrum and bispectrum data up to $k_{\rm max} = 0.52$ h/Mpc, we find significant improvements on all parameters. Most notably, we find a factor $3.65$ (local), $3.54$ (equilateral) and $2.90$ (orthogonal) improvement on the marginalized errors of $f_{\rm NL}$ as compared to only using the pre-reconstructed data. We show that these improvements can be attributed to a combination of reduced data covariance and parameter degeneracy. The results constitute an important step towards more optimal inference of primordial non-Gaussianity from non-linear scales., Comment: 22 pages, 8 figures, 3 tables, codes available at https://github.com/tsfloss/URecon and https://github.com/tsfloss/DensityFieldTools. v2 matches version accepted for JCAP

Cosmic Dawn (CD) and Epoch of Reionization (EoR) are epochs of the Universe which host invaluable information about the cosmology and astrophysics of X-ray heating and hydrogen reionization. Radio interferometric observations of the 21-cm line at high redshifts have the potential to revolutionize our understanding of the universe during this time. However, modeling the evolution of these epochs is particularly challenging due to the complex interplay of many physical processes. This makes it difficult to perform the conventional statistical analysis using the likelihood-based Markov-Chain Monte Carlo (MCMC) methods, which scales poorly with the dimensionality of the parameter space. In this paper, we show how the Simulation-Based Inference (SBI) through Marginal Neural Ratio Estimation (MNRE) provides a step towards evading these issues. We use 21cmFAST to model the 21-cm power spectrum during CD-EoR with a six-dimensional parameter space. With the expected thermal noise from the Square Kilometre Array (SKA), we are able to accurately recover the posterior distribution for the parameters of our model at a significantly lower computational cost than the conventional likelihood-based methods. We further show how the same training dataset can be utilized to investigate the sensitivity of the model parameters over different redshifts. Our results support that such efficient and scalable inference techniques enable us to significantly extend the modeling complexity beyond what is currently achievable with conventional MCMC methods., Comment: 15 pages, 9 figures. Accepted for publication in MNRAS

Diffusion damping of the cosmic microwave background (CMB) power spectrum results from imperfect photon-baryon coupling in the pre-recombination plasma. At redshift $5 \times 10^4 < z < 2 \times 10^6$, the plasma acquires an effective chemical potential, and energy injections from acoustic damping in this era create $\mu$-type spectral distortions of the CMB. These $\mu$ distortions trace the underlying photon density fluctuations, probing the primordial power spectrum in short-wavelength modes $k_\mathrm{S}$ over the range $50 \ \mathrm{Mpc}^{-1} \lesssim k \lesssim 10^4 \ \mathrm{Mpc}^{-1}$. Small-scale power modulated by long-wavelength modes $k_\mathrm{L}$ from squeezed-limit non-Gaussianities introduces cross-correlations between CMB temperature anisotropies and $\mu$ distortions. Under single-field inflation models, $\mu \times T$ correlations measured from an observer in an inertial frame should vanish up to a factor of $(k_\mathrm{L}/k_\mathrm{S})^2 \ll 1$. Thus, any measurable correlation rules out single-field inflation models. We forecast how well the next-generation ground-based CMB experiment CMB-S4 will be able to constrain primordial squeezed-limit non-Gaussianity, parameterized by $f_\mathrm{NL}$, using measurements of $C_{\ell}^{\mu T}$ as well as $C_{\ell}^{\mu E}$ from CMB $E$ modes. Using current experimental specifications and foreground modeling, we expect $\sigma(f_\mathrm{NL}) \lesssim 1000$. This is roughly four times better than the current limit on $f_\mathrm{NL}$ using $\mu \times T$ and $\mu \times E$ correlations from Planck and is comparable to what is achievable with LiteBIRD, demonstrating the power of the CMB-S4 experiment. This measurement is at an effective scale of $k \simeq 740 \ \text{Mpc}^{-1}$ and is thus highly complementary to measurements at larger scales from primary CMB and large-scale structure., Comment: 19 pages, 6 figures

In a sky-averaged 21-cm signal experiment, the uncertainty on the extracted signal depends mainly on the covariance between the foreground and 21-cm signal models. In this paper, we construct these models using the modes of variation obtained from the Singular Value Decomposition of a set of simulated foreground and 21-cm signals. We present a strategy to reduce this overlap between the 21-cm and foreground modes by simultaneously fitting the spectra from multiple different antennas, which can be used in combination with the method of utilizing the time dependence of foregrounds while fitting multiple drift scan spectra. To demonstrate this idea, we consider two different foreground models (i) a simple foreground model, where we assume a constant spectral index over the sky, and (ii) a more realistic foreground model, with a spatial variation of the spectral index. For the simple foreground model, with just a single antenna design, we are able to extract the signal with good accuracy if we simultaneously fit the data from multiple time slices. The 21-cm signal extraction is further improved when we simultaneously fit the data from different antennas as well. This improvement becomes more pronounced while using the more realistic mock observations generated from the detailed foreground model. We find that even if we fit multiple time slices, the recovered signal is biased and inaccurate for a single antenna. However, simultaneously fitting the data from different antennas reduces the bias and the uncertainty by a factor of 2-3 on the extracted 21-cm signal., Comment: 12 pages, 13 figures. Accepted for publication in MNRAS. Accompanying code is available https://github.com/anchal-009/SAVED21cm

Weak gravitational lensing by the intervening large-scale structure (LSS) of the Universe is the leading non-linear effect on the anisotropies of the cosmic microwave background (CMB). The integrated line-of-sight mass that causes the distortion -- known as lensing convergence -- can be reconstructed from the lensed temperature and polarization anisotropies via estimators quadratic in the CMB modes, and its power spectrum has been measured from multiple CMB experiments. Sourced by the non-linear evolution of structure, the bispectrum of the lensing convergence provides additional information on late-time cosmological evolution complementary to the power spectrum. However, when trying to estimate the summary statistics of the reconstructed lensing convergence, a number of noise-biases are introduced, as previous studies have shown for the power spectrum. Here, we explore for the first time the noise-biases in measuring the bispectrum of the reconstructed lensing convergence. We compute the leading noise-biases in the flat-sky limit and compare our analytical results against simulations, finding excellent agreement. Our results are critical for future attempts to reconstruct the lensing convergence bispectrum with real CMB data., Comment: 40 pages (including refs), 19 figures, comments very welcome

Observations of the 21-cm line from primordial hydrogen promise to be one of the best tools to study the early epochs of the Universe: the Dark Ages, the Cosmic Dawn, and the subsequent Epoch of Reionization. In 2018, the EDGES experiment caught the attention of the cosmology community with a potential detection of an absorption feature in the sky-averaged radio spectrum centred at 78 MHz. The feature is deeper than expected, and, if confirmed, would call for new physics. However, different groups have re-analyzed the EDGES data and questioned the reliability of the signal. The Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) is a sky-averaged 21-cm experiment aiming at improving the current observations by tackling the issues faced by current instruments related to residual systematic signals in the data. The novel experimental approach focuses on detecting and jointly explaining these systematics together with the foregrounds and the cosmological signal using Bayesian statistics. To achieve this, REACH features simultaneous observations with two different antennas, an ultra wideband system (redshift range 7.5 to 28), and a receiver calibrator based on in-field measurements. Simulated observations forecast percent-level constraints on astrophysical parameters, potentially opening up a new window to the infant Universe., Comment: Publication date: 21 July 2022

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process., Comment: 57 pages

In the pursuit of primordial non-Gaussianities, we hope to access smaller scales across larger comoving volumes. At low redshift, the search for primordial non-Gaussianities is hindered by gravitational collapse, to which we often associate a scale $k_{\rm NL}$. Beyond these scales, it will be hard to reconstruct the modes sensitive to the primordial distribution. When forecasting future constraints on the amplitude of primordial non-Gaussianity, $f_{\rm NL}$, off-diagonal components are usually neglected in the covariance because these are small compared to the diagonal. We show that the induced non-Gaussian off-diagonal components in the covariance degrade forecast constraints on primordial non-Gaussianity, even when all modes are well within what is usually considered the linear regime. As a testing ground, we examine the effects of these off-diagonal components on the constraining power of the matter bispectrum on $f_{\rm NL}$ as a function of $k_{\rm max}$ and redshift, confirming our results against N-body simulations out to redshift $z=10$. We then consider these effects on the hydrogen bispectrum as observed from a PUMA-like 21-cm intensity mapping survey at redshifts $2

Rayleigh scattering of the cosmic microwave background (CMB) by neutral hydrogen shortly after recombination leaves frequency-dependent imprints on intensity and polarization fluctuations. High signal-to-noise observations of CMB Rayleigh scattering would provide additional insight into the physics of recombination, including greater constraining power for parameters like the primordial helium fraction, the light relic density, and the sum of neutrino masses. However, such a measurement of CMB Rayleigh scattering is challenging due to the presence of astrophysical foregrounds, which are more intense at the high frequencies, where the effects of Rayleigh scattering are most prominent. Here we forecast the detectability of CMB Rayleigh scattering including foreground removal using blind internal linear combination methods for a set of near-future surveys. We show that atmospheric effects for ground-based observatories and astrophysical foregrounds pose a significant hindrance to detecting CMB Rayleigh scattering with experiments planned for this decade, though a high-significance measurement should be possible with a future CMB satellite., Comment: 14 pages, 12 figures

The Laser Interferometer Space Antenna (LISA) has two scientific objectives of cosmological focus: to probe the expansion rate of the universe, and to understand stochastic gravitational-wave backgrounds and their implications for early universe and particle physics, from the MeV to the Planck scale. However, the range of potential cosmological applications of gravitational wave observations extends well beyond these two objectives. This publication presents a summary of the state of the art in LISA cosmology, theory and methods, and identifies new opportunities to use gravitational wave observations by LISA to probe the universe.

Cosmic inflation provides a window to the highest energy densities accessible in nature, far beyond those achievable in any realistic terrestrial experiment. Theoretical insights into the inflationary era and its observational probes may therefore shed unique light on the physical laws underlying our universe. This white paper describes our current theoretical understanding of the inflationary era, with a focus on the statistical properties of primordial fluctuations. In particular, we survey observational targets for three important signatures of inflation: primordial gravitational waves, primordial non-Gaussianity and primordial features. With the requisite advancements in analysis techniques, the tremendous increase in the raw sensitivities of upcoming and planned surveys will translate to leaps in our understanding of the inflationary paradigm and could open new frontiers for cosmology and particle physics. The combination of future theoretical and observational developments therefore offer the potential for a dramatic discovery about the nature of cosmic acceleration in the very early universe and physics on the smallest scales., Comment: Contribution to Snowmass 2021; 103 pages, 8 figures, 378 endorsers; v2: minor changes

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan., Comment: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), 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 velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed 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., Comment: Contribution to Snowmass 2021. Note some text overlap with CMB-HD Astro2020 APC and RFI (arXiv:1906.10134, arXiv:2002.12714). Science case further broadened and updated

We investigate tomography of 21-cm brightness temperature fluctuations during the Dark Ages as a probe for constraining primordial non-Gaussianity. We expand the 21- cm brightness temperature up to cubic order in perturbation theory and improve previous models of the signal by including the effect of the free electron fraction. Using modified standard perturbation theory methods that include baryonic pressure effects we derive an improved secondary bispectrum and for the first time derive the secondary trispectrum of 21-cm brightness temperature fluctuations. We then forecast the amount of information available from the Dark Ages to constrain primordial non-Gaussianity, including the imprints of massive particle exchange during inflation and we determine how much signal is lost due to secondary non-Gaussianity. We find that although secondary non-Gaussianity swamps the primordial signal, primordial non-Gaussianity can still be extracted with signal-to-noise ratios that surpass current and future CMB experiments by several orders of magnitude, depending on the experimental setup. Furthermore, we conclude that for the bi- and trispectra of massive particle exchange marginalizing over other primordial shapes affects signal-to-noise ratios more severely than secondary shapes. Baryonic pressure effects turn out to have a negligible impact on our forecasts, even at scales close to the Jeans scale. The results of this work reinforce the prospects of 21-cm brightness temperature fluctuations from the Dark Ages as the ultimate probe for primordial non-Gaussianity., Comment: 47 pages, 20 figures, code available at https://github.com/tsfloss/21-cm-Dark-Ages-Fisher-Forecast. v2: minor typos corrected. v3: matches published version

We present a reproduction of the Planck 2018 angular power spectra at $\ell > 30$, and associated covariance matrices, for intensity and polarization maps at 100, 143 and 217 GHz. This uses a new, publicly available, pipeline that is part of the PSpipe package. As a test case we use the same input maps, ancillary products, and analysis choices as in the Planck 2018 analysis, and find that we can reproduce the spectra to 0.1$\sigma$ precision, and the covariance matrices to 10%. We show that cosmological parameters estimated from our re-derived products agree with the public Planck products to 0.1$\sigma$, providing an independent cross-check of the Planck team's analysis. Going forward, the publicly-available code can be easily adapted to use alternative input maps, data selections and analysis choices, for future optimal analysis of Planck data with new ground-based Cosmic Microwave Background data., Comment: 25 pages, 18 figures, code available at https://simonsobs.github.io/planck-pr3-web/

We introduce and validate a delensing framework for the Simons Observatory (SO), which will be used to improve constraints on inflationary gravitational waves (IGWs) by reducing the lensing noise in measurements of the $B$-modes in CMB polarization. SO will initially observe CMB by using three small aperture telescopes and one large-aperture telescope. While polarization maps from small-aperture telescopes will be used to constrain IGWs, the internal CMB lensing maps used to delens will be reconstructed from data from the large-aperture telescope. Since lensing maps obtained from the SO data will be noise-dominated on sub-degree scales, the SO lensing framework constructs a template for lensing-induced $B$-modes by combining internal CMB lensing maps with maps of the cosmic infrared background from Planck as well as galaxy density maps from the LSST survey. We construct a likelihood for constraining the tensor-to-scalar ratio $r$ that contains auto- and cross-spectra between observed $B$-modes and the lensing $B$-mode template. We test our delensing analysis pipeline on map-based simulations containing survey non-idealities, but that, for this initial exploration, does not include contamination from Galactic and extragalactic foregrounds. We find that the SO survey masking and inhomogeneous and atmospheric noise have very little impact on the delensing performance, and the $r$ constraint becomes $\sigma(r)\approx 0.0015$ which is close to that obtained from the idealized forecasts in the absence of the Galactic foreground and is nearly a factor of two tighter than without delensing. We also find that uncertainties in the external large-scale structure tracers used in our multi-tracer delensing pipeline lead to bias much smaller than the $1\,\sigma$ statistical uncertainties., Comment: 22 pages, 14 figures

Following the reported detection of an absorption profile associated with the 21~cm sky-averaged signal from the Cosmic Dawn by the EDGES experiment in 2018, a number of experiments have been set up to verify this result. This paper discusses the design process used for global 21~cm experiments, focusing specifically on the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH). This experiment will seek to understand and compensate for systematic errors present using detailed modelling and characterization of the instrumentation. There is detailed the quantitative figures of merit and numerical modelling used to assist the design process of the REACH dipole antenna (one of the 2 antenna designs for REACH Phase I). This design process produced a 2.5:1 frequency bandwidth dipole. The aim of this design was to balance spectral smoothness and low impedance reflections with the ability to describe and understand the antenna response to the sky signal to inform the critically important calibration during observation and data analysis., Comment: 32 pages, 30 figures, to be submitted to the Journal of Astronomical Instrumentation

In this work we analyse in detail the possibility of using small and intermediate-scale gravitational wave anisotropies to constrain the inflationary particle content. First, we develop a phenomenological approach focusing on anisotropies generated by primordial tensor-tensor-scalar and purely gravitational non-Gaussianities. We highlight the quantities that play a key role in determining the detectability of the signal. To amplify the power of anisotropies as a probe of early universe physics, we consider cross-correlations with CMB temperature anisotropies. We assess the size of the signal from inflationary interactions against so-called induced anisotropies. In order to arrive at realistic estimates, we obtain the projected constraints on the non-linear primordial parameter $F_{\rm NL}$ for several upcoming gravitational wave probes in the presence of the astrophysical gravitational wave background. We further illustrate our findings by considering a concrete inflationary realisation and use it to underscore a few subtleties in the phenomenological analysis., Comment: 47 pages, 16 figures

Cross-correlations between Cosmic Microwave Background (CMB) temperature and polarization anisotropies and $\mu$-spectral distortions have been considered to measure (squeezed) primordial scalar bispectra in a range of scales inaccessible to primary CMB bispectra. In this work we address whether it is possible to constrain tensor non-Gaussianities with these cross-correlations. We find that only primordial tensor bispectra with statistical anisotropies leave distinct signatures, while isotropic tensor bispectra leave either vanishing or highly suppressed signatures. We discuss how the angular dependence of squeezed bispectra in terms of the short and long momenta determine the non-zero cross-correlations. We also discuss how these non-vanishing configurations are affected by the way in which primordial bispectra transform under parity. By employing the so-called BipoSH formalism to capture the observational effects of statistical anisotropies, we make Fisher-forecasts to assess the detection prospects from $\mu T$, $\mu E$ and $\mu B$ cross-correlations. Observing statistical anisotropies in squeezed $\langle \gamma \gamma \gamma\rangle$ and $\langle \gamma \gamma \zeta\rangle$ bispectra is going to be challenging as the imprint of tensor perturbations on $\mu$-distortions is subdominant to scalar perturbations, therefore requiring a large, independent amplification of the effect of tensor perturbations in the $\mu$-epoch. In absence of such a mechanism, statistical anisotropies in squeezed $\langle \zeta \zeta \gamma\rangle$ bispectrum are the most relevant sources of $\mu T$, $\mu E$ and $\mu B$ cross-correlations. In particular, we point out that in anisotropic inflationary models where $\langle \zeta \zeta \zeta \rangle$ leaves potentially observable signatures in $\mu T$ and $\mu E$, the detection prospects of $\langle \zeta \zeta \gamma\rangle$ from $\mu B$ are enhanced., Comment: 36 pages, 6 figures. v3: Other typos corrected. Matched published version

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

We study the cosmic variance limit on constraining primordial non-Gaussianity for a variety of theory-motivated shapes. We consider general arguments for 2D and 3D surveys, with a particular emphasis on the CMB. A scale-invariant $N$-point correlator can be measured with a signal-to-noise that naively scales with the square root of the number of observed modes. This intuition generally fails for two reasons. First, the signal-to-noise scaling is reduced due to the blurring of the last scattering surface at short distances. This blurring is caused by the combination of projection and damping, but the loss of signal is not due to exponential decay, as both signal and noise are equally damped. Second, the behavior of the $N$-point correlator in the squeezed and collapsed (for $N>3$) limits can enhance the scaling of the signal-to-noise with the resolution, even with a reduced range of momenta probing these limits. We provide analytic estimates for all $N$-point correlators. We show that blurring affects equilateral-like shapes much more than squeezed ones. We discuss under what conditions the optimistic scalings in the collapsed limit can be exploited. Lastly, we confirm our analytical estimates with numerical calculations of the signal-to-noise for local, orthogonal and equilateral bispectra, and local trispectra. We also show that adding polarization to intensity data enhances the scaling for equilateral-like spectra., Comment: 39 pages, 6 figures. Added a paragraph on lensing effect. Matches journal version

We quantify the calibration requirements for systematic uncertainties for next-generation ground-based observatories targeting the large-angle $B$-mode polarization of the Cosmic Microwave Background, with a focus on the Simons Observatory (SO). We explore uncertainties on gain calibration, bandpass center frequencies, and polarization angles, including the frequency variation of the latter across the bandpass. We find that gain calibration and bandpass center frequencies must be known to percent levels or less to avoid biases on the tensor-to-scalar ratio $r$ on the order of $\Delta r\sim10^{-3}$, in line with previous findings. Polarization angles must be calibrated to the level of a few tenths of a degree, while their frequency variation between the edges of the band must be known to ${\cal O}(10)$ degrees. Given the tightness of these calibration requirements, we explore the level to which residual uncertainties on these systematics would affect the final constraints on $r$ if included in the data model and marginalized over. We find that the additional parameter freedom does not degrade the final constraints on $r$ significantly, broadening the error bar by ${\cal O}(10\%)$ at most. We validate these results by reanalyzing the latest publicly available data from the BICEP2/Keck collaboration within an extended parameter space covering both cosmological, foreground and systematic parameters. Finally, our results are discussed in light of the instrument design and calibration studies carried out within SO., Comment: 41 pages, 18 figures

Anisotropies in the cosmic microwave background (CMB) are primarily generated by Thomson scattering of photons by free electrons. Around recombination, the Thomson scattering probability quickly diminishes as the free electrons combine with protons to form neutral hydrogen off which CMB photons can scatter through Rayleigh scattering. Unlike Thomson scattering, Rayleigh scattering is frequency dependent resulting in the generation of anisotropies with a different spectral dependence. Unfortunately the Rayleigh scattering efficiency rapidly decreases with the expansion of the neutral universe, with the result that only a small percentage of photons are scattered by neutral hydrogen. Although the effect is very small, future CMB missions with higher sensitivity and improved frequency coverage are poised to measure Rayleigh scattering signal. The uncorrelated component of the Rayleigh anisotropies contains unique information on the primordial perturbations that could potentially be leveraged to expand our knowledge of the early universe. In this paper we explore whether measurements of Rayleigh scattering anisotropies can be used to constrain primordial non-Gaussianity (NG) and examine the hints of anomalies found by WMAP and \textit{Planck} satellites. We show that the additional Rayleigh information has the potential to improve primordial NG constraints by $30\%$, or more. Primordial bispectra that are not of the local type benefit the most from these additional scatterings, which we attribute to the different scale dependence of the Rayleigh anisotropies. Unfortunately this different scaling means that Rayleigh measurements can not be used to constrain anomalies or features on large scales. On the other hand, anomalies that may persist to smaller scales, such as the potential power asymmetry seen in WMAP and \textit{Planck}, could be improved by the addition of Rayleigh measurements., Comment: 12 pages, 5 figures

CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5\sigma$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL., Comment: 24 pages, 8 figures, 9 tables, submitted to ApJ. arXiv admin note: text overlap with arXiv:1907.04473

The cosmic microwave background (CMB) has been a treasure trove for cosmology. Over the next decade, current and planned CMB experiments are expected to exhaust nearly all primary CMB information. To further constrain cosmological models, there is a great benefit to measuring signals beyond the primary modes. Rayleigh scattering of the CMB is one source of additional cosmological information. It is caused by the additional scattering of CMB photons by neutral species formed during recombination and exhibits a strong and unique frequency scaling ($\propto \nu^4$). We will show that with sufficient sensitivity across frequency channels, the Rayleigh scattering signal should not only be detectable but can significantly improve constraining power for cosmological parameters, with limited or no additional modifications to planned experiments. We will provide heuristic explanations for why certain cosmological parameters benefit from measurement of the Rayleigh scattering signal, and confirm these intuitions using the Fisher formalism. In particular, observation of Rayleigh scattering allows significant improvements on measurements of $N_{\rm eff}$ and $\sum m_\nu$., Comment: 20 pages, 10 figures

Large-scale Fourier modes of the cosmic density field are of great value for learning about cosmology because of their well-understood relationship to fluctuations in the early universe. However, cosmic variance generally limits the statistical precision that can be achieved when constraining model parameters using these modes as measured in galaxy surveys, and moreover, these modes are sometimes inaccessible due to observational systematics or foregrounds. For some applications, both limitations can be circumvented by reconstructing large-scale modes using the correlations they induce between smaller-scale modes of an observed tracer (such as galaxy positions). In this paper, we further develop a formalism for this reconstruction, using a quadratic estimator similar to the one used for lensing of the cosmic microwave background. We incorporate nonlinearities from gravity, nonlinear biasing, and local-type primordial non-Gaussianity, and verify that the estimator gives the expected results when applied to N-body simulations. We then carry out forecasts for several upcoming surveys, demonstrating that, when reconstructed modes are included alongside directly-observed tracer density modes, constraints on local primordial non-Gaussianity are generically tightened by tens of percents compared to standard single-tracer analyses. In certain cases, these improvements arise from cosmic variance cancellation, with reconstructed modes taking the place of modes of a separate tracer, thus enabling an effective "multitracer" approach with single-tracer observations., Comment: 30 pages plus 14 pages appendices, 19 figures

The search for primordial gravitational waves through the $B$-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary $B$-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper, we revise estimates of the reionization-induced $B$-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-$\alpha$ forest. To derive these $B$-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Ly$\alpha$ data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles $\ell\approx 100$, reionization induces $B$-mode power with $\ell(\ell+1)C_\ell^{BB}/2\pi\approx 4\times 10^{-6}\,\mu$K$^2$. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio $r<10^{-4}$, and can increase by a factor of $\sim 50$ if reionization is sourced by highly clustered sources residing in haloes with mass of $\sim 10^{11}$ M$_\odot$. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments., Comment: Published in JCAP

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., Comment: Response to request for information (RFI) by the Panel of Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Astro2020 Decadal Survey regarding the CMB-HD APC (arXiv:1906.10134). Note some text overlap with original APC. Note also detector count and cost have been reduced by 1/3, and observing time increased by 1/3 compared to original APC; science goals expanded

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project., Comment: 46 pages, 16 figures, 7 tables; response to the request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Astro2020 Decadal Survey regarding PUMA APC submission (arXiv:1907.12559); v2: updated with correct bbl file

In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a "science-first" approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics., Comment: 22 pages, 1 figure, to appear in General Relativity and Gravitation

The next generation of ground-based CMB experiments aim to measure temperature and polarization fluctuations up to $\ell_{\rm max} \approx 5000$ over half of the sky. Combined with Planck data on large scales, this will provide improved constraints on primordial non-Gaussianity. However, the impressive resolution of these experiments will come at a price. Besides signal confusion from galactic foregrounds, extra-galactic foregrounds and late-time gravitational effects, gravitational lensing will introduce large non-Gaussianity that can become the leading contribution to the bispectrum covariance through the connected 4-point function. Here, we compute this effect analytically for the first time on the full sky for both temperature and polarization. We compare our analytical results with those obtained directly from map-based simulations of the CMB sky for several levels of instrumental noise. Of the standard shapes considered in the literature, the local shape is most affected, resulting in a 35\% increase of the estimator standard deviation for an experiment like the Simons Observatory (SO) and a 110\% increase for a cosmic-variance limited experiment, including both temperature and polarization modes up to $\ell_{\rm max} = 3800$. Because of the nature of the lensing 4-point function, the impact on other shapes is reduced while still non negligible for the orthogonal shape. Two possible avenues to reduce the non-Gaussian contribution to the covariance are proposed. First by marginalizing over lensing contributions, such as the ISW-lensing 3pt function in temperature, and second by delensing the CMB. We show the latter method can remove almost all extra covariance, reducing the effect to below $<$5\% for local bispectra. At the same time, delensing would remove signal biases from secondaries induced by lensing, such as ISW-lensing., Comment: 22 pages, 9 figures. Abstract abridged

Upcoming cosmic microwave background (CMB) data can be used to explore harmonic 3-point functions that involve the B-mode component of the CMB polarization signal. We focus on bispectra describing the non-Gaussian correlation of the B-mode field and the CMB temperature anisotropies (T) and/or E-mode polarization, i.e. , , and . Such bispectra probe violations of the tensor consistency relation: the model-independent behavior of cosmological correlation functions that involve a large-wavelength tensor mode (gravitational wave). An observed violation of the tensor consistency relation would exclude a large number of inflation models. We describe a generalization of the Komatsu-Spergel-Wandelt (KSW) bispectrum estimator that allows statistical inference on this type of primordial non-Gaussianity with data of the CMB temperature and polarization anisotropies. The generalized estimator shares its statistical properties with the existing KSW estimator and retains the favorable numerical scaling with angular resolution. In this paper we derive the estimator and present a set of Fisher forecasts. We show how the forecasts scale with various experimental parameters such as lower and upper angular band-limit, relevant for e.g. the upcoming ground-based Simons Observatory experiment and proposed LiteBIRD satellite experiment. We comment on possible contaminants due to secondary cosmological and astrophysical sources., Comment: 40 pages, 7 figures

Many models of the early universe predict that there should be primordial tensor perturbations. These leave an imprint into the temperature and polarisation anisotropies of the cosmic microwave background (CMB). The differential equation describing the primordial tensor perturbations is a second order differential equation and thus has two solutions. Canonically, the decaying solution of this equation in radiation domination is dropped as it diverges at early times and on superhorizon scales while it is then suppressed at late times. Furthermore, if there is an inflationary phase prior to the radiation domination phase, the amplitude of the decaying mode will also be highly suppressed as it enters the radiation phase, thus its effect will be negligible. In this study we remain agnostic to the early universe models describing pre-radiation domination physics and allow this mode to be present and see what effect it has on the CMB anisotropies. We find that the decaying mode, if normalised at the same time on subhorizon scales as the growing mode leaves an imprint on the CMB anisotropies that is identical to the growing mode. Contrary to expectation, on large scales both modes are poorly constrained for a scale invariant spectrum, and the apparent divergence of the decaying mode does not lead to a divergent physical observable. Quantitatively, the decaying mode can be more constrained both from temperature and polarisation anisotropies. We use a model independent, non-parametric, approach to constrain both of these primordial tensor perturbations using the temperature and polarisation anisotropies. We find that both modes are best constrained at the reionisation and recombination bumps and crucially, at the reionisation bump the decaying mode can be distinguished from the growing mode., Comment: 10 pages, comments welcome!

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

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged], Comment: 24 pages, 8 figures, 3 tables; White paper submitted to ESA Voyage 2050

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., Comment: Project White Paper submitted to the 2020 Decadal Survey, 10 pages plus references. arXiv admin note: substantial text overlap with arXiv:1907.04473

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively., Comment: 10 pages + references, 3 figures, 3 tables; project white paper submitted to the Astro2020 decadal survey; further details in updated arXiv:1810.09572

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., Comment: Astro2020 Decadal Project Whitepaper. arXiv admin note: text overlap with arXiv:1808.07445

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

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., Comment: APC White Paper for the Astro2020 Decadal, with updated proposing team

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged)., Comment: 50 pages, 20 figures