En-gauging Naturalness

The discovery of a 125.5 GeV Higgs with standard model-like couplings and naturalness considerations motivate gauge extensions of the MSSM. We analyse two variants of such an extension and carry out a phenomenological study of regions of the parameter space satisfying current direct and indirect constraints, employing state-of-the art two-loop RGE evolution and GMSB boundary conditions. We find that due to the appearance of non-decoupled D-terms it is possible to obtain a 125.5 GeV Higgs with stops below 2 TeV, while the uncolored sparticles could still lie within reach of the LHC. We compare the contributions of the stop sector and the non-decoupled D-terms to the Higgs mass, and study their effect on the Higgs couplings. We further investigate the nature of the next-to lightest supersymmetric particle, in light of the GMSB motivated searches currently being pursued by ATLAS and CMS.


Introduction
The recent discovery of a ∼ 125.5 GeV particle consistent with the properties of the standard model Higgs boson [1,2] and no direct evidence of supersymmetry (SUSY) in the current LHC data are pushing traditional setups of gauge-mediated supersymmetry breaking (GMSB) for the minimal supersymmetric standard model (MSSM) into fine-tuned territory. In the MSSM, to obtain a lightest CP-even scalar of the observed mass requires either heavy stops, thereby introducing a naturalness or fine-tuning problem, or substantial leftright stop mixing, which, being strongly dependent on the trilinear soft term A t , is heavily influenced by the mechanism of supersymmetry-breaking mediation that is invoked. In particular, in GMSB trilinear terms such as A t are vanishing at the supersymmetry breaking scale M , and a large A t can only be generated via renormalisation group evolution. This requires the scale M to be very high, which is also detrimental to the naturalness of the theory. Moreover, in minimal SUSY models, choosing heavy stops results in the sparticle spectrum becoming heavier and beyond the reach of the LHC, and consequently phenomenologically less interesting. The heart of the problem appears to be that in the MSSM, the tree-level Higgs mass is simply too small, the upper bound being the mass of the Z boson. If however we drop minimality from our criteria, natural scenarios of supersymmetry breaking with discovery potential still persist. These usually involve lifting the Higgs potential at tree-level through non-decoupled F or D type terms [3][4][5], the later of which arise when the MSSM is extended by additional gauge groups. The less-studied gauge extensions involving non-decoupled D-terms could enhance the tree-level Higgs mass, resulting in detectable deviations in the Higgs couplings, and further induce suppressions to the scalar masses compared to minimal GMSB. Such models may therefore have a direct impact on phenomenology and, unsurprisingly, have of late found increasing interest [6][7][8][9][10][11][12][13][14][15][16].
In this paper we study a quiver model or gauge extension of the MSSM within the framework of GMSB, to determine whether it is indeed possible to reproduce the observed Higgs boson mass while keeping the stop masses below roughly 2 TeV. In particular we build a tailor-made spectrum generator for our model using the publicly available tool SARAH [17][18][19]. This allows us to perform the renormalization group evolution at two loops and analyze several aspects of the model's phenomenology. Although we do not carry out a thorough "naturalness" or fine-tuning study, it is at least clear that qualitively, having stops lighter than benchmark minimal GMSB certainly improves the relative naturalness of the model. We therefore study the resulting spectra for the model, consistent with experimental results, in particular the Higgs sector. We further investigate possible signatures of this model at the LHC, taking into account the latest results of GMSB motivated searches.
It is useful to summarise the current status of the literature that explores two site or minimal quiver models. Initially a simple two site deconstruction, similar to our model MI (described later in section 3.1), with gauge mediated boundary conditions was proposed in [7], in particular in which both quiver gauge groups G A and G B are U (1) × SU (2) × SU (3) and in which both site A and B gauge couplings unify separately. As in such models sfermion soft masses are typically smaller than gaugino soft masses it unwittingly allows for the foundations of a natural spectrum. However the resulting tree-level Higgs mass is described by the MSSM, such that the observed Higgs mass [1,2] must therefore be generated by large A t or heavy stops. Later in [10] the non decoupled D-terms of [3] in the Higgs sector were included by hand, although they did not contribute significantly and the resulting Higgs mass of their benchmarks, of around 116 GeV, are also now firmly excluded. In [8,9] a model similar to our model MII 1 , was sketched in which it is argued that one could obtain a) a linking field vev v < 10 TeV b) unification of each gauge site separately, as well as c) lighter 3rd generation scalars relative to 1st and 2nd, due to location of generations in different sites. Most of the studies so far mentioned were carried out at tree-level. The formulas for the soft masses in the case of gauge mediation were also given in [12,23,24] where some interest had developed in finding and making precise models where sfermions were lighter than gauginos at the messenger scale. In light of the Higgs discovery [1,2] and naturalness, the split families models of [8,9] have re-emerged [13,14,25].
In [14] a study using the MSSM spectrum generator, SOFTSUSY [26] combined with some private codes found that the model based on [3, 8-10, 12, 13, 23-25] could not simultaneously achieve unification and obtain the necessary enhancement of the Higgs mass. If the necessary enhancement was obtained then the gauge groups hit a Landau pole much before the GUT scale. Some models involving three sites were proposed to alleviate the problem. The study of [14] misses some leading order effects to some one loop RGEs and as such a more comprehensive study and implementation of these models is necessary. Furthermore, simpler (more minimal) variants of this class of model may fair better. The core issue is that the SU (3) A × SU (3) B matter component of the quiver introduces additional flavours to each side respectively, charged also under the electroweak groups, and this additional matter content can cause the Landau pole problem, if one wishes for g A > g B to increase the effect of the Higgs enhancement. Given the importance of the Higgs enhancement, perhaps it is worth removing the copy of SU (3) (and possibly complete unification) in favour of naturalness. Such a setup should anyway test what is the maximum allowed enhancement to the Higgs from the electroweak quiver without worrying about unification or Landau poles.
There are clearly a number of important and unresolved issues which our paper addresses. The most detailed study so far, in [14], although useful does not include some leading one-loop contributions to the RGEs (see appendix A and [27]), it is essentially based on an MSSM spectrum generator [26] such that the Higgs enhancement is added afterwards by hand. Considering the importance not only on Higgs observables but the pressing issue of Naturalness this is not sufficient and a dedicated (and publicly available) spectrum generator and study of this type of model is required. It is of course essential that the HEP community has access to a myriad of spectrum generators, particularly those relevant to the Higgs sector, and more ideally several spectrum generators dedicated to the exact same model. This has been incredibly fruitful in studies based on the MSSM and NMSSM. This paper is an important step in that direction as it supplies a full two loop spectrum generator, with the non decoupled D-term Higgs enhancements included in the 2 nd gen.

