Improved $(g-2)_\mu$ Measurements and Wino/Higgsino Dark Matter

The electroweak (EW) sector of the Minimal Supersymmetric Standard Model (MSSM) can account for a variety of experimental data. In particular, it can explain the persistent 3-4 sigma discrepancy between the experimental result for the anomalous magnetic moment of the muon and its Standard Model (SM) prediction. The lightest supersymmetric particle (LSP), which we take as the lightest neutralino, can furthermore account for the observed Dark Matter (DM) content of the universe via coannihilation with the next-to-LSP (NLSP), while being in agreement with negative results from Direct Detection (DD) experiments. Concerning the unsuccessful searches for EW superparticles at the LHC, owing to relatively small production cross-sections, a comparably light EW sector of the MSSM is in full agreement with the experimental data. The DM relic density can fully be explained by a mixed bino/wino LSP. Here we take the relic density as an upper bound, which opens up the possibility of wino and higgsino DM. We first analyze which mass ranges of neutralinos, charginos and scalar leptons are in agreement with all experimental data, including relevant LHC searches. We find roughly an upper limit of ~ 600 GeV for the LSP and NLSP masses. In a second step we assume that the new result of the Run 1 of the 'MUON G-2' collaboration at Fermilab yields a precision comparable to the existing experimental result with the same central value. We analyze the potential impact of the combination of the Run 1 data with the existing muon g-2 data on the allowed MSSM parameter space. We find that in this case the upper limits on the LSP and NLSP masses are substantially reduced by roughly 100 GeV. We interpret these upper bounds in view of future HL-LHC EW searches as well as future high-energy electron-positron colliders, such as the ILC or CLIC.


Introduction
existing data, such as the DD bounds and the EW searches at the LHC. The derived upper limits on the EW masses are discussed in the context of the upcoming searches at the HL-LHC as well as at possible future e + e − colliders, such as the ILC [36,37] or CLIC [37,38].
2 The electroweak sector of the MSSM In our notation for the MSSM we follow exactly Ref. [15]. Here we retrict ourselves to a very short introduction of the relevant parameters and symbols of the EW sector of the MSSM, consisting of charginos, neutralinos and scalar leptons. The scalar quark sector is assumed to be heavy and not to play a relevant role in our analysis. Throughout this paper we also assume that all parameters are real, i.e. the absence of CP-violation.
The masses and mixings of the neutralinos are determined (besides SM parameters)) by U (1) Y and SU (2) L gaugino masses M 1 and M 2 , the Higgs mixing parameter µ and tan β, the ratio of the two vacuum expectation values (vevs) of the two Higgs doublets of MSSM, tan β = v 2 /v 1 . After diagonalization, the four eigenvalues of the matrix give the four neutralino masses mχ0 1 < mχ0 2 < mχ0 3 < mχ0 4 . The masses and mixings of the charginos are determined (besides SM parameters) by M 2 , µ and tan β. Diagonalizing the mass matrix two chargino-mass eigenvalues mχ± 1 < mχ± 2 can be obtained. For the sleptons, as in Ref. [15], we choose common soft SUSY-breaking parameters for all three generations. The charged slepton mass matrix are determined (besides SM parameters) by the diagonal soft SUSY-breaking parameters m 2 l L and m 2 l R and the trilinear coupling A l (l = e, µ, τ ), where the latter are taken to be zero. Mixing between the "lefthanded" and "right-handed" sleptons is only relevant for scalar taus, where the off-diagonal entry in the mass matrix is given by −m τ µ tan β. Thus, for the first two generations, the mass eigenvalues can be approximated as ml 1 ml L , ml 2 ml R . In general we follow the convention thatl 1 (l 2 ) has the large "left-handed" ("right-handed") component. Besides the symbols equal for all three generations, we also explicitly use the scalar electron, muon and tau masses, mẽ 1,2 , mμ 1,2 and mτ 1,2 . The sneutrino and slepton masses are connected by the usual SU(2) relation.
Overall, the EW sector at the tree level can be described with the help of six parameters: M 1 , M 2 , µ, tan β, ml L , ml R . Throughout our analysis we neglect CP-violation and assume µ, M 1 , M 2 > 0. In Ref. [15] it was shown that choosing these parameters positive covers the relevant parameter space once the (g − 2) µ results are taken into account (see, however, the discussion in Sect. 7).
Following the stronger experimental limits from the LHC [39, 40], we assume that the colored sector of the MSSM is sufficiently heavier than the EW sector, and does not play a role in this analysis. For the Higgs-boson sector we assume that the radiative corrections to the light CP-even Higgs boson (largely originating from the top/stop sector) yield a value in agreement with the experimental data, M h ∼ 125 GeV. This naturally yields stop masses in the TeV range [41,42], in agreement with the above assumption. Concerning the Higgsboson mass scale, as given by the CP-odd Higgs-boson mass, M A , we employ the existing experimental bounds from the LHC. In the combination with other data, this results in a mostly non-relevant impact of the heavy Higgs bosons on our analysis, as will be discussed below.

