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.

The cosmic microwave background (CMB) is sensitive to the recent phase of accelerated cosmic expansion through the late-time integrated Sachs-Wolfe (ISW) effect, which manifests as secondary temperature fluctuations on large angular scales. However, the large cosmic variance from primary CMB fluctuations limits the usefulness of this effect in constraining dark energy or modified gravity. In this paper, we propose a novel method to separate the ISW signal from the primary signal using gravitational lensing, based on the fact that the ISW signal is, to a good approximation, not gravitationally lensed. We forecast how well we can isolate the ISW signal for different experimental configurations, and discuss various applications, including modified gravity, large-scale CMB anomalies, and measurements of local-type primordial non-Gaussianity. Although not within reach of current experiments, the proposed method is a unique way to remove the cosmic variance of the primary signal, allowing for better CMB-based constraints on late-time phenomena than previously thought possible.

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.

It is predicted that the spectrum of radio emission from the whole sky should show a dip arising from the action of the light from the first stars on the hydrogen atoms in the surrounding gas, which causes the 21-cm line to appear in absorption against the cosmic microwave background. Bowman et al. 2018 identified a broad flat-bottomed absorption profile centred at 78 MHz, which could be this feature, although the depth of the profile is much larger than expected. We have examined the modelling process they used and find that their data implies unphysical parameters for the foreground emission and also that their solution is not unique in the sense that we found other simple formulations for the signal that are different in shape but that also fit their data. We argue that this calls into question the interpretation of these data as an unambiguous detection of the cosmological 21-cm absorption signature.

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.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

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. <TTB>, <EEB>, and <TEB>. 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.

In the standard big bang model the universe starts in a radiation dominated era, where the gravitational perturbations are described by second order differential equations, which will generally have two orthogonal set of solutions. One is the so called {\it growing(cosine)} mode and the other is the {\it decaying(sine)} mode, where the nomenclature is derived from their behaviour on super-horizon(sub-horizon) scales. The decaying mode is qualitatively different to the growing mode of adiabatic perturbations as it evolves with time on \emph{super-horizon} scales. The time dependence of this mode on super-horizon scales is analysed in both the synchronous gauge and the Newtonian gauge to understand the true gauge invariant behaviour of these modes. We then explore constraints on the amplitude of this mode on scales between $k \sim 10^{-5}$ Mpc$^{-1}$ and $k \sim 10^{-1}$ Mpc$^{-1}$ using the temperature and polarization anisotropies from the cosmic microwave background, by computing the Fisher information. Binning the primordial power non-parametrically into 100 bins, we find that the decaying modes are constrained at comparable variance as the growing modes on scales smaller than the horizon today using temperature anisotropies. Adding polrisation data makes the decaying mode more constrained. The decaying mode amplitude is thus constrained by $\sim 1/l$ of the growing mode. On super-horizon scales, the growing mode is poorly constrained, while the decaying mode cannot substantially exceed the scale-invariant amplitude. This interpretation differs substantially from the past literature, where the constraints were quoted in gauge-dependent variables, and resulted in illusionary tight super-horizon decaying mode constraints. The results presented here can generally be used to non-parametrically constrain any model of the early universe.

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.

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.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

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.

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

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

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

