Search for dark matter in events with heavy quarks and missing transverse momentum in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp$$\end{document}pp collisions with the ATLAS detector

This article reports on a search for dark matter pair production in association with bottom or top quarks in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20.3 \mathrm {~fb}^{-1}$$\end{document}20.3fb-1 of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp$$\end{document}pp collisions collected at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s} = 8$$\end{document}s=8 TeV by the ATLAS detector at the LHC. Events with large missing transverse momentum are selected when produced in association with high-momentum jets of which one or more are identified as jets containing \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b$$\end{document}b-quarks. Final states with top quarks are selected by requiring a high jet multiplicity and in some cases a single lepton. The data are found to be consistent with the Standard Model expectations and limits are set on the mass scale of effective field theories that describe scalar and tensor interactions between dark matter and Standard Model particles. Limits on the dark-matter–nucleon cross-section for spin-independent and spin-dependent interactions are also provided. These limits are particularly strong for low-mass dark matter. Using a simplified model, constraints are set on the mass of dark matter and of a coloured mediator suitable to explain a possible signal of annihilating dark matter.


Introduction
The existence of dark matter (DM) in the Universe is highly motivated by many astrophysical and cosmological observations [1][2][3][4]. However, its nature remains a mystery. One of the best motivated candidates for a DM particle is a weakly interacting massive particle (WIMP) [5]. At the Large Hadron Collider (LHC), one can search for DM particles (χ ) that are pair produced in pp collisions. These studies are sensitive to low DM masses (m χ ≤ 10 GeV), and therefore provide information complementary to direct DM searches, which are most sensitive to larger DM masses [6][7][8][9].
If the particles that mediate the interactions between DM and Standard Model (SM) particles are too heavy to be produced directly in the experiment, their interactions can be e-mail: atlas.publications@cern.ch described by contact operators in the framework of an effective field theory [10][11][12]. For each operator considered, the reach is expressed in terms of the effective mass scale of the interaction, M * , and of the χ -nucleon cross-section, σ χ −N , as a function of m χ .
Since DM particles do not interact in the detector, the main signature of DM pair production at colliders is large missing transverse momentum. Initial-state radiation (ISR) of jets, photons, Z , or W bosons, was used to tag DM pair production at colliders in several searches at the Tevatron [13] and the LHC [14][15][16][17][18][19][20][21][22].
A new search for DM pair production in association with one b-quark or a pair of heavy quarks (b or t) was proposed in Ref. [23]. The dominant Feynman diagrams for these processes are shown in Fig. 1. To search for these processes, dedicated selections are defined to reconstruct the various production and decay modes of these heavy-quark final states. For final states containing a semileptonic decay of a top quark, the results of the search for a supersymmetric partner of the top quark are used [24].
The analysis presented in this article is particularly sensitive to effective scalar interactions between DM and quarks described by the operator [12] O scalar = q m q M N * q qχχ, (1) where N = 3 for Dirac DM (D1 operator) and N = 2 for complex scalar DM (C1 operator). The quark and DM fields are denoted by q and χ , respectively. The scalar operators are normalized by m q , which mitigates contributions to flavour-changing processes, strongly constrained by flavour physics observables [25,26], through the framework of minimal flavour violation (MFV). The dependence on the quark mass makes final states with bottom and top quarks the most sensitive to these operators. This search is also sensitive to tensor couplings between DM and quarks. The tensor operator (D9), which describes a magnetic moment coupling, is parameterized as [12]: MFV suggests that the D9 operator should have a mass dependence from Yukawa couplings although canonically this is not parametrised as such.
The results are also interpreted in light of a bottom-Flavoured Dark Matter model (b-FDM) [27]. The b-FDM model was proposed to explain the excess of gamma rays from the galactic centre, recently observed by the Fermi Gamma-ray Space Telescope, and interpreted as a signal for DM annihilation [28]. This analysis of the data recorded by the Fermi-LAT collaboration favours DM with a mass of approximately 35 GeV annihilating into b-quarks via a coloured mediator. In this model, a new scalar field, φ, mediates the interactions between DM and quarks as shown in Fig. 2. DM is assumed to be a Dirac fermion that couples to right-handed, down-type quarks. The lightest DM particle, which constitutes cosmic DM, preferentially couples to b-quarks. The collider signature of this model is b-quarks produced in association with missing transverse momentum. This analysis sets constraints on the mass of the mediator and DM particle in the framework of the b-FDM model.

Detector description and physics objects
The ATLAS detector [34] at the LHC covers the pseudorapidity 1 range of |η| < 4.9 and is hermetic in azimuth φ. It con-1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector, and sists of an inner tracking detector surrounded by a superconducting solenoid, electromagnetic and hadronic calorimeters, and an external muon spectrometer incorporating large superconducting toroidal magnets. A three-level trigger system is used to select events for subsequent offline analysis. The data set used in this analysis consists of 20.3 fb −1 of pp collision data recorded at a centre-of-mass energy of √ s = 8 TeV with stable beam conditions [35] during the 2012 LHC run. All subsystems listed above were required to be operational.
This analysis requires the reconstruction of muons, electrons, jets, and missing transverse momentum. Muon candidates are identified from tracks that are well reconstructed inside both the inner detector and the muon spectrometer [36]. To reject cosmic-ray muons, muon candidates are required to be consistent with production at the primary vertex, defined as the vertex with the highest ( p track T ) 2 , where p track T refers to the transverse momentum of each track. Electrons are identified as tracks that are matched to a well-reconstructed cluster in the electromagnetic calorimeter. Electron candidates must satisfy the tight electron shower shape and track selection criteria of Ref. [37]. Both electrons and muons are required to have transverse momenta p T > 20 GeV and |η| < 2.5. Potential ambiguities between overlapping candidate objects are resolved based on their angular separation. If an electron candidate and a jet overlap within R < 0.2, then the object is considered to be an electron and the jet is discarded. If an electron candidate and any jet overlap within 0.2 < R < 0.4, or if an electron candidate and a b-tagged jet overlap within R < 0.2 of each other, then the electron is discarded and the jet is retained.
Jet candidates are reconstructed using the anti-k t clustering algorithm [39] with a radius parameter of 0.4. The inputs to this algorithm are three-dimensional topological clusters [40]. The four-momentum of the jet is defined as the vector sum of the four-momenta of the topological clusters, assuming that each cluster originates from a particle defined to be massless and to come from the interaction point.
To calibrate the reconstructed energy, jets are corrected for the effects of calorimeter response and inhomogeneities using energy-and η-dependent calibration factors based on simulation and validated with extensive test-beam and collision-data studies [40]. In the simulation, this procedure calibrates the jet energies to those of the corresponding jets constructed from stable simulated particles. In-situ measure-Footnote 1 continued the z-axis along the beam line. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam line. The pseudorapidity η is defined in terms of the polar angle θ as η = − ln tan(θ/2). Observables labeled "transverse" are projected into the x-y plane. ments are used to further correct the data to match the energy scale in simulated events. Effects due to additional pp interactions in the same and preceding bunch crossings (pile-up effects) are corrected [41]. Only jets with p T > 20(25) GeV and |η| < 4.5(2.5) are considered in this analysis for final states involving b (t) quarks.
Jets containing particles from the hadronisation of a bquark (b-jets) are tagged using a multivariate algorithm [42,43]. The b-tagging algorithm combines the measurement of several quantities distinguishing heavy quarks from light quarks based on their longer lifetime and heavier mass. These quantities include the distance of closest approach of tracks in the jet to the primary event vertex, the number and position of secondary vertices formed by tracks within the jet, as well as the invariant mass associated with such vertices. The algorithm is trained on Monte Carlo (MC) simulations and its performance is calibrated using data. To optimize the sensitivity of this analysis, a requirement on the output of the b-tagging algorithm which provides a 60 % (70 %) b-jet efficiency operating point is used in signal regions (SR) 1 and 2 (3 and 4) defined below. The corresponding misidentification probability is 15 % (20 %) for c-jets, and less than 1 % for light-quark jets. The aforementioned b-tagging efficiencies and misidentification probabilities were derived in a simulated tt sample with jet transverse momenta of p T > 20 GeV and |η| < 2.5.
The missing transverse momentum, with magnitude E miss T , is defined as the negative vector sum of the transverse momenta of jets, muons, electrons, photons, and topological clusters not assigned to any reconstructed objects [44].

Event selection
Candidate signal events containing at least one high-p T jet and large E miss T are assigned to one of four orthogonal signal regions. The first two signal regions focus on events with DM produced in conjunction with one (SR1) or two (SR2) b-quarks in the final state. SR3 and SR4 target events in which DM is produced in conjunction with a tt pair, where either both top quarks decay hadronically (SR3) or one top quark decays hadronically and the other semileptonically (SR4). SR4 was developed for a top squark search by the ATLAS Collaboration and coincides with the "tNbC_mix" signal region described in Ref. [24]. The four signal regions provide the complementary information needed in case of observation of a signal.
Events assigned to SR1 and SR2 are required to pass a calorimeter-based E miss T trigger with a threshold of 80 GeV. To enrich the sample in pp → χχ +b(b), events are required to have a low jet multiplicity (n jets <5), E miss T > 300 GeV, and the most energetic b-tagged jet must have a p T > 100 GeV. The azimuthal separation between the directions of the jets and the missing transverse momentum is required to be more than 1.0 radian. Events with at least one identified muon or electron are discarded to reject leptonic decays of W and Z bosons. Events satisfying these selection criteria are assigned to SR1 provided that the jet multiplicity does not exceed two. Events are assigned to SR2 when at least three jets are reconstructed in the event and the second most energetic jet has p T > 100 GeV. If there is a second b-tagged jet it has to satisfy p T > 60 GeV.
Events assigned to SR3 are required to pass triggers specifically designed to select hadronic decays of top quark pairs. Such triggers require either five jets with p T ≥ 55 GeV each or four jets with p T ≥ 45 GeV, of which one is tagged as a b-jet. To select pp → χχ + tt events, at least five reconstructed jets are required, of which at least two are b-tagged, and E miss T > 200 GeV. Furthermore, the azimuthal separation between the most energetic b-jet and the missing transverse momentum is required to be at least 1.6 radians. To reduce W/Z leptonic decays and leptonic top quark decays, events with at least one identified muon or electron are discarded. To maximize the rejection of the abundant tt background, the Razor variable R [33] is used. This variable utilizes both transverse and longitudinal information about the event to fully exploit the kinematics of the decay. To separate signal and background, R > 0.75 is required.
To enrich the sample in pp → χχ + tt with one semileptonic decay of the t quark, events assigned to SR4 use singlelepton or E miss T triggers, and require exactly one isolated lepton (electron or muon) with p T > 25 GeV, at least four high-p T jets, where one jet is b-tagged with p T > 60 GeV. Events with E miss T > 270 GeV are selected when the transverse mass 2 formed by the lepton and E miss T , m T ( , E miss T ), exceeds 130 GeV and E miss where the jets are ordered by decreasing p T . The azimuthal angle between the missing transverse momentum and the two most energetic jets is required to be greater than 0.6 radians.
Special variables, such as the asymmetric transverse mass am T2 [29][30][31] and the topness variable [32], are used to reject the dileptonic tt component of the background. Details can be found in Ref. [24]. The diboson background is suppressed by a requirement on the three-jet invariant mass (m j j j < 360 GeV) [24]. A τ veto rejects tt events with hadronically decaying τ leptons in the final state. Addi- 2 Since the longitudinal component of the momentum of the neutrinos is not measured, the measured properties of the W boson candidates are limited to their transverse momentum and transverse mass, defined as is the magnitude of the missing transverse momentum vector, p T is the transverse momentum of the lepton and p x and p y (E miss x and E miss y ) are the magnitude of the x and y components of the lepton momentum (missing transverse momentum) respectively.
represent the transverse momentum of the ith jet (b-tagged jet). The asymmetric transverse mass am T2 [29][30][31], topness [32], m j j j and Razor R [33] are used to reject the abundant top quark background on the angles between the lepton and the various jets are imposed to further reduce the tt background. Table 1 provides an overview of the selections applied in all four signal regions. The product of the detector acceptance A and the reconstruction efficiency for the selections described above varies between 0.1 and 8 % depending on the signal region, operator, and specific channel considered. SR1 and SR2 have the highest efficiencies (A × > 2 %) for the D9 operator, while SR3 and SR4 are most efficient for the D1 and C1 operators (A × > 1 %).
The dominant background for SR1 and SR2 is due to Z → νν events produced in conjunction with one or more jets. This irreducible background is estimated from data using two control regions (CRs). The first CR exploits Z +jets events with Z → μ + μ − , while the second uses γ +jets events for which the production at high transverse momentum ( p γ T > M Z ) mimics that of Z +jets [45]. The γ +jets control region substantially increases the number of events at large missing transverse momentum. The transverse momentum of the dimuon pair or photon is added vectorially to the E miss T of the event to simulate the Z → νν background. Corrections to compensate for the differences in efficiency and acceptance between the Z (νν)+jets and Z (μ + μ − )+jets or γ +jets are derived from data using control regions without b-tagged jets before applying any requirements on the missing transverse momentum. Remaining kinematic selections correspond to the ones described in Table 1. A muon control region is chosen because the energy loss of muons in the detector is comparatively small. The systematic uncertainties introduced by this data-driven procedure on the Z (νν)+jets background are approximately 10 %, mainly from the flavour composition of background processes, kinematic differences between the control and signal regions and relative normalizations of backgrounds.
Production of W/Z +jets with subsequent leptonic decays of W and to a much smaller degree Z is also a substantial source of background for SR1 and SR2 when the resulting charged leptons fail to be identified or if the W or Z bosons decay to τ leptons. These contributions are estimated from Z ( + − )+jets and W ( ν)+jets MC samples generated using ALPGEN2.3 [46] with the CTEQ6L1 [47] parton distribution function (PDF) set. The procedure used for the normalization of this sample is described in reference [48]. These samples are generated with up to five light partons (u, d, s) and one c quark or two heavy quarks (c, b) per event. W + b production is highly suppressed and therefore negligible. A control region enriched in W ( ν)+jets events is selected by adding a lepton requirement to the selection and is used to validate the estimate of this background. The purity of W ( ν)+jets in the control region for SR1 (SR2) is 67 % (47 %). After full selection the contribution of b(c)-quarks to the dominant W ( ν)+jets background is approximately 39 % (38 %) for SR1 and 52 % (37 %) for SR2. The systematic uncertainty on this background is approximately 20 %. Finally, the small contribution from tt is estimated using MC samples and validated in data control regions before applying signal selection requirements. The tt process is selected with very high purity Table 2 Expected background and signal yields for m χ = 10 GeV compared with observed yields in data for the various signal regions. For the b-FDM model, m φ is 600 GeV. The row labeled "total expected background" shows the sum of all background components. The quoted uncertainties include all statistical and systematic effects added in quadrature. The effective mass scale, M * , is set to be 100/40/600 GeV for the D1/C1/D9 operators, approximately corresponding to the expected limit. The probabilities of the background-only hypothesis, p values, are also given. by requiring events with one lepton and large jet multiplicities.
The dominant source of background for SR3 and SR4 is tt events. In SR3, this contribution is estimated from data using a control region not overlapping with SR4 and largely dominated by tt events with one of the two top quarks decaying semileptonically. The five-jets requirement is relaxed to three jets. Additionally, the event is required to contain exactly one lepton with p e(μ) T > 30 (25) GeV and must fulfill E miss T + m T > 25 (30) GeV for the electron (muon) channel. The potential signal contribution to this selection is less than 0.1 %. The uncertainties are small because the SR3 data control region uses a kinematic region similar to the signal region with the lepton veto and jet multiplicity being the main difference. These effects were studied and considered as systematic uncertainties. Dominant uncertainties are related to jets and the top quark momentum distribution. Corrections to compensate for the differences in efficiency and acceptance between hadronic and semi-leptonic top decays are derived from MC samples generated using the POWHEG BOX generator [49] interfaced with JIMMY4.31 [50] with the nextto-leading-order (NLO) PDF set CT10 [51]. The systematic uncertainty on the tt background in SR3 of approximately 7 % is derived by studying corrections for the top quark momentum distribution, and shower modelling by interfacing the same generator with PYTHIA6 [52,53].
In SR4, the tt background is estimated from data using a control region obtained by requiring 60 GeV < m T < 90 GeV and loosening the selection criteria on E miss T , am T2 , and E miss T / H 4 j T . A similar selection, but applying an inverted b-tagging requirement, is used to estimate the W ( ν)+jets background. The uncertainty on the tt background is estimated to be approximately 20 % [24], which is larger than the uncertainty in SR3 due to the limited statistics. These uncertainties are evaluated by varying the renormalisation and factorisation scale of the simulations, comparing alternative PDF sets, and studying the effects of different shower generators and of ISR and final-state radiation.
Additional sources of background, which include singletop, tt + Z /W , and diboson production, are estimated in all signal regions using simulations and NLO cross sections [54,55]. The single-top (s-channel) and W t background is generated using the POWHEG generator. The single-top tchannel is generated with ACERMC3.8 [56] interfaced with PYTHIA6. Associated production of tt and a vector boson (W , Z ) are generated with MADGRAPH5 [57] with up to two additional partons interfaced with PYTHIA6. The crosssections for tt production in association with a W (Z ) boson are determined using the MSTW2008 NLO (CTEQ6.6M) PDF sets. The diboson samples are generated using HER-WIG6.520 [58,59] and JIMMY4.31 with the CTEQ6L1 PDF set. The multijet background is estimated using data-driven methods [60] and is found to be negligible in all signal regions after full selection.
Object reconstruction efficiencies in simulated events are corrected to reproduce the performance measured in data. The systematic uncertainty of the background estimates derived from simulation combines the uncertainties on the efficiency of the b-tagging algorithm, the uncertainties on the determination of the energy scale and resolution of the jet energy and E miss T , the theoretical uncertainty on the various cross-sections, changes in the shapes of distributions used to extrapolate event counts from control regions to the signal region, data driven corrections and the PDF uncertainties. Overall, the systematic uncertainty on the background estimated from simulation is calculated to be between 12 and 18 %, depending on the signal region.
The simulation of the signal samples of pp → χ χ +b(b), pp → χ χ + tt, and b-FDM employs the MADGRAPH5 generator interfaced with PYTHIA6 using the CTEQ6L1 PDF. Samples are generated for operators D1, C1, and D9, assuming M * = 1 TeV and m χ between 10 and 1300 GeV. Samples for the b-FDM model are generated for m χ values between 1 and 1300 GeV and mediator masses, m φ , between 5 and 3000 GeV. The instrumental uncertainties on the simulated signal yields for D1, C1, and D9 operators are between 11 and 15 %, depending on the signal region. The equivalent uncertainties for the b-FDM model range between 6 and 16 % depending on m χ and the mediator mass. The uncertainties from the PDF are computed by comparing the rates obtained with the default PDF set (CTEQ6L1) with those obtained with two alternative sets (MSTW2008LO and NNPDF21LO [61,62]). The uncertainties on the signal acceptance from PDF and scale variations are estimated to be approximately 10 % for the D1, C1, and D9 operators for m χ = 10 GeV and approximately 6 % for b-FDM models.
The validity of the effective field theory assumption depends on the momentum transfer of the process modelled, which should be below the energy scale of the underlying interactions [63]. To account for this, the momentum transfer m(χχ) = Q tr in the events is required to be less than the energy scale probed. Specifically, Q tr must be smaller than the mass M of the heavy mediator. For an ultraviolet completion this implies M * = M/ √ g q g χ . Along with perturbativity of the couplings g q g χ < 4π this leads to the following validity requirements on MC truth level: Q tr < 4π(M 3 * /m q ) 1/2 (D1), Q tr < 4π M * (D9), Q tr < (4π) 2 M 2 * /m q (C1). Table 2 shows the expected background from various sources in the four signal regions as well as the observed yields in data. The expected signal yields for the operators D1, C1, and D9, as well as for the b-FDM model are also shown. The probabilities of the background-only hypothesis, p values, for the signal regions SR1, SR2, SR3, and SR4 are 0.09, 0.29, 0.24, and 0.18, respectively. As no significant excess is observed, limits on the signal yield are set using a profile likelihood ratio test following the C L s prescription [64]. Also given is the 95 % confidence level ( [6,7,65]. The coupling is assumed to be g q g χ = g = 4π Figure 4 shows the 90 % CL exclusion curves for the effective mass scale M * as a function of m χ . The results for the operators D1, C1, and D9 are presented individually for all four signal regions. The best limits on the D1 and C1 operators are obtained using SR4, while SR1 provides the best limits on the D9 operator, as shown in Fig. 4. These limits are then converted into limits on the χ -nucleon cross-section [12]. Figures 5 and 6 show the corresponding 90 % CL exclusion curves for the spin-independent and spindependent χ -nucleon cross-section for the scalar (D1) and tensor (D9) operators as a function of m χ for the strongest results obtained in any signal region. The most stringent limits set by direct detection experiments [6][7][8][9] are also shown. Only m χ where more then 90 % of the events fulfill the effec-   Fig. 6 Upper limits at 90 % CL on the spin-dependent χ-nucleon cross-section (σ SD χ−N ) for the tensor operator D9 (red) as a function of m χ . The yellow and green curves represent the exclusion limits recently set by the COUPP and PICASSO collaborations [8,9,65]. The coupling is assumed to be g q g χ = g = 4π The limits shown are especially strong in the low-mass region where several collaborations [28,[66][67][68] have recently claimed possible observations of DM. The results reported in this article represent the first ATLAS limits on the scalar operator C1 and they significantly improve the sensitivity to χ -nucleon interactions mediated by the scalar operator D1 compared to previous ATLAS results [14,16,18,19]. Figure 7 shows the exclusion curves observed and expected for the b-FDM model as a function of the mediator and DM masses. For each point in (m χ , m φ ), the signal region with the best expected sensitivity is used, with SR1 dominating over the other signal regions. For a DM particle of approximately 35 GeV, as suggested by the interpretation of data recorded by the Fermi-LAT collaboration, mediator masses between approximately 300 and 500 GeV are excluded at 95 % CL.

Conclusions
In summary, this article reports a search for dark-matter pair production in association with bottom or top quarks. The analysis is performed using 20.3 fb −1 of pp collisions collected at √ s = 8 TeV by the ATLAS detector at the LHC. The results are interpreted in the framework of an effective field theory to set stringent limits on scalar and tensor interactions between Standard Model and DM particles. The data are found to be consistent with the Standard Model expectations, and limits are set on the mass scale of effective field theories that describe scalar and tensor interactions between DM and Standard Model particles. The exclusion limits are strongest at low DM masses. The limit on the χ -nucleon cross-section mediated by the D1 operator is improved significantly with respect to previously published ATLAS results by obtaining sensitivities of approximately σ SI χ −N = 10 −42 cm 2 for m χ = 10 GeV. Constraints on b-Flavoured Dark Matter models, suitable to explain a possible signal of annihilating DM, are also presented. The excluded regions depend on m χ and m φ . For m χ = 35 GeV, mediator particles with m φ = 300-500 GeV are excluded.
Acknowledgments We would like to thank Tongyan Lin (University of Chicago) for helpful discussions about the models presented and the interplay between collider DM constraints and direct and indirect DM experiments. We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4.0.