
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 highenergy neutrinos from
dark matter annihilation in the Sun with IceCube, and searches for gamma rays
from annihilation in dwarf galaxies with the FermiLAT. 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 lowmass 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 finetuned.

Dark matter in the cosmological concordance model is parameterised by a
single number, describing the covariantly conserved energy density of a
nonrelativistic fluid. Here we test this assumption in a modelindependent and
conservative way by considering the possibility that, at any point during the
cosmological evolution, dark matter may be converted into a noninteracting
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 largescale 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 Sommerfeldenhanced dark
matter annihilation. We identify a small part of parameter space that both
mitigates the tension between cosmic microwave and largescale structure data
and allows for velocitydependent dark matter selfinteractions strong enough
to address the smallscale 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 modelindependent 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 nonthermalisation. We find that
even the HLLHC cannot possibly produce more than a few events with energy
above 10 GeV involving states from a nonthermalised 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 gammaray
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 spinindependent and
spindependent 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 predominantly with Standard Model gauge
bosons (socalled axionlike 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 longlived 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 freezeout, leading to a highly predictive scenario consistent
with all existing constraints but testable with singlephoton searches at Belle
II in the near future.

Selfinteracting dark matter may have striking astrophysical signatures, such
as observable offsets between galaxies and dark matter in merging galaxy
clusters. Numerical Nbody 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
selfscatterings: deceleration and evaporation. We point out that
selfscatterings 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 selfscattering cross section, and
cannot be predicted from the momentum transfer cross section alone.

Dark matter particles interacting via the exchange of very light spin0
mediators can have large selfinteraction rates and obtain their relic
abundance from thermal freezeout. 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
selfinteractions in astrophysical systems. For the case of a mediator with
purely scalar couplings we point out the highly relevant impact of
lowthreshold direct detection experiments like CRESSTII, which essentially
rule out the simplest realization of this model. These constraints can be
significantly relaxed if the mediator has CPviolating 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 CRESSTIII, SuperCDMS SNOLAB and EDELWEISSIII 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 selfinteractions 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 spin0 simplified models
employed in analyses of Run I data these models are selfconsistent, 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
nontrivial constraints. We emphasise in this context the sensitivity of the
$t\bar t + E_{T, \rm miss}$, mono$Z$ and monoHiggs channels, which yield
stronger limits than monojet searches in large parts of the parameter space.
Constraints from spin0 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 largeradius 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 monojet or
dijet searches.

Coupling dark matter to light new particles is an attractive way to combine
thermal production with strong velocitydependent selfinteractions. 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 swave annihilation these constraints exclude the
entire parameter space where the selfinteractions are large enough to address
the smallscale problems of structure formation.

Weaklycoupled TeVscale 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 "monoX" 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 axialvector mediators, based on a workshop that took
place on September 19/20, 2016 and subsequent discussions. These suggestions
include how to extend the spin1 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 nonrelativistic effective
interactions with stateoftheart haloindependent methods. This approach
makes it possible to analyse direct detection experiments for arbitrary dark
matter interactions and quantify the goodnessoffit 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
haloindependent way using a single tonscale 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 wellmotivated example, we study twomediator 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 freezeout 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
vectormediated DMnucleon scattering leads to relevant constraints despite the
velocitysuppressed cross section, and that indirect detection can be important
if DM annihilations into both mediators are kinematically allowed.

The potential sensitivity to isospinbreaking effects makes LHC searches for
monoW 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 monoW channel to the relative magnitude and sign of the
uptype and downtype quark couplings to dark matter is an artefact of
unitarity violation. We provide three different solutions to this monoW
problem in the context of spin1 simplified models and briefly discuss the
impact that our findings have on the prospects of monoW 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 schannel and tchannel scenarios.
For schannel, spin0 and spin1 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 nearfuture 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 socalled axionlike 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 nonminimal 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 (Sinflation) 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 unitarityviolation 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.

Selfinteractions 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 selfinteraction 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
selfinteractions 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.