Supersymmetric dark matter after LHC run 1.

Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of the lightest neutralino, [Formula: see text], assumed here to be the lightest SUSY particle (LSP) and thus the dark matter (DM) particle, into the range allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly degenerate next-to-lightest supersymmetric particle such as the lighter stau [Formula: see text], stop [Formula: see text] or chargino [Formula: see text], resonant annihilation via direct-channel heavy Higgs bosons H / A, the light Higgs boson h or the Z boson, and enhanced annihilation via a larger Higgsino component of the LSP in the focus-point region. These mechanisms typically select lower-dimensional subspaces in MSSM scenarios such as the CMSSM, NUHM1, NUHM2, and pMSSM10. We analyze how future LHC and direct DM searches can complement each other in the exploration of the different DM mechanisms within these scenarios. We find that the [Formula: see text] coannihilation regions of the CMSSM, NUHM1, NUHM2 can largely be explored at the LHC via searches for [Formula: see text] events and long-lived charged particles, whereas their H / A funnel, focus-point and [Formula: see text] coannihilation regions can largely be explored by the LZ and Darwin DM direct detection experiments. We find that the dominant DM mechanism in our pMSSM10 analysis is [Formula: see text] coannihilation: parts of its parameter space can be explored by the LHC, and a larger portion by future direct DM searches.

Introduction

Fig. 1

The planes in the CMSSM (upper left), the NUHM1 (upper right) and the NUHM2 (lower left), and the plane in the pMSSM10. Regions in which different mechanisms bring the CDM density into the allowed range are shaded as described in the legend and discussed in the text. The red and blue contours are the and 5.99 contours found in global fits to these models, corresponding approximately to the 68 and 95 % CL contours, with the green stars indicating the best-fit points, and the solid purple contours show the current LHC 95 % exclusions from  searches. In the CMSSM, NUHM1, and NUHM2 cases, the dashed purple contours show the prospective 5 discovery reaches for  searches at the LHC with 3000/fb at 14 TeV, corresponding approximately to the 95 % CL exclusion sensitivity with 300/fb at 14 TeV. In the pMSSM10 case, the dashed purple contour shows the 95 % CL exclusion sensitivity of the LHC with 3000/fb assuming , and the dash-dotted lines bound the corresponding sensitivity region assuming

Fig. 3

The planes in the CMSSM (upper left), the NUHM1 (upper right) and the NUHM2 (lower left), showing the regions where the lowest- points have and  s  s: the lifetimes [14] of these points are color-coded, as indicated in the legends. The red and blue contours are for relative to the absolute minimum. Also shown in these panels as solid purple contours are the current LHC 95 % exclusions from  searches in the coannihilation regions, and as dashed purple contours the prospective 5 discovery reaches for  searches at the LHC with 3000/fb at 14 TeV, corresponding approximately to the 95 % CL exclusion sensitivity with 300/fb at 14 TeV. As discussed in the text, the sensitivities of LHC searches for metastable ’s in the coannihilation region are expected to be similar [98]. The lower right panel shows the one-dimensional function in the pMSSM10 for the lifetime of the in the range  s  s

Fig. 4

The planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). The red and blue solid lines are the and 5.99 contours, and the solid (dashed) purple lines are the current and (projected) 95 % exclusion contours for  searches at the LHC (with 300/fb of data at 14 TeV). The solid lines are almost identical with the contours for 5 discovery with 3000/fb

Fig. 8

The planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). The red and blue solid lines are the and 5.99 contours, and the solid purple lines show the projected 95 % exclusion sensitivity of the LUX-Zeplin (LZ) experiment [124]. The green and black lines show the current sensitivities of the XENON100 [122] and LUX [123] experiments, respectively, and the dashed orange line shows the astrophysical neutrino ‘floor’ [126], below which astrophysical neutrino backgrounds dominate (yellow region)

