# Supersymmetry Breaking in a Minimal Anomalous Extension of the MSSM

Authors: A. Lionetto; A. Racioppi
Journal: ISRN High Energy Physics
 (2012)

---

## Abstract

We study a supersymmetry breaking mechanism in the context of a minimal
anomalous extension of the MSSM. The anomaly cancellation mechanism is achieved
through suitable counterterms in the effective action, that is, the Green-Schwarz terms.
We assume that the standard MSSM superpotential is perturbatively realized; that is, all terms allowed by gauge symmetries except for theμ-term which has a nonperturbative origin. The presence of this term is expected in many intersecting D-brane models which can be considered as the ultraviolet completion of our model.
We show how soft supersymmetry breaking terms arise in this framework, and we
study the effect of some phenomenological constraints on this scenario.


---

## Body

## 1. Introduction

The LHC era has begun and the high-energy physics community is analyzing and discussing the first results. One of the key goals of LHC, beside shedding light on the electroweak (EW) symmetry breaking sector of the standard model (SM), is to find some signature of physics beyond the SM. Supersymmetric particles and extra neutral gauge bosonsZ' are widely studied examples of such signatures. A large class of phenomenological and string models aiming to describe the low-energy physics accessible to LHC predict the existence of additional abelian U(1) gauge groups as well as N=1 supersymmetry softly broken roughly at the TeV scale. In particular in string theory the presence of extra anomalous U(1)’s seems ubiquitous. D-brane models in orientifold vacua contain several abelian factors, and they are typically anomalous [1–57]. In [58] we studied a string-inspired extension of the (minimal supersymmetric SM) MSSM with an additional anomalous U(1) (see [59–67] for other anomalous U(1) extensions of the SM and see [68–71] for extensions of the MSSM). The term anomalous refers to the peculiar mechanism of gauge anomaly cancellation [72–74] which does not rely on the fermion charges but rather on the presence of suitable counterterms in the effective action. These terms are usually dubbed as Green-Schwarz (GS) [59–67, 75–77] and generalized Chern-Simons (GCS) [78–85]. They can be considered as the low-energy remnants of the higher-dimensional anomaly cancellation mechanism in string theory. In our model we assumed the usual MSSM superpotential and soft supersymmetry breaking terms allowed by the symmetries (the well-known result [86]). In this paper we address the question of the origin of the latter in the context of a global supersymmetry breaking mechanism. This means that we do not rely on a supergravity origin of the soft terms but rather on a local setup based, for example, on intersecting D-brane constructions in superstring theory in which gravity is essentially decoupled (see, for instance, [87] for, a recent attempt in this direction). Moreover in [58] we made the assumption that all the MSSM superpotential terms were perturbatively realized, that is, allowed by the extra abelian U(1) symmetries. In the following we assume instead that the μ-term is perturbatively forbidden. The origin of this term is rather non-perturbative and can be associated to an exotic instanton contribution which naturally arises from euclidean D-brane in the framework of a type IIA intersecting brane model (see [88] and references therein).The paper is organized as follows: in Section2 we describe the basic setup of the model and we discuss the perturbative and non-perturbative origin of the superpotential terms. We argue how the latter can naturally come from an intersecting D-brane model considered as the ultraviolet (UV) completion of our model. In Section 3 we describe the (global) supersymmetry breaking mechanism that gives mass to all the soft terms. In Section 5 we compute the gauge vector boson masses while in Section 4 we study the scalar potential of the theory in the neutral sector. In Section 6 we describe the neutralino sector while in Section 7 we describe the sfermion mass matrices. In Section 8 we study the phenomenology of our model and the bounds that can be put by some experimental constraints. Finally in Section 9 we draw our conclusions.


## 2. Model Setup

The model is an extension of the MSSM with two extra abelian gauge groups,U(1)A and U(1)B. The first one is anomalous while the second one is anomaly-free. This assumption is quite generic since in models with several anomalous U(1) symmetries there exists a unique linear combination which is anomalous while the other combinations are anomaly-free. The charge assignment for the chiral superfields is shown in Table 1. The vector and matter chiral multiplets undergo the usual gauge transformations
(2.1)V→V+i(Λ-Λ†),Φ→e-iqΛΦ.
The anomaly cancellation of the U(1)A gauge group is achieved by the four-dimensional analogue of the higher-dimensional GS mechanism which involves the Stückelberg superfield S=s+2θψS+θ2FS transforming as a shift
(2.2)S→S-2iMVAΛ,
where MVA is a mass parameter related to the anomalous U(1)A gauge boson mass. It turns out that not all the anomalies can be cancelled in this way. In particular the so-called mixed anomalies between anomalous and non-anomalous U(1)’s require the presence of trilinear GCS counterterms. For further details about the anomaly cancellation mechanism, see Appendix A (see also, for instance, [58] and [80–82]). The effective superpotential of our model at the scale E=MVA is given by
(2.3)W=WMSSM+λe-kSHuHd+mΦ+Φ-,
where WMSSM is given by
(2.4)WMSSM=yuijQiUjcHu-ydijQiDjcHd-yeijLiEjcHd,
which is the usual MSSM superpotential without the μ-term which is forbidden for a generic choice of the charges qHu and qHd. The second term in (2.3) is the only gauge invariant coupling allowed between the Stückelberg superfield and the two Higgs fields. This is the only allowed coupling with matter fields for a field transforming as (2.2). We will argue later about how nonperturbative effects can generate such a term. The last term in (2.3) is a mass term for Φ± which are charged under both U(1)A and U(1)B. These fields have been considered as supersymmetry breaking mediators in the context of anomalous models by Dvali and Pomarol [89]. They play a key role in generating gaugino masses. In the effective lagrangian, beside the usual kinetic terms (they are charged under both U(1)A and U(1)B), the two U(1)B fields Φ± couple to the gauge field strength Waα through the dimension six effective operator
(2.5)ℒg=caΦ+Φ-Λ2WaαWαa,
where a=A,B,Y,2,3, Λ is the cut-off scale of the theory while ca are constants that have to be computed in the UV completion of the theory.Table 1
Charge assignment.
S
U
(
3
)
c
S
U
(
2
)
L
U
(
1
)
Y
U
(
1
)
A
U
(
1
)
B
Q
i
3
2
1
/
6
q
Q
0
U
i
c
3
¯
1
-
2
/
3
q
U
c
0
D
i
c
3
¯
1
1
/
3
q
D
c
0
L
i
1
2
-
1
/
2
q
L
0
E
i
c
1
1
1
q
E
c
0
H
u
1
2
1
/
2
q
H
u
0
H
d
1
2
-
1
/
2
q
H
d
0
Φ
+
1
1
0
1
1
Φ
-
1
1
0
-
1
-
1The non perturbative term in (2.3) is expected to be generated in the effective action of intersecting D-brane models which can be considered as the UV completion of our model. This is the leading order term when the coupling HuHd is not allowed by gauge invariance. In string theory there are many axions related to the GS mechanism of anomaly cancellation which are charged under some Ramond-Ramond (RR) form. For example, in type IIA orientifold model with D6-branes, axion fields are associated to the C3 RR-form (see for a recent review [90]). Instantons charged under this RR-form, such as Euclidean E2-branes wrapping some γ3 3-cycle in the Calabi-Yau (CY) compactification manifold, give a contribution to the holomorphic couplings in the N=1 superpotential. Our analysis does not rely on any concrete intersecting brane model but rather on the generic appearance of such instanton induced terms. The exponential suppression factor of the classical instanton action is
(2.6)e-VolE2/gs,
where VolE2 is the volume of the 3-cycle in the CY wrapped by a E2-brane measured in string units while gs is the string coupling. Such exponential factor is independent of the d=4 gauge coupling, and thus this instanton is usually termed as stringy or exotic instanton (see [88, 91] and references therein). Moreover the instanton contribution can be sizable even in the case gs=1 if VolE2=1 measured in string units.In type IIA orientifold models with intersecting branes, the complexified moduli, whose imaginary part are the generalized axion fields (depending on the cycleγ3i), can be written as
(2.7)Ui=e-φ∫γ3i‍Ω3+i∫γ3i‍C3,
where φ is the dilaton, Ω3 is the CY volume 3-form (which is a complex form), and C3 is the RR-form. The integral of this form is dual to the axion whose shift symmetry is gauged in the GS mechanism. The generic contribution of an E2 instanton is formally given by
(2.8)W~∏i=1n‍Φai,bie-SE2,
where Φai,bi are chiral superfields localized at the intersection of two D6-branes described by open strings while SE2 denotes the instanton classical action:
(2.9)e-SE2=exp[-2πls3(1gs∫γ‍Re(Ω3)-i∫γ‍C3)].
This result can be immediately extended to the supersymmetric case which involves the complete Stückelberg multiplet. The appearance of the exponential suppression factor is dictated by the fact that the superpotential is a holomorphic quantity. Thus the only allowed functional dependence on the string coupling gs=e〈φ〉 and the axionic field is an exponential. Any other dependence can be excluded due to the shift transformation (2.2).


