Next Article in Journal
On Possible Origin of an Artificial Wormhole
Previous Article in Journal
Spinor-Vector Duality and the Swampland
Previous Article in Special Issue
Progress in the Simulation and Modelling of Coherent Radio Pulses from Ultra High-Energy Cosmic Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Dark Matter in Supersymmetry

1
Instituto de Física Teórica (IFT) UAM-CSIC, Campus de Cantoblanco, 28049 Madrid, Spain
2
Departamento de Física Teórica, Campus de Cantoblanco, Universidad Autónoma de Madrid (UAM), 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Universe 2022, 8(8), 427; https://doi.org/10.3390/universe8080427
Submission received: 31 March 2022 / Revised: 27 July 2022 / Accepted: 29 July 2022 / Published: 18 August 2022
(This article belongs to the Special Issue Studying the Universe from Spain)

Abstract

:
Supersymmetry is a well-motivated theory for physics beyond the Standard Model. In particular, supersymmetric models can naturally possess dark matter candidates that can give rise to the measured dark matter content of the universe. We review several models that have been analyzed with regard to dark matter by groups based in Spain in recent years. These models include, in particular, the Minimal Supersymmetric Standard Model (MSSM) and the ‘ μ from ν ’ Supersymmetric Standard Model ( μ ν SSM) in various versions.

1. Introduction

