Trilinear-Augmented Gaugino Mediation

We consider a gaugino-mediated supersymmetry breaking scenario where in addition to the gauginos the Higgs fields couple directly to the field that breaks supersymmetry. This yields non-vanishing trilinear scalar couplings in general, which can lead to large mixing in the stop sector providing a sufficiently large Higgs mass. Using the most recent release of FeynHiggs, we show the implications on the parameter space. Assuming a gravitino LSP, we find allowed points with a neutralino, sneutrino or stau NLSP. We test these points against the results of Run 1 of the LHC, considering in particular searches for heavy stable charged particles.


Introduction
Gaugino mediation [1,2] is a mechanism for mediating supersymmetry breaking in a setup with extra spacetime dimensions, which avoids flavor problems by suppressing the soft sfermion masses at a high-energy scale.The original version of the model also yields suppressed trilinear scalar couplings, which is unfortunate since the measured Higgs mass [3] then requires a unified gaugino mass of m 1/2 3 TeV and thus very heavy sparticles [4].
However, a simple extension of the scenario does allow for non-vanishing trilinears and thus a lighter sparticle spectrum [5].The couplings arise proportional to Yukawa couplings and thus do not lead to problematic flavor violation.We will investigate this possibility in detail in section 2, demonstrating explicitly how the trilinear couplings can be obtained.
In section 3, we study the parameter space of the extended setup.We show that the nonzero trilinears make it possible to reach the observed Higgs mass with sparticle masses that are accessible at the LHC.In gaugino mediation the gravitino can be the lightest supersymmetric particle (LSP) [6], making it a viable dark matter candidate [7]. 1 We assume this scenario, in which case the next-to-lightest sparticle (NLSP) can be a stau, a tau sneutrino or a neutralino [8].We determine the corresponding parts of the parameter space and constrain them by a careful analysis of LHC searches using data of the complete Run 1, in particular searches for long-lived heavy charged particles, extending the analysis in [5].

General setup
The present work considers one out of a class of higher-dimensional models.There are in general D spacetime dimensions, D − 4 of which are compact with volume V D−4 .This size determines the energy scale M c ≡ (1/V D−4 ) 1 D−4 needed to resolve the compact dimensions, referred to as the compactification scale.Fields can either live in the whole D-dimensional space referred to as the bulk or be localized on 3 + 1-dimensional branes that are located at different positions in the extra dimensions.The D-dimensional Lagrangian is [9] L D = L bulk Φ(x, y) + j δ (D−4) (y − y j ) L j Φ(x, y j ), φ j (x) , where j runs over the branes, x are coordinates on the branes, y are coordinates in the bulk, Φ is a bulk field2 and φ j is a field localized on the jth brane.Hats denote bulk fields with canonically normalized kinetic terms in D dimensions.We consider a model with two branes: the MSSM brane, where the visible matter fields are localized, and the hidden brane with a chiral superfield S, which is a singlet under the Standard Model (SM) gauge groups.Supersymmetry (SUSY) is broken by the vacuum expectation value (VEV) F S of the auxiliary field of S. The gauge and Higgs superfields propagate in the bulk.Therefore, they can couple directly to the SUSY-breaking field and obtain soft masses proportional to F S .In contrast, sfermion soft masses are strongly suppressed due to the separation between the MSSM and hidden brane, which avoids unacceptably large flavor-changing neutral currents [1,2].

