Upcoming searches for the stochastic background of inflationary gravitational waves (GWs) offer the exciting possibility to probe the evolution of our Universe prior to Big Bang nucleosynthesis. In this spirit, we explore the sensitivity of future GW observations to a broad class of beyond-the-Standard-Model scenarios that lead to a nonstandard expansion history. We consider a new scalar field whose coherent oscillations dominate the energy density of the Universe at very early times, resulting in a scalar era prior to the standard radiation-dominated era. The imprint of this scalar era on the primordial GW spectrum provides a means to probe well-motivated yet elusive models of particle physics. Our work highlights the complementarity of future GW observatories across the entire range of accessible frequencies.

Scattering and decay processes of thermal bath particles involving heavy leptons can dump hot axions in the primordial plasma around the QCD phase transition. We compute their relic density, parameterized by an effective number $\Delta N_{\rm eff}$ of additional neutrinos. For couplings allowed by current bounds, production via scattering yields $\Delta N_{\rm eff} \lesssim 0.6$ and $\Delta N_{\rm eff} \lesssim 0.2$ for the cases of muon and tau, respectively. Flavor violating tau decays to a lighter lepton plus an axion give $\Delta N_{\rm eff} \lesssim 0.3$. Such values of $\Delta N_{\rm eff}$ can alleviate the tension between the direct local measurement of the Hubble constant $H_0$ and the inferred value from observations of the Cosmic Microwave Background, assuming $\Lambda$CDM. We analyze present cosmological data from the Planck collaboration and baryon acoustic oscillations with priors given in terms of the axion-lepton couplings. For axions coupled to muons, the tension can be alleviated below the 3$\sigma$ level. Future experiments will measure $\Delta N_{\rm eff}$ with higher precision, providing an axion discovery channel and probing the role of hot axions in the $H_0$ tension.

We propose a novel framework where light (sub-GeV) dark matter (DM) is detectable with future MeV gamma-ray telescopes without conflicting with Cosmic Microwave Background (CMB) data. The stable DM particle $\chi$ has a very low thermal relic abundance due to its large pair-annihilation cross section. The DM number density is stored in a heavier, meta-stable partner $\psi$ with suppressed pair-annihilation rates, that does not perturb the CMB, and whose late-time decays $\psi \rightarrow \chi$ fill the universe with $\chi$ DM particles. We provide explicit, model-independent realizations for this framework, and discuss constraints on late-time decays, and thus on parameters of this setup, from CMB, Big Bang Nucleosynthesis, and Large Scale Structure.

If the dark matter is produced in the early universe prior to Big Bang nucleosynthesis, a modified cosmological history can drastically affect the abundance of relic dark matter particles. Here, we assume that an additional species to radiation dominates at early times, causing the expansion rate at a given temperature to be larger than in the standard radiation-dominated case. We demonstrate that, if this is the case, dark matter production via freeze-in (a scenario when dark matter interacts very weakly, and is dumped in the early universe out of equilibrium by decay or scattering processes involving particles in the thermal bath) is dramatically suppressed. We illustrate and quantitatively and analytically study this phenomenon for three different paradigmatic classes of freeze-in scenarios. For the frozen-in dark matter abundance to be as large as observations, couplings between the dark matter and visible-sector particles must be enhanced by several orders of magnitude. This sheds some optimistic prospects for the otherwise dire experimental and observational outlook of detecting dark matter produced by freeze-in.

We study UV complete theories where the Standard Model (SM) gauge group is extended with a new abelian $U(1)$, and the field content is augmented by an arbitrary number of scalar and fermion SM singlets, potentially including dark matter (DM) candidates. Considerations such as classical and quantum gauge invariance of the full theory and S-matrix unitarity, not applicable within a simplified model approach, are shown to have significant phenomenological consequences. The lack of gauge anomalies leads to compact relations among the $U(1)$ fermion charges, and puts a lower bound on the number of dark fermions. Contrary to naive expectations, the DM annihilation to Zh is found to be p-wave suppressed, as hinted by perturbative unitarity of S-matrix, with dramatic implications for DM thermal relic density and indirect searches. Within this framework, the interplay between dark matter, new vector boson and Higgs physics is rather natural and generic.

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.

We study vector portal dark matter models where the mediator couples only to leptons. In spite of the lack of tree-level couplings to colored states, radiative effects generate interactions with quark fields that could give rise to a signal in current and future experiments. We identify such experimental signatures: scattering of nuclei in dark matter direct detection; resonant production of lepton-antilepton pairs at the Large Hadron Collider; and hadronic final states in dark matter indirect searches. Furthermore, radiative effects also generate an irreducible mass mixing between the vector mediator and the $Z$ boson, severely bounded by ElectroWeak Precision Tests. We use current experimental results to put bounds on this class of models, accounting for both radiatively induced and tree-level processes. Remarkably, the former often overwhelm the latter.

