text:
# 1 Introduction

## 1 Introduction

In low energy supersymmetric extensions of the Standard Model (SM), unless sparticle masses are considerably increased, present limits on flavour violating (FV) decays and electric dipole moments (EDMs) respectively allow for a quite small amount of fermion-sfermion misalignment in the flavour basis and constrain the phases in the diagonal elements of sfermion masses, involving the parameters $`\mu `$ and $`A`$, to be rather small. The bounds on the supersymmetric contribution to lepton (L)FV decays and EDMs have the advantage, as compared to the corresponding squark sector ones, of being not biased by the SM contribution - nor by the non-supersymmetric seesaw  contribution . Experimental limits on LFV decays and EDMs are then a direct probe of the flavour and CP pattern of slepton masses - see e.g. Ref.  for a recent collection. From a theoretical perspective, understanding why CP phases and deviations from alignment are so strongly suppressed is one of the major problems of low energy supersymmetry, the CP and flavour problem. On the other hand, precisely because FV decays and EDMs provide strong constraints, they can share some light on the features of the (possibly) supersymmetric extension of the SM.

Indeed, it is well known that even if these CP phases and misalignments were absent - or suppressed enough - from the effective broken supersymmetric theory defined at $`M_{Pl}`$, they would be generated at low energy by the RGE corrections due to other flavour and CP violating sources already present in the theory, in particular from the Yukawa couplings of heavy states like the right-handed neutrinos of the seesaw model  and the Higgs triplets of $`SU(5)`$ grand unified theories (GUT) . Effects of radiative origin have the nice features of being naturally small, exactly calculable once specified a certain theory and - most interestingly - if experiments are sensitive to them, they yield limits on some combination of Yukawa couplings 111Of course, several effects of different origin could be simultaneously present, but there would be no reason for a destructive interference among them to occur..

Recently, it has been stressed that the present (planned) limit on $`\mu e\gamma `$  ($`\tau \mu \gamma `$ ) is sensitive to the lepton-slepton misalignments induced by the radiative corrections in the framework of the supersymmetric seesaw model and that the associated constraints on neutrino Yukawa couplings have a remarkable feedback on neutrino mass model building . However, LFV decays cannot provide any information on the pattern of CP violation of slepton masses, while EDMs are sensitive to both LF and CP violations. Since the present sensitivities to the electron and muon EDMs, $`d_e<10^{27}`$ e cm  and $`d_\mu <10^{18}`$ e cm , could be lowered by planned experiments by up to three to five  and six to eight  orders of magnitude respectively, it is natural to wonder whether lepton EDMs would explore the range associated to LF and CP violations of radiative origin.

In this work we analyze the predicted range for $`d_e`$ and $`d_\mu `$ when the radiative corrections due to the right-handed neutrinos of the seesaw model and, possibly, the Higgs triplets of $`SU(5)`$ GUT provide the main source of CP violation and misalignment in slepton masses. This allows to extract many informations. If experimental searches could explore this range, one could obtain limits on the imaginary part of the relevant combination of neutrino Yukawa couplings. Moreover, the eventual discovery of a lepton EDM in this range might be interpreted as indirectly suggesting the existence of such a kind of fundamental particles which are too heavy to be more manifest. On the other hand, finding $`d_e`$ and/or $`d_\mu `$ above this predicted range would prove the existence of a source of CP (and likely also LF) violation other than these heavy states.

The pure seesaw case has been considered in Ref.  \- see also the related studies  on specific seesaw textures -, where it has been pointed out that threshold effects due to hierarchical right-handed neutrino masses enhance the radiatively-induced $`m(A_{ii})`$, $`i=e,\mu ,\tau `$. By extensively reappraising this framework, we find that for $`\mathrm{tan}\beta 10`$ the amplitude with a LL RR double insertion - proportional to $`\mathrm{tan}^3\beta `$ \- dominates over the one involving $`m(A_{ii})`$ \- insensitive 222In Ref.  a dependence on $`\mathrm{tan}\beta `$ arises due to the particular class of seesaw textures studied. to $`\mathrm{tan}\beta `$ -, the exact ratio depending on the particular choice of supersymmetric mass parameters. We provide general expressions from which it is easy to recognize the model dependencies and which complete previous analyses. Lepton EDMs are then strongly enhanced in models with large $`\mathrm{tan}\beta `$. To see how close experiments are getting to the seesaw induced lepton EDMs range, we compare the upper estimates for the radiatively-induced misalignments and CP phases with the corresponding present and planned experimental limits collected in .

The range for a seesaw-induced $`d_\mu `$ turns out to be quite far from present searches and, in particular, its eventual discovery above $`10^{23}`$ e cm would prove the existence of a source of CP violation other than the neutrino Yukawa couplings. On the contrary, the present sensitivity to $`d_e`$ explores the seesaw-induced range in models with large $`\mathrm{tan}\beta `$ and small R-slepton masses, say $`m_R`$ around $`100200`$ GeV. The planned improvements for $`d_e`$ would allow to test also models with $`m_R`$ up to the TeV region and moderate $`\mathrm{tan}\beta `$. When present or planned experimental limits turns out to overlap with these allowed ranges, bounds on the imaginary part of the relevant combinations of neutrino Yukawa couplings are obtained and plotted in the plane $`(\stackrel{~}{M}_1,m_R)`$, respectively the bino and the average R-slepton masses at low energy. These plots allow to check whether any particular seesaw model is consistent with present data and, if so, which level of experimental sensitivity would test it. In any case, an experimentally interesting contribution requires hierarchical right-handed neutrino masses.

This feature is no more necessary when, in addition to the seesaw, a stage of $`SU(5)`$ grand unification is present above the gauge coupling unification scale. It is well known that the main drawback of minimal $`SU(5)`$ is that the value of the triplet mass, $`M_T`$, required by gauge coupling unification is sizeably below the lower bound on $`M_T`$ derived from proton lifetime (see for instance  for recent reviews and references). In our analysis we therefore keep the triplet mass as a free parameter. We nevertheless exploit the minimal $`SU(5)`$ relations between the doublet and triplet Yukawa couplings since in general they are mildly broken in non-minimal versions of $`SU(5)`$. The simultaneous presence of right-handed neutrinos and heavy triplets turns out to further enhance the amplitude with the LL RR double insertion over the one with $`m(A_{ii})`$, essentially unaffected by triplets. This cannot be derived by the - otherwise elegant - technique based on the allowed invariants .

When right-handed neutrinos and Higgs triplets are simultaneously present, the predicted range for the radiatively-induced $`d_e`$ is already sizeably excluded by the present experimental limit. This in turn is translated into a strong constraint on the imaginary part of the combination of Yukawa couplings which is relevant for the LL RR amplitude. Such a constraint could be hardly evaded even for large values of the triplet mass and unfavorable supersymmetric mass parameters. The radiatively-induced $`d_\mu `$ in the presence of triplets should not exceed $`10^{23}`$ e cm, as was the case for the pure seesaw. However, at difference of the latter case, in the former one $`d_\mu `$ only mildly depends on the spectrum of right-handed neutrinos. Notice that planned searches for $`d_\mu `$ could constrain (depending on the triplet mass mass, of course) the imaginary part of the combination of Yukawa couplings whose absolute value could be independently constrained by the LFV decay $`\tau \mu \gamma `$.

The paper is organized as follows. In Section 2 we introduce our notations, discuss the framework and make some preliminary considerations. Section 3 considers the pure seesaw case by separately analyzing the flavour conserving (FC) and flavour violating (FV) amplitudes contributing to lepton EDMs. In Section 4 the framework of seesaw and $`SU(5)`$ is discussed along the same lines. Concluding remarks are drawn in Section 5. Finally, in Appendix A and B we collect the RGE in the case of the seesaw, without and with a minimal $`SU(5)`$ unification respectively.

## 2 Framework and Method

In this section we recall the expression for the lepton EDMs in the mass insertion approximation  to display the supersymmetric mass parameters that are constrained by the present and projected searches for $`d_e`$ and $`d_\mu `$. We then draw some preliminary considerations to introduce our procedure to calculate the radiative 1-loop contribution to these mass parameters. The relevant RGE can be found in the Appendix.

We adopt here the following conventions for the $`6\times 6`$ slepton mass matrix in the lepton flavour (LF) basis where the charged lepton mass matrix, $`m_{\mathrm{}}`$, is diagonal:

$$\left(\begin{array}{cc}\stackrel{~}{\mathrm{}}_L^{}& \stackrel{~}{\mathrm{}}_R^{}\end{array}\right)\left(\begin{array}{cc}m_{LL}^2& A_e^{}v_d\mu \mathrm{tan}\beta m_{\mathrm{}}\\ A_ev_d\mu ^{}\mathrm{tan}\beta m_{\mathrm{}}& m_{RR}^2\end{array}\right)\left(\begin{array}{c}\stackrel{~}{\mathrm{}}_L\\ \stackrel{~}{\mathrm{}}_R\end{array}\right)$$
(1)
where $`A_e`$ is the $`3\times 3`$ matrix of the trilinear coupling, the $`A`$term. All deviations from alignment in this mass matrix are gathered in the $`\delta `$ matrices, which contain 30 real parameters (including 12 phases) and are defined as:

$`m_{LL}^2=m_L^2(𝕀+\delta ^{LL})`$ $`m_{RR}^2=m_R^2(𝕀+\delta ^{RR})`$
$`A_e^{}v_d\mu \mathrm{tan}\beta m_{\mathrm{}}`$ $`=`$ $`(A_{\mathrm{}\mathrm{}}^{}v_d\mu \mathrm{tan}\beta m_{\mathrm{}})+m_Lm_R\delta ^{LR}`$ (2)
where $`m_L,m_R`$ are average masses for L and R sleptons respectively and $`A_{\mathrm{}\mathrm{}}O(m_{susy}m_{\mathrm{}}/v_d)`$ are the diagonal elements of $`A_e`$, so that $`\delta ^{LR}`$ has only non-diagonal, flavour violating, elements.

The supersymmetric contributions to $`d_i`$, $`(i=e,\mu ,\tau )`$, can be splitted in two parts, involving respectively only flavour conserving (FC) or flavour violating (FV) elements of the slepton mass matrix (1):

