Realization of a spontaneous gauge and supersymmetry breaking vacuum

It is one of the major issues to realize a vacuum which breaks supersymmetry (SUSY) and R-symmetry, in a supersymmetric model. We study the model, where the same sector breaks the gauge symmetry and SUSY. In general, the SUSY breaking model without gauge symmetry has a flat direction at the minimum of F-term scalar potential. When we introduce U(1) gauge symmetry to such a SUSY breaking model, there can appear a runaway direction. Such a runway direction can be lifted by loop effects, and the gauge symmetry breaking and SUSY breaking are realized. The R-symmetry, that is assigned to break SUSY, is also spontaneously broken at the vacuum. This scenario can be extended to non-Abelian gauge theories. We also discuss application to the Pati-Salam model and the SU(5) grand unified theory. We see that non-vanishing gaugino masses are radiatively generated by the R-symmetry breaking and the gauge messenger contribution.


Introduction
It is important to study physics beyond the standard model (SM). Indeed, several types of extensions have been studied. One direction of extensions is to assume larger gauge groups, e.g. U(1) extension and grand unified theories (GUTs) such as SU (5) and SO (10). Another direction is supersymmetric extension such as minimal supersymmetric standard model (MSSM). Supersymmtric gauge-extended models such as supersymmetric GUTs are motivated well by the explanation of the origins of the electroweak (EW) scale and the SM gauge groups. In such models, it is an important key how gauge symmetries and supersymmetry (SUSY) break down. It is also interesting to construct models that both gauge symmetries and SUSY are broken spontaneously by the same sector and their breaking is tightly related with each other [1][2][3][4][5][6][7].
In general, spontaneous SUSY breaking models without gauge symmetries have flat direction at the tree-level potential minimum [8][9][10] like the O'Raifeartaigh model [11]. Such a flat direction could be lifted up by one-loop effects. SUSY breaking models with U(1) gauge symmetry have been studied as well. Then, it is found that the U(1) Dterm potential does not stabilize the flat direction of the F-term scalar potential. However, there can appear a runaway direction along which the D-term potential becomes vanishing, when D-term is non-vanishing at the minimum of F-term scalar potential [12]. Such a runaway direction could be lifted up by one-loop effects and a minimum would appear at non-vanishing finite field value. Then, we could realize both gauge symmetry and SUSY breaking. Note that R-symmetry, that should be broken to realize finite gaugino masses, JHEP11(2017)073 could also be spontaneously broken at the vacuum. Thus, we can evade the vanishing gaugino masses that are often predicted in the gauge mediation models [10,13,14].
In this paper, we study the above scenario, that is, the runaway direction and its lifting to realize both gauge symmetry and SUSY breaking by the same sector. At first we discuss the U(1) model, and then extend it to non-Abelian models. As illustrative models toward realistic GUTs, we discuss the Pati-Salam model [15,16] and the SU(5) GUT [17][18][19][20]. In the Pati-Salam model, the gauge symmetry is SU(4) × SU(2) R × SU(2) L . The gauge symmetry SU(4) × SU(2) R and SUSY are broken at the same time. That can be also realization of gauge messenger models, which can lead a specific spectrum of superpartners [21][22][23]. Similarly, we see that our SUSY breaking scenario can be applied to the SU(5) GUT and compare the result with the one in the Pati-Salam model. This paper is organized as follows. In section 2, we study the SUSY breaking model with U(1) gauge symmetry. We show that there is a runaway direction and it can be lifted by one-loop effects. In section 3, we extend the U(1) model to non-Abelian models, and we apply the above model to the Pati-Salam model and the flipped SU(5) GUT. Section 4 is devoted to conclusion.

SUSY breaking model
In this section, we study flat directions and runaway directions in SUSY breaking models with and without U(1) gauge symmetry. We show that such a runaway direction can be lifted by one-loop effects. Most of the content in this section is review except lifting the runaway direction by one-loop effects.

