Chiara Arina
email: chiara.arina@uclouvain.be

Mihailo Backović
email: mihailo.backovic@uclouvain.be

Eric Conte
email: eric.conte@iphc.cnrs.fr

Benjamin Fuks
email: fuks@lpthe.jussieu.fr

Jun Guo

Jan Heisig
email: heisig@physik.rwth-aachen.de

Benoît Hespel
email: benoit.hespel@uclouvain.be

Michael Krämer
email: mkraemer@physik.rwth-aachen.de

Fabio Maltoni
email: fabio.maltoni@uclouvain.be

Antony Martini
email: antony.martini@uclouvain.be

Kentarou Mawatari
email: kentarou.mawatari@lpsc.in2p3.fr

Mathieu Pellen
email: mpellen@physik.uni-wuerzburg.de

Eleni Vryonidou
email: eleni.vryonidou@uclouvain.be

A comprehensive approach to dark matter studies: exploration of simplified top-philic models

Keywords:

Studies of dark matter lie at the interface of collider physics, astrophysics and cosmology. Constraining models featuring dark matter candidates entails the capability to provide accurate predictions for large sets of observables and compare them to a wide spectrum of data. We present a framework which, starting from a model Lagrangian, allows one to consistently and systematically make predictions, as well as to confront those predictions with a multitude of experimental results. As an application, we consider a class of simplified dark matter models where a scalar mediator couples only to the top quark and a fermionic dark sector (i.e. the simplified top-philic dark matter model). We study in detail the complementarity of relic density, direct/indirect detection and collider searches in constraining the multi-dimensional model parameter space, and efficiently identify regions where individual approaches to dark matter detection provide the most stringent bounds. In the context of collider studies of dark matter, we point out the complementarity of LHC searches in probing different regions of the model parameter space with final states involving top quarks, photons, jets and/or missing energy. Our study of dark matter production at the LHC goes beyond the tree-level approximation and we show examples of how higherorder corrections to dark matter production processes can affect the interpretation of the experimental results.

JHEP11(2016)111

A Mediator width [START_REF] Craig | Searching for signs of the second Higgs doublet[END_REF] 1 Introduction

Evidence for the existence of dark matter (DM), although indirect, is quite convincing [START_REF] Bertone | Particle dark matter: evidence, candidates and constraints[END_REF][START_REF] Bertone | Particle dark matter: observations, models and searches[END_REF][START_REF] Drees | Mini-review of dark matter[END_REF]. Measurements of the cosmic microwave background and baryonic acoustic oscillations predict a dominant dark matter component in the matter budget of the Universe (in the framework of standard cosmology). In addition, detection of gravitational anomalies, such as the flattening of galaxy rotation curves and the presence of gravitational lensing in the absence of visible matter (e.g. the bullet cluster [START_REF] Clowe | A direct empirical proof of the existence of dark matter[END_REF]), strongly favours gravitational interactions of dark matter as plausible explanations.

JHEP11(2016)111

The many hints for dark matter sparked a huge endeavour to detect it and measure its properties, leading to a number of experiments and searches which exploit very different ideas and approaches to dark matter detection. The experiments can be broadly grouped into three categories:

• A wide range of underground nuclear recoil experiments aimed at detecting galactic dark matter scattering off atomic nuclei;

• Searches for dark matter annihilation in the galaxy or nearby dense sources via measurements of, for instance, gamma-rays;

• Collider searches in channels with large missing transverse energy ( / E T ).

However, despite an enormous experimental effort, the detection of the dark matter particles remains elusive. In fact, there is no clear indication that dark matter interacts with ordinary matter via forces other than gravity, and current experimental results are not able to put stringent bounds on the dark matter properties and couplings in a modelindependent way.

As so little is known about the true nature of dark matter, it is a useful strategy to try and constrain viable dark matter scenarios in the most model-independent way (i.e. via simplified models [START_REF] Alves | Simplified models for LHC new physics searches[END_REF]), confronting them with results from collider experiments, direct dark matter searches, astrophysical observations and cosmology. If or when a signal is observed, the aforementioned approach will help us to determine more accurately both the particle properties (mass, couplings, etc.) and astroparticle properties (halo properties, thermal relic density, etc.) of dark matter. Conversely, if searches result only in limits on dark matter parameters, combining constraints from different approaches aids us in excluding specific scenarios and hence narrow down the scope of viable dark matter theories.

Recent collider searches have focused mostly on studies of dark matter in the simplified model framework, where a single dark matter candidate of arbitrary spin couples to visible matter (e.g. quarks) via an s-channel or a t-channel mediator, whose quantum numbers are fixed by assumed local and global symmetries [START_REF] Abercrombie | Dark matter benchmark models for early LHC run-2 searches: report of the ATLAS/CMS dark matter forum[END_REF]. The minimal implementations of simplified dark matter models involve four basic parameters: the mass m X of the dark matter particle, the mass m Y of the mediator, the coupling constant g X of the dark matter to the mediator and the universal coupling g SM of the mediator to the visible sector (the width of the mediator is a derived quantity). Fast and efficient studies of the full simplified model parameter space require parameter scanning technology beyond simple sequential grids, due to the relatively high dimensionality of the parameter space. Past studies of simplified dark matter models have hence been limited to explorations of the parameter space in two-dimensional projections while keeping the remaining parameters fixed (see e.g. the works of refs. [START_REF] Buckley | Scalar simplified models for dark matter[END_REF][START_REF] Haisch | Simplified dark matter top-quark interactions at the LHC[END_REF][START_REF] Heisig | Constraints on Majorana dark matter from the LHC and IceCube[END_REF][START_REF] Bell | Mono-W dark matter signals at the LHC: simplified model analysis[END_REF][START_REF] Brennan | Collide and conquer: constraints on simplified dark matter models using Mono-X collider searches[END_REF][START_REF] Du Pree | Cosmological constraints on dark matter models for collider searches[END_REF][START_REF] Jackson | Gamma-ray lines and one-loop continuum from s-channel dark matter annihilations[END_REF][START_REF] Abdallah | Simplified models for dark matter searches at the LHC[END_REF][START_REF] Godbole | A simplified model for dark matter interacting primarily with gluons[END_REF][START_REF] Xiang | Searches for dark matter signals in simplified models at future hadron colliders[END_REF][START_REF] Difranzo | Simplified models for dark matter interacting with quarks[END_REF][START_REF] Harris | Closing up on dark sectors at colliders: from 14 to 100 TeV[END_REF][START_REF] Khoze | Spectroscopy of scalar mediators to dark matter at the LHC and at 100 TeV[END_REF][START_REF] Harris | Constraining dark sectors at colliders: beyond the effective theory approach[END_REF][START_REF] Zhang | Top quark mediated dark matter[END_REF] and the references therein).

In this paper we illustrate how comprehensive studies of simplified dark matter models can be performed, exploring their full four-dimensional parameter space while taking into account constraints from collider physics, astroparticle physics and cosmology. For concreteness, we focus on a class of simplified models where the dark matter dominantly couples via a scalar mediator to top quarks (i.e. 'top-philic dark matter' scenarios). Yet,

JHEP11(2016)111

the methodology we employ is general and can be applied to other scenarios as well. We provide detailed examinations of the two-dimensional projections of the full parameter space, and we demonstrate that striking features in the structure of the viable parameter space emerge through the combination of all current constraints. We also stress that in addition to collider searches for dark matter in channels with large missing energy, in this study we also consider resonance searches in channels with fully reconstructed final states, which can be useful to constrain the properties of the mediators.

We perform the study of simplified top-philic dark matter models by using a combination of simulation tools, including the MadGraph5 aMC@NLO (MG5 aMC shorthand) event generator [START_REF] Alwall | The automated computation of tree-level and next-to-leading order differential cross sections and their matching to parton shower simulations[END_REF], the FeynRules package [START_REF] Alloul | FeynRules 2.0 -a complete toolbox for tree-level phenomenology[END_REF][START_REF] Degrande | Automatic evaluation of UV and R 2 terms for beyond the Standard Model Lagrangians: a proof-of-principle[END_REF], the MadAnalysis 5 platform [START_REF] Conte | MadAnalysis 5, a user-friendly framework for collider phenomenology[END_REF][START_REF] Conte | Designing and recasting LHC analyses with MadAnalysis 5[END_REF][START_REF] Dumont | Toward a public analysis database for LHC new physics searches using MadAnalysis 5[END_REF], the Delphes 3 detector simulator [28] and the MadDM program [START_REF] Backović | MadDM v.1.0: computation of dark matter relic abundance using MadGraph5[END_REF][START_REF] Backović | Direct detection of dark matter with MadDM v.2.0[END_REF], together with an efficient parameter sampling technology based on the MultiNest algorithm [START_REF] Feroz | Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis[END_REF][START_REF] Feroz | MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics[END_REF]. We explore the full four-dimensional parameter space of the model in the light of existing collider and astroparticle constraints. Our analysis thus also represents a proof of concept for a unified numerical framework for comprehensive dark matter studies at the interface of collider physics, astrophysics and cosmology. This has direct implications for dark matter searches at colliders, as comprehensive phenomenological studies of dark matter models can be used to drive the experimental efforts towards the regions of the parameter space that are not already ruled out by astrophysical and cosmological constraints. In addition, we have also implemented previously unavailable experimental analyses into the MadAnalysis 5 platform, providing an added benefit of our work for future collider studies which go beyond searches for dark matter.

The article is organised as follows. Section 2 describes the details of the simplified topphilic dark matter model under consideration and discusses the constraints on the model parameter space that are implemented in our analysis setup. All cosmology and astrophysics constraints are discussed in section 3. More precisely, the relic density constraints are illustrated in section 3.1. We discuss the direct detection constraints in section 3.2, while constraints from gamma-ray flux measurements are detailed in section 3.3. Collider constraints are investigated in section 4. We study constraints from searches with and without missing transverse energy in section 4.1 and 4.2 respectively. Before concluding in section 6, we provide in section 5 a detailed discussion of the overall combined information coming from all the considered data. We provide more information on the mediator width in appendix A. As a validation of our calculations, we perform a detailed comparison between MadDM and MicrOMEGAs in appendix B. [START_REF] Bertone | Particle dark matter: evidence, candidates and constraints[END_REF], give details on the annihilation cross section of dark matter in the top-philic model in appendix B.2 and present the validation of the CMS t t + / E T and monojet implementation in the MadAnalysis 5 framework in appendix C.

Simplified top-philic dark matter model and its numerical implementation

The simplified top-philic dark matter model that we consider is constructed by supplementing the Standard Model (SM) with a Dirac-type fermionic dark matter candidate X and a

JHEP11(2016)111

scalar mediator Y 0 . The interactions of the two particles are described by the Lagrangian

L Y 0 t,X = -g t y t √ 2 tt + g X XX Y 0 , (2.1)
where the new physics interaction strengths are denoted by g t and g X for the mediator couplings to the Standard Model sector and to dark matter respectively. We have assumed a minimally flavour-violating [START_REF] D'ambrosio | Minimal flavor violation: an effective field theory approach[END_REF] scalar theory where the mediator couples to quarks with a strength proportional to the Standard Model Yukawa couplings, so that we neglect all light quark flavour couplings and only include the coupling of the mediator to the top quark, y t = √ 2m t /v where v = 246 GeV is the Higgs vacuum expectation value and m t is the top quark mass. The above Lagrangian is not invariant under the Standard Model gauge group. In ultraviolet completions where Y 0 is one of the components of an SU(2) L doublet, as in the Two-Higgs-Doublet Model, our simplified model corresponds to scenarios with a high degree of alignment. In this case, the couplings of the gauge bosons to the heavier CP -even scalar are suppressed, e.g. cos(βα) ∼ 0 where tan β consists of the ratio of the vacuum expectation values of the neutral components of the two Higgs doublets and α is the CP -even Higgs mixing angle (see e.g. [START_REF] Craig | Searching for signs of the second Higgs doublet[END_REF][START_REF] Carena | Impersonating the Standard Model Higgs boson: alignment without decoupling[END_REF]). Such setups are in particular common in minimal supersymmetric realisations. In contrast, if the mediator is a gauge singlet, it should mix with the Higgs sector, leading to a more complex phenomenology as in the so-called Higgs portal dark matter models (see e.g. the works of refs. [START_REF] Kim | Singlet fermionic dark matter[END_REF][START_REF] Baek | Search for the Higgs portal to a singlet fermionic dark matter at the LHC[END_REF][START_REF] Lopez-Honorez | Higgs portal, fermionic dark matter and a Standard Model like Higgs at 125 GeV[END_REF][START_REF] Baek | Invisible Higgs decay width vs. dark matter direct detection cross section in Higgs portal dark matter models[END_REF][START_REF] Baek | Beyond the dark matter effective field theory and a simplified model approach at colliders[END_REF][START_REF] Ko | Search for Higgs portal DM at the ILC[END_REF]). On the other hand, it has also recently been shown that the use of simplified models for LHC and future (feasible) collider studies does not break perturbative unitarity [START_REF] Englert | S-channel dark matter simplified models and unitarity[END_REF][START_REF] Barducci | Monojet searches for momentum-dependent dark matter interactions[END_REF].

The model contains four free parameters (two couplings and two masses),

{g t , g X , m X , m Y } , (2.2)
while the width Γ Y is fixed by the remaining model parameters. In addition to the Lagrangian of eq. (2.1), we could also have considered mediator couplings to leptons. They however cannot be well constrained by LHC searches and dark matter direct detection data, and we have excluded them from our model description. We will nonetheless comment on their relevance for relic density predictions and dark matter indirect detection signals in sections 3.1 and 3.3. The Lagrangian of eq. (2.1) induces dimension-five couplings of the mediator to gluons and photons via loop diagrams of top quarks. The loop-induced operators can be relevant in the context of both astrophysical and collider searches for dark matter. The couplings of the mediator to gluons and photons are given, at the leading order (LO), by the effective operators

L Y 0 g = - 1 4 g g (Q 2 ) v G a µν G a,µν Y 0 and L Y 0 γ = - 1 4 g γ (Q 2 ) v F µν F µν Y 0 , (2.3)
with the effective couplings being

g g (Q 2 ) = g t α s 3π 3 2 F S 4m 2 t Q 2 and g γ (Q 2 ) = g t 8α e 9π 3 2 F S 4m 2 t Q 2 .
(2.4)

JHEP11(2016)111

In the above expressions, Q 2 denotes the virtuality of the s-channel resonance, while F S is the one-loop form factor

F S (x) = x 1 + (1 -x) arctan 2 1 √ x -1 , (2.5)
with F S (x) → 2/3 for x 1. Eq. (2.4) contrasts with the Standard Model Higgs case where the effective Higgs-photon coupling receives contributions from vector-boson loopdiagrams that are absent in our simplified dark matter model setup. As a result, the gluon and photon effective couplings to Y 0 are characterised by a larger hierarchy compared to their Higgs counterparts.

The tree-level partial decay widths of the scalar mediator are given by

Γ(Y 0 → t t) = g 2 t 3y 2 t m Y 16π β 3 t Θ(m Y -2m t ) , (2.6)
Γ(Y 0 → X X) = g 2 X m Y 8π β 3 X Θ(m Y -2m X ) , (2.7)
where β t,X = 1 -4m 2 t,X /m 2 Y and Θ(x) is the Heaviside step function, and we ignored the top quark width in the expression for Γ(Y 0 → t t). The loop-induced Y 0 partial widths are

Γ(Y 0 → gg) = g 2 t α 2 s m 3 Y 72π 3 v 2 3 2 F S 4m 2 t m 2 Y 2 , (2.8)
Γ(Y 0 → γγ) = g 2 t α 2 e m 3 Y 81π 3 v 2 3 2 F S 4m 2 t m 2 Y 2 . (2.9)
The Y 0 partial width to photons is by construction always smaller than the partial decay width into a pair of gluons by virtue of α 2 s /α 2 e ∼ 100. In addition to a coupling suppression, other decay processes such as the loop-induced Y 0 decays into Zγ, ZZ and hh final states receive a kinematic suppression. Couplings of Y 0 to ZZ and hh could also appear at tree level in our model, but in the spirit of simplified models, we define them to be vanishing. In the following we hence safely approximate the total decay width for the mediator to be the sum of eqs. (2.6), (2.7) and (2.8).

The total decay width and the branching ratios of the mediator into t t, X X, gg and γγ final states are displayed in figure 1 for different choices of new physics couplings and masses. Light mediators with masses below the top-quark pair or the dark matter pair decay thresholds are narrow states, while above these thresholds, large Γ Y /m Y values are possible in particular for large couplings. For mediators with m Y m t , m X , the dominant decay channel is into a pair of gluons. In contrast, heavy mediators with mass m Y > m t , m X decay predominantly into pairs of top quarks and/or dark matter particles, where the exact details of the partial width values strongly depend on the masses and couplings. The branching ratio of Y 0 to photons is always suppressed, as argued above. We present in appendix A the dependence of the Γ Y /m Y ratio on the g t and g X couplings for different mass choices and on the m Y and m X masses for different coupling choices.

Our top-philic dark matter model can be probed in different ways including astrophysical and collider searches, as listed in table 1. The relative importance of the various searches depends on the hierarchy of the dark matter, mediator and top-quark masses, as well as on the hierarchy between the couplings. Starting with the dark matter relic density, the annihilation cross section is dominated by subprocesses with top-quark final states for m X > m t , and by annihilation into gluons and to a lesser extent photons for light dark matter particles with m X < m t . If the mediator is lighter than the dark matter state, an additional annihilation channel into a pair of mediators can open up. The annihilation mechanisms into top-quarks, gluons/photons and mediators moreover provide an opportunity to indirectly search for dark matter, e.g. in gamma-ray data. The interactions of the dark matter particles with nuclei, relevant for direct detection experiments, proceed via mediator exchanges. The mediator-nucleon coupling is in turn dominated by the scattering off gluons through top-quark loops.

Dark matter production at the LHC proceeds either through the production of the mediator in association with top quarks, or from gluon-fusion through top-quark loops. Searches at the LHC can be classified into two categories regarding the nature of the final states that can contain missing transverse energy / E T or not. Searches involving missing energy may include final state systems containing a top-quark pair and probe in this way the associated production of a top-antitop-mediator system where the mediator subsequently decays into a pair of dark matter particles. Alternatively, the mediator can be produced via gluon fusion through top-quark loops, where the probe of the associated events consists of tagging an extra radiated object. This yields the well-known monojet, mono-Z and mono-Higgs signatures. We do not consider the monophoton channel, as photon emission is forbidden at LO in our simplified model by means of charge conjugation invariance. The second search category is related to final states without any missing energy, i.e. when the mediator decays back into Standard Model particles. This includes decays into top-quarks, leading to final states comprised of four top quarks, into a top-quark pair, as well as into a dijet or a diphoton system via a loop-induced decay. This is, however, relevant only for on-shell (or close to on-shell) mediator production. We proceed with a description of the numerical setup for our calculations. In the following sections, we explore the full four-dimensional model parameter space and present results in terms of two-dimensional projections. We perform the four-dimensional sampling using the MultiNest algorithm [START_REF] Feroz | Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis[END_REF][START_REF] Feroz | MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics[END_REF], where we assume Jeffeys' prior on all the free parameters in order not to favour a particular mass or coupling scale. The choice of prior ranges for the parameters is summarised in table 2, in which we have chosen to limit the coupling values to a maximum of π to ensure perturbativity. We implement the relic density constraints into MultiNest using a Gaussian likelihood profile, while for the direct detection limits we assume a step likelihood function smoothed with half a Gaussian. In addition, the sampling imposes that the model is consistent with values of Γ Y such that the mediator Y 0 decays promptly within the LHC detectors. Table 3 summarises the constraints that we have imposed on the model parameter space.

JHEP11(2016)111

m X > m t Cosmology relic indirect m X < m t Planck, FermiLAT Astrophysics m X > m Y direct m X > 1 GeV LUX, CDMSLite Colliders / E T m Y > 2m X +t t m Y > 2m X +j, +Z, +h no / E T m Y > 2m t 4t m Y > 2m t t t m Y < 2m X , 2m t jj, γγ

JHEP11(2016)111

MultiNest parameter Prior Throughout our study, we assume that X is the dominant dark matter component, namely that it fully accommodates a relic density Ω DM h 2 as measured by the Planck satellite [START_REF] Ade | Planck 2015 results. XIII. Cosmological parameters[END_REF]. Concerning the direct detection of dark matter, we consider the currently most stringent bounds on the spin-independent (SI) nucleon-DM cross section as measured by LUX for dark matter with m X > 8 GeV [START_REF] Akerib | First results from the LUX dark matter experiment at the Sanford Underground Research Facility[END_REF] and by CDMSLite for 1 GeV< m X < 8 GeV [START_REF] Agnese | New results from the search for low-mass weakly interacting massive particles with the CDMS low ionization threshold experiment[END_REF]. In section 3.3, we focus on indirect detection constraints that are imposed on the basis of the gamma-ray measurements achieved by the Fermi-LAT telescope [START_REF] Ackermann | Searching for dark matter annihilation from milky way dwarf spheroidal galaxies with six years of Fermi Large Area Telescope data[END_REF][START_REF] Ackermann | Updated search for spectral lines from galactic dark matter interactions with pass 8 data from the Fermi Large Area Telescope[END_REF]. Those bounds are however not applied at the level of the likelihood function encoded in our MultiNest routine, and we have chosen instead to reprocess the scan results for those parameter points that are consistent with both the relic density and direct detection considerations. For the purpose of the relic density and direct detection cross section calculations, we utilise both the MadDM [START_REF] Backović | MadDM v.1.0: computation of dark matter relic abundance using MadGraph5[END_REF][START_REF] Backović | Direct detection of dark matter with MadDM v.2.0[END_REF] and MicrOMEGAs [START_REF] Bélanger | MicrOMEGAs4.1: two dark matter candidates[END_REF] numerical packages, although we only present the results obtained with MadDM. The consistency checks that we have performed with both codes are detailed in appendix B.1.

log(m X / GeV) 0 → 3 log(m Y / GeV) 0 → 3.7 log(g X ) -4 → log(π) log(g t ) -4 → log(π)
We derive collider constraints on the simplified top-philic dark matter model using the MG5 aMC [START_REF] Alwall | The automated computation of tree-level and next-to-leading order differential cross sections and their matching to parton shower simulations[END_REF] framework and the recast functionalities of MadAnalysis 5 [START_REF] Conte | MadAnalysis 5, a user-friendly framework for collider phenomenology[END_REF][START_REF] Conte | Designing and recasting LHC analyses with MadAnalysis 5[END_REF][START_REF] Dumont | Toward a public analysis database for LHC new physics searches using MadAnalysis 5[END_REF] (where appropriate). We apply the LHC constraints on the top-philic dark matter model with two different procedures. On one side, similarly to what has been performed for the indirect detection bounds, we reprocess the scenarios that accomodate the observed relic density and that are compatible with LUX and CDMSLite data. However, we also study the collider bounds on the parameter space independently of any astrophysics and cosmology consideration and by relaxing the narrow width requirement (allowing Γ Y /m Y to be of O(1)) as well. In order to increase the sensitivity of the LHC searches, we allow

JHEP11(2016)111

for wider coupling ranges of 10 -2 < g X < 2π and 10 -2 < g t < 2π. The collider study without any cosmological and astrophysical constraint therefore includes the cases where the dark matter is not a standard thermal relic (i.e. its relic density is a result of a nonthermal mechanism or a non-standard evolution of the Universe). Details are provided in section 4 and appendix C for what concerns the validation of the CMS analyses that we have implemented in MadAnalysis 5 for this work.

In conclusion to this section, we point out that even though our current work focuses on a dark matter candidate which is a Dirac fermion, a more general implementation of simplified dark matter models in FeynRules [START_REF] Alloul | FeynRules 2.0 -a complete toolbox for tree-level phenomenology[END_REF][START_REF] Degrande | Automatic evaluation of UV and R 2 terms for beyond the Standard Model Lagrangians: a proof-of-principle[END_REF] can also account for pseudoscalar mediators as well as for CP -mixed states and for dark matter particles which are real or complex scalars [START_REF] Backović | Higher-order QCD predictions for dark matter production at the LHC in simplified models with s-channel mediators[END_REF][START_REF] Mattelaer | Dark matter production through loop-induced processes at the LHC: the s-channel mediator case[END_REF][START_REF] Neubert | Higher-order QCD predictions for dark matter production in mono-Z searches at the LHC[END_REF]. The corresponding model files have been used in this work and can be downloaded from the FeynRules model repository [START_REF]Simplified dark matter models webpage[END_REF] that also includes a model where the mediator is a spin-1 state that couples to either a fermionic or a scalar dark matter candidate [START_REF] Backović | Higher-order QCD predictions for dark matter production at the LHC in simplified models with s-channel mediators[END_REF]. All the models allow for the automated calculation of next-toleading-order (NLO) effects and loop-induced leading-order (LO) processes in QCD in the context of LHC predictions.

Cosmological and astrophysical constraints

We begin our analysis of the simplified top-philic dark matter model with a detailed discussion of the cosmological and astrophysical constraints.

Constraints from dark matter relic density

Dark matter annihilation in the early Universe is determined, in the simplified top-philic dark matter model, by a combination of three processes,

X X → t t (I) , X X → gg (II) , and
X X → Y 0 Y 0 (III) ,
where we have omitted the annihilation into photons as it is always suppressed compared to the annihilation into gluons. The analytic expressions for the thermally averaged annihilation cross section in the non-resonant region σv rel corresponding to each of the processes listed above are provided in appendix B.2. The first two processes proceed via an s-channel Y 0 exchange (first two rows of table 1), while the third process consists of a t-channel X exchange (third row of table 1). The resonance structure of the s-channel processes implies that the width of Y 0 potentially plays an important role in the determination of the relic density assuming a dominant annihilation via the processes (I) and (II), while the effects of the Y 0 width are mostly negligible if the annihilation dominantly proceeds via the t-channel X exchange process (III).

According to the hierarchy between the dark matter mass m X , the mediator mass m Y and the top quark mass m t , different situations can occur. Qualitatively, one expects that:

• for m Y m X m t : process (I) is dominant as the tree-level annihilation into a pair of top quarks is kinematically allowed, the annihilation into gluons being loop suppressed, and the one into a pair of mediators kinematically suppressed;

JHEP11(2016)111

• for m X m t , m Y : dark matter annihilates into a pair of gluons as in process (II), since it is the only kinematically allowed channel;

• for m t m X m Y : relic density is determined by process (III) since annihilation into top quarks is kinematically forbidden and the one into gluons occurs away from the resonant pole of m Y ;

• for m X > m t , m Y and m Y < 2m t : similarly to the case above, the dominant annihilation mechanism is process (III), as annihilation into top quarks occurs far from the resonant pole and is suppressed kinematically;

• for m X > m t , m Y and m Y > 2m t : processes (I) and (III) are competitive and the dominant process among the two is determined by the hierarchy between the g t and g X couplings.

Requiring our simplified top-philic dark matter model to result in a dark matter relic density consistent with the most recent Planck measurements [START_REF] Ade | Planck 2015 results. XIII. Cosmological parameters[END_REF] implies strong constraints on the viable regions of the parameter space. As an illustration, we consider the region of the parameter space in which m t m X m Y , where we expect the dominant annihilation mechanism of dark matter to be process (III) and to give rise to a pair of mediators. In this region, the thermally averaged annihilation cross section approximately reads

σv rel ann ∼ g 4 X m 2 X ∼ 10 -9 GeV -2 , (3.1)
so that it is clear that imposing that the relic density predictions agree with Planck data leads to a stringent constraint on the ratio g 2 X /m X . The argument is more involved in parameter space regions where the total mediator width Γ Y plays a role, as the relevant quantity involved in the relic density calculation is in general not σv rel ann but dx σv rel ann (x) where x ≡ m X /T and σv rel ann is a non trivial function of x. This is especially true, for instance, for the Breit-Wigner-type amplitudes that appear in processes (I) and (II).

In order to provide a more detailed quantitative analysis, we have performed a fourdimensional scan the top-philic dark matter model parameter space and examined the effects of imposing relic density constraints on the allowed/ruled out parameter sets. Figure 2 reveals the rich structure of the four-dimensional parameter space allowed by relic density measurements. The bulk of the allowed parameter points lies in the region where m X > m Y , and the annihilation cross section is dominantly driven by process (III). This region of the parameter space has the particularity of not being reachable by traditional monojet, monophoton, mono-Z and mono-Higgs searches at colliders. The decay of the mediator into a pair of dark matter particles is indeed not kinematically allowed, so that any new physics signal will not contain a large amount of missing energy. The model can however be probed at colliders via dijet, diphoton, t t (plus jets) and four-top analyses. We elaborate on this point more in section 4.2. The characteristic mediator width Γ Y in this region tends to be extremely small, with values of at most 10 -4 GeV as shown in the top left panel of figure 2. This is expected as the width is mostly controlled by the decays into JHEP11(2016)111 gluons, and into top quarks in the regions where this decay is kinematically allowed, the decay into a pair of dark matter particles being forbidden.

In the region where m X m t and m Y 2m t , the mediator decay into a t t final state is kinematically allowed and the dark matter annihilation cross section is driven by the X X → Y 0 → t t process. The only other parameter space region that is not ruled out by the relic density data is centered around the resonance region where m Y ∼ 2m X . The extension of the region away from the resonance pole is due to the Y 0 width that can reach O(10) GeV. The resonant region extends to lower m X and m Y values, and is the only allowed region when both m X and m Y are smaller than m t (but with Y 0 decays into a pair of dark matter particles being allowed). This has interesting implications for LHC searches as the low dark matter/mediator mass region is the one where colliders have the best sensitivity, in particular through monojet searches (see for instance section 4.1.2).

-11 -JHEP11(2016)111 The structure of the ruled out parameter space regions shows several other interesting qualities. The most striking feature is that almost the entire region where m Y 2m X does not lead to predictions of a dark matter relic density in agreement with the observations. There are also no allowed points for m X m t , except very close to the resonance line. This region is characterised by a dominant mediator decay into gluons, which results in typical Γ Y /m Y 1, a small total dark matter annihilation cross section, and hence an overproduced dark matter. The upper limit imposed on the size of the couplings (see table 2) is largely responsible for the absence of allowed points in the region. For instance, taking any m X value so that the predicted relic density agrees with the observed value, an increase in m Y will result in a decrease of the annihilation cross section, in turn leading to a higher relic density. The only way (away from the resonance) to restore the correct relic density is then to increase the size of g X and/or g t . However, our results show that even for couplings of O(1), the cross section in this region is too small not to overproduce dark matter.

The region of parameter space between m Y ∼ m X and m Y ∼ 2m X is consistent with the above-mentioned argument. This strip of the ruled out parameter space can be seen as a part of the larger ruled out region for which m Y m X . Tuning m X to be close to m Y /2 and assuming a relatively small Γ Y value is the only way to enhance the dark matter annihilation cross section and not overproduce dark matter.

In addition to projections of the allowed parameter space onto the (m Y , m X ) plane, we have also studied several other projections. Figure 3 shows the projections of our results onto the (g X , g t ) plane (left) and (m X , g t ) plane (right), where we show m X and m Y as a colourmap in the first and second panel respectively. Regardless of the value of m X and m Y in the considered scan range, there are no solutions for g X and g t which satisfy the JHEP11(2016)111 relic density constraint in the region where g t 10 -2 and g X 10 -1 . This finding is consistent with the left lower panel of figure 2 where we have found that a correct relic density favours g X couplings of O(1). Furthermore, we can observe that only couplings of O(1) result in Γ Y O(1) GeV, while in the majority of the allowed (g X , g t ) parameter space regions Γ Y /m Y 1. We find no striking features in the (m X , g t ) projection of the scanned parameter space. The unpopulated regions in the lower left corners are artifacts of the lower limit on the coupling size of 10 -4 .

As a validation, we have cross checked our calculations with the MicrOMEGAs code. The results obtained with MadDM and MicrOMEGAs agree in most of the parameter space, except in the region where g t and m X are small. Some numerical discrepancies are expected to occur in this region, as shown in appendix B.1 and by comparing figures 2 and 19.

As a last remark, allowing the scalar mediator to couple to all quarks and leptons would only have a minor impact on our results. The region dominated by the process (III) will indeed stay unchanged, since it is insensitive to the coupling between Y 0 and the Standard Model fermions. As far as it concerns dark matter annihilation via an s-channel Y 0 exchange, one would have to sum up over all the possible final states kinematically open. This would increase the total annihilation cross-section and decrease Ω DM h 2 , implying that the constraint of having Ω DM h 2 ∼ 0.12 leads to a rescaling of all fermionic couplings towards smaller values with respect to the g t values shown in this work. The major difference would reside in a potentially larger decay width for Y 0 and hence wider "bands" of the resonance regions of the allowed parameter space.

Constraints from direct detection

Simplified models of dark matter which feature couplings to quarks and gluons can also be bounded by results from underground direct dark matter detection experiments. In topphilic dark matter scenarios, dark matter scatters off nucleons via the t-channel exchange of Y 0 , where the scattering off gluons via triangle top loops accounts for the dominant contribution to the DM-nucleon scattering rate.

The spin independent (SI) dark matter-nucleon cross section is given by

σ n SI = 4 π m X m n m X + m n 2 2 27 m n m t g X g t m 2 Y f G 2 , (3.2)
where f G ≡ 1 -q≤3 f q = 0.921 [START_REF] Beringer | Review of particle physics (RPP)[END_REF][START_REF] Alarcon | The chiral representation of the πN scattering amplitude and the pion-nucleon sigma term[END_REF] 1 is the gluon form factor and the sum runs over the light quarks q = u, d, s, where m n ≈ 0.938 GeV is the nucleon mass and m t = 173 GeV is the top quark mass. The expression in eq. (3.2) does not depend on Γ Y , simplifying the constraints which can possibly be derived from direct detection. For instance, considering a scenario in which generic m X and m Y masses are fixed and where the dominant annihilation process is process (I), direct detection directly constraints the product g X g t . Extracting the constraint on this quantity in a generic fashion is much more complicated in the case JHEP11(2016)111 GeV and accommodate the direct detection constraints (cf. table 3).

of dark matter annihilation in the early Universe and at colliders, as the processes involved in dark matter relic density and dark matter production calculations intrinsically depend on a quantity which is proportional to g X g t /Γ Y .

The running of the g X and g t couplings could have an effect on the value of the spin independent DM-nucleon scattering cross section [START_REF] D'eramo | Connecting dark matter UV complete models to direct detection rates via effective field theory[END_REF][START_REF] Vecchi | WIMPs and un-naturalness[END_REF]. However, a proper inclusion of the running couplings would require a careful treatment of the renormalisation group evolution via multiple energy scales which is beyond the scope of our current effort. Instead, we restrict here our calculations to constant g X and g t values. The effect of the running couplings would then be equivalent to a rescaling of g X and g t to different values.

Next, we have repeated the four-dimensional parameter scan from section 3.1 including into the MultiNest likelihood function also bounds stemming from direct detection. Figure 4 shows the results of the scan, projected onto three different planes, where we removed -14 -

JHEP11(2016)111

the points excluded by the 95% confidence limit (CL) bound from LUX and CDMSLite. Direct detection rules out a major portion of (m Y , m X ) space allowed by the relic density constraints (regardless of the coupling value) in the region where m X m Y , where collider bounds are irrelevant. Figure 4 hence serves as a good example for the complementarity among direct detection, relic density and collider bounds. In the (g X , g t ) plane, direct detection does not rule out a well defined-region (top-right panel of figure 4), indicating that for any pair of couplings (g X , g t ) in the range of [10 -4 , π] which are allowed by the relic density constraint, it is always possible to find a pair of (m Y , m X ) values which are not ruled out by direct detection data. In the (m X , g t ) projection, we finally observe that direct detection rules out a well-defined portion of the parameter space. Furthermore, the constraint also rules out small width points for g X 0.1 and m X m t . Direct detection bounds are indeed more sensitive to dark matter masses in the ballpark of 10 to 200 GeV and quickly deteriorate at larger dark matter masses, since the event rate in the detector scales as 1/m2 X . We also see that the direct detection exclusion limit is able to rule out a large portion of the parameter space where Y 0 is light, below 30 GeV, while the sensitivity is quickly lost for heavier masses of the scalar mediator. This can be understood by the 1/m 2 Y dependence of the SI elastic cross section of eq. (3.2). Both mass dependences are illustrated by the lower panel of figure 4.

Constraints from indirect detection

Top-philic dark matter annihilation in the present Universe could result in fluxes of cosmic rays and prompt gamma-rays, which can also be used to infer useful limits on the model parameter space. The annihilation of a X X pair in the galactic halo (or in dense environments of galactic centers) and the subsequent production of a secondary gamma ray flux is dictated by the same processes (I), (II) and (III) that set the relic abundance. These processes give rise to a continuum of secondary photons due to the decay and subsequent QED showering of the pair-produced top quarks, gluons and/or mediators. As already mentioned in section 2, a direct coupling of the mediator to a pair of prompt photons is induced at higher order in perturbation theory via a loop of top quarks. Hence, analogously to process (II), the process X X → γγ exists and yields the production of two monochromatic photons that could be detected in searches for lines in the gamma-ray spectrum. 2 Similarly to the relic density case, measurements of the gamma-ray fluxes can potentially constrain the coupling g X for the t-channel process (III) or the product of couplings g X g t in the case of an s-channel annihilation via the processes (I) and (II). However, it is important to highlight the differences between factors which are constrained by the dark matter relic density and by its indirect detection. The relic density is an integrated result over the thermal history of the Universe. Hence, the width of the resonance is important,

even if |m Y -2m X |
Γ Y (except in the case where m Y 2m X ). Conversely, the characteristic velocity of the dark matter particles today is of the order of v ∼ 10 -3 , implying highly non-relativistic dark matter annihilation. The width of the mediator in an s-channel JHEP11(2016)111 . Dark matter annihilation cross section at present time that is relevant for gamma-ray limits extracted from dwarf spheroidal galaxies measurements (left) and gamma-ray line searches (right). We show a maximal estimate of (σv rel ) tot and (σv rel ) γγ obtained by choosing v rel = 2v ∞ , where v ∞ is the escape velocity for dwarf spheroidal galaxies and the galactic center, respectively. All represented points feature a relic density in agreement with Planck data, Γ Y /m Y ≤ 0.2, Γ Y > 10 -11 GeV and accommodate the direct detection constraints (cf. table 3).

dark matter annihilation process is hence relevant for indirect detection only in the case of

|m Y -2m X | Γ Y .
Searches for gamma-ray signals of dark matter annihilation weakly constrain our simplified top-philic dark matter model. We have investigated results from gamma-ray line searches in the inner galactic region [START_REF] Ackermann | Updated search for spectral lines from galactic dark matter interactions with pass 8 data from the Fermi Large Area Telescope[END_REF], as well as continuum gamma-ray measurements from dwarf spheroidal galaxies [START_REF] Ackermann | Searching for dark matter annihilation from milky way dwarf spheroidal galaxies with six years of Fermi Large Area Telescope data[END_REF] and found no meaningful exclusion of the parameter space once the relic density and direct detection constraints are imposed. The lack of additional useful bounds is expected, as the annihilation of dark matter in the present Universe is p-wave suppressed, i.e. σv rel ∝ v 2 rel for all three annihilation channels (see appendix B.2 for more detail). This contrasts with scenarios in which the mediator is a pseudoscalar state that implies that the p-wave suppression at low dark matter velocity is only present for process (III), so that the gamma-ray constraints should be significantly stronger.

The gamma-ray line searches constrain the velocity-averaged cross section for the direct dark matter annihilation into two photons. Due to its p-wave suppression, this quantity is very sensitive to the choice of the velocity distribution of the dark matter in the galaxy which is subject to large uncertainties (see e.g. ref. [START_REF] Mao | Halo-to-halo similarity and scatter in the velocity distribution of dark matter[END_REF]). We adopt a conservative viewpoint here, evaluating the annihilation cross section at the highest possible velocity v rel = 2v ∞ with v ∞ being the escape dark matter velocity for our galaxy which we take to be v ∞ = 550 km/s [START_REF] Lisanti | The dark matter at the end of the galaxy[END_REF]. The left panel of figure 5 shows the respective result for (σv rel ) γγ . The limits from gamma-ray line searches lie between 2 × 10 -32 cm 3 s -1 (for dark matter masses around 1 GeV) and 4 × 10 -28 cm 3 s -1 (for dark matter masses around 500 GeV).

Photons arising from process (III) when the mediator subsequently decays into two photons provide a box-shaped gamma-ray spectrum, which, for a sufficiently compressed spectrum, can also resemble a line signal [START_REF] Ibarra | Dark matter constraints from box-shaped gamma-ray features[END_REF][START_REF] Ibarra | Gamma-ray boxes from axion-mediated dark matter[END_REF]. However, in analogy to X X → γγ, this process is both p-wave suppressed and coupling (8α 2 e /9α 2 s ) suppressed. Gamma ray

JHEP11(2016)111

line searches hence do not provide any further constraint on the parameter space in this compressed region. Searches for gamma-ray signals in dwarf spheroidal galaxies constrain the total the annihilation cross section at (two times) the escape velocity, the escape velocity of the considered dwarf spheroidal galaxies being typically much smaller and of the order of 10 km/s [START_REF] Walker | A universal mass profile for dwarf spheroidal galaxies[END_REF], which leads to a heavy suppression of the dark matter annihilation cross section. The right panel of figure 5 shows the annihilation cross section evaluated for v ∞ = 50 km/s. The cross sections are much smaller than the constraints which are around 10 -26 cm 3 s -1 , the exact details depending on the dark matter mass and the relevant annihilation processes.

In cases where we would have allowed for leptonic couplings of the scalar mediator Y 0 , our general conclusion about the poor ability of indirect dark matter searches to constrain the model parameter space remains unchanged. Dark matter annihilation into leptonic final states could give rise to additional continuum gamma-ray or positron fluxes, but the overall normalisation of σv rel would not change significantly and remain four to five orders of magnitude below the current bounds. Even under the most aggressive assumptions, all obtained bounds would still be far from being able to constrain a top-philic dark matter model with scalar mediators.

Collider constraints

As discussed in section 2, simplified top-philic dark matter scenarios can be probed at colliders through the production of the mediator either in association with a top-quark pair or through a top-quark loop. Depending on the mass and coupling hierarchy, the mediator decays either into a pair of dark matter particles, which results in signatures including missing transverse energy ( / E T ), or into Standard Model final states. The size of the cross sections associated with these two classes of mediator production mechanisms is depicted in figure 6 where we present their dependence on the mediator and dark matter masses m Y and m X . For the case where the mediator is singly produced, we use the Higgs cross section values that are reported in the Higgs Cross section Working Group documentation [START_REF] Andersen | Handbook of LHC Higgs cross sections: 3. Higgs properties[END_REF] and that are evaluated at the next-to-next-to-leading order (NNLO) accuracy in QCD. For all the other cases, the hard-scattering cross section is convoluted with the NNPDF 2.3 [START_REF] Ball | Parton distributions with LHC data[END_REF] set of parton distribution functions (PDF) within MG5 aMC, the PDFs being accessed via the LHAPDF library [START_REF] Whalley | The Les Houches accord PDFs (LHAPDF) and LHAGLUE[END_REF][START_REF] Buckley | LHAPDF6: parton density access in the LHC precision era[END_REF]. We employ a five-flavournumber scheme, and leading-order (LO) and next-to-leading-order (NLO) PDFs are used where relevant. The renormalisation and factorisation scales are set to half the sum of the transverse mass of all the final-state particles both for LO and NLO calculations, and the scale uncertainty is estimated by varying the two scales independently by a factor of two up and down. Additional details on the calculation of the Y 0 t t cross section are provided in ref. [START_REF] Backović | Higher-order QCD predictions for dark matter production at the LHC in simplified models with s-channel mediators[END_REF] while loop-induced processes are extensively documented in ref. [START_REF] Mattelaer | Dark matter production through loop-induced processes at the LHC: the s-channel mediator case[END_REF].

All the cross sections shown in figure 6 are proportional to g 2 t and we therefore arbitrarily choose g t = 1 as a benchmark. In this case, sizeable cross sections of 10 8 TeV LHC production mode. Requiring an extra hard jet in the final state reduces the cross section by a factor which depends on the missing energy (or the jet transverse momentum p T ) selection, and the production rates are not sensitive to the mediator mass as soon as the latter is smaller than the / E T selection threshold. The cross sections for producing the mediator in association with a Standard Model Higgs or Z boson are further suppressed. In contrast, the cross section related to the production of the mediator in association with a top-quark pair is significant for light mediators, but falls off quickly with the increase in the mediator mass due to phase-space suppression. As a result, a change in the collider energy from 8 to 13 TeV is important for heavy mediators and the cross section can be enhanced by about an order of magnitude. In the right panel of figure 6, we further show first that the cross sections are constant when the dark matter particle pair is produced through the decay of an on-shell mediator, and next that they are considerably suppressed when the mediator is off-shell, especially for the t tX X channel.

Y 0 tt (NLO) 2m t Y 0 j (LO p T Y 0 >300) Y 0 Z (LO p T Y 0 >150) Y 0 h (LO p T Y 0 >150) g t = 1
As already mentioned, the collider searches which provide the most relevant constraints on simplified top-philic dark matter models are based on the production channels shown in figure 6 and can in general be divided into two categories. The first category involves signals with missing transverse energy originating from the production of dark matter particles that do not leave any trace in the detectors and that are accompanied by one or more Standard Model states. The most relevant searches of this type are the production of dark -18 - matter in association with a top-quark pair and the loop-induced production of dark matter in association with a jet, a Z boson or a Higgs boson. This is discussed in section 4.1.

JHEP11(2016)111

The second category of searches relies on Y 0 resonant contributions to Standard Model processes. In our scenario, dijet, diphoton, top-pair and four-top searches are expected to set constraints on the model parameter space. This is discussed in section 4.2. As shown below, missing-energy-based searches and resonance searches are complementary and necessary for the best exploration of the model parameter space at colliders. In the rest of this section, we study collider constraints independently from the cosmological and astrophysical ones, and we dedicate section 5 to their combination. We moreover allow the mediator couplings to be as large as 2π and do not impose any constraint on the mediator width over mass ratio. We summarise the relevant 8 TeV LHC constraints used in this study in table [START_REF] Clowe | A direct empirical proof of the existence of dark matter[END_REF] and give details on the t t + / E T and monojet searches that have been recast in the MadAnalysis 5 framework in appendix C. Dark matter production in association with a top-quark pair (t t + / E T ) has been explored by both the ATLAS [START_REF]Search for dark matter in events with heavy quarks and missing transverse momentum in pp collisions with the ATLAS detector[END_REF] and CMS [START_REF]Search for the production of dark matter in association with top-quark pairs in the single-lepton final state in proton-proton collisions at √ s = 8 TeV[END_REF] collaborations within the 8 TeV LHC dataset, and limits have been derived in particular in the effective field theory approach [START_REF] Cheung | The top window for dark matter[END_REF][START_REF] Lin | Probing dark matter couplings to top and bottom quarks at the LHC[END_REF]. Such analyses could however be used to derive constraints in other theoretical contexts, and we choose to recast the CMS search to constrain the parameters of the simplified top-philic dark matter model under consideration. In this work, we simulate t tX X events at the NLO accuracy in QCD by making use of MG5 aMC. The first study of the genuine NLO effects on the production of a system composed of a pair of top quarks and a pair of dark matter particles has been presented in ref. [START_REF] Backović | Higher-order QCD predictions for dark matter production at the LHC in simplified models with s-channel mediators[END_REF] in which NLO K-factors have been investigated both at the total cross-section and differential distribution level for a series of representative benchmark scenarios. Here, we explore the impact of the NLO corrections on the exclusion limits originating from the t t + / E T channel. In order to examine the reach of the CMS search, we start by performing a twodimensional scan of the mediator and dark matter masses with fixed mediator couplings, similar to figure 7 in ref. [START_REF] Haisch | Simplified dark matter top-quark interactions at the LHC[END_REF]. The same scan is performed at both LO and NLO accuracy concerning the simulation of the hard scattering process, which allows us to determine the impact of the QCD corrections on the exclusion bounds. Before presenting the results for the excluded regions and to facilitate the discussion, we show the dependence of the LO cross section for g t = g X = 4 on the new physics masses and the corresponding K-factors in figure 7. The cross section is the largest in the low mass regions where the mediator can resonantly decay to a pair of dark matter particles, and falls steeply in the off-shell regions. In particular, the region where 2m X < m Y < 2m t is characterised by mediator decays either into a pair of dark matter particles or into a pair of gluons. These two decay rates are related by (see section 2)

Γ(Y 0 → gg) Γ(Y 0 → X X) = g 2 t g 2 X α 2 s 9π 2 β 3 X m 2 Y v 2 ∼ g t g X 2 × 10 -5 m Y GeV , (4.1)
which suggests that the decay rate into a pair of dark matter particles is always significantly higher, except in the case of a large hierarchy between the couplings (g t /g X 100). For m Y > 2m t , the Y 0 → t t decay mode is open, and t t + / E T production turns out to be suppressed by the visible decay channels of the mediator, unless g X > g t . Such a feature has already been illustrated in figure 1. The NLO K-factors related to t t + / E T production (right panel of figure 7) are found to vary from 0.96 to 1.15 in the range of masses examined here, the QCD corrections being more important in the low mass region.

The results for the exclusion regions are shown in figure 8 when LO (left panel) and NLO (right panel) simulations are used; see more details on the recasting procedure in appendix C.1. Setups excluded at the 40%, 68% and 95% confidence level (CL) are marked separately in the figures. As expected from the total cross section results, all excluded points (at the 95% CL) lie in the triangular low-mass region where the mediator resonantly decays into a dark matter particle pair. The exclusion region reaches mediator masses of about 200-250 GeV if close to threshold (m Y ∼ 2m X ). This region is in fact not exactly triangular as for a given mediator mass, not all dark matter masses below m Y /2 are excluded.

- [START_REF] Harris | Constraining dark sectors at colliders: beyond the effective theory approach[END_REF] This is related to the parametric choice of g t = g X = 4 for which the mediator width can become large. In this case, the narrow width approximation is not valid and the t t + / E T cross section acquires a dependence on the dark matter mass even in the resonant region.

Comparing the LO and NLO results, we observe that in the low mass resonant region where the K-factor is small and of about 1.10, the exclusion contours are mildly modified and this small 10% shift in the cross section does not lead to any significant change. For larger mediator masses, the K-factors are ∼ 1 and therefore do not imply a modification of the exclusion regions, if the central prediction at the default choice of scale is considered. However, the inclusion of NLO corrections significantly reduces the theoretical error and thus leads to sharper exclusion bounds as discussed below.

In order to further investigate the effects of the NLO corrections, we select three benchmark scenarios for which we perform a detailed study. These benchmarks are defined in table 5 where they are presented along with the corresponding LO and NLO cross-sections and the CL exclusion obtained with MadAnalysis 5. As discussed in appendix C.1, the most relevant observables for this analysis consist of the / E T , M T ( , / E T ) and M W T 2 for which distributions are shown in figure 9. We normalise the distributions to 100, 10 and 1 for the scenarios I, II and III respectively to ensure that they are all clearly visible in the figure. Moreover, we also indicate the scale uncertainty bands that have been obtained from a scale variation of 0.5µ 0 < µ R,F < 2µ 0 . In agreement with the findings of ref. [START_REF] Backović | Higher-order QCD predictions for dark matter production at the LHC in simplified models with s-channel mediators[END_REF], higher-order corrections have a rather mild effect on the distribution shapes for all key observables. Using NLO predictions however leads to a significant reduction of the scale uncertainties compared to the LO case. In table 5, one can also see that the use of NLO predictions leads to a significant reduction of the uncertainty in the cross section which propagates down to the CLs. NLO predictions therefore allow us to draw more reliable conclusions on whether a parameter point is excluded. 86.5 +8.6%

JHEP11(2016)111

-5.5%

Table 5. Benchmark scenarios used to investigate the impact of the NLO corrections on the t t + / E T CMS search. The LO and NLO cross sections at 8 TeV LHC are shown together with the CL exclusion obtained from MadAnalysis 5. The uncertainties originating from scale variation (0.5µ 0 < µ R,F < 2µ 0 ) are also shown.

Mono-X final states

