The low-energy electronic recoil spectrum in XENON1T provides an intriguing hint for potential new physics. At the same time, observations of horizontal branch stars favor the existence of a small amount of extra cooling compared to the one expected from the Standard Model particle content. In this note, we argue that a hidden photon with a mass of $\sim 2.8$ keV and a kinetic mixing of $\sim 10^{-15}$ allows for a good fit to both of these excesses. In this scenario, the signal detected in XENON1T is due to the absorption of hidden photon dark matter particles, whereas the anomalous cooling of horizontal branch stars arises from resonant production of hidden photons in the stellar interior.

The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis.

We present a simple model of two dark matter species with opposite millicharge that can form electrically neutral bound states via the exchange of a massive dark photon. If bound state formation is suppressed at low temperatures, a sub-dominant fraction of millicharged particles remains at late times, which can give rise to interesting features in the 21 cm absorption profile at cosmic dawn. The dominant neutral component, on the other hand, can have dipole interactions with ordinary matter, leading to non-standard signals in direct detection experiments. We identify the parameter regions predicting a percent-level ionisation fraction and study constraints from laboratory searches for dark matter scattering and dark photon decays.

Studies of axion-like particles (ALPs) commonly focus on a single type of interaction, for example couplings only to photons. Most ALP models however predict correlations between different couplings, which change the phenomenology in important ways. For example, an MeV-scale ALP coupled to Standard Model gauge bosons at high energies will in general interact with photons, $W^\pm$ and $Z$ bosons as well as mesons and nucleons at low energies. We study the implications of such scenarios and point out that astrophysical constraints, in particular from SN1987A, may be substantially relaxed, opening up new regions of parameter space that may be explored with laboratory experiments such as NA62.

We study the cosmology and LHC phenomenology of a consistent strongly interacting dark sector coupled to Standard Model particles through a generic vector mediator. We lay out the requirements for the model to be cosmologically viable, identify annihilations into dark vector mesons as the dominant dark matter freeze-out process and discuss bounds from direct detection. At the LHC the model predicts dark showers, which can give rise to semi-visible jets or displaced vertices. Existing searches for di-jet resonances and for missing energy mostly probe the parameter regions where prompt decays are expected and constrain our model despite not being optimised for dark showers. We also estimate the sensitivity of dedicated analyses for semi-visible jets and emphasize the complementarity of different search strategies.

We present the first $N$-body simulations that adapt the equations of smoothed particle hydrodynamics to capture the effect of dark matter self-interactions which are too frequent to be resolved explicitly. The relevant energy transfer equations are derived, the appropriate thermal conductivity is determined and the effects of different smoothing kernels are studied. We apply our framework to simulate the formation of isothermal cores in isolated dark matter haloes and determine the core growth rate as a function of the self-scattering cross section. Our approach may be combined with explicit simulations of rare scatterings in order to simulate accurately the effects of arbitrary dark matter self-interactions in future cosmological simulations.

We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through-wall experiments, haloscopes, helioscopes, the axion relic density, horizontal branch stars, supernova 1987A, white dwarf cooling, and gamma-ray observations. We carry out both frequentist and Bayesian analyses, with and without the inclusion of white dwarf cooling. We explore the degree of fine-tuning present in different models and identify parameter regions where it is possible for QCD axion models to account for both the dark matter in the Universe and the cooling hints, comparing them to specific DFSZ- and KSVZ-type models. We find the most credible parameter regions, allowing us to set (prior-dependent) upper and lower bounds on the axion mass. Our analysis also suggests that QCD axions in this scenario most probably make up a non-negligible but sub-dominant component of the dark matter in the Universe.

New gauge bosons at the MeV scale with tiny gauge couplings (so-called dark photons) can be responsible for the freeze-in production of dark matter and provide a clear target for present and future experiments. We study the effects of thermal mixing between dark photons and Standard Model gauge bosons and of the resulting plasmon decays on dark matter production before and after the electroweak phase transition. In the parameter regions preferred by the observed dark matter relic abundance, the dark photon is sufficiently long-lived to be probed with fixed-target experiments and light enough to induce direct detection signals. Indeed, current limits from XENON1T already constrain Dirac fermion dark matter in the GeV to TeV range produced via the freeze-in mechanism. We illustrate our findings for the case of a $U(1)_{B-L}$ gauge extension and discuss possible generalisations.

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

