We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We call these flow-based likelihoods (FBL). We analyze the accuracy and precision of the reconstructed likelihoods on mock Gaussian data, and show that simply gauging the quality of samples drawn from the trained model is not a sufficient indicator that the true likelihood has been learned. We nevertheless demonstrate that the likelihood can be reconstructed to a precision equal to that of sampling error due to a finite sample size. We then apply FBLs to mock weak lensing convergence power spectra, a cosmological observable that is significantly non-Gaussian (NG). We find that the FBL captures the NG signatures in the data extremely well, while other commonly-used data-driven likelihoods, such as Gaussian mixture models and independent component analysis, fail to do so. This suggests that works that have found small posterior shifts in NG data with data-driven likelihoods such as these could be underestimating the impact of non-Gaussianity in parameter constraints. By introducing a suite of tests that can capture different levels of NG in the data, we show that the success or failure of traditional data-driven likelihoods can be tied back to the structure of the NG in the data. Unlike other methods, the flexibility of the FBL makes it successful at tackling different types of NG simultaneously. Because of this, and consequently their likely applicability across datasets and domains, we encourage their use for inference when sufficient mock data are available for training.

A promising avenue to measure the total, and potentially individual, mass of neutrinos consists of leveraging cosmological datasets, such as the cosmic microwave background and surveys of the large-scale structure of the universe. In order to obtain unbiased estimates of the neutrino mass, however, many effects ought to be included. Here we forecast, via a Markov Chain Monte Carlo likelihood analysis, whether measurements by two galaxy surveys: DESI and {\it Euclid}, when added to the CMB-S4 experiment, are sensitive to two effects that can alter neutrino-mass measurements. The first is the slight difference in the suppression of matter fluctuations that each neutrino-mass hierarchy generates, at fixed total mass. The second is the growth-induced scale-dependent bias (GISDB) of haloes produced by massive neutrinos. We find that near-future surveys can distinguish hierarchies with the same total mass only at the $1\sigma$ level; thus, while these are poised to deliver a measurement of the sum of neutrino masses, they cannot significantly discern the mass of each individual neutrino in the foreseeable future. We further find that neglecting the GISDB induces up to a $1\sigma$ overestimation of the total neutrino mass, and we show how to absorb this effect via a redshift-dependent parametrization of the scale-independent bias.

Cosmological data provide a powerful tool in the search for physics beyond the Standard Model (SM). An interesting target are light relics, new degrees of freedom which decoupled from the SM while relativistic. Nearly massless relics contribute to the radiation energy budget, and are commonly searched through variations in the effective number $N_{\rm eff}$ of neutrino species. Additionally, relics with masses on the eV scale (meV-10 eV) become non-relativistic before today, and thus behave as matter instead of radiation. This leaves an imprint in the clustering of the large-scale structure of the universe, as light relics have important streaming motions, mirroring the case of massive neutrinos. Here we forecast how well current and upcoming cosmological surveys can probe light massive relics (LiMRs). We consider minimal extensions to the SM by both fermionic and bosonic relic degrees of freedom. By combining current and upcoming cosmic-microwave-background and large-scale-structure surveys, we forecast the significance at which each LiMR, with different masses and temperatures, can be detected. We find that a very large coverage of parameter space will be attainable by upcoming experiments, opening the possibility of exploring uncharted territory for new physics beyond the SM.

Galaxy-galaxy strong gravitational lenses have become a popular probe of dark matter (DM) by providing a window into structure formation on the smallest scales. In particular, the convergence power spectrum of subhalos within lensing galaxies has been suggested as a promising observable to study DM. However, the distances involved in strong-lensing systems are vast, and we expect the relevant volume to contain line-of-sight (LOS) halos that are not associated with the main lens. We develop a formalism to calculate the effect of LOS halos as an effective convergence power spectrum. The multi-lens plane equation couples the angular deflections of consecutive lens planes, but by assuming that the perturbations due to the LOS halos are small, we show that they can be projected onto the main-lens plane as effective subhalos. We test our formalism by simulating lensing systems using the full multi-plane lens equation and find excellent agreement. We show how the relative contribution of LOS halos and subhalos depends on the source and lens redshift, as well as the assumed halo and subhalo mass functions. For a fiducial system with fraction of DM halo mass in substructure $f_{\rm sub}=0.4\%$ for subhalo masses $[10^5-10^8]\rm{M}_{\odot}$, the interloper contribution to the power spectrum is at least several times greater than that of subhalos for source redshifts $z_s\gtrsim0.5$. Furthermore, it is likely that for the SLACS and BELLS lenses the interloper contribution dominates: $f_{\rm sub}\gtrsim2\%$ ($4\%$) is needed for subhalos to dominate in SLACS (BELLS), which is higher than current upper bounds on $f_{\rm sub}$ for our mass range. Since the halo mass function is better understood from first principles, the dominance of interlopers in galaxy-galaxy lenses with high-quality imaging can be seen as a significant advantage when translating this observable into a constraint on DM.

Angular cosmological correlators are infamously difficult to compute due to the highly oscillatory nature of the projection integrals. Motivated by recent development on analytic approaches to cosmological perturbation theory, in this paper we present an efficient method for computing cosmological four-point correlations in angular space, generalizing previous works on lower-point functions. This builds on the FFTLog algorithm that approximates the matter power spectrum as a sum over power-law functions, which makes certain momentum integrals analytically solvable. The computational complexity is drastically reduced for correlators in a "separable" form---we define a suitable notion of separability for cosmological trispectra, and derive formulas for angular correlators of different separability classes. As an application of our formalism, we compute the angular galaxy trispectrum at tree level, with and without primordial non-Gaussianity. This includes effects of redshift space distortion and bias parameters up to cubic order. We also compute the non-Gaussian covariance of the angular matter power spectrum due to the connected four-point function, beyond the Limber approximation. We demonstrate that, in contrast to the standard lore, the Limber approximation can fail for the non-Gaussian covariance computation even for large multipoles.

We present a novel technique for Cosmic Microwave Background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to Generalized Morphological Component Analysis (GMCA), we introduce Hierarchical GMCA (HGMCA), a Bayesian hierarchical graphical model for source separation. We test our method on $N_{\rm side}=256$ simulated sky maps that include dust, synchrotron, free-free and anomalous microwave emission, and show that HGMCA reduces foreground contamination by $25\%$ over GMCA in both the regions included and excluded by the Planck UT78 mask, decreases the error in the measurement of the CMB temperature power spectrum to the $0.02-0.03\%$ level at $\ell>200$ (and $<0.26\%$ for all $\ell$), and reduces correlation to all the foregrounds. We find equivalent or improved performance when compared to state-of-the-art Internal Linear Combination (ILC)-type algorithms on these simulations, suggesting that HGMCA may be a competitive alternative to foreground separation techniques previously applied to observed CMB data. Additionally, we show that our performance does not suffer when we perturb model parameters or alter the CMB realization, which suggests that our algorithm generalizes well beyond our simplified simulations. Our results open a new avenue for constructing CMB maps through Bayesian hierarchical analysis.

The distribution of matter fluctuations in our universe is key for understanding the nature of dark matter and the physics of the early cosmos. Different observables have been able to map this distribution at large scales, corresponding to wavenumbers $k\lesssim 10$ Mpc$^{-1}$, but smaller scales remain much less constrained. The 21-cm line is a promising tracer of early stellar formation, which took place in small haloes (with masses $M\sim 10^6-10^8M_\odot$), formed out of matter overdensities with wavenumbers as large as $k\approx100$ Mpc$^{-1}$. Here we forecast how well both the 21-cm global signal, and its fluctuations, could probe the matter power spectrum during cosmic dawn ($z=12-25$). We find that the long-wavelength modes (with $k\lesssim40$ Mpc$^{-1}$) are highly degenerate with astrophysical parameters, whereas the modes with $k= (40-80)$ Mpc$^{-1}$ are more readily observable. This is further illustrated in terms of the principal components of the matter power spectrum, which peak at $k\sim 50$ Mpc$^{-1}$ both for a typical experiment measuring the 21-cm global signal and its fluctuations. We find that, imposing broad priors on astrophysical parameters, a global-signal experiment can measure the amplitude of the matter power spectrum integrated over $k= (40-80)$ Mpc$^{-1}$ with a precision of tens of percent. A fluctuation experiment, on the other hand, can constrain the power spectrum to a similar accuracy over both the $k=(40-60)$ Mpc$^{-1}$ and $(60-80)$ Mpc$^{-1}$ ranges even without astrophysical priors. The constraints outlined in this work would be able to test the behavior of dark matter at the smallest scales yet measured, for instance probing warm-dark matter masses up to $m_{\rm WDM}=8$ keV for the global signal and $14$ keV for the 21-cm fluctuations. This could shed light on the nature of dark matter beyond the reach of other cosmic probes.

Strong gravitational lensing is a promising way of uncovering the nature of dark matter, by finding perturbations to images that cannot be well accounted for by modeling the lens galaxy without additional structure, be it subhalos (smaller halos within the smooth lens) or line-of-sight (LOS) halos. We present results attempting to infer the presence of substructure from images without requiring an intermediate step in which a smooth model has to be subtracted, using a simple convolutional neural network (CNN). We find that the network is only able to infer the presence of subhalos with $>75\%$ accuracy when they have masses of $\geq 5 \times 10^9$M$_{\odot}$ if they lie within the main lens galaxy. Since less massive foreground LOS halos can have the same effect as higher mass subhalos, the CNN can probe lower masses in the halo mass function. The accuracy does not improve significantly if we add a population of less massive subhalos. With the expectation of experiments such as HST and Euclid yielding thousands of high-quality strong lensing images in the next years, having a way of analyzing images quickly to identify candidates that merit further analysis to determine individual subhalo properties while preventing extensive resources being used for images that would yield null detections could be very useful. By understanding the sensitivity as a function of substructure mass, non-detections could be combined with the information from images with substructure to constrain the cold dark matter scenario, in particular if the sensitivity can be pushed to lower masses.

The forthcoming generation of wide-field galaxy surveys will probe larger volumes and galaxy densities, thus allowing for a much larger signal-to-noise ratio for higher-order clustering statistics, in particular the galaxy bispectrum. Extracting this information, however, is more challenging than using the power spectrum due to more complex theoretical modeling, as well as significant computational cost of evaluating the bispectrum signal and the error budget. To overcome these challenges, several proxy statistics have been proposed in the literature, which partially or fully capture the information in the bispectrum, while being computationally less expensive than the bispectrum. One such statistics are {\it weighted skew-spectra}, which are cross-spectra of the density field and appropriately weighted quadratic fields. Using Fisher forecasts, we show that the information in these skew-spectra is equivalent to that in the bispectrum for parameters that appear as amplitudes in the bispectrum model, such as galaxy bias parameters or the amplitude of primordial non-Gaussianity. We consider three shapes of the primordial bispectrum: local, equilateral and that due to massive particles with spin two during inflation. To obtain constraints that match those from a measurement of the full bispectrum, we find that it is crucial to account for the full covariance matrix of the skew-spectra.

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.

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

We present the strongest constraints to date on anisotropies of cosmic microwave background (CMB) polarization rotation derived from 150 GHz data taken by the BICEP2/Keck Array CMB experiments up to and including the 2014 observing season (BK14). The definition of the polarization angle in BK14 maps has gone through self-calibration in which the overall angle is adjusted to minimize the observed TB and EB power spectra. After this procedure, the QU maps lose sensitivity to a uniform polarization rotation but are still sensitive to anisotropies of polarization rotation. This analysis places constraints on the anisotropies of polarization rotation, which could be generated by CMB photons interacting with axionlike pseudoscalar fields or Faraday rotation induced by primordial magnetic fields. The sensitivity of BK14 maps ($\sim 3\mu$K-arcmin) makes it possible to reconstruct anisotropies of the polarization rotation angle and measure their angular power spectrum much more precisely than previous attempts. Our data are found to be consistent with no polarization rotation anisotropies, improving the upper bound on the amplitude of the rotation angle spectrum by roughly an order of magnitude compared to the previous best constraints. Our results lead to an order of magnitude better constraint on the coupling constant of the Chern-Simons electromagnetic term $g_{a\gamma}\leq 7.2\times 10^{-2}/H_I$ (95% confidence) than the constraint derived from the B-mode spectrum, where $H_I$ is the inflationary Hubble scale. This constraint leads to a limit on the decay constant of $10^{-6}\lesssim f_a/M_{\rm pl}$ at mass range of $10^{-33}< m_a< 10^{-28}$ eV for $r=0.01$, assuming $g_{a\gamma}\sim\alpha/(2\pi f_a)$ with $\alpha$ denoting the fine structure constant. The upper bound on the amplitude of the primordial magnetic fields is 30nG (95% confidence) from the polarization rotation anisotropies.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

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.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

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.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

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.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

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.

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

The next generation of instruments designed to measure the polarization of the cosmic microwave background (CMB) will provide a historic opportunity to open the gravitational wave window to the primordial Universe. Through high sensitivity searches for primordial gravitational waves, and tighter limits on the energy released in processes like phase transitions, the CMB polarization data of the next decade has the potential to transform our understanding of the laws of physics underlying the formation of the Universe.

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

We investigate the bounds on the sum of neutrino masses in a cosmic-acceleration scenario where the equation of state $w(z)$ of dark energy (DE) is constructed in a model-independent way, using a basis of principal components (PCs) that are allowed to cross the phantom barrier $w(z)=-1$. We find that the additional freedom provided to $w(z)$ means the DE can undo changes in the background expansion induced by massive neutrinos at low redshifts. This has two significant consequences: (1) it leads to a substantial increase in the upper bound for the sum of the neutrino masses ($M_{\nu} < 0.33 - 0.55$ eV at the 95\% C.L. depending on the data sets and number of PCs included) compared to studies that choose a specific parametrization for $w(z)$; and (2) it causes $\sim1\sigma$ deviations from $\Lambda$CDM in the luminosity distance and the Hubble expansion rate at higher redshifts ($z \gtrsim 2$), where the contribution of DE is subdominant and there is little constraining data. The second point consequently means that there are also observable deviations in the shear power spectrum and in the matter power spectrum at low redshift, since the clustering of matter throughout cosmic time depends on the expansion rate. This provides a compelling case to pursue high-$z$ BAO and SN measurements as a way of disentangling the effects of neutrinos and dark energy. Finally, we find that the additional freedom given to the dark energy component has the effect of lowering $S_8$ with respect to $\Lambda$CDM.

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

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

Dark matter (DM) could couple to particles in the Standard Model (SM) through a light vector mediator. In the limit of small coupling, this portal could be responsible for producing the observed DM abundance through a mechanism known as freeze-in. Furthermore, the requisite DM-SM couplings provide a concrete benchmark for direct and indirect searches for DM. In this paper, we present updated calculations of the relic abundance for DM produced by freeze-in through a light vector mediator. We identify an additional production channel: the decay of photons that acquire an in-medium plasma mass. These plasmon decays are a dominant channel for DM production for sub-MeV DM masses, and including this channel leads to a significant reduction in the predicted signal strength for DM searches. Accounting for production from both plasmon decays and annihilations of SM fermions, the DM acquires a highly non-thermal phase space distribution which impacts the cosmology at later times; these cosmological effects will be explored in a companion paper.

The BICEP and Keck Array experiments are a suite of small-aperture refracting telescopes observing the microwave sky from the South Pole. They target the degree-scale B-mode polarization signal imprinted in the Cosmic Microwave Background (CMB) by primordial gravitational waves. Such a measurement would shed light on the physics of the very early universe. While BICEP2 observed for the first time a B-mode signal at 150 GHz, higher frequencies from the Planck satellite showed that it could be entirely due to the polarized emission from Galactic dust, though uncertainty remained high. Keck Array has been observing the same region of the sky for several years, with an increased detector count, producing the deepest polarized CMB maps to date. New detectors at 95 GHz were installed in 2014, and at 220 GHz in 2015. These observations enable a better constraint of galactic foreground emissions, as presented here. In 2015, BICEP2 was replaced by BICEP3, a 10 times higher throughput telescope observing at 95 GHz, while Keck Array is now focusing on higher frequencies. In the near future, BICEP Array will replace Keck Array, and will allow unprecedented sensitivity to the gravitational wave signal. High resolution observations from the South Pole Telescope (SPT) will also be used to remove the lensing contribution to B-modes.

We present results from an analysis of all data taken by the BICEP2/Keck CMB polarization experiments up to and including the 2015 observing season. This includes the first Keck Array observations at 220 GHz and additional observations at 95 & 150 GHz. The $Q/U$ maps reach depths of 5.2, 2.9 and 26 $\mu$K$_{cmb}$ arcmin at 95, 150 and 220 GHz respectively over an effective area of $\approx 400$ square degrees. The 220 GHz maps achieve a signal-to-noise on polarized dust emission approximately equal to that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz. We evaluate the joint likelihood of the spectra versus a multicomponent model of lensed-$\Lambda$CDM+$r$+dust+synchrotron+noise. The foreground model has seven parameters, and we impose priors on some of these using external information from Planck and WMAP derived from larger regions of sky. The model is shown to be an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint $r_{0.05}<0.07$ at 95% confidence, which tightens to $r_{0.05}<0.06$ in conjunction with Planck temperature measurements and other data. The lensing signal is detected at $8.8 \sigma$ significance. Running maximum likelihood search on simulations we obtain unbiased results and find that $\sigma(r)=0.020$. These are the strongest constraints to date on primordial gravitational waves.

The epoch of the formation of the first stars, known as the cosmic dawn, has emerged as a new arena in the search for dark matter. In particular, the first claimed 21-cm detection exhibits a deeper global absorption feature than expected, which could be caused by a low baryonic temperature, and has been interpreted as a sign for electromagnetic interactions between baryons and dark matter. This hypothesis has a striking prediction: large temperature anisotropies sourced by the velocity-dependent cooling of the baryons. However, in order to remain consistent with the rest of cosmological observations, only part of the dark matter is allowed to be charged, and thus interactive. Here we compute, for the first time, the 21-cm fluctuations caused by a charged subcomponent of the dark matter, including both the pre- and post-recombination evolution of all fluids. We find that, for the same parameters that can explain the anomalous 21-cm absorption signal, any percent-level fraction of charged dark matter would source novel 21-cm fluctuations with a unique acoustic spectrum, and with an amplitude above any other known effects. These fluctuations are uncorrelated with the usual adiabatic anisotropies, and would be observable at high significance with interferometers such as LOFAR and HERA, thus providing a novel probe of dark matter at cosmic dawn.

Strong gravitational lensing has been identified as a promising astrophysical probe to study the particle nature of dark matter. In this paper we present a detailed study of the power spectrum of the projected mass density (convergence) field of substructure in a Milky Way-sized halo. This power spectrum has been suggested as a key observable that can be extracted from strongly lensed images and yield important clues about the matter distribution within the lens galaxy. We use two different $N$-body simulations from the ETHOS framework: one with cold dark matter and another with self-interacting dark matter and a cutoff in the initial power spectrum. Despite earlier works that identified $ k \gtrsim 100$ kpc$^{-1}$ as the most promising scales to learn about the particle nature of dark matter we find that even at lower wavenumbers - which are actually within reach of observations in the near future - we can gain important information about dark matter. Comparing the amplitude and slope of the power spectrum on scales $0.1 \lesssim k/$kpc$^{-1} \lesssim 10$ from lenses at different redshifts can help us distinguish between cold dark matter and other exotic dark matter scenarios that alter the abundance and central densities of subhalos. Furthermore, by considering the contribution of different mass bins to the power spectrum we find that subhalos in the mass range $10^7 - 10^8$ M$_{\odot}$ are on average the largest contributors to the power spectrum signal on scales $2 \lesssim k/$kpc$^{-1} \lesssim 15$, despite the numerous subhalos with masses $> 10^8$ M$_{\odot}$ in a typical lens galaxy. Finally, by comparing the power spectra obtained from the subhalo catalogs to those from the particle data in the simulation snapshots we find that the seemingly-too-simple halo model is in fact a fairly good approximation to the much more complex array of substructure in the lens.

BICEP Array is a degree-scale Cosmic Microwave Background (CMB) experiment that will search for primordial B-mode polarization while constraining Galactic foregrounds. BICEP Array will be comprised of four receivers to cover a broad frequency range with channels at 30/40, 95, 150 and 220/270 GHz. The first low-frequency receiver will map synchrotron emission at 30 and 40 GHz and will deploy to the South Pole at the end of 2019. In this paper, we give an overview of the BICEP Array science and instrument, with a focus on the detector module. We designed corrugations in the metal frame of the module to suppress unwanted interactions with the antenna-coupled detectors that would otherwise deform the beams of edge pixels. This design reduces the residual beam systematics and temperature-to-polarization leakage due to beam steering and shape mismatch between polarized beam pairs. We report on the simulated performance of single- and wide-band corrugations designed to minimize these effects. Our optimized design alleviates beam differential ellipticity caused by the metal frame to about 7% over 57% bandwidth (25 to 45 GHz), which is close to the level due the bare antenna itself without a metal frame. Initial laboratory measurements are also presented.

Targeting faint polarization patterns arising from Primordial Gravitational Waves in the Cosmic Microwave Background requires excellent observational sensitivity. Optical elements in small aperture experiments such as Bicep3 and Keck Array are designed to optimize throughput and minimize losses from transmission, reflection and scattering at millimeter wavelengths. As aperture size increases, cryostat vacuum windows must withstand larger forces from atmospheric pressure and the solution has often led to a thicker window at the expense of larger transmission loss. We have identified a new candidate material for the fabrication of vacuum windows: with a tensile strength two orders of magnitude larger than previously used materials, woven high-modulus polyethylene could allow for dramatically thinner windows, and therefore significantly reduced losses and higher sensitivity. In these proceedings we investigate the suitability of high-modulus polyethylene windows for ground-based CMB experiments, such as current and future receivers in the Bicep/Keck Array program. This includes characterizing their optical transmission as well as their mechanical behavior under atmospheric pressure. We find that such ultra-thin materials are promising candidates to improve the performance of large-aperture instruments at millimeter wavelengths, and outline a plan for further tests ahead of a possible upcoming field deployment of such a science-grade window.

Bicep Array is a cosmic microwave background (CMB) polarization experiment that will begin observing at the South Pole in early 2019. This experiment replaces the five Bicep2 style receivers that compose the Keck Array with four larger Bicep3 style receivers observing at six frequencies from 30 to 270GHz. The 95GHz and 150GHz receivers will continue to push the already deep Bicep/Keck CMB maps while the 30/40GHz and 220/270GHz receivers will constrain the synchrotron and galactic dust foregrounds respectively. Here we report on the design and performance of the Bicep Array instruments focusing on the mount and cryostat systems.

BICEP Array is the newest multi-frequency instrument in the BICEP/Keck Array program. It is comprised of four 550 mm aperture refractive telescopes observing the polarization of the cosmic microwave background (CMB) at 30/40, 95, 150 and 220/270 GHz with over 30,000 detectors. We present an overview of the receiver, detailing the optics, thermal, mechanical, and magnetic shielding design. BICEP Array follows BICEP3's modular focal plane concept, and upgrades to 6" wafer to reduce fabrication with higher detector count per module. The first receiver at 30/40 GHz is expected to start observing at the South Pole during the 2019-20 season. By the end of the planned BICEP Array program, we project $\sigma(r) \sim 0.003$, assuming current modeling of polarized Galactic foreground and depending on the level of delensing that can be achieved with higher resolution maps from the South Pole Telescope.

BICEP3 is a 520mm aperture on-axis refracting telescope observing the polarization of the cosmic microwave background (CMB) at 95GHz in search of the B-mode signal originating from inflationary gravitational waves. BICEP3's focal plane is populated with modularized tiles of antenna-coupled transition edge sensor (TES) bolometers. BICEP3 was deployed to the South Pole during 2014-15 austral summer and has been operational since. During the 2016-17 austral summer, we implemented changes to optical elements that lead to better noise performance. We discuss this upgrade and show the performance of BICEP3 at its full mapping speed from the 2017 and 2018 observing seasons. BICEP3 achieves an order-of-magnitude improvement in mapping speed compared to a Keck 95GHz receiver. We demonstrate $6.6\mu K\sqrt{s}$ noise performance of the BICEP3 receiver.

Features during inflation and reionization leave corresponding features in the temperature and polarization power spectra that could potentially explain anomalies in the Planck 2015 data but require a joint analysis to disentangle. We study the interplay between these two effects using a model-independent parametrization of the inflationary power spectrum and the ionization history. Preference for a sharp suppression of large scale power is driven by a feature in the temperature power spectrum at multipoles $\ell \sim 20$, whereas preference for a component of high redshift ionization is driven by a sharp excess of polarization power at $\ell \sim 10$ when compared with the lowest multipoles. Marginalizing inflationary freedom does not weaken the preference for $z \gtrsim 10$ ionization, whereas marginalizing reionization freedom slightly enhances the preference for an inflationary feature but can also mask its direct signature in polarization. The inflation and reionization interpretation of these features makes predictions for the polarization spectrum which can be tested in future precision measurements especially at $10\lesssim \ell \lesssim 40$.

Cosmological data can be used to search for---and characterize---light particles in the standard model, if these populate our Universe. In addition to the well-known effect of these light relics in the background cosmology, usually parametrized through a change in the effective number $N_{\rm eff}$ of neutrino species, these particles can become nonrelativistic at later times, affecting the growth of matter fluctuations due to their thermal velocities. An extensively studied example is that of massive neutrinos, which are known to produce a suppression in the matter power spectrum due to their free streaming. Galaxies, as biased traces of matter fluctuations, can therefore provide us with a wealth of information about both known and unknown degrees of freedom in the standard model. To harness this information, however, the galaxy bias has to be determined in the presence of massive relics, which is expected to vary with scale. Here we present the code RelicFast, which efficiently computes the scale-dependent bias induced by relics of different masses, spins, and temperatures, through spherical collapse and the peak-background split. Using this code, we find that, in general, the bias induced by light relics partially compensates the suppression of power, and should be accounted for in any search for relics with galaxy data. In particular, for the case of neutrinos, we find that both the normal and inverted hierarchies present a percent-level step in the Lagrangian bias, with a size scaling linearly with the neutrino-mass sum, in agreement with recent simulations. RelicFast can compute halo bias in under a second, allowing for this effect to be properly included for different cosmologies, and light relics, at little computational cost.

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

Massive spinning particles, if present during inflation, lead to a distinctive bispectrum of primordial perturbations, the shape and amplitude of which depend on the masses and spins of the extra particles. This signal, in turn, leaves an imprint in the statistical distribution of galaxies; in particular, as a non-vanishing galaxy bispectrum, which can be used to probe the masses and spins of these particles. In this paper, we present for the first time a new theoretical template for the bispectrum generated by massive spinning particles, valid for a general triangle configuration. We then proceed to perform a Fisher-matrix forecast to assess the potential of two next-generation spectroscopic galaxy surveys, EUCLID and DESI, to constrain the primordial non-Gaussianity sourced by these extra particles. We model the galaxy bispectrum using tree-level perturbation theory, accounting for redshift-space distortions and the Alcock-Paczynski effect, and forecast constraints on the primordial non-Gaussianity parameters marginalizing over all relevant biases and cosmological parameters. Our results suggest that these surveys would potentially be sensitive to any primordial non-Gaussianity with an amplitude larger than $f_{\rm NL}\approx 1$, for massive particles with spins 2, 3, and 4. Interestingly, if non-Gaussianities are present at that level, these surveys will be able to infer the masses of these spinning particles to within tens of percent. If detected, this would provide a very clear window into the particle content of our Universe during inflation.

We derive new limits on the elastic scattering cross-section between baryons and dark matter using Cosmic Microwave Background data from the Planck satellite and measurements of the Lyman-alpha forest flux power spectrum from the Sloan Digital Sky Survey. Our analysis addresses generic cross sections of the form $\sigma\propto v^n$, where v is the dark matter-baryon relative velocity, allowing for constraints on the cross section independent of specific particle physics models. We include high-$\ell$ polarization data from Planck in our analysis, improving over previous constraints. We apply a more careful treatment of dark matter thermal evolution than previously done, allowing us to extend our constraints down to dark matter masses of $\sim$MeV. We show in this work that cosmological probes are complementary to current direct detection and astrophysical searches.

Strong lensing is a sensitive probe of the small-scale density fluctuations in the Universe. We implement a novel approach to modeling strongly lensed systems using probabilistic cataloging, which is a transdimensional, hierarchical, and Bayesian framework to sample from a metamodel (union of models with different dimensionality) consistent with observed photon count maps. Probabilistic cataloging allows us to robustly characterize modeling covariances within and across lens models with different numbers of subhalos. Unlike traditional cataloging of subhalos, it does not require model subhalos to improve the goodness of fit above the detection threshold. Instead, it allows the exploitation of all information contained in the photon count maps, for instance, when constraining the subhalo mass function. We further show that, by not including these small subhalos in the lens model, fixed-dimensional inference methods can significantly mismodel the data. Using a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass function can be probed even when many subhalos in the sample catalogs are individually below the detection threshold and would be absent in a traditional catalog. With the planned Wide Field Infrared Space Telescope (WFIRST), simultaneous probabilistic cataloging of dark subhalos in high-resolution, deep strong lens images has the potential to constrain the subhalo mass function at even lower masses.

The presence of additional particles during inflation leads to non-Gaussianity in late-time correlators of primordial curvature perturbations. The shape and amplitude of this signal depend on the mass and spin of the extra particles. Constraints on this distinct form of primordial non-Gaussianity, therefore, provide a wealth of information on the particle content during inflation. We investigate the potential of upcoming galaxy surveys in constraining such a signature through its impact on the observed galaxy power spectrum. Primordial non-Gaussianity of various shapes induces a scale-dependent bias on tracers of large-scale structure, such as galaxies. Using this signature we obtain constraints on massive particles during inflation, which can have non-zero spins. In particular, we show that the prospects for constraining particles with spins 0 and 1 are promising, while constraining particles with spin 2 from power spectrum alone seems challenging. We show that the multi-tracer technique can significantly improve the constraints from the power spectrum by at least an order of magnitude. Furthermore, we analyze the effect of non-linearities due to gravitational evolution on the forecasted constraints on the masses of the extra particles and the amplitudes of the imprinted non-Gaussian signal. We find that gravitational evolution affects the constraints by less than a factor of 2.

Studying the smallest self-bound dark matter structure in our Universe can yield important clues about the fundamental particle nature of dark matter. Galaxy-scale strong gravitational lensing provides a unique way to detect and characterize dark matter substructures at cosmological distances from the Milky Way. Within the cold dark matter (CDM) paradigm, the number of low-mass subhalos within lens galaxies is expected to be large, implying that their contribution to the lensing convergence field is approximately Gaussian and could thus be described by their power spectrum. We develop here a general formalism to compute from first principles the substructure convergence power spectrum for different populations of dark matter subhalos. As an example, we apply our framework to two distinct subhalo populations: a truncated Navarro-Frenk-White subhalo population motivated by standard CDM, and a truncated cored subhalo population motivated by self-interacting dark matter (SIDM). We study in detail how the subhalo abundance, mass function, internal density profile, and concentration affect the amplitude and shape of the substructure power spectrum. We determine that the power spectrum is mostly sensitive to a specific combination of the subhalo abundance and moments of the mass function, as well as to the average tidal truncation scale of the largest subhalos included in the analysis. Interestingly, we show that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the subhalos. In particular, the SIDM power spectrum exhibits a characteristic steepening at large wave number absent in the CDM power spectrum, opening the possibility of using this observable, if at all measurable, to discern between these two scenarios.

Cosmology is intrinsically intertwined with questions in fundamental physics. The existence of non-baryonic dark matter requires new physics beyond the Standard Model of elemenatary-particle interactions and Einstein's general relativity, as does the accelerating expansion of the universe. Current tensions between various cosmological measurements may be harbingers of yet more new physics. Progress on understanding dark matter and cosmic acceleration requires long term, high-precision measurements and excellent control of systematics, demanding observational programs that are often outside the discovery/characterization mode that drives many areas of astronomy. We outline potential programs through which the Hubble Space Telescope (HST) could have a major impact on issues in fundamental physics in the coming years. To realize this impact, we suggest the introduction of a "HST Fundamental Physics" observational program that would be subject to a modified proposal and review process.

Through likelihood analyses of both current and future data that constrain both the expansion history of the universe and the clustering of matter fluctuations, we provide falsifiable predictions for three broad classes of models that explain the accelerated expansions of the universe: $\Lambda$CDM, the quintessence scenario and a more general class of smooth dark energy models that can cross the phantom barrier $w(z)=-1$. Our predictions are model independent in the sense that we do not rely on a specific parametrization, but we instead use a principal component (PC) basis function constructed a priori from a noise model of supernovae and Cosmic Microwave Background observations. For the supernovae measurements, we consider two type of surveys: the current JLA and the upcoming WFIRST surveys. We show that WFIRST will be able to improve growth predictions in curved models significantly. The remaining degeneracy between spatial curvature and $w(z)$ could be overcome with improved measurements of $\sigma_8 \Omega_m^{1/2}$, a combination that controls the amplitude of the growth of structure. We also point out that a PC-based Figure of Merit reveals that the usual two-parameter description of $w(z)$ does not exhaust the information that can be extracted from current data (JLA) or future data (WFIRST).

The amplitude of large-scale matter fluctuations inferred from the observed Sunyaev-Zeldovich (SZ) cluster mass function and from weak gravitational lensing studies, when taken at face value, is in tension with measurements of the cosmic microwave background (CMB) and baryon acoustic oscillation (BAO). In this work, we revisit whether this possible discrepancy can be attributed to new interactions in the dark matter sector. Focusing on a cosmological model where dark matter interacts with a dark radiation species until the epoch of matter-radiation equality, we find that measurements of the Lyman-alpha flux power spectrum from the Sloan Digital Sky Survey provides no support to the hypothesis that new dark matter interactions can resolve the possible tension between CMB and large-scale structure (LSS). Indeed, while the addition of dark matter-dark radiation interactions leads to an improvement of $2\Delta\ln\mathcal{L}=12$ with respect to the standard $\Lambda$ cold dark matter ($\Lambda$CDM) model when only CMB, BAO, and LSS data are considered, the inclusion of Lyman-alpha data reduces the improvement of the fit to $2\Delta\ln\mathcal{L}=6$ relative to $\Lambda$CDM. We thus conclude that the statistical evidence for new dark matter interactions (largely driven by the Planck SZ dataset) is marginal at best, and likely caused by systematics in the data. We also perform a Fisher forecast analysis for the reach of a future dataset composed of a CMB-S4 experiment combined with the Large Synoptic Survey Telescope galaxy survey. We find that the constraint on the effective number of fluid-like dark radiation species, $\Delta N_{\rm fluid}$, will be improved by an order of magnitude compared to current bounds.

We explore the relationship between features in the Planck 2015 temperature and polarization data, shifts in the cosmological parameters, and features from inflation. Residuals in the temperature data at low multipole $\ell$, which are responsible for the high $H_0\approx 70$ km s$^{-1}$Mpc$^{-1}$ and low $\sigma_8\Omega_m^{1/2}$ values from $\ell<1000$ in power-law $\Lambda$CDM models, are better fit to inflationary features with a $1.9\sigma$ preference for running of the running of the tilt or a stronger $99\%$ CL local significance preference for a sharp drop in power around $k=0.004$ Mpc$^{-1}$ in generalized slow roll and a lower $H_0\approx 67$ km s$^{-1}$Mpc$^{-1}$. The same in-phase acoustic residuals at $\ell>1000$ that drive the global $H_0$ constraints and appear as a lensing anomaly also favor running parameters which allow even lower $H_0$, but not once lensing reconstruction is considered. Polarization spectra are intrinsically highly sensitive to these parameter shifts, and even more so in the Planck 2015 TE data due to an outlier at $\ell \approx 165$, which disfavors the best fit $H_0$ $\Lambda$CDM solution by more than $2\sigma$, and high $H_0$ value at almost $3\sigma$. Current polarization data also slightly enhance the significance of a sharp suppression of large-scale power but leave room for large improvements in the future with cosmic variance limited $E$-mode measurements.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

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

Correlations between intrinsic galaxy shapes on large-scales arise due to the effect of the tidal field of the large-scale structure. Anisotropic primordial non-Gaussianity induces a distinct scale-dependent imprint in these tidal alignments on large scales. Motivated by the observational finding that the alignment strength of luminous red galaxies depends on how galaxy shapes are measured, we study the use of two different shape estimators as a multi-tracer probe of intrinsic alignments. We show, by means of a Fisher analysis, that this technique promises a significant improvement on anisotropic non-Gaussianity constraints over a single-tracer method. For future weak lensing surveys, the uncertainty in the anisotropic non-Gaussianity parameter, $A_2$, is forecast to be $\sigma(A_2)\approx 50$, $\sim 40\%$ smaller than currently available constraints from the bispectrum of the Cosmic Microwave Background. This corresponds to an improvement of a factor of $4-5$ over the uncertainty from a single-tracer analysis.

Primordial features are one of the most important extensions of the Standard Model of cosmology, providing a wealth of information on the primordial universe, ranging from discrimination between inflation and alternative scenarios, new particle detection, to fine structures in the inflationary potential. We study the prospects of future large-scale structure (LSS) surveys on the detection and constraints of these features. We classify primordial feature models into several classes, and for each class we present a simple template of power spectrum that encodes the essential physics. We study how well the most ambitious LSS surveys proposed to date, including both spectroscopic and photometric surveys, will be able to improve the constraints with respect to the current Planck data. We find that these LSS surveys will significantly improve the experimental sensitivity on features signals that are oscillatory in scales, due to the 3D information. For a broad range of models, these surveys will be able to reduce the errors of the amplitudes of the features by a factor of 5 or more, including several interesting candidates identified in the recent Planck data. Therefore, LSS surveys offer an impressive opportunity for primordial feature discovery in the next decade or two. We also compare the advantages of both types of surveys.

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.

A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.

We present measurements of polarization lensing using the 150 GHz maps which include all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season (BK14). Despite their modest angular resolution ($\sim 0.5^\circ$), the excellent sensitivity ($\sim 3\mu$K-arcmin) of these maps makes it possible to directly reconstruct the lensing potential using only information at larger angular scales ($\ell\leq 700$). From the auto-spectrum of the reconstructed potential we measure an amplitude of the spectrum to be $A^{\phi\phi}_{\rm L}=1.15\pm 0.36$ (Planck $\Lambda$CDM prediction corresponds to $A^{\phi\phi}_{\rm L}=1$), and reject the no-lensing hypothesis at 5.8$\sigma$, which is the highest significance achieved to date using an EB lensing estimator. Taking the cross-spectrum of the reconstructed potential with the Planck 2015 lensing map yields $A^{\phi\phi}_{\rm L}=1.13\pm 0.20$. These direct measurements of $A^{\phi\phi}_{\rm L}$ are consistent with the $\Lambda$CDM cosmology, and with that derived from the previously reported BK14 B-mode auto-spectrum ($A^{\rm BB}_{\rm L}=1.20\pm 0.17$). We perform a series of null tests and consistency checks to show that these results are robust against systematics and are insensitive to analysis choices. These results unambiguously demonstrate that the B-modes previously reported by BICEP / Keck at intermediate angular scales ($150\lesssim\ell\lesssim 350$) are dominated by gravitational lensing. The good agreement between the lensing amplitudes obtained from the lensing reconstruction and B-mode spectrum starts to place constraints on any alternative cosmological sources of B-modes at these angular scales.

We show that intrinsic (not lensing-induced) correlations between galaxy shapes offer a new probe of primordial non-Gaussianity and inflationary physics which is complementary to galaxy number counts. Specifically, intrinsic alignment correlations are sensitive to an anisotropic squeezed limit bispectrum of the primordial perturbations. Such a feature arises in solid inflation, as well as more broadly in the presence of light higher spin fields during inflation (as pointed out recently by Arkani-Hamed and Maldacena). We present a derivation of the all-sky two-point correlations of intrinsic shapes and number counts in the presence of non-Gaussianity with general angular dependence, and show that a quadrupolar (spin-2) anisotropy leads to the analog in galaxy shapes of the well-known scale-dependent bias induced in number counts by isotropic (spin-0) non-Gaussianity. Moreover, in presence of non-zero anisotropic non-Gaussianity, the quadrupole of galaxy shapes becomes sensitive to far superhorizon modes. These effects come about because long-wavelength modes induce a local anisotropy in the initial power spectrum, with which galaxies will correlate. We forecast that future imaging surveys could provide constraints on the amplitude of anisotropic non-Gaussianity that are comparable to those from the Cosmic Microwave Background (CMB). These are complementary as they probe different physical scales. The constraints, however, depend on the sensitivity of galaxy shapes to the initial conditions which we only roughly estimate from observed tidal alignments.

We present results from an analysis of all data taken by the BICEP2 & Keck Array CMB polarization experiments up to and including the 2014 observing season. This includes the first Keck Array observations at 95 GHz. The maps reach a depth of 50 nK deg in Stokes $Q$ and $U$ in the 150 GHz band and 127 nK deg in the 95 GHz band. We take auto- and cross-spectra between these maps and publicly available maps from WMAP and Planck at frequencies from 23 GHz to 353 GHz. An excess over lensed-LCDM is detected at modest significance in the 95x150 $BB$ spectrum, and is consistent with the dust contribution expected from our previous work. No significant evidence for synchrotron emission is found in spectra such as 23x95, or for correlation between the dust and synchrotron sky patterns in spectra such as 23x353. We take the likelihood of all the spectra for a multi-component model including lensed-LCDM, dust, synchrotron and a possible contribution from inflationary gravitational waves (as parametrized by the tensor-to-scalar ratio $r$), using priors on the frequency spectral behaviors of dust and synchrotron emission from previous analyses of WMAP and Planck data in other regions of the sky. This analysis yields an upper limit $r_{0.05}<0.09$ at 95% confidence, which is robust to variations explored in analysis and priors. Combining these $B$-mode results with the (more model-dependent) constraints from Planck analysis of CMB temperature plus BAO and other data, yields a combined limit $r_{0.05}<0.07$ at 95% confidence. These are the strongest constraints to date on inflationary gravitational waves.

The Keck Array is a system of cosmic microwave background (CMB) polarimeters, each similar to the BICEP2 experiment. In this paper we report results from the 2012 and 2013 observing seasons, during which the Keck Array consisted of five receivers all operating in the same (150 GHz) frequency band and observing field as BICEP2. We again find an excess of B-mode power over the lensed-$\Lambda$CDM expectation of $> 5 \sigma$ in the range $30 < \ell < 150$ and confirm that this is not due to systematics using jackknife tests and simulations based on detailed calibration measurements. In map difference and spectral difference tests these new data are shown to be consistent with BICEP2. Finally, we combine the maps from the two experiments to produce final Q and U maps which have a depth of 57 nK deg (3.4 $\mu$K arcmin) over an effective area of 400 deg$^2$ for an equivalent survey weight of 250,000 $\mu$K$^{-2}$. The final BB band powers have noise uncertainty a factor of 2.3 times better than the previous results, and a significance of detection of excess power of $> 6\sigma$.

If observations confirm BICEP2's claim of a tensor-scalar ratio $r\approx 0.2$ on CMB scales, then the inflationary consistency relation $n_{t}=-r/8$ predicts a small negative value for the tensor spectral index $n_t$. We show that future CMB polarization experiments should be able to confirm this prediction at several sigma. We also show how to properly extend the consistency relation to solar system scales, where the primordial gravitational wave density $\Omega_{gw}$ could be measured by proposed experiments such as the Big Bang Observer. This would provide a far more stringent test of the consistency relation and access much more detailed information about the early universe.

We report the results of a joint analysis of data from BICEP2/Keck Array and Planck. BICEP2 and Keck Array have observed the same approximately 400 deg$^2$ patch of sky centered on RA 0h, Dec. $-57.5\deg$. The combined maps reach a depth of 57 nK deg in Stokes $Q$ and $U$ in a band centered at 150 GHz. Planck has observed the full sky in polarization at seven frequencies from 30 to 353 GHz, but much less deeply in any given region (1.2 $\mu$K deg in $Q$ and $U$ at 143 GHz). We detect 150$\times$353 cross-correlation in $B$-modes at high significance. We fit the single- and cross-frequency power spectra at frequencies $\geq 150$ GHz to a lensed-$\Lambda$CDM model that includes dust and a possible contribution from inflationary gravitational waves (as parameterized by the tensor-to-scalar ratio $r$), using a prior on the frequency spectral behavior of polarized dust emission from previous \planck\ analysis of other regions of the sky. We find strong evidence for dust and no statistically significant evidence for tensor modes. We probe various model variations and extensions, including adding a synchrotron component in combination with lower frequency data, and find that these make little difference to the $r$ constraint. Finally we present an alternative analysis which is similar to a map-based cleaning of the dust contribution, and show that this gives similar constraints. The final result is expressed as a likelihood curve for $r$, and yields an upper limit $r_{0.05}<0.12$ at 95% confidence. Marginalizing over dust and $r$, lensing $B$-modes are detected at $7.0\,\sigma$ significance.

We conduct a model-independent analysis of temporal features during inflation in the large-scale CMB temperature power spectrum allowing for the possibility of non-negligible tensor contributions. Of 20 principal components of the inflationary history, the suppression of power at low multipoles beginning with a glitch at multipoles $\ell \sim 20-40$ implies deviations in 2-3 of them with 2-3$\sigma$ deviations in each, with larger values reflecting cases where tensors are allowed.If tensors are absent, the corresponding $E$-mode polarization features follow a similar pattern but are predicted to be up to twice as large. They offer the opportunity to soon double the significance of inflationary features or eliminate them as an explanation of temperature features. The tensor degeneracy with features in the temperature power spectrum is broken not only by $B$ but also by $E$-polarization. A precision measurement of $E$-mode polarization at multipoles from $\ell\sim 20-60$ can potentially provide an independent constraint on tensors that is less subject to dust foreground uncertainties.

In light of the recent BICEP2 B-mode polarization detection, which implies a large inflationary tensor-to-scalar ratio r_{0.05}=0.2^{+0.07}_{-0.05}, we re-examine the evidence for an extra sterile massive neutrino, originally invoked to account for the tension between the cosmic microwave background (CMB) temperature power spectrum and local measurements of the expansion rate H0 and cosmological structure. With only the standard active neutrinos and power-law scalar spectra, this detection is in tension with the upper limit of r<0.11 (95% confidence) from the lack of a corresponding low multipole excess in the temperature anisotropy from gravitational waves. An extra sterile species with the same energy density as is needed to reconcile the CMB data with H0 measurements can also alleviate this new tension. By combining data from the Planck and ACT/SPT temperature spectra, WMAP9 polarization, H_0, baryon acoustic oscillation and local cluster abundance measurements with BICEP2 data, we find the joint evidence for a sterile massive neutrino increases to DeltaNeff=0.98\pm 0.26 for the effective number and ms= 0.52\pm 0.13 eV for the effective mass or 3.8 sigma and 4 sigma evidence respectively. We caution the reader that these results correspond to a joint statistical evidence and, in addition, astrophysical systematic errors in the clusters and H0 measurements, and small-scale CMB data could weaken our conclusions.

The detection of B-modes in the Cosmic Microwave Background (CMB) polarization by the BICEP2 experiment, if interpreted as evidence for a primordial gravitational wave background, has enormous ramifications for cosmology and physics. It is crucial to test this hypothesis with independent measurements. A gravitational wave background leads to B-modes in galaxy shape correlations (shear) both through lensing and tidal alignment effects. Since the systematics and foregrounds of galaxy shapes and CMB polarization are entirely different, a detection of a cross-correlation between the two observables would provide conclusive proof for the existence of a primordial gravitational wave background. We find that upcoming weak lensing surveys will be able to detect the cross-correlation between B-modes of the CMB and galaxy shapes. However, this detection is not sufficient to confirm or falsify the hypothesis of a primordial origin for CMB B-mode polarization.

The recent BICEP2 measurement of primordial gravity waves (r = 0.2^{+0.07}_{-0.05}) appears to be in tension with the upper limit from WMAP (r<0.13 at 95% CL) and Planck (r<0.11 at 95% CL). We carefully quantify the level of tension and show that it is very significant (around 0.1% unlikely) when the observed deficit of large-scale temperature power is taken into account. We show that measurements of TE and EE power spectra in the near future will discriminate between the hypotheses that this tension is either a statistical fluke, or a sign of new physics. We also discuss extensions of the standard cosmological model that relieve the tension, and some novel ways to constrain them.

In this follow-up work to the High Energy Physics Community Summer Study 2013 (HEP CSS 2013, a.k.a. Snowmass), we explore the scientific capabilities of a future Stage-IV Cosmic Microwave Background polarization experiment (CMB-S4) under various assumptions on detector count, resolution, and sky coverage. We use the Fisher matrix technique to calculate the expected uncertainties in cosmological parameters in $\nu \Lambda$CDM that are especially relevant to the physics of fundamental interactions, including neutrino masses, effective number of relativistic species, dark-energy equation of state, dark-matter annihilation, and inflationary parameters. To further chart the landscape of future cosmology probes, we include forecasted results from the Baryon Acoustic Oscillation (BAO) signal as measured by DESI to constrain parameters that would benefit from low redshift information. We find the following best 1-sigma constraints: $\sigma$ constraints: $\sigma(M_{\nu})= 15$ meV, $\sigma(N_{\rm eff})= 0.0156$, Dark energy Figure of Merit = 303, $\sigma(p_{ann})= 0.00588\times3\times10^{-26}$ cm$^3$/s/GeV, $\sigma(\Omega_K)= 0.00074$, $\sigma(n_s)= 0.00110$, $\sigma(\alpha_s)= 0.00145$, and $\sigma(r)= 0.00009$. We also detail the dependences of the parameter constraints on detector count, resolution, and sky coverage.

The intrinsic alignments of galaxies are usually regarded as a contaminant to weak gravitational lensing observables. The alignment of Luminous Red Galaxies, detected unambiguously in observations from the Sloan Digital Sky Survey, can be reproduced by the linear tidal alignment model of Catelan, Kamionkowski & Blandford (2001) on large scales. In this work, we explore the cosmological information encoded in the intrinsic alignments of red galaxies. We make forecasts for the ability of current and future spectroscopic surveys to constrain local primordial non-Gaussianity and Baryon Acoustic Oscillations (BAO) in the cross-correlation function of intrinsic alignments and the galaxy density field. For the Baryon Oscillation Spectroscopic Survey, we find that the BAO signal in the intrinsic alignments is marginally significant with a signal-to-noise ratio of 1.8 and 2.2 with the current LOWZ and CMASS samples of galaxies, respectively, and increasing to 2.3 and 2.7 once the survey is completed. For the Dark Energy Spectroscopic Instrument and for a spectroscopic survey following the EUCLID redshift selection function, we find signal-to-noise ratios of 12 and 15, respectively. Local type primordial non-Gaussianity, parametrized by fNL = 10, is only marginally significant in the intrinsic alignments signal with signal-to-noise ratios < 2 for the three surveys considered.

We derive constraints on elastic scattering between baryons and dark matter using the cosmic microwave background (CMB) data from the Planck satellite and the Lyman-alpha forest data from the Sloan Digital Sky Survey. Elastic scattering allows baryons and dark matter to exchange momentum, affecting the dynamics of linear density perturbations in the early Universe. We derive constraints to scattering cross sections of the form sigma \propto v^n, allowing for a wide range of velocity dependencies with n between -4 and 2. We improve and correct previous estimates where they exist, including velocity-independent cross section as well as dark matter millicharge and electromagnetic dipole moments. Lyman-alpha forest data dominates the constraints for n>-3, where the improvement over CMB data alone can be several orders of magnitude. Dark matter-baryon scattering cannot affect the halo mass function on mass scales M>10^{12} M_{solar}. Our results imply, model-independently, that a baryon in the halo of a galaxy like our own Milky Way, does not scatter from dark matter particles during the age of the galaxy.

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.

We show that dark matter annihilation around the time of recombination can lead to growing ionization fraction perturbations, that track the linear collapse of matter over-densities. This amplifies small scale cosmological perturbations to the free electron density by a significant amount compared to the usual acoustic oscillations. Electron density perturbations distort the CMB, inducing secondary non-gaussianity. We calculate the CMB bispectrum from recombination, that is marginally observable by Planck. Even though electron perturbations can be markedly boosted compared with the Standard Model prediction, the dark matter effect in the CMB bispectrum turns out to be small and will be difficult to disentangle from the standard model in the foreseeable future.

We study the effect of the supersonic baryon--CDM flow, which has recently been shown to have a large effect on structure formation during the dark ages 10 <~ z <~ 1000, on the abundance of luminous, low-mass satellite galaxies around galaxies like the Milky Way. As the supersonic baryon--CDM flow significantly suppresses both the number of halos formed and the amount of baryons accreted onto such halos of masses 10^6 < M_{halo} / M_solar < 10^8 at z >~ 10, a large effect results on the stellar luminosity function before reionization. As halos of these masses are believed to have very little star formation after reionization due to the effects of photo-heating by the ultraviolet background, this effect persists to the present day. We calculate that the number of low-mass 10^6 < M_{halo} / M_solar < 10^8 halos that host luminous satellite galaxies today is typically suppressed by 50 percent, with values ranging up to 90 percent in regions where the initial supersonic velocity is high. We show that this previously-ignored cosmological effect resolves most of the tension between the observed and predicted number of low-mass satellites in the Milky Way, obviating the need for any other mass-dependent star-formation suppression before reionization.

We provide analytic solutions for the power spectrum and bispectrum of curvature fluctuations produced by a step feature in the inflaton potential, valid in the limit that the step is short and sharp. In this limit, the bispectrum is strongly scale dependent and its effective non-linearity attains a large oscillatory amplitude. The perturbations to the curvature power spectrum, on the other hand, remain a small component on top of the usual spectrum of fluctuations generated by slow roll. We utilize our analytic solutions to assess the observability of the predicted non-Gaussian signatures and show that, if present, only very sharp steps on scales larger than ~ 2 Gpc are likely to be able to be detected by Planck. Such features are not only consistent with WMAP7 data, but can also improve its likelihood by 2 Delta ln L ~ 12 for two extra parameters, the step location and height. If this improvement were due to a slow roll violating step as considered here, a bispectrum or corresponding polarization power spectrum detection would provide definitive checks as to its primordial origin.

We develop a fast technique based on the generalized slow roll (GSR) approach for computing the curvature bispectrum of inflationary models with features. We show that all triangle configurations can be expressed in terms of three simple integrals over the inflationary background with typical accuracy of better than ~20%. With a first order GSR approach the typical accuracy can be improved to better than the 5% level. We illustrate this technique with the step potential model that has been invoked to explain the WMAP temperature power spectrum glitches at ell ~ 20-40 and show that the maximum likelihood model falls short of observability by more than a factor of 100 in amplitude. We also explicitly demonstrate that the bispectrum consistency relation with the local slope of the power spectrum is satisfied for these models. In the GSR approach, the bispectrum arises from integrals of nearly the same function of the background slow-roll parameters as the power spectrum but with a stronger weight to the epoch before horizon crossing. Hence this technique enables reverse engineering of models with large bispectrum but small power spectrum features.

The predictions of the inflationary LCDM paradigm match today's high-precision measurements of the cosmic microwave background anisotropy extremely well. The same data put tight limits on other sources of anisotropy. Cosmic strings are a particularly interesting alternate source to constrain. Strings are topological defects, remnants of inflationary-era physics that persist after the big bang. They are formed in a variety of models of inflation, including string theory models such as brane inflation. We assume a "Nambu-Goto" model for strings, approximated by a collection of unconnected segments with zero width, and show that measurements of temperature anisotropy by the South Pole Telescope break a parameter degeneracy in the WMAP data, permitting us to place a strong upper limit on the possible string contribution to the CMB anisotropy: the power sourced by zero-width strings must be <1.75% (95% CL) of the total or the string tension Gmu <1.7x10^{-7}. These limits imply that the best hope for detecting strings in the CMB will come from B-mode polarization measurements at arcminute scales rather than the degree scale measurements pursued for gravitational wave detection.

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

Using a principal component (PC) basis that accommodates order unity features in the slow roll parameters as fine as 1/10 of a decade across more than 2 decades of the inflationary expansion, we test slow roll and single field inflation with the WMAP7 data. Detection of any non-zero component would represent a violation of ordinary slow roll and indicate a feature in the inflaton potential or sound speed. Although one component shows a deviation at the 98% CL, it cannot be considered statistically significant given the 20 components tested. The maximum likelihood PC parameters only improves 2DeltalnL by 17 for the 20 parameters associated with known glitches in the WMAP power spectrum at multipoles l<60. We make model-independent predictions for the matching glitches in the polarization spectrum that would test their inflationary origin. This complete analysis for bandlimited features in the source function of generalized slow roll can be used to constrain parameters of specific models of the inflaton potential without requiring a separate likelihood analysis for each choice. We illustrate its use by placing bounds on the height and width of a step-like feature in the potential proposed to explain the glitch at 20<l<40. Even allowing for the presence of features in the temperature spectrum, single field inflation makes sharp falsifiable predictions for the acoustic peaks in the polarization whose violation would require extra degrees of freedom.

We place functional constraints on the shape of the inflaton potential from the cosmic microwave background through a variant of the generalized slow roll approximation that allows large amplitude, rapidly changing deviations from scale-free conditions. Employing a principal component decomposition of the source function G'~3(V'/V)^2 - 2V''/V and keeping only those measured to better than 10% results in 5 nearly independent Gaussian constraints that maybe used to test any single-field inflationary model where such deviations are expected. The first component implies < 3% variations at the 100 Mpc scale. One component shows a 95% CL preference for deviations around the 300 Mpc scale at the ~10% level but the global significance is reduced considering the 5 components examined. This deviation also requires a change in the cold dark matter density which in a flat LCDM model is disfavored by current supernova and Hubble constant data and can be tested with future polarization or high multipole temperature data. Its impact resembles a local running of the tilt from multipoles 30-800 but is only marginally consistent with a constant running beyond this range. For this analysis, we have implemented a ~40x faster WMAP7 likelihood method which we have made publicly available.

We develop a variant of the generalized slow roll approach for calculating the curvature power spectrum that is well-suited for order unity deviations in power caused by sharp features in the inflaton potential. As an example, we show that predictions for a step function potential, which has been proposed to explain order unity glitches in the CMB temperature power spectrum at multipoles l=20-40, are accurate at the percent level. Our analysis shows that to good approximation there is a single source function that is responsible for observable features and that this function is simply related to the local slope and curvature of the inflaton potential. These properties should make the generalized slow roll approximation useful for inflation-model independent studies of features, both large and small, in the observable power spectra.

B-modes in CMB polarization from patchy reionization arise from two effects: generation of polarization from scattering of quadrupole moments by reionization bubbles, and fluctuations in the screening of E-modes from recombination. The scattering contribution has been studied previously, but the screening contribution has not yet been calculated. We show that on scales smaller than the acoustic scale (l>300), the B-mode power from screening is larger than the B-mode power from scattering. The ratio approaches a constant ~2.5 below the damping scale (l>2000). On degree scales relevant for gravitational waves (l<100), screening B-modes have a white noise tail and are subdominant to the scattering effect. These results are robust to uncertainties in the modeling of patchy reionization.

The angular power spectrum of the cosmic microwave background temperature anisotropy observed by WMAP has an anomalous dip at l~20 and bump at l~40. One explanation for this structure is the presence of features in the primordial curvature power spectrum, possibly caused by a step in the inflationary potential. The detection of these features is only marginally significant from temperature data alone. However, the inflationary feature hypothesis predicts a specific shape for the E-mode polarization power spectrum with a structure similar to that observed in temperature at l~20-40. Measurement of the CMB polarization on few-degree scales can therefore be used as a consistency check of the hypothesis. The Planck satellite has the statistical sensitivity to confirm or rule out the model that best fits the temperature features with 3 sigma significance, assuming all other parameters are known. With a cosmic variance limited experiment, this significance improves to 8 sigma. For tests of inflationary models that can explain both the dip and bump in temperature, the primary source of uncertainty is confusion with polarization features created by a complex reionization history, which at most reduces the significance to 2.5 sigma for Planck and 5-6 sigma for an ideal experiment. Smoothing of the polarization spectrum by a large tensor component only slightly reduces the ability of polarization to test for inflationary features, as does requiring that polarization is consistent with the observed temperature spectrum given the expected low level of TE correlation on few-degree scales. A future polarization satellite would enable a decisive test of the feature hypothesis and provide complementary information about the shape of a possible step in the inflationary potential. (Abridged.)

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

We introduce a new statistical technique for extracting the inhomogeneous reionization signal from future high-sensitivity measurements of the cosmic microwave background temperature and polarization fields. If reionization is inhomogeneous, then the optical depth to recombination will be a function tau(n) of position on the sky. Anisotropies in tau(n) alter the statistics of the observed CMB via several physical mechanisms: screening of the surface of last scattering, generation of new polarization via Thomson scattering from reionization bubbles, and the kinetic Sunyaev-Zel'dovich effect. We construct a quadratic estimator for the modes of the tau field. This estimator separates the patchy reionization signal from the CMB in the form of a noisy map, which can be cross-correlated with other probes of reionization or used as a standalone probe. A future satellite experiment with sufficient sensitivity and resolution to measure the lensed B-mode on most of the sky can constrain key parameters of patchy reionization, such as the duration of the patchy epoch or the mean bubble radius, at the ~10% level.

Anomalies in the large-scale CMB temperature sky measured by WMAP have been suggested as possible evidence for a violation of statistical isotropy on large scales. In any physical model for broken isotropy, there are testable consequences for the CMB polarization field. We develop simulation tools for predicting the polarization field in models that break statistical isotropy locally through a modulation field. We study two different models: dipolar modulation, invoked to explain the asymmetry in power between northern and southern ecliptic hemispheres, and quadrupolar modulation, posited to explain the alignments between the quadrupole and octopole. For the dipolar case, we show that predictions for the correlation between the first 10 multipoles of the temperature and polarization fields can typically be tested at better than the 98% CL. For the quadrupolar case, we show that the polarization quadrupole and octopole should be moderately aligned. Such an alignment is a generic prediction of explanations which involve the temperature field at recombination and thus discriminate against explanations involving foregrounds or local secondary anisotropy. Predicted correlations between temperature and polarization multipoles out to l = 5 provide tests at the ~ 99% CL or stronger for quadrupolar models that make the temperature alignment more than a few percent likely. As predictions of anomaly models, polarization statistics move beyond the a posteriori inferences that currently dominate the field.