2HDM portal for singlet-doublet dark matter

We present an extensive analysis of a model in which the (Majorana) Dark Matter candidate is a mixture between a SU(2) singlet and two SU(2) doublets. This kind of setup takes the name of singlet-doublet model. We will investigate in detail an extension of this model in which the Dark Matter sector interactions with a 2-doublet Higgs sector enforcing the complementarity between Dark Matter phenomenology and searches of extra Higgs bosons.


Introduction
The WIMP paradigm is a compelling solution of the Dark Matter (DM) problem. It relates the achievement of the DM relic density, as measured with unprecedented precision by the PLANCK collaboration [1], to a specific range of values of the thermally averaged pair annihilation cross-section of the DM. An implication of this setup is that the DM should possess sizeable interactions with the Standard Model (SM) particles, making possible a detection at present experimental facilities.
A fermionic DM candidate, if it is a SM singlet, can couple, in pairs, with the Higgs boson only through D > 4 operators. 1 The so called singlet-doublet models [17][18][19][20][21] 2 overcome this problem by enlarging the spectrum of BSM states by two SU (2) doublets, so that the DM is a mixture of their neutral components as well as of the singlet originally introduced. This has, however, the consequence that the DM can interact, in pairs, also with the Z boson, as well as, through the charged component of the extra doublets, with the W boson.
As recently reviewed in [23] (see also e.g. [24,25]), DM interactions mediated exclusively by the Higgs and the Z bosons are disfavored by DM Direct Detection (DD), expecially in the case that the DM is a dirac fermion [19,26,27].
In this work we will investigate in detail whether this problem can be encompassed by extending, with a second dou-blet, the Higgs sector of the theory. A similar investigation has been already presented in [28], but with focus only on the possibility of a light pseudoscalar boson. While including this scenario in our discussion, we will, however, investigate the parameter space of the theory from a more general perspective. We will pinpoint, furthermore, the complementarity with constraints from searches at collider and in low energy processes of extra Higgs bosons.
The paper is organized as follows. We will first introduce, in Sect. 2, our model setup. Section 3 will be devoted to a brief review of the two-doublet extension of the Higgs sector and to the discussion of the theoretical and experimental limits which can impact the viable parameter space for DM. The most salient features of DM phenomenology will be then discussed in Sect. 4. We will finally present and discuss our findings in Sect. 5.

The model
2.1 2HDM and coupling to the SM We will adopt, for our study, a 2HDM model described by the the following potential: where two doublets are defined by: further assuming that all the couplings above are real. Equation 1 further assumes the absence of two additional couplings λ 6 and λ 7 . This is for example obtained by assuming that the Higgs doublets transform as H 1 → H 1 and H 2 → −H 2 under a discrete Z 2 symmetry [29] (see also [30][31][32]). We also introduce, as usual, the β angle defined as v 2 /v 1 = tan β. Imposing CP-conservation in the scalar sector, the spectrum of physical states is constituted by two CP even neutral states, h, identified with the 125 GeV Higgs, and H , the CP-odd Higgs A and finally the charged Higgs H ± . The transition from the interaction basis (H 1 , H 2 ) T to the mass basis (h, H, A, H ± ) depends on two mixing angles, α and β according the following relations: cos β − sin β sin β cos β G + H + , with G 0 , G + being Goldstone bosons, eaten, after EW symmetry breaking, as longitudinal degrees of freedom, by the SM gauge bosons. The quartic couplings of the scalar potential (1) can be expressed as function of the masses of the physical states as: where M ≡ m 12 /(s β c β ).
The SM fermions cannot couple freely with both Higgs doublets since, otherwise, FCNCs would be induced at tree level. Four specific configurations, labelled Type-I, Type-II, Lepton specific and Flipped, avoid this eventuality. In the physical basis for the scalar sector the interaction Lagrangian between the Higgses and the SM fermions reads: where v = v 2 1 + v 2 Table 1. Constraints from 125 GeV Higgs signal strength limit the values of α and β. These bounds will be discussed in more detail below. We just mention that one can automatically comply with them by going to the so-called "alignment" limit, i.e. β − α = π 2 , which makes automatically the couplings of the h state SM like, i.e. ξ u,d,l h = 1, and the other ξ Table 1 Couplings of the Higgses to the SM fermions as a function of the angles α and β and in the alignment limit where (β − α) → π/2 Type-I Type-II Lepton-specific Flipped parameters only dependent on tan β. A further implication of the alignment limit is that the coupling of the CP-even state H with the W and Z bosons is null (the couplings of the A boson with Z Z and W W are null as long CP is conserved in the Higgs sector). In our study we will keep as free as possible the parameters of the Higgs sector. We will then do not strictly impose the alignment limit but rather keep α and β as free parameters and impose on them the relevant constraints.

Coupling to the DM
In the scenario under investigation, the DM arises from the mixture of a SM singlet N and the neutral components of two (Weyl) SU (2) doublets L L and L R with hypercharges equal, respectively, to −1 and +1, defined as: The new fermions are coupled with two Higgs doublets H 1 and H 2 according the following Lagrangian: Notice that, given their quantum numbers, the new fermions could be coupled, through the Higgs, also with SM leptons and, hence, mix with them after EW symmetry breaking. This would imply, in particular, that the DM is not stable. To avoid this possibility we assume the existence of a discrete Z 2 symmetry, which we label as Z DM 2 , under which the new and the SM fermions are, respectively, odd and even.
Similarly to SM fermions, it is appropriate to avoid to couple these new states freely with both Higgs doublets. Several different configurations are nevertheless still possible (see for example [28]). For definiteness we will then focus, along this work, on two specific configurations (we will comment on other possibilities in Appendix B): y L 1 = y 1 , y L 2 = y R 2 = 0, y R 1 = y 2 and :y L 1 = y 1 , y L 2 = y R 1 = 0, y R 2 = y 2 . For the first configuration we will further assume that the SM fermions are coupled as in the standard Type-I 2HDM and globally label as "Type-I" the model defined in this way. The second configuration for the new fermions is, instead, accompanied by coupling of the SM fermions with H 1 , H 2 as in the Type-II 2HDM and we will refer this scenario as "Type-II' model.
After EW symmetry breaking, mixing between N and the neutral components of L L and L R occurs, so that the physical spectrum of the new fermions is represented by three Majorana fermions ψ i=1,3 defined by: and one charged dirac fermion ψ ± with mass m ψ ± ≈ M L . The matrix U diagonalizes a mass matrix of the form: where, for definiteness we have considered the Type-II model. Similarly to [18,20,28], we will adopt, in spite of y 1 , y 2 , the free parameters y, θ defined by: In the mass basis, the relevant interaction Lagrangian for DM phenomenology reads: where: Notice, in particular, that, as expected from its Majorana nature, the vectorial coupling of the DM with the Z boson is null.

Bounds on the potential
The quartic couplings λ i=1,5 should comply with a series of constraints coming from the unitarity and boundness from below of the scalar potential as well as perturbativity (see for example. [33,34] for more detailed discussions). These bounds, can be translated, through Eq. 4, into bounds on the masses of the new Higgs bosons as function of the angles α and β. For completeness we list below the main constraints: -Scalar potential bounded from below: -tree level s-wave unitarity: where: -global minimum in (v 1 , v 2 ) for the scalar potential [35]: where the mass parameters m 11 , m 22 , m 12 should satisfy:

