Paleo-detectors are a proposed experimental technique in which one would search for traces of recoiling nuclei in ancient minerals. Natural minerals on Earth are as old as $\mathcal{O}(1)\,$Gyr and, in many minerals, the damage tracks left by recoiling nuclei are also preserved for timescales long compared to 1 Gyr once created. Thus, even reading out relatively small target samples of order 100 g, paleo-detectors would allow one to search for very rare events thanks to the large exposure, $\varepsilon \sim 100\,{\rm g}\,{\rm Gyr} = 10^5\,{\rm t}\,{\rm yr}$. Here, we explore the potential of paleo-detectors to measure nuclear recoils induced by neutrinos from galactic core collapse supernovae. We find that they would not only allow for a direct measurement of the average core collapse supernova rate in the Milky Way, but would also contain information about the time-dependence of the local supernova rate over the past $\sim$1 Gyr. Since the supernova rate is thought to be directly proportional to the star formation rate, such a measurement would provide a determination of the local star formation history. We investigate the sensitivity of paleo-detectors to both a smooth time evolution and an enhancement of the core collapse supernova rate on relatively short timescales, as would be expected for a starburst period in the local group.

We propose a multi-messenger probe of QCD axion Dark Matter based on observations of black hole-neutron star binary inspirals. It is suggested that a dense Dark Matter spike may grow around intermediate mass black holes ($10^{3}-10^{5} \mathrm{\,M_{\odot}}$). The presence of such a spike produces two unique effects: a distinct phase shift in the gravitational wave strain during the inspiral and an enhancement of the radio emission due to the resonant axion-photon conversion occurring in the neutron star magnetosphere throughout the inspiral and merger. Remarkably, the observation of the gravitational wave signal can be used to infer the Dark Matter density and, consequently, to predict the radio emission. We study the projected reach of the LISA interferometer and next-generation radio telescopes such as the Square Kilometre Array. Given a sufficiently nearby system, such observations will potentially allow for the detection of QCD axion Dark Matter in the mass range $10^{-7}\,\mathrm{eV}$ to $10^{-5}\,\mathrm{eV}$.

Despite strong evidence for the existence of large amounts of dark matter (DM) in our Universe, there is no direct indication of its presence in our own solar system. All estimates of the local DM density, crucial for all direct DM searches, rely on extrapolating estimates on much larger scales. We demonstrate for the first time the possibility of measuring the local DM density with a direct detection experiment. It relies on the assumption that incoming DM particles frequently scatter on terrestrial nuclei prior to detection, inducing an additional time-dependence of the signal. We show that for sub-GeV DM, with a spin-independent DM-proton cross section $\sim10^{-32}\,\text{cm}^2$, future direct detection experiments should be able to reconstruct the local DM density with 20% uncertainty.

A dark matter overdensity around a black hole may significantly alter the dynamics of the black hole's merger with another compact object. We consider here intermediate mass-ratio inspirals of stellar-mass compact objects with intermediate-mass black holes "dressed" with dark matter. We first demonstrate that previous estimates based on a fixed dark-matter dress are unphysical for a range of binaries and dark-matter distributions by showing that the total energy dissipated by the compact object through dynamical friction, as it inspirals through the dense dark matter environment towards the black hole, is larger than the gravitational binding energy of the dark-matter dress itself. We then introduce a new formalism that allows us to self-consistently follow the evolution of the dark-matter dress due to its gravitational interaction with the binary. We show that the dephasing of the gravitational waveform induced by dark matter is smaller than previously thought, but is still potentially detectable with the LISA space interferometer. The gravitational waves from such binaries could provide powerful diagnostics of the particle nature of dark matter.

