Simplified dark matter models with a spin-2 mediator at the LHC

We consider simplified dark matter models where a dark matter candidate couples to the standard model (SM) particles via an $s$-channel spin-2 mediator, and study constraints on the model parameter space from the current LHC data. Our focus lies on the complementarity among different searches, in particular monojet and multijet plus missing energy searches and resonance searches. For universal couplings of the mediator to SM particles, missing-energy searches can give stronger constraints than $WW$, $ZZ$, dijet, dihiggs, $t\bar t$, $b\bar b$ resonance searches in the low-mass region and/or when the coupling of the mediator to dark matter is much larger than its couplings to SM particles. The strongest constraints however come from diphoton and dilepton resonance searches. Only if these modes are suppressed, missing-energy searches can be competitive in constraining dark matter models with a spin-2 mediator.


Introduction
Convincing astrophysical and cosmological observations for the existence of dark matter (DM) provide us one of the strong motivations to consider physics beyond the standard model (SM). The search for DM is thus one of the main pillars of the LHC physics program.
As the nature of DM is known so little, a so-called simplified-model approach [1] has been widely adopted, and concrete simplified DM models have recently been proposed by the LHC DM working group to conduct the systematic DM searches at the LHC Run-II [2]. Following the proposal, the Run-I data as well as the early Run-II data have already been analysed to constrain simplified DM models with s-channel spin-1 and spin-0 mediators, see e.g. [3][4][5][6][7][8][9][10][11][12]. On the other hand, the model with a spin-2 mediator [13,14] has not been fully explored for the LHC yet-it is one of the next-generation simplified DM models [15].
In this article, we consider simplified DM models where a DM candidate couples to the SM particles via an schannel spin-2 mediator, and study constraints on the model parameter space from searches in final states with and without missing energy in the current LHC data. This work follows the DMsimp framework [16][17][18], which provides the DM model files for event generators such as MadGraph5 aMC@NLO [19] as well as for DM tools such as micrOMEGAs [20][21][22] and MadDM [23,24]. The same framework was used previously to study the cases of s-channel spin-1 and spin-0 mediators.
We note that, to keep the analysis of the LHC constraints fully general, we do not impose any astrophysical constraints like relic density or (in)direct detection limits on the DM candidate, as these partly depend on astrophysical assumptions. Moreover, in a full model, the DM may couple to other new particles that are irrelevant for the collider phenomenology discussed here. We refer readers to [13,14] for the astrophysical constraints, and to [25] for a discussion of spectral features in indirect detection.
The article is organised as follows. The simplified model is presented in Section 2, and the production and decays of the spin-2 mediator in Section 3. The re-interpretation of the LHC results is discussed in Section 4. Section 5 contains a summary and conclusions. Supplemental material for recasting is provided in the Appendix.
SM particles via an s-channel spin-2 mediator. The interaction Lagrangian of a spin-2 mediator (Y 2 ) with DM (X) is given by [13] L Y2 where Λ is the scale parameter of the theory, g T X is the coupling parameter, and T X µν is the energy-momentum tensor of a DM field. Here, we consider three types of DM independently; a real scalar (X R ), a Dirac fermion (X D ), and a vector (X V ). The interaction with SM particles is obtained by where i denotes each SM field, i.e. the Higgs doublet (H), quarks (q), leptons ( ), and SU (3) C , SU (2) L and U (1) Y gauge bosons (g, W, B). Following [26,27] we introduce the phenomenological coupling parameters without assuming any UV model. 1 The energy-momentum tensors of the DM are where F µν is the field strength tensor. Those of the SM fields are similar; see e.g. [28] for the explicit formulae.
Complying with the simplified-model idea, it is instructive to consider universal couplings between the spin-2 mediator and the SM particles: With this simplification, the model has only four independent parameters, two masses and two couplings: where we dropped the superscript T for simplicity. Such a universal coupling to SM particles is realised, e.g., in the 1 One may also assign independent coupling parameters for each flavour, especially for heavy flavours [28].  original Randall-Sundrum (RS) model of localised gravity [29]. The parameters are related as where x 1 = 3.83 is the first root of the Bessel function of the first kind, k is the curvature of the warped extra dimension, and M Pl = 2.4×10 18 GeV is the reduced fourdimensional Planck scale. On the other hand, in the socalled bulk RS model [30,31], where the SM particles also propagate in the extra dimension, g T i can take different values depending on the setup.
In [28], the SM sector of the above model was implemented in FeynRules/NloCT [32,33] (based on [34][35][36]), and the Y 2 production and decay rates at next-toleading order (NLO) QCD accuracy were presented. In this work, we include the three DM species (X R , X D , X V ) with the corresponding interactions, and add the model into the DMsimp framework [37] as the simplified DM model with a spin-2 mediator. Regarding LHC phenomenology, let us begin by discussing the spin-2 mediator decays. The partial widths for the decays into a pair of spin-0 (S = X R , h), spin-1/2 (F = X D , q, ) and spin-1 (V = X V , g, γ, Z, W ) DM or SM particles are given by where β i = √ 1 − 4r i with r i = m 2 i /m 2 Y , g γ = g B cos 2 θ W + g W sin 2 θ W and g Z = g B sin 2 θ W + g W cos 2 θ W with the weak-mixing angle θ W , and f (r V ) = 1 + 1 12 κ 2 For gluons and photons, κ H = 0 in f (r V ), while κ H = 1 for vector DM. The factors N ν = 1/2 for neutrinos and N s = 1/2 for two identical particles, and are unity otherwise; N F,V C is the number of colours. We note that B(Y 2 → Zγ) = 0 for g W = g B as the decay rate is proportional to We see that, due to the different overall prefactors, the partial widths become larger in order of scalar, fermion, vector DM. Moreover, the different powers (5, 3, 1) of the velocity factor β i indicate that the decay proceeds mainly via a D, P, and S wave for the scalar, fermion, and vector case, respectively. Figure 1 shows the Y 2 total width scaled by the mass, Γ Y /m Y , and the decay branching ratios for the case that only decays into SM particles are allowed. MadWidth [38] provides the partial decay rates numerically for each parameter point. In Table 1 we provide the explicit values for a few representative mass points. We see that, for a universal coupling g SM , decays into gluons and light quarks, leading to a dijet signature, are completely dom-inant ( 80% depending on m Y ). The diphoton channel has 4-5% branching ratio; other diboson channels (W W and ZZ) as well as tt are important as well when kinematically allowed. Finally, it is important to note that decays into neutrinos have 3-4% branching ratio, leading to missing energy signatures independent of decays to DM. 2 The width is proportional to m 3 Y , and from the upper panel in Fig. 1 we see that for g SM /Λ (3 TeV) −1 , the resonance is always very narrow (Γ Y /m Y < 1%) up to m Y ∼ 1 TeV. Note here, that Λ is simply a scale parameter, not a physical cut-off of the theory.
When decays into DM are allowed, their relative importance depends on g X and the type of DM (scalar, Dirac or vector) as illustrated in Fig. 2; see also Eqs. (10)-(12). Two mass scales are considered: m Y = 100 GeV and 1 TeV, with m X = 10 GeV and g SM = 1. 3 We see that decays into DM can be important and even dominant, but the resonance remains narrow for any choice of Λ 3 TeV for m Y 1 TeV. Another important observation is that for scalar DM (X R ), for g X ∼ g SM the decay into Y 2 → X R X R is practically irrelevant; one needs g X /g SM ≈ 3 for the decay into DM to exceed the one into neutrinos, and g X /g SM ≈ 5-6 to reach the 10% level. For Dirac (X D ) and vector (X V ) DM, the decays into DM and into neutrinos are of comparable magnitude at g X ∼ g SM , both contributing to missing-energy signatures. For attains about 10% (20%). These differences depending on the type of DM will be important later for the collider limits.

Production of the spin-2 mediator
Turning to the production modes, the potentially interesting channels are inclusive Y 2 production (pp → Y 2 ), as well as the production with an extra hard tagging jet (pp → Y 2 j) or an electroweak boson (e.g. pp → Y 2 γ). With the Y 2 decaying into SM particles, the former gives resonant peak signatures (without missing energy). On the other hand, the latter two give the typical monojet or monophoton signatures when the mediator decays invisibly. Moreover, the latter two play a role in the low-mass resonance search in dijet events with initial-state radiation (ISR) as seen later.
The Y 2 production cross sections at NLO QCD accuracy for pp collisions at 13 TeV are depicted in Fig. 3 as a function of the mediator mass. 4 We employ Mad-Graph5 aMC@NLO [19]   the sum of the transverse masses of the final states as a dynamical scale choice. In our simplified model, the cross sections depend solely on g SM /Λ and scale with (g SM /Λ) 2 . The dashed lines showing g SM /Λ = (10 TeV) −1 are therefore an order of magnitude below the corresponding solid lines for g SM /Λ = (3 TeV) −1 . Also noteworthy is the fact that pp → Y 2 is mostly gluon-initiated for the low-mass case [39]; 97%, 83%, and 28% of the LO total rate for m Y = 100 GeV, 1 TeV, and 5 TeV, respectively, stem from gg fusion. Since the radiation of an initial-state photon (Z/W ) can only occur in the quark-initiated process, Y 2 + photon (Z/W ) production is very much suppressed as compared to Y 2 + jet production. This is also the reason that the process has a huge K factor especially in the low-mass region [28]. 5 In the context of DM searches, the monojet signature is expected to give important constraints on the model. The fiducial cross sections for pp → Y 2 j with p j T > 200 GeV and |η j | < 5 are shown in Fig. 3, where one can estimate the monojet cross section by taking into account the Y 2 branching ratio into DM particles (and/or neutrinos) when m Y > 2m X . In Fig. 4 we also plot the fiducial cross sections for pp → j + / E T as a function of the DM mass,  separating the contributions from neutrinos (black lines) and DM (red lines) produced through the spin-2 mediator. For definiteness, we take m Y = 500 GeV, Λ = 3 TeV, g SM = 1 and compare g X = 1, 2 and 10 for Dirac DM. As already seen in Fig. 2, their relative importance depends on g X . For m Y < 2m X , a pair of DM is produced via the off-shell mediator and the cross section is strongly suppressed. Therefore, the neutrino contribution always dominates the monojet signature for the m Y < 2m X region even if g X /g SM = 10. For the other DM types, scalar and vector, the picture is similar, but the relative importance to the neutrino channel is different; see Fig. 2. This is one of the characteristic features of the spin-2 mediator DM model with universal couplings, as compared to the s-channel spin-1 and spin-0 models, whose mediators do not couple to charged leptons and neutrinos in the minimal setup [2].

Constraints from current LHC data 4.1 Searches with missing energy
The ATLAS and CMS experiments have been searching for new physics in a large variety of final states. As mentioned above, in the context of DM searches, the monojet signature is regarded as particularly interesting. In practice, at 13 TeV, the monojet analyses require one hard jet recoiling against / E T , but allow for additional jets from QCD radiation. Therefore one can expect that multijet+ / E T searches are also relevant [41,42]. To work out the current constraints on the spin-2 mediator DM model from these searches, we consider the following early Run-II analyses: In the monojet analysis [5], a simplified DM model with an s-channel spin-1 mediator is considered. Events are required to have at least one hard jet with p T > 250 GeV and |η| < 2.4, and a maximum of four jets with p T > 30 GeV and |η| < 2.8 are allowed. Several inclusive and exclusive signal regions (SRs) are considered with increasing / E T requirements from 250 GeV to 700 GeV. The multijet+ / E T analysis [43] is designed to search for squarks and gluinos in supersymmetric models, where neutralinos lead to missing energy. Several SRs are characterised by minimum jet multiplicity from two to six; / E T > 200 GeV is required for all SRs, while different thresholds are applied on jet momenta and on the azimuthal separation between jets and / E T . To reinterpret the above analyses in the context of our spin-2 mediator simplified DM model, we use Check-MATE2 [44], which is a public recasting tool providing confidence limits from simulated signal events and includes a number of 13 TeV analyses. We generate hadronlevel signal samples by using the tree-level matrix-element plus parton-shower (ME+PS) merging procedure. In practice, we make use of the shower-k T scheme [45], implemented in MadGraph5 aMC@NLO [19] with Pythia6 [46], and generate signal events with parton multiplicity from one to two partons. We impose / E T > 200 GeV and set Q cut = 200 GeV for the merging separation parameter at the parton level; these values are chosen for an efficient event generation without affecting the final results. The event rate is normalised to the pp → Y 2 j NLO cross sections shown in Fig. 3. (Note, however, that NLO corrections may also affect the shapes of the kinematic distributions, as shown for the spin-1 and spin-0 cases in [17]; a detailed study of this aspect will be reported elsewhere.) It turns out that, for an on-shell mediator of given mass, the selection efficiencies are independent of the mass and spin of the invisible decay products. Moreover, contributions from off-shell production are negligible for the scenarios considered here. The efficiencies can thus be evaluated as a function of the mediator mass only; see also Appendix A.1. In the following, we normalise the number of events with NLO cross sections, shown in Fig. 3, and the total branching ratio into invisible final states (DM and neutrino). We note that for a given mediator mass the leading jet for the spin-2 mediator case is harder and more forward than that for the spin-1 case. This is partly because the spin-2 mediator with a parton is produced not only through the qq and qg initial states but also dominantly through the gg initial state, and partly because the spin-2 mediator is also emitted from a gluon as well as from the gggY 2 and qqgY 2 four-point vertices.   Figure 5 shows the ratio of signal events over the number of events excluded at 95% confidence level (CL), S/S 95 , as a function of the mediator mass, for the three types of DM (taking g X = 1 or 2 with Λ = 3 TeV, g SM = 1 and m X = 10 GeV as a benchmark case). As expected from the discussion in the previous section, the scalar DM case is the least constrained, with the / E T coming dominantly from the neutrino channel; for g X = 1 (2), we find the limit m Y 600 (750) GeV from the monojet analysis and m Y 750 (850) GeV from the multijet+ / E T analysis. 6 For Dirac DM the limit increases to m Y 950 (1300) GeV owing to the contribution from Y 2 → X D X D . Finally, for vector DM we have m Y 1100 (1550) GeV. For the monojet analysis, the inclusive SR with the / E T cut of 500, 600, and 700 GeV (denoted as IM5, IM6, and 6 While both analyses have very similar sensitivity, i.e. their expected limits are basically the same, the monojet results have over-and under-fluctuations in some SRs. Therefore the expected and observed limits slightly differ from each other for the monojet analysis. Overall, the multijet+ / E T analysis tends to give the stronger limit. IM7 in [5]) gives the limit for the low (100 ∼ 300 GeV), middle (300 ∼ 450 GeV), and high ( 450 GeV) mass region, respectively. For the multijet+ / E T analysis, the 2jet loose (2jl) SR gives the limit for the mass range of 100 ∼ 300 GeV, while the 2-jet medium (2jm) SR does for 300 GeV. See [43] for the detailed selection criteria.
As the production rate scales as 1/Λ 2 , the upper limit of Λ can be estimated from the plots. For instance, for vector DM with m Y = 100 GeV, Λ should be larger than around 10 TeV for g SM = g X = 1. It should be noted that, due to the K factors of 1.7−1.2 for m Y = 100−2000 GeV (see Fig. 3), these limits are slightly stronger than what would be obtained with LO production rates.
The 95% CL exclusion in the m X vs. m Y plane is shown in Fig. 6. Due to the different threshold behaviours, as seen in Eqs. (10)-(12), the excluded region near m Y = 2m X strongly depends on the type of DM.
We note that we compared the CheckMATE results with those obtained by the equivalent analysis implementations in MadAnalysis 5 [47,48] (recast codes [49,50]) and Rivet 2.5 [51] for a couple of representative mass  choices and found agreement at the level of 20% within all three tools. The monophoton (as well as mono-Z/W ) signature could also be interesting to explore the spin-2 model. However, as seen in Sec. 3.2, the production rate for a pair of DM with a photon is strongly suppressed. We checked that there is no constraint for the above benchmark points from the CMS 13 TeV monophoton analysis (12.9 fb −1 ) [12].
An interesting alternative to the universal coupling g SM is a leptophobic scenario with In this case, the / E T signatures come exclusively from decays into DM, because Y 2 decays into neutrinos are switched off. Moreover, constraints from dilepton resonance searches, which as we will see in the next subsection are quite severe, are evaded. The results for the leptophobic scenario are presented in Figs. 7 and 8 in analogy to Figs. 5 and 6. As expected, the m Y < 2m X region is no longer constrained. Also, for g X = 1, the exclusion becomes considerably weaker for all the DM types; in particular there is no more constraint for scalar DM. For g X = 2, except scalar DM, the mediator decays into DM dominates the neutrino decay mode even for the universal coupling scenario (see Fig. 2), and hence the m Y limits are very similar.

Resonance searches
Direct resonance searches can also be used to explore schannel mediator DM models, see e.g. [65,66] for the spin-1 and spin-0 mediator models, respectively. Results from Run-II data are already available for a large variety of final states (dijet, dilepton, diphoton, W W , ZZ, bb, tt, hh) from ATLAS [52-55, 57, 59, 60] and CMS [9, 56,58,[67][68][69], and give powerful constraints for mediator masses of a few hundred GeV up to several TeV. Lower masses are partly covered by Run-I results. 7 Table 2 lists the current resonance search results which we use to constrain our spin-2 simplified model. The RS massive graviton is considered in the analyses for pairs of electroweak gauge or Higgs bosons [54,57,58,60,63,64] Table 2. Constraints from resonance searches used in this study. The observed 95% CL upper limits on resonant production cross section (σ) times branching ratio (B) (times acceptance (A)) from each analysis are shown in Fig. 12 Constraints on Λ/g SM from resonance searches . Constraints on Λ/gSM from observed 95% CL upper limits of resonance searches at the 13 TeV (solid) and 8 TeV (dashed) LHC as a function of the spin-2 mediator mass. We assume a negligible branching ratio to DM, except for a dotted line, where the vector DM coupling gX /gSM = 10 with mX = 10 GeV is taken into account as a reference. Regions below each line are excluded. Information on the mediator width-to-mass ratio is given by the grey dotted lines.
as one of the new physics hypotheses. For the fermionic and jet final states in [52,53,55,56,59], on the other hand, Z and a model-independent Gaussian-shaped resonance have been studied. Except the dijet and di-b-jet analyses at 13 TeV and the low-mass diphoton analysis at 8 TeV from ATLAS, the limits are provided directly on the cross section in the given channel, and hence we obtain the model constraints by simply using the Y 2 production cross section and the branching ratio discussed in Section 3. For the analyses with different hypotheses from the spin-2 resonance, we assume that the acceptance and efficiency are similar. When limits are given on the fiducial cross section, σ × B × A, we generate LO events normalised by the NLO cross section and apply the fiducial cuts at the parton level by using MadAnalysis5 [70].
We recall that, for a given mediator mass, the Y 2 production cross section depends solely on g SM /Λ, while the branching ratio depends also on the parameters related to DM, i.e. g X and m X , as well as on the type of DM. In the decoupling limit of the dark sector, the constraints on Λ/g SM are the most stringent. When decays to DM are relevant, the branching ratios to SM particles become smaller and hence the constraints are weakened. Figure 9 shows the constraints on Λ/g SM from the observed 95% CL upper limits of the resonance searches listed in Table 2 as a function of the mediator mass, where we assume a negligible branching ratio to DM particles, i.e. g X 1 and/or m Y < 2m X . Although the branching ratio is small, B(Y 2 → γγ) ∼ 4% at high mass, the diphoton resonance searches give the most stringent limit for the whole mass range, resulting in Λ/g SM 100 TeV for m Y 1 TeV. The dilepton channel, also having a branching ratio of about 4%, provides a similarly strong constraint for mediator masses above 200 GeV. The dijet

Constraints on spin-2 mediated dark matter
Vector DM Fig. 10. Summary of the constraints on Λ/gSM from searches with and without missing energy at the 13 TeV LHC as a function of the spin-2 mediator mass, for mX = 10 GeV with gX /gSM = 1 (left) and 10 (right). The labelling of the constraints from resonance searches is the same as in Fig. 9. For gX /gSM = 1 the differences among the different types of DM for the limits from the resonance searches are not visible. For gX /gSM = 10, however, they are quite relevant so only the vector DM case is shown. The figure assumes a universal gSM but is also valid for the leptophobic case when ignoring the lines.
and W W/ZZ resonance searches lead to a constraint of a few tens of TeV on Λ/g SM for around 1 TeV mediator mass. We note again that the limits are obtained based on the NLO production rates which are larger than the LO ones, especially for pp → (Y 2 → jj)γ; see Fig. 3. We also note that, as indicated by grey dotted lines in Fig. 9, the mediator width can be very large at high mass and low Λ/g SM ; as the experimental analyses often assume a narrow width, this region has to be regarded with caution. The weakening of the constraints when Y 2 decays into DM are allowed is demonstrated for the dilepton channel in Fig. 9, depicted by a dotted line, where we assume vector DM and take g X = 10 and m X = 10 GeV. For instance, at m Y = 1 TeV, the dilepton (electron and muon) branching ratio becomes 0.8%, i.e. the dilepton production rate becomes smaller by a factor of five, reducing the limit on Λ/g SM by 1/ √ 5. As seen in Fig. 2, the above assumption gives the largest DM branching ratio within the scenarios we consider. 8 Therefore, the diphoton resonance searches, and for m Y > 200 GeV also the dilepton resonance searches, provide stronger constraints on the universal coupling scenario than the searches with missing energy.
To avoid such severe constraints from resonance searches, it is interesting to consider scenarios beyond the universal coupling case. The dilepton constraints could be avoided, for example, in the leptophobic scenario, g T = 0, as already discussed in the previous subsection. To avoid the diphoton constraints is somewhat more complicated. One possibility would be the gravity-mediated DM model [13,14], where the KK graviton mainly couples to massive particles -DM, Higgs, massive gauge bosons and top quarks-while the couplings to photons, gluons and light quarks are highly suppressed. In such scenarios, the branching ratios and the production cross sections of the spin-2 resonance strongly depend on the setup and can be very different from those in the universal coupling case. In fact associated production of the mediator with a W or Z boson, or mediator production in vector boson fusion may be more relevant than s-channel production in qq or gg fusion. While such setups can in principle be studied easily in the simplified model framework by appropriately choosing the free parameters g T X and g T i in Eq. (3), such an analysis is beyond the scope of this paper. A final caveat is that nonuniversal couplings to gluons and quarks, g T g = g T q , give rise to a unitarity violating behaviour at higher order in QCD [36]. We therefore only consider phenomenological scenarios with g T g = g T q .

Summary
We considered a simplified DM model where the DM candidate couples to the SM particles via an s-channel spin-2 mediator, Y 2 , and studied the constraints from the current LHC data. In particular, we compared the constraints from searches with and without missing energy. For universal couplings of the mediator to SM particles, we found that diphoton resonance searches provide the strongest constraints, Λ/g SM 100 TeV for Y 2 masses up to ∼ 1 TeV. For Λ/g SM = 10 (3) TeV, the exclusion extends up to 4 (beyond 5) TeV in m Y . The dilepton channel provides a similarly strong constraint for mediator masses above 200 GeV. Monojet and multijet+ / E T searches are competitive only if the mediator decays into photons and leptons are heavily suppressed; in this case they could provide complementary constraints to the other resonance searches in particular in the low-mass region below 0.5-1 TeV, depending on g X /g SM .
For m Y < 2m X , / E T signatures arise solely from Y 2 decays into neutrinos, leading to m Y 700 GeV for g X /Λ = g SM /Λ = (3 TeV) −1 , based on 3.2 fb −1 of data at √ s = 13 TeV. For m Y > 2m X , the limit crucially depends on g X and the type of dark matter. The dependence on the DM mass is less pronounced unless one approches the threshold region. For m X = 10 GeV and g X /Λ = g SM /Λ = (3 TeV) −1 , we found m Y 750, 950, and 1100 GeV for scalar, Dirac, and vector DM, respectively. This increases to 850, 1300, and 1550 GeV when doubling g X . We note that the obtained limits are based on the NLO-QCD predictions, which give a larger production rate than at the LO. The K factor depends on the mediator mass and the production channel.
The complementarity among the different searches is illustrated in Fig. 10, where we have rescaled the reach of the jets + / E T searches from 3.2 to 15 fb −1 in order to make a fair comparison. We see that, for the same amount of data, in case of g X g SM the missing energy searches are roughly competitive with the dijet and heavy diboson (W W , ZZ) searches, pushing Λ/g SM beyond 20 TeV. (As mentioned, when the dilepton and diphoton constraints hold, they give even stronger limits.) For g X /g SM = 10 (or g X /ĝ SM = 10), also the resonance constraints strongly depend on the type of DM. Therefore, in the right plot in Fig. 10 only the vector DM case is shown. We see that the jets+ / E T searches give stronger constraints than the dijet and heavy diboson searches up to mediator masses of about 1.2 TeV. The dilepton and diphoton constraints are weakened by about a factor of two but still give the strongest constraints.
We hope our work will be useful to find reasonable benchmark scenarios for spin-2 mediated DM searches at the LHC as well as to construct viable UV-completed models which can give predictions for those parameters. We also note that our study on resonance searches in Sec. 4.2 can be applied not only for spin-2 mediated DM models but also for usual RS-type graviton searches; see also, e.g. [71]. As a final remark we like to point out that in a full model the presence of KK excitations might alter the LHC phenomenology as compared to the simplified model scenarios discussed here. Examples are limits on gauge KK modes providing additional constraints on light gravitons, or KK excitations of the DM fields contributing to / E T signatures. While this goes well beyond the simplified model picture, it is certainly an interesting topic for future studies. A Supplemental material for recasting A.1 Searches with missing energy As mentioned in the main part of the paper, in case of the monojet and the 2-6 jets + / E T searches, the signal comes solely from on-shell mediator production with the Y 2 decaying into neutrinos and/or DM. The signal selection efficiency (more precisely acceptance times efficiency, A × ) depends only on the properties of the mediator, but not on those of the invisible decay products. Figure 11 shows A × for those SRs which, depending on m Y , can be the most sensitive ones in each of the two ATLAS analyses considered in this paper. As a service to the reader and potential user of our work, the complete A × tables for all SRs are available in numerical form at [72].

A.2 Resonance searches
In Fig. 12 we show observed 95% CL upper limits on resonant production cross section times branching ratio (times Observed 95% CL upper limits on resonant production cross section times branching ratio (times acceptance) as a function of the resonance mass from each experimental paper; see Table 2 for more detailed information. Dashed lines denote limits including cut acceptance. For reference, NLO production cross sections of the spin-2 mediator are shown by dotted lines for different values of gSM/Λ. acceptance) as a function of the resonance mass from each experimental paper. The analyses denoted by solid lines present the limit on σ × B, while those by dashed lines provide the limit on σ × B × A (× for bb); see Table 2 for more detailed information.
As indicated in Table 2, the dijet (+ ISR jet/photon) and tt analyses at 13 TeV as well as the ATLAS 8 TeV diphoton analysis provide tables with the numbers corresponding to the lines in the exclusion plots, which is very convenient for our purpose. The other analyses do not provide explicit values, and hence we have to extract these data from the exclusion plots 'by hand', e.g. using WebPlotDigitizer [73], a public software. To avoid that other people have to redo this exercise, our digitised data files are available at [72] and on the new PhenoData database [74]. We encourage the experimental collaborations to provide digitised data together with their plots, in order to make it easier to use their results.
Finally, we notice a caveat regarding the re-interpretation of the low-mass resonance search in dijet plus ISR final states [53]. We found that final state radiation (FSR) may be also important and give rise to a non-trivial structure in the dijet invariant mass spectrum. Technically, simulated event shapes can differ by including FSR or not in the matrix elements, which may affect the parameter fitting procedure for a bump search.