Constraints on mediator-based dark matter and scalar dark energy models using s = 13 TeV pp collision data collected by the ATLAS detector

Constraints on selected mediator-based dark matter models and a scalar dark energy model using up to 37 fb−1 √ s = 13 TeV pp collision data collected by the ATLAS detector at the LHC during 2015–2016 are summarised in this paper. The results of experimental searches in a variety of final states are interpreted in terms of a set of spin1 and spin-0 single-mediator dark matter simplified models and a second set of models involving an extended Higgs sector plus an additional vector or pseudo-scalar mediator. The searches considered in this paper constrain spin-1 leptophobic and leptophilic mediators, spin-0 colour-neutral and colour-charged mediators and vector or pseudo-scalar mediators embedded in extended Higgs sector models. In this case, also √ s = 8 TeV pp collision data are used for the interpretation of the results. The results are also interpreted for the first time in terms of light scalar particles that could contribute to the accelerating expansion of the universe (dark energy).


JHEP05(2019)142
For most of the models in this paper, WIMPs are potentially pair-produced in pp collisions at the Large Hadron Collider (LHC) [74]. These particles, denoted by the symbol χ throughout this paper, are stable over cosmological scales and do not interact with the detector. To identify events with DM, additional particle(s), X = jet, γ, W, Z, h, (t)t, (b)b, need to be produced in association with DM in a pp collision, in order to tag the event and detect the recoiling WIMPs as missing transverse momentum (with magnitude E miss T ). If the DM candidates can be produced at the LHC via an s-channel exchange of a new particle, then this mediator could also decay back into SM final states: resonance searches can therefore also be used to constrain DM models. The interplay of resonance and X + E miss T searches depends on the specific model choice and is further outlined in this paper. In the models under study, some of which are new with respect to previous ATLAS publications, one or more new particles mediate the interaction of DM with the SM particles. The first category considers simplified models mediated by a vector, axial-vector, scalar or pseudo-scalar mediator. In the case of simplified vector and scalar mediators, different types of interactions are explored (baryon-charged, neutral-flavour-changing and coloured interactions). The second category considers less simplified models involving an extended Higgs sector plus an additional mediator, either a vector or a pseudo-scalar particle. The assumptions and choices of the models closely follow the work of the DM Forum/LHC DM Working Group [75][76][77][78]. Analyses focusing on signatures compatible with (unstable) longlived particles decaying in the detector volume are also not considered in this paper [79].
Results from particle physics experiments may be used to elucidate the microscopic nature of DE [80,81]. Hadron collider data considering X + E miss T final states are used to constrain Horndeski models of DE [82] in an effective field theory (EFT) framework [83].
This paper aims to provide an overview of the experimental programme of ATLAS searches [84] for mediator-based DM production performed to date using 13 TeV protonproton collisions delivered by the LHC in 2015 and 2016. The studies presented in this paper use public ATLAS results. Since no significant excess over the expected SM background was found in any of these analyses, the results are used to constrain the available phase space for DM models. Furthermore, DE models are also constrained using these analyses.
The paper is structured as follows. The DM and DE models evaluated in this paper are outlined in section 2, while the data and simulation samples are described in section 3. The data analyses for each different signature are briefly described in section 4, where the complementarity of different final states is also discussed. Finally, the dominant systematic uncertainties affecting the modelling of the signal samples are highlighted in section 5 and the results are presented in section 6, followed by the conclusions (section 7).