Observational constraints on gamma rays produced by the annihilation of weakly interacting massive particles around primordial black holes (PBHs) imply that these two classes of Dark Matter candidates cannot coexist. We show here that the successful detection of one or more PBHs by radio searches (with the Square Kilometer Array) and gravitational waves searches (with LIGO/Virgo and the upcoming Einstein Telescope) would set extraordinarily stringent constraints on virtually all weak-scale extensions of the Standard Model with stable relics, including those predicting a WIMP abundance much smaller than that of Dark Matter. Upcoming PBHs searches have in particular the potential to rule out almost the entire parameter space of popular theories such as the minimal supersymmetric standard model and scalar singlet Dark Matter.

Dark matter substructure can contribute significantly to local dark matter searches and may provide a large uncertainty in the interpretation of those experiments. For direct detection experiments, sub-halos give rise to an additional dark matter component on top of the smooth dark matter distribution of the host halo. In the case of dark matter capture in the Sun, sub-halo encounters temporarily increase the number of captured particles. Even if the encounter happened in the past, the number of dark matter particles captured by the Sun can still be enhanced today compared to expectations from the host halo as those enhancements decay over time. Using results from an analytical model of the sub-halo population of a Milky Way-like galaxy, valid for sub-halo masses between $10^{-5}\,M_\odot$ and $10^{11}\,M_\odot$, we assess the impact of sub-halos on direct dark matter searches in a probabilistic way. We find that the impact on direct detection can be sizable, with a probability of $\sim 10^{-3}$ to find an $\mathcal{O}(1)$ enhancement of the recoil rate. In the case of the capture rate in the Sun, we find that $\mathcal{O}(1)$ enhancements are very unlikely, with probability $\lesssim 10^{-5}$, and are even impossible for some dark matter masses.

We study the prospects for the Large Synoptic Survey Telescope (LSST) to find new dwarf spheroidal galaxies in the Milky Way. Adopting models of Milky-Way halo substructure and phenomenological prescriptions connecting subhalos and satellite galaxies, we obtain surface brightness distributions of $V$-band magnitude that lead us to predict that LSST will discover tens to hundreds of dwarf spheroidal galaxies above its sensitivity. The soon-to-be-discovered dwarfs will be interesting targets for indirect searches of dark matter annihilation yields. We forecast the distribution function of gamma-ray emission from dark matter annihilation in these objects, and discuss the detectability of these signals at both Fermi Large Area Telescope (LAT) and Cherenkov Telescope Array (CTA). By combining information from the predicted dwarf galaxies, we obtain an expected sensitivity to the annihilation cross section $\langle \sigma v \rangle$ of $10^{-26}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 10 GeV with Fermi-LAT) and $5\times 10^{-24}$ cm$^3$ s$^{-1}$ (for dark matter particles of mass 500 GeV with CTA). We find that the current uncertainties in the mass measurement of the Milky-Way halo are relatively minor compared with the Poisson errors associated to drawing the most promising dwarfs from the underlying flux distribution.

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.

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

Even if Dark Matter (DM) is neutral under electromagnetism, it can still interact with the Standard Model (SM) via photon exchange from higher-dimensional operators. Here we classify the general effective operators coupling DM to photons, distinguishing between Dirac/Majorana fermion and complex/real scalar DM. We provide model-independent constraints on these operators from direct and indirect detection. We also constrain various DM-lepton operators, which induce DM-photon interactions via RG running or which typically arise in sensible UV-completions. This provides a simple way to quickly assess constraints on any DM model that interacts mainly via photon exchange or couples to SM leptons.

The EDELWEISS collaboration has performed a search for dark matter particles with masses below the GeV-scale with a 33.4-g germanium cryogenic detector operated in a surface lab. The energy deposits were measured using a neutron-transmutation-doped Ge thermal sensor with a 17.7~eV (RMS) baseline heat energy resolution leading to a 60~eV analysis energy threshold. Despite a moderate lead shielding and the high-background environment, the first sub-GeV spin-independent dark matter limit based on a germanium target has been achieved. The experiment provides the most stringent, nuclear recoil based, above-ground limit on spin-independent interactions above 600~MeV/c$^{2}$. The experiment also provides the most stringent limits on spin-dependent interactions with protons and neutrons below 1.3~GeV/c$^{2}$. Furthermore, the dark matter search results were studied in the context of Strongly Interacting Massive Particles, taking into account Earth-shielding effects, for which new regions of the available parameter space have been excluded. Finally, the dark matter search has also been extended to interactions via the Migdal effect, resulting for the first time in the exclusion of particles with masses between 45 and 150~MeV/c$^{2}$ with spin-independent cross sections ranging from $10^{-29}$ to $10^{-26}$~cm$^2$.

Paleo-detectors are a recently proposed method for the direct detection of Dark Matter (DM). In such detectors, one would search for the persistent damage features left by DM--nucleus interactions in ancient minerals. Initial sensitivity projections have shown that paleo-detectors could probe much of the remaining Weakly Interacting Massive Particle (WIMP) parameter space. In this paper, we improve upon the cut-and-count approach previously used to estimate the sensitivity by performing a full spectral analysis of the background- and DM-induced signal spectra. We consider two scenarios for the systematic errors on the background spectra: i) systematic errors on the normalization only, and ii) systematic errors on the shape of the backgrounds. We find that the projected sensitivity is rather robust to imperfect knowledge of the backgrounds. Finally, we study how well the parameters of the true WIMP model could be reconstructed in the hypothetical case of a WIMP discovery.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

The formation of astrophysical and primordial black holes influences the distribution of dark matter surrounding them. Black holes are thus expected to carry a dark matter `dress' whose properties depend on their formation mechanism and on the properties of the environment. Here we carry out a numerical and analytical study of the merger of dressed black holes, and show that the distribution of dark matter around them dramatically affects the dynamical evolution of the binaries. Although the final impact on the merger rate of primordial black holes is rather small with respect to the case of `naked' black holes, we argue that our analysis places the calculation of this rate on more solid ground, with LIGO-Virgo observations ruling out a dark matter fraction of $10^{-3}$ for primordial black holes of 100 solar masses, and it paves the way to more detailed analyses of environmental effects induced by dark matter on the gravitational wave emission of binary black holes.

The theoretical interpretation of dark matter (DM) experiments is hindered by uncertainties on the dark matter density and velocity distribution inside the Solar System. In order to quantify those uncertainties, we present a parameter that characterizes the deviation of the true velocity distribution from the standard Maxwell-Boltzmann form, and we then determine for different values of this parameter the most aggressive and most conservative limits on the dark matter scattering cross section with nuclei; uncertainties in the local dark matter density can be accounted for trivially. This allows us to bracket, in a model independent way, the impact of astrophysical uncertainties on limits from direct detection experiments and/or neutrino telescopes. We find that current limits assuming the Standard Halo Model are at most a factor of $\sim 2$ weaker than the most aggressive possible constraints. In addition, combining neutrino telescope and direct detection constraints (in a statistically meaningful way), we show that limits on DM in the mass range $\sim 10 - 1000$ GeV cannot be weakened by more than around a factor of 10, for all possible velocity distributions. We finally demonstrate that our approach can also be employed in the event of a DM discovery, allowing us to avoid bias in the reconstruction of the DM properties.

Forecasting the signal discrimination power of dark matter (DM) searches is commonly limited to a set of arbitrary benchmark points. We introduce new methods for benchmark-free forecasting that instead allow an exhaustive exploration and visualization of the phenomenological distinctiveness of DM models, based on standard hypothesis testing. Using this method, we reassess the signal discrimination power of future liquid Xenon and Argon direct DM searches. We quantify the parameter regions where various non-relativistic effective operators, millicharged DM, and magnetic dipole DM can be discriminated, and where upper limits on the DM mass can be found. We find that including an Argon target substantially improves the prospects for reconstructing the DM properties. We also show that only in a small region with DM masses in the range 20-100 GeV and DM-nucleon cross sections a factor of a few below current bounds can near-future Xenon and Argon detectors discriminate both the DM-nucleon interaction and the DM mass simultaneously. In all other regions only one or the other can be obtained.

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta.

We investigate the radiative decay of the cosmic neutrino background, and its impact on the spectrum of the cosmic microwave background (CMB) that is known to be a nearly perfect black body. We derive exact formulae for the decay of a heavier neutrino into a lighter neutrino and a photon, $\nu_j \to \nu_i + \gamma$, and of absorption as its inverse, $\nu_i + \gamma \to \nu_j$, by accounting for the precise form of the neutrino momentum distribution. Our calculations show that if the neutrinos are heavier than $\mathcal O(0.1)$ eV, the exact formulae give results that differ by $\sim$50%, compared with approximate ones where neutrinos are assumed to be at rest. We also find that spectral distortion due to absorption is more important for heavy neutrino masses (by a factor of $\sim$10 going from a neutrino mass of 0.01 eV to 0.1 eV). By analyzing the CMB spectral data measured with COBE-FIRAS, we obtain lower limits on the neutrino lifetime of $\tau_{12} \gtrsim 4 \times 10^{21}$ s (95% C.L.) for the smaller mass splitting and $\tau_{13} \sim \tau_{23} \gtrsim 10^{19}$ s for the larger mass splitting. These represent up to one order of magnitude improvement over previous CMB constraints. With future CMB experiments such as PIXIE, these limits will improve by roughly 4 orders of magnitude. This translates to a projected upper limit on the neutrino magnetic moment (for certain neutrino masses and decay modes) of $\mu_\nu < 3 \times 10^{-11}\, \mu_B$, where $\mu_B$ is the Bohr magneton. Such constraints would make future precision CMB measurements competitive with lab-based constraints on neutrino magnetic moments.

Direct searches for Dark Matter (DM) are continuously improving, probing down to lower and lower DM-nucleon interaction cross sections. For strongly-interacting massive particle (SIMP) Dark Matter, however, the accessible cross section is bounded from above due to the stopping effect of the atmosphere, Earth and detector shielding. We present a careful calculation of the SIMP signal rate, focusing on super-heavy DM ($m_\chi \gtrsim 10^5 \,\,\mathrm{GeV}$) for which the standard nuclear-stopping formalism is applicable, and provide code for implementing this calculation numerically. With recent results from the low-threshold CRESST 2017 surface run, we improve the maximum cross section reach of direct detection searches by a factor of around 5000, for DM masses up to $10^8 \,\,\mathrm{GeV}$. A reanalysis of the longer-exposure, sub-surface CDMS-I results (published in 2002) improves the previous cross section reach by two orders of magnitude, for masses up to $10^{15} \,\,\mathrm{GeV}$. Along with complementary constraints from SIMP capture and annihilation in the Earth and Sun, these improved limits from direct nuclear scattering searches close a number of windows in the SIMP parameter space in the mass range $10^6$ GeV to $10^{13}$ GeV, of particular interest for heavy DM produced gravitationally at the end of inflation.

The analysis of signals in directional dark matter (DM) detectors typically assumes that the directions of nuclear recoils can be measured in the Galactic rest frame. However, this is not possible with all directional detection technologies. In nuclear emulsions, for example, the recoil events must be detected and measured after the exposure time of the experiment. Unless the entire detector is mounted and rotated with the sidereal day, the recoils cannot be reoriented in the Galactic rest frame. We examine the effect of this `time integration' on the primary goals of directional detection, namely: (1) confirming that the recoils are anisotropic; (2) measuring the median recoil direction to confirm their Galactic origin; and (3) probing below the neutrino floor. We show that after time integration the DM recoil distribution retains a preferred direction and is distinct from that of Solar neutrino-induced recoils. Many of the advantages of directional detection are therefore preserved and it is not crucial to mount and rotate the detector. Rejecting isotropic backgrounds requires a factor of 2 more signal events compared with an experiment with event time information, whereas a factor of 1.5-3 more events are needed to measure a median direction in agreement with the expectation for DM. We also find that there is still effectively no neutrino floor in a time-integrated directional experiment. However to reach a cross section an order of magnitude below the floor, a factor of 8 larger exposure is required than with a conventional directional experiment. We also examine how the sensitivity is affected for detectors with only 2D recoil track readout, and/or no head-tail measurement. As for non-time-integrated experiments, 2D readout is not a major disadvantage, though a lack of head-tail sensitivity is.

It was recently pointed out that direct detection signals from at least three different targets may be used to determine whether the Dark Matter (DM) particle is different from its antiparticle. In this work, we examine in detail the feasibility of this test under different conditions, motivated by proposals for future detectors. Specifically, we perform likelihood fits to mock data under the hypotheses that the DM particle is identical to or different from its antiparticle, and determine the significance with which the former can be rejected in favor of the latter. In our analysis, we consider 3 different values of the DM mass ($50$ GeV, $300$ GeV, $1$ TeV) and 4 different experimental ensembles, each consisting of at least 3 different targets -- Xe and Ar plus one among the following: Si, Ge, $\mathrm{CaWO_4}$, or Ge/$\mathrm{CaWO_4}$. For each of these experimental ensembles and each DM mass, the expected discrimination significance is calculated as a function of the DM-nucleon couplings. In the best case scenario, the discrimination significance can reach $\mathcal{O}(3\sigma)$ for three of the four ensembles considered, and $\mathcal{O}(5\sigma)$ for the ensemble including $\mathrm{Si}$, highlighting the need for a variety of experimental targets in order to determine the DM properties. These results show that future direct detection signals could be used to exclude, at a statistically significant level, a Majorana or a real DM particle, giving a critical clue about the identity of the Dark Matter.

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.

Direct detection experiments search for the interactions of Dark Matter (DM) particles with nuclei in terrestrial detectors. But if these interactions are sufficiently strong, DM particles may scatter in the Earth, affecting their distribution in the lab. We present a new analytic calculation of this `Earth-scattering' effect in the regime where DM particles scatter at most once before reaching the detector. We perform the calculation self-consistently, taking into account not only those particles which are scattered away from the detector, but also those particles which are deflected towards the detector. Taking into account a realistic model of the Earth and allowing for a range of DM-nucleon interactions, we present the EarthShadow code, which we make publicly available, for calculating the DM velocity distribution after Earth-scattering. Focusing on low-mass DM, we find that Earth-scattering reduces the direct detection rate at certain detector locations while increasing the rate in others. The Earth's rotation induces a daily modulation in the rate, which we find to be highly sensitive to the detector latitude and to the form of the DM-nucleon interaction. These distinctive signatures would allow us to unambiguously detect DM and perhaps even identify its interactions in regions of the parameter space within the reach of current and future experiments.

Directionally sensitive dark matter (DM) direct detection experiments present the only way to observe the full three-dimensional velocity distribution of the Milky Way halo local to Earth. In this work we compare methods for extracting information about the local DM velocity distribution from a set of recoil directions and energies in a range of hypothetical directional and non-directional experiments. We compare a model independent empirical parameterisation of the velocity distribution based on an angular discretisation with a model dependent approach which assumes knowledge of the functional form of the distribution. The methods are tested under three distinct halo models which cover a range of possible phase space structures for the local velocity distribution: a smooth Maxwellian halo, a tidal stream and a debris flow. In each case we use simulated directional data to attempt to reconstruct the shape and parameters describing each model as well as the DM particle properties. We find that the empirical parametrisation is able to make accurate unbiased reconstructions of the DM mass and cross section as well as capture features in the underlying velocity distribution in certain directions without any assumptions about its true functional form. We also find that by extracting directionally averaged velocity parameters with this method one can discriminate between halo models with different classes of substructure.

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.

Cosmological observations indicate that most of the matter in the Universe is Dark Matter. Dark Matter in the form of Weakly Interacting Massive Particles (WIMPs) can be detected directly, via its elastic scattering off target nuclei. Most current direct detection experiments only measure the energy of the recoiling nuclei. However, directional detection experiments are sensitive to the direction of the nuclear recoil as well. Due to the Sun's motion with respect to the Galactic rest frame, the directional recoil rate has a dipole feature, peaking around the direction of the Solar motion. This provides a powerful tool for demonstrating the Galactic origin of nuclear recoils and hence unambiguously detecting Dark Matter. Furthermore, the directional recoil distribution depends on the WIMP mass, scattering cross section and local velocity distribution. Therefore, with a large number of recoil events it will be possible to study the physics of Dark Matter in terms of particle and astrophysical properties. We review the potential of directional detectors for detecting and characterizing WIMPs.

Dark matter (DM) direct detection experiments which are directionally-sensitive may be the only method of probing the full velocity distribution function (VDF) of the Galactic DM halo. We present an angular basis for the DM VDF which can be used to parametrise the distribution in order to mitigate astrophysical uncertainties in future directional experiments and extract information about the DM halo. This basis consists of discretising the VDF in a series of angular bins, with the VDF being only a function of the DM speed $v$ within each bin. In contrast to other methods, such as spherical harmonic expansions, the use of this basis allows us to guarantee that the resulting VDF is everywhere positive and therefore physical. We present a recipe for calculating the event rates corresponding to the discrete VDF for an arbitrary number of angular bins $N$ and investigate the discretisation error which is introduced in this way. For smooth, Standard Halo Model-like distribution functions, only $N=3$ angular bins are required to achieve an accuracy of around $10-30\%$ in the number of events in each bin. Shortly after confirmation of the DM origin of the signal with around 50 events, this accuracy should be sufficient to allow the discretised velocity distribution to be employed reliably. For more extreme VDFs (such as streams), the discretisation error is typically much larger, but can be improved with increasing $N$. This method paves the way towards an astrophysics-independent analysis framework for the directional detection of dark matter.

The framework of non-relativistic effective field theory (NREFT) aims to generalise the standard analysis of direct detection experiments in terms of spin-dependent (SD) and spin-independent (SI) interactions. We show that a number of NREFT operators lead to distinctive new directional signatures, such as prominent ring-like features in the directional recoil rate, even for relatively low mass WIMPs. We discuss these signatures and how they could affect the interpretation of future results from directional detectors. We demonstrate that considering a range of possible operators introduces a factor of 2 uncertainty in the number of events required to confirm the median recoil direction of the signal. Furthermore, using directional detection, it is possible to distinguish the more general NREFT interactions from the standard SI/SD interactions at the $2\sigma$ level with $\mathcal{O}(100-500)$ events. In particular, we demonstrate that for certain NREFT operators, directional sensitivity provides the only method of distinguishing them from these standard operators, highlighting the importance of directional detectors in probing the particle physics of dark matter.

With positive signals from multiple direct detection experiments it will, in principle, be possible to measure the mass and cross sections of weakly-interacting massive particle (WIMP) dark matter. Recent work has shown that, with a polynomial parameterisation of the WIMP speed distribution, it is possible to make an unbiased measurement of the WIMP mass, without making any astrophysical assumptions. However, direct detection experiments are not sensitive to low-speed WIMPs and, therefore, any model-independent approach will lead to a bias in the cross section. This problem can be solved with the addition of measurements of the flux of neutrinos from the Sun. This is because the flux of neutrinos produced from the annihilation of WIMPs which have been gravitationally captured in the Sun is sensitive to low-speed WIMPs. Using mock data from next-generation direct detection experiments and from the IceCube neutrino telescope, we show that the complementary information from IceCube on low-speed WIMPs breaks the degeneracy between the cross section and the speed distribution. This allows unbiased determinations of the WIMP mass and spin-independent and spin-dependent cross sections to be made, and the speed distribution to be reconstructed. We use two parameterisations of the speed distribution: binned and polynomial. While the polynomial parameterisation can encompass a wider range of speed distributions, this leads to larger uncertainties in the particle physics parameters.

In a recent paper, a new parametrization for the dark matter (DM) speed distribution f(v) was proposed for use in the analysis of data from direct detection experiments. This parametrization involves expressing the logarithm of the speed distribution as a polynomial in the speed v. We present here a more detailed analysis of the properties of this parametrization. We show that the method leads to statistically unbiased mass reconstructions and exact coverage of credible intervals. The method performs well over a wide range of DM masses, even when finite energy resolution and backgrounds are taken into account. We also show how to select the appropriate number of basis functions for the parametrization. Finally, we look at how the speed distribution itself can be reconstructed, and how the method can be used to determine if the data are consistent with some test distribution. In summary, we show that this parametrization performs consistently well over a wide range of input parameters and over large numbers of statistical ensembles and can therefore reliably be used to reconstruct both the DM mass and speed distribution from direct detection data.

The search for weakly-interacting massive particle (WIMP) dark matter is multi-pronged. Ultimately, the WIMP-dark-matter picture will only be confirmed if different classes of experiments see consistent signals and infer the same WIMP properties. In this work, we review the ideas, methods, and status of direct-detection searches. We focus in particular on extracting WIMP physics (WIMP interactions and phase-space distribution) from direct-detection data in the early discovery days when multiple experiments see of order dozens to hundreds of events. To demonstrate the essential complementarity of different direct-detection experiments in this context, we create mock data intended to represent the data from the near-future Generation 2 experiments. We consider both conventional supersymmetry-inspired benchmark points (with spin-independent and -dependent elastic cross sections just below current limits), as well as benchmark points for other classes of models (inelastic and effective-operator paradigms). We also investigate the effect on parameter estimation of loosening or dropping the assumptions about the local WIMP phase-space distribution. We arrive at two main conclusions. Firstly, teasing out WIMP physics with experiments depends critically on having a wide set of detector target materials, spanning a large range of target nuclear masses and spin-dependent sensitivity. It is also highly desirable to obtain data from low-threshold experiments. Secondly, a general reconstruction of the local WIMP velocity distribution, which will only be achieved if there are multiple experiments using different target materials, is critical to obtaining a robust and unbiased estimate of the WIMP mass.

Determining the dark matter (DM) mass is of paramount importance for understanding dark matter. We present a novel parametrization of the DM speed distribution which will allow the DM mass to be accurately measured using data from Weakly Interacting Massive Particle (WIMP) direct detection experiments. Specifically, we parametrize the natural logarithm of the speed distribution as a polynomial in the speed v. We demonstrate, using mock data from upcoming experiments, that by fitting the WIMP mass and interaction cross-section, along with the polynomial coefficients, we can accurately reconstruct both the WIMP mass and speed distribution. This new method is the first demonstration that an accurate, unbiased reconstruction of the WIMP mass is possible without prior assumptions about the distribution function. We anticipate that this technique will be invaluable in the analysis of future experimental data.

Direct detection experiments searching for weakly interacting massive particle (WIMP) dark matter typically use a simplified model of the Galactic halo to derive parameter constraints. However, there is strong evidence that this Standard Halo Model is not a good approximation to our Galaxy. We discuss previous attempts to extract the WIMP mass, cross-section and speed distribution from direct detection data and show that these lead to significant biases in the reconstructed parameter values. We develop and test an alternative model-independent method based on parametrising the momentum distribution of the WIMPs. This allows us to limit the analysis only to those regions of momentum space to which the experiments are sensitive. The method can be applied to a single experiment to extract information about the distribution function, as well as information on the degenerate WIMP mass and interaction cross-section combined in a single parameter. This degeneracy can be broken by including data from additional experiments, meaning that the WIMP mass and speed distribution can be recovered. We test the momentum parametrisation method using mock datasets from proposed ton-scale direct detection experiments, showing that it exhibits improved coverage properties over previous methods, as well as significantly reduced bias. We are also able to accurately reconstruct the shape of the WIMP speed distribution but distinguishing between different underlying distributions remains difficult.