In addition to the constraints that can be derived by means of t t+ / E T probes and that have been discussed in the previous section, mono-X searches can also be relevant for obtaining bounds on our top-philic dark matter model. Monojet [START_REF]Search for dark matter, extra dimensions and unparticles in monojet events in proton-proton collisions at √ s = 8 TeV[END_REF][START_REF]Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector[END_REF][START_REF]Search for dark matter with jets and missing transverse energy at 13 TeV[END_REF], mono-Z [START_REF]Search for dark matter and unparticles produced in association with a Z boson in proton-proton collisions at √ s = 8 TeV[END_REF][START_REF]Search for dark matter in events with a hadronically decaying W or Z boson and missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector[END_REF][START_REF]Search for dark matter in events with a Z boson and missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector[END_REF][START_REF]Search for new physics in the V /jet + MET final state[END_REF][START_REF]Search for dark matter produced in association with a hadronically decaying vector boson in pp collisions at √ s = 13 TeV with the ATLAS detector at the LHC[END_REF] and mono-Higgs [START_REF]Search for dark matter produced in association with a Higgs boson decaying to two bottom quarks in pp collisions at √ s = 8 TeV with the ATLAS detector[END_REF][START_REF]Search for dark matter in events with missing transverse momentum and a Higgs boson decaying to two photons in pp collisions at √ s = 8 TeV with the ATLAS detector[END_REF][START_REF]Search for new phenomena in events with missing transverse momentum and a Higgs boson decaying to two photons in pp collisions at √ s = 13 TeV with the ATLAS detector[END_REF][START_REF]Search for dark matter in association with a Higgs boson decaying to b-quarks in pp collisions at √ s = 13 TeV with the ATLAS detector[END_REF] signals have been searched for during the first run of the LHC, and these search results could be recast to constrain the dark matter model studied in this work. In contrast to tree-level dark matter production in association with a pair of top quarks, the production of a pair of dark matter particles with a jet, a Z-boson or a Higgs boson proceeds via a gluon fusion top-quark loop diagram. Although they have been largely studied by ATLAS and CMS, monophoton analyses cannot be used as charge conjugation invariance forbids the existence of a monophoton signal for the spin-0 mediator scenario.

Monojet. We start by discussing constraints that can be imposed by the CMS 8 TeV monojet analysis [START_REF]Search for dark matter, extra dimensions and unparticles in monojet events in proton-proton collisions at √ s = 8 TeV[END_REF]. For this study, hard-scattering events are generated at the LO accuracy within MG5 aMC, and the matching with parton showers is made with Pythia 6. The results are analysed in MadAnalysis 5 that also takes care of the detector simulation using its interface with Delphes 3. This recasting procedure allows us to exclude any specific parameter space point at any desired confidence level, our exclusion being conservatively derived on the basis of the signal region that drives the strongest bound. This limitation is related to the lack of public information, the statistical model used by CMS for the combination being not available. One can find more details for the recasting procedure in appendix C.2.

Similar to the t t + / E T analysis of the previous section, we perform a two-dimensional scan on the mediator and dark matter masses while fixing both new physics couplings to g t = g X = 4 (as in figure 5 in ref. [START_REF] Haisch | Simplified dark matter top-quark interactions at the LHC[END_REF]). Figure 10 shows our results, where we represent the scenarios excluded at the 40%, 68% and 95% CL. The bulk of the excluded points lie again in the triangular low-mass region where the mediator resonantly decays into a pair of dark matter particles. Except for the small subset of points excluded at the 40% and 68% CL in the region where m Y < 2m X , the extent of the exclusion region is determined by the significant reduction of the monojet cross section below the resonant production threshold already presented in figure 6. The pp → Y 0 j cross section indeed rapidly falls with m Y , reaching levels beyond the sensitivity of the 8 TeV search at m Y ∼ 500 GeV. In addition to the decrease of the Y 0 j production cross section, the opening of the mediator decay mode into a top-antitop system when m Y > 2m t leads to a further reduction of the monojet production rate. In comparison with the t t + / E T case, the monojet search overall appears to be more constraining, especially for higher mediator mass values thanks to the larger monojet cross section.

As shown in ref. [START_REF] Mattelaer | Dark matter production through loop-induced processes at the LHC: the s-channel mediator case[END_REF], the shape of key monojet differential distributions differs in the resonant and in the off-shell parameter space regions. While the total cross section falls dramatically in the off-shell region m Y < 2m X (as shown in figure 6), the / E T and jet transverse momentum distributions tend to be harder for off-shell production. We demonstrate this feature with a detailed investigation of three benchmark points defined in table 6. They consist of two resonant scenarios with different mediator masses and one non-resonant scenario. The monojet production rate is also indicated in the table, and we present normalised distributions relevant for the monojet analysis in figure 11. The off-shell scenario yields harder distributions compared to the resonant cases. This implies that a larger fraction of events features high missing transverse energy ( / E T >250 GeV) and populates the different signal regions of the CMS analysis. As a result, a better sensitivity is found than what one might expect from considering the total cross section alone. This feature leads to the exclusion of dark matter scenarios where m Y < 2m X , as depicted in figure 10. In our simulation of the monojet signal, we have ignored the possible impact of the merging of event samples featuring different final state jet multiplicities. A reliable description of the high transverse momentum spectra of the leading jet typically necessitates the merging of event samples including at least one and two jets in the final state [START_REF] Mattelaer | Dark matter production through loop-induced processes at the LHC: the s-channel mediator case[END_REF]. We have explicitly verified that for both resonant and off-shell scenarios, employing a merged sample does not have a big impact on the / E T distribution and therefore on the resulting exclusion contours. This originates from the analysis selection strategy that requires one single hard jet and rather loose requirements on the second jet, so that the configuration that dominates consists of a single hard jet recoiling against the missing energy. Such a configuration is described similarly by the one-jet and merged samples. We nevertheless stress that the importance of the merging procedure has to be checked on a case-by-case basis as this depends on the analysis, so that higher multiplicity samples might be necessary to accurately describe the relevant distributions.

JHEP11(2016)111

-25 -

JHEP11(2016)111

Mono-Z and mono-Higgs. In addition to the use of monojet processes, we explore the possibility of constraining the parameter space of our model using mono-Z and mono-Higgs production. While the production rates are much smaller than the monojet rate as seen in figure 6, the backgrounds can be also small. Therefore, these search channels can be sensitive to the top-philic simplified dark matter model, as we will see below. Here, instead of employing a full recasting procedure as in the t t + / E T and monojet analyses, we perform parton-level analyses to provide rough estimates of the constraints on our model parameters.

We rely on the CMS search for dark matter production in association with a Z-boson that decays leptonically [START_REF]Search for dark matter and unparticles produced in association with a Z boson in proton-proton collisions at √ s = 8 TeV[END_REF], in which a 95% CL upper limit on the visible cross section of 0.85 fb is obtained once a / E T requirement of at least 150 GeV and the minimal detector selection requirements for the leptons (p T > 20 GeV and |η | < 2.5) are considered. We generate events for this process, and after applying the above fiducial selection requirements we obtain a cross section of 0.30 fb for (m Y , m X ) = (100, 10) GeV and g t = g X = 1. We show in figure 12 the / E T and leading lepton transverse momentum distributions (red lines) without and with applying the selection strategy. While we have not performed a detailed study, simple estimates show good prospects for setting limits on the parameter space of the model using the mono-Z analysis results. Using the upper limit of 0.85 fb, scenarios with couplings close to g t ∼ 2 could be excluded in the resonant region (m Y > 2m X ) with m Y < 100 GeV. For larger mediator masses, the cross section starts to fall due to the reduction of the phase space. In the off-shell region (m Y < 2m X ), the mono-Z cross section suffers from the same drastic decrease seen in figure 6 for the t t + / E T and monojet cases. The same procedure can be repeated to constrain the parameter space of the model using mono-Higgs events on the basis of the results of the ATLAS search for dark matter production in association with a Higgs boson decaying into two bottom quarks [START_REF]Search for dark matter produced in association with a Higgs boson decaying to two bottom quarks in pp collisions at √ s = 8 TeV with the ATLAS detector[END_REF]. This search results in a 95% CL upper limit on the visible cross section of 3.6 fb for a / E T threshold of 150 GeV. In order to estimate a limit, we generate events for (m Y , m X ) = (100, 10) GeV and g t = g X = 1, and require the two b-quarks to have a transverse momentum p b 1 T > 100 GeV and p b 2 T > 25 GeV, a pseudorapidity |η b | < 2.5 and to be separated in the transverse plane by an angular distance ∆R(b 1 , b 2 ) < 1.5. Moreover, we only select events exhibiting at least 150 GeV of missing transverse energy. We show again in figure 12 the / E T and leading b-quark transverse momentum distributions (blue lines) without and with applying the above-mentioned selection requirements. We then include a b-tagging efficiency of 60% and extract an upper limit on the g t coupling by comparing our results to the ATLAS limit. Coupling values of g t > 2 are found to be excluded for m Y > 2m X with m Y < 100 GeV. All other parameter space regions suffer from the same limitations as the mono-Z case.

From our naive parton-level analysis, we have seen that mono-Z and mono-Higgs signals show promising signs of setting constraints on the parameter space of the model and therefore deserve dedicated studies, which will be reported elsewhere (see also ref. [START_REF] Goncalves | Boosting invisible searches via ZH: from the Higgs boson to dark matter simplified models[END_REF]). The sensitivity to such signals will benefit from applying more aggressive / E T thresholds to ensure the reduction of the corresponding backgrounds. As seen in figure 12, we obtain a rather hard / E T distribution [START_REF] Mattelaer | Dark matter production through loop-induced processes at the LHC: the s-channel mediator case[END_REF], especially for mono-Z production. The result implies that an increase in the / E T threshold requirement in future analyses could lead to a sig-

JHEP11(2016)111

Resonance search constraints on top-philic DM The different coloured areas are excluded by the diphoton [START_REF]Search for diphoton resonances in the mass range from 150 to 850 GeV in pp collisions at √ s = 8[END_REF] (orange), t t [START_REF]A search for t t resonances using lepton-plus-jets events in proton-proton collisions at √ s = 8 TeV with the ATLAS detector[END_REF] (magenta) and t tt t [START_REF]Search for Standard Model production of four top quarks in the lepton + jets channel in pp collisions at √ s = 8 TeV[END_REF] (blue) searches. We include information on the mediator width to mass ratios (green curves). We assume a negligible branching ratio to the invisible sector.

nificant improvement of the sensitivity, especially given the the fact that Standard Model backgrounds rapidly fall off with the increase in missing energy.

Constraints from searches without missing transverse energy

Dijet and diphoton resonances. Dijet and diphoton resonance search results could (in principle) be used to constrain the simplified top-philic dark matter model. Due to double-loop suppressions, mediator-induced contributions to dijet and diphoton production are only relevant in the parameter space regions where m Y < 2m X , 2m t (i.e. where the mediator cannot decay into top quarks and/or dark matter particles). The partial mediator decay rate into gluons is then always dominant (as mentioned in section 2) since

Γ(Y 0 → γγ) Γ(Y 0 → gg) ∼ 8 9 α 2 e α 2 s ≈ 10 -3 . (4.2)
All LHC dijet resonance searches focus on the dijet high invariant-mass region, leading to no useful constraints on the top-philic dark matter model. The lowest mediator mass that is probed is ∼ 500 GeV, with a visible cross section restricted to be smaller than 10 pb [START_REF]Search for narrow resonances in dijet final states at √ s = 8 TeV with the novel CMS technique of data scouting[END_REF].

Although the branching ratio of the mediator into a photon pair is very small, the background associated with a diphoton signal is low so that one expects to be able to obtain stringent constraints on the model from the diphoton search results. We focus here on the CMS 8 TeV diphoton search [START_REF]Search for diphoton resonances in the mass range from 150 to 850 GeV in pp collisions at √ s = 8[END_REF] that investigates resonance masses ranging from 150 GeV to 850 GeV and derives limits on the corresponding cross section. For instance, the 95% CL upper bound on the mediator-induced diphoton production cross section σ(pp → Y 0 → γγ)

JHEP11(2016)111

is of 20 fb (4 fb) for a mediator mass of 150 GeV (300 GeV). Making use of the pp → Y 0 cross section values shown in figure 6 and the Y 0 → γγ branching ratio computed from the formulas shown in section 2, we present diphoton constraints on the model in the (m Y , g t ) plane in figure 13. These results assume that the dark matter particle is much heavier than the mediator that can thus not resonantly decay invisibly. The constraints are found to be stringent below the 2m t threshold, where the g t coupling cannot be larger than 0.6.

Top-antitop resonances. For scenarios with mediator masses above the top-antitop threshold (m Y > 2m t ), t t resonance searches [START_REF]A search for t t resonances using lepton-plus-jets events in proton-proton collisions at √ s = 8 TeV with the ATLAS detector[END_REF][START_REF]Search for resonant t t production in proton-proton collisions at √ s = 8 TeV[END_REF] can be used as probes of the model. In our setup, loop-induced resonant mediator contributions can indeed enhance the t t signal, in particular when there is a large coupling hierarchy (g t g X ) or mass hierarchy (2m t < m Y < 2m X ). We derive constraints on our model from the ATLAS 8 TeV t t resonance search [START_REF]A search for t t resonances using lepton-plus-jets events in proton-proton collisions at √ s = 8 TeV with the ATLAS detector[END_REF] that relies on the reconstruction of the invariant mass of the topquark pair to derive a 95% CL exclusion on the existence of a new scalar particle coupling to top quarks. The associated cross section limits range from 3.0 pb for a mass of 400 GeV to 0.03 pb for m Y = 2.5 TeV, assuming that the narrow width approximation is valid with a mediator width being of at most 3% of its mass and that there is no interference between the new physics and Standard Model contributions to the t t signal. Constraints are computed using the NNLO mediator production cross section (see figure 6) and the relevant top-antitop mediator branching ratio derived from the formulas presented in section 2. The latter is in fact very close to one in the relevant region, the mediator decays into dark matter particle pairs being kinematically forbidden and those into gluons and photons loop-suppressed. The results are presented in the (m Y , g t ) plane in figure 13. This shows that scalar mediators with masses ranging from 400 GeV to 600 GeV could be excluded for g t couplings in the [START_REF] Bertone | Particle dark matter: evidence, candidates and constraints[END_REF][START_REF] Clowe | A direct empirical proof of the existence of dark matter[END_REF] range, the exact details depending on m Y and on the fact that the narrow-width approximation must be valid. This demonstrates the ability of the t t channel to probe a significant portion of the m Y > 2m t region of the model parameter space. In the region where 2m t , 2m X < m Y , the partial decay Y 0 → X X reduces the t t signal and therefore limits the sensitivity of the search.

Four-top signals. Scenarios featuring a mediator mass above twice the top-quark mass can be probed via a four-top signal, since the mediator can be produced in association with a pair of top quarks and further decay into a top-antitop system. Theoretically, the Standard Model four-top cross section has been calculated with high precision [START_REF] Bevilacqua | Constraining BSM physics at the LHC: four top final states with NLO accuracy in perturbative QCD[END_REF], but the sensitivity of the 8 TeV LHC run was too low to measure the cross section. Instead, an upper limit on the cross section at a centre-of-mass energy of 8 TeV has been derived [START_REF]Search for Standard Model production of four top quarks in the lepton + jets channel in pp collisions at √ s = 8 TeV[END_REF][START_REF]Search for new physics in events with same-sign dileptons and jets in pp collisions at √ s = 8 TeV[END_REF].

The four-top production rate is constrained to be below 32 fb [START_REF]Search for Standard Model production of four top quarks in the lepton + jets channel in pp collisions at √ s = 8 TeV[END_REF], a value that has to be compared to the Standard Model prediction of about 1.3 fb. Only models with new physics contributions well above the background (see e.g. ref. [START_REF] Beck | Probing top-philic sgluons with LHC run I data[END_REF]) can therefore be constrained by the four-top experimental results.

In our top-philic dark matter model, the new physics contributions to the four-top cross section can be approximated by the t tY 0 cross section, the branching ratio B(Y 0 → t t) ∼ 1. Using the NLO cross section (see figure 6), we derive limits that we represent in the (m Y , g t ) plane in figure 13. A small region of the parameter space with g t > 2.5 and in which the

JHEP11(2016)111

mediator mass lies in the [2m t , ∼ 450 GeV] mass window turns out to be excluded. The weakness of the limit is related to the steeply decreasing cross section for pp → Y 0 t t with the increase in m Y .

The mediator width. In all the above studies where the final state does not contain any missing energy, the mediator width has been assumed narrow. Concerning the diphoton channel, this assumption holds within the entire excluded region as only loop-suppressed gluon and photon mediator decays are allowed. In the region where m Y > 2m t , the width of the mediator rises quickly with its mass, and the width over mass ratio rapidly exceeds the 3% value that has been imposed in the ATLAS t t resonance search [START_REF]A search for t t resonances using lepton-plus-jets events in proton-proton collisions at √ s = 8 TeV with the ATLAS detector[END_REF] as can be seen in figure 13. The reinterpretation of the ATLAS results to a generic t t resonance model should therefore be made carefully, as the limit cannot be necessarily applied to scenarios featuring significantly larger mediator widths. This is shown in figure 13 by a dotted line, and we can also observe that most of the points that would have been excluded by the ATLAS search do not fulfil the requirement of a width below 3% of the mediator mass. In our excluded region of the parameter space, we allow the mediator width to reach 8% of its mass, by the virtue of the experimental resolution on the invariant mass of the t t system. This leads to the exclusion of scenarios with mediator masses up to 600 GeV.

The ATLAS resonance t t study claims that varying the width of the resonance from 10% to 40% for the massive gluon model results in a loss in sensitivity by a factor 2 for a 1 TeV resonance. An extension of the reinterpretation of the ATLAS limits on our simplified top-philic dark matter model to the case of larger resonance widths could then be performed by rescaling the limits by the appropriate correction factor. We have nonetheless found that no additional points are excluded even without rescaling the sensitivity of the search as the ATLAS analysis rapidly loses sensitivity for resonance masses above 600 GeV. Considering model points with a mediator width to mass ratio of at most about 8% therefore provides a realistic exclusion over the entire model parameter space.

Concluding remarks on direct mediator searches. Mediator resonance searches at 8 TeV show good prospects of constraining our simplified top-philic dark matter model, especially in the mediator mass range of 150-345 GeV and 400-600 GeV by means of the diphoton and top-pair searches respectively. So far, the t t resonance searches are strictly applicable to a limited parameter space region of the simplified model, and considering larger widths in the interpretation of the future results would allow for a more straightforward reinterpretation of the limits to a wider range of parameters. Concerning the four-top analysis, it can presently only exclude a restricted part of the parameter space, but future measurements are expected to lead to more competitive bounds.

Finally, the pp → tt + j channel could also be used to probe dark matter models coupling preferably to top quarks. This has been for instance shown in ref. [START_REF] Greiner | Model-independent production of a top-philic resonance at the LHC[END_REF] where a loop-induced production of t tj can in some cases lead to interesting constraints on topphilic models of new physics. In our case, they are nonetheless not expected to give more stringent constraints than the t t resonance searches. One could also consider the pp → t ttj and pp → t tW t processes [START_REF] Greiner | Model-independent production of a top-philic resonance at the LHC[END_REF]. Because of the magnitude of the electroweak couplings, these processes are characterised by smaller cross sections than when four top quarks are -29 -

JHEP11(2016)111

LHC constraints on top-philic dark matter involved, and are hence not likely to set more stringent constraints on the class of models under consideration.

Combined constraints

The final segment of our comprehensive study of top-philic dark matter simplified models is a combined study of astrophysical and collider constraints. We find that in the region where g X , g t ≤ π, the 8 TeV collider results that provide relevant bounds (once the relic density and direct detection constraints are imposed) originate from direct mediator production searches when the mediator further decays into a pair of Standard Model particles. Figure 14 illustrates our results and shows the scenarios that are excluded by resonant diphoton and top-pair searches as well as by the four-top analysis. All points in the plot accommodate the dark matter relic density and direct detection constraints, while the colours indicate points excluded by individual complementary collider bounds. The vast majority of excluded points lie in the region where 2m X > m Y with m Y ∈ [150, 600] GeV. This is the region where the mediator decay into a pair of dark matter particles is kinematically forbidden, ensuring large branching fractions for decays into Standard Model particles. The diphoton resonance search excludes points below the 2m t threshold, while t t results constrain the 400 < m Y < 600 GeV region. The four-top probe is able to exclude a narrow parameter space region close to m Y ∼ 2m t , in agreement with the findings shown in figure 13.

Relaxing the requirements on the relic density, the direct detection and the upper bound on the coupling strengths allows for another meaningful study of combined collider constraints. For this purpose we have performed a joint analysis of collider bounds on the top-philic simplified dark matter model in the scope of a four-dimensional parameter scan -30 -JHEP11(2016)111 In case of resonant t t, four top, and the combined constraints, we only show the 95% CL exclusion as the t t and four top results have not been obtained using a recast LHC analysis.

JHEP11(2016)111

with a flat likelihood function over all dimensions. We have performed the scan by restricting the couplings to be smaller than 2π, as well as by allowing the mediator widths to reach 50% of the mediator mass. Figure 15 shows our results, where the upper left panel shows the model points excluded by the combination of all collider results, and the rest of the panels show the points excluded by individual LHC Run I collider results. We find that the 8 TeV monojet searches exclude model points which lie mainly in and around the triangle bounded by the m Y = 2m X and m Y = 2m t lines, where the characteristic g t which is excluded by the 8 TeV results is of O [START_REF] Bell | Mono-W dark matter signals at the LHC: simplified model analysis[END_REF]. The region in which the excluded points are located is reasonable, as we expect any significant monojet signal in the region where m Y > 2m X . Furthermore, we expect the branching ratio to missing energy to be lower in the region where m Y > 2m t due to the kinematically allowed decays into a pair of top quarks. This in turn leads to a lower signal cross section in all channels with missing energy and hence a lower number of points which can be excluded by monojet searches in the m Y > 2m t region.

The points excluded by the 8 TeV t t + / E T measurements lie in roughly the same region as the points excluded by the monojet search, but with a more defined edge of m Y = 2m t . Conversely, the 8 TeV t t resonance search provides constraints in the region of m Y ∈ [400, 600] GeV and m X 100 GeV, and is able to rule out g t couplings of O(1). The four top searches constrain roughly the same region of the (m Y , m X ) parameter space as the t t searches. However, the characteristic size of the couplings four top searches are able to constrain is significantly larger than the case of t t.

Finally the diphoton resonant search excludes m Y ∈ [150, 2m t ] GeV with 2m X > m Y , ruling out g t couplings larger than 0.6. In the (m Y , m X ) plane, we can observe that the constraints arising from all mediator resonance searches, i.e. the diphoton and t t analyses, are largely complementary to those issued from searches in channels with large missing energy.

Conclusions

We presented a comprehensive analysis of simplified top-philic dark matter models, in the scope of collider physics, astrophysics and cosmology. Our study considered the full four dimensional model parameter space, where we treated the experimental constraints on the model space both separately and in conjunction with each other. The requirement of predicting the measured relic density Ω DM h 2 gives the most stringent constraint on the viable regions of the parameter space. Most of the region where m Y > m X cannot accommodate the observed relic density, except near the resonance m Y ∼ 2m X and for m X > m t . Direct detection data complementary excludes large portions of the parameter space in the m Y < m X region once experimental results from LUX and CDMSLite are accounted for. In the context of dark matter indirect detection, we studied prospects for further model constraints from gamma-ray flux measurements originating from dwarf spheroidal galaxies and the gamma-ray lines issued from the inner galactic region. In the specific model we consider, the dark matter annihilation cross section is p-wave suppressed, leading to indirect detection bounds which are too weak to provide additional constraints on the parameter space.

Collider searches from LHC Run 1 at √ s = 8 TeV can constrain the parameter space beyond the limits obtained from the relic density and direct detection, but apply mostly JHEP11(2016)111

in the limit of coupling values 1. We found that for couplings of π, the resonant t t and diphoton searches are able to exclude a fraction of model points in the regions of m Y ∼ 400 -600 GeV and m Y ∼ 150 -350 GeV respectively, even upon assuming astrophysical and relic density constraints.

In addition to studying collider signatures of the top-philic dark matter simplified model as a complementary way of dark matter detection, we performed a study of collider constraints without assuming relic density and direct detection (as well as extended the parameter range to include coupling values of < 2π and Γ Y ≤ 0.5m Y ). Our results for a four dimensional parameter scan show that (in the scenario where astrophysical and cosmological constraints are not relevant), / E T + j and / E T + t t 8 TeV results provide meaningful bounds on the model parameter space in the 2m X < m Y < 2m t region, but only for g t , g X π. In the m X > m t region, the resonant t t searches are again able to exclude some model points in the m Y ∼ 400 -500 GeV region, while γγ measurements provide constraints in the m Y < 2m t region. We have also explored the prospects of using rarer processes such as four-top production as well as mono-Z and mono-Higgs production to constrain our model. While we have not performed a detailed analysis we have found that these processes show promising signs of further constraining the parameter space of our model and deserve dedicated studies.

For the purposes of our study we have recast the CMS monojet and / E T + t t searches in the framework of MadAnalysis 5, which allows us to reliably extract constraints on our model, and can benefit future collider studies which go beyond our simplified model and even beyond dark matter searches. Another important aspect of our work, is the use of NLO QCD predictions for the / E T + t t process to constrain our model. While we find that K-factors for this process are close to one, the importance of taking higher order effects into account lies in the reduced theoretical uncertainties of the NLO results. We have shown that the uncertainties in the CL estimates significantly reduce with the inclusion of higher order QCD terms which clearly illustrates the importance of higher order corrections on the interpretation of dark matter searches at colliders.

The work presented in this paper also represents a proof-of-concept for unified numerical framework for dark matter studies at the interface of collider physics, astrophysics and cosmology in a generic model. As a part of consistency checks, we have ensured that the scan covers similar regions of the parameter space both in case of MadDM and micrOMEGAs. Figure 18 shows the results for distributions of masses and couplings in the scans, where the blue/red lines refer to MadDM/micrOMEGAs respectively. Similarities in the distributions of figure 18 indicate that parameter scanning was performed consistently between the two codes. We also made an explicit comparison between the projections of the four dimensional parameter scans obtained with MadDM and micrOMEGAs respectively. Figures 19 and20 show several examples. We don't find significant deviations between the results obtained in the two codes except in the region of m Y ∼ 2m X and m X ∼ O(1) GeV. In these regions we expect some discrepancies due to the possible numerical instabilities in integration of the thermally averaged cross section for amplitudes which feature resonances of extremely small widths.

JHEP11(2016)111

JHEP11(2016)111

B.2 Details on the dark matter annihilation cross sections

In this appendix we give the detailed analytic expression of the three annihilation processes described in section 3.1.

The s-channel annihilation cross section X X → t t (process (I)) is given by:

σ(X X → t t) = 3g 2 X g 2 t y 2 t 32πs (s -4m 2 t ) 3/2 s -4m 2 X (m 2 Y -s) 2 + m 2 Y Γ 2 Y . (B.1)
Process (II) denotes the annihilation of dark matter into a pair of gluons via s-channel and is given by:

σ(X X → gg) = g 2 g g 2 X 16πv 2 s 3/2 s -4m 2 X (m 2 Y -s) 2 + m 2 Y Γ 2 Y (B.2)
Finally the process (III), namely X X → Y 0 Y 0 via t-channel is given by:

σ(X X → Y 0 Y 0 ) = g 4 X 64π h(t 0 ) -h(t 1 ) s(s -4m 2 X ) , (B.3)
where t 0,1 are the integration extrema: and the undefined integral h(t) has the form:

t 0,1 = - 1 4 s -4m 2 X ∓ s -m 2 Y 2 , (B.4) -37 - JHEP11 (
h(t) ≡ (m 2 Y -4m 2 X ) 2 m 2 X -u - (m 2 Y -4m 2 X ) 2 m 2 X -t -4t + 6m 2 Y -4m 2 Y (4m 2 X + s) -32m 4 X + 16m 2 X s + s 2 2m 2 Y -s log t -m 2 X m 2 X -u , (B.5)
with t and u Mandelstam variables such that u = 2m 2 X + 2m 2 Y -st. In general the thermally averaged cross section can be approximated in the non relativistic regime by expanding the cross section in powers of the dark matter relative velocity v rel , with s m 2 X (4+v 2 rel ), weighting with the appropriate Maxwell-Boltzmann distribution:

σv rel (x) i = A i + 3 2 B i x + O(x -2 ) (B.6)
where the index i indicated the annihilation process, x ≡ m X /T and T is the temperature of the dark matter gas. In case of s-channel annihilation, along the resonance the thermal average is much more complex and requires the full computation of the integral dx σv rel (x). The approximation given in eq. (B.6) holds in all regions far away from the resonance and is useful to show the dependence on v rel of σv rel (x) for each specific process.

-38 -JHEP11(2016)111 For all processes (I), (II) and (III) the first coefficient is always null, A i = 0. The first non negligible term in the expansion eq. (B.6) is then B:

B t t = 3g 2 X g 2 t y 2 t 16π m 2 X (1 -m 2 t /m 2 X ) 3/2 (m 2 Y -4m 2 X ) 2 , (B.7) B gg = 2g 2 g g 2 X m 4 X 2πv 2 m 2 Y -4m 2 X 2 , (B.8) B Y 0 Y 0 = g 4 X 24π m 2 X (9m 4 X -8m 2 X m 2 Y + 2m 4 Y ) (2m 2 X -m 2 Y ) 4 1 - m 2 Y m 2 X . (B.9)
This is equivalent to say that all three process are p-wave suppressed for dark matter annihilation at present epoch. The case of Dirac dark matter particles communicating with the SM via a pseudoscalar mediator has been described in [START_REF] Arina | Dark matter with pseudoscalar-mediated interactions explains the DAMA signal and the galactic center excess[END_REF], where analytic expressions for σv rel can be found. Similarly to scalar mediator Y 0 the t-channel process is again p-wave suppressed, while the s-channel annihilation is dominated by s-wave.

C Recasting of LHC searches within the MadAnalysis 5 framework

In this appendix, we detail the implementation, within the MadAnalysis 5 framework [START_REF] Conte | MadAnalysis 5, a user-friendly framework for collider phenomenology[END_REF][START_REF] Conte | Designing and recasting LHC analyses with MadAnalysis 5[END_REF][START_REF] Dumont | Toward a public analysis database for LHC new physics searches using MadAnalysis 5[END_REF], of the two dark matter searches that we have investigated in this work. More precisely,

JHEP11(2016)111

this consists of the CMS-B2G-14-004 analysis [START_REF]Search for the production of dark matter in association with top-quark pairs in the single-lepton final state in proton-proton collisions at √ s = 8 TeV[END_REF] that probes final states comprised of a top-antitop system produced in association with a pair of invisible dark matter particles (see section C.1) and the CMS-EXO-12-048 analysis [START_REF]Search for dark matter, extra dimensions and unparticles in monojet events in proton-proton collisions at √ s = 8 TeV[END_REF] related to the production of a pair of dark matter particles together with a hard jet (see section C.2). Both recasting codes have been validated within the version 1.3 of MadAnalysis 5, although the monojet search reimplementation is also compatible with the version 1.2 of the program. The simulation of the detector response is performed with the standard Delphes 3 package that we have run from the MadAnalysis 5 platform. In the monojet case, we have used the standard CMS detector parameterisation that is the shipped with MadAnalysis 5, while in the top-antitop plus missing energy case, we have designed a dedicated detector card. For both setups, jets are reconstructed on the basis of the anti-k T algorithm [START_REF] Cacciari | The anti-k t jet clustering algorithm[END_REF] with a radius parameter set to 0.5, as implemented in FastJet [START_REF] Cacciari | FastJet user manual[END_REF].

The validation of both our reimplementations is based on material provided by CMS. Two UFO models [START_REF] Degrande | UFO -the Universal FeynRules Output[END_REF], one for each of the recast analyses, have been shared so that we have been allowed to generate specific dark matter signals for which CMS has released public cutflow charts and differential distributions. Using MG5 aMC [START_REF] Alwall | The automated computation of tree-level and next-to-leading order differential cross sections and their matching to parton shower simulations[END_REF] (with the leading order set of CTEQ6 parton densities [START_REF] Pumplin | New generation of parton distributions with uncertainties from global QCD analysis[END_REF]) and Pythia 6 [START_REF] Sjöstrand | PYTHIA 6.4 physics and manual[END_REF] (with the Z * 2 tune [START_REF] Field | Min-bias and the underlying event at the LHC[END_REF] for the description of the underlying events) for the simulation of the hard scattering process and of the parton showering and hadronisation, respectively, we have generated signal events that have been analyzed with MadAnalysis 5. Our results have been confronted to the CMS official numbers, which has allowed us to assess the validity of our recasting codes. Our simulation procedure moreover includes the generation of matrix elements containing up to two extra jets that we have merged according to the MLM prescription [START_REF] Mangano | Matching matrix elements and shower evolution for top-quark production in hadronic collisions[END_REF][START_REF] Alwall | QCD radiation in the production of heavy colored particles at the LHC[END_REF], the merging scale being set to 40 GeV.

All Pythia 6, Delphes 3 and MG5 aMC configuration cards can be downloaded from the public analysis database webpage of MadAnalysis 5, http://madanalysis. irmp.ucl.ac.be/wiki/PublicAnalysisDatabase, while the recasting C++ codes associated with the CMS-EXO-12-048 and CMS-B2G-14-004 analyses can be found on In-Spire [START_REF] Guo | MadAnalysis5 implementation of the CMS monojet search[END_REF][START_REF] Fuks | MadAnalysis5 implementation of the CMS search for dark matter production with top quark pairs in the single lepton channel[END_REF].

C.1 The CMS top-antitop plus missing energy CMS-B2G-14-004 search In order to validate our reimplementation of the CMS-B2G-14-004 search in MadAnalysis 5, we focus on a new physics model that features the production a pair of dark matter particle X of mass m X = 1 GeV in association with a top-antitop pair via a four-fermion interaction. The CMS event selection strategy requires a large amount of missing transverse energy, a single isolated lepton and multiple jets, and uses 19.7, fb -1 of proton-proton collision data recorded at a center-of-mass energy of √ s = 8 TeV.

The CMS-B2G-14-004 analysis relies on single electron and muon triggers, with lower p T thresholds of 27 GeV and 24 GeV respectively, and the reconstructed electron (muon) candidate is imposed to be isolated in such a way that the sum of the transverse momenta of all objects lying in a cone of radius R = 0.3 centered on the lepton has to be smaller than 10% (12%) of the lepton p T . Event preselection finally requires that the lepton p T is larger than 30 GeV and pseudorapidity |η| is smaller than 2.5 ( demands the presence of at least three jets of p T > 30 GeV and |η| < 2.4 with one of them being b-tagged, as well as missing energy / E T > 160 GeV. The signal region is defined by selecting events with a large amount of missing transverse energy / E T > 320 for which the transverse mass M T that is constructed from the lepton and the missing energy is larger than 160 GeV. Moreover, the missing transverse momentum and the two leading jets are asked to be well separated in azimuth, ∆Φ j 1,2 , / E T > 1.2, and the M W T 2 variable [107] is enforced to be greater than 200 GeV.

JHEP11(2016)111

In table 7, we confront the cutflow chart that has been obtained with MadAnalysis 5 to the official results of CMS for the benchmark scenario under consideration. For each step of the selection, we have calculated the related efficiency defined as

i = n i n i-1 , (C.1)
where n i and n i-1 mean the event number after and before the considered cut, respectively. The relative difference information given in the table corresponds to the difference between the MadAnalysis 5 and the CMS efficiencies, normalized to the CMS result,

δ rel i = 1 - MA5 i CMS i . (C.2)
An agreement at the percent level has been found all over the selection procedure. Moreover, we compare several (normalized) differential distributions as calculated with Mad-Analysis 5 when all selection steps but the one related to the represented kinematic variable are included with the public CMS results in figure 21. A very good agreement can again be observed.

C.2 The CMS monojet CMS-EXO-12-048 search

The validation of our implementation of the CMS-EXO-12-048 search in MadAnalysis 5 has been achieved on the basis of a benchmark scenario that is inspired by refs. [START_REF] Beltrán | Maverick dark matter at colliders[END_REF][START_REF] Goodman | Constraints on dark matter from colliders[END_REF][START_REF] Goodman | Constraints on light Majorana dark matter from colliders[END_REF][START_REF] Bai | The Tevatron at the frontier of dark matter direct detection[END_REF]. In this context, monojet events arise from the associated production of a pair of invisible Dirac fermions of mass of 1 GeV with at least one hard jet. The interactions of the dark particle with the Standard Model are mediated by a new gauge boson Z of mass and width of 40 TeV and 10 GeV respectively, and all new physics interactions have been assumed to have a vector coupling structure and a strength set equal to 1. Concerning our signal simulation setup, we have imposed that all parton-level jets have a transverse momentum p T larger than 20 GeV and that the leading jet has a p T > 80 GeV.

The CMS monojet search relies on an integrated luminosity of 19.7 fb -1 of protonproton collisions at a center-of-mass energy of √ s = 8 TeV. It focuses on a signal containing a very hard jet with a transverse momentum satisfying p T > 110 GeV and a pseudorapidity smaller than 4.5 in absolute value. A second jet is moreover allowed, provided that its transverse momentum is larger than 30 GeV, its pseudorapidity satisfies |η| < 4.5 and if it is well separated from the first jet by 2.5 radians in azimuth. Events featuring more than two jets (with p T > 30 GeV and |η| < 4.5), isolated electrons or muons with a transverse momentum p T > 10 GeV or hadronically decaying tau leptons with a transverse momentum p T > 20 GeV and a pseudorapidity satisfying |η| < 2.3 are discarded. The analysis then contains seven inclusive signal regions in which the missing energy / E T is required to be above specific thresholds of 250, 300, 350, 400, 450, 500 and 550 GeV respectively.

The selection strategy of the CMS monojet analysis thus consists of six preselection cuts followed by one region-dependent cut, when we ignore the first two requirements of the analysis related to the cleaning of the events from the detector noise that cannot be handled with Delphes 3. For the benchmark scenario under consideration, we compare -42 - the results that have been derived with our MadAnalysis 5 reimplementation with those provided by the CMS collaboration in table 8.

JHEP11(2016)111

JHEP11(2016)111

We have found that all selection steps are properly described by our implementation, with the exception the missing energy selection / E T > 250 GeV for which a disagreement of about 20% has been observed. It is however not uncommon that low missing energy is difficult to simulate with a fast-simulation of the detector based on Delphes 3. We have verified that for missing energy values of interest, the description of the missing energy agree relatively well with CMS, as illustrated in figure 22 where we compare, for a benchmark scenario where the Z mass has been set to 900 GeV, the missing energy distribution as obtained by CMS to the one derived with MadAnalysis 5.

Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

Figure 1 .

 1 Figure1. Ratio of the mediator width to its mass Γ Y /m Y (upper panels) and mediator branching ratios (lower panels) as a function of the mediator mass for different coupling choices and a dark matter mass fixed to m X = 50 GeV (solid lines) and 300 GeV (dashed lines).

Figure 2 .

 2 Figure 2. Results of our four-dimensional parameter scan using MadDM projected onto the (m Y , m X ) plane. The first three panels show the projections with the colourmap representing the values of Γ Y , g X and g t respectively. The right-most panel shows the zoomed-in upper right region of the left-most panel. All represented points feature a relic density in agreement with Planck data, Γ Y /m Y ≤ 0.2 and Γ Y > 10 -11 GeV (cf. table3).

Figure 3 .

 3 Figure 3. Results of our four-dimensional parameter scan using MadDM, projected onto the (g X , g t ) plane (left) and (m X , g t ) plane (right). All represented points feature a relic density in agreement with Planck data, Γ Y /m Y ≤ 0.2 and Γ Y > 10 -11 GeV (cf. table3).

Figure 4 .

 4 Figure 4. Results of our four-dimensional parameter scan using MadDM. The top left panel shows the projection of the scan into the (m Y , m X ) plane with a colourmap representing the values of Γ Y . The top right panel shows the projection of the allowed points into the (g X , g t ) plane with a colourmap given by m X . Finally the lower panel shows a projection onto the (m X , g t ) plane with a colourmap denoting the values of m Y . All represented points feature a relic density in agreement with Planck data, Γ Y /m Y ≤ 0.2, Γ Y > 10 -11 GeV and accommodate the direct detection constraints (cf. table3).

Figure 5

 5 Figure 5. Dark matter annihilation cross section at present time that is relevant for gamma-ray limits extracted from dwarf spheroidal galaxies measurements (left) and gamma-ray line searches (right). We show a maximal estimate of (σv rel ) tot and (σv rel ) γγ obtained by choosing v rel = 2v ∞ , where v ∞ is the escape velocity for dwarf spheroidal galaxies and the galactic center, respectively. All represented points feature a relic density in agreement with Planck data, Γ Y /m Y ≤ 0.2, Γ Y > 10 -11 GeV and accommodate the direct detection constraints (cf. table3).

1 -

 1 10 3 pb are expected for the production of light mediators with m Y 100 GeV at a centre-of-mass energy of 8 TeV (left panel), the dominant mechanism being the loop-induced gg → Y 0

Figure 6 .

 6 Figure 6. Left: total cross sections (with scale uncertainties) for various mediator production channels (with g t = 1) at a centre-of-mass energy of √ s = 8 TeV as a function of the mediator mass. The NNLO cross section for single mediator production σ(Y 0 ) is taken from the Higgs Cross section Working Group report, the Y 0 t t one is computed at NLO accuracy and all other loopinduced processes are evaluated at LO accuracy. The monojet (Y 0 j), mono-Z (Y 0 Z) and mono-Higgs (Y 0 h) cross sections include a transverse momentum cut on the mediator as indicated in the figure.In the lower panel, we show the ratios of the cross sections evaluated at a centre-of-mass energy of √ s = 13 TeV over those at 8 TeV. Right: cross sections for t t + / E T and monojet (with / E T > 150 GeV) production for a mediator mass of m Y = 50 and 100 GeV and at a centre-of-mass energy of √ s = 8 TeV given as a function of the dark matter mass.

4. 1

 1 Constraints from searches with missing transverse energy 4.1.1 The t t + / E T final state

1 Figure 7 .

 17 Figure 7. LO cross sections (left) and corresponding K-factors (right) for pp → t tX X at √ s = 8 TeV as a function of the mediator and dark matter masses. The top and dark matter couplings to the mediator are set to 4.

Figure 9 .

 9 Figure 9. Differential distributions for / E T , M T ( , / E T ) and M W T 2 for the three scenarios of table 5 at LO and NLO. The distributions are normalised to 100, 10 and 1 for Scenarios I, II and III respectively, and the scale uncertainty bands obtained by varying the renormalisation and factorisation scale in the range of 0.5µ 0 < µ R,F < 2µ 0 are also shown.

Figure 11 .Figure 12 .

 1112 Figure 11. Differential distributions for the missing transverse energy / E T and the hardest jet transverse momentum p T (j 1 ) for the three scenarios defined in table6. The distributions are normalised to one.

Figure 13 .

 13 Figure 13. Resonance search constraints from the LHC results at a collision centre-of-mass energy of 8 TeV on the simplified top-philic dark matter model presented in terms of the mediator mass m Y and the g t coupling.The different coloured areas are excluded by the diphoton[START_REF]Search for diphoton resonances in the mass range from 150 to 850 GeV in pp collisions at √ s = 8[END_REF] (orange), t t[START_REF]A search for t t resonances using lepton-plus-jets events in proton-proton collisions at √ s = 8 TeV with the ATLAS detector[END_REF] (magenta) and t tt t[START_REF]Search for Standard Model production of four top quarks in the lepton + jets channel in pp collisions at √ s = 8 TeV[END_REF] (blue) searches. We include information on the mediator width to mass ratios (green curves). We assume a negligible branching ratio to the invisible sector.

Figure 14 .

 14 Figure 14. Results of our four-dimensional parameter scan projected onto the (m Y , m X ) plane once constraints set from the LHC results are imposed. All represented points satisfy the relic density constraint, Γ Y /m Y ≤ 0.2, Γ Y > 10 -11 GeV and the direct detection constraints (cf. table3).

Figure 15 .

 15 Figure 15. Constraints derived from the LHC Run I results on the simplified top-philic dark matter model. The panels show results of a four-dimensional parameter scan, uniform on the linear scale. The upper left panel shows a combination of all relevant collider constraints. The upper right panel shows the points excluded by monojet constraints, while the third panel shows the points excluded by t t + /E T constraints. The resonant t t searches constraints are shown in the fourth panel, while the last two panels show the four top and the diphoton constraints. The results assume couplings smaller than 2π and Γ Y /m Y < 0.5, with no constraints from astrophysics or cosmology being imposed. In case of resonant t t, four top, and the combined constraints, we only show the 95% CL exclusion as the t t and four top results have not been obtained using a recast LHC analysis.

Figure 16 .Figure 17 .

 1617 Figure 16. Γ Y /m Y in the (g t , g X ) plane for different mass choices (expressed in GeV). The colour bar shows the numerical value of the width to mass ratio.

Figure 18 .

 18 Figure 18. Consistency check for the parameter scan. The panels show the distribution of couplings, masses and Y 0 widths and relic densities resulting from a MultiNest parameter scan. The vertical lines in the panel showing the distribution of relic density represent the 5σ Planck bound.

Figure 19 .

 19 Figure 19. Four dimensional parameter scan using micrOMEGAs, projected onto the m X , m Y plane. The first three panels show the projections with the colourmap representing the values of Γ Y , g X and g t respectively. The right-most panel shows the zoomed-in upper right region of the left-most panel. The scan is made under the same assumptions as figure 2.

Figure 20 .

 20 Figure 20. Consistency check for the parameter scan. The panels show projections of the four dimensional parameter scan using micrOMEGAs with the same assumptions as in figure 4.

2 RFigure 21 .

 221 Figure 21. Missig transverse energy, M T , M W T 2 and ∆Φ j 1,2 , / E T spectrum as obtained with MadAnalysis 5 (blue) once all selection steps but the one related to the represented variable are applied, compard to the CMS official results (red).

Figure 22 .

 22 Figure 22. Missing energy spectrum as obtained with MadAnalysis 5 (green dashed line) after the CMS-EXO-12-048 monojet preselection, compared to the CMS official results (red solid line). The last bin is the overflow bin.

Table 1 .

 1 Signatures of our simplified top-philic dark matter model.

Table 2 .

 2 MultiNest parameters and prior ranges for the four free parameters. All priors are uniform over the indicated range. Comment

	Measurement	ΩDMh 2	0.1198 ± 0.0015	Planck 2015 [44]
	Limits	ΓY /mY	< 0.2	Narrow width approximation
		ΓY	> 10 -11 GeV	Ensures prompt decay at colliders
		σ SI n	< σ SI LUX (90% CL)	LUX bound [45] (mX > 8 GeV)
		σ SI

n < σ SI CDMS (95% CL) CDMSlite bound

[START_REF] Agnese | New results from the search for low-mass weakly interacting massive particles with the CDMS low ionization threshold experiment[END_REF]

(1 GeV < mX < 8 GeV)

Table 3 .

 3 Summary of the observables and constraints used in this analysis and encoded into our MultiNest routine. The relic density constraints assume a Gaussian likelihood function, while the direct detection limits use step likelihood functions smoothed with half a Gaussian.

Table 4 .

 4 Summary of the 8 TeV LHC constraints used in this paper.

	Final state Imposed constraint	Reference	Comments
	/ E T + t t	See appendix C.1	CMS [68]	Semileptonic top-antitop decay
	/ E T + j	See appendix C.2	CMS [69]
	/ E T + Z	σ( / E T > 150 GeV) < 0.85 fb CMS [70]	Leptonic Z-boson decay
	/ E T + h	σ( / E T > 150 GeV) < 3.6 fb	ATLAS [71] h → b b decay
	jj	σ(m Y = 500 GeV) < 10 pb	CMS [72]	Only when m Y > 500 GeV
	γγ	σ(m Y = 150 GeV) < 30 fb	CMS [73]	Only when m Y > 150 GeV
	t t	σ(m Y = 400 GeV) < 3 pb	ATLAS [74] Only when m Y > 400 GeV
	t tt t	σ < 32 fb	CMS [75]	Upper limit on the SM cross section

  Figure 8. Constraints on simplified top-philic dark matter scenarios from the CMS 8TeV t t + / E T analysis[START_REF]Search for the production of dark matter in association with top-quark pairs in the single-lepton final state in proton-proton collisions at √ s = 8 TeV[END_REF]. The top and dark matter couplings to the mediator are set to 4 while the mediator and dark matter masses are allowed to vary freely. LO and NLO exclusions are respectively shown in the left and right panels of the figure.

			ttX X constraints (g X =g t =4; LO; 8 TeV LHC )			ttX X constraints (g X =g t =4; NLO; 8 TeV LHC )
		300	All points						300	All points
		250	Excluded points (40% CL) Excluded points (68% CL)				250	Excluded points (40% CL) Excluded points (68% CL)
			Excluded points (95% CL)					Excluded points (95% CL)
		200							200
	X [GeV]	150						X [GeV]	150
	m							m
		100							100
		50							50
		0 0	100	200	300 m Y [GeV]	400	500	600	0 0	100	200	300 m Y [GeV]	400	500	600	JHEP11(2016)111
								-

Table 6 .

 6 Figure 10. Constraints on the simplified top-philic dark matter model from the CMS 8 TeV monojet analysis[START_REF]Search for dark matter, extra dimensions and unparticles in monojet events in proton-proton collisions at √ s = 8 TeV[END_REF]. The top and dark matter couplings to the mediator are set to 4 while the mediator and dark matter masses are allowed to vary freely. Benchmarks used to investigate the differential distributions related to the CMS monojet analysis. The corresponding cross sections for a / E

	Monojet constraints (g X =g t =4; 8 TeV LHC)
	All points
	Excluded points (40% CL)
	Excluded points (68% CL)
	Excluded points (95% CL)
	(m Y , m X )	σ LO [pb]
	(100, 10) GeV	0.605
	(300, 10) GeV	0.194
	(100, 100) GeV 0.00261

T > 150 GeV selection are shown in the second column.

  2.1 for muons). It additionally

		Selection step	CMS	CMS i	MA5	MA5 i	δ rel i
	0	Nominal	224510		224510
	1	Preselection			15468.5 0.069
	2	/ E T > 320 GeV	4220.8		4579.8 0.296
	3	M T > 160 GeV	3390.1 0.803 3648.2 0.797 0.75%
	4 ∆Φ(j 1,2 , / E T ) > 1.2 2963.5 0.874 3124.3 0.856 2.06%
	5	M W T 2 > 200 GeV	2267.6 0.765	2403	0.769 -0.52%
							JHEP11(2016)111
			-40 -

Table 7 .

 7 Comparison of results obtained with our MadAnalysis 5 reimplementation (MA5) and those provided by the CMS collaboration (CMS). The efficiencies are defined in eq. (C.1) and the relative difference between the CMS and the MadAnalysis 5 results δ rel i in eq. (C.2).

Table 8 .

 8 Comparison of results obtained with our MadAnalysis 5 reimplementation (MA5) and those provided by the CMS collaboration (CMS). The efficiencies are defined in eq. (C.1) and the relative difference between the CMS and the MadAnalysis 5 results δ rel i in eq. (C.2).

		1 1						CMS CMS
								MA5 MA5
		0.1 0.1
	Events/25GeV Events/25GeV	0.01 0.01
		0.001 0.001
		0.0001 0.0001
		300 300	400 400	500 500	600 600	700 700	800 800	900 900	1000 1000
					E miss T E miss T	[GeV] [GeV]

The gluon form factor suffers from relatively large uncertainties on the strange quark content of the nucleons[START_REF] Alarcon | The strangeness content of the nucleon from effective field theory and phenomenology[END_REF], which we here omit.

Dark matter annihilation into two prompt photons is always suppressed by a factor 8α 2 e /9α 2 s with respect to annihilation into a pair of jets in the considered class of scenarios.

Acknowledgments

We would like to thank the LHC dark matter forum for motivation to pursue this project and for useful discussions during the April 2016 Amsterdam workshop. We are also grateful to Lana Beck for providing her four top analysis code for cross checks of our results, as well as to Yong Yang from the B2G CMS group for providing useful material to validate the recast of the B2G-14-004 analysis. This work has been supported by the research unit New physics at the LHC of the German research foundation DFG, the Helmholtz Alliance for Astroparticle Physics and the German Federal Ministry of Education and Research BMBF. CA is supported by the ATTRACT -Brains back to Brussels grant from Innoviris. MB is supported by the MOVE-IN Louvain Cofund grant. EC, BF and KM

A Mediator width

As supplementary material, figure 16 shows the relative mediator width Γ Y /m Y in the (g t , g X ) plane for different mass choices. The magnitude of the mediator width depends on the hierarchy among the different decay processes i.e. t t, X X, gg and γγ. Diphoton channel is negligible compared to the others and will not be discussed. Figure 16 shows that the Y 0 resonance can be considered as narrow (i.e. Γ Y /m Y < 0.03) when only the gg decay mode is involved. As soon as t t and X X decay channels are opened, the ratio Γ Y /m Y grows quickly, reaching 20% for g t , g X ∼ 2. The narrow width approximation is valid below couplings of O(1). Figure 17 shows the relative mediator width Γ Y /m Y in the (m Y , m X ) plane for different coupling choices. When g t , g X ≤ 1, Γ Y /m Y never exceeds 10% and the narrow width approximation is reliable for a wide region of the parameter space. In the kinematic regions where Y 0 decays to X and/or top quarks is allowed, increase in either g t or g X quickly leads to Γ Y /m Y ratio above 20%.