Low-Energy Probes of No-Scale SU(5) Super-GUTs

We explore the possible values of the $\mu \to e \gamma$ branching ratio, $\text{BR}(\mu\rightarrow e\gamma)$, and the electron dipole moment (eEDM), $d_e$, in no-scale SU(5) super-GUT models with the boundary conditions that soft supersymmetry-breaking matter scalar masses vanish at some high input scale, $M_{\rm in}$, above the GUT scale, $M_{\rm GUT}$. We take into account the constraints from the cosmological cold dark matter density, $\Omega_{CDM} h^2$, the Higgs mass, $M_h$, and the experimental lower limit on the lifetime for $p \to K^+ \bar \nu$, the dominant proton decay mode in these super-GUT models. Reconciling this limit with $\Omega_{CDM} h^2$ and $M_h$ requires the Higgs field responsible for the charge-2/3 quark masses to be twisted, and possibly also that responsible for the charge-1/3 and charged-lepton masses, with model-dependent soft supersymmetry-breaking masses. We consider six possible models for the super-GUT initial conditions, and two possible choices for quark flavor mixing, contrasting their predictions for proton decay with versions of the models in which mixing effects are neglected. We find that $\tau\left(p\rightarrow K^+ \bar{\nu}\right)$ may be accessible to the upcoming Hyper-Kamiokande experiment, whereas all the models predict $\text{BR}(\mu\rightarrow e\gamma)$ and $d_e$ below the current and prospective future experimental sensitivities or both flavor choices, when the dark matter density, Higgs mass and current proton decay constraints are taken into account. However, there are limited regions with one of the flavor choices in two of the models where $\mu \to e$ conversion on a heavy nucleus may be observable in the future. Our results indicate that there is no supersymmetric flavor problem in the class of no-scale models we consider.


Introduction
Supersymmetry remains an attractive prospective extension of the Standard Model (SM), despite its non-appearance during Runs 1 and 2 of the LHC [1,2]. Indeed, the discovery of a 125-GeV Higgs boson at the LHC [3] has supplemented the traditional arguments for supersymmetry, which include the naturalness of the electroweak scale [4], the unification of the fundamental interactions [5] and the existence of a cold dark matter candidate (if R-parity is conserved) [6]. The minimal supersymmetric extension of the SM (MSSM) predicted the existence of a Higgs boson with mass M h 130 GeV [7], and is a prime example of new physics capable of stabilizing the electroweak vacuum for M h ∼ 125 GeV [8]. Furthermore, global fits in the framework of simple supersymmetric models suggest that the couplings of the lightest supersymmetric Higgs boson should be very similar to those of the Higgs boson in the SM, as is indicated by the ATLAS and CMS experiments [9,10]. When the supersymmetric particle masses are large, which is the case we consider, the Higgs couplings resemble even more closely the couplings predicted by the SM.
However, the continuing absence of supersymmetry at the LHC [1,2] reinforces the need to seek complementary indications of supersymmetry outside colliders. It is in this context that we address the questions of proton decay, contributions to the electron dipole moment and µ flavor violation observables in the SU(5) models based on no-scale supergravity that were introduced in [11]. There, Higgs fields were assigned to twisted chiral supermultiplets with a suitable choices of modular weights in order to obtain the correct mass of the observed Higgs boson and the cold dark matter density, while avoiding proton decay in violation of the current limits.
Contrary to what happens in the Standard Model, where flavor and CP violation are controlled by the Cabibbo-Kobayashi-Maskawa (CKM) matrix, there is no established mechanism for flavor and CP violation in supersymmetry, the so-called supersymmetric flavor problem. Experiments show that many low-energy predictions of CKM mixing must be reproduced in any extension of the SM, which is therefore an important constraint on any supersymmetric model that is studied.
In a previous study of super-GUT no-scale models in [12] we adopted a pragmatic approach to this challenge, using particular Ansätze for Yukawa couplings to study flavor violation constraints in a scenario with maximal sfermion flavor violation at the input scale M in > M GUT . Here we revisit flavor violation and proton decay, considering alternative options for the flavor mixing associated with different embeddings of the MSSM fields in GUT multiplets. 1 In the SM, Yukawa couplings in the up-and down-quark sectors are described by a couple of 3 × 3 complex matrices whose diagonalizations each require two unitary matrices, one acting on left-handed quarks and the other on right-handed quarks. The two left-handed matrices, one in the up-quark sector and the other in the down-quark sector, combine to form the CKM matrix, whereas the right-handed matrices remain unobservable. In supersymmetry, however, the right-handed matrices propagate into the soft-breaking terms and hence become constrained by flavor observables. These observables clearly indicate that off-diagonal elements of the right-handed sfermion mixing matrices should be tiny 2 . Any model of supersymmetric flavor must specify how to reproduce the CKM matrix via the two down-and up-quark left-handed matrices that diagonalize the Yukawa couplings.
One choice is to associate the CKM matrix with the up-quark Yukawa matrix, for which electroweak (EW) precision observables play an important role in constraining how this is propagated into the supersymmetric sector, as was studied for the CMSSM in [15]. Another is to associate the CKM matrix with the down-quark sector, as we considered in [11]. In this case the constraints from flavor observables are more stringent than those from EW observables, particularly for the low tan β values that we use.
We study in this paper six different no-scale super-GUT SU(5) models, some with both electroweak Higgs representations in twisted chiral supermultiplets, and some with only one twisted Higgs supermultiplet. The soft supersymmetry-breaking masses of the MSSM matter sfermions vanish at the input scale M in in all the models, but they have different boundary conditions for other supersymmetry-breaking parameters. Four of the models have M in = 10 16.5 GeV, whereas the other two have M in = 10 18 GeV, in which case there are larger renormalization-group running effects above the GUT scale, M GUT . For each model, we study predictions for proton decay, µ → eγ and the electron EDM, using two possible choices for the flavor embeddings of the quarks and leptons into SU(5) multiplets that illustrate the ambiguity discussed in the previous paragraph. We find that proton decay rates are relatively insensitive to the treatment of flavor mixing, whereas µ → eγ and the electron EDM are more sensitive. In general, the predictions for these flavor observables are below the present experimental limits when the cosmological dark matter density and the proton lifetime are taken into account, though there are limited regions with one of the flavor choices in two of the models where µ → e conversion on a heavy nucleus may be observable in the future. These no-scale super-GUT models have no supersymmetric flavor problem, as also argued in [12].
This paper is organized as follows. In Section 2 we introduce the class of no-scale SU(5) super-GUT models we study, including the specification of different choices for the embedding of MSSM fields in GUT multiplets and the corresponding Ansätze for matter Yukawa coupling matrices, the no-scale boundary conditions on soft supersymmetry breaking at M in , and our treatment of the renormalization-group running down to the electroweak scale. Then in Section 3 we discuss how proton decay, µ → e flavor-violating observables and the electron EDM arise in these models, and review the available experimental information. In Section 4 we introduce the specific no-scale models we study, and analyze their predictions for these observables. We then present our conclusions in Section 5.
2 Model Framework 2.1 Embedding the MSSM in SU (5) In the minimal supersymmetric SU(5) GUT model, the three generations of matter superfields are embedded into three pairs of 5 and 10 representations. There are also two chiral electroweak Higgs superfields H u and H d , whose vacuum expectation values (vevs) break the electroweak SU(2)×U(1) gauge group down spontaneously to U(1) EM . They are embedded in 5 and 5 representations, H and H, which also contain 3 and 3 colored Higgs superfields H C and H C , respectively. The SU(5) GUT gauge group is broken spontaneously down to the Standard Model (SM) gauge group by the vev of a 24 chiral superfield, Σ ≡ √ 2 Σ A T A , where T A (A = 1, . . . , 24) are the generators of SU (5) with Tr(T A T B ) = δ AB /2. The vev of the adjoint is given by Σ = V · diag(2, 2, 2, −3, −3), with V = 4µ Σ /λ ′ . We follow the notation of [16,12,17,18] for the SU(5) superpotential parameters: where we have suppressed all SU (5) indices. Once SU(5) is broken, the GUT gauge bosons acquire masses M X = 5g 5 V , where g 5 is the SU(5) gauge coupling. Doublet-triplet separation within the H and H representations can be achieved by a fine-tuning condition: µ H − 3λV ≪ V , in which case the color-triplet Higgs states have masses M H C = 5λV . We note also that the masses of the color and weak adjoint components of Σ are equal to M Σ = 5λ ′ V /2, while the singlet component of Σ acquires a mass M Σ 24 = λ ′ V /2.
Our notation for the Yukawa couplings of MSSM fields is specified by the following low-energy superpotential: Note that we use a "Left-Right" (LR) notation for Yukawa couplings, which means that the first index of the Yukawa couplings corresponds to the SU(2) doublets, and the second index to the SU(2) singlets. In order to match the GUT theory (1) to the MSSM (2), in particular for the proton decay operators we discuss below, we decompose the second row of the SU(5) superpotential (1) into MSSM component fields, yielding the Yukawa couplings of the MSSM fields in terms of the SU(5) field couplings, as follows: where the superscripts C on Higgs multiplets indicate their color triplet components. We recall that the embedding of the MSSM fields into the SU(5) model is ambiguous, and various Ansätze are possible. In particular, the following SU(5) Yukawa couplings were chosen in [18] where V GCKM is the CKM matrix at the GUT scale. Transforming the fields where the phase factors φ i satisfy the condition 3 Throughout this work,ĥ denotes a diagonalized Yukawa matrix.
so that only two of them are independent 4 . It is well known that the masses of the leptons and down-type quarks of the first two generations are not consistent with unification at the GUT scale, 5 whereas those of the third generation are in reasonable agreement with Yukawa unification. We determine the SU(5) Yukawa couplings by using the following matching conditions for the MSSM couplings after renormalization group (RG) running them from the electroweak scale up to the GUT scale: Thus, the Yukawa couplings of the charge-2/3 quarks are matched directly to the GUT-scale couplings of the 10 representations, up to a numerical factor, as are those of the first two generations of quarks in the5 representations. 6 Recalling that the third-generation Yukawa couplings for b and τ are similar, we match an average of these Yukawa couplings to that of the third generation of5 fermions. Using as input the values for the Yukawa couplings at the EW scale discussed further below, we use Eq. (7) to determine the SU(5) Yukawa couplings at the GUT scale, which we then run up to M in . Note that we also run the Yukawa couplings of the first two generations of charged leptons up to the GUT scale. These are not used as a basis for further running to M in , but are subsequently run back down to the EW scale.
There are ambiguities in the description of flavor mixing in the supersymmetric GUT model. Various options were considered in [12], including the contrasting cases V R = 1 and V R = V GCKM . If we choose V R = 1 in (4), we obtain from Eq. (7) and the embedding (5) the following relations between the MSSM couplings and the diagonal GUT-scale couplings (after running down from M in to M GUT ): Because of the lack of Yukawa coupling unification, we do not relate h E and h D (3,3) to h5 at the GUT scale. We also do not relate h U (3,3) to h 10 at the GUT scale, in order to converge more efficiently to the observed top quark mass. For h E , h D (3,3) and h U (3,3) , the previous values at M GUT are used for running back down to the EW scale, as will become clear when we discuss the RGE boundary conditions below. This is one of three choices for the treatment of flavor that we consider in this paper: • We call choice A the embedding (5) combined with V R = 1 in (4). This is the Ansatz A2 considered in [12].
We consider also the embedding (after shifting only U c i → e −φ i U c i ), Choosing again V R = 1, we obtain once again Eq. (8) for matching when running down from M GUT to the EW scale.
• We call this choice of embedding B, noting that it is equivalent to Ansatz A3 of [12]. 7 At this point A and B are identical. There would be no difference if we had Yukawa i.e., equal to the diagonal SU(5) coupling as in case A. However, since we do not match h E from the 5-plet, we can only"mimic" this condition at the EW scale and, as we see below, the boundary conditions for A and B differ at the EW scale. We emphasize that in the case of perfect unification the choices A and B would make identical predictions for all observables. A and B would not be distinct cases but rather different ways of formulating the same model for specifying the lepton sector in terms of the 5-plet of SU(5) and possibly additional operators. The motivation to consider cases A and B here is to explore the sensitivity to the precise way the couplings in the charged-lepton sector alter flavor observables.
• We also compare our results for τ (p → K +ν ) with these flavor choices to models that ignore the flavor structure by limiting the RG running to diagonal matrix elements. We label this choice NF.
The Yukawa couplings of the MSSM fields entering the dimension-six operators mediating proton decay can be defined from the Yukawa couplings of the SU(5) theory, Eq. (3), as follows: where (V 10 ) kl =V GCKM for A and 1 for B, while V R = 1 for both of the choices A and B.