In all supersymmetric theories, gravitinos, with mass suppressed by the Planck scale, are an obvious candidate for dark matter; but if gravitinos ever reached thermal equilibrium, such dark matter is apparently either too abundant or too hot, and is excluded. However, in theories with an axion, a saxion condensate is generated during an early era of cosmological history and its late decay dilutes dark matter. We show that such dilution allows previously thermalized gravitinos to account for the observed dark matter over very wide ranges of gravitino mass, keV < $m_{3/2}$ < TeV, axion decay constant, $10^9$ GeV < $f_a$ < $10^{16}$ GeV, and saxion mass, 10 MeV < $m_s$ < 100 TeV. Constraints on this parameter space are studied from BBN, supersymmetry breaking, gravitino and axino production from freeze-in and saxion decay, and from axion production from both misalignment and parametric resonance mechanisms. Large allowed regions of $(m_{3/2}, f_a, m_s)$ remain, but differ for DFSZ and KSVZ theories. Superpartner production at colliders may lead to events with displaced vertices and kinks, and may contain saxions decaying to $(WW,ZZ,hh), gg, \gamma \gamma$ or a pair of Standard Model fermions. Freeze-in may lead to a sub-dominant warm component of gravitino dark matter, and saxion decay to axions may lead to dark radiation.

We consider a modification to the standard cosmological history consisting of introducing a new species $\phi$ whose energy density red-shifts with the scale factor $a$ like $\rho_\phi \propto a^{-(4+n)}$. For $n>0$, such a red-shift is faster than radiation, hence the new species dominates the energy budget of the universe at early times while it is completely negligible at late times. If equality with the radiation energy density is achieved at low enough temperatures, dark matter can be produced as a thermal relic during the new cosmological phase. Dark matter freeze-out then occurs at higher temperatures compared to the standard case, implying that reproducing the observed abundance requires significantly larger annihilation rates. Here, we point out a completely new phenomenon, which we refer to as $\textit{relentless}$ dark matter: for large enough $n$, unlike the standard case where annihilation ends shortly after the departure from thermal equilibrium, dark matter particles keep annihilating long after leaving chemical equilibrium, with a significant depletion of the final relic abundance. Relentless annihilation occurs for $n \geq 2$ and $n \geq 4$ for s-wave and p-wave annihilation, respectively, and it thus occurs in well motivated scenarios such as a quintessence with a kination phase. We discuss a few microscopic realizations for the new cosmological component and highlight the phenomenological consequences of our calculations for dark matter searches.

We study direct detection in simplified models of Dark Matter (DM) in which interactions with Standard Model (SM) fermions are mediated by a heavy vector boson. We consider fully general, gauge-invariant couplings between the SM, the mediator and both scalar and fermion DM. We account for the evolution of the couplings between the energy scale of the mediator mass and the nuclear energy scale. This running arises from virtual effects of SM particles and its inclusion is not optional. We compare bounds on the mediator mass from direct detection experiments with and without accounting for the running. In some cases the inclusion of these effects changes the bounds by several orders of magnitude, as a consequence of operator mixing which generates new interactions at low energy. We also highlight the importance of these effects when translating LHC limits on the mediator mass into bounds on the direct detection cross section. For an axial-vector mediator, the running can alter the derived bounds on the spin-dependent DM-nucleon cross section by a factor of two or more. Finally, we provide tools to facilitate the inclusion of these effects in future studies: general approximate expressions for the low energy couplings and a public code runDM to evolve the couplings between arbitrary energy scales.

We explore the phenomenology of a class of models where the dark matter particle can inelastically up-scatter to a heavier excited state via off-diagonal dipolar interactions with the interstellar plasma (gas or free electrons). The heavier particle then rapidly decays back to the dark matter particle plus a quasi-monochromatic photon. For the process to occur at appreciable rates, the mass splitting between the heavier state and the dark matter must be comparable to, or smaller than, the kinetic energy of particles in the plasma. As a result, the predicted photon line falls in the soft X-ray range, or, potentially, at arbitrarily lower energies. We explore experimental constraints from cosmology and particle physics, and present accurate calculations of the dark matter thermal relic density and of the flux of monochromatic X-rays from thermal plasma excitation. We find that the model provides a natural explanation for the observed 3.5 keV line from clusters of galaxies and from the Galactic center, and is consistent with null detections of the line from dwarf galaxies. The unique line shape, which will be resolved by future observations with the Hitomi (formerly Astro-H) satellite, and the predicted unique morphology and target-temperature dependence will enable easy discrimination of this class of models versus other scenarios for the generation of the 3.5 keV line or of any other unidentified line across the electromagnetic spectrum.

We study the effective field theory obtained by extending the Standard Model field content with two singlets: a 750 GeV (pseudo-)scalar and a stable fermion. Accounting for collider productions initiated by both gluon and photon fusion, we investigate where the theory is consistent with both the LHC diphoton excess and bounds from Run 1. We analyze dark matter phenomenology in such regions, including relic density constraints as well as collider, direct, and indirect bounds. Scalar portal dark matter models are very close to limits from direct detection and mono-jet searches if gluon fusion dominates, and not constrained at all otherwise. Pseudo-scalar models are challenged by photon line limits and mono-jet searches in most of the parameter space.

