We propose two effective parameters that fully characterise galactic-scale structure formation at high redshifts ($z\gtrsim5$) for a variety of dark matter (DM) models that have a primordial cutoff in the matter power spectrum. Our description is within the recently proposed ETHOS framework and includes standard thermal Warm DM (WDM) and models with dark acoustic oscillations (DAOs). To define and explore this parameter space, we use high-redshift zoom-in simulations that cover a wide range of non-linear scales from those where DM should behave as CDM ($k\sim10\,h\,{\rm Mpc}^{-1}$), down to those characterised by the onset of galaxy formation ($k\sim500\,h\,{\rm Mpc}^{-1}$). We show that the two physically motivated parameters $h_{\rm peak}$ and $k_{\rm peak}$, the amplitude and scale of the first DAO peak, respectively, are sufficient to parametrize the linear matter power spectrum and classify the DM models as belonging to effective non-linear structure formation regions. These are defined by their relative departure from Cold DM ($k_{\rm peak}\rightarrow\infty$) and WDM ($h_{\rm peak}=0$) according to the non-linear matter power spectrum and halo mass function. We identify a region where the DAOs still leave a distinct signature from WDM down to $z=5$, while a large part of the DAO parameter space is shown to be degenerate with WDM. Our framework can then be used to seamlessly connect a broad class of particle DM models to their structure formation properties at high redshift without the need of additional $N$-body simulations.

Cosmological measurements of the 21-cm line of neutral hydrogen are poised to dramatically enhance our understanding of the early universe. In particular, both the epochs of reionization and cosmic dawn remain largely uncharted, and the 21-cm signal is one of the few probes to reach them. The simplest 21-cm measurement is the global signal (GS), which corresponds to the averaged absorption or emission of 21-cm photons across the entire sky. While bright radio foregrounds swamp the cosmic signal over the entire frequency range observable, presenting a formidable hurdle, they can in principle be subtracted, given enough sensitivity. Here, however, we point out an additional---and irreducible---source of uncertainty for the 21-cm GS: cosmic variance. The cosmic-variance noise arises from the finite volume of the universe accessible to 21-cm experiments. Due to the cosmological redshifting of 21-cm photons, each observed frequency probes our universe during a particular cosmic age, corresponding to a narrow redshift slice. The presence of large 21-cm fluctuations makes the GS within each slice different than the GS averaged over the entire universe. We estimate the size of this cosmic-variance noise, and find that for a standard scenario it has a size of $\sim 0.1$ mK, which is $\sim 10\%$ of the size of the expected instrumental noise of a year-long experiment. Interestingly, cosmic variance can overtake instrumental noise for scenarios with extreme 21-cm fluctuations, such as those suggested to explain the sharpness of the claimed EDGES detection. Moreover, as large-scale 21-cm fluctuations are coherent over long distances, cosmic variance correlates the measurements of the GS at nearby redshifts, leading to off-diagonal uncertainties that have so far been neglected.

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.

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

In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging.

Gravitational lensing has emerged as a powerful probe of the matter distribution on subgalactic scales, which itself may contain important clues about the fundamental origins and properties of dark matter. Broadly speaking, two different approaches have been taken in the literature to map the small-scale structure of the Universe using strong lensing, with one focused on measuring the position and mass of a small number of discrete massive subhalos appearing close in projection to lensed images, and the other focused on detecting the collective effect of all the small-scale structure between the lensed source and the observer. In this paper, we follow the latter approach and perform a detailed study of the sensitivity of galaxy-scale gravitational lenses to the ensemble properties of small-scale structure. As in some previous studies, we adopt the language of the substructure power spectrum to characterize the statistical properties of the small-scale density field. We present a comprehensive theory that treats lenses with extended sources as well as those with time-dependent compact sources (such as quasars) in a unified framework for the first time. Our approach uses mode functions to provide both computational advantages and insights about couplings between the lens and source. The goal of this paper is to develop the theory and gain the intuition necessary to understand how the sensitivity to the substructure power spectrum depends on the source and lens properties, with the eventual aim of identifying the most promising targets for such studies.

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.

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