Soft Supersymmetry Breaking
We write the soft supersymmetry-breaking terms in the Lagrangian in the SU(5) GUT symmetry limit as − a 10 10 10 H + a5 10 5 H + h.c. , 7 If we take V R = V GCKM with this embedding, we obtain Ansatz A4 of [12]. This choice turns out to be problematic for the observables we discuss below, and is not considered further here.
where theλ A are the SU(5) gaugino fields. For convenience, we make no distinction in notation between chiral superfields and their scalar components.
In super-GUT models [20,21,16,17,11,18], the soft supersymmetry-breaking mass parameters are taken to be universal at some input scale, M in , that is greater than the GUT scale, M GUT . The RG running of the couplings and masses then takes place in two stages. We run the 2-loop MSSM beta functions for Yukawa couplings, trilinear terms, soft masses-squared, m H d , m Hu , B, and µ between the electroweak scale, M EW , and M GUT , including three generations of fermions and sfermions, the SU(3)×SU(2)×U(1) gauge bosons and gauginos, and the SU(2)-doublet Higgs bosons and Higgsinos. Then, between M GUT and M in the SU(5) GUT parameters are run also with three generations of fermions and sfermions, SU(5) gauge bosons and gauginos, Higgses and Higgsinos. For the sake of clarity we now specify all the boundary conditions we impose at M in and M GUT .
Our boundary conditions at M in are derived from no-scale supergravity [22,23,24]. We assume a Kähler potential of the form where T is a volume modulus, the φ i are untwisted matter fields and include the SU(5) matter multiplets. The ϕ a are twisted fields, which include H and/orH, and the n a are the modular weights of the twisted fields. We also allow for modular weights in the superpotential, writing where c is an arbitrary constant, and W 2,3 denote bilinear and trilinear terms with modular weights β, α that are in general non-zero and can differ for each superpotential term. When φ, ϕ = 0, the effective potential for T is completely flat at the tree level, with an undetermined vev, and the gravitino mass is undetermined, varying with the value of this volume modulus 8 . We assume here that some Planck-scale dynamics fixes T =T = c, and assume the representative value c = 1/2 in the following. 9 Finally, we assume a universal gauge kinetic function f ab = δ ab , so that at M in there is a universal gaugino mass, m 1/2 . We work with the no-scale framework introduced in [17], where m 0 = 0, but allow for the possibility that the Higgs 5-plets are twisted, in which case either one or both of their soft masses may be non-zero. It was shown in [11] that in models in which matter and both Higgs supermultiplets are untwisted, the minimal SU(5) super-GUT model considered here is unable to provide simultaneously a dark matter relic density and Higgs mass in agreement with experimental values, and at the same time provide a sufficiently long proton lifetime. It was concluded in [11] that either one or both of the Higgs multiplets must be twisted. The bilinear and trilinear soft supersymmetry-breaking terms b Σ , b H , a ′ , a, a 10 , a 5 may also be non-zero. Each gets a contribution from the modular weight in Eq. (13) and an additional contribution that depends on the specific superpotential term and whether the 5-plets are twisted or not. Our boundary conditions at M in are therefore: The parameters p, q = (0, 1) depend whether (H,H) is untwisted (0) or twisted (1). The parameters r F = p, q, p + q, 0, for F = 10, 5, λ, λ ′ , and p S = p + q, 0 for S = H, Σ. The different modular weights, α, β, chosen for the different models are specified in Section 4.1. We take all the n a = 0. Other quantities run up to M in , such as the SU(5) Yukawa couplings, are not reset at M in .

