A universal framework for t-channel dark matter models

We present the DMSimpt model implementation in FeynRules, which aims to offer a unique general framework allowing for all simulations relevant for simplified $t$-channel dark matter models at colliders and for the complementary cosmology calculations. We describe how to match next-to-leading-order QCD fixed-order calculations with parton showers to derive robust bounds and predictions in the context of LHC dark matter searches, and moreover validate two model restrictions (relevant for Dirac and Majorana fermionic dark matter respectively) to exemplify how to evaluate dark matter observables to constrain the model parameter space. More importantly, we emphasise how to achieve these results by using a combination of publicly available automated tools, and discuss how dark matter predictions are sensitive to the model file and software setup. All files, together with illustrative Mathematica notebooks, are available from the URL http://feynrules.irmp.ucl.ac.be/wiki/DMsimpt.


Introduction
Despite of convincing evidence for its existence [1], dark matter still evades direct detection both in dedicated underground nuclear recoil experiments and at colliders. Getting insights on the nature of dark matter and the way in which it interacts with the Standard Model particles therefore consists in one of the hot topics in particle and astroparticle physics today. One potential strategy that could shed light on this matter involves simplified models [2,3] in which the Standard Model is minimally extended in terms of particles and new couplings. This approach allows for the exploration a e-mail: chiara.arina@uclouvain.be b e-mail: fuks@lpthe.jussieu.fr (corresponding author) c e-mail: luca.mantani@uclouvain.be of viable dark matter scenarios in a model-independent way and the comparison of theoretical predictions with results of direct, indirect and collider searches. This however requires the ability of making predictions for large classes of models, both at colliders and for what concerns cosmology.
The FeynRules package [4] offers such a possibility, as from a unique FeynRules implementation of any given dark matter model, it is subsequently possible to generate model files suitable for various high-energy physics tools such as MG5_aMC [5], MadDM [6] or MicrOMEGAs [7]. Following the general strategy for new physics computations outlined in Ref. [8], such a joint usage of various packages has two major advantages in the dark matter context. First, it allows for the straightforward and automatic calculation of the dark matter relic density, as well as of the direct and indirect detection cross sections to verify the cosmological viability of any model. Second, it enables the extraction of the exclusion levels of various searches at colliders through the automated generation of realistic collision events and the recasting of the corresponding LHC analyses. In the latter case, the MG5_aMC framework allows in particular for simulations including next-to-leading order corrections in α s , so that predictions are accurate enough to derive robust constraints when LHC recasting is at stake through, e.g., the MadAnalysis 5 platform [9] that includes, from version 1.8, the propagation of the theoretical uncertainties on the signal predictions up to the derived exclusion levels [10].
In most simplified models for dark matter, the dark matter is assumed to be a single massive particle that interacts weakly with the Standard Model through a mediator particle. In s-channel setups [11][12][13][14], the mediator is colour-singlet and is enforced to couple to a pair of either dark matter particles, or Standard Model particles. Such a configuration generally arises in scenarios in which the dark matter stability is guaranteed by a Z 2 discrete symmetry under which all Standard Model fields and the mediator are even, and the dark matter particle is odd. A comprehensive approach for achieving automatic and straightforward cosmological calculations and collider simulations for s-channel dark matter models has been recently proposed [13], the cornerstone being a unique FeynRules implementation driving any subsequent computation.
The present work is dedicated to a general implementation, in the FeynRules package, of a large set of t-channel dark matter models in which the mediator interacts with one of the Standard Model quarks and dark matter. We have used this FeynRules implementation to generate a UFO library [15]. that can subsequently be imported in programmes like MG5_aMC or MadDM for undertaking various simulations and computations for a large class of tchannel dark matter models. Our implementation in particular allows for collider simulations systematically including next-to-leading-order (NLO) QCD corrections to all new physics processes involving either the dark matter particle, the mediator or both.
Such a possibility requires however a specific treatment of the real emission contributions that feature, in t-channel dark matter models, narrow s-channel resonances. Realemission corrections to a given process (e.g. dark matter pairproduction whose real emission contributions include the production of a system comprised of a dark-matter pair and a jet) may indeed include partonic sub-processes featuring an s-channel resonance corresponding to another Born process (e.g. mediator/dark matter associated production) followed by a tree-level decay (e.g. mediator decay into dark matter and jets). The integration of such contributions over the phase space leads to a growth proportional to an inverse power of the resonance width, so that such contributions could be numerically dominant and apparently spoil the convergence of the perturbative series. We recall that this type of configuration also exists in the Standard Model, in particular in the context of t W production. Real correction to the latter process include diagrams describing the production of an s-channel resonant tt final state, followed by a top decay into a W b system.
Moreover, when all new physics processes allowed by the model are considered as a whole (as each subprocess contributes to the new physics signal), these resonant diagrams could be double-counted and lead to incorrect predictions. They therefore need to be treated consistently. Different strategies to treat these resonances have been recently automated within the MG5_aMC framework [16], hence enabling NLO QCD simulations for the considered t-channel dark matter models in a way that is as easy as for the s-channel case.
In order to illustrate the strength of our approach, we focus on two limiting cases and study their phenomenology at colliders and in cosmology, which allows for the validation of our implementation. We in particular compare the performances of MadDM and MicrOMEGAs and present, for the first time, automated computations for loop-induced processes relevant for dark matter indirect detection. Such a feature, which will be available from the next release of MadDM, greatly eases the phenomenological analysis of tchannel dark matter models. More specifically, we consider the case of a fermionic dark matter particle whose interactions with the Standard Model are mediated by a scalar particle coupling to the right-handed up quark, both for what concerns Dirac and Majorana dark matter. In the following, we coin these two configurations, that have been vastly studied in the literature (as shown e.g. in Refs. [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32]) and that are particularly promising for LHC and dark matter searches (see e.g. Refs. [33,34]), the S3D_uR and S3M_uR models, respectively.
The rest of the paper is structured as follows. In the next section, we present the model conventions, its implementation into FeynRules and the restrictions (i.e. the limiting cases) shipped with our general implementation. In Sect. 3, we detail how to match NLO QCD calculations with parton showers for collider simulations, providing extensive details on how to make use of MG5_aMC in order to ensure a consistent treatment of the resonant contributions appearing at O(α s ). We then present, for the first time, total rate and differential distributions extracted from accurate predictions matching NLO QCD calculations with parton showers, and derive the corresponding constraints from selected LHC searches. In Sect. 4, we briefly outline the dark matter observables relevant for t-channel dark matter models, how to compute them with MadDM, and present the results for the S3D_uR and S3M_uR model restrictions to validate our implementation against known results. We summarise our work in Sect. 5.