We perform a series of cosmological hydrodynamic simulations to investigate the effects of non-gravitational dark matter (DM) interactions on the intergalactic medium (IGM). In particular, we use the ETHOS framework (Cyr-Racine et al. 2016; Vogelsberger et al. 2016) to compare statistics of the Lyman-$\alpha$ forest in cold dark matter (CDM) with an alternative model in which the DM couples strongly with a relativistic species in the early universe. These models are characterised by a cutoff in the linear power spectrum, followed by a series of 'dark acoustic oscillations' (DAOs) on sub-dwarf scales. While the primordial cutoff delays the formation of the first galaxies, structure builds-up more rapidly in the interacting DM model compared to CDM. We show that although DAOs are quickly washed away in the non-linear clustering of DM at $z\lesssim10$, their signature can be imprinted prominently in the Lyman-$\alpha$ flux power spectrum at $z>5$. On scales larger than the cutoff ($k\sim0.08$ s/km for the specific model considered here), the relative difference to CDM is reminiscent of a warm dark matter (WDM) model with a similar initial cutoff; however, the redshift evolution on smaller scales is distinctly different. The appearance and disappearance of DAOs in the Lyman-$\alpha$ flux spectrum provides a powerful way to distinguish interacting DM models from WDM and, indeed, variations in the thermal history of the IGM that may also induce a small-scale cutoff.

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/ .

Of the many proposed extensions to the $\Lambda$CDM paradigm, a model in which neutrinos self-interact until close to the epoch of matter-radiation equality has been shown to provide a good fit to current cosmic microwave background (CMB) data, while at the same time alleviating tensions with late-time measurements of the expansion rate and matter fluctuation amplitude. Interestingly, CMB fits to this model either pick out a specific large value of the neutrino interaction strength, or are consistent with the extremely weak neutrino interaction found in $\Lambda$CDM, resulting in a bimodal posterior distribution for the neutrino self-interaction cross section. In this paper, we explore why current cosmological data select this particular large neutrino self-interaction strength, and by consequence, disfavor intermediate values of the self-interaction cross section. We show how it is the $\ell \gtrsim 1000$ CMB temperature anisotropies, most recently measured by the Planck satellite, that produce this bimodality. We also establish that smaller scale temperature data, and improved polarization data measuring the temperature-polarization cross-correlation, will best constrain the neutrino self-interaction strength. We forecast that the upcoming Simons Observatory should be capable of distinguishing between the models.

New physics in the neutrino sector might be necessary to address anomalies between different neutrino oscillation experiments. Intriguingly, it also offers a possible solution to the discrepant cosmological measurements of $H_0$ and $\sigma_8$. We show here that delaying the onset of neutrino free-streaming until close to the epoch of matter-radiation equality can naturally accommodate a larger value for the Hubble constant $H_0=72.3 \pm 1.4$ km/s/Mpc and a lower value of the matter fluctuations $\sigma_8=0.786\pm 0.020$, while not degrading the fit to the cosmic microwave background (CMB) damping tail. We achieve this by introducing neutrino self-interactions in the presence of a non-vanishing sum of neutrino masses. This strongly interacting neutrino cosmology prefers $N_{\rm eff} = 4.02 \pm 0.29$, which has interesting implications for particle model-building and neutrino oscillation anomalies. We show that the absence of the neutrino free-streaming phase shift on the CMB can be compensated by shifting the value of other cosmological parameters, hence providing an important caveat to the detections made in the literature. Due to their impact on the evolution of the gravitational potential at early times, self-interacting neutrinos and their subsequent decoupling leave a rich structure on the matter power spectrum. In particular, we point out the existence of a novel localized feature appearing on scales entering the horizon at the onset of neutrino free-streaming. While the interacting neutrino cosmology provides a better global fit to current cosmological data, we find that traditional Bayesian analyses penalize the model as compared to the standard cosmological. Our analysis shows that it is possible to find radically different cosmological models that nonetheless provide excellent fits to the data, hence providing an impetus to thoroughly explore alternate cosmological scenarios.

Local measurements of the Hubble parameter are increasingly in tension with the value inferred from a $\Lambda$CDM fit to the cosmic microwave background (CMB) data. In this paper, we construct scenarios in which evolving scalar fields significantly ease this tension by adding energy to the Universe around recombination in a narrow redshift window. We identify solutions of $V \propto \phi^{2 n}$ with simple asymptotic behavior, both oscillatory (rocking) and rolling. These are the first solutions of this kind in which the field evolution and fluctuations are consistently implemented using the equations of motion. Our findings differ qualitatively from those of the existing literature, which rely upon a coarse-grained fluid description. Combining CMB data with low-redshift measurements, the best fit model has $n=2$ and increases the allowed value of $H_0$ from 69.2 km/s/Mpc in $\Lambda$CDM to 72.3 km/s/Mpc at $2\sigma$. Future measurements of the late-time amplitude of matter fluctuations and of the reionization history could help distinguish these models from competing solutions.

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

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

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

For nearly 40 years, dark matter has been widely assumed to be cold and collisionless. Cold dark matter models make fundamental predictions for the behavior of dark matter on small (<10 kpc) scales. These predictions include cuspy density profiles at the centers of dark matter halos and a halo mass function that increases as dN/dM ~ M^-1.9 down to very small masses. We suggest two observational programs relying on extremely large telescopes to critically test these predictions, and thus shed new light on the nature of dark matter. (1) Combining adaptive optics-enabled imaging with deep spectroscopy to measure the three-dimensional motions of stars within a sample of Local Group dwarf galaxies that are the cleanest dark matter laboratories known in the nearby universe. From these observations the inner slope of the dark matter density profile can be determined with an accuracy of 0.20 dex, enabling a central cusp to be distinguished from a core at 5 sigma significance. (2) Diffraction-limited AO imaging and integral field spectroscopy of gravitationally lensed galaxies and quasars to quantify the abundance of dark substructures in the halos of the lens galaxies and along the line of sight. Observations of 50 lensed arcs and 50 multiply-imaged quasars will be sufficient to measure the halo mass function over the range 10^7 < M < 10^10 Msun at cosmological scales, independent of the baryonic and stellar composition of those structures. These two observational probes provide complementary information about the small scale structure, with a joint self-consistent analysis mitigating limitations of either probe. This program will produce the strongest existing constraints on the properties of dark matter on small scales, allowing conclusive tests of alternative warm, fuzzy, and self-interacting dark matter models.

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.

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.

We contrast predictions for the high-redshift galaxy population and reionization history between cold dark matter (CDM) and an alternative self-interacting dark matter model based on the recently developed ETHOS framework that alleviates the small-scale CDM challenges within the Local Group. We perform the highest resolution hydrodynamical cosmological simulations (a 36~Mpc$^3$ volume with gas cell mass of $\sim10^5\mathrm{M}_{\odot}$ and minimum gas softening of $\sim180$~pc) within ETHOS to date -- plus a CDM counterpart -- to quantify the abundance of galaxies at high redshift and their impact on reionization. We find that ETHOS predicts galaxies with higher ultraviolet (UV) luminosities than their CDM counterparts and a faster build-up of the faint end of the UV luminosity function. These effects, however, make the optical depth to reionization less sensitive to the power spectrum cut-off: the ETHOS model differs from the CDM $\tau$ value by only 10 per cent and is consistent with Planck limits if the effective escape fraction of UV photons is 0.1-0.5. We conclude that current observations of high-redshift luminosity functions cannot differentiate between ETHOS and CDM models, but deep JWST surveys of strongly-lensed, inherently faint galaxies have the potential to test non-CDM models that offer attractive solutions to CDM's Local Group problems.

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.

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta.

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.