$$d_i=d_i^{FC}+d_i^{FV}$$
(3)
$`d_i^{FC}`$ $`=`$ $`{\displaystyle \frac{e}{2}}{\displaystyle \frac{\alpha M_1}{4\pi |\mu |^2\mathrm{cos}^2\theta _W}}[m_im(\mu )\mathrm{tan}\beta (I_B+{\displaystyle \frac{1}{2}}I_LI_R+I_2)v_dm(A_{ii}^{})I_B]`$ (4)
$`d_i^{FV}`$ $`=`$ $`{\displaystyle \frac{e}{2}}{\displaystyle \frac{\alpha M_1}{4\pi |\mu |^2\mathrm{cos}^2\theta _W}}[m_Rm_L(m(\delta ^{LL}\delta ^{LR})_{ii}I_{B,L}^{}+m(\delta ^{LR}\delta ^{RR})_{ii}I_{B,R}^{})`$ (5)
$`+\mathrm{tan}\beta (m(\delta ^{LL}\mu \eta _{\mathrm{}}m_{\mathrm{}}\delta ^{RR})_{ii}I_B^{\prime \prime }+m(\delta ^{LR}\mu ^{}\eta _{\mathrm{}}^{}m_{\mathrm{}}\delta ^{LR})_{ii}I_B^{\prime \prime })]`$
where $`\eta _{\mathrm{}}𝕀A_{\mathrm{}\mathrm{}}v_d/(m_{\mathrm{}}\mu \mathrm{tan}\beta )`$, the functions $`I`$ are defined as in  and terms that are less relevant 333E.g. a contribution to $`d_i^{FC}`$ of the form of (5), with $`\delta ^{LR}`$ $`(A_{\mathrm{}\mathrm{}}^{}v_d\mu \mathrm{tan}\beta m_{\mathrm{}})/(m_Lm_R)`$, is in principle present, but it is negligible with respect to (4). or higher order in the $`\delta `$’s matrix elements are omitted. Notice that $`\eta _{\mathrm{}}𝕀`$ for relatively large values of $`\mathrm{tan}\beta `$, favored in mSUGRA and for which the LF and CP violations are most likely to be detected. The FC and FV contributions could result from different seeds of CP violation but could also be correlated in many different ways in models. Anyway, the experimental limit can be put on both because, due to the different nature of the many parameters involved, an eventual cancellation between these contributions appears unnatural. The distinction between the FC and FV contributions is also phenomenologically relevant: some of the $`|\delta |`$’s in the FV terms are already constrained to be smaller than $`O(1)`$ by LFV decays, while $`m(\mu )`$ and $`m(A_{ii})`$ are directly constrained by the EDMs.

Our aim here is to estimate the radiative contribution to the lepton EDMs induced by the seesaw interactions, first alone and subsequently accompanied by a stage of $`SU(5)`$ GUT. Since the present experimental bounds on LFV decays and EDMs already point towards family blind soft terms (i.e. sparticles with the same quantum numbers must have the same soft terms) with very small CP phases, at the cut-off $`\mathrm{\Lambda }=M_{Pl}`$ corresponding to the decoupling of gravitational interactions we assume real and flavour blind soft terms, namely in eq. (1),

$$m_{LL}^2=m_{RR}^2=m_0^2𝕀,A_e=y_ea_0,$$
(6)
with real $`a_0`$, $`m_0`$ and $`\mu `$ term. In the next sections we separately study the FC and FV contributions to $`d_i`$ and obtain explicit approximate expressions for $`m(A_{ii})`$ and the products of $`\delta `$’s in (5). The effects of a more general family independence assumption are important but not crucial and can be easily included in our analysis. By means of these general approximated expressions we will:
i) derive the upper prediction for the radiatively induced leptonic EDMs, stressing the model dependences;
ii) compare it with the experimental limits;
iii) when allowed by the experiment, obtain an upper bound on the imaginary part of the relevant combination of Yukawa couplings.

A couple of preliminary considerations are in order before presenting the results of the next sections.

### 2.1 On the naive scaling relation

In the limit that all slepton masses are family independent, the FV contribution vanishes and the FC one is proportional to the mass of the $`i`$th lepton ($`m(A_{ii})v_dm(a_0)m_i`$) leading, except an accidental cancellation 444For dedicated studies on cancellations between amplitudes, also in more general frameworks, see e.g. . with the $`\mu `$-term amplitude, to the ”naive” scaling relation

$$d_i/d_j=m_i/m_j.$$
(7)
Then, due to the present experimental limit on $`d_e`$, $`d_\mu `$ could not exceed $`d_em_\mu /m_e210^{25}`$ e cm, which roughly corresponds to the planned sensitivity and, if the limit on $`d_e`$ were still to be lowered, next generation experiments would have no chance of measuring $`d_\mu `$.

Such considerations could provide interesting informations because (7) strictly apply only to the FC $`\mu `$-contribution to the EDM, while in general both $`m(A_{ii})`$ and the FV terms may strongly violate it. This is the case for the radiatively induced $`m(A_{ii})`$. Some of the FV contributions are instead naturally proportional to a different lepton mass, $`m_k`$, possibly heavier than $`m_i`$, as discussed in  \- and, before, for the quark sector, in . In particular, it will turn out in the next sections that the FV contribution can even take over the FC one. Hence, a value of $`d_\mu `$ above $`210^{25}`$ e cm is a possibility that deserves experimental tests and, interestingly enough, it would imply the source of lepton EDM being either the FV contribution or a non-universal $`m(A_{ii})`$, so providing a remarkable hint for our understanding of CP violation.

### 2.2 On the relevant combination of Yukawa couplings

Once a theory is specified, it is relatively easy to list the combinations of Yukawa couplings appearing in the slepton mass radiative corrections that may contribute to lepton EDMs. However, the actual calculation is more involved as we now turn to discuss.

As an example, let us look for $`m(A_{ii})`$ by studying the evolution of $`A_e`$ in the case of degenerate right-handed neutrino masses, $`\overline{M}`$, for simplicity. Let us first consider the case of the pure seesaw and define 555Here and in the following Dirac mass terms are always written as $`\overline{f}_Rm_ff_L`$. the hermitian matrices $`Ey_e^{}y_e`$, $`Ny_\nu ^{}y_\nu `$. Notice that $`E`$ is real and diagonal in the LF defining basis and $`N`$ is diagonalized by a unitary matrix similar to the CKM one with only one phase (even for non-degenerate right-handed neutrinos). By solving the RGE for $`A_e`$, eq. (37), linearly in $`t_31/(4\pi )^2`$ $`\mathrm{ln}(M_{Pl}/\overline{M})`$, $`A_e`$ can only be proportional to $`y_eE`$ and $`y_eN`$, whose diagonal elements are real. At $`O(t_3^2)`$ only $`y_eEE,y_eNN,y_eNE`$ and $`y_eEN`$ appear, whose diagonal elements are again real. A potential $`m(A_{ii})`$ shows up at $`O(t_3^4)`$, through $`m(y_eN[N,E]N)_{ii}`$. On the other hand, when also Yukawa interaction of the $`SU(5)`$ triplets are present, eq. (48) shows that, at first order in $`t_T1/(4\pi )^2`$ $`\mathrm{ln}(M_{Pl}/M_T)`$, $`A_e`$ can be proportional to $`y_eE,y_eN`$ but also to $`U^{}y_e`$, where $`Uy_u^{}y_u`$. The latter have real diagonal elements but allow at $`O(t_T^2)`$ for a combination with diagonal imaginary part, namely $`(U^{}y_eN)_{ii}`$.

Accordingly, in order to evaluate the actual coefficients in front of the products of Yukawa coupling matrices which are likely to have phases, we solve the RGE for the soft parameters by a Taylor expansion in the small parameters $`t_{if}=1/(4\pi )^2\mathrm{ln}(Q_i/Q_f)`$ associated to the intervals between the successive decoupling thresholds of the various heavy states. For instance, by integrating the RGE for $`A_e`$ in the case of $`SU(5)`$ plus seesaw, a non vanishing coefficient is obtained for the combination $`m(U^{}y_eN)_{ii}`$ at $`O(t_T^2)`$. Now, lepton FV transitions and CP phases are naturally defined in the LF basis where $`y_e`$ and the Majorana masses $`M_R`$ are diagonal and real. This basis is not invariant under the RGE evolution and one should diagonalize $`y_e`$ and $`M_R`$ again at the lower scale. Therefore, one has to find out the effect of these final rotations.

We adopt the rotating basis method introduced in Ref. , where the RGE are modified to incorporate the fact that the matrices are defined in the LF basis at each scale, so that $`y_e`$ and $`M_R`$ are always diagonal. It is worth to stress that, even if the rotations needed to diagonalize $`y_e`$ are small, their effect could be crucial for $`m(A_{ii})`$ and the four products of $`\delta `$’s in eq. (5). For instance, in the case of $`SU(5)`$ plus seesaw it turns out that this correction exactly cancels the term proportional to $`t_T^2m(U^{}y_eN)_{ii}`$ in $`A_e`$. On the contrary, these rotations can be safely neglected when deriving approximate expressions for the radiatively-induced LFV decays because the latter are only sensitive to absolute values of $`\delta `$’s.

## 3 Lepton EDMs and Seesaw

In this section we consider the predictions for lepton EDMs in the context of the supersymmetric extension of the seesaw model, namely we assume the MSSM supplemented with the seesaw Yukawa interactions and Majorana masses for the right-handed neutrinos as the effective theory valid up to the cut-off $`\mathrm{\Lambda }=M_{Pl}`$ where gravitational interactions decouple. Starting with real and universal boundary conditions at $`M_{Pl}`$, we solve the RGE displayed in Appendix A by expanding in the small parameters defined by the right-handed neutrino thresholds