Generalities
We consider a generic t-channel dark matter simplified model in which the Standard Model (SM) is extended by several incarnations of two extra fields, a dark matter candidate (that we generically denote by X ) and a mediator lying in the fundamental representation of SU (3) c (that we generically denote by Y ). In order to maintain the model as general as possible, we allow for several options for the spin of the new particles and therefore include six new dark matter fieldsS, S,χ, χ ,Ṽ μ and V μ , all lying in the singlet representation of the SM gauge group SU (3) c × SU (2) L × U (1) Y . These fields respectively correspond to a real scalar field, a complex scalar field, a Majorana spinor, a Dirac spinor, a real vector field and a complex vector field.
2,000,001 2,000,003 2,000,005 The most general Lagrangian embedding all the interactions of these fields with the SM can be written, after imposing that electroweak gauge invariance is preserved, as where L SM is the SM Lagrangian and L kin contains gaugeinvariant kinetic and mass terms for all new fields. The fermionic, scalar and vector dark matter Lagrangians read In our notation, Q stands for the SU (2) L doublet of lefthanded quarks and u and d are the up-type and down-type SU (2) L singlets of right-handed quarks respectively. The scalar mediators ϕ Q , ϕ u and ϕ d are chosen to solely interact with the Q, u and d quarks, as for the fermionic mediators ψ Q , ψ u and ψ d (that are thus vector-like). The mediators therefore lie in the same SM representation as their quark partners. In the above expression, we have understood all flavour indices for clarity. The λ Q , λ u and λ d coupling strengths are hence 3×3 matrices in the flavour space, that we moreover consider real and flavour-diagonal for simplicity.
The new physics particles of the simplified model are given in Table 1, together with their representation under the gauge and Poincaré groups, their potential Majorana nature, the adopted particle name in the FeynRules implementation and the adopted Particle Data Group (PDG) identifiers [35]. The conventions for the different coupling parameters are summarised in Table 2, in which they are given together with the name used in the FeynRules implementation and the Les Houches (LH) blocks [36] storing their numerical values when running tools like MG5_aMC or MadDM. By relying on a joint usage of the FeynRules [4], NLOCT [37] and FeynArts [38] packages, we automatically generate a UFO model [15] that can be used by MG5_aMC for both leading order (LO) and NLO computations. This UFO model includes all UV counterterms allowing for the renormalisation of the model with respect to the QCD interactions, as well as all R 2 Feynman rules that are relevant for the numerical evaluation of one-loop integrals in four dimensions. Table 3 List of all restrictions included in the DMSimpt UFO model. In each case, the simplified model contains a single class of massdegenerate mediators (where f = 1, 2, 3 is a flavour index), a specific dark matter candidate and universal and flavour-conserving dark matter couplings λ ϕ andλ ψ

Name DM Mediators Parameters
S3M_uniχ The model is shipped with a large ensemble of restrictions dedicated to specific t-channel simplified models. These are summarised in Table 3 where for each restriction, we specify the active new physics states, all other states being taken decoupled and non-interacting. In other words, each restriction consists in a simplified model in which the SM is extended by a specific class of mediators, and a given dark matter state. In order to reduce the number of free parameters, all (active) mediators are taken mass-degenerate. A given restriction named XYZ can be loaded in MG5_aMC In the model restrictions whose name ends with the uni suffix, all twelve flavours of mediators are considered, their mass and interaction strengths being taken flavourconserving and universal, for F = Q, u and d. In model restrictions of the uR class, only mediators coupling to the right-handed up quark are taken as active, all other couplings being vanishing, whilst in the 3rd class of model restrictions, we only consider the mediator coupling to the third generation of SM quarks, all other couplings being again assumed vanishing.

The S3M/S3D class of models
In S3M-type and S3D-type models, the dark matter is taken to be respectively the Majorana and Dirac stateχ and χ of mass M χ . As written in Sect. 2.1, all mediators are considered degenerate of mass M ϕ , and all new physics interactions are universal and flavour-conserving with a global strength λ ϕ . The generic Lagrangian L F of Eq. (2) therefore simplifies to where X = χ (S3D) orχ (S3M) equivalently refers to Dirac or Majorana dark matter, and f is a generation index. The model is thus defined by three parameters, In the universal S3M_uni and S3D_uni restrictions, the simplified model includes all twelve mediators, whilst in the S3M_3rd and S3D_3rd restrictions, the setup is further simplified and dark matter only couples to the third generation via the four corresponding mediators. In the S3M_uR and S3D_uR restrictions, only a coupling to the right-handed up quark u 1 is considered, through a single mediator. The associated Lagrangians read,

The F3S/F3C class of models
In F3S-type and F3C-type models, the dark matter consists of the real and complex scalar stateS and S of mass M S respectively. As in the previous subsection, all mediators are assumed to be degenerate of mass M ψ , and all new physics interactions are universal and flavour-conserving with a strengthλ ψ . The Lagrangian L S of Eq. (2) therefore simplifies to where X =S (F3S) and S (F3C) in the real and complex case. The model is defined by three parameters, In the universal F3S_uni and F3C_uni restrictions, dark matter couples to all SM quark eigenstates through twelve mediators. In the third generation F3S_3rd and F3C_3rd models, its couplings are restricted to the bottom and top quark ones and the corresponding four mediators, while in the F3S_uR and F3C_uR models, dark matter only couples to the right-handed up quark. The associated Lagrangians are 2.4 The F3V/F3W class of models In the F3V and F3W types of models, the dark matter is a real and complex vector stateṼ μ and V μ of mass M V respectively. All mediators are degenerate of mass M ψ , and all new physics interactions are universal and flavour-conserving with a common strengthλ ψ . The Lagrangian L V of Eq. (2) is simplified to where X =Ṽ (F3V) or V (F3W) in the real and complex vector case. The model is defined by three parameters, In the F3V_uni and F3W_uni restrictions, all twelve mediators are included. In contrast, in the F3V_3rd and F3W_3rd restrictions, only the four mediators relating dark matter to the top and bottom quarks are included, whilst in the F3V_uR and F3W_uR models, the only non-vanishing coupling is the one to the right-handed up-quark. The associated Lagrangians read 3 Matching NLO QCD fixed-order calculations with parton showers

Generalities
In the class of simplified models under consideration, the computation of NLO QCD corrections involve real emission diagrams possibly featuring intermediate s-channel resonances. These should be treated consistently in order not to apparently spoil the convergence of the perturbative series by yielding NLO cross sections much larger than the associated LO ones. This occurs when the cross section related to the production of the resonant state is much larger than the one of the initially considered process. Moreover, we aim at combining events describing all possible new physics processes of a given model at the NLO accuracy in QCD. We will hence consider the production of a pair of dark matter particles ( pp → X X), of any mediators ( pp → Y i Y j ), as well as the associated production of a mediator and a dark matter state ( pp → XY i ). Therefore, the subtraction of all resonant contributions in the real corrections is mandatory to avoid their double-counting when combining the three types of processes. Different strategies dealing with the treatment of these resonances have been recently automated within the MG5_aMC framework [16]. They include diagram removal methods with or without keeping the interferences between the resonant and non-resonant contributions [39], as well as various techniques to subtract the resonant contribution from the full amplitude [40]. In the following, we employ one of such strategies, in which all squared resonant diagram contributions are discarded whilst the interferences of the resonant and non-resonant diagrams are kept. All available methods should however lead to numerically similar results if the resonant process can be consistently defined.
In practice, MG5_aMC has to be run together with the MadSTR plugin 1 that can be activated by starting MG5_aMC as

mg5_aMC --mode=MadSTR
The code is then used to simulate events, at the NLO accuracy in QCD, relevant for all new physics processes allowed by t-channel dark matter models. The considered processes can be classified into three categories, This corresponds to the production of a pair of dark matter particles (generically denoted by X X), the associated pro- Representative LO Feynman diagrams describing the production of a pair of dark matter particles (left) and the associated production of a mediator with a dark matter state (right). The mediator decay into dark matter and a quark is included duction of a mediator and a dark matter particle (generically denoted by XY ), and the production of a pair of mediators (generically denoted by Y Y ). In the latter two cases, the mediator further decays into a SM quark and dark matter. After simulating each process separately, the different contributions are combined, which is only possible if all resonant pieces from the real emission to the three subprocesses are subtracted. For instance, the diagrams associated with the second Born subprocess ( pp → XY → X X j) are included in those related to the real corrections to the first subprocess ( pp → X X). In order to avoid any double counting, we include the resonant component into the Born contribution to XY production, and the non-resonant one into the real corrections to X X production.
We import the DMSimpt UFO model in MG5_aMC to deal with the generation of hard-scattering events at the NLO accuracy for all the processes of Eq. (15), using the MadSTR plugin and convoluting the matrix elements with the NLO set of NNPDF 3.0 parton distribution functions (PDFs) [41] accessed via the LHAPDF 6 library [42]. Mediator decays are handled with the MadSpin [43] and MadWidth [44] programmes, which allows for the factorisation of the production and decay processes in a way retaining both off-shell and spin correlation effects. The resulting partonic events are matched with parton showers as described by the Pythia 8 package [45], following the MC@NLO prescription [46]. We also use Pythia 8 to handle hadronisation. We then reconstruct the hadron-level events by clustering hadrons according to the anti-k T algorithm with a separation parameter set to ΔR = 0.4 [47], as implemented in the FastJet software [48] that we drive from MadAnalysis 5 [9,49]. The latter programme is also used for the generation of the differential distributions studied in Sect. 3.3, and the reinterpretation analysis of the LHC results in Sect. 3.4.

