Long-lived higgsinos as probes of gravitino dark matter at the LHC

We investigate the LHC sensitivity to supersymmetric models with light higgsinos, small R-parity breaking and gravitino dark matter. The limits on decaying gravitino dark matter from gamma-ray searches with the Fermi-LAT put a lower bound on the higgsino-like neutralino NLSP decay length, giving rise to a displaced-vertex collider signature. Using publicly available tools for simulation of signal, background and detector response, we find that higgsinos with masses of 100 – 400 GeV and R-parity violation of ζ ~ 10−8 – 10−9 can show up in the 8 TeV LHC data with 10 – 30 fb−1 of integrated luminosity. We demonstrate that in the case of a signal, the higgsino mass can be determined by reconstruction of the dimuon mass edge.


Introduction
The Large Hadron Collider (LHC) at CERN has brought the field of high energy physics into a whole new era. The observation by both the ATLAS and CMS experiments of a new resonance at ∼ 126 GeV [1,2] is consistent with the long-sought Higgs boson. If this is indeed the scalar particle needed to induce electroweak symmetry breaking in the standard model (SM), the search for new physics does not end there. Such a light scalar receives quadratically divergent corrections to its mass, leading to large fine-tuning within the SM. Furthermore, although the gravitational evidence for dark matter is strong (see e.g. [3]), its particle nature has yet to be determined.
The minimal supersymmetric standard model (MSSM) addresses both these issues but also faces new problems. While unbroken supersymmetry only introduces one new parameter, the higgsino mass parameter µ, soft supersymmetry breaking generically introduces many new parameters along with sources of large flavour-and CP-violation (see e.g. [4]).

JHEP02(2013)133
The latter problems would be alleviated if the superpartners have large (multi-TeV) masses. Indeed, this would be consistent with a relatively large value of ∼ 126 GeV for the Higgs mass and the absence of LHC signals for any other new particles so far.
The disadvantage of such a decoupling solution is that heavy superpartners possibly reintroduce fine-tuning in the electroweak sector. Nevertheless, a light-higgsino scenario (LHS), in which two higgsino-like neutralinos and a higgsino-like chargino are light (of the order 100 GeV) and the other superparticles heavy, is a theoretical possibility. Experimental limits from chargino searches tell us that µ cannot be zero and there is a priori nothing relating its value to the soft breaking parameters. Even if the electroweak scale can be obtained with cancelations of terms with large soft parameters, this would be spoiled again if µ is too large (see e.g. [5]).
In fact, a spectrum with higgsino masses around the electroweak scale and the other sparticles typically in the TeV range has been shown recently to arise in models of hybrid gauge-gravity mediated supersymmetry breaking [6]. Such models are motivated by ways of obtaining the MSSM and grand unification in string theory that feature a hidden sector of exotic states. They thus make use of the mechanism present in gravity-mediated supersymmetry breaking to generate an electroweak-scale µ, as well as of the advantages of gauge mediation by providing a large number of messengers.
In this set-up we have where M 1,2 are the bino-and wino mass parameters, and mass splittings in the higgsino sector of the order m 2 Z/|M1,2|. Because the higgsinos are nearly mass degenerate and the strongly interacting superparticles are out of reach, such a scenario within the usual MSSM is difficult to probe at the LHC [7]. As will be the subject of this work, the prospects change if we allow for lepton number violating, and therefore R-parity violating, couplings.
R-parity conserves baryon and lepton number and is imposed in the usual MSSM in order to forbid proton decay. Its conservation also renders the lightest supersymmetric particle stable, making the lightest neutralino a natural WIMP dark matter candidate. From a theoretical point of view, however, R-parity conservation is not particularly favoured and the stability of the proton can be ensured by requiring either baryon or lepton number violation to be small. In fact, it is a disadvantage of the MSSM LHS scenario that the relic density of neutralinos is too low, due to coannihilations, to provide the dark matter. Here, we will consider an extension of the MSSM with broken R-parity where the dark matter candidate is instead the gravitino.
The possibility of gravitino dark matter has an interesting connection to leptogenesis. While the gravitino is a prediction of the desirable promotion of supersymmetry to a local symmetry, it leads to the gravitino problem [8][9][10][11][12]. The thermal production of gravitinos depends on the reheating temperature and the gravitino and gluino masses [13]. In thermal leptogenesis, the lepton asymmetry is created through the decays of heavy right-handed neutrinos and then transferred to a baryon asymmetry via sphaleron processes. In order to generate enough CP asymmetry, as well as to account for the small neutrino masses generated via the seesaw mechanism, the right-handed neutrinos need to be very heavy JHEP02(2013)133 and therefore high reheating temperatures are required to produce them thermally [14][15][16][17]. This would also lead to gravitinos being produced in great abundance. If the gravitino is not the lightest supersymmetric particle (LSP), its decays would interfere with big bang nucleosynthesis (BBN) [18]. If the gravitino is the LSP and a dark matter candidate, the BBN bounds instead apply to the next-to-lightest supersymmetric particle (NLSP) which is meta-stable. This is where R-parity violation (RPV) comes in.
The requirement of successful baryogenesis puts an upper bound on the amount of RPV that can be allowed, by the condition that the baryon asymmetry is not erased before the electroweak phase transition sets in. A small amount of RPV, however, leads to 1) a gravitino lifetime exceeding the age of the universe, because of the double suppression of gravitino decays by the Planck mass and the small RPV coupling and 2) a sufficiently short lifetime of the NLSP to be consistent with primordial nucleosynthesis. This makes it possible to have a good gravitino dark matter candidate even with the high reheating temperatures needed for leptogenesis, thus solving the gravitino problem [19].
It has been shown that the gravitino can account for the observed dark matter abundance for typical gluino masses and different types of NLSPs [20]. In the LHS, we will have to allow for a very large gluino mass and we discuss in this work how the gravitino dark matter abundance can be accounted for in this case. There are also direct bounds on the RPV couplings from cosmology, which will be of particular relevance for the present work. Decaying dark matter of this kind would lead to a diffuse gamma-ray flux observable by the Fermi-LAT. The non-observation of such an excess gives an upper bound on the RPV and thereby a lower bound on the NLSP decay length [21]. The finite NLSP decay length leads to the prediction of the displaced vertex signatures at the LHC that we study here.
Displaced vertex signatures are extremely powerful in suppressing SM backgrounds and have been studied, in the context of different SUSY models, for macroscopic decay lengths ranging from 1 mm to hundreds of meters [22][23][24][25][26][27]. Unlike e.g. models where neutrino masses are generated by RPV, which lead to decay lengths up to 1 mm, the scenario described above predicts decay lengths that are orders of magnitude larger. This will lead to displacements of decay vertices in the outer layers of the multipurpose LHC detectors ATLAS and CMS, which motivates the muon signature that we consider.
Previous studies of RPV at the LHC have been motivated by the fact that bounds on sparticle masses can be weakened due to RPV, thus providing a possible explanation to why supersymmetry has not been discovered yet even if it is already being produced at the LHC [25,28]. Here, we will show that RPV can also allow for detection of weakly coupled new physics that generically is not probed by the usual LHC searches. We will discuss why also the signal we consider would hide from the LHC searches performed so far, and show that the search strategy that we propose here could reveal new physics already in the data accumulated during the LHC runs with proton collisions at 8 TeV center-of-mass energy.
In this work, we consider the MSSM extended by bilinear R-parity violating couplings. These violate lepton number and R-parity by the introduction of only a small number of free parameters. Baryon number is conserved also at loop level, ensuring the stability of the proton. The bilinear interactions can be rotated under the symmetries of the theory to allow for a description in terms of trilinear Yukawa interactions, which simplifies the JHEP02(2013)133 phenomenological analysis. In this framework, we study the case of a higgsino-like neutralino NLSP and the prospects for discovery with the data from proton collisions at a center-of-mass energy of 8 TeV. All the strongly interacting superpartners are assumed to be out of reach, and the higgsinos would be pair-produced via a virtual Z or W boson. The heavier neutralino and chargino will decay into the NLSP, and the NLSP will travel in the detector before decaying into SM particles, typically a W boson and a charged lepton in the case of a higgsino NLSP. NLSP decays into a gravitino LSP are, due to suppression by the Planck mass, orders of magnitude less probable. A clear signature arises from at least one of the two NLSPs decaying inside the detector, giving rise to two opposite-sign muons from a secondary vertex.
In section 2 we set up the framework of bilinear RPV, derive the branching ratios for the higgsino-like neutralino NLSP and the gravitino LSP and discuss the cosmological bounds on the RPV couplings. In section 3 we describe the LHC signature and analysis tools, and present the result of our detector level study for a few benchmark models. We conclude in section 4.

Decaying dark matter in the light-higgsino scenario
If we abandon the requirement of R-parity, the additional terms in the MSSM superpotential together with the soft terms introduce 99 new free parameters into the model [29]. By allowing only for the bilinear terms, baryon number is conserved, and the number of new parameters is reduced to 9 (c.f. [29,30]). Such a scenario can be realised through the spontaneous breaking of B − L, the difference between baryon and lepton number [19].
Compared to the case studied in [25] where the lightest neutralino was assumed to be bino-like, direct production of higgsino-like neutralinos will have larger cross-sections since they are not suppressed by mixing angles [4]. In this section we derive the branching ratios of relevance for our phenomenological study and introduce the relevant parameters for the study of RPV in the LHS.

Bilinear R-parity breaking
In the extension of the MSSM with bilinear R-parity breaking that we consider, mass mixing terms between lepton and Higgs fields appear in the superpotential, 1 as well as in the scalar potential, induced by supersymmetry breaking, These mixing terms, together with the R-parity conserving superpotential, 3)

JHEP02(2013)133
the scalar mass terms, 4) and the standard SU(3) × SU(2) × U(1) Y gauge interactions define the supersymmetric standard model with bilinear R-parity breaking. Note that the Higgs mass terms m 2 u and m 2 d contain the contributions both from the superpotential (2.3) and the soft supersymmetry breaking terms. For simplicity, we have assumed flavour diagonal mass matrices in (2.4).
As discussed in [21], it is convenient to work in a basis of SU(2) doublets where the mass mixings µ i , B i and m 2 id in eqs. (2.1) and (2.2) are traded for R-parity breaking Yukawa couplings. This can be achieved by field redefinitions: the standard rotation of the superfields H d and l i , followed by a non-supersymmetric rotation involving all scalar SU(2) doublets, where ε is the usual SU(2) matrix, ε = iσ 2 , and i as well as i are functions of B, B i , m 2 id , m 2 li , m 2 u and m 2 d [21]. The virtue of these two rotations are the disappearing mixing terms between the Higgs and the leptons doublets as well as vanishing sneutrino VEVs, which allows to calculate all RPV decays with usual Yukawa-like Feynman diagrams.
The R-parity breaking Yukawa terms contain couplings between gauginos, lepton doublets and Higgs doublets. After electroweak symmetry breaking one obtains new mass mixings between higgsinos, gauginos and leptons, where we have defined: Here g, g and h e ij are the SU(2) and the U(1) Y gauge couplings and the charged lepton Yukawa couplings, respectively. Note that one also obtains couplings of the bino and wino to the lepton doublets and the Higgs doublets [31] where κ i = i sin(−α) + i cos(α). (2.11) In the Higgs decoupling limit, α β − π /2, we have κ i ζ i . The Higgs decoupling limit is satisfied in the models considered in the present work. The details of the following derivation are given in appendix A. The gaugino and higgsino mass terms together with the mixing terms in eq. (2.7) represent the 7 × 7 neutralino mass matrix in the basis of gauginos b, w 3 , higgsinos h 0 u , h 0 d and the three gauge eigenstates of the neutrinos ν i (A.3a), and also the 5 × 5 chargino mass matrix of gaugino, higgsino and the gauge eigenstates of the charged leptons (A.3b). Both mass matrices have to be diagonalized by (bi-) unitary transformation matrices U in order to obtain the mass eigenstates of the neutralinos χ 0 and charginos χ ± , respectively. The currents which couple gauge fields to neutralinos and charginos are modified as well by these transformations and then depend on CKM-type matrix elements of neutral V (χ 0 ,ν) , charged V (χ 0 ,e) , which are functions of the transformation matrices U . Furthermore we have derived the coupling to the supercurrents U (χ 0 ,ν) , as well as the coupling to the lightest Higgs V (ν,χ 0 ) . These R-parity breaking matrix elements are calculated in appendix A and read where the photino matrix element and the photino mass parameter are defined as and m is the largest out of the supersymmetric mass parameters M 1 , M 2 and µ in the neutralino (chargino) mass matrix. Hence the approximated diagonalisation does not depend on the details of the supersymmetric spectrum, but in fact depends only on the ratio between the electroweak scale and the largest supersymmetric parameter. Figure 1: Gravitino decay into photon and neutrino.

Gravitino and neutralino decays
The partial width for a gravitino decaying into a photon and a neutrino (see figure 1) is given by [32] (2.14) Inserting the matrix element (2.12d) one obtains for the gravitino lifetime to leading order in m Z/m [21]: and V are the charged and neutral current matrix elements, whereas is the matrix element for the decay into the Higgs boson. We will evaluate these elements at the leading order using the expressions given in eqs. (2.12a), (2.12b) and (2.12c), respectively. We will assume throughout the rest of this work that the Higgs decoupling limit is satisfied. Figure 3a shows the phase space suppression factors that are important for small neutralino masses. We have assumed that the mass of the lightest Higgs boson is 125 GeV. The total neutralino decay width is given by the sum However, the evaluation of the partial decay widths in the case of a higgsino-like neutralino (see figure 3b) leads to the conclusion that the lifetime of the lightest higgsino can be estimated solely from the decay into a W boson and a charged lepton. Below, we only provide the approximative formula, where we have taken into account that tan β in hybrid gauge-gravity mediation is in general large. In the numerical evaluation, however, we use the full results. Using eq.  arrive at the following expression for the higgsino-like neutralino lifetime where f W is the phase space suppression factor appearing in eq. (2.18) and we have set m χ 0 1 µ. The neutralino lifetime depends on the neutralino mass, the gravitino mass and its lifetime, and additionally on the bino and wino mass parameters. Expanding the higgsino lifetime in µ /m allows us to arrive at a formula which is to first order independent of the higher neutralino mass scale m ∼ M 1 ∼ M 2 : The factor including the gaugino masses M 1 and M 2 depends only on their ratio and is e.g. in the case of unified GUT masses very close to 1 /8.

Cosmological bounds
Both lower and upper bounds on ζ can be derived from cosmology. The lower bound comes from the BBN contraints on the NLSP when the gravitino is the LSP. An upper bound can in principle be derived by demanding that the baryon asymmetry generated by leptogenesis is not washed out before the electroweak phase transition in the early universe [34][35][36][37]. However, the bound from the constraints on decaying dark matter from

JHEP02(2013)133
the Fermi gamma-ray searches is stronger. As we will see, for our analysis the lower bound is not very constraining while the upper bound will be the motivation for our LHC search strategy.
Having derived the decay widths of the gravitino LSP and the higgsino NLSP, we are now ready to estimate the gravitino mass range allowing for gravitino dark matter and successful leptogenesis [20]. This allows us to connect the results from gamma-ray searches with displaced neutralino decays at the LHC.

Big bang nucleosynthesis
To start with, we need to make sure that the decays of higgsino NLSPs do not interfere with BBN. Hence we demand that all higgsinos decay during the first 100 seconds of the universe [11,38,39]. Deriving the neutralino lifetime (2.22) without substituting the R-parity violating parameter ζ for the gravitino mass and lifetime leads to which then characterizes the lower bound on ζ.

Gravitino dark matter mass
The minimal gravitino mass is limited by the requirement that the gravitino abundance does not overclose the universe. Since gravitinos are produced in thermal SQCD scatterings gg →gΨ µ [40], the gravitino mass must increase with increasing gluino mass for a given reheating temperature. As we are interested in models in which the colored particles are inaccessible at the LHC, gluinos will typically be very heavy. For example, in the hybrid gauge-gravity mediation scenario in [41], which gives rise to a Higgs mass close to the tentative LHC result, the gluino mass is close to 4 TeV. In order to still allow for small gravitino masses, we will assume that the hot phase of the universe was created in the decay of the false vacuum of unbroken B − L [42,43]. Since right-handed neutrinos are created from B−L Higgs decays, this scenario allows for gravitino dark matter, leptogenesis and the correct values for the neutrino mass parameters while requiring lower reheating temperatures compared to the thermal leptogenesis case. The lower bound on the gravitino mass obtained in [43] for mg = 1 TeV is m min 3 /2 = 10 GeV. It is possible to scale this bound to other gluino masses using [42] m3 /2 = m min (2.24) Assuming a lower bound on the gluino mass of 2 TeV, the minimal gravitino mass is 40 GeV and therefore a neutralino NLSP with a mass of 100 GeV is viable.

Fermi-LAT bound on the gravitino lifetime
With the help of Fermi-LAT data we are able to restrict the lifetime of gravitinos for a given mass. Using the isotropic diffuse gamma-ray flux one can derive a lower bound of τ3 /2 3×10 28 s [21]. A stronger bound of τ3 /2 6×10 28 s can be derived as a consequence of

JHEP02(2013)133
the non-observation of any gamma-ray lines [44][45][46]. For a LHS with a bino mass of roughly M 1 ∼ 2 TeV this translates via (2.17) into an upper bound on the R-parity violation of ζ 4.70 × 10 −8 and ζ 3.32 × 10 −8 , respectively. 3 Finally we are able to derive a lower bound on the decay length of the lightest higgsino as a function of its mass as well as the mass and lifetime of the gravitino which is well within the reach of the multipurpose detectors at the LHC. Even when the higgsino decay length is larger than the detector dimensions, some higgsinos would, due to the statistical nature of the process, decay inside the detector.
3 Higgsino production and decay at the LHC In the LHS, higgsino-like charginos and neutralinos would be pair produced at the LHC via virtual Z and W bosons. Heavier higgsinos decay into lighter ones, the lightest one being the neutralino NLSP. The mass difference between the lightest chargino and the lightest neutralino is where we have used The mass difference between the next-to-lightest neutralino and the lightest neutralino is Hence, in the case of heavy gauginos the mass difference is rather small and to first order proportional to m Z/m. Therefore standard model products at this stage of the decay chain will be too soft to be detectable. In the presence of RPV, however, the NLSP would travel in the detector and further decay, prominently into a W boson and a lepton, yielding detectable SM objects coming from a displaced vertex. The lifetime of the lightest higgsino (2.25) is as a function of ζ Figure 4: Typical R-parity violating decay chain involving higgsino-like neutralinos at the LHC. The secondary vertices as well as the two possibilities of interesting muon combinations are highlighted. The Z boson decay is invisible, due to the small mass difference between the heavier higgsinos and the lightest higgsino (see eqs. (3.1), (3.3) and table 4). The signature is essentially the same for chargino production, since also in this case the decays into the lightest higgsino lead only to particles with small p T .

Signatures and search strategy
As shown in section 2, neutralinos that are produced in proton collisions at the LHC decay in a secondary vertex into a W boson and a charged lepton in almost 100 % of all decays. Figure 4 shows an example of a cascade decay with muons in the final state. The distance between the collision point and the secondary vertex depends on the decay width of the neutralino (2.19), and hence on the R-parity breaking parameter ζ.
The larger the value of ζ, the larger is the probability that the NLSP decays inside the detector. There are two multipurpose experiments at the LHC, the ATLAS and CMS detectors. Each detector consists of several subdetectors, from the inner detector, for track reconstruction and primary and secondary vertex reconstruction, to the calorimeters and outermost the muon system. Since the models under study here give rise to rather large displacements, we will choose to rely on the identification of muon objects in our search strategy. ATLAS has larger dimensions than CMS, with its muon system stretching out to a radius of about 20 m [47]. In our analysis, we choose to use requirements on the radial and beam-line coordinates corresponding to the CMS detector geometry, which will then be the somewhat more conservative choice. We do not expect our results to depend much on which of the two detectors is considered. The different detector layers of the CMS detector are [48]: • The inner detector or tracker, which stretches out to a radius of r ∼ 110 cm transverse to the beam. Its innermost part, the pixel detector, covers r ≤ 11 cm.
• The electromagnetic calorimeter which measures electron and photon energies and stretches to r 2 m. • The hadronic calorimeter, for measuring strongly interacting particles and identification of jets, which stretches out to r 3 m.
• The system of muon detectors or muon chambers, for identification of muons and measurement of their momentum, out to radius of r ∼ 7.4 m. Table 1 summarizes all possible LHC signatures if the NLSP is a higgsino-like neutralino in the case when both neutralinos decay inside of the tracker volume. The signatures are classified according to the final states in the neutralino decays. Unlike the previously studied case of bilinear RPV in the MSSM in [25], when the spectra contain light squarks and gluinos and a bino-like neutralino NLSP, there are only two types of signatures: leptonic signatures involving only charged leptons in the final state, and semi-leptonic signatures involving at least two charged leptons and jets.
Most supersymmetry searches for such final states at the LHC so far rely on the large production cross sections of the strongly interacting squarks and/or gluinos. The searches that are starting to probe direct electroweak production of charginos and neutralinos have been interpreted in R-parity conserving models with a stable neutralino LSP [49-51]. They rely on leptonic signatures with larger missing transverse energy than what would be present in our scenario in the case of both neutralinos decaying inside the tracker. Searches in the first LHC data for RPV have also been performed [52-59] but because of differences in the scenarios considered and differing signatures they do not apply here. Searches for long-lived neutral particles have been carried out as well but do not apply to our model as they assume either the wrong event topologies [  For smaller values of ζ, one of the neutralinos may decay inside or after the muon system leading to signatures with a larger amount of / E T as shown in table 2. We show in figure 5 how this possibly gives rise to a missing energy signature as ζ decreases. However, this situation would still not be covered by the current direct neutralino-and chargino searches by ATLAS and CMS [49-51] since the lepton reconstruction in these searches requires a track in the pixel detector, often with a certain maximum impact parameter to the primary vertex. For the same reason that we obtain missing energy from one of the neutralinos decaying outside of the detector, it will also be less probable that the other neutralino decays early enough for the leptons to fulfill such requirements on their inner tracks, as will be illustrated below.

JHEP02(2013)133
For very small R-parity violation both neutralinos may escape the detector without being observed. This case cannot be distinguished from the one with R-parity conservation and is for stable higgsinos and heavy coloured particles very challenging to detect [7]. This can be understood just by looking at the typical LHS process shown in figure 4 and trying to imagine possible signatures in the case of the lightest neutralino being stable.
We illustrate in table 3 how often the different situations of none, one or two of the neutralinos, produced in processes like in figure 4, decaying inside the detector occurs depending on the RPV and the higgsino mass. We also check how often the events would give inner tracks, possibly with a small impact parameter to the primary vertex, since this is a common requirement on lepton objects in existing searches. Table 3 illustrates: • For decreasing ζ, we approach a situation indistinguishable from the case of R-parity conserving LHS (see numbers in red).
• The second column in each subtable shows that when we can have large missing transverse energy in the event due to one of the neutralinos decaying outside of the detector, very few events will have the other neutralino decaying sufficiently close to the primary vertex for the event to be seen in existing searches for direct production of charginos and neutralinos (see numbers in blue).
• We also see that for a given ζ, increasing higgsino mass will lead to more decays inside the detector.
The leptonic decays of at least one of the neutralinos inside the detector lead to an opposite-sign dimuon signature, which will be exploited in the present work. We will focus    In the case where only one of the two neutralinos escapes the detector, corresponding to a large value in the second columns in table 3b, we notice an increased missing p T around the neutralino mass. In the case where both neutralino decay outside the detector corresponding to a large value in the first columns in table 3b, the values for small missing p T are increased. This reflects that the neutralinos are mostly back-to-back, due to the absence of other high p T objects in the decay cascade. on events with two opposite sign muons originating either in a secondary vertex in the tracker, far away from the primary interaction point, or having no associated track at all, being identified solely by the muon system. The muon identification and reconstruction process applied is described in [25] and will be briefly summarized in section 3.4.

Benchmark points
In models with hybrid gauge-gravity mediated supersymmetry breaking it is possible to realize a spectrum with higgsino masses around the electroweak scale, Higgs partners around the TeV scale and all other particles at the multi TeV scale [6].
The gravity mediated higgsino parameter µ of the superpotential and hence the lightest neutralinos and charginos can be of order 100 GeV. The existing lower bound on the chargino mass of roughly 95 GeV for degenerate spectra comes from LEP [64,65]. In this study we have chosen four benchmark points for which we have varied the higgsino mass parameter µ in three steps from the lower bound of about 100 GeV to 400 GeV, see  Table 4: Mass spectrum of light higgsinos in our four benchmark models with a higgsino mass parameter µ between 100 GeV and 400 GeV. All masses are given in units of GeV.
The masses of the MSSM Higgs particles are in this case set by the CP-odd Higgs mass parameter m A , which we have taken to be 800 GeV. All other particles are governed by the gauge mediated parameters m 0 and m1 /2 which are chosen to be 3 TeV, putting them out of reach of the LHC. In all our benchmark points the Higgs mass is around 125 GeV, in agreement with the observed Higgs-like resonance at the LHC [1,2]. Furthermore, while the LHCb result of an excess in the search for the rare decay B 0 s → µ + µ − [66] excludes many models with large tan β [67], the LHS is unaffected by this constraint due to the large mass splitting between the µ parameter and the squark masses.
We have used these parameter choices as input values for a full RGE calculation performed with SOFTSUSY. As expected the production cross sections for all supersymmetric particles except the light higgsino states are negligible. The higgsino production cross sections for the four benchmark points are listed in table 5.

Background
The SM processes that dominate the dimuon channel are: In table 6 we give the next-to-leading order (NLO) cross section for the processes that we have simulated for our study. As we will see, these will be efficiently removed by the requirement of a secondary vertex.
In our analysis, we will require the muons to be isolated, which efficiently removes leptons originating in jets, and we further remove possible contributions from displaced b quarks by a sufficiently large cut on tracks in the inner radius. For low background levels, however, other background sources might come into play. These are: • cosmic muons, • pion and kaon decays in flight, • hadronic punch-throughs,     An estimation of such contributions to our background has to be done with real LHC data, and is beyond the scope of this work. We argue here that most of this background, should it contribute, can be removed without significant loss of signal. Cosmic muons can be vetoed against by using the timing information, as discussed in [25], or a cut on back-toback muons. Punch-throughs are also not simulated in Delphes but should in principle be possible to veto since in this case the muon would be associated with a jet. Most of any possible contribution to displaced muons from decay in flight should be removed by our high p T requirement on muons. Pileup was estimated in a partly similar analysis to give a systematic uncertainty in the event selection efficiency of 2 % [63]. The displacement due to pileup is in general much smaller than the secondary vertices we are expecting. Therefore, such a background can be reduced by increasing the minimal impact parameter value required, which in our case of larger displacements would not lead to a large decrease in signal efficiency.

JHEP02(2013)133
In the following we therefore neglect these backgrounds to our displaced muon channel. However, as will be described in section 3.5, we will in our statistical analysis allow some margin for systematic uncertainty in case of vanishing estimated background levels by requiring our predicted signal to amount to a certain number of observed events.

Tools and settings
All Monte Carlo samples were generated using MADGRAPH 4.4.44 [71] interfaced with PYTHIA 6.4.22 [72]. We have chosen parton distribution functions given by CTEQ6L1 [73] and have used a modified version of SOFTSUSY 3.2.4 [74] for the calculation of the R-parity violating decays according to the formulas (2.18). All other signal decays were calculated using SDECAY [75].
The generic detector simulation DELPHES 1.9 [76], tuned to the CMS detector, 4 was used in order to account for effects of event reconstruction at the detector level. The finite radial size of the detector important for studies with secondary vertices has been taken into account for the case of muons following the muon reconstruction procedure described in [25]. We have reconstructed displaced muons with and without inner tracks by using the information at the event generator level to deduce where they come into existence. The generator level particles and the detector level objects are matched including the distance in pseudorapidity and azimuthal angle. In the reconstruction process we have assumed that a muon can be reconstructed as long as it is created before the muon system and that an inner track can be reconstructed as long as the muon is created within the first third of the tracker chamber.
After the muon reconstruction process we are left with two non-overlapping classes of muon-objects: first what we will refer to as chamber muons, that are identified solely by the muon system, and second the tracker muons having in addition the information about their origin from the associated inner track. 5 The simulated muon reconstruction efficiency is close to 80 % in the inner parts of the detector and drops rapidly if the muon is created in the muon system.

Cuts
We focus our search solely on the muon objects as we assume that we can trust the detector simulation results in this case even in the presence of secondary vertices. We aim to reconstruct clear signatures where one of the neutralinos decays into a muon and a W boson, which in turn decays into another muon and a neutrino. Thus, we demand two opposite sign muons that are either chamber muons or tracker muons with a secondary vertex far away from the primary interaction point. We remind the reader that displaced muons are not expected in SM processes giving two isolated muons, and we will show in this section that the requirement of large displacements efficiently removes the SM dimuon background. Furthermore, we require that the invariant mass of the two-muon system is not too small,

JHEP02(2013)133
thus suppressing back-to-back signals. This reduces not only the SM background, but also helps to decrease the background of cosmic muons.
Only DELPHES muon objects with a p T larger than 10 GeV are passed in our additional muon reconstruction processes described above. We demand that the secondary vertex of the muons lies before the muon system, meaning r(Vertex) < 4 m and z(Vertex) < 6 m, where r and z are the radial coordinate perpendicular to the beam and the coordinate parallel to the beam, respectively. We assume that tracks can be reconstructed reliably as long as the secondary vertex lies inside the cylinder defined by r(Vertex) < 40 cm and z(Vertex) < 1.3 m. The cut on the pseudorapidity for muons is taken to be η < 2.5. In addition to the isolation requirements in DELPHES, we only select muons that have no overlap with jet objects in angular coordinates ∆R > 0.1, where and ∆φ and ∆η are the appropriate angle differences between the muon and the jet in the usual detector coordinates. Then we perform a selection cut on the total number of muons (chamber and tracker muons) in the event: • N (muons) ≥ 2.
We define two event classes: • Class 1: the event contains exactly two tracker muons, i.e. N (tracker muons) = 2, with opposite charge.
• Class 2: the event does not fulfill the conditions for Class 1 and contains exactly two chamber muons, i.e. N (chamber muons) = 2, of opposite charge.
As the amount of R-parity breaking decreases, more and more events will fall into the second class. In accordance with the description of the signal presented above, we implement the following cuts on the Class 1 events: • m inv (µ + µ − ) > 5 GeV: we compute the invariant mass of the muon pair and demand that it is larger than 5 GeV.
• d(Vertex) > 5 mm: each of the tracks associated with the two tracker muons should have a vertex which is further than 5 mm away from the primary vertex. This value is approximately one order of magnitude larger than the resolution of the inner tracker [48,77].
• ∆d(Vertex) ij < 1 mm: the distance between the two track vertices should be less than 1 mm, in order to capture events where both muons originate in the same secondary vertex.
If the event fails one of the above cuts and fulfills the criteria for Class 2 events it is classified as Class 2:    • m inv (µ + µ − ) > 5 GeV: also in this case we demand that the invariant mass of the muon pair is larger than 5 GeV.
As expected, all of the SM background events are removed by the cut on the minimal distance of the vertex from the primary interaction point and the requirement that the reconstructed secondary vertices are close to each other (see table 7). The LHS model events, however, survive these cuts to such an extent that a signal is detectable (see table 8).

Mass determination
Rejecting the background-only hypothesis is only one contribution to the degree of belief that new physics has been discovered. Whether the signal hypothesis is a plausible one should be tested in other ways as well. Here, we show that the chosen signature allows for the determination of the neutralino mass via the well-known mass edge method [78,79].  (3.6) Following [80] we fold the phase space function with a Gaussian to model the mass edge: where the endpoint m cut and the height of the triangle σ are the free parameters to be fitted to the dilepton invariant mass m ll distribution to reconstruct the dimuon mass edge. We implemented this mass edge formula in the MINUIT class of the ROOT package. For this method to work, a sufficiently large sample of signal events is needed. In figure 6 we show examples of the mass edge reconstruction for different numbers of observed events in the case of our benchmarks model with µ = 100 GeV. We conclude from figure 6b that a total number of events between 26 and 50 should give an accurate estimate of the higgsino mass.

Result
The signal cross sections after cuts for all our LHS benchmark models are given in 0.028 ± 0.007 0.023 ± 0.009 0.014 ± 0.003 0.0037 ± 0.0008 Table 9: Signal cross sections after cuts for all benchmark models and different values of the RPV parameter ζ, in units of fb. The errors are Poisson errors and the center-of-mass energy 8 TeV.
When dealing with very low background levels, a Gaussian approximation may not be adequate, and one should assume the number of events to be Poisson distributed. Under the null hypothesis of B background events, the probability of observing N or fewer events is then given that the expectation value B is the true mean. We denote the expected number of events predicted by the model with S. To estimate the integrated luminosity 6 needed for a 5 sigma detection, one can require that there is a value for the minimum number of observed events N obs = N + 1 such that 1 − P (N ; B) < 2.9 × 10 −5 %, corresponding to 5 standard deviations in the case of a one-sided Gaussian. In addition, N obs has to satisfy (1 − P (N ; S + B)) being larger than some probability P obs for observation under the hypothesis of our model. 7 In our case, the expectation is B = 0, and in principle any S > 0 would constitute a signal. In a real measurement, however, the estimated background will be known only to a limited precision, and we will require N obs ≥ 5 in order to have some margin to allow for systematic uncertainties. In figure 7 we present results assuming a P obs of 50 %, 90 % and 99 %. We see that the integrated luminosities of ∼ 30 fb −1 expected in the 8 TeV run at the LHC would suffice for discovery of the lightest higgsinos with R-parity violation in the range ζ ∼ 2 × 10 −9 -2 × 10 −8 , and for ζ above 6 × 10 −9 masses of µ = 400 GeV may be reached.
After applying the mass edge method described in section 3.4.3 to one of our benchmark models, we estimate that 30 signal events are sufficient to reconstruct the neutralino mass with a couple of GeV's precision. Figure 8 shows the resulting integrated luminosities at which the higgsino mass could be determined at the LHC running at 8 TeV for our benchmark models.
We use our 8 TeV results to estimate the reach when LHC runs at the design centre-ofmass energy of 14 TeV, by the same statistical analysis applied after a naive scaling of the cross sections after cuts with the factors σ tot 14 /σ tot 8 presented in table 5. Since the background 6 The integrated luminosity is L = σN where σ is the cross section and N the number of events. 7 For a discussion of the statistical measures used for this kind of study, see appendix B of [81].  is assumed to be completely removed by our cuts, the reach at 14 TeV would be significantly improved. Larger higgsino masses can be reached at smaller integrated luminosities, as can be observed in figure 9. Also the luminosity which is needed in order to reconstruct the neutralino masses is reduced, as can be seen in figure 10. These results are approximate and we expect that changes in e.g. the p T cut on the muons may be needed to deal with systematic effects such as increasing pileup.

Conclusion
We have investigated the LHC detection prospects for the light-higgsino scenario, or LHS, in the MSSM extended with bilinear R-parity breaking terms. A spectrum with the higgsinolike neutralinos and chargino light, of the order of 100 GeV, and the other superparticles in the multi-TeV range is consistent with a Higgs mass of ∼ 126 GeV within the MSSM and

JHEP02(2013)133
can be obtained in GUT models. Because the higgsinos are nearly mass degenerate and the strongly interacting superparticles are out of reach, such a scenario within the usual MSSM is difficult to probe at the LHC. The prospects change if we allow for R-parity violation, which leads to a consistent cosmology where leptogenesis and gravitino dark matter can be accounted for without conflict with BBN.
We calculated the R-parity violating decay modes for a higgsino-like NLSP, and thereby diagonalized the full R-parity breaking neutralino matrix for this case. We showed that neither the resulting mass eigenstates nor the transformation matrices differ from the ones derived previously in the bino NLSP case. In fact, the parametrization only depends on a large hierarchy between the fermi scale and the largest neutralino mass.
We were thus able to use the limits on decaying gravitino dark matter from gammaray searches with the Fermi-LAT to put an upper bound on the RPV parameter ζ and thereby a lower bound on the higgsino NLSP decay length. This motivated our study of a dimuon LHC signature with macroscopic and, compared to other scenarios studied in the literature, relatively large displaced vertices. We simulated events and detector response for a few benchmark models, varying the value of the higgsino mass parameter µ and the RPV parameter ζ. We found that this scenario can show up already in the 8 TeV run at the LHC with the 10 -30 fb −1 of data expected at this center-of-mass energy, and that the reach in ζ is improved by up to an order of magnitude compared to the current reach of gamma-ray searches.
We also demonstrated that in the case of a signal, the LHS with RPV hypothesis can be tested further by a mass edge reconstruction. Except for the largest ζ and smallest higgsino masses considered, the luminosity accumulated during the 8 TeV run may not be enough. However, based on a our rough estimation of the 14 TeV reach, the higgsino mass in all our benchmark scenarios can be determined in the 14 TeV run, requiring integrated luminosities in the range 30 -1000 fb −1 .
The complementarity of gamma-ray and LHC signatures for the LHS with RPV may also be used to falsify the model. A future observation of a gamma-ray line consistent with decaying dark matter would fix the gravitino mass and the value of ζ. If also an LHC signal of a higgsino-like neutralino shows up, which may even be possible in the scenario considered here if the stops are light enough to be produced at the higher center-of-mass energy, the absence of any displaced decays would then render the LHS with RPV in conflict with experiments.

A Higgsino and gravitino branching ratios
The couplings of the gauge fields to charged and neutral matter are given by

JHEP02(2013)133
The currents in the gauge eigenstate basis are The upper indices of the electromagnetic currents indicate the transformation properties of the fields in the current under the SU (2) L . Having derived the currents we now have to transform them into the mass-eigenstate basis of the fermions. To this end, we have to diagonalize the mass matrices M N and M C , which in the gauge eigenbasis are 8 8 Note the extra factors of √ 2 in the charged mass matrix compared to [21].