Gravitino vs Neutralino LSP at the LHC

Using the latest LHC data, we analyse and compare the lower limits on the masses of gluinos and the lightest stop in two natural supersymmetric motivated scenarios: one with a neutralino being the lightest supersymmetric particle (LSP) and the other one with gravitino as the LSP and neutralino as the next-to-lightest supersymmetric particle (NLSP). In the second case our analysis applies to neutralinos promptly decaying to very light gravitinos, which are of cosmological interest, and are generic for low, of order O(100) TeV, messenger scale in gauge mediation models. We find that the lower bounds on the gluino and the lightest stop masses are stronger for the gravitino LSP scenarios but, nevertheless, the low messenger scale still ameliorates the fine tuning in the electroweak potential.


Introduction
In this paper we analyse the latest LHC data, to obtain and compare the bounds on the gluino and the lightest stop masses in two general scenarios, one with neutralino as the lightest supersymmetric particle (LSP) and the other one with gravitino as the LSP. Focusing on the bounds for gluino and the lightest stop is motivated by the question of naturalness of the Higgs potential in supersymmetric (SUSY) models. It has been well appreciated already for many decades that those particle masses, together with the higgsino mass, are most crucial for the degree of fine tuning; see Ref. [1,2] for a recent analysis. Our analysis is performed in the framework of three simplified models described in detail in Section 2. The main assumption here is that only the stop, gluino and the neutralino LSP (and the gravitino in the second scenario) are relevant for the production and decay processes leading to the collider signatures under study. It is assumed that the other superpartners decouple because they are either heavy or have small cross sections. It should be stress that such a situation indeed occurs in many explicit models studied in the literature. The three simplified models correspond to the lightest neutralino being pure bino, wino or higgsino. In the scenario with neutralino LSP, the bounds on the the gluino and stop masses have been extensively studied using the LHC data in experimental and theoretical papers, see e.g. [3][4][5][6][7][8][9]. We revisit this class of simplified models taking into account recent LHC results and directly compare with the gravitino LSP scenarios.
For the case with gravitino LSP, we assume that neutralino is the next-to-lightest supersymmetric particle (NLSP) particle decaying promptly, with cτχ0 1 1 mm, and focus on the corresponding LHC signatures 1 [11] (for a similar analysis based on the Tevatron data; see Ref. [10,13]). It follows from this assumption that the obtained bounds apply when gravitinos are very light, with masses m 3/2 1 keV for neutralinos lighter than O(1)TeV (see Section 2 for details). For substantially heavier neutralinos it is not possible to apply our recasting strategy since the neutralino becomes long-lived.
Gravitino LSP is generic for gauge mediated SUSY breaking (GMSB) scenarios. Its mass reads m 3/2 = F √ 3M pl , where F is the supersymmetry breaking F -term (or their combination), so the gravitino mass is model dependent. In general, gaugino masses in the GMSB scenario are given by one-loop contributions, roughly of the order m SUSY ∼ 0.01 F Λ . This means the messenger scale Λ is linearly related to m 3/2 , for a fixed m SUSY , and the mass range of gravitinos considered in this paper is consistent with the region Λ 100-1000 TeV. Very light gravitinos are theoretically interesting for several reasons. One motivation is to relax the apparent fine-tuning in the Higgs sector thanks to low values of the messenger scale. 2 Very light gravitino, m 3/2 O(10) eV, is also motivated by cosmology since it can evade the gravitino problem [15,16] and simultaneously allow for very high reheating temperature, convenient for leptogenesis. For high enough reheating temperature, gravitinos can be in thermal equilibrium in the early universe. In this case, the relic abundance of gravitinos can be approximately estimated in terms of the relativistic degrees of freedom, g * 3/2 ∼ 100, at the gravitino decoupling epoch as (see e.g. [18]) Thus, the relic abundance of gravitinos, that were originally in thermal equilibrium, exceeds the observed abundance of dark matter (DM) for m 3/2 100 eV. On the other hand, fits to the matter fluctuations at small scales (small scale structure formation) put a lower bound on the mass of a warm dark matter particle, m 3/2 O(1) keV, if it fully accounts for the observed DM abundance [17]. Therefore, the light thermal gravitinos cannot be the dominant component of the dark matter. Assuming that DM consists of the light thermal gravitinos and some additional cold dark matter (CDM) constituent, Ω DM = Ω CDM + Ω 3/2 , the constraint from the Ly-α and sets the upper bound, m 3/2 < 16 eV [17]. More recent study [18] sets the tighter upper bound, m 3/2 < 4.7 eV, using the data of the cosmic microwave background (CMB) lensing collected by Planck and of cosmic shear measured by the Canada France Hawaii Lensing Survey, combined with analyses of the primary CMB anisotropies and the baryon acoustic oscillations in galaxy distributions. The light thermal gravitino cannot be the dominant component of the dark matter but remains cosmologically interesting for the reasons mentioned above. Moreover, in presence of very light neutralinos, light gravitinos are even favoured [19]. Its mass range suggested by the mentioned above upper bounds corresponds to the O(100) TeV messenger masses in the gauge mediation models. Very light gravitino on the other hand, in the whole mass range considered in this paper, remains a viable dark matter matter candidate if produced thermally in collisions of the Standard Model particles at sufficiently low reheating temperature [16].

The light gravitino and model setup
Our study is organized around the question of how the presence of gravitino as the lightest supersymmetric particle-to this end at least from a collider point of view-changes exclusion limits at the LHC. We focus on models that have been motivated by the concept of natural SUSY, hence we consider spectra where among light particles we have gauginos/higgsinos, serving as the (next-to-lightest) supersymmetric particle, stops and gluinos. Other SUSY particles are assumed to be heavy, i.e. outside the collider reach and with a negligible impact on decay chains of light particles.
We examine three types of the simplified models: • a stop model, where gluinos are decoupled; • a gluino model, where stops are assumed to be heavy, however the gluino decays are mediated by off-shell stops; • a stop-gluino model where both particles are within or close to kinematic reach of the LHC and no definite mass hierarchy is assumed.
In all the models the strongly interacting particles are accompanied by electroweakinos: either the higgsino doublets, wino triplet or bino singlet -which are assumed to be lighter than the stops and gluinos. Finally, the models can be completed by adding the light gravitino. Generally the mass hierarchy is as follows: These models are parametrised by the masses of participating particles: m g , m t 1 , m χ . Since we assume that only one type of electroweakinos is present at a time, the mixing between gaugino-higgsino states is neglected, which implies that the additional charginos and neutralinos are almost mass degenerate with the lightest neutralino.

The electroweakino sector
The electroweakino sector is characterised by soft SUSY-breaking masses of bino and wino, M 1 and M 2 , respectively, the higgsino mass µ and the ratio of Higgs vacuum expectation values tan β. The neutralino gauge eigenstates ψ 0 = ( B, W 0 , H d , H u ) T , are related to mass eigenstates by means of the neutralino mixing matrix N ij ; see e.g. Ref. [20,21] for details. For our purposes the most important property of the lightest neutralino is its gauginohiggsino composition. While this can be completely general, in our simplified model setup we only have one type of light electroweakinos. In such a case, the mixing will be small and can be neglected. Therefore, we classify our models according to the dominant component of the NLSP in three types, defined in terms of the neutralino mixing components as 3 bino : N 11 ≈ 1, N 1i ≈ 0 (i = 1) wino : For each of these cases a different decay pattern emerges which additionally depends on the NLSP mass and, for higgsinos, on tan β. Taking the bino mass M 1 much smaller than the wino, M 2 , and higgsino, µ, masses we obtain the bino scenario. It means that the relevant electroweakino sector consist of just one lightest neutralino. Because of its vanishing coupling to the gauge bosons the bino mass is currently unconstrained in collider experiments.
The wino scenario is realized when M 2 M 1 , |µ|, there are two charged states together with the neutralino, which we write W + , W − and W 0 , respectively. Since they form an SU(2) triplet in the limit v → 0, they are almost mass degenerate. The charged and neutral states will split after including the electroweak symmetry breaking effects as well as the radiative correction. The size of this mass splitting is in general small and irrelevant from the collider perspective, unless it is so small that the charged state becomes long-lived.
In such a situation, the long-lived charged wino gives rise to so-called disappearing track signature. The ATLAS and CMS collaborations have analysed such a signature and put a strong constraint on the long-lived wino, m W + 460 GeV [22,23]. In this paper, we do not consider such a case and assume that the charged winos are short-lived at the LHC.
For |µ| M 1 , M 2 we obtain the higgsino scenario for which there are two neutral and two charged states coming from four components of the two higgsino doublets, h 0 1 , h 0 2 , h + , h − . Up to the electroweak symmetry breaking and radiative corrections, they are nearly mass degenerate. On the other hand, the mass splitting among the higgsino states is in general larger than that for winos, and we again assume that the higgsinos are short-lived in our analysis.

NLSP decay to gravitino
In models with the gravitino LSP we assume it is light enough so that it can be treated as massless at the LHC. We also assume that the decay of the NLSP neutralino to gravitino is prompt, taking a conservative limit cτχ0 1 1 mm, cf. Ref. [24][25][26][27][28][29]. The partial decay rates of the lightest neutralino into the gravitino are given by [10,30] where N ij is the neutralino mixing matrix and In the left panel of Figure 1 we plot contours of a fixed neutralino lifetime cτ χ 0 1 = 1 mm in the gravitino-neutralino mass plane. The three contours correspond to the lightest neutralino which is predominantly bino (red-solid), wino (blue-dashed) and higgsino (pinkdotted-dashed). The prompt region (cτ χ 0 1 < 1 mm) is located above contours (the top-left part of the plots), allowing fairly light gravitinos with m 3/2 5 eV-1 keV for neutralinos lighter than O(1)TeV. It justifies our assumption that the gravitino can be treated as a massless particle in dealing with its kinematics at colliders and we conveniently fix m 3/2 to 1 eV throughout our analysis. The right panel of Figure 1 recasts the calculation in the Λ-m NLSP plane, where Λ is the messenger scale and assuming 0.01 F Λ = 1 TeV. The branching ratio of the bino-like neutralino can be obtained by substituting N 1i = δ 1i into Eq. (4), which leads to Numerically, the χ 0 1 → γ G mode dominates the Z-boson mode for any bino mass. Therefore the models are mainly constrained by the analysis targeting photon final states as we will see later.
For the wino-like NLSP the branching ratio of W 0 is obtained by taking N 1i = δ 2i in Eq. (4), which gives Compared to the bino-like neutralino, the branching ratio to the photon final state is suppressed by the weak mixing angle squared, s 2 W , and the Z-boson mode is dominant for winos heavier than m χ easy to see that BR( h 0 1 → γ G) = 0 and In the large tan β and heavy h 0 1 limit, these mode will have the equal branching ratio of 50%, though the h 0 1 → Z G mode is generally favoured due to the difference in phase-space and tan β effect.
In Figure 2 we show the branching ratios of different classes of the NLSP. For binos the dominant decay mode is to the photon, regardless of its mass. For light winos photonic decay mode dominate as well, however for the heavier winos the dominant decay mode is to Z boson. Finally higgsinos decay either to the Higgs boson or Z boson, and for higgsinos heavier than 200 GeV either decay mode has a similar share.

Naturalness
One of the motivations for a light gravitino is to relax the apparent fine-tuning in the Higgs sector. In the leading-log approximation, the Higgs mass-squared parameters in a moderate or large tan β regime are roughly given by where Q is the electroweak symmetry breaking scale and Λ is the messenger scale of SUSY breaking. The first, second and third terms come from tree, one-loop and two-loop level, and are sensitive to the higgsino mass (µ), the stop mass (m stop ) and the gluino mass (M 3 ), respectively. In our analysis we take all the mass parameters real and positive for simplicity. It is clear that the small µ is crucial for naturalness. From the above formula it is also evident that the second and third terms can be made not-too-large by taking Λ to be small. The right panel of Fig. 1 shows the region of (m χ 0 1 , Λ) that is consistent with the prompt decay requirement cτ χ 0 1 1 mm. One can see that our prompt decay requirement is consistent with the region Λ 100 − 1000 TeV, where the fine-tuning can be largely relaxed for given µ, m stop and M 3 .
In particular, for a meaningful estimation of the contributions from each of the three terms in Eq. (9) we have to specify the value of µ. Therefore, we will discuss the impact of the gravitino LSP on the fine-tuning problem only for the higgsino-like neutralino case. In the other two cases, we have assumed that higgsinos are irrelevant for the collider signatures which, in practice, means that they are heavier than stops and gluinos. The limits obtained for stops and gluinos depend on this assumption. The first term in Eq. (9) is then the most important one irrespectively of the value of Λ. 4

Recasting LHC analyses
We confront our simplified models with various ATLAS and CMS analyses searching for beyond the Standard Model. For each analysed model we generate a grid of points each described by masses of particles characteristic for the model. The branching ratios are calculated using SDecay [31]. The exclusion of each point is determined using CheckMATE-2.0.26 [32].
The events are generated within the CheckMATE framework using Pythia 8 [33] including parton shower and hadronization. The signal is simulated separately for the coloured particle production and electroweak particle production. The cross sections are rescaled to the next-to-the-leading-log (NLL) accuracy computed with NLL-fast [34] and Resummino [35] for coloured and electroweak particles, respectively. The fast detector simulation is carried out using Delphes 3 [36]. Finally, CheckMATE emulates the various ATLAS and CMS analyses, estimating the signal efficiency and confronting the signal yield with the model-independent 95 % CL upper limit for each signal region. In order to see the impact of individual analysis, we do not combine multiple analyses, but rather show the exclusion limit obtained from single analysis.
The analyses used in this study are summarised in Table 1. The above Monte Carlo chain is applied to the all analyses except for cms_1801_03957 [37]. In cms_1801_03957, the dedicated analysis has been carried out for the production of charginos and neutralinos decaying to a massless gravitino and a Z-boson or Higgs boson. The limit has been derived as a function of Br(h 0 → hG) with Br(h 0 → hG) + Br(h 0 → ZG) = 1. The limit is strongest (mh 750 GeV) for Br(h 0 → hG) = 1, whilst most conservative (mh 650 GeV) for Br(h 0 → hG) ∼ 0.5. Guided by our considerations in the previous section, see Figure 2, we impose the limit mh 650 GeV on our higgsino-gravitino simplified models.
As we shall see in the following sections, apart from the already mentioned cms_1801_03957, three other searches included in the current analysis will define exclusion limits. The first one is atlas_1709_04183 which is the search for direct production of top squarks. The analysis

Results: the LHC limits
Our presentation of the exclusion limits is organised in the following way. For each of the three models: stop, gluino and stop-gluino we dedicate a separate subsection. There we specify to models for which the strongly produced particles are accompanied by different classes of electroweakinos: bino, wino or higgsino, and finally we compare exclusions for the case with the electroweakino or gravitino LSP.

Stop simplified model
We start the analysis of the LHC constraints by looking at the simplified model with stops and electroweakinos. The pattern of stop decays will crucially depend on the nature of electroweakinos. Gluinos are assumed to be heavy, m g 2.5 TeV , which is well above the current limits. In the simplest scenario, when mB mW , mh, the lightest neutralino is predominantly composed of the bino. The only available decay mode is: provided mt 1 > m t + mχ0 1 . In the following analysis we assume that this relation is satisfied. Otherwise,t 1 may decay into bWχ 0 1 , bjjχ 0 1 or cχ 0 1 depending on the mass spectrum and the parameters.
In the wino-like neutralino scenario, the lightest stop will decay though itst L component into the winos. There are two possible decay modes and in the limit of heavy t 1 , we have, stop-wino : For smaller stop mass, the phase-space factor becomes important, which further favours the t 1 → b W + mode. In particular, for m t 1 < m W + m t , the the top-quark decay mode vanishes and Br( t 1 → b W + ) = 100%. In the higgsino scenario there are three competing stop decay modes: where the branching ratios are in the m t 1 mh limit. On the other hand, for lighter stops the t 1 → bh + mode is preferred due to larger phase-space. In particular Br( t 1 → bh + ) = 100% for m t 1 < mh+m t . In our Monte Carlo simulation the branching ratios for different parameter points are obtained by SDecay [31].
The LHC constraints on the stop simplified model are summarized in Fig. 3, without gravitino (left column) and with gravitino (right column). For the model without gravitino the strongest constraint comes from the ATLAS search for direct stop production atlas_1709_04183 which targets stops decaying into top and the neutralino or into a bottom quark and chargino with subsequent decays via (offshell) W 's. The limit is the strongest for the bino LSP and extends up to mt ∼ 1000 GeV which is consistent with the ATLAS results. For the wino and higgsino case the limit is slightly weaker due to different competing decay modes and extends up to mt = 850-900 GeV.
The exclusion limit is drastically changed when a light gravitino is added as the lightest supersymmetric particle. In the NLSP bino case (right top plot in Fig. 3), we find atlas_1802_03158 (GMSB; γ(≥ 1) + E miss T ) [39] is by far the strongest, excludingt 1 up to mt 1 1250 GeV for mB ∼ 200 GeV and mt 1 1400 GeV for mB ∼ 1000 GeV. This analysis have managed to reduce the Standard Model background drastically by requiring at least one energetic photon in conjunction with large E miss T . Unlike the previous case without gravitino, the limit is not weakened for heavier bino, because the photons produced from B decays are more energetic for larger mB in general.
In the middle row we see the limit on the stop-wino-gravitino model. As in the bino-like χ 0 1 scenario with the light gravitino, atlas_1802_03158 is the most sensitive across the plane. The exclusion is driven by the electroweak production of chargino, χ + 1 χ − 1 , and charginoneutralino pairs, χ ± 1 χ 0 1 . In the heavy stop region the lower limit on winos is m W 900 GeV which is slightly lower than the corresponding ATLAS limit. For lighter stops due to the  additional small contribution from stop production the limit goes as far as m W 1050 GeV. Generally we find the limit, m t 1 1050 GeV, which is slightly weaker than the corresponding limit in bino-like χ 0 1 scenario. This is because the branching ratio for χ 0 1 → γ G is suppressed by s 2 W compared to the bino-like χ 0 1 case. We see that atlas_1802_03158 excludes winos below ∼ 900 GeV almost model-independently, as long as a massless gravitino is present in the spectrum.
The right panel shows the limit on the stop-higgsino-gravitino model. As can be seen, there is a strong constraint on the higgsino mass, mh 650 GeV, coming from cms_1801_03957 (electroweak; diboson final states [W, Z, h]) [37]. This limit has been obtained in cms_1801_03957 by a dedicated analysis targeting pair productions of higgsino states decaying into a massless gravitino. For the higgsino heavier than 650 GeV the relevant limit on the stop mass comes from atlas_1709_04183 (t; b-jets + E miss T ) [6], which exclude the stop up to ∼ 750 GeV.

Gluino simplified model
The second type of the analysed simplified model assumes that the only SUSY particles accessible at the LHC are gluinos and electroweakinos. Stops are assumed to be lighter than other squarks but heavier than gluinos, so that the gluino decays via the off-shell top quarks in the three-body decay modes. The pattern of gluino decays will therefore follow the pattern of stop decays discussed in the previous section: Note that for the wino and higgsino cases the above branching ratios assume that m g m χ . For smaller gluino mass, the g → tbW ± mode becomes more favoured by the phase-space. In particular, Br( g → ttW 0 ) = 0 for m g < 2m t + m W . In our Monte Carlo simulation, we calculate the stop and gluino branching ratios using SDecay [31].
The LHC constraints on the gluino simplified model are summarized in Fig. 4, without gravitino (left column) and with gravitino (right column). For the model without gravitino the strongest constraint comes from the ATLAS search for direct stop production atlas_conf_2018_041. It provides stricter limits than atlas_1709_04183 which considers stops produced in gluino decays. However, the latter search only takes into account pre-2018 results.
The limit is the strongest for the bino LSP and extends up to mg ∼ 2200 GeV. For the wino and higgsino case the limit is slightly weaker due to different competing decay modes and extends up to mg ∼ 2000 GeV. In either case the exclusion is up to the LSP mass of ∼ 1000 GeV for the bino and higgsino scenarios, while the limit extends to ∼ 1200 GeV in the wino case.
For the models with the gravitino LSP, the right column of Fig. 4, the shape of excluded region depends on the nature of the NLSP. In the bino case, two searches atlas_conf_2018_041 and atlas_1802_03158 provide comparable limits and gluino with masses up to 2000-2200 GeV are excluded, depending on the NLSP mass. While the atlas_1802_03158 search targets specifically the final states with energetic photons it is using just 36 fb −1 of data. On the other hand atlas_conf_2018_041 does not take the advantage of the additional photons, but includes additional 44 fb −1 of data from year 2017.
The situation is different for the wino NLSP. As we have seen in the stop-wino-gravitino model, the light wino below ∼ 900 GeV is excluded independently of the gluino mass by atlas_1802_03158 exploiting the photon final state. The limit is again due to the electroweak production of chargino, χ + 1 χ − 1 , and chargino-neutralino pairs, χ ± 1 χ 0 1 . However, in the heavy wino region, m W 900 GeV, the most stringent constraint comes from atlas_conf_2018_041 exploiting the hadronic final state with multiple b-jets accompanied by missing energy. This limit only slightly depends on the wino mass, and exclude the gluino up to ∼ 2000 GeV. It is easy to understand why the dedicated GMSB search atlas_1802_03158 gives a weaker exclusion for heavy winos looking at Fig. 2. In the wino scenario, the branching ratio Br( W 0 → γ G) becomes smaller than the Z decay mode, contrary to the bino case.
Similarly to the stop simplified model discussed in the previous section, higgsinos with mass mh 650 GeV, are excluded by cms_1801_03957. For heavier higgsinos, the strongest limit comes again from atlas_conf_2018_041, extending up to 2200 GeV. The limit becomes slightly weaker for higgsinos heavier than ∼ 1500 GeV.

Stop-gluino simplified model
Finally we analyse a model that admits both stops and gluinos. Depending on the mass hierarchy there will be two possible decay chains for gluinos. If m g > m t 1 + m t the gluino will decay as g → t t 1 . Otherwise the gluino will have a three-body decay with the pattern described in the previous section. The limits are summarized in Fig. 5.
In the left column of Fig. 5 we show the limits on the stop-gluino simplified model with electroweakino LSP in the ( g, t 1 ) plane for fixed LSP masses: m B = 400 GeV, m W = 1 TeV and m h = 650 GeV. While there is no experimental limit on the bino [42], the values for winos and higgsinos are motivated by the direct searches, in particular cms_1801_03957. It turns out that regardless of the nature of the LSP the most constraining analysis is atlas_conf_2018_041 and the exclusion is very similar for all types of the electroweakino LSP. The gluinos of mass up to 1800-2100 GeV are excluded for any stop mass, while we do not observe additional exclusion for stops, which is consistent with findings from previous sections.
We now turn to stop-gluino model with the gravitino LSP. For the bino NLSP, as can be seen in top-right panel of Fig. 5, the strongest limit again comes from atlas_1802_03158 (GMSB; γ(≥ 1) + E miss T ) [39] except for the small area where atlas_conf_2018_041 is slightly stronger. The limits consistent with the previous two-dimensional scenarios are obtained when one of the particles is very heavy; m g 2000 GeV for heavy t 1 , and m t 1 1300 GeV for heavy g. The limit is only slightly stronger if both of these particles are light.  search. For very light stop mass, ∼ 1 TeV, atlas_1802_03158 provides the strongest limit and excludes the gluino mass up to 2200 GeV. The lower bound on the stop mass in this scenario is given by 1 TeV, which directly comes from our assumption; m t 1 > m W = 1 TeV. Finally, in the higgsino NLSP case (right-bottom panel of Fig. 5) the most stringent limit is again provided by atlas_conf_2018_041. The limits m g > 2000 GeV and m t 1 > 1300 GeV generally hold. A somewhat stronger exclusion is obtained in the intermediate mass regime for 2000 < m g < 2200 GeV and 1200 < m t 1 < 1500 GeV For the scenarios with the higgsino (N)LSP we also show a contribution to the Higgs mass fine-tuning, ∆ x ≡ 27, over the plane, the largest contribution comes from the gluino, which amounts ∆g 700 due to the gluino mass limit mg 1800 GeV. The gluino contribution to the fine-tuning is particularly large for scenarios with high scale cut-off because of its square dependence, In the light gravitino scenario, the cut-off scale can be lowered, and we take Λ = 100 TeV as an example. Compared to the case without gravitino, the fine-tuning contribution from gluino and stop are significantly relaxed due to much smaller logarithmic factors. One can see that in the region mt 1 1600 GeV and mg 2400 GeV, the fine-tuning contributions from stop and gluino can be less than 30, and a similar contribution will be due to higgsino ∆h ∼ 27 for mh = 650 GeV. We see in general that the LHC constraint on the coloured particle masses is stronger for light gravitino scenario. However, low fine-tuning still favours the light gravitino scenario due to the suppression on the logarithmic factors with smaller cut-off scale.

Conclusions
In this paper we have compared the lower bounds on the gluino and stop masses in two supersymmetric scenarios, one with neutralino as the LSP and the other one with neutralino decaying into gravitino. The analysis of the latest LHC data is based on three simplified models, with bino, wino and higgsino as the (N)LSP and under the assumption that other than those four particles are irrelevant for the collider signatures. For the gravitino LSP, the considered signatures are for prompt neutralino decays, so that obtained bounds apply when gravitinos are very light, in the eV to keV mass range. Such light gravitinos are very interesting cosmologically and typical for gauge mediation models with the messenger mass scale of order O(100 − 1000) TeV. One may expect that low messenger scale ameliorates the fine tuning problem of the Higgs potential.
Generically, the lower limits on the stop and gluino masses obtained with the gravitino LSP are significantly stronger than for the neutralino LSP. This is due to the fact that the gravitino decay as the last step in the chain of decays gives signatures which have very low SM background. One may therefore expect that our conclusions for the comparison of the two scenarios remain qualitatively valid beyond the simplified models studied in this paper.
The fine tuning in the Higgs potential has been discussed in the model with higgsino (N)LSP. In this case all three relevant mass parameters, higgsino, stop and gluino masses, can be fixed at their minimal allowed values following from our analysis. Although with gravitino LSP the lower limits on the gluino and stop masses are stronger than with higgsino LSP for the same higgsino mass, in gauge mediation models the net effect on the fine tuning in the Higgs potential is substantially ameliorated thanks to a very low messenger (cut-off) scale.