Two-Higgs-Doublet-Portal Dark-Matter Models in Light of Direct Search and LHC Data

We explore simple Higgs-portal models of dark matter (DM) with spin 1/2, 3/2, and 1, respectively, applying to them constraints from the LUX and PandaX-II direct detection experiments and from LHC measurements on the 125-GeV Higgs boson. With only one Higgs doublet, we find that the spin-1/2 DM having a purely scalar effective coupling to the doublet is viable only in a narrow range of mass near the Higgs pole, whereas the vector DM is still allowed if its mass is also close to the Higgs pole or exceeds 1.4 TeV, both in line with earlier analyses. Moreover, the spin-3/2 DM is in a roughly similar situation to the spin-1/2 DM, but has surviving parameter space which is even more restricted. We also consider the two-Higgs-doublet extension of each of the preceding models, assuming that the expanded Yukawa sector is that of the two-Higgs-doublet model of type II. We show that in these two-Higgs-doublet-portal models significant portions of the DM mass regions excluded in the simplest scenarios by direct search bounds can be reclaimed due to suppression of the effective DM interactions with nucleons at some ratios of the $CP$-even Higgs bosons' couplings to the up and down quarks. The regained parameter space contains areas which can yield a DM-nucleon scattering cross-section that is far less than its current experimental limit or even goes below the neutrino-background floor.


I. INTRODUCTION
The latest direct searches for weakly interacting massive particle (WIMP) dark matter (DM) by the LUX and PandaX-II Collaborations [1,2] have turned up null results, leading to the strictest upper-limits to date on the cross section of spin-independent elastic WIMP-nucleon scattering in the mass region from about 4 GeV to 100 TeV. For lower WIMP masses down to 0.5 GeV, the existing corresponding limits were set a little earlier in the CRESST [3] and CDMSlite [4] experiments. Ongoing and planned efforts to detect the DM directly [5] will likely improve upon these findings if it still eludes discovery.
The above recent measurements translate into substantial restrictions on WIMP DM models, especially minimal Higgs-portal ones, which are also subject to constraints from continuing quests at the LHC [6,7] for decays of the 125-GeV Higgs boson into final states that would signal new physics beyond the standard model (SM). The combination of restraints from direct detection and LHC data has ruled out in particular the DM mass region below 0.5 TeV or so in the simplest Higgs-portal real-scalar-DM model, with the exception of a narrow range around the Higgs pole [8][9][10]. If the DM is instead a spin-1/2 fermion with a purely scalar effective coupling to the Higgs, only this small region near the Higgs pole remains viable [9,11].
Nevertheless, as previously demonstrated in the Higgs-portal scalar-DM case [8,12,13], by adding another Higgs doublet to the simplest scenario it is possible to decrease the effective interactions of the DM with nucleons sufficiently and thereby to regain at least some of the parameter space disallowed by the direct searches. This motivates us to explore similar ideas in other simple models. Specifically, in this paper we first revisit the minimal Higgs-portal scenarios in which the DM is a fermion of spin 1/2 or 3/2 or a spin-1 boson, to see how the aforesaid restraints impact them. Subsequently, we consider a somewhat expanded version of each of the models by incorporating another Higgs doublet and arranging the new Yukawa sector to be that of the so-called two-Higgs-doublet model (THDM) of type II. We will show that these two-Higgs-doublet-portal DM models, like their scalar-DM counterpart, potentially have ample parameter space that can avoid all the latest direct-detection limits and may even evade future ones. It is worth mentioning here that there have been a number of studies in the past on various Higgs-portal scenarios in which a THDM was supplemented with a SM-gauge-singlet DM candidate having spin 0 [12][13][14][15][16][17], 1/2 [16][17][18], or 1 [17].
The organization of the rest of the paper is as follows. In Sec. II we take another look at the available constraints from DM direct searches and LHC Higgs measurements on the three different minimal Higgs-portal models having WIMP DM candidates with spin 1/2, 3/2, and 1, respectively. In Sec. III, we deal with the two-Higgs-doublet extensions of these models where the enlarged Yukawa sector is that of the type-II THDM. We will address how the extended models can escape some of the restrictions in the presence of sizable breaking of isospin symmetry in the DM effective interactions with nucleons. 1 We give our conclusions in Sec. IV. We collect additional formulas and extra details in a couple of appendices.

A. Spin-1/2 dark matter
In the most economical fermionic scenario [22], the SM is slightly enlarged by the inclusion of a spin-1/2 Dirac field ψ which is a singlet under the SM gauge group and serves as the WIMP DM candidate. In this model, hereafter referred to as SM+ψ, the DM is stable due to an exactly preserved Z 2 symmetry under which ψ → −ψ, while the SM fields are unchanged. Since ψ cannot couple directly to SM members in a renormalizable way, without explicitly introducing other new ingredients one can explore ψ interactions with the SM sector that are induced by effective nonrenormalizable operators. The simplest ones with a Higgs doublet H are the dimension-five combinations ψ(1 ± γ 5 )ψ H † H, which are invariant under the SM gauge group. Assuming that the Lagrangian L ψ for ψ conserves CP symmetry, making the pseudoscalar coupling absent, one can then write [22] where µ ψ and Λ ψ are real constants of dimension mass and Λ ψ absorbs the parameters of the underlying heavy physics. 2 After electroweak symmetry breaking, where h is the physical Higgs and v ≃ 246 GeV the vacuum expectation value (VEV) of H. The DM mass m ψ and the DM-Higgs coupling λ ψh are the only free parameters in L ψ .
The λ ψh terms in Eq. (2) are responsible for the DM relic density. It results fromψψ annihilation into SM particles, which happens mainly via the Higgs-exchange processψψ → h * → X sm . If theψψ pair has a center-of-mass (c.m.) energy √ s > 2m h , the channelψψ → hh needs to be taken into account. Thus, we can express the cross section σ ann of DM annihilation as where the formula for σ ψ ψ → hh is relegated to Appendix A, the sum in the second line is over SM final-states X i,sm = hh, andh refers to a virtual Higgs having a mass mh = √ s.
Subsequently, we can extract λ ψh from the observed DM abundance, as outlined in Appendix A, and then test the result with various constraints.
One of the important restrictions on λ ψh applies in the region m ψ < m h /2, where the invisible decay channel h →ψψ is open and contributes to the Higgs' total width Γ h = Γ sm h +Γ(h →ψψ). From Eq. (2), we obtain the partial rate The experiments at the LHC offer information pertaining to this process. According to the latest combined analysis by the ATLAS and CMS Collaborations on their Higgs measurements [7], the branching fraction of h decay into channels beyond the SM is B exp bsm = 0.00 +0. 16 , which can be interpreted as placing a cap on the invisible decay of h, explicitly B(h → invisible) exp < 0.16. Consequently, we can demand In numerical work, we fix m h = 125.1 GeV, based on the current data [31], and correspondingly the SM Higgs width Γ sm h = 4.08 MeV [32]. Another major constraint on λ ψh is supplied by direct detection experiments, which look for nuclear recoil effects caused by the DM colliding with a nucleon, N, nonrelativistically at momentum transfers small relative to the nucleon mass, m N . In the SM+ψ, this is an elastic transition, ψN → ψN, which is mediated by the Higgs in the t channel. Its cross section is where g N N h is the effective Higgs-nucleon coupling. Numerically, we adopt g N N h = 0.0011 following Ref. [8]. The strongest restraints on σ N el to date for m ψ 5 GeV are provided by LUX [1] and PandaX-II [2].
(2) is the coefficient of a dimension-5 effective operator, the size of λ ψh is also capped by the extent of validity of the effective field theory (EFT) description for the ψ-H interactions. To derive a rough estimate on the minimum of Λ ψ , we entertain the possibility that this operator is induced by a tree-level diagram mediated by a heavy scalar X with mass m X and couplings to ψ and h described by L X ⊃ −g ψψ ψX − g h h 2 X in the UV-complete theory. In addition, we suppose that g ψ ∼ m ψ /v X and g h ∼ λ hX v X , where v X is the VEV of X and λ hX is dimensionless, inspired by the fermionic and scalar couplings in the SM, ignoring potential modifications due to h-X mixing. The EFT will then remain a good approximation and perturbative if 1/|Λ ψ | ∼ 2|λ hX |m ψ /m 2 X < |λ hX |/(2m ψ ) < 2π/m ψ , as the s-channelψψ energy √ s satisfies m 2 X > s > 4m 2 ψ and |λ hX | < 4π for perturbativity. 3 We then have |λ ψh | < 2πv/m ψ . As this follows from the most relaxed requisite on λ hX , it is likely that the EFT description breaks down at a significantly smaller λ ψh . Therefore, alternatively it is reasonable to set |λ hX | < 2, leading to |λ ψh | < v/m ψ . In the m ψ < m h /2 region, this restriction turns out to be far weaker than that from Eq. (5) for the Higgs invisible decay, as will be seen shortly.
To illustrate how the model confronts these different requirements, we display in Fig. 1(a) the values of |λ ψh | derived from the observed DM relic abundance (green solid curve) and compare them to the upper limits on |λ ψh | inferred from Eq. (5) based on LHC data [7] (black dotted curve) and from the validity extent of the EFT approach. For the latter, based on the discussion in the preceding paragraph we draw the magenta band corresponding to the region |λ ψh | ∈ [1, 2π]v/m ψ in which the EFT description may be expected to have broken down. Thus, we can regard the lower boundary of this band as the upper limit for the reliability of the EFT approximation. We determine the range allowed by these constraints to be 54 GeV m ψ 0.8 TeV, which translates into the solid portion of the green curve in Fig. 1(b) for the ψ-N cross-section, σ N el . This green solid curve turns out to be forbidden by the LUX bound, except in a slender zone near the Higgs pole, more precisely 55.8 GeV m ψ 62.3 GeV. Similar results were found in Refs. [9,11]. The ongoing PandaX-II as well as the planned XENON1T, DarkSide G2, and LZ experiments [5] will likely be able to probe the surviving parameter space exhaustively.
It is worth remarking that the λ ψh values in Fig. 1(a) are much bigger than most of their counterparts in the simplest scalar-DM model [8][9][10]. This enlargement is compensation for the suppression of the DM annihilation rate, σ ann v rel , by two powers of the c.m. relative speed v rel ofψ and ψ in the nonrelativistic limit, as can be easily checked. 4 As a consequence, the The corresponding cross-section σ N el of ψ-nucleon elastic scattering (green curve), compared to the measured upper-limits at 90% confidence level from LUX [1] (red dashed curve) and PandaX-II [2] (orange long-dashed curve). Also shown are the sensitivity projections [5] of XENON1T [34] (black dotted curve), DarkSide G2 [35] (purple dash-dot-dotted curve), and LZ [36] (maroon dash-dotted curve), and the WIMP discovery lower-limit due to coherent neutrino scattering backgrounds [37] (brown dashed curve). The dotted portions of the green curve are excluded by the LHC and EFT restrictions in (a). enhanced prediction for σ N el is in conflict with the LUX bound over a wider mass region than in the scalar-DM case.