MSSM DSB
H u H d Figure 1. A picture to represent the quiver module of the electroweak sector for MI as in model and in the self energies for the calculation of the masses. What is particularly exciting and novel regarding our work is that our spectrum generator is on a par with all currently publicly available spectrum generators (such as [26,[28][29][30]) and is the first dedicated spectrum generator to include the non decoupled D-term enhancement to the Higgs sector. Secondly a dedicated study can address the important question of, if by removing the quiver of SU (3), one can then attempt to maximise the Higgs enhancement and achieve a linking field vev, v, below 10 TeV, whilst keeping all gauge couplings perturbative to the Planck scale. In fact we found that this cannot be done. In addition this paper supplies a number of important and new results: it supplies a full derivation of the non decoupled D-term effect to all scalars and not just the Higgs fields, these contributions have so far gone unnoticed in the literature, but may be quite relevent to precisions studies at e + , e − colliders. Further, this implementation includes all RGEs to two loops for all parameters and all anomalous dimensions for all fields, supplying for the first time, the full anatomy of such a setup.
The outline of the paper is as follows: In section 2 we give a first account of our key findings. In section 3 we present the specific realization of a quiver gauge theory that we will be studying in the following sections, discuss its general features and our choice of parameterisation of the soft term boundary conditions. Amusingly, these setups in some sense also model the effects of a truncated theory of gauge fields in an extra dimension and we make a small digression in section 3.5 to discuss this connection. In section 4 we discuss the concrete implementation of our quiver setup, the model's parameter space and the phenomenological constraints it is subject to. Then, in section 5 we present our results on the phenomenology of the model and discuss the features of its particle spectrum. Finally, in section 6 we conclude. The appendix A contains details about the implementation of our framework along with some important relations. This paper is accompanied by a supplementary document containing details with regards to the implementation [27].

Key results
It is useful here to give a first summary of our key results. In this work we will show that • With the setup we adopt it is perfectly possible to obtain a 125.5 GeV Higgs with stops lighter than 2 TeV, gluinos of order 1600 GeV, light electroweakinos and sleptons potentially within the reach of the LHC, all within a GMSB framework (see benchmarks in table 5).
• The correct Higgs mass is obtained through non decoupled D-term contributions in the low-energy lagrangian that lift the Higgs boson mass at tree level, as shown in figure 4.
• These terms also modify the Higgs branching ratios but well within current LHC bounds, and could be probed by the ILC as seen in figure 6.
• The light uncoloured sparticle spectrum is achieved primarily due to the specific supersymmetry breaking mediation mechanism we employ (see figure 7).
These results are obtained by implementing our model into the publicly available package SARAH, which enables us to create a spectrum generator in order to perform the RGE evolution of all model parameters at two loops, from the messenger scale M down to the TeV scale, with GMSB boundary conditions. We have implemented five gauge groups at full one-and two-loop running, plus oneloop self energies, from the GUT or messenger scale, Higgsing and breaking to the diagonal subgroup of 4 gauge groups, while finite shifts and threshold corrections are carefully applied for each degree of freedom.
We finally stress that we have implemented a conservative (precise) formulation of GMSB with full two loop equations for soft masses in the hidden sector. Still at this level of specification a reasonably natural spectrum is obtained, which demonstrates the ease with which much lighter spectra would be obtainable if these high standards were relaxed or some more phenomenological parameterisation adopted. The framework that we have developed quite straightforwardly admits numerous extensions such as inclusion of U (1) kinetic mixing or quivering the SU (3) sector, tasks which are left for future work. All of these remarks will be clarified in the following sections.

An electroweak quiver
In this paper we wish to explore two different quiver models for comparison. The first carries the generic features of non-decoupled D-terms and, in the case of GMSB, suppressed scalar soft masses versus gauginos. The second is a flavourful extension of the first model to achieve lighter stops than the first two generations, which still obeys all anomaly cancellations. A common feature is that we will apply a gauge mediated supersymmetry breaking scenario to both and both are characterised by the scale of supersymmetry breaking, M , Table 1. Gauge superfields of the model. and the vevs of the linking fields 2 , v. In this section we outline these models and their features.

The models and features
Let us consider an electroweak two-site quiver with gauge group The two sites are connected by means of a pair of linking chiral superfieldsL ,L. These superfields will play a crucial role both in the breaking of the enlarged gauge group to the MSSM gauge groups, by obtaining vevs, and in the mediation of supersymmetry breaking effects. Moreover, the setup includes a singlet chiral superfield K, whose role will be clarified shortly, as well as an additional superfield A transforming as the adjoint of SU (2) B that serves the role of giving masses to certain fermionic components of the linking fields upon G A ×G B breaking. Much below the higgsing scale, v, the quiver fields usually decouple and so for phenomenological purposes at low energies the model is essentially MSSM-like with the addition of an effective action for the Higgs potential. It will be useful then to refer to the enlarged gauge groups as regime 1 and the MSSM as regime 2. This paper is based on two models which are as follows: [Model I]: The first model (MI) is a basic example of a quiver model where the MSSM chiral superfields are taken to be charged under site A identically to their charges under the MSSM gauge group and are neutral under site B (see figure 1 and table 2). The superpotential of the MSSM-like matter is given by with i, j, k labelling SU (2) indices, and as this group is pseudo real, the2 φ i is simply ij φ j of the 2 representation φ i of SU (2). The superpotential of the quiver module is given by A model with a similar structure albeit based on a more enlarged gauge group was first introduced in [7]. The general features of this model will be outlined below and unless  2 and table 3). Similar representation assignments have been considered, for example, in [8,9] and then later in [13,14,25] in the framework of models of natural supersymmetry that could potentially further address the flavour problem. The superpotential we use in regime 1 is plus the quiver superpotential Eqn. (3.2). In regime 2 after returning to the MSSM gauge groups, we adopt the MSSM superpotential by Eqn. (3.1).

Gauge symmetry breaking
We now describe certain features of the general setup. The superpotential Eqn. (3.2) gives rise to a scalar potential which when minimized sets a vacuum expectation value for the scalar components of the linking fields of the model. Denoting  in the absence of supersymmetry breaking, we write These break the gauge group G A × G B down to the diagonal subgroup G Diag = SU (2) L × U (1) Y , which are simply the MSSM gauge groups. The symmetry breaking pattern takes the form Now, including soft breaking masses m 2 L , m 2 L for the linking fields, we expand the diagonal complex scalars into real scalar and pseudoscalar components (3.7) The σ's play the role of Goldstone bosons and get eaten by the gauge fields. Minimizing the scalar potential with respect to these fields leads to the tadpole equations, which at tree-level read and a similar expression for φ L4 . The value of the vev v can be obtained by requiring that the tadpoles should vanish. In practice it turns out to be much more convenient to take v as an input parameter and compute the superpotential parameter V 2 from Eqn. (3.8).
As a result of the quiver structure at different renormalisation scales Q the following occur: • Regime 1 is characterised by M ≥ Q > v with the full matter content and gauge groups of the quiver.  • In Regime 2, characterised by v > Q, the vevs of the linking fields break the groups to the diagonal and the MSSM superfields transform under G Diag ≡ G MSSM in the usual way.
• The U (1) gauge bosons B A , B B mix to generate a massless and a massive state B 0 and B m , the massless one being then the U (1) Y boson and the B m being a heavy state. Similarly W i A , W i B mix to form the massless SU (2) L W i 0 gauge bosons as well as three heavy states W i m with the corresponding mixing angles discussed below. The masses of the heavy gauge bosons are simply given by In this setup we have not considered the quiver structure for SU (3). Whilst this is mostly due to practical reasons, given the difficulty of a full and proper implementation of a Higgsed SU (3) as in [10,14], it is also not necessary for our purposes. Indeed naively we sacrifice a GUT completion, but clearly we are expecting our model to be valid only up to the messenger scale in the case of GMSB and anyway it should be quite straightforward to tidy up this setup to restore gauge unification without sacrificing the key results of this work. In this sense our setup is both minimal enough, and yet concrete enough to be a "theoretical simplified model" which captures relevant features of a much larger range of possibilities, for example, it incorporates an extra SU (2) and U (1).