We present results from N-body simulations of self-interacting dark matter (SIDM) subhalos, which could host ultra-faint dwarf spheroidal galaxies, inside a Milky-Way-like main halo. We find that high-concentration subhalos are driven to gravothermal core collapse, while low-concentration subhalos develop large (kpc-sized) low-density cores, with both effects depending sensitively on the satellite's orbit and the self-interaction cross section over mass $\sigma/m$. The overall effect for $\sigma/m \gtrsim 3 \ \rm cm^2/g$ is to increase the range of inner densities, potentially explaining the observed diversity of Milky Way satellites, which include compact systems like Draco and Segue 1 that are dense in dark matter, and less dense, diffuse systems like Sextans and Crater II. We discuss possible ways of distinguishing SIDM models from collisionless dark matter models using the inferred dark matter densities and stellar sizes of the dwarf spheroidal galaxies.

We study the cosmological evolution and phenomenological properties of scalar bosons in the keV to MeV range that have a tiny mixing with the Standard Model Higgs boson. The mixing determines both the abundance of light scalars produced via the freeze-in mechanism and their lifetime. Intriguingly, the parameters required for such scalars to account for all of the dark matter in the present Universe generically predict lifetimes comparable to the sensitivity of present and future indirect detection experiments. In order to accurately determine the relic abundance of light scalars, we calculate freeze-in yields including effects from finite temperatures and quantum statistics and develop a new approach for solving the Boltzmann equation for number-changing processes in the dark sector. We find that light scalars can potentially explain the anomalous x-ray emission at 3.5 keV, while evading constraints from structure formation and predicting potentially observable self-interaction cross sections.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

Scalar singlet dark matter is one of the simplest and most predictive realisations of the WIMP (weakly-interacting massive particle) idea. Although the model is constrained from all directions by the latest experimental data, it still has viable regions of parameter space. Another compelling aspect of scalar singlets is their ability to stabilise the electroweak vacuum. Indeed, models of scalar dark matter are not low-energy effective theories, but can be valid all the way to the Planck scale. Using the GAMBIT framework, we present the first global fit to include both the low-energy experimental constraints and the theoretical constraints from UV physics, considering models with a scalar singlet charged under either a $\mathbb{Z}_2$ or a $\mathbb{Z}_3$ symmetry. We show that if the model is to satisfy all experimental constraints, completely stabilise the electroweak vacuum up to high scales, and also remain perturbative to those scales, one is driven to a relatively small region of parameter space. This region has a Higgs-portal coupling slightly less than 1, a dark matter mass of 1 to 2 TeV and a spin-independent nuclear scattering cross-section around 10$^{-45}$ cm$^2$.

Dark Matter particles are commonly assumed to be weakly interacting massive particles (WIMPs) with a mass in the GeV to TeV range. However, recent interest has shifted towards lighter WIMPs, which are more difficult to probe experimentally. A detection of sub-GeV WIMPs would require the use of small gap materials in sensors. Using recent estimates of the WIMP mass, we identify the relevant target space towards small gap materials (100-10 meV). Dirac Materials, a class of small- or zero-gap materials, emerge as natural candidates for sensors for Dark Matter detection. We propose the use of informatics tools to rapidly assay materials band structures to search for small gap semiconductors and semimetals, rather than focusing on a few preselected compounds. As a specific example of the proposed strategy, we use the organic materials database (omdb.diracmaterials.org) to identify organic candidates for sensors: the narrow band gap semiconductors BNQ-TTF and DEBTTT with gaps of 40 and 38 meV, and the Dirac-line semimetal (BEDT-TTF)$\cdot$Br which exhibits a tiny gap of $\approx$ 50 meV when spin-orbit coupling is included. We outline a novel and powerful approach to search for dark matter detection sensor materials by means of a rapid assay of materials using informatics tools.

Dark matter (DM) particles with mass in the sub-GeV range are an attractive alternative to heavier weakly-interacting massive particles, but direct detection of such light particles is challenging. If however DM-nucleus scattering leads to ionisation of the recoiling atom, the resulting electron may be detected even if the nuclear recoil is unobservable. We demonstrate that including this effect significantly enhances direct detection sensitivity to sub-GeV DM. Existing experiments set world-leading limits, and future experiments may probe the cross sections relevant for thermal freeze-out.

Dark matter in the cosmological concordance model is parameterised by a single number, describing the covariantly conserved energy density of a non-relativistic fluid. Here we test this assumption in a model-independent and conservative way by considering the possibility that, at any point during the cosmological evolution, dark matter may be converted into a non-interacting form of radiation. This scenario encompasses, but is more general than, the cases where dark matter decays or annihilates into these states. We show that observations of the cosmic microwave background allow to strongly constrain this scenario for any conversion time after big bang nucleosynthesis. We discuss in detail, both from a Bayesian and frequentist point of view, in which sense adding large-scale structure observations may even provide a certain preference for a conversion of dark matter to radiation at late times. Finally we apply our general results to a specific particle physics realisation of such a scenario, featuring late kinetic decoupling and Sommerfeld-enhanced dark matter annihilation. We identify a small part of parameter space that both mitigates the tension between cosmic microwave and large-scale structure data and allows for velocity-dependent dark matter self-interactions strong enough to address the small-scale problems of structure formation.

The relative sensitivity of different direct detection experiments depends sensitively on the astrophysical distribution and particle physics nature of dark matter, prohibiting a model-independent comparison. The situation changes fundamentally if two experiments employ the same target material. We show that in this case one can compare measurements of an annual modulation and exclusion bounds on the total rate while making no assumptions on astrophysics and no (or only very general) assumptions on particle physics. In particular, we show that the dark matter interpretation of the DAMA/LIBRA signal can be conclusively tested with COSINUS, a future experiment employing the same target material. We find that if COSINUS excludes a dark matter scattering rate of about $0.01\,\text{kg}^{-1}\,\text{days}^{-1}$ with an energy threshold of $1.8\,$keV and resolution of $0.2\,$keV, it will rule out all explanations of DAMA/LIBRA in terms of dark matter scattering off sodium and/or iodine.

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

Self-interacting dark matter may have striking astrophysical signatures, such as observable offsets between galaxies and dark matter in merging galaxy clusters. Numerical N-body simulations used to predict such observables typically treat the galaxies as collisionless test particles, a questionable assumption given that each galaxy is embedded in its own dark matter halo. To enable a more accurate treatment we develop an effective description of small dark matter haloes taking into account the two major effects due to dark matter self-scatterings: deceleration and evaporation. We point out that self-scatterings can have a sizeable impact on the trajectories of galaxies, diminishing the separation between galaxies and dark matter in merging clusters. This effect depends sensitively on the underlying particle physics, in particular the angular dependence of the self-scattering cross section, and cannot be predicted from the momentum transfer cross section alone.

Dark matter particles interacting via the exchange of very light spin-0 mediators can have large self-interaction rates and obtain their relic abundance from thermal freeze-out. At the same time, these models face strong bounds from direct and indirect probes of dark matter as well as a number of constraints on the properties of the mediator. We investigate whether these constraints can be consistent with having observable effects from dark matter self-interactions in astrophysical systems. For the case of a mediator with purely scalar couplings we point out the highly relevant impact of low-threshold direct detection experiments like CRESST-II, which essentially rule out the simplest realization of this model. These constraints can be significantly relaxed if the mediator has CP-violating couplings, but then the model faces strong constraints from CMB measurements, which can only be avoided in special regions of parameter space.

We explore the potential of future cryogenic direct detection experiments to determine the properties of the mediator that communicates the interactions between dark matter and nuclei. Due to their low thresholds and large exposures, experiments like CRESST-III, SuperCDMS SNOLAB and EDELWEISS-III will have excellent capability to reconstruct mediator masses in the MeV range for a large class of models. Combining the information from several experiments further improves the parameter reconstruction, even when taking into account additional nuisance parameters related to background uncertainties and the dark matter velocity distribution. These observations may offer the intriguing possibility of studying dark matter self-interactions with direct detection experiments.

Weakly-coupled TeV-scale particles may mediate the interactions between normal matter and dark matter. If so, the LHC would produce dark matter through these mediators, leading to the familiar "mono-X" search signatures, but the mediators would also produce signals without missing momentum via the same vertices involved in their production. This document from the LHC Dark Matter Working Group suggests how to compare searches for these two types of signals in case of vector and axial-vector mediators, based on a workshop that took place on September 19/20, 2016 and subsequent discussions. These suggestions include how to extend the spin-1 mediated simplified models already in widespread use to include lepton couplings. This document also provides analytic calculations of the relic density in the simplified models and reports an issue that arose when ATLAS and CMS first began to use preliminary numerical calculations of the dark matter relic density in these models.

The interpretation of dark matter direct detection experiments is complicated by the fact that neither the astrophysical distribution of dark matter nor the properties of its particle physics interactions with nuclei are known in detail. To address both of these issues in a very general way we develop a new framework that combines the full formalism of non-relativistic effective interactions with state-of-the-art halo-independent methods. This approach makes it possible to analyse direct detection experiments for arbitrary dark matter interactions and quantify the goodness-of-fit independent of astrophysical uncertainties. We employ this method in order to demonstrate that the degeneracy between astrophysical uncertainties and particle physics unknowns is not complete. Certain models can be distinguished in a halo-independent way using a single ton-scale experiment based on liquid xenon, while other models are indistinguishable with a single experiment but can be separated using combined information from several target elements.