B. Spin-3/2 dark matter
The WIMP DM in this scenario is described by a Rarita-Schwinger field [38] which is denoted here by a Dirac four-spinor Ψ ν with a vector index ν and satisfies the relation γ ν Ψ ν = 0. 5 In the minimal model, called SM+Ψ, the DM is a SM-gauge singlet, its stability is maintained by an unbroken Z 2 symmetry under which Ψ ν → −Ψ ν , the SM fields being unaffected, and the Higgs-portal interactions arise from dimension-5 operators [26], like in the SM+ψ. The DM Lagrangian, assumed again to be CP invariant, is then where µ Ψ and Λ Ψ are real constants, Λ Ψ encodes the underlying heavy physics, and the DM mass and coupling respectively, are the only free parameters in L Ψ .
Accordingly, we can derive the main quantities relevant to the DM phenomenology in analogy to the spin-1/2 case. Thus, the DM-annihilation cross-section σ ann is given by 6 where β x is defined in Eq. (3), the formula for σ(ΨΨ → hh) is described in Appendix A, and the Higgs' 5 The basic properties of this kind of spin-3/2 fermion, especially in the DM context, have been elaborated in [26,39]. 6 TheΨΨ annihilation rate, σ ann v rel , like its spin-1/2 counterpart, suffers from v 2 rel suppression in the nonrelativistic limit. If L Ψ is not CP -invariant, it can include the pseudoscalar combination Ψ ν γ 5 Ψ ν H † H. In its presence, there is generally an admixture of scalar and pseudoscalar contributions to Ψ-h interactions which can ameliorate the v 2 rel suppression in σ ann v rel .
which is subject to based on the aforementioned LHC Higgs data [7]. For the h-mediated scattering ΨN → ΨN, the cross section is With the formulas above, we arrive at the |λ Ψh | values consistent with the observed relic density and the corresponding σ N el , which are depicted by the green curves in Figs. 2(a) and 2(b), respectively. In Fig. 2(a), the black dotted curve marks the upper bound on |λ Ψh | from the Higgs measurements, and the magenta band represents the zone |λ Ψh | ∈ [1, 2π]v/m Ψ in which the EFT description may be expected to have broken down, analogously to the spin-1/2 case. Consequently, like before the lower boundary of this band may be taken to be the upper limit for the validity of the EFT approximation. The m Ψ range fulfilling these requirements translates into the solid part of the green curve in Fig. 2(b) for the Ψ-nucleon cross-section, σ N el . Evidently, the direct search bounds reduce the viable DM-mass region in the SM+Ψ relative to the SM+ψ. More precisely, only 58.0 GeV m Ψ 61.8 GeV remains viable. This minuscule stretch toward the bottom of the Higgs-resonance dip is almost eliminated by the final LUX limit and will likely be fully excluded by the forthcoming probes of PandaX-II, unless it discovers Ψ.

