Alternatives to the cold, collisionless dark matter (DM) paradigm in which DM behaves as a collisional fluid generically suppress small-scale structure. Herein we use the observed population of Milky Way (MW) satellite galaxies to constrain the collisional nature of DM, focusing on DM-baryon scattering. We first derive analytic upper limits on the velocity-independent DM-baryon scattering cross section by translating the upper bound on the lowest mass of halos inferred to host satellites into a characteristic cutoff scale in the linear matter power spectrum. We then confirm and improve these results through a detailed probabilistic inference of the MW satellite population that marginalizes over relevant astrophysical uncertainties. This yields $95\%$ confidence upper limits on the DM-baryon scattering cross section of $2\times10^{-29}\ \rm{cm}^2$ ($6\times 10^{-27}\ \rm{cm}^2$) for DM particle masses $m_\chi$ of~$10\ \rm{keV}$ ($10\ \rm{GeV}$); these limits scale as $m_\chi^{1/4}$ for $m_\chi \ll 1\ \rm{GeV}$ and $m_\chi$ for~$m_\chi \gg 1\ \rm{GeV}$. This analysis improves upon cosmological bounds derived from cosmic-microwave-background anisotropy measurements by multiple orders of magnitude over a wide range of DM masses, excluding regions of parameter space previously unexplored by other methods, including direct-detection experiments. Our work reveals a mapping between DM-baryon scattering and other alternative DM models, and we discuss the implications of our results for warm and fuzzy DM scenarios.

Baryon acoustic oscillations (BAO) provide a robust standard ruler, and can be used to constrain the expansion history of the Universe at low redshift. Standard BAO analyses return a model-independent measurement of the expansion rate and the angular diameter distance as function of redshift, normalized by the sound horizon at radiation drag. However, this methodology relies on anisotropic distance distortions of a fixed, pre-computed template (obtained in a given fiducial cosmology) in order to fit the observations. Therefore, it may be possible that extensions to the consensus $\Lambda$CDM add contributions to the BAO feature that cannot be captured by the template fitting. We perform mock BAO fits to power spectra computed assuming cosmological models which modify the growth of perturbations prior to recombination in order to test the robustness of the standard BAO analysis. We find no significant bias in the BAO analysis for the models under study ($\Lambda$CDM with a free effective number of relativistic species, early dark energy, and a model with interactions between neutrinos and a fraction of the dark matter), even for cases which do not provide a good fit to Planck measurements of the CMB power spectra. This result supports the use of the standard BAO analysis and its measurements to perform cosmological parameter inference and to constrain exotic models. In addition, we provide a methodology to reproduce our study for different models and surveys, as well as discuss different options to handle eventual biases in the BAO measurements.

We use a semianalytic approach that is calibrated to N-body simulations to study the evolution of self-interacting dark matter cores in galaxies. We demarcate the regime where the temporal evolution of the core density follows a well-defined track set by the initial halo parameters and the cross section. Along this track, the central density reaches a minimum value set by the initial halo density. Further evolution leads to an outward heat transfer, inducing gravothermal core collapse such that the core shrinks as its density increases. We show that the time scale for the core collapse is highly sensitive to the outer radial density profile. Satellite galaxies with significant mass loss due to tidal stripping should have larger central densities and significantly faster core collapse compared to isolated halos. Such a scenario could explain the dense and compact cores of dwarf galaxies in the Local Group like Tucana (isolated from the Milky Way), the classical Milky Way satellite Draco, and some of the ultrafaint satellites. If the ultimate fate of core collapse is black hole formation, then the accelerated time scale provides a new mechanism for creating intermediate mass black holes.

We present the Model-Agnostic Dark Halo Analysis Tool (MADHAT), a numerical tool which implements a Fermi-LAT data-driven, model-independent analysis of gamma-ray emission from dwarf satellite galaxies and dwarf galaxy candidates due to dark matter annihilation, dark matter decay, or other nonstandard or unknown astrophysics. This tool efficiently provides statistical upper bounds on the number of observed photons in excess of the number expected, based on empirical determinations of foregrounds and backgrounds, using a stacked analysis of any selected set of dwarf targets. It also calculates the resulting bounds on the properties of dark matter under any assumptions the user makes regarding dark sector particle physics or astrophysics. As an application, we determine new bounds on Sommerfeld-enhanced dark matter annihilation in a set of eight dwarfs. MADHAT v1.0 includes 58 dwarfs and dwarf candidate targets, and we discuss future planned developments. MADHAT is available and will be maintained at https://github.com/MADHATdm

We calculate the effective $J$-factors, which determine the strength of indirect detection signals from dark matter annihilation, for 25 dwarf spheroidal galaxies (dSphs). We consider several well-motivated assumptions for the relative velocity dependence of the dark matter annihilation cross section: $\sigma_A v$: $s$-wave (velocity independent), $p$-wave ($\sigma_A v \propto v^2$), $d$-wave ($\sigma_A v \propto v^4$), and Sommerfeld-enhancement in the Coulomb limit ($\sigma_A v \propto 1/v$). As a result we provide the largest and most updated sample of J-factors for velocity-dependent annihilation models. For each scenario, we use Fermi-LAT gamma-ray data to constrain the annihilation cross section. Due to the assumptions made in our gamma-ray data analysis, our bounds are comparable to previous bounds on both the $p$-wave and Sommerfeld-enhanced cross sections using dSphs. Our bounds on the $d$-wave cross section are the first such bounds using indirect detection data.

We present the effective $J$-factors for the Milky Way for scenarios in which dark matter annihilation is p-wave or d-wave suppressed. We find that the velocity suppression of dark matter annihilation can have a sizable effect on the morphology of a potential dark matter annihilation signal in the Galactic Center. The gamma-ray flux from the innermost region of the Galactic Center is in particular suppressed. We find that for dark matter density profiles with steep inner slopes, the morphology of the Inner Galaxy gamma-ray emission in p-wave models can be made similar to the morphology in standard s-wave models. This similarity may suggest that model discrimination between s-wave and p-wave is challenging, for example, when fitting the Galactic Center excess. However, we show that it is difficult to simultaneously match s- and p-wave morphologies at both large and small angular scales. The $J$-factors we calculate may be implemented with astrophysical foreground models to self-consistently determine the morphology of the excess with velocity-suppressed dark matter annihilation.

We consider the effect of velocity-dependent dark matter annihilation on the angular distribution of gamma rays produced in dark matter subhalos. We assume that the dark matter potential is spherically symmetric, characterized by a scale radius and scale density, and the velocity distribution is isotropic. We find that the effect of velocity-dependent dark matter annihilation is largely determined by dimensional analysis; the angular size of gamma-ray emission from an individual subhalo is rescaled by a factor which depends on the form of the dark matter distribution, but not on the halo parameters, while the relative normalization of the gamma-ray flux from different mass subhalos is rescaled by a factor which depends on the halo parameters, but not on the form of the dark matter distribution. We apply our results to a Navarro-Frenk-White profile for the case of an individual subhalo and comment on the application of these results to a distribution of subhalos.

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 Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

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.

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

We discuss the theory of pulsar-timing and astrometry probes of a stochastic gravitational-wave background with a recently developed "total-angular-momentum" (TAM) formalism for cosmological perturbations. We review the formalism, emphasizing in particular the features relevant for this work and describe the observables we consider (i.e. the pulsar redshift and stellar angular displacement). Using the TAM approach, we calculate the angular power spectra for the observables and from them derive angular auto- and cross-correlation functions. We provide the full set of power spectra and correlation functions not only for the standard transverse-traceless propagating degrees of freedom in general relativity, but also for the four additional non-Einsteinian polarizations that may arise in alternative-gravity theories. We discuss how pulsar-timing and astrometry surveys can complement and serve as cross checks to one another and comment on the importance of testing the chirality of the gravitational-wave background as a tool to understand the nature of its sources. A simple rederivation of the power spectra from the plane-wave formalism is provided in an Appendix.

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.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

