Collective excitations in condensed matter systems, such as phonons and magnons, have recently been proposed as novel detection channels for light dark matter. We show that excitation of i) optical phonon polaritons in polar materials in an ${\mathcal O}$(1 T) magnetic field (via the axion-photon coupling), and ii) gapped magnons in magnetically ordered materials (via the axion wind coupling to the electron spin), can cover the difficult-to-reach ${\mathcal O}$(1-100) meV mass window of QCD axion dark matter with less than a kilogram-year exposure. Finding materials with a large number of optical phonon or magnon modes that can couple to the axion field is crucial, suggesting a program to search for a range of materials with different resonant energies and excitation selection rules; we outline the rules and discuss a few candidate targets, leaving a more exhaustive search for future work. Ongoing development of single photon, phonon and magnon detectors will provide the key for experimentally realizing the ideas presented here.

Scattering of light dark matter with sub-eV energy deposition can be detected with collective excitations in condensed matter systems. When dark matter has spin-independent couplings to atoms or ions, it has been shown to efficiently excite phonons. Here we show that, if dark matter couples to the electron spin, magnon excitations in materials with magnetic dipole order offer a promising detection path. We derive general formulae for single magnon excitation rates from dark matter scattering, and demonstrate as a proof of principle the projected reach of a yttrium iron garnet target for several dark matter models with spin-dependent interactions. This highlights the complementarity of various collective excitations in probing different dark matter interactions.

Dark matter substructure on small scales is currently weakly constrained, and its study may shed light on the nature of the dark matter. In this work we study the gravitational effects of dark matter substructure on measured pulsar phases in pulsar timing arrays (PTAs). Due to the stability of pulse phases observed over several years, dark matter substructure around the Earth-pulsar system can imprint discernible signatures in gravitational Doppler and Shapiro delays. We compute pulsar phase correlations induced by general dark matter substructure, and project constraints for a few models such as monochromatic primordial black holes (PBHs), and Cold Dark Matter (CDM)-like NFW subhalos. This work extends our previous analysis, which focused on static or single transiting events, to a stochastic analysis of multiple transiting events. We find that stochastic correlations, in a PTA similar to the Square Kilometer Array (SKA), are uniquely powerful to constrain subhalos as light as $\sim 10^{-13}~M_\odot$, with concentrations as low as that predicted by standard CDM.

The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.

Pulsars act as accurate clocks, sensitive to gravitational redshift and acceleration induced by transiting clumps of matter. We study the sensitivity of pulsar timing arrays (PTAs) to single transiting compact objects, focusing on primordial black holes and compact subhalos in the mass range from $10^{-12} M _{\odot}$ to well above $100~M_\odot$. We find that the Square Kilometer Array can constrain such objects to be a subdominant component of the dark matter over this entire mass range, with sensitivity to a dark matter sub-component reaching the sub-percent level over significant parts of this range. We also find that PTAs offer an opportunity to probe substantially less dense objects than lensing because of the large effective radius over which such objects can be observed, and we quantify the subhalo concentration parameters which can be constrained.

We show that polar materials are excellent targets for direct detection of sub-GeV dark matter due to the presence of gapped optical phonons as well as acoustic phonons with high sound speed. We take the example of Gallium Arsenide (GaAs), which has the properties needed for experimental realization, and where many results can be estimated analytically. We find GaAs has excellent reach to dark photon absorption, can completely cover the freeze-in benchmark for scattering via an ultralight dark photon, and is competitive with other proposals to detect sub-MeV dark matter scattering off nuclei.

We present a detailed description of the electromagnetic filter for the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Starting with an initial estimate for the orbital magnetic moment, the higher-order drift process of ExB is configured to balance the gradient-B drift motion of the electron in such a way as to guide the trajectory into the standing voltage potential along the mid-plane of the filter. As a function of drift distance along the length of the filter, the filter zooms in with exponentially increasing precision on the transverse velocity component of the electron kinetic energy. This yields a linear dimension for the total filter length that is exceptionally compact compared to previous techniques for electromagnetic filtering. The parallel velocity component of the electron kinetic energy oscillates in an electrostatic harmonic trap as the electron drifts along the length of the filter. An analysis of the phase-space volume conservation validates the expected behavior of the filter from the adiabatic invariance of the orbital magnetic moment and energy conservation following Liouville's theorem for Hamiltonian systems.

We consider gravitationally bound states of asymmetric dark matter (ADM stars), and the impact of ADM capture on the stability of neutron stars. We derive and interpret the equation of state for ADM with both attractive and repulsive interactions, and solve the Tolman-Oppenheimer-Volkoff equations to find equilibrium sequences and maximum masses of ADM stars. Gravitational wave searches can utilize our solutions to model exotic compact objects (ECOs). Our results for attractive interactions differ substantially from those in the literature, where fermionic ADM with attractive self-interactions was employed to destabilize neutron stars more effectively than non-interacting fermionic ADM. By contrast, we argue that fermionic ADM with an attractive force is no more effective in destabilizing neutron stars than fermionic ADM with no self-interactions.

We propose to achieve the proof-of-principle of the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Each of the technological challenges described in [1,2] will be targeted and hopefully solved by the use of the latest experimental developments and profiting from the low background environment provided by the LNGS underground site. The first phase will focus on the graphene technology for a tritium target and the demonstration of TES microcalorimetry with an energy resolution of better than 0.05 eV for low energy electrons. These technologies will be evaluated using the PTOLEMY prototype, proposed for underground installation, using precision HV controls to step down the kinematic energy of endpoint electrons to match the calorimeter dynamic range and rate capabilities. The second phase will produce a novel implementation of the EM filter that is scalable to the full target size and which demonstrates intrinsic triggering capability for selecting endpoint electrons. Concurrent with the CNB program, we plan to exploit and develop the unique properties of graphene to implement an intermediate program for direct directional detection of MeV dark matter [3,4]. This program will evaluate the radio-purity and scalability of the graphene fabrication process with the goal of using recently identified ultra-high radio-purity CO2 sources. The direct detection of the CNB is a snapshot of early universe dynamics recorded by the thermal relic neutrino yield taken at a time that predates the epochs of Big Bang Nucleosynthesis, the Cosmic Microwave Background and the recession of galaxies (Hubble Expansion). Big Bang neutrinos are believed to have a central role in the evolution of the Universe and a direct measurement with PTOLEMY will unequivocally establish the extent to which these predictions match present-day neutrino densities.

Nuggets---very large stable bound objects arising in the presence of a sufficiently attractive and long-range force and in the absence of a dark Coulomb force---are a smoking gun signature for Asymmetric Dark Matter (ADM). The cosmology of ADM nuggets is both generic and unique: nuggets feature highly exothermic fusion processes, which can impact the shape of the core in galaxies, as well as give rise to rare dark star formation. We find, considering the properties of nuggets in a generic extended nuclear model with both attractive and repulsive forces, that self-interaction constraints place an upper bound on nugget masses at the freeze-out of synthesis in the ballpark of $M_{\rm fo} \lesssim 10^{16}$ GeV. We also show that indirect detection strongly constrains models where the scalar mediator binding the nuggets mixes with the Higgs.

We compute the mass function of bound states of Asymmetric Dark Matter--nuggets--synthesized in the early Universe. We apply our results for the nugget density and binding energy computed from a nuclear model to obtain analytic estimates of the typical nugget size exiting synthesis. We numerically solve the Boltzmann equation for synthesis including two-to-two fusion reactions, estimating the impact of bottlenecks on the mass function exiting synthesis. These results provide the basis for studying the late Universe cosmology of nuggets in a future companion paper.

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 propose two-dimensional materials as targets for direct detection of dark matter. Using graphene as an example, we focus on the case where dark matter scattering deposits sufficient energy on a valence-band electron to eject it from the target. We show that the sensitivity of graphene to dark matter of MeV to GeV mass can be comparable, for similar exposure and background levels, to that of semiconductor targets such as silicon and germanium. Moreover, a two-dimensional target is an excellent directional detector, as the ejected electron retains information about the angular dependence of the incident dark matter particle. This proposal can be implemented by the PTOLEMY experiment, presenting for the first time an opportunity for directional detection of sub-GeV dark matter.

Semiconductors are by now well-established targets for direct detection of MeV to GeV dark matter via scattering off electrons. We show that semiconductor targets can also detect significantly lighter dark matter via an absorption process. When the dark matter mass is above the band gap of the semiconductor (around an eV), absorption proceeds by excitation of an electron into the conduction band. Below the band gap, multi-phonon excitations enable absorption of dark matter in the 0.01 eV to eV mass range. Energetic dark matter particles emitted from the sun can also be probed for masses below an eV. We derive the reach for absorption of a relic kinetically mixed dark photon or pseudoscalar in germanium and silicon, and show that existing direct detection results already probe new parameter space. With only a moderate exposure, low-threshold semiconductor target experiments can exceed current astrophysical and terrestrial constraints on sub-keV bosonic dark matter.

A remarkable prediction of the Standard Model is that, in the absence of corrections lifting the energy density, the Higgs potential becomes negative at large field values. If the Higgs field samples this part of the potential during inflation, the negative energy density may locally destabilize the spacetime. We use numerical simulations of the Einstein equations to study the evolution of inflation-induced Higgs fluctuations as they grow towards the true (negative-energy) minimum. These simulations show that forming a single patch of true vacuum in our past light cone during inflation is incompatible with the existence of our Universe; the boundary of the true vacuum region grows outward in a causally disconnected manner from the crunching interior, which forms a black hole. We also find that these black hole horizons may be arbitrarily elongated---even forming black strings---in violation of the hoop conjecture. By extending the numerical solution of the Fokker-Planck equation to the exponentially suppressed tails of the field distribution at large field values, we derive a rigorous correlation between a future measurement of the tensor-to-scalar ratio and the scale at which the Higgs potential must receive stabilizing corrections in order for the Universe to have survived inflation until today.