Renormalization-Group Running of Parameters
Having specified the theoretical boundary conditions at M in , we now discuss the renormalization-group (RG) running of the model parameters. This involves matching parameters at M GUT , since the fundamental degrees of freedom and hence the RG equations differ above and below this scale, and the phenomenological inputs for the gauge and Yukawa couplings are measured at the electroweak scale. The RG equations are run up and down between the electroweak scale and M in iteratively until a convergent solution is found. We use the following matching and boundary conditions.
Matching boundary conditions at M GUT : There are two sets of boundary conditions at M GUT , one corresponding to RG running from the EW scale to M in , and the other when running back down. We first specify the matching conditions for the gauge couplings when running up from the EW scale. At one-loop level in the DR renormalization scheme [25], we have where g 1 , g 2 , and g 3 , are the U(1), SU(2), and SU(3) gauge couplings, respectively, and Q is a renormalization scale taken in our analysis to be the unification scale: Q = M GUT .
The last terms in Eqs. (16) -(18) represent a possible contribution from the dimensionfive operator W ∆g eff = where W ≡ W A T A denotes the superfields corresponding to the field strengths of the SU(5) gauge vector bosons V ≡ V A T A . Since V /M P ≃ 10 −2 , these terms can be comparable to the one-loop threshold corrections, and their possible presence should be taken into account when discussing gauge-coupling unification [26]. Including the c 5 coupling is essential for our purposes, as it allows us to choose independently the Higgs couplings λ and λ ′ , which we specify at the GUT scale. Eqs. (16 -18) can be combined to give The masses, M H C and M Σ , have implicit dependences on the gauge couplings, including g 5 , making it impossible to write an analytic expression for the matching of the three low-energy gauge couplings, g i , to g 5 . Nevertheless, we can solve for g 5 iteratively.
The matching conditions for the Yukawa couplings were given in Eq. (7). As noted there, we take the average of h E 3,3 and h D3,3 for the third-generation charged-lepton and charge-1/3 quark Yukawa couplings, which are close to the unification expected in SU(5). We adopt a similar approach for the trilinear terms and the soft squared masses. For the embedding A, when matching from M GUT to M in we take for the trilinear couplings and for the soft squared masses For the embedding B, when matching from M GUT to M in we take the same matching conditions for a 10 and m 2 5 as for the embedding A, see Eqs. (22,23), respectively, with and We note that by taking these averages we are effectively generating two inequivalent models at the GUT scale, which in turn produce different values for observable quantities. If one was not required to use the averages in Eq. (7) and Eqs. (21,24,25,26), perfect unification would allow us simply to formulate the SU(5) theory with the quark-sector couplings, and all quark-sector differences between the two models would vanish.
The remaining matching conditions for masses at M GUT when running up from the EW scale are: The matching of the gaugino masses to M 5 when running up to M GUT is chosen to be consistent with the matching of the gaugino masses to M 5 when running down from M GUT to the EW scale as discussed below. Finally, when running from M GUT to M in , m 2 Σ is set equal to its value from the previous iterative run down from M in where it was initially set to 0 as in Eq. (15).
At M in , the soft mass terms are reset according to Eq. (15) and the theory is run down to M GUT , where the matching conditions for the soft squared-mass terms and Yukawa couplings are For the trilinear terms, we use which can be obtained from Eq. (30) by setting g 1 (M GUT ) = g 2 (M GUT ). It is important to note that allowing c 5 = 0 enables us to increase the colored Higgs mass, thereby increasing the proton lifetime [11,18]. The matching conditions for the gaugino masses [26,30,31,18] are Finally, we must match the MSSM µ and B-terms to their SU(5) counterparts [32] with As noted earlier, in the minimal SU(5) GUT model studied here we must tune |µ H − 3λV | to be O(M SUSY ). The parameters µ and B can be determined at the electroweak scale by the minimization of the Higgs potential as in the CMSSM. These are then run up to the scale where Eqs. (37) and (38) are applied. However, the GUT A-and B-terms are specified at the input scale by Eq. (15) and, in general, the condition (38) will not be satisfied. This mismatch can be rectified by adding a Giudice-Masiero (GM) term to the Kähler potential [33]: where we have allowed for the possibility of additional modular weights, γ H and γ Σ . This term induces shifts in both the µ-terms and B-terms [34,11,18]: As a result, there is a shift in ∆ given by Then any mismatch in (38) can be corrected by where we have used µ Σ = λ ′ V /4 and µ H = 3λV . If λ ≫ λ ′ , we can ignore, c H , and use (44) to determine c Σ (for a given value of γ Σ ).
Boundary conditions at M EW : Although the soft supersymmetry-breaking parameters are input at the high scale, M in , some of the phenomenological inputs are set by boundary conditions at the electroweak scale, M EW , namely the ratio of electroweak Higgs vevs, tan β, m f and V CKM . The Higgs vevs are in principle determined by the minimization of the Higgs potential at the weak scale. However, it is common in constrained models to fix these by using the experimental value of M Z and tan β, and solve for µ and the pseudoscalar Higgs mass, or equivalently the MSSM B-term. In very constrained models such as the no-scale models considered here, B is fixed by the high-scale boundary conditions and as a consequence, either tan β is an output rather than an input [35], or a GM term is used to fix the matching conditions for the B-terms. We adopt the latter approach here, and treat tan β as a weak-scale input.
We also use the experimental values of the masses of the six quarks and the three charged leptons, m f . The matching of Yukawa couplings is done in terms of the CKM matrix elements, using experimental input for the CKM matrix at M EW . In general h D and h E can be written as follows where y µ , y τ ) are the diagonalized mass matrices containing the mass eigenvalues for the D-type quarks and charged leptons, respectively. The U matrices aid with the diagonalization of these matrices. When running up to the M GUT scale they should match Eqs. (8) at M GUT for the choices A and B, respectively. Hence, in both cases we start with At M GUT the RG evolution determines the evolution of V CKM into V GCKM , while U D R , U E L and U E R are no longer diagonal. However, since we match the SU(5) fields to the MSSM fields at M GUT with Eq. (7), once the RG program has converged, U R D is in practice equal to 1. We match the Yukawa couplings for the first two generations of charged leptons at M GUT , so that they converge rapidly to satisfy U E R = U E L = 1 for A and U E R = V * CKM and U E L = 1 for B at the EW scale. Any remaining non-diagonality can be absorbed into the embedding of the MSSM fields into SU(5), and does not alter the Yukawa couplings relevant for proton decay. Finally, all the fermion masses are converted appropriately to the DR scheme and then matched to the supersymmetric theory at M Z .

Proton Decay
The most important constraint on the supersymmetric SU(5) GUT model from searches for proton decay comes from the decay mode p → K +ν , for which the current experimental limit is [36] τ p → K +ν > 6.6 × 10 33 yrs .
10 In the way we define the Yukawa couplings, these enter the SM interaction Lagrangian as In this paper we will refer to this limit as the proton life-time limit if not otherwise specified. In the future, the Hyper-Kamiokande (HK) experiment is expected to be sensitive to τ (p → K +ν ) ∼ 5 × 10 34 yrs [37], an improvement by almost an order of magnitude. Since generic amplitudes for dimension-5 proton decay are inversely proportional to sparticle masses (see below), the HK reach for proton decay will provide sensitivity to supersymmetric model parameters ∼ 3 times larger than the current constraints from τ (p → K +ν ).

