Revisiting singlino dark matter of the natural $Z_3$-symmetric NMSSM in the light of LHC

Inspired by the fact that relatively small values of the effective higgsino mass parameter of the $Z_3$-symmetric Next-to-Minimal Supersymmetric Standard Model (NMSSM) could render the scenario `natural', we explore the plausibility of having relatively light neutralinos and charginos (the electroweakinos or the ewinos) in such a scenario with a rather light singlino-like Lightest Supersymmetric Particle (LSP), which is a Dark Matter (DM) candidate, and singlet-dominated scalar excitations. By first confirming the indications in the existing literature that finding simultaneous compliance with results from the Large Hadron Collider (LHC) and those from various DM experiments with such light states is, in general, a difficult ask, we proceed to demonstrate, with the help of a few representative benchmark points, how exactly and to what extent could such a highly motivated `natural' setup with a singlino-like DM candidate still remains plausible.


Introduction
A key ingredient that renders a popular supersymmetry (SUSY) scenario like the phenomenological Minimal or Next-to-Minimal Supersymmetric Standard Model (pMSSM or pNMSSM) 'natural' [1][2][3] is a relatively small SUSY conserving higgsino mass parameter 'µ' in the pMSSM [4][5][6][7] or, similarly, µ eff in the pNMSSM. In both scenarios, this would imply presence of at least two light neutralinos and a similarly light chargino (electroweakinos or ewinos) which are higgsino-like. Though theoretically much motivated, such light ewinos generally derive significant constraints from their null searches at the colliders. These searches target pair or associated productions of such ewinos. A stronger set of bounds emerge in scenarios with significant mass-splits between such states and the lightest neutralino which is the Lightest SUSY Particle (LSP). The LSP is stable when a well-known discrete symmetry called R-parity is conserved and thus, can be a viable candidate for the Dark Matter (DM) [8,9].
In the MSSM, an optimally healthy split between χ 0 2,3 /χ ± 1 and χ 0 1 is not possible when µ M 1 , M 2 for which these states are almost purely higgsinos and hence nearly degenerate, where M 1 and M 2 stand for the soft SUSY-breaking masses of the U (1) and SU (2) gauginos, respectively. However, with M 1 < µ < M 2 , one could find m χ 0 and hence would obtain reasonable mass-splits ∆m (χ 0 2,3 ,χ 0 1 ) and ∆m (χ ± 1 ,χ 0 1 ) leading to hard enough leptons/jets in the cascades. This renders these searches viable and hence yielding constraints. 2 Critical studies as to how strong and robust a constraint the LHC experiments could impose on such relatively light higgsino-like ewinos (and hence on 'µ') have recently been undertaken by various groups [18,19]. Incidentally, the observed SM-like Higgs boson (h SM ) of mass around 125 GeV could at best be the lightest of the MSSM Higgs bosons while all its cousins have to be much heavier (the so-called decoupling limit). Thus, there is only a limited scope for the decay χ 0 2,3 → χ 0 1 h SM to dominate over χ 0 2,3 → χ 0 1 Z and hence weakening of lepton-rich final states is not expected to be very common. Consequently, a notable relaxation on the masses of these light ewinos (and hence on 'µ') is unlikely to be a common occurrence.
In contrast, the situation can get very different in the NMSSM when the coefficient 'κ' of the superpotential term cubic in the singlet chiral superfield gets vanishingly small (the Peccei-Quinn symmetric limit). First, a rather light scalar (h 1 ) and a pseudoscalar (a 1 ) Higgs bosons with m h 1 ,a 1 < m Z , both of which are singlet-dominated, are inevitable [20,21]. The Higgs sector of the NMSSM has been studied in great details in refs. [22][23][24][25][26][27][28][29][30][31][32][33][34][35][36][37]. Second, a rather light singlino-dominated neutralino LSP (mass ranging from sub-GeV to a few tens of a GeV) is naturally present in the spectrum. These two together could easily allow for a much smaller value of µ eff leading to two next-to-LSP neutralino states (χ 0 2,3 ) and the lighter chargino (χ ± 1 ) all of which are higgsino-dominated with masses ∼ µ eff and having prominent decays χ 0 2,3 → χ 0 1 h 1 /a 1 . Also, these facilitate the simultaneous opening up of the decays χ 0 2,3 → χ 0 1 h 2 (h SM ) thus reinforcing the combined branching fractions of χ 0 2,3 to Higgs bosons over the same to Z-boson. Some specific consequences of such possibilities had been studied in the past which include rather light scalars decaying to (i) ττ and leading to soft multilepton final state [38], (ii) bb [39] and (iii) two photons [40,41]. For m a 1 1 GeV, even mesons can be produced out of a boosted pair of light quarks that a 1 /h 1 might decay to. 3 1 There have also been experimental searches involving two soft leptons [12,14], opposite sign di-leptons, as well as final states with b-jets and photons [13], effectively constraining the chargino-neutralino spectra. The implications of these searches for our present study will be discussed in some detail later in this work. 2 For M2 < µ, χ ± 1 becomes wino-dominated and degenerate in mass with a wino-dominated neutralino LSP. This would result in softer leptons/jets in the cascades of χ ± 1 thus eroding experimental sensitivity to multi-lepton/jets final states, in general, and tri-lepton final state, in particular [17], while soft-lepton searches could emerge more relevant [12,14]. 3 Furthermore, as we would appreciate later in this work, one could have a possible situation, without sacrificing much of the essential features of such a scenario, when even the decay χ ± 1 → χ 0 2 W ± * could compete -2 -On the DM front, presence of singlet-like light scalars along with a stable singlinodominated LSP with a critical higgsino admixture (thanks to a not so large µ eff ), would have nontrivial consequences [42]. First, the higgsino admixture could now enable the LSP annihilate efficiently enough in the early Universe yielding DM relic in the right ballpark. Second, the same enhanced interaction of such an LSP could make it sensitive to DM Direct Detection (DMDD) experiments. Third, the light scalars (a 1 and h 1 ) could offer new annihilation 'funnels' that are efficient handles on the DM relic. Furthermore, in the context of DMDD experiments sensitive to spin-independent (SI) scattering, there may appear the socalled blind spots [43][44][45][46] either due to vanishing LSP-Higgs coupling or due to a destructive interference between the contributions from the CP -even Higgs bosons. These could suppress the DMDD-SI cross section to a value still allowed by experiments.
The collider and the DM aspects of such an NMSSM scenario are thus expected to be connected in a rather nontrivial way. It is encouraging to find a few recent works addressing these aspects, focussing mainly on one or the other of them. Their broad scopes are as follows.
• Ref. [15] is the first one to discuss the case of a light singlino-like LSP as the DM candidate with light bino (higgsino)-like neutralino(s) and a higgsino-like chargino as the next heavier sparticle(s) and the combined constraint such a scenario draws from various DM and collider experiments. It points out the roles played by relatively light a 1 /h 1 (i) in obtaining the DM Relic Density (DMRD) in the right ballpark, (ii) in complying with the DMDD constraints on the SI scattering cross section using the blind spot mechanism and (iii) in evading (degrading) the LHC bounds on such light ewinos.
• Ref. [46] undertakes a detailed scan of the 'natural' NMSSM parameter space requiring relatively light higgsino-like states compatible with the relic density bound from Planck experiment [47,48], the bound from DMDD experiments like XENON-1T [49,50] and those from the 13 TeV run (with up to 36 fb −1 of data) of the LHC. In our current context, the most relevant finding is that only a singlino-dominated LSP with a small higgsino admixture (µ eff m χ 0 1 ) might survive the combined constraints if m χ 0 1 90 GeV and, that also, for a compressed spectrum for the LSP and χ ± 1 .
• Ref. [16] is mainly concerned with the impact of recent multi-lepton searches at the LHC on the ewinos of the NMSSM in the presence of light singlet scalars, h 1 and a 1 .
The study chooses to remain agnostic about the detailed bounds in the DM sector except for respecting only the upper bound on the relic density. The discussion on the scenario with a relatively light singlino-like LSP that could annihilate via singlet-Higgs funnel(s) are of particular relevance for our present work.
The former would add to jet activity via the decay χ 0 2 → χ 0 1 h1/a1 and hence could potentially alter bounds obtained from the studies which vetoes extra jets. Otherwise, BR(χ ± 1 → χ 0 1 W ± * ) would remain 100% and hence collider constraints derived solely by studying χ ± 1 -pair production would hold in a robust manner.
-3 -In this work we focus on an Z 3 -symmetric NMSSM scenario with a relatively small µ eff (preferably less than ∼ 300 GeV and not exceeding 500 GeV) that ensures enhanced 'naturalness' and with a singlino-enriched (> 95%) LSP neutralino as the DM candidate with mass around or below the SM Higgs boson funnel, i.e. 62 GeV. The purpose is to find how such a scenario could still be compatible with all pertinent experimental data from both DM and collider fronts. Our study goes beyond what was found in ref. [46] which excludes the possibilities of having a singlino-dominated LSP below ∼ 90 GeV and away from the coannihilation regime. As we would elucidate soon, allowing for some modest bino content in the lighter neutralinos by considering an appropriately small M 1 could provide us with a much lighter and a viable singlino-dominated DM candidate which finds right funnels in various light states like the SM Higgs boson, the Z-boson and even the lighter singlet-like Higgs states of the scenario to. This renders, not only the DM neutralino, but the entire system of lighter neutralinos 'well-tempered' [51]. Constraints imposed by us include the one on DMRD within 10% uncertainty, those from the DMDD experiments like XENON-1T studying the SI [49] and the spin-dependent (SD) [50,52] DM-nucleon scattering cross sections. Our study also takes into account all relevant LHC analyzes that considers up to ∼ 36 fb −1 worth data via use of the package CheckMATE [53,54].
The paper is organized as follows. In section 2 we present the structures and the salient features of the NMSSM ewino and the Higgs sectors along with their interactions that are relevant to the present work. These are followed by a discussion on the nature of the spectrum in our scenario and on the important decay modes of the light ewinos to Higgs bosons. Section 3 contains our results where the impact of experimental bounds on the DM observables is quantitatively assessed leading to our choice of suitable benchmark points with low enough µ eff . A dedicated CheckMATE-based analysis follows to assess the viability of the benchmark points in view of the LHC data. In section 4 we conclude.