C. Vector dark matter
In the minimal model, dubbed SM+V, the only new ingredient beyond the SM is a spin-1 state acting as the WIMP DM candidate [24]. We assume that it is described by a real field V ν which is a singlet under the SM gauge group and odd under an unbroken Z 2 symmetry which does not alter SM members. Unlike its fermionic counterparts, V ν can couple to the Higgs doublet via a dimension-four operator. The DM Lagrangian that respects the SM gauge symmetry is then [24] where V κν = ∂ κ V ν −∂ ν V κ and µ V , λ V , and λ h are real constants. Although the terms in L V are at most of dimension four, it is actually nonrenormalizable and violates unitarity [40]. 7 Hereafter, we make no assumption about the details of the UV completion of L V , implying that V is not necessarily a gauge boson and may even be a composite object in the dark sector. Consequently, we can generally treat µ V and λ V,h as independent parameters. For the calculations below, λ V does not play any essential role, and so we can express where m V = µ 2 V + λ h v 2 1/2 is the V mass. As in the fermionic DM models, λ h has to meet the different requirements on DM annihilation, invisible decay h → VV, and DM-nucleon scattering. From Eq. (14), we obtain the formula for σ(VV → hh) is relegated to Appendix A and the h width We will again demand based on the Higgs data [7]. For the Higgs-mediated DM-nucleon collision VN → VN, the cross section is There are also theoretical considerations relevant to restraining λ h . Since L V in Eq. (14) leads to unitarity violation [40], we need to ensure that it does not occur with the extracted λ h values. As discussed in Appendix B, this implies that we need to have |λ h | < √ 2π m V /v. There is additionally a complementary restraint from the expectation that the theory remains perturbative. As also explained in Appendix B, this translates into the bound |λ h | < 1.
In Fig. 3(a), we present the λ h values fulfilling the relic density requirement (green curve). Also shown are the upper limits from the Higgs invisible decay data (black dashed curve) and from the unitarity and perturbativity considerations (maroon dashed curves). We plot the corresponding V -nucleon cross-section from Eq. (18) in Fig. 3(b), where the dotted sections of the green curve are disallowed by the restrictions in Fig. 3(a). We find that m V values approximately below 54 GeV and from 62.6 GeV to 1.42 TeV are in conflict with LHC and LUX data. Comparable results were obtained in Ref. [9]. However, the graphs also reveal that for m V > 3.9 TeV the effective theory probably is no longer perturbative. Thus, overall the situation is rather similar to that in the simplest scalar-DM model [8][9][10], where the annihilation rate does not suffer from the v 2 rel suppression and consequently the viable parameter space is far greater than in its fermionic counterparts. The corresponding cross-section σ N el of DM-nucleon scattering (green curve) compared to the same data and projections as in Fig. 1(b). The dotted sections of the green curve are disallowed by the constraints in (a).