Inflation may provide unique insight into the physics at the highest available energy scales that cannot be replicated in any realistic terrestrial experiment. Features in the primordial power spectrum are generically predicted in a wide class of models of inflation and its alternatives, and are observationally one of the most overlooked channels for finding evidence for non-minimal inflationary models. Constraints from observations of the cosmic microwave background cover the widest range of feature frequencies, but the most sensitive constraints will come from future large-scale structure surveys that can measure the largest number of linear and quasi-linear modes.

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

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

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

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 investigate the detection prospects for gravitational lensing of three-dimensional maps from upcoming line intensity surveys, focusing in particular on the impact of gravitational nonlinearities on standard quadratic lensing estimators. Using perturbation theory, we show that these nonlinearities can provide a significant contaminant to lensing reconstruction, even for observations at reionization-era redshifts. However, we show how this contamination can be mitigated with the use of a "bias-hardened" estimator. Along the way, we present an estimator for reconstructing long-wavelength density modes, in the spirit of the "tidal reconstruction" technique that has been proposed elsewhere, and discuss the dominant biases on this estimator. After applying bias-hardening, we find that a detection of the lensing potential power spectrum will still be challenging for the first phase of SKA-Low, CHIME, and HIRAX, with gravitational nonlinearities decreasing the signal to noise by a factor of a few compared to forecasts that ignore these effects. On the other hand, cross-correlations between lensing and galaxy clustering or cosmic shear from a large photometric survey look promising, provided that systematics can be sufficiently controlled. We reach similar conclusions for a single-dish survey inspired by CII measurements planned for CCAT-prime, suggesting that lensing is an interesting science target not just for 21cm surveys, but also for intensity maps of other lines.

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

Cosmic microwave background (CMB) constraints on dark matter annihilation are a uniquely powerful tool in the quest to understand the nature of dark matter. Annihilation of dark matter to Standard Model particles between recombination and reionization heats baryons, ionizes neutral hydrogen, and alters the CMB visibility function. Surprisingly, CMB bounds on dark matter annihilation are not expected to improve significantly with the dramatic improvements in sensitivity expected in future cosmological surveys. In this paper, we will present a simple physical description of the origin of the CMB constraints and explain why they are nearly saturated by current observations. The essential feature is that dark matter annihilation primarily affects the ionization fraction which can only increase substantially at times when the universe was neutral. The resulting change to the CMB occurs on large angular scales and leads to a phenomenology similar to that of the optical depth to reionization. We will demonstrate this impact on the CMB both analytically and numerically. Finally, we will discuss the additional impact that changing the ionization fraction has on large scale structure.

The observed dipole anisotropy of the cosmic microwave background (CMB) temperature is much larger than the fluctuations observed on smaller scales and is dominated by the kinematic contribution from the Doppler shifting of the monopole due to our motion with respect to the CMB rest frame. In addition to this kinematic component, there is expected to be an intrinsic contribution with an amplitude about two orders of magnitude smaller. Here we explore a method whereby the intrinsic CMB dipole can be reconstructed through observation of temperature fluctuations on small scales which result from gravitational lensing. Though the experimental requirements pose practical challenges, we show that one can in principle achieve a cosmic variance limited measurement of the primary dipole using the reconstruction method we describe. Since the primary CMB dipole is sensitive to the largest observable scales, such a measurement would have a number of interesting applications for early universe physics, including testing large-scale anomalies, extending the lever-arm for measuring local non-Gaussianity, and constraining isocurvature fluctuations on super-horizon scales.

The Thomson optical depth from reionization is a limiting factor in measuring the amplitude of primordial fluctuations, and hence in measuring physics that affects the low-redshift amplitude, such as the neutrino masses. Current constraints on the optical depth, based on directly measuring large-scale cosmic microwave background (CMB) polarization, are challenging due to foregrounds and systematic effects. Here, we consider an indirect measurement of large-scale polarization, using observed maps of small-scale polarization together with maps of fields that distort the CMB, such as CMB lensing and patchy reionization. We find that very futuristic CMB surveys will be able to reconstruct large-scale polarization, and thus the mean optical depth, using only measurements on small scales.

It is generally expected that heavy fields are present during inflation, which can leave their imprint in late-time cosmological observables. The main signature of these fields is a small amount of distinctly shaped non-Gaussianity, which if detected, would provide a wealth of information about the particle spectrum of the inflationary Universe. Here we investigate to what extent these signatures can be detected or constrained using futuristic 21-cm surveys. We construct model-independent templates that extract the squeezed-limit behavior of the bispectrum, and examine their overlap with standard inflationary shapes and secondary non-Gaussianities. We then use these templates to forecast detection thresholds for different masses and couplings using a 3D reconstruction of modes during the dark ages ($z\sim 30-100$). We consider interactions of several broad classes of models and quantify their detectability as a function of the baseline of a dark ages interferometer. Our analysis shows that there exists the tantalizing possibility of discovering new particles with different masses and interactions with future 21-cm surveys.

Primordial gravitational waves leave a characteristic imprint on the cosmic microwave background (CMB) in the form of $B$-mode polarization. Photons are also deflected by large scale gravitational waves which intervene between the source screen and our telescopes, resulting in curl-type gravitational lensing. Gravitational waves present at the epoch of reionization contribute to both effects, thereby leading to a non-vanishing cross-correlation between $B$-mode polarization and curl lensing of the CMB. Observing such a cross correlation would be very strong evidence that an observation of $B$-mode polarization was due to the presence of large scale gravitational waves, as opposed to astrophysical foregrounds or experimental systematic effects. We study the cross-correlation across a wide range of source redshifts and show that a post-SKA experiment aimed to map out the 21-cm sky between $15 \leq z \leq 30$ could rule out non-zero cross-correlation at high significance for $r \geq 0.01$.

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.

I will discuss to what degree the cosmic microwave background (CMB) can be used to constrain primordial non-Gaussianity involving one tensor and two scalar fluctuations, focusing on the correlation of one $B$-mode polarization fluctuation with two temperature fluctuations (BTT). In the simplest models of inflation, the tensor-scalar-scalar primordial bispectrum is non-vanishing and is of the same order in slow-roll parameters as the scalar-scalar-scalar bispectrum. I will show that constraints from an experiment like CMB-Stage IV using this observable are more than an order of magnitude better than those on the same primordial coupling obtained from temperature measurements alone. I will argue that $B$-mode non-Gaussianity opens up an as-yet-unexplored window into the early Universe, demonstrating that significant information on primordial physics remains to be harvested from CMB anisotropies. BTT presents a measure of both primordial tensors and primordial non-Gaussianity, two of the most sought after signatures of the inflationary paradigm.

Detecting a deviation from a featureless primordial power spectrum of fluctuations would give profound insight into the physics of the primordial Universe. Depending on their nature, primordial features can either provide direct evidence for the inflation scenario or pin down details of the inflation model. Thus far, using the cosmic microwave background (CMB) we have only been able to put stringent constraints on the amplitude of features, but no significant evidence has been found for such signals. Here we explore the limit of the experimental reach in constraining such features using 21 cm tomography at high redshift. A measurement of the 21 cm power spectrum from the Dark Ages is generally considered as the ideal experiment for early Universe physics, with potentially access to a large number of modes. We consider three different categories of theoretically motivated models: the sharp feature models, resonance models, and standard clock models. We study the improvements on bounds on features as a function of the total number of observed modes and identify parameter degeneracies. The detectability depends critically on the amplitude, frequency and scale-location of the features, as well as the angular and redshift resolution of the experiment. We quantify these effects by considering different fiducial models. Our forecast shows that a cosmic variance limited 21 cm experiment measuring fluctuations in the redshift range $30\leq z \leq 100$ with a 0.01-MHz bandwidth and sub-arcminute angular resolution could potentially improve bounds by several orders of magnitude for most features compared to current Planck bounds. At the same time, 21 cm tomography also opens up a unique window into features that are located on very small scales.

We consider the presence of oscillations in the primordial bispectrum, inspired by three different cosmological models; features in the primordial potential, resonant type non-Gaussianities and deviation from the standard Bunch Davies vacuum. In order to put constraints on their bispectra, a logical first step is to put these into factorized form which can be achieved via the recently proposed method of polynomial basis expansion on the tetrahedral domain. We investigate the viability of such an expansion for the oscillatory bispectra and find that one needs an increasing number of orthonormal mode functions to achieve significant correlation between the expansion and the original spectrum as a function of their frequency. To reduce the number of modes required, we propose a basis consisting of Fourier functions orthonormalized on the tetrahedral domain. We show that the use of Fourier mode functions instead of polynomial mode functions can lead to the necessary factorizability with the use of only 1/5 of the total number of modes required to reconstruct the bispectra with polynomial mode functions. Moreover, from an observational perspective, the expansion has unique signatures depending on the orientation of the oscillation due to a resonance effect between the mode functions and the original spectrum. This effect opens the possibility to extract informa- tion about both the frequency of the bispectrum as well as its shape while considering only a limited number of modes. The resonance effect is independent of the phase of the reconstructed bispectrum suggesting Fourier mode extraction could be an efficient way to detect oscillatory bispectra in the data.

We study the degree to which the cosmic microwave background (CMB) can be used to constrain primordial non-Gaussianity involving one tensor and two scalar fluctuations, focusing on the correlation of one polarization $B$ mode with two temperature modes. In the simplest models of inflation, the tensor-scalar-scalar primordial bispectrum is non-vanishing and is of the same order in slow-roll parameters as the scalar-scalar-scalar bispectrum. We calculate the $\langle BTT\rangle$ correlation arising from a primordial tensor-scalar-scalar bispectrum, and show that constraints from an experiment like CMB-Stage IV using this observable are more than an order of magnitude better than those on the same primordial coupling obtained from temperature measurements alone. We argue that $B$-mode non-Gaussianity opens up an as-yet-unexplored window into the early Universe, demonstrating that significant information on primordial physics remains to be harvested from CMB anisotropies.

We develop the tools necessary to assess the statistical significance of resonant features in the CMB correlation functions, combining power spectrum and bispectrum measurements. This significance is typically addressed by running a large number of simulations to derive the probability density function (PDF) of the feature-amplitude in the Gaussian case. Although these simulations are tractable for the power spectrum, for the bispectrum they require significant computational resources. We show that, by assuming that the PDF is given by a multi-variate Gaussian where the covariance is determined by the Fisher matrix of the sine and cosine terms, we can efficiently produce spectra that are statistically close to those derived from full simulations. By drawing a large number of spectra from this PDF, both for the power spectrum and the bispectrum, we can quickly determine the statistical significance of candidate signatures in the CMB, considering both single frequency and multi-frequency estimators. We show that for resonance models, cosmology and foreground parameters have little influence on the estimated amplitude, which allows to simplify the analysis considerably. A more precise likelihood treatment can then be applied to candidate signatures only. We also discuss a modal expansion approach for the power spectrum, aimed at quickly scanning through large families of oscillating models.

We propose optimal estimators for bispectra from excited states. Two common properties of such bispectra are the enhancement in the collinear limit, and the prediction of oscillating features. We review the physics behind excited states and some of the choices made in the literature. We show that the enfolded template is a good template in the collinear limit, but does poorly elsewhere, establishing a strong case for an improved estimator. Although the detailed scale dependence of the bispectra differs depending on various assumptions, generally the predicted bispectra are either effectively 1 or 2-dimensional and a simple Fourier basis suffices for accurate reconstruction. For an optimal CMB data analysis, combining all n-point functions, the choice for the excited state needs to be the same when computing power spectrum, bispectrum and higher order correlation functions. This has not always been the case, which could lead to wrong conclusions. We calculate the bispectrum for different choices previously discussed for the power spectrum, setting up a consistent framework to search for evidence of excited states in the CMB data.

We present the first attempt to use a combination of CMB, LIGO, and PPTA data to constrain both the tilt and the running of primordial tensor power spectrum through constraints on the gravitational wave energy density generated in the early universe. Combining measurements at different cosmological scales highlights how complementary data can be used to test the predictions of early universe models including the inflationary consistency relation. Current data prefers a slightly positive tilt ($n_t = 0.06^{+0.63}_{-0.89}$) and a negative running ($n_{t, {\rm run}} < -0.22$) for the tensor power spectrum spectrum. Interestingly, the addition of direct gravitational wave detector data alone puts strong bounds on the tensor-to-scalar ratio $r < 0.2$ since the large positive tensor tilt preferred by the \textit{Planck} temperature power spectrum is no longer allowed. Adding the recently released BICEP2/KECK and Planck 353 GHz Polarization cross correlation data puts an even stronger bound $r<0.1$ . We comment on possible effects of a large positive tilt on the background expansion and show that depending on the assumptions regarding the UV cutoff ($k_{\rm UV}/k_* = 10^{24}$) of the primordial spectrum of gravitational waves, the strongest bounds on $n_t = 0.07^{+0.52}_{-0.80}$ are derived from this effect.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

We propose an (optimal) estimator for a CMB bispectrum containing logarithmically spaced oscillations. There is tremendous theoretical interest in such bispectra, and they are predicted by a plethora of models, including axion monodromy models of inflation and initial state modifications. The number of resolved logarithmical oscillations in the bispectrum is limited due to the discrete resolution of the multipole bispectrum. We derive a simple relation between the maximum number of resolved oscillations and the frequency. We investigate several ways to factorize the primordial bispectrum, and conclude that a one dimensional expansion in the sum of the momenta $\sum k_i = k_t$ is the most efficient and flexible approach. We compare the expansion to the exact result in multipole space and show for $\omega_{\rm eff}=100$ that $\mathcal{O}(10^3)$ modes are sufficient for an accurate reconstruction. We compute the expected $\sigma_{f_{\rm NL}}$ and find that within an effective field theory (EFT) the overall signal to noise scales as $S/N\propto \omega^{3/2}$. Using only the temperature data we find $S/N\sim\mathcal{O}(1-10^2)$ for the frequency domain set by the EFT.

There exists some tension on large scales between the Planck data and the LCDM concordance model of the Universe, which has been amplified by the recently claimed discovery of non-zero tensor to scalar ratio $r$. At the same time, the current best-fit value of $r$ suggests large field inflation delta phi>M_p, which requires a UV complete description of inflation. A very promising working example that predicts large tensor modes and can be UV completed is axion monodromy inflation. This realization of inflation naturally produces oscillating features, as consequence of a broken shift symmetry. We analyse a combination of Planck, ACT, SPT, WMAP low l polarization and BICEP2 data, and show a long wavelength feature from a periodic potential can alleviate the tension at low multipoles with an improvement delta chi^2 ~2.5-4 per degree of freedom, depending on the level of foreground subtraction. As with an introduction of running, one expects that any scale dependence should lead to a worsened fit at high multipoles. We show that the logarithmic nature of the axion feature in combination with a tilt n_s~1 allows the fit to be identical to a no feature model at the 2 percent level on scales 100 =< l =<3500, and quite remarkable actually slightly improves the fit at scales l >2000. We also consider possible unremoved dust foregrounds and show that including these hardly changes the best-fit parameters. Our analysis suggests an axion decay constant of f/M_p ~O(.01). We discuss how Planck measurements of the TE and EE spectra can further constrain axion monodromy inflation with such a large feature. A measurement of the large scale structure power spectrum is even more promising, as the effect is much bigger since the tensor modes do not affect the large scales. At the same time, a feature could also lead to a lower sigma8, lifting the tension between CMB and SZ constraints on sigma8.