$$t_3=\frac{1}{(4\pi )^2}\mathrm{ln}\frac{\mathrm{\Lambda }}{M_3}t_2=\frac{1}{(4\pi )^2}\mathrm{ln}\frac{M_3}{M_2}t_1=\frac{1}{(4\pi )^2}\mathrm{ln}\frac{M_2}{M_1},$$
(8)
with the ordering $`M_3>M_2>M_1`$. We thus obtain approximate analytic expressions for the radiatively induced $`\delta `$’s and $`m(A_{ii})`$ which depend on $`m_0`$, $`a_0`$ and the Yukawa couplings defined in the LF basis at the scale $`\mathrm{\Lambda }`$. Of course the latter can be immediately translated into the corresponding ones defined at any scale, e.g. $`M_1`$ or $`m_{susy}`$.

Precisely because of the lepton flavour and CP violating Yukawa couplings, the LF basis is continuously rotated and rephased with the RGE evolution and, as already discussed, we handle this by working in a rotating basis. It turns out that the basis transformation introduces negligible corrections in the FC terms (4) but important ones in the FV terms (5). Notice also that the seesaw effects stop at the decoupling threshold of the lightest right-handed neutrino, $`M_1`$, and that the RGE evolution of these effects down to the supersymmetric scales where CP and lepton FV transitions are estimated is generically small and can be neglect for our estimates. We now study separately the FC and the FV contributions.

### 3.1 Flavour conserving contribution

Starting in the LF basis at $`M_{Pl}`$ from $`A_e=y_ea_0`$, the seesaw interactions generate a $`m(A_{ii})`$ in the LF basis at $`M_1`$. At leading order in the $`t`$’s and defining $`y_\nu ^{}P_ay_\nu N_a`$ ($`a=1,2,3`$) and the right-handed neutrinos projectors $`P_1=\mathrm{diag}(1,0,0)`$, $`P_2=\mathrm{diag}(1,1,0)`$, $`P_3=𝕀`$, the latter reads 666Actually, (9) is an approximation by excess and improves as the involved Yukawa couplings become small. However, for our estimates it is reliable up to $`O(1)`$ Yukawa couplings.:

$$m(A_{ii})=8a_0y_{e_i}\left(t_2t_3m(N_2N_3)_{ii}+t_1t_3m(N_1N_3)_{ii}+t_1t_2m(N_1N_2)_{ii}\right),$$
(9)
where the Yukawas in the r.h.s. are evaluated at $`\mathrm{\Lambda }`$. A similar 777As far as the comparison with the expressions in Refs.  is possible, we find agreement with the last papers up to the coefficient of $`m(N_1N_2)_{ii}`$. formula were previously presented in Refs. , but with the various Yukawas involved defined at different scales. Then, one must be more cautious in drawing general conclusions and the authors validate theirs with specific numerical examples. Our approach offers the advantage that the corresponding results become more transparent and allow for an easier estimate of the effects once a pattern of seesaw parameters is assigned. Anyway, the crucial point  is that the more right-handed neutrinos are hierarchical, the more the FC contribution increases. Indeed, it vanishes in the limit that right-handed neutrinos are degenerate, in which case a contribution to $`m(A_{ii})`$ only appears at fourth order, proportional to $`m(y_eN[N,E]N)_{ii}`$. Notice also that the naive scaling relation is generally violated  by (9).

Eq. (9) displays a linear dependence on the unknown parameter $`a_0`$. By defining $`A_{ii}|A_{ii}|e^{i8\varphi _{A_i}}`$, then $`\varphi _{A_i}t_2t_3m(N_2N_3)_{ii}+`$ $`t_1t_3m(N_1N_3)_{ii}+`$ $`t_1t_2m(N_1N_2)_{ii}`$ is completely specified by the seesaw parameters. However, as appears from eq. (4) the experimental limit on $`d_i`$ doesn’t probe directly $`\varphi _{A_i}`$, rather it gives a bound on $`m(A_{ii})v_dm_im(a_i)`$ once supersymmetric masses are fixed (and up to unnatural conspiracies between the various amplitudes). Figs. 1, taken from Ref. , show the present upper limits on $`|m(a_e)|/m_R`$ and the planned ones for $`|m(a_\mu )|/m_R`$ in the plane $`(\stackrel{~}{M}_1,m_R)`$, respectively the bino and R slepton mass at $`m_{susy}`$. Notice that these limits are quite model independent because, apart from $`m_R`$ and $`\stackrel{~}{M}_1`$, there is only a mild dependence on $`m_L`$, which we have fixed as in mSUGRA for definiteness. Indeed, since the $`A`$-term amplitude in $`d_i^{FC}`$ arises from pure bino exchange, it does not involve 888The function $`I_B`$ in eq. (4) contains a factor $`|\mu |^2`$ that cancels the one in the overall coefficient. The exact expressions for $`I_B`$ can be found in , as well as approximate ones suitable for various pattern of supersymmetric masses. $`\mu `$ nor $`\mathrm{tan}\beta `$. In mSUGRA there is an unphysical region in the plane $`(\stackrel{~}{M}_1,m_R)`$ corresponding to $`m_0^2<0`$ and which has been indicated in light grey in the plots. Anyway, the dark grey region and below is also excluded because $`m_R\stackrel{~}{M}_1`$, in contrast with the requirement of neutrality for the LSP.

To find an upper estimate for the seesaw induced $`d_i^{FC}`$, let us evaluate $`|m(a_i)|/m_R`$ from eq. (9) by considering only its first term. This situation is representative because the terms proportional to $`t_1`$ are negligible when the lightest right-handed neutrino has smaller Yukawa couplings, as happens in many models and as one would guess from similarity with the charged fermion sectors. Anyway, the following discussion is trivially adapted to other cases. For definiteness, we adopt this set of reference threshold values: $`M_2=10^{12}`$ GeV, $`M_3=10^{15}`$ GeV, $`\mathrm{\Lambda }=210^{18}`$ GeV. Since, as demanded by perturbativity, $`(N_2N_3)_{ii}O(1)`$,

$$\frac{m(a_i)}{m_R}O\left(0.02\frac{t_2t_3}{210^3}\right)\frac{a_0}{m_R}.$$
(10)
This upper estimate is easily adapted to any given model once the relation between $`a_0`$ and $`m_R`$ is made explicit. To make the dependence more manifest, let us introduce the following two mSUGRA situations: a) $`a_0^2=2m_0^2`$; b) $`a_0^2=\stackrel{~}{M}_{1/2}^2`$, where $`\stackrel{~}{M}_{1/2}`$ is the universal gaugino mass at the gauge coupling unification scale. In the first case $`(a_0/m_R)^22(10.9r^2)`$, with $`r\stackrel{~}{M}_1/m_R`$. The ratio $`a_0/m_R`$ thus displays a mild excursion as it decreases from $`\sqrt{2}`$ down to $`0.45`$ when moving from the vertical left axes, where $`r1`$, to the joining line of the dark grey region, where $`r=1`$. For case b), the situation is opposite and not so mild: $`(a_0/m_R)2.5r`$. In the more natural cases in between, the ratio $`a_0/m_R`$ should thus be rather stable.

Let us now compare the upper estimate (10) with the experimental bound. As appears from eq. (4), the latter improves linearly with the experimental sensitivity to $`d_i`$. Taking $`a_0m_R`$, figs. (1) show that values around $`0.02`$ for $`|m(a_i)|/m_R`$ would require an experimental sensitivity to $`d_e`$ and $`d_\mu `$ respectively at the level of $`10^{28}10^{29}`$ e cm and $`210^{26}210^{27}`$ e cm, the exact value depending on the particular point of the $`(\stackrel{~}{M}_1,m_R)`$ plane. However, to avoid charge and color breaking, the constraint $`a_0/m_R3`$ in general applies. Thus, focusing for instance around $`m_R500`$ GeV, the FC contribution cannot exceed $`10^{28}`$ e cm for $`d_e`$ and $`210^{26}`$ e cm for $`d_\mu `$, even with highly hierarchical right-handed neutrinos. Notice that the former value could be at hand of future experimental searches for $`d_e`$, while the latter is at the very limit of the planned experimental sensitivity to $`d_\mu `$. Allowing for smaller $`m_R`$ values, $`m_R200`$ GeV, the FC seesaw upper estimate comes close to the present limit for $`d_e`$ while it cannot be more than $`210^{25}`$ e cm for $`d_\mu `$.

The result for $`d_\mu `$ is thus a kind of ”negative” one, since its eventual future discovery above $`210^{25}`$ e cm couldn’t be attributed to the radiative FC contribution of the seesaw. Even if the situation for $`d_e^{FC}`$ appears more optimistic, it is worth to underline that the previous upper estimate actually applies to models with at least four $`O(1)`$ neutrino Yukawa couplings and large CP phases, like those in Ref. . However, as we now turn to discuss, the FV contribution, underestimated by previous analyses, could drastically enhance the $`d_i`$ upper estimate.

### 3.2 Flavour violating contributions

In the following we isolate the potentially most important products of $`\delta `$’s and give an estimate of their relative magnitude with respect to the FC amplitude. The relevant approximations for the flavour violating elements, $`ij`$, of the $`\delta `$’s are:

$`m_Lm_R\delta _{ij}^{LR}`$ $`=`$ $`a_0m_i[2t_3𝒩_{ij}+{\displaystyle \underset{a}{}}{\displaystyle \frac{t_a^2}{2}}F_A(a,a)_{ij}+{\displaystyle \underset{a>b}{}}t_at_bF_A(a,b)_{ij}]`$ (11)
$`m_R^2\delta _{ij}^{RR}`$ $`=`$ $`{\displaystyle \underset{a}{}}{\displaystyle \frac{t_a^2}{2}}F_R(a,a)_{ij}+{\displaystyle \underset{a>b}{}}t_at_bF_R(a,b)_{ij}`$ (12)
$`m_L^2\delta _{ij}^{LL}`$ $`=`$ $`(6m_0^2+2a_0^2)t_3𝒩_{ij}+{\displaystyle \underset{a}{}}{\displaystyle \frac{t_a^2}{2}}F_L(a,a)_{ij}+{\displaystyle \underset{a>b}{}}t_at_bF_L(a,b)_{ij}`$ (13)
where $`a,b=1,2,3`$ , the matrix $`𝒩`$ is defined as $`t_3𝒩t_3N_3+t_2N_2+t_1N_1`$ and