We obtain the first cosmological constraints on interactions between dark matter and protons within the formalism of nonrelativistic effective field theory developed for direct detection. For each interaction operator in the effective theory, parametrized by different powers of the relative velocity of the incoming particles, we use the Planck 2015 cosmic microwave background (CMB) temperature, polarization, and lensing anisotropy to set upper limits on the scattering cross section for all dark matter masses above 15 keV. We find that for interactions associated with a stronger dependence on velocity, dark matter and baryons stay thermally coupled for longer, but the interaction strengths are suppressed at the low temperatures relevant for Planck observations and are thus less constrained. At the same time, cross sections with stronger velocity dependencies are more constrained in the limit of small dark matter mass. In all cases, the effect of dark matter-proton scattering is most prominent on small scales in the CMB power spectra and in the matter power spectrum, and we thus expect substantial improvement over the current limits with data from ground-based CMB experiments and galaxy surveys.

We perform a comprehensive analysis of the most common early- and late-Universe solutions to the $H_0$, Ly-$\alpha$, and $S_8$ discrepancies. When considered on their own, massive neutrinos provide a natural solution to the $S_8$ discrepancy at the expense of increasing the $H_0$ tension. If all extensions are considered simultaneously, the best-fit solution has a neutrino mass sum of $\sim 0.4$ eV, a dark energy equation of state close to that of a cosmological constant, and no additional relativistic degrees of freedom. However, the $H_0$ tension, while weakened, remains unresolved. Motivated by this result, we perform a non-parametric reconstruction of the evolution of the dark energy fluid density (allowing for negative energy densities), together with massive neutrinos. When all datasets are included, there exists a residual $\sim1.9\sigma$ tension with $H_0$. If this residual tension remains in the future, it will indicate that it is not possible to solve the $H_0$ tension solely with a modification of the late-Universe dynamics within standard general relativity. However, we do find that it is possible to resolve the tension if either galaxy BAO or JLA supernovae data are omitted. We find that \textit{negative} dark energy densities are favored near redshift $z\sim2.35$ when including the Ly-$\alpha$ BAO measurement (at $\sim 2\sigma$). This behavior may point to a negative curvature, but it is most likely indicative of systematics or at least an underestimated covariance matrix. Quite remarkably, we find that in the extended cosmologies considered in this work, the neutrino mass sum is always close to $0.4$ eV regardless of the choice of external datasets, as long as the $H_0$ tension is solved or significantly decreased.

We present the first cosmological constraint on dark matter scattering with protons in the early Universe for the entire range of dark matter masses between 1 keV and 1 TeV. This constraint is derived from the Planck measurements of the cosmic microwave background (CMB) temperature and polarization anisotropy, and the CMB lensing anisotropy. It improves upon previous CMB constraints by many orders of magnitude, where limits are available, and closes the gap in coverage for low-mass dark matter candidates. We focus on two canonical interaction scenarios: spin-independent and spin-dependent scattering with no velocity dependence. Our results exclude (with 95% confidence) spin-independent interactions with cross sections greater than $5.3 \times 10^{-27}$ cm$^2$ for 1 keV, $3.0 \times 10^{-26}$ cm$^2$ for 1 MeV, $1.7 \times 10^{-25}$ cm$^2$ for 1 GeV, and $1.6 \times 10^{-23}$ cm$^2$ for 1 TeV dark matter mass. Finally, we discuss the implications of this study for dark matter physics and future observations.