The light ewinos and the light Higgs bosons
The superpotential of the Z 3 -symmetric NMSSM is given by where W MSSM | µ = 0 is the MSSM superpotential sans the higgsino mass term (the µ-term), H u , H d and S are the usual MSSM SU (2) Higgs doublets and the NMSSM-specific singlet superfields, respectively while 'λ' and 'κ' are dimensionless coupling constants. The µ-term is generated dynamically from the second term when the singlet scalar field 'S' develops a vacuum expectation value (vev) S =v S (i.e., µ eff = λv S ) thus offering a solution to the puzzling µ-problem [55]. The NMSSM-specific part of the soft SUSY-breaking Lagrangian is given by where m 2 S is the squared soft SUSY-breaking mass for the singlet scalar field 'S' while A λ and A κ are the NMSSM-specific trilinear soft terms having dimensions of mass. In the following -4 -subsections, we briefly discuss the sectors that are directly involved in the present study, i.e., the ewino and the Higgs sectors of the scenario.

The ewino sector
The ewino sector is comprised of the neutralino and the chargino sectors. The neutralino sector is augmented in the NMSSM by the presence of the singlino (S) state, when compared to the same for the MSSM. Thus, the symmetric 5 × 5 neutralino mass matrix, in the basis where g 1 and g 2 stand for the gauge couplings of the U (1) and SU (2) gauge groups, respec- GeV. The above mass-matrix can be diagonalized by a matrix N , i.e., The resulting neutralino mass-eigenstates (χ 0 i , in order of increasing mass as 'i' varies from 1 to 5), in terms of the weak eigenstates (ψ 0 j , with j = 1, . . . , 5), is given by It is possible to find analytic expressions for the masses and the elements of the mixing matrix, N ij , when two of the five states get decoupled. Hence, to start with, for our purposes, we consider the bino and the wino states to be decoupled. This would describe our basic setup fairly robustly with a rather light singlino-like LSP and a relatively small µ eff (thus aiding 'naturalness') leading to two light higgsino-like states. Such a scenario can be realized for λv |µ eff | along with κ/λ 1. The ratios of higgsino to singlino admixtures in a given neutralino (in particular, in the LSP) would remain to be much instrumental in our present analysis. In the above-mentioned situation, these are given by [43,56] where N i3 , N i4 and N i5 denote the two higgsino and the singlino components, respectively, in the i-th mass eigenstate with i = 1, 2, 3 and m χ 0 1 < m χ 0 2 < m χ 0 3 . Subsequently, we note that a relatively small value of M 1 ( µ eff ) could have a nontrivial impact on the combined DM and collider phenomenology of such a scenario with light ewinos. However, this compels one to work with a 4 × 4 neutralino mass-matrix for which analytical expressions for N ij , similar to those in eq. (2.6), would not be much illuminating. On top of that, when M 2 is allowed to become small, the eigenvalue problem seeks solution of a polynomial of degree 5 of which a general solution does not exist. Hence, for smaller values of M 1 and/or M 2 , we adopt a numerical approach. On the other hand, the 2 × 2 chargino mass matrix of the NMSSM is structurally the same as that of the MSSM with µ → µ eff and, in the basis is given by [22] As in the MSSM, this can be diagonalized by two 2 × 2 unitary matrices U and V : (2.9) As noted in the Introduction, to ensure our scenario remains reasonably 'natural', we choose to work with relatively low values of µ eff . This yields two light neutralinos along with a lighter chargino with masses ∼ µ eff , all of which can be dominantly higgsino-like. However, their actual masses and compositions depend much on the extent they mix with the singlino and the bino (for the neutralinos only) and with the wino states. In particular, we are interested in a scenario where, 2κv S µ eff (i.e., for κ λ/2). This could lead to a singlino-dominated LSP. However, it may contain a crucial higgsino admixture thus making it a viable DM candidate. Implications of such an LSP in the context of various DM and collider experiments have recently been studied in the literature [15,16,46], though as parts of more general studies. Apart from the subtle role played by our proposed manoeuvring by allowing for M 1 µ eff , this brings in a fourth relatively light neutralino in the picture. We will further assume the wino-like neutralino to be the heaviest of them all and hence would require M 2 > µ eff , M 1 . This would help avoid stringent collider constraints by restricting heavier ewinos cascading via such wino-like states. In the next subsection, we discuss that such a scenario is necessarily accompanied by light singlet-like scalars which characterize our scenario of interest.

The Higgs sector
The superpotential of eq. (2.1) leads to the following Lagrangian containing soft masses and couplings for the NMSSM Higgs sector: The neutral Higgs fields are parameterized about the real vev's v d , v u and v S for the three neutral fields H 0 d , H 0 u and S, respectively as where "R" and "I" denote, for each field, the CP -even and the CP -odd states, respectively. The CP -even squared mass matrix, M 2 S , in the basis {H dR , H uR , S R }, is given by [22] where g 2 = (g 2 1 + g 2 2 )/2. The squared mass of the singlet-like CP -even eigenstate (up to a mixing with the doublet states) is given by the (3,3) component, i.e., Out of the other two eigenstates, one has to turn out to be the SM-like Higgs boson with mass ∼ 125 GeV, the other one being a relatively heavy, doublet-dominated neutral Higgs boson with its squared mass around µ eff (A λ + κv S )/ sin 2β. Thus, a more realistic basis to work in is with m 2 A = 2µ eff (A λ + κv S )/ sin 2β representing the squared mass of the doublet-like CPodd scalar, as in the MSSM. The mass-squared for the singlet CP -odd scalar (modulo some mixing) is given by the (2,2) element of the above matrix, i.e., The mass eigenstates of the the CP -even (h i ) and the CP -odd (a i ) sectors are given by [43,45] where E (3 × 3) and O (2 × 2) are the matrices that diagonalize the mass-squared matrices for the CP -even scalars in the basis {H 1 , H 2 , S R } and that for the CP -odd scalar of eq. (2.14). Clearly, the scalar masses have rather complex dependencies on as many as six input parameters like λ, κ, A λ , A κ , µ eff and tan β. However, for our scenario of interest with a light singlino-like LSP (m χ 0 1 ≈ mS ∼ 2κv S ) and light singlet-like scalars (given by eqs. (2.13) and (2.15)), one could find the following (approximate) sum-rule [15,57] relating their masses when the singlet-doublet mixing among the scalar (Higgs) states can be safely ignored, i.e., in the decoupling limit (λ, κ → 0) or for a sizable tan β and not too large λ, κ and A λ : This clearly indicates that the masses of the singlino and those for the singlet-like scalar and the pseudoscalar are rather closely tied. The relationship becomes handy in discussions on DM-annihilation via light scalar funnels [15,46].

Interactions among the ewinos and the scalars
Interactions among the ewinos and the Higgs-like scalars states take the central stage in our present study. Their subtle dependence on various NMSSM parameters and their interplay crucially shape the phenomenology on both DM and collider fronts, sometimes in a rather complementary fashion.
To be a little more specific, conformity with the observed value of DMRD would depend not only on a mass-spectrum that offers efficient DM-annihilation mechanisms via funnels and/or coannihilations 4 but also on the strengths of the involved interactions. The latter, in turn, could also control the DM-nucleon interactions that are studied at the DMDD experiments. Hence requiring an efficient DM-annihilation to meet the DMRD observations might imply a strong enough DM-nucleon interaction strength that is ruled out by the DMDD experiments. The converse is also true. This highlights a built-in tension in finding a simultaneous explanation of the two crucial observations in the DM sector alone.
On the collider front, the interactions among the ewinos and the scalars determine the branching fractions of the former to the latter. Such modes include the ones beyond what are being routinely considered in the LHC analyzes in the context of 'simplified scenarios' and result in new final states. These are likely to result in relaxed mass-bounds on the higgsinolike states thus offering enhanced 'naturalness'. Interestingly enough, as we will discuss soon in section 3, these might also help satisfy the constraints in the DM sector. Thus, for decoupled sfermions and a gluino, the interactions that are of paramount importance are those among (i) various neutralinos and the gauge (Z-) boson and (ii) various neutralinos and Higgs bosons, of both CP -even (scalar) and CP -odd (pseudoscalar) types, from both doublet and the singlet sectors.
The neutralino DM interacts with the Z-boson only through its higgsino admixture. This interaction governs the self-annihilation of DM via Z-boson funnel thus controlling the DMRD as well as the DMDD-SD cross section and is given by α Zχ 0 . On the other hand, a doublet-like Higgs scalar has an MSSM-like interaction with a higgsino and a gaugino. In addition, in the Z 3 -symmetric NMSSM, as can be gleaned from eq. (2.1), this also interacts with a higgsino and a singlino while the singlet-like scalar interacts with two higgsinos, both strengths being proportional to 'λ'. The S 3 term in eq. (2.1) further implies that the singlet scalar has an interaction with two singlinos whose strength goes as 'κ'. Thus, if the gaugino (bino and/or wino) admixture in a singlino-dominated LSP can be ignored (which is somewhat ensured by the neutralino mass matrix of the NMSSM), the generic coupling of such an LSP with the CP -even Higgs scalars are given by [43] α h i χ 0 The couplings α a i χ 0 1 χ 0 1 for the CP -odd scalar counterparts would be somewhat similar except for the appearance of an overall factor of imaginary 'i' and that in the rotated basis for the pseudoscalar sector there are only two massive eigenstates.
Furthermore, while the CP -even Higgs states from the doublet and the singlet sectors contribute to both DMRD and DMDD-SI, their CP -odd counterparts could contribute only to DMRD and practically nothing to any DMDD processes [59]. Ref. [43] discusses the issue of the blind spots for DM-nucleon interaction in a few specific and motivated scenarios in the Z 3 -symmetric NMSSM. Among these, the scenario that is germane to our present study is the one (section 6) that discusses blind spots arising from destructive interference between CP -even singlet-like (h 1 ) and doublet-like (h 2 , the observed SM-like Higgs boson) scalar states where m h 1 < m h 2 , while the heavier MSSM-like CP -even Higgs state is virtually decoupled. This yields a DM-nucleon SI cross section well below the threshold of sensitivity of the relevant DMDD-SI experiments.
Up to this point, the relative strengths of all the couplings that matter are essentially governed the ratios presented in eq. (2.6). It is now instructive to note that if the singlinodominated LSP could be infused with a bino/wino component, it would alter the higgsino shares in the same. This can be achieved by allowing bino/wino to mix substantially with higgsinos, given that, at the lowest order, this only can (indirectly) induce some gaugino admixture in an otherwise singlino-dominated LSP. Hence such a regime would reign as long as M 1 (or M 2 , though decreasing it beyond a point could attract severe experimental constraints) is not too far away from µ eff . In certain regions of the NMSSM parameter space, with mS < M 1 < µ eff , this causes the coupling-strength α Zχ 0 weakening to a minimum due to rather involved variations of N ij 's as functions of M 1 . This we will discuss soon in a little more detail. This would then diminish the DMDD-SD cross section thus helping us evade the related experimental bound. Clearly, under such circumstances, eq. (2.6) ceases to hold and improving the same in the presence of an active bino state is unlikely to be illuminating enough, given the complicated structure the situation presents. We thus take a numerical route for the rest of the present study and frequently confront the results with broad-based expectations for checking their basic sanity.
It may further be noted that the higgsino content of the LSP (given by N 2 13 + N 2 14 ) could contribute only partially to the DMDD-SI cross section (for the DM-nucleon scattering process mediated by the doublet CP -even Higgs bosons) while there could be a significant additional contribution from the singlet-like Higgs exchange in such a scattering. However, in the region of parameter space of our interest for which κ ∼ O(10 −2 ), this contribution is expected to be suppressed. An increase in the total higgsino fraction could attract severe experimental constraints from the DMDD-SI experiments. However, its effect may get subdued in the presence of blind spots in the SI processes. In this work we exploit these two simultaneous effects in our favour, by manoeuvring M 1 and/or M 2 , to find compliance with the DMDD data while still obtaining a thermal relic density within the Planck-allowed range.
In the upper panel of figure 1 we present the variations of the quantity |N     left (upper right) plot. One clearly finds that the magnitude of |N 2 13 − N 2 14 | could practically drop to a vanishing level as M 1 decreases. However, as can be gleaned from the plots in the upper panel, for what exact value of M 1 this happens, depends on the input value of M 2 , although it becomes more or less insensitive to M 2 for its larger values. These two plots also reveal that such a phenomenon occurs only for M 1 and 'κ' having no relative sign between them, a situation in which the mixing between the two involved sectors is known to get maximal. It is worth pointing out that even though the relevant null entry in the neutralino -10 -mass matrix (eq. (2.3)) prohibits a direct mixing between the bino and the singlino states, a possible mixing through the higgsino portal could give rise to something that drastic with important phenomenological consequence akin to a blind spot for DMDD-SI scattering, but this time occurring for DMDD-SD scattering. We exploit this effect in our study the results of which are presented in section 3.The onset of discontinuous flat line segments seen at the top left (right) part of the plots on left (right) has its origin in the bino-like LSP with a negative mass-eigenvalue turning instantly to a singlino-dominated one with a positive eigenvalue, for certain particular values of M 1 depending on values of other input parameters.
Plots in the lower panel of figure 1  It worths a mention that the crucial variation is the one that of |N 13 |. While it may not be outright unexpected that lowering of M 1 would immediately result in an enhanced bino admixture in the LSP, at the expense of mostly a decreasing singlino fraction in the same, it is somewhat curious to note that a decreasing M 1 boosts the otherwise subdominant higgsino content of the LSP in the form of |N 13 |. It is possible that a decreasing M 1 , given its healthy connection to the higgsino sector, drags the higgsino along on a collective bid to deplete the singlino content in the LSP. The discontinuity of the curves appearing for certain negative M 1 values in the upper left plot are also efficiently explained by the plots in the lower panels.

The spectrum and the decays
Discussions in the previous susbsections reveal that both the light (singlet-like) Higgs sector and the neutralino sector get simultaneously affected in a rather intricate way as 'κ' turns smaller. This includes non-trivial modifications of the involved couplings among these states via mixings effects in both sectors and resulting mass-splits between the physical states. Together these could alter the phenomenology in an essential manner and experimental analyzes need to take due note of the same.
As has been already pointed out, in the scenario under study, the lightest neutralino (the LSP) is singlino-like whereas the immediately heavier neutralinos, to start with, are higgsino-like. The latter could have enhanced decay branching fractions to singlet-like Higgs bosons, h 1 and a 1 , which can become light enough for suitably small values of 'κ'. Under such a circumstance, the SM-like Higgs boson is the second lightest CP -even Higgs boson (h 2 ∼ h SM ) and this is always the case in our present study. As mentioned in the Introduction, the decay branching fractions of the neutralinos to lighter (singlet-like) Higgs bosons could then compete with (or could even exceed) those for the popularly considered modes like χ 0 2,3 → χ 0 1 Z ( * ) /h 2 (h SM ) and this is likely to relax the existing bounds on the ewino sector.
-11 -A further nontrivial alteration of the decay branching fractions of the higgsino-like neutralinos may take place if one allows for a bino/wino-like neutralino (now χ 0 2 ) sneak in below the formerly higgsino-like states (now χ 0 3,4 ). Thus, more involved cascades could kick in, viz., χ 0 3,4 → χ 0 2 (→ χ 0 1 h i )Z/h i (with i = 1, 2, 3 standing for the two light singlet-like and the SM-like Higgs bosons) thanks to some higgsino admixture in an otherwise gaugino (bino)-dominated χ 0 2 . This would have important bearing on collider phenomenology. In section 3, we shall discuss how such an intermediate state plays a crucial role in finding an all-round compliance with the experimental results pertaining to the DM-sector (as pointed out in section 2.3) as those from the colliders.
In passing, it is to be noted that presence of light Higgs states would not directly affect the decay of the lighter chargino for which the experiments assume BR(χ ± 1 → χ 0 1 W ±( * ) ) to be 100% when the other Higgs states of the NMSSM, along with the sfermions, are all much heavier. Thus, at the first sight, it might appear that bounds imposed on the lighter chargino sector, in particular, by looking for its pair-production, and, consequently, on µ eff (for a higgsino-like lighter chargino) would still hold and need to be respected. However, there are a couple of caveats. First, since the presence of a light singlino state could significantly modify the NMSSM neutralino spectrum through its mixing with the light higgsino states, a reasonable mass-split between χ 0 2,3 and χ ± 1 cannot be ruled out. This could open up competing decay modes of χ ± 1 in the form χ ± 1 → χ 0 2,3 W ±( * ) . While these would still lead to final states with leptons thanks to the presence of W ±( * ) , the same are likely to be contaminated with the decay products of χ 0 2,3 , as noted in the last paragraph. Second, as discussed above in the case for the neutralinos, the competing decay mode in the form of χ ± 1 → χ 0 2 W ±( * ) could again open up for the lighter chargino when we require, as discussed in section 2.3, M 1 to be brought down below µ eff . In both cases, experimental bounds even from the study of chargino-pair production would likely to get relaxed. 5 As for the light Higgs states (h 1 , a 1 ) appearing in the cascades of the lighter neutralinos, those could have significant branching fractions to bb similar to the case of the SM-like Higgs boson. However, in general, constraints derived from neutralino cascades involving such Higgs states are weaker when compared to those obtained with cascades involving Z ( * ) [13]. Thus, enhanced branching fraction for the decay χ 0 2,3 → χ 0 1 h 1 /a 1 (at the expense of BR(χ 0 2,3 → χ 0 1 Z ( * ) )) are expected to relax the existing experimental bounds on the ewino sector thus capable of opening up a more 'natural' region of the NMSSM parameter space.
In this work we confine ourselves to a region of the Z 3 -symmetric NMSSM parameter space for which the LSP is a singlino-dominated (> 95%), the lighter chargino and two neutralinos are higgsino-like with masses 300 GeV, with a further possibility of having an intermediate (gaugino-like) neutralino lighter than the higgsino-like states. In addition, the setup offers singlet-like scalars that are lighter than the SM-like Higgs bosons which could 5 Note, however, that if M2 µ eff , this would present us with a lighter chargino which is wino-like and close in mass with χ 0 2 . Hence the second effect mentioned above would be absent and BR(χ ± 1 → χ 0 1 W ± ) would be 100%. This would thus invite the standard, stronger bound on M2 from null searches for chargino pair-production at the LHC.
-12 - M 2 (TeV) 0.05-0.2 0.001-0.05 1-60 ≤ 300 ≤ 10 ≤ 100 50-500 0.2-1 Table 1. Ranges of various model parameters adopted for scanning the Z 3 -symmetric NMSSM parameter space. All parameters are defined at the scale Q 2 = (2m 2 Q + m 2 U + m 2 D )/4, except for tan β which is defined at m Z (see text for details). even turn out to be lighter than the LSP. Phenomenological possibilities discussed in the previous paragraphs are realized in such a set up. Scan-ranges adopted for various model parameters are summarized in table 1. The soft masses for the SU (3) gaugino (M 3 ), those for the sfermions and the soft trilinear parameters A τ,b,t are all fixed at around 5 TeV while A e,µ is set to zero.

Results
We now present our results for the broad scenario discussed in the previous section which is characterized by a light singlino LSP, accompanied by rather light singlet-like scalars, along with higgsino-dominated lighter chargino and neutralinos that ensure a healthy degree of 'naturalness'. The focus is on if such a scenario can be compatible with recent constraints pertaining to the DM sector (i.e., those involving DMRD, DMDD-SI and DMDD-SD) and those coming from various past and recent collider experiments that include the LEP and the LHC experiments. In particular, it emerges from the recent literature [16,46] that such an allround compliance is not easy to find. As pointed out in the Introduction, ref. [46] concludes that this may be only possible in the coannihilation region marked by a near-degeneracy of the singlino-dominated LSP and the higgsino-dominated chargino (and neutralinos). Our goal is to go beyond this and to find if such a thorough compliance with DM and collider data is possible away from the coannihilation region while still retaining the essential features of the broad scenario.
Results are obtained via a random scan over the parameter space of the Z 3 -symmetric NMSSM using the package NMSSMTools-v5.1.0 [62][63][64]. Experimental constraints (at 2σ level) implemented in NMSSMTools are automatically imposed on our analysis. These include various constraints from the LEP experiments and those on the B-physics observables. Compliance with experimental results on (g − 2) µ is not demanded. In addition, constraints from various Higgs boson searches at LEP and Tevatron and compatibility to the Higgs boson observed at the LHC are considered/checked by using the packages HiggsBounds-v4.3.1 [65,66] and HiggsSignals-v1.4.0 [67,68]. DM-related computations are done using an adapted version of the package micrOMEGAs-v4.3 [69][70][71] that is built-in to NMSSMTools. Finally, we employ the package CheckMATE-v2.0.26 [53,54] to check our benchmark points (that pass all relevant constraints including the DM-related ones) if they are passing all relevant LHC analyzes.

Impact of bounds from the DM sector
Unless otherwise stated, in the present work, bounds from the DM sector would imply strict adherence to a relic density within 10% of the central value of Ωh 2 = 0.119 measured by the Planck experiment [47,48], i.e., 0.107 < Ωh 2 < 0.131. Allowed maximum values for the DMnucleon scattering cross sections are taken (somewhat conservatively, for DM-mass 30 GeV, for which the DMDD-SI bound is the strongest) to be σ SI - 14 -In figure 2 we illustrate various relevant aspects of the regions of Z 3 -symmetric NMSSM parameter space that are simultaneously compatible with all experimental data pertaining to DMRD, DMDD-SI and DMDD-SD. These aspects are as follows.
• The plot on the top is in the m χ 0 1 −m χ ± 1 plane with the value of the combined branching fraction of χ 0 2,3,4 in the decay mode χ 0 2,3,4 → χ 0 1 Z indicated in the palette which could reach a possible maximum value of '3'. Visibly, over the dark patch along the diagonal, the singlino-dominated DM neutralino is nearly mass-degenerate with the lighter chargino (χ ± 1 ) and the next two lighter neutralinos (χ 0 2,3 ) all of which are higgsino-like. Hence coannihilation of the DM neutralino with these states is rather efficient. This renders DMRD in the right experimental ballpark. 6 It is also important to note that due to this degeneracy, the bounds on the ewino sector are also much relaxed over this region [14]. Hence parameter points from this region have a good chance to survive bounds obtained from both DM experiments and the LHC. In fact, ref. [46] pointed out this to be the only region for a singlino-dominated LSP which could exhibit such a simultaneous compliance with data. Note that given M 1 < µ eff is a possibility in our scan, there may be a situation when χ 0 2 becomes bino-dominated while χ 0 3,4 become higssino-like. For such a spectrum, the decay χ 0 2 → χ 0 1 Z may be kinematically disfavoured while χ 0 3,4 → χ 0 1 Z could open up and become relevant. In agreement with ref. [46], our scan also finds strips of DM-allowed points at LSP masses with the SM Higgs and Z-boson funnels, i.e., for m χ 0 1 = m h 2 /2 and at m χ 0 1 = m Z /2, respectively. However, it appears that these strips extend to much higher values of m χ ± 1 (apparently limited only by our choice of the upper limit of µ eff ( 300 GeV)) when compared to what was found in ref. [46]. Also, unlike in ref. [46], the bottom sections of the funnel strips for the SM Higgs boson and the Z-boson are found to be notably populated. We indeed notice that compliance with DMDD-SD data is facilitated with low values of M 1 , as discussed in section 2.3. In this region there is a substantial mass difference between χ ± 1 and χ 0 1 due to which hard enough leptons are expected from decays of χ ± 1 . Furthermore, the 3-body decays χ 0 2 , χ 0 3 → ¯ χ 0 1 (presumably via an off-shell Z boson) could contribute significantly. Thus, this region is expected to get severely constrained from tri-lepton searches at the LHC [12,13,72]. One could as well expect a corresponding 3-body decay of χ 0 2,3 that involves a bottom quark pair. Hence searches involving b-jets in the final states [13] are likely to get sensitive to the said region of parameter space.
Furthermore, we find a DM-allowed region with lighter LSP masses ( 20 GeV) possessing funnels in light singlet scalars (mostly a 1 , and only occasionally, h 1 ). Refs. [15,42,46] had correctly argued on the difficulty in realizing an a 1 funnel. However, as envisaged in ref. [15], we now find a generic region with a singlino-dominated LSP with mass 20 GeV that possesses a 1 funnel for even (an optimally) small 'λ' along with rather large A λ . As predicted, the region indeed yields a rather light h 1 which, ref. [46] argued, would yield too large a DMDD scattering rate to survive the experimental data. Here, it is our specific observation that a suitably low value of M 1 could again do the trick by pushing the DMDD rate down to a safe level. 7 A closer inspection reveals that funnels at work for a singlino-like LSP with mass 40 GeV can be that of h 1 or a 1 or both, simultaneously. In addition, emergence of points only in discrete strips, even though the LSP and the light scalar masses are varying, are due to stringent requirement of having the relic density within a specific band about its observed central value.
Of some interest are the points in darker shades in the funnel strips. These are the points for which the collective branching fractions in the decay modes of χ 0 2,3,4 containing a real Z-boson are tiny. Thus, it may be expected that these could evade some pertinent collider bounds while being still consistent with all DM data, unless m χ ± 1 is too small, as is the case at the bottom these strips. This is since the latter kinematically prohibits the decay(s) of one or more of the participating heavier neutralinos (χ 0 2,3,4 ) to χ 0 1 Z. Nevertheless, three-body decays (via an off-shell Z-boson) into leptonic final states may remain significant, as discussed before. In addition, we find regular (sparse) population of darker points within the strips representing h SM (Z-boson and h 1 /a 1 ) funnel(s) for higher values of m χ ± 1 as well. These result from opening up of new decay modes involving lighter Higgs bosons for χ 0 2,3,4 due to genuine (dynamical) suppressions of the strengths for the χ 0 2,3,4 χ 0 1 Z interaction in the presence of competing χ 0 2,3,4 χ 0 1 h i interactions. Clearly these points need to be subjected to thorough examination to ascertain their viability against LHC data. We undertake this exercise, for relevant final states involving leptons mostly, using CheckMATE in section 3.3 with reference to a few benchmark points picked from all the three funnel regions.
• The plots in the bottom row of figure 2 convey the interplay of M 1 and µ eff keeping m χ 0 1 and the combined branching fraction BR(χ 0 2,3,4 → χ 0 1 Z) in reference. Thus, while the left plot reveals the funnel strips over specific m χ 0 1 ranges (thus corresponding exactly to the plot on the top) having the branching fraction to Z-boson either less (indicated by '+' marks) or greater (indicated by circular blobs) than 1.5, the right plot explicitly displays the same branching fraction with the three specific (funnel) ranges for the associated m χ 0 1 being indicated by three different symbols: ' ' for the SM Higgs funnel, ' ' for the Z-boson funnel and '•' for the singlet-like scalar(s) funnel. These two plots clearly reveal that to achieve a dominant (≥ 1.5) combined branching fraction to every 7 Some such situations are discussed in ref. [16] as specific benchmark points. However, given that the work focuses on the impacts of the LHC data, it remains agnostic as to whether such points would satisfy various DM-related constraints but for the DMRD upper bound. We observe that most of these points possess a rather light h1 (m h 1 ( 20 GeV)) which would make it difficult to survive DMDD bounds unless for suitable M1 < µ eff thus yielding a bino-like χ 0 2 .
-16 - other mode save χ 0 1 Z (thereby evading relevant collider bounds; represented by the '+' symbol) one requires M 1 < µ eff . It is somewhat curious to note that the combined branching fraction to Z-boson could systematically go down to a value of ≈ 1 but does not drop further if M 1 > 250 GeV.
In figure 3 we illustrate how the values of 'λ' (left) and hence N 2 15 (the singlino admixture in the LSP) are distributed in the σ SI From the left plot one can clearly see that low values of 'λ' ( 0.2) are preferred. This is not unexpected for the following reasons.
• First, The DMDD-SI cross section dominantly involves coupling of a DM(LSP)-pair to the singlet-like Higgs bosons (h 1 χ 0 1 χ 0 1 ) which is enhanced for a mixed singlino-higgsino LSP. Given the higgsino admixture in an otherwise singlino-dominated LSP is proportional to 'λ' (for given fixed values of m χ 0 1 and µ eff ), the coupling in context grows with its value and could lead to a large enough SI cross section that is ruled out by the experiments.
• Second, the DMDD-SD cross section, in contrast, involves coupling of a DM(LSP)-pair to a Z-boson. This, on the other hand, depends on the higgsino content of an otherwise singlino-dominated LSP and hence grows as λ 2 (for given fixed values of m χ 0 1 and µ eff ). This could result in a large enough SD cross section which again could be ruled out by relevant experiments.
-17 - The plot on the right first corroborates the correlation between 'λ' and N 2 15 that is explained above, i.e., the smaller is the value of 'λ', the smaller (larger) is the higgsino (singlino) admixture in the singlino-dominated LSP. Furthermore, one finds that the reddish/purple part on the right edge of the plot has an enhanced higgsino fraction in the LSP and hence always gets ruled out by DMDD-SI data. However, DMDD-SD data may still allow such 'λ' values (in the bottom right quadrant) which is due to somewhat smaller sensitivity of SD rates to 'λ', as hinted above. In contrast, the black regions are very special in the sense that these have the LSP which is bino-dominated (when M 1 goes below mS ∼ 2κv S in our scan). The only admixture that is pertinent here is in the form of higgsinos (since bino does not mix directly to singlino at the lowest order) and 'λ' is likely to decouple from DM physics. Thus, a small higgsino admixture in the LSP could suffice to result in a large enough SI and SD cross sections that are ruled out by experiments. Nonetheless, we find a tiny bino-dominated region in the intersection of the two boundaries that separate the DD-allowed regions. For clarity, it may be mentioned that the points that appear in the allowed (bottom left) quadrant comply with all DM data and hence are the same data-points that show up in figure 2.

Benchmark scenarios
In this subsection we briefly discuss our strategy to choose a few representative benchmark points that worth thorough scrutiny against recent LHC data in order to establish their viability. We choose our benchmark points from the scan described earlier by ensuring that these all have a singlino-dominated (> 95%) LSP, have low values of µ eff and satisfy basic experimental constraints mentioned earlier including those from the DM-sector. The scenarios are divided into three categories according to the DM-annihilation funnels at work, i.e., singlet (pseudo)scalar funnel, Z-boson funnel and SM-like Higgs funnel. Note that we have ensured, apart from satisfying the DMRD and the DMDD constraints, our benchmark points also satisfy various other constraints from indirect DM searches [73][74][75] thanks to a small annihilation cross-section at late times ( σv O(10 −29 ) cm 3 s −1 ). Next, we look for if the combined decay branching fraction of the heavier neutralinos to χ 0 1 Z could be on the smaller side so that such points stand higher chance of evading LHC constraints on the lighter ewino sector. Furthermore, we try to ensure that the decay branching fraction for χ ± 1 → χ 0 2 W ± competes or even exceeds that for χ ± 1 → χ 0 1 W ± adopted in the standard paradigm for experimental analyzes. This would further relax the existing bounds in this sector.
In table 2 we present these benchmark points by indicating the relevant input parameters, the resulting spectra, the contents of the LSP and the next-to-lightest neutralinos, various relevant branching fractions along with the values for the DM observables. Finally, we summarize for each of these points, their status in view of recent LHC analyzes obtained via CheckMATE.

Impact of recent LHC results: a CheckMATE -based analysis
In this section we describe the status of the benchmark scenarios presented in table 2 in the light of the LHC results. As can be seen, these scenarios feature a light higgsino-like chargino,  It may be reiterated that when we consider only singlino-and higgsino-like light neutralinos in the presence of a light (pseudo-) scalar Higgs in the spectrum, the following decay channels are of importance: where h ≡ {h 1 , h 2 (h SM )} ('a') represents a CP -even (CP -odd) Higgs boson. In presence of a bino-like state in the spectrum, typically χ 0 2 for our benchmark points, the following additional decay modes can be relevant too.
Depending on the mass-difference between the heavier higgsino-like states and χ 0 1 , on-or off-shell gauge/scalar bosons may appear in the above decays of the light ewinos. Since we mainly focus on the uncompressed region, with rather sizable mass-split between the heavier higgsino-like states and χ 0 1 , on-shell gauge bosons feature in all our benchmark scenarios. Considering final states with leptons, the following final states are going to be relevant.
• Chargino pair production (pp → χ ± 1 χ ∓ 1 ) can lead to 2 + / E T (missing E T or MET). In the presence of a bino-like χ 0 2 , there could be significant number of events with up to four accompanying b-jets, assuming the Higgs boson in the cascade dominantly decays into two b-quarks.
• Chargino-heavier neutralinos associated production (pp → χ ± 1 χ 0 2,3,4 ) can lead to 3 + / E T and + 2b + / E T . As in the previous case, the presence of a bino-like χ 0 2 , either produced in the hard scattering or in the cascade of heavier neutralinos, might lead to final states with an enhanced b-jet multiplicity.
• Finally, heavier neutralino-pair production could lead to up to 2 +2-jets/4 + / E T where, in our case, the pairs of leptons come from the decay of on-shell Z-bosons. The presence of a bino-like χ 0 2 in the cascade, as before, would ensure enhanced b-jet multiplicity in the final state.
We use CheckMATE-v2.0.26 to test our benchmark scenarios against relevant experimental analyzes by the ATLAS and the CMS collaborations (which are already implemented in CheckMATE and have been validated) at the 13 TeV LHC with up to 36 fb −1 worth data. The mono-jet/γ + / E T [76,77] are relevant for pair production of the lightest neutralino, together with an ISR jet or a photon. Searches for two soft leptons [14,78] can be relevant for compressed spectra of light ewinos. While searches in these final states could, in general, put reasonable constraints on the chargino-neutralino spectra, these are not expected to be much constraining in the current context. In the case of mono-jet/mono-photon searches, the insensitivity stems from the small production cross-section of χ 0 1 -pair in the present scenario. Soft lepton searches are insensitive since in our case the heavier ewinos and the lightest neutralino are already well-separated in mass.
-20 -Several other searches for strongly interacting particles have been performed by both the ATLAS and the CMS collaborations. The inclusion of b-tagged jets, together with leptons can be relevant in our present context. However, these searches consider large jet multiplicity (typically ≥ 4 − 6 jets). Generic absence of large jet multiplicity in our situations make them immune to any constraint whatsoever derived from these searches.
The most relevant searches, in our case, involve multi-leptons and b-tagged jets along with / E T , low jet multiplicity [12,72,79]. 8 Out of the multi-lepton analyzes implemented in the CheckMATE version that we employed, the most stringent constraints appear to arise from the 3 + / E T final states, as well as from the ones with an opposite-sign di-lepton pair in the final states [12,13,72]. 9 We use MadGraph5-v2.4.3 [81] to simulate ewino pair/associated-production. Events are generated for pp → χ j χ k , (χ j ∈ {χ 0 i , χ ± 1 }), with up to one additional parton in the final state. These result in 10 (15) distinct production channels when 3 (4) light neutralino states are considered. For each production channel, 0.3 million parton level events are generated. We then use the built-in version of PYTHIA6 [82,83] for showering and hadronization and for decays of unstable particles. We have used the MLM [84, 85] prescription for the matching of jets from matrix elements with those from parton showers, as implemented in MadGraph.
Typically, on merging and matching of partonic jets, the number of simulated events per production channel reduces to around 0.2 million, on an average. Such a volume of generated event-samples is expected to be healthy enough to ensure a stable statistics and hence could be used for reliable estimates in subsequent analyzes. The cross sections for all the processes have been computed at the leading order in MadGraph. A flat K-factor of 1.25 [86] has been multiplied to the cross sections of all relevant ewino-pair production processes to factor in the approximate NLO+NLL contributions. This is expected to help CheckMATE make conservative estimates of the lowest values of the ewino masses that the recent LHC data could allow. Finally, we have used CheckMATE [54] (see also [87][88][89][90][91]) to examine the viability of the benchmark scenarios in the light of 13 TeV LHC results. CheckMATE reports an r-value for each of the benchmark scenarios where r = (S −1.64∆S)/S95 and 'S', ∆S and S95 denote the predicted number of signal events, its Monte Carlo error and the experimental limit on 'S' at 95% confidence level, respectively.
The benchmark scenarios in table 2 are so chosen that they yield r < 1 which, going by the CheckMATE convention, are dubbed 'allowed' by the LHC analyzes employed for the purpose. We are aware of a stricter criteria used in some literature (say, r < 0.67 [16]) for definiteness in such a conclusion. In that sense, our approach is only semi-conservative. Thus, the best that can be said about these points is that most of them are on the verge of being ruled out by the LHC experiments and might soon get to be so with some additional data. However, at present, they are indicative of how low a µ eff could still be viable under different 8 Final states involving 'τ ' leptons have also been considered in the literature [12,80]. However, our benchmark scenarios are not sensitive to the signal regions discussed in those works. 9 Final states with leptons and b-jets have been considered in ref. [13] and certain signal regions discussed there can be relevant for our present study. However, the experimental results have not been implemented in CheckMATE version we used and hence it is beyond the scope of the present work.
-21 -scenarios when the LSP is singlino-dominated. Table 2 reveals that µ eff as low as ∼ 200 GeV cannot yet be ruled out with a reasonable certainty.
Recently the ATLAS collaboration has analyzed 139 fb −1 of data and has derived constraints by studying the pair-production of charginos where they used the di-lepton + / E T data for the purpose [92]. However, the analysis assumes that the chargino decays 100 % of the times to χ 0 1 W ± . Since in the presence of a bino-like χ 0 2 , as demonstrated in our benchmark scenarios, there is a substantial contribution from the decay mode χ ± 1 → χ 0 2 W ± , the constraints derived from the above analysis do not apply directly to our cases.

Conclusions
A low value of µ eff is known to ensure an enhanced degree of 'naturalness' in a Z 3 -symmetric NMSSM scenario. An interesting possibility in such a scenario is a light singlino-dominated LSP DM. These two together form the edifice of a singlino-higgsino LSP as a possible candidate for the DM. Motivated by these, in this work, we have explored in some detail the viability of relatively low values of µ eff with the LSP being singlino-dominated.
We agree with the observations made in the recent literature that for a singlino-dominated LSP it is not easy to meet the relevant constraints from the DM and the collider sectors simultaneously. Compliance has been reported only when the higgsino-like ewinos are nearly degenerate with the singlino-like LSP. This ensures its efficient coannihilation with a degenerate higgsino-like state thus producing a relic at the right (experimentally observed) ballpark. At colliders, this presents a compressed spectrum that results in relaxed bounds on the higgsino-like states which could then be light and still evading generic searches.
We have presented new regions of the target parameter space (with relatively light singlino-like LSP of mass m h SM /2, with a purity level > 95% and with relatively small µ eff ) which exhibit such an overall compliance with experimental data. These regions comprise of theoretically much-motivated regions that offer DM-annihilation funnels in the SM-like Higgs boson, in the Z-boson and in the singlet-like scalars. The higgsino admixture in the LSP DM is anyway necessary to secure their optimal annihilation in order to find compliance with the observed relic density. However, this needs moderation since otherwise the cross section for DM scattering off the nucleon in the DMDD experiments (in particular, DMDD-SD) becomes too large and violates the reported bounds. We have demonstrated that allowing for a smaller value of M 1 and/or M 2 (∼ µ eff ) can be a helpful one-shot manoeuvre that could favorably tweak the dynamics and the kinematics simultaneously. This way it helps achieve the right balance among various relevant interaction strengths and decay branching fractions thus offering simultaneous compliance with data from both DM experiments and the colliders. In a sense, this presents the scope and the requirement of an indispensable and nontrivial tempering of the singlino-like LSP for the purpose. Further studies in the area of tempered neutralinos in the NMSSM are in progress [93].