We examine in depth a recent proposal to utilize superfluid helium for direct detection of sub-MeV mass dark matter. For sub-keV recoil energies, nuclear scattering events in liquid helium primarily deposit energy into long-lived phonon and roton quasiparticle excitations. If the energy thresholds of the detector can be reduced to the meV scale, then dark matter as light as ~MeV can be reached with ordinary nuclear recoils. If, on the other hand, two or more quasiparticle excitations are directly produced in the dark matter interaction, the kinematics of the scattering allows sensitivity to dark matter as light as ~keV at the same energy resolution. We present in detail the theoretical framework for describing excitations in superfluid helium, using it to calculate the rate for the leading dark matter scattering interaction, where an off-shell phonon splits into two or more higher-momentum excitations. We validate our analytic results against the measured and simulated dynamic response of superfluid helium. Finally, we apply this formalism to the case of a kinetically mixed hidden photon in the superfluid, both with and without an external electric field to catalyze the processes.

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 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 apply effective field theory methods to compute bino-nucleon scattering, in the case where tree-level interactions are suppressed and the leading contribution is at loop order via heavy flavor squarks or sleptons. We find that leading log corrections to fixed-order calculations can increase the bino mass reach of direct detection experiments by a factor of two in some models. These effects are particularly large for the bino-sbottom coannihilation region, where bino dark matter as heavy as 5-10 TeV may be detected by near future experiments. For the case of stop- and selectron-loop mediated scattering, an experiment reaching the neutrino background will probe thermal binos as heavy as 500 and 300 GeV, respectively. We present three key examples that illustrate in detail the framework for determining weak scale coefficients, and for mapping onto a low energy theory at hadronic scales, through a sequence of effective theories and renormalization group evolution. For the case of a squark degenerate with the bino, we extend the framework to include a squark degree of freedom at low energies using heavy particle effective theory, thus accounting for large logarithms through a "heavy-light current." Benchmark predictions for scattering cross sections are evaluated, including complete leading order matching onto quark and gluon operators, and a systematic treatment of perturbative and hadronic uncertainties.

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.

Superconducting targets have recently been proposed for the direct detection of dark matter as light as a keV, via elastic scattering off conduction electrons in Cooper pairs. Detecting such light dark matter requires sensitivity to energies as small as the superconducting gap of O(meV). Here we show that these same superconducting devices can detect much lighter DM, of meV to eV mass, via dark matter absorption on a conduction electron, followed by emission of an athermal phonon. We demonstrate the power of this setup for relic kinetically mixed hidden photons, pseudoscalars, and scalars, showing the reach can exceed current astrophysical and terrestrial constraints with only a moderate exposure.

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.

Identifying the true theory of dark matter depends crucially on accurately characterizing interactions of dark matter (DM) with other species. In the context of DM direct detection, we present a study of the prospects for correctly identifying the low-energy effective DM-nucleus scattering operators connected to UV-complete models of DM-quark interactions. We take a census of plausible UV-complete interaction models with different low-energy leading-order DM-nuclear responses. For each model (corresponding to different spin-, momentum-, and velocity-dependent responses), we create a large number of realizations of recoil-energy spectra, and use Bayesian methods to investigate the probability that experiments will be able to select the correct scattering model within a broad set of competing scattering hypotheses. We conclude that agnostic analysis of a strong signal (such as Generation-2 would see if cross sections are just below the current limits) seen on xenon and germanium experiments is likely to correctly identify momentum dependence of the dominant response, ruling out models with either "heavy" or "light" mediators, and enabling downselection of allowed models. However, a unique determination of the correct UV completion will critically depend on the availability of measurements from a wider variety of nuclear targets, including iodine or fluorine. We investigate how model-selection prospects depend on the energy window available for the analysis. In addition, we discuss accuracy of the DM particle mass determination under a wide variety of scattering models, and investigate impact of the specific types of particle-physics uncertainties on prospects for model selection.

We propose and study a new class of superconducting detectors which are sensitive to O(meV) electron recoils from dark matter-electron scattering. Such devices could detect dark matter as light as the warm dark matter limit, mX > keV. We compute the rate of dark matter scattering off of free electrons in a (superconducting) metal, including the relevant Pauli blocking factors. We demonstrate that classes of dark matter consistent with terrestrial and cosmological/astrophysical constraints could be detected by such detectors with a moderate size exposure.

Measurements of the Higgs boson and top quark masses indicate that the Standard Model Higgs potential becomes unstable around $\Lambda_I \sim 10^{11}$ GeV. This instability is cosmologically relevant since quantum fluctuations during inflation can easily destabilize the electroweak vacuum if the Hubble parameter during inflation is larger than $\Lambda_I$ (as preferred by the recent BICEP2 measurement). We perform a careful study of the evolution of the Higgs field during inflation, obtaining different results from those currently in the literature. We consider both tunneling via a Coleman-de Luccia or Hawking-Moss instanton, valid when the scale of inflation is below the instability scale, as well as a statistical treatment via the Fokker-Planck equation appropriate in the opposite regime. We show that a better understanding of the post-inflation evolution of the unstable AdS vacuum regions is crucial for determining the eventual fate of the universe. If these AdS regions devour all of space, a universe like ours is indeed extremely unlikely without new physics to stabilize the Higgs potential; however, if these regions crunch, our universe survives, but inflation must last a few e-folds longer to compensate for the lost AdS regions. Lastly, we examine the effects of generic Planck-suppressed corrections to the Higgs potential, which can be sufficient to stabilize the electroweak vacuum during inflation.

We explore models for the GeV Galactic Center Excess (GCE) observed by the Fermi Telescope, focusing on $\chi \chi \rightarrow f \bar f$ annihilation processes in the $Z_3$ NMSSM. We begin by examining the requirements for a simplified model (parametrized by the couplings and masses of dark matter (DM) and mediator particles) to reproduce the GCE via $\chi \chi \rightarrow f \bar f$, while simultaneously thermally producing the observed relic abundance. We apply the results of our simplified model to the $Z_3$ NMSSM for Singlino/Higgsino(S/H) or Bino/Higgsino(B/H) DM. In the case of S/H DM, we find that the DM must be be very close to a pseudoscalar resonance to be viable, and large $\tan\beta$ and positive values of $\mu$ are preferred for evading direct detection constraints while simultaneously obtaining the observed Higgs mass. In the case of B/H DM, by contrast, the situation is much less tuned: annihilation generally occurs off-resonance, and for large $\tan\beta$ direct detection constraints are easily satisfied by choosing $\mu$ sufficiently large and negative. The B/H model generally has a light, largely MSSM-like pseudoscalar with no accompanying charged Higgs, which could be searched for at the LHC.

We study the decay signatures of Asymmetric Dark Matter (ADM) via higher dimension operators which are responsible for generating the primordial dark matter (DM) asymmetry. Since the signatures are sensitive both to the nature of the higher dimension operator generating the DM asymmetry and to the sign of the baryon or lepton number that the DM carries, indirect detection may provide a window into the nature of the mechanism which generates the DM asymmetry. We consider in particular dimension-6 fermionic operators of the form ${\cal O}_{ADM} = X {\cal O}_{B-L}/M^2$, where ${\cal O}_{B-L} = u^c d^c d^c,~\ell \ell e^c,~q \ell d^c$ (or operators related through a Hermitian conjugate) with the scale $M$ around or just below the GUT scale. We derive constraints on ADM particles both in the natural mass range (around a few GeV), as well as in the range between 100 GeV to 10 TeV. For light ADM, we focus on constraints from both the low energy gamma ray data and proton/anti-proton fluxes. For heavy ADM, we consider $\gamma$-rays and proton/anti-proton fluxes, and we fit $e^+/e^-$ data from AMS-02 and H.E.S.S. (neglecting the Fermi charged particle fluxes which disagree with AMS-02 below 100 GeV). We show that, although the best fit regions from electron/positron measurement are still in tension with other channels on account of the H.E.S.S. measurement at high energies, compared to an ordinary symmetric dark matter scenario, the decay of DM with a primordial asymmetry reduces the tension. Better measurement of the flux at high energy will be necessary to draw a definite conclusion about the viability of decaying DM as source for the signals.

We examine the effect of nuclear response functions, as laid out in [Fitzpatrick et al, arXiv:1203.3542], on dark matter (DM) direct detection in the context of well-motivated UV completions, including electric and magnetic dipoles, anapole, spin-orbit, and pseudoscalar-mediated DM. Together, these encompass five of the six nuclear responses extracted from the non-relativistic effective theory of [Fitzpatrick et al, arXiv:1203.3542] (with the sixth difficult to UV complete), with two of the six combinations corresponding to standard spin-independent and -dependent responses. For constraints from existing direct detection experiments, we find that only the COUPP constraint, due to its heavy iodine target with large angular momentum and an unpaired spin, and its large energy range sensitivity, is substantially modified by the new responses compared to what would be inferred using the standard form factors to model the energy dependence of the response. For heavy targets such as xenon and germanium, the behavior of the new nuclear responses as recoil energy increases can be substantially different than that of the standard responses, but this has almost no impact on the constraints derived from experiments such as LUX, XENON100 and CDMS since the maximum nuclear recoil energy detected in these experiments is relatively low. We simulate mock data for 80 and 250 GeV DM candidates utilizing the new nuclear responses to highlight how they might affect a putative signal, and find the new responses are most important for momentum-suppressed interactions such as the magnetic dipole or pseudoscalar-mediated interaction when the target is relatively heavy (such as xenon and iodine).

We present constraints on decaying and annihilating dark matter (DM) in the 4 keV to 10 GeV mass range, using published results from the satellites HEAO-1, INTEGRAL, COMPTEL, EGRET, and the Fermi Gamma-ray Space Telescope. We derive analytic expressions for the gamma-ray spectra from various DM decay modes, and find lifetime constraints in the range 10^24-10^28 sec, depending on the DM mass and decay mode. We map these constraints onto the parameter space for a variety of models, including a hidden photino that is part of a kinetically mixed hidden sector, a gravitino with R-parity violating decays, a sterile neutrino, DM with a dipole moment, and a dark pion. The indirect constraints on sterile-neutrino and hidden-photino DM are found to be more powerful than other experimental or astrophysical probes in some parts of parameter space. While our focus is on decaying DM, we also present constraints on DM annihilation to electron-positron pairs. We find that if the annihilation is p-wave suppressed, the galactic diffuse constraints are, depending on the DM mass and velocity at recombination, more powerful than the constraints from the Cosmic Microwave Background.

We examine the consistency of light dark matter (DM) elastic scattering in CoGeNT, DAMA, and CDMS-Silicon in light of constraints from XENON, CDMS, LUX, PICASSO and COUPP. We consider a variety of operators that have been employed to reconcile anomalies with constraints, including anapole, magnetic dipole, momentum-dependent, and isospin-violating DM. We find that elastic scattering through these alternative operators does not substantially reduce the tension between the signals and the null constraints for operators where at least two of the three purported signals map onto a common space in the DM mass--scattering cross-section plane. Taking a choice of the scintillation efficiency that lies at the $-1 \sigma$ region of the Manzur et al measurement relieves tension between signals and the LUX constraint---in particular for a magnetic dipole interaction and a xenophobic interaction (though for the latter the signal regions do not substantially overlap). We also find that modest changes in the halo model does not alter this result. We conclude that, even relaxing the assumption about the type of elastic scattering interaction and taking a conservative choice for the scintillation efficiency, LUX and the results from other null experiments remain in tension with a light DM elastic scattering explanation of direct detection anomalies.