We apply our recently developed code to search for resonance features in the Planck CMB temperature data. We search both for log spaced oscillations or linear spaced oscillations and compare our findings with results of our WMAP9 analysis and the Planck team analysis. While there are hints of log spaced resonant features present in the WMAP9 data, the significance of these features weaken with more data. With more accurate small scale measurements, we also find that the best fit frequency has shifted and the amplitude has been reduced. We confirm the presence of a several low frequency peaks, earlier identified by the Planck team, but with a better improvement of fit (delta chi^2 ~ 12). We further investigate this improvement by allowing the lensing potential to vary as well, showing mild correlation between the amplitude of the oscillations and the lensing amplitude. We find that the improvement of the fit increases even more (delta chi^2 ~ 14) for the low frequencies that modify the spectrum in a way that mimics the lensing effect. Since these features were not present in the WMAP data, they are primarily due to better measurements of Planck at small angular scales. For linear spaced oscillations we find a maximum delta chi^2 ~ 13 scanning two orders of magnitude in frequency space, and the biggest improvements are at extremely high frequencies. We recover a best fit frequency very close to the one found in WMAP9, which confirms that the fit improvement is driven by low l. Further comparisons with WMAP9 show Planck contains many more features, both for linear and log space oscillations, but with a smaller improvement of fit. We discuss the improvement as a function of the number of modes and study the effect of the 217 GHz map, which appears to drive most of the improvement for log spaced oscillations. We conclude that none of the detected features are statistically significant.

A small deviation from scale invariance in the form of oscillations in the primordial correlation spectra has been predicted by various cosmological models. In this paper we review a recently developed method to search for these features in the data in a more effective way. By Taylor expanding the small features around the 'background' cosmology, we have shown we are able to improve the search for these features compared to previous analyses. In this short paper we will extend that work by combining this method with a multi nested sampler. We recover our previous findings and are able to do so in 192 CPU hours. We will also briefly discuss the possibility of a long wavelength feature in the data to alleviate some tension between CMB data and the LCDM+r concordance cosmology.

Fluctuations of the 21 cm brightness temperature before the formation of the first stars hold the promise of becoming a high-precision cosmological probe in the future. The growth of over densities is very well described by perturbation theory at that epoch and the signal can in principle be predicted to arbitrary accuracy for given cosmological parameters. Recently, Tseliakhovich and Hirata pointed out a previously neglected and important physical effect, due to the fact that baryons and cold dark matter (CDM) have supersonic relative velocities after recombination. This relative velocity suppresses the growth of matter fluctuations on scales $k \sim 10-10^3$ Mpc$^{-1}$. In addition, the amplitude of the small-scale power spectrum is modulated on the large scales over which the relative velocity varies, corresponding to $k \sim 0.005-1$ Mpc$^{-1}$. In this paper, the effect of the relative velocity on 21 cm brightness temperature fluctuations from redshifts $z \geq 30$ is computed. We show that the 21 cm power spectrum is affected on ${\it most}$ scales. On small scales, the signal is typically suppressed several tens of percent, except for extremely small scales ($k \gtrsim 2000$ Mpc$^{-1}$) for which the fluctuations are boosted by resonant excitation of acoustic waves. On large scales, 21 cm fluctuations are enhanced due to the non-linear dependence of the brightness temperature on the underlying gas density and temperature. The enhancement of the 21 cm power spectrum is of a few percent at $k \sim 0.1$ Mpc$^{-1}$ and up to tens of percent at $k \lesssim 0.005$ Mpc$^{-1}$, for standard $\Lambda$CDM cosmology. In principle this effect allows to probe the small-scale matter power spectrum not only through a measurement of small angular scales but also through its effect on large angular scales.

In this first of two papers, we present a new method for searching for oscillatory features in the primordial power spectrum. A wide variety of models predict these features in one of two different flavors: logarithmically spaced oscillations and linearly spaced oscillations. The proposed method treats the oscillations as perturbations on top of the scale-invariant power spectrum, allowing us to vary all cosmological parameters. This perturbative approach reduces the computational requirements for the search as the transfer functions and their derivatives can be precomputed. We show that the most significant degeneracy in the analysis is between the distance to last scattering and the overall amplitude at low frequencies. For models with logarithmic oscillations, this degeneracy leads to an uncertainty in the phase. For linear spaced oscillations, it affects the frequency of the oscillations. In this first of two papers, we test our code on simulated Planck-like data, and show we are able to recover fiducial input oscillations with an amplitude of a few times order 10^{-2}. We apply the code to WMAP9-year data and confirm the existence of two intriguing resonant frequencies for log spaced oscillations. For linear spaced oscillations we find a single resonance peak. We use numerical simulations to assess the significance of these features and conclude that the data do not provide compelling evidence for the existence of oscillatory features in the primordial spectrum.

We consider the cross-correlation between free electrons and neutral hydrogen during the epoch of reionization. The free electrons are traced by the optical depth to reionization tau, while the neutral hydrogen can be observed through 21 cm photon emission. As expected, this correlation is sensitive to the detailed physics of reionization. Foremost, if reionization occurs through the merger of relatively large halos hosting an ionizing source, the free electrons and neutral hydrogen are anti-correlated for most of the reionization history. A positive contribution to the correlation can occur when the halos that can form an ionizing source are small. A measurement of this sign change in the cross-correlation would directly measure the bubble bias and as, such, the halo mass. We estimate the signal-to-noise of the cross-correlation using the estimator for inhomogeneous reionization tau_{lm} proposed by Dvorkin and Smith (2009). We find that with upcoming radio interferometers and CMB experiments, the cross-correlation is measurable going up to multipoles l ~ 1000. We also derive parameter constraints and conclude that, despite the foregrounds, the cross-correlation proofs a complementary measurement of the EoR parameters to the 21 cm and CMB polarization auto-correlations expected to be observed in the coming decade.

Models of axion inflation are particularly interesting since they provide a natural justification for the flatness of the potential over a super-Planckian distance, namely the approximate shift-symmetry of the inflaton. In addition, most of the observational consequences are directly related to this symmetry and hence are correlated. Large tensor modes can be accompanied by the observable effects of a the shift-symmetric coupling $\phi F\tilde F$ to a gauge field. During inflation this coupling leads to a copious production of gauge quanta and consequently a very distinct modification of the primordial curvature perturbations. In this work we compare these predictions with observations. We find that the leading constraint on the model comes from the CMB power spectrum when considering both WMAP 7-year and ACT data. The bispectrum generated by the non-Gaussian inverse-decay of the gauge field leads to a comparable but slightly weaker constraint. There is also a constraint from mu-distortion using TRIS plus COBE/FIRAS data, but it is much weaker. Finally we comment on a generalization of the model to massive gauge fields. When the mass is generated by some light Higgs field, observably large local non-Gaussianity can be produced.

We use the WMAP 7 data to place constraints on oscillations supplementing an almost scale-invariant primordial power spectrum. Such oscillations are predicted by a variety of models, some of which amount to assuming there is some non-trivial choice of the vacuum state at the onset of inflation. In this paper we will explore data-driven constraints on two distinct models of initial state modifications. In both models the frequency, phase and amplitude are degrees of freedom of the theory for which the theoretical bounds are rather weak: both the amplitude and frequency have allowed values ranging over several orders of magnitude. This requires many computationally expensive evaluations of the model CMB spectra and their goodness-of-fit, even in a Markov Chain Monte Carlo (MCMC), normally the most efficient fitting method for such a problem. To search more efficiently we first run a densely spaced grid, with only 3 varying parameters; the frequency, the amplitude and the baryon density. We obtain the optimal frequency and run an MCMC at the best fit frequency, randomly varying all other relevant parameters. To reduce the computational time of each power spectrum computation, we adjust both comoving momentum integration and spline interpolation (in l) as a function of frequency and amplitude of the primordial power spectrum. Applying this to the WMAP 7 data allows us to improve existing constraints on the presence of oscillations. We confirm earlier findings that certain frequencies can improve the fitting over a model without oscillations. For those frequencies we compute the posterior probability, allowing us to put some constraints on the primordial parameter space of both models.

There exist several models of inflation that produce primordial bispectra that contain a large number of oscillations. In this paper we discuss these models, and aim at finding a method of detecting such bispectra in the data. We explain how the recently proposed method of mode expansion of bispectra might be able to reconstruct these spectra from separable basis functions. Extracting these basis functions from the data might then lead to observational constraints on these models.

We present high-spectral-resolution, optical spectra of the Herbig Be star MWC 147, in which we spectrally resolve several emission lines, including the [O I] lines at 6300 and 6363\deg. Their highly symmetric, double-peaked line profiles indicate that the emission originates in a rotating circumstellar disk. We deconvolve the Doppler-broadened [O I] emission lines to obtain a measure of emission as a function of distance from the central star. The resulting radial surface brightness profiles are in agreement with a disk structure consisting of a flat, inner, gaseous disk and a flared, outer, dust disk. The transition between these components at 2 to 3 AU corresponds to the estimated dust sublimation radius. The width of the double-peaked Mg II line at 4481\deg suggests that the inner disk extends to at least 0.10 AU, close to the corotation radius.

In this short note we translate the best available observational bounds on the CMB bispectrum amplitudes into constraints on a specific scale-invariant New Physics Hypersurface (NPH) model of vacuum state modifications, as first proposed by Danielsson, in general models of single-field inflation. As compared to the power spectrum the bispectrum constraints are less ambiguous and provide an interesting upper bound on the cut-off scale in general models of single-field inflation with a small speed of sound. This upper bound is incompatible with the power spectrum constraint for most of the parameter domain, leaving very little room for minimal cut-off vacuum state modifications in general single-field models with a small speed of sound.

Deviations from the Bunch-Davies vacuum during an inflationary period can leave a testable imprint on the higher-order correlations of the CMB and large scale structures in the Universe. The effect is particularly pronounced if the statistical non-Gaussianity is inherently large, such as in models of inflation with a small speed of sound, e.g. DBI. First reviewing the motivations for a modified vacuum, we calculate the non-Gaussianity for a general action with a small speed of sound. The shape of its bispectrum is found to most resemble the 'orthogonal' or 'local' templates depending on the phase of the Bogolyubov parameter. In particular, for DBI models of inflation the bispectrum can have a profound 'local' template feature, in contrast to previous results. Determining the projection into the observational templates allows us to derive constraints on the absolute value of the Bogolyubov parameter. In the small sound speed limit, the derived constraints are generally stronger than the existing constraint derived from the power spectrum. The bound on the absolute value of the Bogolyubov parameter ranges from the 10^-6 to the 10^-3 level for H/\Lambda_c = 10^-3, depending on the specific details of the model, the sound speed and the phase of the Bogolyubov parameter.

Modifications of the initial-state of the inflaton field can induce a departure from Gaussianity and leave a testable imprint on the higher order correlations of the CMB and large scale structures in the Universe. We focus on the bispectrum statistics of the primordial curvature perturbation and its projection on the CMB. For a canonical single-field action the three-point correlator enhancement is localized, maximizing in the collinear limit, corresponding to enfolded or squashed triangles in comoving momentum space. We show that the available local and equilateral template are very insensitive to this localized enhancement and do not generate noteworthy constraints on initial-state modifications. On the other hand, when considering the addition of a dimension 8 higher order derivative term, we find a dominant rapidly oscillating contribution, which had previously been overlooked and whose significantly enhanced amplitude is independent of the triangle under consideration. Nevertheless, the oscillatory nature of (the sign of) the correlation function implies the signal is nearly orthogonal to currently available observational templates, strongly reducing the sensitivity to the enhancement. Constraints on departures from the standard Bunch-Davies vacuum state can be derived, but also depend on the next-to-leading terms. We emphasize that the construction and application of especially adapted templates could lead to CMB bispectrum constraints on modified initial states already competing with those derived from the power spectrum.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.