SUSY breaking models without gauge symmetry
In this section, we review that a generic SUSY breaking model has a flat direction at the potential minimum [8][9][10] like the O'Raifeartaigh model.
We consider renormalizable superpotential W (φ i ) with i = 1, · · · , n. Here, we use the notation that the chiral superfield φ i and its lowest component are written by the same letter. Then, the F-term scalar potential V F is obtained by assuming canonical Kähler potential. Here W i denotes the first derivative of W (φ i ) by φ i , and we use a similar notation for higher derivatives. We assume that the potential minimum is obtained at φ i = φ (0) i and SUSY is broken there. That is, the stationary condition is satisfied as Since we assume that SUSY is broken, some of W i (φ (0) ) must be non-vanishing. Actually, the fermion field along the direction v i = W i (φ (0) ) is massless, and corresponds to the Nambu-Goldstone fermion caused by the SUSY breaking.

JHEP11(2017)073
The mass squared matrix of scalar fields is written by .

(2.3)
Let us evaluate the mass squared along the direction v i , which is the superpartner direction of the Nambu-Goldstone fermion. Its supersymmetric mass is vanishing because of the stationary condition (2.2). Thus, the mass squared along this direction is obtained as If this mass squared is non-vanishing, it can be negative and the vacuum is not stable. For the vacuum to be stable, the above value should vanish, i.e., Now, let us consider the following direction, z: By use of the above results, we find That is, the F-term scalar potential is flat, i + zWī(φ (0) )), along the above direction in eq. (2.6).
Such a flat direction could be lifted by the one-loop effects [24] In the regime that the magnitude of soft SUSY breaking term is much smaller than the field value, the full potential including loop effects could be written by [25] where Z i denotes the wave-function renormalization of φ i .

SUSY breaking model with U(1) gauge symmetry
Here, we study a SUSY breaking model with U(1) gauge symmetry. We show that the flat direction in the previous section is still flat even including the D-term potential, but there can appear a runaway direction, along which certain fields go to infinity, i.e., φ → ∞, and the potential becomes lower [12,26] (see also ref. [27]). The full scalar potential is given by where q i denotes the U(1) charge of φ i . First, we show useful relations among the superpotential, D-term and their derivatives. The superpotential must be invariant under the U(1) gauge transformation, Its derivative by φ j is written as We also obtain In addition, the stationary condition is written as By use of this, we can obtain the following relation: This relation implies that the D-term is non-vanishing only if at least one charged field has non-vanishing F-term. Otherwise, the D-term vanishes. When D is vanishing, the structure of the full potential is the same as one of V F . The mass squared matrix of the scalar fields φ i is written by . (2.17) It is found that when both D = 0 and W i = 0 are satisfied, the direction v i has the same mass squared as eq. (2.4). Using eq. (2.12), we can show that the direction z in eq. (2.6) is flat when the minimum of i ) = 0. Now, let us study the other case that the minimum of Obviously, we find that V (φ (0) ) is larger than V F (φ (0) ) because of the nonvanishing D-term. We examine the value of D along the direction in eq. (2.6):

JHEP11(2017)073
The third term vanishes because of the relation (2.12). We also find that q i |W i (φ (0) )| 2 is vanishing according to the relations (2.14) and (2.2). This result implies that the potential, V F + V D , is flat along the direction defined in eq. (2.6) even adding non-vanishing D-term.
On the other hand, we can show that there is the following runaway direction: where c i is constant. The D-term along this direction is evaluated as We can choose c i such that they satisfy Then, we find that D = 0 in the limit z ∞ → ∞, that is, V → V F (φ (0) ). Hence, there is, in general, a runaway direction in the model, where common fields contribute to both U(1) gauge symmetry and SUSY breaking. We have shown that there is a runaway direction in generic model when the minimum of the F-term scalar potential corresponds to the non-vanishing D-term. Such a runaway direction would be lifted by loop effects. In the next section, we discuss lifting of the runaway direction by using an explicit model.

A concrete model
In this section, we study a concrete model that causes SUSY breaking and predicts a runaway direction at the tree level [12]. Our model includes five chiral superfields, X 0 , X ± , and φ ± . The superfields X + and φ + (X − and φ − ) have U(1) charge, +1 (−1), while X 0 is neutral. We write the superpotential, We also assign R-symmetry to cause SUSY breaking. The fields, X 0 and X ± , have R-charge 2, while φ ± have vanishing R-charge.
The constants, f , λ and m 1,2 , can be defined as positive real values. Assuming that m 1 m 2 < λf is satisfied, the minimum of the F-term scalar potential is given by where F and F 0 are defined as

JHEP11(2017)073
At this minimum, the F-terms are obtained as (2.25) and the F-term scalar potential is written by Furthermore, the minimum of the F-term scalar potential has the following flat direction: At the minimum of the F-term scalar potential, the D-term is evaluated as (2.28) Unless m 1 = m 2 , the D-term D (0) does not vanish. We can confirm that the value of D-term does not change along the direction of eq. (2.27).
Based on the discussion in section 2.2, there is a runaway direction in this kind of model. We investigate the following direction: We choose c ± such that they satisfy that corresponds to the condition (2.21). For large z, the D-term and the D-term potential behave as On the other hand, the F-terms of φ ± behave as Then, the full scalar potential becomes where C is given by (2.34) Thus, this potential has the runaway direction z → ∞. The minimum of C is obtained as .

(2.35)
Now, let us evaluate loop-effects, assuming |λ| 2 ≫ g 2 . We expand φ ± around the minimum, The mass term of δφ ± in the superpotential is written by That is, the non-vanishing X 0 generates the supersymmetric mass of δφ ± . In addition, we have the following term in the scalar potential, that makes the mass splitting between scalars and fermions of δφ ± . Then we obtain the one-loop potential, Note that the full potential can be written approximately [25,28] where γ X 0 is the anomalous dimension of X 0 . At any rate, the above one-loop correction can lift up the runaway direction. The potential for z can be approximated as Then, the stationary condition, ∂V ∂|z| = 0, is satisfied at Thus, by including the one-loop effects, we can obtain the potential minimum with finite vacuum expectation values (VEVs), where both U(1) gauge symmetry and SUSY break down. X 0 and X ± carry the non-vanishing R-charges, so that the R symmetry is also broken at this vacuum. For simple illustrating estimation, we take the parameters such that m 1 ≫ m 2 . Then, we can approximately evaluate the VEVs: and |X (0)

JHEP11(2017)073
Also, the F-terms are approximated by Substituting sample values, let us evaluate the parameters quantitatively. For instance, fixing the parameters at (F, m 2 , g) = (10 × m 2 1 , 0.5 × m 1 , 0.1), we estimate the SUSY and gauge symmetry breaking scales as (2.47) Note that |z (0) | is approximately evaluated as 240/m 1 at this reference point, so that |X It is important to investigate the masses of the fields in the SUSY breaking sectors. At this reference point, the scalar masses squared normalized by m 2 1 are quantitatively estimated as In addition, there is a massive mode from the real part of z, whose mass is given by the one loop correction in eq. (2.42). The imaginary part of z corresponds to the Goldstone boson of the R symmetry. Note that the superpotential in eq. (2.22) leads only SUSY breaking vacua. Adding the D-term, we also find a SUSY breaking vacuum with vanishing X 0 and X ± at the tree level. At this vacuum, φ ± and the F-terms of X ± develop the VEVs, and the SUSY and the gauge symmetry are broken. This vacuum, however, suffers from tachyonic masses of the sfermions, as discussed in section 3.2. The vacuum we have obtained at the one-loop level is located at the point with non-vanishing X 0 and X ± . There, the D-term is suppressed by |z| 2 and the one-loop correction given by the non-vanishing F-terms can easily stabilize the vacuum. The distance between the two SUSY breaking vacua is enough large for our vacuum to be long-lived, because of the runaway behavior. Thus, we focus on this vacuum and construct some models with the GUT gauge symmetries.
Before the application to the GUT models, let us comment on the theoretical aspects of our SUSY breaking model. Above, we have shown that the runaway direction can be lifted up by one-loop effects in one concrete model. The runway behavior is the generic feature in a certain class of SUSY breaking models with gauge symmetries as explained in the previous section. Similarly, runaway directions in generic models could be stabilized by loop-effects in proper parameter regions. It would be important to discuss conditions on lifting of runaway directions in generic models, but it is beyond our scope.

JHEP11(2017)073
Here, we also give a comment on the R-symmetry. The above model has the Rsymmetry, whose charges are assigned such that X 0 and X ± have the R-charge 2 and φ ± have vanishing charge. At the minimum studied above, the fields X 0 and X ± develop VEVs, and then the R-symmetry is spontaneously broken. Note that the U(1) charges of X 0 and X ± are different. For example, if a VEV of a single field breaks the R-symmetry and U(1) symmetry, a new R-symmetry, which is a linear combination of the R-symmetry and broken U(1) symmetry, would remain. However, in the above model, such a new R-symmetry does not remain. Then, the gauge messenger contribution produces non-vanishing gaugino masses at the one-loop level. We see the predictions in some illustrative models.
So far, we have studied the SUSY breaking model with the U(1) gauge symmetry. We can extend this model to the model with non-Abelian gauge symmetry G. In the next section, we apply the above study to models with non-Abelian gauge symmetry, and discuss the applications to the Pati-Salam Model and the SU(5) GUT.

Non-Abelian gauge models
In this section, we extend the previous discussion on U(1) models to non-Abelian gauge models.

SU(N ) model
Here, we consider the extension of the U(1) model to non-Abelian gauge theory. We replace X + and φ + by chiral matter fields with R representation under non-Abelian gauge symmetry (G), and X − and φ − by chiral matter fields with conjugate representation, R, while X 0 is the singlet. For concreteness, we study the model with G = SU(N ) gauge symmetry, where X + and φ + are the N fundamental representations and X − and φ − are its conjugate representations. Then, we consider the same superpotential as eq. (2.22) with the above replacement of representations. Similar to eqs. (2.23) and (2.25), some components in N +N representations develop VEVs and non-vanishing F-terms. By using SU(N ) rotation, we can fix the VEV directions as Thus, the gauge symmetry SU(N ) is broken to SU (N − 1). Similarly, we obtain nonvanishing F-terms along the following directions: The F-term, W X 0 , is the same as eq. (2.25). The D-terms corresponding to the broken generators are non-vanishing at X 0 = X ± = 0, but the tree-level potential has a runaway direction, which is the same as eq. (2.29). Furthermore, similar to eqs. (2.39), (2.40), and (2.42), the potential including one-loop effects would be written as Here, γ X 0 denotes the anomalous dimension of X 0 , which depends on the coupling and multiplicity N . Then, the minimum is estimated as Note that the gauge symmetry breaking scale is given by Compared to the U(1) model, there are extra fields from the decomposition of φ ± and X ± . The VEVs of X 0 , X ± and φ ± can make the remnant fields massive at the tree level, except for z. Then, we obtain the SU(N-1) gauge theory, effectively. Integrating out the remnant fields at the breaking scale, the mass of the SU(N-1) gaugino is radiatively induced. In addition, the mass squared of extra fields charged under SU(N) would be also generated. In order to check the stability of our vacuum, we need estimate the soft SUSY breaking terms. Below, we study the stabilities in some concrete models.
Similarly, we can construct a model, where SU(N ) × U(1) is broken by fields with N q representation and its conjugate where q is U(1) charge. Also we can construct a model, where SU(N ) × SU(M ) gauge symmetries are broken by the fields with (N, M ) representation and its conjugate, respectively. Such models would be interesting from the viewpoint of phenomenological applications: the flipped SU(5) × U(1) model [18][19][20] and SU(4) × SU(2) L × SU(2) R model [15,16] could correspond to the case. In the next section, we discuss the application including quark and lepton chiral superfields and study the soft SUSY breaking terms in each model. Inclusion of squarks and sleptons, however, makes the potential complicated and in general there are directions, where squarks and sleptons develop their VEVs. If all of the squark and slepton masses squared are positive, such a vacuum would be (meta-)stable. We assume that quarks and leptons have no couplings with X 0 , X ± and φ ± . Then, we estimate soft SUSY breaking terms through the gauge mediation. We give some comments on the stability of our vacuum in each setup.

Illustrative models
Based on the above discussion, we construct illustrative models where gauge symmetry and SUSY are simultaneously broken. In the previous section, we introduce the extension

JHEP11(2017)073
to the model with SU(N) gauge symmetry. In the same manner, we can consider a model with G 1 × G 2 gauge symmetries as well. Here, G A (A = 1, 2) is Abelian or non-Abelian gauge symmetry, and both X ± and φ ± are charged under G 1 × G 2 , while X 0 is the singlet.
In our SUSY breaking model, the VEVs of X ± and φ ± break gauge symmetry. If G 1 × G 2 has a bigger rank than the SM gauge symmetry, we could discuss the simple scenario that the SM gauge symmetry is embedded into G 1 and/or G 2 like the GUT and the SUSY breaking sector also causes the GUT breaking. Since the dynamics of SUSY breaking and GUT breaking is explicitly given in this kind of model, the soft SUSY breaking terms for the supersymmetric SM fields are explicitly predicted according to the gauge mediation. Thus, in this subsection, we evaluate soft SUSY breaking terms from the gauge mediation. We neglect D-term contributions in the study below.
Let us assume that one of the SM gauge groups (G SM a ) is given by the part of G 1 × G 2 , the gaugino mass of G SM a is generated at the gauge symmetry breaking scale µ as [29] M a (µ) = α a (µ) 4π ∆b a F X X . (3.8) Here, ∆b a denotes the difference between the beta-function coefficients of G SM a and of G 1 × G 2 . 1 For instance, if the SM SU(3) comes from G 1 , ∆b 3 is given by ∆b 3 1 are the beta-function coefficients of SU(3) and G 1 , respectively. Here, we assume that chiral superfields integrated out at µ obtain the masses from the non-vanishing VEV, X. F X is the F-term of the superfield developing the VEV.
When the MSSM chiral superfield, Q I , is charged under G 1 × G 2 , the non-vanishing A-term and B-term are generated as follows [29]: Here, c A I and c a I are the second Casimir operators of G A and G SM a . The SUSY breaking trilinear coupling corresponding to the Yukawa coupling y IJK , y IJK A IJK Q I Q J Q K , and the SUSY breaking bilinear coupling corresponding to the µ-term, µ H BH u H d , are given by A IJK = A I + A J + A K and B = A Hu + A H d . Note that H u and H d denote the SU(2) L -doublet Higgs fields in the MSSM.
It is a critical feature of this model that non-vanishing A-terms and B-term are generated at the one-loop level. In order to realize 125 GeV Higgs mass, a sizable A-term involving top squark is favorable. Besides, a proper value of the B-term is also necessary to cause the EW symmetry breaking. Then, this feature would be appropriate to construct a realistic supersymmetric model.
Next, we estimate the scalar mass squared in our model. As discussed in ref. [22], there are one-loop corrections to the scalar masses squared in this kind of supersymmetric model. In our model, φ ± and their F-terms also develop non-vanishing VEVs, and the

JHEP11(2017)073
VEVs drive the masses squared negative according to the one-loop level [22]. We estimate the one-loop corrections as where M 2 1 is given by (3.11) M 2 1 is vanishing in the limit that |φ (0) ± | and Fφ ± go to zero. In our model, |φ ± is also suppressed by |z (0) | as shown in eq. (2.47), so that the one-loop corrections are expected to be small. For instance, M 2 1 is estimated as (2.5 × 10 −5 )/|z (0) | 2 , at the reference point in eq. (2.47). On the other hand, the two-loop contributions to the masses squared are estimated as is a subgroup of G 1 . Here, F X /X is dominantly given by X ± and estimated as 1/|z (0) |. Thus, the two-loop contributions could dominate over the one-loop, as far as the gauge couplings are not too small. At the reference point in eq. (2.47), the minimum size of the gauge coupling is about 0.04 for the two-loop correction to be dominant compared to the one-loop. Note that the one-loop contribution is suppressed more significantly, if F is assumed to be much larger than m 2 1,2 . Even if the two-loop contributions dominate the masses squared, the beta-function coefficient of G A may give a negative contribution to the masses squared, as shown in eq. (3.12). In such a case, we would conclude that the vacuum is not stable, when only the gauge mediation is dominant. We need additional contributions to sfermion masses, e.g. gravity mediation, unless large RG corrections are expected. We will give a comment on the extra contributions in section 3.3.
Below, we especially introduce two different models: the Pati-Salam model [15,16] and the SU(5) × U(1) GUT, namely the flipped SU(5) GUT [18][19][20]. In each model, we show the soft SUSY breaking terms and discuss the phenomenological impacts. We also give a short discussion about the conventional SU(5) model [17]. Concerned with the soft SUSY breaking terms, we investigate the one-loop corrections for the gaugino and the A-terms and especially the two-loop corrections for the mass squared. The one-loop corrections may be dominant, depending on the parameters. The one-loop, however, gives negative mass squared, so that we discuss the possibility that the two-loop corrections to the mass squared compensate the tachyonic mass in each model.
In our model, all fields from X ± and φ ± can gain the masses around the SUSY breaking scale. Then, ∆b a are evaluated as follows, assuming that the chiral superfields in the SUSY breaking sector are integrated out at µ: (∆b Y , ∆b 2 , ∆b 3 ) = (−10/3, 0, −1). (3.14) These values lead vanishing wino mass and relatively small gluino mass, according to eq. (3.8). Following eq. (3.9) and eq. (3.12), the A-terms and masses squared are also evaluated. We see that non-vanishing A-terms are generated, if Q I is charged under SU(4) × SU(2) R . In the mass squared, the signs of b ′ A and b a play a crucial role in avoiding the tachyonic masses. In our setup, X ± and φ ± largely contribute to the beta-function coefficients of SU(4) × SU(2) R : b ′ SU(4) = 4 and b ′ SU(2) R = −8. Also, b 3 = −5 is led by this matter content, so that the masses squared of right-handed squarks tend to be negative. The soft-SUSY breaking terms relevant to down-type and up-type squarks are obtained as follows: Here, Q L , u R , and d R denote the SU(2) L -doublet, SU(2) L -singlet up-type, and down-type quark superfields respectively. In these descriptions, the gauge coupling of SU(4) is the same as the one of the SM SU(3). In addition, α R (µ) denotes the gauge coupling of SU(2) R symmetry, and satisfies the following relation at the breaking scale;

JHEP11(2017)073
where α Y denotes the U(1) Y gauge coupling. As we see, the sizable α R (µ) gives the negative contributions to m 2 d R and m 2 u R . Depending on the breaking scale, α R becomes compatible with α 3 and makes m 2 d R and m 2 u R negative. This means that up-type and down-type squarks become tachyonic at the low scale even if the two-loop contributions are dominant, as far as large positive RG corrections are not expected. In this model, the gluino mass is relatively light, so that the RG correction is relatively small.
In the mass squared for right-handed slepton, there is also a negative contribution from SU(2) R : The SU(4) gauge interaction, however, compensates for the negative contribution, so that m 2 e R can become larger than m 2 d R and m 2 u R . Note that the mass squared for left-handed lepton is also positive, because of no SU(2) R contribution.
We conclude that this application of our SUSY breaking scenario to the Pati-Salam model works well to cause both SUSY breaking and GUT breaking. The R-symmetry is spontaneously broken, so that finite gaugino masses are generated by the gauge mediation. This model may, however, suffer from the tachyonic squark masses, if the gauge mediation contribution is dominant in the soft SUSY breaking terms. If the breaking scale is lower than 10 10 GeV, all masses squared can be positive because of small α R . Otherwise, we need other sizable mediation effects such as gravity mediation and anomaly mediation, to lead a realistic supersymmetric SM model. The vanishing wino mass also requires such effects.

Flipped SU(5) GUT
Next, we consider another application of our SUSY breaking scenario to the GUT model: G 1 × G 2 ≡ SU(5) × U(1) X . If U(1) Y is given by the linear combination of U(1) X and the subgroup of SU(5), the GUT model could correspond to the flipped SU(5) GUT [18][19][20]. In the flipped SU(5) GUT, we consider the charge assignment of SU(5) × U(1) X for X ± and φ ± as follows: Higgs fields should be GUT-scale. In our setup, we can write down the following terms: X + X + H, X + φ + H, X − X − H, X − φ − H, and so on. 2 Then, we expect that the colored Higgs fields can obtain the masses around the GUT scale, and mediate the SUSY breaking to the visible sector. Note that we may have to assign

JHEP11(2017)073
R-symmetry to the visible sector and we need some mechanisms to generate the low-scale µ term, that is the supersymmetric mass term of the Higgs doublets. This issue is beyond our scope, and we estimate the soft SUSY breaking terms assuming that the colored Higgs fields are also integrated out at the GUT breaking scale and mediate the SUSY breaking effect.
The threshold corrections, that correspond to the coefficients of the gaugino masses, are given by (∆b 1 , ∆b 2 , ∆b 3 ) = (1, −3, 1). (3.22) Note that ∆b 3 is relatively small, and it is vanishing if the colored Higgs fields do not contribute to the soft SUSY breaking term. In this GUT model, the beta-function coefficients are not so large: b ′ SU(5) = 2 and b ′ U(1) X = −8. The coefficient, b 3 , that appear in the mass squared for squark, is estimated as b 3 = −1. The soft SUSY breaking terms concerned with the squark and slepton masses are estimated as follows:

JHEP11(2017)073
Let us comment on the not-flipped SU(5) GUT case [17]. In this case, the gauge symmetry consists of two symmetry: G 1 × G 2 ≡ SU(5) × U(1) ′ . U(1) Y comes from the subgroup of SU(5) and we could, for instance, consider the following charge assignment for the SUSY breaking sector: This setup, however, leads very large negative b ′ SU (5) , because of many adjoint chiral superfields: b ′ SU(5) = −12. This large value leads Landau pole just above the breaking scale. Besides, we face the big issue concerned with the masses of the colored Higgs fields. In the SU(5) GUT, we introduce two terms, W H = µ H HH + λ Σ ΣHH, where Σ is the adjoint field to break the SU(5) gauge symmetry. We have to allow the fine-tuning between µ H and λ Σ Σ , but in principal we obtain the large hierarchy between the EW Higgs doublet and the colored Higgs fields. Now, we can expect that either X ± or φ ± plays a role of Σ and realizes the hierarchy. The U(1) ′ symmetry, however, forbids either µ H or λ Σ , so that it is impossible to gain the hierarchy in this setup. Besides, the VEVs of X ± are expected to be large, and then X ± should be identical to Σ in W H . This setup, however, causes the bilinear term of the scalar components of H and H, according to the non-vanishing F-terms of X ± . Therefore, it is difficult to realize the realistic EW symmetry breaking vacuum.