Relevant constraints
The experimental result for a µ := (g − 2) µ /2 is dominated by the measurements made at the Brookhaven National Laboratory (BNL) [43], resulting in a world average of [44] a exp µ = 11659209.1(5.4)(3.3) × 10 −10 , where the first uncertainty is statistical and the second systematic. The SM prediction of a µ is given by [45] (based on Refs. [1,2,[46][47][48][49][50][51][52][53][54][55][56][57][58][59][60][61][62][63] ) 2 , Comparing this with the current experimental measurement in Eq. (1) results in a deviation of ∆a old µ = (28.1 ± 7.6) × 10 −10 , corresponding to a 3.7 σ discrepancy. This "current" result will be used below with a hard cut at 2 σ uncertainty. Efforts to improve the experimental result at Fermilab by the "MUON G-2" collaboration [4] and at J-PARC [64] aim to reduce the experimental uncertainty by a factor of four compared to the BNL measurement. For the second step in our analysis we consider the upcoming Run 1 result from the Fermilab experiment [3]. The Run 1 data is expected to have roughly the same experimental uncertainty as the current result in Eq. (1). We furthermore assume that the Run 1 data yields the same central value as the current result. Consequently, we anticipate that the experimental uncertainty shrinks by 1/ √ 2, yielding a future value of ∆a fut µ = (28.1 ± 6.2) × 10 −10 , corresponding to a 4.5 σ discrepancy. Thus, the combination of Run 1 data with the existing experimental (g − 2) µ data has the potential to (nearly) establish the "discovery" of BSM physics. This "anticipated future" result will be used below with a hard cut at 2 σ uncertainty.
Recently a new lattice calculation for the leading order hadronic vacuuum polarization (LO HVP) contribution to a SM µ [65] has been reported, which, however, was not used in the new theory world average, Eq. (2) [45]. Consequently, we also do not take this result into account, see also the discussions in Ref. [15,[65][66][67][68][69]. On the other hand, we are also aware that our conclusions would change substantially if the result presented in [65] turned out to be correct.

Other constraints
All other experimental constraints are taken into account exactly as in Ref. [15]. These comprise • Vacuum stability constraints: All points are check to possess a stable and correct EW vacuum, e.g. avoiding charge and color breaking minima. This check is performed with the public code Evade [79,80].
• Constraints from the LHC: All relevant EW SUSY searches are taken into account, mostly via CheckMATE [81][82][83], where many analysis had to be implemented newly [15]. As in Ref. [15], the constraints coming from "compressed spectra" searches [84], corresponding to very low splittings between mχ± 1 , mχ0 2 , ml 1 and mχ0 1 are applied directly on our parameter space. In addition to the searches described in Ref. [15], we take into account the latest constraints from the disappearing track searches at the LHC [85,86]. These are particularly important for wino DM scenario where the mass gap betweenχ ± 1 andχ 0 1 can be ∼ a few hundred MeV. The long-livedχ ± 1 (lifetime ∼ O(ns)) decays into final states involving aχ 0 1 and a soft pion which can not be reconstructed within the detector. Thus, the signal involves a charged track fromχ ± 1 that produces hits only in the innermost layers of the detector with no subsequent hits at larger radii.
• Dark matter relic density constraints: We use the latest result from Planck [5].
As stressed above, we take the relic density as an upper limit (evaluated from the central value plus 2 σ. The relic density in the MSSM is evaluated with MicrOMEGAs [87][88][89][90]. An additional DM component could be, e.g., a SUSY axion [91], which would then bring the total DM density into agreement with the Planck measurement of Ω CDM h 2 = 0.120 ± 0.001 [5]. In the case of wino DM, because of the extremely small mass splitting, the effect of "Sommerfeld enhancement" [92] can be very important. For wino DM providing the full amount of DM it shifts the allowed range of mχ0 1 from ∼ 2.0 TeV to about ∼ 2.9 TeV. Since here we are interested in the case that the wino DM only gives a fraction of the whole DM relic density, see Eq. (7), we can safely neglect the Sommerfeld enhancement.
The upper limit on mχ0 1 is given, as will be shown below, by the (g −2) µ constraint, but not by the DM relic density. Allowing higher masses here (as would be the case if the Sommerfeld enhancement had been taken into account) could thus not lead to a larger allowed parameter space. On the other hand, for a point with a relic density fulfilling Eq. (7) the Sommerfeld enhancement would only lower the "true" DM density, which still fulfills Eq. (7).
• Direct detection constraints of Dark matter: We employ the constraint on the spin-independent DM scattering cross-section σ SI p from XENON1T [8] experiment, evaluating the theoretical prediction for σ SI p using MicrOMEGAs [87][88][89][90]. A combination with other DD experiments would yield only very slightly stronger limits, with a negligible impact on our results. For parameter points with Ωχh 2 ≤ 0.118 (2 σ lower limit from Planck [5]), we scale the cross-section with a factor of (Ωχh 2 /0.118) to account for the fact thatχ 0 1 provides only a fraction of the total DM relic density of the universe. Here the effect of neglecting the Sommerfeld enhancement leads to a more conservative allowed region of parameter space.

Parameter scan
We scan the relevant MSSM parameter space to obtain lower and upper limits on the relevant neutralino, chargino and slepton masses. In order to achieve a "correct" DM relic density, see Eq. (7), by the lightest neutralino,χ 0 1 , some mechanism such as a specific co-annihilation or pole annihilation has to be active in the early universe. At the same time mχ0 1 must not be too high, such that the EW sector can provide the contribution required to bring the theory prediction of a µ into agreement with the experimental measurement, see Sect. 3. The combination of these two requirements yields the following possibilities. (The cases present a certain choice of favored possibilities, upon which one can expand, as will briefly discussed in Sect. 7.)

(B) Wino DM
This scenario is characterized by a small value of M 2 . Such a scenario is also naturally realized in Anomaly Mediation SUSY breaking (see e.g. Ref. [100] and references therein). We scan the following parameters: The choice of M 2 M 1 , µ leads (at tree-level) to a very degenerate spectrum with mχ± 1 − mχ0 1 = O(1 eV). However, this spectrum does not correspond to the on-shell (OS) masses of all six charginos and neutralinos. Since only three (soft SUSY-breaking) mass parameters are available (M 1 , M 2 and µ), only three out of the six masses can be renormalized OS. The (one-loop) shifts for the three remaining masses are obtained via the CCN i renormalization scheme [101] with i ∈ {2, 3, 4}. In a CCN i scheme thẽ χ ± 1 ,χ ± 2 andχ 0 i are chosen OS. This automatically yields a good renormalization for M 2 and µ. Theχ 0 i has to be chosen, parameter point by parameter point, such that also M 1 is renormalized well (which excludes the CCN 1 for M 2 M 1 , µ). For each point we choose theχ 0 i to be renormalized OS such that the maximum shift of all three shifted neutralino masses is minimized. We have explicitly checked for each scanned point that the such chosen CCN i indeed yields reasonably small shifts for three shifted neutralino masses, in particular for mχ0 1 . 4 Only this transition to OS masses yields a mass splitting between mχ± 1 and mχ0 1 that allows then for the decayχ ± 1 →χ 0 1 π ± .
(C) Mixed bino/wino DM This scenario has been analyzed in Ref. [15]. It can in principle be realized in three different versions corresponding to the coannihilation mechanism (see Ref. [15] for a detailed discussion). However, a larger wino component is found only forχ ± 1coannihilation. The scan parameters are chosen as, Here we choose one soft SUSY-breaking parameter for all sleptons together. While this choice should not have a relevant effect in theχ ± 1 -coannihilation case, this have an impact in the next case. In our scans we will see that the chosen lower and upper limits are not reached by the points that meet all the experimental constraints. This ensures that the chosen intervals indeed cover all the relevant parameter space.

(D) Bino DM
This scenario covers the coannihilation with sleptons, and has also been analyzed in Ref. [15]. In this scenario "accidentally" the wino component of theχ 0 1 can be nonnegligible. However, this is not a distinctive feature of this scenario. We cover the two distinct cases that either the SU(2) doublet sleptons, or the singlet sleptons are close in mass to the LSP. (D1) Case-L: SU(2) doublet 4 We thank C. Schappacher for evaluating the mass shift for our wino DM points (following Ref. [102]).
In all scans we choose flat priors of the parameter space and generate O(10 7 ) points. The mass parameters of the colored sector have been set to high values, such that the resulting SUSY particle masses are outside the reach of the LHC, the light CP-even Higgsboson is in agreement with the LHC measurements (see, e.g., Refs. [41,42]), where the concrete values are not relevant for our analysis. M A has also been set to be above the TeV scale. Consequently, we do not include explicitly the possibility of A-pole annihilation, with M A ∼ 2mχ0 1 . As we will discuss below the combination of direct heavy Higgs-boson searches with the other experimental requirements constrain this possibility substantially (see, however, also Sect. 6). Similarly, we do not consider h-or Z-pole annihilation, as such a light neutralino sector likely overshoots the (g − 2) µ contribution (see, however, the discussion in Sect. 6).

Analysis flow
The data sample is generated by scanning randomly over the input parameter range mentioned above, using a flat prior for all parameters. We use SuSpect [103] as spectrum and SLHA file generator, which yields DR values for the chargino and neutralino masses. While in most of the cases the difference between DR and on-shell (OS) masses is not phenomenologically relevant for our analysis, this is different for wino DM. Because of the very small mass difference of the DR values for mχ0 1 and mχ± 1 we take into account the transition to OS masses, as discussed above. In the next step the points are required to satisfy theχ ± 1 mass limit from LEP [104]. The SLHA output files from SuSpect are then passed as input to GM2Calc and MicrOMEGAs for the calculation of (g − 2) µ and the DM observables, respectively. The parameter points that satisfy the current (g − 2) µ constraint, Eq. (3), the DM relic density, Eq. (7), the direct detection constraints and the vacuum stability constraints as checked with Evade are then taken to the final step to be checked against the latest LHC constraints implemented in CheckMATE. The branching ratios of the relevant SUSY particles are computed using SDECAY [105] and given as input to CheckMATE.

Higgsino DM
We start our discussion with the case of higgsino DM, as discussed in Sect. 4.1. We follow the analysis flow as described in Sect. 4.2, where the vacuum stability constraints did not affect any of the LHC allowed points. We also remark that we have checked that the (possible) one-loop corrections to the chargino/neutralino masses do not change the results in a relevant way and were thus omitted (see, however, Sect. 5.2). In the following we denote the points surviving certain constraints with different colors: • grey (round): all scan points (i.e. points excluded by (g − 2) µ ).
• green (round): all points that are in agreement with (g − 2) µ , taking into account the current or anticipated future limits, see Eqs.
• blue (triangle): points that additionally obey the upper limit of the DM relic density, see Eq. (7) (i.e. points excluded by the DD constraints). In some plots below all points that pass the (g − 2) µ constraint are also in agreement with the DM relic density constraint, resulting in only blue (but no green) points to be visible.
• cyan (diamond): points that additionally pass the DD constraints, see Sect. 3 (i.e. points excluded by the LHC constraints).
• red (star): points that additionally pass the LHC constraints, see Sect. 3 (i.e. points that pass all constraints).
In Fig. 1 we show our results in the mχ0 1 -mχ± 1 plane for the current (left) and future (right) (g − 2) µ constraint, see Eqs.
(3) and (4), respectively. By definition the points are clustered in the diagonal of the plane, and mχ0 2 ≈ mχ± 1 . Starting with the (g − 2) µ constraint (green points), they all pass also the relic density constraint (dark blue) and thus are visible as dark blue points in the plot. Overall, one can observe a clear upper limits from (g − 2) µ of about 660 GeV for the current limits and about 600 GeV from the anticipated future accuracy. As mentioned above, the DM relic density constraint does not yield changes in the allowed parameter space. Applying the DD limits, on the other hand, forces the points to have even smaller mass differences betweenχ 0 1 andχ ± 1 . It has also an important impact on the upper limit, which is reduced to about 500 (480) GeV for the current (future) (g − 2) µ bounds. Applying the LHC constraints, corresponding to the "surviving" red points (stars), does not yield a further reduction from above for the current (g − 2) µ constraint, whereas for the anticipated future accuracy yields a reduction of ∼ 30 GeV. The LHC constraints (see also the discussion below) also cut always (as anticipated) points in the lower mass range, resulting in a lower limit of ∼ 120 GeV for mχ0 1 ≈ mχ± 1 . The LHC constraint which is most effective in this parameter plane is the one designed for compressed spectra, as demonstrated in Fig. 2, showing our scan points in the mχ0 1 -∆m(= mχ0 2 − mχ0 1 ) parameter plane. The color coding is as in Fig. 1. The black line indicates the bound from compressed spectra searches [84]. One can clearly observe that the compressed spectra bound sits exactly in the preferred higgsino DM region, i.e. at ∆m(= mχ0 2 − mχ0 1 ) = O(10 GeV). Other LHC constraint that is effective in this case is the bound from slepton pair production leading to dilepton and E / T in the final state [106], as will be discussed in more detail below. The bounds from the disappearing track searches [86] turn out to be ineffective in this case because of the very short lifetime of theχ ± 1 [107]. From Figs. 1, 2 one can conclude that the experimental data set an upper as well as a lower bound, yielding a clear search target for the upcoming LHC runs, and in particular for future e + e − colliders, as will be discussed in Sect. 6. In particular, this collider target gets (potentially) sharpened by the improvement in the (g − 2) µ measurements. The impact of the DD experiments is demonstrated in Fig. 3. We show the mχ0 1 -σ SI p plane for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding of the points (from yellow to dark blue) denotes M 2 /µ, whereas in red we show the points fulfilling (g − 2) µ , relic density, DD and the LHC constraints. The black line indicates the current DD limits, here taken for sake of simplicity from XENON1T [8], as discussed in Sect. 3. It can be seen that a slight downward shift of this limit, e.g. due to additional DD experimental limits from LUX [6] or PANDAX [7], would not change our results in a strong way, but only slightly reduce the upper limit on mχ0 1 . The scanned parameter space extends from large σ SI p values, given for the smallest scanned M 2 /µ values to the smallest ones, reached for the largest scanned M 2 /µ, i.e. the σ SI p constraints are particularly strong for small M 2 /µ, which can be understood as follows. The most important contribution to DM scattering comes from the exchange of a light CP-even Higgs boson in the t-channel. The corresponding hχ 0 1χ 0 1 coupling at tree level is given by [108] where we have assumed µ > 0. Thus, the coupling becomes large for µ ∼ M 2 or µ ∼ M 1 . Therefore, the XENON1T DD bound pushes the allowed parameter space into the almost pure higgsino-LSP region, with negligible bino and wino component. The impact of the compressed spectra searches is visible in the lower mχ0 1 region. Given both CDM constraints ) plane for the higgsino DM scenario for current (left) and anticipated future limits (right) from (g − 2) µ . For the color coding: see text. and the LHC constraints, shown in red, the smallest M 2 /µ value we find is 2.2 for both current and anticipated future (g − 2) µ bound. This result depends mildly on the assumed (g − 2) µ constraint, as this cuts away the largest mχ0 1 values. All of the points will be conclusively probed by the future DD experiment XENONnT [109].
The distribution of ml 1 (where it should be kept in mind that we have chosen the same masses for all three generations, see Sect. 2) is presented in the mχ0 1 -ml 1 plane in Fig. 4, with the same color coding as in Fig. 1. The (g − 2) µ constraint places important constraints in this mass plane, since both types of masses enter into the contributing SUSY diagrams, see Sect. 3. The constraint is satisfied in a roughly triangular region with its tip around (mχ0 1 , ml 1 ) ∼ (650 GeV, 700 GeV) in the case of current (g − 2) µ constraints, and around ∼ (600 GeV, 600 GeV) in the case of the anticipated future limits, i.e. the impact of the anticipated improved limits is clearly visible as an upper limit for both masses. Since no specific other requirement is placed on the slepton sector in the higgsino DM case the slepton masses are distributed over the (g − 2) µ allowed region. The DM relic density constraint, as discussed above, does not yield any further bounds on the allowed parameter space. The inclusion of the DM DD bounds, as visible by the cyan and red points, only cuts away the very largest slepton masses (for a given LSP mass).
The LHC constraints cut out all points with mχ0 1 < ∼ 125 GeV, as well as a triangular region with the tip around (mχ0 1 , ml 1 ) ∼ (340 GeV, 450 GeV). The first "cut" is due to the searches for compressed spectra. The second cut is mostly a result of the constraint coming from slepton pair production searches leading to dilepton and E / T in the final state [106]. The bound obtained by recasting the experimental search in CheckMATE is substantially weaker than the original limit from ATLAS. That limit is obtained for a "simplified model" with BR(l 1 ,l 2 → lχ 0 1 ) = 100%, an assumption which is not stritly valid in our parameter space. The small mass gap amongχ 0 1 ,χ 0 2 andχ ± 1 allows significant BR of the sleptons to final states involvingχ ± 1 andχ 0 2 . This reduces the number of signal leptons and hence weakens the exclusion limit. Overall we can place an upper limit on the light slepton mass of about ∼ 1200 GeV and 1050 GeV for the current and the anticipated future accuracy of (g − 2) µ , respectively. Since larger values of slepton masses are reached for lower values of mχ0 1 , the impact of (g − 2) µ is relatively weaker than in the case of chargino/neutralino masses.
We finish our analysis of the higgsino DM case with the mχ0 1 -tan β plane presented in Fig. 5 with the same color coding as in Fig. 1. The (g − 2) µ constraint is fulfilled in a triangular region with the largest neutralino masses allowed for the largest tan β values (where we stopped our scan at tan β = 60), following the analytic dependence of the (g − 2) µ contributions in Sect. 3, a µ ∝ tan β/m 2 EW (where we denote with m EW an overall EW mass scale. In agreement with the previous plots, the largest values for the lightest neutralino masses are ∼ 650 GeV (∼ 600 GeV) for the current (anticipated future) (g − 2) µ constraint. The DM relic density does not give any additional constraint. The points allowed by the DM DD limits (cyan and red) yield the observed reduction to about ∼ 500 GeV. The LHC constraints cut out all points at low mχ0 1 , but nearly independent of tan β. As observed before, they yield a small further reduction in the case of the anticipated future (g − 2) µ accuracy.
In Fig. 5 we also show as black lines the current bound from LHC searches for heavy  neutral Higgs bosons [110] in the channel pp → H/A → τ τ in the M 125 h (χ) benchmark scenario (based on the search data published in Ref. [111] using 139 fb −1 .) 5 . In this scenario light charginos and neutralinos are present, suppressing the τ τ decay mode and thus yielding relatively weak limits in the M A -tan β plane (see, e.g., Fig. 5 in [110]). The black lines correspond to mχ0 1 = M A /2, i.e. roughly to the requirement for A-pole annihilation, where points above the black lines are experimentally excluded. It can be observed that all points allowed by (g − 2) µ are above the exclusion curve. This renders the effects of A-pole annihilation in this scenario effectively irrelevant (and justifies our choice to fix M A above the TeV scale).

Wino DM
The next case under investigation is the wino DM case, as discussed in Sect. 4.1. We follow the analysis flow as described in Sect. 4.2 and denote the points surviving certain constraints with different colors as defined in Sect. 5.1. The vacuum stability test had no effect on the points passing all other constraints. ) plane for the wino DM scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding is as in Fig. 1.
(3) and (4), respectively. We display the results for ∆m rather than for mχ± 1 , since the mass difference is very small, and the various features are more easily visible in this plane. It should be remembered that we have applied the one-loop shift to, in particular, mχ0 1 that allows for the decayχ ± 1 →χ 0 1 π ± , see the discussion in Sect. 4.1. This results in the lower bound on ∆m of ∼ 0.15 GeV. As in the -τ plane for current (left) and anticipated future limits (right) from (g − 2) µ , where τ is the lifetime of the chargino decaying to π ±χ0 1 . The current limit from CMS [86] is shown as the black solid line. The color coding is as in Fig. 6. higgsino scenario, see Sect. 5.1, all points that pass the (g − 2) µ constraint (current and anticipated future) also pass the relic density constraint, shown as blue triangles in Fig. 6. The highest allowed chargino masses are bounded from above by the (g −2) µ constraint. The overall allowed parameter space, shown as red stars, is furthermore bounded from "above" by the DD limits and from "below" by the LHC constraints. The DD limits cut away larger mass differences, which can be understood as follows. The hχ 0 1χ 0 1 coupling for a wino-likeχ 0 1 is given by [108] in the limit of ||µ|−M 2 | M Z and a decoupled CP-odd Higgs boson (assuming also that the h-exchange dominates over the H contribution in the (spin independent) DD bounds). This coupling becomes large for µ ∼ M 2 . On the other hand, the tree level mass splitting between the wino-like statesχ ± 1 andχ 0 1 generated (mainly by the mixing of the lighter chargino with the charged higgsino) is given as [117] ∆m(= mχ± for The masss splitting increases for smaller µ values and thus coincides with larger DD cross sections, as discussed with Eq. (14). All limits together yield maximum ∆m ∼ 2(0.2) GeV for mχ± 1 ∼ 100(600) GeV for the current (g − 2) µ constraint. The upper limit is reduced to ∼ 500 GeV for the future anticipated (g − 2) µ constraint.
The relevant LHC constraint is further analyzed in Fig. 7, where we show the plane mχ± 1 -τχ± 1 . τχ± 1 denotes the lifetime of the chargino decaying to π ±χ0 1 . Overlaid as black line is the bound from (CMS) charged disappearing track analysis [86]. One can observe that this constraint, cutting out parameter points between τχ± 1 ∼ 0.01 ns and ∼ 1 ns is responsible for the main LHC exclusion. It should be noted that in the "red star area" also some points appear cyan, i.e. excluded by (other) LHC searches, where the most relevant channels are pair production of sleptons leading to two leptons and E / T in the final state [106]. However, these channels are not strong enough to exclude more of the mχ± 1 -τχ± 1 plane than the disappearing track search. It can be expected in the (near) future that improved DD bounds, cutting the allowed parameter space from smallχ ± 1 lifetime and improved disappearing track searches, cutting from large lifetimes, may substantially shrink the allowed parameter space, sharpening the upper limit on mχ± 1 and the prospects for future collider searches. The two bounds together have the potential to firmly rule out the case of wino DM in the MSSM, will be discussed below. The impact of the DD experiments is demonstrated in Fig. 8. We show the mχ0 1 -σ SI p plane for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding of the points (from yellow to dark green) denotes µ/M 2 , whereas in red we show the points fulfilling (g − 2) µ , relic density, DD and the LHC constraints. The solid black line indicates the current DD limits, here taken for sake of simplicity from XENON1T [8], as discussed in Sect. 3. It can be seen that a slight downward shift of this limit, e.g. due to additional DD experimental limits from LUX [6] or PANDAX [7], would not change our results in a relevant way. However, moderately improved limits may have a strong impact, as discussed above. The scanned parameter space extends from large σ SI p values, given for the smallest scanned µ/M 2 values to the smallest ones, reached for the largest scanned µ/M 2 , i.e. the σ SI p constraints are particularly strong for small µ/M 2 . Given both CDM constraints and the LHC constraints, shown in red, the smallest µ/M 2 value we find is 1.5 for both the current and anticipated future (g − 2) µ bound. As mentioned above, the DD bound can become relevantly stronger with future experiments. We show as dashed line the projected limit of XENONnT [109], and the dot-dashed line indicates the neutrino floor [118]. One can see that the XENONnT result will either firmly exlude or detect a wino DM candidate, possibly in conjunction with improved disappearing track searches at the LHC, as discussed above.
The distribution of the lighter slepton mass (where it should be kept in mind that we have chosen the same masses for all three generations, see Sect. 2) is presented in the mχ0 1ml 1 plane in Fig. 9, with the same color coding as in Fig. 6. The (g − 2) µ constraint places important constraints in this mass plane, since both types of masses enter into the contributing SUSY diagrams, see Sect. 3 (and obviously all points pass the DM relic density upper limit). The (g − 2) µ constraint is satisfied in a triangular region with its tip around (mχ0 1 , ml 1 ) ∼ (700 GeV, 700 GeV) in the case of current (g − 2) µ constraints, and around ∼ (600 GeV, 700 GeV) in the case of the anticipated future limits. The highest slepton masses reached are about 1500(1200) GeV, respectively; i.e. the impact of the anticipated improved limits is clearly visible as an upper limit. After including the DD and LHC constraints the upper limits for the LSP are slightly reduced by ∼ 100 GeV, whereas the upper limits on sleptons are not affected. Concerning the LHC searches, the constraints from slepton pair production searches rule out some parameter points from the low ml 1 region. Larger slepton masses are excluded by the same search if the "second" slepton turns out to be relatively light. Points excluded by compressed spectra searches [84], depending on the lightest chargino mass, can be found all over the plane.
We finish our analysis of the wino DM case with the mχ0 1 -tan β plane presented in Fig. 10 with the same color coding as in Fig. 6. The (g − 2) µ constraint is fulfilled in a triangular region with largest neutralino masses allowed for the largest tan β values (where we stopped our scan at tan β = 60), following the analytic dependence of the (g − 2) µ contributions in Sect. 3, a µ ∝ tan β/m 2 EW (where we denote with m EW an overall EW mass scale. In agreement with the previous plots, the largest values for the lightest neutralino masses are ∼ 600 GeV (∼ 500 GeV) for the current (anticipated future) (g − 2) µ constraint. The points allowed by the DM constraints (blue/cyan) are distributed all over the allowed region. The LHC constraints also cut out points distributed all over the allowed triangle, in agreement with the previous discussion.
In Fig. 10 we also show as black lines the current bound from LHC searches for heavy neutral Higgs bosons [110,111], see the discussion in Sect. 5.1. Points above the black lines are experimentally excluded. There are a few points passing the current (g − 2) µ constraint below the black A-pole line, reaching up to mχ0 1 ∼ 220 GeV, for which the A-pole annihilation could provide the correct DM relic density. For the anticipated future accuracy in (g − 2) µ this mechanism would effectively be absent, making the A-pole annihilation in this scenario marginal.

Bino/wino and bino DM
In this section we analyze the case of bino/wino and bino DM, as defined in Sect. 4.1. The three cases defined there correspond exactly to the set of analyses in Ref. [15]. However, we now apply the DM relic density as an upper bound ("DM upper bound"), whereas in Ref. [15] the LSP was required to give the full amount of CDM ("DM full"). Since overall the results are similar to the ones found in Ref. [15], we will keep the discussion brief, but try to highlight the differences w.r.t. Ref. [15]. For all scenarios we find that the vacuum stability bounds have no impact on the final mass limits found.

Bino/wino DM withχ ±
1 -coannihilation In Fig. 11 we show our results in the mχ0 1 -mχ± 1 plane for the current (left) and future (right) (g −2) µ constraint, see Eqs. (3) and (4), respectively. The color coding is defined in Sect. 5.1. By definition ofχ ± 1 -coannihilation the points are clustered in the diagonal of the plane. Overall we observe here exactly the same pattern of points as in the case of "DM full" [15]. After taking into account all constraints we find upper limits of ∼ 600(500) GeV in the case of the current (future) (g −2) µ limits. Thus, the experimental data set about the same upper as well as lower bounds as in Ref. [15] (modulo differences due to point density artefacts). Consequently, the search targets for the upcoming LHC runs particular for future e + e − colliders remain about the same. plane for the bino-winoχ ± 1coannihilation scenario for current (left) and anticipated future limits (right) from (g − 2) µ . For the color coding: see text.
The impact of the DD experiments is demonstrated in Fig. 12. We show the mχ0 1 -σ SI p plane for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding of the points (from yellow to dark green) denotes µ/M 1 , whereas in red we show the points fulfilling all constraints including the LHC ones. As in Fig. 3 the solid black line indicates the current DD limits from XENON1T [8]. Also here the results are in very good agreement with Ref. [15] (again modulo point density artefacts). The scanned parameter space extends from large σ SI p values, given for the smallest scanned µ/M 1 values to the smallest ones, reached for the largest scanned µ/M 1 , i.e. the σ SI p constraints are particularly strong for small µ/M 1 . Given in red the points fulfilling (g − 2) µ , both CDM constraints and the LHC constraints, the smallest µ/M 1 value we find is 1.67 for the current and 1.78 for the anticipated future (g − 2) µ bound. The dashed line indicates the projected XENONnT limit [109], and the dot-dashed line indicates the neutrino floor [118]. One can see that XENONnT will not be able to fully test the chargino co-annihilation scenario, with some points that pass all constraints (red) being even below the neutrino floor. -σ SI p plane for bino-winoχ ± 1 -coannihilation scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding of the points denotes µ/M 1 and the black lines indicates the DD limits (see text). In red we show the points fulfilling (g − 2) µ , the relic density, DD and the LHC constraints.
The distribution of the lighter slepton mass (where it should be kept in mind that we have chosen the same masses for all three generations, see Sect. 2) is presented in the mχ0 1 -ml 1 plane in Fig. 13, with the same color coding as in Fig. 11. The (g − 2) µ constraint places important constraints in this mass plane, since both types of masses enter into the contributing SUSY diagrams, see Sect. 3. The constraint is satisfied in a triangular region with its tip around (mχ0 1 , ml 1 ) ∼ (700 GeV, 800 GeV) in the case of current (g − 2) µ constraints, and around ∼ (600 GeV, 700 GeV) in the case of the anticipated future limits, i.e. the impact of the anticipated improved limits is clearly visible as an upper limit. These results remain unchanged w.r.t. the "DM full" case [15]. The points fulfilling the DM relic density constraint (blue/cyan/red) are distributed all over the (g − 2) µ allowed range. They extend to somewhat higher slepton masses as compared to Ref. [15], due to the less restrictive DM upper bound. 6 The DD limits cut away the largest values, in particular for ml 1 , which can be understood as follows. Large ml 1 values, with correspondingly large sneutrino masses, require smaller µ to satisfy the (g − 2) µ constraint. This in turn puts them in tension with the DD bounds, see Fig. 12. -ml 1 plane for the bino-winoχ ± 1 -coannihilation scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding is as in Fig. 11.
The LHC constraints cut out lower slepton masses, following the same pattern as in Ref. [15]. They cut away masses up to ml 1 < ∼ 450 GeV, as well as part of the very low mχ0 1 points nearly independent of ml 1 . Here the latter "cut" is due to the searches for compressed spectra withχ ± 1 ,χ 0 2 decaying via off-shell gauge bosons [84]. The first "cut" is mostly a result of the searches for slepton pair production with a decay to two leptons plus missing energy [106]. As was demonstrated and discussed in detail in Ref. [15] for this limit it is crucial to employ a proper re-cast of the LHC searches, rather than a naive application of the published bounds, Overall we can place an upper limit on the light slepton mass of about ∼ 1050 GeV and ∼ 950 GeV for the current and the anticipated future accuracy of (g − 2) µ , respectively.
We finish our analysis of theχ ± 1 -coannihilation case with the mχ0 1 -tan β plane presented in Fig. 14 with the same color coding as in Fig. 11. For this plane we find that the results are in full agreement with the "DM full" case as analyzed in Ref. [15]. The (g − 2) µ constraint is fulfilled in a triangular region with largest neutralino masses allowed for the largest tan β values (tan β = 60). In agreement with the previous plots, the largest values for the lightest neutralino masses are ∼ 600 GeV (∼ 500 GeV) for the current (anticipated future) (g − 2) µ (a) (b) Figure 14: The results of our parameter scan in the mχ0 1 -tan β plane in the bino-winoχ ± 1 -coannihilation scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding is as in Fig. 11. The black line indicates the current exclusion bounds for heavy MSSM Higgs bosons at the LHC (see text).
constraint. The points allowed by the DM constraints (blue/cyan) are distributed all over the allowed region. The LHC constraints cut out points at low mχ0 1 , but nearly independent on tan β.
As in the previous scenarios, in Fig. 14 we also show as black lines the current bound from LHC searches for heavy neutral Higgs bosons [110] in the channel pp → H/A → τ τ in the M 125 h (χ) benchmark scenario. As before the black lines correspond to mχ0 1 = M A /2, i.e. roughly to the requirement for A-pole annihilation, where points above the black lines are experimentally excluded. The improved limits of the experimental analysis based on 139fb −1 [111], but now with the relaxed DM relic density bound, still allow parameter points that cannot be regarded as excluded. They are found for mχ0 1 < ∼ 240(200) GeV and tan β < ∼ 12 in the case of the current (anticipated future) (g − 2) µ constraints. To analyze these points a dedicated analysis of the A-pole annihilation would be required. This dedicated analysis we leave for future work.

Bino DM withl ± -coannihilation case-L
We now turn to the case of bino DM withl ± -coannihilation. As discussed in in Sect. 4.1 we distinguish two cases, depending which of the two slepton soft SUSY-breaking parameters is set to be close to mχ0 1 . We start with the Case-L, where we chose ml L ∼ M 1 , i.e. the left-handed charged sleptons as well as the sneutrinos are close in mass to the LSP. We find that all six sleptons are close in mass and differ by less than ∼ 50 GeV.
In Fig. 15 we show the results of our scan in the mχ0 1 -mμ 1 plane, as before for the current  (g − 2) µ constraint (left) and the anticipated future constraint (right). The color coding of the points is the same as in Fig. 1, see the description in the beginning of Sect. 5.1. By definition of the scenario, the points are located along the diagonal of the plane. The result is qualitatively similar to the analysis in Ref. [15], where the relic density constraint was used as a direct measurement, and not only as an upper limit. However, one can observe that a slightly larger mass difference between theμ 1 and theχ 0 1 is allowed in comparison with Ref. [15]. Taking all bounds into account we find upper limits on the LSP mass of ∼ 500 GeV and ∼ 450 GeV for the current and anticipated future (g − 2) µ accuracy, respectively. The limits on mμ 1 are about ∼ 50 − 100 GeV larger.
In Fig. 16 we show the results in the mχ0 1 -mχ± 1 plane with the same color coding as in Fig. 15. The results are again qualitatively in agreement with Ref. [15]. However, in particular for the future (g −2) µ constraint, slightly larger values of mχ± 1 are allowed. Overall we find for the light charginos mass an upper limit of ∼ 3200 GeV (∼ 2900 GeV) for the current (anticipated future (g − 2) µ ) constraint. Clearly visible are the LHC constraints forχ 0 2 -χ ± 1 pair production leading to three leptons and E / T in the final state [119]. At very low values of both mχ0 1 and mχ± 1 the compressed specta searches [84] cut away points up to The results for thel ± -coannihilation Case-L in the mχ0 1 -tan β plane are presented in Fig. 17 7 . The overall picture is similar to theχ ± 1 -coannhiliation case shown above in Fig. 14. As above they are also qualitatively and quantitatively similar to the results of Case-L with 7 Clearly visible in the two plots are different point densities due to indepedent samplings that were the relic density taken as a direct measurement as presented in Ref. [15]. Larger LSP masses are allowed for larger tan β values. On the other hand the combination of small mχ0 1 and large tan β leads to a too large contribution to a SUSY µ and is thus excluded. As in the other scenarios we also show the limits from H/A searches at the LHC as a solid line. In this case for the current and anticipated future (g − 2) µ limit substantially more points passing the (g − 2) µ constraint "survive" below the black line(s), i.e. they are potential candidates for A-pole annihilation. The masses reach up to ∼ 300 GeV and ∼ 250 GeV, respectively. These limits are slightly higher compared to the case with the relic density measurement taken as a direct measurement [15].

Bino DM withl ± -coannihilation case-R
We now turn to our fifth scenario, bino DM withl ± -coannihilation Case-R, where in the scan we require the "right-handed" sleptons to be close in mass with the LSP. It should be kept in mind that in our notation we do not mass-order the sleptons: for negligible mixing as it is given for selectrons and smuons the "left-handed" ("right-handed") slepton corresponds tol 1 (l 2 ). As it will be seen below, in this scenario all relevant mass scales are required to be relatively light by the (g − 2) µ constraint 8 .
We start in Fig. 18 with the mχ0 1 -mμ 2 plane with the same color coding as in Fig. 11. By definition of the scenario the points are concentrated on the diagonal. The current (future) (g − 2) µ bound yields upper limits on the LSP of ∼ 700(600) GeV, as well as an upper limit on mμ 2 (which is close in mass to theẽ 2 andτ 2 ) of ∼ 800(700) GeV. These limits remain unchanged by the inclusion of the DM relic density bound. The limits agrees well with the ones found in the case with the relic density measurement taken as a direct measurement [15]. Including the DD and LHC constraints, these limits reduce to ∼ 500 (380) GeV for the LSP for the current (future) (g − 2) µ bounds, and correspondingly to ∼ 550 (440) GeV for mμ 2 . The LHC constraints cut out some, but not all lower-mass points, where the searches for compressed slepton/neutralino spectra [84] are most relevant. Due to the larger splitting the "right-handed" stau turns out to be the NLSP in this scenario, where the uppoer bounds are found at ∼ 500 (380) GeV for the current (anticipated future) (g − 2) µ bounds.
The distribution of the heavier slepton is displayed in the mχ0 1 -mμ 1 plane in Fig. 19. As before, the results agree well with the ones found in Ref. [15]. Although the "left-handed" sleptons are allowed to be much heavier, the (g − 2) µ constraint imposes an upper limit of ∼ 950 (800) GeV in the case for the current (future) (g − 2) µ precision. Taking into account the CDM and LHC constraints we find upper limits for mμ 1 of ∼ 850 GeV and ∼ 800 GeV for the current and anticipated future (g − 2) µ constraint, respectively.
In Fig. 20 we show the results in the mχ0 1 -mχ± 1 plane with the same color coding as in Fig. 18. The results are again in qualitative agreement with the case of the DM relic density taken as a direct measurement [15]. As in the Case-L the (g − 2) µ limits on mχ0 1 become slightly stronger for larger chargino masses. The upper limits on the chargino mass, subsequently joined. However, the ranges covered by the different searches (or colors) are clearly visible and not affected by the different point densities. 8 The scan for case-R turned out to be computationally extremely expensive, resulting in a relatively low density of LHC tested points. Consequently, all upper bounds on particle masses should be taken with a relatively large uncertainty.   Figure 19: The results of our parameter scan in the mχ0 1 -mμ 1 plane for thel-coannihilation case-R scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding is as in Fig. 18. however, are substantially stronger as in the Case-L. Taking all constraints into account, they are reached at ∼ 1540 GeV for the current and ∼ 1350 GeV for the anticipated future precision in a µ . plane for thel-coannihilation case-R scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding is as in Fig. 18.
We finish our analysis of thel ± -coannihilation Case-R with the results in the mχ0 1 -tan β plane, presented in Fig. 21. The overall picture is similar to the previous cases shown above in Figs. 14, 17, and also qualitatively in agreement with the results found in Ref. [15]. Larger LSP masses are allowed for larger tan β values. On the other hand the combination of small mχ0 1 and very large tan β values, tan β > ∼ 40 leads to stau masses below the LSP mass, which we exclude for the CDM constraints. The LHC searches mainly affect parameter points with tan β < ∼ 20. Larger tan β values induce a larger mixing in the third slepton generation, enhancing the probability for charginos to decay via staus and thus evading the LHC constraints. As before we also show the limits from H/A searches at the LHC, where we set (as above) mχ0 1 = M A /2, i.e. roughly to the requirement for A-pole annihilation, where points above the black lines are experimentally excluded. Comparing Case-R and Case-L, here for the current (g − 2) µ limit substantially less points are passing the current (g − 2) µ constraint below the black line, i.e. are potential candidates for A-pole annihilation. The masses reach only up to ∼ 200 GeV. With the anticipated future (g − 2) µ limit hardly any point survives, leaving A-pole annihilation as a quite remote possibility with a strict upper bound on mχ0 1 . -tan β plane in thel-coannihilation case-R scenario for current (left) and anticipated future limits (right) from (g − 2) µ . The color coding is as in Fig. 18. The black line indicates the current exclusion bounds for heavy MSSM Higgs bosons at the LHC (see text).

Lowest and highest mass points
In this section we present some sample spectra for the five cases discussed in the previous subsections. For each case, higgsino DM, wino DM, bino/wino DM withχ ± 1 -coannihilation, l ± -coannihilation Case-L and Case-R, we present three parameter points that are in agreement with all constraints (red points): the lowest LSP mass, the highest LSP with current (g − 2) µ constraints, as well as the highest LSP mass with the anticipated future (g − 2) µ constraint. They will be labeled as "H1, H2, H3", "W1, W2, W3", "C1, C2, C3", "L1, L2, L3", "R1, R2, R3" for higgsino DM, wino DM, bino/wino DM withχ ± 1 -coannihilation, l ± -coannihilation Case-L and Case-R, respectively. While the points are obtained from "random sampling", nevertheless they give an idea of the mass spectra realized in the various scenarios. In particular, the highest mass points give a clear indication on the upper limits of the NLSP mass. It should be noted that the σ SI p values given below are scaled with a factor of (Ωχh 2 /0.118) to account for the lower relic density.
In Tab. 1 we show the 3 parameter points ("H1, H2, H3") from higgsino DM scenario, which are defined by the six scan parameters: M 1 , M 2 , µ, tan β and the two slepton mass parameters, ml L and ml R (corresponding roughly to mẽ 1 ,μ 1 and mẽ 2 ,μ 2 , respectively). Together with the masses and relevant BRs we also show the values of the DM observables and a SUSY µ . By definition of the scenario we have µ < M 1 , M 2 and mχ0 1 ∼ mχ0 2 ∼ mχ± 1 with theχ ± 1 as NLSP. We find in all three points M 2 M 1 , which can be understood from the (g − 2) µ constraint. For (relatively) small µ the chargino-sneutrino contribution, Eq. (5), is ∝ 1/M 2 and thus becomes large for smaller M 2 . Conversely, the neutralino-slepton contri-bution, Eq. (6), is ∝ M 1 and thus becomes larger for larger M 1 . As anticipated, the relic density is found to be quite low. For all of the three points, the decays ofχ ± 1 andχ 0 2 to sleptons are not kinematically accessible. Therefore they mainly decay via off-shell sleptons and/or gauge-bosons to a pair of SM particles and the LSP. The sleptons represented byl 1 andl 2 decay to chargino/neutralino and the corresponding SM particle. Thel 2 is mostly "right-handed" and thus decays preferably to a bino and an electron. However, this mode is kinematically open only for H3, leading to the observedl 2 decay pattern for H1 and H2.   Table 1: The masses (in GeV) and relevant BRs (%) of three points from the higgsino DM scenario, corresponding to the lowest LSP mass, the highest LSP mass with current (g − 2) µ constraints, as well as the highest LSP mass with the anticipated future (g − 2) µ constraint. Herel (l) refers toẽ (e) and µ (µ) together. ν denotes ν e , ν µ and ν τ together. Only BRs above 0.1 % are shown. The values of (g − 2) µ and DM observables are also shown. σ SI p is scaled with a factor of (Ωχh 2 /0.118) and given in the units of pb.
In Tab. 2 we show the 3 parameter points ("W1, W2, W3") from wino DM scenario, defined in terms of the same set of input parameters as the higgsino DM scenario. Together with the masses and relevant BRs we also show the values of the DM observables and a SUSY µ . By definition of the scenario we have M 2 < M 1 , µ and mχ0 1 ∼ mχ± 1 , i.e. theχ ± 1 being the NLSP decaying dominantly asχ ± 1 →χ 0 1 π ± , but also the decay to the LSP and eν e or µν µ occurs (either with a soft pion or a soft charged lepton). As anticipated, the relic density is found to be quite low. The second lightest neutralino is found substantially heavier than the LSP, with a variety of decay channels to be open, diluting possible signals. The sleptons represented byl 1 andl 2 decay to chargino/neutralino and the corresponding SM particle with relevant BRs to charged leptons, which will be crucial for their detection. Thel 2 is mostly "right-handed" and thus decays preferably to a bino and an electron (i.e. the neutralino in the decay is determined by the mass ordering of µ and M 1 ).  In Tab. 3 we show the 3 parameter points ("C1, C2, C3") fromχ ± 1 -coannihilation scenario, with the same definitions and notations as for the previous tables. We find that theτ is the NLSP in these three cases, and the combined contribution fromτ -coannihilation together withχ ± 1 -coannihilation brings the relic density to the ballpark value. For all of the three points, the decays ofχ ± 1 andχ 0 2 to first two generations of sleptons are not kinematically accessible or strongly phase space suppressed. Therefore they decay with a very large BR to third generation charged sleptons and sneutrinos. This makes them effectively invisible to the LHC searches looking for electrons and muons in the signal. LHC analyses designed to specifically look for τ -rich final states can prove beneficial to constrain these points, which are much less powerful, as discussed above. We refrain from showing results forτ decays, since the corresponding dedicated searches turn out to be weaker (and thus not effective) than other applicable searches.  Table 3: The masses (in GeV) and relevant BRs (%) of three points fromχ ± 1 -coannihilation scenario corresponding to the lowest LSP mass, the highest LSP mass with current (g − 2) µ constraints, as well as the highest LSP mass with the anticipated future (g − 2) µ constraint. The notation, definitions and units are as in Tab. 2. Only BRs above 0.1 % are shown.
In Tab. 4 we show three parameter points ("L1, L2, L3") taken froml ± -coannihilation scenario Case-L, defined in the same way as in theχ ± 1 -coannihilation case. The character of theχ 0 2 andχ ± 1 depend on the mass ordering of µ and M 2 and are thus somewhat random. This leads to a large variation in the various BRs of these two particles, possibly diluting any signal involving a certain SM particle, such as a charged lepton. The selectrons or smuons, on the other hand, decay to the LSP and the corresponding SM lepton, which tend to be soft in this case, making way for compressed spectra searches at the LHC.  Table 4: The masses (in GeV) and relevant BRs (%) of three points froml ± -coannihilation scenario Case-L corresponding to the lowest LSP mass, the highest LSP mass with current (g − 2) µ constraints, as well as the highest LSP mass with the anticipated future (g −2) µ constraint. The notation, definitions and units are as in Tab. 1. Only BRs above 0.1 % are shown.
The masses, BRs and values of the (g − 2) µ and DM observables of the parameter point for the Case-R ("R1, R2, R3") are shown in Tab. 5, defined in the same way as in theχ ± 1coannihilation case. As in the case-L the character of theχ 0 2 andχ ± 1 depend on the mass ordering of µ and M 2 and are thus somewhat random. In particular for the two high-mass points R2 and R3 the two mass parameters are relatively close, leading to larger mixings for these two states. This in turn leads to a larger variation in the various BRs of these two particles, possibly diluting any signal involving a certain SM particle, such as a charged lepton. The selectrons or smuons, on the other hand, decay preferably to the LSP and the corresponding SM lepton. As in case-L these tend to be soft, making way for compressed spectra searches at the LHC.  Table 5: The masses (in GeV) and relevant BRs (%) of three points froml ± -coannihilation scenario Case-R corresponding to the lowest LSP mass, the highest LSP mass with current (g − 2) µ constraints, as well as the highest LSP mass with the anticipated future (g −2) µ constraint. The notation, definitions and units are as in Tab. 1. Only BRs above 0.1 % are shown.

Prospects for future colliders
In this section we briefly discuss the prospects of the direct detection of the (relatively light) EW particles at the approved HL-LHC, the hypothetical upgrade to the HE-LHC, the potential future FCC-hh, and at a possible future e + e − collider such as ILC [36,37] or CLIC [37,38]. We concentrate on the compressed spectra searches, relevant for higgsino DM, wino DM and bino/wino DM withχ ± 1 -coannihilation. Results for the future prospects for slepton coannihilation (although with the relic DM density taken as a direct measurement) can be found in Ref. [15]. -∆m plane with anticipated limits from compressed spectra searches at the HL-LHC, the HE-LHC, FCC-hh and ILC500, ILC1000, CLIC380, CLIC1500, CLIC3000 (see text), original taken from Ref. [120]. Not included are disappearing track searches. Shown in blue, dark blue, turquoise are the points surviving all current constraints in the case of higgsino DM, wino DM and bino/wino DM withχ ± 1 -coannihilation, respectively.
In Fig. 22 we present our results in the mχ± 1 -∆m plane (with ∆m := mχ± 1 − mχ0 1 ), which was presented (in its original form, i.e. down to ∆m = 0.7 GeV) in Ref. [120] for the higgsino DM case, but also directly valid for the wino DM case [121] (see the discussion below). Shown are the following anticipated limits for compressed spectra searches 9 : • HL-LHC with 3 ab −1 at √ s = 14 TeV: solid blue and solid red: di-lepton searches (with soft leptons) [124].
[124], but rescaled by 1.5 to take into account the higher energy [125].
Not shown in Fig. 22 are searches for disappearing tracks, which are relevant for very small mass differences, see Figs. 6 and 7. They cut out the lower edge of the dark blue points, wino DM. The life-time limits are expected to improve at the HL-LHC by a factor of ∼ 2 [125], only slightly cutting into the still allowed region.
In the higgsino DM case (blue) it can be seen that the HL-LHC can cover part of the allowed parameter space, but that a full coverage can only be reached at a high-energy e + e − collider with √ s < ∼ 1000 GeV (i.e. ILC1000 or CLIC1500). The wino DM case (dark blue), has larger production cross sections at pp colliders. However, the ∆m is so small in this scenario that it largely escapes the HL-LHC searches (and the different production cross section turns out to be irrelevant). It can be expected that ∆m > ∼ 0.7 GeV can be covered with the monojet searches at the FCC-hh. Very low, but still allowed mass differences can be covered by the disappearing track searches. However, as for the higgsino DM case, also here a high-energy e + e − collider will be necessary to cover the whole allowed parameter space. While the currently allowed points would require CLIC1500, a parameter space reduced by the HL-LHC disappearing track searches, resulting e.g. in mχ0 1 < ∼ 500 GeV could be covered by the ILC1000.
The more complicated case for the future collider analysis is given by the bino/wino parameter points (turquoise), since the limits shown not only assume a small mass difference betweenχ 0 1 andχ ± 1 , but also pp production cross sections as for the higgsino case. For bino/wino DM these production cross sections turn out to be larger (as for the pure wino case), i.e. displaying the wino/bino points in Fig. 22 should be regarded as conservative with respect to the pp based limits. Consequently, it is expected that the HE-LHC or "latest" the FCC-hh would cover this scenario entirely. On the other hand, the e + e − limits should be directly applicable, and large parts of the parameter space will be covered by the ILC1000, and the entire parameter space by CLIC1500.
As discussed in Sect. 4.1 we have not considered explicitly the possibility of Z or h pole annihilation to find agreement of the relic DM density with the other experimental measurements. However, it should be noted that in this context an LSP with M ∼ mχ0 would yield a detectable cross-section e + e − →χ 0 1χ 0 1 γ in any future high-energy e + e − collider. Furthermore, in the case of higgsino or wino DM, this scenario automatically yields other clearly detectable EW-SUSY production cross-sections at future e + e − colliders. For bino/wino DM this would depend on the values of M 2 and/or µ. We leave this possibility for future studies.
On the other hand, the possibility of A-pole annihilation was discussed for all five scenarios. While it appears a rather remote possibility (particularly for higgsino and wino DM), it cannot be fully excluded by our analysis. However, even in the "worst" case of l ± -coannihilation Case-L an upper limit on mχ0 1 of ∼ 300 GeV can be set (about ∼ 40 GeV as as in the case with the relic DM density taken as a direct measurement [15]. While not as low as in the case of Z or h-pole annihilation, this would still offer good prospects for future e + e − colliders. We leave also this possibility for future studies.

Conclusions
The electroweak (EW) sector of the MSSM, consisting of charginos, neutralinos and scalar leptons can account for a variety of experimental data. Concerning the CDM relic abundance, where the MSSM offers a natural candidate, the lightest neutralino,χ 0 1 , while satisfying the (negative) bounds from DD experiments. Concerning the LHC searches, because of comparatively small EW production cross-sections, a relatively light EW sector of the MSSM is also in agreement with the latest experimental exclusion limits. Most importantly, the EW sector of the MSSM can account for the long-standing 3 − 4 σ discrepancy of (g − 2) µ . Improved experimental results are expected soon [3] by the publication of the Run 1 data of the "MUON G-2" experiment.
In this paper we assume that theχ 0 1 provides MSSM DM candidate, where we take the DM relic abundance measurements [5] as an upper limit. We analyzed several SUSY scenarios, scanning the EW sector (chargino, neutralino and slepton masses as well as tan β), taking into account all relevant experimental data: the current limit for (g − 2) µ , the relic density bound (as an upper limit), the DD experimental bounds, as well as the LHC searches for EW SUSY particles. Concerning the latter we included all relevant existing data, mostly relying on re-casting via CheckMATE, where several channels were newly implemented into CheckMATE [15]. Concretely, we analyzed five scenarios, depending on the hierarchy between M 1 , M 2 and µ as well as the mechanism that brings the relic density in agreement with the experimental data. For µ < M 1 , M 2 we have higgsino DM, for M 2 < M 1 , µ wino DM, whereas for M 1 < M 2 , µ one can have mixed bino/wino DM ifχ ± 1 -coannihilation is responsible for the correct relic density (i.e. M 1 < ∼ M 2 ), or one can have bino DM ifl ± -coannihilation yields the correct relic density, with the mass of the "left-handed" ("right-handed") slepton close to mχ0 1 , Case-L (Case-R). Our scans naturally also consist "mixed" cases. Higgsino and wino DM, together with the (g − 2) µ constraint can only be fulfilled if the relic density is substantially smaller than the measured density. The other three scenarios can easily accommodate the relic DM density as direct measurement, which was analyzed in Ref. [15]. In the present analysis these three scenarios were extended to the case where the relic density is used only as an upper bound.
We find in all five cases a clear upper limit on mχ0 1 . These are ∼ 500 GeV for higgsino DM, ∼ 600 GeV for wino DM, while forχ ± 1 -coannihilation we find ∼ 600 GeV, forl ± -coannihilation Case-L ∼ 500 GeV and for Case-R values up to ∼ 500 GeV are allowed. Similarly, upper limits to masses of the coannihilating SUSY particles are found as, mχ0 2 ∼ mχ± 1 < ∼ 510 GeV for higgsino DM, mχ± 1 < ∼ 600 GeV for wino DM, mχ± 1 < ∼ 650 GeV for bino/wino DM withχ ± 1 -coannihilation and ml L < ∼ 600 GeV for bino DM case-L and ml R < ∼ 600 GeV for bino DM case-R. For the latter, in thel ± -coannihilation case-R, the upper limit on the lighterτ is even lower, mτ 2 < ∼ 500 GeV. The current (g − 2) µ constraint also yields limits on the rest of the EW spectrum, although much loser bounds are found. These upper bounds set clear collider targets for the HL-LHC and future e + e − colliders.
In a second step we assumed that the new result of the Run 1 of the "MUON G-2" collaboration at Fermilab yields a precision comparable to the existing experimental result with the same central value. We analyzed the potential impact of the combination of the Run 1 data with the existing result on the allowed MSSM parameter space. We find that the upper limits on the LSP masses are decreased to about ∼ 480 GeV for higgsino DM, ∼ 500 GeV for wino DM, as well as ∼ 500 GeV forχ ± 1 -coannihilation, ∼ 450 GeV forl ± -coannihilation Case-L and ∼ 380 GeV in Case-R, sharpening the collider targets substantially. Similarly, the upper limits on the NLSP masses go down to about mχ0 2 ∼ mχ± 1 < ∼ 490 GeV for higgsino DM, mχ± 1 < ∼ 600 GeV for wino DM, mχ± 1 ∼ 550 GeV for bino/wino DM withχ ± 1 -coannihilation, as well as ml L < ∼ 550 GeV for bino DM case-L and ml R < ∼ 440 GeV, mτ 2 < ∼ 380 GeV for bino DM case-R.
For the three cases with a small mass difference between the lightest chargino and the LSP (higgsino DM, wino DM and bino/wino DM withχ ± 1 -coannihilation) we have also briefly analyzed the prospects for future collider searches, specially targeting these small mass differences. The results for the points surviving all (current) constraints have been displayed in the plane of the Next-to-LSP (in these cases mχ± 1 ) vs. ∆m = mχ± 1 − mχ0 1 , overlaid with the anticipated limits from future collider searches from HL-LHC, HE-LHC, FCC-hh, ILC and CLIC [120,121]. (These limits are directly applicable to higgsino and wino DM, and likely to be overly conservative for bino/wino DM withχ ± 1 -coannihilation.) While parts of the parameter spaces can be covered by future pp machines, a full exploration of the parameter space requires a future e + e − colliders. Taking into account limits from future e + e − machines, a center-of-mass energy of √ s = 1 TeV will be sufficient to conclusively explore higgsino DM, wino DM and bino/wino DM withχ ± 1 -coannihilation. While we have attempted to cover nearly the full set of possibilities that the EW spectrum of the MSSM presents, while being in agreement with all the various experimental constraints, our studies can be extended/completed in the following ways. One can analyze the cases of: (i) complex parameters in the chargino/neutralino sector (then also taking EDM constraints into account); (ii) different soft SUSY-breaking parameters in the three generations of sleptons, and/or between the left-and right-handed entries in the case of χ ± 1 -coannihilation; (iii) A-pole annihilation, in particular in the case ofl ± -coannihilation for very low mχ0 1 and tan β values; (iv) h-and Z-pole annihilation, which could be realized for sufficiently heavy sleptons. On the other hand, one can also restrict oneself further by assuming some GUT relations between, in particular, M 1 and M 2 . We leave these analyses for future work.
In this paper we have analyzed in particular the impact of (g − 2) µ measurements on the EW SUSY spectrum. The current measurement sets clear upper limits on many EW SUSY particle masses. On the other hand we have also clearly demonstrated the potential of the upcoming measurements of the "MUON G-2" collaboration, which have a strong potential of sharpening the future collider experiment prospects. We are eagerly awaiting the new "MUON G-2" result to illuminate further the possibility of relatively light EW BSM particles.