## 3. Supersymmetry Breaking

TheD-term contribution of the U(1)A vector multiplet VA relevant to supersymmetry breaking is given, in the limit of vanishing kinetic mixing δYA,δAB=0, by the following lagrangian:
(3.1)ℒ=12DADA+∑i‍gAqiϕi†DAϕi+ξDA,
where the sum is extended to all the scalars charged under the U(1)A. There is no D-term contribution related to the U(1)B except that of ϕ± since all the MSSM chiral fields are uncharged under U(1)B (see Table 1). The last term in (3.1) is a tree-level field-dependent Fayet-Iliopoulos (FI) term which comes from the supersymmetrized Stückelberg lagrangian
(3.2)ℒaxion=14(S+S†+2MVAVA)2|θ2θ-2+⋯=MVA(S+S†)VA|θ2θ-2+⋯=MVAαDA+⋯,
where in the last line α denotes the real part of the lowest component of the Stückelberg chiral multiplet s=α+iφ. The fields α and φ are called the saxion and the axion, respectively (with a slight abuse of notation with respect to the previous section where we denoted the dilaton with φ). We assume that the real part α gets an expectation value. This gives a contribution to the gauge coupling constants which can be absorbed in the following redefinition:
(3.3)116ga2τa=116g~a2τa-12baa〈α〉,
where the gauge factors τa take the values 1,1,1,1/2,1/2, and the baa constants are given in (A.2). The tree-level FI term is then given by
(3.4)ξ=MVA〈α〉.
Moreover in the following we assume that 1-loop FI terms are absent (see the discussion in [92]). The FI term induces a mass term for the scalars. This can be seen by solving the equations of motion for DA:
(3.5)DA+∑i‍gAqiϕi†ϕi+ξ=0,
where the index i runs over all chiral superfields. The D-term contribution to the scalar potential is given by
(3.6)V(ϕi,ϕi†)=12(ξ+gA∑i‍qi|ϕi|2)2.
The quadratic part gives the scalar mass term
(3.7)∑i‍ξgAqi|ϕi|2=∑i‍mi2|ϕi|2,
where we have defined
(3.8)mi2=ξgAqi=〈α〉gAMVAqi=qimξ2,
with
(3.9)mξ2=〈α〉gAMVA=gAξ.
The typical scale for the mass mξ is of the order of few hundreds of GeV if MVA~〈α〉~1 TeV and gA~0.1. It is interesting to note that in this scenario a low sub-TeV supersymmetry breaking scale mξ is due to the Stückelberg mechanism which gives mass to VA. This is the most important difference with the scenario proposed in [89], where the scale mξ is dynamically generated by some dynamics in a strong coupling regime.Mass terms for the gauginos, that is,λaλa, are generated by the dimension six effective operator (2.5) in the broken phase where ϕ± get vacuum expectation value (vev). The contribution coming from this mechanism is
(3.10)Ma=ca〈F+ϕ-〉+〈F-ϕ+〉Λ2=cam(v+2+v-2)2Λ2,
where v±/2=〈ϕ±〉 and where in the right-hand side we have used the F-term equations of motion for F±(3.11)F±=-∂W*∂ϕ±*=-mϕ∓*,
having assumed m real without any loss in generality. We assume ca=c for each a. This is an assumption of universality as a boundary condition at the cut-off scale Λ which does not affect in a crucial way our analysis. In Section 4 we study the scalar potential of our model and we derive the conditions for having a vev for ϕ± different from zero. Since we are breaking supersymmetry in the global limit in which the Planck mass MP→∞, the F-term induced contribution to the scalar masses
(3.12)mi2~〈F±〉MP2
vanishes leaving (3.8) as the leading contribution.The requirement of gauge invariance of the superpotential implies the following constraints on theU(1)A charges
(3.13)qUc=-qQ-qHu,qDc=-qQ-qHd,qEc=-qL-qHd,(3.14)k=qHu+qHd2MVA.
As we said at the beginning of this section, we assume that the net kinetic mixing between U(1)Y and U(1)A vanishes (we postpone the discussion about the kinetic mixing between U(1)A and U(1)B to the next section). There are two contributions for the U(1)Y-U(1)A kinetic mixing: the 1-loop mixing δYA and bYA coming from the GS coupling SWYWA (see (A.1)). The following conditions imply a bound on the charges
(3.15)δYA=0⇒∑f‍qfYf=0,bYA=0⇒∑f‍qf2Yf=0,
where the sum is extended over all the chiral fermions in the theory. Constraints (3.15), can be solved in terms of qQ and qL. By using conditions (3.13), we get
(3.16)qL=14(3qHu-4qHd),qQ=-112(5qHu-2qHd).
The positive squared mass condition for the sfermions
(3.17)mf~2=gAqfMVA〈α〉>0,
implies qf>0 for all the sfermions having assumed without loss of generality 〈α〉>0. Using the constraints (3.13) and (3.16), we get the allowed parameter space
(3.18)qHu<0,52qHu<qHd<34qHu.


## 4. Scalar Potential