A reliable comparison of different dark matter (DM) searches requires models that satisfy certain consistency requirements like gauge invariance and perturbative unitarity. As a well-motivated example, we study two-mediator DM (2MDM). The model is based on a spontaneously broken $U(1)'$ gauge symmetry and contains a Majorana DM particle as well as two $s$-channel mediators, one vector (the $Z'$) and one scalar (the dark Higgs). We perform a global scan over the parameters of the model assuming that the DM relic density is obtained by thermal freeze-out in the early Universe and imposing a large set of constraints: direct and indirect DM searches, monojet, dijet and dilepton searches at colliders, Higgs observables, electroweak precision tests and perturbative unitarity. We conclude that thermal DM is only allowed either close to an $s$-channel resonance or if at least one mediator is lighter than the DM particle. In these cases a thermal DM abundance can be obtained although DM couplings to the Standard Model are tiny. Interestingly, we find that vector-mediated DM-nucleon scattering leads to relevant constraints despite the velocity-suppressed cross section, and that indirect detection can be important if DM annihilations into both mediators are kinematically allowed.

We show that simplified models used to describe the interactions of dark matter with Standard Model particles do not in general respect gauge invariance and that perturbative unitarity may be violated in large regions of the parameter space. The modifications necessary to cure these inconsistencies may imply a much richer phenomenology and lead to stringent constraints on the model. We illustrate these observations by considering the simplified model of a fermionic dark matter particle and a vector mediator. Imposing gauge invariance then leads to strong constraints from dilepton resonance searches and electroweak precision tests. Furthermore, the new states required to restore perturbative unitarity can mix with Standard Model states and mediate interactions between the dark and the visible sector, leading to new experimental signatures such as invisible Higgs decays. The resulting constraints are typically stronger than the 'classic' constraints on DM simplified models such as monojet searches and make it difficult to avoid thermal overproduction of dark matter.

A gauge singlet scalar with non-minimal coupling to gravity can drive inflation and later freeze out to become cold dark matter. We explore this idea by revisiting inflation in the singlet direction (S-inflation) and Higgs Portal Dark Matter in light of the Higgs discovery, limits from LUX and observations by Planck. We show that large regions of parameter space remain viable, so that successful inflation is possible and the dark matter relic abundance can be reproduced. Moreover, the scalar singlet can stabilise the electroweak vacuum and at the same time overcome the problem of unitarity-violation during inflation encountered by Higgs Inflation, provided the singlet is a real scalar. The 2-$\sigma$ Planck upper bound on $n_s$ imposes that the singlet mass is below 2 TeV, so that almost the entire allowed parameter range can be probed by XENON1T.

Self-interactions of dark matter particles can potentially lead to an observable separation between the dark matter halo and the stars of a galaxy moving through a region of large dark matter density. Such a separation has recently been observed in a galaxy falling into the core of the galaxy cluster Abell 3827. We estimated the DM self-interaction cross section needed to reproduce the observed effects and find that the sensitivity of Abell 3827 has been significantly overestimated in a previous study. Our corrected estimate is $\tilde{\sigma}/m_\text{DM} \sim 3\:\text{cm}^2\:\text{g}^{-1}$ when self-interactions result in an effective drag force and $\sigma/m_\text{DM} \sim 1.5\:\text{cm}^2\:\text{g}^{-1}$ for the case of contact interactions, in some tension with previous upper bounds.

Uncertainty in the local dark matter velocity distribution is a key difficulty in the analysis of data from direct detection experiments. Here we propose a new approach for dealing with this uncertainty, which does not involve any assumptions about the structure of the dark matter halo. Given a dark matter model, our method yields the velocity distribution which best describes a set of direct detection data as a finite sum of streams with optimised speeds and densities. The method is conceptually simple and numerically very efficient. We give an explicit example in which the method is applied to determining the ratio of proton to neutron couplings of dark matter from a hypothetical set of future data.

When a dark matter halo moves through a background of dark matter particles, self-interactions can lead to both deceleration and evaporation of the halo and thus shift its centroid relative to the collisionless stars and galaxies. We study the magnitude and time evolution of this shift for two classes of dark matter self-interactions, viz. frequent self-interactions with small momentum transfer (e.g. due to long-range interactions) and rare self-interactions with large momentum transfer (e.g. contact interactions), and find important differences between the two cases. We find that neither effect can be strong enough to completely separate the dark matter halo from the galaxies, if we impose conservative bounds on the self-interaction cross-section. The majority of both populations remain bound to the same gravitational potential and the peaks of their distributions are therefore always coincident. Consequently any apparent separation is mainly due to particles which are leaving the gravitational potential, so will be largest shortly after the collision but not observable in evolved systems. Nevertheless the fraction of collisions with large momentum transfer is an important characteristic of self-interactions, which can potentially be extracted from observational data and provide an important clue as to the nature of dark matter.

The CDMS-II collaboration has reported 3 events in a Si detector, which are consistent with being nuclear recoils due to scattering of Galactic dark matter particles with a mass of about 8.6 GeV and a cross-section on neutrons of about 2 x 10^-41 cm^2. While a previous result from the XENON10 experiment has supposedly ruled out such particles as dark matter, we find by reanalysing the XENON10 data that this is not the case. Some tension remains however with the upper limit placed by the XENON100 experiment, independently of astrophysical uncertainties concerning the Galactic dark matter distribution. We explore possible ways of ameliorating this tension by altering the properties of dark matter interactions. Nevertheless, even with standard couplings, light dark matter is consistent with both CDMS and XENON10/100.

The latest results from LHC searches for jets in association with missing transverse energy place strong bounds on the scattering cross section of dark matter. For the case of spin-dependent or momentum suppressed interactions these limits seem to be superior to the bounds from direct detection experiments. In this article, we show that loop contributions can significantly alter this conclusion and boost direct detection bounds, whenever they induce spin-independent interactions. This effect is most striking for tensor and pseudotensor interactions, which induce magnetic and electric dipole moments at loop level. For axialvector and anapole interactions a relevant contribution to direct detection signals arises from loop-induced Yukawa-like couplings between dark matter and quarks. We furthermore compare the resulting bounds to additional constraints on these effective operators arising from indirect searches and relic density requirements.

Improved limits as well as tentative claims for dark matter annihilation into gamma-ray lines have been presented recently. We study the direct detection cross section induced from dark matter annihilation into two photons in a model-independent fashion, assuming no additional couplings between dark matter and nuclei. We find a striking non-standard recoil spectrum due to different destructively interfering contributions to the dark matter nucleus scattering cross section. While in the case of s-wave annihilation the current sensitivity of direct detection experiments is insufficient to compete with indirect detection searches, for p-wave annihilation the constraints from direct searches are comparable. This will allow to test dark matter scenarios with p-wave annihilation that predict a large di-photon annihilation cross section in the next generation of experiments.

We study the impact of the assumed velocity distribution of galactic dark matter particles on the interpretation of results from nuclear recoil detectors. By converting experimental data to variables that make the astrophysical unknowns explicit, different experiments can be compared without implicit assumptions concerning the dark matter halo. We extend this framework to include the annual modulation signal, as well as multiple target elements. Recent results from DAMA, CoGeNT and CRESST-II can be brought into agreement if the velocity distribution is very anisotropic and thus allows a large modulation fraction. However constraints from CDMS and XENON cannot be evaded by appealing to such astrophysical uncertainties alone.

DAMA observes an annual modulation in their event rate, as might be expected from dark matter scatterings, while CoGeNT has reported evidence for a similar modulation. The simplest interpretation of these findings in terms of dark matter-nucleus scatterings is excluded by other direct detection experiments. We consider the robustness of these exclusions with respect to assumptions regarding the scattering and find that isospin-violating inelastic dark matter helps alleviate this tension and allows marginal compatibility between experiments. Isospin-violation can significantly weaken the XENON constraints, while inelasticity enhances the annual modulation fraction of the signal, bringing the CoGeNT and CDMS results into better agreement.

We provide a new way of constraining the relative scintillation efficiency L_eff for liquid xenon. Using a simple estimate for the electronic and nuclear stopping powers together with an analysis of recombination processes we predict both the ionization and the scintillation yields. Using presently available data for the ionization yield, we can use the correlation between these two quantities to constrain L_eff from below. Moreover, we argue that more reliable data on the ionization yield would allow to verify our assumptions on the atomic cross sections and to predict the value of L_eff. We conclude that the relative scintillation efficiency should not decrease at low nuclear recoil energies, which has important consequences for the robustness of exclusion limits for low WIMP masses in liquid xenon Dark Matter searches.