The density of cold dark matter (CDM) in the Universe is now very tightly constrained, in particular by measurements of the cosmic microwave background radiation, which yield [1] and are consistent with other, less precise, determinations. This determination of the CDM density at the percent level imposes a corresponding constraint on the parameters of any model that provides the dominant fraction of the CDM density. This is, in particular, true for supersymmetric (SUSY) models with conserved R-parity in which the CDM is provided by the stable lightest SUSY particle (LSP) [2, 3]. In a series of recent papers incorporating the data from LHC run 1 and elsewhere, we have implemented the dark matter (DM) density constraint in global analyses of the parameter spaces of different variants of the minimal SUSY extension of the Standard Model (MSSM), assuming that the LSP is the lightest neutralino . The models studied included the constrained MSSM (CMSSM) with universal soft SUSY-breaking parameters ( and , in standard notation) at the GUT scale [4], the NUHM1(2) in which universality is relaxed for both together (each separately) of the soft SUSY-breaking contributions to the masses-squared of the Higgs multiplets [4, 5], and a version of the pMSSM10 [6], in which 10 of the effective Lagrangian parameters (3 gaugino masses , 2 squark masses , a common slepton mass , a common trilinear coupling , the Higgs mixing parameter , the pseudoscalar Higgs mass , and the ratio of Higgs vevs ) are treated as independent inputs specified at the electroweak scale.

Reproducing correctly the cosmological CDM density requires, in general, some special choice of the SUSY model parameters, which may be some particular combination of sparticle masses and/or couplings. Examples of the former include hypersurfaces in the SUSY parameter space where the LSP is almost degenerate in mass with some next-to-lightest SUSY particle (NLSP), such as the lighter stau [7-14], stop [15-21] or chargino [22-26], or where is almost half the mass of a boson such as a heavy Higgs H / A [27-31], a light Higgs h or Z [32, 33], in which case rapid direct-channel annihilation may bring the CDM density into the allowed range. Examples of special coupling combinations include the focus-point region [34-39], where the LSP acquires a significant Higgsino component.

We have commented in our previous work on the relevances of these DM mechanisms for our global analyses. Here we discuss systematically which DM mechanisms are dominant in which subspaces of the CMSSM [40-51], NUHM1 [52-55], NUHM2 [54-57], and pMSSM10 (see, for example, [58-72]) parameter spaces, what the corresponding experimental signatures are, and how one might discover SUSY in each of these different DM regions.

Our analysis of the possible detectability of supersymmetry in the CMSSM, NUHM1, NUHM2, and pMSSM10, depending on the dominant DM mechanisms, is summarized in Table 1.

Table 1

Summary of the detectability of supersymmetry in the CMSSM, NUHM1, NUHM2, and pMSSM10 models at the LHC in searches for  events, long-lived charged particles (LL) and heavy A / H Higgs bosons, and in direct DM search experiments, depending on the dominant mechanism for bringing the DM density into the cosmological range. The symbols , () and indicate good prospects, interesting possibilities and poorer prospects, respectively. The symbol – indicates that a DM mechanism is not important for the corresponding model. The LHC information is drawn largely from Figs. 1, 3 and 4, and the direct DM search information from Fig. 8

DM mechanism	Exp’t	Models
		CMSSM	NUHM1	NUHM2	pMSSM10
τ~_1	LHC	✓  /E_T , ✓ LL	(✓  /E_T , ✓ LL)	(✓  /E_T , ✓ LL)	(✓  /E_T ), × LL
Coann.	DM	(✓)	(✓)	×	×
χ~1±	LHC	–	×	×	(✓  /E_T )
Coann.	DM	–	✓	✓	(✓)
t~_1	LHC	–	–	✓  /E_T	–
Coann.	DM	–	–	×	–
A / H	LHC	✓ A / H	(✓ A / H)	(✓ A / H)	–
Funnel	DM	✓	✓	(✓)	–
Focus	LHC	(✓  /E_T )	–	–	–
Point	DM	✓	–	–	–
h, Z	LHC	–	–	–	(✓  /E_T )
Funnels	DM	–	–	–	(✓)

Measures of mass degeneracy

We first introduce measures on the MSSM parameters that quantify the relevant mass degeneracies and define each of the above-mentioned subspaces in the CMSSM, NUHM1 and NUHM2 [5, 73]:In each case we also indicate the color coding we use in the subsequent figures. The measures (1) that we use are empirical, but we have verified extensively that CMSSM, NUHM1, and NUHM2 points that satisfy the DM density constraint do fulfill at least one of these conditions, and that they indeed correspond to the dominant DM mechanisms (in the sense of giving the largest fractions of final states, generally  %) [5, 73]. We have found that there are some ‘hybrid’ regions where the dominant mechanism requires two of these conditions simultaneously. In particular, there are regions where the dominant DM mechanism is or , processes involving both stau coannihilation and annihilation via the A / H funnel, which we color purple. There are also regions where the chargino coannihilation condition is satisfied simultaneously with the stau coannihilation or A / H funnel condition. However, a dedicated study using MicrOMEGAs [74] shows that chargino coannihilation is the dominant DM mechanism in these regions, and that hybrid processes dependent on the and some other degeneracy conditions being valid simultaneously are unimportant, so we color these regions the same green as the other chargino coannihilation regions.

The above DM mechanism conditions need to be modified for our analysis of the pMSSM10. First, as we shall see later, funnels due to annihilations via direct-channel h and Z resonances can be important [32, 33], so for the pMSSM10 we add to (1) the supplementary criteria:

Second, we find that chargino coannihilation dominates in the pMSSM10 also when the second condition in (1) is relaxed: we use later the condition in our subsequent analysis, which reproduces better the domains of dominance by coannihilation.1 Finally, we recall that the focusing property of the RGEs is not relevant in the pMSSM10. However, the LSP annihilation rate may still be enhanced when the fifth measure in (1) is satisfied, due to a larger Higgsino component in the LSP, though we find that the dominant DM mechanism in the pMSSM10 generally does not involve this property. We use the same cyan color to identify regions where this condition is satisfied, though it is not due to focus-point behavior.

Our discussion here of DM mechanisms is based on our previously published global likelihood analyses of the CMSSM and NUHM1 [4], the NUHM2 [5] and the pMSSM10 [6]. The reader wishing to know details of our treatments of the various experimental, phenomenological, theoretical and cosmological constraints, as well as our strategies for scanning the parameter spaces of these models is referred to [4–6, 75–77]. However, we note here a couple of important points.

An important contribution to the global function comes from , which we evaluate using FeynHiggs 2.10.0 [78]. This combines Feynman-diagrammatic results with a resummation of leading and subleading logarithmic corrections from the stop/top sector. As documented in [78], this is a significant improvement on earlier versions of FeynHiggs, and differences with other earlier codes and their implications are discussed in [79]. FeynHiggs 2.10.0 also provides an estimate of the uncertainty in the calculation of . This is in general well approximated by (thus assigning a uncertainty at the 95 % CL), which we use together with the current experimental value of to calculate the contribution to the global function. Another important contribution to our calculation of the global function is made by the experimental measurement of  [80, 81], which differs by from many theoretical calculations. This is in tension with LHC constraints on the CMSSM, NUHM1, and NUHM2 [4, 5], but not the pMSSM10 [6].

Dominant dark matter mechanisms

In this section we discuss the various mechanisms that play dominant rôles in bringing the relic density into the experimentally measured interval in our four models. We display in Fig. 1 planes for the CMSSM (upper left) [4], the NUHM1 (upper right) [4] and the NUHM2 (lower left) [5], while for the pMSSM10 we show the plane (lower right), where denotes the mass of the squarks of the first two generations, which we assume to be common.2 The and 5.99 contours that we found in global fits to these models, corresponding approximately to the 68 and 95 % CL contours, are shown as solid red and blue lines, respectively. Here and elsewhere, the green stars indicate the best-fit points, whose exact locations in some parameter planes are poorly determined and do not carry much useful information, in general, as the minima are quite shallow. Also shown, as solid purple lines, is the current 95 % CL CMSSM exclusion from  searches at the LHC.3 The dashed purple contours in the CMSSM, NUHM1, and NUHM2 cases show the prospective 5 discovery reaches for  searches at the LHC with 3000/fb at 14 TeV, corresponding approximately to the 95 % CL exclusion sensitivity with 300/fb at 14 TeV. In the pMSSM10 case the dashed purple contour shows the 95 % CL exclusion sensitivity of the LHC with 3000/fb assuming , and the dash-dotted lines bound the corresponding sensitivity region assuming .

All the parameter planes we show are color-coded as listed in (1, 2), with regions where none of these processes are dominant left uncolored. We see in the upper left panel of Fig. 1 that three DM mechanisms dominate in the CMSSM: coannihilation at low , the H / A funnel at larger and , and the focus point at larger and smaller where the neutralino becomes a ‘well-tempered’ mixture of bino and Higgsino [82]. There is also a hybrid region extending up to . In the case of the NUHM1 shown in the upper right panel of Fig. 1, there is an analogous coannihilation region. However, it has a much larger hybrid region, which has an extension to low for . On the other hand, coannihilation dominates in a large region with . Here , so that the and are nearly degenerate in mass and the has mainly a Higgsino composition. A similar coannihilation region is visible in the NUHM2 in the lower left panel of Fig. 1, where we also see a more extensive coannihilation region extending to ,4 with a relatively small hybrid region. This is the only case where we see a region of dominance by coannihilation, in islands around .

In contrast, as shown in the lower right panel of Fig. 1, we found in our version of the pMSSM10 [6] that the dominant DM mechanism is usually coannihilation, this time with a Bino-like LSP and .5 We also note in Fig. 1 bands at low where rapid annihilations via the h and Z funnels are dominant. Not shown in Fig. 1 are scatterings of points with where coannihilation can also be important, and of points with where the fifth condition in (1) comes into play.6 We see later that this condition and the coannihilation mechanism dominate in specific regions of other projections of the pMSSM10 parameter space, as does ‘bulk’ annihilation where none of the conditions (1) and ( 2) are satisfied.

One may also consider the possibility that the LSP provides only a fraction of the CDM. A complete discussion of this possibility is beyond the scope of this paper, but we note that in some regions, e.g., those dominated by or coannihilation, lowering the CDM density requires reducing the NLSP–LSP mass difference. This would also have the effect of reducing correspondingly the maximal LSP mass, which would favor sparticle detection at the LHC. On the other hand, direct detection would be more difficult if only a small fraction of the galactic halos were composed of LSPs. For some recent discussions, see [83-85].

Before concluding this section, we note that the marginalization procedure we use to produce the two-dimensional planes we display does not display directly all the physics information that would be contained in a complete set of two-dimensional slices with fixed parameters such as and used in many analyses, cf., [54, 55, 79]. On the other hand, these carry very little if any statistical information about the global parameter spaces. The color coding of the dominant annihilation channels used here contains the information of interest to us in this paper.

The LHC sensitivity

In this section we discuss the prospective reaches of future LHC searches (see for example, [86, 87]) and their impacts in the contexts of the various preferred DM mechanisms. We will see that in many cases the preferred DM mechanism can be directly probed via the appropriate LHC searches, as summarized in Table 1.

searches

Looking now at the physics reach of the LHC with  searches (we base our analysis on Aad et al. [88, 89]), we see that in the CMSSM the preferred stau coannihilation region lies just outside the current LHC 95 % CL exclusion region (solid purple line). On the other hand, the A / H-funnel region allowed in the CMSSM at the 95 % CL lies well outside this region. However, this is no longer the case in the NUHM1 and particularly the NUHM2, where portions of the A / H-funnel region lie much closer to the LHC 95 % CL exclusion. This is possible mainly because the Higgs mass constraint is less restrictive in these models. In the CMSSM, the future LHC sensitivity estimated in Fig. 1 of [90] (dashed purple line) covers the region where coannihilation is dominant, and a part of the hybrid region. It also covers slices of the H / A funnel region and of the focus-point region. Similar features are seen in the NUHM1, except that the focus-point region is not visible in this case, but the LHC  search has no sensitivity in the coannihilation region. In the case of the NUHM2, the  search is sensitive to only part of the coannihilation, and none of the coannihilation region, but it does cover part of the H / A funnel region and all the coannihilation region. In the pMSSM10 case, parts of the coannihilation region and the low- band are accessible to future LHC searches,7 though the future LHC  searches are less sensitive if (dashed line) than if (dotted line). Table 1 summarizes the observability of events that we estimate in the different models considered, depending on the dominant DM mechanism in each case.

The possibility of a long-lived charged sparticle

In some circumstances, a charged sparticle such as the or the may be the NLSP, and be only slightly more massive than the LSP, so that it can in principle be long-lived. As we now show, the presence of a long-lived is indeed a distinctive prospective signature in the CMSSM, NUHM1, and NUHM2.

Figure 2 displays on log-linear scales the one-dimensional profile likelihoods as functions of the mass differences in the coannihilation regions in the CMSSM (upper left), the NUHM1 (upper right) and the NUHM2 (lower left). Also shown (in the lower right panel) is the corresponding distribution in the pMSSM10, although in this case the coannihilation region is disfavored: in this case we see that the likelihood function increases sharply for , rising to for very small mass differences.

Fig. 2

The one-dimensional profile likelihood functions in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right) for . In the CMSSM, NUHM1, and NUHM2, low values of are found for points with , whereas in the pMSSM10 rises to at small

In the CMSSM, NUHM1, and NUHM2 panels the one-dimensional profile likelihood function is quite a flat function of , and there are points with that have very small values of this mass difference . Hence, it is possible that the might live long enough ( ns) to appear at the LHC as a long-lived (LL) charged particle, or even long enough ( s) to affect Big Bang nucleosynthesis [91-97]. However, we emphasize that the mass differences required to realize these possibilities () require quite special parameter sets. One should presumably require  s in order to avoid destroying the success of Big Bang nucleosynthesis, a constraint that we impose in the following figures.

If , corresponding to  ns, the would live long enough to be detectable at the LHC as a LL charged particle [98]. Figure 3 displays the regions of the plane in the CMSSM (upper left panel), the NUHM1 (upper right panel) and the NUHM2 (lower left panel) where the lowest- point has  s  s: the lifetimes of these points are color-coded, as indicated in the legends. The contours for relative to the absolute minimum of our data set are shown as solid red and blue lines, respectively.

On the other hand, the lower right panel of Fig. 3 displays the one-dimensional function in the pMSSM10 for the lifetime of the in the range  s  s. We see that throughout the displayed range, indicating that a long-lived is not expected in the region of the pMSSM10 parameter space that is favored by present data. This is because in our analysis the measurement and the DM constraint favor light sfermions and higgsinos, whereas long-lived charginos typically require sfermion and higgsino masses larger than 10 TeV [99].8

The sensitivity of the LHC to a long-lived has been compared in [98] to the sensitivity to  events, and found to be comparable within the uncertainties. We therefore assume that the projected sensitivity of the LHC to  events is a good approximation to its sensitivity to parameter sets in the coannihilation region with . However, while the reach in the parameter space is similar, the long-lived stau would constitute a spectacular additional signature, and would give complementary information to the direct searches for colored sparticles. The purple contours in the CMSSM, NUHM1, and NUHM2 panels of Fig. 3 again show the present and prospective reaches of the LHC for such events. We infer that long-lived charged particles could be detectable at the LHC throughout the lower parts of the allowed regions in the CMSSM, NUHM1 and NUHM2. Table 1 also summarizes the observability of long-lived charged sparticles in the different models considered.

Squark and gluino searches

Figure 4 displays the planes for the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). In each panel we show the and 5.99 contours as red and blue solid lines, respectively. The current 95 % CL exclusions from ATLAS  searches are shown as solid purple lines, and the estimated reaches of  searches for 95 % exclusion with 300/fb of data at 14 TeV [100] (very similar to the reaches for 5 discovery with 3000/fb) are shown as dashed purple lines. The CMSSM panel shows again that the coannihilation region is within the LHC reach in this model. However, in the NUHM1 and the NUHM2 only portions of the coannihilation regions are accessible at the LHC, along with small pieces of the H / A funnel regions. In the case of the NUHM2 the small coannihilation regions are also well within the LHC reach.

The pMSSM10 panel shows a completely different picture: coannihilation dominates throughout the plane, as discussed at the end of Sect. 3, and the likelihood function is very flat beyond the current LHC  exclusion. The LHC at 14 TeV will explore a large part of the plane, but a (more) complete exploration would be a task for a higher-energy collider [90].

Charginos and neutralinos

The differences between the dominant DM mechanisms in the pMSSM10 and the other models studied are highlighted in Fig. 5, which displays the planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). The diagonal dashed lines indicate where . As usual, the red and blue solid lines are the and 5.99 contours. In the pMSSM10 case, the region preferred at the 68 % CL is a narrow strip where is small, whereas in the other models much of the 68 % CL region is in a narrow strip where . This reflects the fact that in the CMSSM, NUHM1, and NUHM2 universal boundary conditions are imposed on the gaugino masses at the GUT scale.

Fig. 5

The planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). The red and blue solid lines are the and 5.99 contours. The solid (dashed) orange lines are the current and projected 3000/fb 95 % CL exclusion sensitivities for searches [101], the green dashed lines the projected 3000/fb 95 % CL exclusion sensitivity for a search [102] (both from [103]), and the magenta dashed line is the projected 3000/fb 95 % CL exclusion sensitivity for searches (from [104])

We see that coannihilation dominates over most of the 95 % CL region in this projection of the pMSSM10 parameter space, though not in the 68 % CL region, which has small and is where coannihilation dominates.9 On the other hand, in the CMSSM, NUHM1, and NUHM2, the H / A funnel dominates most of the 95 % CL regions in the planes, and also the 68 % CL region in the CMSSM, whereas coannihilation dominates the 68 % CL region in the NUHM1 and the NUHM2, as shown in Fig. 1. At the 95 % CL there are also small coannihilation regions in the CMSSM, NUHM1, and NUHM2 where is small, and the CMSSM and NUHM1 also have small focus-point regions.

Figure 5 also displays the present and prospective future sensitivities of various LHC searches in the planes. The solid (dashed) orange lines are the current and projected 3000/fb 95 % CL exclusion sensitivities for searches [101],10 the green dashed lines show the projected 3000/fb 95 % CL exclusion sensitivities for the search [102], both taken from [103], and the magenta dashed line is the projected 3000/fb 95 % CL exclusion sensitivity for searches, taken from [104]. We see that these searches have very limited sensitivities to the CMSSM, NUHM1, and NUHM2, but could explore significant parts of the coannihilation region in the pMSSM10. However, they would largely miss the coannihilation region, i.e., the dominant pMSSM10 DM mechanism cannot be explored by direct searches at the LHC.

The lighter stop squark

Next we study the differences in the impacts of the dominant DM mechanisms on the pMSSM10 and the other models in the planes shown in Fig. 6. In each of these planes, we indicate by dashed lines where and where . In the CMSSM, as shown in the upper left panel, we see that the coannihilation region (which contains all of the parameter space that is allowed at the 68 % CL) is well separated from the H / A funnel region, and that only a small part of the displayed portion of the plane is in the hybrid region. In this model we find that at the 95 % CL, and we do not find a coannihilation region, but we do see a focus-point region and a small coannihilation region. The situation in the NUHM1 (upper right panel of Fig. 6) exhibits significant differences. The coannihilation region (which again dominates the 68 % CL region) and the H / A funnel region still dominate the displayed portion of the plane, but there is a larger hybrid region, the focus-point region has disappeared and the coannihilation region has remained small, but has moved to larger . We also note the appearance of a small coannihilation ‘island’ at the 95 % CL in this model. In the case of the NUHM2 (lower left panel), the 68 % CL region is dominated by coannihilation, whereas its extension to the 95 % CL is dominated by the H / A funnel, with small areas of coannihilation. In this case there is a much more prominent coannihilation strip at the 95 % CL. Finally, in the pMSSM10 (lower right panel), we see two patches with small connected by a narrow ‘isthmus’, and a ‘continental’ region at large , where the 68 % CL is located. As indicated by the green and pink shadings, the dominant DM mechanisms in the ‘islands’ are and coannihilation, rather than coannihilation. We also note the reappearance of the h and Z funnel bands at low . In this model, the lower-mass ‘island’ and part of the higher-mass ‘island’ can be explored by future LHC searches for  [6], since . This search channel is less powerful in the CMSSM, NUHM1, and NUHM2, where in general, particularly in the coannihilation regions that are favored at the 68 % CL, as seen in Fig. 5.

Fig. 6

The planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). The red and blue solid lines are the and 5.99 contours. The diagonal black dashed lines correspond to and . In each of the CMSSM, NUHM1, and NUHM2 panels, the solid purple line is the current 95 % CL limit from the search in [107], and the dashed purple line is the 3000/fb projection from [108]. The solid and dashed purple lines for the pMSSM10, obtained using [109, 110], respectively, are the LHC run 1 95 % CL limit and the projected 3000/fb 95 % CL exclusion sensitivity with 3000/fb for a search, assuming a 100 % branching ratio and

Figure 6 also displays as purple lines the sensitivities of the most relevant present (solid) and prospective 3000/fb (dashed) searches, namely those for in the CMSSM, NUHM1 and NUHM2 cases, and for followed by + soft particles in the pMSSM10 case. We see that the current search does not impact the CMSSM, NUHM1 or NUHM2. In the case of the pMSSM10, the solid purple line is the current 95 % CL limit from the search in [109], assuming a 100 % branching ratio. This analysis is sensitive to stop topologies with b quarks in the final state where the decay products of the subsequent decay are undetected, and it was used in [6] to constrain compressed stop spectra. However, it is not sensitive to decays involving on-shell W bosons or . We conclude that future searches have the potential to explore parts of the coannihilation regions of the CMSSM, NUHM1, and NUHM2, and of the coannihilation region in the pMSSM10 case,11 but no DM channel can be fully explored by LHC searches.

The heavy Higgs bosons

We now study the differences between the dominant DM mechanisms in the pMSSM10 and the other models in the planes shown in Fig. 7. In the case of the CMSSM, the regions allowed at the 95 % CL and preferred at the 68 % CL (blue and red contours, respectively) are generally at considerably larger than the LHC bound (shown as a solid purple line) [111].12 The coannihilation mechanism dominates in a region around for , the H / A funnel dominates for , and there is an intermediate hybrid region. On the other hand, coannihilation dominates in an arc at larger . In the NUHM1, the hybrid and coannihilation regions are greatly expanded, and values of closer to the LHC bound are allowed at the 95 % CL. In the NUHM2, on the other hand, essentially all values of consistent with the LHC bound are allowed at the 95 % CL, the coannihilation mechanism dominates over most of the plane, leaving coannihilation to dominate for and . Finally, we see that in the pMSSM10 coannihilation dominates the 68 % CL region, that there is also a region at where coannihilation may be important, and that there are intermediate uncolored regions where neither of these mechanisms dominate. The LHC bound on is again saturated for .

Fig. 7

The planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right). The red and blue solid lines are the and 5.99 contours, and the solid purple line is the current LHC 95 % CL exclusion in the scenario

We have also estimated (not shown) the prospective LHC 95 % CL exclusion sensitivity in the H / A plane with 300/fb of data for the scenario, scaling the current limit (using results from [114-120]), and comparing with the estimated limits in [121], where good overall agreement was found. We estimate that could be explored for , reducing to for . This would cover much of the H / A funnel and hybrid regions in the CMSSM, portions of these regions in the NUHM1, and parts of the coannihilation regions in the NUHM2 and the pMSSM10, though not the regions of coannihilation in these models. Table 1 also summarizes the observability of the heavy Higgs bosons A / H in the different scenarios considered.

Direct dark matter searches

We now turn to the capabilities of direct DM search experiments to cast light on the various DM mechanisms. Figure 8 displays the planes for the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 (lower left) and the pMSSM10 (lower right), where is the cross section for spin-independent scattering on a proton. Our computation of follows the procedure described in [4], and we have once again adopted for the -nucleon term the value MeV. As previously, the and 5.99 contours are shown as red and blue lines. The sensitivities of the XENON100 [122] and LUX [123] experiments are shown as green and black lines, respectively, and the prospective sensitivity of the LUX-Zeplin (LZ) experiment [124] is shown as a purple line: the projected sensitivity of the XENON1T experiment [125] lies between the current LUX bound and the future LZ sensitivity. Also shown, as a dashed orange line, is the neutrino ‘floor’, below which astrophysical neutrino backgrounds would dominate any DM signal [126] (yellow region).

In the CMSSM case, we see that the current XENON100 and LUX data already put strong pressure on models where the focus-point or coannihilation mechanism dominates. There are borderline regions that are formally excluded by the data considered in isolation, but become permitted at the 95 % CL in a global fit including other observables, and also due to the uncertainties in the calculation of that have been included in the evaluation of the global function [4]. We also see that the coannihilation region and most of the H / A funnel region would be accessible to the planned LZ experiment. However, much of the coannihilation region lies below the LZ sensitivity, though it could be accessible to a 20-tonne DM experiment such as Darwin [127]. However, this region can be covered in the complementary direct searches at the LHC, as discussed previously. In the case of the NUHM1, the coannihilation and H / A funnel regions still lie largely within the LZ range, and the pure H / A funnel region is also within reach of LZ. However, the coannihilation and hybrid regions extend below the reach of LZ and Darwin, even below the neutrino ‘floor’. Again, this region could be partially covered by the complementary LHC direct searches. Similar qualitative conclusions apply to the NUHM2, with the additional observation that the small coannihilation regions lie below the LZ sensitivity and straddle the neutrino ‘floor’.