Dimension-5 Proton Decay Operators
In [38] a complete analysis of proton decay operators in supersymmetric SU(5) theories was given, including in particular the explicit forms of the Wilson Coefficients (WCs) C 5L and C 5R entering into the dimension-five Lagrangian generated by integrating out the colored Higgs multiplets [39]: where i, j, k and ℓ are flavor indices, and where a, b, c are color indices. Normalizing these operators at the GUT scale, M GUT , and matching the Yukawa matrices using Eq. (7), we find The Yukawa matrices appearing in Eq. (49) are different for the different embeddings, as seen in Eq. (10). This is because each of the terms in the superpotential Eq. (3) that are relevant for proton decay depend on V R and the choices of the h 10 and h 5 Yukawa matrices in Eqs. (4). The leading-order RG evolutions of the C ijkl 5L and C ijkl 5R between M GUT and the supersymmetry breaking scale are given by [38] where Λ is the renormalization scale. Below the supersymmetry-breaking scale, we use the RGEs given in Ref. [40]. We write the effective Lagrangian for p → K +ν i decay in the following form: The operators C LL (usdν k ) and C LL (udsν k ) are mediated by Wino exchange, and C RL (usdν τ ) and C RL (udsν τ ) are mediated by higgsino exchange (see Eqs. (23) and (27) of [40]). At the EW scale, the operators entering into the proton decay amplitudes are C 221i 5L and C 331i 5L , i = 1, 2, 3, which contribute to C LL (usdν k ) and C LL (udsν k ), and C * 3311 5R and C * 3312 5R , which contribute to C RL (usdν τ ) and C RL (udsν τ ).
However, due to the off-diagonal nature of the Yukawa matrices, the evolution from M GUT down to M EW induces contributions from some other operators. Consider as an example C 3312 5L , whose leading-order RG terms are The terms in Eq. (53) involving h D and h E are not diagonal, and generate contributions to the β functions of the operators mentioned above. In particular We are therefore required to run C 3332 5L between the weak and GUT scales, using the initial condition set by Eq. (10), even though the corresponding operator does not contribute directly to the effective Lagrangian (52) defined at the EW scale. We note, on the other hand, that the combinations h † U h U and h U h † U appearing in Eq. (50) for C 5L and Eq. (51) for C 5R , respectively, remain diagonal as in the case considered in [40] (see Eq. (22) of that reference). 11 The dimension-6 operator coefficients C LL (ud p d q ν k ) and C RL (ud p d q ν τ ), p, q = 1, 2, are related to the dimension-five WCs C 221i 5L , C 331i 5L , C * 3311

5R
, and C * 3312 5R (which were obtained by integrating out the colored Higgs multiplets in Eq. (49)) via CKM mixing angle factors and loop integrals: where 11 We have omitted contributions that are proportional y u and y c , given their smallness in comparison to y t , and we have omitted terms proportional to C 1312 5L .
Here m t R , m τ R , m Q j , and m L k are the masses of the right-handed stop, the right-handed stau, left-handed squarks, and left-handed sleptons, respectively, α i ≡ g 2 i /(4π), and The loop integrals (56) yield the dimension-6 operator coefficients at the supersymmetrybreaking scale, and they must then be run down to the EW scale. The corresponding RGEs are given in [40], where many other details of the calculation are provided. Finally, as also given in [40], the partial decay width for p → K +ν i decay is where The proton decay rates (58) depend on the Yukawa coupling matrices through the various WCs, and hence on our choice of diagonalization scheme. As an illustration of this sensitivity, in Fig. 1 ) is about 10% larger for choice NF (i.e., when off-diagonal flavor-violating effects are ignored) mainly because of the treatment of the Yukawa couplings. When off-diagonal terms are considered in the Yukawa couplings, off-diagonal terms appear also in the soft masses-squared, trilinear terms, etc., which affect the running of the gauge couplings, with the largest effect being that on g 2 3 . Note that the off-diagonal terms in h d affect not only y d and y s as shown in Fig. 2 (see below) but also the CKM matrix elements at the GUT scale. At the electroweak scale, the difference in the WCs is about the same, (roughly 10% between choices A/B and NF) though the magnitudes of the coefficients are about 3-4 times larger.
The sensitivities to the Yukawa couplings of the charge-1/3 quarks d, s and b are also significant. We can understand this by considering the one-loop β function of h d , which is given by where  to evolving only the diagonal elements, we see from the hierarchy of the elements of the Yukawa couplings that the lightest eigenvalue, corresponding to y d , will be affected mainly by |h 11 We see that, due to the term proportional to y 2 t in β , which in turn produces a higher value of y 2 s at M GUT than when the running of off-diagonal Yukawa couplings is neglected, because no information on the evolution of h 23 d is considered in that case. A comparison of the squared Yukawa couplings, y 2 d , y 2 s , and y 2 b as functions of m 1/2 for the set of inputs used in Fig. 1 is shown in Fig. 2. We see that while the differences between A and B do not manifest themselves in any of the down-quark Yukawa couplings, they do differ from the NF choice for the first two generations. The flavor-violating contributions are negligible for the bottom quark because m b ≫ m d,s , and the three choices considered give results that are nearly identical. The fact that the difference between y 2 d and y 2 s is larger for choices A and B than for the NF choice is a reflection of the larger magnitudes of the off-diagonal Yukawa couplings.

Hadronic Uncertainties
In addition to the WCs, the proton decay amplitudes in Eq. (59) depend on hadronic matrix elements. As discussed in detail in [18], in order to apply the limit in Eq. (46), one needs to know not only the central values of the matrix elements but also their uncertainties. The relevant systematic uncertainties of the form factors were taken into account for the first time in [41]. The total uncertainties found in K final states were 20%-40%, whereas they were 30%-40% for π final states, which were reduced to 10% -15% in [42]. The uncertainties in all of the matrix elements in Eq. (46) must be taken into account in order to determine the region of parameter space for which τ (p → K +ν ) > 6.6 × 10 33 yrs. For the matrix elements contributing to the relevant amplitude A(p → K +ν τ ) in Eq. (59), Ref. [42] found where we have quoted the total error obtained by combining the statistical and systematic errors in quadrature. We note that the matrix elements K + |(us) L d L |p and K + |(ud) L s L |p are the most relevant, since the C LL coefficients dominate over C RL .