$`F_A(a,b)`$ $`=`$ $`15\{E,N_a\}5\{E,N_b\}+12\{N_a,N_b\}+4[N_b,u_E^{(b)}]`$ (14)
$`+`$ $`2((N_{b}^{}{}_{(d)}{}^{}+D_e)N_a+(N_{a}^{}{}_{(d)}{}^{}+D_e)N_b)+4(\stackrel{~}{D}_\nu ^{(a)}/a_0+D_\nu ^{(a)})N_b`$
$`+`$ $`[N_a,E]+4[N_b,N_a]+7[E,N_b],`$
$`F_R(a,b)`$ $`=`$ $`8[(6m_0^2+4a_0^2)y_eN_ay_e^{}(6m_0^2+2a_0^2)y_eN_by_e^{}]`$ (15)
$`F_L(a,b)`$ $`=`$ $`2(6m_0^2+2a_0^2)(\{3N_b+E,N_a\}+2D_\nu ^{(a)}N_b+[N_a,u_E^{(b)}])+2m_{}^{2}{}_{H_u}{}^{}N_b`$ (16)
$`+`$ $`4a_0^2(\{3N_b+E,N_a\}+\{E,N_b\})+2(G_L+4a_0\stackrel{~}{D}_\nu ^{(a)})N_b`$
with the Yukawas in the r.h.s. evaluated at the scale $`\mathrm{\Lambda }`$. The subscript $`(d)`$ indicates to take only the diagonal elements of the matrix, $`u_E^{(a)}`$ is defined through $`[u_E^{(a)},E]=\{E,3E+N_a+D_e\}`$ and the definition of all the other quantities can be found in Appendix A.

By means of the above expressions, the predictions for the imaginary part of the various product of $`\delta `$’s can be studied. An imaginary part in the products $`(\delta ^{LL}\delta ^{LR})_{ii}`$, $`(\delta ^{LR}\delta ^{RR})_{ii}`$, …, only arises at third order in at least two different $`t`$’s.

#### 3.2.1 The LL RR contribution

Let us firstly discuss the LL RR double insertion, $`m(\delta ^{LL}\mu \eta _{\mathrm{}}m_{\mathrm{}}\delta ^{RR})_{ii}`$, which turns out to be the quantitatively most interesting one. Since in general $`\eta _{\mathrm{}}𝕀`$ and the phase of $`\mu `$ \- if any - is extremely small so that it can be safely neglected in the discussion, $`m(\delta ^{LL}\mu \eta _{\mathrm{}}m_{\mathrm{}}\delta ^{RR})_{ii}\mu m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{ii}`$. From eqs. (12, 13) it appears that $`\delta _{ij}^{RR}`$, $`\delta _{ij}^{LL}`$ are respectively of second and first order in the $`t`$’s. Thus, the lowest order is the cubic and the LL RR contribution reads:

$$m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{ii}=8m_i\frac{(6m_0^2+2a_0^2)(6m_0^2+3a_0^2)}{m_L^2m_R^2}\underset{a>b}{}t_at_b(t_a+t_b)m(N_aEN_b)_{ii}.$$
(17)
Notice that $`m(N_aEN_b)_{ii}=\mathrm{tan}^2\beta m(N_am_{\mathrm{}}^2N_b)_{ii}/m_t^2`$, where $`m_t`$ is the top mass. Due to the hierarchy in $`m_{\mathrm{}}`$, a potentially important effect could only come from $`m(N_{a}^{}{}_{i3}{}^{}N_{b}^{}{}_{3i}{}^{})`$. Notice also that, on the contrary of the FC contribution, (17) doesn’t strongly depend on $`a_0`$.

The relative importance between this amplitude and the FC one can be easily appreciated by neglecting the terms proportional to $`t_1`$, which, as already mentioned, is justified when the lightest right-handed neutrino has the smallest Yukawa couplings:

$$\frac{d_i^{FV_{LLRR}}}{d_i^{FC}}=\frac{\mu \mathrm{tan}^3\beta }{a_0}\frac{(6m_0^2+2a_0^2)(6m_0^2+3a_0^2)}{m_R^2m_L^2}(t_2+t_3)\frac{m(N_3m_{\mathrm{}}^2N_2)_{ii}}{m_t^2m(N_3N_2)_{ii}}\frac{I_B^{\prime \prime }}{I_B}.$$
(18)
For realistic values, the ratio of the two loop functions is slightly smaller than one. To obtain a rule of thumb, let us take $`m_0a_0m_{L,R}`$ and the reference ratio $`\mathrm{\Lambda }/M_2=210^6`$. Then, unless ad hoc fine-tunings in the structure of $`N`$, $`m(N_3m_{\mathrm{}}^2N_2)_{ii}m_\tau ^2m(N_3N_2)_{ii}`$, and one finds

$$\frac{d_i^{FV_{LLRR}}}{d_i^{FC}}0.5\frac{\mathrm{tan}^3\beta }{10^3}\frac{\mu }{\mu _{ew}}\frac{t_2+t_3}{0.1}(\mathrm{rule}\mathrm{of}\mathrm{thumb}),$$
(19)
where $`|\mu _{ew}|^20.5m_R^2+20\stackrel{~}{M}_1^2`$ is the value of $`\mu `$ accounting for radiative electroweak breaking in mSUGRA. Despite being of third order in the $`t`$’s, the FV amplitude can take over the FC one thanks to its $`\mathrm{tan}^3\beta `$ dependence. The precise value of this ratio is displayed in figs. 2 in the plane $`(\stackrel{~}{M}_1,m_R)`$. Cases a) and b) are separately displayed so that the behavior for any situation in between can be easily extrapolated. The relevant generalizations are also reminded.

The plots show that the rule of thumb is quite reliable and allow to extract, for each point of the plane, the value of $`\mathrm{tan}\beta `$ for which the FV amplitude takes over the FC one. In both cases a) and b), this happens for $`\mathrm{tan}\beta 10`$ \- the only exception being the region with $`r1`$ for case a), where $`\mathrm{tan}\beta 20`$ is required.

#### 3.2.2 The other contributions

The amplitudes with the other products of $`\delta `$’s are less important than the LL RR one. Consider for instance $`m(\delta ^{LR}\delta ^{RR})_{ii}`$. The flavour violating elements in $`\delta ^{RR}`$ and $`\delta ^{LR}`$ are respectively of second and first order, so that the product appears at third order:

$$m(\delta ^{LR}\delta ^{RR})_{ii}=8\frac{m_ia_0(12m_0^2+6a_0^2)}{m_R^3m_L}\underset{a>b}{}t_at_b(t_a+t_b)m(N_aEN_b)_{ii}.$$
(20)
It vanishes in the limit $`a_00`$, as the FC contribution. Neglecting the terms proportional to $`t_1`$, the ratio of the LR RR amplitude and the FC one reads

$$\frac{d_i^{FV_{LRRR}}}{d_i^{FC}}=\mathrm{tan}^2\beta \frac{12m_0^2+6a_0^2}{m_R^2}(t_2+t_3)\frac{m(N_3m_{\mathrm{}}^2N_2)_{ii}}{m_t^2m(N_3N_2)_{ii}}\frac{I_{B(R)}^{}}{I_B}.$$
(21)
It is easy to check that this ratio is smaller than one (unless ad hoc fine-tunings in the structure of $`N`$): the ratio of the two loop functions is slightly smaller than one for realistic values of the supersymmetric parameters and, taking for instance $`m_0=a_0`$ and $`\mathrm{\Lambda }/M_2=210^6`$, (21) should not exceed $`10^4\mathrm{tan}^2\beta `$.

For $`m(\delta ^{LL}\delta ^{LR})_{ii}`$, no imaginary part can arise at second order in $`t`$’s because both $`\delta ^{LL}`$ and $`\delta ^{LR}`$ are proportional to $`t_3𝒩`$. At third order there are many contributions and it is lengthy but straightforward to check that they are proportional to at least two different $`t`$’s. This contribution is also proportional to $`a_0`$ and could be comparable to the FC one but is in any case smaller than the LL RR one. The expression for the double LR insertion is also quite involved. It is proportional to $`a_0^2`$ and, being also suppressed by a factor $`m_{\mathrm{}}^2/m_L^2`$ with respect to the LL RR one, it can be safely neglected.

### 3.3 Predicted range and constrains on Yukawas

For values of $`\mathrm{tan}\beta 10`$, for which the EDMs are enhanced and thus most likely to be observed, the FV amplitude with the LL RR double insertion is generically dominant with respect to all other amplitudes. Then, when $`\mathrm{tan}\beta 10`$, all the considerations made previously for $`d_i^{FC}`$ actually apply to $`d_i`$ when strengthened by the rule of thumb factor (19). To give an example, if $`t_2+t_30.1`$ and $`\mu \mu _{ew}`$, $`d_e`$ cannot exceed $`0.510^{28(27)}\mathrm{tan}^3\beta /10^3`$ e cm and $`d_\mu `$ $`10^{26(25)}\mathrm{tan}^3\beta /10^3`$ e cm when $`m_R500(200)`$ GeV. Planned experimental sensitivities to $`d_\mu `$ could then test the seesaw radiative contribution for models with small $`m_R`$ and/or large $`\mathrm{tan}\beta `$, in which case the factor (19) could be up to $`5060`$ so that $`d_\mu `$ should not exceed $`O(10^{23})`$ e cm. On the contrary, the range of the seesaw induced $`d_e`$ already overlaps with the present experimental limit for values of $`m_R`$ up to $`1/2`$ TeV when $`\mathrm{tan}\beta 30`$.

Barring unnatural cancellations, planned searches for $`d_e`$ could thus test each term of the sum in (17), namely each 999Needless to say, this remains true even if the FC were the dominant amplitude, in which case the bound would be stronger.

$$t_at_b(t_a+t_b)m(N_a\frac{m_{\mathrm{}}^2}{m_\tau ^2}N_b)_{11}(a>b).$$
(22)
The effect of an eventual two orders of magnitude improvement for $`d_e`$ on the upper limit on $`m(N_am_{\mathrm{}}^2N_b)_{11}/m_\tau ^2`$ is displayed in fig. 3 by taking as reference values $`\mathrm{tan}\beta =30`$ and $`t_at_b(t_a+t_b)=210^4`$. It turns out that planned limits could be severe enough to test models with hierarchical neutrino Yukawa couplings. This cannot be done by present limits. For any given seesaw model, it is straightforward to extrapolate from the plot the level of experimental sensitivity required to test it. Notice also that for $`d_\mu `$, limits on $`m(N_am_{\mathrm{}}^2N_b)_{22}/m_\tau ^2`$ as strong as the present ones on $`m(N_am_{\mathrm{}}^2N_b)_{11}/m_\tau ^2`$ would require a sensitivity to $`d_\mu `$ at the level of $`210^{25}`$ e cm.

## 4 Lepton EDMs, Seesaw and $`SU(5)`$ Triplets

We now add to the supersymmetric seesaw model a stage of a minimal $`SU(5)`$ GUT above the gauge couplings unification scale, $`M_{GUT}210^{16}`$ GeV. Namely, we include the contribution of the higgs triplets Yukawa interactions to the RGE evolution of slepton masses from $`M_{Pl}`$ down to their threshold decoupling scale $`M_T`$, which is very likely to be bigger than $`M_3`$. Notations are defined in Appendix B.

It is not restrictive to work (at any scale) in the basis where $`(y_d^T=)y_e`$ and $`M_R`$ are real and diagonal and $`y_u=V^Td_u\varphi _uV`$, where $`d_u`$ are the (real and positive) eigenvalues of $`y_u`$, $`V`$ is the CKM matrix in the standard parameterization (more on this later) and $`\varphi _u`$ is a diagonal $`SU(3)`$ matrix. The RGE for the radiatively induced misalignments are written in eqs. (B.1) to (50). At first order in $`t_{1,2,3}`$ defined in eq. (8) and

$$t_T\frac{1}{(4\pi )^2}\mathrm{ln}\frac{M_{Pl}}{M_T},$$
(23)
their solutions at the scale $`M_1`$ reads:

$`m_R^2\delta ^{RR}`$ $`=`$ $`(6m_0^2+2a_0^2)(2t_1y_e^2+3t_TU^{})`$
$`m_L^2\delta ^{LL}`$ $`=`$ $`(6m_0^2+2a_0^2)((3t_T+t_1)y_e^2+t_3𝒩)`$ (24)
$`m_Lm_R\delta ^{RL}`$ $`=`$ $`a_0(6t_TU^{}m_{\mathrm{}}+2m_{\mathrm{}}t_3𝒩)`$
where the matrix $`𝒩`$ is defined as $`t_3𝒩t_3N_3+t_2N_2+t_1N_1`$. It is understood that all the quantities in the r.h.s. of (24) are evaluated at $`\mathrm{\Lambda }`$. The small effects of the subsequent evolution from $`M_1`$ down to $`m_{susy}`$ can be neglected in the following discussion. We explicitly write only the first order terms 101010Of course eqs. (24) overestimate the misalignment. However, for our estimates here they are reliable up to $`y_ty_{\nu _3}1`$. in the $`t`$’s because, contrarily to the seesaw case, they already produce a potential imaginary part for the FV contribution, eq. (5). Instead, for the diagonal part of $`A_e`$ this is not the case, as we now turn to discuss.

### 4.1 Flavour conserving contribution

As anticipated in the simplified discussion of Section 2, a potential candidate for $`m(A_{ii})`$, proportional to $`m(U^{}y_eN)_{ii}`$, could show up at $`O(t_T^2)`$. It is lengthy but straightforward to see that such term is exactly canceled by the effect of rotating the basis. For the same reason, in the general case with different thresholds $`M_{1,2,3}`$, $`M_T`$, the overall coefficient of $`m(U^{}y_e𝒩)_{ii}`$ is zero. Therefore, the second order contribution to $`m(A_{ii})`$ is just

$`m(A_{ii})=8a_0y_{e_i}[({\displaystyle \frac{3}{2}}t_T+t_3)t_2m(N_2N_3)_{ii}+({\displaystyle \frac{3}{2}}t_T+t_3)t_1m(N_1N_3)_{ii}`$
$`+t_2t_1m(N_1N_2)_{ii}],`$ (25)
namely the pure seesaw one discussed in the previous section, eq. (10), with the substitution $`t_3(3/2t_T+t_3)`$, due to the fact that above $`M_T`$ also the triplets circulate in the loop renormalizing the wave functions. $`t_T`$ is naturally expected to be small (triplets will not decouple much below $`10^{16}`$ GeV) so that eventual higher order contributions involving $`t_T`$ are expected to be negligible with respect to (25). As a result, a stage of $`SU(5)`$-like grand unification, cannot enhance by much the FC contribution with respect to the pure seesaw case.

### 4.2 Flavour violating contribution

On the contrary, products of two $`\delta `$’s have an imaginary part proportional to $`t_3t_T`$ $`m(U^{}y_e𝒩)_{ii}`$ and the FV contribution to $`d_i`$ is potentially bigger than the FC one. Most interestingly, the predicted range for the radiatively induced $`d_e`$ turns out to have been already sizeably excluded by the present experimental bounds. Planned searches for $`d_\mu `$ would also get close to test the range corresponding to radiatively induced misalignments.

Let us consider in turn the predictions for the imaginary part of the products of $`\delta `$’s, eq. (5), from their expressions given in eq. (24). As before, the most important contribution comes out from the LL RR double insertion. Since $`\eta _{\mathrm{}}𝕀`$ and the phase of $`\mu `$ \- if any - is experimentally small enough to be safely neglected in the present discussion, one has at the lowest relevant order in the $`t`$’s

$$m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{ii}=3t_Tt_3\frac{(6m_0^2+2a_0^2)^2}{m_R^2m_L^2}m(𝒩m_{\mathrm{}}U^{})_{ii}.$$
(26)
This FV contribution is potentially much bigger than the FC one because, apart from Yukawas and numerical coefficients, it is enhanced by a factor $`(m_\tau \mu \mathrm{tan}\beta )/(m_ia_0)`$.

The other FV combinations in (5) are:

$`m(\delta ^{LR}m_{\mathrm{}}\delta ^{LR})_{ii}`$ $`=`$ $`12t_Tt_3{\displaystyle \frac{a_0^2}{m_R^2}}{\displaystyle \frac{m(𝒩m_{\mathrm{}}^3U^{})_{ii}}{m_L^2}}`$ (27)
$`m(\delta ^{LR}\delta ^{RR})_{ii}`$ $`=`$ $`6t_Tt_3{\displaystyle \frac{a_0(6m_0^2+2a_0^2)}{m_R^3}}{\displaystyle \frac{m(𝒩m_{\mathrm{}}U^{})_{ii}}{m_L}}`$ (28)
$`m(\delta ^{LL}\delta ^{LR})_{ii}`$ $`=`$ $`{\displaystyle \frac{m_R^2}{m_L^2}}m(\delta ^{LR}\delta ^{RR})_{ii}.`$ (29)
The only contribution which does not vanish in the limit $`a_00`$ is the LL RR one. To compare the amplitudes, (26), (27) and (28), (29) have to be multiplied respectively by $`\mu \mathrm{tan}\beta `$ and $`m_Lm_R`$ and also by the appropriate loop functions, which have the same sign and in general are of the same order of magnitude (for more details see ). Then, if $`\mu >0`$ the four FV amplitudes have the same sign. However, the double LR insertion is always negligible with respect to the LL RR one because of the suppression factor $`m_{\mathrm{}}^2/m_L^2`$. For the other amplitudes (28),(29) the suppression factor with respect to (26) is $`m_{R,L}/(\mu \mathrm{tan}\beta )`$ (actually smaller due to the numerical coefficients). Then, even in the case of $`\mu <0`$, a reduction of the LL RR amplitude due to accidental cancellations seems unrealistic. Of course, also the contributions due to different thresholds in the right-handed neutrino spectrum are present, in exact analogy to what has been discussed in the previous section.

It is instructive to focus on the magnitude and dependencies of the combination $`m(𝒩m_{\mathrm{}}U^{})_{ii}`$. For $`d_e`$ and $`d_\mu `$, neglecting subleading terms proportional to $`y_c^2`$, $`y_u^2`$ and defining $`V_{td}|V_{td}|e^{i\varphi _{td}}`$:

$`m(𝒩m_{\mathrm{}}U^{})_{22}`$ $``$ $`m_\tau y_t^2V_{ts}\left(m(𝒩_{23}){\displaystyle \frac{m_e}{m_\tau }}|V_{td}|m(e^{i\varphi _{td}}𝒩_{12})\right)`$ (30)
$`m(𝒩m_{\mathrm{}}U^{})_{11}`$ $``$ $`m_\tau y_t^2|V_{td}|\left(m(e^{i\varphi _{td}}𝒩_{13})+{\displaystyle \frac{m_\mu }{m_\tau }}V_{ts}m(e^{i\varphi _{td}}𝒩_{12})\right)`$ (31)
where we exploited the fact that $`V_{ts}`$ is real in the standard parameterization. The latter is convenient to stress that the CP phases involved in the above combinations could be naturally large - as is indeed the case for $`\varphi _{td}`$ \- but, of course, any other choice must give equivalent results 111111Had we exploited the freedom of parameterizing $`V`$ in such a way that $`V_{ti}`$ are real numbers, then the CKM phase of $`V`$ would have been hidden in the redefinition of $`𝒩`$.. The contribution proportional to $`m(𝒩_{12})`$ has important suppression factors. Moreover $`|𝒩_{12}|`$ is independently constrained to be quite small from the present limits on $`\mu e\gamma `$. A plot of the present upper limit on $`|C_{12}|`$, with $`C(4\pi )^2t_3𝒩`$, in the plane $`(\stackrel{~}{M}_1,m_R)`$ can be found in Ref. . The limit were derived for the seesaw but also applies without significant modifications to the case of $`SU(5)`$ plus seesaw. As a result, once fixed $`M_T`$, experimental searches for $`d_e`$ and $`d_\mu `$ represent a test for $`m(e^{i\varphi _{td}}𝒩_{13})`$ and $`m(𝒩_{23})`$ respectively. Although present searches for $`\tau \mu \gamma `$ ($`\tau e\gamma `$) are not able by now to interestingly constrain $`|C_{23}|`$ ($`|C_{13}|`$), eventual experimental improvements would have an impact on $`d_\mu `$ ($`d_e`$) .

Notice also that the naive scaling relation is violated according to

$$\frac{d_e}{d_\mu }=\frac{|V_{td}|}{V_{ts}}\frac{m(e^{i\varphi _{td}}𝒩_{13})}{m(𝒩_{23})}$$
(32)
and that the combinations of Yukawas relevant for $`d_i`$ are independent on the phases of the diagonal $`SU(3)`$ matrix, $`\varphi _u`$. On the contrary, the latter affects (see for instance Ref. ) the proton decay lifetime due to d=5 operators, whose most important decay mode in the case of minimal $`SU(5)`$ is $`pK^+\overline{\nu }`$.

### 4.3 Predicted range and constraints on Yukawas

To understand how close to the experimental sensitivity is the radiatively induced EDM range, we plot in figs. 5 and 8 the upper estimate for the most important products of $`\delta `$’s. We consider a degenerate spectrum of right-handed neutrinos to pick out just the effect of the triplets, substitute the $`M_{Pl}`$-values for $`y_t`$ ($`.7`$) and the relevant CKM elements and choose as reference values $`\mathrm{\Lambda }/M_T=10^2`$ and $`\mathrm{\Lambda }/M_3=210^3`$. Then, the upper estimate follows by requiring perturbativity, $`m(N_{23})1`$, $`m(e^{i\varphi _{td}}N_{13})1`$. Solid and dashed lines refers to cases a) and b) respectively. We don’t show case b) for $`m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{ii}/m_\tau `$, because the predicted value is essentially flat. The upper estimate for $`m(\delta ^{LL}\delta ^{LR})_{ii}`$ is not shown, being closely related to that on $`m(\delta ^{LR}\delta ^{RR})_{ii}`$ (see eq. (29)). The bounds on $`m(\delta ^{LR}m_{\mathrm{}}\delta ^{LR})_{ii}/m_\tau `$ are not displayed because they are too small to be of potential interest.

For an easy comparison, we have taken from the sleptonarium  the experimental limits on the same quantities, figs. 4 and 7 . The experimental limits on $`m(\delta ^{LR}\delta ^{RR})_{ii}`$ are close to those on $`m(\delta ^{LL}\delta ^{LR})_{ii}`$ and mildly depend on $`m_L`$, which has been fixed as in mSUGRA in the plots. On the contrary, the limits on $`m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{ii}/m_\tau `$ are proportional to $`1/(\mu \mathrm{tan}\beta )`$. For definiteness, $`\mathrm{tan}\beta =10`$ and $`\mu =\mu _{ew}`$ have been assumed. The experimental bounds shown in figs. 4 and 7 correspond to the present bound $`d_e<10^{27}`$ e cm and to the planned sensitivity $`d_\mu <10^{24}`$ e cm. Since the experimental bounds are proportional to the bound on $`d_i`$, it is straightforward to extrapolate the sensitivity to $`d_\mu `$ and $`d_e`$ required to test the radiatively induced $`d_e`$ and $`d_\mu `$.

Let us firstly discuss the $`d_e`$ range. For our reference values, the upper estimate for $`m(\delta ^{LR}\delta ^{RR})_{11}`$ is $`O(10^6)`$ and the present bound on $`d_e`$ already constrains it to be smaller in a large region of the plane. Although we already know that they are not dominant, it is worth to discuss the LR RR amplitude and the similar LL LR one because, as already mentioned, they are quite model independent. Allowing for higher triplet masses, however, these amplitudes shift below the present experimental sensitivity. This is not the case for the LL RR amplitude. The maximum value allowed for $`m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{11}/m_\tau `$ is displayed in the right panel of fig. 5 for case a). In case b) it is $`10^3`$ everywhere. Then, the present experimental bound on $`d_e`$ has already explored the radiative range for roughly $`23`$ orders of magnitude.

This can be translated into an upper limit on $`m(e^{i\varphi _{td}}𝒩_{13})`$, as in fig. 6. Notice that this limit is indeed very strong and can be hardly evaded even looking for less favorable parameters than those taken as reference. To check this, keep e.g. $`t_3`$ fixed and try to worsen the limit: a reduction of $`\mu `$ with respect to $`\mu _{ew}`$ is unlikely to reduce it more than one order of magnitude; a suppression by another factor $`10`$ would require a $`\mathrm{\Lambda }/M_T1.6`$ ; on the other hand, for values of $`\mathrm{tan}\beta `$ larger than 10 the limit linearly improves. As a result, in the framework of the seesaw accompanied by $`SU(5)`$, the limit on $`m(e^{i\varphi _{td}}𝒩_{13})`$ can be considered robust.

Let us now turn to discuss $`d_\mu `$. By comparing figs. 7 and 8, it turns out that the upper bound on the LR RR insertion would require a sensitivity to $`d_\mu `$ at $`O(10^{26})`$ e cm, at the very limit of planned experimental improvements. Quantitatively, the upper estimate for $`m(\delta ^{LL}m_{\mathrm{}}\delta ^{RR})_{22}/m_\tau `$ which, for case b) is $`410^3`$, is more promising. The major part of the plane in fig. 8 could be tested with $`d_\mu `$ at the level of $`10^{24}`$$`10^{25}`$ e cm and in general $`d_\mu `$ should not exceed $`O(10^{23})`$ e cm. As a result, the eventual presence of triplets doesn’t enhance by much the range for $`d_\mu `$ with respect to the pure seesaw case. Nevertheless, here too the possibility of constraining $`m(𝒩_{23})`$ can be envisaged, as shown in fig. 9. Notice that, due to (32), a limit on $`m(𝒩_{23})`$ comparable to the present one on $`m(e^{i\varphi _{td}}𝒩_{13})`$ would require to improve the $`d_\mu `$ sensitivity down to $`510^{27}`$ e cm.

## 5 Conclusions

Planned experiments might significantly strengthen the limit on $`d_e`$  and $`d_\mu `$ . Their eventual discovery could be interpreted as an indirect manifestation of supersymmetry but could not reveal which source of CP (and possibly flavour) violation is actually responsible for the measured effect. Clearly, all sources in principle able to give the lepton EDM even at an higher level would be automatically constrained while those which fail in giving the lepton EDM at the desired level would be automatically excluded from the list of possible candidates.

In this work, we have estimated the ranges for the lepton EDMs induced by the Yukawa interactions of the heavy neutrinos, both alone and with the simultaneous presence of the heavy $`SU(5)`$ triplets. It turns out that the FV LLRR amplitude is in general larger or comparable to the FC one. So, EDMs are enhanced for large values of $`\mathrm{tan}\beta `$ and do not strongly depend on $`a_0`$.

The pure seesaw, even with large $`\mathrm{tan}\beta `$ and very hierarchical right-handed neutrinos, cannot account for $`d_\mu `$ above $`10^{23}`$ e cm. Its eventual discovery above this level would then signal the presence of some source of CP and LF violation other than the neutrino Yukawa couplings. The heavy triplets Yukawa couplings would be excluded from the list of possible sources because their additional presence do not significantly enhance the predicted range for $`d_\mu `$. Notice however that in the latter case a hierarchical right-handed neutrino spectrum is no more essential to end up with $`d_\mu `$ at an interesting level for planned searches. From the theoretical point of view, finding $`d_\mu `$ above $`10^{23}`$ e cm would indeed have a remarkable impact.

Interestingly enough, the present experimental sensitivity to $`d_e`$ is already testing the simultaneous presence of triplets and right-handed neutrinos. Correspondingly, constrains on $`m(e^{i\varphi _{td}}𝒩_{13})`$ have been derived which are significant even for quite large values of the triplet mass and unfavorable supersymmetric masses. Without the triplets, the radiatively-induced $`d_e`$ is close to the present experimental sensitivity only in models with large $`\mathrm{tan}\beta `$ and small slepton masses. Therefore, an experimental improvement would eventually provide interesting limits on the imaginary part of the relevant combination of neutrino Yukawa couplings and right-handed neutrino masses.

In the present discussion we have been looking for results as general as possible. Indeed, the specification of any particular seesaw model has been avoided and the attention has rather focused on the dependencies on the supersymmetric masses and heavy thresholds. Although some relevant seesaw models deserve a dedicated analysis , figs. 3, 6 and 9 are suitable for a quick check of the status of any given seesaw model with respect to the present and planned experimental limits on lepton EDMs.

## Acknowledgements

We thank C.A. Savoy for useful discussions and collaboration in the early stage of this work. I.M. acknowledge the CNRS and the “A. Della Riccia” Foundation for support and the SPhT, CEA-Saclay, for kind hospitality. We also acknowledge M. Peskin and Y. Farzan for pointing out a wrong numerical coefficient in the previous version of eq. (15); the impact of this correction is quantitatively negligible and the subsequent analysis and results stay unaffected.

## Appendix A Appendix: Seesaw

In the basis where charged fermion and right handed Majorana neutrino masses are diagonal

$$Wu^{cT}y_uQH_u+d^{cT}y_dQH_d+e^{cT}y_eLH_d+\nu ^{cT}y_\nu LH_u+\frac{1}{2}\nu ^{cT}M_R\nu ^c$$
(33)
where $`Q=\left(ud\right)^T`$, $`L=\left(\nu e\right)^T`$ and $`H_{d\left(u\right)}^0=v_{d\left(u\right)}`$. Soft scalar masses are defined as

$`_{soft}`$ $``$ $`\stackrel{~}{u}_R^{}m_u^2\stackrel{~}{u}_R+\stackrel{~}{d}_R^{}m_d^2\stackrel{~}{d}_R+\stackrel{~}{Q}^{}m_Q^2\stackrel{~}{Q}+\stackrel{~}{e}_R^{}m_e^2\stackrel{~}{e}_R+\stackrel{~}{L}^{}m_L^2\stackrel{~}{L}+\stackrel{~}{\nu }_R^{}m_\nu ^2\stackrel{~}{\nu }_R`$ (34)
$`+`$ $`(\stackrel{~}{u}_R^{}A_u\stackrel{~}{u}_Lv_u+\stackrel{~}{d}_R^{}A_d\stackrel{~}{d}_Lv_d+\stackrel{~}{e}_R^{}A_e\stackrel{~}{e}_Lv_d+\stackrel{~}{\nu }_R^{}A_\nu \stackrel{~}{\nu }_Lv_u+h.c.)`$
Let us also introduce the following notations:

$`y_x^{}y_xXy_xy_x^{}\stackrel{~}{X}\left(x=e,u,d\right)`$
$`P_ay_\nu y_\nu ^{\left(a\right)}y_\nu ^{}P_ay_\nu N_aP_ay_\nu y_\nu ^{}P_a\stackrel{~}{N}_a\left(a=1,2,3\right)`$
where $`P_2`$, $`P_1`$ project out $`M_3`$ and $`M_{3,2}`$ respectively, $`P_2=\mathrm{diag}(1,1,0)`$, $`P_1=\mathrm{diag}(1,0,0)`$, and $`P_3=𝕀`$.

### A.1 Running

Defining $`t\frac{1}{\left(4\pi \right)^2}\mathrm{ln}Q`$, the running of the Yukawa coupling constants is governed by:

$`{\displaystyle \frac{dy_\nu ^{\left(a\right)}}{dt}}`$ $`=`$ $`y_\nu ^{\left(a\right)}\left[3N_a+E+D_\nu ^{\left(a\right)}\right]`$
$`{\displaystyle \frac{dy_e}{dt}}`$ $`=`$ $`y_e\left[3E+N_a+D_e\right]`$
$`{\displaystyle \frac{dy_u}{dt}}`$ $`=`$ $`y_u\left[3U+D+D_u^{\left(a\right)}\right]`$ (35)
$`{\displaystyle \frac{dy_d}{dt}}`$ $`=`$ $`y_d\left[3D+U+D_d\right]`$
$`{\displaystyle \frac{d\left(P_aM_RP_a\right)}{dt}}`$ $`=`$ $`2\left[P_aM_R\stackrel{~}{N}_a^T+\stackrel{~}{N}_aM_RP_a\right]`$
where

$`D_\nu ^{\left(a\right)}`$ $`=`$ $`\left[Tr\left(3U+N_a\right)\left(3g_2^2+{\displaystyle \frac{3}{5}}g_1^2\right)\right]𝕀`$
$`D_e`$ $`=`$ $`\left[Tr\left(3D+E\right)\left(3g_2^2+{\displaystyle \frac{9}{5}}g_1^2\right)\right]𝕀`$
$`D_u^{\left(a\right)}`$ $`=`$ $`\left[Tr\left(3U+N_a\right)\left({\displaystyle \frac{16}{3}}g_3^2+3g_2^2+{\displaystyle \frac{13}{15}}g_1^2\right)\right]𝕀`$ (36)
$`D_d`$ $`=`$ $`\left[Tr\left(3D+E\right)\left({\displaystyle \frac{16}{3}}g_3^2+3g_2^2+{\displaystyle \frac{7}{15}}g_1^2\right)\right]𝕀.`$
For the trilinear couplings, defining $`P_aA_\nu A_\nu ^{\left(a\right)}`$

$`{\displaystyle \frac{dA_\nu ^{\left(a\right)}}{dt}}`$ $`=`$ $`4\stackrel{~}{N}_aA_\nu ^{\left(a\right)}+5A_\nu ^{\left(a\right)}N_a+2y_\nu ^{\left(a\right)}y_e^{}A_e+A_\nu ^{\left(a\right)}E+D_\nu ^{\left(a\right)}A_\nu ^{\left(a\right)}+2\stackrel{~}{D}_\nu ^{\left(a\right)}y_\nu ^{\left(a\right)}`$
$`{\displaystyle \frac{dA_e}{dt}}`$ $`=`$ $`4\stackrel{~}{E}A_e+5A_eE+2y_ey_{\nu }^{\left(a\right)}{}_{}{}^{}A_\nu ^{\left(a\right)}+A_eN_a+D_eA_e+2\stackrel{~}{D}_ey_e`$ (37)
$`{\displaystyle \frac{dA_u}{dt}}`$ $`=`$ $`4\stackrel{~}{U}A_u+5A_uU+2y_uy_d^{}A_d+A_uD+D_u^{\left(a\right)}A_u+2\stackrel{~}{D}_u^{\left(a\right)}y_u`$
$`{\displaystyle \frac{dA_d}{dt}}`$ $`=`$ $`4\stackrel{~}{D}A_d+5A_dD+2y_dy_u^{}A_u+A_dU+D_dA_d+2\stackrel{~}{D}_dy_d`$
where

$`\stackrel{~}{D}_\nu ^{\left(a\right)}`$ $`=`$ $`\left[Tr\left(3y_u^{}A_u+y_{\nu }^{\left(a\right)}{}_{}{}^{}A_\nu ^{\left(a\right)}\right)\left(3g_2^2\stackrel{~}{M}_2+{\displaystyle \frac{3}{5}}g_1^2\stackrel{~}{M}_1\right)\right]𝕀`$
$`\stackrel{~}{D}_e`$ $`=`$ $`\left[Tr\left(3y_d^{}A_d+y_e^{}A_e\right)\left(3g_2^2\stackrel{~}{M}_2+{\displaystyle \frac{9}{5}}g_1^2\stackrel{~}{M}_1\right)\right]𝕀`$ (38)
$`\stackrel{~}{D}_u^{\left(a\right)}`$ $`=`$ $`\left[Tr\left(3y_u^{}A_u+y_{\nu }^{\left(a\right)}{}_{}{}^{}A_\nu ^{\left(a\right)}\right)\left({\displaystyle \frac{16}{3}}g_3^2\stackrel{~}{M}_3+3g_2^2\stackrel{~}{M}_2+{\displaystyle \frac{13}{15}}g_1^2\stackrel{~}{M}_1\right)\right]𝕀`$
$`\stackrel{~}{D}_d`$ $`=`$ $`\left[Tr\left(3y_d^{}A_d+y_e^{}A_e\right)\left({\displaystyle \frac{16}{3}}g_3^2\stackrel{~}{M}_3+3g_2^2\stackrel{~}{M}_2+{\displaystyle \frac{7}{15}}g_1^2\stackrel{~}{M}_1\right)\right]𝕀.`$
For soft scalars, defining $`P_am_\nu ^2P_am_{\nu }^{2}{}_{}{}^{\left(a\right)}`$

$`{\displaystyle \frac{dm_L^2}{dt}}`$ $`=`$ $`\{m_L^2,E+N_a\}+2\left(y_e^{}m_e^2y_e+m_{H_d}^2E+A_e^{}A_e\right)+2\left(y_{\nu }^{\left(a\right)}{}_{}{}^{}m_{\nu }^{2}{}_{}{}^{\left(a\right)}y_\nu ^{\left(a\right)}+m_{H_u}^2N_a+A_{\nu }^{\left(a\right)}{}_{}{}^{}A_\nu ^{\left(a\right)}\right)+G_L`$
$`{\displaystyle \frac{dm_e^2}{dt}}`$ $`=`$ $`2\{m_e^2,\stackrel{~}{E}\}+4\left(y_em_L^2y_e^{}+m_{H_d}^2\stackrel{~}{E}+A_eA_e^{}\right)+G_e`$
$`{\displaystyle \frac{dm_{\nu }^{2}{}_{}{}^{\left(a\right)}}{dt}}`$ $`=`$ $`2\{m_{\nu }^{2}{}_{}{}^{\left(a\right)},\stackrel{~}{N}_a\}+4\left(y_\nu ^{\left(a\right)}m_L^2y_{\nu }^{\left(a\right)}{}_{}{}^{}+m_{H_u}^2\stackrel{~}{N}_a+A_\nu ^{\left(a\right)}A_{\nu }^{\left(a\right)}{}_{}{}^{}\right)`$ (39)
$`{\displaystyle \frac{dm_Q^2}{dt}}`$ $`=`$ $`\{m_Q^2,U+D\}+2\left(y_u^{}m_u^2y_u+m_{H_u}^2U+A_u^{}A_u\right)+2\left(y_d^{}m_d^2y_d+m_{H_d}^2D+A_d^{}A_d\right)+G_Q`$
$`{\displaystyle \frac{dm_u^2}{dt}}`$ $`=`$ $`2\{m_u^2,\stackrel{~}{U}\}+4\left(y_um_Q^2y_u^{}+m_{H_u}^2\stackrel{~}{U}+A_uA_u^{}\right)+G_u`$
$`{\displaystyle \frac{dm_d^2}{dt}}`$ $`=`$ $`2\{m_d^2,\stackrel{~}{D}\}+4\left(y_dm_Q^2y_d^{}+m_{H_d}^2\stackrel{~}{D}+A_dA_d^{}\right)+G_d`$
where

$`G_L=\left({\displaystyle \frac{6}{5}}g_1^2\stackrel{~}{M}_1^2+6g_2^2\stackrel{~}{M}_2^2\right)𝕀G_e=\left({\displaystyle \frac{24}{5}}g_1^2\stackrel{~}{M}_1^2\right)𝕀`$
$`G_Q=\left({\displaystyle \frac{2}{15}}g_1^2\stackrel{~}{M}_1^2+6g_2^2\stackrel{~}{M}_2^2+{\displaystyle \frac{32}{3}}g_3^2\stackrel{~}{M}_3^2\right)𝕀`$ (40)
$`G_u=\left({\displaystyle \frac{32}{15}}g_1^2\stackrel{~}{M}_1^2+{\displaystyle \frac{32}{3}}g_3^2\stackrel{~}{M}_3^2\right)𝕀G_d=\left({\displaystyle \frac{8}{15}}g_1^2\stackrel{~}{M}_1^2+{\displaystyle \frac{32}{3}}g_3^2\stackrel{~}{M}_3^2\right)𝕀.`$
Finally,

