arxiv: v3 [hep-ph] 8 Jun 2017

Transcription

1 Prepared for submission o JHEP CERN-TH , DESY Simplified dark maer models wih wo Higgs doubles: I. Pseudoscalar mediaors arxiv: v3 [hep-ph] 8 Jun 2017 Marin Bauer, 1 Ulrich Haisch 2,3 and Felix Kahlhoefer 4 1 Insiu für Theoreische Physik, Universiä Heidelberg, Philosophenweg 16, Heidelberg, Germany 2 Rudolf Peierls Cenre for Theoreical Physics, Universiy of Oxford, OX1 3NP Oxford, Unied Kingdom 3 CERN, Theoreical Physics Deparmen, CH-1211 Geneva 23, Swizerland 4 DESY, Nokesraße 85, D Hamburg, Germany bauer@hphys.uni-heidelberg.de, ulrich.haisch@physics.ox.ac.uk, felix.kahlhoefer@desy.de Absrac: We sudy a new class of renormalisable simplified models for dark maer searches a he LHC ha are based on wo Higgs double models wih an addiional pseudoscalar mediaor. In conras o he spin-0 simplified models employed in analyses of Run I daa hese models are self-consisen, uniary and bounds from Higgs physics ypically pose no consrains. Predicions for various missing ransverse energy (E T,miss ) searches are discussed and he reach of he 13 TeV LHC is explored. I is found ha he proposed models provide a rich specrum of complemenary observables ha lead o non-rivial consrains. We emphasise in his conex he sensiiviy of he + E T,miss, mono-z and mono-higgs channels, which yield sronger limis han mono-je searches in large pars of he parameer space. Consrains from spin-0 resonance searches, elecroweak precision measuremens and flavour observables are also derived and shown o provide furher imporan handles o consrain and o es he considered dark maer models.

2 Conens 1 Inroducion 2 2 THDM plus pseudoscalar exensions Scalar poenial Yukawa secor 5 3 Anaomy of he parameer space Alignmen/decoupling limi Yukawa assignmens Di-au searches Di-op searches Flavour consrains EW precision consrains Perurbaiviy and uniariy 9 4 Parial decay widhs and branching raios Ligher pseudoscalar a Ligher scalar h Heavier scalar H Heavier pseudoscalar A Charged scalar H ± 15 5 Anaomy of mono-x signaures Mono-je channel /b b + E T,miss channels Mono-Z channel Mono-Higgs channel Mono-W channel 20 6 Numerical resuls Signal generaion Background esimaes Inerference effecs Summary plos LHC Run II reach 27 7 Conclusions 30 1

3 1 Inroducion Simplified models of dark maer (DM) and a single mediaor overcome many of he shorcomings of DM effecive field heories, bu remain general enough o represen a large class of popular heories of DM (see he reviews [1 3] for a complee lis of references). In paricular, including conribuions from on-shell producion of he mediaors allows o capure he full kinemaics of DM producion a colliders, making meaningful comparisons wih bounds from direc and indirec deecion experimens possible. In simplified DM models he ineracions beween he mediaors and he Sandard Model (SM) fermions are usually wrien as gauge or Yukawa couplings of mass dimension four. In many cases hese ineracions are however only apparenly renormalisable, because in a full SU(2) L U(1) Y gauge-invarian heory hey in fac arise from higher-dimensional operaors or hey signal he presence of addiional paricles or couplings ha are needed o resore gauge invariance [4 9]. These feaures can lead o parameer regions which are heoreically inaccessible or o misleading/unphysical predicions ofen relaed o uniariy violaion. Models in which he mediaors mix wih he SM bosons avoid such inconsisencies. The exising LEP and LHC measuremens of he Z-boson and Higgs-boson properies however severely resric he corresponding mixing angles, and as a resul classic E T,miss searches like mono-jes are ypically no he leading collider consrains in his class of simplified DM models [6, 10, 11]. In his aricle, we sudy a new class of simplified DM models for spin-0 mediaors based on wo Higgs double models (THDMs), which are an essenial ingredien of many wellmoivaed heories beyond he SM. In conras o iner THDMs, where he DM paricle is he lighes neural componen of he second Higgs double and is sabilised by an adhoc Z 2 symmery [12 15], our focus is on he case where he DM candidae is a SM single fermion. To couple he DM paricle o he SM, we inroduce a new spin-0 mediaor, which mixes dominanly wih he scalar or pseudoscalar parners of he SM Higgs. In his way consrains from Higgs signal srengh measuremens [16] can be saisfied and one obains a framework in which all operaors are gauge invarian and renormalisable. In wha follows we will explore he phenomenology of pseudoscalar mediaors, while scalar porals will be discussed in deail in an accompanying paper [17] (see also [18]). Pseudoscalar mediaors have he obvious advanage of avoiding consrains from DM direcdeecion experimens, so ha he observed DM relic abundance can be reproduced in large regions of parameer space and LHC searches are paricularly relevan o es hese models. Similar invesigaions of THDM plus pseudoscalar simplified DM models have been presened in [19 21]. Whenever indicaed we will highligh he similariies and differences beween hese and our work. The mono-x phenomenology of he considered simplified pseudoscalar models urns ou o be surprisingly rich. We examine he consrains from searches for j + E T,miss [22, 23], /b b + E T,miss [24 27], Z + E T,miss [28 30], h + E T,miss [31 34] and W + E T,miss [35, 36] and presen projecions for he 13 TeV LHC. In paricular, we provide benchmark scenarios ha are consisen wih bounds from elecroweak (EW) precision, flavour and Higgs observables including invisible decays [37, 38]. For he simplified pseudoscalar model recommended 2

4 g g a g g H a Z g g A a h Figure 1. Assored diagrams ha give rise o a + E T,miss (lef), Z + E T,miss (middle) and h + E T,miss (righ) signal in he simplified pseudoscalar model considered in our work. The exchanged spin-0 paricles are of scalar (H) or pseudoscalar (a, A) ype. Furher Feynman graphs ha conribue o he differen mono-x channels can be found in Figures 7 o 11. by he ATLAS/CMS DM Forum (DMF) [3] consrains from mono-je searches dominae hroughou he parameer space [39], whereas for he model considered here + E T,miss, mono-z and mono-higgs searches yield compeiive bounds and ofen provide he leading consrains. See Figure 1 for an illusraion of he various E T,miss processes ha are of mos ineres in our simplified model. This complemenariy of differen searches is he resul of he consisen reamen of he scalar secor, inducing gauge and rilinear scalar couplings of he mediaor beyond he ones presen in he DMF pseudoscalar model. I is paricularly appealing ha he Z +E T,miss and h+e T,miss signaures are sronges in he heoreically bes moivaed region of parameer space, where he couplings of he ligh Higgs are SM-like. In his region of parameer space, couplings of he new scalar saes o SM gauge bosons are srongly suppressed and play no role in he phenomenology, leading o gluon-fusion dominaed producion and a very predicive paern of branching raios. In consequence, a complemenary search sraegy can be advised, wih he exciing possibiliy o observe DM simulaneously in a number of differen channels, some of which are no limied by sysemaic errors and can be improved by saisics even beyond 300 fb 1 of luminosiy. The imporance of di-op resonance searches [40, 41] o probe neural spin-0 saes wih masses above he hreshold is also sressed, and i is poined ou ha for model realisaions wih a ligh scalar parner of he SM Higgs, di-au resonance searches should provide relevan consrains in he near fuure. We finally commen on he impac of boom-quark (b b) iniiaed producion. This paper is srucured as follows. In Secion 2 we describe he class of simplified DM models ha we will sudy hroughou our work, while Secion 3 conains a comprehensive review of he non-e T,miss consrains ha have o be saisfied in order o make a given model realisaion phenomenologically viable. The parial decay widhs and he branching raios of he spin-0 paricles arising in he considered simplified DM models are sudied in Secion 4. The mos imporan feaures of he resuling E T,miss phenomenology are described in Secion 5. In Secion 6 we finally presen he numerical resuls of our analyses providing summary plos of he mono-x consrains for several benchmark scenarios. The resul-oriened reader migh wan o skip direcly o his secion. Our conclusions and a brief oulook are given in Secion 7. 3