Flavor Violation
The embedding of the MSSM in SU (5), as in either Eq. (5) or (9), can make an important difference. In particular, the different embeddings for the SU(5) Yukawa matrices h 10 and h 5 lead to different effective mass matrices for h D , h E and h U , as we have seen in Section 2.
Minimal SU(5) corresponds to the embedding (9) (without the phases), where h D and h E are the transposes of each other. When h D involves the CKM matrix, h E inevitably leads to large right-handed currents, enhancing the branching ratio BR(µ → eγ), which can be written as where we use the notation in [12] for the amplitude of the decay µ → eγ. The experimental upper limit BR EXP (µ → eγ) ≤ 4.2 × 10 −13 [43] imposes the constraints These limits on the coefficients a µeγL and a µeγR constrain the amount of flavor violation mediated by charginos and neutralinos in the MSSM. We note that care must be taken in an analysis in terms of mass-insertion operators in the presence of off-diagonal entries in all the soft supersymmetry-breaking sectors, because there are correlations among the elements of the matrices and some cancellations may occur.
In order to understand the order of magnitude of possible contributions to a µeγR and a µeγL that are consistent with the limits in Eq. (64), we consider simplified formulae for the neutralino contributions. There are significant contributions coming from chargino exchange, but these are suppressed relative to the neutralino exchange contributions. Fig. 3 displays the diagrams making the most important contributions to a µeγR . First the neutralino exchange diagram is shown in the mass-eigenstate basis, and then we identify four main contributions in the interaction basis. The contributions from a (I) µeγR , which requires a mass insertion outside the loop, can be approximated as As we will see when we consider specific models in Section 4, (m 2 E ) 12 is similar in both the cases A and B. Similarly the contributions a Here N 11 and N 13 are mixing elements of the neutralinos, with N 11 ≈ 1 when the lightest neutralino is mainly bino, and N 31 characterizes the mixing between the the HiggsinoH 0 d and the bino. These diagram factors are also proportional to (m 2 E ) 12 , a common factor between cases A and B. Although the diagram a (IIa) also depends on (a E ) 22 , this quantity is also similar in cases A and B for the models we consider below. In contrast, the diagram corresponding to a (IIc) µeγR , which may be approximated by is proportional to (a E ) 21 . This mixing term is very different in cases A and B and can lead to differences in the total value of a µeγR by an order of magnitude or more, as we see below. There are similar contributions to a µeγL , but they are mediated by (m 2 L ) 12 instead of (m 2 E ) 12 , and hence suppressed for these models, as we see below in Section 4. The reason why the right-handed contribution to BR(µ → eγ), which is encoded in |a µeγR | and associated with m 2 E , is significantly larger than m 2 L is that m 2 E is matched to m 2 10 at M GUT , and m 2 L to m 2 5 . Both m 2 10 and m 2 5 start at zero at M in (see Eq. (28)), but they evolve differently: Consequently, (m 2 10 ) ii is typically twice as large as (m 2 5 ) ii at M GUT (see, e.g., Fig. 2 of [12]). We conclude this discussion by noting that the future MEG II experiment is expected to be sensitive to BR(µ → eγ) = 6 × 10 −14 [44].

µ → eee
There are other proposals for future experiments that are sensitive to muon flavor violation, e.g., to the µ → eee mode. The current experimental limit on this mode is provided by the SINDRUM experiment: BR(µ → eee) < 1.0 × 10 −12 [45], and the Mu3e experiment aims at a sensitivity of ∼ 10 −16 in the future [46].
The µ → eγ dipole processes shown in Fig. 3 give the dominant contributions to µ → eee in many supersymmetric models [47]. 12 In this case BR(µ → eee) is related to BR(µ → eγ) by [48] BR(µ → eee) This relation indicates that currently µ → eee gives a much weaker limit on lepton flavor violation than µ → eγ, but will offer a better sensitivity in the future.

µ → e conversion
Another promising process is µ → e conversion on a nucleus. The tightest current experimental bound on the µ → e conversion rate is provided for gold nuclei by the SINDRUM II collaboration: BR(µ + Au → e + Au) < 7 × 10 −13 [49]. . We infer that µ → e conversion processes may be more promising than µ → eγ and µ → eee in the future.

Electric dipole moments (EDMs)
The new limit on the electron EDM, |d e | < 1.1 × 10 −29 e.cm [54] could in principle constrain parts of the parameter space that would otherwise be allowed if no flavor-violating terms in the soft terms were considered. At the one-loop level, there are supersymmetric contributions to the electron EDM mediated by charginos and neutralinos. A general expression is given by [55] where In these expressions, the sfermion mass matrix, in the basis where Yukawa couplings are diagonal, for each family is given by and m f , f = d, e, u, are the masses of the fermions corresponding to the sfermions F = D, E, U. Note that the indices 1L and 1R of (K F ) k,1L and (K * F ) k,1R in Eq. (71) for F = E correspond to e L and e R , respectively, which are the external lines in Fig. 4. The index k corresponds to the mass eigenstates from k = 1, . . . , 6, where the sfermion mass eigenstates are defined by such that M 2 F = K F M 2 K † F is a diagonal matrix. In the models considered here, the lightest neutralino,χ 0 1 , typically gives the dominant contribution to the electron EDM in Eq. (70). This term is proportional to N 1i K Ek1 , so the most important contribution to d e comes from and where the m 2 e k are the slepton mass eigenstates. This contribution is depicted in the left Feynman diagram of Fig. 4 in the flavor basis, and we see that, in the absence of off-diagonal and imaginary terms in a E , the EDM is zero. However, once the CKM matrix is introduced to seed flavor violation, as in the flavor choices A and B discussed earlier, imaginary parts appear in the soft squared-mass matrices (m 2 L ) 1j , (a E ) ij and (m 2 E ) 1j . We note that the function B in Eq. (70) varies slowly over the range ∼ 0.2 to ∼ 0.3 for all of the spectra we consider and for all the indices k. Therefore the individual contributions in the terms of Eq. (70) depend mainly on the combination The imaginary part above can be easily understood in terms of the second diagram of Fig. 4, since where b, c = 1, 2, 3. The imaginary parts of each of the contributions to the sum above can be written as We find that there are important contributions from the terms involving Re(a E ) 33 but, depending on the model, contributions containing (a E ) 11 can dominate for models with the flavor choice A, and contributions containing (a E ) 21 and (a E ) 31 can also be important. In Eq. (78) with b = c = 3, we find that in the models considered in Section 4 Then the dominant term in Eqs. (77) and (78) contains Re(a E ) 33 and reduces to In the models considered below, this contribution is similar in both of the choices A and B. However, for choice B, the contribution from (a E ) 31 can also be important, as we will see in Section 4. It is relatively easy to understand how the contribution from (a E ) 11 can dominate in models with flavor choice A relative to choice B. As seen in Eqs. (77) and (78), when b = c = 1 (and noting that the imaginary parts of (m 2 E ) 11 and (m 2 L ) 11 both vanish), the term containing (a E ) 11 reduces to As we will see in Section 4 below, Im {(a E ) 11 } is typically four orders of magnitude larger in choice A than in choice B. This can be traced to the matching condition in Eq. (29). In fact, the contribution containing (a E ) 11 can be even larger than that containing (a E ) 33  When two or more contributions are important we can still use the formula above for each slepton, taking the signs of the Im {η E1k } into account. Overall, therefore, we find that while the contributions from b = c = 3 for choices A and B are similar, the contribution from b = c = 1 is much greater in choice A, and we expect the EDM to be larger in choice A than in choice B.

Models
In the continued absence of supersymmetry at the LHC, the allowed parameter space in constrained supersymmetric models has been pushed to ever higher mass scales [40,17,11,56,57,18,58,59]. For this reason, also in order to obtain a Higgs boson with mass consistent with the experimental value, M h ≃ 125 GeV [3], and a sufficiently long proton decay lifetime [60], supersymmetric mass scales in the range from 1 to 5 TeV are favored. Then, the requirement that the relic dark matter density agree with Planck results, Ω χ h 2 ≃ 0.12 [61], imposes significant constraints on models and their parameters, as do the upper limits on dark matter scattering on matter [62]. It was found in the context of no-scale supergravity models that one or both of the MSSM Higgs fields must be twisted [11], i.e., they must acquire masses different from the universal masses for squarks and sleptons, which vanish at the input scale in no-scale models. 13 Models with universal input scalar masses suffer from tension between the Higgs mass measurement, proton decay limits and the cosmological relic density. With all fields untwisted it was possible to find parameters with a sufficiently large Higgs mass and acceptable relic density or long proton lifetime, but not both [11].
As discussed earlier, the trilinear and bilinear soft terms depend on the nature of the twisted Higgs fields, and on the assignments of the modular weights that appear in the superpotential (see Eq. (15)). We outline here the sample model classes that we use for our analysis, which are adapted from some studied previously in [11]. The models are distinguished by the parameters p and q that take values 0 or 1 depending on whether the H and H fields are twisted or not, as well as the choices of modular weights. Here, we take p = 1, and allow q to take values 0 or 1. Once these are specified, the models have six free continuous parameters and one sign: We recall that in the absence of the dimension-five coupling, c 5 , we cannot choose independently the two GUT couplings, λ and λ ′ . In this case, typically the colored Higgs mass is low and proton decay is rapid. However, when c 5 = 0, the colored Higgs mass is sufficiently large for small λ ′ . As in previous work [11,17,18], we fix λ ′ = 10 −5 in all of the models considered here in order to ensure a sufficiently large colored Higgs mass, M H C and hence a sufficiently long nucleon lifetime. The lifetime for the dominant proton decay mode, τ (p → K +ν ), increases with λ, for which we adopt either λ = 0.6 or λ = 1.0. We take tan β = 7 in all models except M1 where tan β = 6 14 , and we choose sign(µ) > 0 in all models. We consider two values of M in : M in = 10 16.5 GeV, for which there is little RG running above M GUT , and M in = 10 18 GeV, for which the RG running is more important. We illustrate the effects of the choice of flavor structure using a subset of the models considered previously in [11]. As noted above, because of the restrictive nature of the untwisted no-scale boundary conditions, we require that either one or both of the Higgs five-plets are twisted in order to obtain simultaneously the correct relic density, Ω χ h 2 = 0.12 [61], and Higgs mass, M h = 125 GeV [3], as well as τ (p → K +ν ) consistent with the lower limit given in [36]. In Models M1 -M4 below, both Higgs multiplets are twisted, whereas for models M5 and M6, only H is twisted. The dark matter, M h and τ (p → K +ν ) constraints can all be reconciled in these models. Ref. [11] also considered models in which only H is twisted. However, we find using FeynHiggs 2.16.0 [68] a drop in the calculated Higgs mass of ∼ 2 GeV, relative to previous versions, making it difficult to reconcile an acceptable relic density with M h ≃ 125 GeV and the τ (p → K +ν ) constraint, and do not consider further such models. The models considered here are as follows: (M1) In this model, we set M in = 10 16.5 GeV, tan β = 6, λ = 0.6, with p = q = 1, and we take all modular weights α F = β S = 0. In this case, This model is similar to that considered in the left panel of Fig. 3 in [11].
(M3) In this model, we consider M in = 10 18 GeV, tan β = 7, λ = 0.6. We again take p = q = 1, with the same modular weights as adopted in M2. This model is similar to that considered in the right panel of Fig. 4 of [11].
(M4) In this model, we consider M in = 10 18 GeV, tan β = 7, λ = 1. We again take p = q = 1, with the same modular weights as adopted in M2. This model is the same as that considered in the right panel of Fig. 4 of [11].
For each of the models M1 -M6, we compute the proton decay lifetime, BR(µ → eγ), and the induced electron EDM, comparing the flavor choices A and B, and also comparing the predictions of the NF scenario for the proton lifetime.

Both Higgs fields in twisted sectors
Since both Higgs five-plets are twisted in models M1 and M2, we must use p = q = 1 in Eq. (15), yielding m H = mH = m 1 = m 2 = m 3/2 at M in , whereas all the other scalar masses vanish there. Once the modular weights α and β appearing in Eq. (15) are specified, all of the bi-and tri-linear terms are fixed relative to m 3/2 , so the models are fully specified. In model M1, we take all modular weights to vanish, yielding the non-zero A-terms A 10 = A5 = m 3/2 and A λ = 2m 3/2 , as well as a non-zero B-term for B H = 2m 3/2 . In models M2 -M4, we take all A-and B-terms to vanish. In what follows, we display our results in (m 1/2 , m 1 ) planes.
Model M1 In this model we fix M in = 10 16.5 GeV, so there is little super-GUT running between M in and M GUT , tan β = 6 and µ > 0. The chosen values of the couplings of the adjoint Higgs supermultiplets are λ = 0.6 and λ ′ = 0.00001. We show in the upper left panel of Fig. 5 the (m 1/2 , m 1 ) plane for this model, where we recall that m 1 = m 3/2 in this model. There is no EW symmetry breaking (EWSB) in the triangular region shaded pink in the upper left corner, i.e., the solution for the MSSM µ parameter has µ 2 < 0. The dark blue shaded strip just below the no-EWSB region corresponds to the focus point [69], with the relic density taking values in the range 0.06 < Ω χ h 2 < 0.2. This is wider than the range determined by Planck [61], but we show an extended range in order to make it more visible on the scale of this figure Proton lifetime: Also shown in Fig. 5 are predictions for the proton lifetime, τ (p → K +ν ). For each case considered, we show 3 sets of 3 contours each, corresponding to the current lower limit on τ (p → K +ν ). The central contour in each set uses the central values for the hadronic matrix elements given in Eq. (62), and the outer contours to either side correspond to the ±1σ variations in these matrix elements indicated there, keeping the masses m s (2 GeV) and m c (2 GeV) fixed at their central values. The sets of solid and dashed blue contours correspond to model choices A and B, respectively, and we see that the choice between these flavor embeddings has very little effect the proton lifetime. Along the dark matter strip, the lower limit on τ (p → K +ν ) (assuming central values of the hadronic matrix elements) corresponds to m 1/2 5(7)    For choice A BR(µ → eγ) is always below the current experimental upper limit of 4.2 × 10 −13 [43]. In the region of greatest interest along the blue relic density strip, the branching ratio may exceed 10 −17 , but a small portion at low (m 1/2 , m 1 ) where it reaches 10 −16 is excluded by the proton decay limit. Moreover, BR(µ → eγ) decreases significantly below the strip and at larger masses. The low values for the branching ratio arise primarily from the choice of an embedding in which h E is diagonal at the EW scale (see Eq. (45)).
In contrast, for choice B the lepton Yukawa couplings are not diagonal at the EW scale and we see in the right panel of Fig. 6 that BR(µ → eγ) is significantly larger, with values above 10 −16 becoming consistent with τ (p → K +ν ), M h and the relic dark matter density. Indeed, BR(µ → eγ) is larger than 10 −18 even at very large gaugino masses > 10 TeV. We note that in this case that the dependence of BR(µ → eγ) on m 1/2 is much stronger than that on m 1 . Nevertheless, there is a stretch of the focus-point strip with 4.5 TeV m 1/2 6 TeV, compatible with the present limit on τ (p → K +ν ) and the Higgs mass, where µ → e conversion may be accessible to the PRISM experiment [52]. 16 In order to understand this behavior, we analyze a benchmark point in model M1 lying on the relic density strip with m 1/2 = 6000 GeV, which corresponds to m 1 = 9070 GeV, and a Higgs mass of M h = 125.2 ± 0.9 GeV according to FeynHiggs 2.16.0. We show in Table 1 the relevant mass parameter values in model M1 that are used to extract the approximate values for τ (p → K +ν ), BR(µ → eγ) and the electron EDM. As one can see, there is essentially no difference in the selectron masses between cases A and B and only a 2% difference in the the smuon masses. As we discussed earlier, (m 2 L ) 12 ≪ (m 2 E ) 12 so that a µeγL is suppressed. We also see that in all models (M1-M4), (m 2 E ) 12 is within a factor of two and (a E ) 22 is nearly identical between cases A and B. The difference seen in Fig. 6 between cases A and B is a result of a (IIc) µeγR (see Eq. (67)), which is proportional to (a E ) 21 and is more than a factor of 10 3 times larger in case B due to the choice of U E R = V * CKM as opposed to U E R =1 in case A. We note that the predictions for τ (p → K +ν ) are similar for flavor choices A and B, beyond the current limit but well within the projected reach of Hyper-Kamiokande [37]. While the predictions for BR(µ → eγ) and the electron EDM differ for flavor choices A and B, they lie significantly below the current limits and prospective experimental sensitivities.
Similar behavior is found for model M2, shown in the lower two panels of Fig. 6. Along the relic density strip (now at lower m 1 relative to M1), the branching ratio ranges from 10 −20 to a few ×10 −19 for flavor choice A. As we saw for model M1, the branching ratio is considerably larger for flavor choice B and may be as large O(10 −16 ) while remaining consistent with proton decay limits. A representative benchmark point along the relic density strip at m 1/2 = 6000 GeV for M2 is also given in Table 1, with m 1 = 5950 GeV. For this point M h = 123.6 ± 0.7 GeV, which is consistent within the uncertainties with the experimental value. In this case, we again see that the dominant difference in the branching ratio between choices A and B is due to (a E ) 12 . The predictions for τ (p → K +ν ) are again similar for flavor choices A and B, and are consistent with the current limit within the current matrix element uncertainties. As in the case of model M1, the predictions for BR(µ → eγ) and the electron EDM again differ for flavor choices A and B, while lying significantly below the current limits.  In Table 2, we show representative benchmark points for models M3 and M4 along the relic density strip with m 1/2 = 6000 GeV. For M3, m 1 = 7850 GeV, giving M h = 124.5 ± 0.7 GeV and for M4, m 1 = 9780 GeV, with M h = 124.4 ± 0.7 GeV. We again see that the large increase in (a E ) 12 in choice B relative to A accounts for the increase in BR(µ → eγ). In both cases, τ (p → K +ν ) should be within reach of the Hyper-Kamiokande experiment [37], but the predictions for BR(µ → eγ) and the electron EDM are below the projected experimental sensitivities for both flavor choices.
Electron EDM: In the upper panels of Fig. 8 we show the values of the electron EDM (eEDM), d e , calculated using SUSY FLAVOR [70] in model M1, presented in the corresponding (m 1/2 , m 1 ) plane used in Figs. 5 and 6. In the absence of flavor effects and in the absence of complex phases in the supersymmetric parameters (which we do not consider here), the EDM would be zero. Once the CKM matrix is introduced as a seed of flavor and CP violation, the CKM phase propagates in all of the spectra, generating a non-zero eEDM. The values of d e displayed in Fig. 8 are for the flavor choices A (left) and B (right).
We see that the eEDM is generally larger for choice A reaching ∼ 10 −32 e.cm in the portion of the dark matter strip that is consistent with τ (p → K +ν ) and M h . We also see that the eEDM is roughly a factor of 10 larger for case A than it is for case B. As one can see from Table 1, the contribution from Eq. (80) in both cases A and B are similar. However since Im(a E ) 11 is about four orders of magnitude larger for case A relative to case B, the contribution from Eq. (81) boosts the eEDM in case A. For M1 and case A, indeed the most important contributions come from (a E ) 11 and (a E ) 33 , where Due to the overall sign in Eq. (75), this is a positive contribution to the eEDM. However, for the choice B the contribution containing (a E ) 11 is negligible due to the smallness of Im {(a E ) 11 }, but the contribution containing (a E ) 31 becomes important, and we have instead 17 Although all the other contributions in the cases A and B above are small, we keep them in Eq. (77) and get with Eq. (82) The reader should keep in mind that the approximation of Eq. (82) should give the right order of magnitude, but the exact numerical factor is difficult to obtain with this approximation, due to the detailed structure of the complete 6 × 6 diagonalization matrices K E .    This range of eEDM values is well below the experimental limit, and with flavor choice B the eEDM remains below 10 −33 e.cm along all the dark matter strip. Indeed, the eEDM falls precipitously as the relic density strip is approached, changing sign as it passes through zero in the thin cross-hatched region, where its magnitude is below 10 −34 e.cm. We see in (85) that Im (K E ) k,1L (K E ) * k,1R 31 > 0 in choice B, and so the opposite signs in (85) would explain the change in sign in d e with respect to choice A if in the exact diagonalization we had as is the case in the full numerical calculation, causing d e to become negative. The value of d e is reduced in choice B, with respect to A, due to a cancellation. We find that the values of Im (K E ) k,1L (K E ) * k,1R 33 at the benchmark point are similar in the two flavor choices, being equal to −6.1 × 10 −13 and −6.2 × 10 −13 for choices A and B, respectively. This accounts for the cancellation in B, but not in A, and we find that the sign of the eEDM at our benchmark point is indeed opposite in choices A and B.
The lower panels of Fig. 8 show the values of the eEDM in model M2 with the flavor choices A (left) and B (right), presented in the corresponding (m 1/2 , m 1 ) planes displayed in Figs. 5 and 6. We see again that larger values of the eEDM are found with flavor choice A than with choice B: < 10 −33 e.cm compared with 10 −34 e.cm. For M2, the contribution from (a E ) 11 dominates over the contribution from (a E ) 33 : whereas for choice B Thus this analytic approximation accounts for an enhancement by a factor of roughly 4.5, the full numerical ratio between A and B being about 22.
In Fig. 9 we show the values of the electron EDM in models M3 (upper panels) and M4 (lower panels), on the corresponding (m 1/2 , m 1 ) planes displayed in Figs. 5 and 7, for flavor choices A (left panels) and B (right panels). We see that the eEDM is generally larger in model M3 than in model M4, and larger with flavor choice A than with flavor choice B. However, along the dark matter strips the τ (p → K +ν ) constraint generally imposes d e < 10 −33 e.cm, except in the case of model M4 with choice A, for which d e may reach a few ×10 −33 e.cm. Note that M4, for both cases A and B, all the elements (a E ) 11 , (a E ) 31 and (a E ) 33 are important.

Models in which onlyH is twisted
As mentioned above, unlessH is twisted we find no solutions for which the relic density, Higgs mass and proton lifetime are consistent with experiment. In this Section, we consider models in which onlyH twisted, i.e., we leave H untwisted so that m 0 = m 1 = 0, whereas m 2 = m 3/2 . As in the previous Section, we consider two choices for the modular weights, one in which the modular weights take values such that all tri-and bi-linear terms vanish, and another in which the weights all vanish, leaving some of the tri-and bi-linear terms non-zero. These are labelled models M5 and M6, respectively.
Model M5 In this model we fix M in = 10 18 GeV, tan β = 7, and µ > 0. The chosen values of the couplings of the adjoint Higgs supermultiplets are λ = 1 and λ ′ = 0.00001. We show in the left panel of Fig. 10 the (m 1/2 , m 2 ) plane for this model, where we recall that m 2 = m 3/2 when onlyH is twisted. There is no EW symmetry breaking (EWSB) in the triangular region shaded pink in the upper left corner, i.e., the solution for the MSSM µ parameter has µ 2 < 0. As in Fig. 5, the red dot-dashed curves show contours of the Higgs mass as calculated using FeynHiggs 2.16.0 [68], and there is a dark blue shaded strip just below the no-EWSB region, corresponding to the focus point [69], where the relic density taking values in the range 0.06 < Ω χ h 2 < 0.2. In addition to this strip, there is a band at lower m 2 , which corresponds to a funnel where rapid annihilation via direct-channel H/A poles when m χ ≃ M H/A /2 brings the relic density into this range. This band actually consists of two unresolved narrow strips with m χ > and m χ < M H/A /2, between which the relic density takes lower values. We note that this funnel strip ends when M h < 124 GeV. Beyond this endpoint, the suppression in the annihilation cross-section due to the large value of m 1/2 is strong enough that the relic density always exceeds the observed value, i.e., Ω χ h 2 > 0.12, even on the H/A poles.

Model M6
The corresponding results for model M6 are shown in the right panel of Fig. 10. In this case, A 10 = A λ = B H = m 2 , and the other bi-and tri-linear terms vanish at M in . The dependence of the A-terms on m 2 induces a weak dependence of M h on m 2 , as is readily seen by comparing the Higgs mass contours in the two panels of Fig. 10. Importantly, in this case we do not find the focus-point strip along the boundary of the no-EWSB region. Indeed, in this case m 2 A < 0, where m A is the pseudoscalar Higgs mass, everywhere in the pink shaded region. There is, nevertheless, a funnel strip at lower m 2 , which extends only as far as m 1/2 ∼ 3 TeV, where M h ∼ 123 GeV in this case.
Proton lifetime: The proton lifetime limits for M5 are weaker than those in M6. In both cases, there is little flavor dependence and the proton lifetimes for cases A and B are nearly identical and also similar to the NF case. In the case of M5, we only see two sets of lines, as the −1σ variations in the hadronic matrix elements, which increase the proton lifetime for fixed supersymmetric model parameters, push the contour for τ p = 6.6 × 10 33 yrs to low values of m 1/2 < 2 TeV, below its displayed range and where M h is too small. There is considerable parameter space in model M5 where the relic density Ω χ h 2 ≈ 0.12 and the proton lifetime constraint is satisfied. In contrast, in model M6 only the portion of the funnel strip between m 1/2 = 2.5 and 3 TeV satisfies the proton decay constraint when the  Flavor violation: Values of the branching ratio BR(µ → eγ) for models M5 and M6 are shown in Fig. 11 (upper and lower panels, respectively) for flavor choices A and B (left and right panels, respectively). As was the case in model M1, for flavor choice A, the branching ratio in model M5 exceeds 10 −17 only at very low m 1/2 , in this case for m 1/2 < 1.4 TeV, for which M h < 121 GeV. For the portion of the focus-point strip with M h > 123 GeV, BR(µ → eγ) < 10 −18 . Furthermore, the branching ratio is over an order of magnitude smaller in the funnel strip than it is in the focus-point strip. In contrast, for choice B, there is little difference in BR(µ → eγ) between the two relic density strips. We again see that overall the branching ratio for choice B is significantly larger than for choice A. We provide in Table 3 the parameters of a benchmark point in model M5 lying on the relic density strip with m 1/2 = 6000 TeV, m 2 = 8385 GeV and M h = 124.4 ± 0.7 GeV. Once again, we see that (a E ) 21 is significantly larger for choice B than for choice A, leading to the increased branching ratio for µ → eγ.
In the case of model M6, BR(µ → eγ) is generally smaller than in M5 for flavor choice A, though it does exceed 10 −16 for m 1/2 < 1.2 TeV. For flavor choice B, the branching ratio exceeds 10 −15 for m 1/2 < 2 TeV, in the unshaded the region. At the tip of the funnel strip, the branching ratio exceeds 10 −16 . In a portion of this strip that is compatible with the present limit on τ (p → K +ν ) and the Higgs mass, µ → e conversion may be accessible to the PRISM experiment. [52] In this case, because of the lack of a focus-point strip, we provide also in Table 3 the parameters of a benchmark point on the funnel strip with m 1/2 = 3000 TeV and m 2 = 4470 GeV, corresponding to M h = 123.2 ± 0.8 GeV. The differences between the branching ratios in choices A and B can again be attributed to the increase in (a E ) 21 for choice B seen in the Table. Electron EDM: Predictions for the electron EDM in models M5 (upper panels) and M6 (lower panels) are shown in Fig. 12, again with flavor choice A in the left panels and flavor choice B in the right panels. Predictions are everywhere significantly below the present experimental sensitivity. Overall, we see that the predicted values are somewhat smaller in model M6 than in model M5, and somewhat larger with choice A than with choice B. In the most favorable case, namely model M5 with flavor choice A, the electron EDM varies between 10 −32 e.cm and 10 −34 e.cm along the focus-point strip, and between 10 −32 e.cm and 10 −33 e.cm along the rapid-annihilation strip. In the least favorable case, namely model M6 with flavor choice B, the electron EDM is below 10 −35 e.cm along all the rapid-annihilation  11 [GeV] 5.9 × 10 −2 e −i 4.7×10 −6 5.9 × 10 −2 e i 10 −9 (a E ) 22 Table 3: Benchmark points in model M5 with m 1/2 = 6000 GeV and in model M6 with m 1/2 = 3000 GeV. For M5, m 2 = 8385 GeV, and for M6, m 2 = 4470 GeV. We list values of the parameters relevant for BR(µ → eγ) and the electron EDM obtained with flavor choices A and B, as well as the corresponding predictions for τ (p → K +ν ) , BR(µ → eγ) and the electron EDM.
strip. In the case of M6, the approximation for |d e | we use in Eq. (77) gives the correct order of magnitude for choice A, but falls short for choice B by an about an order of magnitude, as the contributions of other elements in the matrix (a E ) bc must be taken into account when determining the total value of d e . As for the previous benchmarks, in both models M5 and M6 τ (p → K +ν ) is within reach of Hyper-Kamiokande [37], whereas BR(µ → eγ) and the electron EDM lie below the prospective future experimental reaches.

Overview and Conclusions
We have studied in this paper the phenomenological scope for SU(5) super-GUTs, in which variants of no-scale boundary conditions are imposed on the soft supersymmetry-breaking parameters at some input scale M in > M GUT . Specifically, the soft supersymmetry-breaking scalar masses for the squarks and sleptons vanish at M in , whereas those for the 5 and 5 Higgs supermultiplets depend whether they have twisted boundary conditions at M in , as seen in (15), where other details of the boundary conditions such as modular weights can be found.
In addition to these input conditions, the low-energy phenomenology of such models depends on the magnitude of the hierarchy between M GUT and M in , for which we consider the illustrative values of 10 16.5 and 10 18 GeV. We consider the constraints on such no-scale SU(5) models that are imposed by the cosmological density of cold dark matter, τ (p → K +ν ) and M h . We find that the Higgs field responsible for the charge-2/3 quark masses must be twisted, while twisting the Higgs responsible for the charge-1/3 and charged-lepton masses is optional.
Within this general framework, we have considered six specific choices for the input boundary conditions. In addition to M in , modular weights and GUT Higgs trilinear couplings -see (1) -these include possible dimension-5 effects on GUT unification. We emphasize also that the super-GUT running between M in and M GUT is sensitive to the way in which the MSSM matter fields are embedded into GUT supermultiplets, and specifically the underlying origin of CKM flavor mixing. For each of these six no-scale models, we have considered two choices for flavor mixing, which yield predictions for BR(µ → eγ) and the electron EDM that are quite different, but less so for τ (p → K +ν ). We also contrast their predictions for τ (p → K +ν ) with those made when neglecting off-diagonal entries in the Yukawa coupling matrices. Though the differences in τ (p → K +ν ) between the two flavor choices are small, the differences from when the mixing is neglected may be larger than the uncertainties associated with hadronic matrix elements in some cases. We note that the ranges of sparticle mass parameters favored by the dark matter density and M h (as calculated using FeynHiggs 2.16.0) are generally beyond the current τ (p → K +ν ) limit as well as the reach of the LHC.
As can be seen in the various panels of Figs. 6, 7 and 11, the predictions for BR(µ → eγ) are strongly dependent on the flavor choice as well as the choice of no-scale model. However, in all cases except portions of the dark matter strips in models M1 and M6 with flavor choice B, the value of BR(µ → eγ) lies significantly below the current and projected experimental sensitivities. The electron EDM is also below the current and projected experimental sensitivities, as can be seen in Figs. 8, 9 and 12. On the other hand, there are significant regions of parameter space for all models where τ (p → K +ν ) is within reach of the Hyper-Kamiokande experiment. As seen in the Tables, this is in particular the case for all the benchmark points highlighted there.
These examples demonstrate explicitly that there is no supersymmetric flavor problem in no-scale models, the reasons being that the no-scale boundary condition that every soft supersymmetry-breaking matter scalar mass vanishes at the input scale M in is flavoruniversal, and that the leading-order renormalization by gauge interactions is also flavoruniversal. Nevertheless, τ (p → K +ν ) may well be within reach.