Direct searches for WIMPs are sensitive to physics well below the weak scale. In the absence of light mediators, it is fruitful to apply an Effective Field Theory (EFT) approach accounting only for dark matter (DM) interactions with Standard Model (SM) fields. We consider a singlet fermion WIMP and effective operators up to dimension 6 which are generated at the mass scale of particles mediating DM interactions with the SM. We perform a one-loop Renormalization Group Evolution (RGE) analysis, evolving these effective operators from the mediators mass scale to the nuclear scales probed by direct searches. We apply our results to models with DM velocity-suppressed interactions, DM couplings only to heavy quarks, leptophilic DM and Higgs portal, which without our analysis would not get constrained from direct detection bounds. Remarkably, a large parameter space region for these models is found to be excluded as a consequence of spin-independent couplings induced by SM loops. In addition to these examples, we stress that more general renormalizable models for singlet fermion WIMP can be matched onto our EFT framework, and the subsequent model-independent RGE can be used to compute direct detection rates. Our results allow us to properly connect the different energy scales involved in constraining WIMP models, and to combine information from direct detection with other complementary searches, such as collider and indirect detection.

We consider an effective field theory for a gauge singlet Dirac dark matter (DM) particle interacting with the Standard Model (SM) fields via effective operators suppressed by the scale $\Lambda \gtrsim 1$ TeV. We perform a systematic analysis of the leading loop contributions to spin-independent (SI) DM--nucleon scattering using renormalization group evolution between $\Lambda$ and the low-energy scale probed by direct detection experiments. We find that electroweak interactions induce operator mixings such that operators that are naively velocity-suppressed and spin-dependent can actually contribute to SI scattering. This allows us to put novel constraints on Wilson coefficients that were so far poorly bounded by direct detection. Constraints from current searches are comparable to LHC bounds, and will significantly improve in the near future. Interestingly, the loop contribution we find is maximally isospin violating even if the underlying theory is isospin conserving.

Hints in the Fermi data for a 130 GeV gamma line from the galactic center have ignited interest in potential gamma line signatures of dark matter. Explanations of this line based on dark matter annihilation face a parametric tension since they often rely on large enhancements of loop-suppressed cross sections. In this paper, we pursue an alternative possibility that dark matter gamma lines could arise from "semi-annihilation" among multiple dark sector states. The semi-annihilation reaction with a single final state photon is typically enhanced relative to ordinary annihilation into photon pairs. Semi-annihilation allows for a wide range of dark matter masses compared to the fixed mass value required by annihilation, opening the possibility to explain potential dark matter signatures at higher energies. The most striking prediction of semi-annihilation is the presence of multiple gamma lines, with as many as order N^3 lines possible for N dark sector states, allowing for dark sector spectroscopy. A smoking gun signature arises in the simplest case of degenerate dark matter, where a strong semi-annihilation line at 130 GeV would be accompanied by a weaker annihilation line at 173 GeV. As a proof of principle, we construct two explicit models of dark matter semi-annihilation, one based on non-Abelian vector dark matter and the other based on retrofitting Rayleigh dark matter.

Pure singlets are typically disfavored as dark matter candidates, since they generically have a thermal relic abundance larger than the observed value. In this paper, we propose a new dark matter mechanism called "assimilation", which takes advantage of the baryon asymmetry of the universe to generate the correct relic abundance of singlet dark matter. Through assimilation, dark matter itself is efficiently destroyed, but dark matter number is stored in new quasi-stable heavy states which carry the baryon asymmetry. The subsequent annihilation and late-time decay of these heavy states yields (symmetric) dark matter as well as (asymmetric) standard model baryons. We study in detail the case of pure bino dark matter by augmenting the minimal supersymmetric standard model with vector-like chiral multiplets. In the parameter range where this mechanism is effective, the LHC can discover long-lived charged particles which were responsible for assimilating dark matter.

The semi-annihilation reaction takes the schematic form psi_i psi_j -> psi_k phi, where psi_i are stable dark matter particles and phi is an unstable state. Such reactions are allowed when dark matter is stabilized by a larger symmetry than just Z_2. They lead to non-trivial dark matter dynamics in the early universe, and the thermal production of the relic particles can be completely controlled by semi-annihilations. This process might also take place today in the Milky Way, enriching the (semi-)annihilation final state spectrum observed in indirect detection experiments.

We show that the thermal relic abundance of dark matter can be affected by a new type of reaction: semi-annihilation. Semi-annihilation takes the schematic form X_i X_j -> X_k phi, where X_i are stable dark matter particles and phi is an unstable state. Such reactions are generically present when dark matter is composed of more than one species with "flavor" and/or "baryon" symmetries. We give a complete set of coupled Boltzmann equations in the presence of semi-annihilations, and study two toy models featuring this process. Semi-annihilation leads to non-trivial dark matter dynamics in the early universe, often dominating over ordinary annihilation in determining the relic abundance. This process also has important implications for indirect detection experiments, by enriching the final state spectrum from dark matter (semi-)annihilation in the Milky Way.