We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

We provide expressions for the nonperturbative matching of the effective field theory describing dark matter interactions with quarks and gluons to the effective theory of nonrelativistic dark matter interacting with nonrelativistic nucleons. We give the leading and subleading order expressions in chiral counting. In general, a single partonic operator already matches onto several nonrelativistic operators at leading order in chiral counting. Thus, keeping only one operator at the time in the nonrelativistic effective theory does not properly describe the scattering in direct detection. Moreover, the matching of the axial--axial partonic level operator, as well as the matching of the operators coupling DM to the QCD anomaly term, naively include momentum suppressed terms. However, these are still of leading chiral order due to pion poles and can be numerically important. We illustrate the impact of these effects with several examples.

We present the effective field theory for dark matter interactions with the visible sector that is valid at scales of O(1 GeV). Starting with an effective theory describing the interactions of fermionic and scalar dark matter with quarks, gluons and photons via higher dimension operators that would arise from dimension-five and dimension-six operators above electroweak scale, we perform a nonperturbative matching onto a heavy baryon chiral perturbation theory that describes dark matter interactions with light mesons and nucleons. This is then used to obtain the coefficients of the nuclear response functions using a chiral effective theory description of nuclear forces. Our results consistently keep the leading contributions in chiral counting for each of the initial Wilson coefficients.

We provide a Mathematica package, DirectDM, that takes as input the Wilson coefficients of the relativistic effective theory describing the interactions of dark matter with quarks, gluons and photons, and matches it onto an effective theory describing the interactions of dark matter with neutrons and protons. The nonperturbative matching is performed at leading order in a chiral expansion. The one-loop QCD and QED renormalization-group evolution from the electroweak scale down to the hadronic scale, as well as finite corrections at the heavy quark thresholds are taken into account. We also provide an interface with the package DMFormFactor so that, starting from the relativistic effective theory, one can directly obtain the event rates for direct detection experiments.

We explore the relationship between astrophysical gamma-ray signals and LHC signatures for a class of phenomenologically successful secluded dark matter models, motivated by recent evidence for a ~130 GeV gamma-ray line. We consider in detail scenarios in which interactions between the dark sector and the standard model are mediated by a vev-less scalar field \phi, transforming as an N-plet (N > 3) under SU(2)_L. Since some of the component fields of \phi carry large electric charges, loop induced dark matter annihilation to \gamma \gamma and \gamma Z can be enhanced without the need for non-perturbatively large couplings, and without overproduction of continuum gamma-rays from other final states. We discuss prospects for other experimental tests, including dark matter-nucleon scattering and production of \phi at the LHC, where future searches for anomalous charged tracks may be sensitive. The first LHC hints could come from the Higgs sector, where loop corrections involving \phi lead to significantly modified h to \gamma \gamma and h to \gamma Z branching ratios.

We present halo-independent methods to analyze the results of dark matter direct detection experiments assuming inelastic scattering. We focus on the annual modulation signal reported by DAMA/LIBRA and present three different halo-independent tests. First, we compare it to the upper limit on the unmodulated rate from XENON100 using (a) the trivial requirement that the amplitude of the annual modulation has to be smaller than the bound on the unmodulated rate, and (b) a bound on the annual modulation amplitude based on an expansion in the Earth's velocity. The third test uses the special predictions of the signal shape for inelastic scattering and allows for an internal consistency check of the data without referring to any astrophysics. We conclude that a strong conflict between DAMA/LIBRA and XENON100 in the framework of spin-independent inelastic scattering can be established independently of the local properties of the dark matter halo.

We show how constraints on the time integrated event rate from a given dark matter (DM) direct detection experiment can be used to set a stringent constraint on the amplitude of the annual modulation signal in another experiment. The method requires only very mild assumptions about the properties of the local DM distribution: that it is temporally stable on the scale of months and spatially homogeneous on the ecliptic. We apply the method to the annual modulation signal in DAMA/LIBRA, which we compare to the bounds derived from the constraints on the time-averaged rates from XENON10, XENON100, CDMS and SIMPLE. Assuming a DM mass of 10 GeV, we show that a DM interpretation of the DAMA/LIBRA signal is excluded at 6.3sigma (4.6sigma) for isospin conserving (violating) spin-independent interactions, and at 4.9sigma for spin-dependent interactions on protons.

We derive constraints on the annual modulation signal in Dark Matter (DM) direct detection experiments in terms of the unmodulated event rate. A general bound independent of the details of DM distribution follows from the assumption that the motion of the earth around the sun is the only source of time variation. The bound is valid for a very general class of particle physics models and also holds in the presence of an unknown unmodulated background. More stringent bounds are obtained, if modest assumptions on symmetry properties of the DM halo are adopted. We illustrate the bounds by applying them to the annual modulation signals reported by the DAMA and CoGeNT experiments in the framework of spin-independent elastic scattering. While the DAMA signal satisfies our bounds, severe restrictions on the DM mass can be set for CoGeNT.

We perform a global fit to data from Dark Matter (DM) direct detection experiments, including the recent CDMS-II results. We discuss possible interpretations of the DAMA annual modulation signal in terms of spin-independent and spin-dependent DM-nucleus interactions, both for elastic and inelastic scattering. We find that in the spin-dependent inelastic scattering off protons a good fit to all data is obtained. We present a simple toy model realizing such a scenario. In all the remaining cases the DAMA allowed regions are disfavored by other experiments or suffer from severe fine tuning of DM parameters with respect to the galactic escape velocity. Finally, we also entertain the possibility that the two events observed in CDMS-II are an actual signal of elastic DM scattering, and we compare the resulting CDMS-II allowed regions to the exclusion limits from other experiments. In this arXiv version of the manuscript we also provide in appendix A the updated fits including recent CoGeNT results.

We consider the hypothesis that Dark Matter (DM) has tree-level interactions only with leptons. Such a framework, where DM recoils against electrons bound in atoms, has been proposed as an explanation for the annually modulated scintillation signal in DAMA/LIBRA versus the absence of a signal for nuclear recoils in experiments like CDMS or XENON10. However, even in such a leptophilic DM scenario there are loop induced DM-hadron interactions, where the photons emitted from virtual leptons couple to the charge of a nucleus. Using a general effective field theory approach we show that, if such an interaction is induced at one or two-loop level, then DM-nucleus scattering dominates over DM-electron scattering. This is because the latter is suppressed by the bound state wave function. One obtains a situation similar to standard DM-nucleus scattering analyses with considerable tension between the results of DAMA and CDMS/XENON10. This conclusion does not apply in the case of pseudoscalar or axial vector coupling between DM and leptons, where the loop diagrams vanish. In this case the explanation of the DAMA signal in terms of DM-electron scattering is strongly disfavored by the spectral shape of the signal. Furthermore, if DM can annihilate into neutrinos or tau leptons, the required cross sections are excluded by many orders of magnitude using the Super-Kamiokande bound on neutrinos from DM annihilations in the Sun.

We propose a scenario where Dark Matter (DM) annihilates into an intermediate state which travels a distance $\lambda \equiv v/\Gamma$ on the order of galactic scales and then decays to Standard Model (SM) particles. The long lifetime disperses the production zone of the SM particles away from the galactic center and hence, relaxes constraints from gamma ray observations on canonical annihilation scenarios. We utilize this set up to explain the electron and positron excesses observed recently by PAMELA, ATIC, and FERMI. While an explanation in terms of usual DM annihilations seems to conflict with gamma ray observations, we show that within the proposed scenario, the PAMELA/ATIC/FERMI results are consistent with the gamma ray data. The distinction from decay scenarios is discsussed and we comment on the prospects for DM production at LHC. The typical decay length $\lambda \gtrsim 10$ kpc of the intermediate state can have its origin from a dimension six operator suppressed by a scale $\Lambda \sim 10^{13}$ GeV, which is roughly the seesaw scale for neutrino masses.

We explain the PAMELA positron excess and the PPB-BETS/ATIC e+ + e- data using a simple two component dark matter model (2DM). The two particle species in the dark matter sector are assumed to be in thermal equilibrium in the early universe. While one particle is stable and is the present day dark matter, the second one is metastable and decays after the universe is 10^-8 s old. In this model it is simple to accommodate the large boost factors required to explain the PAMELA positron excess without the need for large spikes in the local dark matter density. We provide the constraints on the parameters of the model and comment on possible signals at future colliders.

Recently, Csaki, Kaloper and Terning (CKT) suggested that axion-photon oscillation in the intergalactic medium can explain the observed dimming of distant type Ia supernovae. This mechanism works only if the initial axion flux is much smaller than the initial photon flux. We study several mechanisms for such initial axion production. The direct axion production in the supernovae, the photon-axion oscillation in the magnetic field of supernovae and in the magnetic field of their host galaxies are addressed. We find it likely that the initial axion flux is very small and therefore does not pose a problem to the CKT mechanism.