$`{\displaystyle \frac{dm_{H_d}^2}{dt}}`$ $`=`$ $`6Tr\left(y_dm_Q^2y_d^{}+y_d^{}m_d^2y_d+m_{H_d}^2D+A_d^{}A_d\right)+2Tr\left(y_em_L^2y_e^{}+y_e^{}m_e^2y_e+m_{H_d}^2E+A_e^{}A_e\right)+G_L`$
$`{\displaystyle \frac{dm_{H_u}^2}{dt}}`$ $`=`$ $`6Tr\left(y_um_Q^2y_u^{}+y_u^{}m_u^2y_u+m_{H_u}^2U+A_u^{}A_u\right)`$
$`+`$ $`2Tr\left(y_\nu ^{\left(a\right)}m_L^2y_{\nu }^{\left(a\right)}{}_{}{}^{}+y_{\nu }^{\left(a\right)}{}_{}{}^{}m_{\nu }^{2}{}_{}{}^{\left(a\right)}y_\nu ^{\left(a\right)}+m_{H_u}^2N_a+A_{\nu }^{\left(a\right)}{}_{}{}^{}A_\nu ^{\left(a\right)}\right)+G_L`$
and we define $`m_{}^{2}{}_{H_{d\left(u\right)}}{}^{}dm_{H_{d\left(u\right)}}^2/dt`$.

For energy scales between $`\mathrm{\Lambda }`$ and $`M_3`$, one has to take $`a=3`$; below $`M_3`$ and above $`M_2`$, $`a=2`$; while below $`M_2`$ and above $`M_1`$, $`a=1`$.

## Appendix B Appendix: $`SU(5)`$ \+ See-saw

We adopt the following notation to write matter and Higgs superfields:

$$\psi _{10}=\frac{1}{\sqrt{2}}\left(\begin{array}{ccccc}0& u_3^c& u_2^c& u^1& d^1\\ u_3^c& 0& u_1^c& u^2& d^2\\ u_2^c& u_1^c& 0& u^3& d^3\\ u^1& u^2& u^3& 0& e^c\\ d^1& d^2& d^3& e^c& 0\end{array}\right)\varphi _{\overline{5}}=\left(\begin{array}{c}d_1^c\\ d_2^c\\ d_3^c\\ e\\ \nu \end{array}\right)\eta _1=\nu ^c$$
(42)
$$\overline{H}=\left(\begin{array}{c}H_{3d}^{}{}_{1}{}^{}\\ H_{3d}^{}{}_{2}{}^{}\\ H_{3d}^{}{}_{3}{}^{}\\ H_{2d}^{}{}_{}{}^{}\\ H_{2d}^{}{}_{}{}^{0}\end{array}\right)H=\left(\begin{array}{c}H_{3u}^{}{}_{}{}^{1}\\ H_{3u}^{}{}_{}{}^{2}\\ H_{3u}^{}{}_{}{}^{3}\\ H_{2u}^{}{}_{}{}^{+}\\ H_{2u}^{}{}_{}{}^{0}\end{array}\right)$$
(43)
where $`<H_{2d\left(2u\right)}^0>v_{d\left(u\right)}`$. The superpotential

$$W\frac{1}{4}\psi ^{AB}y_u\psi ^{CD}H^Eϵ_{ABCDE}+\sqrt{2}\psi ^{AB}y_e\varphi _A\overline{H}_B+\eta _1y_\nu \varphi _AH^A+\frac{1}{2}\eta _1M_R\eta _1$$
(44)
where $`A,B,\mathrm{}=1,\mathrm{},5`$ and flavour indices are understood, gives rise to (33) with $`y_d=y_e^T`$ and $`y_u=y_u^T`$. For the soft breaking part of the Lagrangian

$`_{soft}`$ $`=`$ $`\stackrel{~}{\psi }^{}m_\psi ^2\stackrel{~}{\psi }+\stackrel{~}{\varphi }^{}m_\varphi ^2\stackrel{~}{\varphi }+\stackrel{~}{\eta }^{}m_\eta ^2\stackrel{~}{\eta }+m_h^2h^{}h+m_{\overline{h}}^2\overline{h}^{}\overline{h}+({\displaystyle \frac{1}{2}}M_5\lambda _5\lambda _5+\mathrm{h}.\mathrm{c}.)`$ (45)
$`+`$ $`\left(\stackrel{~}{u}^{cT}A_u\stackrel{~}{u}v_u+\stackrel{~}{d}^{cT}A_e^T\stackrel{~}{d}v_d+\stackrel{~}{e}^{cT}A_e\stackrel{~}{e}v_d+\stackrel{~}{\nu }^{cT}A_\nu \stackrel{~}{\nu }v_u\right)+\mathrm{h}.\mathrm{c}.`$
where the scalar fields are $`\stackrel{~}{\psi }=(\stackrel{~}{u}^c,\stackrel{~}{u},\stackrel{~}{d},\stackrel{~}{e}^c)`$, $`\stackrel{~}{\varphi }=(\stackrel{~}{d}^c,\stackrel{~}{e},\stackrel{~}{\nu })`$ and $`\stackrel{~}{\eta }=\stackrel{~}{\nu }^c`$ and gauginos are denoted with $`\lambda _5`$. In this way, in the scalar lepton mass matrix, $`m_{LL}^2=m_\varphi ^2`$, $`m_{RR}^2=m_\psi ^2`$.

### B.1 Running

Setting $`t\frac{1}{\left(4\pi \right)^2}\mathrm{ln}Q`$, $`dg_5/dt=3g_5^3`$, $`dM_5/dt=6g_5^2M_5`$ and

$`{\displaystyle \frac{dy_e}{dt}}`$ $`=`$ $`6y_eE+3\stackrel{~}{U}y_e+y_eN+G_ey_e`$
$`{\displaystyle \frac{dy_u}{dt}}`$ $`=`$ $`6y_uU+2\stackrel{~}{E}y_u+2y_uD+G_uy_u`$ (46)
$`{\displaystyle \frac{dy_\nu }{dt}}`$ $`=`$ $`6y_\nu N+4y_\nu E+G_\nu y_\nu `$
where $`G_e=84/5g_5^2+4Tr\left(E\right)`$, $`G_u=96/5g_5^2+Tr\left(3U+N\right)`$, $`G_\nu =48/5g_5^2+Tr\left(3U+N\right)`$. For scalar masses

$`{\displaystyle \frac{dm_\psi ^2}{dt}}`$ $`=`$ $`\{m_\psi ^2,2\stackrel{~}{E}+3\stackrel{~}{U}\}+4\left(y_em_\varphi ^2y_e^{}+m_{\overline{h}}^2\stackrel{~}{E}+A_eA_e^{}\right)+6\left(y_um_\psi ^2y_u^{}+m_h^2\stackrel{~}{U}+A_uA_u^{}\right)+G_\psi `$
$`{\displaystyle \frac{dm_\varphi ^2}{dt}}`$ $`=`$ $`\{m_\varphi ^2,4E+N\}+8\left(y_e^{}m_\psi ^2y_e+m_{\overline{h}}^2D+A_e^{}A_e\right)+2\left(y_\nu ^{}m_\eta ^2y_\nu +m_h^2N+A_\nu ^{}A_\nu \right)+G_\varphi `$
$`{\displaystyle \frac{dm_\eta ^2}{dt}}`$ $`=`$ $`5\{m_\eta ^2,\stackrel{~}{N}\}+10\left(y_\nu m_\varphi ^2y_\nu ^{}+m_h^2\stackrel{~}{N}+A_\nu A_\nu ^{}\right)`$ (47)
where $`G_\psi =144/5g_5^2M_5^2𝕀`$, $`G_\varphi =96/5g_5^2M_5^2𝕀`$. For trilinear couplings

$`{\displaystyle \frac{dA_e}{dt}}`$ $`=`$ $`10\stackrel{~}{E}A_e+3\stackrel{~}{U}A_e+8A_eE+6A_uy_u^{}y_e+A_eN+2y_ey_\nu ^{}A_\nu +G_eA_e+2\stackrel{~}{G}_ey_e`$ (48)
$`{\displaystyle \frac{dA_u}{dt}}`$ $`=`$ $`9\stackrel{~}{U}A_u+2\stackrel{~}{E}A_u+9A_uU+4A_ey_e^{}y_u+2A_uD+4y_uy_d^{}A_d+G_uA_u+2\stackrel{~}{G}_uy_u`$ (49)
$`{\displaystyle \frac{dA_\nu }{dt}}`$ $`=`$ $`7\stackrel{~}{N}A_\nu +11A_\nu N+4A_\nu E+8y_\nu y_e^{}A_e+G_\nu A_\nu +2\stackrel{~}{G}_\nu y_\nu `$ (50)
where $`\stackrel{~}{G}_e=84/5g_5^2M_5+4Tr\left(y_e^{}A_e\right)`$, $`\stackrel{~}{G}_u=96/5g_5^2M_5+Tr\left(3y_u^{}A_u+y_\nu ^{}A_\nu \right)`$, $`\stackrel{~}{G}_\nu =48/5g_5^2M_5+Tr\left(3y_u^{}A_u+y_\nu ^{}A_\nu \right)`$. Finally, for scalar higgses

$`{\displaystyle \frac{dm_h^2}{dt}}`$ $`=`$ $`Tr\left(6\stackrel{~}{U}+2\stackrel{~}{N}\right)m_h^2+6Tr\left(y_um_\psi ^2y_u^{}+y_u^{}m_\psi ^2y_u+A_uA_u^{}\right)`$ (51)
$`+`$ $`2Tr\left(y_\nu m_\psi ^2y_\nu ^{}+y_\nu ^Tm_\eta ^2y_\nu ^{}+A_\nu A_\nu ^{}\right){\displaystyle \frac{96}{5}}g_5^2M_5^2`$
$`{\displaystyle \frac{dm_{\overline{h}}^2}{dt}}`$ $`=`$ $`8Tr\left(\stackrel{~}{E}\right)m_{\overline{h}}^2+8Tr\left(y_em_\varphi ^2y_e^{}+y_e^{}m_\psi ^2y_e+A_eA_e^{}\right){\displaystyle \frac{96}{5}}g_5^2M_5^2.`$