Extended Gauge Mediation in the NMSSM with Displaced LHC Signals

We analyze models of extended Gauge Mediation in the context of the NMSSM, concentrating on supersymmetric spectra with light gluinos, low fine-tuning and decays of the lightest neutralino leading to displaced vertices. While the minimal scenario has rather heavy gluinos as a result of restrictions from the Higgs sector, we propose two new models in which the gluino can be as light as allowed by direct searches at the LHC, with a mass of about 1.7 TeV and 2.0 TeV, respectively. Both models have a tuning of a few permille, and lead to an interesting phenomenology due to a light singlet sector. A singlet state at around 98 GeV can account for the LEP excess, while the singlino has a mass of the order of 100 GeV and decays to b-jets and the gravitino, with decay lengths of a few cm.


Introduction
Gauge mediated supersymmetry breaking (GMSB) is one of the most elegant ways to explain the absence of sizable contributions to flavour violating processes in supersymmetric (SUSY) models. However, its simplest realizations are inconsistent with the measured Higgs mass of 125 GeV, unless one is willing to accept a very heavy (and thus unnatural) SUSY spectrum. In general there are two ways to make GMSB viable with sparticles light enough to be in the LHC reach. The first option to boost the Higgs mass is to maximize the loop corrections with sizable stop mixing, which requires direct Higgs-messenger couplings [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]. The second a mbadziak@fuw.edu.pl b nishita.desai@umontpellier.fr c cyril.hugonie@umontpellier.fr d robert.ziegler@cern.ch (less explored) option is to introduce additional treelevel contributions to the Higgs mass. Such corrections may originate from the mixing of the Higgs with a lighter singlet [22], and can be realized in a simple model proposed by Delgado, Giudice and Slavich (DGS) [23] that combines the NMSSM and GMSB with direct singletmessenger couplings. Indeed it has been shown in Ref. [24] that the most interesting region in parameter space features a light SUSY spectrum and a singlet around 90 ÷100 GeV, which gives a sizable tree-level contribution to the SM-like Higgs mass through mixing. Besides this "push-up" of the Higgs mass, the general framework of NMSSM+GMSB with a light singlet has other advantages. First, it provides a natural solution to the µ-B µ problem [23], as µ and B µ are generated dynamically through the vacuum expectation value (VEV) of the singlet. Second, the Next-to-Lightest SUSY particle (NLSP) is typically the singlino, whose decay into the gravitino is suppressed, leading to novel displaced signatures at colliders [25].
In the present paper, we focus on GMSB + NMSSM models with direct matter-messenger couplings, with a special emphasis on displaced signatures. This is motivated by recent LHC results for heavy Higgs searches in the τ τ channel, which has excluded the most interesting part of the parameter space of the DGS model. We investigate two simple extensions of the DGS model leading to a relatively light sparticle spectrum, which can be probed in the second run of the LHC, in particular associated with displaced decay signatures.
The rest of the article is organized as follows: In section 2 we analyse the general features of SUSY models with light singlets. In section 3 we review the DGS model and investigate the impact of recent LHC results on its parameter space. In sections 4 and 5 we introduce two new models with singlet-messenger couplings, and summarize and conclude in section 6. In an Appendix we provide complete expressions for soft terms in the respective models.

EWSB in the NMSSM with a Light Singlet
Let us start with a brief discussion of Electroweak Symmetry Breaking (EWSB) in the Z 3 invariant NMSSM, in which the µ-term in the MSSM superpotential is replaced by the following singlet couplings and the NMSSM specific soft terms are given by In order to generate a sufficiently large VEV for the singlet S ≡ s, the singlet soft mass m 2 S must be negative or at least much smaller than A 2 κ . This statement can be quantified in the limit s v, where v is the electroweak scale, in which one finds The approximate condition for proper EWSB in a global minimum of the potential reads In Minimal Gauge Mediation soft terms for the singlet sector do not arise at leading order which prevents successful EWSB. However, sufficiently large soft terms can be generated if there are direct couplings of the singlet to the messenger sector, as originally proposed in Ref. [26] and worked out in detail in Ref. [23]. It was further demonstrated in Ref. [24] that this model allows to realize the interesting "push-up" scenario where the SM-like Higgs mass gets a large positive contribution from mixing with a lighter singlet-like Higgs, thus allowing to lower the overall scale of the SUSY spectrum that drives the radiative corrections to the Higgs mass. In the following we revisit this model (DGS) in the light of updated experimental constraints, and compare it to the phenomenology of two similar models that feature additional singlet-messenger couplings. We are particularly interested in SUSY spectra with experimental signatures that may be probed at the LHC using displaced vertices, as studied in Ref. [25].
Before discussing these models, we note that in their most interesting regions of the parameter space, where the SUSY spectrum can be relatively light thanks to the push-up effect of Higgs-singlet mixing, correct EWSB generically implies large values of tan β. As a consequence, the most stringent constraints on these models often come from LHC searches for heavy MSSM-like Higgs bosons H/A decaying to τ τ , because the crosssection for H/A production grows with tan 2 β. Since this constraint turns out so important, it is instructive to discuss the origin of large tan β in the push-up region of NMSSM models with Gauge Mediation. Neglecting m 2 H d , which is usually a good approximation in the models under consideration, the expression for tan β is relatively simple: The soft parameters can in turn be related to the physical Higgs boson masses using the following approximate expressions for CP-odd Higgs masses: and the lightest (singlet-like) CP-even Higgs mass: We also note that for µ > 0, which we always assume in this analysis 1 , A κ must be negative to avoid tachyons in the singlet sector. In the limit w 1, corresponding to the case where EWSB is driven by a large negative m 2 S , one finds indeed large values of tan β Using this relationship between tan β and m a2 at face value with m h1 ≈ 100 GeV, the recent ATLAS constraint from Ref. [27], implies a limit tan β 25 or equivalently m a2 1.2 TeV, assuming no heavy Higgs decays to SUSY particles and neglecting threshold corrections to bottom quark Yukawa couplings. In typical GMSB models the value of m a2 is correlated to other sparticle masses, including squarks and gluino, so that a stronger bound on m a2 typically results in a stronger bound on coloured sparticles.
For smaller w the value of tan β slightly decreases. In the limit w ≈ 1/3 one finds Neglecting the terms proportional to A λ , tan β is smaller by a factor two with respect to the limit w 1, so 1 Also µ < 0 leads to viable spectra, but we will not consider these scenarios in the following, since the gluino is always quite heavy with a mass above 3.2 TeV.
that the LHC constraints on the a 2 are expected to be relaxed. Taking into account non-zero A λ , one can suppress (enhance) tan β when A λ is positive (negative). We also note that A λ may get positive contributions via RG running from a negative top trilinear term A t . Finally we provide an approximate expression for the Z-boson mass that is convenient to assess the finetuning. In the large tan β limit one has with z defined in Eq. (3). For the fine-tuning ∆ we use the Barbieri-Giudice measure [28] ∆ ≡ max where the maximum is taken over all UV parameters λ i . For a thorough discussion of the tuning measure in extended GMSB models see e.g. Ref. [10].

The DGS Model
The field content of the DGS model [23] consists of the NMSSM fields (the MSSM fields plus a gauge singlet S), in addition to two copies of messengers in 5 + 5 representations of SU (5). The superpotential is given by the NMSSM (see Appendix for our conventions), the spurion-messenger couplings of ordinary gauge mediation and new singlet-messenger couplings. Apart from the NMSSM part we have where Φ u + Φ d and Φ T + Φ T denote the doublet and triplet components in 5 + 5, respectively, and X denotes the SUSY breaking spurion that takes the VEV X = M + F θ 2 . This superpotential gives rise to soft SUSY breaking terms that can be found in the Appendix. These are determined by six parameters: the messenger scale M , the NMSSM couplings λ and κ, the DGS couplings ξ D and ξ T and the effective scale of soft SUSY breaking termsm ≡ 1/(16π 2 )F/M . One of these parameters (following DGS we choose κ) can be eliminated by requiring correct EWSB. We also impose a unification condition for ξ D and ξ T that allows to eliminate one additional parameter, In the following analysis we will always assume this relation, but we have checked that the more general case of independent ξ D and ξ T leads to similar phenomenology.
For heavy singlet-like scalars, all sparticles must be very heavy in order to satisfy the Higgs mass constraint in this model, and fall outside the discovery reach of the LHC [23]. It was found in Ref. [24] that a significantly lighter SUSY spectrum is possible in the presence of a light singlet that pushes up the Higgs mass via Higgs-singlet mixing. Having the light singlet mass fixed around 90÷100 GeV (where the constraints on the Higgs-singlet mixing from LEP [29] are weakest) implies ξ ∼ 10 −2 , while the Higgs-singlet mixing maximizing the push-up effect on the Higgs mass requires λ ∼ 10 −2 . For such small values of the couplings one has m 2 S ≈ −16m 2 ξ 2 T g 2 3 , so that the condition in Eq. (4) for correct EWSB is always fulfilled. As a result, the DGS model with a light singlet has w 1, so that tan β is approximately given by Eq. (8). An interesting prediction of the DGS model with a light singlet is a singlino NLSP with mass about mÑ 1 ≈ 100 GeV that decays to gravitino LSP and a 1 or h 1 , which typically decay further to bb. An estimate for the decay length up to O(1) factors is given by Therefore the singlino is long-lived and results in displaced vertices within the LHC detectors for sufficiently low values of the messenger scale M 10 7 GeV, corresponding to an NLSP decay length cτÑ 1 O(10) m. However, in order to see such a long-lived singlino at the LHC, it must have been produced from the decay of a heavy SUSY particle like gluino or squark, since the direct production of a singlino is negligible due its small couplings. Therefore displaced signatures require a sufficiently light colored spectrum, which might be in conflict with existing LHC searches.
In order to assess this issue, we have updated the analysis of the DGS model in Ref. [24] using the latest LHC constraints. For this analysis we have used the public code NMSSMTools [30,31], which computes the SUSY and Higgs spectrum and checks the latest LHC constraints in the Higgs sector (LHC constraints on SUSY sector will be dealt with later with CheckMATE2). We have performed large MCMC scans of the parameter space of the DGS model and sorted the results in a 2D histogram in the plane of the gluino mass vs. the decay length of the singlino NLSP. In each bin of this histogram, we have kept the point in parameter space minimising the fine-tuning as defined in Eq. (11). In Fig. 1 we present a map of fine-tuning in the plane of the gluino mass and the decay length of the singlino NLSP. We see that for NLSP decay lengths corresponding to a displaced vertex, the gluino mass is pushed far beyond 3 TeV. Thus, the DGS model does not predict any displaced signatures that could be observed at the LHC. We also note that even for larger NLSP decay length, for which the NLSP is stable from the LHC perspective, the gluino must be heavy and may be beyond the LHC discovery reach. Numerically, we find a lower bound on the gluino mass of 1.8 TeV, and a lower bound on the tuning of ∆ ≥ 600.
These negative conclusions can be traced back to recent heavy Higgs searches in the τ τ channel [27], which exclude previously viable points with light gluinos. This happens because successful EWSB in the DGS model with a light singlet requires rather large values of tan β 30, which pushes the heavy MSSM-like Higgs bosons to values above 1.3 TeV [27], and thus requires larger SUSY scales.
Therefore, in order to have a light sparticle spectrum that can be tested at the LHC, a model of NMSSM + GMSB is desirable that can provide sufficiently heavy MSSM-like Higgs bosons. In the DGS model the gluino mass is correlated with m a2 since the scale of both parameters is set bym. In order to avoid stringent lower bounds on the gluino mass, one should look for a model in which the correlation between gluino mass and m a2 is broken by new contributions to the UV soft masses and/or w ≈ 1/3 (for which tan β is generically smaller so the LHC constraints on m a2 are weaker). In the following sections we present two models that satisfy these requirements and therefore simultaneously allow both for a gluino much lighter than in the DGS model and sufficiently small singlino decay lengths to have displaced vertices at the LHC.

The U Model
In this new model, we take one copy of messengers and add all possible linear couplings of the messenger Φ u : This gives rise to soft SUSY breaking terms that can be found in the Appendix, and are now determined by six parameters:m, M, λ, κ, λ t and λ S d . Note that one might impose the condition λ S d y t = λ t λ that would result from explicit messenger-Higgs mixing [1]. However, we found that for small values of λ 10 −2 (as required to avoid experimental constraints on the Higgs-singlet mixing), λ S d would be too small to allow for m h1 to be in the preferred window between about 90 and 100  Fig. 2 The same as in Fig. 1 but for the U model and without applying constraints from LHC sparticle searches. Also shown are the benchmark points P1-P3 of Table 1.
GeV in order to to significantly enhance the Higgs mass. Thus, in the following we assume that λ S d and λ t are independent parameters, so that there is one additional parameter as compared to the DGS model.
In Fig. 2 we present a map of fine-tuning in the plane of the gluino mass and the decay length of the singlino NLSP, obtained analogously to Fig. 1. One can clearly distinguish two different regions here, in the upper part the "DGS-like" region, that is characterised by a large NLSP decay length cτÑ 1 100 cm, small fine-tuning ∆ 10 3 , a light singlet spectrum (m a1 ≈ 30 ÷ 40 GeV, m h1 ≈ 90 GeV, mS ≈ 100 GeV) and an essentially con-stant input value λ S d ∼ 10 −2 . Representative for the DGS-like region are the benchmark points P1 and P2 in Table 1. A closer look reveals that this region actually falls into two sub-regions, which can be separated by the value of the gluino mass as visible in Fig. 2. For gluino masses above 1.8 TeV, the value of |λ t | is essentially constant, ∼ 0.3, while it increases up to ∼ 0.7 towards lighter masses. Moreover, these parts are distinguished by the value of tanβ, which is about 10 in the left part and about 20 in the heavy gluino part. The latter is represented by P1 while P2 exemplifies the light gluino part of the DGS region.
The other region with smaller NLSP decay lengths cτÑ 1 100 cm is represented by benchmark P3 in Table 1. This region features much larger tuning ∆ 10 3 , essentially constant tan β ≈ 10, a heavier singlet spectrum (m a1 , mÑ We have therefore taken into account the LHC limits from direct SUSY searches with CheckMATE2 [32]. To briefly summarise the workflow, CheckMATE2 uses Pythia8 [33] to generate all accessible 2 → 2 processes followed by a detector simulation with Delphes3 [34]. Variables and cuts used in experimental analyses are then implemented as closely as possible to "recast" the analysis, and the expected number of signal events that pass the cuts are validated against published benchmarks and cut flows. This validated analysis can then be used to test New Physics models against published upper limits. In the presence of multiple signal regions that may potentially be sensitive to model predictions, CheckMATE2 selects only the most sensitive region with respect to the expected background. We summarize the results in Fig. 3 in the same parameter space as in Fig. 2, where red (green) points are excluded (allowed) by current LHC constraints. This plot shows that direct LHC searches still allow for a gluino as light as about 1.7 TeV, for essentially any NLSP decay length. The bound is significantly weaker than in typical simplified models presented by the experimental collaborations. The main reason for the relaxed constraints is that the wino is lighter than the lefthanded sleptons of the first two generations, and therefore dominantly decays to the lightest stau (which has some non-negligible left-handed component), resulting in τ 's instead of leptons in the final state. This feature strongly relaxes the mass limits both for direct produc-  Table 1.
tion of winos, as well as production of gluinos decaying predominantly to winos. We found that the most constraining searches for this model are the ATLAS searches with two same-sign leptons or tri-leptons [35], and jets and missing energy [36].
Let us now discuss how such a light gluino in the U model can be compatible with the constraints from the SM-like Higgs mass and experimental searches for heavy Higgses. First of all, the points with a light gluino feature large λ t , which implies large A t , so that the Higgs mass is enhanced not only by the push-up effect but also by the loop contribution from stop mixing. However, we recall that in the DGS model the lower bound on the SUSY scale does not arise from the Higgs mass constraint, but from direct searches for heavy MSSM-like Higgs bosons. In the U model instead, the pseudoscalar Higgs mass m a2 is enhanced by large λ t , as a result of the contribution to the soft Higgs mass parameter m 2 Hu ∼ −9y 2 t λ 2 tm 2 , cf. (A.14). Moreover, larger values of λ t are correlated with smaller tan β, as shown in Fig. 4. This is because larger λ t results in smaller w, which follows from and Eq. (3). For λ t 0.7, w approaches 1/3 for which tan β is twice as small as in the limit w 1, as ex-  plained 2 in Section 2. Therefore, the combined effect of increased m a2 and decreased tan β renders the LHC searches for heavy Higgs bosons with τ τ final states essentially insensitive to the U model with large λ t and light gluinos.
2 For large λt, tan β is suppressed even more due to the RG effect of At which makes A λ positive at the EW scale and cannot be neglected in Eq. (9). This effect is partially compensated by the fact that λt increases m a2 that enters in the numerator of the formula for tan β in Eq. (9).
We however notice that the points with light gluinos are more fine-tuned than DGS-like solutions with smaller λ t and gluinos above 1.8 TeV. This is because the tuning associated with λ t dominates the total tuning in almost all points of the U model. Using the approximate expression for M 2 Z in Eq. (10), one can easily show that the λ t contribution to the tuning ∆ λt is approximately given by where the first term is the contribution to m 2 Hu , and the second is the contribution to the effective µ-term.
While λ 2 /κ 2 ∼ 10 3 is essentially constant in both regions, the smallness of λ S d ∼ 10 −2 in the DGS-like region implies that the second term in Eq. (17) is small, and the tuning is essentially determined by λ 2 tm 2 , with a minimal value of ∆ ≈ 450. Instead the other region features large λ S d ∼ 0.2, which together with large λ t leads to a significant enhancement of the tuning with respect to the DGS-like region, by about a factor 100. This estimate explains the large increase in the tuning when crossing between the two different regions in Fig. 2. It also accounts for the approximately constant tuning over the two DGS-like sub-regions, where a decrease in the gluino mass (i.e.m) is compensated by an increase in λ t , such that the total tuning remains essentially constant, up to O(1) factors.
We conclude this section with a discussion of the three benchmark points P1-P3 in Table 1, which represent sample spectra of the U model with a light SUSY spectrum compatible with all LHC constraints, with a singlino decay length roughly below two meters. The point P1 features a tuning of ∆ ≈ 500, with a gluino mass around 2 TeV and NLSP singlino decay length about 2 m. This point represents the DGS-like region in the upper right corner with large tan β. The point P2 has larger λ t ≈ 0.5, which allows for a lighter gluino, although the tuning slightly increases as discussed above. The larger λ t implies smaller tan β ≈ 11, while the light singlet spectrum is quite similar to P1, with a NLSP singlino decay length that is still relatively large, cτÑ 1 ≈ 1.9 m. For both P1 and P2, the singlino NLSP mainly decays to the gravitino and the lightest pseudoscalar a 1 (with mass ∼ 30 GeV) which in turns decays promptly to bb. Characteristic for these benchmarks is a light Higgs state h 1 with mass around 90 GeV, which is mainly a SM singlet with a small doublet admixture corresponding to a mixing angle cos 2 θ ≈ 0.84. This state might account for the 2σ excess observed at LEP [29,37] in the bb channel, and could also explain the tentative hint for a light Higgs decaying into diphotons at CMS [38], although the signal strength of the light singlet state in the diphoton channel is about a factor of three smaller than the current central value of the CMS excess.
On the other hand, NLSP decay lengths below 1 m can only be obtained for the price of increased finetuning, in the region where λ t ≈ 0.7 is constant. This is exemplified by the point P3, which features cτÑ 1 ≈ 10 cm and a tuning of at least ∆ ≈ 10 4 . These points have a heavier Higgs spectrum, and in particular the lightest pseudoscalar is now heavier than 250 GeV, while the NLSP singlino is above 200 GeV. In contrast to P1 and P2, here the singlino mainly decays to gravitino and h 1 (with mass ∼ 110 GeV), since a 1 is too heavy.

The DGSU Model
Motivated by the previous analysis, we finally consider a model that combines the virtues of the DGS and the U model, thus allowing for light gluinos, small singlino decay lengths and low tuning. In this scenario, dubbed the "DGSU model", we take two messenger copies (i.e. N = 2) and introduce the following couplings: We impose the DGS unification condition at the GUT scale In contrast to the U model we can 3 now impose the Higgs-messenger mixing condition (for simplicity at the messenger scale) The superpotential couplings lead to additional contributions that can be found in the Appendix, and are given in terms of the six independent parametersm, M , λ, κ, ξ and λ t (the same number as in the U model).
In Fig. 5 we present a map of fine-tuning in the plane of the gluino mass and the decay length of the singlino NLSP, obtained analogously to Figs. 1 and 2. However, comparing to the exclusion limits from direct LHC searches shown in Fig. 6, we see that they are stronger than in the U model and the allowed points can have gluino masses only slightly below 2 TeV. The reason for this is that in the DGSU model there are two messengers (in contrast to one in the U model). Since the minimal gauge mediation contribution to gaugino masses is proportional to the number of messengers N while the corresponding contribution to sfermion masses scales as √ N , cf. (A.3), this makes sfermions lighter for a given gaugino masses. As a result, the squarks of the first two generations are comparable to the gluino mass and their production cross-section is non-negligible. For the same reason left-handed sleptons of the first two generations can now be lighter than wino, which results in more leptons in the final state (instead of τ 's as in the U model). The latter feature of the DGSU model also explains why the constraints for larger NLSP decay lengths are weaker, since the larger decay length is obtained for larger messenger scale, which in turn results in heavier sleptons (as compared to gauginos) due to their stronger renormalization by electroweak gauginos. Nevertheless, a gluino mass of 2 TeV is still viable for a decay length O(1) m, while decay lengths roughly between 20 and 100 cm implies a limit for the gluino mass of almost 2.5 TeV. We should also emphasize that this limit is particularly strong not only due to gluino/squark production, but also due to direct production of sleptons and winos, whose mass is correlated with the gluino mass. Therefore also direct searches for direct electroweak production set important constraints [39], besides analyses using jets and missing energy.
Note that for decay lengths of O(10) cm, there is a small strip of allowed points (represented by benchmark point P6), which however correspond to larger tuning, cf. Fig. 5. This is a consequence of large values of |λ t |, which also explains the relaxed LHC constraints. For large values of |λ t |, there is a large positive contribution to stop masses at the messenger scale, cf. (A.14), which leads to heavier stop masses at low scales, so that gluino decays to stops are kinematically forbidden.
We notice that the DGSU model shares many features with the DGS model, or rather the DGS-like region of the U model. In particular, all points (represented by benchmarks P4-P6 in Table 1) have a welldefined singlet sector with a singlet-like scalar h 1 around 90-100 GeV, a singlet-like pseudoscalar a 1 around 20-30 GeV and a singlino NLSP between 90-100 GeV. The singlino NLSP decay length can be rather short, of the order of O(10) cm, and it decays mainly to gravitino and a 1 , which decays promptly to bb. The input parameters are essentially constant throughout the whole region, with λ ∼ 10 −3 and ξ ∼ 10 −2 . As in the DGS-like region of the U model one can further distinguish two sub-regions, which are characterised by different values of |λ t | (that is essentially constant ∼ 0.4 for gluino masses above ∼ 2.2 TeV, and below this value starts increasing towards lighter gluino masses, up to ∼ 0.8) and tan β, which is correlated with λ t similar as in Figure 4. These two parameters essentially control the mass of the MSSM-like pseudoscalar a 2 and the tuning. The mass of the former grows for larger λ t (a feature inherited from the U model), and thus is no longer correlated with the gluino mass so that H/A → τ τ searches do not constrain this model. Moreover, large values of λ t lead to a large contribution to the Higgs mass from stop mixing, which implies that the overall SUSY scale can be lowered, and gluinos and stops can be quite light. The DGSU model is thus a perfect realization of the pNMSSM scenario analysed in Ref. [25], and motivates the combined searches using displaced and prompt signatures advocated in that article.
Finally, the DGSU model is also much less finetuned than the points with light gluinos in the U model. This can be again understood from the fact that λ t (which is sizable) dominates the total tuning, and contributes as in Eq. (17) with λ S d replaced by ξ. Similar to the DGS-like region of the U model, this contribution is small, so that the tuning is controlled by λ 2

Summary and Conclusions
In this article we have analyzed models of extended Gauge Mediation in the context of the NMSSM. The simplest scenario (DGS) requires rather heavy gluinos, which are constrained not by direct SUSY searches but rather a two-fold restriction coming from the Higgs sec-tor. On the one hand, the mass of the SM-like Higgs requires sizable loop corrections, on the other hand direct LHC Higgs searches in the τ τ channel set stringent lower limits on the mass of the heavy Higgses, which in turn is tied to the SUSY scale. We have proposed two new models with singlet/Higgs-messenger couplings (the U model and the DGSU model), which ease these constraints and thus allow for gluinos as light as allowed by direct searches. First, the presence of a direct coupling of stops to the messenger fields generate sizable A t at the messenger scale, which in turn allows for much lighter stops (and thus also other sparticles) consistent with the measured SM-like Higgs mass. Second, the same coupling also generates new contributions to the soft Higgs mass, which allows to both increase the MSSM-like Higgs doublet mass and reduce tan β, such that the LHC searches for heavy Higgs are satisfied without raising the overall SUSY scale. To identify the lower bound on the gluino mass in these models, we have recasted the existing LHC searches and found that gluinos can be as light as about 1.7 TeV (U model) and 2.0 TeV (DGSU model). The tuning in both models is rather low, and can be as small as 2 (U model) and 3 (DGSU model). The phenomenology is quite different in the two models: in the U model there are essentially two distinct regions, one with a DGS-like spectrum featuring a light singlet sector, large singlino decay lengths roughly above a meter (but small enough to see a displaced vertex in the LHC detectors) and low fine-tuning. The other region is characterized by a heavy singlet sector, less displaced singlino decays with cτÑ 1 ≈ 1 ÷ 10 cm and large tuning. The DGSU model combines the most interesting phenomenological features of both regions: light gluinos in the reach of LHC, a light singlet sector with a 98 GeV state (that might account for the LEP excess and improves the fit to the CMS data hinting at a new light state decaying to γγ), and displaced singlino decays into bb + MET with decay lengths as small as a few cm. In Table 1 we have collected benchmark points representative for the relevant parameter regions and the two models.
In this Appendix we provide further details about the models discussed in the main text, including the full superpotential, its motivation by symmetries, and the complete list of soft terms (obtained using the results of Ref. [10]).

0) NMSSM + Gauge Mediation
In order to fix notation, we define the NMSSM and minimal gauge mediation by the superpotential where X denotes the SUSY breaking spurion that takes the vev X = M +F θ 2 , N is the number of messengers in complete 5+5 representations of SU(5), and Φ u +Φ d and Φ T +Φ T denote the doublet and triplet components in 5 + 5, respectively. The spurion vev will induce soft masses and A-terms defined as Without introducing additional couplings, the A-terms vanish at the messenger scale and gaugino masses and sfermion masses are given by the usual minimal gauge mediation expressions:  particular the soft mass of the singlet vanishes at the messenger scale,m 2 S = 0. Together with the vanishing A-terms, this makes it difficult to trigger EWSB in minimal gauge mediation and motivates the introduction of additional interactions among messengers and NMSSM fields.

1) DGS Model
The following singlet-messenger couplings are added to the superpotential of Eq. (A.1): Notice that two copies of messengers are introduced (i.e. N = 2) in order to avoid that S has the same quantum numbers of X, which would lead to tadpoles terms that destabilize the hierarchy. The superpotential in Eq. (A.1) and Eq. (A.4) is the most general one allowed by a U (1) Z ×Z 3 symmetry with quantum numbers as in Table 3. Note that one can impose a unifi-  The new couplings leads to additional contributions to soft terms, on top of the soft terms from minimal GM in Eq. (A.3). Now A-terms for the singlet are generated at one loop: