Warm dark matter has recently become increasingly constrained by observational inferences about the low-mass end of the subhalo mass function, which would be suppressed by dark matter free streaming in the early Universe. In this work, we point out that a constraint can be placed on ultralight bosonic dark matter (often referred to as "fuzzy dark matter") based on similar considerations. Recent limits on warm dark matter from strong gravitational lensing of quasars and from fluctuations in stellar streams separately translate to a lower limit of $\sim 2.1 \times 10^{-21}$ eV on the mass of an ultralight boson comprising all dark matter. These limits are complementary to constraints on ultralight dark matter from the Lyman-$\alpha$ forest and are subject to a completely different set of assumptions and systematic uncertainties. Taken together, these probes strongly suggest that dark matter with a mass $\sim 10^{-22}$ eV is not a viable way to reconcile differences between cold dark matter simulations and observations of structure on small scales.

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

The twin Higgs mechanism is a solution to the little hierarchy problem in which the top partner is neutral under the Standard Model (SM) gauge group. The simplest mirror twin Higgs (MTH) model -- where a $\mathbf{Z}_2$ symmetry copies each SM particle -- has too many relativistic degrees of freedom to be consistent with cosmological observations. We demonstrate that MTH models can have an observationally viable cosmology if the twin mass spectrum leads to twin neutrino decoupling before the SM and twin QCD phase transitions. Our solution requires the twin photon to have a mass of $\sim 20$ MeV and kinetically mix with the SM photon to mediate entropy transfer from the twin sector to the SM. This twin photon can be robustly discovered or excluded by future experiments. Additionally, the residual twin degrees of freedom present in the early Universe in this scenario would be detectable by future observations of the cosmic microwave background.

Self-interacting dark matter provides a promising alternative for the cold dark matter paradigm to solve potential small-scale galaxy formation problems. Nearly all self-interacting dark matter simulations so far have considered only elastic collisions. Here we present simulations of a galactic halo within a generic inelastic model using a novel numerical implementation in the Arepo code to study arbitrary multi-state inelastic dark matter scenarios. For this model we find that inelastic self-interactions can: (i) create larger subhalo density cores compared to elastic models for the same cross section normalisation; (ii) lower the abundance of satellites without the need for a power spectrum cutoff; (iii) reduce the total halo mass by about 10%; (iv) inject the energy equivalent of O(100) million Type II supernovae in galactic haloes through level de-excitation; (v) avoid the gravothermal catastrophe due to removal of particles from halo centers. We conclude that a ~5 times larger elastic cross section is required to achieve the same central density reduction as the inelastic model. This implies that well-established constraints on self-interacting cross sections have to be revised if inelastic collisions are the dominant mode. In this case significantly smaller cross sections can achieve the same core density reduction thereby increasing the parameter space of allowed models considerably.

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

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.

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.

Dark matter could be a thermal relic comprised of strongly interacting massive particles (SIMPs), where $3 \rightarrow 2$ interactions set the relic abundance. Such interactions generically arise in theories of chiral symmetry breaking via the Wess-Zumino-Witten term. In this work, we show that an axion-like particle can successfully maintain kinetic equilibrium between the dark matter and the visible sector, allowing the requisite entropy transfer that is crucial for SIMPs to be a cold dark matter candidate. Constraints on this scenario arise from beam dump and collider experiments, from the cosmic microwave background, and from supernovae. We find a viable parameter space when the axion-like particle is close in mass to the SIMP dark matter, with strong-scale masses of order a few hundred MeV. Many planned experiments are set to probe the parameter space in the near future.

If a component of the dark matter has dissipative interactions, it could collapse to form a thin dark disk in our Galaxy that is coplanar with the baryonic disk. It has been suggested that dark disks could explain a variety of observed phenomena, including periodic comet impacts. Using the first data release from the $Gaia$ space observatory, we search for a dark disk via its effect on stellar kinematics in the Milky Way. Our new limits disfavor the presence of a thin dark matter disk, and we present updated measurements on the total matter density in the Solar neighborhood.

The recent discovery of gravitational waves from merging black holes has generated interest in primordial black holes as a possible component of the dark matter. In this paper, we show that pulsar timing may soon have sufficient data to constrain $1$-$1000\,M_{\odot}$ primordial black holes via the non-detection of a third-order Shapiro time delay as the black holes move around the Galactic halo. We present the results of a Monte Carlo simulation which suggests that future data from known pulsars may be capable of constraining the PBH density more stringently than other existing methods in the mass range ~1-30$\,M_{\odot}$. We find that timing new pulsars discovered using the proposed Square Kilometre Array may constrain primordial black holes in this mass range to comprise less than ~1-10% of the dark matter.

We present a new method for interferometric imaging that is ideal for the large fields of view and compact arrays common in 21 cm cosmology. We first demonstrate the method with simulations for two very different low frequency interferometers, the Murchison Widefield Array (MWA) and the MIT Epoch of Reionization (MITEoR) Experiment. We then apply the method to the MITEoR data set collected in July 2013 to obtain the first northern sky map from 128 MHz to 175 MHz at about 2 degree resolution, and find an overall spectral index of -2.73+/-0.11. The success of this imaging method bodes well for upcoming compact redundant low-frequency arrays such as HERA. Both the MITEoR interferometric data and the 150 MHz sky map are publicly available at http://space.mit.edu/home/tegmark/omniscope.html.

We show that a two-excitation process in superfluid helium, combined with sensitivity to meV energy depositions, can probe dark matter down to the ~keV warm dark matter mass limit. This mass reach is three orders of magnitude below what can be probed with ordinary nuclear recoils in helium at the same energy resolution. For dark matter lighter than $\sim 100$ keV, the kinematics of the process requires the two athermal excitations to have nearly equal and opposite momentum, potentially providing a built-in coincidence mechanism for controlling backgrounds.

We compute the non-Gaussian contribution to the covariance of the matter power spectrum at one-loop order in Standard Perturbation Theory (SPT), and using the framework of the effective field theory (EFT) of large scale structure (LSS). The complete one-loop contributions are evaluated for the first time, including the leading EFT corrections that involve seven independent operators, of which four appear in the power spectrum and bispectrum. We compare the non-Gaussian part of the one-loop covariance computed with both SPT and EFT of LSS to two separate simulations. In one simulation, we find that the one-loop prediction from SPT reproduces the simulation well to $k_i + k_j \sim$ 0.25 h/Mpc, while in the other simulation we find a substantial improvement of EFT of LSS (with one free parameter) over SPT, more than doubling the range of $k$ where the theory accurately reproduces the simulation. The disagreement between these two simulations points to unaccounted for systematics, highlighting the need for improved numerical and analytic understanding of the covariance.

Pulsar timing arrays (PTAs) are placing increasingly stringent constraints on the strain amplitude of continuous gravitational waves emitted by supermassive black hole binaries on subparsec scales. In this paper, we incorporate independent information about the dynamical masses $M_{bh}$ of supermassive black holes in specific galaxies at known distances and use this additional information to further constrain whether or not those galaxies could host a detectable supermassive black hole binary. We estimate the strain amplitudes from individual binaries as a function of binary mass ratio for two samples of nearby galaxies: (1) those with direct dynamical measurements of $M_{bh}$ in the literature, and (2) the 116 most massive early-type galaxies (and thus likely hosts of the most massive black holes) within 108 Mpc from the MASSIVE Survey. Our exploratory analysis shows that the current PTA upper limits on continuous waves (as a function of angular position in the sky) can already constrain the mass ratios of hypothetical black hole binaries in many galaxies in our samples. The constraints are stronger for galaxies with larger $M_{bh}$ and at smaller distances. For the black holes with $M_{bh} \gtrsim 5\times 10^9 M_\odot$ at the centers of NGC 1600, NGC 4889, NGC 4486 (M87) and NGC 4649 (M60), any binary companion in orbit within the PTA frequency bands would have to have a mass ratio of a few percent or less.

We compute the connected four point correlation function (the trispectrum in Fourier space) of cosmological density perturbations at one-loop order in Standard Perturbation Theory (SPT) and the Effective Field Theory of Large Scale Structure (EFT of LSS). This paper is a companion to our earlier work on the non-Gaussian covariance of the matter power spectrum, which corresponds to a particular wavenumber configuration of the trispectrum. In the present calculation, we highlight and clarify some of the subtle aspects of the EFT framework that arise at third order in perturbation theory for general wavenumber configurations of the trispectrum. We consistently incorporate vorticity and non-locality in time into the EFT counterterms and lay out a complete basis of building blocks for the stress tensor. We show predictions for the one-loop SPT trispectrum and the EFT contributions, focusing on configurations which have particular relevance for using LSS to constrain primordial non-Gaussianity.

Self-interacting dark matter scenarios have recently attracted much attention, as a possible means to alleviate the tension between N-body simulations and observations of the dark matter distribution on galactic and sub-galactic scales. The presence of internal structure for the dark matter --- for example, a nearly-degenerate state in the spectrum that could decay, or be collisionally excited or de-excited --- has also been proposed as a possible means to address these discrepancies. Such internal structure can be a source of interesting signatures in direct and indirect dark matter searches, for example providing a novel explanation for the 3.5 keV line recently observed in galaxies and galaxy clusters. We analyze a simple model of dark matter self-scattering including a nearly-degenerate excited state, and develop an accurate analytic approximation for the elastic and inelastic s-wave cross sections, which is valid outside the perturbative regime provided the particle velocity is sufficiently low (this condition is also required for the s-wave to dominate over higher partial waves). We anticipate our results will be useful in incorporating inelastic self-scattering into N-body simulations, in order to study the quantitative impact of nearly-degenerate states in the dark matter spectrum on galactic structure and dynamics, and in computing the indirect signatures of multi-state dark matter.

We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.

Recent measurements by the {\it Planck} experiment of the power spectrum of temperature anisotropies in the cosmic microwave background radiation (CMB) reveal a deficit of power in low multipoles compared to the predictions from best-fit $\Lambda$CDM cosmology. The low-$\ell$ anomaly may be explained by the presence of primordial isocurvature perturbations in addition to the usual adiabatic spectrum, and hence may provide the first robust evidence that early-universe inflation involved more than one scalar field. In this paper we explore the production of isocurvature perturbations in nonminimally coupled two-field inflation. We find that this class of models readily produces enough power in the isocurvature modes to account for the \emph{Planck} low-$\ell$ anomaly, while also providing excellent agreement with the other {\it Planck} results.

Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.