5 2 THDM plus pseudoscalar exensions In his secion we describe he srucure of he simplified DM model wih a pseudoscalar mediaor. We sar wih he scalar poenial and hen consider he Yukawa secor. In boh cases we will poin ou which are he new parameers corresponding o he ineracions in quesion. 2.1 Scalar poenial The ree-level THDM scalar poenial ha we will consider hroughou his paper is given by he following expression (see for example [42, 43]) ( ) V H = µ 1 H 1 H 1 + µ 2 H 2 H 2 + µ 3 H 1 H ( 2 + h.c. + λ 1 H 1 H ) 2 ( 1 + λ2 H 2 H ) 2 2 ( + λ 3 H 1 H )( 1 H 2 H ) ( 2 + λ4 H 1 H )( 2 H 2 H ) ( 1 + [λ 5 H 1 H ) ] h.c.. (2.1) Here we have imposed a Z 2 symmery under which H 1 H 1 and H 2 H 2 o suppress flavour-changing neural currens (FCNCs), bu allowed for his discree symmery o be sofly broken by he erm µ 3 H 1 H 2 + h.c. The vacuum expecaion values (VEVs) of he Higgs doubles are given by H i = (0, v i / 2) T wih v = v1 2 + v GeV and we define an β = v 2 /v 1. To avoid possible issues wih elecric dipole momens, we assume ha he mass-squared erms µ j, he quaric couplings λ k and he VEVs are all real and as a resul he scalar poenial as given in (2.1) is CP conserving. The hree physical neural Higgses ha emerge from V H are in such a case boh mass and CP eigensaes. The mos economic way o couple fermionic DM o he SM hrough pseudoscalar exchange is by mixing a CP-odd mediaor P wih he CP-odd Higgs ha arises from (2.1). This can be achieved by considering he following ineracion erms V P = 1 2 m2 P P 2 + P ( ) ( ) ib P H 1 H 2 + h.c. + P 2 λ P 1 H 1 H 1 + λ P 2 H 2 H 2, (2.2) where m P and b P are parameers wih dimensions of mass. We assume ha V P does no break CP and hus ake b P o be real in he following. In his case P does no develop a VEV and remains a pure CP eigensae. Neverheless, his erm does lead o a sof breaking of he Z 2 symmery. Noice ha compared o [19 21] which include only he rilinear poral coupling b P, we also allow for quaric poral ineracions proporional o λ P 1 and λ P 2. A quaric self-coupling of he form P 4 has insead no been included in (2.2), as i does no lead o any relevan effec in he observables sudied in our paper. The ineracions in he scalar poenial (2.1) mix he neural CP-even weak eigensaes and we denoe he corresponding mixing angle by α. The poral coupling b P appearing in (2.2) insead mixes he wo neural CP-odd weak eigensaes wih θ represening he associaed mixing angle. The resuling CP-even mass eigensaes will be denoed by h and H, while in he CP-odd secor he saes will be called A and a, where a denoes he exra degree of freedom no presen in THDMs. The scalar specrum also conains wo charged mass eigensaes H ± of idenical mass. 4

6 Diagonalising he mass-squared marices of he scalar saes leads o relaions beween he fundamenal parameers enering V H and V P. These relaions allow o rade he parameers m P, µ 1, µ 2, µ 3, b P, λ 1, λ 2, λ 4, λ 5 for sines and cosines of mixing angles, VEVs and he masses of he physical Higgses. This procedure ensures in addiion ha he scalar poenial is posiive definie and ha he vacuum soluion is an absolue minimum. In he broken EW phase he physics of (2.1) and (2.2) is hence fully capured by he angles α, β, θ, he EW VEV v, he quaric couplings λ 3, λ P 1, λ P 2 and he masses M h, M H, M A, M H ±, M a. We will use hese parameers as inpu in our analysis. 2.2 Yukawa secor The couplings beween he scalars and he SM fermions are resriced by he sringen experimenal limis on flavour observables. A necessary and sufficien condiion o avoid FCNCs associaed o neural Higgs ree-level exchange is ha no more han one of he Higgs doubles couples o fermions of a given charge [44, 45]. This so-called naural flavour conservaion hypohesis is auomaically enforced by he aforemenioned Z 2 symmery acing on he doubles, if he righ-handed fermion singles ransform accordingly. The Yukawa couplings are explicily given by L Y = ( QY i u Hi u R + QY d i H id R + LY ) l i H il R + h.c.. (2.3) i=1,2 Here Yf i are Yukawa marices acing on he hree fermion generaions and we have suppressed flavour indices, Q and L are lef-handed quark and lepon doubles, while u R, d R and l R are righ-handed up-ype quark, down-ype quark and charged lepon singles, respecively. Finally, Hi = ɛhi wih ɛ denoing he wo-dimensional anisymmeric ensor. The naural flavour conservaion hypohesis can be saisfied by four discree assignmens, where by convenion up-ype quarks are always aken o couple o H 2 : Yu 1 = Yd 1 = Y l 1 = 0, (ype I), Yu 1 = Yd 2 = Y l 2 = 0, (ype II), Yu 1 = Yd 1 = Y l 2 = 0, (ype III), Yu 1 = Yd 2 = Y l 1 = 0, (ype IV). (2.4) The dependence of our resuls on he choice of he Yukawa secor will be discussed in some deail in he nex secion. Taking DM o be a Dirac fermion χ a separae Z 2 symmery under which χ χ can be used o forbid a coupling of he form L H 1 χ R + h.c. A he level of renormalisable operaors his leaves L χ = iy χ P χγ 5 χ, (2.5) as he only possibiliy o couple he pseudoscalar mediaor P o DM. In order o no violae CP we require he dark secor Yukawa coupling y χ o be real. The parameer y χ and he DM mass m χ are furher inpu parameers in our analysis. 5

