We perform a global fit within the pseudo-Nambu-Goldstone dark matter (DM) model emerging from an additional complex scalar singlet with a softly broken global U(1) symmetry. Leading to a momentum-suppressed DM-nucleon cross section at tree level, the model provides a natural explanation for the null results from direct detection experiments. Our global fit combines constraints from perturbative unitarity, DM relic abundance, Higgs invisible decay, electroweak precision observables and latest Higgs searches at colliders. The results are presented in both frequentist and Bayesian statistical frameworks. Furthermore, post-processing our samples, we include the likelihood from gamma-ray observations of $\mathit{Fermi}$-LAT dwarf spheroidal galaxies and compute the one-loop DM-nucleon cross section. We find two favoured regions characterised by their dominant annihilation channel: the Higgs funnel and annihilation into Higgs pairs. Both are compatible with current $\mathit{Fermi}$-LAT observations, and furthermore, can fit the slight excess observed in four dwarfs in a mass range between about 30-300 GeV. While the former region is hard to probe experimentally, the latter can partly be tested by current observations of cosmic-ray antiprotons as well as future gamma-ray observations.

The current paradigm for dark matter direct detection is to assume that the dark sector is solely composed of a single particle species. In this short paper, we make the observation that dark matter comprising both a light and a heavy component that modulate out of phase leads to interesting phenomenology in annual modulation experiments. For an illustrative example, we use the recently released DAMA/LIBRA phase-2 results with a lower energy threshold. Immediately after, it was argued that a one-component spin-independent dark matter explanation of the observed annual modulation is strongly disfavored or excluded unless isospin-violating couplings are invoked. We show that a simple two-component extension can reproduce the observed spectrum without the need to invoke fine-tuned couplings. Using the publicly available DAMA/LIBRA data, we perform a fit of the DAMA/LIBRA energy spectrum of the annual modulation amplitude to a scenario with two dark matter components. We also take into account how gravitational focusing affects the phases of the light and a heavy components differently, which leads to nontrivial effects in the total time-dependent rate. Our results show that there exists a unique solution in agreement with the data in the simplest case of isospin-conserving couplings with equal cross sections. The distinctive features found in this work are crucial for a dark matter interpretation of any observed annual modulation.

We study the case of multi-component dark matter, in particular how direct detection signals are modified in the presence of several stable weakly-interacting-massive particles. Assuming a positive signal in a future direct detection experiment, stemming from two dark matter components, we study the region in parameter space where it is possible to distinguish a one from a two-component dark matter spectrum. First, we leave as free parameters the two dark matter masses and show that the two hypotheses can be significantly discriminated for a range of dark matter masses with their splitting being the critical factor. We then investigate how including the effects of different interaction strengths, local densities or velocity dispersions for the two components modifies these conclusions. We also consider the case of isospin-violating couplings. In all scenarios, we show results for various types of nuclei both for elastic spin-independent and spin-dependent interactions. Finally, assuming that the two-component hypothesis is confirmed, we quantify the accuracy with which the parameters can be extracted and discuss the different degeneracies that occur. This includes studying the case in which only a single experiment observes a signal, and also the scenario of having two signals from two different experiments, in which case the ratios of the couplings to neutrons and protons may also be extracted.

Dark matter searches in gamma ray final states often make use of the fact that photons can be produced from final state muons. Modern Monte Carlo generators and DM codes include the effects of final state radiation from muons produced in the dark matter annihilation process itself, but neglect the O(1%) radiative correction that arises from the subsequent muon decay. After implementing this correction we demonstrate the effect that it can have on dark matter phenomenology by considering the case of dark matter annihilation to four muons via scalar mediator production. We first show that the AMS-02 positron excess can no longer easily be made consistent with this final state once the Fermi-LAT dwarf limits are calculated with the inclusion of radiative muon decays, and we next show that the Fermi-LAT galactic centre gamma excess can be improved with this final state after inclusion of the same effect. We provide code and tables for the implementation of this effect in the popular dark matter code micrOMEGAs, providing a solution for any model producing final state muons.

We present LUXCalc, a new utility for calculating likelihoods and deriving WIMP-nucleon coupling limits from the recent results of the LUX direct search dark matter experiment. After a brief review of WIMP-nucleon scattering, we derive LUX limits on the spin-dependent WIMP-nucleon couplings over a broad range of WIMP masses, under standard assumptions on the relevant astrophysical parameters. We find that, under these and other common assumptions, LUX excludes the entire spin-dependent parameter space consistent with a dark matter interpretation of DAMA's anomalous signal, the first time a single experiment has been able to do so. We also revisit the case of spin-independent couplings, and demonstrate good agreement between our results and the published LUX results. Finally, we derive constraints on the parameters of an effective dark matter theory in which a spin-1 mediator interacts with a fermionic WIMP and Standard Model fermions via axial-vector couplings. A detailed appendix describes the use of LUXCalc with standard codes to place constraints on generic dark matter theories.