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.

Thermal dark matter at the MeV scale faces stringent bounds from a variety of cosmological probes. Here we perform a detailed evaluation of BBN bounds on the annihilation cross section of dark matter with a mass $1\,\text{MeV} \lesssim m_\chi \lesssim 1\,\text{GeV}$. For $p-wave suppressed annihilations, constraints from BBN turn out to be significantly stronger than the ones from CMB observations, and are competitive with the strongest bounds from other indirect searches. We furthermore update the lower bound from BBN on the mass of thermal dark matter using improved determinations of primordial abundances. While being of similar strength as the corresponding bound from CMB, it is significantly more robust to changes in the particle physics model.

Meta-stable dark sector particles decaying into electrons or photons may non-trivially change the Hubble rate, lead to entropy injection into the thermal bath of Standard Model particles and may also photodisintegrate light nuclei formed in the early universe. We study generic constraints from Big Bang Nucleosynthesis on such a setup, with a particular emphasis on MeV-scale particles which are neither fully relativistic nor non-relativistic during all times relevant for Big Bang Nucleosynthesis. We apply our results to a simple model of self-interacting dark matter with a light scalar mediator. This setup turns out to be severely constrained by these considerations in combination with direct dark matter searches and will be fully tested with the next generation of low-threshold direct detection experiments.

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$.

Light vector mediators can naturally induce velocity-dependent dark matter self-interactions while at the same time allowing for the correct dark matter relic abundance via thermal freeze-out. If these mediators subsequently decay into Standard Model states such as electrons or photons however, this is robustly excluded by constraints from the Cosmic Microwave Background. We study to what extent this conclusion can be circumvented if the vector mediator is stable and hence contributes to the dark matter density while annihilating into lighter degrees of freedom. We find viable parts of parameter space which lead to the desired self-interaction cross section of dark matter to address the small-scale problems of the collisionless cold dark matter paradigm while being compatible with bounds from the Cosmic Microwave Background and Big Bang Nucleosynthesis observations.

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 study constraints from Big Bang Nucleosynthesis on inert particles in a dark sector which contribute to the Hubble rate and therefore change the predictions of the primordial nuclear abundances. We pay special attention to the case of MeV-scale particles decaying into dark radiation, which are neither fully relativistic nor non-relativistic during all temperatures relevant to Big Bang Nucleosynthesis. As an application we discuss the implications of our general results for models of self-interacting dark matter with light mediators.

The theoretical interpretation of dark matter direct detection experiments is hindered by uncertainties of the microphysics governing the dark matter-nucleon interaction, and of the dark matter density and velocity distribution inside the Solar System. These uncertainties are especially relevant when confronting a detection claim to the null results from other experiments, since seemingly conflicting experimental results may be reconciled when relaxing the assumptions about the form of the interaction and/or the velocity distribution. We present in this paper a halo-independent method to calculate the maximum number of events in a direct detection experiment given a set of null search results, allowing for the first time the scattering to be mediated by an arbitrary combination of various interactions (concretely we consider up to 64). We illustrate this method to examine the compatibility of the dark matter interpretation of the three events detected by the silicon detectors in the CDMS-II experiment with the null results from XENON1T and PICO-60.

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.

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.

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

Antibaryons are produced in our Galaxy in collisions of high energy cosmic rays with the interstellar medium and in old supernova remnants, and possibly, in exotic sources such as primordial black hole evaporation or dark matter annihilations and decays. The search for signals from exotic sources in antiproton data is hampered by large backgrounds from spallation which, within theoretical errors, can solely account for the current data. Due to the higher energy threshold for antideuteron production, which translates into a suppression of the low energy flux from spallations, antideuteron searches have been proposed as a probe for exotic sources. We perform in this paper a comprehensive analysis of the antideuteron fluxes at the Earth expected from known and hypothetical sources in our Galaxy, and we calculate their maximal values consistent with current antiproton data from AMS-02. We find that supernova remnants generate a negligible flux, whereas primordial black hole evaporation and dark matter annihilations or decays may dominate the total flux at low energies. On the other hand, we find that the detection of cosmic antideuterons would require, for the scenarios studied in this paper and assuming optimistic values of the coalescence momentum and solar modulation, an increase of the experimental sensitivity compared to ongoing and planned instruments by at least a factor of 2. Finally, we briefly comment on the prospects for antihelium-3 detection.

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.

We examine the dark matter interpretation of the modulation signal reported by the DAMA experiment from the perspective of effective field theories displaying Galilean invariance. We consider the most general effective coupling leading to the elastic scattering of a dark matter particle with spin 0 or 1/2 off a nucleon, and we analyze the compatibility of the DAMA signal with the null results from other direct detection experiments, as well as with the non-observation of a high energy neutrino flux in the direction of the Sun from dark matter annihilation. To this end, we develop a novel semi-analytical approach for comparing experimental results in the high-dimensional parameter space of the non-relativistic effective theory. Assuming the standard halo model, we find a strong tension between the dark matter interpretation of the DAMA modulation signal and the null result experiments. We also list possible ways-out of this conclusion.

Recent years have seen increased theoretical and experimental effort towards the first-ever detection of cosmic-ray antideuterons, in particular as an indirect signature of dark matter annihilation or decay. In contrast to indirect dark matter searches using positrons, antiprotons, or gamma-rays, which suffer from relatively high and uncertain astrophysical backgrounds, searches with antideuterons benefit from very suppressed conventional backgrounds, offering a potential breakthrough in unexplored phase space for dark matter. This article is based on the first dedicated cosmic-ray antideuteron workshop, which was held at UCLA in June 2014. It reviews broad classes of dark matter candidates that result in detectable cosmic-ray antideuteron fluxes, as well as the status and prospects of current experimental searches. The coalescence model of antideuteron production and the influence of antideuteron measurements at particle colliders are discussed. This is followed by a review of the modeling of antideuteron propagation through the magnetic fields, plasma currents, and molecular material of our Galaxy, the solar system, the Earth's geomagnetic field, and the atmosphere. Finally, the three ongoing or planned experiments that are sensitive to cosmic-ray antideuterons, BESS, AMS-02, and GAPS, are detailed. As cosmic-ray antideuteron detection is a rare event search, multiple experiments with orthogonal techniques and backgrounds are essential. Many theoretical and experimental groups have contributed to these studies over the last decade, this review aims to provide the first coherent discussion of the relevant dark matter theories that antideuterons probe, the challenges to predictions and interpretations of antideuteron signals, and the experimental efforts toward cosmic antideuteron detection.

We present a comprehensive study of a model where the dark matter is composed of a singlet real scalar that couples to the Standard Model predominantly via a Yukawa interaction with a light quark and a colored vector-like fermion. A distinctive feature of this scenario is that thermal freeze-out in the early universe may be driven by annihilation both into gluon pairs at one-loop ($gg$) and by virtual internal Bremsstrahlung of a gluon ($q \bar{q} g$). Such a dark matter candidate may also be tested through direct and indirect detection and at the LHC; viable candidates have either a mass nearly degenerate with that of the fermionic mediator or a mass above about 2 TeV.

We propose a method that allows to place an upper limit on the dark matter elastic scattering cross section with nucleons which is independent of the velocity distribution. Our approach combines null results from direct detection experiments with indirect searches at neutrino telescopes, and goes beyond previous attempts to remove astrophysical uncertainties in that it directly constrains the particle physics properties of the dark matter. The resulting halo-independent upper limits on the scattering cross section of dark matter are remarkably strong and reach $\sigma_{\text{SI}}^{p} \lesssim 10^{-43} \, (10^{-42})~\text{cm}^2$ and $\sigma_{\text{SD}}^{p} \lesssim 10^{-37} \, (3\times 10^{-37 })~\text{cm}^2$, for dark matter particles of $m_{\text{DM}}\sim 1~\text{TeV}$ annihilating into $W^+W^-$ ($b\bar b$), assuming $\rho_\text{loc}=0.3\text{ GeV}/\text{cm}^3$.

The precise measurement of cosmic-ray antiparticles serves as important means for identifying the nature of dark matter. Recent years showed that identifying the nature of dark matter with cosmic-ray positrons and higher energy antiprotons is difficult, and has lead to a significantly increased interest in cosmic-ray antideuteron searches. Antideuterons may also be generated in dark matter annihilations or decays, offering a potential breakthrough in unexplored phase space for dark matter. Low-energy antideuterons are an important approach because the flux from dark matter interactions exceeds the background flux by more than two orders of magnitude in the low-energy range for a wide variety of models. This review is based on the "dbar14 - dedicated cosmic-ray antideuteron workshop", which brought together theorists and experimentalists in the field to discuss the current status, perspectives, and challenges for cosmic-ray antideuteron searches and discusses the motivation for antideuteron searches, the theoretical and experimental uncertainties of antideuteron production and propagation in our Galaxy, as well as give an experimental cosmic-ray antideuteron search status update. This report is a condensed summary of the article "Review of the theoretical and experimental status of dark matter identification with cosmic-ray antideuteron" (arXiv:1505.07785).

We analyze the direct detection signals of a toy model consisting of a Dirac dark matter particle which couples to one Standard Model fermion via a scalar mediator. For all scenarios, the dark matter particle scatters off nucleons via one loop-induced electromagnetic and electroweak moments, as well as via the one-loop exchange of a Higgs boson. Besides, and depending on the details of the model, the scattering can also be mediated at tree level via the exchange of the scalar mediator or at one loop via gluon-gluon interactions. We show that, for thermally produced dark matter particles, the current limits from the LUX experiment on these scenarios are remarkably strong, even for dark matter coupling only to leptons. We also discuss future prospects for XENON1T and DARWIN and we argue that multi-ton xenon detectors will be able to probe practically the whole parameter space of the model consistent with thermal production and perturbativity. We also discuss briefly the implications of our results for the dark matter interpretation of the Galactic GeV excess.

The search for sharp features in the gamma-ray spectrum is a promising approach to identify a signal from dark matter annihilations over the astrophysical backgrounds. In this paper we investigate the generation of gamma-ray lines and internal bremsstrahlung signals in a toy model where the dark matter particle is a real scalar that couples to a lepton and an exotic fermion via a Yukawa coupling. We show that the Fermi-LAT and H.E.S.S. searches for line-like spectral features severely constrain regions of the parameter space where the scalar dark matter is thermally produced. Finally, we also discuss the complementarity of the searches for sharp spectral features with other indirect dark matter searches, as well as with direct and collider searches.

Cosmic-ray anti-nuclei provide a promising discovery channel for the indirect detection of particle dark matter. Hadron showers produced by the pair-annihilation or decay of Galactic dark matter generate anti-nucleons which can in turn form light anti-nuclei. Previous studies have only focused on the spectrum and flux of low energy antideuterons which, although very rarely, are occasionally also produced by cosmic-ray spallation. Heavier elements ($A\geq3$) have instead entirely negligible astrophysical background and a primary yield from dark matter which could be detectable by future experiments. Using a Monte Carlo event generator and an event-by-event phase space analysis, we compute, for the first time, the production spectrum of \antihe and \antiT for dark matter annihilating or decaying to $b\bar{b}$ and ${W^+}{W^-}$ final states. We then employ a semi-analytic model of interstellar and heliospheric propagation to calculate the \antihe flux as well as to provide tools to relate the anti-helium spectrum corresponding to an arbitrary antideuteron spectrum. Finally, we discuss prospects for current and future experiments, including GAPS and AMS-02.

Dark matter particles captured in the Sun would annihilate producing a neutrino flux that could be detected at the Earth. In some channels, however, the neutrino flux lies in the MeV range and is thus undetectable at IceCube, namely when the dark matter particles annihilate into $e^+e^-$, $\mu^+\mu^-$ or light quarks. On the other hand, the same interaction that mediates the annihilations into light fermions also leads, via higher order effects, to the production of weak gauge bosons (and in the case of quarks also gluons) that generate a high energy neutrino flux potentially observable at IceCube. We consider in this paper tree level annihilations into a fermion-antifermion pair with the associated emission of one gauge boson and one loop annihilations into two gauge bosons, and we calculate the limits on the scattering cross section of dark matter particles with protons in scenarios where the dark matter particle couples to electrons, muons or light quarks from the non-observation of an excess of neutrino events in the direction of the Sun. We find that the limits on the spin-dependent scattering cross section are, for some scenarios, stronger than the limits from direct detection experiments.

We investigate the prospects to observe a high energy neutrino signal from dark matter annihilations in the Sun in scenarios where the dark matter is a Majorana fermion that couples to a quark and a colored scalar via a Yukawa coupling. In this minimal scenario, the dark matter capture and annihilation in the Sun can be studied in a single framework. We find that, for small and moderate mass splitting between the dark matter and the colored scalar, the two-to-three annihilation $q \bar q g$ plays a central role in the calculation of the number of captured dark matter particles. On the other hand, the two-to-three annihilation into $q \bar q Z$ gives, despite its small branching fraction, the largest contribution to the neutrino flux at the Earth at the highest energies. We calculate the limits on the model parameters using IceCube observations of the Sun and we discuss their interplay with the requirement of equilibrium of captures and annihilations in the Sun and with the requirement of thermal dark matter production. We also compare the limits from IceCube to the limits from direct detection, antiproton measurements and collider searches.

The search for cosmic antideuterons has been proposed as a promising method to indirectly detect dark matter, due to the very small background flux from spallations expected at the energies relevant to experiments. The antideuteron flux from dark matter annihilations or decays is, however, severely constrained by the non-observation of an excess in the antiproton-to-proton fraction measured by PAMELA. In this paper we calculate, for representative dark matter annihilation and decay channels, upper limits on the number of antideuteron events at AMS-02 and GAPS from requiring that the associated antiproton flux is in agreement with the PAMELA data. To this end, we first analyze in detail the formation of antideuterons in the coalescence model using an event-by-event Monte Carlo simulation and using data from various high energy experiments. We find that the resulting coalescence momentum shows a dependence on the underlying process and on the center of mass energy involved. Then, we calculate, using a diffusion model, the flux of antideuterons at the Earth from dark matter annihilations or decays. Our results indicate that, despite the various sources of uncertainty, the observation of an antideuteron flux at AMS-02 or GAPS from dark matter annihilations or decays will be challenging.

We investigate in this paper the antideuteron flux produced in high energy collisions of cosmic rays with the interstellar matter. We employ the Monte Carlo generator DPMJET-III together with the coalescence model to simulate, in an event-by-event basis, the antideuteron production in cosmic ray collisions in our Galaxy. Then, we use a diffusion model to calculate the expected flux. We find a secondary antideuteron flux at the Earth which is approximately a factor of two smaller than in previous calculations.