Trilinear couplings
The supersymmetric part of the MSSM Lagrangian contains both bulk fields and fields constrained to the visible brane, where W is the visible-sector superpotential, Ŵ the field strength superfield and K the Kähler potential.Using the notation of equation ( 1), we have j = 1 for the visible brane and will accordingly use j = 2 for the hidden brane.On this brane, the gauge and Higgs superfields interact with the hidden-sector field S, where h, a, b u,d , c u,d and d are dimensionless couplings.The dots refer to terms containing only hidden-sector fields.Setting b u,d = 0 reduces the present case to the one considered in [2].
Setting also a = c u,d = d = 0, i.e., not placing the Higgs fields in the bulk, reduces our case to the one in [1].Note that the localizations of S and the sfermions forbid terms like S ū Ĥu Q and SQ † Q, which would directly yield trilinear couplings and sfermion soft masses.Interactions between the bulk fields and the hidden-sector field are non-renormalizable, so L D describes an effective theory valid up to some fundamental scale M .To obtain the 4-dimensional effective theory valid below the compactification scale, we integrate over the extra dimensions and keep only the zero modes of the bulk fields, which are constant in the extra dimensions.The integration yields a volume factor V D−4 in the kinetic terms of the bulk fields, so we define fields with canonical kinetic terms in 4D by Φ ≡ V D−4 Φ.Thus, the part of the effective 4D Lagrangian describing the interactions of S with the visible sector is The first term generates gaugino masses [1,2].We assume a unified gauge theory above the compactification scale, so that there is a unified gaugino mass m 1/2 .The remaining terms produce the Bµ-term, soft Higgs masses m 2 Hu and m 2 H d , and a contribution to the µ-term [2].The terms proportional to b u and b d , which were not included in the original versions of gaugino mediation [1,2], contribute to the soft Higgs masses and Bµ as well.Most importantly, however, they yield trilinear scalar couplings [5].This can be seen by absorbing them via the field redefinitions H u,d ≡ H u,d 1 + b u,d S M , from the general expressions for soft SUSY-breaking terms in the supergravity formalism, see e.g.[10,11], or by integrating out the Higgs auxiliary fields.We find it instructive to show the latter calculation for our particular case.
First, the part of the Lagrangian (2) that contains the Higgs supermultiplets' auxiliary fields where φ X denotes the scalar component of the superfield X. Adding the D-terms from equation (4) and employing the equations of motion ∂L/∂F † H u,d = 0 yields where we have omitted terms that do not contribute to SUSY-breaking trilinears. 3The solutions are thus omitting irrelevant higher-order terms in φ S and replacing the extra dimensions' volume by the compactification scale in the last step.Substituting F Hu and F H d into the Lagrangian (5) and replacing F S by its VEV finally gives rise to the desired trilinear terms, Consequently, we obtain trilinear scalar couplings proportional to the SUSY-breaking VEV and the Yukawa matrices, with Due to the proportionality of trilinear matrices and Yukawa matrices in the relations (9), these matrices are simultaneously diagonalized when changing to the super-CKM basis.Although the running to low energies leads to deviations from the exact proportionality, they are small enough to suppress flavor-changing neutral currents below the experimental upper limits.
Interestingly, the proportionality factors A u0 for the up-type squarks and A d0 for the downtype squarks and charged sleptons are different in general, in contrast to other simple setups for SUSY breaking like the Constrained MSSM or non-universal Higgs mass (NUHM) scenarios [12].In the following we will restrict ourselves to the simplest possibility A u0 = A d0 ≡ A 0 .

Constraints from naïve dimensional analysis
We will now estimate an upper limit on the trilinears, arguing that the couplings between the hidden-sector brane field S and the bulk fields can be constrained by naïve dimensional analysis (NDA) [9].This discussion generalizes results of [8], where the specific case of a 6-dimensional model was considered, to an arbitrary number of dimensions.
We write the Lagrangian (4) in terms of dimensionless fields Hu,d and S defined by where is the factor suppressing one-loop diagrams in D dimensions, and C is a group theory factor depending on the unified theory valid above M c .The volume factor V D−4 ensures canonical kinetic terms in 4D for the zero modes of the bulk fields.In this way, we obtain for the part of the Lagrangian coupling S to the Higgs fields According to NDA, the theory is weakly coupled below the cutoff scale M , if all couplings inside the curly brackets in equation ( 12) are smaller than one.This implies the constraints Combined with equation (10), they translate into the upper bound on the trilinears.For comparison, the NDA constraint on the gaugino mass is [6] Consequently, the ratio of the upper bounds is simply If the limit on m 1/2 is saturated, it is thus possible for the trilinear couplings to be somewhat larger than the gaugino mass, but not by orders of magnitude.

