Gauge mediation at the LHC: status and prospects

We show that the predictivity of general gauge mediation (GGM) with TeV-scale stops is greatly increased once the Higgs mass constraint is imposed. The most notable results are a strong lower bound on the mass of the gluino and right-handed squarks, and an upper bound on the Higgsino mass. If the μ-parameter is positive, the wino mass is also bounded from above. These constraints relax significantly for high messenger scales and as such long-lived NLSPs are favored in GGM. We identify a small set of most promising topologies for the neutralino/sneutrino NLSP scenarios and estimate the impact of the current bounds and the sensitivity of the high luminosity LHC. The stau, stop and sbottom NLSP scenarios can be robustly excluded at the high luminosity LHC.


Introduction
With the early phase of run II of the LHC under way, it is a particularly good time to review what run I has taught us about weak-scale supersymmetry, and to look ahead about what we can expect to learn from the next data set. The most significant lesson from run I was without any doubt the discovery of a SM-like Higgs with a mass of 125 GeV [1,2]. This has profound implications on the spectrum of the superpartners, as it tends to favor a somewhat higher supersymmetry breaking scale. In the light of this observation, it is perhaps not so surprising that no concrete evidence for supersymmetry has emerged from the 8 TeV data.
Considering the constraints from precision flavor experiments, manifestly flavor-safe forms of supersymmetry breaking are highly motivated. In this category gauge mediated supersymmetry breaking (GMSB) [3][4][5][6][7][8][9][10][11] is perhaps the most elegant and certainly the JHEP01(2017)135 most widely studied paradigm. In previous work [12] we obtained a detailed quantitative understanding of the role of the Higgs mass constraint in the framework of General Gauge Mediation (GGM), which encompasses a wide class of GMSB models [13,14]. In particular, we showed that it is highly non-trivial within GGM to simultaneously achieve electroweak symmetry breaking and a sufficiently large A-term to accommodate a 125 GeV Higgs with LHC-accessible stops. Although it is possible to relax this tension in various extensions of GGM , it is also interesting to take it at face value. In this case requiring m h = 125 GeV induces strong correlations in spectrum of the superpartners, as well as interesting lower and/or upper bounds on some of the sparticle masses. As a result, the GGM framework is surprisingly predictive, even though it contains as many as 8 independent parameters.
In the present paper we make use of the understanding of the GGM parameter space developed in [12] to survey the phenomenology and to investigate the existing constraints and discovery prospects for GGM at the LHC. While our previous paper was primarily intended for an audience of theorists, we hope this paper will be of some interest to experimentalists as well. A general lesson is that accounting for the Higgs mass constraint systematically in top-down SUSY scenarios can drastically modify our expectations for the "first signatures" of SUSY at hadron colliders. This makes full phenomenological surveys of UV motivated scenarios complementary to the bottom up approach based on the full coverage of SUSY decay topologies in terms of simplified models [39,40]. (For simplified models specific to gauge mediation, see [41][42][43][44][45][46][47][48][49].) For instance we will see how the GGM parameter space, with the Higgs mass accounted for, features strong correlations in the spectrum which single out a small set of particularly promising decay topologies.
While the full MSSM parameter space is very rich and complicated, we can get a sense of the broader picture by dissecting its phenomenology in terms of the parameters controlling the colored production. These are the squark masses, here parametrized by the geometric mean of the stop masses M s = √ mt 1 mt 2 , and the gluino soft mass M 3 . (This (over)simplified picture is somewhat analogous to the m 0 and M 1/2 parametrization in the mSUGRA framework and relies on the fact that first and second generation squarks masses are tightly correlated with the stop masses. ) It is possible that the Higgs is as heavy as 125 GeV mostly because the scalars, and in particular the stops, are heavier than 5 TeV. In this case the squarks are clearly inaccessible at the LHC, while the gluino may or may not be accessible. The case with accessible gluinos is often referred to as (mini-)split supersymmetry [50][51][52][53][54][55], and may produce interesting signatures in the form of prompt or delayed gluino decays. If on the other hand, the gauginos are as heavy as the squarks, there is little hope for hints of SUSY at the LHC. The latter scenario includes minimal gauge mediation (MGM), as recently discussed in [56].
In this paper we focus on the collider phenomenology of the remaining part of the parameter space, where the stops are below 5 TeV. This part of parameter space is complementary to the more standard split-SUSY and heavy SUSY scenarios, and was somehow neglected in previous studies of the GGM phenomenology [57][58][59][60][61], essentially because it is difficult to study by taking different 2D slices in terms of the UV parameters. We estimate the current constraints and the reach of the high luminosity LHC (HL-LHC) by mapping JHEP01(2017)135 the viable GGM parameter space into a set of representative simplified models and identify the most relevant simplified topologies. For this purpose we focus on the scenario where the NLSP is long-lived on detector time-scales, which is, as we will see, the most generic case once the Higgs mass is imposed.
The remainder of this paper is organized as follows: in section 2 we first summarize the main results of [12] and then extend these results with a discussion of the expected production cross sections and the various NLSP types. In section 3 we focus on the phenomenology of neutral NLSP. We present the most promising decay topologies with their corresponding branching fractions and estimate the current bounds and the projected bounds at HL-LHC. In section 4 we discuss the phenomenology of the spectra with a charged or colored NLSP and we summarize our results in section 5. Appendix A contains a brief review of some aspects of the NLSP decay in gauge mediation. We reserve some additional results on squeezed spectra for appendix B.

GGM in a nutshell
The defining feature of gauge mediation models is that supersymmetry breaking is communicated via a messenger sector to the MSSM by heavy "messenger" states which are charged under the SM gauge interactions. Throughout this paper, we will refer to the mass scale of the heavy messengers as the "messenger scale" (M mess ). As a consequence, the MSSM fields feel SUSY-breaking only through SM gauge interactions. Since the SM gauge interactions are manifestly flavor blind, the necessary flavor alignment is achieved trivially.
Without specifying the nature of the messengers sector, it was shown in [13] that all gauge mediation models can be captured in a single, very predictive equivalence class. This framework goes by the name of "General Gauge Mediation" (GGM) [13,14] and includes all models where the SM gauge interactions are the only source of communication between the supersymmetry breaking sector and the MSSM. While this is a very general condition, it nevertheless yields a number of strong predictions regarding the SUSY-breaking parameters: • The Higgsino mass parameter µ is set by hand, as a parameter independent from the rest of the soft spectrum.
• The trilinear scalar couplings (A-terms) as well as the B µ term vanish at the messenger scale M mess .
• All soft masses are flavor universal at the messenger scale and satisfy the sum rules 1 also at the messenger scale.
1 These sum rules are typically broken in extensions of GGM with extra yukawa-like interactions involving the MSSM Higgs fields  and also if D-tadpoles are present in the messenger sector [13,62].

JHEP01(2017)135
Most of these relations receive important corrections under renormalization group (RG) running from the messenger scale down to the weak scale, but they nevertheless leave a strong imprint on the low energy spectrum. The full GGM parameter space can then be described in terms of 7 parameters plus the messenger scale (M mess ), which sets the length of the RG-flow. The parametrization we choose here is where we take real gaugino masses and µ, but allow for both positive and negative values all the soft masses and the µ-term. In [14] it was shown that the full parameter space in (2.2) can be spanned by explicit models with weakly coupled messengers. In this paper we do not restrict ourselves to a particular model, but instead deal with the complete parameter space in (2.2). We will organize our presentation around the following, model independent phenomenological features of GGM: • the gravitino is always the LSP 2 • the NLSP decays to the gravitino LPS and a SM state (which depend on the nature of the NLSP) with a decay width which is suppressed by M mess . In appendix A we briefly review the main features of the NLSP decay to the gravitino and a SM state.
In this paper we focus on NLSP masses around or below the TeV scale (within the reach of LHC) and we take two benchmark datasets for high and low messenger scale gauge mediation, with respectively M mess = 10 15 GeV and M mess = 10 7 GeV. In the former case, we assume that Planck-suppressed contributions from gravity mediation are small compared to those from the gauge mediation sector. (For a discussion of gauge-gravity hybrid models, see for instance [64,65].) As we will see, the phenomenology of these two benchmarks differs greatly.

Features of the GGM spectrum
While the parameter space in (2.2) is of course an enormous reduction from the 100+ parameters which characterize the general MSSM, it is still challenging to fully survey this 8 dimensional parameter space. In earlier efforts, this problem has been partially addressed by taking lower dimensional slices [57][58][59][60][61], of which the "Minimal Gauge Mediation" [66] slice is the most well known. While these studies capture a number of generic features, they are insufficient to get a complete picture of the surviving parameter space after the Higgs mass constraint is imposed. In [12] we addressed this deficiency by obtaining both a numerical and a semi-analytic solution of the full 8 dimensional parameter space of GGM, including the Higgs mass constraint. 3 Motivated by the Higgs mass, we restrict ourselves to the regime where tan β 10. In this regime the spectrum is fairly insensitive to tan β,

JHEP01(2017)135
which we fix to tan β = 20. Also solving for m Z and m h then reduces the parameter space from 8 to 5 dimensions. Of these 5 remaining parameters M mess and M 1 can be scanned coarsely without sacrificing a smooth interpolation, such that a manageable 3 dimensional volume remains to be mapped out carefully. The most important features of the solution are most easily understood in terms of a set of approximate relations between the soft parameters at the weak scale [12]. In the large tan β limit these relations are These relations are direct consequences of flavor universality, the electroweak symmetry breaking (EWSB) conditions and the GGM UV sum-rules in (2.1), and are therefore independent of the messenger scale. A number of phenomenological features immediately follow: • 1st/2nd generation Q and U squarks are always heavier of the 3rd generation squarks, although the mass splitting can be small if the left-handed slepton and the Higgsino are both light.
• D squarks are always heavier than the quadratic mean of m Q 3 and m U 3 .
• The right-handed sleptons E are always heavier than the left-handed sleptons L if m Q 3 < m U 3 . Reversely, for a relatively light Higgsino, the m U 3 m Q 3 implies that the left-handed sleptons are heavier than the right-handed sleptons.
• If both the left-handed slepton and the Higgsino are light, so must be the pseudoscalar Higgs A 0 .
In addition to these simple, M mess independent relations, there are a number of features which do depend on the choice of M mess in a fundamental way. Ultimately, all these features can be traced back to the Higgs mass constraint: for stop masses in the few TeV range, the Higgs in the MSSM can only be sufficiently heavy if the top A-term (A t ) is large. Since GGM predicts A t = 0 at the messenger scale, the A-term must be generated in the RG-running, which requires a heavy gluino and/or high messenger scale [72]. In [12] we showed how this arrangement has the following indirect effects: • The RG-equation for m 2 Hu depends strongly on both A t and the stop masses, where the net result is that the radiative corrections tend to drive m 2 Hu upwards in the regime of interest. In combination with the GGM sum-rule in (2.1), this produces a tension JHEP01(2017)135 between proper EWSB and the absence of slepton tachyons. In practice this results in a strong lower bound on the stop masses, as shown in figure 1, where we projected the full parameter space on the plane of the IR stop soft masses (m Q 3 vs m U 3 ). The bound is asymmetric in the stop mass plane due to the m E relation in (2.3).
• Since the horizontal boundary in figure 1 is due to the conflict between achieving EWSB and avoiding a slepton tachyon, both the Higgsino and the sleptons are constrained to be light nearby the boundary. 4 For m Q 3 < m U 3 (m Q 3 > m U 3 ) the lightest slepton is left(right)-handed, again due to (2.3). The sleptons, being scalars, decouple relatively quickly, however the Higgsino must remain light in a substantial part of the parameter space, as shown in figure 2.
• The milder lower bound on m Q 3 has a different origin. Since the gluino must be heavy to ensure a sufficiently large A t , as shown in figure 2, the small tree-level stop masses receive a large gluino-induced threshold correction [73], which tends to drive the stop tachyonic. In this part of parameter space there is effectively a little hierarchy problem between the stop and the gluino. This effect is especially pronounced for M mess = 10 7 GeV, since the gluino must be much heavier in this case.
• There is an interesting correlation between the sign of µ and the gaugino spectrum, which greatly affects the phenomenology. In particular, we find that M 2 is bounded from above (below) for µ > 0 (µ < 0). (These bounds can be traced back to a combination of the Higgs mass constraint, the EWSB conditions and the requirement that  for µ < 0 the gluino mass tends to be lighter than for µ > 0.
• Finally, if both wino and the Higgsino are light, they induce a small but positive one loop contribution to m h [73]. The impact of this correction is especially important for M mess = 10 7 GeV, since a large A-term is more difficult to achieve in this case. As discussed in the previous bullet point, this spectrum is only possible if µ > 0. As a consequence there is no parameter space for LHC-accessible squarks with µ < 0 and M mess = 10 7 GeV.

The LEP bound
Before discussing the different GGM scenarios for the LHC and their corresponding limits, we take a moment to stress the importance of the limits set already by LEP. Due to the nature of the experiments, searches for SUSY at LHC are necessarily less inclusive and more complex than the corresponding searches at LEP. For the GGM parameter space, the relevant limits are those on the chargino's, as the sleptons are almost always heavier than the kinematic reach of LEP. Even though the LEP limit on the lightest chargino is more inclusive than most LHC searches, there is still some dependence on the full chargino/neutralino spectrum. To account for this, the LEP experiments analysed a wino-like and a Higgsino-like benchmark model and found limits between roughly 90 and 105 GeV [74,75], depending on the model and theχ ± 1 -χ 0 1 mass splitting. Throughout this paper we impose the conservative limit of mχ± 1 90 GeV. We find that the constraints on the parameter space are not significantly modified by a ∼ 10 GeV shift in this limit.
Surprisingly, these LEP limits provide a fairly strong constraint on the GGM parameter space, despite their very limited mass reach. This can best be seen from figure 1, where the LEP limits on the charginos exclude stop masses far outside the reach of LEP, and sometimes even outside the reach of the LHC. The reason lays in the upper bounds GGM enforces on the Higgsino and wino masses, as discussed in the previous section and as depicted in figure 2. This leads to a particularly strong constraint for M mess = 10 7 GeV, as in this case both the wino and Higgsino must be light in much of the parameter space, to accommodate the Higgs mass constraint. For a thorough proof of this effect we again refer to [12].

SUSY cross section
The constraints on the GGM spectrum summarized in the discussion above have a strong impact on the dominant production channels in GGM and consequently on the asymptotic reach of HL-LHC. In fig 2 we show the minimal cross section at the 14 TeV LHC, separately for colored and electroweak production (indicated as σ c and σ ew in the plot), where we marginalized over the full parameter space. We computed the NLO cross section for colored sparticles with prospino2 [76], while for the electroweakinos we compute LO cross section with pythia 8 [77,78]. Since the gluino mass is above 2.5 TeV in most of the parameter space, the colored cross section is predominantly controlled by the squarks and is mostly independent of the sign of µ and M mess . However comparing the two signs of µ for M mess = 10 15 GeV, we see that the colored cross section is somewhat larger for µ < 0 because of the lighter gluino mass with respect to µ > 0. The strong correlation between the squark masses of different families encoded in (2.3) implies that m Q 3 and m U 3 provide a good parametrization of the colored cross section. This is clearly borne out in the left-hand column of figure 2. For M mess = 10 15 GeV, the lower bound on m U 3 implies that the cross section is largest for small m Q 3 , where it is dominated by the left-handed squarksq L ,t L andb L . For M mess = 10 7 GeV this continues to be true, however the strong lower bounds on both m Q 3 and m U 3 result in a very suppressed cross section for colored production after the LEP bounds are imposed. As expected, the minimal cross section decreases as the squarks become heavier. For µ < 0 and m Q 3 ∼ m U 3 4.5 TeV the minimal colored cross section however increases again, since gluino mass is bounded from above by the Higgs mass constraint: in this part of the parameter space, a larger gluino mass would lead to a large enough A-term to overshoot the Higgs mass. To some extent this is an artifact of our parametrization where we fixed tan β = 20, and this effect can be compensated by lowering tan β as the stops become heavier. In this case the minimal colored cross section would continue to decrease as expected.
In contrast with the colored cross section, the electroweak cross section does depend strongly on the sign of µ and M mess (left-hand column of figure 2). For µ < 0 the wino is decoupled and the minimal electroweak cross section tracks the maximal Higgsino mass (for example for µ ∼ 1 TeV we see that σ ew ∼ 1 fb). More interestingly, for µ > 0 the wino and the Higgsino can both be light and both contribute to the electroweak cross section. This plays a particularly important role for the collider phenomenology of M mess = 10 7 GeV, where colored production is suppressed.

The NLSP in GGM
After giving an overview of the dominant production channels in the stop mass plane, it is time to isolate the characteristic decay topologies of GGM. This task is generally challenging since searches at the LHC depend strongly on multiple features of the spectrum. A natural organizing principle in gauge mediation is given by the fact that the gravitino is always the LSP, such that the nature of the NLSP is of particular importance for the collider phenomenology.

JHEP01(2017)135
While the NLSP is guaranteed to decay to the gravitino plus a standard model state, the time scale of this decay may vary over many orders of magnitude. In particular, the NLSP decay tends to occur outside the detector unless the messenger scale is low (i.e. M mess 10 7 GeV), in which case the decay may occur prompt, displaced or outside the detector. In this paper we focus on the phenomenology of long-lived NLSPs: once the Higgs mass is imposed, this is the most generic scenario, under some fairly broad conditions on the nature of the SUSY-breaking sector (see appendix A).
We assume R-parity conservation, such that the NLSP is produced in pairs, either through direct production or through a cascade decay. The associated signatures strongly depend on the SM quantum numbers of the NLSP itself, specifically whether the NLSP is charged under the strong and/or electromagnetic force: • A neutral NLSP (neutralino or sneutrino) escapes the detector without leaving any track and is only seen as missing transverse energy. Its direct production at a hadron collider can generically only be bounded by mono-jet searches [79,80]. (Although sometimes more efficient search strategies can be designed if the NSLP is a member of a quasi-degenerate electroweak multiplet, see below.) The strongest bounds therefore almost always come from cascade decays of other, heavier superpartners.
• On the other hand, a charged and/or colored NLSP (stop, sbottom, stau or gluino) produces a spectacular signature in the form of a pair of highly ionizing tracks, and searches for this final state are nearly background free [81,82]. Extra activity in the event from a cascade is therefore not needed to further reduce the background, although a cascade decay from a colored state may still greatly increase the cross section in the case of a stau NLSP.
In figure 3 we show how requiring a particular NLSP type correlates with the stop masses. This correlation, together with the total SUSY cross section shown in figure 2, already provides a rough idea of the reach of LHC at 14 TeV for a given NLSP. Before discussing this in more detail, we stress that the projection of the possible NLSP types on the m Q 3 vs m U 3 plane in figure 3 only shows whether a given NLSP type is allowed for certain values of the squark masses, which roughly correspond to the colored cross section (see figure 2). Figure 3 should therefore not be interpreted as a measure of the relative abundance of a given NLSP in the full parameter space. Since certain UV completions may very well favor one NLSP type over another, we will treat every allowed NLSP in GGM on the same footing in this paper.
The neutralino NLSP is possible for any value of the stop masses. In general the lightest neutralino is a mixture between the neutral components of the bino, wino and Higgsino, but in most of the parameter space these masses are relatively split, such that the NLSP can usually be thought of as either mostly bino, wino or Higgsino. Depending on the relative hierarchy among M 1 , M 2 and µ a neutralino NLSP can be: • Bino-like: if |M 1 | < |µ|, |M 2 |. The bino NLSP plays a special role in our discussion, as M 1 only enters in the RG-equations only through terms suppressed by ∼ g 2 1 . As a result, M 1 is an approximate flat direction in the parameter space which does not JHEP01(2017)135 significantly affect the rest of the spectrum as long as M 1 is in the TeV-range. Since there is no direct limit on a bino NLSP (not even from LEP if µ is heavy enough [83]) a bino NLSP is allowed everywhere in the parameter space. Note however that near the lower bound on m U 3 , µ and M 2 (if µ > 0) are bounded from above, such that NLSP is necessarily somewhat mixed, even if |M 1 | < |µ|, |M 2 |.
• Higgsino-like: if |µ| < |M 1,2 |. A Higgsino-like NLSP is possible for all values of the stop masses. Its mass is bounded from above as was shown in figure 2, and bounded from below by the LEP bound of mC 1 > 103.5 GeV. The Higgsino multiplet consist of two neutral states and one charged state which are typically split between 1 and 10 GeV for M 1 in the TeV range. The charged state decays promptly to the lightest JHEP01(2017)135 neutral state but the standard model particles in this decay are generally too soft to benefit from dedicated search strategies at LHC [84,85].
• Wino-like: if |M 2 | < |µ|, |M 1 |. As was mentioned earlier and shown in figure 2, for µ < 0 |M 2 | is bounded from below and the wino can only be the NLSP if the stop masses are both above 3 TeV and M 2 is larger than 2 TeV. Conversely, for µ > 0 the wino may be the NLSP in nearly all of the stop mass plane. Also in this case there is a strong upper bound on M 2 for stops close to the boundary of the viable region (see figure 2). This bound is particularly strong for M mess = 10 7 GeV. The wino multiplet consists of a charged (χ ± 1 ) and a neutral state (χ 0 1 ), typically split by ∼ a few 100 MeV to a few GeV, depending on whether the Higgsino is nearby in the spectrum. χ + 1 decays to χ 0 1 by emitting a pion, which is too soft to be a useful observable, as was the case for the Higgsino NLSP. An important exception can arise if the splitting is smaller than ∼ 200 MeV, in which case the decay is displaced and χ ± 1 is seen as a disappearing track [86,87].
For the Higgsino and wino NLSP scenarios we always assume a decoupled bino for simplicity. If the bino mass were instead in the proximity of the Higgsino or wino NLSP mass, the splitting between the components of the multiplet would typically increase.
A slepton NLSP can be either mostly left-handed (L) or right-handed (E). The third generation sleptons are typically slightly lighter than the first two generations, such that a slepton NSLP is usually a right-handed stau or a left-handed snutau. 5 From the GGM sum-rules (2.1) we thus expect aτ R NLSP for m Q 3 > m U 3 andν τ NLSP for m Q 3 < m U 3 , which is very clearly confirmed by our numerical solution in figure 3.
is the NLSP, the final state consists out of a pair of highly ionizing charged tracks [82,89]. In the region where this configuration is possible, the stop and the gluino are however fairly heavy, and the dominant contribution to the cross section comes from electroweak states. In the region where the left-handed slepton is the lightest slepton (m Q 3 m U 3 ) the NLSP is always the neutralν τ since the charged states of the L doublet are always pushed upwards by a small but irreducible D-term contribution of the form The splitting between the slepton and the sneutrino is usually around 10 − 30 GeV and it decreases with increasing sneutrino mass. Direct production of slepton pairs will therefore lead to soft leptons in association with missing energy. The snutau NLSP behaves mostly like a neutralino NLSP, and the only bound on its direct production comes from LEP [75]. However, since some fine tuning is involved in order to get the sneutrino light, it is very atypical to get it within the reach of LEP experiments (practically speaking its mass never goes below ∼ 150 GeV in our dataset). Fortunately, for a snutau NLSP the squarks may be light enough to make this scenario testable at LHC.

JHEP01(2017)135
A stop or sbottom NLSP does not occur for M mess = 10 7 GeV, barring very finetuned regions in the parameter space. They can however be realized for M mess = 10 15 GeV if m Q 3 m U 3 , but even in this case its mass is still always smaller than roughly 1.5 TeV (again see 3). Since the stop/sbottom NLSP must always be left-handed, the second colored state in the SU(2) multiplet is always nearby in mass.
Finally, the gluino can always be the NLSP in the (mini-)split regime where the scalars are heavy. (In this case no large A-term and therefore no large M 3 is needed, since the Higgs mass comes predominantly from the mass of the stops.) If the gluino is the NLSP, the squarks are never accessible at the LHC. Since this case only arises in a very small part of our parameter space and is studied elsewhere already [50][51][52][53][54][55], we do not discuss it in any more detail.

Neutral NLSPs
In this section we present the most relevant simplified topologies for the neutral NLSPs in GGM, and estimate the current constraints and the future reach of the High Luminosity LHC (HL-LHC) with 3000 fb −1 . For this purpose we take the neutral NLSPs to be longlived on detector length scales, which is preferred for M mess 10 7 GeV under some fairly broad assumptions about the nature of the SUSY-breaking sector (see appendix A). For M mess 10 7 GeV, the neutral NLSP may decay either prompt, displaced or outside the detector. A visible NLSP decay results in two additional SM bosons in the final state, which typically strengthens the limits on a neutralino NLSP. (For a discussion of prompt decays of neutralino NLSP's we refer to [43,46,48] and to [44,90] for an analysis of displaced decays.) On the other hand, a sneutrino NLSP decays fully invisibly to a neutrino and a gravitino, and in this case the collider phenomenology therefore does not depend on the NLSP lifetime.
The collider phenomenology depends on whether the NLSP is a sneutrino or a neutralino with a dominant wino, Higgsino or bino component. Moreover many features of the spectrum drastically change with the choice of M mess and the sign(µ), as we discussed in the previous section. We therefore organize our discussion around the NLSP type, M mess and sign(µ). For each case we discuss the most important simplified topologies and present the mass range for each particle in the SUSY spectrum which is compatible with the Higgs mass constraint, the GGM boundary conditions, proper EWSB and LEP constraints. We subsequently overlay our estimates for the constraints from the existing LHC data, as well as an estimate of the future sensitivity of the HL-LHC.
It is always possible to trivially relax all collider constraints by simply decoupling the whole superpartner spectrum, and for concreteness we therefore require M s ≡ √ mt 1 mt 2 < 4 TeV. This allows us to focus on the "low" M s phenomenology of gauge mediation, where at least some of the scalar superpartners are within the reach of LHC and where the gluino does not sensibly contribute to colored sparticle production. For M s > 4 TeV the squarks are decoupled but the gluino is not, and the phenomenology strongly resembles that of (mini)-split supersymmetry [50][51][52][53][54][55]. The most important collider bounds for long-lived neutral NLSPs come from SUSY cascades initiated by colored or electroweak superpartners. In our treatment of the collider bounds we neglect the efficiency drop of SUSY JHEP01(2017)135 searches for squeezed spectra (i.e. when the NLSP mass is nearby the mass of the sparticle responsible for the production). Since we effectively marginalize over the NLSP mass, including squeezing would imply that the limits would be effectively dominated by the squeezed spectra. Such spectra only occur in a small part of the parameter space, and this would therefore result in an overly pessimistic picture of the constraints on the parameter space as a whole. For completeness, we give an example of how much squeezed spectra can affect the constraint on the GGM parameter space with Higgsino NLSP in appendix B.

Wino NLSP
The structure of the GGM spectrum with a wino NLSP is summarized in figure 4, where the various bands indicate the viable interval for each mass parameter. The most important feature is that a wino NLSP within the LHC reach is only possible if µ > 0. (See also figure 2 and figure 3.) In the remainder of this section we therefore only focus on µ > 0.
Having fixed M 1 to be decoupled, the splitting between charged (χ ± 1 ) and neutral (χ 0 1 ) components of the wino multiplet depends on the Higgsino mass and may be sufficiently small for the disappearing track searches to become relevant [86,87]. If the splitting betweenχ + 1 andχ 0 1 is too large for the decay to be displaced, the decay product is typically too soft to provide a good enough handle to separate the SUSY signal from the SM background. Bothχ + 1 andχ 0 1 can then effectively be treated as MET. In this case the bounds on the wino NLSP scenario thus come exclusively from cascade decays, as we will detail in the next subsection.
After the collider bounds are imposed, figure 4 reveals two particularly interesting correlations in the spectrum. Firstly, for HL-LHC there are very mild upper bounds on JHEP01(2017)135 the wino and gluino masses. This is a result of the anti-correlation between M 2 and M 3 , discussed in section 2.2, which indicates that a lower bound on M 2 implies an upper bound on M 3 and vica versa. A second interesting feature is that the left-handed slepton and A 0 are both forced to be heavy, which holds regardless of the collider constraints. This is again a consequence of the Higgs mass constraint and the EWSB conditions [12]. The LHC further constrains these particles indirectly through the limit on the wino mass. Finally, in [12] it was shown that for µ > 0, µ increases if M 2 increases. This correlation results in a strong, indirect lower bound on M 2 from the LHC bound on µ (see next subsection).
In summary, we see that once future HL-LHC bounds are accounted for, the only states below one TeV are electroweakinos, controlled by µ, M 2 and the right-handed sleptons.

Constraints and simplified topologies
We first discuss the bounds on direct wino production. The disappearing track search is relevant if mC 1 − mÑ 1 200 MeV [86,87], which corresponds to µ 500 GeV. The electroweak production of one or two disappearing tracks is therefore always dominated by direct wino production. The bound on disappearing tracks from the 8 TeV (HL-LHC) data on a directly produced wino ranges between 100 GeV (300 GeV) for mC 1 − mÑ 1 ∼ 200 MeV up to 500 GeV (1300 GeV) for mC 1 − mÑ 1 ∼ 140 MeV. For the HL-LHC reach we obtained a naive bound by rescaling the projected bound obtained in [91,92] to different values of mC 1 −mÑ 1 by assuming the same expected sensitivity as in [86]. When imposing the bound on disappearing tracks we are neglecting a possible enhancement of the signal strength coming from left-handed squarks decaying to the lightest chargino either directly or through a SUSY cascade involving the Higgsino. Recasting and optimizing disappearing tracks searches for SUSY spectra with a non-negligible colored production would not qualitatively modify our conclusions regarding the currently excluded GGM spectra in figure 4 but it could provide one of the most sensitive probes of the wino NLSP scenario at the HL-LHC.
If the splitting mC 1 − mÑ 1 200 MeV, the only bound on direct wino production comes from mono-jet and VBF+MET searches, which are very weak even at the HL-LHC [91,93,94]. The sparticles higher up in spectrum thus dominate the collider limits, and it is therefore most convenient to think about the phenomenology in terms of simplified topologies. Both the gluino and the right-handed squarks are decoupled as far as the current LHC data is concerned. Therefore the most promising production channels are the Higgsino, the left-handed squarks and the left-handed stop/sbottom. The sleptons tend to be relatively heavy, except very near the boundary of the viable region. Given that their production cross section is also very small, their direct production typically does not play a role in the constraints.
The four most relevant simplified topologies are shown in figure 5: in the first topology all colored states are decoupled, such that the production cross section is provided by the Higgsino (figure 5a). The decay to the wino NLSP occurs through the emission of gauge and Higgs bosons, of which W Z+MET final state is the most constraining. In figure 6 we show the effective rate for this process, where we combined all production channels, weighted by the branching ratio to W Z+MET. Also indicated on figure 6 is the current limit, as well as the projected limit for the HL-LHC, which we rescaled from [95] and [96,97] respectively. Throughout this paper, the various branching ratios were computed with SUSY-HIT [98].
The second simplified topology (figure 5b) parametrizes the case where the left-handed squarks are accessible. The branching ratio to jets+MET is 100%, but we rescale the bound in [99] to account for the fact that only theq L are accessible. 6 In this rescaling we made use of the NLO squark production cross section as computed by prospino2 [76], where we account for the (mild) dependence of the squark cross section on the gluino mass (see figure 2). For gluino masses in the middle of the allowed range, this roughly amounts to mq L 1100 GeV (mq L 1050 GeV) for M mess = 10 15 GeV (M mess = 10 7 GeV). At the HL-LHC, theũ R andd R may be accessible as well, although they are most likely not degenerate withq L . In this case we estimate the limits by requiring that the total squark cross section is smaller than the projected limit in [96]. A priori there is a second squark-initiated topology, where the Higgsino is in between the squarks and the NLSP. However because of the small coupling of the lowest generation squarks to the Higgsino, the branching ratio of the cascade decay via the Higgsino is negligible, and we do not need to consider this case separately.
The final production mode is through the left-handed stop and sbottom. In this case the presence of the Higgsino in the middle of the spectrum does make a qualitative difference. We first consider the case where the Higgsino decouples, as in figure 5c. Even though the A-term is large, the stop mixing is still relatively small in the region where the lightest stop and sbottom are accessible with the current data. The reason is again the lower bound on the right-handed squarks, which enforces a sizable mass splitting between the stop gauge eigenstates. The left-handed stop and sbottom are therefore degenerate to good approximation, and as such their production cross sections are approximately equal. To establish a bound, we can therefore treat them as a single state with branching ratios which are the average of those of the stop and the sbottom. The branching ratios to do not depend strongly on the mass parameters and are roughly Br(bb + MET) ≈ 0.31, Br(tb + MET) ≈ 6 In principle there are a number of important caveats to this approach which should be mentioned [100].
First, the efficiencies of the searches drop dramatically for squark masses below ∼ 700 GeV. For the righthanded squarks this region is already excluded by the constraints in figure 1, while for the left-handed squarks this region is excluded by constraints on the left-handed stop and/or sbottom. Secondly, the difference in the parton luminosity functions implies that the cross sections of theũ andd are significantly higher than that forc ands. However since the squark spectra in gauge mediation are flavor degenerate in the first two generations, this effect is unimportant.   point a set of analogous simplified topologies with the right-handed squarks will be come relevant, in addition to those discussed in this section.

Higgsino NLSP
For the Higgsino NLSP, summarized in figure 8, the gluino and the right-handed squarks are also too heavy to be significantly constrained by the current data. However in contrast with the wino NLSP, it is possible to probe the µ < 0 branch with current data. On the other hand, the splitting between the components of the Higgsino multiplet is always large enough to allow for prompt decays of the two heaviest states of the multiplet. The visible decay products of this decay are generally still too soft to be efficiently discriminated from background, such that the strongest limits come from cascade decays of heavier sparticles. We discuss the most relevant topologies in the next section.
As for the wino NLSP, there are a number of GGM-specific correlations in the spectrum. Firstly, the HL-LHC can also set an interesting indirect lower bound on the Higgsino mass of ∼ 300 GeV for the case of µ < 0, as is apparent from the summary plot in figure 8. This bound is sensiblity stronger than the direct bound from monojet+MET [91,93], which remains applicable for µ > 0. The indirect bound is a consequence of the bound on the gluino and a particular correlation in the GGM parameter space which is unique to µ < 0 [12]. (It is also worth noting that only for µ < 0 the HL-LHC projection of the reach on gluino production (mg 2900 GeV [96]) plays a role in constraining GGM.) Secondly, figure 8 reveals the same anti-correlation between M 2 and M 3 as was present for the wino NLSP. A third example of an important GGM-specific correlation is apparent in the lower bounds on the left-handed sleptons. For µ < 0, M 2 is always in the multi-TeV range, and keeping the sleptons light therefore requires fine tuning. For µ > 0, m L is correlated  with µ such that the LEP bound on the chargino results in an indirect lower bound on the sleptons. This correlation is strongest for high messenger scales. For more details, we again refer to [12]. The lower bound on the left-handed slepton masses subsequently induces a lower bound on m A 0 (see eq. (2.3)).
In summary, we see that the current data do not lead to a strong lower bound on the SUSY spectrum, but the HL-LHC will be able to push every SUSY particle above 1 TeV except the Higgsino and the right-handed sleptons, regardless of the choice for M mess and sign(µ). In addition, for µ > 0 the wino can be below 1 TeV. If M mess is also low, the left-handed sleptons and the CP-odd Higgs could be relatively light as well. This feature is most pronounced for the Higgsino and sneutrino NLSP's due to the relation in eq. (2.3). The pseudo-scalar and slepton masses also tend to be lower for low messenger scales, which can be understood from combining the RG-running with the EWSB conditions, see eq. (2.5) in [12].

Simplified topologies
As for the wino NLSP the bounds from monojet+MET searches are very weak, and are expected to probe Higgsino only up to ∼ 200 GeV at HL-LHC [91,93]. The remaining simplified topologies for the wino NLSP also apply for the Higgsino NLSP, with the roles of the wino and Higgsino interchanged. (See figures 9a, 9b, 9c and 9d.) Forq L -Higgsino topology the analysis is identical as for the corresponding case with a wino NLSP and we assume the same limits. The decay in thet L /b L -Higgsino topology is dominated by the top yukawa and the final states are therefore predominantly top-rich. In particular, for tan β = 20 the branching ratio to tt+MET is roughly 90% for botht L andb L . Rescaling the (projected) limits in [99] and [103]   For the wino-Higgsino topology, which is relevant only for µ > 0, the branching ratios are modified with respect to the Higgsino-wino topology in figure 5a. This leads to a different effective rate for W Z+MET final state, 7 as shown in figure 10. Similarly, the branching ratios for thet L /b L -wino-Higgsino are modified with respect to its analogue for with a wino NLSP (see figure 11). Finally, for the Higgsino NLSP there is aq L -wino-Higgsino topology (figure 9e), whose analogue was not needed for the wino NLSP. In this topology, the cascade decay via the wino always dominates, unless the wino is very close in mass to the squarks (see figure 12). Also here, since the wino is always heavy for µ < 0, topologies 9a, 9d and 9e only apply if µ > 0. As for the wino NLSP, we do not attempt a proper recasting of the existing limits in terms of the cascade topologies of the colored sparticles, and instead we take the same squarks limits as for the wino NLSP. Especially for the squarks, the limits are expected to be very sensitive to the wino mass, as was shown in [104] for a bino NLSP. However when we marginalize over the full parameter space, as in figure 8, our crude approximation should give nevertheless a reasonable idea of the overall strength and relative importance of the various limits.

Bino NLSP
The bino NLSP scenario, summarized in figure 13, is considerably simpler than the wino and Higgsino NSLP: while the wino and Higgsino masses strongly correlate with the rest of the spectrum, the bino mass has essentially no impact. This also implies that the parameter space for the bino NLSP is substantially larger than for the Higgsino and the wino NLSP, since it is always possible to put the bino at the bottom of the spectrum. More importantly, the bottom of the spectrum then consists out of single neutralino, rather than a quasi degenerate multiplet of several electroweakinos. This greatly simplifies the analysis of the branching ratios in most of the simplified topologies. For this reason simplified 7 It is possible that the sleptons sit below the wino, which enhances the sensitivity, since a cascade decay with additional leptons in the final state becomes possible. This however happens only in a relatively small part of the parameter space, and we do not consider this possibility separately.  topologies with bino (N)LSP have been a popular method to parametrize the experimental limits, and for most topologies considered here a direct limit is available.

JHEP01(2017)135
Again neglecting the sleptons, the relevant simplified topology for electroweak production simply consists of the wino and/or Higgsino in addition to the bino NLSP. ATLAS has set a limit directly on this topology in the M 2 -µ plane [95], where the strongest limit again comes from the W Z+MET final state, as shown in figure 14. In figure 14 we also estimate the reach for the HL-LHC, where we rescale the projection for the wino-bino simplified topology in [96] to also allow for the possibility of an accessible Higgsino.
For colored production the topologies are the same as those for the Higgsino NLSP, upon replacing the Higgsino multiplet with the bino. In addition to the trivialq L -bino topology, theq L -wino-bino topology is relevant as well. To estimate the current and projected bounds, we follow the same procedure as for wino and Higgsino NLSP. For thẽ q L -wino-bino model the limit in principle varies significantly with the wino mass [104], but when we marginalize over the parameter space in figure 13 our simple minded procedure still provides a reasonable idea of the impact of the bound on the GGM parameter space.   [95] and estimated reach (red) on the wino-Higgsino-bino simplified topology with M 1 = 50 GeV. The projected limit was rescaled from [96].

JHEP01(2017)135
The limits on thet L /b L -bino topology can be read off directly from the stop and sbottom searches. The strongest (projected) limit is mt L ,b L 900 GeV [101,102] (mt L ,b L 1500 GeV [96]), where we again took the left-handed stop and sbottom to be quasi degenerate. In a significant part of the parameter space it is moreover possible that both the wino and the Higgsino are in between the stop/sbottom and the bino NLSP. In this case a very complicated, multi-step cascade is possible. Finally, for the HL-LHC reach of the gluino, we employ the same estimate as for the Higgsino NLSP.

Sneutrino NLSP
The sneutrino NLSP behaves qualitatively similar to the neutralino NLSPs, except that the decay chains tend to be lepton and/or tau rich, which can result in a very rich and interesting phenomenology [105][106][107][108][109]. Although the snutau is typically slightly lighter than the remaining sneutrino's, we can take all three flavors as degenerate for the purpose of this discussion. The mass splitting betweenτ L andν τ is typically 5 GeV or less, which means that any leptons produced in theτ L → W * ν τ decay tend to be fairly soft. This is especially so for directτ + Lτ − L production, which is usually produced on or near threshold. Direct production of the sneutrino NLSP is therefore notoriously difficult to constrain at the LHC, and the presence of heavier states with larger cross sections is therefore essential to probe this scenario efficiently.
The sneutrino NLSP can only occur when m Q 3 < m U 3 , as shown in figure 3. Moreover for M mess = 10 15 GeV with µ < 0 a large amount of tuning is needed to achieve a sneutrino NLSP and EWSB simultaneously, and a specialized study is therefore required to make quantitatively accurate statements about this scenario. Since this concerns a very small part of the GGM parameter space, we do not attempt this here. For M mess = 10 15 GeV with µ > 0, a sneutrino NLSP is more likely, but even in this case the remainder of the spectrum is very constrained by the EWSB conditions and the GGM sum rules: aside from the m Q 3 < m U 3 requirement, we find that the wino mass is always in the (multi-)TeV range, as shown in figure 15. This means that it is irrelevant for direct production and for on-shell cascade decays, such that the relevant simplified topologies are just those shown in figure 16.
For M mess = 10 7 GeV the spectrum is even more restricted, and the only relevant topology is the Higgsino-sneutrino topology in figure 16a. For this topology, the most promising channel is 2τ +MET, which has branching fraction 1/4. However even in this channel the SM background from W W production is very large, and there is currently no limit from the LHC once the LEP bound onν is accounted for [110,111]. (It is however possible that mild bounds could be obtained by accounting for possible soft leptons in theτ L decay.) In particular for low scale gauge mediation we thus expect that it will be very difficult for the LHC to significantly constrain the sneutrino NLSP scenario in this way. On the other hand, the CP-odd Higgs tends to be light if the Higgsino is light (see eq. (2.3)). As for the Higgsino NLSP, this effect is more pronounced for low messenger scales and therefore provides an interesting, complementary constraint of m A 0 500 GeV from current data [112,113]. For the HL-LHC we obtain m A 0 900 GeV by rescaling the projection in [114].
For M mess = 10 15 GeV with µ > 0, the left-handed squarks, stop and sbottom can be accessible with current data, which could provide a much needed boost to the signal cross section. An important difference with the neutralino NSLP is that the squarks, stop and sbottom cannot directly decay to the NLSP, but must go through an intermediate electroweakino. This intermediate electroweakino can be either an on-shell Higgsino (figures 16d and 16e) or an off-shell wino and/or Higgsino (figures 16b and 16c). The on-shell cascade takes priority whenever it is available.
Since the wino is heavy, the on-shell decays must go through the Higgsino and these decay chains are therefore characterized by (nu)tau-rich final states with missing energy JHEP01(2017)135 in addition to tops, results or jets. A particularly interesting mode ist LtL to a pair of on-shell Higgsinos, which predominantly decays in ttτ ± τ ∓ +MET and ttτ ± τ ± +MET with roughly 40% branching ratio each. This decay topology is currently not explicitly covered, although it should be possible to regain some sensitivity by recasting the multi-lepton searches or the 2b + τ + τ − + MET search [115]. Since the latter search is inclusive as far as the number of jets is concerned, we can obtain mt L 500 GeV as a naive estimate of the bound by rescaling the observed bound to account for the branching ratio 2t + τ + τ − , with fully hadronic top decays. The off-shell, three-body squark decay in figure 16c prefers an off-shell wino and therefore tends to be lepton flavor democratic. For the stop/sbottom initiated three-body decay in figure 16b, the branching ratios are sensitive to the relative masses of the wino and the Higgsino. In this paper we do not attempt to extract an approximate bound for this case.
In summary, the phenomenology oft L /b L /q L is very rich and complicated: some channels are currently not or poorly covered, while others require detailed recasting [116]. To get very rough, qualitative sense of what impact of such a program would be on the GGM parameter space, we imposed the same (projected) squark bounds as for the Higgsino NLSP in figure 15, in addition our naive rescaling of the 2b + τ + τ − + MET search. The true bounds are likely to be stronger due to the additional leptons and/or tau's in most of the channels.

Charged/colored NLSP
We again assume the NLSP to be long-lived, something which is preferred for M mess 10 7 GeV. For M mess 10 7 GeV prompt or displaced decays are possible, however as we will see, this case does not allow for a stop/sbottom NLSP. A promptly decaying stau NLSP is possible, and in this case the limits are very weak since colored production is not accessible for low M mess . (See for [45,48] for a discussion of electroweak production with a promptly decaying stau NLSP.) For a discussion of the displaced decay of a stau NLSP we refer to [117,118]. For a long-lived stau, stop and sbottom NLSP's the situation is considerably simpler than for a neutral NLSP. While there can be some efficiency loss for cascades with particularly large hierarchies or with heavy NLSPs [49], generally the details of the spectrum are much less important than for searches which rely on MET. The searches for heavy stable charged particles (HSCP's) can therefore be interpreted as a rough bound on total inclusive sparticle cross section. The current bounds and projected limits are summarized in table 1. To estimate the current and future bounds on the parameter space we compute the total SUSY cross section by adding up the individual pair production cross sections of each sparticle and compare this with the corresponding limit in table 1.

Slepton NLSP
The stau NLSP scenario can be realized only if m Q 3 > m U 3 , which implies that colored production is very suppressed due to the strong lower bound on the right-handed squarks (see figure 2). The cross section is therefore dominated by electroweak production of JHEP01(2017)135  The present bounds from the LHC are already severely constraining the parameter space with a stau NLSP, pushing the mass of the lightest chargino χ ± 1 to be heavier than ∼ 900 GeV for M mess = 10 15 GeV. For M mess = 10 7 GeV both M 2 and µ are both forced to be light (see fig 2), which increases the electroweak cross section relative to M mess = 10 15 GeV. As a result, the stau NLSP for M mess = 10 7 GeV is already excluded by the present LHC data. The remainder of the stau NLSP parameter space for M mess = 10 15 GeV can be probed at the HL-LHC.

Stop/sbottom NLSP
Stop or sbottom can only be the NLSP if M mess = 10 15 GeV, in which case one expects the NLSP to be long-lived. Conversely to the stau NLSP case, direct production of NLSP pairs dominates the SUSY cross section almost everywhere in the parameter space. The present and projected constraints are shown in figure 18: the stop NLSP scenario is already very tightly constrained with the existing data, and will be fully probed at the HL-LHC.

Summary and conclusions
Based on the solution obtained in [12], we presented a full characterization of the phenomenology of general gauge mediation with stops below 4 TeV. The Higgs mass constraint is hereby crucial, and dramatically reduces the freedom inherent in the GGM parameter space. In particular: • The gluino is almost always decoupled at LHC, except for a small region of parameter space with µ < 0. The right-handed squarks are currently not accessible, but for high messenger scales they could be within the reach of the high luminosity LHC.
• The Higgsino mass is bounded from above. For µ > 0 the wino mass is also bounded from above while for µ < 0, the wino mass is bounded from below.
• The nature of the slepton NLSP is very tightly correlated with the nature of the lightest squarks, which determine the maximal colored production. While for a sneutrino NLSP the lightest squarks are left-handed and the colored production can be sizable, the lightest squarks are right-handed and strongly bounded from below if the NLSP is a stau. In the latter case only electroweak states dominate the SUSY cross section.
• Higher messenger scales enlarge the allowed parameter space. As a result, long NLSP life-times are favored in most models.
These features have a strong impact on the collider phenomenology: we find that both the stop/sbottom and stau NLSP scenarios are already very strongly constrained by JHEP01(2017)135 the existing data. For low messenger scale these NLSP types are already excluded, while for higher messenger scale both scenarios will be probed completely at HL-LHC. For the neutral NLSP's (wino, Higgsino, bino and sneutrino) there is a relatively small number of simplified topologies (figures 5, 9 and 16) which span nearly all of the accessible parameter space. We employed those handful of topologies to estimate the limits and asymptotic reach for each NLSP type in figures 4, 8, 13 and 15. Figure 19 gives an indication of the relative impact of the various datasets on the full GGM parameter space, where we marginalized over all NLSP types. In this figure the upper bound on the Higgsino mass is apparent by the presence of a meaningful limit from LEP, even for stop masses which are far outside of the reach of the experiment. For M mess = 10 7 GeV, the limits on the electroweakinos remain more important than the squark limits, even at the HL-LHC. The HL-LHC reach on the squark masses is therefore greater than one would expect from direct squark searches alone. This is no longer the case for high messenger scales (M mess = 10 15 GeV), where the direct squark bounds determine the reach of the HL-LHC.
Our analysis reveals a number of simplified topologies which, to the best of our knowledge, are currently not yet studied in detail. In particular, most attention so far has been devoted towards topologies with a bino NLSP, which are well-motivated benchmarks both because of their simplicity and because of their ubiquity in minimal gauge mediation. However now that the SUSY program at the LHC has matured, it may be worthwhile to investigate the topologies with Higgsino, wino or sneutrino NLSP somewhat more thoroughly. As we have shown in section 3, their phenomenology can be very complex and interesting. It would therefore be useful to analyze these topologies in sufficient detail to establish accurate limits and/or to identify potential blind spots or possible improvements to the existing analysis strategies. For example, for the Higgsino/wino NLSP, the stop/sbottom can each decay to both a t+MET and b+MET, which implies that stop/sbottom pair production can give rise to a tb+MET signature with a sizable branching ratio. Some of the latest CMS limits account now for this effect [120,121], and in particular for [121] the sensitivity depends rather strongly on the branching ratios. Other interesting examples are inclusive searches for disappearing tracks optimized for colored production, or more JHEP01(2017)135 generally simplified topologies with a wino and/or Higgsino and left-handed squarks (both 3th and lowest generations), and similar topologies with a sneutrino NLSP.

A NLSP life-time
A generic feature of gauge mediation scenarios is the presence of a light gravitino LSP whose couplings to the other MSSM particles are determined by the universal 2-body decay [122] Γ( X → XG) = The GGM phenomenology is then determined by the nature of the NLSP which decays to the gravitino LSP and its Standard Model partner. The NLSP life time in gauge mediation is then a function of the NLSP mass itself and the gravitino mass. The latter is directly related to the vacuum energy √ F 0 via super-Higgs mechanism m 3/2 = F 0 √ 3M P . Given the soft spectrum at the messenger scale we can estimate the SUSY-breaking scale F felt by the messengers by assuming a standard gauge mediation mechanism (with gaugino masses generated at 1-loop and scalar squared masses at 2-loops). Putting everything together we can write the decay length as where the factor k accounts for the fact that the vacuum energy can in general differ from the SUSY-breaking scale felt by the messenger sector (i.e F = kF 0 following the notation in [123]). In calculable models of SUSY-breaking typically k 1, which we assume throughout this paper. However one should keep in mind that there are no general results putting an upper bound on k.
Up to O(1) effects we can estimate the SUSY-breaking scale in terms of the gluino soft mass, which is typically the largest soft mass in our setup. We then get F ≈ 4π α 3 M mess mg and rewrite the decay length of the NLSP as  This formula shows that our benchmark of M mess = 10 15 GeV has long-lived NLSPs as a robust prediction independently on the details of the spectrum and on the particular UV completion. For M mess = 10 7 GeV we are instead in intermediate regime where O(1) effects become important and the NLSP decay length will generically depend on the details of the model. In order to simplify the phenomenological discussion, we assume the NLSPs to be always long-lived, also for M mess = 10 7 GeV.

B The effect of compressed spectra
We briefly discuss here how the LHC reach on GGM scenarios is affected by accounted for compressed spectra, where we consider the Higgsino NLSP as a example. The wino NLSP is has similar features, while for the bino NLSP the NLSP mass uncorrelated with the remainder of the spectrum. Our results are for the Higgsino NLSP shown in figure 20, which should be compared with figure 8 where compression was neglected. The lightest shading in figure 20, corresponding to the current bounds from the LHC, gets sensibly reduced with respect to the one in figure 8. The Higgsino NLSP can be compressed below the left-handed colored states, reducing the exclusion power of jet+MET, 2b + MET and 2t + MET searches. However it is important to keep in mind that, since we marginalize over the NLSP masses, the importance of the spectrum configurations which minimize the LHC bounds gets magnified in figure 20. The exclusions presented here therefore correspond to the exceptions, rather than to "typical" spectra.
Notice however that we still get a lower bound on the scale of left-handed colored states in figure 20. That is because the presence of multiple colored states close in mass JHEP01(2017)135 (the lightest stop, the lightest sbottom and the left-handed 1st and 2nd generation squarks) makes it difficult to squeeze the Higgsino NLSP below all of them making, such that some LHC searches remain effective. This is a common feature of SUSY scenario in which the production mechanism is dominated by left-handed states.
Comparing the darker shaded bands between figure 20 and figure 8, we see that compression effects are completely unimportant for HL-LHC and the asymptotic reach on GGM spectra at high-luminosity is essentially unchanged after accounting for compressed spectra. The reason is that the right-handed colored states become accessible and the upper bound on the mass of the Higgsino NLSP gives a minimum mass splitting between the right-handed states and the Higgsino itself.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.