Theoretical framework
All DM results presented in this paper are interpreted in the framework of simplified DM models [75,78,[85][86][87][88], where a new particle (or set of particles) mediates the interaction of DM with the SM particles. These DM simplified models, which overcome some of the shortcomings of previous EFT-based DM models [87,[89][90][91][92][93], can be classified according to content and properties of the particles that mediate the interaction between DM and -2 -

JHEP05(2019)142
Short description Acronym Symbol J P Charge Signatures Results (Section 4) Section: Vector/axial-vector mediator  Table 1. Summary of the mediator-based simplified models considered in this paper, along with the associated model acronym (2 nd column, defined in the text) and mediator symbol (3 rd column) used throughout. The 4 th and 5 th columns indicate the quantum numbers of the mediator. The '×' indicates the cases where no other charge than the new mediator's interaction is present. The 6 th column indicates the final-state signatures (details in section 4) and the 7 th column gives the reference to the interpretation.

Baryon-number-charged interaction
The baryon-number-charged mediator simplified model [75,106] (VBC) considers a spin-1 vector mediator. It also assumes that the charge of the U(1) symmetry coincides with the baryon number and it is spontaneously broken by a baryonic Higgs scalar. The DM candidate in this model is a stable baryonic state and it is neutral under the SM gauge symmetry. While the model can provide an ISR signature through s-channel Z B -mediator production subsequently decaying into a pair of DM candidates as for the V/AV models described in the previous section (figure 1a), it can also exhibit a distinctive h + E miss T signature [106], as shown in figure 2a. The model has five parameters [106], whose values are chosen to enhance the cross-section for h + E miss T final states relative to traditional ISR signatures. The mixing angle between the baryonic and the SM Higgs bosons, θ, is fixed to sin θ = 0.3 in order to comply with the current Higgs boson coupling measurements. The coupling of the mediator Z B with the quarks, g q , and the DM, g χ , are set to 1/3 and 1, respectively. The coupling of the mediator with the Higgs boson, g Z B , is set to the ratio of the mediator mass to the vacuum expectation value (VEV) of the baryonic Higgs boson: m Z B /v B . The mediator is naturally leptophobic, thus evading the current constraints coming from the dilepton resonance searches. Different mediator and DM masses are investigated.

Flavour-changing interaction
The flavour-changing vector mediator model (VFC) [109] permits the interaction of the DM candidate with the top quark. A spin-1 colour-neutral mediator Z VFC enables a flavourchanging interaction of the DM with ordinary matter, for instance between the top quark and the up quark. For simplicity, the mediator is allowed to couple only to the right-handed component of the top-quark field [109]. This model predicts flavour-changing neutral current (FCNC) processes which are suppressed in the SM. A representative diagram of the on-shell production of the new mediator Z VFC is shown in figure 2b. The mediator can either decay invisibly, leading to a final state involving a single top quark and large missing transverse momentum, or decay visibly, producing a distinctive final state containing two top quarks with the same electric charge (tt/tt). A similar topology arises from the t-channel exchange of the Z VFC mediator, as depicted in figures 2c and 2d. The model is fully predictive once the four main parameters are specified [110]: the mass of the mediator m Z VFC , the DM mass m χ , and the couplings of the mediator to the DM particles and to the quarks, g χ and g ut , respectively. In the context of the analyses described in this paper, the mass of the DM candidate m χ has negligible impact on the kinematics, provided that m Z VFC > 2m χ , and it is fixed to 1 GeV. This reduces the number of dimensions of the parameter space to three. The sensitivity of the experimental analyses to this model is explored in three scenarios that investigate different parameter planes as a function of m Z VFC , g ut and the invisible branching ratio of the Z VFC mediator, B(χχ).

Colour-neutral interaction
In the scalar or pseudo-scalar simplified models (S/PS) a new spin-0 gauge particle mediates the interaction, at tree level, with a DM particle [75,113]. The mediator is considered to be either a scalar (φ) or a pseudo-scalar (a). This model has four parameters [25]: the mass of -6 -JHEP05(2019)142 the mediator m φ or m a ; the DM mass; the DM-mediator coupling, g χ ; and the coupling of the mediator with the SM fermions. The latter is composed of a flavour-universal term, g q , which is a free parameter of the model and multiplies the SM-Yukawa coupling for each of the fermions [113]. This particular form of interaction, common to all models with spin-0 mediators evaluated in this paper, is typically referred to as the minimal flavour violation (MFV) ansatz and by construction, it relaxes the severe restrictions on the coupling of new spin-0 colour-neutral particles to the SM fermions imposed by flavour measurements [136][137][138]. Furthermore, it implies that these mediators are sizeably produced through loopinduced gluon fusion or in association with heavy-flavour quarks (see figure 3). According to whether the mediator decays into a pair of DM or SM particles, different final states are sensitive to these models. Due to the Yukawa-like structure of the couplings, visible final states with two or four top quarks are particularly important signatures. Final states involving a single top quark and E miss T may also play an important role in constraining these models [139][140][141][142][143][144]. Despite the absence of a dedicated parameter that regulates the relative importance of up-type and down-type quark couplings (otherwise present in UV completions of these models as in section 2.3.2), it is also important to study final states involving bottom quarks separately, since these become a relevant signature if the up-type couplings are suppressed.

Colour-charged interaction
The scalar colour-charged interaction model (SCC) assumes that the scalar mediator couples to left-or right-handed quarks and it is a colour triplet. The DM particle(s) is produced via a t-channel exchange of the mediator which leads to a different phenomenology from that of colour-neutral interactions. These models have a strong connection with the minimal supersymmetric Standard Model (MSSM) [145,146] with a neutralino DM and firstand second-generation squarks with universal masses. They share with it the same crosssections and phenomenology when the mediator is pair-produced via strong interaction. Nevertheless, additional production diagrams are taken into account in this scenario, since values assumed for the couplings of the mediator to quarks and DM differ from those of the MSSM.
As in the case of the MSSM, it is reasonable to decouple the first two generations from the third, considering the different mass scales. For this purpose, three different models are considered: 2 1. In the SCC q model, the mediator, η q , couples to the left-handed quarks of the first and second generations and is a SU (2) singlet under the SM. The mediator decays into a quark-DM pair, so that the strongest sensitivity for these models is provided by searches involving jets and missing transverse momentum. The three model parameters are the mediator mass, the DM mass, and the flavour-universal coupling to quarks and DM, λ q .  Figure 4. Schematic representation of the dominant production and decay modes for the SCC models.
2. In the SCC b model, the mediator, η b , couples to the right-handed bottom quark. Following previous publications [25,147], the specific realisation of this model is obtained within the framework of "flavoured" DM, where the DM candidate is the lightest component of a flavour triplet [131]. With these assumptions, the mediator always decays into a b-quark-DM pair. Of the three parameters of the model, the mediator and DM masses and the coupling, λ b , only the first two are varied, while the last one is set to the value predicting a DM relic density compatible with astrophysical observations [136]. Representative diagrams for these models are presented in figures 4b and 4c.
3. In the SCC t model, the mediator, η t , consists of a SU(2) L -singlet field that couples to right-handed quarks, and is produced by down-type quark-anti-quark fusion, and it decays into a top quark and a DM particle. The representative diagram is shown in figure 4d. This specific realisation of the model [109], which gives rise to a characteristic signature composed of a single top quark and an invisible particle, can be related to the MSSM if an additional R-parity violating interaction of the top squark with the down-type quarks is assumed. The coupling strength of the mediator to DM and top quarks, denoted by λ t , and the coupling strength to light-flavour down-type quarks, g ds , are free parameters of the model.

Extended Higgs sector dark matter models
The third category of models aims to extend the simplified DM mediator models by involving an extended two-Higgs-doublet sector (2HDM) [148][149][150][151][152][153][154][155][156], together with an additional mediator to DM, either a vector or a pseudo-scalar. This embeds the simplified models in a UV-complete and renormalisable framework and allows the investigation of a broad phenomenology predicted by these types of models. In both models, the 2HDM sector has a CP-conserving potential and a softly broken Z 2 symmetry [157], and the alignment limit is assumed, so that the lightest CP-even state, h, of the Higgs sector can be identified with the SM Higgs boson. Figure 5. Schematic representation of the dominant production and decay modes for the 2HDM+Z V model.

Two-Higgs-doublet models with a vector mediator
The first two-Higgs-doublet model [158], denoted for brevity 2HDM+Z V in the following, is based on a type-II 2HDM [157,159] [160], and therefore these signatures are not considered further in the context of this interpretation. The model has six parameters [160]: tan β, the ratio of the vacuum expectation values of the two Higgs doublets, is set to unity; m χ , the DM mass, is set to 100 GeV; and g Z , the coupling of the new Z V U(1) gauge symmetry, is set to 0.8. The masses m h and m H = m H ± of the two CP-even and charged Higgs bosons are set to 125 GeV, and 300 GeV, respectively, while m A , the mass of the CP-odd Higgs partner and m Z V are free parameters and varied in the interpretation.

Two-Higgs-doublet models with a pseudo-scalar mediator
The second 2HDM model [152], 2HDM+a, includes an additional pseudo-scalar mediator, a. In this case also, the 2HDM coupling structure is chosen to be of type-II, although many of the interpretations in this paper hold for a type-I case too. The additional pseudo-scalar mediator of the model couples the DM particles to the SM and mixes with the pseudoscalar partner of the SM Higgs boson. The physics of the model is fully captured by 14 parameters: the masses of the CP-even (h and H), CP-odd (a and A) and charged (H ± ) bosons; the mass of the DM particle (m χ ); the three quartic couplings between the scalar -9 -

JHEP05(2019)142
doublets and the a boson (λ P1 , λ P2 and λ 3 ); and the coupling between the a boson and the DM, y χ ; the electroweak VEV, v; the ratio of the VEVs of the two Higgs doublets, tan β; and the mixing angles of the CP-even and CP-odd weak eigenstates, α and θ, respectively. The coupling y χ = 1 is chosen, having a negligible effect on the kinematics in the final states of interest. The alignment and decoupling limit (cos(β − α) = 0) is assumed, thus h is the SM Higgs boson and v = 246 GeV. The quartic coupling λ 3 = 3 is chosen to ensure the stability of the Higgs potential for our choice of the masses of the heavy Higgs bosons which are themselves fixed to the same value (m A = m H ± = m H ) to simplify the phenomenology and evade the constraints from electroweak precision measurements [152]. The other quartic couplings are also set to 3 in order to maximise the trilinear couplings between the CP-odd and the CP-even neutral states. This model is characterised by a rich phenomenology. The production of the lightest pseudo-scalar is dominated by loop-induced gluon fusion, followed by associated production with heavy-flavour quarks or associated production with a Higgs or Z boson (figures 6a-6c). Furthermore, according to the Higgs sector's mass hierarchy, Higgs and Z bosons can be produced in the resonant decay of the heavier bosons into the lightest pseudo-scalar (see for example figures 6d-6f). The pseudo-scalar mediator can subsequently decay into either a pair of DM particles or a pair of SM particles (mostly top quarks if kinematically allowed), giving rise to very diverse signatures. The four-top-quark signature [161] is particularly interesting in this model if the neutral Higgs partner masses are kept above the tt decay threshold, since, when kinematically allowed, all heavy neutral bosons can contribute to this final state, as depicted in the diagram of figure 6c. Four benchmark scenarios [78] that are consistent with bounds from electroweak precision, flavour and Higgs observables are chosen to investigate the sensitivity to this model as a function of relevant parameters: m a , m A , tan β, sin θ and m χ .

EFT model of scalar dark energy
The Horndeski theories [94] introduce a dark energy scalar which couples to gravity and provide a useful framework for constraining the cosmological constant problem and the source of the acceleration of the expansion of the universe. The model considered in this paper is an EFT implementation of these theories [82]. In this model, the dark energy field is assumed to couple to matter universally. The model contains two classes of effective operators: operators which are invariant under shift-symmetry φ DE → φ DE + constant, where φ DE denotes the DE scalar field, and operators which break this symmetry. Shiftsymmetric operators contain derivative interactions of φ DE with the SM particles, while operators that break the shift-symmetry contain direct interactions of φ DE with the SM. In the former case the DE scalar is pair-produced and does not decay within the volume of collider experiments, thereby resulting in E miss T in the detector, while the latter case includes Yukawa-type interactions φ DEψ ψ, which allow the scalar to decay into SM fermions, thereby changing the expected signatures. The interactions arising from the shift-symmetry breaking operators are tightly constrained [162] and are not evaluated here.
There are nine shift-symmetric Lagrangian effective operators in the model, each suppressed by powers of a characteristic energy scale M according to the operator's dimen- Figure 6. Schematic representation of the dominant production and decay modes for the 2HDM+a model. sionality: where d is the operator's dimension and c i are the Wilson coefficients. Operators L 1 -L 5 correspond to interactions of the DE field with SM fields. The leading, i.e. least suppressed, operators of dimension eight are where T µν is the energy-momentum tensor corresponding to the SM Lagrangian. The L 1 operator corresponds to a derivative coupling of the DE field to the conformal anomaly, T ν ν (= mψψ for a Dirac field), and is therefore proportional to the mass of the SM fermions to which DE couples. Signatures which probe DE production in association with tt are therefore the most sensitive to this type of coupling and are used here. The L 2 operator involves derivatives of the SM fields and is therefore proportional to their momenta. Final states involving large momentum transfers, such as the jet+E miss T signature, offer the highest sensitivity to this type of coupling. The L 1 and L 2 operators are referred to as (kinetically dependent) conformal [163] and disformal, respectively. Operators L 3 -L 5 correspond to higher-order versions of L 1 and L 2 . The operator L 6 corresponds to a generalised kinetic term for the DE scalar and operators L 7 -L 9 correspond to the non-trivial Galilean terms [164]. In this paper, only L 1 and L 2 are considered. Due to the absence of terms allowing the decay of the DE scalars into SM particles, the DE particles (φ DE ) are stable and they escape the detector producing a missing-momentum signature.
The validity of the EFT approach in the context of collider data [165][166][167] is assessed with the procedure described in ref. [75], imposing the condition √ŝ < g * M , where √ŝ is the centre-of-mass energy of the hard interaction and g * is the effective coupling associated with the UV completion of the EFT.
Representative Feynman diagrams corresponding to the L 1 and L 2 operators for the tt + E miss T and mono-jet signatures are shown in figure 7.

Dataset and signal simulation
This paper interprets analyses of pp collision data recorded at a centre-of-mass energy of √ s = 13 TeV by the ATLAS detector during 2015 and 2016. Unless otherwise specified, the integrated luminosity of the dataset, after requiring that all detector subsystems were operational during data recording, amounts to 36.1 ± 0.8 fb −1 .
Monte Carlo (MC) simulated event samples were used to aid in the estimation of the background from SM processes and to model the DM and DE signals. Simulated events were processed either through a detector simulation [168] based on Geant4 [169] or through a fast simulation [168] with a parameterisation of the calorimeter response and Geant4 for the other parts of the detector [170]. Either of these ATLAS detector simulations were used for background processes (details in the specific analysis references) and most of the signal processes, as detailed in the following.
Two sets of samples were used for the modelling of the signal processes considered in this paper. One set of samples is based on signal events processed through the ATLAS detector simulation, referred to as "reconstructed" samples. The second set of samples consists of signal events composed of particle-level objects, defined according to the guiding principles outlined in ref. [171], and not including any resolution effect due to the ATLAS detector. These are referred to as "particle-level" samples. Particle-level samples were used to define a rescaling procedure specifically designed to broaden the range of signal models and parameter choices considered in the interpretation of the results. The procedure allows the use of less extensive computational resources that would be needed to provide a JHEP05(2019)142 full detector simulation for the large set of considered signals, while providing a complete picture of the current experimental coverage for these models. The rescaling procedure calculated a set of correction weights for a reference model as the ratio of the acceptance for a baseline signal sample to the acceptance of the signal sample of interest. Both of these acceptances are derived in a particle-level simulation. These weights were then applied to the reconstructed baseline signal sample of the reference model, assuming similar detector effects for the two models. The same procedure was used in some cases to rescale between signal samples of the same reference model but for different parameter choices which affect the kinematics of the final state. Closure-tests were performed to determine the reliability of this procedure and assign specific systematic uncertainties when needed. Further details about the rescaling used in the V/AV, VFC and the 2HDM+a signal samples are given in appendix A.
The generation settings for signal models considering a spin-1 mediator are summarised in table 4. For each model the table indicates the Universal FeynRules Output (UFO) model [172] implementation, the matrix element (ME) generator, the parton shower (PS), and the cross-section normalisation, at QCD leading-order or next-to-leading order accuracy (LO and NLO, respectively). Following the notation of the previous section, the simplified models are indicated with Z V/A , while the baryon-charged and flavour-changing interactions are indicated as Z B and Z VFC , respectively. The 2HDM model with an additional vector mediator is indicated as 2HDM+Z V . When relevant for the generations settings, each separate final state considered in this paper is indicated for each model.
The generation settings for signal models considering a spin-0 mediator are summarised in table 5. Following the notation of the previous section, the colour-neutral (colourcharged) simplified models are indicated with φ/a (η). The 2HDM with additional pseudoscalar mediator is indicated as 2HDM+a.
Finally, the generation settings for the DE model are also indicated in table 6.

Experimental signatures
Dark matter searches are an important component of the ATLAS physics programme. Several final-state signatures are targeted to maximise the discovery potential. This section presents summaries of the different searches for DM and is not intended to be exhaustive. More details are available in the reference papers. Table 2 summarises the DM searches for invisible final states, while table 3 summarises the searches for visible final states. These tables also provide an overview of the models (table 1) which are constrained by each of these signatures and which of these intepretations have not been presented elsewhere. Electrons, muons, photons and jets are reconstructed by combining the signals from the different components of the ATLAS detector 3 [173][174][175][176][177]. Leptons ( ) in the following JHEP05(2019)142 refers to electrons or muons. In several analyses, events with identified leptons are rejected from the signal region selection. This is referred to as a lepton veto. The analyses may implement different lepton and photon selection criteria for particle identification [173,175,176,178], isolation [173,175,179], and kinematic requirements (p T , η). Small-R and large-R jets are reconstructed from energy deposits in the calorimeters using the anti-k t jet algorithm [180,181] and using a radius parameter of R = 0.4 and R = 1.0, respectively. Reclustered large-R jets are reconstructed from small-R jets using a radius parameter of either R = 0.8 or R = 1.2. Multivariate algorithms are used to identify small-R jets with p T > 20 GeV containing b-hadrons (b-jets) [182,183]. This is referred to as b-tagging. For large-R jets, b-tagging is applied to their associated track-jets, which are constructed from tracks reconstructed in the inner detector using the anti-k t jet algorithm with R = 0.2. The missing transverse momentum p miss T (with magnitude E miss T ) is calculated from the negative vector sum of transverse momenta (p T ) of electrons, muons and jet candidates and an additional soft term [184] which includes activity in the tracking system originating from the primary vertex but not associated with any reconstructed particle. Some analyses may also consider photons in the E miss T reconstructions.

Searches for invisible final states
Searches for WIMP candidates at the LHC are characterised by the requirement of large E miss T since WIMPs escape detection. Therefore, final states with additional visible particles are considered in the selection of the events. These additional particles may come from initial-state radiation or from associated production. Several signatures that are listed in the following are exploited and optimised to enhance the sensitivity to different DM models.
Jet +E miss T . The jet+E miss T analysis [26], commonly referred to as the mono-jet analysis, is characterised by the presence of an energetic jet and large E miss T . The analysis selects events with E miss T > 250 GeV, at least one jet with p T > 250 GeV, and at most three additional jets with p T > 30 GeV. Events are required to pass a lepton veto. To reduce the contribution from multi-jet background where large E miss T can originate from jet energy under-measurement, a minimum separation in the azimuthal angle between each selected jet and the E miss T direction is also required: ∆φ(jet, p miss T ) > 0.4. The W +jets, Z+jets, and top-quark-related backgrounds are constrained using MC event samples normalised to data in selected control regions containing leptons. In the case of W +jets and Z+jets events, MC predictions are reweighted to account for higher-order QCD and electroweak corrections as a function of the vector-boson p T [185]. The normalisation factors for these backgrounds are extracted simultaneously using a binned likelihood fit of the E miss T distributions in all control and signal regions that includes systematic uncertainties. The remaining SM backgrounds from diboson processes are determined using MC simulated samples, while the multi-jet background contribution is extracted from data. plane, φ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). Angular distance is measured in units of ∆R ≡ (∆η) 2 + (∆φ) 2 . The rapidity is defined at y = 1/2 ln[(E + pz)/(E − pz)], where E is the energy and pz is the component of its momentum along the z-axis. The rapidity difference between two jets is defined as y * = (y1 − y2)/2. h(inv). Searches for invisible Higgs boson decays have been performed using several production and decay channels at a centre-of-mass energy of √ s = 8 TeV [29]. Results of searches in the vector-boson fusion (VBF) production channel and in associated production of a Higgs boson with a W/Z boson are statistically combined with the measured production and decay rates of the Higgs boson in the γγ, ZZ, W W , Zγ, bb, τ τ , and µµ channels to set an upper limit on the Higgs boson's invisible branching ratio of 0.23 at 95% confidence level (CL). This combined limit is used in the results quoted in section 6. Among the direct searches, the VBF production of Higgs bosons decaying into invisible particles [186] is the most sensitive one, setting an upper limit on the invisible branching ratio of 0.28. The VBF+E miss T analysis requires E miss T > 150 GeV and two jets with p T > 35 GeV. Three orthogonal signal regions are defined by varying the threshold imposed on the leading jet p T and the invariant mass of the two jets. Additional requirements on the angular separation of the two jets are applied to enhance the sensitivity to VBF production. In particular, the two leading jets are required to be well separated in pseudorapidity. Lepton and bjet vetoes are applied to reduce contamination from W +jets and top-quark backgrounds, respectively. Dedicated control regions with one and two leptons in the final state are used to constrain the contributions from dominant Z/W +jets backgrounds, through a simultaneous fit together with the signal region. The multi-jet background is estimated using a data-driven technique. Searches for Zh(inv) and V h(inv) [20,24,187] have been performed at centre-of-mass energy √ s = 13 TeV. Constraints using a VBF+E miss T analysis are also available using √ s = 13 TeV pp collision data [188,189]. However, the 8 TeV combination gives more stringent limits, thus it is used here.
γ + E miss T . Events in the γ + E miss T analysis [21] are required to pass the lepton veto and to have a photon with E T > 150 GeV. Events with more than one jet (p T > 30 GeV) or with a jet fulfilling ∆φ(jet, p miss T ) < 0.4 are rejected. Three exclusive signal regions with E miss T ranges between 150 GeV, 225 GeV, 300 GeV and above are defined. The W γ, Zγ, and γ+jets backgrounds are normalised in control regions using a simultaneous likelihood fit of all E miss T regions, but with independent normalisation factors for each region. The backgrounds due to photons from the misidentification of electrons or jets in processes such as W/Z+jets, diboson, and multi-jet events are estimated using data-driven techniques.
The event selection criteria in this analysis [24] require large E miss T and a pair of high-p T leptons. Two opposite-sign, same-flavour leptons satisfying p T > 30 GeV and p T > 20 GeV are required. The lepton pair is required to have an invariant mass between 76 GeV and 106 GeV to be consistent with originating from a Z boson. Events with an additional lepton with p T > 7 GeV or a b-jet with p T > 20 GeV are vetoed. To target events consistent with a boosted Z boson produced in the direction opposite to p miss T , additional requirements are imposed on the azimuthal angle between the dilepton system and p miss T and on the angular distance between leptons. A single inclusive E miss T signal region is defined with E miss T > 90 GeV for each of the ee and µµ channels. The dominant background in this analysis, ZZ production, is estimated from MC simulation. The W Z background is normalised to data in a three-lepton control region. The contributions from Z+jets and non-resonant-backgrounds are estimated using data-driven techniques. A statistical combination of the two decay channels is used for the final results. . This analysis [20] selects events with E miss T > 150 GeV and a hadronically decaying W or Z boson candidate. The vector-boson candidate is defined with one large-R jet with p T > 250 GeV in a boosted topology (E miss T > 250 GeV) or with two small-R jets with p T > 20 GeV in a resolved topology. In both cases, a lepton veto is applied. Additional requirements are applied to the invariant mass of the boson candidate. Several signal regions are defined according to the b-jet multiplicity. Similarly, several control regions are defined according to lepton and b-jet multiplicity. The normalisations of the tt and W/Z+jets background processes are constrained using a simultaneous fit of all control and signal regions of the E miss T distribution. The subdominant contribution from diboson and single-top-quark production is obtained from simulation. Multi-jet contributions are estimated with a data-driven technique. ), between the missing momentum direction and the direction of any selected jet in the event. Two sets of signal regions are defined targeting moderatemomentum (resolved) and high-momentum (boosted) Higgs boson candidates. In each case, the regions are further split according to whether there are one or two b-jets. The resolved regime, defined in three exclusive E miss T regions between 150 GeV and 500 GeV, selects a Higgs boson candidate reconstructed from the two leading b-tagged small-R jets (or from a b-tagged and a non-b-tagged small-R jet) with p T > 20 GeV. In the boosted regime, defined by E miss T > 500 GeV, the leading large-R jet with p T > 200 GeV is the Higgs boson candidate. The b-jet multiplicity is defined by the number of b-tagged track-jets associated with the large-R jet. Backgrounds involving the production of W/Z bosons in association with heavy-flavour quarks or top-quark pairs are normalised in dedicated control regions distinct from the signal regions by requiring one or two leptons. A simultaneous binned likelihood fit to the invariant mass of the Higgs boson candidate is performed in all signal and control regions. The multi-jet background is obtained with a data-driven technique. Other subdominant backgrounds are estimated from simulation.
The h(γγ) + E miss T events in this analysis [22] are selected by requiring at least two photons with p T > 25 GeV. The two leading photons are chosen to reconstruct the Higgs candidate, which is required to satisfy 105 GeV < m γγ < 160 GeV. The leading (sub-leading) photon is also required to have E γ T /m γγ > 0.35 (0.25). Events with leptons are vetoed. Events with p T (γγ) > 90 GeV and E miss T / E T > 7 GeV 1/2 in ref.
[22] are used for the interpretation of DM models, where E T is the scalar sum of the transverse momentum of all reconstructed objects in the event. The backgrounds are extracted by fitting an analytic function to the diphoton invariant mass distribution. In the case of the non-resonant background, the normalisation and shape are obtained by fitting the invariant mass distribution in data to an exponential function. The SM Higgs boson background shape is modelled with a double-sided Crystal Ball function and fitted to the MC simulation. The t+E miss T analysis [30] searches for events with one top quark and relatively large E miss T . Two signal regions are defined depending on the decay channel of the top quark. The leptonic channel selects events with a positively charged lepton with p T > 30 GeV, E miss T > 50 GeV, and transverse mass of the lepton and the E miss T , m W T , larger than 260 GeV. One b-jet with p T > 30 GeV is additionally required. The hadronic channel is optimised to select events with a top quark produced with a large boost. Events are selected with E miss T > 200 GeV and one large-R jet with p T > 250 GeV with one b-tagged track-jet associated with it. Events failing the lepton veto are rejected. Dedicated control regions are defined to constrain the tt and W/Z+jets backgrounds from data. The multi-jet background is estimated from data, whereas other remaining backgrounds are taken from simulation. All signal and control regions for the two decay channels are fitted simultaneously to extract the background normalisation. In the case of the hadronic channel, the transverse mass of the large-R jet and the E miss T are the discriminating variables, while for the leptonic channel, the E miss T distribution is used to discriminate signal from background.
selects events with two energetic jets (p T > 160 GeV), at least one of which is b-tagged, E miss T > 650 GeV, and additional total hadronic energy restricted to be less than 100 GeV. This last requirement rejects top-quark background. The dominant background for this analysis, Z+jets events, is constrained from data in a dedicated control region, which is fitted together with the signal region. The bb +E miss T analysis [25] instead exploits a selection with at least two b-jets and a considerably lower E miss T requirement, E miss T > 180 GeV. The azimuthal separations between the b-jets and p miss T are exploited to enhance the separation between the signal and the irreducible background in this channel (Z(νν)+bb), which is constrained from data in a dedicated control region. The results are extracted by fitting an observable that relies on the pseudorapidity difference between the two b-jets: cos θ * bb = |tanh (∆η bb /2)|.
tt + E miss T . The tt +E miss T analysis [25, 27, 28] is split into three channels according to the decays of the W bosons from the top-quark decays: 0-lepton, where both W bosons decay hadronically, 1-lepton, where one of the two W bosons decays leptonically and 2-leptons where both W bosons decay leptonically. The analyses targeting the 0-lepton channel exploit two sets of signal topologies: spin-0 DM models [25], used for the DM interpretations presented in this paper, and top-squark decays into a top quark and a neutralino [28], used for the DE interpretation in this paper. Both analyses are characterised by a set of signal regions which select events with at least four energetic jets, at least two of which are b-tagged, and relatively high E miss T . Requirements on the invariant mass of reclustered large-R jets are imposed to identify events where a W boson or a top quark are boosted. The dominant backgrounds (Z+jets, top-quark processes and tt + Z) are constrained in dedicated control regions. The three signal regions used for the DE interpretation are statistically combined, while the two signal regions in the DM analysis are not. The analysis targeting the 1-lepton channel selects events with at least four energetic jets, at least one of which is b-tagged, one isolated lepton and high E miss T . The events are also required to have at least one hadronic top candidate with invariant mass loosely compatible with the mass of the top quark. Requirements Analytic function fit of mγγ. , large-R jet)) in the leptonic (hadronic) channel. are used to suppress semileptonic and dileptonic tt events, respectively. Requirements on the azimuthal angle between the lepton and p miss T and on ∆φ(jets, p miss T ) are also exploited to further suppress the background contamination of the signal regions. All top-quark background processes are estimated in dedicated control regions. Finally, the analysis targeting the 2-lepton channel selects events with two opposite-sign leptons which are inconsistent with being produced in the decay of a Z boson. At least one b-jet is also required in the selections. The  visible particles, thus requiring the presence of reconstructed objects such as jets or leptons, covering a variety of kinematic regions. In some of the analyses described below, further identification techniques are employed to select final states with top quarks.