7 3 Anaomy of he parameer space In his secion we examine he anaomy of he parameer space of he model inroduced above and discuss a number of imporan simplificaions. We briefly explain he alignmen/decoupling limi and describe he dependence of he predicions on he choice of Yukawa secor. The consrains on he mixing angles, quaric couplings and Higgs masses from spin-0 resonance searches, flavour physics, EW precision measuremens, perurbaiviy and uniariy are also elucidaed. 3.1 Alignmen/decoupling limi Afer EW symmery breaking he kineic erms of he Higgs fields H i lead o ineracions beween he CP-even mass eigensaes and he massive EW gauge bosons. These ineracions ake he form L ( sin (β α) h + cos (β α) H) ( 2MW 2 W µ + W µ + M 2 ) Z v v Z µz µ. (3.1) In order o simplify he furher analysis, we concenrae on he well-moivaed alignmen/decoupling limi of he THDM where α = β π/2. In his case sin (β α) = 1 meaning ha he field h has SM-like EW gauge boson couplings. I can herefore be idenified wih he boson of mass M h 125 GeV discovered a he LHC and he consrains from he Run I combinaion of he ATLAS and CMS measuremens of he Higgs boson producion and decay raes o SM final saes [16] are readily fulfilled. Noice ha in he alignmen/decoupling limi he scalar H does no inerac wih W -boson or Z-boson pairs a ree level because in his limi one has cos (β α) = Yukawa assignmens Working in he alignmen/decoupling limi he fermion-scalar ineracions mos relevan for he furher discussion are given by L y [ h + ξf M H i ξm f 2 f=b,τ y [ f f h + ξf M H + i ξm f 2 ( cos θ A sin θ a ) γ5 ] ( cos θ A sin θ a ) γ5 ] f y 2 V b ξ M H + R b L + y b 2 V b ξ M b H+ L b R + h.c. ) iy χ (sin θ A + cos θ a χγ 5 χ, (3.2) where y f = 2m f /v denoe he SM Yukawa couplings and V ij are he elemens of he Cabibbo-Kobayashi-Maskawa (CKM) marix. The couplings ξf M encode he dependence on 6

8 he choice of Yukawa secor (2.4). In erms of an β one has ξ I = ξ I b = ξi τ = co β, (ype I), ξ II = co β, ξb II = ξii τ = an β, (ype II), ξ III ξ IV = ξ III b = co β, ξ III τ = an β, (ype III), = ξ IV τ = co β, ξ IV b = an β, (ype IV). (3.3) Since he producion of he pseudoscalar mediaor a as well as pp h, H, A is driven by op-quark loops ha ener he gluon-fusion (gg) channel a he LHC (see for insance [46] for a discussion in he conex of E T,miss searches) we will in he following focus on he region of small an β. In his limi he couplings of H, A, a o down-ype quarks and charged lepons in (3.2) are srongly Yukawa suppressed irrespecively of he chosen Yukawa assignmen (3.3). As a resul exising bounds on he neural scalar masses from flavour observables such as B s µ + µ ha are known o receive an β enhanced correcions [49] are wihin experimenal limis [50] even for a ligh scalar specrum. 3.3 Di-au searches In order o undersand wheher he exising LHC searches for heavy neural Higgses in fermionic final saes such as f f = τ + τ, b b pose consrains on he low an β region of our simplified model, i is imporan o realise ha while he pseudoscalars A and a couple boh o DM, he heavy scalar H does no, as can be seen from (3.2). If he channels A/a χ χ are open, he discovery poenial for H f f is herefore generically larger han ha for he corresponding pseudoscalar modes. In fac, he consrains from pp H f f are mos sringen for model realisaions wih M H < 2m and M a > max (M H M Z, M H /2), so ha he decays H, H aa and H az are kinemaically forbidden and in consequence H is forced o decay o ligh SM fermions (see Secion 4.3). The ypical resricions ha resul from LHC searches for heavy scalars can be illusraed by considering M H = 300 GeV and employing he 95% confidence level (CL) limi σ (pp H) BR (H τ + τ ) < 0.4 pb [47, 48] ha is based on 13 fb 1 of 13 TeV daa. Using he nex-o-nex-o-nex-o-leading order resuls [51] for inclusive H producion in gluon fusion, we hen find ha he curren di-au searches only exclude a narrow sliver of parameers in he M a an β plane wih 0.55 an β 0.65 and M a 210 GeV in he case of a Yukawa secor of ype II. A reducion of he quoed upper limi on he producion cross secion imes branching raio o 0.2 pb would however improve he range of excluded an β values o 0.3 an β 1.2. As we will see in Secion 6.4, such a consrain would be very valuable because probing models wih an β = O(1) and M H M a 300 GeV urns ou o be difficul by oher means. 3.4 Di-op searches Heavy scalar and pseudoscalar bosons decaying dominanly ino op-quark pairs can be searched for by sudying he resuling invarian mass specra m. In conras o di-op searches for spin-1 or spin-2 saes, a peak in he m disribuion ha one generically 7

9 expecs in he narrow-widh approximaion (NWA) is however no he only signaure of a spin-0 resonance in his case. Indeed, he gg H/A signal will inerfere wih he QCD background which a he LHC is mainly generaed by he gluon-fusion channel gg. The signal-background inerference will depend on he CP naure of he inermediae spin-0 boson, is mass and is oal decay widh. The observed inerference paern can be eiher consrucive or desrucive, leading o a raher complex signaure wih a peak-dip srucure in he m specrum [52, 53]. The pp H/A channel provides hence an ineresing bu challenging opporuniy for hadron colliders o search for addiional spin-0 bosons (see for insance [54, 55] for recen phenomenological discussions). The firs LHC analysis ha akes ino accoun inerference effecs beween he signal process gg H/A and he SM background gg is he ATLAS search [40]. I is based on 20.3 fb 1 of 8 TeV LHC daa and considers he m specrum in final saes wih a single charged lepon (elecron or muon), large E T,miss and a leas four jes. The search resuls are inerpreed in he conex of a pure THDM of ype II for wo differen mass poins and employ he alignmen/decoupling limi, i.e. sin (β α) = 1. For a neural scalar H (pseudoscalar A) wih a mass of 500 GeV, he ATLAS analysis excludes he parameer values an β < 0.45 (an β < 0.85) a he 95% CL, while for he 750 GeV mass poin no meaningful consrain on an β can be se. Recasing hese limis ino bounds on he parameer space of spin-0 simplified DM models is sraighforward [41] and we will analyse he resuling resricions on our model in Secion Flavour consrains Indirec consrains on he charged Higgs-boson mass M H ± arise from Z b b [56 58], B X s γ [59 61] and B q B q mixing [62 65] since he laer processes receive correcions from he H + R b L +h.c. and H + L b R +h.c. erms in (3.2). We find ha B X s γ provides he sronges indirec consrain on M H ± for small an β values in models of ype I and III a presen, while B s B s oscillaions represen he leading indirec consrain in he oher wo cases. For M H ± = 750 GeV we obain he bound an β 0.8 from a combinaion of B-meson physics observables irrespecive of he choice of he Yukawa secor. A modelindependen lower limi of an β 0.3 can also be obained from he requiremen ha he op-quark Yukawa coupling remains perurbaive [43]. The laes LHC search limis on he charged Higgs mass in he pp bh ± (H ± b) channel [66, 67] are saisfied for an β 0.2 if M H ± = 750 GeV is assumed, and herefore provide no relevan consrain. 3.6 EW precision consrains A scalar poenial wih wo doubles such as he one inroduced in (2.1) leads o addiional Higgs ineracions compared o he SM, which can violae he cusodial symmery presen in he SM Higgs secor. I can be shown [68 72] ha he ree-level poenial V H is cusodially invarian for M A = M H ± or M H = M H ±. Only in hese wo cases can H or A have a sizeable mass spliing from he res of he Higgses wihou being in conflic wih EW precision measuremens, mos imporanly ρ. Since he poenial (2.2) mixes he pseudoscalar degree of freedom in H i wih P, in he heory described by V H + V P here are however addiional sources of cusodial symmery breaking compared o he case of he pure THDM. 8

10 In he alignmen/decoupling limi and aking M A = M H ±, we find ha he exended scalar secor gives rise o he following one-loop correcion ρ = 1 M 2 [ ] H ± (4π) 2 v f(m H, M a, M H ±) + f(m a, M H, M H ±) sin 2 θ, (3.4) wih m 4 ( 1 m 2 f(m 1, m 2, m 3 ) = 2 m 2 ) ( ) 3 m 2 ( ( ) m 2 3 m 2 1 m2) 2 m 2 1 m 2 ln 1 3 m 2. (3.5) 3 Noice ha ρ 0 in he limi sin θ 0 in which he wo CP-odd weak eigensaes are also mass eigensaes or if he scalar mass specrum is fully degenerae. In he alignmen/decoupling limi wih M H = M H ±, he cusodial symmery is insead no broken by V H + V P and as a resul one has ρ = 0 a he one-loop level. From he above discussion i follows ha only cases wih M A = M H ± are subjec o he consrains from he EW precision measuremens, while scenarios wih M H = M H ± are no. In order o derive he resuling consrains in he former case, we employ he 95% CL bound ρ [ 1.2, 2.4] 10 3, (3.6) which corresponds o he value exraced in [73] from a simulaneous deerminaion of he Peskin-Takeuchi parameers S, T and U. The fac ha (3.4) is proporional o he produc of mass differences M H ± M H and M H ± M a as well as sin 2 θ implies ha he exising EW precision daa allow o se sringen bounds on sin θ if he relevan mass spliings in he scalar secor are sizeable. Taking for insance M H ± = 750 GeV and M a = 65 GeV, we find ha for M H = 500 GeV (M H = 300 GeV) he inequaliy sin θ < 0.35 (sin θ < 0.25) has o be saisfied in order o be compaible wih (3.6). We will see in Secion 4.3 ha he resricions on sin θ can have a visible impac on he decay paern of he scalar H, which in urn affecs he mono-z phenomenology discussed in Secion Perurbaiviy and uniariy Perurbaiviy [74, 75] and uniariy [76 79] also pu resricions on he scalar masses and he magniudes and signs of he quaric couplings. In our numerical analysis we will resric our aenion o he parameer space ha saisfies M H, M A, M a M H ± = O(1 TeV) and always keep λ 3, λ P 1 and λ P 2 of O(1) or below. For such inpu parameer choices all consrains discussed in his secion are saisfied if an β is no oo far below 1. We also only consider parameers for which he oal decay widhs of H and A are sufficienly small so ha he NWA applies, i.e. Γ i M i /3 for i = H, A. This requiremen ses an upper limi on he mass of he charged Higgs boson ha is ofen sronger han bounds from perurbaiviy. 4 Parial decay widhs and branching raios This secion is devoed o he discussion of he parial decay widhs and he branching raios of he spin-0 paricles arising in he simplified DM model inroduced in Secion 2. For concreeness we will focus on he alignmen/decoupling limi of he heory. We will 9

11 furhermore pay special aenion o he parameer space wih a ligh DM paricle, small values of an β and scalar specra where he new pseudoscalar a and he scalar h are he lighes degrees of freedom. 4.1 Ligher pseudoscalar a As a resul of CP conservaion he field a has no couplings of he form aw + W, azz and ahh. In conras he ahz verex is allowed by CP symmery bu vanishes in he alignmen/decoupling limi. A ree level he pseudoscalar a can hus only decay ino DM paricles and SM fermions. The corresponding parial decay widhs are given by Γ (a χ χ) = y2 χ 8π M aβ χ/a cos 2 θ, Γ ( a f f ) = N c f ( ) ξ M 2 f 8π m 2 f v 2 M aβ f/a sin 2 θ, (4.1) where β i/a = 1 τ i/a is he velociy of he paricle i in he res frame of he final-sae pair and we have defined τ i/a = 4m 2 i /M a 2. Furhermore Nc f = 3 (1) denoes he relevan colour facor for quarks (lepons) and he explici expressions for he couplings ξf M can be found in (3.3). A he loop level he pseudoscalar a can also decay o gauge bosons. The larges parial decay widh is he one o gluon pairs. I akes he form Γ (a gg) = α2 s 32π 3 v 2 M a 3 q=,b,c ξq M f(τ q/a ) 2 sin 2 θ, (4.2) wih ( ) 1 f(τ) = τ arcan 2. (4.3) τ 1 For small an β and non-zero values of sin θ he couplings of a o DM and op quarks dominae over all oher couplings. As a resul, he decay paern of a is in general very simple. This is illusraed in he panels of Figure 2 for wo differen choices of parameer ses. The lef panel shows he branching raio of a for a very ligh DM paricle wih m χ = 1 GeV. One observes ha below he hreshold one has BR (a χ χ) = 100% while for M a > 2m boh decays o DM and op-quarks pairs are relevan. In fac, sufficienly far above he hreshold one obains BR (a χ χ) /BR (a ) 0.7 y 2 χ an 2 β/ an 2 θ independen of he specific realisaion of he Yukawa secor. In he righ panel we presen our resuls for a DM sae of m χ = 100 GeV. In his case we see ha below he χ χ hreshold he pseudoscalar a decays dominanly ino boom-quark pairs bu ha also he branching raios o aus and gluons exceed he percen level. Compared o he lef plo one also observes ha in he righ plo he raio BR (a χ χ) /BR (a ) is significanly larger for M a > 2m due o he differen choice of sin θ. 4.2 Ligher scalar h For sufficienly heavy pseudoscalars a he decay paern of h resembles ha of he SM Higgs boson in he alignmen/decoupling limi. For M a < M h /2 on he oher hand decays o wo 10

12 sinθ = 1/ 2, m χ = 1 GeV sinθ = 1/2, m χ = 100 GeV χχ bb χχ BR (a X ) BR (a X ) cc gg τ + τ Figure 2. Branching raios of he ligher pseudoscalar a as a funcion of is mass for wo differen choices of sin θ and m χ as indicaed in he headline of he plos. The oher relevan parameers have been se o an β = 1, M H = M A = M H ± = 750 GeV and y χ = 1. Noice ha for his specific an β value he branching raios of he pseudoscalar a do no depend on he choice of Yukawa secor. on-shell a mediaors are possible. The corresponding parial decay widh reads wih Γ (h aa) = 1 32π g2 haa M hβ a/h, (4.4) g haa = 1 [ (M 2 M h v h 2MH 2 + 4MH 2 2M 2 ± a 2λ 3 v 2) sin 2 θ 2 ( λ P 1 cos 2 β + λ P 2 sin 2 β ) ] v 2 cos 2 θ. (4.5) Noice ha he haa coupling conains erms proporional o boh sin 2 θ and cos 2 θ. These conribuions resul from he rilinear and quaric couplings in he scalar poenial (2.2), respecively. In our THDM plus pseudoscalar exension, h aa decays are even possible in he limi θ 0, which is no he case in he simplified model considered in [19 21]. Since he oal decay widh of he SM Higgs is only abou 4 MeV, hree-body decays of h ino final saes wih a single a can also be relevan in he mass range M h /2 < M a M h. Phenomenologically he mos imporan hree-body decay is he one where a is accompanied by a pair of DM paricles bu decays o an a and SM fermions are also possible. The corresponding parial decay widhs are given by y2 χ Γ (h aχ χ) = 32π 3 g2 haa M hβ χ/a g(τ a/h ) cos 2 θ, Γ ( h af f ) = N c f ( ) ξ M 2 f m 2 f 32π 3 v 2 g2 haa M hβ f/a g(τ a/h ) sin 2 θ, (4.6) 11

13 m χ = 1 GeV m χ = 20 GeV BR (h X ) aχχ bb W + W - gg τ + τ - cc BR (h X ) aχχ bb W + W - gg τ + τ - cc 10-2 ZZ ZZ abb Figure 3. Branching raios of he ligher scalar h as a funcion of he pseudoscalar mass M a for wo differen choices of m χ as indicaed in he headline of he plos. The oher relevan parameers have been se o an β = 1, M H = M A = M H ± = 750 GeV, sin θ = 1/ 2, λ 3 = λ P 1 = λ P 2 = 0 and y χ = 1. wih [80] g(τ) = 1 ( τ )] [4 8 (τ 4) ln 4 5τ 4 [ ( ) ( )] τ arcan τ 1 2 arcan. τ 1 τ 1 (4.7) In Figure 3 we show he branching raios of he SM Higgs h for wo differen values of he DM mass. We observe ha for a ligh pseudoscalar mediaor a one has in boh cases BR (h aχ χ) = 100%. In fac, he oal decay widh of he ligher scalar h exceeds 3 GeV for masses M a 70 GeV. Such large values of Γ h are in conflic wih he model-independen upper limis on he oal decay widh of he Higgs as measured by boh ATLAS and CMS in LHC Run I [81, 82]. Noice ha since he pseudoscalar a decays wih 100% o DM pairs for he considered values of m χ one has BR (h aχ χ) = BR (h 2χ2 χ). This implies ha for ligh DM he simplified model presened in Secion 2 is subjec o he consrains arising from invisible decays of he Higgs boson [37, 38]. We will analyse he resuling resricions on he parameer space in Secion 6.4. The righ panel finally illusraes ha in cases where m χ is close o a quarer of he SM Higgs mass also decays such as h ab b wih a χ χ can have branching raios of a few percen (or more) for a narrow range of M a values. Noice ha for he choice an β = 1 used in he figure he resul for BR ( h ab b ) does no depend on he paricular Yukawa assignmen. 4.3 Heavier scalar H In he alignmen/decoupling limi of he pseudoscalar exensions of he THDM model he heavier scalar H does no couple o W + W and ZZ pairs. In addiion he Hhh verex 12

14 vanishes. Under he assumpion ha M H > M a and aking A o be sufficienly heavy, he scalar H can hence decay only o SM fermions or he aa and az final sae a ree level. The corresponding parial decay widhs are wih Γ ( H f f ) = N c f ( ) ξ M 2 f 8π m 2 f v 2 M H β 3 f/h, Γ (H aa) = 1 32π g2 Haa M H β a/h, Γ (H az) = 1 λ 3/2 (M H, M a, M Z ) 16π MH 3 sin 2 θ, v2 g Haa = 1 [ co (2β) ( 2Mh 2 M H v 4M H 2 + 4MH 2 2λ ± 3 v 2) sin 2 θ ] + sin (2β) (λ P 1 λ P 2 ) v 2 cos 2 θ, denoing he Haa coupling. We have furhermore inroduced (4.8) (4.9) λ(m 1, m 2, m 3 ) = ( m 2 1 m 2 2 m 2 3) 2 4m 2 2 m 2 3, (4.10) which characerises he wo-body phase space for hree massive paricles. Noice ha he appearance of λ P 1 and λ P 2 in he parial decay widh Γ (H aa) indicaes again a qualiaive difference beween he scalar ineracions considered in [19 21] and he more general poenial (2.2). A he one-loop level he heavier scalar H can in addiion decay o gluons and oher gauge bosons, bu he associaed branching raios are very suppressed and hus have no impac on our numerical resuls. The dominan branching raios of H as a funcion of M a are displayed in Figure 4 for wo parameer ses. In he lef panel he case of a scalar H wih sin θ = 1/ 2 and M H = 750 GeV is shown. One observes ha for M a 350 GeV he decay mode H az has he larges branching raio, while for heavier a he H channel represens he leading decay. Noice ha for model realisaions where he decay channel H az dominaes, ineresing mono-z signaures can be expeced [20, 21]. We will come back o his poin in Secion 5.3. The decay paern of H is however srongly dependen on he mass of H since for M H < M H ± he mixing angle θ is consrained o be small by EW precision measuremens (see Secion 3.6). This behaviour is easy o undersand from (4.8) which in he limi of small sin θ, an β = O(1) and large M H imply ha Γ (H ) m 2 /(M H an 2 β), Γ (H aa) v 4 /MH 3 (λ P 1 λ P 2 ) 2 and Γ (H az) M H sin 2 θ. For M H > 2m he decay mode H can hence dominae over he whole M a range of ineres. This feaure is illusraed on he righ-hand side of he figure for sin θ = 0.35 and M H = 500 GeV. One also sees from his panel ha he branching raio of H aa can be relevan as i does no end o zero in he sin θ 0 limi if he combinaion λ P 1 λ P 2 of quaric couplings is non-zero. For an β 2 and λ P 1 λ P 2 1, BR (H aa) can even be he larges branching raio for M a < M H /2. This happens because he erms proporional o sin 2 θ and cos 2 θ in (4.9) boh give a sizeable conribuion o he Haa coupling, while he H coupling is suppressed by 1/ an 2 β. 13

15 sinθ = 1/ 2, M H = 750 GeV sinθ = 0.35, M H = 500 GeV az 0.50 BR (H X ) BR (H X ) aa az Figure 4. Branching raios of he heavier scalar H as a funcion of M a for wo differen choices of sin θ and M H as indicaed in he headline of he plos. The oher used inpu parameers are an β = 1, M A = M H ± = 750 GeV, λ 3 = λ P 2 = 0 and λ P 1 = Heavier pseudoscalar A For M A > M a and assuming ha decays o H are kinemaically inaccessible, he pseudoscalar A can only decay o DM, SM fermions and he ah final sae a ree level. In he alignmen/decoupling limi he corresponding parial decay widhs ake he form wih Γ (A χ χ) = y2 χ 8π M Aβ χ/a sin 2 θ, Γ ( A f f ) = N c f ( ) ξ M 2 f 8π m 2 f v 2 M Aβ f/a cos 2 θ, Γ (A ah) = 1 λ 1/2 (M A, M a, M h ) gaah 2 16π M, A g Aah = 1 [ Mh 2 M A v 2M H 2 MA 2 + 4MH 2 M 2 ± a 2λ 3 v ( λ P 1 cos 2 β + λ P 2 sin 2 β ) v 2 ] sin θ cos θ, (4.11) (4.12) denoing he Aah coupling, and he analyic expression for he wo-body phase-space funcion λ(m 1, m 2, m 3 ) can be found in (4.10). Like in he case of H, loop-induced decays of he heavier pseudoscalar A can be negleced for all pracical purposes. In Figure 5 we presen our resuls for he branching raios of he pseudoscalar A as a funcion of M a for wo differen parameer choices. The lef panel illusraes he case M A = 750 GeV and one sees ha for such an A he branching raios are all above 10% and he hierarchy BR (A ah) > BR (A ) > BR (A χ χ) is observed for M a 200 GeV. 14

16 sinθ = 1/ 2, M A = 750 GeV sinθ = 1/2, M A = 500 GeV 0.50 ah 0.50 ah BR (A X ) χχ BR (A X ) χχ Figure 5. Branching raios of he heavier pseudoscalar A as a funcion of M a for wo differen choices of M A and sin θ as indicaed in he headline of he plos. The oher parameer choices are an β = 1, M H = M H ± = 750 GeV, λ 3 = λ P 1 = λ P 2 = 0, y χ = 1 and m χ = 1 GeV. As shown on he righ-hand side of he figure, his hierarchy no only remains inac bu is even more pronounced for a moderaely heavy A unil he hreshold M a = M A M h is reached. For larger M a values only decays o χ χ and final saes maer and he raio of heir branching raios is approximaely given by BR (A χ χ) /BR (A ) 0.9 yχ 2 an 2 β an 2 θ irrespecive of he paricular Yukawa assignmen. Noice ha a sizeable A ah branching raio is a generic predicion in he THDM plus pseudoscalar exensions wih small an β, since he charged Higgs has o be quie heavy in his case in order o avoid he bounds from B X s γ and/or B s -meson mixing. Since a χ χ is ypically he dominan decay mode of he ligher pseudoscalar a, appreciable mono-higgs signals are hence a firm predicion in a cerain region of parameer space of our simplified model. This poin will be furher explained in Secion Charged scalar H ± Since in he alignmen/decoupling limi he H + hw + verex vanishes, he parial decay widhs of he charged scalar H + ha are relevan in he small an β regime read Γ ( H + b ) = N c V b 2( ξ M ) 2 m 2 ( ) 2 8π v 2 M H ± 1 m2 MH 2, ± Γ ( H + HW +) = 1 λ 3/2 (M H ±, M H, M W ) 16π MH 3 v 2, ± Γ ( H + AW +) = 1 λ 3/2 (M H ±, M A, M W ) 16π MH 3 v 2 cos 2 θ, ± (4.13) Γ ( H + aw +) = 1 λ 3/2 (M H ±, M a, M W ) 16π MH 3 v 2 sin 2 θ, ± 15

17 sinθ = 0.35, M H = 500 GeV sinθ = 1/ 2, M A = 500 GeV BR (H + X ) aw + b HW + BR (H + X ) aw + AW + b Figure 6. Branching raios of he charged scalar H + as a funcion of M a for wo differen ses of inpu parameers as indicaed in he headline of he plos. In he lef (righ) panel in addiion an β = 1 and M A = M H ± = 750 GeV (M H = M H ± = 750 GeV) is used. where in he case of H + b we have negleced erms of O(m 2 b /M H 2 ) in he expression ± for he parial decay widh. Noice ha in THDMs of ype II and III also he decay H + τ + ν τ can be imporan if an β 1. The resul for Γ (H + τ + ν τ ) can be obained from he expression given above for Γ ( H + b ) by obvious replacemens. The main branching raios of he charged Higgs H + are displayed in Figure 6. On he lef-hand side of he figure he case of sin θ = 0.35 and M H = 500 GeV is displayed and one observes ha BR ( H + b ) > BR (H + HW + ) > BR (H + aw + ) for he shown values of M a. Noice ha for scenarios wih M H < M H ± he hierarchy BR (H + HW + ) > BR (H + aw + ) is a raher model-independen predicion since in such cases EW precision measuremens require sin θ o be small and Γ (H + aw + ) /Γ (H + HW + ) sin 2 θ. The same is no rue for he hierarchy beween BR ( H + b ) and BR (H + HW + ) which depends sensiively on he choice of an β since Γ ( H + b ) /Γ (H + HW + ) 1/ an 2 β. I follows ha for values of an β > 1 he H + HW + channel can also be he dominan decay mode. In model realisaions wih M A < M H ± here are no consrains from ρ on sin θ and in urn he H + aw + branching raio can dominae for sufficienly large mixing in he pseudoscalar secor. This feaure is illusraed by he righ panel in he figure using sin θ = 1/ 2 and M A = 500 GeV. For his choice of inpu parameers we find ha BR (H + aw + ) > BR ( H + b ) for masses M a 300 GeV. Since he pseudoscalar a predominanly decays via a χ χ i follows ha THDM plus pseudoscalar exensions wih M A < M H ± can lead o a resonan mono-w signal. We will discuss he LHC prospecs for he deecion of such a E T,miss signaure in Secion

18 5 Anaomy of mono-x signaures In his secion we will discuss he mos imporan feaures of he mono-x phenomenology of he pseudoscalar exensions of he THDM. We examine he mono-je, he + E T,miss, he mono-z and he mono-higgs signaure. The b b + E T,miss and mono-w channel are also briefly considered. Our numerical analysis of he mono-x signals is posponed o Secion Mono-je channel A firs possibiliy o search for pseudoscalar ineracions of he form (3.2) consiss in looking for a mono-je signal, where he mediaors ha pair produce DM are radiaed from heavyquark loops [39, 46, 83 89]. Represenaive examples of he possible one-loop Feynman diagrams are shown in Figure 7. For m a > 2m χ and M A M a only graphs involving he exchange of he ligh pseudoscalar a will conribue o he j + E T,miss signal. As a resul he normalised kinemaic disribuions of he mono-je signal in he pseudoscalar exensions of he THDM are idenical o hose of he DMF pseudoscalar model. Working in he NWA and assuming ha an β is small, he raio of he fiducial cross secions in he wo models is hus approximaely given by he simple expression ( ) σ (pp j + E T,miss ) yχ sin θ 2. (5.1) σ (pp j + E T,miss ) DMF g χ g q an β Here g χ (g q ) denoes he DM-mediaor (universal quark-mediaor) coupling in he corresponding DMF spin-0 simplified model. Noice ha he above relaion is largely independen of he choice of Yukawa secor as long as an β = O(1) since boom-quark loops have only an effec of a few percen on he j + E T,miss disribuions (see for insance [90] for a relaed discussion in he conex of Higgs physics). Using he approximaion (5.1) i is sraighforward o recas exising mono-je resuls on he DMF pseudoscalar model such as hose given in [23] ino he THDM plus pseudoscalar model space. The numerical resuls presened in he nex secion however do no employ any approximaion since hey are based on a calculaion of he j + E T,miss cross secions including boh op-quark and boom-quark loops as well as he exchange of boh a and A mediaors. 5.2 /b b + E T,miss channels A second channel ha is known o be a sensiive probe of op-philic pseudoscalars wih large invisible decay widhs is associaed producion of DM and pairs [39, 86, 89, 91 94]. Figure 8 displays examples of ree-level diagrams ha give rise o a + E T,miss signaure in he pseudoscalar exensions of he THDM model. In he case ha A is again much heavier han a, he signal srengh for + E T,miss in our simplified model can be obained from he predicion in he DMF pseudoscalar scenario from a rescaling relaion analogous o he one shown in (5.1). Using such a simple recasing procedure we find ha he mos recen ATLAS [24] and CMS searches for + E T,miss [25] ha are based on 13.2 fb 1 and 2.2 fb 1 of 13 TeV LHC daa, respecively, only allow o se very weak bounds on an β. For insance for M a = 100 GeV, y χ = 1 and m χ = 1 GeV a 17

19 g g g a g g g a Figure 7. Examples of diagrams ha give rise o a j + E T,miss signaure hrough he exchange of a ligher pseudoscalar a. Graphs involving a heavier pseudoscalar A also conribue o he signal in he pseudoscalar exensions of he THDM bu are no shown explicily. lower limi of an β 0.2 is obained. The + E T,miss consrains on he parameer space of he pseudoscalar exensions of he THDM are however expeced o improve noably a forhcoming LHC runs. The numerical resuls ha will be presened in Secion 6.4 are based on he search sraegy developed recenly in [94] which employs a shape fi o he difference in pseudorapidiy of he wo charged lepons in he di-leponic channel of + E T,miss. Besides +E T,miss also b b+e T,miss producion [91, 92] has been advocaed as a sensiive probe of spin-0 poral couplings o heavy quarks. Recasing he mos recen 13 TeV LHC b b + E T,miss searches [26, 27] by means of a simple rescaling similar o (5.1) we find ha no relevan bound on he parameer space of our simplified model can be derived unless he ab b coupling is significanly enhanced. From (3.3) we see ha such an enhancemen can only arise in THDMs of ype II and IV, while i is no possible for he oher Yukawa assignmens. Since in he limi of large an β also direc searches for he ligh pseudoscalar a in final saes conaining boom quarks or charged lepons are relevan (and naively even provide he leading consrains) we do no consider he b b+e T,miss channel in wha follows, resricing our numerical analysis o he parameer space wih small an β. 5.3 Mono-Z channel A mono-x signal ha is srongly suppressed in he case of he spin-0 DMF models [88] bu will urn ou o be relevan in our simplified DM scenario is he mono-z channel [21]. A sample of one-loop diagrams ha lead o such a signaure are displayed in Figure 9. Noice ha he lef diagram in he figure allows for resonan Z + χ χ producion hrough a HaZ verex for a sufficienly heavy scalar H. Unlike he graph on he righ-hand side i has no counerpar in he spin-0 DMF simplified models. As firs emphasised in [20] he appearance of he conribuion wih virual H and a exchange no only enhances he mono-z cross secion compared o he spin-0 DMF models, bu also leads o quie differen kinemaics in Z + χ χ producion. In fac, for masses M H > M a + M Z he prediced E T,miss specrum urns ou o be peaked a E max T,miss λ1/2 (M H, M a, M Z ) 2M H, (5.2) 18

20 g a g a g g Figure 8. Two possible diagrams ha give rise o a + E T,miss signal. Graphs wih boh an exchange of an a and A conribue in he THDM plus pseudoscalar exensions bu only he former are displayed. where he wo-body phase-space funcion λ(m 1, m 2, m 3 ) has been defined in (4.10). Denoing he lower experimenal requiremen on E T,miss in a given mono-z search by ET,miss cu he laer resul can be used o derive a simple bound on M H for which a significan fracion of he oal cross secion will pass he cu. We obain he inequaliy M H M a + MZ 2 + ( ET,miss) cu 2. (5.3) Given ha in he laes mono-z analyses [28 30] selecion cus of ET,miss cu 100 GeV are imposed i follows ha he scalar H has o have a mass of M H 500 GeV if one wans o be sensiive o pseudoscalars a wih masses up o he hreshold M a 350 GeV. Our deailed Mone Carlo (MC) simulaions of he Z + E T,miss signal in Secion 6.4 however reveals ha he above kinemaical argumen alone is insufficien o undersand he shape of he mono-z exclusion in he M a an β plane in all insances. The reason for his is wofold. Firs, in cases where sin θ is small H az is ofen no he dominan H decay mode and as a resul he Z + E T,miss measuremens lose already sensiiviy for masses M a below he bound implied by he esimae (5.3). Second, Z +χ χ producion in gg az and gg AZ is also possible hrough box diagrams, and he inerference beween riangle and box graphs urns ou o be very relevan in models ha have a ligh scalar H or pseudoscalar A wih a mass below he hreshold. We add ha for an β > O(10) also resonan mono-z producion via b b az and b b AZ can be relevan in models of ype II and IV. In he conex of he pure THDM such effecs have been sudied for insance in [95]. 5.4 Mono-Higgs channel In cerain regions of parameer space anoher possible smoking gun signaure of he pseudoscalar exensions of he THDM urns ou o be mono-higgs producion. As illusraed in Figure 10 his signal can arise from wo differen ypes of one-loop diagrams. For M A > M a + M h he riangle graph wih an Aah verex depiced on he lef-hand side allows for resonan mono-higgs producion and hus dominaes over he conribuion of he box diagram displayed on he righ. In consequence he mono-higgs producion cross secions in he THDM plus pseudoscalar exensions can exceed by far he small spin-0 DMF model raes for he h + E T,miss signal [88]. 19

21 g g H a Z g Z g a Figure 9. Represenaive Feynman diagrams ha lead o a Z + E T,miss signal in he pseudoscalar exensions of he THDM. In he case of riangle diagram (lef) only he shown graph conribues, while in he case of he box diagram (righ) insead of an a also an A exchange is possible. Like in he case of he mono-z signal he presence of riangle diagrams wih a rilinear scalar coupling also leads o a peak in he E T,miss disribuion of h + χ χ producion if he inermediae heavy pseudoscalar A can be resonanly produced. The peak posiion in he mono-higgs case is obained from [20] E max T,miss λ1/2 (M A, M a, M h ) 2M A. (5.4) I follows ha in order for evens o pass he E T,miss cu necessary for a background suppression in mono-higgs searches, he relaion M A M a + Mh 2 + ( ET,miss) cu 2, (5.5) has o be fulfilled. A lesson o learn from (5.5) is ha mono-higgs searches in he h b b channel [31, 32] are less suied o consrain he parameer space of our simplified model han hose ha focus on h γγ [33, 34], because he minimal E T,miss requiremens in he former analyses are always sricer han hose in he laer. To give a relevan numerical example le us consider ET,miss cu 100 GeV, which represens a ypical E T,miss cu imposed in he mos recen h + χ χ (h γγ) searches. From (5.5) one sees ha in such a case mono-higgs analyses are very sensiive o masses up o M a 330 GeV for M A 500 GeV. Like in he mono-z case he above kinemaical argumen however allows only for a qualiaive undersanding of he numerical resuls for he pp h+χ χ (h γγ) exclusions, since inerference effecs can be imporan in scenarios wih a pseudoscalar A of mass M A < 2m. Noice ha if M a > M A + M h he role of A and a is inerchanged and he h + E T,miss signal can receive large correcions from resonan a exchanges, as we will see explicily in Secion 6.4. Finally in ype II and IV models resonan mono-higgs producion from b b iniial saes can also be imporan if an β is sufficienly large. 5.5 Mono-W channel The las E T,miss signal ha we consider is he mono-w channel [35, 36]. Two represenaive Feynman graphs ha lead o a resonan W +E T,miss signaure in he pseudoscalar exension of he THDM are shown in Figure 11. These diagrams describe he single producion of a 20

22 g g A a h g h g a Figure 10. Sample diagrams in he THDM wih an exra pseudoscalar ha induce a h + E T,miss signal in he alignmen/decoupling limi. Graphs in which he role of a and A is inerchanged can also provide a relevan conribuion. charged Higgs H ± via he annihilaion of ligh quarks followed by H ± aw ± (a χ χ). One way o assess he prospecs for deecing a mono-w signaure consiss in comparing he producion cross secions of H ± o ha of H and A. Using for insance an β = 1, we find σ (pp H + ) 1.0 fb for M H ± = 500 GeV and σ (pp H + ) 0.2 fb for M H ± = 750 GeV a he 13 TeV LHC. The corresponding cross secions in he case of he heavy neural spin-0 resonances read σ (pp H) 1.4 pb and σ (pp A) 3.1 pb and σ (pp H) 0.2 pb and σ (pp A) 0.3 pb, respecively. These numbers srongly sugges ha an observaion of a mono-w signal is compared o ha of a mono-z or mono-higgs signaure much less probable. We hus do no consider he W + E T,miss channel any furher. Le us finally add ha besides a simple mono-w signaure also W + E T,miss and W b + E T,miss signals can appear in he DM model inroduced in Secion 2. For he relevan charged Higgs producion cross secions we find a 13 TeV he resuls σ ( g b H + ) 0.17 pb (σ ( g b H + ) 0.04 pb) and σ (gg H + b) 0.10 pb (σ (gg H + b) 0.02 pb) using an β = 1 and M H ± = 500 GeV (M H ± = 750 GeV). Given he small H ± producion cross secion in gb and gg fusion, we expec ha searches for a W +E T,miss or a W b+e T,miss signal will in pracice provide no relevan consrain in he small an β regime. 6 Numerical resuls The numerical resuls of our mono-x analyses are presened in his secion. Afer a brief descripion of he signal generaion and he background esimaes, we firs sudy he impac of inerference effecs beween he a and A conribuions o he j +χ χ and +χ χ channels. Then he consrains on he parameer space of he THDM plus pseudoscalar exensions are derived for several well-moivaed benchmark scenarios. In he case of he mono-z and mono-higgs searches we also discuss he LHC Run II reach in some deail. 6.1 Signal generaion The saring poin of our MC simulaions is a UFO implemenaion [96] of he simplified model as described in Secion 2. This implemenaion has been obained by means of he FeynRules 2 [97] and NLOCT [98] packages. The generaion of he j + E T,miss, 21

23 u H + W + s H W d a c a Figure 11. Examples of diagrams ha lead o a W + E T,miss signaure hrough he exchange of a charged Higgs H ± and a ligher pseudoscalar a in he THDM plus pseudoscalar exension. Z + E T,miss (Z l + l ) and h + E T,miss (h γγ) signal samples is performed a leading order (LO) wih MadGraph5_aMC@NLO [99] using PYTHIA 8.2 [100] for showering and NNPDF2.3 [101] as paron disribuion funcions. The whole MC chain is seered wih CheckMATE 2 [102] which iself employs FasJe [103] o reconsruc hadronic jes and Delphes 3 [104] as a fas-deecor simulaion. The resuls of he CheckMATE 2 analyses have been validaed agains MadAnalysis 5 [105, 106]. The selecion requiremens imposed in our analyses resemble hose used in he recen LHC mono-je [22], mono-z [28] and mono-higgs [34] search, respecively. For wha concerns our + E T,miss ( lbν) recas we rely on he resuls of he sensiiviy sudy [94]. In his analysis he DM signal has been simulaed a nex-o-leading order (NLO) wih MadGraph5_aMC@NLO and PYTHIA 8.2 using a FxFx NLO je maching prescripion [107] and he final-sae op quarks have been decayed wih MadSpin [108]. 6.2 Background esimaes For he j+e T,miss, Z +E T,miss (Z l + l ) and h+e T,miss (h γγ) recass our background esimaes rely on he background predicions obained in he 13 TeV LHC analyses [22], [28] and [34], respecively. The given background numbers correspond o 3.2 fb 1, 13.3 fb 1, 2.3 fb 1 and we exrapolae hem o 40 fb 1 of inegraed luminosiy o be able o assess he near-erm reach of he differen mono-x channels. Our exrapolaions assume ha while he relaive sysemaic uncerainies remain he same, he relaive saisical errors scale as 1/ L wih luminosiy L. Depending on he signal region he relaive sysemaic uncerainies amoun o around 4% o 9% in he case of he mono-je search, abou 7% for he mono-z analysis and approximaely 20% for he mono-higgs channel. Since he j + E T,miss search is already sysemaics limied a 40 fb 1 is consraining power will depend sensiively on he assumpion abou he fuure sysemaic uncerainy on he associaed SM background. This should be kep in mind when comparing he differen exclusions presened below, because a beer undersanding of he backgrounds can have a visible impac on he obained resuls. Since he + E T,miss ( lbν) search will sill be saisically limied for 40 fb 1, we base our forecas in his case on a daa se of 300 fb 1 assuming ha he relevan SM background is known o 20%. In he mono-z and mono-higgs cases we will presen below, besides 40 fb 1 projecions, resuls for 100 fb 1 and 300 fb 1 of daa. From hese resuls one can assess if he exising Z + E T,miss and h + E T,miss search sraegies will a some poin become sysemaics limied in LHC Run II. 22