Phenomenology of the model
Let us now explore the parameter space of gaugino mediation extended by trilinear couplings.As explained in section 2, the model contains the five free parameters m 1/2 , m 2 Hu , m 2 H d , A 0 , and Bµ.The soft squark and slepton masses are negligibly small.This is a realization of the NUHM2 scenario [13] with the restriction m 0 = 0.These input parameters are boundary conditions at the compactification scale, which we identify with the scale of gauge coupling unification, M c 10 16 GeV.As usual, we trade Bµ for tan β and use the measured Z mass to determine the absolute value of µ.We choose µ to be positive and restrict ourselves to negative values for A 0 ; changing the sign of both parameters would lead to a similar phenomenology.
One of the most important model restrictions is the Higgs mass required to match the value measured at the LHC, see section 3.1 for details.The allowed parameter space accommodates various choices of the lightest sparticle of the MSSM, discussed in section 3.2.It comprises the lightest neutralino, the tau sneutrino and the lighter stau.As the latter two are not phenomenologically viable dark matter candidates we assume here that the LSP is a non-MSSM sparticle with very weak interactions. 4In the framework of supergravity, this could be the gravitino.In this case the lightest sparticle of the MSSM is the NLSP.Gaugino mediation allows for gravitino masses m 3/2 10 GeV [6], in which case the NLSP becomes stable on collider time-scales and the collider signature of the considered model vitally depends on the choice of NLSP.While a neutralino or sneutrino NLSP provides a signature containing missing transverse momentum, detector-stable staus provide a distinct signature of heavy stable charged particles (HSCPs), for which the LHC sensitivity is very high.LHC constraints for the respective signatures are discussed in section 3.3.Bounds from color or charge breaking minima of the scalar potential are briefly discussed in section 3.4.In section 3.5 we comment on the cosmological constraints on the model.

Higgs mass
One of the most important constraints on the parameter space is the experimentally observed Higgs mass of 125.09 ± 0.24 GeV [3].The theoretical uncertainty of the Higgs mass prediction in the MSSM is on the order of ∼ 2 GeV [15,16].As the theoretical error is large compared to the experimental one, we do not consider the latter.Furthermore, we assume that the lightest CPeven Higgs of the MSSM plays the role of the observed Higgs.Hence, we consider points with a theoretically predicted mass of the lightest CP-even Higgs in the rage 123 GeV m h 127 GeV to be consistent with observations.
In order to compute the Higgs mass we proceed as follows.First we use SPheno 3.3.8[17,18] for the calculation of the sparticle masses and low-energy Lagrangian parameters.The output from SPheno is then used as input to FeynHiggs 2.12.2 [15,16,[19][20][21][22][23], which we use to more accurately calculate the lightest Higgs pole mass.Both programs incorporate two-loop diagrams in the calculation of m h .However, FeynHiggs 2.12.2 includes a more complete treatment of the calculation, including momentum dependent two-loop QCD contributions [16], leading threeloop contributions [15] and additionally, by combining an effective field theory approach with the fixed-order calculation, it incorporates up to NNLL contributions resummed to all orders [19].This treatment can significantly reduce the theoretical uncertainties with respect to the pure fixed-order calculation, in particular for large M susy ≡ √ m t1 m t2 [19,24].
The result for the Higgs mass 5 is shown in figure 1, where the left panel shows the contour for which m h = 125.09GeV in the A 0 -m 1/2 plane.The darker and lighter shaded regions around it denote the ±1 and ±2 GeV bands respectively.As mentioned above, we use the Higgs mass as computed by FeynHiggs, represented by the blue curve and bands on the plot.The right panel shows the Higgs mass dependence on tan β, m 2  Hu and m 2 H d for a fixed value of m 1/2 and three choices of A 0 .
For tan β = 10 and vanishing A 0 , very large values of m 1/2 on the order of 6 TeV are needed to achieve a suitable Higgs mass of 125 GeV.With growing negative A 0 , the required m 1/2 drops to a minimum around m 1/2 2 TeV, beyond which the Higgs mass rises again.This minimum corresponds to the maximal mixing scenario, where [26] for a detailed discussion.This result shows that only with a non-zero trilinear coupling A 0 , a Higgs mass of around 125 GeV can be obtained with m 1/2 such as to obtain a sufficiently light spectrum to be observable in upcoming collider experiments.See further discussion in section 3.3.
The ±1 and ±2 GeV bands span a large range, reflecting the relatively large uncertainty in the 5 We used the most recent results available in [25] FeynHiggs, m 1/2 = 3 TeV required value of m 1/2 between 3 and 8 TeV.However, this uncertainty band shrinks significantly for large negative A 0 .The dependence on tan β is shown in the right panel of figure 1.Both very small and very large values of tan β cause the Higgs mass to drop drastically, making it hard to achieve the correct Higgs mass even for very large m 1/2 .Note that for large tan β and large negative A 0 , the spectrum acquires tachyonic states.Therefore, not all curves extend to tan β = 50.
The influence of the Higgs soft masses m 2 Hu and m 2 H d on the Higgs mass is small throughout the explored parameter space.The most significant effect arises for large tan β, cf. the solid and dashed curves in the right panel of figure 1.
The Higgs mass contour as computed by SPheno, presented by the red curve and shaded bands in the left panel of figure 1, is included for comparison. 6The required Higgs mass is reached with considerably smaller m 1/2 for a given A 0 , as the SPheno result for m h is typically around 3 GeV larger than the one from FeynHiggs.In particular for large M susy , NNLL resummation can yield important corrections that significantly contribute to the difference between the results obtained by the two codes, see e.g.[16,19,24] for details.

Particle spectrum
The phenomenology of the model regarding collider searches, astrophysics and cosmology strongly depends on the nature of the NLSP.As mentioned above, we compute the sparticle spectrum with SPheno.In the considered parameter space, we encounter three possible candidates for the NLSP: the neutralino, the sneutrino, or the lighter stau, which can be predominantly left-or right handed.Figure 2  The results summarized in figure 2 demonstrate the relationship between the Higgs soft masses and the NLSP.As the ratio r ≡ (m 2 Hu − m 2 H d )/m 2 1/2 becomes more negative, the NLSP can shift from the stau, to the neutralino and finally to the sneutrino, depending on the value of tan β and A 0 .If tan β is relatively large and A 0 is large and negative, only a stau NLSP is possible.Interestingly, the stau NLSP is also observed to shift through regions of right-chirality, large mixing and left-chirality with decreasing r (cf. the gray solid curve in the plots of figure 3, showing the stau mixing angle).In addition, figure 2 depicts the NLSP sensitivity to the value of tan β, showing that the stau NLSP region grows with tan β.In fact, for tan β 30, the entire region contains only a stau NLSP.We also find that some of the regions of interest contain unphysical tachyonic spectra, meaning negative soft-masses squared.This occurs when A 0 has a large negative value compared to m 1/2 , and becomes more frequent with increasing tan β.
We would like to explain some of this behavior in a rough analytical manner, beginning with the chirality switch of the stau.This can be understood from analyzing the one-loop RGE's for the third generation leptonic soft masses [11] where For m 2 H d m 1/2 , we can neglect the gaugino masses in the above formula, and the running will depend mostly on the Σ parameter.From equation (18), one sees that for very large m 2 H d , this value is negative, and will therefore lower the value of the left-chiral soft mass term but increase the size of the right-chiral term.Therefore, the NSLP will become more left-chiral with increasing m 2 H d .For larger values of m 2 Hu , the absolute value of the Σ term is smaller, and the progression from right-to left-chirality happens at larger values of m 2 H d .Regions where the sneutrino becomes the LSP are also determined by equations (17).Again, these regions occur in the limit m 2 H d m 1/2 , so we can make the same approximation and assume that the stau is mostly left-chiral.When the stau is mostly left-chiral, it is a delicate matter which of the two particles becomes the NLSP.The sneutrino mass is completely determined by equation (17a), as there are no right-chiral neutrinos in the MSSM, whereas there is mixing in the stau sector.The off-diagonal elements in the stau mixing matrix, which are A 0 and tan β dependent, push the eigenvalue down.However, the diagonal elements, which are predominantly dependent on the soft masses m 2 L 3 and m 2 ē3 , but also depend on the "hyperfine splitting" arising from EWSB, increase the eigenvalues.In figure 3 we show the masses of the staus, the tau sneutrino and the neutralino for the two slices denoted by the black dotted lines in the lower panels of figure 2. It reveals the small mass difference between τ 1 and ν τ for large m 2 H d /m 2 1/2 .The tan β and A 0 dependence can be understood by first noting that the neutralino mass is pushed up with tan β, and larger values of A 0 push the third generation leptonic soft masses down by increasing χ τ .This explains the shrinking neutralino region seen in the lower panels of figure 2. Large values of A 0 also increase mixing in the stau sector, pushing down the smallest eigenvalue of the stau mass matrix, implying the sneutrino LSP region should also shrink with larger tan β.Hu .All panels have m 1/2 = 2 TeV.In the white region below, we run into a tachyonic region, i.e., negative soft masses squared.The red dashed curve indicates the transition from a predominantly right-to left-handed stau NLSP, i.e., the contour sin 2 θτ = 1/2.The black dotted lines in the lower plots denote the slices in parameter space that are considered in figure 3.

Heavy stable charged particles
The lighter stau is the NLSP for a large part of the considered parameter space in our model.In order to determine the 95% CL exclusion limits from collider searches for HSCPs, we first compute the total cross section for the production of sparticles with Pythia 6 [27].For points with σ tot 8 TEV > 1/L int 8 TEV , i.e. for an expected total signal of more than one event we perform a Monte Carlo simulation of the signal at the 8 TeV LHC with the MadGraph5_aMC@NLO event generator [28].We generate 10 k events for each point in the model parameter space, taking into account all possible sparticle production channels.The decay, showering and hadronization is performed with Pythia 6 [27].We do not perform a detector simulation.Instead we determine the signal efficiencies with the method introduced in Ref. [29], which allows for the direct analysis of the hadron-level events on the basis of the kinematic properties of isolated HSCP candidates.In order to identify isolated HSCP candidates we first impose the isolation criteria and where the sums include all charged and visible particles, respectively, in a cone of ∆R = ∆η 2 + ∆φ 2 < 0.3 around the direction of the HSCP candidate, p T i denotes their transverse momenta and E i their energy.Muons are not considered as visible particles as their energy deposition in the calorimeter is small.|p| is the magnitude of the three-momentum of the HSCP candidate.The HSCP candidate itself is not included in either sum.
We compute the signal efficiency by averaging the probabilities for events to pass the on-and off-line selection criteria [29], where the sum runs over all N generated events i.For events containing one or two HSCP candidates the probabilities are given by on/off, i = P on/off (k 1 i ) or P (2) respectively, where k 1,2 i are the kinematical vectors of the HSCP candidates in the ith event.k = (η, p T , β) contains the candidate's pseudo-rapidity, η, transverse momentum, p T , and velocity, β.
The CMS analysis [29] requires a minimum reconstructed mass, m rec , for the candidate.The probabilities P on/off (k) are provided for four distinct mass cuts m rec > 0, 100, 200, 300 GeV , which we here consider to be four different signal regions.Due to detector resolution effects, the reconstructed mass is typically m rec 0.6 m HSCP [29].Hence, we set the efficiencies to zero if 0.6 m HSCP is below the respective mass cut of the signal region.This prescription is also used in Ref. [30], where it is validated by reproducing the efficiencies and cross section upper limits for the gauge mediated supersymmetry breaking model from the full CMS detector simulation [29] with a relative error below 5%.25)-( 26).
The grey dotted curves show the contours of the lighter stau mass m τ 1 .For tan β = 50 and −A0 2.3 TeV the HSCP limit (dot-dashed curve) extents into the region of a tachyonic spectrum, in this region this limit is only an extrapolation.
The resulting limits are shown in figure 4, projected onto the A 0 -m 1/2 plane for two slices in parameter space, where m 2 Hu = m 2 H d = 0, and tan β = 10 (left panel) and tan β = 50 (right panel).Both choices are characterized by a stau NLSP in the entire considered parameter plane.The considered CMS search for HSCPs at the 8 TeV LHC excludes the region below the red dotdashed line (red shaded region) at 95% CL.The exclusion reach depends strongly on the overall sparticle mass spectrum, which is indicated by drawing several contours for the mass of the stau NLSP.The exclusion limits turn out to cut at around m τ 1 400 GeV with a mild dependence on the other parameters.This translates into a limit on m 1/2 between 1 and 2 TeV for tan β = 10 in the considered region of A 0 , but can be much larger for large tan β, as shown in the right panel.The existing limit only touches the −2 GeV band regarding the Higgs mass, and leaves most of the parameter space that provides a Higgs mass in the range 123 GeV m h 127 GeV unchallenged.
The 13 TeV LHC runs have pursued searches for heavy stable charged particles, and (preliminary) results from an integrated luminosity of 2.5 fb −1 [31] (12.9 fb −1 [32]) have been released.For the 13 TeV searches, no on-/off-line probabilities (as in Ref. [29]), have been provided, such that these searches cannot be easily reinterpreted.We do, however, expect to obtain a meaningful estimate of the 13 TeV sensitivity as described in the following.The signal efficiencies for the 8 TeV LHC increase with increasing m τ 1 for the tested points, and are only mildly dependent on the other parameters within the considered model.In particular, we found that the efficiency is always above 0.5 for m τ 1 > 350 GeV, and above 0.6 for m τ 1 > 500 GeV in our scan.Assuming a similar detector performance, the efficiencies at the 13 TeV LHC for a certain stau mass will to first approximation be the same as for the 8 TeV efficiency, for a mass that is smaller by a factor of 8/13.Hence, for the 13 TeV LHC we assume an efficiency of 0.5, which is expected to provide a mostly conservative estimate for stau masses above 600 GeV.Furthermore, as for m rec > 200 GeV the signal region is typically background-free [29,31] we require 3 signal events in the signal region supporting a 95% CL exclusion limit.In this way we estimated the projected sensitivity for 300 fb −1 at 13 TeV, for which we computed the production cross sections with Pythia 6 [27], see the red dot-dot-dashed curves in figure 4. The projected exclusion reach cuts into a larger portion of the parameter space providing the correct Higgs mass.In particular, the maximal mixing scenario for moderate values for tan β can be tested.With 300 fb −1 , stau masses up to around 1 TeV could be tested.
Note that performing the same estimate for the analysis at 2.5 fb −1 (12.9 fb −1 ) provides an estimated limit very close to (slightly above) the 8 TeV limit, which we do not show in figure 4 for the sake of better readability.

Missing energy signatures
As discussed in section 3.2, a high enough m 2 H d relative to m 2 Hu and m 1/2 results in a neutralino or even sneutrino NLSP.If present in collision events, neutral NLSPs lead to a missing transverse energy (MET) signature at the LHC.
In order to test the compatibility with current LHC results, we perform a Monte Carlo simulation with the MadGraph5_aMC@NLO event generator [28] for the 8 TeV LHC.We generate 20k events.The decay, showering and hadronization is performed by Pythia 6 [27].The results are used as input to CheckMate 1 [33], 7 allowing us to simultaneously test the signal against various LHC searches for missing transverse energy.
We test our model against all ATLAS analyses implemented in CheckMate 1 .These analyses search for final states containing a significant amount of missing transverse energy, in addition to jets or leptons.The signal is compared to experimental limits in the respective signal regions of the analysis at 95% CL.The most sensitive region from all the analyses is used to conclude whether the model can be excluded or not.Among the points that provide a Higgs mass m h > 123 GeV, we tested the lighter part of the spectrum, i.e., m 1/2 ≤ 3 TeV for various slices in parameter space regarding tan β, A 0 and Higgs soft masses.We found that even for the lightest spectra the signal falls below the exclusion limits by at least an order of magnitude.Since the spectrum becomes heavier for larger values of m 1/2 , we expect no sensitivity of searches for MET in the region m h > 123 GeV.The analysis which most frequently has the largest signal region is the search for direct stop pair production in final states with two leptons [38].

Charge and color breaking
In addition to the collider constraints, we investigate whether and in which regions of parameter space the current model is limited by charge-and color-breaking minima of the scalar potential.The MSSM contains 26 scalars, most of which carry electric or color charge.Hence, there is a danger of introducing charge and color breaking (CCB), depending on the VEVs where the scalar potential has its minimum.Due to the large number of scalars in the theory, the scalar potential is very complex, limiting an analytical approach to only considering certain rays in field space.It is common to investigate directions in field space where the VEVs of the Higgses and τL/R or tL/R have the same value, and to neglect the D-term of the potential, which is a gauge interaction and positive for non-zero values of the scalar fields, as well as loop corrections.Based on criteria for CCB as found in [59][60][61], we use the same condition as [5] for the stop trilinear coupling, namely By analogy, we take the bound on the stau trilinear to be For large tan β, one can derive an upper bound on the product µ tan β requiring the standard electroweak vacuum to be stable or metastable with a lifetime larger than the age of the universe [62][63][64][65].We use [65], where tan and These bounds are superimposed in figure 4, where we show the most constraining bound from equations ( 25)- (26).For tan β = 10, the region below the purple dashed line violates equation (25), while for tan β = 50 it violates equation (26).For large negative A 0 , the CCB bound cuts into the part of the parameter space that provides the correct Higgs mass.Note that we impose these bounds as a first estimate, indicating the region where CCB constraints might exclude points in the parameter space.It has been shown [66,67] that these bounds are useful, but not entirely reliable in determining vacuum stability when more sophisticated analyses are performed.We leave a detailed numerical analysis of the vacuum stability utilizing Vevacious [68] for future work.

Cosmological constraints
Scenarios with long-lived NLSPs are subject to constraints from big bang nucleosynthesis (BBN) because the presence and late decays of the NLSPs can change the primordial abundances of light elements [69][70][71].In our case, the NLSP decays comparatively early on BBN timescales due to the relatively heavy sparticle spectrum.For example, for gaugino mediation with a stau NLSP a lower bound of m τ 400 GeV was found in [72], which roughly coincides with the lower limit from HSCP searches.Therefore we do not perform a detailed analysis here.
Another constraint we did not include is the non-thermal production of gravitino dark matter by NLSP decays, which may not exceed the observed dark matter density.This is interesting from a theoretical point of view because it leads to an upper bound on the sparticle masses but less relevant for phenomenology, since the constraint becomes relevant only for very large values of m 1/2 [72], which are far beyond the reach of the LHC.

Conclusions
We have considered phenomenological constraints on the gaugino mediation model of supersymmetry breaking.First, we verified that the model allows for soft trilinear scalar interaction terms.These terms were originally assumed to vanish in gaugino mediation and play a crucial role in achieving a Higgs mass in agreement with the observed value of 125 GeV.The trilinear matrices are proportional to the Yukawa coupling matrices, thus avoiding flavor problems.The proportionality factor can be different for up-and down-type sfermions.
Second, we explored the phenomenological consequences of non-vanishing trilinears.The first constraint we discussed is the experimentally observed Higgs mass, calculating the low-energy parameters and the sparticle spectrum with SPheno and the Higgs mass with FeynHiggs.We determined the parameter space regions where the Higgs mass lies within the LHC limits.Large negative trilinears are required to obtain an acceptable Higgs mass if the SUSY scale is to be kept near the reach of the LHC.We also observe that FeynHiggs 2.12.2 -incorporating important NNLL contributions -predict a Higgs mass around 3 GeV lower compared to the SPheno calculation in the parameter regions considered.
We also considered the phenomenological implications of the non-universal soft Higgs masses.We found that these parameters mainly affect which sparticle becomes the NLSP (we assume a gravitino LSP and that the lightest MSSM sparticle is the NLSP).Values of the ratio r ≡ (m 2 Hu − m 2 H d )/m 2 1/2 near zero correspond to a stau NLSP.As r is pushed to larger negative values, the NLSP can become the neutralino and eventually the tau sneutrino.This behavior also depends on A 0 and tan β.For sufficiently large |A 0 | and tan β, the composition of the stau NLSP changes from mainly τ R to mainly τ L as r becomes large and negative, passing through regions with large mixing.
Proceeding to investigate the LHC sensitivity of the scenario, we found that for a neutral NLSP, the viable part of parameter space is not challenged by missing energy searches.However, for a stau NLSP, the corresponding searches for heavy stable charged particles become sensitive and cut into the region where 123 GeV m h 127 GeV.The projection for an integrated luminosity of 300 fb −1 reaches a large portion of this part of parameter space, especially in the maximal-mixing scenario.
Finally, we indicate in which regions of parameter space the model might be limited by charge-and color-breaking minima of the scalar potential by using (semi-)analytic estimates for the CCB conditions.It turns out that only a small part of the allowed Higgs mass region is in conflict with these CCB bounds.

F
e y n H i g g s S P h e n o m h = (125.09,±1, ±2) GeV tan β = 10 , m 2 Hu = m 2 H d

Figure 2 :
Figure 2: Regions characterized by a stau (red), neutralino (blue) and sneutrino (green) NLSP in the m 2 H d /m 2 1/2 -A0/m 1/2 plane for four choices of tan β and m 2Hu .All panels have m 1/2 = 2 TeV.In the white region below, we run into a tachyonic region, i.e., negative soft masses squared.The red dashed curve indicates the transition from a predominantly right-to left-handed stau NLSP, i.e., the contour sin 2 θτ = 1/2.The black dotted lines in the lower plots denote the slices in parameter space that are considered in figure3.

Figure 3 : 1 (
Figure 3: Sparticle masses m τ 1 (red long-dashed curve), m τ 2 (orange short-dashed curve), m ντ (green dot-dashed curve) and m χ 0 1 (blue dotted curve) as a function of m 2 H d /m 2 1/2 for two choices of tan β and m 2 Hu .The stau mixing angle is indicated through the gray solid line showing sin 2 θτ , labelled on the right axis.

Figure 4 :
Figure 4: Contours of m h = 125.09GeV computed by FeynHiggs (blue solid curve) in the A0-m 1/2 plane, as well as constraints from searches for heavy stable charged particles (HSCP) at the 8 TeV LHC (red shaded region below the red dot-dashed curve).Projections for the 13 TeV LHC at 300 fb −1 are indicated by the red dot-dotdashed curve.The purple dashed line represents the strongest of the CCB constraints from equations (25)-(26).The grey dotted curves show the contours of the lighter stau mass m τ 1 .For tan β = 50 and −A0 2.3 TeV the HSCP limit (dot-dashed curve) extents into the region of a tachyonic spectrum, in this region this limit is only an extrapolation.
shows several projections of the parameter space in the plane m 2 H d /m 2 1/2 -A 0 /m 1/2 .We have rescaled m H d and A 0 by m 1/2 as the nature of the NLSP is almost independent of the overall mass scale that is governed mostly by m 1/2 .In other words, for fixed ratios A 0 /m 1/2 , m 2 H d /m 2 1/2 and m 2 Hu /m 2 1/2 , the sparticle spectrum is mainly shifted with m 1/2 and the shown projections remain approximately unchanged.