We investigate the hypothesis that Coulomb-type interactions between dark matter (DM) and baryons explain the anomalously low 21cm brightness-temperature minimum at redshift z ~ 17 that was recently measured by the EDGES experiment. In particular, we reassess the validity of the scenario where a small fraction of the total DM is millicharged, focusing on newly derived constraints from Planck 2015 cosmic microwave background (CMB) data. Crucially, the CMB power spectrum is sensitive to DM-baryon scattering if the fraction of interacting DM is larger than (or comparable to) the fractional uncertainty in the baryon energy density. Meanwhile, there is a mass-dependent lower limit on the fraction for which the required interaction to cool the baryons sufficiently is so strong that it drives the interacting-DM temperature to the baryon temperature prior to their decoupling from the CMB. If this occurs as early as recombination, the cooling saturates. We precisely determine the viable parameter space for millicharged DM, and find that only a fraction (m_chi/MeV) 0.0115% <~ f <~ 0.4% of the entire DM content, and only for DM-particle masses between 0.5 MeV - 35 MeV, can be charged at the level needed to marginally explain the anomaly, without violating limits from SLAC, CMB, Big-Bang nucleosynthesis (BBN), or stellar and SN1987A cooling. In reality, though, we demonstrate that at least moderate fine tuning is required to both agree with the measured absorption profile and overcome various astrophysical sources of heating. Finally, we point out that a ~0.4% millicharged DM component which is tightly coupled to the baryons at recombination may resolve the current 2-sigma tension between the BBN and CMB determinations of the baryon energy density. Future CMB-S4 measurements will be able to probe this scenario directly.

We perform an improved cosmic microwave background (CMB) analysis to search for dark matter-proton scattering with a momentum-transfer cross section of the form $\sigma_0 v^n$ for $n=-2$ and $n=-4$. In particular, we present a new and robust prescription for incorporating the relative bulk velocity between the dark matter and baryon fluids into the standard linear Boltzmann calculation. Using an iterative procedure, we self-consistently include the effects of the bulk velocities in a cosmology in which dark matter interacts with baryons. With this prescription, we derive CMB bounds on the cross section, excluding $\sigma_0 > 2.3 \times 10^{-33}~\mathrm{cm}^2$ for $n=-2$ and $\sigma_0 > 1.7 \times 10^{-41}~\mathrm{cm}^2$ for $n=-4$ at $95\%$ confidence, for dark matter masses below 10 MeV. Furthermore, we investigate how these constraints change when only a subcomponent of dark matter is interacting. We show that Planck limits vanish if $\lesssim 0.4\%$ of dark matter is tightly coupled to baryons. We discuss the implications of our results for present and future cosmological observations.

We forecast constraints on dark matter (DM) scattering with baryons in the early Universe with upcoming and future cosmic microwave background (CMB) experiments, for DM particle masses down to 15 keV. In terms of the upper limit on the interaction cross section for a velocity-independent spin-independent elastic scattering, compared to current Planck results, we find a factor of $\sim$6 improvement with CMB-Stage 3, a factor of $\sim$26 with CMB-Stage 4, and a factor of $\sim$200 with a cosmic-variance limited experiment. Once the instrumental noise reaches the proximity of 1 $\mu$K-arcmin, the constraints are entirely driven by the lensing measurements. The constraints benefit from a wide survey, and show gradual improvement for instrumental noise levels from 10 $\mu$K-arcmin to 1 $\mu$K-arcmin and resolution from 5 arcmin to 1 arcmin. We further study degeneracies between DM interactions and various other signatures of new physics targeted by the CMB experiments. In the primary temperature and polarization only, we find moderate degeneracy between the effects of DM scattering, signals from massive neutrinos, and from the effective number of relativistic degrees of freedom. The degeneracy is almost entirely broken once the lensing convergence spectrum is included into the analyses. We discuss the implications of our findings in context of planned and upcoming CMB measurements and other cosmological probes of dark-sector and neutrino physics.

We present a general, model-independent formalism for determining bounds on the production of photons in dwarf spheroidal galaxies via dark matter annihilation, applicable to any set of assumptions about dark matter particle physics or astrophysics. As an illustration, we analyze gamma-ray data from the Fermi Large Area Telescope to constrain a variety of nonstandard dark matter models, several of which have not previously been studied in the context of dwarf galaxy searches.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

