We give a brief review of the current constraints and prospects for detection of higgsino dark matter in low-scale supersymmetry. In the first part we argue, after performing a survey of all potential dark matter particles in the MSSM, that the (nearly) pure higgsino is the only candidate emerging virtually unscathed from the wealth of observational data of recent years. In doing so by virtue of its gauge quantum numbers and electroweak symmetry breaking only, it maintains at the same time a relatively high degree of model-independence. In the second part we properly review the prospects for detection of a higgsino-like neutralino in direct underground dark matter searches, collider searches, and indirect astrophysical signals. We provide estimates for the typical scale of the superpartners and fine tuning in the context of traditional scenarios where the breaking of supersymmetry is mediated at about the scale of Grand Unification and where strong expectations for a timely detection of higgsinos in underground detectors are closely related to the measured 125GeV mass of the Higgs boson at the LHC.
1. Introduction
From the particle physics point of view, the simplest, most popular, and arguably most robust mechanism leading to the correct amount of cold dark matter (DM) in the early Universe is thermal freeze-out (see, e.g., [1–4]). Briefly stated, one assumes that the DM consists of one or more matter species that were originally in thermal equilibrium with the Standard Model (SM) after the Big Bang and that, as the Universe expanded and cooled down, “froze” out of equilibrium when their number density became too low for annihilation and creation processes to take place.
As is well known, in the context of the freeze-out mechanism the measurement of the relic abundance provided by WMAP and Planck, ΩPLh2=0.1188±0.0010 [5, 6], implies a rather specific value for the thermally averaged annihilation cross section of the DM into SM particles: 〈σv〉≈3×10-26cm3/s≈1pb. Nevertheless, the thermal mechanism fails to provide any additional information on the nature of the DM itself since a cross section of that size can result from a discouraging wide range of DM mass values, spin quantum numbers, and DM-SM coupling strengths. Thus, in lack of more information, one has almost always to resort to some theoretical assumptions in order to narrow the search for DM down.
Since the 1990s, expectations about the scale of the new physics beyond the SM (BSM) have been driven by the theorists’ discomfort with the hierarchy problem. This is the well-known fact that, in a low-energy effective theory that includes one or more light fundamental scalars (as likely is the SM with a Higgs boson), one expects enormous quantum corrections to the scalar’s mass from the physics in the UV (the Planck scale, in the absence of anything else). Given the broad separation between the characteristic energies in play, this means that in order to get electroweak symmetry breaking (EWSB) one should fine tune the fundamental (unknown) Lagrangian parameters at the level, again in the absence of anything lighter than the Planck scale, of one part in ~1028, unless, of course, additional degrees of freedom were present, preferably close to the Higgs mass itself (say ~100-1000GeV).
Remarkably, simply on dimensional grounds, if one of these expected TeV-scale BSM particles were to be the DM, its coupling to the SM extracted from the freeze-out mechanism would be of the size of the electroweak coupling constant, g≈(16πmDM2·1pb)1/4≈0.1-1. This fascinating coincidence, which, in light of its singling out specifically weakly interacting massive particles, or WIMPs, is known as the “WIMP miracle,” maintains its attractiveness to these days, even if the LHC has failed to discover new particles below the scale of approximately 2TeV [7, 8].
Arguably the most complete and well motivated of the known BSM theories still remains low-scale supersymmetry (SUSY) (see, e.g., [9], for a popular review). From the theoretical point of view, not only does SUSY provide possibly the most elegant solution to the hierarchy problem (if one allows for the possibility that, given the current LHC bounds, the theory might have to be amended to regain full naturalness), but it also leads to a more precise UV unification of the gauge couplings than in the SM alone; it provides a solid rationale for the measured value of the Higgs boson and top quark masses and, by extension, for radiative EWSB. From the phenomenological point of view, the Minimal Supersymmetric Standard Model (MSSM) contains all the necessary ingredients for successful baryogenesis and provides a framework for cosmic inflation. It thus makes sense that, of all possible candidates for WIMPs, through the years a lot of attention was dedicated to the particles of the MSSM.
In this review we give a compact summary of the subject of DM in the traditional MSSM. After briefly surveying the particles with the potential of providing a good DM candidate, we argue that the nearly pure higgsino neutralino survives to these days as perhaps the only one that is not in substantial tension with any phenomenological constraint. Interestingly, it does so in a relatively model-independent way, without the need of resorting to narrow or secluded regions of the parameter space. We will thus review the higgsino’s prospects for detection in direct underground DM searches, indirect searches for DM in gamma-ray, and neutrino telescopes and at the LHC. Incidentally we will show that, in those models where SUSY breaking is transmitted to the visible sector at the scale of Grand Unification (GUT), the detection prospects of higgsino DM become tightly bound to the typical mass of the sfermions in the spectrum and, as a direct consequence, to the size of the Higgs boson mass.
In recent months several comprehensive reviews on the status of WIMP dark matter have appeared in the literature [10–13], one of which, coauthored by one of us, dedicated a full chapter to the MSSM neutralino with particular attention to the detection prospects of a ~1TeV higgsino. While that work is broader in scope, casting light on the experimental opportunities provided by neutralinos in the context of the wider picture of thermal DM models, DM constraints, and existing experimental anomalies, we concentrate here instead on the specific physical characteristics of higgsinos, underlining what we believe makes them currently stand out as the most interesting elements in the DM panorama of the MSSM. In this we are not dissimilar, perhaps, to recently appearing studies in the same tone [14, 15].
The structure of the review is as follows. In Section 2 we recall the particles of the MSSM that can provide a good DM candidate, classifying them according to their transformation properties under the SM gauge symmetry group. In Section 3 we single out the higgsino as the most promising candidate of the list and review its detection prospects in different and complementary experimental venues. We dedicate an additional subsection to the calculation of typical fine tuning and expectations for the scale of the superpartners in models constrained at the GUT scale. We summarize the main treated points and conclude in Section 4.
2. Dark Matter in the MSSM
One of the features making the MSSM very attractive from a phenomenological point of view is that its gauge symmetry structure originates directly from the supersymmetrization of the SM itself. As such, the fundamental gauge symmetry is SU(3) × SU(2) × U(1), and the dimensionless couplings are of the strong, electroweak, or SM Yukawa type.
One of the consequences is that a potentially viable DM particle is also expected to interact with SM-like strength. Since cosmological observations have long excluded the possibility of DM particles being charged under color [16] and, on the other hand, the DM is by definition “dark,” or practically electrically neutral [3, 17], one is led to conclude that all viable DM candidates in the MSSM must be classifiable on the basis only of the SU(2) representation they belong to. Moreover, the available representations are limited to those that can be found in the SM: SU(2) singlets, doublets, and the adjoint.
Before we proceed to briefly review these three groups individually, we remind the reader that in order to make the lightest SUSY particle (LSP) stable on cosmological time scales, one introduces in the MSSM an additional discrete symmetry, R-parity [18–22], under which only the superpartners of the SM fermions, gauge bosons, and any Higgs scalar field are odd. The origin of R-parity is still an active subject of research, and addressing the issue goes beyond the scope of the present review. We just point out that R-parity violation is strongly constrained phenomenologically, by the proton decay rate and electroweak precision measurements [23].
The only particles of the MSSM that are electrically and color-neutral are the neutrinos, their scalar superpartners, called sneutrinos, and, finally, the neutralinos. Neutralinos, χi=1,..,4, are Majorana fermion mass eigenstates emerging, after EWSB, from the diagonalization of the mass matrix of four electrically and color-neutral SUSY states (see [24–27] for early studies and [3] for a comprehensive, classic review). Two of these particles are gauginos, fermionic superpartners of the SM gauge bosons. The bino, B~, in particular, is the partner of the U(1) gauge boson, while the neutral wino, W~, is the partner of the SU(2) gauge boson W3. The other two states are neutral higgsinos, H~u and H~d, which belong to a vector-like pair of Higgs doublet superfields. If the lightest neutralino, hereafter, indicated simply with χ, is the LSP it can be the DM particle.
At the tree level, the neutralino mass matrix takes the following well-known form:(1)Mχ=M10-g′vd2g′vu20M2gvd2-gvu2-g′vd2gvd20-μg′vu2-gvu2-μ0,where g and g′ are SU(2) and U(1) gauge couplings, respectively, vu and vd are the vacuum expectation values (vev) of the neutral components of the scalar Higgs doublets, M1 and M2 are the soft SUSY-breaking bare masses of the bino and wino, respectively, and μ is the vector-like mass parameter of the Higgs doublet superfields.
In the remainder of this section we give an overview of the mentioned DM candidates of the MSSM, highlighting the strongest phenomenological constraints that can be applied in each case. We will not, however, discuss the neutrinos. It has been long known [28, 29] that the SM neutrinos do not provide, on their own, a viable candidate for cold DM. Their mass is O(<eV), so that they are relativistic at the time of decoupling and therefore incur strong constraints from structure formation [30–32]. On the other hand, heavy right-handed neutrinos, whose existence might be postulated on the ground of the observed neutrino masses and could provide a naturally expected extension of the traditional MSSM, also do not provide a good candidate for DM because they are not protected by R-parity and therefore not stable over cosmological scales in most scenarios.
2.1. SU(2) Singlets
(Nearly) Pure Bino. The first SU(2) singlet DM candidate we present is the bino. Because of EWSB, a pure bino state does not exist in the MSSM, but the lightest neutralino behaves like a pure bino to a very good approximation, after the diagonalization of Mχ, if |M1|≪M2,μ.
The interactions of the bino-like neutralino with the SM fields are easily found by directly supersymmetrizing the SM gauge-fermion-fermion interaction and applying the R-parity conservation constraint. The resulting vertex takes the form bino-sfermion-fermion, L⊃-XLf~Lχ¯PLf-XRf~Rχ¯PRf, where tree-level couplings, XL,R=2g′YL,R, are expressed in terms of the hypercharge assignment YL,R of the fermion Weyl spinors.
The pair-annihilation of bino-like neutralinos in the early Universe proceeds at the leading order through the t-channel diagram shown in Figure 1(a). The region of the MSSM parameter space where Ωh2≈0.12 is obtained in this way is historically known as the bulk [33, 34]. One can calculate the thermal cross section for binos, given approximately by the following [35]:(2)σvB~≈∑f~g′4Yf~42πmχ2mf~4+mχ4mf~2+mχ24TFmχ,in terms of the neutralino (bino) mass, mχ, sfermions’ mass mf~, hypercharge Yf~, and freeze-out temperature TF, which parameterizes the dependence on velocity of the p-wave cross section, and is set here approximately at TF≈(0.04-0.05)mχ.
(a) The dominant early-Universe annihilation channel for a nearly pure bino-like neutralino. (b), (c) Examples of annihilation and coannihilation tree-level channels into gauge bosons for a predominantly higgsino-like neutralino.
The bulk has been long known to be strongly constrained by direct SUSY searches at colliders. To give a semiquantitative estimate of these constraints, let us assume that only selectrons and smuons belong to the light SUSY spectrum, a reasonable ansatz in light of the strong LHC bounds on particles with color [36–38]. Assuming all four left- and right-handed slepton states have the same mass, and inserting Yf~L=-1/2, Yf~R=-1 in (2) one finds that the cross section is typically much smaller than ~1pb, except in the range mχ<mf~≲100GeV. A charged slepton mass of this size has been long excluded by direct searches at LEP [39].
If, instead of selectrons and smuons, the light sfermions happen to be status, the parameter space opens up a little, mτ~1≲150GeV for mχ≈50GeV, due to the nonnegligible mixing between left and right chiral slepton states, which introduces an s-wave component to the annihilation cross section (see, e.g., [40]). Nevertheless, LHC bounds on electroweak production [41], implying mτ~1≳109GeV, are by now becoming strongly constraining for these scenarios too, which will be probed even more deeply soon [42]. Finally, as we have mentioned, SUSY parameter space where bulk sfermions are charged under color is strongly excluded by LHC direct searches.
A way to evade the strong collider bounds is provided, if the bino-like neutralino and some other sparticles (sfermions f~ or other gauginos) are nearly degenerate in mass, by the mechanism of coannihilation [43–45]. In this case the cross section of (2) should be replaced with an effective quantity that takes into account the thermal average of all annihilations and coannihilations of the kind χχ,χf~,f~f~→SMSM, some of which are likely to be much more efficient than χχ→SMSM alone.
However, without any guidance from the theory in the UV, coannihilation of the bino with other sparticles can only be achieved in narrow slices of the parameter space, which require some tuning of the initial parameters to engineer the desired coincidence of neutralino and sfermion mass. And in models that are instead defined in terms of a limited number of free parameters in the UV, like the CMSSM [46], in which slepton or stop coannihilation with the bino can occur naturally for particular choices of the initial conditions, the preferred regions of the parameter space are incurring increasingly strong limits from direct LHC searches [13, 47–50]. Besides, with gaugino universality at the GUT scale, it is a struggle to fully accommodate the measured value of the Higgs mass at the LHC [47, 51] (this problem is resolved if the gluino mass is a free parameter, e.g., [52]). Thus, even if coannihilation of the bino with other sparticles can still lead to viable regions of the parameter space in the most generic parametrizations of the MSSM [53], it is also perhaps not exceedingly attractive from a natural point of view.
R Sneutrino. The second SU(2) singlet DM candidate of the MSSM is the scalar “right-handed” sneutrino. The right-handed sneutrino does not properly belong to the MSSM, which in its original formulation features massless neutrinos, but naturally emerges in SM extensions with right-handed neutrinos, which can give rise to the neutrino mass via small Yukawa couplings (if the right-handed neutrino is Dirac), or through the see-saw mechanism (if the right-handed neutrino is Majorana, see, e.g., [54] and references therein).
The phenomenology of right-handed sneutrinos as DM, however interesting, is very model-dependent. In traditional see-saw models with large-scale Majorana mass the right-handed sneutrino is too heavy to be the DM. On the other hand, for a sneutrino of the “Dirac” type, or, in alternative, Majorana but such that the bare mass is of the order of the superpartners’ mass [55, 56], the only really model-independent vertex with the SM involves a very small Yukawa coupling L⊃-yνRe¯LH~u±ν~R-yνRν¯LH~u0ν~R. Thus, the induced t-channel processes similar to Figure 1(a), with sneutrinos (charginos) in place of neutralinos (sfermions), and a tiny coupling constant are not strong enough to get the correct Ωh2.
On the other hand, the correct relic density can certainly be obtained thanks to the mixing with the left-handed sneutrino, and SUSY breaking can generate A-terms of the order of the SUSY scale, which provide large couplings to the SM Higgs boson. The phenomenology of these cases can be very rich and exceeds the scope of this review. We direct the reader to the vast literature on sneutrino DM for further details (see, e.g., [57–60], for early studies and bounds, and [61] for a recent LHC analysis).
2.2. SU(2) Doublets
We have seen that singlet DM candidates in the MSSM are accompanied by some uncomfortable features: they are either strongly constrained by collider bounds, are only viable in fine-tuned regions of the parameter space, or present a phenomenology that is highly model-dependent. We therefore move on to reviewing the next set of candidates, the SU(2) doublets.
(Nearly) Pure Higgsino. The most popular SU(2) doublet DM candidate, and the one that appears to us most attractive from a phenomenological point of view, is the higgsino, which is the main subject of this review. As was the case for the bino, there is no pure higgsino state after EWSB, but one obtains an almost pure higgsino-like neutralino by diagonalizing Mχ in (1) in the limit |μ|≪M1,M2.
As supersymmetry assigns a Weyl spinor to each complex state in the scalar Higgs doublets one counts four physical higgsino states, which, after EWSB, give rise to two Majorana neutralinos, χ1 (or χ) and χ2, and a Dirac chargino, χ±. When |μ|≪M1≈M2, the tree-level mass splitting between the two higgsino-like neutralinos is of approximately the size of mZ2/M1,2 [9], and the splitting between the higgsino-like chargino and the lightest neutralino is approximately half of that. Moreover, radiative corrections also induce a nonnegligible and irreducible mass splitting (~100s MeV) between the charged and neutral states (see, e.g., [62, 63]).
To correctly compute the thermally averaged effective cross section that yields the DM relic abundance, one must take into account all possible annihilations and coannihilations of higgsino states. For mχ above the W threshold the dominant final state is into W and Z bosons (Figures 1(b) and 1(c) give examples of possible diagrams for this processes), to which higgsino-like neutralinos and charginos couple through the electroweak charged and neutral currents [3]:(3)L⊃-g2Wμ+χ¯γμχ--g4cosθWZμχ¯1γμχ2+h.c.-g2cosθWZμχ¯+γμ1-2sin2θWχ-.
The effective cross section can be obtained at the leading order in the limit of all four states being degenerate (see, e.g., [35]):(4)σvH~eff≈21g4+3g2g′2+11g′2512πmχ2.
For heavy, very pure higgsinos, one should include in the calculation of 〈σv〉H~(eff) corrections due to the Sommerfeld enhancement, a well-known nonperturbative effect originating from the fact that if a DM particle is much heavier than the electroweak gauge bosons and relatively slow, the weak force becomes effectively long-range and the impact of the nonrelativistic potential on the interaction cross section becomes significant [64, 65]. However, in the case of the higgsino the splitting between its charged and neutral components is almost always large enough to effectively wash out substantial nonperturbative effects originating from the resummation of ladder diagrams [66–68], so that in a first approximation (4) provides a fairly accurate estimate of 〈σv〉H~(eff).
One can see that the cross section is typically much larger than ~1pb, unless mχ≈1TeV (the precise numerical value is more about 1.1TeV, as we shall see). Thus, ~1TeV higgsino is on its own a good candidate for the DM in the Universe [69], while a higgsino much lighter than 1TeV requires one to assume the existence of an additional DM component (e.g., axion [70, 71]), needed to get Ωh2≈0.12.
As we shall see in the next sections, ~1TeV higgsino is generally associated with a large SUSY-breaking scale and for this reason it is not currently very constrained from a phenomenological point of view. However, its characteristic properties can give us hope for a timely detection in direct and indirect DM searches and even, if mχ≪1TeV, in collider searches.
L Sneutrino. We conclude this subsection by reviewing the properties of the only other SU(2) doublet DM candidate in the MSSM: the “left-handed” sneutrino, scalar superpartner of the SM left-handed neutrino.
The left-handed sneutrino is a complex scalar field with SU(2) × U(1) quantum numbers equal to the higgsino’s. Like the higgsino, it has charged and neutral current couplings to the W and Z bosons, L~-ig/2Wμ+ν~L∗∂μe~L-+Wμ-e~L+∂μν~L-ig/(2cosθW)Zμν~L∗∂μν~L. The mass splitting of the charged and neutral components of the SU(2) doublet is, however, much larger for sneutrinos/sleptons than for higgsinos, being generated through hypercharge D-term contributions [9]: me~L2-mν~L2≈-mW2cos2β, where tanβ≡vu/vd. Thus, one should resist the temptation of interpreting (4) as an accurate estimate of the effective cross section for sneutrinos too, since the coannihilation of charged and neutral states becomes somewhat less efficient. It turns out [60] that the mass required to produce 〈σv〉ν~L(eff)≈1pb is about mν~L≈600-700GeV. Sneutrinos lighter than that imply the existence of an additional component of DM.
A very important constraint on left-handed sneutrinos as DM arises because they, unlike the Majorana higgsino-like neutralinos, are not their own antiparticle, so that their elastic scattering with nuclei in direct detection experiments proceeds also through t-channel exchange of a Z boson. By virtue of the sneutrino’s neutral current coupling, the spin-independent cross section is approximately given by a Fermi-like contact interaction, σpSI≈μred2GF2/8π≈10-3pb=10-39cm2, where reduced mass μred≈mp for mν~L≫mp. Cross sections of this size have been long excluded in underground detector searches [72, 73].
2.3. SU(2) Adjoint Triplet
(Nearly) Pure Wino. The only SU(2) triplet DM candidate in the MSSM is the wino-like neutralino, dominated by the fermionic superpartner of the W3 weak gauge boson. The wino belongs to the adjoint representation of the gauge group (hypercharge Y=0) and the wino-like neutralino emerges, after EWSB, from the diagonalization of (1) in the limit |M2|≪M1,μ. One finds a Majorana neutralino, χ, and a Dirac chargino, χ±, mass-degenerate at the tree level. In the context of UV complete models of SUSY-breaking, spectra with a light wino can arise, for example, in scenarios where SUSY breaking is transmitted via anomaly mediation [93, 94].
If the wino LSP is heavier than the electroweak gauge bosons, its dominant final state channel for annihilation (and coannihilation with charginos) in the early Universe is into W (but not Z) boson final states, to which it couples as L~-gWμ±χ¯γμχ∓. The thermal annihilation cross section is dominated by coannihilations of the three wino states, similarly to what happens for the doublet higgsinos. Annihilation into fermion–antifermion final states through a t-channel sfermion exchange, reminiscent of the bino bulk mechanism, has been instead long excluded by LEP limits on the charged slepton masses.
Unlike higgsinos, in the wino case mass splitting between the charged and neutral fermion component of the SU(2) multiplet is generated exclusively by radiative corrections, ΔMW~=(g2/4π)mWsin2(θW/2)≈166MeV [95]. Note that the mass splitting is typically much smaller than for higgsinos, so that one cannot neglect the effects of the Sommerfeld resummation on the calculation of the thermal cross section. When one includes the Sommerfeld enhancement numerically, the correct relic density is obtained for mχ≈2.7-2.8TeV [66–68]. For a lighter mass, winos do not saturate the relic abundance.
The Sommerfeld enhancement induces more dramatic modifications of the effective DM annihilation cross section when the average kinetic energy of the WIMP corresponds to speeds of the order of 10-3c, as in the present-day Universe. This fact has led to the derivation of powerful indirect astrophysical constraints on the annihilation cross section of wino-like neutralinos [91, 96–99]. By taking into account the effects of Sommerfeld-enhanced contributions to the annihilation of winos into monochromatic gamma rays, as well as bounds on the present-day cross section to W+W- from diffuse gamma radiation from the Galactic Center and Dwarf Spheroidal satellite galaxies (dSphs), measured in terrestrial and space telescopes HESS [89, 100] and Fermi-LAT/MAGIC [88], and from cosmic-ray (CR) antiproton data at AMS-02 [90, 91], one can derive strong independent constraints (albeit affected by significant systematic uncertainties) which steeply raise the stakes on the wino as a viable DM particle, especially in scenarios where it saturates the relic abundance.
2.4. Mixed Cases
The four neutralinos of the MSSM are all Majorana fermions that, after EWSB, remain neutral under U(1)em and color. In the absence of a well-separated hierarchy among M1, M2, and μ, the lightest mass eigenstate will be an admixture of the SU(2) gauge multiplets discussed in Sections 2.1–2.3 but, unlike those cases, it will present properties that differ significantly from a pure gauge eigenstate.
When |M1|≈|μ| the neutralino is in a highly mixed bino/higgsino state. Mixed neutralinos of this kind (sometimes also called “well-tempered” [35]), originally observed in mSUGRA parameter space [101–103] but that can arise under different boundary conditions (e.g., [104, 105]), enjoyed some popularity, especially before the advent of the LHC, because they can easily lead to Ωh2≈0.12 for values of the μ parameter as low as few hundreds GeV, which are favored to solve the hierarchy problem. However, the rapid progress made in the bounds on the spin-independent cross section of the neutralino scattering off nuclei in direct WIMP detection searches, combined with a failure to directly observe scalar fermions and heavy Higgs bosons at the LHC, have rendered scenarios where the lightest neutralino is a rich admixture of gaugino and higgsino much less appealing if not excluded altogether (see, e.g. [106], for a very recent update of the constraints on bino-higgsino, and [99] for wino-higgsino scenarios).
To briefly set the issue on quantitative grounds, let us estimate the strength of the coupling with which neutralino admixtures of higgsino and gaugino contribute to the spin-independent cross section. We recall that, in the limit of the squarks and heavy Higgs bosons being much heavier than mh=125GeV, which has become a reasonable assumption after the first two runs of the LHC, the main interaction between the neutralino and heavy nuclei in underground detectors proceeds as in Figure 2, via t-channel exchange of the 125GeV Higgs boson and an effective coupling to gluons through the heavy quark loops.
The main interaction between the neutralino and heavy nuclei in underground detectors in the limit of squarks and heavy Higgs bosons being much heavier than mh=125GeV and in general outside of LHC reach.
As the neutralino LSP-Higgs-neutralino LSP tree-level vertex directly stems from applying the gauge covariant derivative on the Higgs doublets, it is nonzero only for a gaugino/higgsino admixture. For tanβ sufficiently large to ensure a predominantly SM-like Higgs boson (tanβ>3-4 is a condition often fulfilled, for instance, in scenarios where EWSB is obtained radiatively via the renormalization group evolution of soft SUSY-breaking parameters constrained at some high scale, as it prevents certain soft masses from running tachyonic at the low scale.), the coupling to the nucleon can thus be expressed entirely in terms of the higgsino fraction (or purity), fh, which depends on the elements of the unitary matrix, N, diagonalizing (1).
If diag[mχ1,mχ2,mχ3,mχ4]=NMχN†, one can define fh≡|N13|2+|N14|2 and express the coupling of interest as L~(gfh1-fh/4)χ¯χh. Note, incidentally, that deriving an explicit form for the elements of matrix N in terms of bare masses M1, M2, and μ is not a trivial task even at the tree level, and useful formulas in this regard can be found in several papers, for example, [107–110]. By simple inspection of (1), however, one can infer a rough approximation for the higgsino fraction in the limit of nearly pure higgsinos, |μ|≪M2≈M1:(5)1-fh≈mW2M1,2-μ2.Equation (5) becomes quite accurate for fh≳0.999.
The spin-independent cross section of the neutralino with protons (nucleons), σpSI=4μred2/πAp2, can be parameterized for moderate-to-large tanβ simply as follows [3]:(6)Apfh≈aefffTG9mpvgfh1-fhmh2,in terms of the gluon fractional content of the proton, fTG (we use the default value for micrOMEGAs v4.3.1 [111], fTG=0.92), and a phenomenological fudge factor, aeff≈0.9-1, which takes into account the dependence of Ap on twist-two operators [112] and higher-order loop corrections [113].
We show in Figure 3 a plot of σpSI as a function of purity fh for a mχ=1TeV neutralino (to a first approximation the DM mass affects the cross section only through the reduced mass leading to μred≈mp). One can see that, for admixtures dominated by the higgsino fraction, the most recent XENON-1T 90% CL upper bound [75] on σpSI enforces fh>98%, so that viable DM candidates ought to be very close to a pure higgsino state.
The neutralino-proton spin-independent cross section, σpSI, for a typical case of predominantly higgsino-like neutralino DM with mχ=1.0TeV as a function of higgsino purity fhiggsino (≡fh).
Since the purity of well-tempered higgsino-dominated neutralinos stays well below 90% in those models attempting to provide a satisfactory solution to the hierarchy problem while saturating the relic abundance [35], we conclude that, barring increasingly narrow corners of the parameter space [106], these scenarios have become very hard to rescue or justify in light of the most recent direct detection bounds.
To conclude this subsection, we finally recall that, in cases where |M1|<|μ|≲1-2TeV, one obtains scenarios where the mixed neutralino is predominantly bino-like, but also acquires couplings that originate from its admixture with higgsino states, so that additional mechanisms for obtaining 〈σv〉≈1pb with respect to Section 2.1 are possible.
These mechanisms, often called funnels, involve resonant or close-to-resonant s-channel annihilation of two neutralino LSPs via a nearly on-shell mediator which could be the Z boson (if mχ≈mZ/2) [25], the SM Higgs boson (if mχ=60-65GeV) [114], or one of the heavy Higgs bosons of the MSSM [33].
Note that the Z-funnel parameter space is strongly constrained by the LHC. The coupling of the lightest neutralino to the Z boson is due exclusively to the isospin neutral current, cf. Section 2.2, which means that in mixed bino-higgsino scenarios it is directly proportional to the higgsino fraction. As a consequence, fh cannot take excessively small values or, in other words, μ cannot be much larger than M1≈mZ/2. The relative proximity of a mostly higgsino-like chargino and a mostly bino-like neutralino subjects this region of the parameter space to strong bounds from direct LHC multi-lepton searches [115].
Light and heavy Higgs boson funnels are less constrained from direct LHC SUSY searches than the Z funnel, since the direct coupling to the lightest neutralino is dependent on fh and the mediator can be quite heavy. However, there exist complementary observables which can constrain these regions, like the branching ratio BR Bs→μ+μ- [116] and direct searches for heavy Higgs bosons in the ττ channel [117]. Moreover, as was the case for the coannihilations of the bino, most phenomenological scenarios require ad hoc arrangement of the parameters to obtain the right ratio of neutralino to scalar mass, although this is not necessarily the case for some parameter-space regions of GUT-constrained scenarios like the CMSSM, in which the renormalization group evolution (RGE) of soft masses from a handful of free parameters can lead more naturally to the right mass coincidence (see, e.g., [118, 119] for early studies).
3. Phenomenology of Higgsino Dark Matter
The discussion of Section 2 has led us to conclude that the sole DM candidate of the MSSM emerging almost unscathed from the wealth of observational data of recent years is the nearly pure higgsino. We therefore dedicate this section to the analysis of the prospects for detection of a higgsino-like neutralino in direct DM detection searches, collider searches, and indirect astrophysical signals, and spend a few words on alternative strategies in other experimental venues. We will also give some predictions for the scale of the superpartner particles in traditional models and briefly discuss the issue of fine tuning.
3.1. Prospects for Detection in Direct and Indirect Searches
We begin in Figure 4(a), where we plot the rescaled spin-independent neutralino-nucleon cross section versus neutralino mass for a nearly pure higgsino under CMSSM/mSUGRA boundary conditions [46] (We remind the reader that this means scanning simultaneously over 4 free parameters: m0, the universal soft SUSY-breaking scalar mass at the GUT scale; m1/2, the universal GUT-scale gaugino mass; A0, the universal GUT-scale soft trilinear coupling; and tanβ, the ratio of the Higgs doublets’ vevs. We scan them in this study over broad ranges: m0,m1/2∈[0.1TeV,30TeV], A0∈[-30TeV,30TeV], tanβ∈[1,62]. Additionally, one chooses the sign of μ, which we set here to positive, as its sign does not much affect the region of parameter space with nearly pure higgsino DM (see, e.g., [47, 116]). Note that the chosen input mass ranges encompass the parameter space region shown in Figure 4 in its entirety. In it one finds m1/2≲0.6m0, with 5TeV≲m0≲25TeV, 2.5TeV≲m1/2≲15TeV due to the Higgs mass measurement, see discussion below.). The color code depicts the higgsino DM relic abundance. For the points of the parameter space corresponding to Ωh2 below the Planck measurement [6], ΩPLh2≈0.12, we directly rescale σpSI by ξ=Ωh2/ΩPLh2, assuming implicitly that the fraction of higgsino DM we measure locally today traces closely its early time large-scale freeze-out value. Solid tilted lines show recent direct upper bounds from the PandaX-II [74] (maroon) and XENON1T [75] (blue) underground experiments. The latter is not much more constraining than an earlier bound from the now decommissioned LUX [76]. Dot-dashed lines show the projected reach of several upcoming and planned experiments.
(a) Spin-independent neutralino-nucleon cross section σpSI rescaled by the relic abundance, as a function of neutralino mass mχ, for a nearly pure higgsino with CMSSM/mSUGRA boundary conditions subject to mh≈125GeV and LHC Higgs bounds. Solid lines show the 90% CL upper bounds from PandaX-II [74] (maroon) and XENON1T [75] (LUX [76]) (blue). Dot-dashed lines show the projected reach for DEAP-3600 [77] (orange), XENON1T/nT [78] (blue), DarkSide G2 [79] (maroon), LZ [80] (black), DARWIN (purple) [81]. Thin solid black line shows the current lower bound on mass from direct searches at the LHC [82, 83]. (b) Rescaled spin-dependent neutralino nucleon cross section σpSD as a function of neutralino mass mχ, for the a nearly pure higgsino in the CMSSM/mSUGRA. Solid lines show the 90% CL indirect upper bounds from IceCube [84] (green) and Antares [85] (red). Dashed lines show projections for LZ [86] (violet), XENON1T [78] (purple), Pico-500 [87] (blue), and DARWIN [81] (black).
We also show in Figure 4(a) as a thin black line the current lower bound on mass from direct searches for compressed electroweakinos in final states with two low-momentum leptons at the LHC ([82, 83], following a proposal and case studies by [120, 121]), which is sensitive to higgsino DM for mass splitting mχ2-mχ1=3-30GeV. One should also be aware of the estimated putative reach of the ILC in testing higgsinos [122], which we do not show in the plot for lack of space. It extends to approximately 240GeV(480GeV), independently of mass splitting, if the beam energy is set to s=(500GeV)2 (s=10002GeV2).
In Figure 4(b) we show the rescaled spin-dependent neutralino-proton elastic scattering cross section, ξσpSD, versus neutralino mass. We show with solid lines existing indirect upper bounds from observations of neutrinos from the Sun in the neutrino telescopes IceCube [84] (green) and Antares [85] (red), interpreted for a predominantly W+W- annihilation final state, which give a good approximation for the nearly pure higgsino case [53, 123]. Dashed lines of different colors give various projections for the future direct reach in σpSD of underground detectors.
The relic density and DM observables are here calculated with micrOMEGAs v4.3.1 [111]. The supersymmetric spectrum is calculated with SPheno v4.0.3 [124, 125], and all model points are subject to LHC Higgs constraints from HiggsSignals/HiggsBounds [126–129] and to the Higgs mass measurement [130]. The Higgs mass is calculated, like the SUSY spectrum, with the latest version of SPheno, which yields, in the regime where soft SUSY-breaking masses are well above ~1TeV, a value in excellent agreement with other numerical packages, SusyHD [131] and FlexibleSUSY [132]. The calculated value is subject to an overall estimated theory uncertainty of approximately 2GeV [133], which we take into account in Figure 4. Note that when the SUSY spectrum lies in the several TeV regime or above, all electroweak precision and flavor observables, including the anomalous magnetic moment of the muon, are expected to roughly maintain their SM value.
We have chosen to show in Figure 4 the higgsino parameter space under CMSSM boundary conditions, which provide a reasonable ansatz for models with scalar universality inspired by supergravity, and more generally cast in a lean framework scenarios in which supersymmetry breaking is transmitted to the visible sector at some high scale (the GUT scale) and EWSB is obtained radiatively around the minima of the MSSM scalar potential. In models defined in this way one observes, for a higgsino-like neutralino, strong correlation between the Higgs boson mass and the allowed minimum value of σpSI. We show this in Figure 5, where we plot the lower bound on σpSI as a function of Higgs mass for a higgsino LSP of arbitrary mass. The correlation between minimum cross section and Higgs mass translates in Figure 4(a) into a lower bound on σpSI when mχ≈1TeV.
Lower bound on σpSI as a function of Higgs mass for a higgsino LSP of arbitrary mass in generic models where the breaking of supersymmetry is transmitted at the GUT scale and the physical spectrum and EWSB are obtained after RGE to the low scale.
To qualitatively understand what is happening, let us recall from Section 2.4 that in order to push down σpSI for a predominantly higgsino-like neutralino one must increase purity fh or, in other words, raise the wino and bino masses, cf. (5). Very heavy winos/binos at the GUT scale feed through the RGE on the low-scale value of the soft SUSY-breaking up-type Higgs doublet mass, which carry SU(2) isospin and hypercharge and also tend to push down the right-handed stop mass. This happens even in scenarios where the gluino mass is not universal and can be found relatively close to the higgsino, like those analyzed in [134].
In order to keep the Higgs doublet soft mass under control, so as to obtain a higgsino-like LSP after EWSB, and avoid tachyonic physical states, numerical scans are in this situation driven to large negative A0 and/or larger soft scalar mass. Both solutions have the net effect of pushing up the Higgs boson mass and give rise to the behavior we observe in Figure 5(The attractiveness, from the phenomenological point of view, of a lower bound on the neutralino scattering cross section determined by the Higgs mass measurement was pointed out early on in Bayesian analyses of the CMSSM/NUHM [13, 47, 116]. The exact minimal cross section depends strongly on the calculation of the Higgs mass itself, and on how it translates into mass predictions for the sparticles. In SPheno v4.0.3, mh≈125GeV leads to less optimistic expectations for the mean SUSY scale than in the versions of SOFTSUSY [135] or FeynHiggs [136] used in [47, 116]. Hence the parameter space in Figure 4(a) extends to lower σpSI values than in those studies.).
There is no apparent lower bound on the scattering cross section if we relax the requirement of radiative EWSB from boundary conditions generated at the GUT scale. This is the case, for example, in models where the typical mass of scalar particles is by several orders of magnitude decoupled from the electroweak vev (see, e.g., [137–139]), and one does not expect to infer strict relations between the mechanism of SUSY-breaking and EWSB. The relic density alone determines then the mass of the higgsino-like DM, and purity fh can be extremely close to 1. We generically indicate with a black arrow in Figure 4(a) the parameter space for higgsino DM in those models, which can extend well below the neutrino background floor [63, 140].
This highly inaccessible part of the higgsino parameter space proves particularly tricky to probe. For underabundant higgsinos, μ≪1TeV, interesting venues for detections can be provided, for very small mass splitting, mχ±-mχ≈150MeV, by future collider searches for disappearing tracks [141, 142]. If there is a sizable CP violating phase, future electron dipole moment experiments might be sensitive to parameter space with purity in excess of 99.99% [63]. And possibly new venues for detection are given by the cooling curve of white dwarfs [15]. Additional opportunities for the future detection of higgsino-like compressed spectra, in particular for long-lived particles with a relatively short lifetime, can arise then in electron-proton colliders [143].
We finally show in Figure 6 the status of indirect detection bounds and projections in gamma-ray searches in space and terrestrial telescopes for ~1TeV higgsino DM under CMSSM/mSUGRA boundary conditions (we implicitly assume that the chances for detection maximize if higgsinos saturate the relic abundance). In Figure 6(a), solid black line shows the most recent 90% CL upper bound on the present-day σv from the statistical combination of Fermi-LAT and MAGIC observations of dSphs [88], and the magenta line draws the recent bound from 10-year observation of the Galactic Center at HESS [89] under the Einasto profile assumption. We adopt the bounds in the W+W- final state interpretation, which give a good approximation for the ~1TeV higgsino.
(a) Indirect detection bounds and projections in gamma-ray searches in space and terrestrial telescopes for ~1TeV higgsino DM under CMSSM/mSUGRA boundary conditions. Solid black line shows 90% CL upper bounds on the present-day annihilation cross section to W+W- from the statistical combination of Fermi-LAT and MAGIC observations of dSphs [88]; solid magenta line shows the recent bound from 10-year observation of the Galactic Center at HESS [89] under the Einasto profile assumption; solid green line shows the upper bound from antiproton cosmic-ray (CR) data at AMS-02 [90] according to [91] for the NFW profile; and dashed blue line shows the projected reach of CTA 500h under the Einasto profile assumption [53]. (b) In magenta, the current 95% CL upper bound on the annihilation cross section (times velocity) to gamma-ray lines, σγγv, from HESS [92] under the Einasto profile assumption, compared to the cross section of our ~1TeV higgsino points.
For the W+W- final state we show in solid green the determination by [91] of the 95% CL upper bound on σv from antiproton CR data at AMS-02 [90], under the NFW profile assumption. Note that the bound is subject to uncertainties related to the choice of diffusion model for CR propagation in the Galaxy. Some of these choices can in fact weaken it [91], and push it up to approximately the level of the HESS limit. Finally, dashed blue line shows the projected statistical reach of CTA 500h, under the Einasto profile assumption [53, 144]. Note that including the systematic uncertainty from diffuse astrophysical radiation will most likely weaken the extent of the projected reach [123, 145]. Also note in Figure 6(a) that some model points are characterized by σv significantly above the thermal relic expectation, due to the presence of the heavy pseudoscalar Higgs mass at mA≈2mχ [47, 53]. Regions of the parameter space that allow for this serendipitous coincidence thus see their indirect detection prospects improve significantly.
We show in Figure 6(b), as a magenta solid line, the current 95% CL upper bound on the annihilation cross section (times velocity) to gamma-ray lines from the final 254h data at HESS [92] under the Einasto profile assumption. The line is compared to the cross section of our ~1TeV higgsino points, which lie well below the limit.
3.2. The Soft SUSY Scale and Fine Tuning
We conclude with a few words about the expected scale of the supersymmetric particles associated with higgsino DM. In truth, little is known in this regard, as the issue is highly model-dependent and there is not one only way of inferring the scale of SUSY breaking.
Of course, expressions similar to (5)-(6) can give us a lower bound on the scale of the electroweak gauginos for every given upcoming new constraint on σpSI, but to be precise one should then take into account the rich parametric dependence of the full formulas. Equivalently, the Higgs mass measurement tells us that in all likelihood stops and gluinos sit well above the LHC reach, but little more than that is known, as expectations depend strongly on parameters like tanβ and the trilinear coupling At.
Thus, without pretence of presenting any universally valid result, but to just show an example of a model where the measurement of the Higgs mass actually does provide predictions for the maximally allowed typical scale of the superpartners, we present in Figure 7(a) the distribution of the mean stop mass, MSUSY=(mt~1mt~2)1/2, under CMSSM/mSUGRA boundary conditions in the (mχ,ξσpSI) plane with higgsino DM. One can see that by approximately the next round of XENON-1T data we will be starting to probe the 10TeV range of the superpartners if the DM is entirely composed of higgsinos. Note also that, for higgsino mass mχ≲140GeV, the LHC is already excluding, with direct soft-lepton bounds on electroweakinos, the parameter space corresponding to MSUSY≲8-10TeV.
(a) A plot of MSUSY=(mt~1mt~2)1/2 in the (mχ,ξσpSI) plane with higgsino DM under CMSSM/mSUGRA boundary conditions. (b) EWSB fine tuning for points with higgsino DM in the (mχ,ξσpSI) plane.
Finally, like all BSM models developed at least in part to deal with the hierarchy problem, after the first two runs of the LHC models with higgsino DM have become marred by a certain amount of EWSB fine tuning. The severity of this issue depends, of course, on the specific features of each model: how EWSB is obtained and the relation to the mass of the Higgs boson. In the context of the CMSSM, the fine tuning associated with higgsino DM is shown in Figure 7(b), where we plot in the (mχ,ξσpSI) plane the size of the usual Barbieri-Giudice measure [146, 147] (following the prescription of [148]) (We remind the reader that the Barbieri-Giudice measure is generally defined as maxpi|∂logMZ2/∂logpi|, where pi are the model’s input parameters at the typical scale of the messengers for SUSY breaking. In the CMSSM these are the GUT-defined parameters m0, m1/2, A0, B0, μ0.). No point shows EWSB fine tuning of less than a part in 100, as direct consequence of the Higgs mass measurement, and one can observe the well-known fact that higgsino points favored by expectations of naturalness correspond to mχ<1TeV and lead to Ωh2≪0.12. For the specific case of the ~1TeV higgsino, a failure to observe a signal in, say, the next round of XENON-1T data will imply a fine tuning greater than one part in 103, with rapid increase with each successive milestone exclusion (There exist ways of embedding the MSSM in UV completions that can lead to lower fine tuning for higgsino DM, see, e.g., [134, 149].).
However, we emphasize that a large fine tuning is by no means exclusive to the CMSSM, to higgsino DM, or even to SUSY in general (see, e.g., [150] for fine tuning in a non-SUSY scenario). As a matter of fact, the majority of phenomenological DM models found in the literature do not even attempt to construct a UV completion that could directly relate their free parameters to the physics of the high scale. It is very possible that once a discovery is finally made many of the suspended questions will start to find their answers. Higgsinos appear to be just in the perfect position to usher, in case of their eventual discovery, a new era of understanding.
4. Summary and Conclusions
The appealing theoretical features of the MSSM have made it, through the years, a natural favorite among the theoretical frameworks incorporating a possible DM particle. In this review, we have given a summary of the current status of phenomenological constraints on the DM candidates of the MSSM and have highlighted the growing consensus that, although available parameter space remains open for most DM aspirant particles, only one of them, the higgsino-like neutralino, is almost entirely free of tension from the increasing amount of observational data.
Much of what makes higgsinos very attractive is the fact that the current constraints are not evaded with specific arrangements of some model parameters, but rather as a consequence only of the higgsino isospin quantum numbers, which lead to a fairly large mass to produce Ωh2 in agreement with observations, and of the mass splittings among its neutral and charged components, which stem directly from EWSB. As these are not exotic features, one reasonably expects that the higgsino parameter space will not remain unexplored indefinitely.
We have thus reviewed the excellent prospects for detection of higgsinos in the traditional experimental venues of direct DM detection in underground searches, indirect detection from astrophysical observations, and collider accelerators, all of which show reasons for optimism. The prospects are particularly enticing in supergravity-inspired scenarios with radiative EWSB, where the overall consistency of the theoretical picture requires a lower bound on the spin-independent cross section for higgsinos, determined indirectly but convincingly by the measured value of the Higgs boson mass.
For those models that might instead be characterized by very large scales for the superpartners (in agreement with the 125GeV Higgs mass when tanβ is close to 1), the prospects for detection are more tricky to assess, but not without hope. We have drawn the reader’s attention to a few references that promoted alternative venues for the explorations of this more fleeting scenarios. Promising venues are given by the experimental determination of dipole moments, disappearing track signatures in colliders, and the measurement of cooling curves in white dwarfs and neutron stars.
Overall, we hope this might serve as an agile but comprehensive report on the consistency of the higgsino DM picture, and on the multiple opportunities that arise for its observation in the not so distant future.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
The authors would like to thank Luc Darmé for his comments on the manuscript and discussions. The use of the CIS computer cluster at the National Centre for Nuclear Research in Warsaw is gratefully acknowledged.
KolbE. W.TurnerM. S.The early universe1990691547GondoloP.GelminiG.Cosmic abundances of stable particles: improved analysis1991360114517910.1016/0550-3213(91)90438-4JungmanG.KamionkowskiM.GriestK.Supersymmetric dark matter19962675-619537310.1016/0370-1573(95)00058-52-s2.0-0030103124DodelsonS.2003Amsterdam, the NetherlandsAcademic Presshttp://www.slac.stanford.edu/spires/find/books/www?cl=QB981:D62:2003WMAP CollaborationKomatsuE.Seven-year wilkinson microwave anisotropy probe (WMAP) observations: cosmological interpretation201119218Planck CollaborationAdeP. A. R.Planck 2015 results. XIII. Cosmological parameters2016549A13https://indico.cern.ch/event/653848/contributions/2719894/attachments/1576222/2489420/CouncilTalk.pdfhttps://indico.cern.ch/event/653848/contributions/2719895/attachments/1576169/2489824/CMSCERN-LHC25.pdfMartinS. P.A Supersymmetry primer19981819810.1142/9789812839657_0001GelminiG. B.Light weakly interacting massive particles201780808220110.1088/1361-6633/aa6e5cArcadiG.DutraM.GhoshP.LindnerM.MambriniY.PierreM.ProfumoS.QueirozF. S.The waning of the WIMP? A review of models, searches, and constraints201878320310.1140/epjc/s10052-018-5662-yPlehnT.Yet another introduction to dark matter2017RoszkowskiL.SessoloE. M.TrojanowskiS.WIMP dark matter candidates and searches—current status and future prospects201881606620110.1088/1361-6633/aab913BaerH.BargerV.SerceH.SUSY under siege from direct and indirect WIMP detection experiments2016941110.1103/PhysRevD.94.115019KrallR.ReeceM.Last electroweak WIMP standing: pseudo-dirac higgsino status and compact stars as future probes201842410.1088/1674-1137/42/4/043105StarkmanG. D.GouldA.EsmailzadehR.DimopoulosS.Opening the window on strongly interacting dark matter19904112359436032-s2.0-000040826010.1103/PhysRevD.41.3594SmithP.BennettJ.A search for heavy stable particles1979149352553310.1016/0550-3213(79)90006-3FarrarG. R.FayetP.Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry197876557557910.1016/0370-2693(78)90858-4DimopoulosS.GeorgiH.Softly broken supersymmetry and SU(5)1981193115016210.1016/0550-3213(81)90522-8WeinbergS.Supersymmetry at ordinary energies. 1. Masses and conservation laws19822628710.1103/physrevd.26.287SakaiN.YanagidaT.Proton decay in a class of supersymmetric grand unified models1982197353354210.1016/0550-3213(82)90457-6DimopoulosS.RabyS.WilzeckF.Proton decay in supersymmetric models1982112213313610.1016/0370-2693(82)90313-6BarbierR.R-parity violating supersymmetry20054201202EllisJ. R.HagelinJ. S.NanopoulosD. V.OliveK. A.SrednickiM.Supersymmetric relics from the big bang1984238245347610.1016/0550-3213(84)90461-92-s2.0-0000476398GriestK.Cross sections, relic abundance, and detection rates for neutralino dark matter1988382357Erratum: Physical Review D, vol. 39, pp. 3802, 198910.1103/PhysRevD.39.38022-s2.0-33744577642GriestK.Calculations of rates for direct detection of neutralino dark matter19886166666692-s2.0-1754438018110.1103/PhysRevLett.61.666GriestK.KamionkowskiM.TurnerM. S.Supersymmetric dark matter above the w mass199041123565358210.1103/PhysRevD.41.3565TremaineS.GunnJ. E.Dynamical role of light neutral leptons in cosmology19794264074102-s2.0-424331742510.1103/PhysRevLett.42.407WhiteS. D.FrenkC. S.DavisM.Clustering in a neutrino-dominated universe1983274L1L510.1086/184139AbazajianK.Linear cosmological structure limits on warm dark matter200673610.1103/PhysRevD.73.063513De PutterR.MenaO.GiusarmaE.HoS.CuestaA.SeoH.RossA. J.WhiteM.BizyaevD.BrewingtonH.KirkbyD.MalanushenkoE.MalanushenkoV.OravetzD.PanK.PercivalW. J.RossN. P.SchneiderD. P.SheldenA.SimmonsA.SneddenS.New neutrino mass bounds from sloan digital sky survey III data release 8 photometric luminous galaxies20127611210.1088/0004-637X/761/1/12LukashV. N.MikheevaE. V.MalinovskyA. M.Formation of the large-scale structure of the Universe20115410983100510.3367/UFNe.0181.201110a.1017DreesM.NojiriM. M.The Neutralino relic density in minimal N = 1 supergravity19934737640810.1103/PhysRevD.47.376BaerH.BrhlikM.Cosmological relic density from minimal supergravity with implications for collider physics199653259760510.1103/PhysRevD.53.597Arkani-HamedN.DelgadoA.GiudiceG. F.The well-tempered neutralino20067411-210813010.1016/j.nuclphysb.2006.02.010CMS CollaborationSirunyanA. M.Search for new phenomena with the MT2 variable in the all-hadronic final state produced in proton-proton collisions at s=13 TeV20177710710ATLAS CollaborationAaboudM.Search for top-squark pair production in final states with one lepton, jets, and missing transverse momentum using 36 fb−1 of s=13 TeV pp collision data with the ATLAS detector2017ATLAS CollaborationAaboudM.Search for squarks and gluinos in final states with jets and missing transverse momentum using 36 fb−1 of s=13 TeV pp collision data with the ATLAS detector201897Particle Data Group CollaborationPatrignaniC.Review of particle physics20164010FukushimaK.KelsoC.KumarJ.SandickP.YamamotoT.MSSM dark matter and a light slepton sector: the incredible bulk201490910.1103/PhysRevD.90.095007ATLAS CollaborationAadG.Search for the electroweak production of supersymmetric particles in s=8 TeV pp collisions with the ATLAS detector2016935ATLAS CollaborationProspect for a search for direct stau production in events with at least two hadronic taus and missing transverse momentum at theHigh Luminosity LHC with the ATLAS Detector2016Geneva, SwitzerlandCERNhttp://cds.cern.ch/record/2220805GriestK.SeckelD.Three exceptions in the calculation of relic abundances19914310319132032-s2.0-000063354510.1103/PhysRevD.43.3191EllisJ.FalkT.OliveK. A.Neutralino-stau coannihilation and the cosmological upper limit on the mass of the lightest supersymmetric particle19984443-436737210.1016/S0370-2693(98)01392-6EllisJ.FalkT.OliveK. A.SrednickiM.Calculations of neutralino–stau coannihilation channels and the cosmologically relevant region of MSSM parameter space2000132-3181213Erratum: Astroparticle Physics, vol. 15, article 413, 200110.1016/S0927-6505(99)00104-8KaneG. L.KoldaC. F.RoszkowskiL. G.WellsJ. D.Study of constrained minimal supersymmetry199449116173621010.1103/physrevd.49.61732-s2.0-33744734726RoszkowskiL.SessoloE. M.WilliamsA. J.What next for the CMSSM and the NUHM: improved prospects for superpartner and dark matter detection201420140806710.1007/JHEP08(2014)067BechtleP.Killing the cMSSM softly201676296HanC.HikasaK.WuL.YangJ. M.ZhangY.Status of CMSSM in light of current LHC Run-2 and LUX data201776947047610.1016/j.physletb.2017.04.026GAMBIT CollaborationAthronP.Global fits of GUT-scale SUSY models with GAMBIT20177712824EllisJ.EvansJ. L.LuoF.OliveK. A.ZhengJ.Stop coannihilation in the CMSSM and SubGUT models201878510.1140/epjc/s10052-018-5831-zAkulaS.NathP.Gluino-driven radiative breaking, Higgs boson mass, muon g-2, and the Higgs diphoton decay in supergravity unification2013871110.1103/PhysRevD.87.115022RoszkowskiL.SessoloE. M.WilliamsA. J.Prospects for dark matter searches in the pMSSM20152015201410.1007/JHEP02(2015)014MohapatraR. N.Theories of neutrino masses and mixings1999Arkani-HamedN.HallL.MurayamaH.SmithD.WeinerN.Small neutrino masses from supersymmetry breaking2001641110.1103/PhysRevD.64.115011BorzumatiF.NomuraY.Low-scale seesaw mechanisms for light neutrinos200164510.1103/PhysRevD.64.053005Tucker-SmithD.WeinerN.Inelastic dark matter200164Tucker-SmithD.WeinerN.The Status of inelastic dark matter200572AsakaT.IshiwataK.MoroiT.Right-handed sneutrino as cold dark matter200673510.1103/PhysRevD.73.051301ArinaC.FornengoN.Sneutrino cold dark matter, a new analysis: relic abundance and detection rates200720071102910.1088/1126-6708/2007/11/029ArinaC.CatalanM. E.KramlS.KulkarniS.LaaU.Constraints on sneutrino dark matter from LHC Run 120152015514210.1007/JHEP05(2015)142DreesM.NojiriM. M.RoyD. P.YamadaY.Light Higgsino dark matter199756276290Erratum: Physical Review D, vol. 64, article 039901, 200110.1103/PhysRevD.64.039901NagataN.ShiraiS.Higgsino dark matter in high-scale supersymmetry20152015102910.1007/JHEP01(2015)029HisanoJ.MatsumotoS.NojiriM. M.Unitarity and higher-order corrections in neutralino dark matter annihilation into two photons200367710.1103/PhysRevD.67.075014HisanoJ.MatsumotoS.NojiriM. M.SaitoO.Nonperturbative effect on dark matter annihilation and gamma ray signature from the galactic center200571610.1103/PhysRevD.71.063528HisanoJ.MatsumotS.NagaiM.SaitoO.SenamiM.Non-perturbative effect on thermal relic abundance of dark matter20076461343810.1016/j.physletb.2007.01.012CirelliM.StrumiaA.TamburiniM.Cosmology and astrophysics of minimal dark matter20077871-215217510.1016/j.nuclphysb.2007.07.023HryczukA.IengoR.UllioP.Relic densities including Sommerfeld enhancements in the MSSM20112011306910.1007/JHEP03(2011)069Zbl1301.81346ProfumoS.YagunaC. E.Statistical analysis of supersymmetric dark matter in the minimal supersymmetric standard model after WMAP200470910.1103/PhysRevD.70.095004BaerH.LessaA.RajagopalanS.SreethawongW.Mixed axion/neutralino cold dark matter in supersymmetric models201111060603110.1088/1475-7516/2011/06/031BaerH.LessaA.SreethawongW.Coupled Boltzmann calculation of mixed axion/neutralino cold dark matter production in the early universe201212010103610.1088/1475-7516/2012/01/036FalkT.OliveK. A.SrednickiM.Heavy sneutrinos as dark matter1994339324825110.1016/0370-2693(94)90639-4HallL. J.MoroiT.MurayamaH.Sneutrino cold dark matter with lepton-number violation19984243-430531210.1016/S0370-2693(98)00196-8PandaX-II CollaborationTanA.Dark matter results from first 98.7 days of data from the PandaX-II experiment201611712XENON CollaborationAprileE.First dark matter search results from the XENON1T experiment201711918LUX CollaborationAkeribD. S.Results from a search for dark matter in the complete LUX exposure20171182DEAP CollaborationAmaudruzP. A.DEAP-3600 dark matter search2016273–275340346XENON CollaborationAprileE.Physics reach of the XENON1T dark matter experiment2016160404027AalsethC. E.The darkside multiton detector for the direct dark matter search20152015854136210.1155/2015/541362LUX CollaborationLZ CollaborationSzydagisM.The present and future of searching for dark matter with LUX and LZ220Proceedings of the 38th International Conference on High Energy PhysicsAugust 2016Chicago, IL, USA10.22323/1.282.0220DARWIN CollaborationAalbersJ.DARWIN: towards the ultimate dark matter detector20161611017ATLAS CollaborationAaboudM.Search for electroweak production of supersymmetric states in scenarios with compressed mass spectra at s=13 TeV with the ATLAS detector2018975CMS CollaborationSearch for new physics in events with two low momentum opposite-sign leptons and missing transverse energy at s=13TeV201716-048Geneva, SwitzerlandCERNIceCube CollaborationAartsenM. G.Search for annihilating dark matter in the Sun with 3 years of IceCube data2017773146ANTARES CollaborationAdrian-MartinezS.Limits on dark matter annihilation in the sun using the ANTARES neutrino telescope2016759697410.1016/j.physletb.2016.05.019LZ CollaborationAkeribD. S.LUX-ZEPLIN (LZ) conceptual design report2015https://indico.cern.ch/event/432527/contributions/1071434/attachments/1320962/1980949/PICO_Dark_Matter_Searches_ICHEP_2016.pdfFermi-LAT CollaborationMAGIC CollaborationAhnenM. L.Limits to dark matter annihilation cross-section from a combined analysis of MAGIC and Fermi-LAT observations of dwarf satellite galaxies201620160203910.1088/1475-7516/2016/02/039H.E.S.S. CollaborationAbdallahH.Search for dark matter annihilations towards the inner Galactic halo from 10 years of observations with H.E.S.S201611711AMS CollaborationAguilarM.Antiproton flux, antiproton-to-proton flux ratio, and properties of elementary particle fluxes in primary cosmic rays measured with the alpha magnetic spectrometer on the international space station20161179CuocoA.HeisigJ.KorsmeierM.KrämerM.Constraining heavy dark matter with cosmic-ray antiprotons201818040400410.1088/1475-7516/2018/04/004H.E.S.S. CollaborationRinchiusoL.MoulinE.VianaA.Van EldikC.VehJ.Dark matter gamma-ray line searches toward the Galactic Center halo with H.E.S.S I893Proceedings of the 35th International Cosmic Ray Conference (ICRC '17)2017RandallL.SundrumR.Out of this world supersymmetry breaking199955779118GiudiceG. F.RattazziR.LutyM. A.MurayamaH.Gaugino mass without singlets199819981202710.1088/1126-6708/1998/12/027CirelliM.FornengoN.StrumiaA.Minimal dark matter20067531-21781942-s2.0-3374779797710.1016/j.nuclphysb.2006.07.012CohenT.LisantiM.PierceA.SlatyerT. R.Wino dark matter under siege201313101006110.1088/1475-7516/2013/10/061FanJ.ReeceM.In wino veritas? Indirect searches shed light on neutralino dark matter201320131012410.1007/JHEP10(2013)124HryczukA.CholisI.IengoR.TavakoliM.UllioP.Indirect detection analysis: wino dark matter case study201414070703110.1088/1475-7516/2014/07/031BenekeM.BharuchaA.HryczukA.RecksiegelS.Ruiz-FemeníaP.The last refuge of mixed wino-Higgsino dark matter20172017110.1007/JHEP01(2017)002Zbl1373.85006H.E.S.S. CollaborationAbramowskiA.Search for photon-linelike signatures from dark matter annihilations with H.E.S.S2013110ChanK. L.ChattopadhyayU.NathP.“Naturalness, weak scale supersymmetry and the prospect for the observation of supersymmetry at the Tevatron and at the CERN LHC199858910.1103/PhysRevD.58.096004FengJ. L.MatchevK. T.MoroiT.Focus points and naturalness in supersymmetry200061710.1103/PhysRevD.61.075005FengJ. L.MatchevK. T.WilczekF.Neutralino dark matter in focus point supersymmetry2000482438839910.1016/S0370-2693(00)00512-8BaerH.MustafayevA.ParkE.TataX.Target dark matter detection rates in models with a well-tempered neutralino20072007070101710.1088/1475-7516/2007/01/017BaerH.MustafayevA.ParkE.TataX.Collider signals and neutralino dark matter detection in relic-density-consistent models without universality200820080505810.1088/1126-6708/2008/05/058BadziakM.OlechowskiM.SzczerbiakP.Is well-tempered neutralino in MSSM still alive after 2016 LUX results?201777022623510.1016/j.physletb.2017.04.059El KheishenM. M.ShafikA. A.AboshoushaA. A.Analytic formulas for the neutralino masses and the neutralino mixing matrix1992454345434810.1103/physrevd.45.43452-s2.0-0001683203ChoiS.KalinowskiJ.Moortgat–PickG.ZerwasP.Analysis of the neutralino system in supersymmetric theories200122356357910.1007/s100520100808ChoiS. Y.DreesM.GaissmaierB.Systematic study of the impact of CP violating phases of the MSSM on leptonic high-energy observables200470110.1103/PhysRevD.70.014010BeylinV.KuksaV.PasechnikR.VereshkovG.Diagonalization of the neutralino mass matrix and boson–neutralino interaction200856339540510.1140/epjc/s10052-008-0660-0BelangerG.BoudjemaF.PukhovA.SemenovA.micrOMEGAs_3: a program for calculating dark matter observables2014185396098510.1016/j.cpc.2013.10.016DreesM.NojiriM. M.Neutralino-nucleon scattering reexamined19934883483350110.1103/PhysRevD.48.3483HisanoJ.MatsumotoS.NojiriM. M.SaitoO.Direct detection of the Wino and Higgsino-like neutralino dark matter at one-loop level200571110.1103/PhysRevD.71.015007EllisJ.RoszkowskiL.LalakZ.Higgs effects on the relic supersymmetric particle density19902453-454555510.1016/0370-2693(90)90688-3CalibbiL.LindertJ. M.OtaT.TakanishiY.LHC tests of light neutralino dark matter without light sfermions201420141110610.1007/JHEP11(2014)106KowalskaK.RoszkowskiL.SessoloE. M.Two ultimate tests of constrained supersymmetry2013201367810.1007/JHEP06(2013)078ArbeyA.BattagliaM.MahmoudiF.Supersymmetric heavy Higgs bosons at the LHC201388110.1103/PhysRevD.88.015007LahanasA. B.NanopoulosD. V.SpanosV. C.Neutralino relic density in a universe with a nonvanishing cosmological constant200062210.1103/PhysRevD.62.023515EllisJ.FalkT.GanisG.OliveK. A.SrednickiM.The CMSSM parameter space at large tanβ20015101–423624610.1016/S0370-2693(01)00541-XGiudiceG. F.HanT.WangK.WangL.-T.Nearly degenerate gauginos and dark matter at the LHC20108110.1103/PhysRevD.81.115011115011SchwallerP.ZuritaJ.Compressed electroweakino spectra at the LHC201420143, article 06010.1007/JHEP03(2014)060FujiiK.The potential of the ILC for discovering new particles2017Cabrera-CatalanM. E.AndoS.WenigerC.ZandanelF.Indirect and direct detection prospect for TeV dark matter in the nine parameter MSSM201592310.1103/PhysRevD.92.035018PorodW.SPheno, a program for calculating supersymmetric spectra, SUSY particle decays and SUSY particle production at e+e− colliders2003153227531510.1016/S0010-4655(03)00222-4PorodW.StaubF.SPheno 3.1: extensions including flavour, CP-phases and models beyond the MSSM2012183112458246910.1016/j.cpc.2012.05.0212-s2.0-84864488086BechtleP.HeinemeyerS.StålO.StefaniakT.WeigleinG.HiggsSignals: confronting arbitrary higgs sectors with measurements at the tevatron and the LHC2014742, article 271110.1140/epjc/s10052-013-2711-42-s2.0-84897617176BechtleP.BreinO.HeinemeyerS.WeigleinG.WilliamsK. E.HiggsBounds: confronting arbitrary higgs sectors with exclusion bounds from LEP and the tevatron2010181113816710.1016/j.cpc.2009.09.003Zbl1205.82001BechtleP.BreinO.HeinemeyerS.WeigleinG.WilliamsK. E.HiggsBounds 2.0.0: confronting neutral and charged Higgs sector predictions with exclusion bounds from LEP and the Tevatron2011182122605263110.1016/j.cpc.2011.07.0152-s2.0-80052322410BechtleP.BreinO.HeinemeyerS.StålO.StefaniakT.WeigleinG.WilliamsK. E.HiggsBounds-4: improved tests of extended Higgs sectors against exclusion bounds from LEP, the Tevatron and the LHC20147431322-s2.0-8489459982310.1140/epjc/s10052-013-2693-2ATLASCMS CollaborationAadG.Combined measurement of the higgs boson mass in pp collisions at s=7 and 8 TeV with the ATLAS and CMS experiments201511410.1063/1.4742082VegaJ. P.VilladoroG.SusyHD: higgs mass determination in supersymmetry20152015715910.1007/JHEP07(2015)159AthronP.ParkJ.SteudtnerT.StöckingerD.VoigtA.Precise higgs mass calculations in (non-)minimal supersymmetry at both high and low scales20172017107910.1007/JHEP01(2017)079StaubF.PorodW.Improved predictions for intermediate and heavy supersymmetry in the MSSM and beyond201777510.1140/epjc/s10052-017-4893-7KowalskaK.RoszkowskiL.SessoloE. M.TrojanowskiS.Low fine tuning in the MSSM with higgsino dark matter and unification constraints20142014416610.1007/JHEP04(2014)166AllanachB. C.SOFTSUSY: a program for calculating supersymmetric spectra2002143330533110.1016/s0010-4655(01)00460-xZbl1009.815882-s2.0-0036498238HahnT.HeinemeyerS.HollikW.RzehakH.WeigleinG.High-precision predictions for the light CP-even higgs boson mass of the minimal supersymmetric standard model20141121410.1103/physrevlett.112.141801HallL. J.NomuraY.Spread supersymmetry201201082MR2949296FoxP. J.KribsG. D.MartinA.Split dirac supersymmetry: an ultraviolet completion of higgsino dark matter201490710.1103/PhysRevD.90.075006BenakliK.DarméL.GoodsellM. D.(O)Mega split201520151110.1007/JHEP11(2015)100HillR. J.SolonM. P.WIMP-nucleon scattering with heavy WIMP effective theory2014112MahbubaniR.SchwallerP.ZuritaJ.Closing the window for compressed dark sectors with disappearing charged tracks20172017611910.1007/JHEP10(2017)061FukudaH.NagataN.OtonoH.ShiraiS.Higgsino dark matter or not: role of disappearing track searches at the LHC and future colliders201878130631110.1016/j.physletb.2018.03.088MR3814469CurtinD.DeshpandeK.FischerO.ZuritaJ.Physics opportunities for long-lived particles at electron-proton colliders2017CTA CollaborationCarrJ.Prospects for indirect dark matter searches with the cherenkov telescope array (CTA)1203Proceedings of the 34th International Cosmic Ray Conference (ICRC '15)2016The Hague, The NetherlandsSilverwoodH.WenigerC.ScottP.BertoneG.A realistic assessment of the CTA sensitivity to dark matter annihilation201520150305510.1088/1475-7516/2015/03/055EllisJ.EnqvistK.NanopoulosD.ZwirnerF.Observables in low-energy superstring models198611576910.1142/S0217732386000105BarbieriR.GiudiceG. F.Upper bounds on supersymmetric particle masses19883061637610.1016/0550-3213(88)90171-xRossG. G.Schmidt-HobergK.StaubF.Revisiting fine-tuning in the MSSM2017201732110.1007/JHEP03(2017)021RossG. G.Schmidt-HobergK.StaubF.On the MSSM Higgsino mass and fine tuning201675911011410.1016/j.physletb.2016.05.053Zbl1367.81134BarnardJ.MurnaneD.WhiteM.WilliamsA. G.Constraining fine tuning in composite higgs models with partially composite leptons201720174910.1007/JHEP09(2017)049