Explore the vast range of titles available.

We introduce DRAKE, a numerical precision tool for predicting the dark matter relic abundance also in situations where the standard assumption of kinetic equilibrium during the freeze-out process may not be satisfied. DRAKE comes with a set of three dedicated Boltzmann equation solvers that implement, respectively, the traditionally adopted equation for the dark matter number density, fluid-like equations that couple the evolution of number density and velocity dispersion, and a full numerical evolution of the phase-space distribution. We review the general motivation for these approaches and, for illustration, highlight three concrete classes of models where kinetic and chemical decoupling are intertwined in a way that quantitatively impacts the relic density: (i) dark matter annihilation via a narrow resonance, (ii) Sommerfeld-enhanced annihilation and (iii) ‘forbidden’ annihilation to final states that are kinematically inaccessible at threshold. We discuss all these cases in some detail, demonstrating that the commonly adopted, traditional treatment can result in an estimate of the relic density that is wrong by up to an order of magnitude. The public release of DRAKE, along with several examples of how to calculate the relic density in concrete models, is provided at drake.hepforge.org

A bstract The irreducible upscattering of cold dark matter by cosmic rays opens up the intriguing possibility of detecting even light dark matter in conventional direct detection experiments or underground neutrino detectors. The mechanism also significantly enhances sensitivity to models with very large nuclear scattering rates, where the atmosphere and rock overburden efficiently stop standard non-relativistic dark matter particles before they could reach the detector. In this article, we demonstrate that cosmic-ray upscattering essentially closes the window for strongly interacting dark matter in the (sub-)GeV mass range. Arriving at this conclusion crucially requires a detailed treatment of both nuclear form factors and inelastic dark matter-nucleus scattering, as well as including the full momentum-transfer dependence of scattering amplitudes. We illustrate the latter point by considering three generic situations where such a momentum-dependence is particularly relevant, namely for interactions dominated by the exchange of light vector or scalar mediators, respectively, and for dark matter particles of finite size. As a final concrete example, we apply our analysis to a putative hexaquark state, which has been suggested as a viable baryonic dark matter candidate. Once again, we find that the updated constraints derived in this work close a significant part of otherwise unconstrained parameter space.

A bstract Traditional direct searches for dark matter, looking for nuclear recoils in deep underground detectors, are challenged by an almost complete loss of sensitivity for light dark matter particles. Consequently, there is a significant effort in the community to devise new methods and experiments to overcome these difficulties, constantly pushing the limits of the lowest dark matter mass that can be probed this way. From a model-building perspective, the scattering of sub-GeV dark matter on nucleons essentially must proceed via new light mediator particles, given that collider searches place extremely stringent bounds on contact-type interactions. Here we present an updated compilation of relevant limits for the case of a scalar mediator, including a new estimate of the near-future sensitivity of the NA62 experiment as well as a detailed evaluation of the model-specific limits from Big Bang nucleosynthesis. We also derive updated and more general limits on DM particles upscattered by cosmic rays, applicable to arbitrary energy- and momentum dependences of the scattering cross section. Finally we stress that dark matter self-interactions, when evaluated beyond the common s -wave approximation, place stringent limits independently of the dark matter production mechanism. These are, for the relevant parameter space, generically comparable to those that apply in the commonly studied freeze-out case. We conclude that the combination of existing (or expected) constraints from accelerators and astrophysics, combined with cosmological requirements, puts robust limits on the maximally possible nuclear scattering rate. In most regions of parameter space these are at least competitive with the best projected limits from currently planned direct detection experiments.

A bstract Relic density calculations of dark matter freezing out from the primordial plasma have reached a high level of sophistication, with several numerical tools readily available that match the observationally required accuracy. Dark matter production via the freeze-in mechanism, on the other hand, is sensitive to much higher temperatures than in the freeze-out case, implying both technical and computational difficulties when aiming for the same level of precision. We revisit the formulation of freeze-in production in a way that facilitates the inclusion of in-medium corrections like plasma effects and the spin statistics of relativistic quantum gases, as well as the temperature dependence of dark matter production rates induced by the electroweak and strong phase transitions, and we discuss in detail the additional complications arising in the presence of s -channel resonances. We illustrate our approach in the context of Higgs portal models, and provide the most accurate calculation to date of the freeze-in abundance of Scalar Singlet dark matter. We explore in particular the case of small reheating temperatures, for which the couplings implied by the freeze-in mechanism may be testable at the LHC. Together with this article we present a major update 6.3 of DarkSUSY with the added capability of performing general freeze-in calculations, including all complications mentioned above.

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework GAMBIT. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from Planck, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.

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 that 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 stand-alone 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.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a 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 ∼ 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.

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 .

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the constrained minimal supersymmetric standard model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

A bstract It is well known that the annihilation of Majorana dark matter into fermions is helicity suppressed. Here, we point out that the underlying mechanism is a subtle combination of two distinct effects, and we present a comprehensive analysis of how the suppression can be partially or fully lifted by the internal bremsstrahlung of an additional boson in the final state. As a concrete illustration, we compute analytically the full amplitudes and annihilation rates of supersymmetric neutralinos to final states that contain any combination of two standard model fermions, plus one electroweak gauge boson or one of the five physical Higgs bosons that appear in the minimal supersymmetric standard model. We classify the various ways in which these three-body rates can be large compared to the two-body rates, identifying cases that have not been pointed out before. In our analysis, we put special emphasis on how to avoid the double counting of identical kinematic situations that appear for two-body and three-body final states, in particular on how to correctly treat differential rates and the spectrum of the resulting stable particles that is relevant for indirect dark matter searches. We find that both the total annihilation rates and the yields can be significantly enhanced when taking into account the corrections computed here, in particular for models with somewhat small annihilation rates at tree-level which otherwise would not be testable with indirect dark matter searches. Even more importantly, however, we find that the resulting annihilation spectra of positrons, neutrinos, gamma-rays and antiprotons differ in general substantially from the model-independent spectra that are commonly adopted, for these final states, when constraining particle dark matter with indirect detection experiments.