Tachyonic mass
We have studied two examples towards constructing realistic models. Indeed, by the mechanism in section 2, we can break the gauge symmetry and SUSY in realistic GUT gauge theories. However, only pure gauge mediation may lead to tachyonic squark and/or slepton masses, especially in the Pati-Salam model. That implies that such vacua are not stable or even meta-stable. In order to stabilize the vacuum, we need another contribution, e.g. gravity mediation, in such a case. For example, we can assume the additional term in Kähler potential, where Q I denotes quark and lepton superfields, such that squarks and sleptons have positive masses squared. Phenomenological aspects of models depend strongly on c 0 and c ± . On the other hand, if we have no additional corrections on the gaugino masses and A-terms except the pure gauge mediation, these can be predictions of our models. Alternatively, we may assume that the anomaly mediation [30,31] is comparable with the gauge mediation discussed above. The pure anomaly mediation leads to tachyonic slepton masses, although squark masses squared are positive.
In the Pati-Salam model, the vanishing wino mass also requires such additional contributions. A proper combination of the gauge mediation, the gravity mediation and anomaly mediation would lead realistic mass spectrum of the SUSY particles in a certain GUT breaking model. Such a study is challenging and we would study it elsewhere.

JHEP11(2017)073 4 Conclusion
It is one of important issues to understand the vacuum structure of our universe. If SUSY really exists in our nature, our vacuum spontaneously breaks the symmetry, so that it is a major issue to construct a SUSY breaking model.
When SUSY breaking is triggered by F-terms of chiral superfields, it is known that the symmetry breaking is accompanied by flat directions in the field space. The flat directions should be stabilized at the non-vanishing VEV to realize the R-symmetry breaking. Besides, it is a big issue to induce non-vanishing gaugino masses in the gauge-mediation models, even if the R-symmetry is broken at our vacuum. Thus, it is not trivial to find the realistic SUSY breaking vacua and construct the SUSY model that predicts massive superpartners of the SM particles.
In this paper, we consider a supersymmetric model with U(1) gauge symmetry and R-symmetry. In this model, both of the gauge symmetry and SUSY are broken by the same fields. We find flat directions triggered by the SUSY breaking, and the D-term of the U(1) gauge symmetry is not vanishing along the flat directions. In this kind of model, it is known that there are also runaway directions at the tree-level [12]. We suggest that such runaway directions can be lifted by the one-loop effect, and the SUSY breaking vacuum can be realized. The gauge symmetry breaking is also caused by the SUSY breaking dynamics, and the R-symmetry also spontaneously breaks down. In such a case, the gauge messenger field can mediate the SUSY breaking effect and can induce non-vanishing gaugino masses.
We can extend this U(1) model to non-Abelian theory. It is quite interesting to apply this mechanism to the GUTs, e.g. the Pati-Salam model and the flipped SU(5). This simple setup may, however, cause the problem that squarks and sleptons develop VEVs according to the one-loop and two-loop corrections. We estimate the soft SUSY breaking terms concerned with sfermions through the gauge mediation. In the Pati-Salam model, the SU(2) R contributions to the mass squared are negative even at the two-loop level, so that especially the squark masses become tachyonic depending on the size of gauge coupling, i.e. the breaking scale. On the other hand, we find that all masses squared can be positive in the flipped SU(5), taking into account the two-loop corrections. We need study in more detail, taking into account how to realize the realistic Yukawa couplings in the MSSM. In the case that the negative mass squared is derived, we propose another contribution, e.g. gravity mediation and anomaly mediation. In particular, such additional contributions are required by the vanishing wino mass in the Pati-Salam model. Those contributions may drastically change the mass spectrum, and phenomenology may depend on details of mediations. Further study on the GUT with our SUSY breaking model will be given in the near future.