Searches for visible final states
Dijet. For this analysis [190] events with at least two small-R jets are selected if the p T of the leading (sub-leading) jet is greater than 440 (60) GeV. The dijet selection requires a rapidity difference |y * | < 0.6 and the invariant mass of the dijet system to be m jj > 1.1 TeV. The background estimation is obtained by fitting the falling m jj distribution. Bin widths are chosen to approximate the m jj resolution, and thus are wider for higher masses. A sliding-window fitting technique is used, where restricted regions of the spectrum are fitted with a functional form. The background is constructed bin-by-bin by performing a likelihood fit to the data in each window and using the fit value in the central bin for the background estimate. The values from the full set of windows are then combined to create the background estimate for the full mass range. Model-independent limits on the visible crosssection for a hypothetical signal that produces a Gaussian contribution to the m jj distribution (for several signal widths) are provided for this analysis (see appendix A of ref. [191]). This analysis was performed in data corresponding to an integrated luminosity of 37.0 fb −1 .
Dijet angular. A dijet selection can also be exploited to search for deviations from the SM expectation in angular distributions, characteristic of wider resonances where the nominal dijet search would lose sensitivity. A dijet angular analysis [190] is performed in events with two jets following the p T requirements of the dijet search, but relaxing the |y * | requirement to 1.7. Due to different kinematics in this loosened selection, the mass of the dijet pair is required to be m jj > 2.5 TeV. The analysis makes use of the variable -20 -JHEP05(2019)142 4 constructed so that, in the limit of massless parton scattering and when only the t-channel scattering contributes to the partonic cross-section, the angular distribution dN/dχ jj is approximately independent of χ jj . MC events from multi-jet production are modelled at LO in QCD, and reweighted to NLO predictions from NLOJET++ [192,193] using mass-and angle-dependent correction factors. Additional electroweak mass-and angle-dependent correction factors are applied. The data are compared with a SM template in different m jj ranges, and different χ jj bins.
Trigger-object-level dijet. For the dijet analysis described before, the high p T threshold imposed on the leading jet is limited by the trigger selection driven by the bandwidth available for single-jet triggers, thus it only targets m jj > 1.5 TeV. The limitation from the high-level trigger selection is overcome by recording only high-level trigger jet information, rather than the full detector readout, to a dedicated data stream, reducing the storage needs per event. This strategy allows to record all events passing the single-jet level-one (L1) trigger (with lower threshold than in the high-level trigger) with minimal bandwidth increase. The dataset collected corresponds to an integrated luminosity of 29.3 fb −1 . This trigger-object-level dijet analysis (TLA dijet) [194] selects events with at least two triggerlevel jets with p T > 85 GeV. Two selection criteria are used: |y * | < 0.6 in the mass range 700 GeV < m jj < 1.8 TeV and |y * | < 0.3 for 450 GeV < m jj < 700 GeV. The leading trigger-level jet is required to have p T > 185 GeV and p T > 220 GeV for the |y * | < 0.3 and |y * | < 0.6 selections, respectively, to ensure full efficiency for the L1 triggers. The search is then interpreted in terms of resonances with a mass between 450 GeV and 1.8 TeV. The background strategy used in the dijet search is also used here.
Resolved dijet + ISR. Another alternative strategy to search for low-mass dijet resonances is to select events with a pair of jets recoiling against a photon or an additional jet from ISR. The resolved dijet + ISR analyses [195] select events with a high-p T ISR object (γ or jet), used to trigger the event, and a relatively low mass dijet resonance. Dijet+γ events contain at least one photon with p γ > 150 GeV and at least two jets with p T > 25 GeV. The two leading jets must satisfy |y * | < 0.8, which allows to probe of dijet invariant masses between 170 GeV and 1.5 TeV. The three-jet selection requires at least one jet with p T > 430 GeV as well as two additional jets with p T > 25 GeV. The leading jet is chosen as the ISR candidate, and the second-and third-highest-p T jets are required to satisfy |y * | < 0.6. This selection probes a mass range between about 300 GeV and 600 GeV. The background contribution is estimated by fitting the m jj distribution. This analysis was performed in 13 TeV collision data corresponding to an integrated luminosity of 15.5 fb −1 .
Boosted dijet + ISR. In the case of a dijet+ISR selection, if the associated ISR photon or jet has large transverse momentum, the dijet resonance candidate is reconstructed as a large-R jet [196] of radius 1.0 with mass m. To enhance the sensitivity to quark pair decays, jet substructure techniques are used to discriminate between a two-particle jet from a decay JHEP05(2019)142 of a boosted resonance and a single-particle jet [197]. Events are required to have a large-R jet, the resonance candidate, and at least one ISR object candidate. The azimuthal angular separation between the resonance candidate and the ISR object is required to satisfy ∆φ > π/2. A p T > 2m requirement ensures sufficient collimation of the resonance candidate. In the ISR jet (photon) channel, the large-R jet satisfies p T > 450 (200) GeV and the ISR jet (photon) has p T > 450 (155) GeV. A data-driven technique is used to model the expected background in the signal region via a transfer factor that extrapolates from a control region with inverted jet substructure requirements.
Dibjet. The dibjet search [198] targets dijet resonances with one or two jets identified as b-jets. Two different analyses cover both the low and high invariant mass regions. Events in the high invariant mass region are selected with at least two jets, one of which has p T > 430 GeV and passes an inclusive jet trigger. The rapidity difference is required to be |y * | < 0.8. This analysis covers the region with m jj > 1.2 TeV. The low invariant mass region uses a trigger targeting events with two jets containing b-hadrons, which provides access to lower dibjet invariant masses (m jj ) compared to the single jet trigger: 570 GeV < m jj < 1.5 TeV. The rapidity difference requirement is tightened to |y * | < 0.6. In this case, only the two-b-jets selection is considered. Because the double b-jet trigger was not available during the full data-taking period, the total integrated luminosity used for the low-mass analysis corresponds to 24.3 fb −1 of 13 TeV collision data. A background estimation strategy similar to that of the dijet analysis is used in these analyses.
Dilepton. The dilepton analysis [199] selects events with at least two same-flavour leptons. The pair of electrons (muons) with highest E T (p T ) are chosen as the candidate decay products of the resonance. Only the muon channel candidates are required to have opposite charge, due to higher charge misidentification for high-E T electrons and the p T misreconstruction associated with wrongly measured charge in muons. Background processes with two prompt leptons are modelled using MC samples. The Z/γ * → background is smoothed for 120 GeV < m < 1 TeV. This is done by fitting the MC spectrum and the resulting fitted function is used to set the expected event yields in that mass range. A data-driven method is employed to estimate backgrounds with at least one misidentified lepton. The m distribution is explored between 80 GeV and 6 TeV.
Same-sign tt. Events in the same-sign tt analysis [110] are selected with exactly two leptons with positive charge and at least one b-jet. Events are required to satisfy H T > 750 GeV, where H T is defined as the scalar sum of the p T of all selected objects, including jets. Additionally, requirements on E miss T and the azimuthal separation between the two leptons are imposed. Signal regions for the different lepton flavours (ee, eµ and µµ) are treated separately. Irreducible SM backgrounds are determined using MC simulation samples. Backgrounds from fake leptons are estimated using data-driven techniques. Dijet [190] V/AV 2 jets, mjj, y * . Sliding-window fit of the mjj distribution. Section 6.1.1 Dijet angular [190] V/AV 2 jets, mjj, y * . Binned likelihood fit of χjj. Section 6.1.1 TLA dijet [194] V/AV 2 trigger-level jets, mjj, y * . Sliding-window fit of the mjj distribution. Section 6.1.1 Resolved dijet+ISR [195] V/AV 3 jets (or 2 jets and 1 photon), mjj, y * . Parametric function fit of the mjj distribution. Section 6.1.1 Boosted dijet+ISR [196] V/AV( * ) 1 large-R jet, 1 jet or photon, mJ . Data-driven extrapolation from control region via transfer factor. Section 6.1.1 Dibjet [198] V/AV 2 jets (1 and 2 b-jets), mjj, y * . Sliding-window parametric fit of the mjj distribution. Section 6.1.1 Dilepton [199] V/AV 2 e or 2 µ. Z/γ * → from fitted MC spectrum.  close by. Events are classified as boosted or resolved depending on their hadronic activity.
In the boosted selection, events contain one large-R jet passing top-tagging requirements.
In the resolved selection, events have at least four small-R jets and fail the boosted selection.
The tt invariant mass, m tt , is reconstructed from the decay products of the two top-quark candidates in the event. The b-jet multiplicity is used for further categorisation. The SM tt production is estimated using MC samples and fixed-order theory calculations. The multi-jet and W +jets background contribution is estimated using data-driven techniques.
tttt. The tttt analysis [201] searches for events characterised by a single lepton and high jet multiplicity. Events are categorised by their jet multiplicity, which is defined using three p T thresholds: 40 GeV, 60 GeV, and 80 GeV. Events are further separated into five bins corresponding to the b-jet multiplicity. The tt+jets and W/Z+jets background production is estimated using a combined fit to these jet and b-jet multiplicity bins. The normalisation of these backgrounds is extrapolated from lower to higher jet multiplicity, while the b-jet multiplicity shape is taken from a parameterised extrapolation from data (simulation) for the tt (W/Z+jets) background. As in the case of the searches for invisible final states, these analyses found no significant deviation from the expected SM backgrounds. Therefore, exclusion limits are placed on the allowed phase space of the corresponding signal models, as discussed in section 6.

Complementarity and combination of signatures
It can be seen from tables 2 and 3 that several analyses are often sensitive to the same model. In cases like the X +E miss T final-states searches, X originating from initial-state radiation or -23 -JHEP05(2019)142 associated production, a variety of final states are evaluated: X = jet, γ, W, Z, h, t(t), b(b). Since the mediator couples DM to SM particles, it is also possible to reinterpret results from resonance searches targeting the mediator directly. The complementarity depends on the choice of model as well as coupling values. For the V/AV model, this paper considers jet/γ/V + E miss T , dijet, dibjet and tt resonance final states. All results for this model are new or have been updated from previous publications, except for the jet/γ + E miss T interpretations. For the VBC and 2HDM+Z V models, this paper considers only h + E miss T , which dominates the sensitivity for the chosen parameter values. All possible final states (same-sign tt and the t + E miss T ) are taken into account for the VFC model. Only invisible final states, tt/bb/jet + E miss T , are considered for the S/PS model. The SCC q , SCC b and SCC t models are each addressed with a specific signature: jet + E miss T , b + E miss T , t + E miss T , respectively, and all results were presented in each specific analysis paper. Various final states are evaluated in order to place the first constraints on the 2HDM+a model by ATLAS searches: Z/h + E miss T , tt/bb+ E miss T , h(inv) and tttt. Finally, the constraints on DE models are set using jet+E miss analysis has a larger reach in mediator masses, however its sensitivity is limited at lower masses by the threshold requirement of the E miss T trigger used to record the events for this analysis. The h(γγ) + E miss T analysis covers a lower mass region owing to its selection based on a photon trigger. For the combination, the luminosity, experimental, and signal modelling uncertainties were taken to be correlated between the two channels. In the h(γγ) + E miss T analysis the systematics uncertainties are not significantly constrained by the fit. This is mainly due to the use of a single signal region and no control regions. In the case of h(bb) + E miss T the systematics uncertainties are constrained due to the use of dedicated control regions. It is observed, however, that the results from the combination and the individual h(bb)+E miss T results are very similar. While the h(bb) + E miss T channel dominates the sensitivity, searches in different decay channels are of interest in probing different kinematic regions defined by different analyses strategies.
For this paper, the tt+E miss T exclusion limits discussed in section 6 are combined based on the best expected exclusion for each signal model, unless separate contours are shown. The combination and comparison of this channel is also a novel result of this paper.