III. TWO-HIGGS-DOUBLET-PORTAL FERMIONIC AND VECTOR DM MODELS
Here we explore extensions of the minimal models in the last section by adding in each case another Higgs doublet. Furthermore, we suppose that the SM fermions in the extended scenarios have the Yukawa interactions of the two-Higgs-doublet model (THDM) of type II, where the down-type fermions get mass from only one of the Higgs doublets, H 1 , and the up-type fermions from the other doublet, H 2 . Accordingly, the Yukawa Lagrangian is [42,43] where j, l = 1, 2, 3 are summed over, Q j,L (L j,L ) stands for left-handed quark (lepton) doublets, U l,R and D l,R (E l,R ) are right-handed quark (charged lepton) fields,H 1,2 = iτ 2 H * 1,2 with τ 2 being the second Pauli matrix, and λ u,d,ℓ represent 3×3 Yukawa-coupling matrices. This Lagrangian respects another discrete symmetry, Z ′ 2 , under which H 2 → −H 2 and U R → −U R , while all the other fields are not affected. Thus, Z ′ 2 prohibits the combinations Q LH1 U R , Q L H 2 D R , L L H 2 E R , and their Hermitian conjugates from entering L Y .
In the scalar sector, the renormalizable potential V H is that of the THDM II, Although dimension-4 combinations with an odd number of H ( †) 2 cannot appear due to Z ′ 2 , in V H we have allowed the m 2 12 terms which softly break Z ′ 2 and are important in relaxing the caps on the Higgs masses [43]. The Hermiticity of V H implies that m 2 11,22 and λ 1,2,3,4 are real parameters. With V H chosen to be CP invariant, m 2 12 and λ 5 are also real constants. To see how V H describes the couplings among the physical states in the Higgs doublets, we first decompose them as where v 1,2 denote the VEVs of H 1,2 , respectively, and are linked to the electroweak mass scale v ≃ 246 GeV by v 1 = v cos β and v 2 = v sin β. The components h + r , h 0 r , and I 0 r are connected to the physical Higgs bosons h, H, A, and H + by where X is any angle or combination of angles and w ± and z are, respectively, the would-be Goldstone bosons that will be eaten by the W ± and Z bosons. We can then express the terms in V H after electroweak symmetry breaking that are relevant to our purposes as where the λs are linked to the physical Higgs masses [43] and the relations are listed in Ref. [8].
Since h and H couple directly to the W and Z bosons, we need to take into account DM annihilation into W + W − and ZZ. The pertinent interactions are given by The presence of the extra Higgs doublet not only offers another portal between the dark and visible sectors, but also can produce modifications to the effective coupling between a Higgs boson H and a nucleon N which is defined by and plays a crucial role in DM-nucleon collisions. The potential changes spring from the quark-Higgs terms in Eq. (19) where the sum is over all quarks, q = u, d, s, c, b, t, and which are generally different from the SM values k h q = 1 and k H q = 0. Relating the above quarkand hadronic-level quantities, one arrives numerically at [8,12,14] Setting k h u,d = 1 in the H = h formulas, we reproduce the SM values g sm pph,nnh ≃ 0.0011 quoted in the last section. However, if k H u,d are away from their SM expectations, g ppH and g nnH can be very dissimilar, manifesting considerable isospin-violation. Especially, if k H u,d have opposite signs, it may be possible to reduce g N N H such that the DM effective interactions with nucleons become weak enough to evade the experimental constraints.
This suggests that to address DM-nucleon collisions in models with a THDM-II portal it is more appropriate to work with either the DM-proton or -neutron cross-section (σ p el or σ n el , respectively) rather than the DM-nucleon one (σ N el ) under the assumption of isospin conservation. For comparison with experiment, in case the DM effective interactions with nucleons violate isospin, the computed σ p,n el can be converted to σ N el by means of [20] σ N el i where each sum is over isotopes of the element in the target material with which the DM interacts dominantly, η i (A i ) stand for the fractional abundances (the nucleon numbers) of the isotopes, In what follows, we select h to be the 125-GeV Higgs boson. Accordingly, k h d,u,V in Eq. (27) need to be compatible with LHC measurements on the h couplings to SM fermions and electroweak bosons. The modification to the h → XX interaction due to new physics can be parameterized by κ X in κ 2 X = Γ h→XX /Γ sm h→XX . Assuming that |κ W,Z | ≤ 1 and the total width of h can be altered by decay modes beyond the SM, the ATLAS and CMS Collaborations have carried out simultaneous fits to their Higgs data to extract [7] κ W = 0.90 ± 0.09 , κ t = 1.43 +0.23 −0.22 , |κ b | = 0.57 ± 0.16 , |κ γ | = 0.90 +0.10 −0.09 , κ Z = 1.00 −0.08 , where [7] κ 2 γ = 0.07 κ 2 t + 1.59 κ 2 W − 0.66 κ t κ W . In the context of the THDM II, we expect these numbers to obey the relations k h V = κ W = κ Z , k h u = κ t ≃ κ g , and k h d = κ b = κ τ within one sigma, but the κ t,g (κ b,τ ) numbers above overlap only at the two-sigma level. Pending improvement in the precision of these parameters from future data and following Ref. [8], based on Eq. (30) we can then impose There are other constraints that we need to consider as well. The extra Higgs particles in the THDM generally modify the so-called oblique electroweak parameters S and T encoding the impact of new physics coupled to the standard SU(2) L gauge boson [44], and so the new scalars must also conform to the empirical requisites on these quantities. To ensure this, we use the pertinent results of Ref. [45] and the S and T data from Ref. [46].
Theoretically, the parameters of the scalar potential V H in Eq. (20) need to fulfill a number of conditions. To keep the theory perturbative, each of the quartic couplings in V H cannot be too big. Another requirement is the stability of V H , implying that it has to be bounded from below. It is also essential to check that the (tree level) amplitudes for scalar-scalar scattering do not violate unitarity constraints. We summarize the expressions pertaining to these conditions in Appendix B.
In the rest of this section, we treat in turn the THDM II+ψ, THDM II+Ψ, and THDM II+V, which are respectively the type-II two-Higgs-doublet extensions of the minimal models of the previous section. Thus, the DM stability is maintained in each case by the exactly conserved Z 2 symmetry as before, under which the DM is odd and all the other fields are even. In addition, we demand that, besides the scalar potential, the DM sector be CP invariant. We will demonstrate that in the presence of the second doublet it is possible to have substantial weakening of the constraints from DM direct detection experiments or perhaps even to evade them in the future.

