Dark Matter and the Higgs in Natural SUSY
Aria Basirnia, Sebastian Macaluso and David Shih
NHETC, Department of Physics and Astronomy, Rutgers University,
Piscataway, NJ 08854
Null results from dark matter (DM) direct detection experiments and the 125 GeV Higgs both pose serious challenges to minimal supersymmetry. In this paper, we propose a simple extension of the MSSM that economically solves both problems: a “dark sector” consisting of a singlet and a pair of doublets. Loops of the dark sector fields help lift the Higgs mass to 125 GeV consistent with naturalness, while the lightest fermion in the dark sector can be viable thermal relic DM, provided that it is mostly singlet. The DM relic abundance is controlled by wave annihilation to tops and Higgsinos, leading to a tight relation between the relic abundance and the spindependent direct detection cross section. As a result, the model will be fully probed by the next generation of direct detection experiments. Finally we discuss the discovery potential at LHC Run II.
1 Introduction and Summary
The MSSM paradigm is under siege from both the LHC and dark matter (DM) direct detection. The Higgs mass at treelevel in the MSSM is famously bounded by , and relying on radiative corrections from stops and other particles in the MSSM forces the stops to be either at least TeV or their terms to be multiTeV (for recent reviews and original references, see e.g. [1, 2, 3]). Together with the null direct search results at the LHC, this puts the finetuning in the MSSM at the percent level or worse. Meanwhile, to evade stringent direct and indirect detection bounds, thermal relic neutralino DM in the MSSM must rely on increasingly contrived numerical accidents (welltempering, blind spots, funnels, coannihilations) or an increasingly heavy SUSY scale (e.g. TeV Higgsinos or TeV winos) (see e.g. [4, 5, 6] for recent comprehensive studies). The latter constitutes a DM version of the little hierarchy problem, whereby the WIMP miracle’s preference for TeVscale DM (as opposed to 100 GeV scale DM) is in tension with naturalness.
This strongly motivates looking beyond the MSSM for both the source of the Higgs mass and dark matter. Although it is logically possible that different sectors are independently responsible for the Higgs mass and dark matter, it is interesting to contemplate more elegant and economical models where a single sector generates both. In this paper, we will study such a model. We will show how to achieve a 125 GeV Higgs and thermal relic WIMP DM consistent with all existing constraints, while greatly ameliorating the finetuning, by just adding a pair of doublets , and a singlet to the MSSM. With a “DM parity” that keeps the lightest state in the dark sector stable, together with matter parity from the MSSM, the most general renormalizable superpotential for this “dark sector” is:
(1.1) 
Although it would be interesting to also consider phases, we will focus on real couplings in this paper for simplicity. Then without loss of further generality, we can take and to be positive.
The idea of extending the Standard Model (SM) with a “singletdoublet DM” sector has been studied previously in [7, 8, 9, 10, 11, 12, 13, 14], motivated by minimality and by the fact that it is a simple generalization of the wellstudied bino/Higgsino system of the MSSM. The idea of lifting the Higgs mass with loops of vectorlike matter has also been wellstudied [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. But to our knowledge, the two ideas have never been combined before.^{1}^{1}1A related idea [33, 34] is to use vectorlike matter to boost the Higgs mass while simultaneously opening up new annihilation channels for binolike DM. Combining these two ideas leads to some important differences with previous works.
First, unlike in previous works on lifting the Higgs mass, our dark sector cannot be truly vectorlike. The scalar soft masssquareds of the dark sector must be positive in order to lift the Higgs mass, making our DM the lightest fermion in the dark sector. It cannot be a Dirac fermion, otherwise it would be ruled out by many orders of magnitude by mediated spinindependent (SI) direct detection. Instead, we make the dark sector fermions Majorana (as shown in (1.1)) by having only one singlet and not a vectorlike pair of them. This only has a minor effect on the contribution to the Higgs mass in this model, which we fully take into account. We will find that a GeV Higgs can be achieved with the finetuning coming from the DM being only 10%, provided that .
Second, we differ from the singletdoublet DM models in that we are supersymmetrizing everything.^{2}^{2}2Actually, in [8] they also added singlets and doublets to the MSSM. However, they considered soft masses purely from GMSB (whereas we are agnostic) and therefore they never have mostlysinglet fermionic DM. Moreover they fix whereas we have them as free parameters. Finally, they do not calculate the contribution to the Higgs mass from the dark sector. A priori, the parameter space of the entire model (MSSM+dark sector) is vast, but most of the soft parameters do not play a significant role in the analysis. As seen in (1.1), our dark sector only couples directly to the Higgs sector and the EW gauge sector of the MSSM. We will keep the Higgsinos light ( GeV), since they contribute to the finetuning of the EW scale at tree level. As a result, DM annihilation to light Higgsinos through superpartners in the dark sector plays a major role in determining the relic abundance of the DM. Meanwhile, it does not change our analysis qualitatively to decouple all other MSSM superpartners (effectively at the TeV scale). This is further motivated by the null results from the LHC.
We will further simplify the analysis of the model by focusing on the regime where the dark matter is mostlysinglet, i.e. and . As we will argue in much more detail in section 5, this regime is absolutely necessary in order to evade direct detection bounds while raising the Higgs mass without finetuning. A key part of the argument, which distinguishes this from the bino/Higgsino system in the MSSM, is that must be in order to lift the Higgs mass without finetuning. This eliminates both the welltempered regime and the mostlydoublet regime vis a vis DM direct detection. The mostlydoublet regime is further unpromising because (by analogy with pure Higgsino DM in the MSSM) it would require a DM mass in excess of 1 TeV, and this would greatly exacerbate the finetuning problem, since the rest of the dark sector would have to be even heavier. This leaves the mostlysinglet regime, where the analysis of the model greatly simplifies, and we are able to understand all the features of the model with simple analytic formulas. A cartoon spectrum of the model that describes these hierarchies qualitatively is shown in fig. 1.
In this work, we will assume the simplest DM scenario, namely that is a thermal relic comprising all of the DM. In the mostlysinglet limit with , we will show that the thermal relic abundance is controlled by just two DM annihilation channels: wave (through channel exchange) and wave Higgsinos (through channel superpartner exchange). Assuming for simplicity, we find:
(1.2) 
where is a common soft mass for the dark sector scalars. As noted above, the second term coming from Higgsinos is a major difference from the nonsupersymmetric singletdoublet DM models that have been studied previously. Having more annihilation channels increases , making it possible to have smaller effective couplings between the DM and the SM. This opens up more parameter space that is not ruled out by direct detection experiments and yet still has the correct thermal relic abundance, as compared to the nonSUSY singletdoublet models.
Interestingly, the DM mass drops out of the annihilation cross section (1.2) in the mostlysinglet limit. The WIMP miracle becomes one for the mediator scale, not the WIMP mass! With , and GeV, mediator scales of TeV are implied by the thermal relic constraint. Meanwhile the DM can be much lighter than this, alleviating the DM little hierarchy problem. It is also interesting to contrast this with the mostly bino limit of the bino/Higgsino system in the MSSM. There the annihilation cross section is not large enough, being suppressed by instead of . Our model (and singletdoublet DM more generally) gets around this in the mostlysinglet regime with Yukawa couplings that are free parameters, not fixed to be by supersymmetry.
Meanwhile, DM direct detection in these models is completely controlled by the effective couplings of the DM to the Higgs and respectively:
(1.3) 
As is wellknown, () controls the SI (SD) direct detection cross section. For direct detection, as we will review, the current best bounds for our DM mass range of interest ( GeV) come from LUX [35] and IceCube [36]. We will convert the official experimental results, which are phrased in terms of the DMnucleon cross section, into limits on and . Furthermore, in the mostlysinglet limit, we will obtain simple analytic expressions for and . We will see that can be naturally small enough for mostlysinglet DM, due to suppression from the heavier doublet scale, as well as a mild blindspot cancellation:
(1.4) 
provided that . We should emphasize here that the Higgs mass depends not just on but also on the effective Yukawa couplings between the Higgs and the other dark sector particles. So even dialing does not qualitatively affect the Higgs mass calculation. Meanwhile is given in the mostlysinglet limit by:
(1.5) 
According to our discussion above, after fixing the thermal relic density constraint , is essentially fixed to lie within a narrow range which depends primarily on the Higgsino mass . Therefore imposing the relic density constraint essentially fixes the SD cross section. Fortunately, this value is not ruled out yet, but the next generation of DM experiments (e.g. Xenon1T [37], LZ [38]) should completely rule out or discover this model.
Although direct detection is controlled by and , the other facets of the model (relic abundance, Higgs mass) depend on more than just these couplings, so our model does not fit completely into the framework of  and portal DM. For instance, we mentioned above that the Higgsino cross section arises entirely from channel superpartner exchange. Also, we find that DM annihilation to dibosons is suppressed more than would be the case in and portal models, in part due to channel exchange of doublet fermions. Similar comments apply to the effective operator formalism: our DM is generally not light enough compared to the mediator scale (the doublet mass) for the annihilation to be accurately captured by effective operators. Evidently, the complete model (1.1) is required for an accurate analysis. This illustrates the shortcomings and limitations of both simplified models and effective operator approaches to dark matter.
We have focused primarily on the standard direct detection searches in this work, because other indirect probes of our dark sector are far less sensitive. For example, the Fermi experiment and others have searched for energetic photons produced through DM annihilating at the centers of dwarf galaxies. For DM masses above GeV, Fermi does not constrain any point with the right relicabundance [39], assuming (as is the case for us) that the relic abundance is determined by wave annihilation. Meanwhile, searches at colliders and electroweak precision tests (EWPT) could have put constraints on our model. However as we will discuss further in section 8.2, LHC bounds [40, 41, 42, 43] on and from monojets+MET and monophoton+MET are orders of magnitude weaker than direct detection for the range of DM masses that we are interested in. We will briefly discuss mono(,,)+MET and show how it could probe the low end of DM masses ( GeV) in our model, with 300/fb at LHC Run II. Finally, limits from Higgs and invisible width do not apply to the mass range of DM that we consider in this work, and we checked that contributions to the and parameters are well within the acceptable range, in agreement with previous studies of these variables in closelyrelated models [11, 19].
In this paper, we will analyze the model using a combination of simple, approximate analytic expressions valid in the mostlysinglet regime, and more precise numerical methods that take into account the full suite of one (and even two) loop threshold corrections. The analytic approach, while being reasonably accurate, is primarily a source of intuition and understanding. The numerical approach is meant to be more accurate and to provide us with the quantitative results. Clearly, having both numerics and analytics is a vital source of crosschecks, giving us greater confidence in our results.
Our numerical methods are based on publicly available codes. Our starting point was the powerful SARAH 4.5.8 framework [44] for automated analysis of general models. Once we properly defined our model, SARAH automatically generated source code for SPheno 3.3.7 [45, 46] and for micrOMEGAs 4.1.8 [47].^{3}^{3}3We are extremely grateful to Florian Staub for his time and patience in helping us set up the SARAH model and link it to these other codes. The former calculates the spectrum while the latter calculates the DM relic abundance and direct detection cross sections. In our numerical calculations, all MSSM soft masses as well as gauginos are taken to be at 1 TeV, and the terms are set to zero. As noted earlier, since appears at tree level in finetuning of the electroweak scale we treat it differently. We pick in our numerical calculations which corresponds roughly to finetuning. We also consider to see the effect of on our analysis. Finally, to saturate the tree level contribution to the Higgs mass, we take the other Higgses to be heavy and in the decoupling limit, and we take .
The outline of our paper is as follows. In section 2 we introduce the model. Then in section 3, we derive direct detection limits from LUX and IceCube on the effective couplings and . We will emphasize that these results are general and are not limited to the model we consider in this work. In section 4 we compute the oneloop corrections to the Higgs mass from the new particles in our model, and we discuss finetuning. We argue in section 5 that the mostly singlet case is the only viable scenario. In the mostly singlet limit, we provide analytic expressions for dark matter annihilation in the early universe for our model in section 6. In section 7 we put everything together to show the viable parameter space that satisfies all direct detection constraints while having the right relic abundance and Higgs mass. Here we demonstrate quantitatively that requiring to be all of the DM essentially fixes (and hence ) to a unique value which is not yet ruled out by direct detection, but will be fully within reach of the next generation of experiments. We conclude by studying the collider signatures for LHC Run II and the UV behavior of the model, and giving suggestions on future directions on section 8. Technical details and validations are reserved for three appendices. In appendix A we review the derivation of the direct detection cross sections from effective DM nucleon couplings. We validate our numerical and analytical calculations of the Higgs mass in appendix B. Finally we provide analytical cross sections for DM production at LHC II in appendix C.
2 The Model
We begin by describing the model in more detail. We add to the MSSM a “dark sector” consisting of a vectorlike pair of doublets , and a gauge singlet .^{4}^{4}4To keep gauge coupling unification as in MSSM, we can assume and are part of complete and multiplets of . We take their colored partners to be heavy and decoupled for simplicity. The dark sector is equipped with an unbroken parity symmetry under which all new fields are odd and all MSSM fields are even. This makes the lightest new state stable and a DM candidate. Finally, we assume MSSM matter parity, under which all the dark sector fields have the same charge; otherwise there will be additional, potentially dangerous terms.^{5}^{5}5The assumption of matter parity implies another stable particle – either the LSP in the MSSM, or the gravitino. Either way, we assume the parameters are such that this will add a negligible additional component to the thermal relic density. This would be the case, for instance, if the LSP is a light Higgsino. The transformation properties of the dark sector under the gauge and global symmetries is summarized in tab. 1.
The most generic superpotential consistent with these symmetries is:
(2.1) 
The superpotential has four new parameters in addition to the MSSM: . There is one physical complex phase, but as discussed in the introduction, we will take these parameters to be real in this paper. In this case, there is still a physical sign. We will take , and to be positive and put the sign into .
For the soft SUSYbreaking Lagrangian, for simplicity we take the minimal case with equal soft masssquareds and no  or terms:
(2.2) 
(We denote the scalar components of the dark sector superfields with lowercase letters.) Allowing different soft masses for the different fields will not change most of the discussion in this paper, only the contributions to Higgs mass.
As we want this new sector to increase the lightest Higgs mass analogous to the MSSM stops, we assume that . This implies that the DM candidate is a fermion. Furthermore it is Majorana, thanks to the fact that we have included only one singlet in the theory. Had we started with a Dirac pair of and and defined the mass term as , our dark matter would have had a vectorlike coupling to the . In that case it would have been impossible to hide it from SI direct detection experiments while keeping the interesting features of our model.
After EWSB, neutral fields in the dark sector mix through the Yukawa couplings in (2.1). The fermion mass matrix of the neutral states is:
(2.3) 
where we have introduced and , with and as usual. We take large in this paper to saturate the upper bound on the tree level Higgs mass. The mass matrix is diagonalized by . The spectrum of the model consists of three Majorana fermions with masses and a Dirac charged fermion with mass . The dark matter candidate is then .
We note that the fermionic part of our dark sector is analogous to BinoHiggsino DM in the MSSM (with everything else decoupled), except that in the BinoHiggsino system, we effectively have , whereas here and are general. In fact, as discussed in the introduction, here we will be primarily interested in , .
After rotating to the mass eigenbasis, DM and DMHiggs couplings are generated:
(2.4) 
where is a 4component Majorana fermion and and are given by:
(2.5) 
and
(2.6) 
As is wellknown [9, 10, 11, 12, 13], and play an important role in the analysis of singletdoublet DM: they are entirely responsible for SI and SD direct detection, respectively. In the next section, we will review the current direct detection constraints on and .
3 DM Direct Detection through the and Portals
In the DM mass range of interest (100 GeV 1 TeV), the LUX experiment currently sets the best bound on SI elastic WIMPnucleon scattering [35]. Meanwhile, the best limits for SD elastic WIMPproton (WIMPneutron) scattering come from IceCube [36] and LUX [35]. The IceCube limits depend on an assumption of DM wave annihilation in the sun exclusively to a single SM final state. As we will show in section 6, our DM annihilates in the wave to both and Higgsinos. Annihilation to Higgsinos could weaken the limits somewhat if the Higgsinos are stable, but that depends in detail on the other parameters of the model (such as , and ). Here we consider the simplest case where annihilation is only to ; this will provide the “worst case scenario” where the SD bound from IceCube is strongest. In section 7 we will also take into account annihilation to stable Higgsinos.
In this section, we will recast these constraints in terms of the couplings and . The discussion here can be viewed as an update of the nice treatment in [4] with the latest experimental results (in particular LUX). It is worth emphasizing that these bounds on and are quite model independent. Any WIMP DM that couples to SM mainly through Higgs and (including MSSM neutralinos) should satisfy these bounds.
To convert the results of these experiments into bounds on and , we first translate and into the couplings appearing in the effective Lagrangian for direct detection:
(3.1) 
In Higgs and portal DM models, the SI (SD) terms arise from Higgs () exchange, as shown in fig. 2. The coefficients of the effective operators are given in terms of and as:
(3.2) 
with being the Yukawa coupling and for downtype (uptype) quarks.
Then we use standard formulas that relate the DMnucleon cross sections to (see appendix A for our conventions and parameter choices). The result, assuming is given by:
(3.3) 
In principle, and are slightly different but the difference is negligible, so we only take to represent both.
The resulting limits on and are shown in fig. 3.^{6}^{6}6We agree with the limits from [4] after taking into account a factor of 2 in both and from 4component vs. 2component notation. We also agree with limits on operators from [48] modulo a factor of 4 between Dirac and Majorana fermions and a factor of a few difference between [48] and the latest LUX bounds. We do not agree with the limits on (related to our via ) reported in fig. 3 of [49]. Their limit on , derived from essentially the same LUX results, is over an order of magnitude weaker than ours. Amusingly, we note that although the constraint on the SI cross section is stronger than the SD crosssection, translated constraints on and are of the same order of magnitude. This is because the Higgsnucleon effective Yukawa coupling () is much weaker than the nucleon effective coupling (). Recall that the Higgsnucleon coupling is mainly due to Higgsgluongluon loopinduced interaction with heavy quarks running in the loop
(3.4) 
where is the number of heavy quarks and comes from QCD beta function at one loop. The second equality can be calculated using QCD scale anomaly that relates the QCD beta function to nucleon mass (see [50] for the original references).
4 Higgs mass and FineTuning
In this section we will describe our calculation of the Higgs mass in the model and its implication for the finetuning of the EW scale. As described in the introduction, we used SARAH 4.5.8 [44] and SPheno 3.3.7 [45, 46] to include all the loop corrections (contributions up to two loops both from the MSSM and the dark sector [51, 52]). Here we will describe an analytic treatment of the dominant oneloop contributions from and . This will serve as a valuable source of intuition, as well as a validation of the full twoloop numerical calculation (for more details on the validation, see appendix B).
The oneloop Higgs mass was previously computed in the literature using the ColemanWeinberg potential in closelyrelated vectorlike extensions of the MSSM [17, 18, 19, 21]. However, there are some key differences with our case that necessitate a fresh look. First, as noted above, in these past works, the vectorlike extension was Dirac, while ours is Majorana (the difference between and ). This leads to small differences in the formula for the Higgs mass. Second, previous works presented analytic formulas for the oneloop Higgs mass only in the simplified limit with common fermion masses (). Motivated by the DM side of the story, we will need the Higgs mass in a rather different regime, the mostlysinglet regime where .
Other effects that we will ignore in our discussion here, but that are taken into account in the full numerical SARAHSPheno calculation, include corrections, twoloop corrections, and the effective terms due to . The effects of are about a 1020% correction to , which amounts to a 2 GeV shift in . That matters for our calculations quantitively but not qualitatively. The values we consider in this paper motivated by naturalness are small enough that has a negligible effect on the Higgs mass. Finally, we are interested in moderatelylarge values of (e.g. ) but for simplicity we will present the limit here. The corrections due to also do not make a qualitative difference. (In particular, there are no blindspot cancellations here.)
With all of these simplifying assumptions, the result of our oneloop ColemanWeinberg calculation is:
(4.1) 
with
(4.2) 
where and . A plot of is shown in fig. 4 (left). We see that asymptotes to a finite value as or . In these limits (corresponding to mostlydoublet and mostlysinglet DM respectively), the dependence on the DM mass drops out, and is controlled by the ratio of the soft mass to the heavier mediator scale ( or respectively).
To raise the Higgs to 125 GeV in this paper, we rely on a combination of the extra vectorlike matter and MSSM stops. For stops at 1 TeV, which satisfy the current experimental bounds and imply about a % tuning of the EW VEV, the MSSM contribution to the Higgs mass is about 110 GeV (for a recent review see e.g. [3]). Therefore the target for from the dark sector is:
(4.3) 
This selects out a contour in the plane as shown in fig. 4 (left), according to the value of .
This has the following implications for the finetuning of the EW scale. Just as the dark sector lifts the physical Higgs mass analogous to stops in the MSSM, it also contributes to the finetuning of the EW scale through the renormalization of . Following [53, 54], we define the measure of finetuning to be:
(4.4) 
where is the running of due to the new fields
(4.5) 
Optimistically we take . We can combine this with (4.1) and (4.3) as follows. For a given value of and a given point in the plane, we can solve (4.3) for the soft mass . Then substituting this into (4.4), we get a value for . Regions of are shown in fig. 4 (right) for different representative values of . We see that we need to have any viable parameter space at all for a natural SUSY 125 GeV Higgs. This is not surprising, since from (4.1), we see that plays the role that plays for the MSSM stops. Of course, corrections we have neglected such as the terms and twoloop effects will modify this quantitatively. However, we will see that the same qualitative implications for finetuning and will persist in our final plots.
5 The need for mostlysinglet DM
In section 2, we derived formulas for and in terms of the parameters of the model, while in section 3 we showed that direct detection limits on and are at the level. Finally, in section 4, we argued that we need in order to have any viable parameter space for a natural SUSY Higgs at 125 GeV. Here we will combine these facts and show that the DM must be mostly singlet in order to be consistent with all the constraints.
Basically there are three possibilities: the welltempered regime where (recall our convention is that and are positive), the mostlydoublet regime where and , , and the mostlysinglet regime where and , . Keeping in mind that we need and large for a natural Higgs mass, the challenge is to decrease and to the level. In fact, alone is enough to rule out all but the mostlysinglet case. We will comment on the implications for in sections 6 and 7.
Examining the formula for (2.5), we see that for and , we have . (In particular, there is a cancellation in the denominator, leaving it ). This rules out the welltempered case.
The mostlydoublet case is ruled out separately by two independent considerations. First, from fig. 4, we see that in order to be natural and mostlydoublet, we must have the DM mass below GeV. However, we know by analogy with pure Higgsinos in the MSSM that the thermal relic density constraint requires TeV. (The mostlydoublet DM in this model has additional annihilation modes due to and , so will be even larger.) So the mostlydoublet scenario is not promising for naturalness.
Also, from direct detection, we are basically forced into the mostlysinglet regime. In order to lower by two orders of magnitude, we must either (a) raise or to increase the denominator of (2.5), or (b) cancel the two terms in the numerator of (2.5).

Increasing the denominator of (2.5) necessitates either or . In the former, corresponding to mostlysinglet DM, we see that and we can achieve the required level of suppression for TeV for and . Meanwhile for the latter, corresponding to mostlydoublet DM, we see that and therefore much larger TeV is required for and . The latter is greatly disfavored by naturalness (it would likely be as finetuned as 10 TeV stops in the MSSM).

Cancelling the two terms in the numerator requires
(5.1) This is the blind spot. Since and we are in the large limit, the RHS is generally much greater than one for any reasonable value of . Therefore we must be in the mostlysinglet DM regime to realize the blind spot.
We conclude that several different constraints independently point at mostlysinglet DM as the only viable possibility.
For later reference we exhibit and in the mostlysinglet limit
(5.2) 
Here we have taken and holding fixed and all the other mass scales. In fig. 5, we exhibit the amount of blind spot cancellation that is required by the SI bounds, for a typical choice of parameters that will lead to a viable relic density. We show this behavior by varying keeping other parameters fixed. We can see that we need only a very mild cancellation to satisfy the constraint on . Most of the suppression of is coming from large , which as we will see in the next section is fixed by the thermal relic abundance constraint.
In the same mostlysinglet regime, we also exhibit :
(5.3) 
As noted in the previous subsection, the Higgs mass in this limit to leading order does not depend at all on the DM mass . So the Higgs mass constraint to leading order in the mostlysinglet regime becomes a constraint on and . For example, according to (5.3), in order to achieve GeV for , we need .
6 DM annihilation in the mostlysinglet regime
An attractive feature of WIMP dark matter is its potential to naturally explain the observed relic abundance via the thermal freezeout mechanism. Following the usual procedure (see e.g. the classic review [55]), we have
(6.1) 
The integral over takes into account annihilation after freezeout, and parametrizes the freezeout temperature. is the thermallyaveraged DM annihilation cross section , summed over all final states and . This is usually expanded in the small velocity limit:
(6.2) 
where is a kinematic phase space factor. At the time of freezeout, the DM relative velocity is typically . Therefore, the annihilation cross section is generally controlled by the wave contributions , unless they are suppressed for some reason.
In our model, the dark matter has many interactions and annihilation channels that should all be considered in full generality. As described in the introduction, for numerical calculations we use micrOMEGAs 4.1.8 [47] source code generated by SARAH 4.5.8 [44] to accurately take these into account. However in the mostly singlet limit that we are interested in, the cross sections simplify and we can have an analytic understanding of the behaviour of our model. We will assume that DM is lighter than all MSSM superpartners except possibly the Higgsinos, which are forced to be light by naturalness. In this case, the freezeout process happens only through annihilation to SM particles and the Higgsinos. Including the Higgsinos in the story is a major difference from simplifiedmodelanalyses of singletdoublet dark matter, which generally just add the singlet and doublets to the SM. As we will see, the Higgsinos can be a major part of the DM annihilation in the early universe.
The full cross sections are too complicated to print here. Instead, we will expand in the mostlysinglet limit , , with the further assumption that . This suffices for our purposes and results in relatively simple expressions. (One exception is the treelevel, wave cross section in the next subsection, for which we can write down an extremely simple exact expression in terms of .)
6.1 DM annihilation to fermions
The fermions have wave contributions
(6.3) 
where , and was defined in section 4. In the second line, we have summed over the various Higgsino final states including both neutralinos and charginos, assuming a pure MSSM Higgsino (i.e. decoupled). The fermion coefficients are always subdominant (suppressed by both and ), so we have not included them here.
The fermion cross sections are all suppressed by the square of the fermion mass, so and Higgsinos are the dominant channels. This is the famous wave helicity suppression of DM annihilation to fermion pairs.
Although and Higgsinos are parametrically similar, their diagrammatic origin is entirely different. The former (latter) arise from channel (channel superpartner) exchange. As a result, the Higgsinos are suppressed by the soft mass . For , we saw in section 5 that we need for GeV, so the suppression is not large. Also, is constrained to be GeV by naturalness. So all in all, the Higgsino contribution ends up generally of the same order or smaller than .
The fact that the SM fermions all arise from channel diagrams means that they have a simple exact expression beyond the small approximation:
(6.4) 
In other words, controls both the SD direct detection cross section and the annihilation to . Therefore, we expect to see a fairly direct correlation between the SD direct detection limits and the relic density constraint.
6.2 DM annihilation to bosons
Meanwhile the diboson cross sections are all wave to leading order:
(6.5) 
Here we took for simplicity; we checked that the corrections are irrelevant. The wave contributions are suppressed by so they are always subdominant to the wave contributions shown here.
Clearly, the diboson cross sections exhibit some interesting features. They are nonvanishing even in the limit, so they can be understood as a consequence of symmetry. These treelevel cross sections arise entirely due to the longitudinal components of the and bosons, which by the Goldstone equivalence theorem are also equivalent to the charged and neutral Goldstones and respectively. Under a rotation, and , while under an rotation, . This explains both relations in (6.5).
Comparing and Higgsinos to the total diboson cross section, we see that parametrically the latter can be larger than the former, for sufficiently large . However the cross over point is generally at very large and . For instance, for and , we find the cross over to be in the range TeV for GeV. This is well beyond the naturalnessmotivated part of the parameter space that we are focusing on in this paper. Therefore we conclude that the total is always dominated by and Higgsinos, and dibosons are always a subdominant part of it.
6.3 Total annihilation cross section
We have shown analytically that the relic density is dominated by wave annihilation to and Higgsinos (assuming of course that the DM is above the respective thresholds):
(6.6) 
A plot comparing our analytics to micrOMEGAs is shown in fig. 6 for fixed choices of the parameters; we see there is excellent agreement across the entire range of relevant DM masses. We confirm that the dibosons are never more than % of the relic density across the entire parameter range of interest. Higgsinos and are comparable for GeV, while for GeV, dominates, as expected from the dependence of the Higgsino cross section (6.3).
One very interesting consequence of (6.6) is that in the limit of large , the DM mass drops out of the annihilation cross section. Furthermore, we have seen that we need for the blind spot, for the Higgs mass, and for naturalness. Thus the WIMP miracle transforms from being a constraint on the WIMP mass to being a constraint on the mediator scale ! This helps to relieve the “WIMP little hierarchy problem”, whereby the preference of the thermal relic constraint for TeVscale WIMPs is in tension with naturalness. Comparing with (5.2), we also expect that the relic density constraint will essentially fix to a unique value. We will confirm this in the next section with our full numerical scans and discuss its implications for SD direct detection.
7 Putting it all together
7.1 Plots in the  plane
Having described the various individual components of the analysis of the model (direct detection, the Higgs mass, and the relic abundance), we will now combine them and describe how the different constraints interact to produce the viable parameter space of the model.
In fig. 7 we show contour plots for numerical scans over the () plane for fixed values of and . We choose four sets of benchmark parameters: large coupling (, ) and small coupling (, ); and large ( GeV) and small ( GeV).
We see the impact of the direct detection limits on the parameter space of the model. The LUX SI and SD limits are strongest almost everywhere except a tiny sliver for large in the , case where IceCube has an impact. (Note that the LUX limits assume the singletdoublet sector comprises all of the DM, regardless of whether it is thermal or not.) The SD (SI) limits primarily cover the lighter (heavier) DM mass region. The heavier DM region is ruled out because we are holding fixed , so as one increases the blind spot cancellation shown in (5.2) becomes less effective.
For every point in the plane we numerically solved (using SPheno) the GeV constraint for the common soft mass ; these contours are shown in fig. 7 along with their corresponding tuning. These contours are mostly vertical; as discussed in section 4, the soft mass and depend primarily on since the dependence on drops out to leading order at large .
Finally, we used micrOMEGAs to numerically solve the thermal relic density constraint [56]; this fixes as a function of and these contours are shown in green for various choices of the parameters. Note the rapid increase in across the top and Higgsino thresholds. Here new wave annihilation channels open up, and so larger values of are needed to maintain the overall annihilation rate at the thermal relic value. This effect is more pronounced for larger values of and for larger values of . Indeed, in section 6 we saw that the annihilation cross sections to and Higgsinos are enhanced for greater , and the Higgsino cross section in particular is proportional to .
Since larger decreases direct detection cross sections, increasing and also increases the viable parameter space for thermal relic DM. The Higgsino channels in particular allow the model to survive the direct detection limits over a wider range of parameter space than would have been the case for nonsupersymmetric singletdoublet DM. Fig. 7 also shows that larger is better for finetuning, confirming our discussion in section 4. The only potential drawback of the larger coupling choice is (as we will discuss in section 8.1) that the former has a lower Landau pole ( TeV vs TeV).
Away from the top and Higgsino thresholds, we see that the relic density contours are mostly vertical, meaning that the relic density constraint becomes a constraint primarily on , once the other parameters (, , , ) are fixed, i.e. the WIMP DM mass drops out to leading order. This confirms our analytics in the previous section.
7.2 Projecting onto the thermal relic contour
Finally let us impose the relic density constraint and see how various parameters vary along the green contours in fig. 7. In fig. 8 we show the finetuning for the points with the correct relic abundance. It is remarkable that there are allowed regions of the parameter space with , making this model much less tuned than the MSSM.
In fig. 9 we show for the points of the parameter space that satisfy (including both allowed and excluded points from direct detection limits). We see that varying we can move toward the blind spot and satisfy the SI direct detection bounds.
Similarly, in fig. 10 we show for the points of the parameter space with (including both allowed and excluded points from direct detection limits). We can see that for contours of constant have an approximately constant . This confirms the discussion based on analytics in section 6.3. Indeed, using (6.6) with the parameter choices here, we find that for , ranges from . These values are clearly illustrated in fig. 10.
We conclude that (and consequently ) is basically fixed by the relic density constraint. Requiring to be all the dark matter leads to a nearly unique prediction for the SD cross section! Fortunately, as shown in fig. 10, these values of are still allowed by the current direct detection experiments, IceCube in particular.^{7}^{7}7It is quite crucial that our DM annihilates almost exclusively to and Higgsinos. The IceCube bound on is by a factor of a few weaker than the cross section, and it saves the model from being already ruled out (stable Higgsinos would not contribute to DM detection in IceCube). With factor of 10100 improvements in cross section expected from Xenon1T [37] and LZ [38], the next generation of DM direct detection experiments will be sensitive to essentially the entire parameter space of this model (assuming is a thermal relic and is all the DM). A discovery might be right around the corner!
8 Outlook
In this section we briefly discuss the UV behavior of the model (in particular the Landau poles) and the potential sensitivity from LHC Run II. Finally we conclude with some thoughts on future directions.
8.1 UV considerations
So far we have been exploring our model at the EW scale and have identified the interesting parts of the parameter space around . Here we want to examine the UV consequences of such large Yukawa couplings and comment on possible solutions to the Landau pole problem.
Let’s focus on the most important couplings, and neglect the effect of the other couplings in finding the scale of Landau poles. Starting from one loop beta functions above the scale of the new fields (including spectator color triplets for unification) we have
(8.1) 
Solving the RGE’s of our model numerically, we can find the lowest scale at which one of the couplings hits its Landau pole. In fig. 11 we show this scale as a function of the Yukawa couplings at TeV.
Note that for (as we have considered in this work), the Landau poles are above TeV. Now we might ask: how can we understand physics above the Landau pole scale, or how can we postpone it to higher energies e.g. the GUT scale? One idea is to use nonAbelian gauge interactions for the new sector to reduce the beta functions of the Yukawa couplings: if we include multiple copies of , , and charge them under a nonAbelian gauge group, the corresponding gauge coupling appears with negative sign in the beta function of (see e.g. [32] for a recent implementation of this idea). As is Majorana, we need to be in a real representation of the new gauge group. A simple example is when the gauge group is and is in the fundamental representation. Another possibility might be to match our model to the magnetic side of a Seiberg duality and interpret physics above the scale of the Landau pole by the electric theory. It will be interesting to explore these ideas further in the future.
8.2 LHC Phenomenology
In addition to direct detection experiments, DM models are also probed by the LHC. In principle, monojet+MET [40, 41] and monophoton+MET [42, 43] searches for direct DM production could be sensitive to our model. Since quarks and gluons only talk to through channel diagrams involving ’s and Higgses, these searches constrain the same and couplings as direct detection. However, these constraints are weaker by several orders of magnitude than those from direct detection under the assumption that our DM candidate is all of the relic density, for the mass range we consider. See e.g. [57, 58] for a recent discussion in terms of simplified DM models.
Instead, let us briefly consider mono(, , )+MET. This can occur in our model through production of and and subsequent decay of the (mostlydoublet) and . A full treatment including estimation of SM backgrounds, detector acceptances, etc. is beyond the scope of this work. Here we will just present the raw production cross sections in our model.
Diagrams contributing to monoHiggs/+MET are shown in fig. 12 (mono+MET is the same as monoHiggs with the final state Higgs replaced by ). Note that we have included the oneloop gluon fusion diagram.^{8}^{8}8We thank Matt Reece for bringing this to our attention. Because of the large, Yukawas , in this model, this contribution can be as much as 60% of the total cross section. We calculated the gluon fusion contribution analytically, and the tree level contributions both analytically and with MadGraph5 [59] using the model file generated by SARAH and the spectrum files generated by SPheno. More details on the analytics are given in appendix C. In both cases, we used the NNPDF2.3 [60] PDF set. Fig. 13 shows the sum of tree level contributions and gluon fusion along the contour. We see that LHC13 will ultimately be able to probe the small mass region. Of course, if is all of the dark matter, then direct detection experiments will discover the model first. In that case, the LHC will only be useful as a postdiscovery confirmation of the model. However, since the LHC is producing directly, it does not depend on the relic density. Therefore if our dark sector is only one component of , the direct detection limits could be greatly relaxed while the LHC would remain sensitive.
8.3 Future directions
The work presented in this paper is a simple realization of a general idea: economically extending the MSSM with a single sector that provides both thermal WIMP dark matter and the 125 GeV Higgs mass. Here we took this sector to be a singlet and a pair of doublets, but one could easily imagine many other possibilities. For instance, very popular ideas for lifting the Higgs mass include the NMSSM (see e.g. [61] for a review and original references) and nondecoupling Dterms [62, 63]. While dark matter in the NMSSM is a wellstudied topic, it would be very interesting to try to connect nondecoupling Dterms to dark matter.
Even within the context of our specific singletdoublet model, there are many interesting open questions. In this work we made some simplifying assumptions in our analysis, and it would be interesting to explore the consequences of relaxing these assumptions. For example, we took all model parameters to be real, but in general there is one physical CPviolating phase. The effect of this phase on direct detection and annihilation cross sections can qualitatively change the model’s behavior. Furthermore, we took all the soft masssquareds to be positive to increase the Higgs mass. One might wonder how the phenomenology of the model would change if one of the soft masses is negative and the DM is a scalar instead of a fermion. We also assumed negligible terms in the dark sector. By analogy to stops, having substantial terms can help in raising the Higgs mass, see e.g. [19]. This could allow for smaller , and open up more of the parameter space. Additionally, we focused on dark matter above GeV. It could be interesting to study the phenomenology of the model for lighter dark matter masses. In particular the annihilations through the Higgs and resonances could be large enough while still having suppressed direct detection signals. Finally, one could relax the assumption that is thermal and is all of the DM, and consider nonthermal relics or multicomponent DM scenarios. All of these directions will become especially wellmotivated if nothing is discovered at the next round of direct detection experiments, as discussed in section 7.
There are also many interesting modelbuilding directions in the UV. For example, enlarging the dark sector to accommodate a nonAbelian gauge symmetry could have potentially interesting consequences. As noted in section 8.1, this may help postpone the Landau pole of the Yukawa couplings, and the new gauge interactions could play an important role in the dynamics of the dark sector. Additionally we have two supersymmetric masses and at the electroweak scale. Perhaps the same dynamics that generates in the MSSM is responsible for generating these masses as well.
Acknowledgements:
We thank Kyu Jung Bae, Matt Buckley, John Paul Chou, Eva Halkiadakis, Ryuichiro Kitano, Stephen Martin, Natsumi Nagata, Matt Reece, Alessandro Strumia and Scott Thomas for helpful discussions. We are especially grateful to Florian Staub for his indispensable help setting up the model in SARAH. We are also grateful to Natsumi Nagata and Florian Staub for comments on the manuscript. This work was supported by DOE grant DESC0013678.
Appendix A Connecting model parameters to DD cross sections
In this appendix, we will review how to relate the SI and SD DMnucleon cross sections to the couplings and appearing in the effective Lagrangian (3.1). To check our results we verify that by calculating SI and SD cross sections analytically, we get the same result as the one we get from micrOMEGAs.
Following [55], the SI and SD cross sections are