Finally, we see that whereas the region of the pMSSM10 parameter space that is favored at the 68 % CL lies within reach of the LZ experiment, as is the case for much of the coannihilation region, there are models in the and coannihilation regions, as well as in the h and Z funnels and uncolored regions where none of these mechanisms dominate, in which cancellations in the spin-independent matrix element may drive below the neutrino ‘floor’. It should be kept in mind here (see the discussion in [6]) that these very low values of are due to cancellations (for some discussions of these cancellations, see [128-132]) between different contributions to the matrix element for spin-independent scattering on protons. In [6] it was shown that similar cancellations hold when the cross section for spin-independent scattering on neutrons is considered, instead of the proton case shown in Fig. 8.

Table 1 also summarizes the observability of DM particles in direct searches in the different scenarios considered. We see a degree of complementarity between the LHC and direct DM searches.

We have focused in this article on the prospects for direct searches for DM scattering. A complementary probe of the properties of supersymmetric DM is through indirect detection, searching for the traces of DM annihilation in the Galaxy. A number of recent works have focused on this. For example, [133] has demonstrated that Fermi-LAT satellite limits on -ray emission in dwarf spheroidal galaxies [134] do not currently affect the parameter space of the pMSSM, although they may do so in the future.13 Constraints from IceCube and the HESS telescope have been investigated in [136, 137]. IceCube limits [138] do not currently affect the pMSSM10 parameter space, while HESS bounds [139] are primarily on pure wino states, which must have masses greater than a TeV [140, 141] in order to provide a thermal relic. Since the mass of the lightest neutralino is low in our models due to the incorporation of the constraint, this does not affect our fits. Accordingly we do not include these data sets in this work, but we plan on implementing likelihoods for these searches in future analyses.

Summary and conclusions

We have analyzed in this paper the mechanisms that play dominant roles in bringing the relic neutralino density into the range allowed by cosmology in the CMSSM, the NUHM1, the NUHM2 and the pMSSM10. We have delineated the regions of the parameter spaces of these models where dominant roles are played by , or coannihilation, or funnel regions where the neutralino annihilates rapidly via the heavy Higgs bosons H / A, and also regions where the neutralino has a significant Higgsino component. In the CMSSM, the NUHM1 and the NUHM2 we find that different mechanisms operate in different regions of the parameter spaces, with relatively small hybrid regions where two mechanisms contribute. In the pMSSM10, coannihilation dominates in most of the parameter space, with some contributions from other processes in specific ragions.

Our assessments of the observability of supersymmetry within different models, depending on the dominant DM mechanisms, are summarized in Table 1. Within the CMSSM, the NUHM1 and the NUHM2,  searches at the LHC can explore significant portions of the coannihilation regions. These regions also offer the possibility that the may be relatively long-lived, and detectable at the LHC as a long-lived charged particle. There are regions of the NUHM2 parameter space where coannihilation dominates, which can also be explored by  searches at the LHC. The coannihilation and focus-point regions of these models can be explored by the LUX-Zeplin direct DM search experiment. Much of the H / A funnel regions in these models may be explored via LHC searches for the heavy Higgs bosons, and also via direct DM searches with the LUX-Zeplin and Darwin experiments. On the other hand, the coannihilation regions seem likely to lie beyond the reaches of the direct DM searches.

Within our analysis of the pMSSM10, coannihilation is the dominant DM mechanism in most of the parameter space, though this might change with a different set of independent pMSSM parameters. In addition, there are regions where coannihilation or direct-channel annihilation via a h and Z funnel may dominate. Parts of the pMSSM10 model space can be explored at the LHC via  , H / A and other searches, and parts by direct DM searches. However, we find no long-lived particle signature in the region of the pMSSM10 parameter space that is currently favored statistically. Overall, large parts of the pMSSM10 parameter space could escape the LHC searches considered here, but a large fraction would be accessible to future DM experiments.

Our analysis shows that the LHC and direct matter searches offer significant prospects for discovering SUSY if it is responsible for the cosmological CDM, and in many cases the mode of discovery can reveal the nature of the dominant mechanism responsible for determining the CDM density. We look forward with interest to learning what the LHC and direct searches will be able to tell us about SUSY DM.