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.

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 all conversion times 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 show under rather general assumptions that hidden sectors that never reach thermal equilibrium in the early Universe are also inaccessible for the LHC. In other words, any particle that can be produced at the LHC must either have been in thermal equilibrium with the Standard Model at some point or must be produced via the decays of another hidden sector particle that has been in thermal equilibrium. To reach this conclusion, we parametrise the cross section connecting the Standard Model to the hidden sector in a very general way and use methods from linear programming to calculate the largest possible number of LHC events compatible with the requirement of non-thermalisation. We find that even the HL-LHC cannot possibly produce more than a few events with energy above 10 GeV involving states from a non-thermalised hidden sector.

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.

Light pseudoscalars interacting pre-dominantly with Standard Model gauge bosons (so-called axion-like particles or ALPs) occur frequently in extensions of the Standard Model. In this work we review and update existing constraints on ALPs in the keV to GeV mass region from colliders, beam dump experiments and astrophysics. We furthermore provide a detailed calculation of the expected sensitivity of Belle II, which can search for visibly and invisibly decaying ALPs, as well as long-lived ALPs. The Belle II sensitivity is found to be substantially better than previously estimated, covering wide ranges of relevant parameter space. In particular, Belle II can explore an interesting class of dark matter models, in which ALPs mediate the interactions between the Standard Model and dark matter. In these models, the relic abundance can be set via resonant freeze-out, leading to a highly predictive scenario consistent with all existing constraints but testable with single-photon searches at Belle II in the near future.

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.

We study a new class of renormalisable simplified models for dark matter searches at the LHC that are based on two Higgs doublet models with an additional pseudoscalar mediator. In contrast to the spin-0 simplified models employed in analyses of Run I data these models are self-consistent, unitary and bounds from Higgs physics typically pose no constraints. Predictions for various missing transverse energy ($E_{T, \rm miss}$) searches are discussed and the reach of the 13 TeV LHC is explored. It is found that the proposed models provide a rich spectrum of complementary observables that lead to non-trivial constraints. We emphasise in this context the sensitivity of the $t\bar t + E_{T, \rm miss}$, mono-$Z$ and mono-Higgs channels, which yield stronger limits than mono-jet searches in large parts of the parameter space. Constraints from spin-0 resonance searches, electroweak precision measurements and flavour observables are also derived and shown to provide further important handles to constraint and to test the considered dark matter models.

We discuss a novel signature of dark matter production at the LHC resulting from the emission of an additional Higgs boson in the dark sector. The presence of such a dark Higgs boson is motivated simultaneously by the need to generate the masses of the particles in the dark sector and the possibility to relax constraints from the dark matter relic abundance by opening up a new annihilation channel. If the dark Higgs boson decays into Standard Model states via a small mixing with the Standard Model Higgs boson, one obtains characteristic large-radius jets in association with missing transverse momentum that can be used to efficiently discriminate signal from backgrounds. We present the sensitivities achievable in LHC searches for dark Higgs bosons with already collected data and demonstrate that such searches can probe large regions of parameter space that are inaccessible to conventional mono-jet or di-jet searches.

Coupling dark matter to light new particles is an attractive way to combine thermal production with strong velocity-dependent self-interactions. Here we point out that in such models the dark matter annihilation rate is generically enhanced by the Sommerfeld effect, and we derive the resulting constraints from the Cosmic Microwave Background and other indirect detection probes. For the frequently studied case of s-wave annihilation these constraints exclude the entire parameter space where the self-interactions are large enough to address the small-scale problems of structure formation.

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.

This review discusses both experimental and theoretical aspects of searches for dark matter at the LHC. An overview of the various experimental search channels is given, followed by a summary of the different theoretical approaches for predicting dark matter signals. A special emphasis is placed on the interplay between LHC dark matter searches and other kinds of dark matter experiments, as well as among different types of LHC searches.

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 analyse a combination of ATLAS and CMS searches for dijet resonances at run I and run II, presenting the results in a way that can be easily applied to a generic Z' model. As an illustrative example, we consider a simple model of a Z' coupling to quarks and dark matter. We first study a benchmark case with fixed couplings and then focus on the assumption that the Z' is responsible for setting the dark matter relic abundance. Dijet constraints place significant bounds on this scenario, allowing us to narrow down the allowed range of dark matter masses for given Z' mass and width.

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.

The potential sensitivity to isospin-breaking effects makes LHC searches for mono-W signatures promising probes of the coupling structure between the Standard Model and dark matter. It has been shown, however, that the strong sensitivity of the mono-W channel to the relative magnitude and sign of the up-type and down-type quark couplings to dark matter is an artefact of unitarity violation. We provide three different solutions to this mono-W problem in the context of spin-1 simplified models and briefly discuss the impact that our findings have on the prospects of mono-W searches at future LHC runs.

This document outlines a set of simplified models for dark matter and its interactions with Standard Model particles. It is intended to summarize the main characteristics that these simplified models have when applied to dark matter searches at the LHC, and to provide a number of useful expressions for reference. The list of models includes both s-channel and t-channel scenarios. For s-channel, spin-0 and spin-1 mediation is discussed, and also realizations where the Higgs particle provides a portal between the dark and visible sectors. The guiding principles underpinning the proposed simplified models are spelled out, and some suggestions for implementation are presented.

This document summarises the proposal of the LHC Dark Matter Working Group on how to present LHC results on $s$-channel simplified dark matter models and to compare them to direct (indirect) detection experiments.

We show that contrary to naive expectations the recently observed diphoton excess can be explained by a vector resonance, which decays to a photon and a light scalar $s$, followed by a decay of the scalar into two photons: $Z' \rightarrow \gamma s \rightarrow 3 \gamma$. As the two photons from the scalar decay are highly boosted, the experimental signature is an apparent diphoton final state. In fact all the necessary ingredients are naturally present in $Z'$ models: Additional fermions with electroweak quantum numbers are required in order to render the theory anomaly free and naturally induce the required effective couplings, while the hidden Higgs which gives mass to the $Z'$ can be very light. In particular no new coloured states are required in this framework. We also show that in such a setup the width of the resonance can be rather large, while all couplings remain perturbative.

With their high beam energy and intensity, existing and near-future proton beam dumps provide an excellent opportunity to search for new very weakly coupled particles in the MeV to GeV mass range. One particularly interesting example is a so-called axion-like particle (ALP), i.e. a pseudoscalar coupled to two photons. The challenge in proton beam dumps is to reliably calculate the production of the new particles from the interactions of two composite objects, the proton and the target atoms. In this work we argue that Primakoff production of ALPs proceeds in a momentum range where production rates and angular distributions can be determined to sufficient precision using simple electromagnetic form factors. Reanalysing past proton beam dump experiments for this production channel, we derive novel constraints on the parameter space for ALPs. We show that the NA62 experiment at CERN could probe unexplored parameter space by running in 'dump mode' for a few days and discuss opportunities for future experiments such as SHiP.

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.