Systematic uncertainties
Systematic uncertainties for both the background and signal models are considered in each of the analyses presented in section 4. These uncertainties, as well as statistical uncertainties, depend on the event selection, the phase space covered by a given analysis, and its background estimation strategy. The systematic uncertainties include experimental and theoretical uncertainties. Experimental uncertainties may include uncertainties in the ab- The signal modelling is subject to experimental and theoretical uncertainties. The experimental uncertainties are the same as for the background processes. Theoretical uncertainties affecting the production cross-section (normalisation) and the acceptance are considered separately. The strategy used to estimate systematic uncertainties for those signal models studied in this paper which are not discussed in previous publications is outlined below.
The results for the 2HDM+a and DE signal models include theoretical systematic uncertainties due to parton distribution functions (PDFs), evaluated following the PDF4LHC recommendations [202]. The choice of different PDFs results in up to 2% uncertainty in the acceptance and up to 10% uncertainty in the cross-section. Uncertainties related to the choice of renormalisation and factorisation scales are derived by varying independently such scales by a factor of 2.0 and 0.5 relative to the nominal values used for the MC generation. They account for an uncertainty in the acceptance below 5% for the different analyses. Uncertainties in initial-and final-state radiation due to the parton shower modelling are estimated by generating MC samples with alternative underlying event and multi-parton interaction parameter choices resulting in uncertainties between 5% and 15% in the signal acceptance, typically increasing at higher mediator masses. In the very large jet multiplicity phase space of the tttt analysis they reach values of 50%.
In some cases, additional uncertainties are included to account for non-closure effects of the rescaling procedure described in section 3. These uncertainties include a contribution from the statistical uncertainty associated with the acceptance ratios determined from the baseline signal sample. For the h(γγ) + E miss T (Z( ) + E miss T ) analysis this translates in up to 7% (8%) uncertainty in the final 2HDM+a signal yields.

Interpretation of the results
This section summarises the exclusion limits set by ATLAS published searches briefly outlined in section 4, on the various signal models described in section 2 (following the notation in table 1). The analyses and corresponding signal regions are referred to by their analysis labels defined in tables 2 and 3. The observed and expected 95% confidence level (CL) exclusion limits are obtained from the signal region or combination of regions of each contributing analysis using the CL s [203] method. The signal contamination in CRs across different analyses are kept minimal through dedicated selection requirements and hence are not explicitly taken into account in the following results. In section 6.3.1 a statistical combination of the h + E miss T final states is used to derive the results.  The V/AV simplified model is strongly constrained by searches for a high-mass resonance decaying into a pair of fermions and searches for associated production of DM particles with an ISR object.
As presented in figure 10 for the case of an axial-vector mediator, each resonance search analysis is sensitive to complementary regions of the mass-coupling parameter space. Couplings above the exclusion line are excluded, as long as the width predicted by the model is smaller than the maximal ratio of width to mediator mass (Γ/m Z ) to which the analysis is sensitive. This limitation arises where the background model is estimated via a sliding-window fit of the m jj distribution. Specifically, the TLA dijet analysis assuming |y * | < 0.6 is sensitive up to around Γ/m Z = 7%, the TLA dijet analysis requiring |y * | < 0.  Figure 11. Regions in a (mediator-mass, DM-mass) plane excluded at 95% CL by dijet, dilepton and E miss T +X searches, for leptophobic (a) or leptophilic (b) vector mediator simplified models described in section 2.1.1. The exclusions are computed for a DM coupling g χ , quark coupling g q , universal to all flavours, and lepton coupling g as indicated in each case. Dashed curves labelled "thermal relic" correspond to combinations of DM and mediator mass values that are consistent with a DM density of Ωh 2 = 0.12 and a standard thermal history as computed in MadDM [77,204]. Above the curve in (a) annihilation processes described by the simplified model deplete Ωh 2 to below 0.12. In (b), this occurs between the two dashed curves. The dotted line indicates the kinematic threshold where the mediator can decay on-shell into DM. The boosted dijet+ISR analysis has the best reach for low masses, excluding Z A mediator masses between 100 GeV and 220 GeV. Two new interpretations, for the dibjet and tt resonance analyses, are presented for these models. The dibjet (tt resonance) analysis places constraints on Z A mediators with masses between 500 GeV and 2.5 (2) TeV, in the same region of sensitivity of the dijet, TLA dijet, and boosted dijet+ISR analyses.
To illustrate the complementarity of the searches [75,77], three different coupling scenarios are also considered in the interpretation of the results: In particular, the lower lepton coupling value is set to highlight the dilepton search sensitivity even for very small values of this parameter.
The exclusions from the resonance searches (dijet, dibjet, dilepton) in the (m Z V/A , m χ ) plane are derived from the limits placed on resonances reconstructed with a Gaussian shape, while the limits from the E miss T +X and tt resonance analyses are derived using a mixture of simulated signal samples and rescaling procedures as described in section 3. For each scenario in figures 11 and 12, dashed curves labelled "thermal relic" correspond to combinations of DM and mediator mass values that are consistent with a DM density of Ωh 2 = 0.12 and a standard thermal history, as computed in MadDM [77,204]. Between the two curves, annihilation processes described by the simplified model deplete the relic density to below the thermal value (except for figure 11a, where this occurs above the dashed curve). A dotted line indicates the kinematic threshold where the mediator can decay on-shell into DM particles. Excluded regions that are in tension with the perturbative unitary considerations of ref. [205] are indicated by shading in figure 12.
The sensitivity reach of the various experimental signatures for the leptophobic vectormediator scenario as a function of the DM and mediator masses is summarised in figure 11a. Since the chosen universal quark coupling is relatively high in comparison with other benchmarks, the strongest limits are obtained from the resonance searches. These analyses are sensitive to mediator masses between 200 GeV and 2.5 TeV with little dependence on the DM mass. Opening of the Z V → χχ decay channel, significantly reduces the sensitivity observed at high mediator masses and for 200 < m Z V < 450 GeV when m Z V > 2m χ . The boosted dijet + ISR search is not reinterpreted here but has sensitivity in this region. The lower limit on the mass is determined by the trigger requirements of the resolved dijet + ISR analysis. For m Z V < 2m χ , masses up to 2.9 TeV are excluded by the resolved dijet + -29 -