Searches for Dark Matter (DM) is one of the main objectives in today’s particle and astroparticle physics. Searches at the Large Hadron Collider (LHC) (or other collider experiments) are complementary to the searches in “direct detection” (DD) experiments. Among the Beyond the Standard Model (BSM) theories that predict a viable DM particle, the Minimal Supersymmetric Standard Model (MSSM) [1,2,3,4] is one of the leading candidates. Supersymmetry (SUSY) predicts two scalar partners for all Standard Model (SM) fermions as well as fermionic partners to all SM bosons. Furthermore, contrary to the SM case, the MSSM requires two Higgs doublets.
This results in five physical Higgs bosons instead of the single Higgs boson in the SM: the light and heavy charge conjugation parity symmetry ( CP )-even Higgs bosons, h and H, the CP -odd Higgs boson, A and the charged Higgs bosons, H ± . The neutral SUSY partners of the (neutral) Higgs and electroweak (EW) gauge bosons give rise to the four neutralinos, χ ˜ 1 , 2 , 3 , 4 0 .
The corresponding charged SUSY partners are the charginos, χ ˜ 1 , 2 ± . The SUSY partners of the SM leptons and quarks are the scalar leptons and quarks (sleptons and squarks), respectively. The lightest SUSY particle (LSP) is naturally the lightest neutralino, χ ˜ 1 0 . It can make up the full DM content of the universe [5,6], or—depending on its nature—only a fraction of it. In the latter case, e.g., a SUSY axion [7] could be an additional DM component, bringing the total DM density into agreement with the experimental measurement.
Another particularly well motivated SUSY model is the ‘ μ from ν ’ Supersymmetric Standard Model ( μ ν SSM ) [8,9], (see also [10] for a recent review of the μ ν SSM and [11] for an vacuum structure analysis of the μ ν SSM . Other information about the model can be found in http://dark.ft.uam.es/mununiverse (accessed on 28 July 2022)). Beyond the well-known appealing features of commonly studied SUSY models, in the μ ν SSM, the tiny neutrino masses and their mixings can be accommodated via an EW seesaw mechanism, where it is required that the matter content is enlarged with regard to the SM by right-handed neutrinos [8,12,13,14,15]. Their superpartners, the “right-handed” (This name is used also for the scalar particles in order to indicate that they are the superpartners of right-handed fermions) scalar neutrinos (sneutrinos) are gauge singlet scalar fields.
If the right-handed sneutrinos acquire vacuum expectation values (vevs), the so-called μ -term of the MSSM can effectively be generated, analogous to the Z 3 symmetric Next-to-MSSM (NMSSM) [16,17]. Consequently, the μ ν SSM also offers a solution to the so-called μ -problem [18]. By construction, the μ ν SSM does not permit a consistent assignment of conserved R-parity charges. It is thus a R-parity violating (RPV) model and therefore has no stable LSP.
Therefore, compared to the (N)MSSM, the collider constraints from the LHC are substantially weaker [19,20,21,22,23,24,25,26]. Regarding the DM content of the universe, it is possible to accommodate the measured relic abundance by means of the decaying—but long-lived—gravitino [27,28,29,30,31] (For an analysis of gravitino DM in the context of the bilinear RPV model, see Ref. [32]) or axino [33], interestingly producing gamma-rays potentially detectable in gamma-ray telescopes. Multicomponent DM scenarios with the axino/gravitino as the LSP and the gravitino/axino as the next-to-LSP (NLSP) were also discussed in the μ ν SSM [31,33]. Concerning cosmology, baryon asymmetry might be realized in the μ ν SSM through EW baryogenesis [34].
In this review, we give an overview of DM analyses in various SUSY models that have been obtained in the recent years by groups (to a relevant degree) based in Spain.(In the past, other DM candidates were proposed, such as the right-handed sneutrino in the NMSSM extended with right-handed neutrino superfields [35,36]. For a recent analysis, see Ref. [37]) The models comprise the pMSSM11 (the phenomenological MSSM with 11 free parameters); the electroweak (EW) MSSM, where the colored sector is assumed to be heavy and does not enter into the low-energy phenomenology; the μ ν SSM with gravitino and/or axino as DM candidates (also the possibility of a sterile right-handed neutrino will be presented); and finally the recently proposed U μ ν SSM  [38,39].
The latter is a U ( 1 ) extension of the μ ν SSM , where the presence of weakly interacting massive particles (WIMPs) as DM candidates dictated by the anomaly cancellation conditions was proven [40]. This a remarkable result, given that the U μ ν SSM is a RPV scenario.

2. Results in the pMSSM11

In this section, we review the results of the DM analysis in the pMSSM11 as obtained in Ref. [41]. The main idea is to investigate phenomenological models that have soft SUSY-breaking parameters that are not constrained by any universality condition (at the Grand Unified Theory (GUT) scale), though subject to milder constraints emanating, in particular, from upper limits on SUSY contributions to flavor-changing processes. These phenomenological MSSM (pMSSM) [42] models that have been studied in the literature contain up to 19 free parameters.

2.1. The Parameter Space

Here, we employ a model with 11 free parameters, given at the EW scale. These are three independent gaugino masses, M 1 , 2 , 3 ; a common scalar mass parameter for the first- and second-generation squarks, m q ˜ ; a distinct mass parameter for the third-generation squarks, m q ˜ 3 ; a common mass parameter m ˜ for the sleptons of the first- and second-generation; a distinct mass parameter for the scalar taus, m τ ˜ ; a single trilinear mixing parameter A, which is taken to be universal at the electroweak scale; the Higgsino mass parameter μ ; the mass of the pseudoscalar Higgs, M A and tan β = v 2 / v 1 ; and the ratio of the Higgs vevs.
The renormalization scale M SUSY given by the geometric mean of the masses of the two scalar tops, M SUSY m t ˜ 1 m t ˜ 2 , which is also the scale at which the EW symmetry breaking conditions are imposed. The ranges of the sampled parameter space are summarized in Table 1.

2.2. The Analysis Framework

A global likelihood analysis of the pMSSM11 was performed, including constraints from direct searches for SUSY particles at the LHC, measurements of the Higgs boson mass and signal strengths, LHC searches for SUSY Higgs bosons, precision electroweak observables, flavor constraints from B- and K-physics observables, the cosmological constraint on the overall cold DM (CDM) density, and upper limits on spin-independent and -dependent LSP-nuclear scattering. Furthermore, ( g 2 ) μ is included as an additional constraint.
The calculation of the observables contributing to the likelihood is performed with the MasterCode tool [41,43,44,45,46,47,48,49,50,51,52]. This tool combines consistently and interfaces various private and public codes employing the SUSY Les Houches Accord (SLHA) [53]. The analysis within MasterCode uses the following codes: SoftSusy 3.3.9 [54] for the spectrum, FeynWZ [55,56,57] for the electroweak precision observables,FeynHiggs 2.11.3 [58,59,60,61,62,63,64,65,66,67] for the Higgs sector and ( g 2 ) μ , SuFla [68,69] and SuperIso [70,71,72] for the flavor physics observables, Micromegas-3.2 [73] for the DM relic density, SSARD [74] for the spin-independent and -dependent elastic scattering cross-sections σ p SI and σ p SD .
The uncertainties in the cross-sections are derived from a straightforward propagation of errors in the input quantities that determine the cross-section. The dominant uncertainties are discussed below in more detail, SDECAY 1.3b [75] for calculating sparticle branching ratios and HiggsSignals 1.4.0 [76,77,78] and HiggsBounds 4.3.1 [79,80,81,82,83] for calculating constraints on the SUSY Higgs sector. The experimental values used for the analysis are given in [41,51] (and references there in).

2.3. DM Results in the pMSSM11

The first set of results for the χ ˜ 1 0 , which is assumed to yield the full amount of CDM, is presented in Figure 1. The upper plot shows the profile likelihood functions in one dimension for m χ ˜ 1 0 in the pMSSM11 with (without) the ( g 2 ) μ constraint in blue (green), as well as dashed (solid) for (not) applying the constraints from LHC Run II. Including all constraints, one can see that a clear preference is found for a relatively low value of m χ ˜ 1 0 . At the 2 σ level, it is restricted to be m χ ˜ 1 0 < 500 GeV .
Comparing the blue and the green lines, it becomes apparent that this result strongly relies on the ( g 2 ) μ bound (see Section 3 for a detailed discussion). The lower plot of Figure 1 shows the triangular presentations of the composition of the χ ˜ 1 0 in the fit with LHC 13-TeV and with the ( g 2 ) μ constraint. The color coding indicates the Δ χ 2 . The best-fit point is marked with a green star. One can observe that a small Wino fraction N 12 2 < 0.1 is strongly favored, while the relative proportions of the Bino fraction N 11 2 and the Higgsino fraction N 13 2 + N 14 2 are relatively unconstrained at the 95% confidence level (CL). The best-fit point is found as a pure bino LSP.
In Figure 2, we show the preferred parameter regions in the m χ ˜ 1 0 σ p SI plane (upper plot) and the m χ ˜ 1 0 σ p SD plane (lower plot). σ p SI and σ p SD denote the spin-independent and spin-dependent DM-nucleon cross sections, respectively. In the upper plot, the upper limits established by the PandaX-II [84,85], XENON1T [86] and LUX [87] Collaborations are shown as blue, magenta and green contours, respectively. The combined limit (with green shading above) is indicated by a black line.
The future projected 90% CL exclusion sensitivities of the XENON1T/nT [88] and LUX-Zeplin (LZ) [89] experiments are shown as dashed blue and solid purple lines, respectively. The background neutrino ‘floor’ (with a yellow shading below) is shown as a dashed orange line. Although the DAMA/LIBRA collaboration observed an annual modulation in the detection rate [90], it is difficult to reconcile it with the negative results from the other experiments. It is true that most of them use different detection techniques than DAMA/LIBRA; however, ANAIS-112 using the same technique and target material, has not reported of hints of the presence of modulation after three-year exposure analysis [91].
The solid lines in red, blue and green indicate the 1 , 2 , 3 σ preferred parameter regions of the global fit, respectively, with the best-fit point shown as a green star. The color coding within the 2 σ area correspond to the DM relic density mechanism. We see that m χ ˜ 1 0 100 GeV , with upper limit m χ ˜ 1 0 < 550 at the 95% CL, in agreement with Figure 1. One furthermore sees that, at the best-fit point, the nominal prediction for σ p SI is at the level of the sensitivities projected for the planned XENON1T/nT and LUX-Zeplin (LZ) experiments (solid purple line). The ranges of the nominal values of σ p SI at the 68% and 95% CL slightly extend below the neutrino floor.
Large values of σ p SI are found in the chargino coannihilation region (indicated by green shading), with other DM mechanisms, including squark coannihilation, yielding large values of σ p SI for m χ ˜ 1 0 > 1 TeV . This and the other DM mechanisms indicated, however, also allow much smaller values of σ p SI .
In the lower plot of Figure 2, the spin-dependent cross sections are analyzed. Here, the upper limit as reported by the PICO Collaboration [92] is shown as a purple contour (with a green shading). For this case, the neutrino floor for σ p SD is taken over from [93]. Furthermore, we show the indicative upper limits from IceCube [94] and SuperKamiokande [95] searches for energetic solar neutrinos obtained, which assume a predominant annihilation of the LSPs into τ + τ . However, these are subject to larger theoretical uncertainties and will not be further discussed (see, however, Ref. [41]).
The preferred regions of the parameter space are indicated as in the upper figure. Similar to the case of σ p SI , one can observe that the permitted values of m χ ˜ 1 0 range from ∼100 GeV to ∼550 GeV. The uncertainties in the σ p SD calculation are substantially smaller than those found for σ p SI . One can see that the regions of the 68 and 95% CL in the nominal σ p SD calculations are found below the upper limit of PICO [92] (as indicated by the solid purple line). Concerning the best-fit point the nominal predictions is found ∼3 orders of magnitude below the limit given currently by PICO.
Overall, the global analysis of the pMSSM11 shows that DD experiments, in particular via the spin-independent searches, can cover a relevant part of the parameter space. However, even the 1 σ preferred regions reach below the neutrino floor and may thus escape the discussed DD experiments. In this case either novel techniques for the DM DD experiements are needed, or collider experiments will be necessary to cover the full parameter space (see also the discussion in Section 3).

3. Results in the EW-MSSM

In this section, we review recent DM analyses in the EW-MSSM [96]. In this analysis all relevant constraints on the EW sector of the MSSM are considered, which, in particular, includes the anomalous magnetic moment of the muon, ( g 2 ) μ .

3.1. The EW-MSSM

The EW-MSSM is characterized by a light EW SUSY sector, where all other particles are assumed to be heavy. The masses and mixings of the neutralinos are defined by (on top of SM parameters) the S U ( 2 ) L and U ( 1 ) Y gaugino masses, M 2 and M 1 , the Higgsino parameter μ , as well as tan β : = v 2 / v 1 , the ratio of the two vevs of the two Higgs doublets. After the diagonalizing the mass matrix the four eigenvalues yield the four neutralino masses m χ ˜ 1 0 < m χ ˜ 2 0 < m χ ˜ 3 0 < m χ ˜ 4 0 . Similarly, the chargino masses and mixings are given (on top of SM parameters) by μ , M 2 and tan β . The diagonalizing the mass matrix gives the two chargino-mass eigenvalues m χ ˜ 1 ± < m χ ˜ 2 ± .
For the sleptons a common soft SUSY-breaking parameters for all three generations was chosen. The mass matrices of the charged sleptons are given (on top of SM parameters) by the soft SUSY-breaking diagonal parameters m l ˜ L 2 and m l ˜ R 2 , as well as A l ( l = e , μ , τ ), the trilinear Slepton-Higgs coupling A l , where, however, the later are set to zero. The mixing between the “right-handed” and “left-handed” scalar leptons is only relevant for scalar taus, since the off-diagonal entry of the mass matrix is dominated by m τ tan β μ .
Therefore, in the first two generations, the masses can be approximated as m l ˜ 1 m l ˜ L , m l ˜ 2 m l ˜ R (where small D-terms are assumed). In general, we follow the convention that l ˜ 2 ( l ˜ 1 ) has the large “right-handed” (“left-handed”) component, i.e., no mass ordering is imposed. In addition m l ˜ 1 and m l ˜ 2 , the symbols are equal for all three generations, we also use m e ˜ 1 , 2 , m μ ˜ 1 , 2 and m τ ˜ 1 , 2 for the masses of the scalar electrons, muons and taus. The charged slepton and the sneutrino masses are connected by the usual S U ( 2 ) relation.
Overall, at the tree level, the EW sector is described in our analysis by six parameters: M 2 , M 1 , μ , m l ˜ L , m l ˜ R and tan β . Here, we furthermore assume M 1 , M 2 , μ to be positive. It was shown in [97] that these positive parameter choices cover the relevant parameter space when the ( g 2 ) μ limits are considered (see the discussion below).
As mentioned above, for strongly interacting particles, we make the assumption that the sector of colored particles in the MSSM is substantially heavier than the EW sector and, therefore, does not play a relevant role in our analysis.Concerning the Higgs-boson sector, we make the assumption that the higher-order corrections to M h , the light CP -even Higgs boson mass, largely originating from the stop/top sector, give a value of M h ∼125 GeV , in agreement with the experimental data. This implies scalar top masses that are naturally in the TeV range [41,98], which are thus in agreement with the LHC bounds. Concerning the heavy Higgs-boson mass scale, M A , it was shown in [97,99,100] that A-pole annihilation is largely excluded. For simplicity, we therefore assume M A to be sufficiently large so as to not play any relevant role in our analysis.

3.2. The Relevant Constraints

The SM prediction of a μ is given by [101]. The comparison with the combined experimental new world average, based on Refs. [102,103] yields a deviation of Δ a μ = ( 25.1 ± 5.9 ) × 10 10 , corresponding to 4.2 σ . This result is used as a hard 2 σ cut on our results. The prediction of ( g 2 ) μ in the MSSM is calculated using GM2Calc [104], implementing two-loop corrections from [105,106,107] (see also [108,109]). Vacuum stability constraints are considered with the public code Evade [110,111]. All relevant SUSY searches for EW particles are considered, mostly via CheckMATE [112,113,114] (see Ref. [97] for details on many analyses newly implemented by our group).
For the DM relic density constraints we use the latest result from Planck [115], either as a direct measurement, or as an upper bound. The relic density in the MSSM is evaluated with MicrOMEGAs [73]. For the DD DM constraints, we use the results for the spin-independent DM scattering cross-section σ p SI from XENON-1T [86] experiment. The theoretical predictions are evaluated using the public code MicrOMEGAs.

3.3. Five Viable DM Scenarios

Five different scenarios were analyzed, classified by the mechanism that brings the LSP relic density into agreement with the measured values. The scenarios differ by the NLSP, or equivalently by the mass hierarchies between the mass scales determining the neutralino, chargino and slepton masses. These mass scales are the gaugino soft-SUSY breaking parameters M 1 and M 2 , the Higgs mixing parameter μ and the slepton soft SUSY-breaking parameters m l ˜ L and m l ˜ R , see [96,97,99,100] for a detailed description. The five scenarios can be summarized as follows,
(i)
higgsino DM ( μ < M 1 , M 2 , m l ˜ L , m l ˜ R ), DM relic density is only an upper bound (the rull relic density implies m χ ˜ 1 0 1 TeV and ( g 2 ) μ cannot be fulfilled), m ( N ) LSP < 500 GeV with m NLSP m LSP 5 GeV ;
(ii)
wino DM ( M 2 < M 1 , μ , m l ˜ L , m l ˜ R ), DM relic density is only an upper bound, (the rull relic density implies m χ ˜ 1 0 3 TeV and ( g 2 ) μ cannot be fulfilled), m ( N ) LSP < 600 GeV with m NLSP m LSP 0.3 GeV ;
(iii)
bino/wino DM with χ ˜ 1 ± -coannihilation ( M 1 < M 2 ), DM relic density can be fulfilled, m ( N ) LSP < 650 ( 700 ) GeV ;
(iv)
bino DM with l ˜ ± -coannihilation case-L ( M 1 < m l ˜ L ), DM relic density can be fulfilled, m ( N ) LSP < 650 ( 700 ) GeV ; and
(v)
bino DM with l ˜ ± -coannihilation case-R ( M 1 < m l ˜ R ), DM relic density can be fulfilled, m ( N ) LSP < 650 ( 700 ) GeV .

3.4. DM Results in the EW-MSSM

In this section, we review our results for the DM DD prospects in the five scenarios [96]. We take into account the projections for the exclusion reach of XENON-nT [116] and of the LZ experiment [117] (which effectively agree with each other). We also include the projections of the DarkSide [118] and Argo [119] experiments, which can go down to even lower cross sections, as well as the neutrino floor (NF) [120].
The results are summarized in Figure 3, where we show the m χ ˜ 1 0 σ p SI planes for higgsino DM (upper left), wino DM (upper right), bino DM case-L (middle left) and case-R (middle right) and bino/wino DM with χ ˜ 1 ± -coannihilation (lower plot). The color code indicates the DM relic density, where the red points are in full agreement with the Planck measurement. The black dashed, blue dashed, blue dot-dashed and black dot-dashed lines indicate the prospects for LZ/Nenon-nT, DarkSide, Argo and the NF, respectively.
One can observe that in for higgsino and wino DM all points will be covered by the next round of DD experiments, LZ and/or Xenon-nT are sufficient to cover the whole parameter space. The situation is different for bino DM case-L/R and bino/wino DM. In these three cases cases a large part of the allowed points cannot be probed by LZ/Xenon-nT. The Argon-based experiments can cover a substantially larger part of the allowed parameter space, where he parameter points giving the full DM relic density are mostly covered by LZ/Xenon-nT; however, DarkSide/Argo might be needed in the case of χ ˜ 1 ± -coannihilation, ascan be seen in the lower plot.
However, in all three scenarios, some parameter points are allowed even below the NF, which makes them unaccessible to current DD techniques (see [96] for a short discussion on future directional detection techniques). The allowed parameter spaces below the NF are relatively restricted in the LSP mass, which is bound to be m χ ˜ 1 0 < 400 GeV .
At the High Luminosity LHC (HL-LHC), these points may still remain elusive due to the small mass splitting between the NLSP and the LSP (see the discussion in Ref. [96]). On the other hand, such small mass splittings are do not pose a problem at e + e colliders. The direct production of EW particles at e + e colliders requires a sufficiently high center-of-mass energy, s . Consequently, we focus here on a proposals for linear e + e colliders, the ILC [121,122], which can reach energies up to 1 TeV , which we denote as ILC1000. We evaluate the cross-sections for the various LSP and NLSP pair production modes for s = 1 TeV .
At the ILC1000 an integrated luminosity of 8 a b 1 is foreseen [123,124]. The cross-section predictions are based on tree-level results, obtained as in Refs. [125,126]. Here, we do not attempt a rigorous experimental analysis, but follow analyses [127,128,129] that indicate that to a good approximation final states with the sum of the masses smaller than the center-of-mass energy can be detected.
In Figure 4, we show the LSP and NLSP pair production cross sections for an e + e collider at s = 1000 GeV as a function of the two (identical) final state masses. The upper plot shows σ ( e + e χ ˜ 1 0 χ ˜ 1 0 ( + γ ) ) production (Our tree level calculation does not include the photon radiation, which appears only starting from the one-loop level. However, such an ISR photon is crucial to detect this process due to the invisible final state. We take our tree-level cross section as a rough approximation of the cross section, including the ISR photon (see also [125]), and use the notation “ ( + γ ) ”.) in green and σ ( e + e χ ˜ 1 + χ ˜ 1 ) in blue. The open circles are the points below the anticipated XENON-nT/LZ limit, whereas the solid circles are the points below the neutrino floor. All points are within the reach of the ILC1000.
The cross sections range roughly from ∼100 f b for low masses to ∼100 f b for larger masses, with only a very few points have smaller cross sections. Overall, assuming an integrated luminosity of 8 ab 1 , this corresponds to ∼80,000–800,000 events. Consequently, in contrast to the HL-LHC, the e + e colliders show a clear and conclusive complementarity to the future DD experiments. The χ ˜ 1 ± -coannihilation scenario will be fully covered by either DD experiments or by searches at the ILC1000.
The lower plots of Figure 4 show the LSP and NLSP production cross section in the l ˜ ± -coannihilation scenario for case-L (left) and case-R (right). The green points show again σ ( e + e χ ˜ 1 0 χ ˜ 1 0 ( + γ ) ) , whereas the violet points left and right show σ ( e + e μ ˜ 1 μ ˜ 1 ) (case-L) and σ ( e + e μ ˜ 2 μ ˜ 2 ) (case-R), respectively. Open and full circles denote, as above, the points below the anticipated XENON-nT/LZ limit and the neutrino floor.
The visible spread in the χ ˜ 1 0 χ ˜ 1 0 ( + γ ) production for case-L with regard to the case-R is a result of the more complex structure of the e ± - e ˜ L - χ ˜ 1 0 coupling as compared to the e ± - e ˜ R - χ ˜ 1 0 coupling, dominating the t-channel exchange diagram, respectively. In both cases, we see that, as for the χ ˜ 1 ± -coannihilation case, all points result in particles that can be pair produced at the ILC1000 (except the very highest mass points in case-L). The cross sections range between 100 to 10 fb for χ ˜ 1 0 χ ˜ 1 0 ( + γ ) , and between 20 to 1 fb for smuon pair production. Even for the smallest production cross section, this corresponds to ∼8000 events in the foreseen 1000 GeV ILC run. Furthermore, these two cases can conclusively be probed in the conjunction of DD experiments and an e + e collider at s = 1000 GeV , in contrast to the HL-LHC, where the prospects are less clear (see the discussion in [96]).

4. Results in the μ ν SSM

In the μ ν SSM [8], the presence in the superpotential of the RPV couplings λ i ν ^ i c H ^ d H ^ u and κ i j k ν ^ i c ν ^ j c ν ^ k c . where ν ^ i c denote the right-handed neutrino superfields, solves dynamically both the μ - and ν -problems. After the successful EW symmetry breaking, the right-handed sneutrinos ν ˜ i R develop vevs of the order of TeV, and therefore the λ couplings generate an effective mass term for Higgsinos, μ = λ i ν ˜ i R , solving the μ -problem. In addition, the κ couplings, generate effective Majorana masses for right-handed neutrinos, M i j = 2 κ i j k ν ˜ k R , solving the ν -problem through an EW-scale seesaw, i.e., it is possible to accommodate the correct neutrino masses and mixing angles with Y i j ν < 10 6 [8,12,13,14,15].
As of these small neutrino Yukawa couplings, the RPV is small in the μ ν SSM . Note that, in the limit Y i j ν 0   ν ^ i c can be identified as pure singlet superfields without lepton number, and, as a consequence, R-parity is conserved similarly to the NMSSM. Thus, Y i j ν are the parameters determining RPV in the μ ν SSM . Since the LSP is not stable, it decays into SM particles with a low decay width related to the smallness of neutrino masses.
This connection between neutrino masses and RPV is relevant to have the gravitino or the axino as DM candidates, with large lifetimes due to the fact that their interactions are further suppressed by the small RPV parameters. In what follows, we briefly analyze the gravitino and/or axino LSP as DM in the μ ν SSM , as well as discuss another DM candidate from the neutrino sector.

4.1. Gravitino

If R-parity is conserved the gravitino LSP is stable and therefore a DM candidate in the framework of supergravity [130,131,132,133] (see also Ref. [134] and references therein). Nevertheless, in the case of RPV the gravitino LSP can also be a (decaying) DM candidate [135,136] and can be detected through the observation of gamma-ray lines (and a smooth spectral signature).
Alternatively, superheavy gravitino DM from Starobinsky supergravity was recently proposed [137], which could also give rise to RPV decays. Decaying gravitino DM was studied in the μ ν SSM for the case of the Fermi Large Area Telescope (Fermi-LAT) [138] in Refs. [27,28,29,30,31]. The prospects for detecting μ ν SSM gravitino DM in future gamma-ray missions, such as enhanced ASTROGAM (e-ASTROGAM) [139] and All-sky Medium Gamma-ray Observatory (AMEGO) [140], were analyzed in Ref. [31].
The gravitino interacts with photon and photino; however, photino and left-handed neutrinos are mixed, and—as a consequence—the gravitino LSP decays into photon and neutrino as shown in Figure 5.
The signals are gamma-ray lines with energies half of the gravitino mass m 3 / 2 . The gravitino decay width is given by [135,136]:
Γ ( Ψ 3 / 2 γ ν i ) m 3 / 2 3 32 π M P 2 | U γ ˜ ν | 2 ,
where Γ ( Ψ 3 / 2 γ ν i ) denotes a sum of the partial decay widths into ν i and ν ¯ i , M P 2.43 × 10 18 GeV is the reduced Planck mass, and the mixing parameter U γ ˜ ν determines the photino content of the neutrino,
U γ ˜ ν 2 = i = 1 3 N i 1 cos θ W + N i 2 sin θ W 2 .
Here, N i 1 ( N i 2 ) is the bino (wino) component of the i-th neutrino, and θ W is the weak mixing angle. One can easily estimate the value of | U γ ˜ ν | in the μ ν SSM [27], with the result
U γ ˜ ν g v i L M 1 ,
where v i L are the vevs of the left-handed sneutrinos. For typical electroweak-scale values for the bino mass M 1 and v i L < 10 4 GeV in the μ ν SSM since their values are determined by the small neutrino Yukawas, one obtains
10 8 | U γ ˜ ν | 10 6 .
This was confirmed in Refs. [27,30] by performing scans in the low-energy parameters of the μ ν SSM in order to reproduce the observed neutrino masses and mixing angles. Relaxing some of the assumptions, such as an approximate GUT relation for gaugino masses and/or TeV scales for them, the lower bound can even be smaller: 10 10 | U γ ˜ ν | 10 6 . As we can see from Equation (1), the gravitino decay is suppressed both, by the small RPV mixing parameter | U γ ˜ ν | and by the scale of the gravitational interaction, making its lifetime much longer than the age of the Universe τ 3 / 2 t t o d a y 10 17 s, with
τ 3 / 2 = Γ 1 ( ψ 3 / 2 γ ν i ) 3.8 × 10 33 s 10 8 | U γ ˜ ν | 2 0.1 GeV m 3 / 2 3 .
There are also three-body decays producing a smooth spectral signature. These processes were included in the analysis of Ref. [30]. The results imply that gravitino mass must be smaller than about 17 GeV and its lifetime larger than 4 × 10 25 s.

4.2. Axino

The Axino LSP is another decaying DM candidate in the μ ν SSM . Similar to the gravitino, it has an interaction term with photon and photino. The interaction is suppressed by the small RPV parameters as in the gravitino case and by the Peccei-Quinn scale instead of the gravitational one. As a consequence, the axino can have a lifetime greater than the age of the Universe and produces a signal with a gamma-ray line with energy half of the axino mass.
This framework was studied in Ref. [33] for Kim–Shifman–Vainshtein–Zakharov (KSVZ) [141,142] and Dine–Fischler–Srednicki–Zhitnitsky (DFSZ) [143,144] axion scenarios. The result implies that Fermi-LAT constraints impose that the axino mass must be smaller than about 3 GeV. Furthermore, a significant region of the parameter space lies in the ballpark of the proposed e-ASTROGAM. This would allow the exploration of masses in the range 2 MeV–3 GeV, as well as lifetimes of about 2 × 10 26 8 × 10 30 s.

4.3. Right-Handed Neutrino

Since the number of right-handed neutrinos is a free parameter in the μ ν SSM , some of them might behave as sterile neutrinos and be viable candidates for warm DM [145]. This is similar to the (non-SUSY) neutrino minimal standard model ( ν MSM) [146]. In the case of the μ ν SSM , we need some of the right-handed neutrinos to have small couplings, Y ν 10 13 and κ 10 8 , in order to obtain obtain keV masses and lifetimes long enough to be DM candidates [147].

4.4. Multicomponent DM

The framework of a DM scenario made of gravitinos and axinos was studied in Refs. [31,33]. If the axino is the LSP, a gravitino NLSP can live enough as to contribute to the relic abundance. Then, both particles can produce a signal detectable by future MeV-GeV gamma-ray telescopes, such as e-ASTROGAM. In addition, a well-tempered mixture of both particles can be found in a particular region producing a double-line signal as a smoking gun. Similar results are obtained with a gravitino LSP and an axino NLSP. Other candidates could also contribute to the total amount of relic abundance, such as the sterile neutrino, the axion, etc.
Summarizing, the μ ν SSM is an interesting framework where the problem of DM can be solved through different candidates. The potential tensions between the standard Λ CDM model and cosmological observations (see, e.g., Ref. [148]) might be relaxed if several of these candidates contribute to DM.

5. Results in the U μ ν SSM

We review, in this section, a specific WIMP DM realization in the framework of the U μ ν SSM [40]. The U μ ν SSM [38] is a U ( 1 ) extension of the μ ν SSM, where baryon-number-violating operators as well as explicit mass terms are forbidden, and the potential domain wall problem is avoided. In order to ensure an anomaly free theory, states charged under the new gauge symmetry are introduced: exotic quarks K i and additional singlets ξ ^ α under the SM gauge group, such that we consider the following relevant superpotential [38]:
W = Y i j e H ^ d L ^ i e ^ j c + Y i j d H ^ d Q ^ i d ^ j c Y i j u H u ^ Q ^ i u ^ j c Y i j ν H u ^ L ^ i ν ^ j c + λ i ν ^ i c H u ^ H ^ d + k i α β ν ^ i c ξ ^ α ξ ^ β + Y i j K ν ^ i c K ^ j K ^ j c
where the summation convention is implied on repeated indexes, with i , j , k = 1 , 2 , 3 , the usual family indexes of the SM and α , β = 1 , 2 in the simplest construction. Our convention for the contraction of two S U ( 2 ) doublets is, e.g., H ^ u H ^ d ϵ a b H ^ u a H ^ d b , a , b = 1 , 2 and ϵ a b the totally antisymmetric tensor with ϵ 12 = 1 .
Masses for these new states are generated dynamically once the right-handed sneutrino acquires a vev, simultaneously generating the μ -term and masses for right-handed neutrinos.

5.1. The U μ ν SSM and WIMP DM

A discrete Z 2 symmetry is present in the superpotential term of Equation (6) containing the superfields of type ξ ^ α , under which they have a charge 1 and the rest of the particle content a charge +1. This symmetry is not an extra requirement but arises from the charge assignment of the model and is a consequence of the gauge anomaly cancellation conditions. Such symmetry remains intact after the spontaneous symmetry breaking of the extra U ( 1 ) by the VEVs of right-handed sneutrinos.
As of this Z 2 symmetry, the superfields of type ξ ^ α can only appear in pairs in the Lagrangian. As a consequence, it is straightforward to realize that vanishing vevs for their scalar components, ξ α = 0 , is a solution of the minimization equations. Thus, the Z 2 symmetry is exact, and therefore we can have either the bosonic or the fermionic components of ξ ^ α as WIMP DM, without introducing R-parity.
Since the hierarchy m ξ α R > m ξ α I > m ξ ˜ α can be naturally satisfied [38], we use, in our analysis, the lightest of the fermionic components of the superfields, say ξ ˜ 1 , as the DM particle. The heaviest state ξ ˜ 2 can decay for example to ξ ˜ 2 ξ ˜ 1 q ¯ q and therefore does not play any role in the phenomenology of interest here. In what follows, we denote our DM candidate by ξ ˜ ξ ˜ 1 .

5.2. Constraints on the Parameter Space

This kind of DM interacts with the SM particle content via exchange of a new massive gauge boson Z , right sneutrinos, SM-like Higgs via scalar mixing, as well as DM exchange (see Ref. [38]). In this setup, SI (SD) DM-nucleon scatterings are mediated by Higgs via scalar mixing ( Z ), by interactions with light quarks within nucleons (see Figure 6). Therefore, DM direct detection experiments can probe regions of our parameter space.
We also point out that the exotic quarks offer an additional channel for SI scatterings by interacting directly with the gluons present in the nucleons and with DM by right sneutrino mediation. As the presence of these exotics is required by the anomaly cancellation conditions, their contribution is a rather general prediction of the U μ ν SSM. Although it turns out to be significant only in specific corners of the parameter space of our scan range, it offers nevertheless the possibility of testing a part of the parameters in the future in the case of low values of the scalar mixing.
Additional constraints on this scenario are imposed by Z LHC searches, which can exclude masses m Z 1 5 TeV depending on the value of the U ( 1 ) gauge coupling (see Ref. [38]), as well as R-hadron searches, which provide a lower bound on the masses of exotic quarks of the order of the TeV scale. Concerning LHC signals of the DM particle itself, the direct production is quite suppressed because it is a SM singlet. However, in regions of the parameter space where the singlets ξ ˜ are lighter than m Z / 2 , they could be produced in Z decays. The decay Z ξ ˜ 2 ξ ˜ 2 with subsequent decay ξ ˜ 2 ξ ˜ 1 + produces two pairs of collimated leptons, which can give striking signatures [149]. Other decay modes, such as ξ ˜ 2 ξ ˜ 1 q q ¯ are likely unobservable, as are the decays Z ξ ˜ 1 ξ ˜ 1 .

5.3. Analysis

We analyzed the possibility of reproducing the observed DM relic abundance via the freeze-out mechanism in this setup, performing a numerical study of the viable parameter space respecting all constraints. Results from Xenon1T experiment already exclude a subdominant portion of the allowed parameter space. We identified two main regions allowed by Xenon1T, at large DM masses m ξ ˜ 2–3 TeV and at smaller DM masses 200 GeV m ξ ˜ 2–3 TeV (see Figure 7).
For the case of large masses, the dynamical generation of the DM mass implies relatively large couplings with the right sneutrino. For this regime, new bosonic (right-handed sneutrinos ν ˜ R , heavy gauge boson Z ) and fermionic (exotic quarks K , neutralinos χ ˜ i and charginos χ ˜ ± ) states are the most frequent particles present in the final states of DM annihilation and therefore essential to reproduce the correct relic abundance.
For the highest masses, m ξ ˜ 10 5 GeV, the viable part of the parameter space requires couplings that are typically on the edge of perturbative unitarity. This part of the parameter space offers optimistic detection prospects as the Darwin experiment should probe the majority of the viable parameters in the following year and the remaining part should be accessible with an increased exposure.
For the region with lower masses, achieving the correct relic abundance is less frequent as most of the annihilation channels mentioned previously are kinematicaly forbidden after imposing constraints on the new states. This regime typically relies on s-channel ν ˜ R or Z resonances with SM particles in the final states, such as quarks. Relatively low couplings are typically required for such masses and therefore the direct detection prospects are less optimistic. Nevertheless, a substantial part of the parameter space will be accessible by the Darwin experiment.
Interestingly, as many annihilation channels are usually required to achieve the correct relic abundance, non-velocity suppressed DM annihilation within large astrophysical structures could offer complementary detection prospects by indirect gamma-ray searches with the upcoming Cherenkov Telescope Array (CTA) [153].

6. Conclusions

We reviewed DM analyses in various SUSY models that have been obtained in the recent years by groups (to a relevant part) based in Spain. The first model was the pMSSM11, the phenomenological MSSM with 11 free parameters, based on Ref. [41]. We found that this model gives a perfect description of all relevant experimental data, including LHC searches, the anomalous magnetic moment of the muon, the DM relic abundance and the DD limits. A large part of the parameter space will be accessible to future spin-independent DD experiments. However, even at the 1 σ level, points below the neutrino floor were found, making them difficult to access in future DD experiments.
The second model analyzed was the electroweak (EW)-MSSM [96], where the colored and the Higgs sector are assumed to be heavy, so that they do not play a role in the phenomenological analysis. Five different scenarios, depending on the mechanism that brings the DM relic abundance in agreement with the Planck limit, were identified. The wino and higgsino DM case can be covered by the next round of DD experiments. In the other three scenarios, two cases of bino DM and mixed bino/wino DM, large parts of the allowed parameter space can be covered by future Argon-based DD experiments. However, all three scenarios contain parameter points below the neutrino floor. For these points, it was shown that possible future e + e colliders with s < 1 TeV can conclusively test these scenarios (while this may not be possible at the HL-LHC).
The third model analyzed was the μ ν SSM [8], based on Ref. [27]. We showed first that the gravitino LSP is an interesting decaying DM candidate, with a lifetime larger than the age of the Universe. In order to avoid too large gamma-ray fluxes that are incompatible with Fermi-LAT observations, gravitino masses must be smaller than about 17 GeV and lifetimes larger than 4 × 10 25 s [29,30].
We also applied this constraint to the case of axino DM, obtaining that its mass must be smaller than about 3 GeV. In addition, a significant region of the parameter space of axino DM lies in the ballpark of future gamma-ray missions, such as the proposed e-ASTROGAM, allowing the exploration of masses and lifetimes in the ranges 2 MeV–3 GeV and 2 × 10 26 –8 × 10 30 s, respectively, [33]. The possibility of a multicomponent DM scenario made of gravitinos and axinos was also discussed [31,33].
We found that an axino or the gravitino can produce a signal detectable by future MeV-GeV gamma-ray telescopes, such as e-ASTROGAM. In addition, there is a parameter region where a well-tempered mixture of both particles is obtained, with a double-line signal arising as a smoking gun. Finally, we suggested that other DM candidates could be available in the μ ν SSM in certain regions of the parameter space, such as a sterile neutrino [147].
The last model analyzed was the U μ ν SSM [38], based on Ref. [40]. This model is a U ( 1 ) extension of the μ ν SSM, where the gauge anomaly-cancellation conditions impose the presence of exotic quark superfields in the spectrum and allow the presence of several singlet superfields under the SM gauge group, in addition to the right-handed neutrino superfields. The gauge structure implies an additional discrete Z 2 symmetry in the superpotential, ensuring the stability of a singlet, which behaves as WIMP DM without invoking R-parity.
We analyzed this novel possibility in detail, using the fermionic component of the singlet as the DM candidate. In particular, we computed its amount of relic density via Z , Higgs-right sneutrino and DM mediated annihilations and its potential signals in DM direct detection experiments. The constraints on the parameter space due to Z direct searches at the LHC were imposed in the analysis as well as those from the hadronization inside the detector of the exotic quarks. Large regions of the parameter space were found to be in the reach of the upcoming Darwin experiment.

Author Contributions

All authors contributed equally to this article. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the Spanish MICINN’s Consolider-Ingenio 2010 Program under Grant MultiDark No. CSD2009-00064, continued by the “Spanish Red Consolider MultiDark” FPA2017-90566-REDC, for long-time support of DM related investigations in Spain. This work is partially supported by the Spanish Research Agency through the grants IFT Centro de Excelencia Severo Ochoa No CEX2020-001007-S, PID2019-110058GB-C21 (S.H.) and PGC2018-095161-B-I00 (C.M.). They are funded by MCIN/AEI/10.13039/501100011033 and the last grant also by ERDF A way of making Europe. We thank J.A. Aguilar Saavedra, A. Albert, E. Bagnaschi, E. Bloom, M. Borsato, O. Buchmüller, M. Chakraborti, E. Charles, Y. Choi, M. Citron, J.C. Costa, A. J. Cuesta, A. De Roeck, M. Dolan, J.R. Ellis, H. Flächer, M. Fornasa, G.A. Gómez-Vargas, M. Grefe, D.E. López-Fogliani, M. Lucio, D. Martínez Santos, M. N. Mazziota, A. Morselli, K. Olive, A.D. Perez, M. Pierre, F. Prada, A. Richards, R. Ruiz de Austri, I. Saha, K. Sakurai, C. Schappacher, V. Spanos, I. Suárez Fernández, G. Weiglein, C. Weniger, G. Yepes and F. Zandanel with whom the results reviewed here were obtained.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nille, H. Supersymmetry, supergravity and particle physics. Phys. Rept. 1984, 110, 1–162. [Google Scholar] [CrossRef]
  2. Barbieri, R. Looking beyond the standard model: The supersymmetric option. Riv. Nuovo Cim. 1988, 11, 1–45. [Google Scholar] [CrossRef]
  3. Haber, H.; Kane, G. The search for supersymmetry: Probing physics beyond the standard model. Phys. Rept. 1985, 117, 175–263. [Google Scholar] [CrossRef]
  4. Gunion, J.; Haber, H. Higgs bosons in supersymmetric models (I). Nucl. Phys. 1986, 272, 1–76. [Google Scholar]
  5. Goldberg, H. Constraint on the photino mass from cosmology. Phys. Rev. Lett. 1983, 50, 1419. [Google Scholar] [CrossRef]
  6. Ellis, J.; Hagelin, J.; Nanopoulos, D.; Olive, K.; Srednicki, M. Supersymmetric relics from the big bang. Nucl. Phys. 1984, 238, 453–476. [Google Scholar] [CrossRef]
  7. Bae, K.J.; Baer, H.; Chun, E.J. Mainly axion cold dark matter from natural supersymmetry. Phys. Rev. 2014, 89, 31701. [Google Scholar] [CrossRef]
  8. Lopez-Fogliani, D.E.; Muñoz, C. Proposal for a supersymmetric standard model. Phys. Rev. Lett. 2006, 97, 041801. [Google Scholar] [CrossRef]
  9. Muñoz, C. Phenomenology of a new supersymmetric standard model: The μνSSM. AIP Conf. Proc. 2010, 1200, 413–416. [Google Scholar]
  10. Lopez-Fogliani, D.E.; Muñoz, C. Searching for supersymmetry: The μνSSM. Eur. Phys. J. 2020, 229, 3263–3301. [Google Scholar] [CrossRef]
  11. Biekötter, T.; Heinemeyer, S.; Weiglein, G. Vacuum (meta-) stability in the μνSSM. Eur. Phys. J. C 2022, 82, 301. [Google Scholar] [CrossRef]
  12. Escudero, N.; Lopez-Fogliani, D.E.; Muñoz, C.; de Austri, R.R. Analysis of the parameter space and spectrum of the μνSSM. J. High Energy Phys. 2008, 12, 99. [Google Scholar] [CrossRef]
  13. Ghosh, P.; Roy, S. Neutrino masses and mixing, lightest neutralino decays and a solution to the µ problem in supersymmetry. J. High Energy Phys. 2009, 4, 69. [Google Scholar] [CrossRef]
  14. Fidalgo, J.; Lopez-Fogliani, D.E.; Muñoz, C.; de Austri, R.R. Neutrino physics and spontaneous CP violation in the μνSSM. J. High Energy Phys. 2009, 8, 105. [Google Scholar] [CrossRef]
  15. Ghosh, P.; Dey, P.; Mukhopadhyaya, B.; Roy, S. Radiative contribution to neutrino masses and mixing in μνSSM. J. High Energy Phys. 2010, 5, 1–50. [Google Scholar] [CrossRef]
  16. Ellis, J.R.; Gunion, J.F.; Haber, H.E.; Roszkowski, L.; Zwirner, F. Higgs bosons in a nonminimal supersymmetric model. Phys. Rev. D 1989, 9, 844. [Google Scholar] [CrossRef] [PubMed]
  17. Miller, D.J.; Nevzorov, R.; Zerwas, P.M. The Higgs sector of the next-to-minimal supersymmetric standard model. Nucl. Phys. B 2004, 681, 3–30. [Google Scholar] [CrossRef]
  18. Kim, J.E.; Nilles, H. The µ-problem and the strong CP-problem. Phys. Lett. B 1984, 138, 150–154. [Google Scholar] [CrossRef]
  19. Ghosh, P.; Lopez-Fogliani, D.E.; Mitsou, V.A.; Muñoz, C.; de Austri, R.R. Probing the μνSSM with light scalars, pseudoscalars and neutralinos from the decay of a SM-like Higgs boson at the LHC. J. High Energy Phys. 2014, 11, 1–57. [Google Scholar] [CrossRef]
  20. Ghosh, P.; Lara, I.; Lopez-Fogliani, D.E.; Muñoz, C.; de Austri, R.R. Searching for left sneutrino LSP at the LHC. Int. J. Mod. Phys. A 2018, 33, 1850110. [Google Scholar] [CrossRef]
  21. Lara, I.; López-Fogliani, D.E.; Muñoz, C.; Nagata, N.; Otono, H.; Austri, R.R.D. Looking for the left sneutrino LSP with displaced-vertex searches. Phys. Rev. D 2018, 98, 75004. [Google Scholar] [CrossRef]
  22. Lara, I.; López-Fogliani, D.E.; Muñoz, C. Electroweak superpartners scrutinized at the LHC in events with multi-leptons. Phys. Lett. B 2019, 790, 176–183. [Google Scholar] [CrossRef]
  23. Kpatcha, E.; Lara, I.; López-Fogliani, D.E.; Muñoz, C.; Nagata, N.; Otono, H.; de Austri, R.R. Sampling the μνSSM for displaced decays of the tau left sneutrino LSP at the LHC. Eur. Phys. J. C 2019, 79, 1–18. [Google Scholar] [CrossRef]
  24. Kpatcha, E.; López-Fogliani, D.E.; Muñoz, C.; de Austri, R.R. Impact of Higgs physics on the parameter space of the μνSSM. Eur. Phys. J. C 2020, 80, 1–43. [Google Scholar] [CrossRef]
  25. Kpatcha, E.; Lara, I.; López-Fogliani, D.E.; Muñoz, C.; Nagata, N. Explaining muon g-2 data in the μνSSM. Eur. Phys. J. C 2021, 81, 1–22. [Google Scholar] [CrossRef]
  26. Heinemeyer, S.; Kpatcha, E.; Lara, I.; López-Fogliani, D.E.; Muñoz, C.; Nagata, N. The new (g-2)μ result and the μνSSM. Eur. Phys. J. C 2021, 81, 802. [Google Scholar] [CrossRef]
  27. Choi, K.-Y.; López-Fogliani, D.E.; Muñoz, C.; de Austri, R.R. Gamma-ray detection from gravitino dark matter decay in the μνSM. J. Cosmol. Astropart. Phys. 2010, 3, 028. [Google Scholar] [CrossRef]
  28. Gómez-Vargas, G.A.; Fornasa, M.; Zandanel, F.; Cuesta, A.J.; Muñoz, C.; Prada, F.; Yepes, G. CLUES on Fermi-LAT prospects for the extragalactic detection of μνSSM gravitino dark matter. J. Cosmol. Astropart. Phys. 2012, 2, 1. [Google Scholar] [CrossRef]
  29. Albert, A.; Gómez-Vargas, G.A.; Grefe, M.; Muñoz, C.; Weniger, C.; Bloom, E.; Charles, E.; Mazziota, M.N.; Morselli, A. Search for 100 MeV to 10 GeV γ -ray lines in the Fermi-LAT data and implications for gravitino dark matter in the μνSSM. J. Cosmol. Astropart. Phys. 2014, 10, 23. [Google Scholar] [CrossRef]
  30. Gómez-Vargas, G.A.; López-Fogliani, D.E.; Muñoz, C.; de Austri, R.R. Search for sharp and smooth spectral signatures of μνSSM gravitino dark matter with Fermi-LAT. J. Cosmol. Astropart. Phys. 2017, 3, 47. [Google Scholar] [CrossRef]
  31. Gómez-Vargas, G.A.; López-Fogliani, D.E.; Muñoz, C.; Perez, A.D. MeV-GeV γ-ray telescopes probing gravitino LSP with coexisting axino NLSP as dark matter in the μνSSM. Astropart. Phys. 2021, 125, 102506. [Google Scholar] [CrossRef]
  32. Restrepo, D.; Taoso, M.; Valle, J.; Zapata, O. Gravitino dark matter and neutrino masses with bilinear R-parity violation. Phys. Rev. D 2012, 85, 023523. [Google Scholar] [CrossRef]
  33. Gómez-Vargas, G.A.; López-Fogliani, D.E.; Muñoz, C.; Perez, A.D. MeV-GeV γ-ray telescopes probing axino LSP/gravitino NLSP as dark matter in the μνSSM. J. Cosmol. Astropart. Phys. 2020, 1, 58. [Google Scholar] [CrossRef]
  34. Chung, D.J.H.; Long, A.J. Electroweak phase transition in the μνSSM. Phys. Rev. D 2010, 81, 123531. [Google Scholar] [CrossRef]
  35. Cerdeno, D.G.; Muñoz, C.; Seto, O. Right-handed sneutrino as thermal dark matter. Phys. Rev. D 2009, 79, 023510. [Google Scholar] [CrossRef]
  36. Cerdeno, D.G.; Peiro, M.; Robles, S. Fits to the Fermi-LAT GeV excess with right-handed sneutrino dark matter: Implications for direct and indirect dark matter searches and the LHC. Phys. Rev. D 2015, 91, 123530. [Google Scholar] [CrossRef]
  37. López-Fogliani, D.E.; Perez, A.D.; de Austri, R.R. Dark matter candidates in the NMSSM with RH neutrino superfields. J. Cosmol. Astropart. Phys. 2021, 4, 67. [Google Scholar] [CrossRef]
  38. Aguilar-Saavedra, J.A.; Lara, I.; López-Fogliani, D.E.; Muñoz, C. U(1) extensions of the μνSSM. Eur. Phys. J. C 2021, 81, 443. [Google Scholar]
  39. Aguilar-Saavedra, J.A.; Lara, I.; López-Fogliani, D.E.; Muñoz, C. Exotic diboson Z’ decays in the U μν SSM. Eur. Phys. J. C 2021, 81, 805. [Google Scholar] [CrossRef]
  40. Aguilar-Saavedra, J.A.; López-Fogliani, D.E.; Muñoz, C.; Pierre, M. WIMP dark matter in the UμνSSM. J. Cosmol. Astropart. Phys. 2022, 5, 4. [Google Scholar] [CrossRef]
  41. Bagnaschi, E.; Sakurai, K.; Borsato, M.; Buchmueller, O.; Citron, M.; Costa, J.C.; De Roeck, A.; Dolan, M.J.; Ellis, J.R.; Flacher, S.; et al. Likelihood analysis of the pMSSM11 in light of LHC 13-TeV data. Eur. Phys. J. C 2018, 78, 56. [Google Scholar] [CrossRef] [PubMed]
  42. AbdusSalam, S.S.; Allanach, B.C.; Dreiner, H.K.; Ellis, J.; Ellwanger, U.; Gunion, J.; Heinmeyer, S.; Kraemer, M.; Mangano, M.L.; Olive, K.A.; et al. Benchmark models, planes, lines and points for future SUSY searches at the LHC. Eur. Phys. J. C 2011, 71, 1835. [Google Scholar] [CrossRef]
  43. Buchmueller, O.; Cavanaugh, R.; De Roeck, A.; Heinmeyer, S.; Isidori, G.; Ronga, F.J.; Weber, A.M.; Weiglein, G. Prediction for the lightest Higgs boson mass in the CMSSM using indirect experimental constraints. Phys. Lett. B 2007, 657, 87–94. [Google Scholar] [CrossRef]
  44. Buchmueller, O.; Cavanaugh, R.; De Roeck, A.; Heinmeyer, S.; Isidori, G.; Ronga, F.J.; Weber, A.M.; Weiglein, G. Supersymmetry in light of 1/fb of LHC data. Eur. Phys. J. 2012, 72, 1878. [Google Scholar] [CrossRef]
  45. Buchmueller, O.; Cavanaugh, R.; De Roeck, A.; Heinmeyer, S.; Isidori, G.; Ronga, F.J.; Weber, A.M.; Weiglein, G. The CMSSM and NUHM1 in light of 7 TeV LHC, Bsμ+μ and XENON100 data. Eur. Phys. J. 2012, 72, 2243. [Google Scholar] [CrossRef]
  46. Buchmueller, O.; Cavanaugh, R.; De Roeck, A.; Heinmeyer, S.; Isidori, G.; Ronga, F.J.; Weber, A.M.; Weiglein, G. Implications of improved Higgs mass calculations for supersymmetric models. Eur. Phys. J. 2014, 74, 2809. [Google Scholar] [CrossRef] [PubMed]
  47. Buchmueller, O.; Cavanaugh, R.; De Roeck, A.; Dolan, M.J.; Ellis, J.R.; Flacher, H.; Heinemeyer, S.; Isidori, G.; Marrouche, J.; Martinez Santos, D.; et al. The CMSSM and NUHM1 after LHC Run 1. Eur. Phys. J. C 2014, 74, 2922. [Google Scholar] [CrossRef]
  48. Buchmueller, O.; Cavanaugh, R.; Citron, M.; De Roeck, A.; Dolan, M.J.; Ellis, J.R.; Flacher, H.; Heinemeyer, S.; Malik, J.; Marrouche, J.; et al. The NUHM2 after LHC Run 1. Eur. Phys. J. C 2014, 12, 3212. [Google Scholar]
  49. De Vries, K.J.; Bagnaschi, E.A.; Buchmueller, O.; Cavanaugh, R.; Citron, M.; De Roeck, A.; Dolan, M.J.; Ellis, J.R.; Flacher, H.; Heinemeyer, S.; et al. The pMSSM10 after LHC run 1. Eur. Phys. J. C 2015, 75, 422. [Google Scholar]
  50. Bagnaschi, E.A.; Buchmueller, O.; Cavanaugh, R.; Citron, M.; De Roeck, A.; Dolan, M.J.; Ellis, J.R.; Flacher, H.; Heinemeyer, S.; Isidori, G.; et al. Supersymmetric dark matter after LHC run 1. Eur. Phys. J. C 2015, 75, 500. [Google Scholar] [CrossRef]
  51. Bagnaschi, E.; Costa, J.C.; Sakurai, K.; Borsato, M.; Buchmueller, O.; Cavanaugh, R.; Chobanova, V.; Citron, M.; De Roeck, A.; Dolan, M.J.; et al. Likelihood analysis of supersymmetric SU(5) GUTs. Eur. Phys. J. C 2017, 77, 104. [Google Scholar] [CrossRef] [PubMed]
  52. Bagnaschi, E.; Borsato, M.; Sakurai, K.M.; Buchmueller, O.; Cavanaugh, R.; Chobanova, V.; Citron, M.; Costa, J.C.; De Roeck, A.; Dolan, M.J.; et al. Likelihood analysis of the minimal AMSB model. Eur. Phys. J. C 2017, 77, 268. [Google Scholar] [CrossRef] [PubMed]
  53. Skands, P.; Allanach, B.C.; Baer, H.; Balazs, C.; Belanger, G.; Boudjema, F.; Djouadi, A.; Godbole, R.; Guasch, J.; Heinemeyer, S.; et al. SUSY Les Houches accord: Interfacing SUSY spectrum calculators, decay packages, and event generators. J. High Energy Phys. 2004, 407, 036. [Google Scholar] [CrossRef]
  54. Allanach, B.C. SOFTSUSY: A program for calculating supersymmetric spectra. Comput. Phys. Commun. 2002, 143, 305–331. [Google Scholar] [CrossRef]
  55. Heinemeyer, S.; Hollik, W.; Stockinger, D.; Weber, A.M.; Weiglein, G. Precise prediction for MW in the MSSM. J. High Energy Phys. 2006, 0608, 052. [Google Scholar] [CrossRef]
  56. Heinemeyer, S.; Hollik, W.; Stockinger, D.; Weber, A.M.; Weiglein, G. Z pole observables in the MSSM. J. High Energy Phys. 2008, 0804, 039. [Google Scholar] [CrossRef]
  57. Heinemeyer, S.; Hollik, W.; Weiglein, G.; Zeune, L. Implications of LHC search results on the W boson mass prediction in the MSSM. J. High Energy Phys. 2013, 1312, 84. [Google Scholar] [CrossRef]
  58. Heinemeyer, S.; Hollik, W.; Weiglein, G. FeynHiggs: A program for the calculation of the masses of the neutral -even Higgs bosons in the MSSM. Comput. Phys. Commun. 2000, 124, 76–89. [Google Scholar] [CrossRef]
  59. Heinemeyer, S.; Hollik, W.; Weiglein, G. The masses of the neutral CP-even Higgs bosons in the MSSM: Accurate analysis at the two-loop level. Eur. Phys. J. 1999, 9, 343–366. [Google Scholar] [CrossRef]
  60. Degrassi, G.; Heinemeyer, S.; Hollik, W.; Weiglein, P.S.G. Towards high-precision predictions for the MSSM Higgs sector. Eur. Phys. J. 2003, 28, 133–143. [Google Scholar] [CrossRef]
  61. Frank, M.; Hahn, T.; Heinemeyer, S.; Hollik, W.; Rzehak, H.; Weiglein, G. LHC physics. JHEP 2007, 47, 0702. [Google Scholar]
  62. Hahn, T.; Heinemeyer, S.; Hollik, W.; Rzehak, H.; Weiglein, G. FeynHiggs: A program for the calculation of MSSM Higgs-boson observables—Version 2.6.5. Comput. Phys. Commun. 2009, 180, 1426–1427. [Google Scholar] [CrossRef]
  63. Hahn, T.; Heinemeyer, S.; Hollik, W.; Rzehak, H.; Weiglein, G. High-Precision Predictions for the Light CP-Even Higgs Boson Mass of the Minimal Supersymmetric Standard Model. Phys. Rev. Lett. 2014, 112, 141801. [Google Scholar] [CrossRef]
  64. Bahl, H.; Hollik, W. Precise prediction for the light MSSM Higgs-boson mass combining effective field theory and fixed-order calculations. Eur. Phys. J. 2016, 76, 499. [Google Scholar] [CrossRef]
  65. Bahl, H.; Heinemeyer, S.; Hollik, W.; Weiglein, G. Reconciling EFT and hybrid calculations of the light MSSM Higgs-boson mass. Eur. Phys. J. 2018, 78, 57. [Google Scholar] [CrossRef]
  66. Bahl, H.; Hahn, T.; Heinemeyer, S.; Hollik, W.; Paßehr, S.; Rzehak, H.; Weiglein, G. Precision calculations in the MSSM Higgs-boson sector with FeynHiggs 2.14. Comput. Phys. Commun. 2020, 249, 107099. [Google Scholar] [CrossRef]
  67. Bahl, H.; Heinemeyer, S.; Hollik, W.; Weiglein, G. Theoretical uncertainties in the MSSM Higgs boson mass calculation. Eur. Phys. J. 2020, 80, 497. [Google Scholar] [CrossRef]
  68. Isidori, G.; Paradisi, P. Hints of large tan β in flavour physics. Phys. Lett. B 2006, 639, 499–507. [Google Scholar] [CrossRef]
  69. Isidori, G.; Mescia, F.; Paradisi, P.; Temes, D. Flavor physics at large tanβ with a binolike lightest supersymmetric particle. Phys. Rev. D 2007, 2007 75, 115019. [Google Scholar] [CrossRef]
  70. Mahmoudi, F. SuperIso: A program for calculating the isospin asymmetry of BK*γ in the MSSM. Comput. Phys. Commun. 2008, 178, 745–754. [Google Scholar] [CrossRef]
  71. Mahmoudi, F. SuperIso v3.0: A program for calculating flavor physics observables in 2HDM and supersymmetry. Comput. Phys. Commun. 2009, 180, 1579. [Google Scholar] [CrossRef]
  72. Eriksson, D.; Mahmoudi, F.; Stål, O. Charged Higgs bosons in minimal supersymmetry: Updated constraints and experimental prospects. JHEP 2008, 035. [Google Scholar] [CrossRef]
  73. Belanger, G.; Boudjema, F.; Pukhov, A.; Semenov, A. micrOMEGAs_3: A program for calculating dark matter observables. Comput. Phys. Commun. 2014, 185, 960–985. [Google Scholar] [CrossRef]
  74. Information about this code is available from K.A. Olive: It contains important contributions from Evans, J.; Falk, T.; Ferstl, A.; Ganis, G.; Luo, F.; Mustafayev, A.; McDonald, J.; Luo, F.; Olive, K.A.; Sandick, P.; Santoso, Y.; Savage, C.; Spanos, V.; Srednicki, M. Available online: https://link.springer.com/article/10.1140/epjc/s10052-018-5633-3 (accessed on 28 July 2022).
  75. Muhlleitner, M.; Djouadi, A.; Mambrini, Y. SDECAY: A Fortran code for the decays of the supersymmetric particles in the MSSM. Comput. Phys. Commun. 2005, 168, 46–70. [Google Scholar] [CrossRef]
  76. Bechtle, P.; Heinemeyer, S.; Staal, O.; Stefaniak, T.; Weiglein, G. HiggsSignals: Confronting arbitrary Higgs sectors with measurements at the Tevatron and the LHC. Eur. Phys. J. C 2014, 74, 2711. [Google Scholar] [CrossRef]
  77. Bechtle, P.; Heinemeyer, S.; Stål, O.; Stefaniak, T.; Weiglein, G. Probing the Standard Model with Higgs signal rates from the Tevatron, the LHC and a future ILC. JHEP 2014, 1411, 39. [Google Scholar] [CrossRef]
  78. Bechtle, P.; Heinemeyer, S.; Staal, O.; Stefaniak, T.; Weiglein, G. HiggsSignals-2: Probing new physics with precision Higgs measurements in the LHC 13 TeV era. Phys. J. C 2021, 81, 145. [Google Scholar] [CrossRef]
  79. Bechtle, P.; Brein, O.; Heinemeyer, S.; Weiglein, G.; Williams, K.E. HiggsBounds: Confronting arbitrary Higgs sectors with exclusion bounds from LEP and the Tevatron. Comput. Phys. Commun. 2010, 181, 138–167. [Google Scholar] [CrossRef]
  80. Bechtle, P.; Brein, O.; Heinemeyer, S.; Weiglein, G.; Williams, K.E. HiggsBounds 2.0.0: Confronting neutral and charged Higgs sector predictions with exclusion bounds from LEP and the Tevatron. Comput. Phys. Commun. 2011, 181, 2605–2631. [Google Scholar] [CrossRef]
  81. Bechtle, P.; Brein, O.; Heinemeyer, S.; Weiglein, G.; Williams, K.E. HiggsBounds-4: Improved tests of extended Higgs sectors against exclusion bounds from LEP, the Tevatron and the LHC. Eur. Phys. J. C 2014, 74, 2693. [Google Scholar] [CrossRef]
  82. Bechtle, P.; Brein, O.; Heinemeyer, S.; Weiglein, G.; Williams, K.E. Applying exclusion likelihoods from LHC searches to extended Higgs sectors. Eur. Phys. J. C 2015, 75, 421. [Google Scholar] [CrossRef]
  83. Bechtle, P.; Dercks, D.; Heinemeyer, S.; Klingl, T.; Stefaniak, T.; Weiglein, G.; Wittbrodt, J. HiggsBounds-5: Testing Higgs sectors in the LHC 13 TeV Era. Eur. Phys. J. C 2020, 12, 1211. [Google Scholar] [CrossRef]
  84. [PandaX-II Collaboration]. Dark Matter Results from 54-Ton-Day Exposure of PandaX-II Experiment. Phys. Rev. Lett. 2017, 119, 181302. [Google Scholar] [CrossRef] [PubMed]
  85. [PandaX-II Collaboration]. Dark Matter Results from First 98.7 Days of Data from the PandaX-II Experiment. Phys. Rev. Lett. 2016, 117, 121303. [Google Scholar] [CrossRef] [PubMed]
  86. [XENON Collaboration]. First Dark Matter Search Results from the XENON1T Experiment. Phys. Rev. Lett. 2017, 119, 181301. [Google Scholar] [CrossRef] [PubMed]
  87. [LUX Collaboration]. Results from a Search for Dark Matter in the Complete LUX Exposure. Phys. Rev. Lett. 2017, 118, 021303. [Google Scholar] [CrossRef] [PubMed]
  88. [XENON Collaboration]. Physics reach of the XENON1T dark matter experiment. J. Cosmol. Astropart. Phys. 2016, 1604, 27. [Google Scholar]
  89. Mount, B.J.; Hans, S.; Rosero, R.; Yeh, M.; Chan, C.; Gaitskell, R.J.; Huang, D.Q.; Makkinje, J.; Malling, D.C.; Pangilinan, M.; et al. LUX-ZEPLIN (LZ) Technical Design Report. arXiv 2017, arXiv:1703.09144. [Google Scholar]
  90. Bernabei, R.; Belli, P.; Caracciolo, V.; Cerulli, R.; Merlo, V.; Cappella, F.; d’Angelo, A.; Incicchitti, A.; Dai, C.J.; Ma, X.H.; et al. The dark matter: DAMA/LIBRA and its perspectives. arXiv 2021, arXiv:2110.04734. [Google Scholar]
  91. Amaré, J.; Cebrián, S.; Cintas, D.; Coarasa, I.; García, E.; Martínez, M.; Oliván, M.A.; Ortigoza, Y.; de Solórzano, A.O.; Puimedón, J.; et al. Annual modulation results from three-year exposure of ANAIS-112. Phys. Rev. D 2021, 103, 102005. [Google Scholar] [CrossRef]
  92. [PICO Collaboration]. Dark matter search results from the PICO-60 CF3I bubble chamber. Phys. Rev. D 2016, 93, 052014. [Google Scholar] [CrossRef]
  93. Ng, K.C.Y.; Beacom, J.F.; Peter, A.H.G.; Rott, C. Solar atmospheric neutrinos: A new neutrino floor for dark matter searches. Phys. Rev. D 2017, 96, 103006. [Google Scholar] [CrossRef]
  94. [IceCube Collaboration]. Search for annihilating dark matter in the Sun with 3 years of IceCube data. Eur. Phys. J. C 2017, 77, 146. [Google Scholar] [CrossRef]
  95. [Super-Kamiokande Collaboration]. Search for Neutrinos from Annihilation of Captured Low-Mass Dark Matter Particles in the Sun by Super-Kamiokande. Phys. Rev. Lett. 2015, 114, 141301. [Google Scholar] [CrossRef]
  96. Chakraborti, M.; Heinemeyer, S.; Saha, I.; Schappacher, C. (g - 2)μ and SUSY Dark Matter: Direct Detection and Collider Search Complementarity. Eur. Phys. J. C 2022, 82, 483. [Google Scholar] [CrossRef]
  97. Chakraborti, M.; Heinemeyer, S.; Saha, I. Improved (g - 2)μ measurements and supersymmetry. Eur. Phys. J. 2020, 80, 984. [Google Scholar]
  98. Slavich, P.; Heinemeyer, S.; Bagnaschi, E.; Bahl, H.; Goodsell, M.; Haber, H.E.; Hahn, T.; Harlander, R.; Hollik, W.; Lee, G.; et al. Higgs-mass predictions in the MSSM and beyond. Eur. Phys. J. 2021, 81, 450. [Google Scholar] [CrossRef]
  99. Chakraborti, M.; Heinemeyer, S.; Saha, I. Improved (g - 2)μ measurements and wino/higgsino dark matter. Eur. Phys. J. 2021, 8, 1069. [Google Scholar] [CrossRef]
  100. Chakraborti, M.; Heinemeyer, S.; Saha, I. The new “MUON G-2” result and supersymmetry. Eur. Phys. J. 2021, 81, 1114. [Google Scholar] [CrossRef]
  101. Aoyama, T.; Aoyama, T.; Asmussen, N.; Benayoun, M.; Bijnens, J.; Blum, T.; Bruno, M.; Caprini, I.; Calame, C.C.; Cè, M.; et al. The anomalous magnetic moment of the muon in the Standard Model. Phys. Rept. 2020, 887, 1–166. [Google Scholar] [CrossRef]
  102. [Muon g-2 Collaboration]. Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm. Phys. Rev. Lett. 2021, 126, 141801. [Google Scholar] [CrossRef] [PubMed]
  103. [Muon g-2 Collaboration]. Final report of the E821 muon anomalous magnetic moment measurement at BNL. Phys. Rev. 2006, 73, 72003. [Google Scholar]
  104. Athron, P.; Bach, M.; Fargnoli, H.G.; Gnendiger, C.; Greifenhagen, R.; Park, J.H.; Pabehr, S.; Stockinger, D.; Stockinger-Kim, H. GM2Calc: Precise MSSM prediction for (g - 2) of the muon. Eur. Phys. J. 2016, 76, 62. [Google Scholar]
  105. Von Weitershausen, P.; Schafer, M.; Stöckinger-Kim, H.; Stöckinger, D. Photonic supersymmetric two-loop corrections to the muon magnetic moment. Phys. Rev. 2010, 81, 093004. [Google Scholar] [CrossRef]
  106. Fargnoli, H.; Gnendiger, C.; Pabehr, S.; Stockinger, D.; Stockinger-Kim, H. Two-loop corrections to the muon magnetic moment from fermion/sfermion loops in the MSSM: Detailed results. J. High Energy Phys. 2014, 1402, 70. [Google Scholar] [CrossRef]
  107. Bach, M.; Park, J.H.; Stockinger, D.; Stockinger-Kim, H. Large muon (g − 2) with TeV-scale SUSY masses for tan β→∞. J. High Energy Phys. 2015, 1510, 26. [Google Scholar] [CrossRef]
  108. Heinemeyer, S.; Stöckinger, D.; Weiglein, G. Two-loop SUSY corrections to the anomalous magnetic moment of the muon. Nucl. Phys. 2004, 690, 62–80. [Google Scholar] [CrossRef]
  109. Heinemeyer, S.; Stöckinger, D.; Weiglein, G. Electroweak and supersymmetric two-loop corrections to (g - 2)μ. Nucl. Phys. 2004, 699, 103–123. [Google Scholar] [CrossRef]
  110. Hollik, W.G.; Weiglein, G.; Wittbrodt, J. Impact of vacuum stability constraints on the phenomenology of supersymmetric models. J. High Energy Phys. 2019, 3, 109. [Google Scholar] [CrossRef]
  111. Ferreira, P.M.; Mühlleitner, M.; Santos, R.; Weiglein, G.; Wittbrodt, J. Vacuum instabilities in the N2HDM. J. High Energy Phys. 2019, 9, 6. [Google Scholar] [CrossRef]
  112. Drees, M.; Dreiner, H.; Schmeier, D.; Tattersall, J.; Kim, J.S. CheckMATE: Confronting your favourite new physics model with LHC data. Comput. Phys. Commun. 2015, 187, 227–265. [Google Scholar] [CrossRef]
  113. Kim, J.S.; Schmeier, D.; Tattersall, J.; Rolbiecki, K. A framework to create customised LHC analyses within CheckMATE. Comput. Phys. Commun. 2015, 196, 535–562. [Google Scholar] [CrossRef]
  114. Dercks, D.; Desai, N.; Kim, J.S.; Rolbiecki, K.; Tattersall, J.; Weber, T. CheckMATE 2: From the model to the limit. Comput. Phys. Commun. 2017, 221, 383–418. [Google Scholar] [CrossRef]
  115. [Planck Collaboration]. Planck 2018 results-VI. Cosmological parameters. Astron. Astrophys. 2020, 641, A6. [Google Scholar] [CrossRef]
  116. [XENON Collaboration]. Projected WIMP sensitivity of the XENONnT dark matter experiment. J. Cosmol. Astropart. Phys. 2020, 11, 31. [Google Scholar]
  117. Akerib, D.S.; Akerlof, C.W.; Alsum, S.K.; Araújo, H.M.; Arthurs, M.; Bai, X.; Bailey, A.J.; Balajthy, J.; Balashov, S.; Bauer, D.; et al. Projected WIMP sensitivity of the LUX-ZEPLIN dark matter experiment. Phys. Rev. 2020, 101, 052002. [Google Scholar] [CrossRef]
  118. Aalseth, C.E.; Acerbi, F.; Agnes, P.; Albuquerque, I.F.M.; Alexander, T.; Alici, A.; Alton, A.K.; Antonioli, P.; Arcelli, S.; Ardito, R.; et al. DarkSide-20k: A 20 tonne two-phase LAr TPC for direct dark matter detection at LNGS. Eur. Phys. J. Plus 2018, 133, 131. [Google Scholar] [CrossRef]
  119. Billard, J.; Boulay, M.; Cebrián, S.; Covi, L.; Fiorillo, G.; Green, A.M.; Kopp, J.; Majorovits, B.; Palladino, K.; Petricca, F.; et al. Direct detection of dark matter—APPEC committee report. arXiv 2021, arXiv:2104.07634. [Google Scholar] [CrossRef]
  120. Ruppin, F.; Billard, J.; Figueroa-Feliciano, E.; Strigari, L. Complementarity of dark matter detectors in light of the neutrino background. Phys. Rev. 2014, 90, 083510. [Google Scholar] [CrossRef]
  121. Baer, H.; Barklow, T.; Fuji, K.; Gao, Y.; Hoang, A.; Kanemura, S.; List, J.; Logan, H.E.; Nomerotski, A.; Perelstein, M. The International Linear Collider Technical Design Report-Volume 2: Physics. CERN: Geneva, Switzerland, 2013. arXiv 2013, arXiv:1306.6352. [Google Scholar]
  122. Moortgat-Pick, G.; Baer, H.; Battaglia, M.; Belanger, G.; Fujii, K.; Kalinowski, J.; Heinemeyer, S.; Kiyo, Y.; Olive, K.; Simon, F.; et al. Physics at e+e Linear Collider. Eur. Phys. J. 2015, 75, 371. [Google Scholar] [CrossRef] [PubMed]
  123. Barklow, T.; Brau, J.; Fujii, K.; Gao, J.; List, J.; Walker, N.; Yokoya, K. ILC Operating Scenarios. arXiv 2015, arXiv:1506.07830. [Google Scholar]
  124. Fujii, K.; Grojean, C.; Peskin, M.E.; Barklow, T.; Gao, Y.; Kanemura, S.; Kim, H.; List, J.; Nojiri, M.; Perelstein, M.; et al. Physics Case for the 250 GeV Stage of the International Linear Collider. arXiv 2017, arXiv:1710.07621. [Google Scholar]
  125. Heinemeyer, S.; Schappacher, C. Chargino and neutralino production at e+e colliders in the complex MSSM: A full one-loop analysis. Eur. Phys. J. 2017, 77, 649. [Google Scholar] [CrossRef]
  126. Heinemeyer, S.; Schappacher, C. Slepton production at e+e colliders in the complex MSSM: A full one-loop analysis. Eur. Phys. J. 2018, 78, 536. [Google Scholar] [CrossRef]
  127. Berggren, M. Simplified SUSY at the ILC. arXiv 2021, arXiv:1308.1461. [Google Scholar]
  128. De Vera, M.T.N.P.; Berggren, M.; List, J. Chargino production at the ILC. arXiv 2019, arXiv:2002.01239. [Google Scholar]
  129. Berggren, M. What pp SUSY limits mean for future e+e- colliders. arXiv 2022, arXiv:2003.12391. [Google Scholar]
  130. Pagels, H.; Primack, J.R. Supersymmetry, Cosmology, and New Physics at Teraelectronvolt Energies. Phys. Rev. Lett. 1982, 48, 223. [Google Scholar] [CrossRef]
  131. Weinberg, S. Cosmological Constraints on the Scale of Supersymmetry Breaking. Phys. Rev. Lett. 1982, 48, 1303. [Google Scholar] [CrossRef]
  132. Khlopov, M.Y.; Linde, A. Is it easy to save the gravitino? Phys. Lett. B. 1984, 138, 265. [Google Scholar] [CrossRef]
  133. Ellis, J.R.; Kim, J.E.; Nanopoulos, D.V. Cosmological gravitino regeneration and decay. Phys. Rev. Lett. 1984, 145, 181–186. [Google Scholar] [CrossRef]
  134. Ellis, J.R.; Olive, K.A.; Santoso, Y.; Spanos, V.C. Gravitino dark matter in the CMSSM. Phys. Rev. Lett. B 2004, 588, 7–16. [Google Scholar] [CrossRef]
  135. Borgani, S.; Masiero, A.; Yamaguchi, M. Light gravitinos as mixed dark matter. Phys. Rev. Lett. B 1996, 386, 189–197. [Google Scholar] [CrossRef]
  136. Takayama, F.; Yamaguchi, M. Gravitino dark matter without R-parity. Phys. Rev. Lett. B 2000, 485, 388–392. [Google Scholar] [CrossRef]
  137. Addazi, A.; Khlopov, M.Y. Dark matter and inflation in R+ζR2 supergravity. Phys. Lett. B 2017, 766, 17. [Google Scholar] [CrossRef]
  138. Fermi-LAT Collaboration; Atwood, W.B.; Abdo, A.A.; Ackermann, M.; Althouse, W.; Anderson, B.; Axelsson, M.; Baldini, L.; Ballet, J.; Band, D.L.; et al. The Large Area Telescope on the Fermi Gamma-ray Space Telescope Mission. Astrophys. J. 2009, 697, 1071. [Google Scholar] [CrossRef]
  139. De Angelis, A.; Tatischeff, V.; Grenier, I.A.; McEnery, J.; Mallamaci, M.; Tavani, M.; Oberlack, U.; Hanlon, L.; Walter, R.; Argan, A.; et al. Science with e-ASTROGAM (A space mission for MeV-GeV gamma-ray astrophysics). arXiv 2021, arXiv:1711.01265. [Google Scholar]
  140. Caputo, R.; Meyer, M.; Sanchez-Conde, M. AMEGO: Dark Matter Prospects. In Proceedings of the 35th International Cosmic Ray Conference (PoS ICRC), Hamme, Belgium, 12–20 July 2017; p. 910. [Google Scholar]
  141. Shifman, M.A.; Vainshtein, A.; Zakharov, V.I. Can confinement ensure natural CP invariance of strong interactions? Nucl. Phys. B 1980, 166, 493–506. [Google Scholar] [CrossRef]
  142. Brandenburg, A.; Steffen, F.D. Axino dark matter from thermal production. J. Cosmol. Astropart. Phys. 2004, 8, 8. [Google Scholar] [CrossRef]
  143. Dine, M.; Fischler, W.; Srednicki, M. A simple solution to the strong CP problem with a harmless axion. Phys. Lett. B 1981, 104, 199–202. [Google Scholar] [CrossRef]
  144. Zhitnitsky, A.R. On possible supression of the axion hadron interactions. Sov. J. Nucl. Phys. 1980, 31, 260. [Google Scholar]
  145. Boyarsky, A.; Drewes, M.; Lasserre, T.; Mertens, S.; Ruchayskiy, O. Sterile neutrino dark matter. Prog. Part. Nucl. Phys. 2019, 104, 1–45. [Google Scholar] [CrossRef]
  146. Asaka, T.; Shaposhnikov, M. The νMSM, dark matter and baryon asymmetry of the universe. Phys. Lett. B 2005, 620, 17–26. [Google Scholar] [CrossRef]
  147. Knees, P.; Lopez-Fogliani, D.E.; Muñoz, C. Phenomenological implications of sterile neutrinos in the μνSSM and dark matter. To be published in Astropart. Phys. arXiv 2022, arXiv:2207.10689. [Google Scholar]
  148. Bringmann, T.; Kahlhoefer, F.; Schmidt-Hoberg, K.; Walia, P. Converting nonrelativistic dark matter to radiation. Phys. Rev. D 2018, 98, 023543. [Google Scholar] [CrossRef]
  149. Aguilar-Saavedra, J.A.; Casas, J.; Quilis, J.; de Austri, R.R. Multilepton dark matter signals. J. High Energy Phys. 2020, 4, 69. [Google Scholar] [CrossRef]
  150. [XENON Collaboration]. Dark Matter Search Results from a One Ton-Year Exposure of XENON1T. Phys. Rev. Lett. 2018, 121, 111302. [Google Scholar] [CrossRef]
  151. Aalbers, J.; Agostini, F.; Alfonsi, M.; Amaro, F.D.; Amsler, C.; Aprile, E.; Arazi, L.; Arneodo, F.; Barrow, P.; Baudis, L.; et al. DARWIN: Towards the ultimate dark matter detector. J. Cosmol. Astropart. Phys. 2016, 11, 17. [Google Scholar] [CrossRef]
  152. Billard, J.; Figueroa-Feliciano, E.; Strigari, L. Implication of neutrino backgrounds on the reach of next generation dark matter direct detection experiments. Phys. Rev. 2014, 89, 023524. [Google Scholar] [CrossRef]
  153. CTA Consortium. Science with the Cherenkov Telescope Array. arXiv 2017, arXiv:1709.07997. [Google Scholar]
Figure 1. Upper plot: one-dimensional profile likelihood functions for the χ ˜ 1 0 mass in the pMSSM11 with (blue) and without the ( g 2 ) μ constraint (green) and with (solid) and without (dashed) applying the constraints from LHC Run II. Lower plot: triangular presentations of the composition of the χ ˜ 1 0 in the fit with LHC 13-TeV and with the ( g 2 ) μ constraint. The best-fit point is marked with a green star. Taken from [41].
Figure 1. Upper plot: one-dimensional profile likelihood functions for the χ ˜ 1 0 mass in the pMSSM11 with (blue) and without the ( g 2 ) μ constraint (green) and with (solid) and without (dashed) applying the constraints from LHC Run II. Lower plot: triangular presentations of the composition of the χ ˜ 1 0 in the fit with LHC 13-TeV and with the ( g 2 ) μ constraint. The best-fit point is marked with a green star. Taken from [41].
Universe 08 00427 g001
Figure 2. Upper plot: m χ ˜ 1 0 σ p SD plane; for the various lines and shaded areas: see text. Lower plot: m χ ˜ 1 0 σ p SD plane; for the various lines and shaded areas: see text. In both plots the red, blue and green solid lines indicate the 1 , 2 , 3 σ preferred parameter regions, respectively, with the best-fit point shown as a green star. The color coding within the 2 σ area correspond to the DM relic density mechanism. Taken from [41].
Figure 2. Upper plot: m χ ˜ 1 0 σ p SD plane; for the various lines and shaded areas: see text. Lower plot: m χ ˜ 1 0 σ p SD plane; for the various lines and shaded areas: see text. In both plots the red, blue and green solid lines indicate the 1 , 2 , 3 σ preferred parameter regions, respectively, with the best-fit point shown as a green star. The color coding within the 2 σ area correspond to the DM relic density mechanism. Taken from [41].
Universe 08 00427 g002
Figure 3. The results of our parameter scan in the five DM scenarios in the m χ ˜ 1 0 σ p SI plane. The color code indicates the DM relic density. Red points are in full agreement with the Planck measurement.
Figure 3. The results of our parameter scan in the five DM scenarios in the m χ ˜ 1 0 σ p SI plane. The color code indicates the DM relic density. Red points are in full agreement with the Planck measurement.
Universe 08 00427 g003
Figure 4. Cross section predictions at an e + e collider with s = 1000 GeV as a function of the sum of two final state masses. Upper plot: χ ˜ 1 ± coannihilation scenario; lower left plot: l ˜ ± coannihilation case−L; lower right plot: l ˜ ± coannihilation case−R. The color code indicates the final state, open circles are below the anticipated XENON−nT/LZ reach, and full circles are below the neutrino floor.
Figure 4. Cross section predictions at an e + e collider with s = 1000 GeV as a function of the sum of two final state masses. Upper plot: χ ˜ 1 ± coannihilation scenario; lower left plot: l ˜ ± coannihilation case−L; lower right plot: l ˜ ± coannihilation case−R. The color code indicates the final state, open circles are below the anticipated XENON−nT/LZ reach, and full circles are below the neutrino floor.
Universe 08 00427 g004
Figure 5. Tree-level diagram for the two-body decay of a gravitino into a photon and a neutrino, via photino-neutrino mixing.
Figure 5. Tree-level diagram for the two-body decay of a gravitino into a photon and a neutrino, via photino-neutrino mixing.
Universe 08 00427 g005
Figure 6. Diagrams contributing to DM scattering with nucleons mediated by (left panel) SM-like Higgs h and right sneutrinos ν ˜ R via mixing, (middle panel) right-handed sneutrinos ν ˜ R and induced by loop effects involving heavy exotic quarks K , (right panel) Z gauge boson.
Figure 6. Diagrams contributing to DM scattering with nucleons mediated by (left panel) SM-like Higgs h and right sneutrinos ν ˜ R via mixing, (middle panel) right-handed sneutrinos ν ˜ R and induced by loop effects involving heavy exotic quarks K , (right panel) Z gauge boson.
Universe 08 00427 g006
Figure 7. Spin-independent DM-nucleon cross section σ SI versus DM mass m ξ ˜ in the parameter space of our model for the relevant two intervals of right-handed sneutrino VEVs (left panel) 1 TeV < v R < 1 TeV and (right panel) v R > 10 TeV. Blue dots are excluded by the Xenon1T experiment [150]. Green dots will be probed by the upcoming Darwin experiment [151]. Red (black) dots correspond to points above (below) the neutrino floor [152].
Figure 7. Spin-independent DM-nucleon cross section σ SI versus DM mass m ξ ˜ in the parameter space of our model for the relevant two intervals of right-handed sneutrino VEVs (left panel) 1 TeV < v R < 1 TeV and (right panel) v R > 10 TeV. Blue dots are excluded by the Xenon1T experiment [150]. Green dots will be probed by the upcoming Darwin experiment [151]. Red (black) dots correspond to points above (below) the neutrino floor [152].
Universe 08 00427 g007
Table 1. The ranges of the sampled pMSSM11 parameters.
Table 1. The ranges of the sampled pMSSM11 parameters.
ParameterRange
M 1 (−4, 4) TeV
M 2 (0, 4) TeV
M 3 (−4, 4) TeV
m q ˜ (0, 4) TeV
m q ˜ 3 (0, 4) TeV
m ˜ (0, 2) TeV
m τ ˜ (0, 2) TeV
M A (0, 4) TeV
A(−5, 5) TeV
μ (−5, 5) TeV
tan β (1, 60)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Heinemeyer, S.; Muñoz, C. Dark Matter in Supersymmetry. Universe 2022, 8, 427. https://doi.org/10.3390/universe8080427

AMA Style

Heinemeyer S, Muñoz C. Dark Matter in Supersymmetry. Universe. 2022; 8(8):427. https://doi.org/10.3390/universe8080427

Chicago/Turabian Style

Heinemeyer, Sven, and Carlos Muñoz. 2022. "Dark Matter in Supersymmetry" Universe 8, no. 8: 427. https://doi.org/10.3390/universe8080427

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop