Constraining sleptons at the LHC in a supersymmetric low-scale seesaw scenario

We consider a scenario inspired by natural supersymmetry, where neutrino data is explained within a low-scale seesaw scenario. We extend the Minimal Supersymmetric Standard Model by adding light right-handed neutrinos and their superpartners, the R-sneutrinos, and consider the lightest neutralinos to be higgsino-like. We consider the possibilities of having either an R-sneutrino or a higgsino as lightest supersymmetric particle. Assuming that squarks and gauginos are heavy, we systematically evaluate the bounds on slepton masses due to existing LHC data.


Introduction
The discovery of the Higgs boson in the 8 TeV run of the LHC [1,2] marks one of the most important milestones in particle physics.Its mass is already known rather precisely: m h = 125.09± 0.21 (stat.)±0.11 (syst.)GeV [3], and the signal strength of various LHC searches has been found consistent with the SM predictions.While this completes the Standard Model (SM) particle-wise, several questions still remain open, for example: (i) Is it possible to include the SM in a grand unified theory where all gauge forces unify?(ii) Is there a particle physics explanation of the observed dark matter relic density?(iii) What causes the hierarchy in the fermion mass spectrum and why are neutrinos so much lighter than the other fermions?What causes the observed mixing patterns in the fermion sector? (iv) What stabilizes the Higgs mass at the electroweak scale?
Supersymmetric models address several of these questions and thus the search for supersymmetry (SUSY) is among the main priorities of the LHC collaborations.Up to now no significant sign for physics beyond SM has been found.The combination of the Higgs discovery with the (yet) unsuccessful searches has led to the introduction of a model class called 'natural SUSY' [4,5,6,7,8,9,10,11,12,13,14,15].Here, the basic idea is to give electroweak-scale masses only to those SUSY particles giving a sizeable contribution to the mass of the Higgs boson, such that a too large tuning of parameters is avoided.All other particle masses are taken at the multi-TeV scale.In particular, a n.cerna@pucp.edu.peb thomas.faber@physik.uni-wuerzburg.dec jones.j@pucp.edu.ped porod@physik.uni-wuerzburg.demasses of the order of a few hundred GeV up to about one TeV are assigned to the higgsinos (the partners of the Higgs bosons), the lightest stop (the partner of the topquark) and, if the latter is mainly a left-stop, also to the light sbottom.In addition the gluino and the heavier stop masses should also be close to at most a few TeV.
Neutrino oscillation experiments confirm that at least two neutrinos have a non-zero mass.The exact mass generation mechanism for these particles is unknown, and both the SM and the MSSM remain agnostic on this topic.Although many ways to generate neutrino mass exist, perhaps the most popular one is the seesaw mechanism [16,17,18,19,20,21].The main problem with the usual seesaw mechanisms lies on the difficulty in testing its validity.In general, if Yukawa couplings are sizeable, the seesaw relations require Majorana neutrino masses to be very large, such that the new heavy states cannot be produced at colliders.In contrast, if one requires the masses to be light, then the Yukawas need to be small, making production cross-sections and decay rates to vanish.A possible way out of this dilemma lies on what is called the inverse seesaw [22], which is based on having specific structures on the mass matrix (generally motivated by symmetry arguments) to generate small neutrino masses.This, at the same time, allows Yukawa couplings to be large, and sterile masses to be light.
We consider here a supersymmetric model where neutrino data are explained via a minimal inverse seesaw scenario where the gauge-singlet neutrinos have masses in the range O (keV) to O (100 GeV).We explore this with a parametrization built for the standard seesaw, and go to the limit where the inverse seesaw emerges, such that Yukawas and mixings become sizeable.Although non-SUSY versions of this scenario can solve the dark matter and matter-antimatter asymmetry problems [23,24,25], we shall make no claim on these issues in our model.
In view of the naturalness arguments, we also assume that the higgsinos have masses of O (100 GeV), whereas the gaugino masses lie at the multi-TeV scale (see [26] for an example of such a scenario).In addition, we assume all squarks are heavy enough such that LHC bounds are avoided, and play no role in the phenomenology within this work 1 .In contrast we allow for fairly light sleptons and investigate the extent to which current LHC data can constrain such scenarios.For further studies in other regions of the parameter space, see [28,29,30].
This paper is organized as follows: in the next section we present the model.Section 3 summarizes the numerical tools used and gives an overview of the LHC analysis used for these investigations.In Section 4 we present our findings for the two generic scenarios which differ in the nature of the lightest supersymmetric particle (LSP): a higgsino LSP and a sneutrino LSP.In Section 5 we draw our conclusions.Appendices A and B give the complete formulae for the neutrino and sneutrino masses.

The Model
We add three sterile neutrino superfields νR,k , and assume conserved R-parity.With this, the superpotential reads as The corresponding soft SUSY breaking terms are given by For the neutrino sector we use a Casa-Ibarra-like parametrization [31,32], the details of which can be found in Appendix A. In this work, for simplicity, we shall use a non-trivial R matrix which will enhance the Yukawa couplings of the two heaviest neutrinos, allowing us to write: Here, m 3 is the mass of the largest light neutrino mass, M i are the masses of the heaviest neutrinos, and a = e, µ, τ .The parameters ρ 56 and γ 56 are the real and imaginary components of a complex mixing angle, appearing in the R matrix.The Z NH a factors [33] depend on the PMNS mixing matrix and ratios of light neutrino masses, and are in general of O (1).The only exception is Z NH e , which is slightly suppressed due to the small s 13 .
The Yukawas can be significantly enhanced by taking a large γ 56 .Furthermore, we can see that, if M 5 = M 6 , the two Yukawa couplings have the same size.From here, it is straightforward to redefine the sterile states, and demonstrate that the resulting mass matrix has the same structure as the one of the inverse seesaw.
For the sneutrino sector, we have written the full sneutrino mass matrix in Appendix B. For simplicity, we neglect terms proportional to Y ν , and take vanishing B ν and T ν 2 .In this case, we do not need to split the sneutrino fields into scalar and pseudoscalar components, and can work with the νL and νc R states.We can then approximately write the sneutrino mass matrix as: such that we can assume that three νi states shall be dominantly νL , and other three states shall be dominantly νR .Thus, we refer to them as L-sneutrinos and R-sneutrinos, respectively.In the following, we take m 2 L and m 2 νR flavour diagonal, so the only source of sneutrino mixing comes from M R , which is very small.Thus, we denote the Lsneutrinos through their interaction eigenstates (ν eL , νµL , ντL ), while R-sneutrinos are denoted as ν1,2,3 .Notice that in our results we use the full formulae shown in Appendices A and B.
As usual, the model contains neutralinos and charginos.As mentioned before, we assume that the gaugino mass parameters are much larger than the higgsino mass parameter |µ|.Therefore, the lightest states are two higgsino-like neutralinos χ0 1,2 and a higgsino-like chargino χ− which are nearly mass degenerate, see e.g.[34] for a discussion of the resulting spectrum.
The best way to probe this model and to distinguish it from the MSSM is by discovering and studying the R-(s)neutrino properties.However, it is clear that their direct production at the LHC is not a very efficient process, as the cross sections are proportional to the corresponding Yukawa couplings.A better way to generate them is through cascade decays of heavier particles, such as Lsleptons or higgsinos.Thus, in this work, we always consider L-sneutrinos heavier than R-sneutrinos.
In Figure 1 we present the cross sections for ẽL,R /ν L production at tree-level.Notice that these cross sections are the same for the μL,R /ν µL and τL,R /ν τ L flavours.It is well-known that QCD corrections shift these to larger values [35], so we apply an overall K-factor of 1.17.Note that the sum of the processes pp → ẽL ν * eL and pp → ẽ * L νeL has by far the largest cross section, followed by ẽL ẽ * L and νeL ν * eL pair production.In the following, we focus on the resulting signal from the decay of these states, as they explain the main features of our results.Nevertheless, in the numerical analysis we have included all possible processes, such that the available data is fully exploited.
Note that left-right mixing in the stau sector is large, such that both τ − 1 and τ − 2 have a relatively large τ − L component.This means that in the following we need to study the decays of both states.
The final states, and thus the signal, depend on the nature of the LSP, which can be either an R-sneutrino or a neutralino.Moreover, in case of an R-sneutrino LSP, we also have a different phenomenology depending on whether the higgsinos are lighter or heavier than the L-sleptons.In addition, in the case of higgsinos being lighter than the L-sleptons, we shall also have a dependence on the size of the small gaugino admixture to the physical charginos and neutralinos.
Before we give details of each scenario, we first review the relevant part of the interaction Lagrangian of lL and lR with charginos and neutralinos: where, for simplicity, we have neglected generation indices as well as left-right mixing.This is a very good approximation for the sneutrinos, the first two slepton generations, and for the staus in case of small to medium values of tan β.
The neutralino mixing matrix N is in the basis b, w3 , Hd , Hu , and in our model we have Knowing these couplings is very convenient at the time of understanding the different branching ratios.For example, if we want to compare the μL decays into higgsinos, for very large values of M 1 , M 2 , we find: where we see that decays into charginos will be subdominant if the mixing with gauginos is small enough.

Sneutrino LSP and light higgsinos (m
This scenario is characterized by subsequent two-body decays.The heavier L-sleptons decay into states involving χ0 1,2 or χ± 1 , which then decay into states involving Rsneutrinos.The decay chains have several branches, with the dominant branching ratio for L-sleptons depending on the size of the couplings and the respective elements of neutralino and chargino mixing matrices.
In the following, for each slepton, we compare two scenarios.In the first one, we set M 1 = M 2 = 2 TeV, such that there is a small but non-negligible gaugino admixture on the neutralino and chargino states.On the second one, both gauginos are "decoupled" from the model by setting their masses at 1 PeV.
For definiteness, we shall set m L = 600 GeV, µ = 120 GeV and tan β = 6.In Figure 2, we show the most important branching ratios of each slepton as a function of M 1 = M 2 , which shall now be discussed.Notice we neglect to comment those cases where two contributions interfere destructively, as this effect is not of our interest.
For the smuon μ− L , the 2 TeV gaugino scenario leads to primarily μ− L → µ − χ0 (75%), followed by μ− L → ν L χ− (25%).The latter decay is due to a gauge coupling, and its branching ratio vanishes when the gauginos decouple.In contrast, μ− L → µ − χ0 is due to a combination of Y µ and gauge contributions, and its branching ratio rises to unity in the decoupling regime.The muon sneutrino νµL follows a very similar pattern, with νµL → µ − χ+ dominating (70% at 2 TeV, and then 100% in the decoupling scenario), followed by νµL → ν L χ0 , which decreases in front of rising gaugino masses.
The case for the selectron ẽ− L is very similar to the one for μ− L in the 2 TeV case, replacing µ by e, and with very similar branching ratios.However, in the gaugino decoupling scenario, we find that the most relevant decay is ẽ− L → ν h χ− (90%), followed by ẽ− L → e − χ0 (9%) and ẽ− L → ν L χ− (1%).The reason for this is that the first decay proceeds through a Y ν coupling, which is larger than the Y e coupling that governs the second decay.Again, the electron sneutrino νeL decays are similar to the νµL for 2 TeV gauginos, and in the decoupling case they change to νeL → ν h χ0 (60%), νeL → e − χ+ (30%) and νeL → ν L χ0 (10%).
Stau decays are somewhat unique, as the mass eigenstates have large components of both τ − L and τ − R .We find that the lightest stau, τ − 1 , which is mostly τ − R , decays in equal proportions through τ − 1 → τ − χ0 and τ − 1 → ν L χ− (50%).In contrast, the heaviest stau, mainly τL , decays through τ − 2 → τ − χ0 (97%) and τ − 2 → ν L χ+ (3%).Ignoring interference effects, the branching ratios of both τ1 and τ2 are independent of the mass of the gauginos.The reason for this is that the stau states can always couple with higgsinos through the Y τ coupling, which is relatively large.For τ2 , the difference in the values of the branching ratios is due to τ − 2 → ν L χ+ being somewhat suppressed due to the need of LR mixing.
The tau sneutrino ντL follows a different pattern, as here there is no large mixing with any νR , such that gaugino couplings can play an important role again.For the 2 TeV gaugino scenario, the dominating decay is ντL → τ − χ+ (90%), followed by ντL → ν L χ0 (10%).The former increases to (100%) in the decoupled case, as in the νµL scenario.
In all scenarios, the charginos decay into a charged lepton, and a light sneutrino: χ− 1 → − ν1,2,3 .The charged lepton is usually a muon or a tau, due to the Z NH a factors in the Yukawa couplings, Eqs (2.3), with the branching ratio into an electron being below 10%.Moreover, due to their couplings, the branching ratio of the decay into ν1 is very suppressed, so charginos decay mostly into ν2,3 .In principle, these should decay further through 3-body processes into additional leptons and ν1 .However, we find ν2,3 to be long-lived, and escape the detector.Thus, charginos contribute to our signal with a charged lepton and missing energy.Neutralinos follow a similar trend, but decay into a light neutrino instead of a charged lepton.Thus, neutralinos can be considered missing energy.
Finally, one has to take into account the decays of the heavy neutrinos, which are independent of the SUSY scenario considered.The heaviest neutrinos form a pseudo-Dirac pair, and shall decay promptly [36].We shall concentrate on decays involving at least one charged lepton: ν 5,6 → − qq or ν 5,6 → −( ) + ν , with off-shell mediators.

Sneutrino LSP and heavy higgsinos
The situation changes drastically once the µ parameter is larger than the slepton mass.In that case, the previous decays are not possible, and one either needs to consider alternative two-body channels, or new three-body decays.We show the available branching ratios in Figure 3, as a function of the slepton mass, where we have fixed µ = 400 GeV.
For ẽ− L and μ− L , we find that the dominant decay is ˜ − L → ν L W − * , with the virtual W − giving jet pairs or a charged lepton plus a light neutrino.As usual, decays with quark final states have larger branching ratios.
Another possibility is to decay directly to an R-sneutrino and a real W − ( ˜ − L → ν2,3 W − ).This process depends on the small LR mixing in the sneutrino sector, so it is proportional to Y ν .We find that the branching ratio for μ− L is generally smaller than 20%.As mentioned earlier, the Z NH e factor in Y ν is slightly suppressed with respect to Z NH µ , so for ẽ− L the branching ratio is smaller.The stau sector has a slightly different phenomenology, due to the large left-right mixing.In particular, this leads the predictions of this scenario to depend strongly on tan β.The mixing splits the states, such that τ − 2 → τ − 1 Z 0 * decay is allowed.The inclusion of this new channel modifies the other ντL W − * and ν2,3 W − branching ratios.
The τ − 1 , on the other hand, for small m L, has similar ντL W − * decays, with non negligible νR W − * contributions.The reason for this is that the mixing-induced mass shift implies that τ − 1 is close in mass to the ντL , leading to a strong kinematical suppression.As a consequence the two-body decay into νR W − clearly dominates once kinematically allowed, despite the fact that there is only a small left-right mixing in the sneutrino sector.Thus, with the exception of τ − 1 , charged slepton decay shall usually produce one additional ν L .These states shall be accompanied by jets more than 50% of the time.We find that all νL flavours have the same behaviour.For low masses, the decays are governed by off-shell Z 0and h 0 -boson exchange, and a νR emission.For larger masses, the branching ratios are dominated by two-body decays into a νR and an on-shell Z 0 or h 0 boson, if kinematically allowed.When both bosons are accessible, the decays into the light Higgs have larger branching ratios.

Higgsino LSP (µ < m νR )
For completeness, we also study the case where the higgsino is the LSP.This is motivated by the fact that it has not been considered so far in the literature.
In this scenario the sleptons and sneutrinos have twobody decays only as described above for the case of M 1 = M 2 = 2 TeV.We note that, at the one-loop level, a mass splitting between the lightest neutralino and the chargino   is induced via the photon-loop yielding the contribution [34] This implies that the chargino will always have a sufficiently large decay width such that it decays inside the detector.However, due to the small mass differences, the decay products of the lightest chargino and the second heaviest neutralino are so soft that they mainly contribute to the missing transverse momentum.Note that, in this case, νi are hardly produced in the decays of the sleptons and heavier L-sneutrinos.

Set-up
For this investigation we have used a series of public programs: As a first step we have used SARAH [37,38,39,40,41] in the SUSY/BSM toolbox 2.0.1 [42,43] to implement the model into the event generator MadGraph5 aMC@NLO 2.5.2 [44].For each set of parameters, we use SPheno 3.3.8[45,46] to numerically calculate the mass spectrum and branching ratios.After generating the hard scattering in MadGraph, the showering and hadronization is carried out internally with PYTHIA 8.233 [47], which uses the CTEQ6L1 PDF set [48].We also use PYTHIA for the heavy neutrino decays.As default we generate 25000 events for every production process.The generated events are then fed into CheckMATE 2.0.7 [49,50], which uses Delphes 3.4.0as detector simulation [51].Given a specific experimental search, CheckMATE compares the number of events passing each signal region with the observed S95 limit obtained by the experiment via the parameter with S being the number of events in the considered signal region, ∆S the error from the Monte Carlo and S 95 obs is the experimentally observed 95% confidence limit on the signal [49,52].In our work, we indicate CheckMATE to compare our predicted signal with all of the available experimental searches.Notice that Eq. (3.1) does not capture all theoretical uncertainties, such as missing higher order corrections in the production and decays of the various particles.Moreover, there are effects due to variations of the input parameters.For example, in case of an νR -LSP, the charginos will decay into either a µ or τ plus one of the R-sneutrinos which escapes detection.The ratio of µ over τ depends on Y ν and varies with the choice of neutrino mixing angles.We therefore follow the basic idea presented in ref. [52] to capture such uncertainties: we do not take the r c = 1 value as sharp boundary but assume that all points with r c ≥ 1.2 (r c ≤ 0.8) are excluded (allowed) whereas for the range in between one would need a more detailed investigation.
For each point, we have also checked that the Z and Higgs invisible width respect experimental bounds, and that the heavy neutrino mixing is small enough to avoid direct detection [53].We have also checked that LFV processes such as µ → eγ, do not exceed the current constraints [54].In this scenario, the SUSY contribution to LFV is very small, either due to the heavy gaugino masses, or due to the small Yukawa couplings.Thus, the non-SUSY part dominates, and as one can see in [33], for the current choice of M 5 , M 6 and γ 56 , it does not saturate the bounds.

Higgsino LSP
We study first the case of a higgsino LSP.In these scenarios the sleptons decay directly into either a lepton and missing energy, or invisibly.The latter case occurs in case of l → ν χ− because the decay products of the charginos are very soft.
For the following investigation we have fixed M 1 = M 2 = 1 TeV implying that (i) the ẽ decays are mainly via  the small gaugino admixtures in the chargino and neutralinos and (ii) there will be practically no right-handed neutrinos produced in the slepton decays.From this point of view we are effectively in the usual MSSM with a higgsino LSP.However, to our knowledge the bounds on the slepton mass parameters due to the LHC data have not been presented in the literature.
In this scenario, the most important analysis is the search for two same sign leptons in combination with large missing transverse energy, carried out in [55].This leads to bounds on the m Ẽ − m L plane, which are shown in Figure 4 for the case µ = 120 GeV and tan β = 10.On this Figure, one can see that m L < 400 GeV is excluded, independent of m Ẽ .This constraint increases up to 500 GeV if, in addition, light ˜ − R are present.In contrast, even ˜ R with a mass of 200 GeV cannot be excluded, which can be seen in the Figure in the limit of heavy ˜ − L .We understand this is due to insufficient LHC data having been analysed.However, this might change in the near future, once the full 2016 data set has been investigated by ATLAS and CMS.
The structure for m L > ∼ 600 GeV and m Ẽ < ∼ 250 GeV can be understood from the interplay of different signal regions defined in [55].These regions differ mainly in the required bound on the 'stransverse' mass m T 2 [56,57]: m T 2 ≥ 90, 120, 150 GeV, corresponding to the signal regions 2LASF, 2LBSF and 2LCSF, respectively.Taking, for example, m L = 625 GeV, one finds that for m Ẽ = 200, 225, 250 and 275 GeV, the signal region 2LASF, 2LBSF, 2LBSF and 2LCSF is the most important one, respectively.
Last but not least, we remind that if the mass difference between the sleptons and the higgsinos gets too small, then the average value of the transverse moment of the lepton could be below 20 GeV.This can be a problem, as the p T cut for the leading (subleading) lepton in this search is of 25 GeV (20 GeV).In this situation, one cannot carry out any exclusions, as the final states are not energetic enough to pass the triggers.We note that the results hardly depend on the value of tan β, whose main effect is to enlarge the mass splitting of the staus for growing values.More important is the size of µ as this affects the kinematics, e.g.larger values of |µ| imply softer leptons for fixed slepton mass parameters.This is demonstrated in Fig. 5, where we display the r cvalue for different values of µ as a function of m L = m Ẽ .As we have mentioned previously, scenarios with r c values below 0.8 are allowed, the ones with r c > 1.2 are excluded whereas for those in between (gray band) no conclusive statement can be made.The structure close to the maxima of the different curves is again due to the interplay of the different signal regions.

R-Sneutrino LSP
As we have seen in Sections 2.1 and 2.2, on the R-sneutrino LSP scenario, different decays occur depending on the size of µ with respect to m L. Thus, in order to study this situation appropriately, we first need to understand the constraints on chargino pair production.
As mentioned previously, after production, each chargino decays into a ν2,3 and a charged lepton.Thus, the main constraints arise from the search for two leptons plus missing transverse energy at 13 TeV [55].For very small values of µ, additional constraints arise from the measurement of the W + W − cross section at 8 TeV, with subsequent decays of the W into leptons [58].In this case, the charginos would contribute more than what is allowed by the experimental uncertainty to the W + W − signal regions.
In Figure 6, we show the exclusion region in the µ − m νR parameter space, based on χ+ χ− production.From the plot, we see that, for vanishing m νR , the bound on m χ± can be as large as 375 GeV.Note that the ν2,3 , for m νR = 0, have nearly the same masses as the right-handed neutrinos.For relatively small values of µ, one finds that R-sneutrino masses lighter than µ − 75 GeV are ruled out, with the allowed region increasing for µ 250 GeV.For completeness, we note that this exclusion does not depend on the sign of µ or the value tan β.
In Section 2.1, we analyzed two scenarios for the gauginos, one where M 1 = M 2 = 2 TeV, and another "decoupled" scenario, where we set M 1 = M 2 = 1 PeV.Given our results for chargino production, we explore two additional possibilities.On the first one ("varying µ"), we set µ = m νR + 25 GeV, such that we always have m νR < µ < m L. On the second one ("fixed µ"), we set µ = 400 GeV, such that one also needs to take into account the m νR < m L < µ case.Thus, four different exclusion plots will be generated.In all of these, we merge exclusions from 8 and 13 TeV data.
We show the constraints on the varying µ scenario in Figure 7. Here, the relevant analysis is again [55], which searches for events with 2-3 leptons and missing energy.We find that, for both choices of gaugino mass, we can rule out values of m L as large as 575 GeV.In addition, for lighter slepton masses, it is possible to rule out Rsneutrino masses as heavy as 175-225 GeV, depending on the amount of gaugino admixture.
The exclusion for "decoupled" gauginos is stronger, which can be understood from Figure 2.For ˜ L , all possible decays shall lead at least to one charged lepton, for all values of gaugino mass (recall that, in this scenario, charginos decay always into final states with visible charged leptons).However, for ν L , one finds that it is possible to have only missing energy on the final state, due to ν L → ν L χ0 decay.This decay channel is suppressed in the "decoupled" scenario, meaning that it is much more likely to have energetic charged leptons on the final state, which strengthens the multi-lepton signal.
The constraints on the fixed µ case are shown in Figure 8, again, for different values of gaugino mass.Here it is very interesting to note that there are very weak con-      straints when m L < µ.The reason for this can be found in Figure 3.We see that in most of the cases, we have the charged slepton decaying into ν L and light fermions.Given the proximity in the slepton masses, most light fermions end up being very soft, and elude detection.On the other hand, when the two resulting ν L decay, the decay products shall involve either an on-shell h 0 or Z 0 .This is again problematic for detection, as fermions coming from these states are generally avoided in new physics searches by suitable cuts to suppress SM background.This leaves us sensitive only to the very low m L region, where L-sneutrino three-body decays are allowed.
For large values of m L, we return to the µ < m L scenario.Here, again, we have [55] giving the relevant constraints.The bounds reach m L as large as 575 GeV for vanishing m νR .Morover, for smaller charged slepton masses, we can bound m νR up to 250 GeV.This is all consistent with our results for the varying µ scenario.

Conclusions
In this work, the MSSM was extended by three righthanded neutrino superfields, with active neutrino masses being provided through the seesaw mechanism.In addition, driven by naturalness arguments, the µ term was kept relatively small, such that the lightest neutralinos were higgsino-like.We considered LHC data on this model, and explored how much existing data constrain such scenarios.
Two possibilities were considered for the nature of the LSP.On the first one, this was a higgsino-like neutralino.In this case, one requires a non-vanishing gaugino admixture in order not to have too long-lived charginos.The sleptons would decay into SM particles and neutralinos.We carried out an analysis considering m L = m Ẽ , for fixed neutralino mass, and found that only m L could be bounded.For µ = 120 GeV, we can rule out at least m L < 400 GeV for all values of m Ẽ , and m L < 500 GeV for m Ẽ = 200 GeV.
On the second possibility we considered, the righthanded neutrino superpartner, the R-sneutrino, was taken as the LSP.This provided a very complex scenario, depending on the relative size of the neutralino / chargino mass with respect to the slepton mass.The phenomenology also depended on the amount of gaugino component within the neutralinos.We found that, as long as µ < m L, we can exclude slepton masses to a maximum of 575 GeV, for vanishing m νR .For lower values of slepton mass, the R-sneutrino masses can be excluded up to about 200 GeV.In case m L < µ, constraints became very weak, as final states were either too soft, or excluded from signal regions.
2. Moreover, U and V are the chargino mixing matrices, in the basis w± , H± , such that in our model we have |U 11 |, |V 11 | |U 12 |, |V 12 |.In addition, we know that Y l , Y ν g , g, the only exception being Y τ , which can become O (g ) in case of very large tan β.

1 M 1 Fig. 2 .
Fig.2.Branching ratios for sleptons as a function of gaugino mass M1 = M2, for µ < mL Decays for charged sleptons (sneutrinos) are shown on the left (right) column.The last panels describe the colours for the branching ratios shown in each column.In case of neutralinos, the sum over the two lightest states is shown.

Fig. 3 .
Fig.3.Branching ratios for sleptons as a function of slepton soft mass mL, for mL < µ 400 GeV.Decays for charged sleptons (sneutrinos) are shown on the left (right) column.The last panels describe the colours for the branching ratios shown in each column.

Fig. 4 .
Fig.4.Constraints on combinations of m Ẽ and m L due to slepton/sneutrino production in case of a higgsino LSP with M1 = M2 = 1 TeV, µ = 120 GeV and tan β = 10.Red points are excluded, blue ones are allowed and in case of the green ones no conclusive statement can be drawn, within the known theoretical and experimental uncertainties.

Fig. 5 .
Fig. 5. r c as a function for tan β = 10, M1 = M2 = 1 TeV and µ = 120, 150, 200 and 250 GeV, respectively.The grey band (0.8 ≤ r c ≤ 1.2) gives the region where one cannot draw a conclusion whether the point is allowed or not, values below are allowed and those above are excluded.