We present updated constraints on the free-streaming nature of cosmological neutrinos from cosmic microwave background (CMB) power spectra, baryonic acoustic oscillation data, and local measurements of the Hubble constant. Specifically, we consider a Fermi-like four-fermion interaction between massless neutrinos, characterized by an effective coupling constant $ G_{\rm eff}$, and resulting in a neutrino opacity $\dot{\tau}_\nu\propto G_{\rm eff}^2 T_\nu^5$. Using a conservative prior on the parameter $\log_{10}\left(G_{\rm eff} {\rm MeV}^2\right)$, we find a bimodal posterior distribution. The first of these modes is consistent with the standard $\Lambda$CDM cosmology and corresponds to neutrinos decoupling at redshift $z_{\nu,{\rm dec}} > 1.3\times10^5$. The other mode of the posterior, dubbed the "interacting neutrino mode", corresponds to neutrino decoupling occurring within a narrow redshift window centered around $z_{\nu,{\rm dec}}\sim8300$. This mode is characterized by a high value of the effective neutrino coupling constant, together with a lower value of the scalar spectral index and amplitude of fluctuations, and a higher value of the Hubble parameter. Using both a maximum likelihood analysis and the ratio of the two mode's Bayesian evidence, we find the interacting neutrino mode to be statistically disfavored compared to the standard $\Lambda$CDM cosmology. Interestingly, the addition of CMB polarization and direct Hubble constant measurements significantly raises the statistical significance of this secondary mode, indicating that new physics in the neutrino sector could help explain the difference between local measurements of $H_0$, and those inferred from CMB data. A robust consequence of our results is that neutrinos must be free streaming long before the epoch of matter-radiation equality.

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.

We formulate an effective theory of structure formation (ETHOS) that enables cosmological structure formation to be computed in almost any microphysical model of dark matter physics. This framework maps the detailed microphysical theories of particle dark matter interactions into the physical effective parameters that shape the linear matter power spectrum and the self-interaction transfer cross section of non-relativistic dark matter. These are the input to structure formation simulations, which follow the evolution of the cosmological and galactic dark matter distributions. Models with similar effective parameters in ETHOS but with different dark particle physics would nevertheless result in similar dark matter distributions. We present a general method to map an ultraviolet complete or effective field theory of low energy dark matter physics into parameters that affect the linear matter power spectrum and carry out this mapping for several representative particle models. We further propose a simple but useful choice for characterizing the dark matter self-interaction transfer cross section that parametrizes self-scattering in structure formation simulations. Taken together, these effective parameters in ETHOS allow the classification of dark matter theories according to their structure formation properties rather than their intrinsic particle properties, paving the way for future simulations to span the space of viable dark matter physics relevant for structure formation.

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.

Recently it was shown that dark matter with mass of order the weak scale can be charged under a new long-range force, decoupled from the Standard Model, with only weak constraints from early Universe cosmology. Here we consider the implications of an additional charged particle $C$ that is light enough to lead to significant dissipative dynamics on galactic times scales. We highlight several novel features of this model, which can be relevant even when the $C$ particle constitutes only a small fraction of the number density (and energy density). We assume a small asymmetric abundance of the $C$ particle whose charge is compensated by a heavy $X$ particle so that the relic abundance of dark matter consists mostly of symmetric $X$ and $\bar{X}$, with a small asymmetric component made up of $X$ and $C$. As the universe cools, it undergoes asymmetric recombination binding the free $C$s into $(XC)$ dark atoms efficiently. Even with a tiny asymmetric component, the presence of $C$ particles catalyzes tight coupling between the heavy dark matter $X$ and the dark photon plasma that can lead to a significant suppression of the matter power spectrum on small scales and lead to some of the strongest bounds on such dark matter theories. We find a viable parameter space where structure formation constraints are satisfied and significant dissipative dynamics can occur in galactic haloes but show a large region is excluded. Our model shows that subdominant components in the dark sector can dramatically affect structure formation.

We revisit constraints on dark matter that is charged under a $U(1)$ gauge group in the dark sector, decoupled from Standard Model forces. We find that the strongest constraints in the literature are subject to a number of mitigating factors. For instance, the naive dark matter thermalization timescale in halos is corrected by saturation effects that slow down isotropization for modest ellipticities. The weakened bounds uncover interesting parameter space, making models with weak-scale charged dark matter viable, even with electromagnetic strength interaction. This also leads to the intriguing possibility that dark matter self-interactions within small dwarf galaxies are extremely large, a relatively unexplored regime in current simulations. Such strong interactions suppress heat transfer over scales larger than the dark matter mean free path, inducing a dynamical cutoff length scale above which the system appears to have only feeble interactions. These effects must be taken into account to assess the viability of darkly-charged dark matter. Future analyses and measurements should probe a promising region of parameter space for this model.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies.

We present the first simulations within an effective theory of structure formation (ETHOS), which includes the effect of interactions between dark matter and dark radiation on the linear initial power spectrum and dark matter self-interactions during non-linear structure formation. We simulate a Milky Way-like halo in four different dark matter models and the cold dark matter case. Our highest resolution simulation has a particle mass of $2.8\times 10^4\,{\rm M}_\odot$ and a softening length of $72.4\,{\rm pc}$. We demonstrate that all alternative models have only a negligible impact on large scale structure formation. On galactic scales, however, the models significantly affect the structure and abundance of subhaloes due to the combined effects of small scale primordial damping in the power spectrum and late time self-interactions. We derive an analytic mapping from the primordial damping scale in the power spectrum to the cutoff scale in the halo mass function and the kinetic decoupling temperature. We demonstrate that certain models within this extended effective framework that can alleviate the too-big-to-fail and missing satellite problems simultaneously, and possibly the core-cusp problem. The primordial power spectrum cutoff of our models naturally creates a diversity in the circular velocity profiles, which is larger than that found for cold dark matter simulations. We show that the parameter space of models can be constrained by contrasting model predictions to astrophysical observations. For example, some models may be challenged by the missing satellite problem if baryonic processes were to be included and even over-solve the too-big-to-fail problem; thus ruling them out.

We perform a detailed study of the weak interactions of standard model neutrinos with the primordial plasma and their effect on the resonant production of sterile neutrino dark matter. Motivated by issues in cosmological structure formation on small scales, and reported X-ray signals that could be due to sterile neutrino decay, we consider $7$ keV-scale sterile neutrinos. Oscillation-driven production of such sterile neutrinos occurs at temperatures $T \gtrsim 100$ MeV, where we study two significant effects of weakly charged species in the primordial plasma: (1) the redistribution of an input lepton asymmetry; (2) the opacity for active neutrinos. We calculate the redistribution analytically above and below the quark-hadron transition, and match with lattice QCD calculations through the transition. We estimate opacities due to tree level processes involving leptons and quarks above the quark-hadron transition, and the most important mesons below the transition. We report final sterile neutrino dark matter phase space densities that are significantly influenced by these effects, and yet relatively robust to remaining uncertainties in the nature of the quark-hadron transition. We also provide transfer functions for cosmological density fluctuations with cutoffs at $k \simeq 10 \ h \ {\rm Mpc}^{-1}$, that are relevant to galactic structure formation.

In the standard model neutrinos are assumed to have streamed across the Universe since they last scattered at the weak decoupling epoch when the temperature of the standard-model plasma was ~MeV. The shear stress of free-streaming neutrinos imprints itself gravitationally on the Cosmic Microwave Background (CMB) and makes the CMB a sensitive probe of neutrino scattering. Yet, the presence of nonstandard physics in the neutrino sector may alter this standard chronology and delay neutrino free-streaming until a much later epoch. We use observations of the CMB to constrain the strength of neutrino self-interactions G_eff and put limits on new physics in the neutrino sector from the early Universe. Recent measurements of the CMB at large multipoles made by the Planck satellite and high-l experiments are critical for probing this physics. Within the context of conventional LambdaCDM parameters cosmological data are compatible with G_eff < 1/(56 MeV)^2 and neutrino free-streaming might be delayed until their temperature has cooled to as low as ~25 eV. Intriguingly, we also find an alternative cosmology compatible with cosmological data in which neutrinos scatter off each other until z~10^4 with a preferred interaction strength in a narrow region around $G_{\rm eff} \simeq 1/({\rm 10 \, MeV})^{2} \simeq 8.6\times10^8 G_{\rm F}$, where $G_{\rm F}$ is the Fermi constant. This distinct self-interacting neutrino cosmology is characterized by somewhat lower values of both the scalar spectral index and the amplitude of primordial fluctuations. While we phrase our discussion here in terms of a specific scenario in which a late onset of neutrino free-streaming could occur, our constraints on the neutrino visibility function are very general.