Supersymmetry breaking and soft breaking terms for gauge mediation
In principle the models above may be combined with any supersymmetry breaking scenario, for example mSUGRA, or some more phenomenological parameterisation. The model contains a large number of soft terms. The soft breaking scalar potential reads a, b are flavour indices and i, j, k SU (2) indices, which are lowered with ij . The soft terms for the fermions are In what follows, the RGE evolution of all parameters in this scalar potential will be accounted for at two loops.
In this work we focus on gauge mediation and take the highest scale of the RGE evolution to be M , the characteristic mass scale of supersymmetry breaking which in perturbative models is the messenger scale. As is typical of these perturbative gauge mediated supersymmetry breaking scenarios, we model the supersymmetry breaking sector with a set of messenger fields coupled to a spurion. Such a messenger sector can (and should) be generalised [12], although we will conform to this standard paradigm. The superpotential we use is of the form where X is a spurion with X = M + θ 2 F and Φ,Φ are representative of fundamental and antifundamental messenger fields respectively, charged under SU (3) c and SU (2) B , U (1) B , but not under the A-site electroweak group. This leads to a scale Λ = F/M which may differ for each gauge group so we can in general write Λ 1,2,3 for the three gauge groups. The messenger fields and spurion are integrated out at M to generate the soft terms, the explicit equations for which are supplied below. 3 Here we describe the gauge mediation parameterisation of the above soft terms.

The trilinear, bilinear and linear terms
The trilinear T-terms (or A-terms) are taken to be zero at the messenger scale, including those corresponding to Y K and Y A . The linear soft term L V 2 for the singlet φ K , is also taken to be vanishing at the messenger scale. In GMSB, the bilinear term B µ is expected to be zero at the supersymmetry breaking scale and should be generated by RG running.
In what follows we will use the standard SUSY Les Houches Accord GMSB conventions [26,[28][29][30][31][32][33] according to which tadpole equations are solved for µ and B µ and tan β is given as an input. Above the scale of G A × G B breaking, the β-function for B µ reads at one-loop which results in a large B µ if the g Ai 's are relatively large. The equations are similar below the quiver breaking scale. Note that there are also two-loop contributions in both regimes.

Gaugino soft masses
where x = F/M 2 , Λ = F/M and r refers to the corresponding gauge group. N = (n 5plets + 3n 10plets ) is the messenger index and the function g(x) is the standard function appearing in GMSB gaugino soft masses.
As the messenger sector is not charged under U (1) A and SU (2) A , the corresponding gauginos are taken to be massless at the messenger scale: One might imagine that such a feature could be detrimental to the low energy spectrum. However, the mass matrices of the gauginos are rather complicated including supersymmetric Dirac masses as well as the above soft terms, so this turns out not to be the case: the mass eigenstates result from a combination of the site A and site B gauginos 4 . The Majorana soft masses of the broken theory can be found by identifying the masses of the relevant components of the mixing matrices, at the threshold scale O(v) (see also appendix A.1). As the A-site gauginos do not obtain significant soft masses until the scale v, the RGEs of matter charged under site A will not feel these effects until a scale Q < v. This turns to be advantageous for naturalness as now the threshold scale T scale = v acts as a cutoff to the leading RGE logarithm. Such an effect will become especially important for an SU (3) A × SU (3) B quiver, as then the A-site gaugino would only influence the RGEs between M ew and T scale and have essentially no effect on the RGEs of site A matter above this scale.

Site A scalar soft masses
The scalar soft masses depend on the site under which the corresponding superfields are actually charged, relative to those of the messenger fields. We will always take the messenger fields to be on site B. Fields charged under G A get soft masses at two loops from mediation along the quiver The form of the function s(x, y), associated with gauge mediation along a two site quiver, can be found in [12,23,24] and is given in both analytical and graphical form in appendix A.5. By inspecting figure 10 we can see how the mediation of supersymmetry breaking along the quiver has the effect of reducing the site A scalar soft masses with respect to usual gauge-mediated supersymmetry breaking. In particular, s(x, y) has the limit where f (x) is the usual GMSB formula. This formula interpolates between y → ∞ of GMSB and the suppressed scalar regime as y → 0, where scalars get their leading soft mass at three loops from additional contributions, which arise anyway from RG evolution.

Site B scalar soft masses
Scalar fields charged under G B and SU (3) c receive standard GMSB soft masses according to The quadratic Casimir invariants are C 1 (Y B ) = 3/5Y 2 B for fields with charge Y B under U (1) B whereas C 2 (2) = 3/4 for the linking fields L and C 2 (3) = 2 for the SU (2) B adjoint A field. All coloured scalars in our setup are SU (3) triplets, for which C 3 (3) = 4/3. Note that in the case of site A coloured scalars, the full soft mass is given by the sum of Eqn. (3.16) and the third term of Eqn. (3.18).

The singlet scalar K soft mass
The superfield K is a gauge singlet and its scalar soft mass is vanishing at the messenger scale, It evolves a positive value through Eqn. (A.29) and so does not pose a phenomenological issue, but if one did wish to assign a tree level soft mass, two approaches are possible: it may be interesting to consider that it is not a singlet under some other group or that it couples directly to messenger fields through a term of the form KΦΦ. In this later case it can develop a one-loop soft mass.

Linking scalar soft masses
The linking fields formally get their soft masses from applying Eqn. (3.18) to describe the soft terms for m 2 L and m 2 L . We will however not be using this formula in order to compute the linking field soft masses. Instead, we will promote them to free parameters of the model. The reasons for this choice will be clarified in the following section. To be noted is that in this setup the two linking field masses are equal, m 2 L = m 2 L .

Linking to the MSSM
Below the quiver breaking scale the gauge group and particle content of the model are those of the minimal supersymmetric standard model with gauge groups The gauge couplings between the unbroken and the broken theory are matched as If one of the two gauge couplings is strong, the other should be weak. Then at low energies the diagonal or MSSM gauge coupling will be of the order of whichever is weaker. This is a key feature which allows for these models to lift the tree-level Higgs mass whilst being consistent with perturbative unification. The various gauge couplings of the model are simply related through two rotation angles 5 In appendix A.4 we present additional comments on threshold effects that enter the coupling constants and other parameters calculation below the quiver breaking scale. The angles θ 1 , θ 2 are free parameters of our setup. Varying these amounts to changing the relative strengths between each site and we typically choose the A-sites to be stronger to enhance additional contributions to the Higgs mass, as we will explain below. 5 Here we should stress an important notation subtlety. In all relations applying to the messenger scale as well as in all RGE expressions, the U (1) coupling constants are taken to be SU (5) GUT-normalized, so for example g1 = g1,GUT = 5/3g , with g being the usual Standard Model hypercharge coupling constant.
In all other relations, the GUT normalization is dropped and g1 identical to g . This is done in order to follow the SARAH package conventions.
One of the most interesting features of this class of models is that non-decoupled D-terms may arise [3,34] in the low energy Lagrangian. The real uneaten scalar components of the linking fields appear in both the A and B site scalar D-term potential and when integrated out generate an effective action which includes the terms where It is particularly informative to see how in this class of models, these terms can work to lift the Higgs mass without large radiative corrections. In the MSSM, the one-loop Higgs mass in the limit m A 0 m Z can be written as [35] where M 2 S = mt 1 mt 2 for mt 1 , mt 2 as defined in appendix A.1, and v ew is the electroweak Higgs vev, such that the upper limit on the tree level Higgs mass (m h,0 ) is set by the Z boson mass m z . This expression assumes that the left and right-handed soft masses of the stops are equal. Note that X t = A t − µ cot β, and for convenience the sfermion mixing matrices are provided explicitly in Eqn. (A.12) of appendix A.2. In our case however, there may in principle be a sizeable shift to the Higgs mass at tree level m h,0 , which takes the precise form, Arguably this enhancement is favoured over that of the NMSSM for a simple reason: in the NMSSM typically m 2 h,0 = m 2 z cos 2β + λ 2 v 2 ew sin 2β (3.26) where λ is the coupling between the Higgs singlet and doublet fields appearing in the superpotential term λSH u H d . This creates a tension between wanting a large tan β to enhance the first term, but a small tan β for the second, forcing one to accept very large values of λ. As a result λ ends up non perturbative before the GUT scale. It is the above observation that forms the basis for the construction of natural spectra in the class of models that we examine: the Higgs mass can now be substantially increased already at tree-level when these new contributions become large. Of course this enhancement is completely independent of the method by which supersymmetry breaking effects are transmitted to the MSSM. We have simply chosen GMSB in this paper on the one hand to demonstrate that it is still a natural candidate for supersymmetry breaking mediation, and on the other hand because in our electroweak GMSB quiver the sleptons can be naturally lighter than their coloured counterparts. This potential enhancement of the tree level Higgs mass is only significant in certains areas of the model's full parameter space. Concretely, for this contribution to be sizeable we must have g 2 Ai /g 2 Bi ≥ 1 and However, this mechanism introduces some additional fine-tuning to the theory, since the Higgs mass now receives an additional quadratically divergent contribution at one loop, induced by the linking fields and cut off by m 2 L . This additional fine-tuning should be kept under control in order to not counterbalance the improved naturalness of the model with respect to traditional mGMSB. An estimate of the maximal size of m L can be found following the arguments of [5]: demanding less than 10% additional fine tuning approximately bounds Requiring, for example, ∆ = 0.2 in Eqn. (3.23) allows m 2 L in the 10 6 − 10 8 GeV 2 range and sets an upper value v < 10 5 GeV so that the additional D-terms do not decouple. So in summary we would ideally want v < 10 5 GeV and m L < 10 TeV. Note that v is also bounded from below both by electroweak precision tests and direct searches for new gauge bosons.
The non-decoupled D terms also appear in the tadpole equations modifying the vacuum structure, as well as the Higgs mixing matrices, which may be found in appendix A.2. Additional soft mass terms for all scalars appear in regime 2 of the model at effective one loop, from integrating out the heavy gauge and linking fields [10] which are also implemented into the model and importantly the soft mass parameters are matched across the threshold scale. 6

An extra dimensional digression
Quiver models are naturally related to extra dimensional setups through deconstruction. For early ideas on the topic we refer the interested reader to [36]. The contemporary formulation of the topic was initiated in [6,[37][38][39], whereas for recent work relating to N = 1 see for instance [12,24,40]. It should therefore be expected that these non decoupled D-terms, Eqn. (3.23), have a natural interpretation in terms of extra dimensional models. We swiftly sketch and motivate this relationship which certainly warrants further study on its own. The quiver construction may be related to N = 1 super Yang-Mills (or N = 2 [41][42][43]) in five dimensions [44], which contains a vector multiplet and chiral adjoint V + Φ. Suppose we compactify on four flat dimensions times a small interval of length R. The scalar component of Φ = (Σ + iA 5 ) and in analogy to the quiver, A 5 plays the role of the Goldstone bosons and are eaten to generate the Kaluza-Klein masses such that we may identify 1/R ∼ v of Eqn. (3.7). To obtain the non decoupled D-terms we write the lagrangian in the off-shell formulation The ellipses denote not just the rest of the bulk action but also any terms generated on bulks and boundaries that may be of use such as the boundary terms To see how this action may generate a non decoupled D-term we define Then, there may be bulk or boundary terms of the usual form Integrating out the auxiliary scalar field D gives rise to the D-term scalar potential. The field Σ, the real uneaten scalar degrees of freedom, corresponds to the real uneaten degrees of freedom in the linking fields L,L of the quiver. It is this field Σ, when integrated out which generates the non decoupled D-term Eqn. (3.23). In such a scenario, the m v of the quiver is related to the Kaluza-Klein mass scale m kk which is O(1/R), the effective length scale of the extra dimension. As such, for these terms to be of relevance π/R ≤ m soft . We hope to return to this topic in a further publication, but for now we effectively model this feature with quiver models as they are a more controlled environment which are more amenable to spectrum generators. It is certainly interesting to speculate that as our model has a v of O(10 4 ) GeV, that this corresponds to an "effective" extra dimensional length scale of roughly O(10 −18 ) cm.

Tools and Observables
In order to study the low-energy phenomenology of our setup in a consistent manner, it is necessary to perform the RGE evolution of all couplings and mass parameters from the highest energy scale of the theory down to the TeV scale, properly imposing all boundary conditions. Here we describe the construction of a tailor-made spectrum generator for the quiver model. We further discuss the parameter space we adopt as well as the constraints it is subject to.

Implementing a quiver framework for phenomenological studies
In order to perform the RGE evolution of the models' parameters and masses and compute the resulting low-energy particle spectra, we implemented the two model variants into the publicly available Mathematica package SARAH 3.3 [17][18][19]. SARAH is a "spectrum generator generator", which includes a library of models that may communicate with HEP tools that are widely used in most phenomenological studies [30,33,[45][46][47][48][49][50][51][52][53][54]. In particular, SARAH performs the task of generating Fortran routines compatible with the SPheno spectrum generator [18].
In order to implement our model, we have used the possiblity offered by the package to implement and link two different "regimes". These regimes correspond to those introduced in section 3.2, each being characterised by a set of gauge groups, a particle content and a superpotential that need to be specified. Regime 1 includes, for both models MI and MII, the full G A × G B × SU (3) c gauge group along with the full quiver particle content, while the superpotential is given by Eqns. 2) for MII. In both cases in regime 2 we have the MSSM, which we supplement with an effective action to account for and study the effects of the non-decoupled D-terms, which are of crucial importance, as discussed in section 3.4. These terms are properly included in all loop calculations, self energies, branching ratios and vertices of regime two. We would have preferred to implement our model using a single regime such that these terms would be automatically generated, and as there are other terms for other fields, however we found it was not practical to integrate out so many fields in full. Moreover it was preferable to include an additional regime with the MSSM gauge and matter configuration to ease communication of SPheno with packages such as HiggsBounds, used to check the compatibility of the model with experimental constraints as discussed below.
In addition to these ingredients, we need to specify on the one hand boundary conditions for all soft parameters at the messenger scale M and on the other hand matching conditions for the parameters of the two regimes at the G A × G B breaking scale which we typically take to coincide with v, the linking field vev. The boundary conditions are applied according to the discussion and relations given in section 3.3, while the matching conditions follow the lines described in section 3.4.
The renormalization group equations for both regimes are then calculated by SARAH at two loops and appropriate Fortran routines are generated that can then be taken over by SPheno to perform the numerical analysis. Concretely, the implementation includes full one-and two-loop RGEs for five gauge groups, mixing matrices for all fields in the quiver and the MSSM including associating goldstones with massive gauge bosons and gauge fixings, full two-loop RGEs for the vev of the linking fields themselves, two-loop RGEs for Bµ, loop-level solutions for the tadpole equations, one-and two-loop anomalous dimensions for all fields and two-loop RGEs for all soft breaking parameters, linear, bilinear and trilinear. The MSSM particle masses are computed at one loop, however the full two-loop corrections to the Higgs mass are implemented in SPheno following the calculation in Refs. [55][56][57][58].
All in all, the RGE evolution is described by three energy scales Following standard practice, the highest energy scale of the theory is taken to be the messenger scale M , where all boundary conditions resulting from the quiver structure, including exact formulae for GMSB soft masses have been implemented and imposed. Again as usual, the running ends at the electroweak scale, which is used as an input scale for the MSSM parameters. The intermediate mass scale T scale is associated with the quiver breaking scale as it separates the two regimes, and at this scale appropriate matching boundary conditions are applied including finite shifts that result from integrating out the heavy fields of the theory (see appendix A.4). We choose it to be equal to the vev of the linking fields, T scale = v. All soft terms of regime 2 are matched to regime 1 as described in section 3.4 and for the soft masses of the winos and binos in appendix A.1. Note that the gluino finite shifts are also accounted for as given in appendix A.4.
It is important to stress that the high-scale boundary conditions themselves may be seen as being separable from the model (the matter content, gauge groups and superpotential) and may be changed with ease, if one wished to explore, for example, different supersymmetry breaking scenarios.
The implementation detailed above, i.e. the construction of a tailor-made spectrum generator for our model, allows us to create a model file for SPheno, which further permits us to study a quiver model in a complete manner, as we can study the influence of RG effects of all the gauge groups and matter content in the highest regime to the low energy spectrum, at the two-loop order. On the practical level, this model can also serve as a first step for the implementation of more complete or complicated setups, such as a model including an additional SU (3). It would further be trivial to change the representation assignments of the Higgs fields in order to study chiral non-decoupled D-terms [59] or flavour models, with the same precision. However, we point out that the implementation of the electroweak only quiver studied here leads to an interesting phenomenology in its own right, allowing for naturally (although moderately) heavier squarks relative to sleptons.

Parameter space and constraints
It is now useful to describe the process through which we choose our parameter space and the regions we will study in the next section. We moreover describe some preliminary findings that could be of interest for model-building purposes.

Choosing the parameter space
The electroweak quiver we consider can be described by a basic set of six parameters where M is the messenger scale, the SUSY breaking scale Λ = F/M , F being the SUSY breaking F-term, tan β is the Higgs vev ratio and θ 1 , θ 2 are the mixing angles between sites A and B for U (1) and SU (2) respectively.
As an initial step, we performed extended scans over large regions of parameter space for MI, imposing the full set of GMSB boundary conditions described in section 3.3, including the exact relations Eqn. (3.18) for the linking field soft masses. We focused in particular on regions where v 40 TeV, where according to the discussion of section 3.4 the non-decoupled D-terms should be most efficient in lifting the tree-level Higgs boson mass. A first finding of these searches is that we could not find viable points when the linking field vev was much below 10 TeV, as often here either the A-site couplings become non-perturbative before the messenger scale, the electroweak vacuum becomes unstable or the RGE code simply wouldn't converge for the numerical precision requirements imposed (a relative error of 0.5%). Where none of these issues occur, the values of Λ are typically low and close to v, such that the linking fields soft masses m L are too low for the D-terms to have an important impact on the Higgs mass.
Motivated by the perturbativity issues, we implemented MII where we expect that by removing some of the matter fields from site A the RGE running will be reduced [14]. We found that although the situation does improve, it still seems to be quite difficult to achieve substantial contributions to the Higgs tree-level mass from the non-decoupled D-terms due to the fact that m 2 L is again driven too low. These results lead us to slightly enlarge our parameter space by promoting the linking field squared soft mass m 2 L to be a free parameter instead of being given by Eqn. (3.18). This is interesting from a theoretical perspective, as it motivates pursuing models that might provide the additional contribution to m 2 L needed in order to achieve a substantial D-term contribution to the Higgs mass. For example, it would be interesting to study whether this can be realized in extensions to the model including an additional SU (3) or U (1) kinetic mixing. Within the scope of our work, the choice to make m 2 L free can be seen as a phenomenological parametrization of the linking field soft masses along the lines of similar choices made in many supersymmetry breaking mediation schemes.
With this small modification, we find that it is indeed perfectly possible to achieve the required D-term size in order to reproduce the observed Higgs mass while keeping the stop masses well below 2 TeV. Furthermore as expected, the mediation of SUSY breaking along the quiver acts as a suppression mechanism for the uncoloured sparticle masses, yielding electroweakinos and sleptons lying roughly in the range [400, 1000] GeV, which is on the boundary of being within the LHC reach [60,61]. At this point, due to the differing bounds on coloured and non-coloured sparticles at the LHC, we introduce a second modification to the original setup that consists of dissociating the scale Λ 3 from Λ 1,2 . Note however that this is a minor modification as the two scales will not differ by orders of magnitude but only by O(1) multiplicative factors. We will see that this setup allows for a rich phenomenology with interesting features.

Constraints
We carried out extensive scans of the parameter space described in the previous paragraph within generous intervals. We are interested in areas of parameter space which are characterised by low values of v and moderate splittings between v and m L , such that the additional D-terms do not decouple from the low-energy theory and the uncoloured scalars are light. In what follows, we will therefore present results that concern a subregion of the parameter space that meets a series of requirements.
First, we wish to obtain a Higgs mass lying in the range [122.5, 128.5] GeV. This interval envelops on the one hand the experimental uncertainty in the Higgs mass measurement [1,2,62], while being sufficiently generous to account for uncertainties in the theoretical mass spectrum determination [63,64] 7 . For naturalness reasons, we require this value for the Higgs mass to be achieved for stop masses as low as possible. The stop mass is governed by Λ 3 , which also controls the masses of the lower generation squarks and the gluino. Strong exclusion limits on these masses arise from ATLAS and CMS null searches for jets plus missing energy, e.g. mg > 1600 GeV for mq 1,2 > 2000 GeV [66,67]. We are guided by these bounds in choosing a lower limit for Λ 3 . Note that in this work, we choose not to quantify the amount of fine-tuning for each point of the models we study, which constitutes a work in its own right involving numerous subtleties (see for example the recent discussion in [68]). It is however at least clear that qualitively, having stops lighter than benchmark minimal GMSB improves the relative naturalness of the model, and this motivates our choice of upper limit on Λ 3 . At the same time, according to the comments made in section 3.4, we should avoid reintroducing excessive fine-tuning via the non-decoupled D-terms. Moreover, in order for the setup to be realistic, we must satisfy the condition m 2 L < m 2 v , but not approach the limit m 2 L m 2 v where the quiver-induced D-terms decouple. These requirements lead us to choose m 2 L within the range [10 7 , 10 8 ] GeV 2 . Also note that as mentioned above, for very low values of v SPheno faces convergence issues. The parameter Λ 1,2 is mainly subject to constraints from searches for charged sleptons and charginos at LEP, i.e. 92 GeV for charginos degenerate with the lightest neutralino, and 103.5 GeV otherwise [69]. Lower limits on sleptons staus and sneutrinos of 68 and 51 GeV respectively were also obtained at LEP [70][71][72]. Finally, given that the non-decoupled D-terms contribute a shift to the Higgs mass as m Z does, i.e. with a factor cos 2β (Eqn. (3.25)), as opposed to the factor sin 2β in the NMSSM (Eqn. (3.26)), we explore a rather standard MSSM-like range for tan β.
From our numerical analysis we find that this set of requirements is satisfied by adopting the following parameter value ranges   Table 4. Low-energy observable constraints imposed in our analysis.
We have moreover chosen signµ = +1, a low value for the messenger index, N = 1, and a fixed common value for the A and K field Yukawa couplings Y A = Y K = 0.8. Apart from the theoretical and experimental constraints so far mentioned, the lowenergy spectrum is subject to further bounds. In the Higgs sector, in addition to obtaining the lightest Higgs boson mass within the observed region, we must ensure that its properties and decay modes comply to current LHC observations. As an example, it is known that enhancing the Higgs mass through non-decoupled D-terms enhances simultaneously the Higgs boson couplings to down-type quarks [5,15,59]. In order to test whether the Higgs sector is compatible with the constraints coming from the LHC and the TeVatron, we have linked SPheno to HiggsBounds-4.0.0 [51]. Taking our analysis a step further, we have also linked SPheno to HiggsSignals-1.0.0 [73], which allows us to test in particular whether the lightest Higgs boson properties are in agreement with all relevant existing mass and signal strength measurements from the LHC and TeVatron.
Finally, we use the in-built functionalities of SPheno in order to apply a set of necessary low-energy constraints, all of which are taken at 3σ: the SUSY contributions to the muon anomalous magnetic moment δa µ and the branching ratios BR(B s → X s γ) and BR(B s → µ + µ − ) and, due to the presence of relatively light sfermions in our spectra, the ρ parameter. The allowed ranges used for these observables are shown in table 4, where theoretical uncertainties and experimental errors are added in quadrature.

Results
Having described our model, and how it is implemented in SARAH, we turn to study the low scale spectrum, which we find has several interesting features. Examples of complementary representative points, one for MI and two for MII, are given in  Table 5. Mass spectra of three example points for MI and MII, along with the associated input parameters as defined in section 4.2. Note that mf 12,L/R and mf 12,L/R are the masses of the lower and third generation left/right-handed squarks and sleptons, and that the sneutrino and sbottom masses can be inferred via mν i ∼ mf i,L , mb L ∼ mt L .
In table 5 we observe that particles of the electroweak sector can be substantially lighter than those of the coloured sector. This arises due to the quiver structure of the model, as explained in section 3.3, which provides a suppression factor s(x, y) for the non-coloured scalar masses, for details see appendix A.5. The suppression is further enhanced by the fact that we have chosen to study the range of parameter space where Λ 1,2 < Λ 3 . Therefore it is possible for the masses of electroweakinos, sleptons and the heavy Higgs bosons to lie well below 1 TeV. One observes in table 5 that this results in the next-to-lightest supersymmetric particles (NLSPs), being either a neutralino, as in MI and MIIb, or stau, as in MIIa. A sneutrino NLSP is also possible as will be discussed in detail later in this section. As we only consider a single SU (3) gauge group at the high scale, the masses of the coloured sparticles do not experience this suppression. This means that the coloured sector lies in general between 1.5 and 2.5 TeV, the stops being the lightest squarks. However, in MII, a splitting is generated between the left and right stop soft masses, for reasons discussed below, as shown in point MIIb of table 5. In the following we will study the spectra of these models in terms of their compatibility with the current experimental constraints described in the previous section and the prospects for detecting signs of TeV-scale sparticles in the near future.

The Higgs mass and couplings
As the non-decoupled D-terms lift the tree-level Higgs mass, as in Eqn. (3.25), here we investigate the range of stop masses in these models for which m h lies in the desired range, and how the stop contribution compares to that of the non-decoupled D-terms. The same non-decoupled D-terms can affect the couplings of the Higgs, so we further investigate these couplings in light of current and future experimental measurements.

The Higgs mass
As mentioned in section 4.2, we have chosen Λ 3 such that the masses of gluino and the first and second generation squarks lie above the LHC exclusion limits. In MI, this translates into the stop masses being close to 2 TeV, which means that the shift in the tree-level Higgs mass required in order to obtain m h ∼ 125.5 GeV is small, and the required soft mass of the linking field remains below 10 TeV. The situation is fairly similar in MII, however we find that there is a slight tendency for a splitting to arise between the left and right-handed stop soft masses, due to the RGEs driving the left-handed mass downwards, and the right-handed mass upwards. This can be understood in terms of the differences between the RGE equations for the two models, where in MII above the quiver-breaking scale the Higgs soft masses are only affected by the third generation squarks, whereas for MI the Higgs soft mass RGE equations contain all generations. This results in a larger splitting between the up and down type Higgs soft masses which further generates a larger splitting between the left and right handed stop. The distribution of the masses of the light and heavy stops i.e. mt 1 and mt 2 as defined in appendix A.1 for the two models are displayed in figure 3. Here the allowed points are shown in yellow and those points excluded by the various constraints described in section 4.2 in grey. We clearly observe that for MII a larger splitting between the stops is possible, and the lighter stop may be as light as 400 GeV, as seen in the benchmark point MIIb in table 5.
In figure 4 we have plotted the Higgs mass as a function of tan β for the two variants of our model. Here the bright red points respect mt 1 < 2 TeV and all constraints imposed, the pale red points only comply with the low-energy constraints and the grey points are excluded. The full two-loop corrections to the Higgs mass are implemented in SPheno following the calculation in Refs. [55][56][57][58]. We conclude that a Higgs mass within the limits ∼ 125.5 ± 3.0 GeV is achievable in both MI and MII. Note that the larger range in m h for MII can be explained by the larger range in stop masses. Indeed, when the left and right handed stop soft masses are not equal, as shown in figure 3 for MII, the simplified expression for the one-loop Higgs mass given in Eqn. (3.24) is no longer valid. An additional correction must be added to Eqn. (3.24) of the form [76] which for the case m 2 in MII induces an enhancement to the Higgs mass of around 1-2 GeV. Note that the sfermion mixing matrix is defined in appendix A.2. As the bright red points correspond to mt 1 < 2 TeV, this further demonstrates that the effect of the non-decoupled D-term seems to reduce the fine tuning by allowing lighter stops than in standard GMSB.
To make the distinction between the non-decoupled D-term and radiative, i.e. stop sector, contributions to the Higgs mass clearer, we compare the tree-level result m h,0 to the full two-loop result m h,2 in figure 5. Here the bright blue points respect mt 1 < 2 TeV and all constraints imposed, the pale blue points only satisfy low-energy restrictions and the grey points are excluded. As opposed to the mGMSB result, where m h,0 is bounded by m Z , here we observe that a shift of up to 10 GeV is possible for both MI and MII, while keeping m L < 10 TeV. This in turn means that the contribution of the radiative corrections required to achieve m h ∼ 125.5 GeV is diminished, rendering the model more natural. Interestingly, the splitting of the stops observed in MII results in a distinct difference between the two plots in figure 5, which can be understood from Eqn. (5.1). Despite the fact that the range in Λ 3 is the same for both MI and MII, the stop splitting enhances the size of the radiative corrections, resulting in a smaller shift in the tree level Higgs mass required to obtain a value of m h in agreement with experiment.

The Higgs couplings
Since the discovery of the Higgs boson, not only has its mass been used to discriminate between supersymmetric models but also its couplings, see e.g. [5,77,78]. The deviation of these couplings from the SM can be parameterised via the set of ratios r i , for i = b, γ, g etc, where The r i are further related to the signal strengths µ i normally quoted by ATLAS and CMS, see e.g. Refs. [79,80]. Note that the errors on the measured signal strengths are still too large to make detailed interpretations about the potential underlying SUSY model, and at present the µ i are all SM compatible, and therefore we will not tackle a precise calculation of the various signal strengths in this work. As mentioned, we however make sure that the lightest Higgs boson signal strengths are in agreement with the existing LHC and TeVatron measurements within 3σ, employing the HiggsSignals code. Recent studies by both ATLAS and CMS [60,61] have found that with a luminosity of 300 fb 3), a precision of 1.8% is quoted in Ref. [81]. In our model, the non-decoupled D-terms result in a tree level contribution to the Higgs coupling to down-type fermions. The ratio r b (= r τ at tree level) takes the form where α is the angle between the two Higgs doublets in the MSSM, and is defined by [14] tan 2α = m 2 A 0 cos 2β + m 2 h,0 Here m A 0 is the pseudoscalar Higgs mass of the MSSM and m h,0 is the tree-level Higgs mass given in Eqn. (3.25). We plot the Higgs mass as a function of r b in figure 6, where again the bright red points respect mt 1 < 2 TeV and all constraints imposed, the pale red points only satisfy low energy constraints (i.e. they do not comply with our requirements in the Higgs sector) and the grey points are excluded. We find that only a ∼ 2% change in r b/τ is required for MI, and a ∼ 4% change for MII in order to obtain a Higgs mass of 125.5 GeV, with mt 1 < 2 TeV. Note that in our model, the enhanced coupling to down-type fermions results in a suppression of the signal strength µ γ [5,59], which was not favoured by initial measurements at the LHC [82]. However, as data has collected, the results for µ γ appear more and more SM-like [79,80]. As in this work the tree-level shift in the Higgs mass only needs to be under 10 GeV, we consider small values of ∆ 1,2 , for which the deviation in the coupling of the Higgs to down-type sfermions are well within the current LHC bounds (see e.g. Ref. [79]) as shown in figure 6. Such deviations should start to become detectable at the √ s = 500 GeV ILC.

Sparticle searches at the LHC
As mentioned in section 4.2, the choice of Λ 3 ensures that the masses of the gluino and lower generation squarks approximately respect the limits from direct searches at the LHC. On the other hand, as the scale of the electroweak sector is set by Λ 1,2 , by allowing Λ 1,2 < Λ 3 we explore a range of parameter space for which the electroweak sector has a greater chance of being observed at the LHC. Further, as discussed previously, the quiver structure means that at the high scale the masses of uncoloured scalar particles lying on site A are suppressed, which can in particular result in light higgsinos or sleptons compared to minimal GMSB.
The phenomenology of the model depends decisively on the nature of the NLSP, as this decides which SM particle is present in the final state along with the gravitinoG. ATLAS and CMS have recently made much progress on constraining gauge mediated models, where they study final states containing missing transverse energy (E miss T ) due to the gravitino (G) escaping the detector. Bino-like NLSPs decay viaχ 1 →Gγ, such that the signature is γγ + E miss T , along with additional jets depending on whether the production process is gg orχ 0 1χ 0 1 [83]. When higgsino-like, the NLSP instead decays to a Higgs which can be detected via b jets, and a mixed higgsinos-bino NLSP can be searched for via a γbb + E miss T signature [84]. For stau or sneutrino NLSPs the τ or ν must be searched for in the final state. where the colours indicate the type of NLSP as shown in the legend. The LEP exclusion limit for the case of theτ andν NLSP is clearly marked, whereas the limit for mχ0 1 is given by the y-axis. The grey points are excluded by experimental constraints as described in the text In order to determine which experimental searches are relevant for these models, in figure 7 we examine the region of the mg-m NLSP plane accessed by our scans, indicating the type of NLSP for each point, which we find may be the neutralino, stau or sneutrino. The LEP exclusion limits (see section 4.2) for both theτ andν NLSP are clearly marked, whereas the limit for mχ0 1 is given by the y-axis. We find that for MII there are allowed points for which either the sneutrino or stau are the NLSP, however for MI no such points were found. In MI, the generations of sfermions are treated equally, such that the staus lie close to the other sleptons. On the other hand in the case of MII, as mentioned earlier, above symmetry breaking the third generation sfermions are on site A, whereas the lower generation ones on site B. This has the result that, as in the stop sector, the left-handed stau soft mass may be lower than the right-handed one, such that a sneutrino NLSP is possible. Therefore, although points in both models were found where the NLSP is the stau or sneutrino, only in MII do points survive the demanding constraints imposed on the Higgs mass and couplings due to measurements at the LHC, as illustrated in figure 8.
As the lightest neutralinoχ 0 1 appears to be the favoured candidate for the NLSP, it is interesting to explore its composition as this will enlighten us as to which decay modes are preferred. Therefore in figure 8, we show the lightest neutralino mass, mχ0 1 , as a function of M 1 . From this plot one can deduce whetherχ 0 1 is higgsino or bino-like respectively, depending on the higgsino and bino masses, approximately given by µ and M 1 respectively. The ubiquitous blue points indicate µ > M 1 whereas the more rare red points, which are even absent for MI, show µ < M 1 . The grey points are excluded by Higgs and low-energy constraints, and the horizontal line demarcates the LEP-excluded region for mχ0 1 . In both MI and MII, the higgsino is rarely lighter than the bino, such that the NLSP is mostly bino-like or mixed bino-higgsino, while both can easily be below 300 GeV. This is interesting in light of the fact that experiments are sensitive to the nature of the neutralino NLSP, and the searches would therefore involve photons and/or Higgs bosons and missing transverse energy. Note that the feature of µ being light results in the model being less fine-tuned.
The experimental search strategy is not only dependent on the SM particle in the final state, but also on the decay length of the NLSP, cτ . This can be approximated by [85], where we take F = ΛM . In the region of parameter space considered in this paper, the NLSP decays within the detector. For the case of the neutralino NLSP decaying to a photon and gravitino, which is the prevalent case in both MI and MII, the excellent time measurement of the electromagnetic calorimeter in both ATLAS and CMS means that the time of arrival of the photon can be measured. If the NLSP decays immediately, i.e. if cτ < 10 −4 m, then the decay is characterised as prompt, but otherwise it is non-prompt and it may be possible to deduce its decay length [86,87]. We therefore show the decay length of the NLSP in figure 9, as a function of the mass m NLSP . From this figure we can confirm that the lightest neutralino NLSP may undergo both promp or non-prompt decays to the photon and the gravitino, although in MI fewer points survive for which the neutralino decays promptly. The most important channels for these models at the LHC are therefore searches for photons and missing transverse energy, where the photons may be prompt or non-prompt. Here the dominant production would be electroweak, as our choice of Λ 3 is such that the gluon and squark pair production is suppressed. Studies so far by CMS have concentrated on strong production of the bino-like NLSP [87,88], whereas ATLAS has considered the diphoton and missing transverse energy final state from direct electroweakino production, for the case of both promptly decaying [83] and long-lived neutralinos [86]. However, the bounds obtained by ATLAS are not directly applicable here, as they are presented for a specific point SPS8 [89], where the neutralino NLSP decays predominantly to the photon and gravitino which is not necessarily the case in our model, especially due to the fact that the higgsino is often light. Therefore the bounds on final states including Higgs bosons, studied in Refs. [90,91], must be taken into account. In order to constrain MII, one must further consider the stau and sneutrino NLSP, such that final states involving τ s and missing transverse energy are of interest [91,92]. It would be of great interest to combine all these excluded cross-sections to extract precise exclusion bounds (along the lines of e.g. Ref. [93]), but this is beyond the scope of this paper. It nonetheless seems that for gauge mediated models an interesting region of parameter space is starting to be probed, and we eagerly await further results.

Conclusions and discussion
In this paper we have examined phenomenological aspects of a minimal gauge extension of the MSSM containing two copies of the electroweak gauge group. Using state-of-the art publicly available HEP tools, we have computed the two-loop RG equations for all parameters of two variants MI and MII of this basic setup, characterised by different assignments for the representation of the MSSM chiral superfields. Although the model may be amenable to any set of soft term boundary conditions, we have chosen to work within the framework of gauge mediatied supersymmetry breaking and we performed the RG evolution from the messenger scale down to the electroweak scale in order to compute the sparticle spectrum. We calculated the corresponding sparticle masses at one loop, the predicted Higgs boson mass at two loops and further investigated the predictions of the models MI and MII for the most relevant experimental observables.
As the extended gauge structure results in non-decoupled D-terms which increase the tree-level Higgs mass, the resulting spectrum can be more natural than in minimal GMSB. We further found that in order to be in agreement with Higgs mass constraints, while keeping the stop masses below 2 TeV, one must generate sufficiently large ∆'s (Eqn. (3.23)). This requires the linking soft mass m L to be O(3−10) times higher than expected from the exact GMSB boundary conditions, which indicates a useful direction in which to extend this work on a theoretical level.
We also found that both variants MI and MII of the model would have interesting phenomenological consequences at colliders, since they could be probed either indirectly through Higgs couplings measurements or via direct sparticle production. The Higgs couplings to down-type fermions deviate from the SM due to the ∆'s by r b 6% for MI and r b 4% for MII, and although at present this is well within experimental limits, such deviations should be measurable at a linear collider. As we have focused on the region of parameter space where the coloured sector is ∼ 2 TeV in order to evade bounds on squark and gluino production from the LHC, while the electroweak sector is kept below 1 TeV, the most promising production channel at the LHC is the direct production of electroweakinos. As the predominant NLSP is the bino-like neutralino, diphoton and missing transverse energy searches offer the most promising search perspectives, though for MII, the NLSP is not limited to the bino such that finals states containing τ or h and missing transverse energy are also relevant. As the LHC exclusions are presented in terms of specific models, we are therefore keen to reinterpret these in order to understand how these bounds translate in the case of our model.
There are a number of ways in which this work may be extended. A first step, as mentioned earlier, would be to determine whether larger ∆'s can be realized without making m L a free parameter by including U (1) kinetic mixing or, more ambitiously, an additional SU (3). This could be achieved by means of the tools we have developed with the help of the publicly available package SARAH. It would further be of interest to study related models of flavour, or models with chiral non-decoupled D-terms at the same level of precision. By moving Higgs fields or generations onto different quiver sites, such models are relatively straightforward to implement in our setup. It would also be ideal to construct single regime models, in particular for cases where the phenomenology of additional light states may become relevant. Furthermore, not only are there SO(N ) and Sp(N ) gauge extensions, but even more general quiver constructions, such as those with 3 or more sites [12,94], may also be implemented in full. The study of these models and their GUT completion is also a noble task from the perspective of string phenomenology which has so far been rather neglected. A Some further comments on the implementation In this section we include some useful comments on the implementation of the model.

A.1 Fermion mixing and soft term matching
In regime 1 there are many fermions that may mix together. The uncharged fermions are the diagonal linking fermions, the light and massive bino-type fermions, the uncharged light and massive wino type fermions and finally the diagonal fermion of the A adjoint superfield as well as the K singlet fermion: The lightest two of these states become the MSSM bino and uncharged wino. The mass matrices may be found in the .pdf of the model file in SARAH.
The charged fermions that mix together are the off-diagonal linking fermions, the charged wino-type light and massive gauginos and the off-diagonal A superfield fermions: the lightest two of these become the MSSM charged winos. The rest of the states, both scalar and fermion, of the linking fields K, A, L,L are integrated out at the threshold scale between the first and second regime.

A.2 MSSM Higgs and sfermion mixing matrices
The non-decoupled D-terms of section 3.4 appear in the tadpole equations as well as the Higgs mixing matrices. For the real components (φ d , φ u ) , (φ d , φ u ) we get where g 1 = g and for convenience we use g 2 12 = g 2 1 + g 2 2 + g 2 1 ∆ 2 1 + g 2 2 ∆ 2 2 .
while for for pseudo-scalar Higgses (σ d , σ u ) , (σ d , σ u ) the relevant expressions are The mass matrix for the charged Higgses where we have used the abbreviationĝ 2 12 = − g 2 1 − g 2 1 ∆ 2 1 + g 2 2 + g 2 2 ∆ 2 2 , and in all the above expressions, the ξ-terms are gauge-dependent contributions (and we work in Feynman gauge throughout this paper).
For completeness, we also include the mixing matrix Mf of a generic sfermionf which may be a squark or charged slepton. This matrix takes the form: for s w = sin θ W where θ W is the Weinberg weak mixing angle, and we make use of the abbreviationM 2 Z ≡ m 2 Z cos 2β. The off-diagonal element X f is defined in terms of the trilinear coupling A f via where cot β applies for the up-type quarks, f = u, c, t, and tan β applies for the down-type fermions, f = d, s, b, e, µ, τ . Note that m f , Q f and I f 3 are the mass, charge and isospin projection of the fermion f , respectively. On diagonalization of this matrix one obtains the light and heavy sfermion masses mf 1 and mf 2 .

A.3 Renormalisation group equations
We evolved the model down from the messenger scale M , to a threshold scale T scale , which is associated with the masses of the linking field states O(m 2 v ). The two-loop renormalisation group equations were used in both regimes 1 and 2, along with one-loop finite energy corrections and two-loop anomalous dimensions.
The beta functions of the gauge couplings of the first regime at one loop are Let us also track the top Yukawas using the "only third family approximation", (A. 16) In the second regime these become In the first regime we find the trilinear A t coupling to be whereas in the MSSM Let us also look at how the gauginos obtain soft mases. The one-loop beta functions for the B-site gaugino soft masses are given by For the A-site gauginos they are given by Even though the A-site gaugino masses are vanishing at the messenger scale M the two loop-contributions which typically depend on all the other gaugino soft mases, feedback into the one-loop contributions. Finally the supersymmetric Dirac masses associated with the quiver structure will lift their mass eigenstates. The two-loop equations are given by The soft masses for the quiver module run too: where by σ 1,3 we denote the soft mass combination The vev of the linking fields runs as well: Further equations may be found in the pdf for this model, including all anomalous dimensions and beta functions. At the electroweak scale one finds and which are after all used to minimise the electroweak tadpole equations.

A.4 Threshold effects
We integrate out various states at the threshold between the two regimes. These include the fermions discussed above as well as all linking scalars and scalars of K and A. To implement this correctly, with two-loop RGEs, we edited by hand the SPheno code to properly account for the finite shifts and the mass orderings of the particles integrated out, given by to account for our particular matter content and mg → mg 1 ± The numbers are associated to a gauge boson, weyl fermion, complex and real scalar respectively. T (R) is the index (half the Dynkin index I(R)), such that T ( ) = 1/2, T (Adj) = N c ). It is divided by the dimension of the representation D(R) as each shift is for the component of the field and not the full multiplet, in SARAH. The massive gauge fields are integrated out either on their own or by including them with the finite shifts of the real eaten goldstone modes.

A.5 Soft mass function
In the SARAH package we included, in the form of a Fortran function, the generalisation of the usual mGMSB formula f (x) to the case of a two site quiver model: s(x, y), as pictured in figure 10. The analytic expression is given by s(x, y) = 1 2x 2 s 0 + s 1 + s 2 y 2 + s 3 + s 4 + s 5 + (x → −x) , (A. 38) where s 0 = 2(1 + x) log(1 + x) − 2Li 2 x 1 + x + 1 2 Li 2 2x 1 + x , 1 + x y 2 , The function h is given by the integral The dilogarithm is defined as Li 2 (x) = − So as not to introduce IR divergent pieces it is best to first evaluate terms with massless propagators. In that case the the function h simplifies to h(0, b) = 1 + Li 2 (1 − b) and has a symmetry h(b, 0)=h(0, b). For four massive poles, the analytic expression for h is used in SARAH For illustration, in figure 10 we depict s(x, y) as a function of the parameter x for some indicative, fixed values of y.

A.6 Generalising non decoupled D-terms
Previously non decoupled D-terms have been used to explore both vector-like and Chiral Dterms for the MSSM Higgses [3,59]. It is actually the case that this effective action effects all fields charged under the relevant symmetries and as a result, there will be effective terms for squarks and sleptons too. This point has so far not been mentioned in the literature. We therefore supply a more general derivation for the two site quiver, whose main result is Eqn. (A.52), although we only include the Higgs contributions in our study. It may also be extend to the three site case. For two abelian or non-abelian gauge groups G A × G B that break to the diagonal, one may write canonical kinetic terms for Chiral superfields charged under only site A or site B: After breaking to the diagonal, there is a massless and massive vector multiplet to parameterise the soft breaking fermion χ, the real uneaten scalar, and the Majorana soft mass for λ, respectively. In terms of standard current multiplets satisfying D 2 J = 0, the Kähler potential may be written to leading order in V H as where J a/b are the currents that contain all the fields charged under site a or site b: with the leading term being the scalar current J A = i φ † i T A φ i . The diagonal gauge group coupling is g d = g SM . The effective lagrangian after integrating out V H will then be of the form where O is the effective super operator