EWPT
The presence of extra Higgs bosons affects the values of the Electroweak Precision Observables (EWPO), possibly making them to deviate from the SM expectations. One option to alleviate these tensions would be represented by imposing specific relations for the masses of the new states; for example deviations of the T parameter can be avoided imposing a custodial symmetry [36,37]. As already pointed, for our DM analysis, we will nevertheless try to keep the parameters of the Higgs sector as free as possible. To identify the viable paramter space we have then computed the S, T, U parameters using the functions provided, for example in [38][39][40][41], as function of the masses m H , m A , m H ± of the new Higgs bosons as well as function of α and β (we remind we are not imposing the alignment limit) and determining the excluded model configurations through the following χ 2 [42][43][44]: with X = (S, T, U ) while X SM , σ and V represent, respectively, the SM expectations of the Electroweak Precision Observables, the corresponding Standard deviation and covariant matrix. Their updated values have been provided here [45] and are also reported here for convenience: We have imposed to each model point to do not induce a deviation, for the EWPO, beyond 3σ from their best fit values. The EWPO can be as well affected by the fermionic sector of the theory [39,40,[46][47][48]. The corresponding contributions to the STU parameters depend on masses and on the Yukawa couplings for the new fermions. In particular in the limit in which the latter parameters vanish also the contributions to EWPO from the new fermions become negligible.
On general grounds one should include in Eq. 23 also the contribution from the new fermionic sector to the EWPO. However, in agreement with the findings of [20] (for the computation of the EWPO we have adapted the expressions provided in [48]), the new fermionic sector mostly affects the T parameter. Given this, we have imposed a priori that the deviation ΔT of the T parameter induced by the new fermions is neglible and then used Eq. 23 to constrain the parameters of the scalar sector. 3 3 One could have considered, as done for example in [27], the possibility of cancellations between the contributions to the EWPO from the new The deviation ΔT induced by the DM sector scales with the combination y 2 1 v 2 a − y 2 2 v 2 b 2 , with a, b set appropriately to 1, 2, depending on whether the Type-I or the Type-II models are considered, and then vanishes in the limit tan θ → 1 (tan β → 1 would be also required for the Type-II model).
We have verified that ΔT remains sensitively below the experimental exclusion unless y > 1 and | tan θ | > 10 (in the case of the Type-II model the contribution to the T parameter increases also with tan β. As will be clarified below, moderate-high values of tan β are also disfavoured by LHC constraints as well as by DM phenomenology and, hence, have not been considered). As already mentioned, for the study of the DM phenomenology we have chosen assignations of the parameters, in the new fermionic sector, automatically complying with the bounds from EWPT.

Collider searches of the new Higgs bosons
H and A bosons can be resonantly produced at colliders through gluon fusion, 4 [73,74]. Among these, the most effective bounds come from searches of τ + τ − pairs, which exclude, for the Type-II model, moderate-high values of tan β, above 10. Searches of diboson final states and of the process A → Zh(H ) (as well as the process with inverse mass ordering) can, in addition, constrain deviations from the alignment limit.
The strongest limits for the latter come, however, from the Higgs signal strenght. For our analysis we have adopted the limit determined in [27] combining the constraints on the Higgs signal strenght into γ γ , W + W − , Z Z,bb and τ + τ − from the LHC run I [75,76] (see also [77]). 5 The excluded regions in the bidimensional plane (cos(β − α), tan β) are sketched in Fig. 1. As evident sizable deviations from the alignment limit are allowed in the case of the Type-I 2HDM. This is not the case, on the contrary, for the Type-II 2HDM, ad exception of a very narrow window corresponding to the so called "wrong-sign" Yukawa regime [81][82][83], i.e. the case in which the couplings of the state h with the down type (Footnote 3 continued) scalar and fermionic sectors. This would allow for higher values of the yukawa couplings, provided very strong constraints in the assignation of m A , m H , m H ± . We have not considered this scenario here. 4 For 2HDM realizations with ξ d H,A enhanced with tan β alsobb fusion can play a relevant role [27]. The Type-II model considered in this work features this property. However, as discussed in the text, limits from searches of resonances decaying into τ + τ − final states put strong limits on the size of tan β so that one can realiably assume gluon fusion as the most relevant production process. 5 Further stronger constraints would be obtained by considering Run II data [78][79][80]. These have not a relevant impact in our analysis. Similarly to [28], we will include in our analysis also the case in which m A < m h . In such a case one should take into account possible limits on the process h → A A. This can be generically constrained through the Higgs signal strenght (i.e. one generically imposes that the branching fraction of this decay does not exceed the allowed value of the invisible branching fraction of the Higgs. This bound is effective in the whole range of masses for which the decay h → A A is kinematically allowed) as well as through dedicated searches [84] (limits are effective only for some range of masses). More contrived are instead the prospects for signals not related to Higgs decays, see anyway [85,86].
Concerning the charged Higgs boson we have first of all a limit m H ± 80 GeV from LEP [87]. Moving to LHC constraints, these come, for m H ± < m t , from searches of top decays t → H ± b with H ± [88][89][90][91] decaying into τ ν τ or cs. The corresponding limits have been reformulated in [92] for different realizations of the 2HDM. Values of the masses of the charged Higgs for which the decay is kinematically allowed are excluded for tan β 10 in the Type-I 2HDM and irrespective of the value of tan β for the Type-II scenario. For m H ± > m t searches rely on direct production of the charged Higgs in association with a top and a bottom quark, followed by the decay of the former in τ ν [90,[93][94][95] or tb [94]. Associated limits are not competitive with the others discussed in here and will be then neglected. The charged Higgs feels, indirectly, also limits from searches of the neutral Higgs bosons since the conditions on the quartic couplings 18-19 impose relations between the masses of the new Higgs bosons (see e.g [92]).

Limits from flavour
While it is possible to avoid that the couplings of the extra Higgs bosons with SM fermions induce FCNC at the tree level, they can impact flavour violating transitions at the loop level. The strongest limits come from processes associated to b → s transitions. Their rates are mostly sensitive to m H ± and tan β. Experimental limits are formulated in terms of these parameters. The most stringent come from the B → X s γ process [96] and are particularly severe in the case of Type-II model, excluding m H ± < 570 GeV irrespective of tan β [97]. Much better is, instead, the situation of the Type-I model where we have just the approximate lower bound tan β 2. A similar exclusion, for both Type-I and Type-II models is also provided by the processes B s → μ + μ − and

Scanning the parameter space
The bounds discussed in the previous subsection have been applied to a parameter scan over the following ranges 6 As already said, we will include in our analysis a very light pseudoscalar A. For this reason we have considered a minimal value of 20 GeV for its mass in the scan. More precisely the scan consisted into random extracting, with flat priors, the model parameters within the ranges reported and keeping the sets of model points satisfying the following constraints: -the constraints 18-21 on the quartic couplings λ 1,5 determined, according 4, as function of m H , m A , m H ± , α, β; -the bounds from Higgs signal strengths, imposing that the values of α, β for each model point lie in the white regions of 1; -bounds from EWPT, according the procedure depicted in sec. III B, and the flavour constraints, illustrated in sec. III D. -compatibility with null results from searches of new Higgs bosons at the LHC. To this purpose we have determined the production cross-sections of the Higgs bosons through the same procedure already employed in [27].
We have then adopted the cross-section provided by the LHC Higgs Cross Section Working Group [98], which have been determined through the package SusHi [99].
In the case of the Type-II model it has been convenient, to our purposes, to adopt the cross-sections computed for the hMSSM [100,101], given the similarity of setups. For what concerns the Type-I model, as discussed in [27], given the fact that the gluon fusion is the dominant production channel and that all the couplings of the H/A bosons with SM fermions have the same dependence on tan β, the production cross sections could be simply computed by rescaling by a 1/ tan 2 β factor the hMSSM cross sections computed for tan β = 1.
This procedure has been performed for both Type-I and Type-II 2HDM individually. For each model a sample of 10,0000 viable solutions has been collected. Among the viable model points we have identified some benchmark assignations for m H , m A , m H ± , α, β which will be used for the DM analysis. We have investigated, in particular, the possibility of obtaining a light mass of the pseudoscalar A, possibly lighter than the SM like Higgs h, while keeping the masses of the other Higgs states above few hundreds GeV. This possibility is easily realized in the Type-I model provided that m A 60 GeV, in order to kinematically forbid the h → A A decay. On the contrary, no viable model points, with a light pseudoscalar, have been found in the Type-II model. This is due to the strong bound on the mass of the charged Higgs from b → s transitions (see Sec. III D) in combination with the strong limits on cos(β − α) from the Higgs signal strengths, which disfavour a sizeable hierarchy between the masses of the different Higgs states. The combination of these two constraints limited the range of viable pseudoscalar masses to m A 500 GeV. 7 7 As evidenced by the right panel of Fig. 1, and already discussed in the text, a sizeable deviation from the alignment limit in the Type-II model is actually still allowed in a narrow strip for cos(β − α) > 0 region, dubbed wrong sign Yukawa region. In this region the hierarchy between the mass of the pseudoscalar and the charged Higgs can be augmented. We have found that in this region the lower bound on m A , We then illustrate in the next section the main constraints coming from DM phenomenology.

Relic density
According to the WIMP paradigm, the DM has sizable enough interactions with the SM particles to be in thermal equilibrium in the Early stages of the history of the Universe. At a later stage the interaction rate of the DM fell below the Hubble expansion rate causing the freeze-out of the DM at temperatures of the order of 1 20 − 1 30 the DM mass. Assuming standard cosmological history, the DM relic density, Ω ψ 1 h 2 , is determined by a single particle physics input, i.e. the DM thermally averaged pair annihilation cross-section. The relation between these two quantities is given by [102]: where T 0 and T f represent, respectively, the present time and freeze-out temperatures while g 1/2 * is a function of the relavistic degrees of freedom at the temperature T [102]. σ eff v is the effective annihilation cross-section [103]: including coannihilation effects from the additional neutral and charged states belonging to the Dark Matter sector. Coannihilation effects are expected to be important in the case M L M N , corresponding to a DM with sizeable or even dominant doublet component, implying that at least the charged fermion ψ ± is very close in mass to it.
We remind that the precise experimental determination Ωh 2 ≈ 0.12 [1] is matched by a value of σ eff v of the order of 10 −26 cm 3 s −1 .
In the model considered in this work, a huge variety of processes contributes to σ eff v . For what regards DM pair annihilations, the most commonly considered are the annihilation into SM fermion pairs, W W , Z Z and Zh, originated by s-channel exchange of the h, H, A bosons as well as, in the case of annihilation into gauge boson pairs, tchannel exchange of the DM and the other new fermions. In this work we will put particular attention also to the case (Footnote 7 continued) indirectly induced by the one on m H ± from Br(B → X s γ ), reduces to m A 300 GeV. All the possible DM annihilation channels have been included in our numerical study, performed through the package micrOMEGAs [104]. The DM relic density has been numerically determined with great accuracy by computing all the relevant annihilation processes, including coannihilations, adopting the Feynman rules, correspoding to the model under investigation, determined by the package Feynrules [105]. For a better insight, we report, nevertheless, in appendix A, the expressions of some phenomenologically interesting channels, by making use of the velocity expansion, σ v = a + bv 2 ψ . We have reported in Fig. 2 the isocontours corresponding to the correct DM relic density, Ωh 2 ≈ 0.12 [1] Fig. 2. Far from the "poles" and the pure doublet regime, the isocontours of the correct relic density are determined by the size of the couplings of the DM with the Higgs bosons. For y = 1 the correct relic density can be obtained for a singlet-like DM thanks, only in the Type-II model, to the tan β enhancement of the annihilation channels intobb and τ + τ − , and, when kinematically accessible, to the annihilation channels with new Higgses in the final state. Once the value of y decreases, the isocontours tend to move towards the "well tempered" [106,107] regime M N ∼ M L . We notice in particular that, for y = 0.2, the lines of the correct relic density, for all the four benchmarks, lie almost exactly in the A similar outcome would be also expected for lower values of y. For this reason we can focus, without loss of generality, to values y ≥ 0.2. As a remark we notice that the Type-II benchmark reported in the top right panel of Fig. 2 is in tension with the constraints on the extensions of the Higgs sector discussed in the previous section. It should be then interpreted just as illustration.

Direct detection
In the considered scenario, the DM features both spin independent (SI) and spin dependent (SD) interactions with nuclei. The former are induced, at three level, by t-channel exchange of the CP-even h and H states. The corresponding cross-section reads (for definiteness we report the case of scattering on protons): where f p q are nucleon form factors (notice that the form factors corresponding to heavy quark are expressed in terms of the gluon form factor [108,109] The strongest limits on spin independent interactions come, at the moment, from the XENON1T experiment [110] (slightly weaker limits have been provided, as well, by the LUX [111] and PANDAX [112]). In our analysis we will consider, in addition, the sensitivity prospectes of the LZ [113], XENONnT [114] (given their similarity, we will represent the sensitivity prospects of these experiments as a single line) and DARWIN [115].
Spin dependent interactions are originated, instead, by interactions of the DM with the Z boson, giving the following cross-section: where Δ p u,d,s are again suitable structure functions 9 while g A u and g A d represent the axial coupling of the Z boson with, respectively, up-type and down-type quarks.
For our phenomenological analysis we have adopted the most recent limits provided by the LUX collaboration [116] (slightly weaker limits have been provided by PANDAX.) As well known, Direct Detection constraints are mostly associated to SI interactions, because of the coherent enhanchment (A 2 ) occuring when these are evaluated at the nuclear, rather than nucleon, level. However the different couplings of the DM with two mediators, the h and H , can induce so called "blind spots" [28,[117][118][119], due to a possible destructive interference between diagrams with h and H exchange or to cancellations of the couplings of the DM with one or both the CP-even Higgs states. An analytical insight can be obtained by inspecting expressions of the couplings of the DM with the h, H states as function of the DM mass and of the (M N , M L , y, θ, β) parameters [28] (for simplicity we consider the alignment limit). These are for the Type-I model: while for the Type-II model: cos 2β cos 2θ ) × m ψ 1 sin 2β cos 2θ + M L cos 2β sin 2θ (31) As evident, in the Type-I model, both the DM couplings with the CP-even bosons become zero in case that m ψ 1 + M L sin 2θ = 0 [18,20]. This is actually a generic blind spot occurring in the case that all the new fermions couple to a same Higgs doublet. The case of the Type-II model is less trivial and the blind spot is obtained either by the condition for destructive interference: or in the case that the coupling of the DM with the SMlike Higgs is null, so m ψ 1 (1 + cos 2β cos 2θ ) + M L sin 2β sin 2θ = 0, while the other CP-even Higgs is heavy enough so that the corresponding contribution to the DD crosssection is not in tension with experimental limits. An illustration of the occurrence of the blind spots, as function of tan θ and some example assignations of the parameters, is shown in Fig. 3 for the Type-I (upper panel) and Type-II (lower panel) model. The plot represents indeed the ratio between the computed scattering cross-section and the limit on the cross-section, for the value of the DM mass corresponding to the considered parameters, as provided by XENON1T. The blind spot corresponds to a sharp decrease, from values eventually greater than 1, implying an experimentally excluded configuration, of this ratio to values much below 1, corresponding to an experimentally viable configuration. In the case of the Type-I model, in agreement with the previous discussion, the position of the blind spot is poorly sensitive to the assignation of m H and y. Since we have considered a value of M N sensitively lower than M L ,  Fig. 3, as well as throughout the rest of the paper, we have also included one-loop corrections to SI cross-section arising from interaction of the DM with the Z and W bosons [120][121][122] and, eventually, with the pseudoscalar boson A [123][124][125] if this is light enough. As evidenced by the figure, even in presence of these corrections, it is possible to achieve a cross-section well below experimental limits for specific values of the angle θ .

Indirect detection
As evidenced by the expressions provided in Appendix A, the DM annihilation channels into SM fermions and gauge bosons pairs feature s-wave, i.e. velocity independent, annihilation cross-sections. Thermal DM production can be thus tested, in our framework, also through Indirect Detection (ID). The most prominent signals come from γ -rays originating mainly from annihilations into W W , Z Z,tt,bb and τ + τ − . Particularly interesting would be, in this context, the scenario of a light pseudoscalar since it would allow for a fit of the γ -ray Galactic Center (GC) excess [28,126,127]. DM interpretations of γ -ray signals are, nevertheless, challenged by the exclusion limits from absence of evidences in Dwarf Spheroidal Galaxies (DSph) [128] as well as, since recently, searches in the Milky-Way Halo away from the GC [129] (these exclude, in particular, thebb interpretation of the GC excess). In this work we will not attempt to provide a DM interpretation of the GC and focus, more conservatively, on the constraints on the DM annihilation cross-section.
At the moment Indirect Detection constraints can probe thermal DM production up to DM masses of around 100 GeV and are rarely competitive with respect to Direct Detection constraints. In order to simplify the presentation of our results we will report Indirect Detection limits only when they are effectively complementary to other experimental searches while and omit them in the other cases.

Invisible decays of the Higgs and of the Z
In the setup under consideration the DM is coupled, in pairs, both to the Higgs and to the Z boson. In the case it is lighter then m h /2, m Z /2, an invisible decay channel for the latter becomes accessible. This possibility is however experimentaly disfavored [130]. We have hence imposed in our analysis that, when the processes are kinematically allowed, the invisible branching ratio of the Higgs fullfills the upper bound

LHC searches of the new fermions
The new fermionic sector featured by our model might be probed also by collider searches. A detailed study goes, however, beyond the scopes of this paper. We will then illustrate qualitatively the most relevant searches and their potential impact on the parameter space of the model.
In the scenario under consideration at least some of new fermions feature sizeable couplings with the SM gauge bosons. They can be then produced, in proton collisions at LHC, through Drell-Yann processes of the type: with the charged fermions subsequently decaying into a W (either on or off shell) and a neutral fermion (typically the lightest one, i.e. the DM candidate) while the heavy neutral fermions features different possible decay channels, i.e. into a lighter one and a Z, or possibly, Higgs boson (in both cases the particles can be off-shell), or the charged fermion ψ ± and a W boson. In the case of production processes involving neutral fermions, the different possible signals depend on the size of the elements of the mixing matrix U , determining the size of the hypercharge and doublet components of the ψ i states.
For what concerns detection prospects, the cleanest signatures are provided by production processes of the heavy fermions, other than the DM ψ 1 , with subsequent decays involving W/Z bosons, leading to events with missing energy accompagned by 2-4 charged leptons (see e.g. [20] for a classification of the possible processes). Among these processes the strongest constraints apply to pp → ψ i=2,3 ψ ± → W ± Z ψ 1 ψ 1 , which leads to a 3 lepton plus missing energy signature. The corresponding limits from CMS and ATLAS searches (see [131,132] for the most recent) should be appropriately recasted in order to be applied for the scenario under consideration. Such recasting has been performed, for example, in [133] for the √ s = 8 TeV data sets [134,135]. It has been found there that these searches constrain the regions m ψ 2 , m ψ 3 , m ψ ± 270 GeV and m ψ 1 75 GeV. Notice that these constraints mostly apply to the case M N < M L . In the opposite case the lightest neutral fermions are mostly doublet like and are very close in mass with the charged state ψ ± . This kind of configuration is more complicated to probe at the collider since it would correspond to the production of long-lived particles, leading to displaced vertices or being eventually detector stable.
The new fermions might also be produced in the decays of resonantly produced new Higgs bosons. In such a case it would be possible to the test directly the production of DM pairs by looking for events with missing energy accompagned by initial state radiation, i.e. monojet. Recasting the corresponding limits [136,137], goes beyond the scopes of this work and will be left to future study (see e.g. [138] for an example of this recasting on a simplified model).
We finally remind that the mass of the charged new fermions m ψ ± M L is also constrained by LEP [139]. For this reason in all the plots shown below, only regions with M L 100 GeV are visualized.

Sommerfeld enhancement for relic density and Indirect Detection
A complete treatment of the DM phenomenology would require the inclusion of the Sommerfeld enhancement for the DM annihilation cross-section. The Sommerfeld enhancement consists into an enhancement of the DM annihilation cross-section at low relative velocity, due to the presence of an attractive long-range potential, typically induced by the multiple exchange of light (with respect to the DM) mediators. In the case in which the DM is, at least partially, charged under the EW interactions, a sizeable Sommerfeld enhancement might be produced by loop effects from exchange of electroweak bosons, when the DM mass is of the order or above the TeV scale. The associated corrections to the DM relic density have been determined in the limit of a pure SU (2) nplet DM in [140][141][142] and subsequently revised and corrected in [143,144]. In case of doublet like DM these corrections have been found to be rather small, hence corresponding to a slight (non-neglible though) increase of the value of the mass (in case determined mostly by M L ), with respect to the one quoted in Sect. 4.1, corresponding to the correct relic density. The computation of the Sommerfeld enhancement in the case of mixed (singled-doublet in our scenario) DM is typically extremely complicated as consequence of its strong dependence on the mass splitting between the DM and the lightest charged and neutral states of the new fermionic spectrum. This kind of computation has been tackled mostly in SUSY inspired scenarios, see e.g. [145][146][147]. More recently more general scenarios have been considered though [144]. The detailed determination of the Sommerfeld correction in the model under study is not in the scopes of this paper. Given the already limited impact on the relic density in the pure SU (2) limit we expect anyway a marginal variation to our results from the inclusion of Sommerfeld corrections. Besides the gauge bosons, the scenario under study features two potential light force carriers, being the SM-like Higgs boson h and the light pseudoscalar A. The Sommerfeld enhanchement associated to h has been considered in the case of the Higgs portal model [148] while studies for a generic pseudoscalar mediator can be found e.g. in [149,150]. Similarly to the case of the gauge bosons, Sommerfeld enhancement induced by exchange of h is expected to be relevant for DM masses well above the TeV scale, marginally considered in our study. Less trivial would be instead the case of the pseudoscalar. As evidenced by Fig. 6, nevertheless, we have considered the case of a light pseudoscalar (we remind that in our setup it cannot be lighter than around 60 GeV in order to comply with collider bounds) in relation to comparatively light DM masses, not above 300 GeV; consequently, Sommerfeld enhancement can be safely neglected.
Being a low velocity effect, Sommerfeld enhancement can impact Indirect Detection to a sensitively larger extent with respect to relic density. Weak constrains are, however, expected in the case of annihilation into scalar and pseudoscalar final states, being them p-wave processes, and, hence, they have been neglected. This might not the case, instead, for annihilations into gauge boson final states. Sommerfeld enhancement can give rise, in particular, to strong γ -ray signals, see e.g [142,[151][152][153][154]. While this effect is particularly pronounced for DM belonging to high SU (2) multiplets, namely n ≥ 3, it is modest, instead, in the case of a SU (2) doublet [140,155,156] unless one considers DM masses of around 7 TeV, well outside the range adopted in In each panel the black curve represents the isocontour of the correct relic density while the blue region is excluded by current constraints from XENON1T. The magenta and purple regions will be ruled out in case of null signals from, respectively, XENONnT/LZ and DARWIN this study. 10 We hence expect that Sommerfeld enhancement does not sensitively affect ID prospects of our scenario.

Results and discussion
We have now all the elements for examining in detail the constraints on the model under consideration. We show first of all in Figs. 4 and 5 the interplay between relic density and Direct Detection, in the bidimensional plane (M L , M N ), for, respectively, Type-I and Type-II scenarios. We have, again, focussed on some specific assignations of the other parameters of the theory and assumed, for simplicity, degenerate masses for the new bosons as well as the alignment limit. As we will see, constraints from Direct Detection are extremely strong, hence we focussed on M N , M L > 100 GeV in order to avoid the regions of maximal sensitivity for these experiments (M L < 100 GeV would be in any case forbidden by LEP limits on production of new charged particles). In each plot the parameter space corresponding to the correct relic density, represented by the black isocontours, is compared with the 10 Sizable Sommerfeld corrections for DM masses of the order of the TeV could be achieved for mass splittings well below the one typically expected from radiative corrections, as could occur in some specific SUSY realizations [156].

Fig. 5
The same as Fig. 4 but for Type-II model (abbreviated Typ-II) and different parameter assignations excluded region (blue) by current limits from Direct Detection, essentially determined by XENON1T [110], as well as the projected sensitivities from XENONnT/LZ [113,114] (magenta, given the similar sensitivity we are assuming the same projected excluded region for both experiments) and DARWIN [115] (purple).
All the panels of Figs. 4 and 5 are grouped in pairs differing by the sign of tan θ , in order to highlight the effect of the blind spots described in Sect. 4.2, since they occur only for negative tan θ .
In the case of the Type-I model the two pairs of benchmarks are mostly differentiated by the assignment of tan θ and tan β. As evidenced by the upper panels of Fig. 4, the moderate value of tan β corresponds to a suppression of the interactions of the DM with the SM fermions, affecting mostly the DD bounds, which are substantially weaker with respect to the case tan β = 2. As already discussed in Sect. 4.2 the blind spots for the Type-I model can be described by the simple equation M N + M L sin 2θ = 0. For tan θ = −10(−6) a blind spot is achieved for M L /M N = 5(3). In the first case the blind spot affects projected excluded regions by XENONnT/LZ and DARWIN while it leaves the current excluded region mostly unchanged. Viable relic density is achieved, for y = 1 and m ψ 1 300 GeV, without relying on blind spots, thanks to the annihilations into the Higgs states. For the pair of benchmarks reported on the lower panels of Fig. 4, on the contrary, the blind spot weakens also the current exclusion bounds and its presence is necessary to achieve the correct DM relic density over a sizeable portion of the parameter space. This viable parameter space would be nev-ertheless almost completely ruled out by a negative signal at XENONnT/LZ. A similar analysis, for two pairs of Type-II benchmarks, is shown in Fig. 5. Contrary to the Type-I model we have chosen a heavier mass of the Higgs mediators, 700 GeV in all panels, in order to comply with the observational constraints discussed in the previous sections. The two pairs of benchmarks are distinguished by the assignations of y, namely 0.2 and 1 for, respectively the upper and the lower panels as well as by the assignations of tan β and tan θ . For the chosen parameter assignations, the blind spot correspond, in good approximation, to the condition g hψ 1 ψ 1 = 0, fulfilled for M L /M N ≈ 4.1, 3.5 for, respectively, the upper and lower panels. As evidenced by the plots, the Type-II model is more constrained. As shown by the upper panels of Fig. 5, viable relic density is achieved, without relying on blind spots, by lowering the value of y to 0.2. On the contrary, by taking y = 1 and a moderate value of tan β, for positive tan θ (implying absence of any blind spot) the whole range of considered values of M N , M L is already excluded by present bounds from DD. By considering negative tan θ a narrow viable region around the s-channel resonance opens. It will be anyway ruled out in absence of signals at next future Direct Detection experiments.
The Type-I model offers another attractive possibility to evade Direct Detection constraints consisting into a light CPodd boson A. In such a case, indeed, it is possible to achieve a sizable s-wave dominated annihilation cross-section of the DM into SM fermions, without strong additional contribution to the scattering cross-section since interactions mediated by a pseudoscalar are momentum suppressed (at least at the tree level). The DM annihilation cross-section can be also enhanced by the presence of h A, A A and Z A final states. Moreover the presence of a velocity independent cross-section would allow ID as complementary probe and possibly fit the GC excess [28,127]. Sizeable constraints from absence γ -ray signals from DSph would be present though.
We have then shown in Fig. 6 On the contrary, the isocontours describing the limit from DSph feature only the cusp associated to the m A resonance, since only the s-channel exchange of the pseudoscalar Higgs leads to a s-wave dominated cross-section (notice that the s-channel exchange of the Z -boson gives a s-wave but helicity suppressed contribution, see Appendix A for details), hence sensitive to ID. Regions of parameters space complying with all observational constraints are nevertheless present. These regions will be, however, fully probed and possibly ruled out by forthcoming Direct Detection experiments.

Conclusions
We have performed an extensive analysis of the DM phenomenology of a model with singlet-doublet Dark Matter coupled with a two doublet Higgs sector. We have considered two scenarios, globally dubbed as Type-I and Type-II model, for the couplings of SM and new fermions with the Higgs doublets. In all cases the most competitive constraints come from limits from Direct Detection. In the case of the Type-II model these can be evaded only by invoking parameter assignations inducing blind-spots in the couplings responsible for Direct Detection. In the case of the Type-I model it is instead possible to evade Direct Detection constraints even without relying on blind spots. The Type-I model presents the additional interesting possibility of a light pseudoscalar Higgs boson.
For all the considered scenarios, next future Direct Detection facilities will fully probe the viable region for thermal DM relic density.
As evident, the s-wave term receives contributions only from s-channel exchange of the pseudoscalar Higgs A and of the Z boson with the latter being, however, helicity suppressed and, hence, relevant only for DM masses close to the mass of the top-quark.
In the regime m ψ 1 < m H,A,H ± the other relevant annihilation channels are in the W W , Z Z and Zh final states. Their cross-sections are given by: where the trilinear couplings used above are given by [33]: The s-wave contributions to the W W and Z Z cross-sections are mostly determined by t-channel exchange of the new neutral and charged fermions, s-channel exchange of the Higgs states is present only in the velocity dependent term. The annihilation cross-section into Zh receives an additional contribution, with respect to the "minimal" singlet-doublet model, from s-channel exchange of the pseudoscalar Higgs A.
In the case that the pseudoscalar Higgs A is lighter than the DM the Z A, h A and A A final states can contribute in a relevant way to the DM relic density. The corresponding cross-sections are given by: and: where: Despite the velocity suppression (v 2 ψ ) the A A and Z A channels can provide not negligible contribution because of the sizeable trilinear scalar couplings.

B Couplings of DM with the Higgs bosons
As already shown in Sect. 3.2 the generic Lagrangian describing the interactions of the new fermions with the two Higgs doublets is (for simplicity we omit this time the mass terms): As already mentioned it is not appropriate to couple each new fermion arbitrarily with both Higgs doublets. This leaves four possible set of assignations of the couplings y L ,R i . Throughout this work we have adopted two of these assignations which, combined with two specific assignations of the couplings of H 1,2 with the SM fermions, have been used to define the models dubbed as Type-I and Type-II. From a bottom-up perspective the choice of the coupling combinations for the new and the SM fermions is substantially free. Consequently the choice of focussing on just two combinations was mostly dictated by definiteness and simplicity.
In this appendix we will briefly discuss about potential changes induced by the choice of different combinations of the couplings as the ones considered in the main text. Leaving to future work the case of the couplings of the Higgs bosons with the SM fermions we will compare here the four possible assignations for the couplings of the SM singlet N and the SU(2) doublets L L ,R with the two Higgs doublets H 1 , H 2 while considering, for all four cases, the same configuration, the Type-II, for the couplings of the Higgs doublets with the SM fermions.
As evidenced in the main text the viable regions for DM phenomenology are mostly determined by the interplay between relic density and (SI) Direct Detection. Concerning the latter, limits are in general very strong, especially when the coupling y is of order 1. It is hence important to identify for which assignations of the model parameters the blind spots occur. As already pointed in the text, and firstly shown in [28], in the case in which the new fermions are coupled only to one Higgs doublet, irrespective whether this is H 1 or H 2 , a blind spot occurs when both couplings g hψ 1 ψ 1 and g H ψ 1 ψ 1 become zero. In all cases this occurs when the condition M N + M L sin 2θ = 0 is fulfilled. In this regard, one of the omitted configurations has been actually accounted in our study. In the case in which the new fermions are selectively Concerning the relic density an analytical description would be similarly difficult because of the different annihilation channels for the DM.
A partial insight is provided by Fig. 7. Similarly to what done in Fig. 2, the isocontours of the correct relic density have been drawn for the three values 0.2,0.5 and 1 of the coupling y while all the other model parameters have been kept fixed. The four panels of the plot correspond to the four possible configurations for the couplings of the new fermionic sector. The assignations adopted for the Type-I and Type-II models studied in this work have been called Case I and Case II respectively (notice that this time the assignation for the coupling of the SM fermions with the Higgs doublets is the same though). The panel dubbed Case III describes the scenarios in which the new fermions couple only with the H 2 doublet; Case IV represent finally the "opposite" coupling configurations as Case II. As evident, while the four panels are qualitatively different, by comparing the upper panels (analogous to the case of study of this work) and the lower panels of Fig. 7, no new feature seems to emerge in the latter. Despite a definitive assessment would require a more systematic study, there is no evidence that potentially relevant cases have been missed. of pp collisions at