We demonstrate that self-interacting dark matter models with interactions mediated by light particles can have significant deviations in the matter power-spectrum and detailed structure of galactic halos when compared to a standard cold dark matter scenario. While these deviations can take the form of suppression of small scale structure that are in some ways similar to that of warm dark matter, the self-interacting models have a much wider range of possible phenomenology. A long-range force in the dark matter can introduce multiple scales to the initial power spectrum, in the form of dark acoustic oscillations and an exponential cut-off in the power spectrum. Using simulations we show that the impact of these scales can remain observationally relevant up to the present day. Furthermore, the self-interaction can continue to modify the small-scale structure of the dark matter halos, reducing their central densities and creating a dark matter core. The resulting phenomenology is unique to this type of models.

If all or a fraction of the dark matter (DM) were coupled to a bath of dark radiation (DR) in the early Universe we expect the combined DM-DR system to give rise to acoustic oscillations of the dark matter until it decouples from the DR. Much like the standard baryon acoustic oscillations, these dark acoustic oscillations (DAO) imprint a characteristic scale, the sound horizon of dark matter, on the matter power spectrum. We compute in detail how the microphysics of the DM-DR interaction affects the clustering of matter in the Universe and show that the DAO physics also gives rise to unique signatures in the temperature and polarization spectra of the cosmic microwave background (CMB). We use cosmological data from the CMB, baryon acoustic oscillations (BAO), and large-scale structure to constrain the possible fraction of interacting DM as well as the strength of its interaction with DR. Like nearly all knowledge we have gleaned about dark matter since inferring its existence this constraint rests on the betrayal by gravity of the location of otherwise invisible dark matter. Although our results can be straightforwardly applied to a broad class of models that couple dark matter particles to various light relativistic species, in order to make quantitative predictions, we model the interacting component as dark atoms coupled to a bath of dark photons. We find that linear cosmological data and CMB lensing put strong constraints on existence of DAO features in the CMB and the large-scale structure of the Universe. Interestingly, we find that at most ~5% of all DM can be very strongly interacting with DR. We show that our results are surprisingly constraining for the recently proposed Double-disk DM model, a novel example of how large-scale precision cosmological data can be used to constrain galactic physics and sub-galactic structure.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

While, to ensure successful cosmology, dark matter (DM) must kinematically decouple from the standard model plasma very early in the history of the Universe, it can remain coupled to a bath of "dark radiation" until a relatively late epoch. One minimal theory that realizes such a scenario is the Atomic Dark Matter model, in which two fermions oppositely charged under a new U(1) dark force are initially coupled to a thermal bath of "dark photons" but eventually recombine into neutral atom-like bound states and begin forming gravitationally-bound structures. As dark atoms have (dark) atom-sized geometric cross sections, this model also provides an example of self-interacting DM with a velocity-dependent cross section. Delayed kinetic decoupling in this scenario predicts novel DM properties on small scales but retains the success of cold DM on larger scales. We calculate the atomic physics necessary to capture the thermal history of this dark sector and show significant improvements over the standard atomic hydrogen calculation are needed. We solve the Boltzmann equations that govern the evolution of cosmological fluctuations in this model and find in detail the impact of the atomic DM scenario on the matter power spectrum and the cosmic microwave background (CMB). This scenario imprints a new length scale, the Dark-Acoustic-Oscillation (DAO) scale, on the matter density field. This DAO scale shapes the small-scale matter power spectrum and determines the minimal DM halo mass at late times which may be many orders of magnitude larger than in a typical WIMP scenario. This model necessarily includes an extra dark radiation component, which may be favoured by current CMB experiments, and we quantify CMB signatures that distinguish an atomic DM scenario from a standard $\Lambda$CDM model containing extra free-streaming particles. [Abridged]

Many models of inflation predict oscillatory features in the bispectrum of primordial fluctuations. Since it has been shown that primordial non-Gaussianity can lead to a scale-dependent halo bias, we investigate the effect of oscillations in the three-point function on the clustering of dark-matter halos. Interestingly, we find that features in the inflaton potential such as oscillations or sharp steps get imprinted in the mass dependence of the non-Gaussian halo bias. In this paper, we focus on models displaying a sharp feature in the inflaton potential as well as Resonant non-Gaussianity. In both cases, we find a strong scale dependence for the non-Gaussian halo bias with a slope similar to that of the local model. In the resonant case, we find that the non-Gaussian bias oscillates with halo mass, a novel feature that is unique to this type of models. In the case of a sharp feature in the inflaton potential, we find that the clustering of halos is enhanced at the mass scale corresponding to the Fourier mode that exited the horizon when the inflaton was crossing the feature in the potential. Both of these are new effects that open the possibility of characterizing the inflationary potential with large-scale-structure surveys. We briefly discuss the prospects for detecting these non-Gaussian effects.

Prior to recombination photons, electrons, and atomic nuclei rapidly scattered and behaved, almost, like a single tightly-coupled photon-baryon plasma. We investigate here the accuracy of the tight-coupling approximation commonly used to numerically evolve the baryon and photon perturbation equations at early times. By solving the exact perturbations equations with a stiff solver starting deep in the radiation-dominated epoch we find the level of inaccuracy introduced by resorting to the standard first-order tight-coupling approximation. We develop a new second-order approximation in the inverse Thomson opacity expansion and show that it closely tracks the full solution, at essentially no extra numerical cost. We find the bias on estimates of cosmological parameters introduced by the first-order approximation is, for most parameters, negligible. Finally, we show that our second-order approximation can be used to reduce the time needed to compute cosmic microwave background angular spectra by as much as ~17%.

Reheating is an important part of inflationary cosmology. It describes the production of Standard Matter particles after the phase of accelerated expansion. We give a review of the reheating process, focusing on an in-depth discussion of the preheating stage which is characterized by exponential particle production due to a parametric resonance or tachyonic instability. We give a brief overview of the thermalization process after preheating and end with a survey of some applications to supersymmetric theories and to other issues in cosmology such as baryogenesis, dark matter and metric preheating.

It has recently been argued that the PAMELA, ATIC and PPB-BETS data showing an anomalous excess of high-energy cosmic ray positrons and electrons might be explained by dark matter annihilating in the Galactic halo with a cross section resonantly enhanced compared to its value in the primeval plasma. We find that with a very large annihilation cross section the flash of energetic photons and electron-positron pairs expected from dark-matter annihilation in the first protohalos that form at redshift z~40 is likely substantial and observable. As a consequence, bounds on the allowed energy injection into the primordial gas and the energy density of the diffuse gamma-ray background give rise to limits on the low-velocity dark matter cross section that can be difficult to reconcile with this interpretation of the PAMELA, ATIC and PPB-BETS results.

We study how the resonant decay of moduli fields arising in the Minimal Supersymmetric Standard Model (MSSM) could affect large scale curvature perturbations in the early universe. It has been known for some time that the presence of entropy perturbations in a multi-component system can act as seeds for the curvature perturbations on all scales. These entropy perturbations could be amplified exponentially if one of the moduli decays via stochastic resonance, affecting the curvature power spectrum in the process. By imposing the COBE normalization on this power spectrum, one could put constraints on the masses and couplings of the underlying particle physics model without having to rely on collider experiments. We discuss in detail the case of the MSSM but this method could be applied to other theories beyond the Standard Model.