Simulating an S3D_uR dark matter signal
In the following, we illustrate how NLO predictions matched with parton showers can be achieved in the S3D_uR class of model. In the X/Y notations of Eq. (15), we thus have, Fig. 2 Same as Fig. 1 but for mediator pair-production (and decay) in the QCD (left) and t-channel dark matter exchange (right) channels Events originating from dark matter pair production at NLO ( pp → X X; see e.g. the left panel of Fig. 1 for a representative LO Feynman diagram) are generated by starting MG5_aMC with the MadSTR plugin switched on. The usual generate and output commands available from the MG5_aMC command line interface [5] are then cast, after having imported the restricted UFO model, In order for the restricted model to be dealt with consistently, it is crucial to explicitly forbid any decoupled particle to run into any virtual loop. This is implemented at the level of the generate command, in which we manually exclude all fermionic mediators, all scalar mediators not coupling to the right-handed up quark u R and all irrelevant dark matter states of the model. The UFO conventions for the particle names follow the FeynRules ones introduced in Table 1 (we recall that the MG5_aMC command line interface is case insensitive), additionally including an integer number for the generation indices. On run time, the MadSTR plugin takes care of identifying and treating any potentially resonant contribution. In our case, the squared resonant contributions are discarded, whilst the interferences of the resonance with the non-resonant continuum are kept.
Events describing the associated production of a mediator with a dark matter particle ( pp → XY ; see e.g. the right panel of Fig. 1 for a representative diagram including the mediator decay process) are generated in a similar fashion, import model DMSimpt-S3D_uR --modelname define dm = xd xdd efine yy1 = ys3u1 ys3u1g enerate p p > dm yy1 / yf3qu1 yf3qu2 \ yf3qu3 yf3qd1 yf3qd2 yf3qd3 yf3u1 yf3u2 \ yf3u3 yf3d1 yf3d2 yf3d3 ys3qu1 ys3qu2 \ ys3qu3 ys3qd1 ys3qd2 ys3qd3 ys3u2 ys3u3 \ ys3d1 ys3d2 ys3d3 xs xm xv [QCD] output where we make use of the dm and yy1 multiparticle labels to guarantee that all potential particle/antiparticle combinations are accounted for.
Mediator pair production is generally dominated by QCD contributions that are independent of the dark matter mass and couplings, as illustrated by the first Feynman diagram of Fig. 2 that also includes the mediator decay process. However, if λ ϕ is large enough, t-channel dark matter exchanges, as depicted by the second diagram of Fig. 2, could significantly contribute. NLO QCD corrections to the strong contributions to mediator pair-production of O(α s ) can be straightforwardly calculated by typing in the MG5_aMC command line interface, import model DMSimpt-S3D_uR --modelname define yy1 = ys3u1 ys3u1g enerate p p > yy1 yy1 / yf3qu1 yf3qu2 \ yf3qu3 yf3qd1 yf3qd2 yf3qd3 yf3u1 yf3u2 \ yf3u3 yf3d1 yf3d2 yf3d3 ys3qu1 ys3qu2 \ ys3qu3 ys3qd1 ys3qd2 ys3qd3 ys3u2 ys3u3 \ ys3d1 ys3d2 ys3d3 xs xm xv [QCD] output without using the MadSTR plugin. Making use of the coupling order information of the model, MG5_aMC automatically restricts the process to its pure QCD contribution, neglecting any t-channel dark matter exchange at the Born level and any real emission or virtual contribution depending on λ ϕ . The considered Born contribution is thus of O(α 2 s ) whilst the NLO component is of O(α 3 s ) and free of any resonance.
The pure t-channel contribution can be evaluated, at the NLO accuracy in QCD, by typing.
Care must be taken when dealing with the mixed-order interferences of the QCD diagrams with the t-channel ones. The version 3.x.y of MG5_aMC being incompatible with MadSTR and the version 2.6.x of the code being unable to handle mixed orders, there is not any publicly available and user-friendly option. One possible way to cure this issue would be to include in the UFO model all UV counterterms and R 2 rational terms necessary for mixed-order NLO calculations in QCD, QED and in the new physics λ coupling, and to implement in MadFKS [50] all necessary subtraction terms. This however goes beyond the scope of this work.
We therefore adopt the strategy of simulating the interferences at LO, and reweight the events by a K -factor assumed to approximate the effect of the QCD corrections. We multiply the interference event weights by the geometric mean of the pure QCD and pure t-channel K -factors, those two NLO to LO ratios being defined differentially. In other words, each distribution will be reweighted bin by bin. Event simulation for the interferences is then performed by typing, in the MG5_aMC command line interface, import model DMSimpt-S3D_uR --modelname define yy1 = ys3u1 ys3u1g enerate p p > yy1 yy1 DMTˆ2==2 / yf3qu1 \ yf3qu2 yf3qu3 yf3qd1 yf3qd2 yf3qd3 yf3u1 \ yf3u2 yf3u3 yf3d1 yf3d2 yf3d3 ys3qu1 \ ys3qu2 ys3qu3 ys3qd1 ys3qd2 ys3qd3 ys3u2 \ ys3u3 ys3d1 ys3d2 ys3d3 xs xm xv output 3.3 Total and differential cross sections for S3D_uR dark matter In order to illustrate how all subprocesses of Eq. (15) could impact a dark matter signal at the LHC, we consider two benchmark scenarios representative of the S3D_uR model. We fix the dark matter and mediator masses to and the new physics coupling to In the first scenario S1, the spectrum is more compressed although there is enough phase space for the light mediator to decay into a dark matter particle and a hard jet. In the second scenario S2, the mediator is heavier, its mass being fixed to a more realistic value with respect to current squark mass limits [51,52]. Whilst present supersymmetry bounds on the strongly-interacting superpartners are usually stricter, they are not directly applicable to our setup by virtue of the different nature of the dark matter and mediator particles. We therefore ignore them for now and address this point in Sect. 3.4.
In Table 4, we present total cross sections for the various processes of Eq. (15) and for the two considered benchmarks, both at LO (upper panel) and NLO (lower panel) accuracy, Table 4 Total cross sections at LO (upper) and NLO (lower), in fb, for the subprocesses of Eq. (15) and the benchmark scenarios defined in eqs. (16) and (17). Our predictions are given together with the scale and parton density uncertainties and for a setup in which the pp → X X process simulation includes a transverse momentum ( p T ) cut of 100 GeV on the leading jet at the matrix-element-generator level. For each of the subprocesses, the NLO K -factor defined as the ratio of the NLO predictions to the LO one is large. This emphasises the relevance of using rates that are NLO-accurate to avoid underestimating signal yields. Our LO and NLO predictions also include theoretical scale uncertainties originating from missing higher-order corrections and those associated with the parton density fit. We estimate the former by a nine-point independent scale variation in which the renormalisation and factorisation scales are varied by a factor of 2 up and down with respect to a central scale set to the average transverse mass of the final-state particles. Except for dark-matter pair-production ( pp → X X) that is insensitive to α s at the lowest order (see e.g. the left panel of Fig. 1), LO predictions are affected by large scale uncertainties that are significantly reduced when NLO corrections are included. This consists in the second major benefit of higher-order calculations: the reduction of the theoretical systematics. The second source of theoretical uncertainties, the PDF errors, yields a similar effect at LO and NLO as the same parton density set has been used. Those errors are reasonably small as our benchmark scenarios feature masses leading to a moderate Bjorken-x regime. Additionally, we have verified that the QCD contribution to mediator pair production agrees with the expectation for squark pair-production in a supersymmetric simplified scenario in which all superpartners but a single squark are decoupled [16].
In the right panel of the table, we investigate the impact of the QCD and t-channel contributions to the production of a pair of mediators. The adopted scenarios, with their large coupling choice of Eq. (17), ensure that the t-channel contribution is relevant and cannot be neglected. On the contrary, for slightly smaller coupling choices, only QCD production would remain, as the t-channel amplitude squared is proportional to λ 4 ϕ and the interference between the QCD and t-channel mode to λ 2 ϕ . For λ ϕ = 1, the relative importance of the QCD and t-channel modes turns out to differ for the two benchmark points under consideration. In the case of the S1 setup, QCD contributions dominate, as expected for such a small mediator mass M ϕ = 500 GeV. In contrast, for scenarios like the S2 scenario in which the mediator is much heavier, the QCD production mode is suppressed by virtue of the steeply falling production rate with M ϕ , so that the t-channel contribution dominates. For both cases, the two contributions are however of a similar order of magnitude and their destructive interferences are large.
The two channels are sensitive to different initial partonic luminosities. QCD production is mostly induced by gluon fusion (at 80% and 60% in the S1 and S2 cases respectively) and t-channel production by quark-antiquark scattering). They feature different PDF uncertainties, as well as a different dependence on the scales.
In Figs. 3 and 4, we present a few properties of the new physics signal induced by the benchmark scenarios S1 and S2 respectively. We focus on observables that are relevant for dark matter searches at the LHC in the monojet channel and consider the description of the missing transverse energy and jet activity. We show differential distributions for the transverse momentum of the leading jet p T ( j 1 ) (left panel), the missing transverse energy E miss T (central panel) and the total hadronic activity H T (right panel) defined as the scalar sum of the p T of all reconstructed jets. For each observable, we present predictions at LO (dashed lines) and NLO (solid lines) for the individual contributions of the processes shown in Eq. (15), as well as for their sum (red). We hence distinguish the QCD-induced (teal) and t-channel-mediated components (blue) of the mediator pair-production channel (Y Y ), the dark-matter pair production mode (X X, purple) and the associated production of a mediator with dark matter (XY , green). The shaded areas around the NLO results correspond to the uncertainty bands obtained as described above, i.e. from a nine-point variation of the unphysical factorisation and renormalisation scales.
Our results show that the dark matter pair-production channel, despite a large production cross section (see Table 4), mainly yields events featuring a small amount of missing energy and not so much hadronic activity, even if at the matrix-element level, our simulation includes a selection cut of 100 GeV on the p T of the leading (parton-level) jet. After matching the fixed-order NLO predictions with parton showers, the emissions originating from this hard parton are often not reclustered back so that a small ensemble  of softer jets are finally reconstructed from the initial hard radiation. Consequently, the dark matter particles turn out to be mostly produced back-to-back, which leads to a small amount of missing energy, and in association with a small hadronic activity. As a consequence, the efficiency of a typical monojet selection is expected to be quite reduced as one generally requires a substantial amount of missing energy and hadronic activity. We refer to Sect. 3.4 for more details.
The leading relevant process from the cross sections presented in Table 4 therefore consists in the associated production of a dark matter particle with a heavier mediator pp → XY . The Y → Xq decay of the mediator leads to the production of a second dark matter state together with a hard parton, which guarantees a much larger missing transverse energy and hadronic activity than for the pp → X X channel. The results, depicted by the green curves in Figs. 3 and 4 for the S1 and S2 scenarios respectively, confirm this, the corresponding E miss T , p T ( j 1 ) and H T distributions being depleted in the low-energy regime. In the high-energy tails, the XY contributions moreover almost match entirely the total new physics signal (red curves), the pp → Y Y contributions being only expected to take over in the very hard part of the phase space (not represented in Fig. 3 and that is drastically phase-space-suppressed, and thus not shown, for the S2 setup in Fig. 4). The XY distributions are indeed steeply falling with the energy scale, compared with the Y Y ones, so that the XY component of the signal is only dominant for moderate observable values of a few hundreds of GeV. The relative importance of the XY process can however be tamed down by reducing the magnitude of the λ ϕ coupling, on which the normalisation of the distributions depends quadratically as the amplitude of the corresponding partonic process qg → XY is linear in both the new physics and strong coupling constants (see the right diagram in Fig. 1). Finally, mediator pair production ( pp → Y Y ) only dominates, as said above, in the harder part of the spectra where all other contributions are kinematically suppressed. In that regime, the decays of the two heavy mediators into Xq systems guarantee an amount of missing energy and hadronic activity greater than for XY production, despite the global rates being reduced by the large mediator mass. In our analysis, we distinguish the QCD production mode whose matrix element is proportional to α 2 s and independent of λ ϕ , and the t-channel one that depends on λ 4 ϕ . Whereas for the adopted λ ϕ = 1 benchmark value, the two production channels contribute equivalently, the t-channel impact can be reduced by fixing λ ϕ to a smaller value. For λ ϕ ∼ g W (g W being the weak coupling constant), we obtain the widely studied supersymmetric limiting scenario in which all superpartners except the right-handed up squark and a bino-like neutralino are decoupled, with a difference on the dark matter nature that is here a Dirac fermion. In this case, only the QCD production of two mediators matters, and the t-channel contribution to Y Y production and XY production can be ignored (at least for the considered mediator masses).

Collider constraints on S3D_uR dark matter
As visible from the results of Sect. 3.3, most of the monojet signal arises, in the considered S1 and S2 benchmark scenarios, from the production of heavy mediators (by pairs or in association with dark matter) that then decay into dark matter and jets. In the following, we reinterpret the results of typical LHC dark matter searches probing final states fea-turing a large amount of missing transverse energy (carried away by the dark matter particles) and an important hadronic activity. We recast two of such ATLAS analyses for which reimplementations within the MadAnalysis 5 Public Analysis Database (PAD) of recasted LHC analyses [58] 2 exist. Starting from Monte Carlo simulations of the X X + XY +Y Y dark matter signal as described in Sect. 3.2, we make use of MadAnalysis 5 to automatically simulate the response of the ATLAS detector by means of appropriate tunes of the Delphes 3 programme [59]. We then assess the sensitivity of the considered analyses to the S1 and S2 signals by using the CL s method [60].
We consider two ATLAS analyses of 36.2 fb −1 of LHC data targeting the production of missing energy recoiling against at least one hard jet and a subleading hadronic activity. We recast the ATLAS-EXOT-2016-27 analysis [53 -55] in which the selection imposes that the dark matter system is produced together with 2 to 4 extra hard jets with quite stringent kinematic requirements. The analysis includes an ensemble of signal regions that are distinguished by different inclusive and exclusive constraints on the missing transverse energy. In the ATLAS-SUSY-2016-07 analysis [56,57], a larger number of jets N j is allowed (N j ≥ 2) and the properties of those jets are less constrained. The analysis includes several signal regions that mainly differ by the minimum number of required jets and a constraint on the effective mass M eff defined by   Table 5 but for the ATLAS-SUSY-2016-07 analysis [56,57]. We indicate here the jet multiplicity requirement and the effective mass M eff threshold defining the most sensitive signal region Our results are presented in Tables 5 and 6 for the ATLAS-EXOT-2016-27 and ATLAS-SUSY-2016-07 analysis, respectively. In each table, we show the confidence level (CL) exclusion obtained when the analysis signal regions are populated by all the X X, XY and Y Y contributions to the signal. The impact of the individual channels is also reported, the Y Y component being further decomposed into its QCD and tchannel part. Our results include theoretical scale uncertainties, which we have extracted by propagating the uncertainties on the total cross sections down to the CL s exclusions that we have computed both at LO and NLO. In our recasting procedure, we conservatively make use of the most sensitive signal region of each analysis, to derive the exclusion levels, as the statistical model used by the ATLAS collaboration for the combination of the various regions is not publicly available. The definition of these regions is provided in the tables, that hence include the required E miss T range for the ATLAS-EXOT-2016-27 analysis, and the thresholds on N j and M eff for the ATLAS-SUSY-2016-07 analysis.
It turns out that both the S1 and S2 scenarios are excluded at the 95% CL by both analyses, even after accounting for the uncertainties on the total rates. However, such a conclusion can only be drawn when more precise NLO simulations are employed and after summing over the X X, XY and Y Y contributions. As already detailed in Sect. 3.3, we have found that dark matter pair production plays no role in the exclusion.
The associated production of a mediator and a dark matter particle (XY ) has the largest impact on the ATLAS-EXOT-2016-27 exclusion, the analysis excluding the S1 model by solely using this component of the signal. This stems from an exclusive region in which Such a range corresponds to a phase-space region containing a significant fraction of the pp → XY events (see Fig. 3). The sensitivity to the S2 scenario, featuring a much heavier mediator (m Y = 1 TeV), is found to be slightly below 2σ when using NLO simulations. In contrast, LO predictions lead to too conservative conclusions, with a sensitivity barely reaching the 1σ level.
The pp → Y Y production cross section being smaller, the sensitivity of the ATLAS-EXOT-2016-27 analysis to this channel is expected to be reduced, although the finalstate objects that are typically reconstructed are significantly harder due to the production of two heavy mediators. In the S1 scenario, this effect is irrelevant as the mediator is light enough (m Y = 500 GeV) to be copiously pair-produced. The subsequent signal is hence excluded by the same signal region as the one defined in Eq. (19). Such a statement can however only be made after using NLO simulations (the LO rates being too small to reach a 95% CL exclusion) and when including not only the QCD-induced production mode, but also the dark matter t-channel exchange one. In the case of the S2 scenario, the signal regions are not populated enough to exclude the model. However, the yields are sufficiently large for driving an exclusion by considering both the Y Y and XY contributions, again provided NLO simulations are used.
We derive similar conclusions from the results obtained by recasting the ATLAS-SUSY-2016-07 analysis. This anal-ysis, that involves more complex cuts, better depicts the NLO impact on the shapes of the differential distributions. The corresponding modifications at the differential level indeed often lead to consider different most sensitive regions at LO and NLO.
With the examples worked out in this section, we have demonstrated the importance of relying on new physics precision simulations including NLO QCD predictions matched with parton showers. The correspondingly more precisely known total and differential cross sections allow for more robust conclusions on the sensitivity to the signal. The differences at the level of the distributions especially play a significant role in modifying the way in controlling how the different signal regions of the LHC analyses are populated. Moreover, it is crucial to consider all the components of a given signal, as their joint contribution may be sufficient to claim an exclusion, in contrast to the individual contributions taken separately.

Generalities
The studied t-channel simplified models are very peculiar as far as their dark matter phenomenology is concerned. While tree-level cross sections can be negligible, if not zero, NLO corrections or loop-induced processes might set up the stage. This is the case for any considered model restriction involving Majorana or scalar dark matter [23][24][25]34,61,62], while it is more model dependent for Dirac dark matter [26]. In the following, we focus on the fermionic dark matter case. In the early universe, the relic abundance is set by the annihilation ofχχ (Majorana) or χχ (Dirac) pairs (see the left diagram in Fig. 5), unless the mediator ϕ and the dark matter are within 20% in mass (r ≡ M ϕ /M χ 1.2). In this case, coannihilations [63] should be included as they dominate over a wide range of the parameter space (see the central and right diagrams in Fig. 5). Moreover, our analysis does not include Sommerfeld enhancement effects [64,65], as they are known not to alter the relic density predictions by more than 15% and only affect specific parts of the parameter space [25,[66][67][68].
We first consider, as in Sect. 3, a model restriction in which Dirac dark matter solely couples to the right-handed up quark (S3D_uR). In this model, both the dark matter spinindependent (SI) and spin-dependent (SD) elastic scattering cross sections off nucleons feature sizeable tree-level contributions stemming from s-channel mediator exchanges. When coannihilations are negligible (r 1.2), indirect detection rates stem from χχ annihilations into pairs of right-handed up quarks, which proceeds via s-wave t-channel mediator exchanges. The associated velocity-averaged cross section is about 3×10 −26 cm 3 /s for large λ ϕ values, and thus in the ball-park of the reach of the Fermi-LAT gamma-ray searches from dwarf spheroidal galaxies [69]. For illustrative purposes, we discuss, in the following, relic density and direct detection predictions. We refer to Refs. [26,30,32] for more comprehensive studies.
The Majorana dark matter restriction (S3M_uR) is similar to a supersymmetric model with bino-like neutralino dark matter and a right-handed up squark mediator. In this configuration, predictions for direct and indirect detection observables are dictated by NLO QCD corrections and loopinduced processes respectively. The direct detection SI elastic scattering cross section is negligible at tree level because of the Majorana nature of the dark matter, for which vectorial currents vanish. NLO QCD contributions at one loop, that include diagrams involving quarks and scalar mediators, therefore dominate and drive the scattering of dark matter off the nucleon constituents [27,28,[70][71][72][73]. The SD elastic scattering cross section is, on the contrary, dominantly dominated by tree-level contributions, and can be of the same order as the current experimental sensitivity.
Present dayχχ → uū annihilations in dense astrophysical environments are p-wave suppressed, as the tree-level s-wave contribution is proportional to the up-quark mass that vanishes in the chiral limit. There however exist two processes that could make Majorana dark matter detectable: virtual internal bremsstrahlung (VIB) in which the quark pair is produced together with a photon emitted by the internal t-channel propagator, and loop-induced annihilations into a photon pair or into a photon and a Z -boson. VIB yields a large correction to the tree-level annihilation cross section, uplifting the p-wave suppression by even a few order of magnitude, and provides a sharp spectral feature at the highest end of the gamma-ray spectrum (see, e.g., Refs. [23][24][25]61,67,68,[74][75][76]).
On the other hand, annihilations into photons have been known since a long time as the smoking gun to detect dark matter, as they produce monochromatic photons pinpointing the dark matter mass (see, e.g., Refs. [77][78][79][80][81][82]). Whilst these two processes are of higher order, the astrophysical background for a sharp gamma-ray spectral feature is very low. This yields a very good experimental sensitivity and annihilation cross sections well below the canonical 10 −26 cm 3 /s value can be probed for a wide range of dark matter masses [83]. Moreover, line searches by the HESS satellite [84,85] are sensitive to very heavy dark matter, with masses of tens of TeV, well above the sensitivity range of the LHC. This thus exhibits a nice complementarity with colliders.
The Majorana dark matter phenomenology briefly sketched here holds for scalar dark matter too. NLO processes even become relevant at freeze-out, the tree-level annihilation cross section being d-wave suppressed [61]. Similarly, any tchannel dark matter model restriction in which the dark mat- Fig. 5 Representative LO Feynman diagrams entering the relic density computation. We consider dark matter annihilations into quarks (left), as well as mediator annihilations (centre) and coannihilations (right) that are relevant for mediator and dark matter mass splittings of about 10-20% ter couples only to the third generation requires to account for QCD corrections already for the relic density predictions [62]. For instance, for all restrictions of the 3rd type, loop-induced dark matter annihilations into gluons turn out to be dominant and set the relic density below the b-quark threshold [86]. These corrections are typically not automatically included in available public software such as MadDM and MicrOMEGAs, and must be implemented following, e.g., Refs. [61,86].

Analysis setup and validation procedure
In our dark matter analysis, we impose the relic abundance for dark matter to match the value measured by the Planck satellite in 2018 [87]. The direct detection predictions are confronted with the exclusion bounds at 90% CL of the XENON1T [88] and of PICO-60 [89] experiments for the SI and SD cases respectively, and we display projections for the neutrino floor in our SI scattering results [90]. Loop-induced gamma-ray line predictions are compared with the Fermi-LAT [83] and HESS [84,85] line searches from the galactic centre, as recasted in Ref. [19] for an Einasto dark matter density profile [91] at 95% CL. We also show the projected sensitivity of the Cherenkov Telescope Array (CTA) [92] as obtained in Ref. [19] at 95% CL.
In order to compute the relevant observables with MadDM in the S3D_uR model restriction, we type in the command line interface of the programme, import model DMSimp_t-S3D_uR --modelname define darkmatter xd define coannihilator ys3u1 generate relic_density add direct_detection output my_project launch In the case of Majorana dark matter, the model name should be changed to DMSimp_t-S3M_uR, and the dark matter candidate name to xm (see Table 1). Scans can be performed by using standard MadDM syntax [6], and details for indirect detection calculations are provided in the next subsections. Representative Feynman diagrams contributing to the Fig. 6 Value of the λϕ coupling, as a function of the dark matter mass M χ , leading to a relic density of Ωh 2 = 0.12 in the S3D_uR model (top). We present MicrOMEGAs predictions for r = 1.1 (red) and 2 (blue). We moreover display the relative difference between MicrOMEGAs and MadDM for the relic density (dashed), the SI (solid green) and SD (solid magenta) scattering cross section off protons for r = 1.1 (centre) and 2 (bottom) thermally-averaged annihilation cross section σ v fo , assuming a standard dark matter freeze-out (fo), are depicted in Fig. 5.
To achieve our calculations in MadDM, we produce a LO UFO library in which all quarks are massive (unless stated otherwise). It differs from the NLO UFO library described in Sect. 2.1 in which all quarks except the top quark are massless. We moreover have additionally generated CalcHep model files [93], which is necessary to validate MadDM predictions obtained with the DMSimpt model against known results [23,24,61] derived with MicrOMEGAs. All model files are available from the URL http://feynrules.irmp.ucl.ac. be/wiki/DMsimpt.

Dark matter observables in the S3D_uR model
In Fig. 6, we validate our S3D_uR model implementation by numerically comparing relic density and direct detection predictions obtained with MadDM and MicrOMEGAs. In the upper panel, we derive the λ ϕ coupling value required to obtain the correct relic density as a function of the dark matter mass, for two choices of the mediator and dark matter mass ratio r .
In the more compressed scenario with r = f.1 (red curve), coannihilations and mediator annihilations are important, especially for small M χ . For M χ 200 GeV, the relic density is indeed mostly independent of λ ϕ , σ v fo being driven by pure QCD processes involving pairs of mediators annihilating into quarks and gluons (second diagram in Fig. 5). To properly evaluate these QCD processes, we include the running of the strong coupling in MadDM, 3 as implemented by default in MG5_aMC [94]. The number of quarks included in the loops depends on the running scale (and can be at most 5), and the QCD beta function can be evaluated at 1, 2 (default) and 3 loops. As far as the dark matter mass increases, processes involving both χ and ϕ become relevant so that λ ϕ has to be sizeable to obtain the right relic density. Already for M χ ∼ 250 GeV, X X annihilations (first diagram in Fig. 5) and XY coannihilations (third diagram in Fig. 5) contribute to the total scattering cross section σ v fo by about 45% and 30% respectively, the reminder being due to mediator-pair annihilations (Y Y ). As in the previous section, X = χ,χ and Y = ϕ, ϕ † in our notations. For dark matter masses larger than 500 GeV, X X annihilations dominate, the XY and Y Y processes contributing only to less than about 20% to the relic density.
In the r = 2 case (blue line), Y is too heavy relatively to dark matter to be relevant at freeze-out. Only X X annihilations contribute, and λ ϕ has to be sizeable and larger than for r = 1.1 for any a given M χ value. This large value compensates the smaller σ v fo cross section stemming from a smaller number of relevant processes than in the r = 1.1 scenario where coannihilations and mediator annihilations play a role.
In order to quantify the numerical differences between MicrOMEGAs and MadDM predictions for a given dark matter observable O, we define the quantity and focus, in Fig. 6, on the relic density (dashed), and the SI (solid green) and SD (solid magenta) direct detection cross sections. We present the dependence of ΔR on the dark matter mass both for the r = 1.1 (middle panel) and r = 2 (lower panel) scenarios. Predictions for both the SI and SD dark matter scattering cross section off protons 4 are found in perfect agreement, the discrepancy between MadDM and MicrOMEGAs being of at most a few percents for the probed dark matter mass range. Relic density predictions are also found to agree quite well, except for dark matter masses in the 100-200 GeV range for the r = 1.1 configuration. In this parameter space region, we get a discrepancy 3 This feature is now standard in the publicly available version of MadDM at the URL https://launchpad.net/maddm. 4 Similar results are obtained in the neutron case.

Fig. 7
Value of the λϕ coupling, as a function of the dark matter mass M χ , leading to a relic density of Ωh 2 = 0.12 in the S3M_uR model (top). We present MicrOMEGAs predictions for r = 1.1 (red) and 2 (blue). We moreover display the relative difference between MicrOMEGAs and MadDM for the relic density (dashed) for r = 1.1 (centre) and 2 (bottom). In the last case, the limits in which all quarks are massless are also presented (solid) reaching 5% to 15% due to the different treatment of the QCD sector in both codes (the relic density being driven by QCD-induced mediator annihilations). MicrOMEGAs indeed includes running quark masses, in addition to the strong coupling running [95].

Dark matter observables in the S3M_uR model
In this section, we focus on Majorana dark matter and derive, in Fig. 7, the values of the λ ϕ coupling that are needed to obtain the correct relic density for r = 1.1 (red) and r = 2 (blue) configurations. Comparing with the Dirac dark matter case, larger couplings are generally required as a consequence of the Majorana nature of dark matter, with the exception of setups featuring dark matter masses below 500 GeV where the relic density is driven by mediator annihilations and coannihilations. Another remarkable difference with Dirac dark matter, in the r = 1.1 scenario, is that there is no phenomenologically-viable solution for M χ 200 GeV.
In the middle and lower panel of the figure, we show that MadDM and MicrOMEGAs predictions agree quite well, as ΔR 5% for both scenarios (dashed lines), except for light dark matter where more important differences stem from the different treatment of the QCD sector. We addi- and r = 2 (blue), with a λ ϕ coupling yielding the right relic density. We compare MicrOMEGAs predictions (solid lines) with the analytical results of Ref. [28]. We additionally include the exclusion bounds extracted from current XENON1T results [88] (dashed black line) and the neutrino floor [90] (dashed yellow line) tionally assess the impact of the quark masses that induce a 10% shift (including running quark mass effects) relatively to the values obtained with MicrOMEGAs (for a massive quark setup).
In Fig. 8, we estimate the NLO SI dark matter elastic scattering cross section off protons for the r = 1.1 (red) and r = 2 (blue) scenarios. For each dark matter mass value, we fix the λ ϕ coupling to reproduce the relic density as observed by the Planck collaboration. We compare the NLO predictions obtained with MicrOMEGAs (solid curves) with the total SI cross-section (markers) given by the sum of the LO contribution obtained with MadDM with the analytically available NLO corrections from Ref. [28]. The analytic expression used here read where M p is the proton mass and the form factors f p and a p are functions of the Wilson coefficients describing the effective interactions with the proton components. These form factors are listed, for the various models, in Ref. [28]. In the case of our S3M_uR model, the SD form factor is given by where Δu p = 0.842 is the spin fraction of the up-quark in the proton. Concerning MicrOMEGAs the NLO SI contribution is automatically included following Ref. [70]. An excellent agreement is found. Confronting those results to the exclusion limits of XENON 1T [88], half of the viable parameter space (M χ 150 GeV) is excluded for both spectrum compression options. Moreover, this shows that most of the currently viable parameter space can be explored by nextgeneration dark matter experiments, as the corresponding SI scattering cross sections are larger than the expectation of the neutrino background [90].
Predictions for the LO SD elastic dark matter cross section off protons are shown in Fig. 9, in which we demonstrate the agreement between the numerical results of MadDM (solid lines) and the analytic expressions of Ref. [28] (markers). Confronting those predictions with the exclusion limits obtained from the PICO-60 experiment [89], it turns out that dark matter masses smaller than 500 GeV in the r = 2 configuration, and lying in the [250, 1200] GeV range in the r = 1.1 case, are disfavoured. Whereas the running of the λ ϕ coupling from the electroweak scale down to the GeV scale is known to largely enhance direct detection cross section predictions [34], this effect is not included neither in MadDM nor in MicrOMEGAs. This goes beyond the scope of this work.
We finally consider dark matter indirect detection in Fig. 10, in which we present predictions for the present time dark matter annihilation cross section in theχχ → uūγ and Predictions for the VIB (qqγ ) and photon-pair production (γ γ ) modes have been automatically computed with MadDM, and compared with the analytic expressions of Refs. [23] and [61] for the VIB and γ γ processes respectively. We additionally show constraints (dashed) from the Fermi-LAT dwarf spheroidal galaxies measurements and from the HESS experiment [19], as well as the expected sensitivity of the CTA experiment [19] χχ → γ γ channels for both considered benchmark scenarios. In both cases, the VIB cross section is larger than the diphoton one, although the latter loop-induced rate is of a similar order of magnitude as the former three-body one. 5 For more split spectra, or equivalently for larger r values, the loop-induced contributions are however known to dominate [24]. Confronting our predictions with the experimental results that are very sensitive to sharp features and lines in the gamma-ray spectra, we observe that there is no sensitivity to the two considered benchmark scenarios. This holds both for current exclusions extracted from the Fermi-LAT dwarf spheroidal galaxy measurements or HESS data, and for the expectation of the future CTA line search from the galactic centre.
In our predictions, we have compared the results of MadDM obtained by using the NLO DMSimpt UFO library of Sect. 2.1 (solid lines) with the analytical approximated expressions of Refs. [23] and [61] for the VIB and diphoton channels, respectively, 5 The extra factor of two accounts for the photon multiplicity.
As discussed in the manual [6], MadDM can automatically handle 2 → 3 generic annihilation processes by typing in, import model DMSimp_t-S3M_uR --modelname define darkmatter xm define coannihilator ys3u1 generate indirect_detection u u˜a output my_project launch The results of Fig. 10 however represents the first validation of a fully automated loop-induced process computation for a dark matter observable. 6 This feature will be available from the future version of MadDM. 7

Summary and conclusion
In this work, we have introduced the DMSimpt framework for dark matter t-channel models, available from the URL http://feynrules.irmp.ucl.ac.be/wiki/DMsimpt. This consists in a unique FeynRules implementation that allows for the calculation, through the various high-energy physics tools interfaced to FeynRules, of a large set of dark matter observables at colliders and in cosmology. The model is shipped with several restrictions relevant for simplified models featuring dark matter and coloured mediators of different spins. We have extensively studied two of those restrictions in which the dark matter is either a Dirac or a Majorana fermion, and the mediator is a scalar state coupling to the right-handed up quark.
We have generated a UFO model including ingredients for the automatic computation of collider observables matching NLO QCD predictions with parton showers. By a joint use of the MG5_aMC, Pythia 8, MadAnalysis 5, FastJet and Delphes 3 programmes, we have investigated the impact of the NLO corrections on various observables relevant for typical dark matter searches at the LHC through monojet probes, and shown how this could affect the sensitivity of the corresponding experimental searches. Our results emphasise the benefits of using NLO simulations to get more realistic predictions for total and differential cross sections including smaller theoretical systematics. We have moreover demonstrated how considering all new physics signals predicted by a given scenario as a whole is necessary for a better assessment of the LHC sensitivity to new phenomena. At the NLO accuracy, this however requires a specific treatment of s-channel resonant contributions usually appearing in the real emission contributions in order to avoid their double counting. Such a task can be automatically achieved within the MG5_aMC framework.
We have then made use of the MadDM programme for the automatic calculation of the dark matter relic density, spinindependent and spin-dependent scattering cross sections off nucleons and indirect detection rates. We have validated our predictions through a comparison with MicrOMEGAs (using a CalcHep model file generated from our general DMSimpt FeynRules implementation) and existing analytical calculations. Our predictions include both NLO corrections and the contributions of loop-induced processes as they could be dominant in specific model configurations, in particular for what concerns Majorana dark matter spinindependent direct detection (that is strongly affected by higher orders) and indirect detection (driven by loop-induced and virtual internal bremsstrahlung subprocesses). While MicrOMEGAs can account for corrections to direct detection, MadDM can automatically evaluate VIB processes. We have moreover pioneered the first automatic calculation of a loop-induced contribution to the production of gamma-ray lines by dark matter annihilations at present time.
In conclusion, our work presents, for the first time, a unified framework to undertake precision dark matter calculations in cosmology and at colliders for a large class of tchannel dark matter models.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theory article: there is no data attached with it.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .