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.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

We argue that extensions of the Standard Model (SM) with a strongly first-order electroweak phase transition generically predict significant deviations of the Higgs couplings to gluons, photons, and Z bosons from their SM values. Precise experimental measurements of the Higgs couplings at the LHC and at the proposed next-generation facilities will allow for a robust test of the phase transition dynamics. To illustrate this point, in this paper we focus on the scenario in which loops of a new scalar field are responsible for the first-order phase transition, and study a selection of benchmark models with various SM gauge quantum numbers of the new scalar. We find that the current LHC measurement of the Higgs coupling to gluons already excludes the possibility of a first-order phase transition induced by a scalar in a sextet, or larger, representation of the SU(3)_c. Future LHC experiments (including HL-LHC) will be able to definitively probe the case when the new scalar is a color triplet. If the new scalar is not colored, an electron-positron Higgs factory, such as the proposed ILC or TLEP, would be required to test the nature of the phase transition. The extremely precise measurement of the Higgsstrahlung cross section possible at such machines will allow for a comprehensive and definitive probe of the possibility of a first-order electroweak phase transition in all models we considered, including the case when the new scalar is a pure gauge singlet.

In a recent paper arXiv:1107.5048, we discussed the correlation between the elastic neutralino-nucleon scattering cross section, constrained by dark matter direct detection experiments, and fine-tuning at tree-level in the electroweak symmetry breaking sector of the Minimal Supersymmetric Standard Model (MSSM). Here, we show that the correlation persists in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), and its variant, lambda-SUSY. Both models are strongly motivated by the recent discovery of a 125 GeV Higgs-like particle. We also discuss the implications of the recently published bound on the direct detection cross section from 225 live days of XENON100 experiment. In both the MSSM and the NMSSM, most of the parameter space with fine-tuning less than 10% is inconsistent with the XENON100 bound. In lambda-SUSY, on the other hand, large regions of completely natural electroweak symmetry breaking are still allowed, primarily due to a parametric suppression of fine-tuning with large \lambda. The upcoming XENON1T experiment will be able to probe most of the parameter space with less than 1% fine-tuning in all three models.

Experiments at electron-positron colliders can search for dark matter particle pair-production in association with a photon. We estimate the sensitivity of this search at the proposed International Linear Collider (ILC), under a variety of run scenarios. We employ the effective operator formalism to provide a quasi-model-independent theoretical description of the signal, and present the reach of the ILC in terms of the scale \Lambda suppressing the dark matter-electron coupling operator. We find that at the 250 GeV center-of-mass energy, the ILC can probe \Lambda up to 1-1.2 TeV, a factor of 2.5-3 above the best current bounds from LEP-2. With 1 TeV energy and polarized beams, the reach can be extended to 3-4 TeV. The ILC can discover this signature even if annihilation to electrons provides only a small fraction of the total dark matter annihilation rate in the early universe. We also argue that large regions of parameter space allowed by current LHC and direct detection bounds will be accessible at the ILC.

We study theoretical implications of direct dark matter searches in the minimal supersymmetric standard model (MSSM). We assume that no accidental cancellations occur in the spin-independent elastic neutralino-quark scattering cross section, but do not impose any relations among the weak-scale MSSM parameters. We show that direct detection cross section below 10^-44 cm^2 requires the lightest supersymmetric particle (LSP) neutralino to be close to either a pure gaugino or pure Higgsino limit, with smaller cross sections correlated with smaller admixture of the subdominant components. The current XENON100 bound rules out essentially all models in which the lightest neutralino has the Higgsino fraction between 0.2 and 0.8. Furthermore, smaller direct detection cross sections correlate with stronger fine-tuning in the electroweak symmetry breaking sector. In the gaugino LSP scenario, the current XENON100 bound already implies some fine-tuning: for example, at least 10% tuning is required if the LSP mass is above 70 GeV. In both gaugino and Higgsino LSP scenarios, the direct dark matter searches currently being conducted and designed should lead to a discovery if no accidental cancellations or fine-tuning at a level below 1% is present.

The introduction of T parity dramatically improves the consistency of Little Higgs models with precision electroweak data, and renders the lightest T-odd particle (LTP) stable. In the Littlest Higgs model with T parity, the LTP is typically the T-odd heavy photon, which is weakly interacting and can play the role of dark matter. We analyze the relic abundance of the heavy photon, including its coannihilations with other T-odd particles, and map out the regions of the parameter space where it can account for the observed dark matter. We evaluate the prospects for direct and indirect discovery of the heavy photon dark matter. The direct detection rates are quite low and a substantial improvement in experimental sensitivity would be required for observation. A substantial flux of energetic gamma rays is produced in the annihilation of the heavy photons in the galactic halo. This flux can be observed by the GLAST telescope, and, if the distribution of dark matter in the halo is favorable, by ground-based telescope arrays such as VERITAS and HESS.

Energetic antiprotons in cosmic rays can serve as an important indirect signature of dark matter. Conventionally, the antiproton flux from dark matter decays or annihilations is calculated by solving the transport equation with a space-independent diffusion coefficient within the diffusion zone of the galaxy, and assuming free propagation outside this zone. Antiproton sources outside of the diffusion zone are ignored. In reality, it is far more likely that the diffusion coefficient increases smoothly with distance from the disk, and the outlying part of the dark matter halo ignored in the conventional approach can be significant, containing as much as 90% of the galactic dark matter by mass in some models. We extend the conventional approach to address these issues. We obtain analytic approximations and numerical solutions for antiproton flux assuming that the diffusion coefficient increases exponentially with the distance from the disk, and including contributions from dark matter annihilations/decays in essentially the full dark matter halo. We find that the antiproton flux predicted in this model deviates from the conventional calculation for the same dark matter parameters by up to about 25%.

If the positron fraction and combined electron-positron flux excesses recently observed by PAMELA, Fermi and HESS have a dark matter origin, final state radiation (FSR) photons from dark matter annihilation into lepton-rich final states may be detected with observations of satellite dwarf galaxies of the Milky Way by ground-based atmospheric Cherenkov telescopes (ACTs). We find that current and near-future ACTs have excellent potential for such detection, although a discovery cannot be guaranteed due to large uncertainties in the distribution of dark matter within the dwarfs. We find that models predicting dark matter annihilation into two-lepton final states and those favoring four-lepton final states (as in, for example, "axion portal" models) can be reliably distinguished using the FSR photon spectrum once measured, and the dark matter particle mass can also be accurately determined.

If the positron fraction and combined electron-positron flux excesses recently observed by PAMELA, FERMI and HESS are due to dark matter annihilation into lepton-rich final states, the accompanying final state radiation (FSR) photons may be detected by ground-based atmospheric Cherenkov telescopes (ACTs). Satellite dwarf galaxies in the vicinity of the Milky Way are particularly promising targets for this search. We find that current and near-future ACTs have an excellent potential for discovering the FSR photons from dwarfs, although a discovery cannot be guaranteed due to large uncertainties in the fluxes resulting from lack of precise knowledge of dark matter distribution within the dwarfs. We also investigate the possibility of discriminating between different dark matter models based on the measured FSR photon spectrum. For typical parameters, we find that the ACTs can reliably distinguish models predicting dark matter annihilation into two-lepton final states from those favoring four-lepton final states (as in, for example, "axion portal" models). In addition, we find that the dark matter particle mass can also be determined from the FSR spectrum.

Energetic positrons produced in annihilation or decay of dark matter particles in the Milky Way can serve as an important indirect signature of dark matter. Computing the positron flux expected in a given dark matter model involves solving transport equations, which account for interaction of positrons with matter and galactic magnetic fields. Existing calculations solve the equations inside the diffusion zone, where galactic magnetic fields confine positrons, and assume vanishing positron density on the boundaries of this zone. However, in many models, a substantial fraction of the dark matter halo lies outside the diffusion zone. Positrons produced there can then enter the diffusion zone and get trapped, potentially reaching the Earth and increasing the expected flux. We calculate this enhancement for a variety of models. We also evaluate the expected enhancement of the flux of energetic photons produced by the inverse Compton scattering of the extra positrons on starlight and cosmic microwave background. We find maximal flux enhancements of order 20% in both cases.

If dark matter consists of weakly interacting massive particles (WIMPs), annihilation of WIMPs in the galactic center may lead to an observable enhancement of high energy gamma ray fluxes. We predict the shape and normalization of the component of the flux due to final state radiation by charged particles produced in WIMP annihilation events. The prediction is made without any assumptions about the microscopic theory responsible for WIMPs, and depends only mildly on the unknown distribution of the total WIMP annihilation cross section among the possible final states. In particular, if the WIMPs annihilate into a pair of charged fermions (leptons or quarks), the photon spectrum possesses a sharp edge feature, dropping abruptly at a photon energy equal to the WIMP mass. If such a feature is observed, it would provide strong evidence for the WIMP-related nature of the flux enhancement, as well as a measurement of the WIMP mass. We discuss the prospects for observational discovery of this feature at ground-based and space-based gamma ray telescopes.

Assuming that cosmological dark matter consists of weakly interacting massive particles, we use the recent precise measurement of cosmological parameters to predict the guaranteed rates of production of such particles in association with photons at electron-positron colliders. Our approach is based on general physical principles such as detailed balancing and soft/collinear factorization. It leads to predictions that are valid across a broad range of models containing WIMPs, including supersymmetry, universal extra dimensions, and many others. We also discuss the discovery prospects for the predicted experimental signatures.

Models with low-scale breaking of global symmetries in the neutrino sector provide an alternative to the seesaw mechanism for understanding why neutrinos are light. Such models can easily incorporate light sterile neutrinos required by the LSND experiment. Furthermore, the constraints on the sterile neutrino properties from nucleosynthesis and large scale structure can be removed due to the non-conventional cosmological evolution of neutrino masses and densities. We present explicit, fully realistic supersymmetric models, and discuss the characteristic signatures predicted in the angular distributions of the cosmic microwave background.

Recent experimental data indicates that the fine structure constant alpha may be varying on cosmological time scales. We consider the possibility that such a variation could be induced by a second order phase transition which occurs at late times (z ~ 1 - 3) and involves a change in the vacuum expectation value (vev) of a scalar with milli-eV mass. Such light scalars are natural in supersymmetric theories with low SUSY breaking scale. If the vev of this scalar contributes to masses of electrically charged fields, the low-energy value of alpha changes during the phase transition. The observational predictions of this scenario include isotope-dependent deviations from Newtonian gravity at sub-millimeter distances, and (if the phase transition is a sharp event on cosmological time scales) the presence of a well-defined step-like feature in the alpha(z) plot. The relation between the fractional changes in alpha and the QCD confinement scale is highly model dependent, and even in grand unified theories the change in alpha does not need to be accompanied by a large shift in nucleon masses.

The possibility that the energy density of the Universe is dominated by a network of low-tension domain walls provides an alternative to the commonly discussed cosmological constant and scalar-field quintessence models of dark energy. We quantify the lower bound on the number density of the domain walls that follows from the observed near-isotropy of the cosmic microwave background radiation. This bound can be satisfied by a strongly frustrated domain wall network. No fine-tuning of the parameters of the underlying field theory model is required. We briefly outline the observational consequences of this model.

We examine the particle production via preheating at the end of inflation in supersymmetric theories. The inflaton and matter scalars are now necessarily complex fields, and their relevant interactions are restricted by holomorphy. In general this leads to major changes both in the inflaton dynamics and in the efficiency of the preheating process. In addition, supersymmetric models generically contain multiple isolated vacua, raising the possibility of non-thermal production of dangerous topological defects. Because of these effects, the success of leptogenesis or WIMPZILLA production via preheating depends much more sensitively on the detailed parameters in the inflaton sector than previously thought.

The existence of bulk sterile neutrinos in theories with large extra dimensions can naturally explain small 4-dimensional Dirac masses for the active neutrinos. We study a model with 3 bulk neutrinos and derive various constraints on the size of the extra dimensions from neutrino oscillation experiments. Our analysis includes recent solar and atmospheric data from SNO and Super-Kamiokande, respectively, as well as various reactor and accelerator experiments. In addition, we comment on possible extensions of the model that could accommodate the LSND results, using natural parameters only. Our most conservative bound, obtained from the atmospheric data in the hierarchical mass scheme, constrains the largest extra dimension to have a radius R < 0.82 microns. Thus, in the context of the model studied here, future gravitational experiments are unlikely to observe the effects of extra dimensions.

Recently there has been a lot of interest in models in which gravity becomes strong at the TeV scale. The observed weakness of gravitational interactions is then explained by the existence of extra compact dimensions of space, which are accessible to gravity but not to Standard Model particles. In this letter we consider graviton emission into these extra dimensions from a hot supernova core. The phenomenology of SN1987A places strong constraints on this energy loss mechanism, allowing us to derive a bound on the fundamental Planck scale. For the case of two extra dimensions we obtain a very strong bound of M > 50 TeV, which corresponds to a radius R < 0.3 microns. While there are a lot of sources of uncertainty associated with this bound, we find that pushing it down to the few-TeV range, which could in principle be probed by ground-based experiments, is disfavored. For three or more extra dimensions the SN1987A constraints do not exclude a TeV gravitational scale.

We compute the annihilation cross-section of two neutralinos at rest into two photons, which is of importance for the indirect detection of neutralino dark matter in the galactic halo through a quasi-monochromatic gamma-ray line. We include all diagrams to one-loop level in the minimal supersymmetric extension of the standard model. We use the helicity formalism, the background-field gauge, and an efficient loop-integral reduction method. We confirm the result recently obtained by Bergstrom and Ullio in a different gauge, which disagrees with other published calculations.