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.

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 bosons 𝑍′ 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 π‘ˆ(1) gauge groups as well as 𝑁=1 supersymmetry softly broken roughly at the TeV scale. In particular in string theory the presence of extra anomalous π‘ˆ(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 π‘ˆ(1) (see [59–67] for other anomalous π‘ˆ(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 π‘ˆ(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 Section 2 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, π‘ˆ(1)𝐴 and π‘ˆ(1)𝐡. The first one is anomalous while the second one is anomaly-free. This assumption is quite generic since in models with several anomalous π‘ˆ(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) The anomaly cancellation of the π‘ˆ(1)𝐴 gauge group is achieved by the four-dimensional analogue of the higher-dimensional GS mechanism which involves the StΓΌckelberg superfield 𝑆=𝑠+2πœƒπœ“π‘†+πœƒ2𝐹𝑆 transforming as a shift π‘†β†’π‘†βˆ’2𝑖𝑀𝑉𝐴Λ,(2.2) where 𝑀𝑉𝐴 is a mass parameter related to the anomalous π‘ˆ(1)𝐴 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 π‘ˆ(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 𝐸=𝑀𝑉𝐴 is given by π‘Š=π‘ŠMSSM+πœ†π‘’βˆ’π‘˜π‘†π»π‘’π»π‘‘+π‘šΞ¦+Ξ¦βˆ’,(2.3) where π‘ŠMSSM is given by π‘ŠMSSM=π‘¦π‘’π‘–π‘—π‘„π‘–π‘ˆπ‘π‘—π»π‘’βˆ’π‘¦π‘‘π‘–π‘—π‘„π‘–π·π‘π‘—π»π‘‘βˆ’π‘¦π‘’π‘–π‘—πΏπ‘–πΈπ‘π‘—π»π‘‘,(2.4) which is the usual MSSM superpotential without the πœ‡-term which is forbidden for a generic choice of the charges π‘žπ»π‘’ and π‘žπ»π‘‘. 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 π‘ˆ(1)𝐴 and π‘ˆ(1)𝐡. 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 π‘ˆ(1)𝐴 and π‘ˆ(1)𝐡), the two π‘ˆ(1)𝐡 fields Φ± couple to the gauge field strength π‘Šπ›Όπ‘Ž through the dimension six effective operator ℒ𝑔=π‘π‘ŽΞ¦+Ξ¦βˆ’Ξ›2π‘Šπ›Όπ‘Žπ‘Šπ›Όπ‘Ž,(2.5) where π‘Ž=𝐴,𝐡,π‘Œ,2,3, Ξ› is the cut-off scale of the theory while π‘π‘Ž are constants that have to be computed in the UV completion of the theory.

The 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 𝐻𝑒𝐻𝑑 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 𝐢3 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 𝑁=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 π‘’βˆ’Vol𝐸2/𝑔𝑠,(2.6) where Vol𝐸2 is the volume of the 3-cycle in the CY wrapped by a 𝐸2-brane measured in string units while 𝑔𝑠 is the string coupling. Such exponential factor is independent of the 𝑑=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 𝑔𝑠=1 if Vol𝐸2=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 𝛾𝑖3), can be written as π‘ˆπ‘–=π‘’βˆ’πœ‘ξ€œπ›Ύπ‘–3Ξ©3ξ€œ+𝑖𝛾𝑖3𝐢3,(2.7) where πœ‘ is the dilaton, Ξ©3 is the CY volume 3-form (which is a complex form), and 𝐢3 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 𝐸2 instanton is formally given by π‘ŠβˆΌπ‘›ξ‘π‘–=1Ξ¦π‘Žπ‘–,π‘π‘–π‘’βˆ’π‘†πΈ2,(2.8) where Ξ¦π‘Žπ‘–,𝑏𝑖 are chiral superfields localized at the intersection of two D6-branes described by open strings while 𝑆𝐸2 denotes the instanton classical action: π‘’βˆ’π‘†πΈ2ξƒ¬βˆ’=exp2πœ‹π‘™3𝑠1π‘”π‘ ξ€œπ›Ύξ€·Ξ©Re3ξ€Έξ€œβˆ’π‘–π›ΎπΆ3ξ‚Άξƒ­.(2.9) 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 𝑔𝑠=π‘’βŸ¨πœ‘βŸ© and the axionic field is an exponential. Any other dependence can be excluded due to the shift transformation (2.2).

3. Supersymmetry Breaking

The D-term contribution of the π‘ˆ(1)𝐴 vector multiplet 𝑉𝐴 relevant to supersymmetry breaking is given, in the limit of vanishing kinetic mixing π›Ώπ‘Œπ΄,𝛿𝐴𝐡=0, by the following lagrangian: 1β„’=2𝐷𝐴𝐷𝐴+ξ“π‘–π‘”π΄π‘žπ‘–πœ™β€ π‘–π·π΄πœ™π‘–+πœ‰π·π΄,(3.1) where the sum is extended to all the scalars charged under the π‘ˆ(1)𝐴. There is no D-term contribution related to the π‘ˆ(1)𝐡 except that of πœ™Β± since all the MSSM chiral fields are uncharged under π‘ˆ(1)𝐡 (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 β„’axion=14𝑆+𝑆†+2𝑀𝑉𝐴𝑉𝐴2|||πœƒ2πœƒ2+β‹―=𝑀𝑉𝐴𝑆+𝑆†𝑉𝐴||πœƒ2πœƒ2+β‹―=𝑀𝑉𝐴𝛼𝐷𝐴+β‹―,(3.2) where in the last line 𝛼 denotes the real part of the lowest component of the StΓΌckelberg chiral multiplet 𝑠=𝛼+π‘–πœ‘. 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: 116𝑔2π‘Žπœπ‘Ž=116̃𝑔2π‘Žπœπ‘Žβˆ’12π‘π‘Žπ‘ŽβŸ¨π›ΌβŸ©,(3.3) where the gauge factors πœπ‘Ž take the values 1,1,1,1/2,1/2, and the π‘π‘Žπ‘Ž constants are given in (A.2). The tree-level FI term is then given by πœ‰=π‘€π‘‰π΄βŸ¨π›ΌβŸ©.(3.4) 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 𝐷𝐴: 𝐷𝐴+ξ“π‘–π‘”π΄π‘žπ‘–πœ™β€ π‘–πœ™π‘–+πœ‰=0,(3.5) where the index 𝑖 runs over all chiral superfields. The D-term contribution to the scalar potential is given by π‘‰ξ‚€πœ™π‘–,πœ™β€ π‘–ξ‚=12ξƒ©πœ‰+π‘”π΄ξ“π‘–π‘žπ‘–||πœ™π‘–||2ξƒͺ2.(3.6) The quadratic part gives the scalar mass term ξ“π‘–πœ‰π‘”π΄π‘žπ‘–||πœ™π‘–||2=ξ“π‘–π‘š2𝑖||πœ™π‘–||2,(3.7) where we have defined π‘š2𝑖=πœ‰π‘”π΄π‘žπ‘–=βŸ¨π›ΌβŸ©π‘”π΄π‘€π‘‰π΄π‘žπ‘–=π‘žπ‘–π‘š2πœ‰,(3.8) with π‘š2πœ‰=βŸ¨π›ΌβŸ©π‘”π΄π‘€π‘‰π΄=π‘”π΄πœ‰.(3.9) The typical scale for the mass π‘šπœ‰ is of the order of few hundreds of GeV if π‘€π‘‰π΄βˆΌβŸ¨π›ΌβŸ©βˆΌ1 TeV and π‘”π΄βˆΌ0.1. It is interesting to note that in this scenario a low sub-TeV supersymmetry breaking scale π‘šπœ‰ is due to the StΓΌckelberg mechanism which gives mass to 𝑉𝐴. This is the most important difference with the scenario proposed in [89], where the scale π‘šπœ‰ is dynamically generated by some dynamics in a strong coupling regime.

Mass terms for the gauginos, that is, πœ†π‘Žπœ†π‘Ž, 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 π‘€π‘Ž=π‘π‘Žξ«πΉ+πœ™βˆ’ξ¬+ξ«πΉβˆ’πœ™+Λ2=π‘π‘Žπ‘šξ€·π‘£2++𝑣2βˆ’ξ€Έ2Ξ›2,(3.10) where 𝑣±/√2=βŸ¨πœ™Β±βŸ© and where in the right-hand side we have used the F-term equations of motion for 𝐹±𝐹±=βˆ’πœ•π‘Šβˆ—πœ•πœ™Β±βˆ—=βˆ’π‘šπœ™βˆ“βˆ—,(3.11) having assumed π‘š real without any loss in generality. We assume π‘π‘Ž=𝑐 for each π‘Ž. 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 π‘€π‘ƒβ†’βˆž, the F-term induced contribution to the scalar masses π‘š2π‘–βˆΌξ«πΉΒ±ξ¬π‘€2𝑃(3.12) vanishes leaving (3.8) as the leading contribution.

The requirement of gauge invariance of the superpotential implies the following constraints on the π‘ˆ(1)𝐴 charges π‘žπ‘ˆπ‘=βˆ’π‘žπ‘„βˆ’π‘žπ»π‘’,π‘žπ·π‘=βˆ’π‘žπ‘„βˆ’π‘žπ»π‘‘,π‘žπΈπ‘=βˆ’π‘žπΏβˆ’π‘žπ»π‘‘,π‘ž(3.13)π‘˜=𝐻𝑒+π‘žπ»π‘‘2𝑀𝑉𝐴.(3.14) As we said at the beginning of this section, we assume that the net kinetic mixing between π‘ˆ(1)π‘Œ and π‘ˆ(1)𝐴 vanishes (we postpone the discussion about the kinetic mixing between π‘ˆ(1)𝐴 and π‘ˆ(1)𝐡 to the next section). There are two contributions for the π‘ˆ(1)π‘Œβˆ’π‘ˆ(1)𝐴 kinetic mixing: the 1-loop mixing π›Ώπ‘Œπ΄ and π‘π‘Œπ΄ coming from the GS coupling π‘†π‘Šπ‘Œπ‘Šπ΄ (see (A.1)). The following conditions imply a bound on the charges π›Ώπ‘Œπ΄ξ“=0β‡’π‘“π‘žπ‘“π‘Œπ‘“π‘=0,π‘Œπ΄ξ“=0β‡’π‘“π‘ž2π‘“π‘Œπ‘“=0,(3.15) where the sum is extended over all the chiral fermions in the theory. Constraints (3.15), can be solved in terms of π‘žπ‘„ and π‘žπΏ. By using conditions (3.13), we get π‘žπΏ=14ξ€·3π‘žπ»π‘’βˆ’4π‘žπ»π‘‘ξ€Έ,π‘žπ‘„1=βˆ’ξ€·125π‘žπ»π‘’βˆ’2π‘žπ»π‘‘ξ€Έ.(3.16) The positive squared mass condition for the sfermions π‘š2𝑓=π‘”π΄π‘žπ‘“π‘€π‘‰π΄βŸ¨π›ΌβŸ©>0,(3.17) implies π‘žπ‘“>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 π‘žπ»π‘’5<0,2π‘žπ»π‘’<π‘žπ»π‘‘<34π‘žπ»π‘’.(3.18)

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: π‘Šinst||πœƒ2=πœ†π‘’βˆ’π‘˜π‘†π»π‘’π»π‘‘||πœƒ2=πœ†π‘’βˆ’π‘˜π‘ β„Žπ‘’πΉπ‘‘+πœ†π‘’βˆ’π‘˜π‘ πΉπ‘’β„Žπ‘‘βˆ’πœ†π‘˜π‘’βˆ’π‘˜π‘ πΉπ‘†β„Žπ‘’β„Žπ‘‘+√2πœ†π‘’βˆ’π‘˜π‘ π‘˜ξ‚€β„Žπ‘’πœ“π‘†ξ‚β„Žπ‘‘+β„Žπ‘‘πœ“π‘†ξ‚β„Žπ‘’ξ‚βˆ’πœ†π‘’βˆ’π‘˜π‘ β„Žπ‘’β„Žπ‘‘πœ“π‘†πœ“π‘†,(4.1) where 𝐹𝑒,𝑑 are the F-terms of 𝐻𝑒,𝑑. Solving the F-terms equations for 𝐻𝑒 and 𝐻𝑑, we get the following contributions for the instanton-induced term in the scalar potential: 𝑉inst=2πœ†2π‘’βˆ’2π‘˜π›Όβ„Žβ€ π‘’β„Žπ‘’+2πœ†2π‘’βˆ’2π‘˜π›Όβ„Žβ€ π‘‘β„Žπ‘‘+πœ†π‘˜π‘’βˆ’π‘˜π›Όξ€·π‘’βˆ’π‘–π‘˜πœ‘πΉπ‘†β„Žπ‘’β„Žπ‘‘ξ€Έ.+h.c.(4.2) 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 𝐹𝑆. By supersymmetry the saxion field 𝛼, being part of the Sẗuckelberg multiplet, has a tree-level mass 𝑀𝑉𝐴. Thus if we want to consider a frozen dynamics for 𝛼 at the TeV scale we have to assume a mass parameter for the anomalous π‘ˆ(1)𝐴 just slightly above the TeV scale; that is, 𝑀𝑉𝐴>1 TeV. In this way the effective instanton-induced potential at a scale 𝐸≃1 TeV is thus given by 𝑉inst=2πœ†2π‘’βˆ’2π‘˜βŸ¨π›ΌβŸ©β„Žβ€ π‘’β„Žπ‘’+2πœ†2π‘’βˆ’2π‘˜βŸ¨π›ΌβŸ©β„Žβ€ π‘‘β„Žπ‘‘+πœ†π‘˜π‘’βˆ’π‘˜βŸ¨π›ΌβŸ©ξ€·βŸ¨πΉπ‘†βŸ©π‘’βˆ’π‘–π‘˜πœ‘β„Žπ‘’β„Žπ‘‘ξ€Έ.+h.c.(4.3) The first two terms are πœ‡-terms while the third one is a b-term. The complete effective scalar potential is given by ξ‚€||πœ‡||𝑉=2+π‘š2β„Žπ‘’||β„Žξ‚ξ‚€0𝑒||2+||β„Ž+𝑒||2+ξ‚€||πœ‡||2+π‘š2β„Žπ‘‘||β„Žξ‚ξ‚€0𝑑||2+||β„Žβˆ’π‘‘||2+ξ‚€|π‘š|2+π‘š2πœ™+||πœ™+||2+ξ‚€|π‘š|2+π‘š2πœ™βˆ’ξ‚||πœ™βˆ’||2+ξ€Ίπ‘π‘’βˆ’π‘–π‘˜πœ‘ξ€·β„Ž+π‘’β„Žβˆ’π‘‘βˆ’β„Ž0π‘’β„Ž0𝑑+1+h.c.8𝑔22+𝑔2π‘Œξ€Έξ€·|β„Ž0𝑒|2+|β„Ž+𝑒|2βˆ’|β„Ž0𝑑|2βˆ’|β„Žβˆ’π‘‘|2ξ€Έ2+12𝑔22||β„Ž+π‘’β„Žπ‘‘0βˆ—+β„Ž0π‘’β„Žπ‘‘βˆ’βˆ—||2+12𝑔2π΄ξ‚ƒπ‘žπ»π‘’ξ‚€||β„Ž0𝑒||2+||β„Ž+𝑒||2+π‘žπ»π‘‘ξ‚€||β„Ž0𝑑||2+||β„Žβˆ’π‘‘||2+||πœ™+||2βˆ’||πœ™βˆ’||2ξ‚„2+12𝑔2𝐡||πœ™+||2βˆ’||πœ™βˆ’||2ξ‚„2,(4.4) where βˆšπœ‡=2πœ†π‘’βˆ’π‘˜βŸ¨π›ΌβŸ©,(4.5)𝑏=πœ†π‘˜π‘’βˆ’π‘˜βŸ¨π›ΌβŸ©βŸ¨πΉπ‘†βŸ©.(4.6) 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 π‘ˆ(1)𝐴 term: π‘š2β„Žπ‘’=π‘žπ»π‘’π‘š2πœ‰,π‘š(4.7)2β„Žπ‘‘=π‘žπ»π‘‘π‘š2πœ‰,π‘š(4.8)2πœ™+=π‘š2πœ‰,π‘š(4.9)2πœ™βˆ’=βˆ’π‘š2πœ‰,(4.10) with π‘š2πœ‰ given by (3.8). The scalar potential depends on the following new parameters: βŸ¨π›ΌβŸ©, βŸ¨πΉπ‘†βŸ©, πœ†, π‘š, 𝑔𝐴,𝐡, π‘žπ»π‘’,𝑑, 𝑀𝑉𝐴.

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 𝑉0=ξ‚€||πœ‡||2+π‘š2β„Žπ‘’ξ‚||β„Ž0𝑒||2+ξ‚€||πœ‡||2+π‘š2β„Žπ‘‘ξ‚||β„Ž0𝑑||2βˆ’ξ€·π‘π‘’βˆ’π‘–π‘˜πœ‘β„Ž0π‘’β„Ž0𝑑+ξ‚€+h.c.|π‘š|2+π‘š2πœ™+||πœ™+||2+ξ‚€|π‘š|2+π‘š2πœ™βˆ’ξ‚||πœ™βˆ’||2+18𝑔22+𝑔2π‘Œξ€Έξ‚€||β„Ž0𝑒||2βˆ’||β„Ž0𝑑||22+12𝑔2π΄ξ‚€π‘žπ»π‘’||β„Ž0𝑒||2+π‘žπ»π‘‘||β„Ž0𝑑||2+||πœ™+||2βˆ’||πœ™βˆ’||22+12𝑔2𝐡||πœ™+||2βˆ’||πœ™βˆ’||2ξ‚„2.(4.11) 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 πœ•π‘‰0/πœ•π‘§π‘–=0 where the scalar field 𝑧𝑖 runs over {πœ‘,β„Ž0𝑒,β„Ž0𝑑,πœ™+,πœ™βˆ’}. The conditions for having a nontrivial minimum boil down to the same condition of the MSSM 𝑏2>ξ‚€||πœ‡||2+π‘š2β„Žπ‘’||πœ‡||2+π‘š2β„Žπ‘‘ξ‚.(4.12) Moreover in order to generate a mass term for the gauginos (see (3.10)), the condition π‘£βˆ’β‰ 0 must hold since 𝑣+=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: |π‘š|2+π‘š2πœ™βˆ’<0,(4.13) The minimum is attained at πœ‘=πœ™+=0. Actually since the potential for the axion πœ‘ is periodic, the minimum condition holds for πœ‘=2π‘›πœ‹/π‘˜ with π‘›βˆˆβ„€. All these minima are physically equivalent, and thus we arbitrarily choose 𝑛=0. The remaining three conditions imply the following constraints on the parameters: π‘š2β„Žπ‘‘+πœ‡2βˆ’π‘π‘‘π›½+18𝑔2π‘Œ+𝑔22𝑣2𝑐2𝛽+12𝑔2π΄π‘žπ»π‘‘ξ‚ƒπ‘£2ξ‚€π‘žπ»π‘‘π‘2𝛽+π‘žπ»π‘’π‘ 2π›½ξ‚βˆ’π‘£2βˆ’ξ‚„π‘š=0,2β„Žπ‘’+πœ‡2βˆ’π‘π‘‘π›½βˆ’1βˆ’18𝑔2π‘Œ+𝑔22𝑣2𝑐2𝛽+12𝑔2π΄π‘žπ»π‘’ξ‚ƒπ‘£2ξ‚€π‘žπ»π‘‘π‘2𝛽+π‘žπ»π‘’π‘ 2π›½ξ‚βˆ’π‘£2βˆ’ξ‚„ξ€·π‘”=0,2𝐴+𝑔2𝐡𝑣2βˆ’βˆ’π‘”2𝐴𝑣2ξ‚€π‘žπ»π‘‘π‘2𝛽+π‘žπ»π‘’π‘ 2𝛽+2|π‘š|2+π‘š2πœ™βˆ’ξ‚=0,(4.14) where we have defined in order to keep a compact notation 𝑐𝛽=cos𝛽,𝑠𝛽=sin𝛽,𝑑𝛽=tan𝛽,𝑐2𝛽=cos(2𝛽),𝑠2𝛽=sin(2𝛽),(4.15) and as usual as tan𝛽=𝑣𝑒/𝑣𝑑.

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) we can neglect all the 𝑔𝐴𝑣 terms. With this approximation, the minima equations read ξ‚π‘š2β„Žπ‘‘+πœ‡2βˆ’π‘π‘‘π›½+18𝑔2π‘Œ+𝑔22𝑣2𝑐2𝛽=0,(4.17)ξ‚π‘š2β„Žπ‘’+πœ‡2βˆ’π‘π‘‘π›½βˆ’1βˆ’18𝑔2π‘Œ+𝑔22𝑣2𝑐2𝛽𝑔=0,(4.18)2𝐴+𝑔2𝐡𝑣2βˆ’ξ‚€+2|π‘š|2+π‘š2πœ™βˆ’ξ‚=0,(4.19) where we have defined ξ‚π‘š2β„Žπ‘‘=π‘š2β„Žπ‘‘βˆ’12𝑔2π΄π‘žπ»π‘‘π‘£2βˆ’,ξ‚π‘š2β„Žπ‘’=π‘š2β„Žπ‘’βˆ’12𝑔2π΄π‘žπ»π‘’π‘£2βˆ’.(4.20) Equations (4.17) and (4.18) have the same functional form as in the MSSM case. Moreover π‘£βˆ’ 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 ξ‚€||πœ‡||𝑉≃2+π‘š2β„Žπ‘’||β„Žξ‚ξ‚€0𝑒||2+||β„Ž+𝑒||2+ξ‚€|πœ‡|2+π‘š2β„Žπ‘‘||β„Žξ‚ξ‚€0𝑑||2+||β„Žβˆ’π‘‘||2+ξ€Ίπ‘π‘’βˆ’π‘–π‘˜πœ‘ξ€·β„Ž+π‘’β„Žβˆ’π‘‘βˆ’β„Ž0π‘’β„Ž0𝑑+1+h.c.8𝑔22+𝑔2π‘Œξ€Έξ‚€||β„Ž0𝑒||2+||β„Ž+𝑒||2βˆ’||β„Ž0𝑑||2βˆ’||β„Žβˆ’π‘‘||22+12𝑔22||β„Ž+π‘’β„Žπ‘‘0βˆ—+β„Ž0π‘’β„Žπ‘‘βˆ’βˆ—||2+12𝑔2π΄ξ‚ƒπ‘žπ»π‘’ξ‚€||β„Ž0𝑒||2+||β„Ž+𝑒||2+π‘žπ»π‘‘ξ‚€||β„Ž0𝑑||2+||β„Žβˆ’π‘‘||2ξ‚βˆ’12𝑣2βˆ’ξ‚„2,(4.21) neglecting further constant terms in π‘£βˆ’. Close to the minima, the relevant term in the last line of (4.21) is the double product of the Higgs part with the 𝑣2βˆ’ term. Hence by using (4.16) we finally get π‘‰β„Ž,πœ‘β‰ƒξ‚€||πœ‡||2+ξ‚π‘š2β„Žπ‘’||β„Žξ‚ξ‚€0𝑒||2+||β„Ž+𝑒||2+ξ‚€||πœ‡||2+ξ‚π‘š2β„Žπ‘‘||β„Žξ‚ξ‚€0𝑑||2+||β„Žβˆ’π‘‘||2+ξ€Ίπ‘π‘’βˆ’π‘–π‘˜πœ‘ξ€·β„Ž+π‘’β„Žβˆ’π‘‘βˆ’β„Ž0π‘’β„Ž0𝑑+1+h.c.8𝑔22+𝑔2π‘Œξ€Έξ‚€||β„Ž0𝑒||2+||β„Ž+𝑒||2βˆ’||β„Ž0𝑑||2βˆ’||β„Žβˆ’π‘‘||22+12𝑔22||β„Ž+π‘’β„Žπ‘‘0βˆ—+β„Ž0π‘’β„Žπ‘‘βˆ’βˆ—||2.(4.22) 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 is 𝑑𝛽≳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.) ξ‚π‘š2β„Žπ‘’<ξ‚π‘š2β„Žπ‘‘. By using (4.20) we get π‘”π΄π‘žπ»π‘’ξ‚€π‘€π‘‰π΄1βŸ¨π›ΌβŸ©βˆ’2𝑔𝐴𝑣2βˆ’ξ‚<π‘”π΄π‘žπ»π‘‘ξ‚€π‘€π‘‰π΄1βŸ¨π›ΌβŸ©βˆ’2𝑔𝐴𝑣2βˆ’ξ‚.(4.23) By assuming 𝑀𝑉𝐴>1 TeV, π‘£βˆ’ in the TeV range, π‘”π΄βˆΌπ‘‚(0.1), the term between brackets is positive and we get the following constraint: π‘žπ»π‘’<π‘žπ»π‘‘(4.24) for the π‘ˆ(1)𝐴 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 𝑀2πœ™+=2|π‘š|2.(4.25) The same holds for the imaginary part of πœ™βˆ’ which becomes the longitudinal mode of the gauge vector 𝑍2. The mass matrix for the real scalar fields {πœ‘,Im(β„Ž0𝑒),Im(β„Ž0𝑑)} is given by ℳ𝑆(Im)=βŽ›βŽœβŽœβŽπ‘π‘‘π›½β‹―β‹―π‘π‘π‘‘π›½βˆ’1β‹―βˆ’π‘π‘˜π‘£π‘ π›½βˆ’π‘π‘˜π‘£π‘π›½π‘π‘˜2𝑣2π‘π›½π‘ π›½βŽžβŽŸβŽŸβŽ .(4.26) The determinant of this matrix is zero. Two eigenvalues are zero which correspond to the Goldstone modes of 𝑍0 and 𝑍1. The physical massive state is an axi-Higgs state with mass given by 𝑀2𝐴0=2𝑏𝑠2π›½βŽ‘βŽ’βŽ’βŽ£11βˆ’ξ€·π‘ž16𝐻𝑒+π‘žπ»π‘‘ξ€Έ2𝑣2𝑀2𝑉𝐴𝑠22π›½βŽ€βŽ₯βŽ₯⎦,(4.27) where we used the relation (3.14). The mass matrix for the real scalar fields {Re(β„Ž0𝑒), Re(β„Ž0𝑑), πœ™βˆ’π‘…β‰‘Re(πœ™βˆ’)} reads as ℳ𝑆(Re)=βŽ›βŽœβŽœβŽœβŽξ‚€14𝑔2πΈπ‘Š+𝑔2π΄π‘žπ»π‘‘2𝑣2𝑐2𝛽+𝑏𝑑𝛽1β‹―β‹―βˆ’π‘βˆ’4𝑔2πΈπ‘Šβˆ’π‘”2π΄π‘žπ»π‘‘π‘žπ»π‘’ξ‚π‘£2𝑐𝛽𝑠𝛽14𝑔2πΈπ‘Š+𝑔2π΄π‘žπ»π‘’2𝑣2𝑠2𝛽+π‘π‘‘π›½βˆ’1β‹―βˆ’π‘”2π΄π‘žπ»π‘‘π‘£π‘£βˆ’π‘π›½βˆ’π‘”2π΄π‘žπ»π‘’π‘£π‘£βˆ’π‘ π›½ξ€·π‘”2𝐴+𝑔2𝐡𝑣2βˆ’βŽžβŽŸβŽŸβŽŸβŽ ,(4.28) where 𝑔2πΈπ‘Š=(𝑔2π‘Œ+𝑔22). 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: 𝑀2β„Ž0,𝐻0≃12βŽ›βŽœβŽœβŽ2𝑏𝑠2π›½βˆ“ξƒŽξ‚΅2𝑏𝑠2π›½βˆ’14𝑔2π‘Œ+𝑔22𝑣2ξ‚Ά2𝑔+2𝑏2π‘Œ+𝑔22𝑣2𝑠2π›½βŽžβŽŸβŽŸβŽ ,𝑀2πœ™βˆ’π‘…β‰ƒξ€·π‘”2𝐴+𝑔2𝐡𝑣2βˆ’.(4.29)

The charged sector is unchanged with respect to the MSSM, so 𝑀2𝐻±=2𝑏𝑠2𝛽+𝑀2π‘Š.(4.30)

As in the standard MSSM case, the mass of the lightest Higgs π‘€β„Ž0 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 𝐷𝐴-term corrections 𝑀2β„Ž0<14𝑔2π‘Œ+𝑔22𝑣2𝑐22𝛽+14𝑔2𝐴𝑣2ξ€Ίπ‘žπ»π‘‘+π‘žπ»π‘’+ξ€·π‘žπ»π‘‘βˆ’π‘žπ»π‘’ξ€Έπ‘2𝛽2,(4.31) where the first term is the MSSM bound. In principle, for arbitrary high values of π‘”π΄π‘žπ»π‘‘, π‘”π΄π‘žπ»π‘’, we get an increasing upper bound. However, as in the standard MSSM case, 𝑀2β„Ž0 undergoes to relatively drastic quantum corrections [93]. Hence in Section 8 we consider tree-level masses for all the particles except for β„Ž0 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 ξ«πœ™Β±ξ¬=π‘£Β±βˆš2,ξ«β„Ž0𝑒,𝑑=𝑣𝑒,π‘‘βˆš2.(5.1) The neutral vector square mass matrix in the base (𝑉𝐡,𝑉𝐴,π‘‰π‘Œ,𝑉32) is ℳ𝑉=βŽ›βŽœβŽœβŽœβŽœβŽœβŽπ‘”2𝐡𝑣2πœ™π‘”β‹―β‹―β‹―π΄π‘”π΅π‘£2πœ™π‘”2𝐴𝑐2π›½π‘žπ»π‘‘2+𝑠2π›½π‘žπ»π‘’2𝑣2+𝑣2πœ™ξ‚„+𝑀2𝑉𝐴01β‹―β‹―2π‘”π΄π‘”π‘Œπ‘žπ»(𝛽)𝑣214𝑔2π‘Œπ‘£2β‹―10βˆ’2𝑔𝐴𝑔2π‘žπ»(𝛽)𝑣2βˆ’14π‘”π‘Œπ‘”2𝑣214𝑔22𝑣2⎞⎟⎟⎟⎟⎟⎠,(5.2) where 𝑣2πœ™=𝑣2++𝑣2βˆ’,𝑣2=𝑣2𝑒+𝑣2𝑑,π‘žπ»ξ‚€π‘ (𝛽)=2π›½π‘žπ»π‘’βˆ’π‘2π›½π‘žπ»π‘‘ξ‚.(5.3) By taking 𝑀𝑉𝐴>1 TeV (see Section 4), 𝑉𝐴 can be considered as decoupled from the low-energy gauge sector (namely, 𝐸≲1 TeV), and we can ignore with very good approximation any mixing term (The kinetic mixing between π‘ˆ(1)𝐴 and π‘ˆ(1)𝐡 deserves some comment, in particular if we relax the 𝑀𝑉𝐴>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(π‘žπ΄π‘žπ΅)β‰ 0 such a mixing arises at the 1-loop level. In such a case it can be assumed that the two π‘ˆ(1)'s are in the kinetic diagonalized basis with Tr(π‘žπ΄π‘žπ΅)=0 thanks to some additional heavy chiral multiplet charged under both π‘ˆ(1)𝐴 and π‘ˆ(1)𝐡. These multiplets generate a counterterm in the effective theory that cancels against 𝛿𝐴𝐡 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 𝑉𝐴. From now on we will apply this approximation.

Since 𝑉𝐡 is a hidden gauge boson, it is decoupled from the SM sector. The charged vector sector is unchanged with respect to the MSSM, so π‘ŠΒ±πœ‡=𝑉21πœ‡βˆ“π‘–π‘‰22πœ‡βˆš2,𝑀2π‘Š=14𝑔22𝑣2.(5.4)

6. Neutralinos

In comparison with the standard MSSM, we now have five new neutral fermionic fields: πœ“π‘†, πœ†π΄, πœ†π΅, ξ‚πœ™Β±. However under the assumption 𝑀𝑉𝐴>1 TeV, πœ“π‘† and πœ†π΄ are not in the low-energy sector because of the 𝑀𝑉𝐴 mass term (we stress that the πœ“π‘†βˆ’πœ†π΄ sector presents a different parameters’ choice with respect to [97–99], where we realized a scenario in which the mixing between πœ“π‘† and πœ†π΄ was suppressed). Thus we have β„’neutralinomass1=βˆ’2ξ€·πœ“0ξ€Έπ‘‡β„³ξ‚π‘πœ“0+h.c.,(6.1) where ξ€·πœ“0𝑇=ξ‚€πœ†π΅,ξ‚πœ™βˆ’,ξ‚πœ™+,πœ†π‘Œ,πœ†02,ξ‚β„Ž0𝑑,ξ‚β„Ž0𝑒.(6.2) In this basis the neutralino mass matrix ℳ𝑁 is written as ℳ𝑁=βŽ›βŽœβŽœβŽœβŽœβŽœβŽœβŽœβŽπ‘€π΅β‹―β‹―β‹―β‹―β‹―β‹―βˆ’π‘”π΅π‘£βˆ’0β‹―β‹―β‹―β‹―β‹―0βˆ’π‘š0β‹―β‹―β‹―β‹―000𝑀1β‹―β‹―β‹―0000𝑀2𝑔⋯⋯000βˆ’1𝑣𝑑2𝑔2𝑣𝑑2𝑔0β‹―0001𝑣𝑒2βˆ’π‘”2𝑣𝑒2βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽŸβŽ βˆ’πœ‡0,(6.3) where πœ‡ is given in (4.5). We recall that gaugino masses arise from the Dvali-Pomarol term (2.5).

ℳ𝑁 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 πœ†π΅ 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 π‘£βˆ’ different from zero, while 𝑣+=0.

Finally there are also corrections coming from the anomalous axino couplings: F-term couplings of the type π‘π‘Žπ‘ŽβŸ¨πΉπ‘†βŸ©πœ†π‘Žπœ†π‘Ž, D-term couplings of the type π‘π‘Žπ‘Žπœ“π‘†πœ†π‘ŽβŸ¨π·π‘ŽβŸ©, and corrections coming from the superpotential term π‘’βˆ’π‘˜π‘†π»π‘’π»π‘‘+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 Section 8 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 π‘š3/2βˆΌπ‘‚β€‰(TeV) in the limit π‘€π‘ƒβ†’βˆž.

7. Sfermion Masses

The sfermion masses receive several contributions. As we have seen in Section 3, 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 π·π‘Œ and 𝐷2 term correction from the Higgs sector. Moreover there are 𝐷𝐴 term corrections from the Higgs and πœ™βˆ’ sector. As an aside, the appearance of such terms in the low-energy action, given our assumption 𝑀𝑉𝐴>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 π‘š2Μƒπ‘’πΏβ‰ƒπ‘š2̃𝑐𝐿=π‘š2𝑄+ξ‚€13𝑔2π‘Œβˆ’π‘”22Δ𝑣28+π‘žπ‘„ξ‚π‘š2𝐷𝐴,π‘š2Μƒπ‘’π‘…β‰ƒπ‘š2̃𝑐𝑅=π‘š2ξ‚π‘ˆπ‘βˆ’π‘”2π‘ŒΞ”π‘£26+π‘žπ‘ˆπ‘ξ‚π‘š2𝐷𝐴,π‘š2ξ‚π‘‘πΏβ‰ƒπ‘š2̃𝑠𝐿=π‘š2𝑄+ξ‚€13𝑔2π‘Œ+𝑔22Δ𝑣28+π‘žπ‘„ξ‚π‘š2𝐷𝐴,π‘š2ξ‚π‘‘π‘…β‰ƒπ‘š2̃𝑠𝑅=π‘š2𝐷𝑐+𝑔2π‘ŒΞ”π‘£212+π‘žπ·π‘ξ‚π‘š2𝐷𝐴,π‘š2Μƒπœˆπ‘’=π‘š2Μƒπœˆπœ‡=π‘š2ξ‚πΏβˆ’ξ€·π‘”2π‘Œ+𝑔22Δ𝑣28+π‘žπΏξ‚π‘š2𝐷𝐴,π‘š2Μƒπ‘’πΏβ‰ƒπ‘š2ξ‚πœ‡πΏ=π‘š2ξ‚πΏβˆ’ξ€·π‘”2π‘Œβˆ’π‘”22Δ𝑣28+π‘žπΏξ‚π‘š2𝐷𝐴,π‘š2Μƒπ‘’π‘…β‰ƒπ‘š2ξ‚πœ‡π‘…=π‘š2𝐸𝑐+𝑔2π‘ŒΞ”π‘£24+π‘žπΈπ‘ξ‚π‘š2𝐷𝐴.(7.1) The first terms on the right-hand side π‘š2𝑄,ξ‚π‘ˆπ‘,𝐷𝑐,𝐿,𝐸𝑐 are the corresponding soft masses (3.8), the second terms are the π·π‘Œ,2 contributions with Δ𝑣2=𝑣2π‘’βˆ’π‘£2𝑑=βˆ’π‘£2𝑐2𝛽, while the last terms are the 𝐷𝐴 corrections given by ξ‚π‘š2𝐷𝐴=12ξ€·π‘žπ»π‘’π‘£2𝑒+π‘žπ»π‘‘π‘£2π‘‘βˆ’π‘£2βˆ’ξ€Έ,(7.2) There is an approximated degeneracy between the sfermions with the same charges.

The mass matrix for the third family sfermions is parametrized as β„³2𝑓=βŽ›βŽœβŽœβŽπ‘€ξ‚π‘“2𝐿𝐿𝑀𝑓2𝐿𝑅𝑀𝑓2𝐿𝑅𝑀𝑓2π‘…π‘…βŽžβŽŸβŽŸβŽ ,(7.3) where the off-diagonal terms are generated by F-term corrections proportional to the Yukawa couplings. The stop mass matrix elements are 𝑀̃𝑑2𝐿𝐿=π‘š2𝑑+π‘š2𝑄+ξ‚€13𝑔2π‘Œβˆ’π‘”22Δ𝑣28+π‘žπ‘„ξ‚π‘š2𝐷𝐴,𝑀̃𝑑2𝑅𝑅=π‘š2𝑑+π‘š2ξ‚π‘ˆπ‘βˆ’π‘”2π‘ŒΞ”π‘£26+π‘žπ‘ˆπ‘ξ‚π‘š2𝐷𝐴,𝑀̃𝑑2𝐿𝑅=βˆ’πœ‡π‘šπ‘‘π‘‘π›½βˆ’1.(7.4) The sbottom mass matrix elements are 𝑀̃𝑏2𝐿𝐿=π‘š2𝑏+π‘š2𝑄+ξ‚€13𝑔2π‘Œ+𝑔22Δ𝑣28+π‘žπ‘„ξ‚π‘š2𝐷𝐴,𝑀̃𝑏2𝑅𝑅=π‘š2𝑏+π‘š2𝐷𝑐+𝑔2π‘ŒΞ”π‘£212+π‘žπ·π‘ξ‚π‘š2𝐷𝐴,𝑀̃𝑏2𝐿𝑅=βˆ’πœ‡π‘šπ‘π‘‘π›½.(7.5) The stau mass matrix elements are π‘€Μƒπœ2𝐿𝐿=π‘š2𝜏+π‘š2ξ‚πΏβˆ’ξ€·π‘”2π‘Œβˆ’π‘”22Δ𝑣28+π‘žπΏξ‚π‘š2𝐷𝐴,π‘€Μƒπœ2𝑅𝑅=π‘š2𝜏+π‘š2𝐸𝑐+𝑔2π‘ŒΞ”π‘£24+π‘žπΈπ‘ξ‚π‘š2𝐷𝐴,π‘€Μƒπœ2𝐿𝑅=βˆ’πœ‡π‘šπœπ‘‘π›½.(7.6) The tau sneutrino mass is π‘š2Μƒπœˆπœ=π‘š2ξ‚πΏβˆ’ξ€·π‘”2π‘Œ+𝑔22Δ𝑣28+π‘žπΏξ‚π‘š2𝐷𝐴,(7.7) where π‘šπ‘‘, π‘šπ‘, and π‘šπœ 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 𝐷𝐴 term correction. Furthermore we stress that there is a mass degeneracy between the three sneutrinos Μƒπœˆπ‘’,πœ‡,𝜏 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 anomalous π‘ˆ(1)𝐴 just slightly above the TeV scale, we fix 𝑀𝑉𝐴=10 TeV. The mass scale in the gaugino sector Ξ› is set to be 𝑂(𝑀𝑉𝐴).

8.1. Charge Bounds

The model parameter space can in principle be constrained by precision EW measurements [101]. However, since 𝑀𝑉𝐴=10 TeV, every value of π‘”π΄π‘žπ»π‘’ and π‘”π΄π‘žπ»π‘‘ is allowed by EW precision data if |π‘”π΄π‘žπ»π‘’|,|π‘”π΄π‘žπ»π‘‘|≲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 (π‘”π΄π‘žπ»π‘’, π‘”π΄π‘žπ»π‘‘).

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 Higgs π‘ˆ(1)𝐴 charges corresponds to the yellow spot in Figure 1: 𝑔𝐴=0.1,π‘€π‘‰π΄π‘ž=10TeV,𝐻𝑑1=βˆ’3,π‘žπ»π‘’ξ‚€2=βˆ’5.(8.1) In order to fix the remaining parameters (βŸ¨π›ΌβŸ©, βŸ¨πΉπ‘†βŸ©, πœ†, π‘š, 𝑔𝐡) we assume 𝑣≃246 GeV, and then we choose some benchmark value for 𝑔𝐡 and π‘£βˆ’ in the π‘ˆ(1)𝐡 sector (we recall that 𝑣+=0  (see Section 4)): (A)𝑔𝐡=0.4,π‘£βˆ’=5TeV,(B)𝑔𝐡=0.1,π‘£βˆ’=4TeV.(8.2) The next step is to solve the minima conditions (4.14) determining βŸ¨πΉπ‘†βŸ©, πœ†, π‘š as function of βŸ¨π›ΌβŸ©. In the limit in which 𝑣2β‰ͺπ‘€π‘‰π΄βŸ¨π›ΌβŸ©,𝑣2βˆ’, we get πœ†2≃18𝑒2βŸ¨π›ΌβŸ©π‘”π΄(π‘žπ»π‘‘+π‘žπ»π‘’)/𝑀𝑉𝐴𝑔𝐴𝑔𝐴𝑣2βˆ’βˆ’2βŸ¨π›ΌβŸ©π‘€π‘‰π΄ξ€·π‘žξ€Έξ€·sec(2𝛽)π»π‘‘βˆ’π‘žπ»π‘’ξ€Έ+π‘žπ»π‘‘+π‘žπ»π‘’,ξ€Έξ€»βŸ¨πΉπ‘†βŸ©β‰ƒβˆ’π‘’βŸ¨π›ΌβŸ©π‘”π΄(π‘žπ»π‘‘+π‘žπ»π‘’)/𝑀𝑉𝐴𝑀𝑉𝐴tan(2𝛽)4ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έπœ†ξ€·π‘žπ»π‘‘βˆ’π‘žπ»π‘’ξ€Έξ€·2βŸ¨π›ΌβŸ©π‘€π‘‰π΄βˆ’π‘”π΄π‘£2βˆ’ξ€Έ,|π‘š|2β‰ƒπ‘”π΄βŸ¨π›ΌβŸ©π‘€π‘‰π΄βˆ’12𝑔2𝐴+𝑔2𝐡𝑣2βˆ’ξ€».(8.3) In Appendix B we report the exact formulae. Thus the only remaining free parameters are 𝑑𝛽 and βŸ¨π›ΌβŸ©, and we perform the following analysis of the mass spectrum as a function of 𝑑𝛽 and βŸ¨π›ΌβŸ©. A lower bound on βŸ¨π›ΌβŸ© as a function of π‘£βˆ’ can be obtained, given the approximation (4.16); from (4.19), βŸ¨π›ΌβŸ©β‰ƒ|π‘š|2𝑔+1/22𝐴+𝑔2𝐡𝑣2βˆ’π‘”π΄π‘€π‘‰π΄,(8.4) where we used the relation π‘š2πœ™βˆ’=βˆ’π‘š2πœ‰=βˆ’βŸ¨π›ΌβŸ©π‘”π΄π‘€π‘‰π΄.(8.5) Thus the lower bound on βŸ¨π›ΌβŸ© is obtained simply by setting |π‘š|=0: βŸ¨π›ΌβŸ©π‘šβ‰ƒξ€·π‘”1/22𝐴+𝑔2𝐡𝑣2βˆ’π‘”π΄π‘€π‘‰π΄.(8.6) The condition βŸ¨π›ΌβŸ©>βŸ¨π›ΌβŸ©π‘š must hold since otherwise we would have a massless scalar field in the spectrum (see (4.25)). Another lower bound, βŸ¨π›ΌβŸ©π‘, can be obtained from the condition (4.12), by solving the minima conditions (4.14) and by substituting the corresponding βŸ¨πΉπ‘†βŸ©, πœ† and π‘š values (8.3). The resulting lower bound can be expressed as ξ€ΊβŸ¨π›ΌβŸ©>maxβŸ¨π›ΌβŸ©π‘š,βŸ¨π›ΌβŸ©π‘ξ€».(8.7) No upper bound can be imposed; hence we decide to perform our analysis by considering βŸ¨π›ΌβŸ©β‰²100 TeV.

The parameters πœ† and βŸ¨πΉπ‘†βŸ© are of a particular phenomenological importance since they appear in the πœ‡ and 𝑏 terms (see (4.5) and (4.6)). In the case A, πœ‡ is in the range (900,6000)  GeV and βˆšπ‘ is in the range (50,1200)  GeV while in the case B, πœ‡ is in the range (500,6000)  GeV and βˆšπ‘ 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 a 𝑑𝛽 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 𝑑𝛽 value, we have 𝑀𝑍1𝑀≃10TeV,𝑍2≃2TeV,(8.8) where with 𝑍1 we denote the 𝑉𝐴-like vector.(B)As in the previous case, we just give the 𝑍1,2 masses 𝑀𝑍1𝑀≃10TeV,𝑍2≃400GeV,(8.9) where as in the previous case 𝑍1 is 𝑉𝐴-like.

We will not give the exact values of the 𝑍0 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 𝑍′ 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 𝑑𝛽 for case A (B) and βˆšΞ›/𝑐=(5)10TeV. The red region is the one in which 𝑀2β„Ž0|1-loop∈[124,126]GeV where the β„Ž0 mass is computed considering 1-loop corrections. Since it turns out that the top squarks have small mixing angle and considering the limit 𝑀𝐴0≫𝑀𝑍0, we have [93] 𝑀2β„Ž0|||1-loop≃𝑀2β„Ž0|||tree+34πœ‹2𝑠2𝛽𝑦2π‘‘π‘š2π‘‘ξƒ©π‘šΜƒπ‘‘ln1π‘šΜƒπ‘‘2π‘š2𝑑ξƒͺ≃𝑀2β„Ž0|||tree+32πœ‹2π‘š4𝑑𝑣2ξƒ©π‘šΜƒπ‘‘ln1π‘šΜƒπ‘‘2π‘š2𝑑ξƒͺ,(8.10) where π‘€β„Ž0|tree is the tree-level β„Ž0 mass and we used π‘šπ‘‘=𝑦𝑑𝑣𝑒/2=𝑦𝑑𝑣𝑠𝛽/2. There is an approximated inverse correlation between βŸ¨π›ΌβŸ© and 𝑑𝛽 in the β„Ž0 mass allowed region because the 1-loop correction in (8.10) increases for increasing values of βŸ¨π›ΌβŸ© or 𝑑𝛽. The β„Ž0 mass allowed region is almost the same for case A and B because of two reasons: (i)the mixing with πœ™Β± is suppressed,(ii)the parameters ξ‚π‘šβ„Žπ‘’,ξ‚π‘šβ„Žπ‘‘,πœ‡,𝑏 in the scalar potential (4.22) are ruled by the square mass parameters π‘”π΄βŸ¨π›ΌβŸ©π‘€π‘‰π΄ and (π‘”π΄π‘£βˆ’)2, and the first one turns out to be dominant.

The magenta region satisfies a milder constraint on the light Higgs boson: 𝑀2β„Ž0|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 (Figure 2) 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 βˆšΞ›/𝑐=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, Λ𝑐=5, βŸ¨π›ΌβŸ©=3 TeV and 𝑑𝛽=50 so that π‘€β„Ž0|1-loop≃121.6 GeV; (ii)case A, Λ𝑐=10, βŸ¨π›ΌβŸ©=5 TeV and 𝑑𝛽=10 so that π‘€β„Ž0|1-loop≃124.7 GeV; (iii)case B, Λ𝑐=10, βŸ¨π›ΌβŸ©=4.5 TeV and 𝑑𝛽=50 so that π‘€β„Ž0|1-loop≃125.1 GeV;(iv)case B, Λ𝑐=10, βŸ¨π›ΌβŸ©=50 TeV and 𝑑𝛽=2.5 so that π‘€β„Ž0|1-loop≃130.1 GeV, and we give the full mass spectrum in Figures 4 and 5.

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 πœ†π΅ while in case B is almost a pure πœ†π΅. (ii)An approximated mass degeneracy of 𝐻0, 𝐴0, and 𝐻± holds, and their masses satisfy the bounds of [96, 109]. (iii)The lightest sleptons is a sneutrino, except for 𝑑𝛽=50 when it is Μƒπœ1(iv)The lightest squark is ̃𝑒𝐿, except for 𝑑𝛽=50 when it is ̃𝑏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)𝐢1(2) is close in mass with 𝑁MSSM1(4). 𝐢2 and 𝑁MSSM4 are heavier than all sfermions. (vii)The gluino is close in mass to 𝐢1 and 𝑁MSSM1 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 between ̃𝑒𝑅 and ̃𝑒𝑅 because using the charge constraints (3.16) and (8.1) we get π‘žπΈπ‘=3 and π‘žπ‘ˆπ‘β‰ƒ2.9. (ix)π‘šπœ™βˆ’π‘…<π‘šπœ™+ except for point (i).

Case B points deserve some more comments.

πœ™+ and 𝑁new2,3 are out of the plot of point (iv) because they are heavier than 6 TeV. 𝑍2 is among the lightest not SM particle, so it can decay only into SM particles, because of energy and R-parity conservation. So 𝑍2 is long lived, because SM particles are coupled to 𝑍2 only through the suppressed 𝑉𝐴,𝐡 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 Figure 5 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 𝑆 to the Higgs sector. This term is related to the generation of a suitable πœ‡ and 𝑏 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 Section 8, 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 𝑑𝛽 come from the combined requirement of 𝑀2β„Ž0|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 Figure 5 we explicitly showed two benchmark mass spectrums for our model with βŸ¨π›ΌβŸ© and 𝑑𝛽 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 𝐻0, 𝐴0 and 𝐻±, and ̃𝑒𝑅 and ̃𝑒𝑅 the lightest sleptons is a sneutrino except for 𝑑𝛽=50 when it is stau, the lightest squark is a ̃𝑒𝐿 except for 𝑑𝛽=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].

Appendices

A. Anomalous Lagrangians

The Lagrangian involved in the anomaly cancellation procedure is ℒ𝑆=14𝑆+𝑆†+2𝑀𝑉𝐴𝑉𝐴2|||πœƒ2πœƒ2ξ“βˆ’2ξƒ―ξƒ¬π‘Žπ‘”2π‘Žπ‘π‘Žπ‘Žξ€·π‘Šπ‘†Trπ‘Žπ‘Šπ‘Žξ€Έ+π‘”π‘Œπ‘”π΄π‘π‘Œπ΄π‘†π‘Šπ‘Œπ‘Šπ΄ξƒ­πœƒ2ξƒ°,+h.c.(A.1) where the index π‘Ž=𝐴,𝐡,π‘Œ,2,3 runs over the π‘ˆ(1)𝐴, π‘ˆ(1)𝐡, π‘ˆ(1)π‘Œ, π‘†π‘ˆ(2), and π‘†π‘ˆ(3) gauge groups respectively, and the constants π‘π‘Žπ‘ are fixed by the anomaly cancellation.

Since we have only one anomalous π‘ˆ(1), we can avoid the use of GCS terms, distributing the anomalies only on the π‘ˆ(1)𝐴 vertices. So we have 𝑏𝐴𝐴𝑔=βˆ’π΄π’œπ΄π΄96πœ‹2𝑀𝑉𝐴,π‘π‘Œπ‘Œπ‘”=βˆ’π΄π’œπ‘Œπ‘Œ32πœ‹2𝑀𝑉𝐴,𝑏22𝑔=βˆ’π΄π’œ2216πœ‹2𝑀𝑉𝐴,𝑏33𝑔=βˆ’π΄π’œ3316πœ‹2𝑀𝑉𝐴,π‘π‘Œπ΄π‘”=βˆ’π΄π’œπ‘Œπ΄32πœ‹2𝑀𝑉𝐴,(A.2) where the π’œβ€™s are the corresponding anomalies: π’œπ΄π΄=βˆ’10π‘žπ»π‘‘3βˆ’9π‘žπ»π‘‘2ξ€·π‘žπΏ+3π‘žπ‘„ξ€Έβˆ’9π‘žπ»π‘‘ξ€·π‘žπΏ2+3π‘žπ‘„2ξ€Έβˆ’7π‘žπ»π‘’3βˆ’27π‘žπ»π‘’2π‘žπ‘„βˆ’27π‘žπ»π‘’π‘žπ‘„2+3π‘žπΏ3,π’œπ‘Œπ‘Œ1=βˆ’2ξ€·7π‘žπ»π‘‘+7π‘žπ»π‘’+3π‘žπΏ+9π‘žπ‘„ξ€Έ,π’œ22=12ξ€·π‘žπ»π‘‘+π‘žπ»π‘’+3π‘žπΏ+9π‘žπ‘„ξ€Έ,π’œ333=βˆ’2ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έ,π’œπ‘Œπ΄=5π‘žπ»π‘‘2+6π‘žπ»π‘‘ξ€·π‘žπΏ+π‘žπ‘„ξ€Έβˆ’π‘žπ»π‘’ξ€·5π‘žπ»π‘’+12π‘žπ‘„ξ€Έ,(A.3) where we used the constraints (3.13). Imposing the conditions (3.16) we get π’œπ΄π΄=1ξ€·64βˆ’1168π‘žπ»π‘‘3+1776π‘žπ»π‘‘2π‘žπ»π‘’βˆ’996π‘žπ»π‘‘π‘žπ»π‘’2+53π‘žπ»π‘’3ξ€Έ,π’œ(A.4)π‘Œπ‘Œ=βˆ’114ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έ,π’œ(A.5)221=βˆ’4ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έ,π’œ(A.6)333=βˆ’2ξ€·π‘žπ»π‘‘+π‘žH𝑒,π’œ(A.7)π‘Œπ΄=0.(A.8) We recall that (A.8) is not a consequence of (3.16), but rather (3.16) is a consequence of imposing (A.8) in order to cancel the π‘ˆ(1)π‘Œβˆ’π‘ˆ(1)𝐴 kinetic mixing.

B. Exact Fixed Parameters

In this Appendix section we give the exact values for the βŸ¨πΉπ‘†βŸ©, πœ†, π‘š parameters determined in Section 8.2. Solving the minima conditions (4.14), we get πœ†2=𝑒2βŸ¨π›ΌβŸ©π‘”π΄(π‘žπ»π‘‘+π‘žπ»π‘’)/𝑀𝑉𝐴×32βˆ’π‘”π΄ξ€·π‘žsec(2𝛽)π»π‘‘βˆ’π‘žπ»π‘’ξ€Έξ€·8βŸ¨π›ΌβŸ©π‘€π‘‰π΄+𝑔𝐴𝑣2ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έ(cos(4𝛽)+3)βˆ’4𝑔𝐴𝑣2βˆ’ξ€Έβˆ’8βŸ¨π›ΌβŸ©π‘”π΄π‘€π‘‰π΄ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έβˆ’2𝑣2ξ€·2𝑔2π΄ξ€·π‘žπ»π‘‘2+π‘žπ»π‘’2ξ€Έ+𝑔2π‘Œ+𝑔22ξ€Έ+4𝑔2𝐴𝑣2βˆ’ξ€·π‘žπ»π‘‘+π‘žπ»π‘’,ξ€Έξ€»βŸ¨πΉπ‘†βŸ©=βˆ’π‘’βŸ¨π›ΌβŸ©π‘”π΄(π‘žπ»π‘‘+π‘žπ»π‘’)/𝑀𝑉𝐴×𝑀𝑉𝐴tan(2𝛽)8π‘”π΄ξ€·π‘žπ»π‘‘+π‘žπ»π‘’ξ€Έπœ†Γ—ξ€Ίπ‘”π΄ξ€·π‘žπ»π‘‘βˆ’π‘žπ»π‘’ξ€Έξ€·4βŸ¨π›ΌβŸ©π‘€π‘‰π΄+𝑔𝐴𝑣2ξ€·π‘žπ»d+π‘žπ»π‘’ξ€Έβˆ’2𝑔𝐴𝑣2βˆ’ξ€Έ+𝑣2𝑔cos(2𝛽)2π΄ξ€·π‘žπ»π‘‘βˆ’π‘žπ»π‘’ξ€Έ2+𝑔2π‘Œ+𝑔22,|π‘š|2=π‘”π΄βŸ¨π›ΌβŸ©π‘€π‘‰π΄βˆ’12𝑔2𝐴+𝑔2𝐡𝑣2βˆ’+𝑔2𝐴𝑣2ξ‚€π‘žπ»π‘‘π‘2𝛽+π‘žπ»π‘’π‘ 2𝛽.(B.1)

Acknowledgments

A. Lionetto acknowledges M. Bianchi, E. Kiritsis, and R. Richter for useful discussions and comments. A. Racioppi acknowledges M. Raidal for discussions and the ESF JD164 contract for financial support.