JHEP05(2019)142
ISR, dijet TLA and dijet searches. Compared to the dijet searches, the tt resonance analysis is particularly sensitive to the reduction in effective cross-section related to changes of the branching ratio, as can be inferred from the coupling reach of figure 10. Conversely, the sensitivity of the E miss T +X signatures is highest in the region m Z V > 2m χ , up to mediator masses of 1.5 TeV and provides unique coverage for masses below 500 GeV. The sensitivity of these analyses is strongly decreased for m Z V < 2m χ , where the DM particles are produced off-shell, with a consequent strong reduction of the production cross-section. For this reason, only the jet+E miss T and γ + E miss T analyses can probe the off-shell regime for this benchmark scenario, and only in the case of very low mediator and DM masses. It is important to highlight that if the value chosen for g q were reduced, the relative interplay between the dijet and E miss T +X searches would change, as exemplified by the change of the dijet limit in the different coupling scenarios described in the following.
The experimental limits for the leptophilic vector-mediator model are summarised in figure 11b. In this case, the mediator decay rates into quarks are reduced in favour of a higher branching ratio to DM particles, reducing the sensitivity of dijet searches to this scenario, whereas the leptonic branching ratio allows dilepton searches to impose constraints on a wide range of mediator masses. The gap in sensitivity of dijet searches around mediator masses of about 1.7 TeV is due to statistical fluctuation in the dijet spectrum. This benchmark highlights the complementarity among dijet, dilepton, and E miss T +X final states. In this case, dibjet and tt resonance searches are not included in the final result. The resonance searches exclude mediator masses between 150 GeV and 2 TeV (except for a small gap around 1.7 TeV), if m Z V < 2m χ and between 150 and 350 GeV for all DM masses. Complementarily, the E miss T +X searches exclude mediator masses up to 1 TeV for m Z V > 2m χ . Similar considerations can be made for the axial-vector mediator models, presented in figure 12, with the exception that in the presence of non-vanishing coupling to leptons (leptophilic scenario), the dilepton resonance search becomes by far the most sensitive analysis for this model, excluding the mass range 150 GeV < m Z A < 2.8 TeV for any DM mass. Also in this case, the sensitivity of this analysis increases when the Z A → χχ is kinematically forbidden and becomes independent of the DM mass above threshold. For m Z A < 2m χ , masses up to 3.5 TeV are excluded.
Collider experiments provide an approach to DM searches which is complementary to direct and indirect detection experiments [76]. It is therefore interesting and informative, though model-dependent, to compare the V/AV limits with the results from other DM searches. Figure 13 shows the translation of the V/AV model limits into limits on the spin-dependent χ-proton and spin-independent χ-nucleon scattering cross-sections as a function of the DM mass. The direct detection experiments dominate the sensitivity by a few orders of magnitude for DM masses above 10 GeV, thanks to coherence effects, the spin-independent interaction cross-section with heavy nuclei is enhanced by A 2 , where A is the number of nucleons in a nucleus. However, with the assumed coupling strengths, the analyses presented in this paper complement direct detection limits in the low DM mass range where the direct DM search experiments have less sensitivity due to the very low energy recoils that such low-mass DM particles would induce. The lower edge of the sensitivity contour for all analyses in figure 13a (figure 13b) is driven by the high-mass  figure 11a (12b), as the scattering cross-section limit is inversely proportional to the mediator mass reach (raised to the fourth power). Conversely, the upper edge of the scattering cross-section contour for the dijet and dilepton analyses is driven by their low-mass sensitivity limit due to the trigger requirements employed in these analyses. Further details of this comparison are discussed in appendix B.

Baryon-charged interaction
In the context of the VBC model, the results from the h(γγ)+E miss

Neutral flavour-changing interaction
For the VFC models, expected and observed limits from the t + E miss T and same-sign tt analyses are derived for each independent subprocess leading to the two signatures, as described in section 2.1.3 and schematically summarised in figures 2b-2d. These individual results are converted into limits for the complete VFC model following the rescaling procedure described in section A.2. These results were obtained subsequently to the original analyses publications.
The sensitivity of the experimental analyses to this model is explored in three scenarios that investigate different interpretation variables (where B is the invisible branching ratio of the mediator): Scenario 1 (m Z VFC , g ut ) interpretation plane, assuming g χ = 0.5 or g χ = 1.0.
The first scenario, presented in figure 15, directly identifies the constraints on the main parameters of the model. The two different g χ coupling values highlight the different contributions of the invisible final state (t + E miss T ), which probes g ut values down to 0.7 for 1 TeV mediators and can exclude couplings down to 0.13 for 1.5 TeV mediators when the DM coupling is set to unity. In this scenario the visible final state (tt) constrains couplings down to 0.3 for mediator masses up to 3 TeV and it is relatively independent of g χ . This is due to the fact that the sensitivity of this final state is dominated by the t-channel exchange of the mediator and therefore it is sensitive to g χ only through the total width of the mediator. In this result, only mediator masses above 1 TeV are under study. However, mediator masses down to 100 GeV are excluded by previous publications [207] for a coupling assumption of g q ∼ 3 · 10 −2 .
Limits are quite similar for the two DM coupling values for the same-sign tt analysis. The sensitivity of this final state is dominated by the t-channel exchange of Z VFC especially for the values of g ut which are probed. This process is sensitive to g χ only through the total with of the mediator involved in the propagator, and therefore does not change drastically when g χ is varied.
The second and third scenarios, presented in figure 16, further highlight the complementarity between the visible and invisible final states as a function of the couplings and the invisible branching ratio of the mediator. Couplings of the Z VFC mediator to SM fermions can be excluded down to 0.14 for any value of g χ or B(χχ).  Figure 15. Exclusion regions in the (m Z VFC , g ut ) plane from the t+E miss T and same-sign tt analyses for the VFC model. The observed exclusion is indicated for each of the two analyses by the filled area. The mass of the DM particle is set to 1 GeV and the DM coupling, g χ , is set to (a) 0.5 or (b) 1.  Figure 16. Exclusion limits from the t + E miss T and same-sign tt analyses for the VFC model as a function of the SM coupling g ut and (a) the DM branching ratio or (b) the DM coupling g χ . The observed exclusion is indicated for each of the two analyses by the filled area. The mass of the DM particle is set to 1 GeV and the mass of the Z VFC boson is set to 1 TeV. The dark shaded area corresponds to an invisible partial width of the mediator above 20%.

Colour-neutral interaction
The most stringent limits on S/PS models are obtained from tt +E miss T final states, which are studied in three channels assuming fully-hadronic, semileptonic and fully-leptonic top pair decays, respectively. The fully leptonic channel excludes scalar-mediator models with unitary couplings g χ = g q = g = 1 up to mediator masses of 45 GeV, setting in this mass range the strongest upper limits on the ratio of the signal production cross-section to the nominal cross-section (signal strength or σ/σ(g = 1.0)), as shown in figure 17a. In the case -34 -

JHEP05(2019)142
of pseudo-scalar mediator models (figure 17b), similar sensitivity is obtained by all channels and mediator masses in the range 15-25 GeV are excluded. In all cases, a DM mass of 1 GeV is assumed, but the results are valid for all DM mass choices for which the mediator's decay into a pair of DM particles is kinematically allowed (m φ/a > 2m χ ). Pseudo-scalar mediator models can also be constrained by jet+E miss T final states, where the mediator is produced through loop-induced gluon fusion. Although the limits obtained by this signature are not competitive with the tt +E miss T final state, except in the mass range above 300 GeV, they provide a complementary constraint, which would become particularly important in case of a discovery. For the scalar model, the jet+E miss T final-states cross-section is instead too small to be probed. Ditop resonance searches in final states with two or four tops can also constrain this parameter space for m φ/a > 2m t . However, tt production through a spin-0 resonance presents a strong interference pattern with SM top pair production [208], which needs to be treated with care. On the other hand, four-top final states are characterised by relatively low event yields with the currently available integrated luminosity. For these reasons the study of these final states is not considered here. Finally, bb +E miss T final states are also used to set constraints on these simplified models, resulting in upper limits on the signal strength between 200 and 300 for mediator masses below 100 GeV. These results quantify the sensitivity to these models if up-type couplings are suppressed.

Colour-charged interaction
The strongest exclusion limits on colour-charged mediators η q that couple to first-and second-generations left-handed quarks are set by the jet+E miss T analysis. Assuming a unitary coupling, η q mediator masses up to 1.7 TeV are excluded for m χ = 50 GeV. Furthermore, η q mediator masses below 600 GeV are excluded for all DM masses such that the decay η q → qχ is kinematically allowed. The strongest exclusion limits on colour-charged mediators η b that couple to third-generation right-handed b-quarks are set by the b + E miss T analysis. Assuming a coupling set to the value that yields a relic density value consistent with astrophysical observations, masses up to 1.4 TeV are excluded for m χ = 1 GeV. Finally, t+E miss T final states are used to constrain the colour-charged mediator's η t coupling to right-handed top quarks. Mediator masses up to 3.4 TeV are excluded, assuming a 10 GeV DM particle mass and setting the coupling strengths of this model to: λ t = 0.4 and g s = 0.2.  figure 18. The statistical combination of the two analyses is also presented. Masses of the pseudo-scalar A in the range 200-600 GeV are excluded for m Z V = 1.5 TeV. The limit in sensitivity is driven by the fact that the A → χχ branching ratio decreases with increasing m A due to decay channels involving top quarks or other heavy bosons of the extended Higgs sector becoming accessible (tt, HZ and W ± H ∓ ). At higher m Z V the loss in branching ratio is combined with the smaller production cross-section so that the reach of the analysis is limited to smaller pseudo-scalar masses. For m A < 2m t and m A > 2m χ , there are no more competing decay channels and the reach of the analysis does not depend on m A any longer. This creates the turnover in the exclusion contour for m Z V = 2.5 TeV.
The two h+E miss T decay signatures are highly complementary at low Z V masses, as can be observed in the enlarged inset in the figure, while the h(bb) + E miss T analysis dominates the sensitivity at high Z V masses. Due to this complementarity, the gain obtained by the statistical combination of the two signatures is limited to the low mass region for this model.

Two-Higgs-doublet models with a pseudo-scalar mediator
As highlighted in section 2.3.2, the 2HDM+a model is characterised by a rich phenomenology. Constraints on this model from ATLAS searches are presented in this paper. Four different benchmark scenarios are used to evaluate the sensitivity to this model achieved by the Z/h+E miss T , tt/bb+E miss T , h(inv), and tttt analyses. These four benchmark scenarios [78] are consistent with bounds from electroweak precision, flavour and Higgs observables and are chosen to highlight the complementarity of the various final states. These scenarios represent two-dimensional and one-dimensional scans of a five-dimensional parameter space, used to present the exclusion limits.  Figure 19. Regions in the (a) (m a , m A ) and (b) (m a , tan β) planes excluded by data at 95% CL by X + E miss T and tttt analyses, following the parameter choices of scenarios 1 and 2 of the 2HDM+a model. The dashed grey regions at the top of (a) and the bottom of (b) indicate the region where the width of any of the Higgs bosons exceeds 20% of its mass. The exclusion limits presented above conservatively neglect the contribution from bb-initiated production, which might be sizeable for tan β ≥ 3 for the Z + E miss T channel and, to a lesser extent, for the h + E miss T one.

JHEP05(2019)142
In all cases, the masses of the heavy pseudo-scalar, heavy scalar, and charged bosons are kept equal (m A = m H = m H ± ). As visible in the results presented in figure 19a, the exclusion sensitivity is vastly dominated by the h(bb) + E miss In the context of 2HDM models, it is customary to investigate the sensitivity in terms of the tan β parameter. This is achieved in the second scenario presented in figure 19b. Although the exclusion reach is dominated also in this case by the h(bb) + E miss T and Z( ) + E miss T analyses, two additional signatures contribute at tan β ∼ 0.5: the tt +E miss T and the tttt analyses. The sensitivity of the former analysis is driven by the production cross-section of the a mediator in association with a top-quark pair and it decreases when the decay of the light pseudo-scalar into a top-quark pair is kinematically allowed and competes with a → χχ. On the other hand, the sensitivity of the latter analysis is fairly independent of m a due to the contribution to the total four-top production cross-section from the heavy bosons H/A → tt, both of which have masses fixed to 600 GeV in this scenario. In the case of the h(bb) + E miss T analysis, the exclusion was not investigated below tan β = 0.5. Given the non-trivial dependency of the width on tan β in this channel, it is not possible to extrapolate beyond the area explored. Figures 20a and 20b present the exclusion limits dependence on the mixing angle, sin θ, for a low-mass and high-mass a hypothesis, as evaluated in the third scenario. The limits are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section of the model. For scenario 3a (figure 20a), the lowest cross-section values are excluded by the Z( ) + E miss T and h(bb) + E miss T analyses. The sensitivity of both Z + E miss T analyses monotonically improves as a function of sin θ, as the cross-section of the non-resonant and resonant production diagrams, in figures 6d and 6e respectively, increases with sin θ. Conversely, the same production diagrams for the h + E miss T signatures have very different dependence on the mixing angle [152] in the two m a regimes explored here. The contribution of each diagram is also affected by the different h(bb)+E miss T and h(γγ)+E miss T analysis selections. For this scenario, both analyses have maximum of sensitivity around sin θ ∼ 0.5. The three heavy-flavour signatures, bb +E miss T , tt +E miss T and tttt, are presented for different tan β assumptions. A value of tan β = 50 is studied for bb +E miss T , with the aim of probing the parameter space where the coupling of the a mediator to down-type quarks is enhanced. However, the tt +E miss  [78,209,210] where the predicted relic density is depleted by the resonant enhancement of the processes χχ → A/a → SM. small mixing angles. Scenario 3b, presented in figure 20b, sets the mass of the light pseudoscalar so that the a → tt decay is kinematically allowed, which introduces an additional sin θ dependence to the X + E miss T analyses interpreted in this scenario. For this reason, the highest sensitivity for each analysis is found to be broadly around (or slightly below) the maximal mixing condition (θ = π/4), except for the tttt and h + E miss T signatures. The tttt signature shows a constant sensitivity as a function of sin θ (with an increase for very high values) due to the mass assumptions of this scenario (m a = 350 GeV and m A/H = 1 TeV) which cause the tttt production cross-section to be completely dominated by the tt + a(tt) process. The h + E miss T signatures have a complex dependence on the mixing angle. This is due to the different contributions of resonant and non-resonant processes to the final selection in the two analyses. In this case it is possible to observe that the h(bb) + E miss T analysis presents a maximum in sensitivity around the maximal mixing condition. The h(γγ) + E miss T analysis instead shows a local sensitivity minimum around sin θ ∼ 0.55. Finally, figure 21 presents the reach of the various experimental searches in a cosmological perspective, following the prescription of the fourth benchmark scenario. In this case, the observed exclusion limits in terms of the ratio of the excluded cross-section to the nominal cross-section of the model are investigated as a function of the DM mass, which is the parameter with the strongest impact on the relic density predicted by the 2HDM+a model.  [78,209,210] where the predicted relic density is depleted by the resonant enhancement of the processes χχ → A/a → SM. The plateau around and above m χ ∼ 200 GeV is determined by the increase in annihilation cross-section of the DM particles close to threshold for χχ → ha → SM and χχ → tt. For DM masses around m a /2 or m χ > 170 GeV the model predicts a relic density which is equal to or below the thermal value, Ωh 2 = 0.12. As the DM mass increases further, annihilation via single s-channel diagrams is more and more suppressed and the observed DM relic density can again be reproduced. At low values of m a this happens around m χ ∼ 10 TeV and is outside the range in figure 21. For all X + E miss T signatures considered, the sensitivity is independent of the DM mass as long as the lightest pseudo-scalar mediator, whose mass is fixed at 250 GeV in this scenario, is allowed to decay into a χχ pair. The Z( ) + E miss T analysis excludes this parameter space. For higher DM masses, the sensitivity of all analyses quickly decreases and no exclusion is observed. For m χ > m a /2 all parameter choices that fulfil or deplete the relic density value are still unconstrained.

Scalar dark energy model
The results of the jet+E miss  figure 22, it is assumed that the EFT approximation is valid for events where Q tr < g * M . For events failing this requirement, the iterative rescaling procedure detailed in ref. [75] is applied. The tt +E miss T analysis yields the most stringent constraints on the L 1 operator (figure 22a), as expected from the fact that the interaction described by L 1 is proportional to the masses of the SM fermions to which the DE scalar couples. The limits are obtained from the search channel (fully hadronic, semileptonic or fully leptonic top pair decays) that provides the smallest expected CL s value. The fully-hadronic and semileptonic channels contribute the most and similarly to the final sensitivity of the analysis, which excludes a suppression scale of about 200 GeV for g * π 2 . The tt +E miss T search is not yet sensitive to weakly coupled models, due to the high momentum transfers involved in the production of the top quarks, which are close to the exclusion limit.
The jet+E miss T analysis yields the most stringent constraints on the L 2 operator (figure 22b), due to the fact that this interaction is proportional to the momenta of the particles involved, excluding up to M 1.2 TeV for g * π. Due to the absence of heavy particles in the final state, the region of EFT validity for the jet+E miss T search is larger, with the constraints extending to lower values of the effective coupling.
These results improve upon the constraints on the disformal operator from astrophysical probes and non-collider experiments by several orders of magnitude [211] and also represent a significant improvement on the limits obtained by a similar reinterpretation of ATLAS and CMS results that made use of a smaller dataset at √ s = 8 TeV [82]. Finally, a Horndeski model for dark energy is studied in the context of ATLAS searches. This model introduces a dark energy scalar which couples to gravity. Limits on the two Lagrangian effective operators, L 1 and L 2 , are set by the tt +E miss T and jet+E miss T analyses, respectively. The suppression scale is excluded up to 200 GeV for g * = π 2 for the L 1 operator. For the L 2 operator, suppression scales up to 1.2 TeV for g * = π are excluded. These results are the first interpretation of a dark energy model by a collider experiment.
In this paper, many interpretations in the context of DM and DE models were added with respect to previous publications. This allowed to restrict very significantly the available parameter space of spin-0, spin-1 and 2HDM-based mediator-DM models as well as EFT DE models.

A Signal models generation details
The model implementations, settings and parameter scans used in this paper follow the prescriptions of the DM Forum/LHC DM Working Group [75][76][77][78]. and all generation settings used for signal models in this paper are summarised in tables 4, 5 and 6.

A.1 V/AV models
For all V/AV models, reconstructed samples were produced only for a specific reference scenario (either a vector or an axial-vector leptophobic mediator model). Rescaling factors for the acceptance (w A ) and the cross-section (w σ ) were calculated to match the acceptance and cross-section of each of the other scenarios to the reference. The acceptance weights were calculated for each (m Z , m χ ) mass hypothesis as the ratio of the particle-level acceptance for each of the NLO benchmark models considered (A NLO truth ) to the particle-level acceptance of the analysis for the reference NLO scenario in a fiducial region (A ref truth ): .
The cross-section weights were calculated for each (m Z , m χ ) mass hypothesis in a similar way, as the ratio of the reference cross-section at NLO to each cross-section of the four NLO benchmark models. The acceptance rescaling weights were found to be consistent with unity for the Z (χχ) + j and Z (χχ) + γ signatures. A few specific exceptions apply to this treatment. In case of the Z (χχ) + j signature, the cross-section rescaling factors were calculated from LO samples (DMSimp [116,212] generated with MG5 aMC@NLO 2.4.3 (LO) [213]) and applied to the samples described in table 4. In the specific case of the Z (χχ)+V signature, the baseline samples were generated at LO and rescaled at particle level to match the NLO samples described in the  the topcolour-assisted technicolour samples of [200] rescaled at particle level to match the DMSimp models described in table 4. The correction weights between the two samples were calculated from the bin-by-bin ratio of the invariant mass distributions of the tt system for the two samples at particle level. An additional uncertainty is assigned to account for this procedure as described in section 5.

A.2 VFC model
The VFC model is studied in two final states: pp → tZ VFC → tχχ and pp → tt(j) (via Z VFC ). A complete set of models with the full ATLAS detector simulation was generated as a function of m(Z VFC ) and assuming minimal width and unitary couplings, following the generation settings summarised in table 4. In order to assess the experimental constraints on all the model parameters, g ut , g χ and subsequently Γ(Z VFC ), a rescaling procedure is applied.
In the case of the tχχ final state, for each point in the parameter space, rescaling factors were calculated at particle level to match the acceptance and cross-section to those of the simulated reference model.
Three different matrix-element amplitudes contribute to the same-sign top-quark signature (tt(j)) relevant for this model: (i) prompt tt production in figure 2c, (ii) figure 2b, (iii) off-shell mediator in figure 2d. The relative contributions of the three amplitudes depend on the model parameters, but not the kinematic properties of each process. The three subprocesses, which were generated separately with full detector simulation, are combined according to the following formula to model the signal kinematics for any choice of parameter values (ξ): The functions α, β, γ can be computed with MadGraph as the ratio of the desired crosssection to the baseline cross-section:

A.3 2HDM + a models with heavy-flavour final states
The χχ + tt/bb signature of the 2HDM+a model can be successfully described [78] as the superposition of the associated production of two heavy-flavour quarks with either the light or the heavy pseudo-scalar mediator, which subsequently decays into a χχ pair. When the masses of the two pseudo-scalar mediators are sufficiently different, the contributions of the two processes can be factorised, and the 2HDM+a model can be described in terms of two sets of colour-neutral pseudo-scalar simplified models, each corresponding to the desired choice for m a and m A .

JHEP05(2019)142
The acceptance A of the analysis for each point of interest in the 2HDM+a parameter space is therefore derived as: where A simp is the acceptance of the analysis for the colour-neutral pseudo-scalar simplified model for a certain mass choice of the A(a)-boson, and σ a computed in fully reconstructed samples and σ A are the production cross-sections for pp → tta(→ χχ) and pp → ttA(→ χχ), respectively. 5 This rescaling is valid in the on-shell region, m a (m A ) > 2m χ [78].

B Comparison with direct and indirect searches
Searches for weakly-interactive massive particles (WIMPs) [57] represent the current paradigm for searches for particle dark matter (DM). Within this paradigm, understanding the nature of DM requires a measurement of its interaction cross-section with Standard Model particles. This can be achieved using three complementary methods [227], schematically depicted in figure 23 and briefly outlined in the following.
Direct searches. These searches aim to measure the elastic scattering of DM with nuclei in low background underground detectors such as CRESST-III [6], LUX [7], PICO [8], DEAP [9], PandaX [10], XENON [11,12] and SuperCDMS [13,14]. These direct detection experiments ultimately measure the strength of the interactions between WIMPs and the partons composing protons and neutrons and are sensitive to the properties of the DM halo around Earth.
Indirect searches. These searches aim to measure the annihilations or decays of DM particles in astrophysical systems, by means of neutrino detectors such as Su-perKamiokande [15] or IceCube [16] or by means of either ground or space telescopes, for example the H.E.S.S. Cherenkov telescope [17,18], AMS [228] and Fermi-LAT [19]. This measurement closely relates to the processes that determine the abundance of DM in the early universe and the interpretation of the results depends on the DM distribution in the universe as well as the SM particles into which the DM preferentially annihilates or decays.
Collider searches. These searches aim to discover DM particles and to measure the DM production cross-section through collisions of high-energy particles. The most stringent results to date on WIMPs are provided by the ATLAS [20-30], CMS [31,32,34,35,229] and LHCb [230,231] experiments at the LHC. Sub-GeV DM candidates are also constrained by the MiniBooNE experiment at Fermilab [232]. The interpretation of these results closely depends on the underlying mechanisms that couple DM to SM particles. In the simplified model framework considered in this paper this underlying mechanism is assumed to be the production of new mediator(s) state(s) which subsequently decay into DM. The present understanding of the DM puzzle is encompassed in the summary and comparison of the experimental results of these three approaches. Likewise, the discovery of DM as an elementary particle will require determination of its interaction cross-section with SM particles via each of these methods. It is convenient and customary to compare these different approaches in terms of spin-dependent (spin-independent) χ-nucleon scattering cross-sections as a function of the DM mass. In this paper, the ATLAS exclusion limits are converted into bounds on the χ-nucleon scattering cross-sections for the following models: • Vector and axial-vector neutral (V/AV) mediator models (two of the benchmark coupling scenarios, see section 6.1.1 for details).
For each model, the translation procedure to convert and compare these limits is well defined and described in ref. [76]. The interpretation in the spin-dependent (SD) and spin-independent (SI) DM-nucleon cross-sections, σ SD and σ SI , respectively, depends on the mediator mass and the couplings assumptions. Each comparison is valid solely in the context of the specific model and coupling assumptions. The ATLAS limits are always shown at 95% confidence level (CL) and the direct detection limits at 90% CL. ATLAS exclusion limits for pseudo-scalar colour-neutral mediator models should be compared with indirect search experiments in terms of the DM annihilation cross-section σv rel , as the rate in direct searches experiments is suppressed by additional velocitydependent terms entering the cross-section. However, the observed exclusion limits for pseudo-scalar mediator models with a unitary coupling assumption are limited to a very narrow mass range, due to a small data excess in the analysis (figure 17b). Therefore, this comparison is deferred to the results with the full Run-2 dataset.
The observed limits for the V/AV, VBC and scalar mediator models are compared with limits from direct search experiments in figures 24-26. The excluded regions are indicated by shaded areas inside the contours. Each combined contour summarises the ATLAS results for each considered model, obtained by using the best expected limit for -50 -

JHEP05(2019)142
each parameter point in the figure. When the contour does not close inside the plotted area, the exclusion of smaller scattering cross-sections does not imply that larger scattering cross-sections (beyond the vertical axis range) are also excluded.
The spin-dependent WIMP-neutron (WIMP-proton) scattering cross-section in the context of the leptophobic Z A mediator model is shown in the upper (lower) panel of figure 24. The difference between the WIMP-proton and WIMP-neutron cross-sections is negligible. The ATLAS exclusion curves are therefore identical in the two panels. The collider searches for this specific leptophobic axial-vector model complement the reach of the direct searches and extend beyond it, being particularly sensitive in the low-DM-mass parameter space, where the LUX and PICO experiments have less sensitivity due to the very low-energy recoils that such low-mass dark matter particles would induce. As in the case of the interpretation of the results in terms of mediator and DM masses (section 6.1.1), if the values chosen for the couplings are reduced, the relative interplay between direct and collider searches changes. This is exemplified by the change of lepton and quark couplings in the leptophilic Z A mediator model shown in figure 25, where the reach of resonant dijet final states is greatly reduced in favour of dilepton searches (differently for the two scenarios) and limited to mediator masses above 200 GeV. The sensitivity of the E miss T + X searches is the same for the two models in figures 24 and 25b. This is a coincidental result of two opposite effects [76]: the fact that the scattering cross-section limit is inversely proportional to the mediator mass reach (raised to the fourth power), which is higher in the leptophobic mediator model (figure 12a), and the fact that the σ SD limit is proportional to g 2 q , which is lower in the leptophilic mediator model. The spin-independent WIMP-nucleon scattering cross-section results for leptophobic, leptophilic, or baryon-charged vector mediator Z and scalar colour-neutral mediator φ are compared with the most stringent direct detection limits to date from the LUX, CRESST-II, XENON1T, SuperCDMS and PandaX experiments in figure 26. One contour for each model is presented in the figure and it includes the combination, based on the best expected limit for each parameter point, of every analysis considered for each model and presented in section 6. The excluded regions are indicated by shaded areas inside the contour. As before, when the contour does not close inside the plotted area, the exclusion of smaller scattering cross-sections does not imply that larger scattering cross-sections (beyond the vertical axis range) are also excluded. The collider searches in this case complement the reach of the direct searches for m χ 5 GeV. By comparing the exclusion reach of the ATLAS searches for each of the four models considered in figure 26, it is possible to gauge the importance of the production mechanism assumptions for the collider limits, which represent a complementary and not exclusive approach to DM searches with respect to direct and indirect searches.
-51 -  Figure 24. A comparison of the inferred limits with the constraints from direct detection experiments on (a) the spin-dependent WIMP-neutron or (b) WIMP-proton scattering cross-section in the context of the Z -like simplified model with axial-vector couplings. The results from this analysis, excluding the region to the left of the contour, are compared with limits from direct detection experiments. LHC limits are shown at 95% CL and direct detection limits at 90% CL. The comparison is valid solely in the context of this model, assuming a mediator width fixed by the dark matter mass and coupling values g q = 0.1, g = 0.1, and g χ = 1. LHC searches and direct detection experiments exclude the shaded areas. Exclusions of smaller scattering cross-sections do not imply that larger scattering cross-sections are also excluded. The resonance and E miss T +X exclusion region represents the union of exclusions from all analyses of that type.  Figure 25. A comparison of the inferred limits with the constraints from direct detection experiments on the spin-independent WIMP-nucleon (spin-dependent WIMP-neutron) scattering crosssection in the context of (a) the Z -like simplified model with leptophilic vector or (b) axial-vector couplings. The results from this analysis, excluding the region to the left of the contour, are compared with limits from the direct detection experiments. LHC limits are shown at 95% CL and direct detection limits at 90% CL. The comparison is valid solely in the context of this model, assuming a mediator width fixed by the dark matter mass and the coupling values highlighted in each figure. LHC searches and direct detection experiments exclude the shaded areas. Exclusions of smaller scattering cross-sections do not imply that larger scattering cross-sections are also excluded. The resonance and E miss T +X exclusion region represents the union of exclusions from all analyses of that type.  Figure 26. A comparison of the inferred limits with the constraints from direct detection experiments on the spin-independent WIMP-nucleon scattering cross-section. The results from ATLAS analyses, excluding the shaded regions, are compared with limits from direct detection experiments. LHC limits are shown at 95% CL and direct detection limits at 90% CL. The comparison is valid solely in the context of this model, assuming a mediator width fixed by the dark matter mass and coupling values g q = 0.25, g = 0 or g q = 0.1, g = 0.01 for the neutral-mediator model and coupling g q = 0.33 for the baryon-charged mediator. The coupling to the DM particle g χ , is set to unity in all cases. LHC searches and direct detection experiments exclude the shaded areas. Exclusions of smaller scattering cross-sections do not imply that larger scattering cross-sections are also excluded. The single dijet and E miss T +X exclusion region represents the union of exclusions from all analyses of that type.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

JHEP05(2019)142
[23] ATLAS collaboration, Search for dark matter produced in association with a Higgs boson decaying to bb using 36 fb −1 of pp collisions at √ s = 13 TeV with the ATLAS detector, -56 -