We review theories of Asymmetric Dark Matter (ADM), their cosmological implications and detection. While there are many models of ADM in the literature, our review of existing models will center on highlighting the few common features and important mechanisms for generation and transfer of the matter-anti-matter asymmetry between dark and visible sectors. We also survey ADM hidden sectors, the calculation of the relic abundance for ADM, and how the DM asymmetry may be erased at late times through oscillations. We consider cosmological constraints on ADM from the cosmic microwave background, neutron stars, the Sun, and brown and white dwarves. Lastly, we review indirect and direct detection methods for ADM, collider signatures, and constraints.

As part of the Snowmass process, the Cosmic Frontier WIMP Direct Detection subgroup (CF1) has drawn on input from the Cosmic Frontier and the broader Particle Physics community to produce this document. The charge to CF1 was (a) to summarize the current status and projected sensitivity of WIMP direct detection experiments worldwide, (b) motivate WIMP dark matter searches over a broad parameter space by examining a spectrum of WIMP models, (c) establish a community consensus on the type of experimental program required to explore that parameter space, and (d) identify the common infrastructure required to practically meet those goals.

Dark matter (DM) self-interactions have important implications for the formation and evolution of structure, from dwarf galaxies to clusters of galaxies. We study the dynamics of self-interacting DM via a light mediator, focusing on the quantum resonant regime where the scattering cross section has a non-trivial velocity dependence. While there are long-standing indications that observations of small scale structure in the Universe are not in accord with the predictions of collisionless DM, theoretical study and simulations of DM self-interactions have focused on parameter regimes with simple analytic solutions for the scattering cross section, with constant or classical velocity (and no angular) dependence. We devise a method that allows us to explore the velocity and angular dependence of self-scattering more broadly, in the strongly-coupled resonant and classical regimes where many partial modes are necessary for the achieving the result. We map out the entire parameter space of DM self-interactions --- and implications for structure observations --- as a function of the coupling and the DM and mediator masses. We derive a new analytic formula for describing resonant s-wave scattering. Finally, we show that DM self-interactions can be correlated with observations of Sommerfeld enhancements in DM annihilation through indirect detection experiments.

We examine cosmological, astrophysical and collider constraints on thermal dark matter (DM) with mass mX in the range 1 MeV to 10 GeV. Cosmic microwave background (CMB) observations, which severely constrain light symmetric DM, can be evaded if the DM relic density is sufficiently asymmetric. CMB constraints require the present anti-DM to DM ratio to be less than 2*10^{-6} (10^{-1}) for DM mass mX = 1 MeV (10 GeV) with ionizing efficiency factor f ~ 1. We determine the minimum annihilation cross section for achieving these asymmetries subject to the relic density constraint; these cross sections are larger than the usual thermal annihilation cross section. On account of collider constraints, such annihilation cross sections can only be obtained by invoking light mediators. These light mediators can give rise to significant DM self-interactions, and we derive a lower bound on the mediator mass from elliptical DM halo shape constraints. We find that halo shapes require a mediator with mass mphi > 4 * 10^{-2} MeV (40 MeV) for mX = 1 MeV (10 GeV). We map all of these constraints to the parameter space of DM-electron and DM-nucleon scattering cross sections for direct detection. For DM-electron scattering, a significant fraction of the parameter space is already ruled out by beam-dump and supernova cooling constraints.

A dark force can impact the cosmological history of dark matter (DM), both explaining observed cores in dwarf galaxies and setting the DM relic density through annihilation to dark force bosons. For GeV - TeV DM mass, DM self-scattering in dwarf galaxy halos exhibits quantum mechanical resonances, analogous to a Sommerfeld enhancement for annihilation. We show that a simple model of DM with a dark force can accommodate all astrophysical bounds on self-interactions in halos and explain the observed relic density, through a single coupling constant.

Evidence continues to grow in the MiniBooNE (MB) antineutrino mode supporting a low-energy excess compatible with the MB neutrino mode and possibly also confirming the results of the LSND experiment. At least one sterile neutrino is required to explain the anomalies consistent with the observations of other experiments. At the same time, there is a strong tension between the positive signals of LSND and MB and the null results of nu_e and nu_mu disappearance experiments. We explore a scenario, first proposed in \cite{Nelson:2010hz}, where the presence of an additional heavy sterile neutrino (with mass well above an eV) can alleviate tension between LSND, MB and the null results of disappearance experiments. We compare and contrast this 3+1+1 scenario with the more standard 3+1 scenario and carry out global fits to all oscillation data including new 2011 MB anti-nu data. We find that the tension can be somewhat alleviated and that a phenomenologically viable window for the heavy neutrino, consistent with rare decays and BBN constraints, can be found if the fifth neutrino has a mass of order 0.3 - 10 GeV. We also find, however, that the 2011 MB anti-nu data exacerbates the tension with null experiments in both the 3+1 and 3+1+1 models when the lowest energy bins are included, resulting in little improvement in the global fit. We also discuss the implications of an additional neutrino for the reactor and gallium anomalies, and show that an oscillation explanation of the anomalies is disfavored by cosmological considerations, direct searches, and precision electroweak tests.

We propose a novel framework in which the observed baryon and dark matter abundances are simultaneously generated via the Affleck-Dine mechanism. In its simplest realization, Affleck-Dine cogenesis is accomplished by a single superpotential operator and its A-term counterpart. These operators explicitly break B-L and X, the dark matter number, to the diagonal B-L+X. In the early universe these operators stabilize supersymmetric flat directions carrying non-zero B-L and X, and impart the requisite CP violation for asymmetry generation. Because B-L+X is preserved, the resulting B-L and X asymmetries are equal and opposite, though this precise relation may be relaxed if B-L and X are violated separately by additional operators. Our dark matter candidate is stabilized by R-parity and acquires an asymmetric abundance due to its non-zero X number. For a dark matter mass of order a few GeV, one naturally obtains the observed ratio of energy densities today, Omega_{DM}/Omega_B \sim 5. These theories typically predict macroscopic lifetimes for the lightest observable supersymmetric particle as it decays to the dark matter.

We study the dynamics of dark matter (DM) particle-antiparticle oscillations within the context of asymmetric DM. Oscillations arise due to small DM number-violating Majorana-type mass terms, and can lead to recoupling of annihilation after freeze-out and washout of the DM density. We derive the density matrix equations for DM oscillations and freeze-out from first principles using nonequilibrium field theory, and our results are qualitatively different than in previous studies. DM dynamics exhibits particle-vs-antiparticle "flavor" effects, depending on the interaction type, analogous to neutrino oscillations in a medium. "Flavor-sensitive" DM interactions include scattering or annihilation through a new vector boson, while "flavor-blind" interactions include scattering or s-channel annihilation through a new scalar boson, or annihilation to pairs of bosons. In particular, we find that flavor-sensitive annihilation does not recouple when coherent oscillations begin, and that flavor-blind scattering does not lead to decoherence.

We consider possibly observable effects of asymmetric dark matter (ADM) in neutron stars. Since dark matter does not self-annihilate in the ADM scenario, dark matter accumulates in neutron stars, eventually reaching the Chandrasekhar limit and forming a black hole. We focus on the case of scalar ADM, where the constraints from Bose-Einstein condensation and subsequent black hole formation are most severe due to the absence of Fermi degeneracy pressure. We also note that in some portions of this constrained parameter space, non-trivial effects from Hawking radiation can modify our limits. We find that for scalar ADM with mass between 100 keV and 10^5 GeV, the constraint from pulsars in globular clusters on the scattering cross-section with neutrons ranges from \sigma_n < 10^{-45} cm^2 to 10^{-52} cm}^2. In particular, for scalar ADM with mass between 1 GeV and 1 TeV (in the case where black hole evaporation due to Hawking radiation is unimportant), the constraint on the scattering cross-section is below what is reachable with ton scale direct detection experiments.

One of the open questions of modern cosmology is the nature and properties of the Dark Matter halo and its substructures. In this work we study the gravitational effect of dark matter substructures on pulsar timing observations. Since millisecond pulsars are stable and accurate emitters, they have been proposed as plausible astrophysical tools to probe the gravitational effects of dark matter structures. We study this effect on pulsar timing through Shapiro time delay (or Integrated Sachs-Wolfe (ISW) effect) and Doppler effects statistically, showing that the latter dominates the signal. For this task, we relate the power spectrum of pulsar frequency change to the matter power spectrum on small scales, which we compute using the stable clustering hypothesis. We compare this power spectrum with the reach of current and future observations of pulsar timing designed for gravitational wave (GW) detection. Our results show that while current observations are unable to detect these signals, the sensitivity of the upcoming Square Kilometer Array (SKA) is only a factor of few weaker than our optimistic predictions.

We consider current observational constraints on the electromagnetic charge of dark matter. The velocity dependence of the scattering cross-section through the photon gives rise to qualitatively different constraints than standard dark matter scattering through massive force carriers. In particular, recombination epoch observations of dark matter density perturbations require that $\epsilon$, the ratio of the dark matter to electronic charge, is less than $10^{-6}$ for $m_X = 1 GeV$, rising to $\epsilon < 10^{-4}$ for $m_X = 10 TeV$. Though naively one would expect that dark matter carrying a charge well below this constraint could still give rise to large scattering in current direct detection experiments, we show that charged dark matter particles that could be detected with upcoming experiments are expected to be evacuated from the Galactic disk by the Galactic magnetic fields and supernova shock waves, and hence will not give rise to a signal. Thus dark matter with a small charge is likely not a source of a signal in current or upcoming dark matter direct detection experiments.

We present a unified framework for the study of late time cosmic acceleration. Using methods of effective field theory, we show that existing proposals for late time acceleration can be subsumed in a single framework, rather than many compartmentalized theories. We construct the most general action consistent with symmetry principles, derive the background and perturbation evolution equations, and demonstrate that for special choices of our parameters we can reproduce results already existing in the literature. Lastly, we lay the foundation for future work placing phenomenological constraints on the parameters of the effective theory. Although in this paper we focus on late time acceleration, our construction also generalizes the effective field theory of inflation to the scalar-tensor and multi-field case for perturbatively constructed backgrounds.

PAMELA's observation that the cosmic ray positron fraction increases rapidly with energy implies the presence of primary sources of energetic electron-positron pairs. Of particular interest is the possibility that dark matter annihilations in the halo of the Milky Way provide this anomalous flux of antimatter. The recent measurement of the cosmic ray electron spectrum by the Fermi Gamma Ray Space Telescope, however, can be used to constrain the nature of any such dark matter particle. In particular, it has been argued that in order to accommodate the observations of Fermi and provide the PAMELA positron excess, annihilating dark matter particles must be as massive as ~1 TeV or heavier. In this article, we revisit Fermi's electron spectrum measurement within the context of annihilating dark matter, focusing on masses in the range of 100-1000 GeV, and considering effects such as variations in the astrophysical backgrounds from the presence of local cosmic ray accelerators, and the finite energy resolution of the Fermi Gamma Ray Space Telescope. When these factors are taken into account, we find that dark matter particles as light as ~300 GeV can be capable of generating the positron fraction observed by PAMELA.

We present a model of weak scale Dark Matter (DM) where the thermal DM density is set by the lepton asymmetry due to the presence of higher dimension lepton violating operators. In these models there is generically a separation between the annihilation cross-section responsible for the relic abundance (through lepton violating operators) and the annihilation cross-section that is relevant for the indirect detection of DM (through lepton preserving operators). Due to this separation, there is a perceived boost in the annihilation cross-section in the galaxy today relative to that derived for canonical thermal freeze-out. This results in a natural explanation for the observed cosmic ray electron and positron excesses, without resorting to a Sommerfeld enhancement. Generating the indirect signals also sets the magnitude of the direct detection cross-section which implies a signal for the next generation of experiments. More generically these models motivate continued searches for DM with apparently non-thermal annihilation cross-sections. The DM may also play a role in radiatively generating Majorana neutrino masses.

The WMAP haze is an excess in the 22 to 93 GHz frequency bands of WMAP extending about 10 degrees from the galactic center. We show that synchrotron emission from electron-positron pairs injected into the interstellar medium by the galactic population of pulsars with energies in the 1 to 100 GeV range can explain the frequency spectrum of the WMAP haze and the drop in the average haze power with latitude. The same spectrum of high energy electron-positron pairs from pulsars, which gives rise to the haze, may also generate the observed excesses in AMS, HEAT and PAMELA. We discuss the spatial morphology of the pulsar synchrotron signal and its deviation from spherical symmetry, which may provide an avenue to determine the pulsar contribution to the haze.

In this letter, we study the possibility that Kaluza-Klein dark matter in a model with one universal extra dimension is responsible for the recent observations of the PAMELA and ATIC experiments. In this model, the dark matter particles annihilate largely to charged leptons, which enables them to produce a spectrum of cosmic ray electrons and positrons consistent with the PAMELA and ATIC measurements. To normalize to the observed signal, however, large boost factors (~10^3) are required. Despite these large boost factors and significant annihilation to hadronic modes (35%), we find that the constraints from cosmic ray antiproton measurements can be satisfied. Relic abundance considerations in this model force us to consider a rather specific range of masses (approximately 600-900 GeV) which is very similar to the range required to generate the ATIC spectral feature. The results presented here can also be used as a benchmark for model-independent constraints on dark matter annihilation to hadronic modes.

In this letter, we suggest that a nearby clump of 600-1000 GeV neutralinos may be responsible for the excesses recently observed in the cosmic ray positron and electron spectra by the PAMELA and ATIC experiments. Although neutralino dark matter annihilating throughout the halo of the Milky Way is predicted to produce a softer spectrum than is observed, and violate constraints from cosmic ray antiproton measurements, a large nearby (within 1-2 kiloparsecs of the Solar System) clump of annihilating neutralinos can lead to a spectrum which is consistent with PAMELA and ATIC, while also producing an acceptable antiproton flux. Furthermore, the presence of a large dark matter clump can potentially accommodate the very large annihilation rate required to produce the PAMELA and ATIC signals.

We study whether spin-independent scattering of weakly-interacting massive particles (WIMPs) with nuclei can account for the annual modulation signal reported by DAMA. We consider both elastic and inelastic scattering processes. We find that there is a region of WIMP parameter space which can simultaneously accommodate DAMA and the null results of CDMS, CRESST, and XENON. This region corresponds to an ordinary, elastically-scattering WIMP with a standard Maxwell-Boltzmann distribution, a mass 3 GeV < m_{DM} <8 GeV, and a spin-independent cross section with nucleons 3 \times 10^{-41} cm^2 < \sigma_p^{SI} < 5 \times 10^{-39} cm^2. This new region of parameter space depends crucially on the recently discovered effect of channeling on the energy threshold for WIMP detection in the DAMA experiment; without the inclusion of this effect, the DAMA allowed region is essentially closed by null experiments. Such low-mass WIMPs arise in many theories of Beyond the Standard Model physics, from minimal extensions of the MSSM to solutions of the baryon-dark matter coincidence problem. We find that inelastic scattering channels do not open up a significant parameter region consistent with all experimental results. Future experiments with low energy thresholds for detecting nuclear recoils, such as CDMSII-Si and those utilizing ultra-low energy germanium detectors, will be able to probe the DAMA region of parameter space.

Warped string compactifications, characterized by non-singular behavior of the metric in the infrared (IR), feature departures from the usual anti-de Sitter warped extra dimensions. We study the implications of the smooth IR cutoff for Randall-Sundrum (RS) type models. We find that the phenomenology of the KK gravitons (including their masses and couplings) depends sensitively on the precise shape of the warp factor in the IR. In particular, we analyze the warped deformed conifold and find that the spectrum differs significantly from that of RS, and present a simple prescription (a mass gap ansatz) which can be used to study the phenomenology of IR modifications to 5-d warped extra dimensions.

In several classes of particle physics models -- ranging from the classical Majoron models, to the more recent scenarios of late neutrino masses or Mass-Varying Neutrinos -- one or more of the neutrinos are postulated to couple to a new light scalar field. As a result of this coupling, neutrinos in the early universe instead of streaming freely could form a self-coupled fluid, with potentially observable signatures in the Cosmic Microwave Background and the large scale structure of the universe. We re-examine the constraints on this scenario from the presently available cosmological data and investigate the sensitivity expected from the Planck satellite. In the first case, we find that the sensitivity strongly depends on which piece of data is used. The SDSS Main sample data, combined with WMAP and other data, disfavors the scenario of three coupled neutrinos at about the 3.5$\sigma$ confidence level, but also favors a high number of freely streaming neutrinos, with the best fit at 5.2. If the matter power spectrum is instead taken from the SDSS Large Red Galaxy sample, best fit point has 2.5 freely streaming neutrinos, but the scenario with three coupled neutrinos becomes allowed at $2\sigma$. In contrast, Planck alone will exclude even a single self-coupled neutrino at the $4.2\sigma$ confidence level, and will determine the total radiation at CMB epoch to $\Delta N_\nu^{eff} = ^{+0.5}_{-0.3}$ ($1\sigma$ errors). We investigate the robustness of this result with respect to the details of Planck's detector. This sensitivity to neutrino free-streaming implies that Planck will be capable of probing a large region of the Mass-Varying Neutrino parameter space. Planck may also be sensitive to a scale of neutrino mass generation as high as 1 TeV.

WIMPs with electroweak scale masses (neutralinos, etc.) remain in kinetic equilibrium with other particle species until temperatures approximately in the range of 10 MeV to 1 GeV, leading to the formation of dark matter substructure with masses as small as $10^{-4} M_{\odot}$ to $10^{-12} M_{\odot}$. However, if dark matter consists of particles with MeV scale masses, as motivated by the observation of 511 keV emission from the Galactic Bulge, such particles are naturally expected to remain in kinetic equilibrium with the cosmic neutrino background until considerably later times. This would lead to a strong suppression of small scale structure with masses below about $10^7 M_{\odot}$ to $10^4 M_{\odot}$. This cutoff scale has important implications for present and future searches for faint Local Group satellite galaxies and for the missing satellites problem.

In concordance cosmology, dark matter density perturbations generated by inflation lead to nonlinear, virialized minihalos, into which baryons collapse at redshift $z \sim 20$. We survey here novel baryon evolution produced by a modification of the power spectrum from white noise density perturbations at scales below $k \sim 10 h {Mpc}^{-1}$ (the smallest scales currently measured with the Lyman-$\alpha$ forest). Exotic dark matter dynamics, such as would arise from scalar dark matter with a late phase transition (similar to an axion, but with lower mass), create such an amplification of small scale power. The dark matter produced in such a phase transition collapses into minihalos, with a size given by the dark matter mass within the horizon at the phase transition. If the mass of the initial minihalos is larger than $\sim 10^{-3} M_\odot$, the modified power spectrum is found to cause widespread baryon collapse earlier than standard $\Lambda$CDM, leading to earlier gas heating. It also results in higher spin temperature of the baryons in the 21 cm line relative to $\Lambda$CDM at redshifts $z > 20$ if the mass of the minihalo is larger than $1 M_\odot$. It is estimated that experiments probing 21 cm radiation at high redshift will contribute a significant constraint on dark matter models of this type for initial minihalos larger than $\sim 10 M_\odot$. Early experiments reaching to $z\approx 15$ will constrain minihalos down to $\sim 10^3 M_\odot$.

We model the formation, evolution and astrophysical effects of dark compact Scalar Miniclusters (``ScaMs''). These objects arise when a scalar field, with an axion-like or Higgs-like potential, undergoes a second order phase transition below the QCD scale. Such a scalar field may couple too weakly to the standard model to be detectable directly through particle interactions, but may still be detectable by gravitational effects, such as lensing and baryon accretion by large, gravitationally bound miniclusters. The masses of these objects are shown to be constrained by the Ly$\alpha$ power spectrum to be less than $\sim 10^4 M_\odot$, but they may be as light as classical axion miniclusters, of the order of $10^{-12} M_\odot$. We simulate the formation and nonlinear gravitational collapse of these objects around matter-radiation equality using an N-body code, estimate their gravitational lensing properties, and assess the feasibility of studying them using current and future lensing experiments. Future MACHO-type variability surveys of many background sources can reveal either high-amplification, strong lensing events, or measure density profiles directly via weak-lensing variability, depending on ScaM parameters and survey depth. However, ScaMs, due to their low internal densities, are unlikely to be responsible for apparent MACHO events already detected in the Galactic halo. A simple estimate is made of parameters that would give rise to early structure formation; in principle, early stellar collapse could be triggered by ScaMs as early as recombination, and significantly affect cosmic reionization.

We show that axion phenomenology may be significantly different than conventionally assumed in theories which exhibit late phase transitions (below the QCD scale). In such theories one can find multiple pseudoscalars with axion-like couplings to matter, including a string scale axion, whose decay constant far exceeds the conventional cosmological bound. Such theories have several dark matter candidates.