A. THDM II+ψ
In this scenario, the Lagrangian for the DM is [17,18] where µ ψ and Λ 1ψ,2ψ are real constants of dimension mass and Λ 1ψ,2ψ contain the parameters of the underlying heavy new physics. The Z ′ 2 symmetry prevents the combinations ψψH † 1 H 2 and ψψH † 2 H 1 from appearing in L ′ ψ . After electroweak symmetry breaking, we can express the relevant terms in L ′ ψ involving the physical bosons as where There is noψψA term under the assumed CP conservation. Since µ ψ and Λ 1ψ,2ψ are free parameters, so are m ψ and λ ψh,ψH . The couplings ofψψ to a pair of Higgs bosons can then be expressed in terms of λ ψh,ψH as If both the h and H couplings to ψ are nonzero, the DM-nucleon scattering ψN → ψN proceeds via tree-level diagrams mediated by h and H, leading to the cross section for momentum transfers small relative to m h,H . Given that the Higgs-nucleon coupling g N N H , for N = p or n and H = h or H, depends on k H u,d according to Eq. (28), it may be possible to suppress g N N H sufficiently with a suitable choice of k H d /k H u to make σ N el evade its experimental limit [12], at least for some of the m ψ values. In addition, the λ h,H terms in Eq. (36) may (partially) cancel each other to lower σ N el as well. These are appealing features of the two-Higgsdoublet scenario that the one-doublet case does not possess. In evaluating model predictions for DM-nucleon reactions later on, we work exclusively with the DM-proton cross-section, σ p el , and then convert it to σ N el with the aid of Eq. (29) for comparison with measurements. As there are countless different possibilities in which h and H may act as portals between the DM and other particles, for definiteness and simplicity hereafter we concentrate on two scenarios in which the 125-GeV Higgs boson h is lighter than the other Higgs bosons, m h < m H,A,H ± . Moreover, we assume particularly that either H or h has a vanishing coupling to the DM, λ ψH = 0 or λ ψh = 0, respectively. Accordingly, either h or H alone serves as the portal, and hence now we have f n /f p = g nnH /g ppH , after ignoring the n-p mass difference.
In the h-portal scenario (λ ψH = 0), the cross section of DM annihilation is where σ ψ ψ → h * → X sm is equal to that in Eq. (3), except the couplings of h to fermions and gauge bosons are multiplied by the appropriate k h u,d,V factors mentioned earlier, the sum is over s 1 s 2 = hh, hH, HH, AA, H + H − with only kinematically permitted channels contributing, and the formulas for σ ψ ψ → s 1 s 2 are listed in to Appendix A. After extracting λ ψh from the relic density data and evaluating g pph with the α and β choices consistent with Eq. (27), we can predict the cross section of the h-mediated transition ψp → ψp, As in Sec. II A, for 2m ψ < m h the invisible channel h →ψψ is open, its rate already written down in Eq. (4). Then B h →ψψ must be consistent with the LHC information on the Higgs invisible decay, and so for this m ψ range we again impose the bound in Eq. (5).
In the H-portal case (λ ψh = 0), the DM-annihilation cross-section is where withH being a virtual H with mass mH = √ s. Given that H is not yet discovered, no empirical restraint on H →ψψ exists. For the ψ-proton scattering via H exchange, the cross section is In applying Eq. (29), we set f n /f p = g nnH /g ppH .
To illustrate the viable parameter space in these λ ψH(ψh) = 0 scenarios, we put together in Table I (II) sample sets of input parameters (the second to seventh columns) which are compatible with Eq. (31) and the other requirements described in the two paragraphs following it. The eighth to fifteenth columns of the tables contain the resulting values of several quantities.  With these input numbers, we show in Figs. 4(a) and 4(b) the λ ψh and λ ψH regions evaluated from the observed relic density. One observes that the |λ ψH | values extracted from the relic density data tend to be bigger than their λ ψh counterparts. This is because the H-portal annihilation rate is relatively more suppressed due to m H > m h . In Fig. 4(a), we also display the upper bound on |λ ψh | inferred from Eq. (5) for the h invisible decay (black dotted curve). Like in the minimal model of Sec. II A, the limited extent of the reliability of the EFT approximation for the ψ-Higgs operators in L ′ ψ implies that in each H-portal instance we also need to ensure |λ ψH | < v/m ψ beyond which the EFT framework may be expected to break down. This condition is represented by the lower sides of the magenta bands in Figs. 4(a) and 4(b). We exhibit the corresponding predictions for σ p el in Fig. 4(c), where the dotted parts of the green and blue curves are excluded by the constraints in Figs. 4(a) and 4(b), respectively.
To test the model with direct search results, which are typically reported in terms of the DM-nucleon cross-section σ N el , we have converted the calculated σ p el in Fig. 4 to the (green and blue) σ N el curves in Fig. 5(a) by means of Eq. (29) with the f n /f p values from the tables, assuming that the target material in the detector is xenon. Since the DarkSide G2 experiment will employ an argon target [35], in Fig. 5(b) we plot the corresponding predictions for σ N el assuming an argon target instead. These graphs reveal some visible differences, especially the xenon curves for Sets 2 and 4 which are significantly lower than their argon counterparts. The differences are not unexpected because the f n /f p numbers in these instances are not far from the xenophobic extreme, f n /f p ≃ −0.7.  Table I. Also plotted are upper limits inferred from LHC data on the Higgs invisible decay (black dotted curve) and from the validity of the EFT approximation (lower sides of magenta bands). (b) The same as (a), except with λ ψh = 0 and input numbers from Sets 3 (blue solid curve) and 4 (blue dashed curve) in Table II, Table I and Sets 3 and 4 (blue curves) in Tables II for (a) xenon and (b) argon targets, compared to the same data and projections as in Fig. 1(b). The dotted portions of the green and blue curves are excluded as in Fig. 4.
Also depicted in Fig. 5 are the same data and projections as in Fig. 1(b). It is obvious that, in stark contrast to the SM+ψ, the THDM II+ψ accommodates a good amount of parameter space which can evade the current direct search restrictions very well. Particularly in these examples, over wide stretches of m ψ the model prediction can also escape future direct detection and even hide below the neutrino floor. In the h-portal cases, the LHC Higgs invisible decay data and the EFT validity limit rule out m ψ < 54 GeV and m ψ > 0.9 TeV, respectively. In the H-portal ones, the LHC Higgs invisible decay restraint does not apply, but the EFT validity limit disallows m ψ values below 165 GeV and above 0.7 TeV.

B. THDM II+Ψ
In this model, the Lagrangian for the spin-3/2 DM is where µ Ψ and Λ 1Ψ,2Ψ are real constants of dimension mass and the latter two contain the parameters of the underlying heavy new physics. After electroweak symmetry breaking, we can then express the relevant terms in L ′ Ψ involving the physical bosons as where m Ψ and the λs are the same in form as their counterparts in Eqs. (34) and (35), but with ψ in the subscripts replaced by Ψ.
It follows that the DM-annihilation cross-section in the h-portal scenario is where σ Ψ Ψ → h * → X sm is equal to that in Eq. (9), except the h couplings of to SM particles are scaled by the suitable k h u,d,V factors, the sum is again over s 1 s 2 = hh, hH, HH, AA, H + H − , and the formulas for σ Ψ Ψ → s 1 s 2 are collected in Appendix A. As in Sec. II B, for 2m Ψ < m h the invisible channel h →ΨΨ is open, its rate given by Eq. (10), and so it must fulfill the condition in Eq. (11). The Ψ-proton cross-section is the same as that in Eq. (38) but with ψ in the subscripts replaced by Ψ.
In the H-portal scenario where The cross section of H-mediated Ψ-proton scattering is equal to that in Eq. (41), but with ψ in the subscripts replaced by Ψ.
Similarly to the THDM II+ψ, we use the numbers from Tables I and II for our examples. We present the results in Figs. 6 and 7. These instances indicate that in the THDM II+Ψ the situation is roughly similar to that in the THDM II+ψ, but the parameter space in the former is less able than in the latter to escape the different restrictions. Specifically for these examples, we estimate the viable zones in the h-and H-portal cases, respectively, to be 56 GeV m Ψ 420 GeV and 200 GeV m Ψ 450 GeV, which translate into predictions for σ N el which are currently below its experimental limits.

C. THDM II+V
In this scenario, the Lagrangian for the vector DM is [17] L where, analogously to the fermionic models, We again look at separate possibilities with either λ H = 0 or λ h = 0.
Thus, the cross section of the h-portal DM-annihilation is where σ(VV → h * → X sm ) equals that in Eq. (15), but with the h couplings to SM members multiplied by the proper k h u,d,V factors. The expressions for σ(VV → s 1 s 2 ) can be found in Appendix A. The invisible channel h → VV has the rate already given in Eq. (16) and hence is also constrained by Eq. (17). For the h-mediated DM-nucleon collision VN → VN, the cross section is In the H-portal scenario, the cross section of DM annihilation is where The cross section of the H-mediated VN → VN scattering has the formula in Eq. (50), but with h in the subscripts replaced by H.
Employing these formulas with the input numbers from Tables I and II for our examples, we arrive at the green and blue curves in Figs. 8 and 9. In Fig. 8(a), we also draw the upper bound on |λ h | from the Higgs data (black dotted curve). As in the SM+V, theoretical considerations concerning unitarity and perturbativity lead to the constraints |λ H | < √ 2π m V /v and |λ H | < 1, respectively, which are represented by the magenta dashed lines in Fig. 8. These instances illustrate that the situation in this model is roughly similar to that in the THDM II plus real scalar DM investigated in Ref. [8]. Specifically, in our h-portal examples the LHC constraint on h → VV rules out m V < 54 GeV and the effective theory is likely to be no longer perturbative for m V > 5 TeV, whereas in the H-portal cases the unitarity and perturbativity conditions disfavor m V < 110 GeV and m V > 3 TeV. Thus, more generally, in great contrast to the THDM II+(ψ, Ψ) as well as the SM+V, the THDM II+V has plentiful parameter space that is far away from the current direct search limits and can even escape future ones.   Table I and Sets 3 and 4 (blue curves) in Tables II for (a) xenon and (b) argon targets, compared to the same data and projections as in Fig. 1(b). The dotted portions of the green and blue curves are disallowed as in Fig. 8.

IV. CONCLUSIONS
The most recent limits from LUX and PandaX-II on DM-nucleon interactions and the available LHC data on the 125-GeV Higgs boson's couplings have led to strong restrictions on the simplest Higgs-portal dark matter models. Taking these constraints into account, we have revisited the minimal models with fermionic DM having spin 1/2 or 3/2 and a purely scalar effective coupling to the Higgs doublet. Realizing also that the EFT description for this coupling has limitations, we have found that these minimal fermionic models are ruled out except in narrow regions of the DM mass in the neighborhood of the Higgs resonance point at m h /2. On the other hand, the simplest Higgs-portal vector-DM model is viable not only around m V = m h /2, but also above m V ∼ 1.4 TeV, although it may lose perturbativity if m V > 3.9 TeV. Slightly expanding each of these models with the addition of another Higgs doublet and assuming the extended Yukawa sector to be that of the type-II two-Higgs-doublet model, we can significantly relax the restraints from direct search and LHC data, even in the fermionic DM scenarios, and recover sizable parts of the regions excluded in the minimal models. This is due to suppression of the effective interactions of the portal CP -even Higgs bosons with nucleons at some values of the ratios of the Higgs couplings to the up and down quarks, rendering the interactions considerably isospinviolating. Sizable portions of the revived parameter space can yield a DM-nucleon scattering cross-section that is much smaller than its current experimental bound or even falls under the neutrino-background floor. Nevertheless, there are also areas in the parameter space of these two-Higgs-doublet-portal DM models that are still within the discovery reach of future quests such as XENON1T and LZ.
is the VEV of X, in analogy to the scalar and weak-boson couplings in the SM, ignoring potential modifications due to h-X mixing. The EFT will then remain a reliable approach and perturbative if 2|λ h | ∼ 4|λ hX |m 2 V /m 2 X < |λ hX | < 4π, as the s-channel energy √ s satisfies m 2 X > s > 4m 2 V . Since it is likely that the EFT description starts to break down at a lower |λ hX |, as was also suggested in Sec. II for the fermionic cases, it is more reasonable to select |λ h | < 1 instead. Similarly, we impose |λ H | < 1 in the H-portal THDM II+V scenarios.