The key ingredient in our model is the instanton-induced term in (2.3) which couples the Stückelberg field to the Higgs fields. The θ2 component of this superpotential term gives the following contribution to the Lagrangian:
(4.1)Winst|θ2=λe-kSHuHd|θ2=λe-kshuFd+λe-ksFuhd-λke-ksFShuhd+2λe-ksk(huψSh~d+hdψSh~u)-λe-kshuhdψSψS,
where Fu,d are the F-terms of Hu,d. Solving the F-terms equations for Hu and Hd, we get the following contributions for the instanton-induced term in the scalar potential:
(4.2)Vinst=2λ2e-2kαhu†hu+2λ2e-2kαhd†hd+λke-kα(e-ikφFShuhd+h.c.).
In the following we assume that α gets a vev different from zero and that the mass of this field is much higher than Λ so that its dynamics is not described by the low-energy effective action. From the point of view of the UV completion (e.g., a type IIA intersecting brane model), this amounts to saying that the closed string modulus related to α is stabilized. Moreover we made the assumption that the same dynamics that stabilizes α also fixes FS. By supersymmetry the saxion field α, being part of the Stu¨ckelberg multiplet, has a tree-level mass MVA. Thus if we want to consider a frozen dynamics for α at the TeV scale we have to assume a mass parameter for the anomalous U(1)A just slightly above the TeV scale; that is, MVA>1 TeV. In this way the effective instanton-induced potential at a scale E≃1 TeV is thus given by
(4.3)Vinst=2λ2e-2k〈α〉hu†hu+2λ2e-2k〈α〉hd†hd+λke-k〈α〉(〈FS〉e-ikφhuhd+h.c.).
The first two terms are μ-terms while the third one is a b-term. The complete effective scalar potential is given by
(4.4)V=(|μ|2+mhu2)(|hu0|2+|hu+|2)+(|μ|2+mhd2)(|hd0|2+|hd-|2)+(|m|2+mϕ+2)|ϕ+|2+(|m|2+mϕ-2)|ϕ-|2+[be-ikφ(hu+hd--hu0hd0)+h.c.]+18(g22+gY2)(|hu0|2+|hu+|2-|hd0|2-|hd-|2)2+12g22|hu+hd0*+hu0hd-*|2+12gA2[qHu(|hu0|2+|hu+|2)+qHd(|hd0|2+|hd-|2)+|ϕ+|2-|ϕ-|2]2+12gB2[|ϕ+|2-|ϕ-|2]2,
where
(4.5)μ=2λe-k〈α〉,(4.6)b=λke-k〈α〉〈FS〉.
These relations give a solution of the well-known μ-problem since both terms have a common origin (see the analysis in Section 8.2). The soft squared masses are generated by the FI U(1)A term:
(4.7)mhu2=qHumξ2,(4.8)mhd2=qHdmξ2,(4.9)mϕ+2=mξ2,(4.10)mϕ-2=-mξ2,
with mξ2 given by (3.8). The scalar potential depends on the following new parameters: 〈α〉, 〈FS〉, λ, m, gA,B, qHu,d, MVA.In order to have a vacuum preserving the electromagnetism, the charged field vevs must vanish. Thus we are left with the problem of finding a minimum for the neutral scalar potential(4.11)V0=(|μ|2+mhu2)|hu0|2+(|μ|2+mhd2)|hd0|2-(be-ikφhu0hd0+h.c.)+(|m|2+mϕ+2)|ϕ+|2+(|m|2+mϕ-2)|ϕ-|2+18(g22+gY2)(|hu0|2-|hd0|2)2+12gA2(qHu|hu0|2+qHd|hd0|2+|ϕ+|2-|ϕ-|2)2+12gB2[|ϕ+|2-|ϕ-|2]2.
Since there are no D-flat directions along which the quartic part vanishes, the potential is always bounded from below. To find the minimum we solve ∂V0/∂zi=0 where the scalar field zi runs over {φ,hu0,hd0,ϕ+,ϕ-}. The conditions for having a nontrivial minimum boil down to the same condition of the MSSM
(4.12)b2>(|μ|2+mhu2)(|μ|2+mhd2).
Moreover in order to generate a mass term for the gauginos (see (3.10)), the condition v-≠0 must hold since v+=0 due to the positive sign of the coefficient of the ϕ+ quadratic term in  (4.9). This implies the following condition for the coefficient of the ϕ- quadratic term:
(4.13)|m|2+mϕ-2<0,
The minimum is attained at φ=ϕ+=0. Actually since the potential for the axion φ is periodic, the minimum condition holds for φ=2nπ/k with n∈ℤ. All these minima are physically equivalent, and thus we arbitrarily choose n=0. The remaining three conditions imply the following constraints on the parameters:
(4.14)mhd2+μ2-btβ+18(gY2+g22)v2c2β+12gA2qHd[v2(qHdcβ2+qHusβ2)-v-2]=0,mhu2+μ2-btβ-1-18(gY2+g22)v2c2β+12gA2qHu[v2(qHdcβ2+qHusβ2)-v-2]=0,(gA2+gB2)v-2-gA2v2(qHdcβ2+qHusβ2)+2(|m|2+mϕ-2)=0,
where we have defined in order to keep a compact notation
(4.15)cβ=cosβ,sβ=sinβ,tβ=tanβ,c2β=cos(2β),s2β=sin(2β),
and as usual as tanβ=vu/vd.In the previous discussion we treated the scalar potential in an exact way. In the following we want to introduce some useful approximation in order to compute the mass eigenstates. Let us go back to the minima equations (4.14). Supposing
(4.16)v≪v-,
we can neglect all the gAv terms. With this approximation, the minima equations read
(4.17)m~hd2+μ2-btβ+18(gY2+g22)v2c2β=0,(4.18)m~hu2+μ2-btβ-1-18(gY2+g22)v2c2β=0,(4.19)(gA2+gB2)v-2+2(|m|2+mϕ-2)=0,
where we have defined
(4.20)m~hd2=mhd2-12gA2qHdv-2,m~hu2=mhu2-12gA2qHuv-2.
Equations (4.17) and (4.18) have the same functional form as in the MSSM case. Moreover v- does not depend on any parameter of the visible sector. Within this approximation the dynamics of the fields ϕ± is decoupled from that of the Higgs sector, and thus the Higgs potential can be studied by fixing ϕ± at their vevs. We get
(4.21)V≃(|μ|2+mhu2)(|hu0|2+|hu+|2)+(|μ|2+mhd2)(|hd0|2+|hd-|2)+[be-ikφ(hu+hd--hu0hd0)+h.c.]+18(g22+gY2)(|hu0|2+|hu+|2-|hd0|2-|hd-|2)2+12g22|hu+hd0*+hu0hd-*|2+12gA2[qHu(|hu0|2+|hu+|2)+qHd(|hd0|2+|hd-|2)-12v-2]2,
neglecting further constant terms in v-. Close to the minima, the relevant term in the last line of (4.21) is the double product of the Higgs part with the v-2 term. Hence by using (4.16) we finally get
(4.22)Vh,φ≃(|μ|2+m~hu2)(|hu0|2+|hu+|2)+(|μ|2+m~hd2)(|hd0|2+|hd-|2)+[be-ikφ(hu+hd--hu0hd0)+h.c.]+18(g22+gY2)(|hu0|2+|hu+|2-|hd0|2-|hd-|2)2+12g22|hu+hd0*+hu0hd-*|2.
This potential has the same form (except for the contribution of the exponential term in φ) of the MSSM potential, and the corresponding minima equations are exactly given in (4.17) and (4.18). Thus all the well-known MSSM results apply here [93].In particular one of the constraints istβ≳1.2 [93] which implies (The presence of the extra field φ does not affect this result since the minima conditions are the same as the MSSM.) m~hu2<m~hd2. By using (4.20) we get
(4.23)gAqHu(MVA〈α〉-12gAv-2)<gAqHd(MVA〈α〉-12gAv-2).
By assuming MVA>1 TeV, v- in the TeV range, gA~O(0.1), the term between brackets is positive and we get the following constraint:
(4.24)qHu<qHd
for the U(1)A Higgs charges.


### 4.1. Higgs Mass Matrices

We discuss the mass eigenvalues starting from the exact form of the scalar potential (4.4), switching to the approximated expression (4.22) when needed. In the neutral sector the singlet scalar ϕ+ does not mix with any other scalar, so it is a mass eigenstate with square mass
(4.25)Mϕ+2=2|m|2.
The same holds for the imaginary part of ϕ- which becomes the longitudinal mode of the gauge vector Z2. The mass matrix for the real scalar fields {φ,Im(hu0),Im(hd0)} is given by
(4.26)ℳS(Im)=(btβ⋯⋯bbtβ-1⋯-bkvsβ-bkvcβbk2v2cβsβ).
The determinant of this matrix is zero. Two eigenvalues are zero which correspond to the Goldstone modes of Z0 and Z1. The physical massive state is an axi-Higgs state with mass given by
(4.27)MA02=2bs2β[1-116(qHu+qHd)2v2MVA2s2β2],
where we used the relation (3.14). The mass matrix for the real scalar fields {Re(hu0), Re(hd0), ϕR-≡Re(ϕ-)} reads as
(4.28)ℳS(Re)=((14gEW2+gA2qHd2)v2cβ2+btβ⋯⋯-b-(14gEW2-gA2qHdqHu)v2cβsβ(14gEW2+gA2qHu2)v2sβ2+btβ-1⋯-gA2qHdvv-cβ-gA2qHuvv-sβ(gA2+gB2)v-2),
where gEW2=(gY2+g22). The matrix can be diagonalized exactly, but the results are cumbersome and difficult to read. It is much more convenient starting from the approximated potential (4.22) neglecting the mixing between Higgses and ϕ-. In this case we can apply the MSSM equations and get the following mass eigenvalues:
(4.29)Mh0,H02≃12(2bs2β∓(2bs2β-14(gY2+g22)v2)2+2b(gY2+g22)v2s2β),MϕR-2≃(gA2+gB2)v-2.The charged sector is unchanged with respect to the MSSM, so(4.30)MH±2=2bs2β+MW2.As in the standard MSSM case, the mass of the lightest HiggsMh0 has a theoretical bound. It is a well-known problem in the MSSM that the upper bound [94, 95] is not compatible with the LEP bound [96]. In our case the bound is increased due to the presence of DA-term corrections
(4.31)Mh02<14(gY2+g22)v2c2β2+14gA2v2[qHd+qHu+(qHd-qHu)c2β]2,
where the first term is the MSSM bound. In principle, for arbitrary high values of gAqHd, gAqHu, we get an increasing upper bound. However, as in the standard MSSM case, Mh02 undergoes to relatively drastic quantum corrections [93]. Hence in Section 8 we consider tree-level masses for all the particles except for h0 for which we use the 1-loop corrected expression (see (8.10)).


## 4.1. Higgs Mass Matrices

We discuss the mass eigenvalues starting from the exact form of the scalar potential (4.4), switching to the approximated expression (4.22) when needed. In the neutral sector the singlet scalar ϕ+ does not mix with any other scalar, so it is a mass eigenstate with square mass
(4.25)Mϕ+2=2|m|2.
The same holds for the imaginary part of ϕ- which becomes the longitudinal mode of the gauge vector Z2. The mass matrix for the real scalar fields {φ,Im(hu0),Im(hd0)} is given by
(4.26)ℳS(Im)=(btβ⋯⋯bbtβ-1⋯-bkvsβ-bkvcβbk2v2cβsβ).
The determinant of this matrix is zero. Two eigenvalues are zero which correspond to the Goldstone modes of Z0 and Z1. The physical massive state is an axi-Higgs state with mass given by
(4.27)MA02=2bs2β[1-116(qHu+qHd)2v2MVA2s2β2],
where we used the relation (3.14). The mass matrix for the real scalar fields {Re(hu0), Re(hd0), ϕR-≡Re(ϕ-)} reads as
(4.28)ℳS(Re)=((14gEW2+gA2qHd2)v2cβ2+btβ⋯⋯-b-(14gEW2-gA2qHdqHu)v2cβsβ(14gEW2+gA2qHu2)v2sβ2+btβ-1⋯-gA2qHdvv-cβ-gA2qHuvv-sβ(gA2+gB2)v-2),
where gEW2=(gY2+g22). The matrix can be diagonalized exactly, but the results are cumbersome and difficult to read. It is much more convenient starting from the approximated potential (4.22) neglecting the mixing between Higgses and ϕ-. In this case we can apply the MSSM equations and get the following mass eigenvalues:
(4.29)Mh0,H02≃12(2bs2β∓(2bs2β-14(gY2+g22)v2)2+2b(gY2+g22)v2s2β),MϕR-2≃(gA2+gB2)v-2.The charged sector is unchanged with respect to the MSSM, so(4.30)MH±2=2bs2β+MW2.As in the standard MSSM case, the mass of the lightest HiggsMh0 has a theoretical bound. It is a well-known problem in the MSSM that the upper bound [94, 95] is not compatible with the LEP bound [96]. In our case the bound is increased due to the presence of DA-term corrections
(4.31)Mh02<14(gY2+g22)v2c2β2+14gA2v2[qHd+qHu+(qHd-qHu)c2β]2,
where the first term is the MSSM bound. In principle, for arbitrary high values of gAqHd, gAqHu, we get an increasing upper bound. However, as in the standard MSSM case, Mh02 undergoes to relatively drastic quantum corrections [93]. Hence in Section 8 we consider tree-level masses for all the particles except for h0 for which we use the 1-loop corrected expression (see (8.10)).


## 5. Vector Mass Matrix

We now discuss the vector mass matrix. All the neutral scalars could in principle take a vev different from zero; hence we assume(5.1)〈ϕ±〉=v±2,〈hu,d0〉=vu,d2.
The neutral vector square mass matrix in the base (VB,VA,VY,V23) is
(5.2)ℳV=(gB2vϕ2⋯⋯⋯gAgBvϕ2gA2[(cβ2qHd2+sβ2qHu2)v2+vϕ2]+MVA2⋯⋯012gAgYqH(β)v214gY2v2⋯0-12gAg2qH(β)v2-14gYg2v214g22v2),
where
(5.3)vϕ2=v+2+v-2,v2=vu2+vd2,qH(β)=(sβ2qHu-cβ2qHd).
By taking MVA>1 TeV (see Section 4), VA can be considered as decoupled from the low-energy gauge sector (namely, E≲1 TeV), and we can ignore with very good approximation any mixing term (The kinetic mixing between U(1)A and U(1)B deserves some comment, in particular if we relax the MVA>1 TeV assumption. Actually the presence of this mixing turns out to be irrelevant for the phenomenology of the visible sector. Anyway one has to take into account that for Tr(qAqB)≠0 such a mixing arises at the 1-loop level. In such a case it can be assumed that the two U(1)'s are in the kinetic diagonalized basis with Tr(qAqB)=0 thanks to some additional heavy chiral multiplet charged under both U(1)A and U(1)B. These multiplets generate a counterterm in the effective theory that cancels against δAB making the net kinetic mixing term equal to zero. This mechanism is analogous to the anomaly cancellation one where the GS mechanism can be generated by an anomaly-free theory with some heavy chiral fermion integrated out of the mass spectrum [80–82].) involving VA. From now on we will apply this approximation.SinceVB is a hidden gauge boson, it is decoupled from the SM sector. The charged vector sector is unchanged with respect to the MSSM, so
(5.4)Wμ±=V21μ∓iV22μ2,MW2=14g22v2.


## 6. Neutralinos

In comparison with the standard MSSM, we now have five new neutral fermionic fields:ψS, λA, λB, ϕ~±. However under the assumption MVA>1 TeV, ψS and λA are not in the low-energy sector because of the MVA mass term (we stress that the ψS-λA sector presents a different parameters’ choice with respect to [97–99],  where we realized a scenario in which the mixing between ψS and λA was suppressed). Thus we have
(6.1)ℒneutralinomass=-12(ψ0)TℳN~ψ0+h.c.,
where
(6.2)(ψ0)T=(λB,ϕ~-,ϕ~+,λY,λ20,h~d0,h~u0).
In this basis the neutralino mass matrix ℳN~ is written as
(6.3)ℳN~=(MB⋯⋯⋯⋯⋯⋯-gBv-0⋯⋯⋯⋯⋯0-m0⋯⋯⋯⋯000M1⋯⋯⋯0000M2⋯⋯000-g1vd2g2vd20⋯000g1vu2-g2vu2-μ0),
where μ is given in (4.5). We recall that gaugino masses arise from the Dvali-Pomarol term (2.5).ℳ
N
~ factorizes in a 4×4 MSSM block in the lower right corner, and in a 3×3 new sector block in the upper left corner. The new sector block is given by the λB and ϕ~± contributions. This last block has a MSSM-like structure that can be easily understood just considering the superpotential (2.3) and the gaugino masses (3.10) and by recalling that ϕ- gets a vev v- different from zero, while v+=0.Finally there are also corrections coming from the anomalous axino couplings:F-term couplings of the type baa〈FS〉λaλa, D-term couplings of the type baaψSλa〈Da〉, and corrections coming from the superpotential term e-kSHuHd+h.c. However such corrections are always subdominant, and thus we neglect them with very good approximation.We assume the lightest supersymmetric particle (LSP) in our model comes from the neutralino sector. In Section8 we show the parameter regions in which this holds true. In order to ensure that the neutralino is the LSP, we keep fixed the gravitino mass m3/2~O (TeV) in the limit MP→∞.


## 7. Sfermion Masses

The sfermion masses receive several contributions. As we have seen in Section3, the leading contribution comes from the induced soft masses (3.8). But there are further contributions. We have MSSM-like contributions: F-term corrections proportional to the Yukawa couplings and DY and D2 term correction from the Higgs sector. Moreover there are DA term corrections from the Higgs and ϕ- sector. As an aside, the appearance of such terms in the low-energy action, given our assumption MVA>1 TeV, can be understood in terms of quantum corrections to, Kähler potential [100]. Considering the first two families, we neglect the corresponding Yukawa couplings (the so-called third family approximation). In this approximation the mass eigenvalues are given by
(7.1)mu~L2≃mc~L2=mQ~2+(13gY2-g22)Δv28+qQm~DA2,mu~R2≃mc~R2=mU~c2-gY2Δv26+qUcm~DA2,md~L2≃ms~L2=mQ~2+(13gY2+g22)Δv28+qQm~DA2,md~R2≃ms~R2=mD~c2+gY2Δv212+qDcm~DA2,mν~e2=mν~μ2=mL~2-(gY2+g22)Δv28+qLm~DA2,me~L2≃mμ~L2=mL~2-(gY2-g22)Δv28+qLm~DA2,me~R2≃mμ~R2=mE~c2+gY2Δv24+qEcm~DA2.
The first terms on the right-hand side mQ~,U~c,D~c,L~,E~c2 are the corresponding soft masses (3.8), the second terms are the DY,2 contributions with Δv2=vu2-vd2=-v2c2β, while the last terms are the DA corrections given by
(7.2)m~DA2=12(qHuvu2+qHdvd2-v-2),
There is an approximated degeneracy between the sfermions with the same charges.The mass matrix for the third family sfermions is parametrized as(7.3)ℳf~2=(MLLf~2MLRf~2MLRf~2MRRf~2),
where the off-diagonal terms are generated by F-term corrections proportional to the Yukawa couplings. The stop mass matrix elements are
(7.4)MLLt~2=mt2+mQ~2+(13gY2-g22)Δv28+qQm~DA2,MRRt~2=mt2+mU~c2-gY2Δv26+qUcm~DA2,MLRt~2=-μmttβ-1.
The sbottom mass matrix elements are
(7.5)MLLb~2=mb2+mQ~2+(13gY2+g22)Δv28+qQm~DA2,MRRb~2=mb2+mD~c2+gY2Δv212+qDcm~DA2,MLRb~2=-μmbtβ.
The stau mass matrix elements are
(7.6)MLLτ~2=mτ2+mL~2-(gY2-g22)Δv28+qLm~DA2,MRRτ~2=mτ2+mE~c2+gY2Δv24+qEcm~DA2,MLRτ~2=-μmτtβ.
The tau sneutrino mass is
(7.7)mν~τ2=mL~2-(gY2+g22)Δv28+qLm~DA2,
where mt, mb, and mτ are the masses of the corresponding standard fermions (i.e., further F-term contributions proportional to the Yukawa couplings). The structure of the diagonal terms of (7.3) is the same as in (7.1): soft masses, MSSM D-term contribution, and DA term correction. Furthermore we stress that there is a mass degeneracy between the three sneutrinos ν~e,μ,τ since the soft masses (3.8) are flavor blind.


## 8. Phenomenology

In the following we derive the phenomenological consequences of our scenario. Following our assumption of having a mass parameter for the anomalousU(1)A just slightly above the TeV scale, we fix MVA=10 TeV. The mass scale in the gaugino sector Λ is set to be O(MVA).


### 8.1. Charge Bounds

The model parameter space can in principle be constrained by precision EW measurements [101]. However, since MVA=10 TeV, every value of gAqHu and gAqHd is allowed by EW precision data if |gAqHu|,|gAqHd|≲0.1. So the only relevant constraints are (3.18) and (4.24) that are plotted, respectively, with a red and a blue region, in Figure 1 in the plane (gAqHu, gAqHd).Figure 1
Higgs couplings bounds. The yellow spot represents our charge choice.


### 8.2. Free Parameters

Here we discuss which parameters remain free in our model after all the constraints discussed in the previous sections are imposed. Our choice for the HiggsU(1)A charges corresponds to the yellow spot in Figure 1:
(8.1)gA=0.1,MVA=10TeV,qHd=-(13),qHu=-(25).
In order to fix the remaining parameters (〈α〉, 〈FS〉, λ, m, gB) we assume v≃246 GeV, and then we choose some benchmark value for gB and v- in the U(1)B sector (we recall that v+=0  (see Section 4)):
(8.2)(A)gB=0.4,v-=5TeV,(B)gB=0.1,v-=4TeV.
The next step is to solve the minima conditions (4.14) determining 〈FS〉, λ, m as function of 〈α〉. In the limit in which v2≪MVA〈α〉,v-2, we get
(8.3)λ2≃18e2〈α〉gA(qHd+qHu)/MVA[gA(gAv-2-2〈α〉MVA)(sec(2β)(qHd-qHu)+qHd+qHu)],〈FS〉≃-e〈α〉gA(qHd+qHu)/MVAMVAtan(2β)4(qHd+qHu)λ(qHd-qHu)(2〈α〉MVA-gAv-2),|m|2≃gA〈α〉MVA-12[(gA2+gB2)v-2].
In Appendix B  we report the exact formulae. Thus the only remaining free parameters are tβ and 〈α〉, and we perform the following analysis of the mass spectrum as a function of tβ and 〈α〉. A lower bound on 〈α〉 as a function of v- can be obtained, given the approximation (4.16); from (4.19),
(8.4)〈α〉≃|m|2+1/2(gA2+gB2)v-2gAMVA,
where we used the relation
(8.5)mϕ-2=-mξ2=-〈α〉gAMVA.
Thus the lower bound on 〈α〉 is obtained simply by setting |m|=0:
(8.6)〈α〉m≃1/2(gA2+gB2)v-2gAMVA.
The condition 〈α〉>〈α〉m must hold since otherwise we would have a massless scalar field in the spectrum (see (4.25)). Another lower bound, 〈α〉b, can be obtained from the condition (4.12), by solving the minima conditions (4.14) and by substituting the corresponding 〈FS〉, λ and m values (8.3). The resulting lower bound can be expressed as
(8.7)〈α〉>max[〈α〉m,〈α〉b].
No upper bound can be imposed; hence we decide to perform our analysis by considering 〈α〉≲100 TeV.The parametersλ and 〈FS〉 are of a particular phenomenological importance since they appear in the μ and b terms (see (4.5) and (4.6)). In the case A, μ is in the range (900,6000)  GeV and b is in the range (50,1200)  GeV while in the case B, μ is in the range (500,6000)  GeV and b is in the range (25,1200)  GeV. These values are in the right range to solve the μ-problem.


### 8.3. Mass Spectrum

(A)
With such choice the gauge vector sector is completely fixed up to atβ dependence. Anyway even such a dependence can be safely ignored with a very good approximation in the new gauge sector since the mixing is strongly suppressed. So for each tβ value, we have
(8.8)MZ1≃10TeV,MZ2≃2TeV,
where with Z1 we denote the VA-like vector.
(B)
As in the previous case, we just give theZ1,2 masses
(8.9)MZ1≃10TeV,MZ2≃400GeV,
where as in the previous case Z1 is VA-like.We will not give the exact values of theZ0 mass. It is enough for our purposes to know that they are compatible with the bounds of  Section 8.1. Both case A and B are compatible with CDF bounds about Z' direct production [102].Recent LHC data have restricted the most probable range for the Higgs particle mass to be[115.5,131] GeV (ATLAS) [103] and [114.5,127] (CMS) [104]. Moreover, there are hints observed by both CMS and ATLAS of an excess of events that might correspond to decays of a Higgs particle with a mass in a range close to 125  GeV. So, in Figures 2 and 3 we give region plots showing the allowed values of 〈α〉 and tβ for case A (B) and Λ/c=(5)10TeV. The red region is the one in which Mh02|1-loop∈[124,126]GeV where the h0 mass is computed considering 1-loop corrections. Since it turns out that the top squarks have small mixing angle and considering the limit MA0≫MZ0, we have [93]
(8.10)Mh02|1-loop≃Mh02|tree+34π2sβ2yt2mt2ln(mt~1mt~2mt2)≃Mh02|tree+32π2mt4v2ln(mt~1mt~2mt2),
where Mh0|tree is the tree-level h0 mass and we used mt=ytvu/2=ytvsβ/2. There is an approximated inverse correlation between 〈α〉 and tβ in the h0 mass allowed region because the 1-loop correction in (8.10) increases for increasing values of 〈α〉 or tβ. The h0 mass allowed region is almost the same for case A and B because of two reasons:  (i)
the mixing withϕ± is suppressed,(ii)
the parametersm~hu,m~hd,μ,b in the scalar potential (4.22) are ruled by the square mass parameters gA〈α〉MVA and (gAv-)2, and the first one turns out to be dominant.Allowed〈α〉 and tanβ values for case A (up) and case B (down), Λc=5 (left) and Λc=10 (right). The red region is the one in which Mh02|1-loop∈[124,126]GeV, the magenta region is the one in which Mh02|1-loop∈[114.5,131]GeV and the blue region satisfies all the mass bounds on the sparticles (from PDG) and requires a neutralino LSP. The yellow dots are our benchmark points.
(a)
(b)
(c)
(d)Allowed〈α〉 and tanβ values for case A (up) and case B (down), Λc=5 (left) and Λc=10 (right). The red region is the one in which Mh02|1-loop∈[124,126]GeV, the magenta region is the one in which Mh02|1-loop∈[114.5,131]GeV and the blue region satisfies all the mass bounds on the sparticles (from preliminary LHC data) and requires a neutralino LSP. The yellow dots are our benchmark points.
(a)
(b)
(c)
(d)The magenta region satisfies a milder constraint on the light Higgs boson:Mh02|1-loop∈[114.5,131]  GeV. In order to be more conservative, we imposed the joint constraints of ATLAS and CMS.The blue region satisfies all the mass bounds on the sparticles and requires a neutralino LSP. We considered two possibilities: one more optimistic (Figure2) using the PDG bounds [105, 106] and one more conservative (Figure 3) using recent LHC data [107, 108]. The combination of the gluino mass bound with a neutralino LSP is a strong constraint that reduces drastically the allowed parameter space. In some cases there is not even a blue region, which means that we cannot satisfy simultaneously all the mass bounds and have a neutralino LSP, so they are completely ruled out. When the gluino mass bound is from PDG, case A is allowed; otherwise it is completely ruled out and only case B for Λ/c=10TeV presents allowed regions. We notice that case A favors low 〈α〉 values, while case B favors high 〈α〉 values. For every allowed case we choose a benchmark point (yellow spots in Figures 2 and 3):(i)
case A,Λc=5, 〈α〉=3 TeV and tβ=50 so that Mh0|1-loop≃121.6 GeV;(ii)
case A,Λc=10, 〈α〉=5 TeV and tβ=10 so that Mh0|1-loop≃124.7 GeV;(iii)
case B,Λc=10, 〈α〉=4.5 TeV and tβ=50 so that Mh0|1-loop≃125.1 GeV;(iv)
case B,Λc=10, 〈α〉=50 TeV and tβ=2.5 so that Mh0|1-loop≃130.1 GeV,and we give the full mass spectrum in Figures 4 and 5.Mass spectrum, case A,Λc=5, 〈α〉=0.3 TeV and tβ=50 (a), Λc=10, 〈α〉=0.5 TeV and tβ=10 (b).
(a)
(b)Mass spectrum, case B,Λc=10, 〈α〉=0.45 TeV and tβ=50 (a), 〈α〉=8 TeV and tβ=2.5 (b).
(a)
(b)All the benchmark points share some common features.(i)
The LSP is the lightest neutralino of the new sector: in case A it is a combination ofϕ~± and λB while in case B is almost a pure λB.(ii)
An approximated mass degeneracy ofH0, A0, and H± holds, and their masses satisfy the bounds of [96, 109].(iii)
The lightest sleptons is a sneutrino, except fortβ=50 when it is τ~1(iv)
The lightest squark isu~L, except for tβ=50 when it is b~1(v)
The first and second family left-handed squarks/sleptons are likely to be lighter than their right-handed counterparts. This is at odds with the usual MSSM cases [93].(vi)
C~1(2) is close in mass with N~1(4)MSSM. C~2 and N~4MSSM are heavier than all sfermions.(vii)
The gluino is close in mass toC~1 and N~1MSSM which are gaugino-like. Moreover it is lighter than all the squarks except for point (i). So it turns out to be long lived, specially in case B where the approximated mass degeneracy involves also the LSP. Long-lived gluinos bind with SM quarks and gluons from the vacuum during the hadronisation process and produce R-hadrons. R-hadrons are among the most interesting searches at LHC. Anyway we will come back to this point with a more detailed study in a forthcoming paper.(viii)
There is an approximated mass degeneracy betweene~R and u~R because using the charge constraints (3.16) and (8.1) we get qEc=3 and qUc≃2.9.(ix)
mϕR-<mϕ+ except for point (i).Case B points deserve some more comments.ϕ
+ and N~2,3new are out of the plot of point (iv) because they are heavier than 6 TeV. Z2 is among the lightest not SM particle, so it can decay only into SM particles, because of energy and R-parity conservation. So Z2 is long lived, because SM particles are coupled to Z2 only through the suppressed VA,B mixing or through the Higgs scalars which present a tiny mixing with ϕ±.It is not an easy task to compare the resulting spectrum we get for our model with those related to the rich zoology of supersymmetry breaking scenarios. It is worth to stress anyway that the two representative spectrums showed in Figure5 which encode the key features of our scenarios listed above are not reproduced in any of the benchmark points showed in [110, 111].


## 8.1. Charge Bounds

The model parameter space can in principle be constrained by precision EW measurements [101]. However, since MVA=10 TeV, every value of gAqHu and gAqHd is allowed by EW precision data if |gAqHu|,|gAqHd|≲0.1. So the only relevant constraints are (3.18) and (4.24) that are plotted, respectively, with a red and a blue region, in Figure 1 in the plane (gAqHu, gAqHd).Figure 1
Higgs couplings bounds. The yellow spot represents our charge choice.


## 8.2. Free Parameters

Here we discuss which parameters remain free in our model after all the constraints discussed in the previous sections are imposed. Our choice for the HiggsU(1)A charges corresponds to the yellow spot in Figure 1:
(8.1)gA=0.1,MVA=10TeV,qHd=-(13),qHu=-(25).
In order to fix the remaining parameters (〈α〉, 〈FS〉, λ, m, gB) we assume v≃246 GeV, and then we choose some benchmark value for gB and v- in the U(1)B sector (we recall that v+=0  (see Section 4)):
(8.2)(A)gB=0.4,v-=5TeV,(B)gB=0.1,v-=4TeV.
The next step is to solve the minima conditions (4.14) determining 〈FS〉, λ, m as function of 〈α〉. In the limit in which v2≪MVA〈α〉,v-2, we get
(8.3)λ2≃18e2〈α〉gA(qHd+qHu)/MVA[gA(gAv-2-2〈α〉MVA)(sec(2β)(qHd-qHu)+qHd+qHu)],〈FS〉≃-e〈α〉gA(qHd+qHu)/MVAMVAtan(2β)4(qHd+qHu)λ(qHd-qHu)(2〈α〉MVA-gAv-2),|m|2≃gA〈α〉MVA-12[(gA2+gB2)v-2].
In Appendix B  we report the exact formulae. Thus the only remaining free parameters are tβ and 〈α〉, and we perform the following analysis of the mass spectrum as a function of tβ and 〈α〉. A lower bound on 〈α〉 as a function of v- can be obtained, given the approximation (4.16); from (4.19),
(8.4)〈α〉≃|m|2+1/2(gA2+gB2)v-2gAMVA,
where we used the relation
(8.5)mϕ-2=-mξ2=-〈α〉gAMVA.
Thus the lower bound on 〈α〉 is obtained simply by setting |m|=0:
(8.6)〈α〉m≃1/2(gA2+gB2)v-2gAMVA.
The condition 〈α〉>〈α〉m must hold since otherwise we would have a massless scalar field in the spectrum (see (4.25)). Another lower bound, 〈α〉b, can be obtained from the condition (4.12), by solving the minima conditions (4.14) and by substituting the corresponding 〈FS〉, λ and m values (8.3). The resulting lower bound can be expressed as
(8.7)〈α〉>max[〈α〉m,〈α〉b].
No upper bound can be imposed; hence we decide to perform our analysis by considering 〈α〉≲100 TeV.The parametersλ and 〈FS〉 are of a particular phenomenological importance since they appear in the μ and b terms (see (4.5) and (4.6)). In the case A, μ is in the range (900,6000)  GeV and b is in the range (50,1200)  GeV while in the case B, μ is in the range (500,6000)  GeV and b is in the range (25,1200)  GeV. These values are in the right range to solve the μ-problem.


## 8.3. Mass Spectrum

(A)
With such choice the gauge vector sector is completely fixed up to atβ dependence. Anyway even such a dependence can be safely ignored with a very good approximation in the new gauge sector since the mixing is strongly suppressed. So for each tβ value, we have
(8.8)MZ1≃10TeV,MZ2≃2TeV,
where with Z1 we denote the VA-like vector.
(B)
As in the previous case, we just give theZ1,2 masses
(8.9)MZ1≃10TeV,MZ2≃400GeV,
where as in the previous case Z1 is VA-like.We will not give the exact values of theZ0 mass. It is enough for our purposes to know that they are compatible with the bounds of  Section 8.1. Both case A and B are compatible with CDF bounds about Z' direct production [102].Recent LHC data have restricted the most probable range for the Higgs particle mass to be[115.5,131] GeV (ATLAS) [103] and [114.5,127] (CMS) [104]. Moreover, there are hints observed by both CMS and ATLAS of an excess of events that might correspond to decays of a Higgs particle with a mass in a range close to 125  GeV. So, in Figures 2 and 3 we give region plots showing the allowed values of 〈α〉 and tβ for case A (B) and Λ/c=(5)10TeV. The red region is the one in which Mh02|1-loop∈[124,126]GeV where the h0 mass is computed considering 1-loop corrections. Since it turns out that the top squarks have small mixing angle and considering the limit MA0≫MZ0, we have [93]
(8.10)Mh02|1-loop≃Mh02|tree+34π2sβ2yt2mt2ln(mt~1mt~2mt2)≃Mh02|tree+32π2mt4v2ln(mt~1mt~2mt2),
where Mh0|tree is the tree-level h0 mass and we used mt=ytvu/2=ytvsβ/2. There is an approximated inverse correlation between 〈α〉 and tβ in the h0 mass allowed region because the 1-loop correction in (8.10) increases for increasing values of 〈α〉 or tβ. The h0 mass allowed region is almost the same for case A and B because of two reasons:  (i)
the mixing withϕ± is suppressed,(ii)
the parametersm~hu,m~hd,μ,b in the scalar potential (4.22) are ruled by the square mass parameters gA〈α〉MVA and (gAv-)2, and the first one turns out to be dominant.Allowed〈α〉 and tanβ values for case A (up) and case B (down), Λc=5 (left) and Λc=10 (right). The red region is the one in which Mh02|1-loop∈[124,126]GeV, the magenta region is the one in which Mh02|1-loop∈[114.5,131]GeV and the blue region satisfies all the mass bounds on the sparticles (from PDG) and requires a neutralino LSP. The yellow dots are our benchmark points.
(a)
(b)
(c)
(d)Allowed〈α〉 and tanβ values for case A (up) and case B (down), Λc=5 (left) and Λc=10 (right). The red region is the one in which Mh02|1-loop∈[124,126]GeV, the magenta region is the one in which Mh02|1-loop∈[114.5,131]GeV and the blue region satisfies all the mass bounds on the sparticles (from preliminary LHC data) and requires a neutralino LSP. The yellow dots are our benchmark points.
(a)
(b)
(c)
(d)The magenta region satisfies a milder constraint on the light Higgs boson:Mh02|1-loop∈[114.5,131]  GeV. In order to be more conservative, we imposed the joint constraints of ATLAS and CMS.The blue region satisfies all the mass bounds on the sparticles and requires a neutralino LSP. We considered two possibilities: one more optimistic (Figure2) using the PDG bounds [105, 106] and one more conservative (Figure 3) using recent LHC data [107, 108]. The combination of the gluino mass bound with a neutralino LSP is a strong constraint that reduces drastically the allowed parameter space. In some cases there is not even a blue region, which means that we cannot satisfy simultaneously all the mass bounds and have a neutralino LSP, so they are completely ruled out. When the gluino mass bound is from PDG, case A is allowed; otherwise it is completely ruled out and only case B for Λ/c=10TeV presents allowed regions. We notice that case A favors low 〈α〉 values, while case B favors high 〈α〉 values. For every allowed case we choose a benchmark point (yellow spots in Figures 2 and 3):(i)
case A,Λc=5, 〈α〉=3 TeV and tβ=50 so that Mh0|1-loop≃121.6 GeV;(ii)
case A,Λc=10, 〈α〉=5 TeV and tβ=10 so that Mh0|1-loop≃124.7 GeV;(iii)
case B,Λc=10, 〈α〉=4.5 TeV and tβ=50 so that Mh0|1-loop≃125.1 GeV;(iv)
case B,Λc=10, 〈α〉=50 TeV and tβ=2.5 so that Mh0|1-loop≃130.1 GeV,and we give the full mass spectrum in Figures 4 and 5.Mass spectrum, case A,Λc=5, 〈α〉=0.3 TeV and tβ=50 (a), Λc=10, 〈α〉=0.5 TeV and tβ=10 (b).
(a)
(b)Mass spectrum, case B,Λc=10, 〈α〉=0.45 TeV and tβ=50 (a), 〈α〉=8 TeV and tβ=2.5 (b).
(a)
(b)All the benchmark points share some common features.(i)
The LSP is the lightest neutralino of the new sector: in case A it is a combination ofϕ~± and λB while in case B is almost a pure λB.(ii)
An approximated mass degeneracy ofH0, A0, and H± holds, and their masses satisfy the bounds of [96, 109].(iii)
The lightest sleptons is a sneutrino, except fortβ=50 when it is τ~1(iv)
The lightest squark isu~L, except for tβ=50 when it is b~1(v)
The first and second family left-handed squarks/sleptons are likely to be lighter than their right-handed counterparts. This is at odds with the usual MSSM cases [93].(vi)
C~1(2) is close in mass with N~1(4)MSSM. C~2 and N~4MSSM are heavier than all sfermions.(vii)
The gluino is close in mass toC~1 and N~1MSSM which are gaugino-like. Moreover it is lighter than all the squarks except for point (i). So it turns out to be long lived, specially in case B where the approximated mass degeneracy involves also the LSP. Long-lived gluinos bind with SM quarks and gluons from the vacuum during the hadronisation process and produce R-hadrons. R-hadrons are among the most interesting searches at LHC. Anyway we will come back to this point with a more detailed study in a forthcoming paper.(viii)
There is an approximated mass degeneracy betweene~R and u~R because using the charge constraints (3.16) and (8.1) we get qEc=3 and qUc≃2.9.(ix)
mϕR-<mϕ+ except for point (i).Case B points deserve some more comments.ϕ
+ and N~2,3new are out of the plot of point (iv) because they are heavier than 6 TeV. Z2 is among the lightest not SM particle, so it can decay only into SM particles, because of energy and R-parity conservation. So Z2 is long lived, because SM particles are coupled to Z2 only through the suppressed VA,B mixing or through the Higgs scalars which present a tiny mixing with ϕ±.It is not an easy task to compare the resulting spectrum we get for our model with those related to the rich zoology of supersymmetry breaking scenarios. It is worth to stress anyway that the two representative spectrums showed in Figure5 which encode the key features of our scenarios listed above are not reproduced in any of the benchmark points showed in [110, 111].


## 9. Conclusions

In this paper we presented a viable mechanism to generate soft supersymmetry breaking terms in the framework of a minimal supersymmetric anomalous extension of the SM. The crucial ingredient is a non perturbative term in the superpotential (2.3) which couples the Stückelberg field S to the Higgs sector. This term is related to the generation of a suitable μ and b terms (see (4.5) and (4.6)) in the low-energy effective action when the Stückelberg gets vev. We argued about the origin of this term from an exotic instanton in an intersecting D-brane setup. We computed the spectrum of our model as a function of the saxion vev 〈α〉 and for different choices of the remaining free parameters. We checked our results against known phenomenological bounds, namely, current lower bounds on the mass of the scalar and fermionic superpartners. We analyzed a scenario in which the anomalous sector is the source of the soft supersymmetry breaking terms while the corresponding vector and Stückelberg multiplets are not present in the low-energy effective action. For what concerns the non anomalous sector, we took into account two different cases (dubbed case A and case B).As we stated in Section8, by applying some phenomenological constraints we were able to derive some bounds on the saxion vev 〈α〉, which is the relevant parameter setting the mass scale of the scalars. The strongest constraints on 〈α〉 and tβ come from the combined requirement of Mh02|1-loop∈[124,126]GeV or ([114.5,131]GeV), a neutralino LSP, and that all mass bounds (specially the gluino one) are fulfilled. In Figure 2 (pre-LHC bounds) and 8.3 (preliminary LHC bounds), we summarize the allowed regions for 〈α〉. In the first case, by requiring a phenomenological appealing neutralino LSP, we get an allowed 〈α〉 of few TeV up to 10 TeV for the A and B scenarios respectively. In the second case (preliminary LHC bounds), we get that only the B scenario is allowed with 〈α〉≳5 TeV. These results can be seen as a bound that a concrete D-brane model has to satisfy. We deserve this analysis for future work.In Figure5 we explicitly showed two benchmark mass spectrums for our model with 〈α〉 and tβ which fulfill the above bounds. The cases shared different peculiar features: the LSP is the lightest neutralino of the new sector, there is a near mass degeneracy between H0, A0 and H±, and e~R and u~R the lightest sleptons is a sneutrino except for tβ=50 when it is stau, the lightest squark is a u~L except for tβ=50, when it is a sbottom, the first and second family left-handed squarks/sleptons are typically lighter than their right-handed counterparts. Moreover in case B the gluino is long lived and can produce R-hadrons. It turns out that these features are not reproduced in any of the widely studied benchmark points presented in [110, 111].



---

*Source: 903106-2012-08-16.xml*