For models in which dark matter annihilation is Sommerfeld-enhanced, the annihilation cross section increases at low relative velocities. Dwarf spheroidal galaxies (dSphs) have low characteristic dark matter particle velocities and are thus ideal candidates to study such models. In this paper we model the dark matter phase space of dSphs as isotropic and spherically-symmetric, and determine the $J$-factors for several of the most important targets for indirect dark matter searches. For Navarro-Frenk-White density profiles, we quantify the scatter in the $J$-factor arising from the astrophysical uncertainty in the dark matter potential. We show that, in Sommerfeld-enhanced models, the ordering of the most promising dSphs may be different relative to the standard case of velocity-independent cross sections. This result can have important implications for derived upper limits on the annihilation cross section, or on possible signals, from dSphs.

Slow-roll inflation can become eternal if the quantum variance of the inflaton field around its slowly rolling classical trajectory is converted into a distribution of classical spacetimes inflating at different rates, and if the variance is large enough compared to the rate of classical rolling that the probability of an increased rate of expansion is sufficiently high. Both of these criteria depend sensitively on whether and how perturbation modes of the inflaton interact and decohere. Decoherence is inevitable as a result of gravitationally-sourced interactions whose strength are proportional to the slow-roll parameters. However, the weakness of these interactions means that decoherence is typically delayed until several Hubble times after modes grow beyond the Hubble scale. We present perturbative evidence that decoherence of long-wavelength inflaton modes indeed leads to an ensemble of classical spacetimes with differing cosmological evolutions. We introduce the notion of per-branch observables---expectation values with respect to the different decohered branches of the wave function---and show that the evolution of modes on individual branches varies from branch to branch. Thus single-field slow-roll inflation fulfills the quantum-mechanical criteria required for the validity of the standard picture of eternal inflation. For a given potential, the delayed decoherence can lead to slight quantitative adjustments to the regime in which the inflaton undergoes eternal inflation.

Many models currently exist which attempt to interpret the excess of gamma rays emanating from the Galactic Center in terms of annihilating or decaying dark matter. These models typically exhibit a variety of complicated cascade mechanisms for photon production, leading to a non-trivial kinematics which obscures the physics of the underlying dark sector. In this paper, by contrast, we observe that the spectrum of the gamma-ray excess may actually exhibit an intriguing "energy-duality" invariance under $E_\gamma \rightarrow E_\ast^2/E_\gamma$ for some $E_\ast$. As we shall discuss, such an energy duality points back to a remarkably simple alternative kinematics which in turn is realized naturally within the Dynamical Dark Matter framework. Observation of this energy duality could therefore provide considerable information about the properties of the dark sector from which the Galactic-Center gamma-ray excess might arise, and highlights the importance of acquiring more complete data for the Galactic-Center excess in the energy range around 1 GeV.

Identifying signatures of dark matter at indirect-detection experiments is generally more challenging for scenarios involving non-minimal dark sectors such as Dynamical Dark Matter (DDM) than for scenarios involving a single dark particle. This additional difficulty arises because the partitioning of the total dark-matter abundance across an ensemble of different constituent particles with different masses tends to "smear" the injection spectra of photons and other cosmic-ray particles that are produced via dark-matter annihilation or decay. As a result, the imprints of the dark sector on these cosmic-ray flux spectra typically take the form of continuum features rather than sharp peaks or lines. In this paper, however, we identify an unambiguous signature of non-minimal dark sectors such as DDM which can overcome these issues and potentially be observed at gamma-ray telescopes operating in the MeV range. We discuss the specific situations in which this signature can arise, and demonstrate that this signature can be exploited in order to significantly enhance our ability to resolve the unique spectral features of DDM and other non-minimal dark sectors at future gamma-ray facilities.

We study the atomic physics and the astrophysical implications of a model in which the dark matter is the analog of hydrogen in a secluded sector. The self interactions between dark matter particles include both elastic scatterings as well as inelastic processes due to a hyperfine transition. The self-interaction cross sections are computed by numerically solving the coupled Schr\"{o}dinger equations for this system. We show that these self interactions exhibit the right velocity dependence to explain the low dark matter density cores seen in small galaxies while being consistent with all constraints from observations of clusters of galaxies. For a viable solution, the dark hydrogen mass has to be in 10--100 GeV range and the dark fine-structure constant has to be larger than 0.02. Precisely for this range of parameters, we show that significant cooling losses may occur due to inelastic excitations to the hyperfine state and subsequent decays, with implications for the evolution of low-mass halos and the early growth of supermassive black holes. Cooling from excitations to higher $n$ levels of dark hydrogen and subsequent decays is possible at the cluster scale, with a strong dependence on halo mass. Finally, we show that the minimum halo mass is in the range of $10^{3.5}$ to $10^7 M_\odot$ for the viable regions of parameter space, significantly larger than the typical predictions for weakly-interacting dark matter models. This pattern of observables in cosmological structure formation is unique to this model, making it possible to rule in or rule out hidden sector hydrogen as a viable dark matter model.

We argue that, under certain plausible assumptions, de Sitter space settles into a quiescent vacuum in which there are no dynamical quantum fluctuations. Such fluctuations require either an evolving microstate, or time-dependent histories of out-of-equilibrium recording devices, which we argue are absent in stationary states. For a massive scalar field in a fixed de Sitter background, the cosmic no-hair theorem implies that the state of the patch approaches the vacuum, where there are no fluctuations. We argue that an analogous conclusion holds whenever a patch of de Sitter is embedded in a larger theory with an infinite-dimensional Hilbert space, including semiclassical quantum gravity with false vacua or complementarity in theories with at least one Minkowski vacuum. This reasoning provides an escape from the Boltzmann brain problem in such theories. It also implies that vacuum states do not uptunnel to higher-energy vacua and that perturbations do not decohere while slow-roll inflation occurs, suggesting that eternal inflation is much less common than often supposed. On the other hand, if a de Sitter patch is a closed system with a finite-dimensional Hilbert space, there will be Poincare recurrences and dynamical Boltzmann fluctuations into lower-entropy states. Our analysis does not alter the conventional understanding of the origin of density fluctuations from primordial inflation, since reheating naturally generates a high-entropy environment and leads to decoherence, nor does it affect the existence of non-dynamical vacuum fluctuations such as those that give rise to the Casimir effect.

The astrophysics community is considering plans for a variety of gamma-ray telescopes (including ACT and GRIPS) in the energy range 1--100 MeV, which can fill in the so-called "MeV gap" in current sensitivity. We investigate the utility of such detectors for the study of low-mass dark matter annihilation or decay. For annihilating (decaying) dark matter with a mass below about 140 MeV (280 MeV) and couplings to first generation quarks, the final states will be dominated by photons or neutral pions, producing striking signals in gamma-ray telescopes. We determine the sensitivity of future detectors to the kinematically allowed final states. In particular, we find that planned detectors can improve on current sensitivity to this class of models by up to a few orders of magnitude.

We consider the prospects for the indirect detection of low mass dark matter which couples dominantly to quarks. If the center of mass energy is below about 280 MeV, the kinematically allowed final states will be dominated by photons and neutral pions, producing striking signatures at gamma ray telescopes. In fact, an array of new instruments have been proposed, which would greatly improve sensitivity to photons in this energy range. We find that planned instruments can improve on current sensitivity to dark matter models of this type by up to a few orders of magnitude.

There is strong evidence in favor of the idea that dark matter is self interacting, with the cross section-to-mass ratio $\sigma / m \sim 1\,\mathrm{cm^2/g} \sim 1\,\mathrm{barn/GeV}$. We show that viable models of dark matter with this large cross section are straightforwardly realized with non-Abelian hidden sectors. In the simplest of such models, the hidden sector is a pure gauge theory, and the dark matter is composed of hidden glueballs with a mass around $100\,\mathrm{MeV}$. Alternatively, the hidden sector may be a supersymmetric pure gauge theory with a $\sim 10\,\mathrm{TeV}$ gluino thermal relic. In this case, the dark matter is largely composed of glueballinos that strongly self interact through the exchange of light glueballs. We present a unified framework that realizes both of these possibilities in anomaly-mediated supersymmetry breaking, where, depending on a few model parameters, the dark matter may be composed of hidden glueballinos, hidden glueballs, or a mixture of the two. These models provide simple examples of multicomponent dark matter, have interesting implications for particle physics and cosmology, and include cases where a subdominant component of dark matter may be extremely strongly self interacting, with interesting astrophysical consequences.

We consider a simple supersymmetric hidden sector: pure SU(N) gauge theory. Dark matter is made up of hidden glueballinos with mass $m_X$ and hidden glueballs with mass near the confinement scale $\Lambda$. For $m_X \sim 1\,\text{TeV}$ and $\Lambda \sim 100\,\text{MeV}$, the glueballinos freeze out with the correct relic density and self-interact through glueball exchange to resolve small-scale structure puzzles. An immediate consequence is that the glueballino spectrum has a hyperfine splitting of order $\Lambda^2 / m_X \sim 10\,\text{keV}$. We show that the radiative decays of the excited state can explain the observed 3.5 keV X-ray line signal from clusters of galaxies, Andromeda, and the Milky Way.

Many modern cosmological scenarios feature large volumes of spacetime in a de Sitter vacuum phase. Such models are said to be faced with a "Boltzmann Brain problem" - the overwhelming majority of observers with fixed local conditions are random fluctuations in the de Sitter vacuum, rather than arising via thermodynamically sensible evolution from a low-entropy past. We argue that this worry can be straightforwardly avoided in the Many-Worlds (Everett) approach to quantum mechanics, as long as the underlying Hilbert space is infinite-dimensional. In that case, de Sitter settles into a truly stationary quantum vacuum state. While there would be a nonzero probability for observing Boltzmann-Brain-like fluctuations in such a state, "observation" refers to a specific kind of dynamical process that does not occur in the vacuum (which is, after all, time-independent). Observers are necessarily out-of-equilibrium physical systems, which are absent in the vacuum. Hence, the fact that projection operators corresponding to states with observers in them do not annihilate the vacuum does not imply that such observers actually come into existence. The Boltzmann Brain problem is therefore much less generic than has been supposed.

The standard $\Lambda$CDM model provides an excellent fit to current cosmological observations but suffers from a potentially serious Boltzmann Brain problem. If the universe enters a de Sitter vacuum phase that is truly eternal, there will be a finite temperature in empty space and corresponding thermal fluctuations. Among these fluctuations will be intelligent observers, as well as configurations that reproduce any local region of the current universe to arbitrary precision. We discuss the possibility that the escape from this unacceptable situation may be found in known physics: vacuum instability induced by the Higgs field. Avoiding Boltzmann Brains in a measure-independent way requires a decay timescale of order the current age of the universe, which can be achieved if the top quark pole mass is approximately 178 GeV. Otherwise we must invoke new physics or a particular cosmological measure before we can consider $\Lambda$CDM to be an empirical success.

The decay and annihilation cross sections of dark matter particles may depend on the value of a chameleonic scalar field that both evolves cosmologically and takes different values depending on the local matter density. This possibility introduces a separation between the physics relevant for freeze-out and that responsible for dynamics and detection in the late universe. We investigate how such dark sector interactions might be implemented in a particle physics Lagrangian and consider how current and upcoming observations and experiments bound such dark matter candidates. A specific simple model allows for an increase in the annihilation cross section by a factor of $10^6$ between freeze-out and today, while more complicated models should also allow for scattering cross sections near the astrophysical bounds.