text:
# Oscillating neutrinos and 𝜇→{𝑒,𝛾}

## 1 Introduction

In the pure Standard Model, flavour is exactly conserved in the leptonic sector since one can always choose a basis in which the (charged) lepton Yukawa matrix, $`𝐘_𝐞`$, and gauge interactions are flavour-diagonal. If neutrinos are massive and mixed, as suggested by the observation of atmospheric and solar fluxes , this is no longer true and there exists a source of lepton flavour violation (LFV), in analogy with the Kobayashi–Maskawa mechanism in the quark sector. Unfortunately, due to the smallness of the neutrinos masses, the predicted branching ratios for these processes are so tiny that they are completely unobservable, namely BR$`(\mu e\gamma )<10^{50}`$ .

In a supersymmetric (SUSY) framework the situation is completely different. Besides the previous mechanism, supersymmetry provides new direct sources of flavour violation in the leptonic sector, namely the possible presence of off-diagonal soft terms in the slepton mass matrices $`\left(m_L^2\right)_{ij}`$, $`\left(m_{e_R}^2\right)_{ij}`$, and trilinear couplings $`A_{ij}^e`$ . There are in the literature very strong bounds on these matrix elements coming from requiring branching ratios for LFV processes to be below the experimental rates . The strongest bounds come precisely from BR($`\mu e\gamma `$).

The most drastic way to avoid these dangerous off-diagonal terms is to impose perfect universality of the $`\left(m_L^2\right)_{ij}`$, $`\left(m_{e_R}^2\right)_{ij}`$, $`A_{ij}^e`$ matrices, i.e. to take them proportional to the unit matrix. Then, of course, they maintain this form in any basis, in particular in the basis where all the gauge interactions and $`𝐘_𝐞`$ are flavour-diagonal, so that we are left just with the previous tiny Kobayashi–Maskawa LFV effects. In principle, there is no theoretical reason to impose the universality constraint, and that is the so-called supersymmetric flavour problem. However, there are some particular supersymmetric scenarios, such as minimal supergravity, dilaton-dominated SUSY breaking or gauge-mediated SUSY breaking, that give rise to the desired universal soft terms (at least with much accuracy). In any case, the universality assumption is obviously the most conservative supposition that one can make about the form of the soft breaking terms when analyzing flavour violation effects; and this will be our starting point.

It turns out, however, that even under this extremely conservative assumption, if neutrinos are massive, radiative corrections may generate off-diagonal soft terms. The reason is that the flavour changing operators giving rise to the (non-diagonal) neutrino matrices will normally contribute to the renormalization group equations (RGEs) of the $`\left(m_L^2\right)_{ij}`$, $`\left(m_{e_R}^2\right)_{ij}`$ and $`A_{ij}^e`$ matrices, inducing off-diagonal entries.

The most interesting example of this occurs when neutrino masses are produced by a see-saw mechanism . Then, above a certain scale $`M`$, there are new degrees of freedom, the right-handed neutrinos, which have conventional Yukawa couplings with a matricial Yukawa coupling $`𝐘_\nu `$. Normally, $`𝐘_𝐞`$ and $`𝐘_\nu `$ cannot be diagonalized simultaneously. Working in the flavour basis where $`𝐘_𝐞`$ is diagonal, the off-diagonal entries of $`𝐘_\nu `$ (more precisely, the off-diagonal entries of $`𝐘_\nu ^+𝐘_\nu `$) drive off-diagonal entries in the previous $`\left(m_L^2\right)_{ij}`$, and $`A_{ij}^e`$ matrices through the RG running .

Actually, the supersymmetric see-saw is a extremely interesting scenario for a number of reasons. First, the see-saw mechanism is probably the most convincing and economical mechanism to naturally produce tiny neutrino masses (certainly it is the most popular one). Second, a non-supersymmetric see-saw suffers from a serious hierarchy problem. Perhaps this point has not been sufficiently acknowledged in the literature and we would like to emphasize it. It is common lore that the Standard Model (SM) presents a generic gauge hierarchy problem. Namely, the Higgs mass acquires quadratically-divergent radiative corrections, which, under the usual interpretation of the renormalization process, are naturally of the size of the scale at which new relevant physics enters. In a see-saw scenario the problem becomes more acute. The presence of very massive fermions (the right-handed neutrinos) coupled to the Higgs field produces logarithmically-divergent radiative corrections to the Higgs mass, which are proportional to $`M^2`$. This can be immediately noticed from the usual expression of the 1-loop effective potential, where the neutrino mass eigenstates which are essentially right-handed neutrinos contribute as $`\frac{1}{64\pi ^2}M^2Y_\nu ^2H^2(\mathrm{log}\frac{M^2}{\mu ^2}\frac{3}{2})+\mathrm{}`$ ($`\mu `$ is the renormalization scale)  . In the supersymmetric framework, this problem is automatically cured by the presence of the right-handed sneutrinos with a similar mass. Consequently, it is fair to say that supersymmetry is the most natural arena to implement the see-saw mechanism.

This is precisely the scenario that we will consider in this paper, computing BR($`\mu e\gamma `$) in the most general case compatible with the observational (atmospheric and solar) indications about neutrino masses and mixing angles. This subject has been addressed in previous works, especially in refs.  (other related works can be found in ref.). Let us briefly comment on the differences between the treatment of these works and the present one. The authors of ref. adopt a “top-down” point of view, i.e. they start with some sets of particular textures for $`𝐘_𝐞`$ and $`𝐘_\nu `$, coming from more or less well motivated symmetries, e.g. family symmetries, which can produce sensible low-energy neutrino mass matrices, $`_\nu `$, at the end of the day. Then, they study the corresponding predictions for BR($`\mu e\gamma `$). On the other hand, the authors of ref. adopt a “bottom-up” point of view, i.e. they start with the experimental data about $`_\nu `$ at low energy, finding out a form of $`𝐘_\nu `$ able to reproduce these data, and then computing the corresponding BR($`\mu e\gamma `$). Both approaches are interesting, but not general, as they are related to particular choices of $`𝐘_\nu `$. As a consequence, their results are not conclusive about what can we expect for BR($`\mu e\gamma `$) based on actual neutrino data. In our case, we also follow a “bottom-up” approach, i.e. we start with the existing experimental information about $`_\nu `$. Then we find the most general textures for $`𝐘_\nu `$ and $``$ (the Majorana mass matrix of the right-handed neutrinos) consistent with that information. In our opinion, this is by itself an interesting issue. The textures turn out to depend on some free parameters, but very few (as we shall see, the textures used in ref. correspond to a particular choice for these parameters). Next, we compute BR($`\mu e\gamma `$) and compare it with present and forthcoming experimental bounds. The approach can be considered as a generalization of that by Hisano et al. in ref.. As a matter of fact, we identify the dominant contributions, which were not considered in ref. as a consequence of the particular texture used there.

The paper is organized as follows. In sect. 2 we determine the most general form of $`𝐘_\nu `$ and $`𝐘_\nu ^+𝐘_\nu `$, compatible with all the phenomenological requirements. In several subsections we specialize the general formulas previously obtained for relevant physical scenarios. In sect. 3 we review the way in which the RG-running induces off-diagonal terms, and the relevance of the latter for $`l_il_j,\gamma `$ processes. In sect. 4 we apply the information of sects. 2, 3 to compute the predictions on BR($`l_il_j,\gamma `$), with particular focus on BR($`\mu e,\gamma `$). These predictions can be understood analytically using some approximations, for the sake of clarity in the discussion. Again, in several subsections we present and discuss the results for different physical scenarios. In sect. 5 we offer a summary of the predictions on BR($`l_il_j,\gamma `$) obtained in the previous sections. In sect. 6 we present our concluding remarks. Finally, we include an Appendix with the most important formulas that we have used for the full computations, specially RGEs and expressions for BR($`l_il_j,\gamma `$).

## 2 General textures reproducing experimental data

The supersymmetric version of the see-saw mechanism has a superpotential

$`W=W_0{\displaystyle \frac{1}{2}}\nu _{R}^{c}{}_{}{}^{T}\nu _R^c+\nu _{R}^{c}{}_{}{}^{T}𝐘_\nu LH_2,`$ (1)
where $`W_0`$ is the observable superpotential, except for neutrino masses, of the preferred version of the supersymmetric SM, e.g. the MSSM. The extra terms involve three additional neutrino chiral fields (one per generation; indices are suppressed) not charged under the SM group: $`\nu _{R}^{}{}_{i}{}^{}`$ ($`i=e,\mu ,\tau `$). $`𝐘_\nu `$ is the matrix of neutrino Yukawa couplings, $`L_i`$ ($`i=e,\mu ,\tau `$) are the left-handed lepton doublets and $`H_2`$ is the hypercharge $`+1/2`$ Higgs doublet. The Dirac mass matrix is given by $`𝐦_𝐃=𝐘_\nu H_2^0`$. Finally, $``$ is a $`3\times 3`$ Majorana mass matrix which does not break the SM gauge symmetry. It is natural to assume that the overall scale of $``$, which we will denote by $`M`$, is much larger than the electroweak scale or any soft mass. Below $`M`$ the theory is governed by an effective superpotential

$`W_{eff}=W_0+{\displaystyle \frac{1}{2}}(𝐘_\nu LH_2)^T^1(𝐘_\nu LH_2),`$ (2)
obtained by integrating out the heavy neutrino fields in (1). The corresponding effective Lagrangian contains a mass term for the left-handed neutrinos:

$`\delta ={\displaystyle \frac{1}{2}}\nu ^T_\nu \nu +\mathrm{h}.\mathrm{c}.,`$ (3)
with

$`_\nu =𝐦_{𝐃}^{}{}_{}{}^{T}^1𝐦_𝐃=𝐘_{\nu }^{}{}_{}{}^{T}^1𝐘_\nu H_2^0^2,`$ (4)
suppressed with respect to the typical fermion masses by the inverse power of the large scale $`M`$. It is convenient to extract the Higgs VEV by defining the $`\kappa `$ matrix as

$`\kappa =_\nu /H_2^0^2=𝐘_{\nu }^{}{}_{}{}^{T}^1𝐘_\nu ,`$ (5)
where $`H_2^0^2=v_2^2=v^2\mathrm{sin}^2\beta `$ and $`v=174`$ GeV. The experimental data about neutrino masses and mixings are referred to the $`_\nu `$ matrix, or equivalently $`\kappa `$, evaluated at low energy (electroweak scale). In this sense, it should be noted that eqs.(4, 5) are not defined at low energy but at the “Majorana scale”, $`M`$. Therefore, in order to compare to the experiment one has still to run $`\kappa `$ down to low energy through the corresponding RGE222 For hierarchical neutrinos, with which we will be mainly concerned, this running does not substantially affect the texture of $`\kappa `$, and thus the mixing angles. It does modify, however, the overall size of $`\kappa `$, which has to be taken into account.. Alternatively, following the bottom-up spirit, one can start with the physical $`\kappa `$ matrix at low energy and then run it upwards to the $`M`$ scale, where eqs.(4, 5) are defined. This has been the procedure we have followed in all the numerical calculations to be presented in the next sections.

Working in the flavour basis in which the charged-lepton Yukawa matrix, $`𝐘_𝐞`$, and gauge interactions are flavour-diagonal, the $`\kappa `$ matrix is diagonalized by the MNS  matrix $`U`$ according to

$$U^T\kappa U=\mathrm{d}iag(\kappa _1,\kappa _2,\kappa _3)D_\kappa ,$$
(6)
where $`U`$ is a unitary matrix that relates flavour to mass eigenstates

$`\left(\begin{array}{c}\nu _e\\ \nu _\mu \\ \nu _\tau \end{array}\right)=U\left(\begin{array}{c}\nu _1\\ \nu _2\\ \nu _3\end{array}\right).`$ (7)
It is possible, and sometimes convenient, to choose $`\kappa _i0`$. Then, $`U`$ can be written as

$`U=V\mathrm{d}iag(e^{i\varphi /2},e^{i\varphi ^{}/2},1),`$ (8)
where $`\varphi `$ and $`\varphi ^{}`$ are CP violating phases (if different from $`0`$ or $`\pi `$) and $`V`$ has the ordinary form of a CKM matrix

$$V=\left(\begin{array}{ccc}c_{13}c_{12}& c_{13}s_{12}& s_{13}e^{i\delta }\\ c_{23}s_{12}s_{23}s_{13}c_{12}e^{i\delta }& c_{23}c_{12}s_{23}s_{13}s_{12}e^{i\delta }& s_{23}c_{13}\\ s_{23}s_{12}c_{23}s_{13}c_{12}e^{i\delta }& s_{23}c_{12}c_{23}s_{13}s_{12}e^{i\delta }& c_{23}c_{13}\end{array}\right).$$
(9)
On the other hand, one can always choose to work in a basis of right neutrinos where $``$ is diagonal

$$=\mathrm{d}iag(_1,_2,_3)D_{},$$
(10)
with $`_i0`$. Then, from eqs.(5, 6)

$`D_\kappa =U^T𝐘_{\nu }^{}{}_{}{}^{T}D_^1𝐘_\nu U=U^T𝐘_{\nu }^{}{}_{}{}^{T}D_\sqrt{^1}D_\sqrt{^1}𝐘_\nu U`$ (11)
where, in an obvious notation, $`D_\sqrt{A}+\sqrt{D_A}`$. Multiplying both members of the eq.(11) by $`D_{\sqrt{\kappa ^1}}`$ from the left and from the right, we get

$`\mathrm{𝟏}=\left[D_\sqrt{^1}𝐘_\nu UD_{\sqrt{\kappa ^1}}\right]^T\left[D_\sqrt{^1}𝐘_\nu UD_{\sqrt{\kappa ^1}}\right]`$ (12)
whose solution is $`D_\sqrt{^1}𝐘_\nu UD_{\sqrt{\kappa ^1}}=R`$, with $`R`$ any orthogonal matrix ($`R`$ can be complex provided $`R^TR=\mathrm{𝟏}`$). Hence, in order to reproduce the physical, low-energy, parameters, i.e. the light neutrino masses (contained in $`D_\kappa `$) and mixing angles and CP phases (contained in $`U`$), the most general $`𝐘_\nu `$ matrix is given by

$`𝐘_\nu =D_{\sqrt{}}RD_\sqrt{\kappa }U^+`$ (13)
So, besides the physical and measurable low-energy parameters contained in $`D_\kappa `$ and $`U`$, $`𝐘_\nu `$ depends on the three (unknown) positive mass eigenvalues of the righthanded neutrinos and on the three (unknown) complex parameters defining $`R`$. We will see, however, that in practical cases the number of relevant free parameters becomes drastically reduced. Let us stress that eq.(13), which is central for us, has to be understood at the $`M`$ scale and in the above-defined basis, i.e. the flavour basis in which the charged-lepton Yukawa matrix, $`𝐘_𝐞`$, the right-handed Majorana mass matrix, $``$, and the gauge interactions are flavour-diagonal. The extension of eq.(13) to other choices of basis is straightforward, taking into account the transformation properties of the $`𝐘_𝐞`$, $`𝐘_\nu `$, $``$ matrices under changes of basis. For example, if one works in a basis of $`\nu _R`$ where the $``$ matrix is non-diagonal, then the most general $`𝐘_\nu `$ matrix reads $`𝐘_\nu =U_M^{}D_{\sqrt{}}RD_\sqrt{\kappa }U^+`$, where $`U_M`$ is such that $`=U_M^{}D_{}U_M^+`$. Let us also mention that eq.(13) is valid as well for the non-supersymmetric see-saw mechanism.

As commented in sect. 1, in ref. a particular choice for $`𝐘_\nu `$ was taken, which corresponds to take $`R=\mathrm{𝟏}`$ in eq.(13). This is equivalent to assume that there exists a basis of $`L_i`$ and $`\nu _{R}^{}{}_{i}{}^{}`$ in which $`𝐘_\nu `$ and $``$ are simultaneously diagonal (though not $`𝐘_𝐞`$). To see this notice that if $`R=\mathrm{𝟏}`$, then $`𝐘_\nu `$, as given by eq.(13), can be made diagonal by rotating $`L_i`$ with the $`U`$-matrix. This hypothesis may be consistent with certain models, namely when all the leptonic flavour violation can be attributed to the sector of charged leptons, but clearly it is not general, and it is inconsistent with other scenarios.

Another special choice of $`R`$ occurs when the $`𝐘_\nu `$ matrix given by eq.(13) has the form $`𝐘_\nu =WD_Y`$, where $`D_Y`$ is a diagonal matrix and $`W`$ is a unitary matrix. From the see-saw formula (5), $`W`$ is related to $`D_{}`$ and $`D_Y`$ through $`D_\sqrt{^1}=W^{}(D_{Y^1}\kappa D_{Y^1})W^+`$. Then, by rotating $`\nu _{R}^{}{}_{i}{}^{}`$ with the $`W`$-matrix, $`𝐘_𝐞`$ and $`𝐘_\nu `$ get simultaneously diagonal, while $``$ gets non-diagonal, more precisely $`=D_Y\kappa ^1D_Y`$. In other words, all the leptonic flavour violation can be attributed to the sector of right-handed neutrinos, in contrast to the previous situation. It is important to remark that both are very special situations, not generic at all.

On the other hand, as we will see in the next section, the rate for $`l_il_j,\gamma `$ processes just depends on the $`𝐘_\nu ^+𝐘_\nu `$ matrix (to be precise, on the matrix element $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$) rather than on $`𝐘_\nu `$, since this is the quantity that enters the RGEs of the off-diagonal soft parameters. From eq.(13) $`𝐘_\nu ^+𝐘_\nu `$ has the general form

$`𝐘_\nu ^+𝐘_\nu =UD_\sqrt{\kappa }R^+D_{}RD_\sqrt{\kappa }U^+`$ (14)
Notice that $`𝐘_\nu ^+𝐘_\nu `$, and therefore eq.(14), do not depend on the $`\nu _R`$-basis, and thus on the fact that $``$ is diagonal or not.

In order to illustrate the use of eqs.(13, 14), let us consider some interesting scenarios that often appear in the literature

### 2.1 $`\nu _L`$’s and $`\nu _R`$’s completely hierarchical

In this case $`\kappa _1\kappa _2\kappa _3`$ and $`_1_2_3`$, so that $`D_\kappa `$ and $`D_{}`$ can be approximated by

$$D_\kappa \mathrm{d}iag(0,\kappa _2,\kappa _3),D_{}\mathrm{d}iag(0,0,_3)$$
(15)
where $`(v_2^2\kappa _2)^2(\mathrm{\Delta }m_\nu ^2)_{sol}`$, $`(v_2^2\kappa _3)^2(\mathrm{\Delta }m_\nu ^2)_{atm}`$. Then, from eq.(13),

$`(𝐘_\nu )_{ij}\sqrt{_3}\delta _{i3}R_{3l}\sqrt{\kappa _l}U_{lj}^+`$ (16)
So, for a certain set of low-energy physical quantities (given by $`\kappa _l`$ and $`U`$), $`𝐘_\nu `$ just depends on two (generally complex) parameters, $`\sqrt{_3}R_{32}`$ and $`\sqrt{_3}R_{33}`$. Let us remark that, although the $`(𝐘_\nu )_{1j}`$ and $`(𝐘_\nu )_{2j}`$ entries are proportional to $`\sqrt{_1}`$ and $`\sqrt{_2}`$ respectively, and thus suppressed, they are normally relevant if one wishes to reconstruct the neutrino mass matrix, $`\kappa `$, from $`𝐘_\nu `$, $`D_{}`$ through the see-saw equation (5). Still, they are not relevant for the $`𝐘_\nu ^+𝐘_\nu `$ entries, and thus for the predictions on BR($`l_il_j,\gamma `$), which is our main goal.

It is interesting to notice that this scenario also includes the possibility of an “inverse hierarchy”, i.e. $`\mathrm{d}iag(_3,0,0)`$, or any other ordering. The reason is that one can permute the entries of $`D_{}`$ by redefining the $`\nu _R`$’s through an appropriate rotation. In particular, one can pass from an inverse hierarchy to an ordinary one through a rotation in the 1–3 “plane”: $`\nu _R=R_p\nu _R^{}`$ with $`R_p=\{(0,0,1),(0,1,0),(1,0,0)\}`$. Then, $`𝐘_\nu 𝐘_{\nu }^{}{}_{}{}^{}=D_{\sqrt{}}^{}R_p^TRD_\sqrt{\kappa }U^+`$ with $`D_{}^{}=R_p^TD_{}R_p`$. Hence, the general formula (13) and its particularization (16) remain identical with the modification $`RR^{}=R_p^TR`$. Since $`R`$ runs over all possible orthogonal matrices, so $`R^{}`$ does. So, both starting points (inverse or ordinary hierarchy) lead to the same complete results.

In many models one can impose an extra condition coming from the unification of Yukawa couplings. In particular, in simple $`SO(10)`$ models, usually one of the $`𝐘_\nu ^+𝐘_\nu `$ eigenvalues (the largest one) coincides with the top Yukawa coupling $`|Y_t|^2`$ at the unification scale, $`M_X`$ (in other more elaborate models, it coincides with $`3|Y_t|^2`$) . We will denote $`|Y_0|^2`$ this maximum eigenvalue. Now, from eq.(16), the largest eigenvalue of $`𝐘_\nu ^+𝐘_\nu `$ is $`|Y_0|^2_3\left(|R_{32}|^2\kappa _2+|R_{33}|^2\kappa _3\right)`$, which (after running up to $`M_X`$) would be identified with $`|Y_t(M_X)|^2`$. This eliminates one real parameter, leaving just $`R_{33}/R_{32}`$, plus one phase, e.g. $`\mathrm{p}hase(R_{33})`$, as the independent parameters. Actually, as mentioned above, the $`l_il_j,\gamma `$ processes just depends on the $`𝐘_\nu ^+𝐘_\nu `$ matrix, which does not depend on the previous phase. Therefore, one is just left with one (complex) parameter, $`R_{33}/R_{32}`$.

A convenient way to visualize this is to take the following parametrization333The parametrization (17) does not count all the possible forms of $`R`$, as it does not include “reflections”. However, for the scenario considered in this subsection it is general enough, since the uncounted forms of $`R`$ do not produce different textures of $`𝐘_\nu `$. of $`R`$

$`R=\left(\begin{array}{ccc}\widehat{c}_2\widehat{c}_3& \widehat{c}_1\widehat{s}_3\widehat{s}_1\widehat{s}_2\widehat{c}_3& \widehat{s}_1\widehat{s}_3\widehat{c}_1\widehat{s}_2\widehat{c}_3\\ \widehat{c}_2\widehat{s}_3& \widehat{c}_1\widehat{c}_3\widehat{s}_1\widehat{s}_2\widehat{s}_3& \widehat{s}_1\widehat{c}_3\widehat{c}_1\widehat{s}_2\widehat{s}_3\\ \widehat{s}_2& \widehat{s}_1\widehat{c}_2& \widehat{c}_1\widehat{c}_2\end{array}\right),`$ (17)
where $`\widehat{\theta }_1`$, $`\widehat{\theta }_2`$, $`\widehat{\theta }_3`$ are arbitrary complex angles. Then, from eq.(16)

$`𝐘_\nu `$ $``$ $`{\displaystyle \frac{|Y_0|}{\sqrt{(|\widehat{s}_1|^2\kappa _2+|\widehat{c}_1|^2\kappa _3}}}{\displaystyle \frac{\widehat{c}_2}{|\widehat{c}_2|}}`$ (18)
$`\times `$ $`\left(\begin{array}{ccc}\widehat{0}& 0& 0\\ 0& 0& 0\\ \widehat{s}_1\sqrt{\kappa _2}U_{12}^{}+\widehat{c}_1\sqrt{\kappa _3}U_{13}^{}& \widehat{s}_1\sqrt{\kappa _2}U_{22}^{}+\widehat{c}_1\sqrt{\kappa _3}U_{23}^{}& \widehat{s}_1\sqrt{\kappa _2}U_{32}^{}+\widehat{c}_1\sqrt{\kappa _3}U_{33}^{}\end{array}\right)`$
which depends just on the angle $`\widehat{\theta }_1`$ and the phase $`\frac{\widehat{c}_2}{|\widehat{c}_2|}`$; $`𝐘_\nu ^+𝐘_\nu `$ depends just on $`\widehat{\theta }_1`$. The previous equation should be understood at the $``$ scale; the prefactor $`|Y_0|`$ represents the square root of the maximal eigenvalue of $`𝐘_\nu ^+𝐘_\nu `$ at the $``$-scale.

In the typical case in which everything is real, very used in the literature, the $`\widehat{\theta }_i`$ angles are real, and hence $`𝐘_\nu `$ and $`𝐘_\nu ^+𝐘_\nu `$ depend on a single real parameter, $`\widehat{\theta }_1`$. Therefore, in spite of the additional degrees of freedom introduced by the lack of knowledge about $`R`$, this scenario is quite predictive.

In the previous discussion, it has been assumed that at least one of the two entries, $`R_{32}`$ or $`R_{33}`$, is different from zero. However, it may happen that $`R_{32}=R_{33}=0`$. Then, $`\kappa _1,_1,_2`$ cannot be neglected anymore in the general equation eq.(13). It is worth-noticing that going to the basis where the $`D_{}`$ matrix is ‘inverse hierarchical’, i.e. $`\nu _RR_p\nu _R`$, $`D_{}\mathrm{d}iag(_3,_2,_1)`$, this possibility translates into $`R_{12}=R_{13}=0`$, which includes the particularly simple case $`R=\mathrm{𝟏}`$.

### 2.2 $`\nu _L`$’s hierarchical and $`\nu _R`$’s degenerate

In this case, $`D_\kappa `$ and $`D_{}`$ have the form

$$D_\kappa \mathrm{d}iag(0,\kappa _2,\kappa _3),D_{}\mathrm{d}iag(,,)$$
(19)
where, again, $`(v_2^2\kappa _2)^2(\mathrm{\Delta }m^2)_{sol}`$, $`(v_2^2\kappa _3)^2(\mathrm{\Delta }m^2)_{atm}`$. Then, from eq.(13),

$`(𝐘_\nu )_{ij}\sqrt{}R_{il}\sqrt{\kappa _l}U_{lj}^+`$ (20)
So, in principle, for a certain set of low-energy physical quantities (given by $`\kappa _q`$ and $`U`$), $`𝐘_\nu `$ depends on one real parameter ($``$) and three complex ones defining the $`R`$ matrix. Notice, however, that in the real case, where $`R^+=R^T`$, $`𝐘_\nu ^+𝐘_\nu `$ takes the simple form

$`𝐘_\nu ^+𝐘_\nu =UD_\kappa U^+,`$ (21)
which just contains one free parameter, $``$. This can be fixed by a unification condition, similarly to the previous case. Here $`|Y_0|^2\kappa _3|Y_t|^2`$ (where the identification must be understood at the $`M_X`$ scale). Therefore, this is also a very predictive scenario. Notice that here $`𝐘_\nu ^+𝐘_\nu `$ does not depend on the form of $`R`$, so in this case the choice $`R=\mathrm{𝟏}`$ is indeed general enough to study $`l_il_j,\gamma `$ processes.

### 2.3 $`\nu _L`$’s quasi-degenerate

In this case $`D_\kappa \mathrm{d}iag(\kappa _1,\kappa _2,\kappa _3)`$, with $`\kappa _1\kappa _2\kappa _3\kappa `$. Then, it is logical to assume that $``$ has degenerate eigenvalues, otherwise a big conspiracy would be needed between $`𝐘_\nu `$ and $``$ so that $`\kappa _1\kappa _2\kappa _3`$ in eq.(5). Hence

$$D_\kappa \mathrm{d}iag(\kappa ,\kappa ,\kappa ),D_{}\mathrm{d}iag(,,).$$
(22)
Consequently, eq.(20) and eq.(21) (for $`R`$ real) hold in this case too. Since $`D_\kappa \kappa \mathrm{𝟏}`$, from eq.(21) we notice

$`𝐘_\nu ^+𝐘_\nu =UD_\kappa U^+|Y_0|^2\mathrm{𝟏},`$ (23)
where $`|Y_0|^2=\kappa _3`$ is the largest eigenvalue of the $`𝐘_\nu ^+𝐘_\nu `$ matrix. This means that for quasi-degenerate neutrinos, with $`R`$ real, we expect small off-diagonal entries in $`𝐘_\nu ^+𝐘_\nu `$ and thus suppressed $`l_il_j,\gamma `$ rates. We will be more quantitative about the size of those entries in subsect. 4.3.

The scenario of partially degenerate neutrinos will also be addressed in that section.

## 3 RG-induced lepton flavour violating soft terms and $`l_il_j,\gamma `$

In addition to the supersymmetric Lagrangian, which contains the gauge interactions plus the part derived from the superpotential (1), one has to consider the soft breaking terms (gaugino and scalar masses, and trilinear and bilinear scalar terms) coming from the (unknown) supersymmetry breaking mechanism. In a self-explanatory notation, they have the form

$`_{\mathrm{soft}}=\left(m_L^2\right)_{ij}\overline{L}_iL_j+\left(m_{e_R}^2\right)_{ij}\overline{e}_{Ri}e_{Rj}+(A_{e}^{}{}_{ij}{}^{}e_{R}^{c}{}_{i}{}^{}H_1L_j+\mathrm{h}.\mathrm{c}.)+\mathrm{etc}.,`$ (24)
where we have written explicitly just the soft breaking terms in the leptonic sector, namely scalar masses and trilinear scalar terms. All the fields in the previous equation denote just the corresponding scalar components. As explained in sect. 1, concerning flavour violation the most conservative starting point for $`_{\mathrm{soft}}`$ is the assumption of universality, which corresponds to take

$`\left(m_L^2\right)_{ij}=m_0^2\mathbf{1},\left(m_{e_R}^2\right)_{ij}=m_0^2\mathbf{1},A_{e}^{}{}_{ij}{}^{}=A_0𝐘_{𝐞}^{}{}_{ij}{}^{},`$ (25)
so that working in the $`L_i`$ and $`e_{R}^{}{}_{i}{}^{}`$ basis where $`𝐘_𝐞`$ is diagonal, the soft terms do not contain off-diagonal (lepton flavour violating) entries.

As mentioned in sect. 1, this is only strictly true at the scale where universality is imposed, e.g. $`M_X`$ in GUT models. Below that scale, the RGEs of the soft terms, which contain non-diagonal contributions proportional to $`𝐘_\nu ^+𝐘_\nu `$, induce off-diagonal soft terms444We will pay attention just to the contribution that appears in the MSSM, i.e. the one $`𝐘_\nu ^+𝐘_\nu `$. This is a minimal possibility. In more complicated models, there can be additional sources of flavour violation in the RGEs, leading to larger predictions for BR($`l_il_j,\gamma `$).. These contributions are decoupled at the characteristic scale of the right-handed neutrinos, $`M`$. As noted e.g. in ref. , and as it is clear from eqs.(6771) in the Appendix, this mechanism is not efficient to generate off-diagonal entries in the $`\left(m_{e_R}^2\right)`$ matrix, since the corresponding RGE does not contain $`𝐘_\nu ^+𝐘_\nu `$ terms. More precisely, in the leading-log approximation555We use the leading-log approximation through the text in order to make the results easily understandable. Nevertheless, the numerical results, to be exposed below, have been obtained by integrating the full set of RGEs., the off-diagonal soft terms at low-energy are given by

$`\left(m_L^2\right)_{ij}`$ $``$ $`{\displaystyle \frac{1}{8\pi ^2}}(3m_0^2+A_0^2)(𝐘_\nu ^+𝐘_\nu )_{ij}\mathrm{log}{\displaystyle \frac{M_X}{M}},`$
$`\left(m_{e_R}^2\right)_{ij}`$ $``$ $`0,`$
$`(A_e)_{ij}`$ $``$ $`{\displaystyle \frac{3}{8\pi ^2}}A_0Y_{l_i}(𝐘_\nu ^+𝐘_\nu )_{ij}\mathrm{log}{\displaystyle \frac{M_X}{M}},`$ (26)
where $`ij`$ and $`Y_{l_i}`$ is the Yukawa coupling of the charged lepton $`l_i`$.

The amplitude for a $`l_il_j,\gamma `$ process has the general form

$`T=ϵ^\alpha \overline{l_j}m_{l}^{}{}_{i}{}^{}i\sigma _{\alpha \beta }q^\beta \left(A_LP_L+A_RP_R\right)l_i,`$ (27)
where $`q`$ is the momentum of the photon, $`P_{RL}=(1\pm \gamma _5)/2`$ and $`A_L`$ ($`A_R`$) is the coefficient of the amplitude when the incoming $`l_i`$ lepton is left (right), and thus the $`l_j`$ lepton is right (left). The expressions for $`A_L`$, $`A_R`$ can be found in the literature  (see Appendix for more details). The corresponding branching ratio is given by

$`\mathrm{BR}(l_il_j,\gamma )={\displaystyle \frac{12\pi ^2}{G_F^2}}\left(|A_L|^2+|A_R|^2\right)`$ (28)
Notice that the $`l_il_j,\gamma `$ process violates the chiral symmetries $`l_{R}^{}{}_{k}{}^{}e^{i\alpha _k}l_{R}^{}{}_{k}{}^{}`$. Under the assumption of initial universality, this symmetry is violated just by the Yukawa coupling of the corresponding lepton. Thus, $`A_L`$ ($`A_R`$) is proportional to the mass of the righthanded lepton involved, i.e. $`l_{j}^{}{}_{R}{}^{}`$ ($`l_{i}^{}{}_{R}{}^{}`$). Since $`m_{l}^{}{}_{j}{}^{2}m_{l}^{}{}_{i}{}^{2}`$, then $`|A_L|^2|A_R|^2`$. This can also be checked by inspection of the diagrams that contribute to $`l_il_j,\gamma `$, taking into account eq.(3). It is interesting to note that the angular distribution of $`l_j`$ with respect to the polarization direction of the incoming lepton ($`P_{l_i}`$) depends on the relative magnitudes of the coefficients $`A_L`$, $`A_R`$. The relation $`|A_L|^2|A_R|^2`$ predicted by the see-saw mechanism, and which does not hold in general in other models of generation of neutrino masses, gives a characteristic ($`1P_{l_i}\mathrm{cos}\theta `$) distribution that could be measured by future experiments .

We have taken into account all contributions to BR($`l_il_j,\gamma `$), using the general expressions given in the literature, in particular from ref. , as explained in the Appendix. These have not been obtained by using the mass-insertion approximation, but by diagonalizing all the mass matrices involved in the task, i.e. those of (left and right) sleptons, charginos and neutralinos. The diagrams have the form shown in Fig. 1. The precise form of BR($`l_il_j,\gamma `$) that we have used in our computations is a rather cumbersome expression, given in the Appendix. However, for the sake of the physical discussion it is interesting to think in the mass-insertion approximation to identify the dominant contributions. As discussed in ref. , these correspond to the mass-insertion diagrams enhanced by $`\mathrm{tan}\beta `$ factors. All of them are proportional to $`m_{L}^{2}{}_{ij}{}^{}`$, and have the generic form shown in Fig. 2. So, in all cases

$`\mathrm{BR}(l_il_j,\gamma ){\displaystyle \frac{12\pi ^2}{G_F^2}}|A_R|^2{\displaystyle \frac{\alpha ^3}{G_F^2}}{\displaystyle \frac{|m_{L}^{2}{}_{ij}{}^{}|^2}{m_S^8}}\mathrm{tan}^2\beta ,`$ (29)
where $`m_S^2`$ represents supersymmetric leptonic scalar masses.

From eq.(29) it is clear that BR($`l_il_j,\gamma `$) depends crucially on the quantity $`m_{L}^{2}{}_{ij}{}^{}`$, which at low energy is given by the integration of its RGE, as explained above. In the scenario examined in this paper (i.e. supersymmetric see-saw), and in the leading-log approximation, $`m_{L}^{2}{}_{ij}{}^{}(𝐘_\nu ^+𝐘_\nu )_{ij}`$, as given by eq.(3), so

$`\mathrm{BR}(l_il_j,\gamma ){\displaystyle \frac{\alpha ^3}{G_F^2m_S^8}}\left|{\displaystyle \frac{1}{8\pi ^2}}(3m_0^2+A_0^2)\mathrm{log}{\displaystyle \frac{M_X}{M}}\right|^2\left|(𝐘_\nu ^+𝐘_\nu )_{ij}\right|^2\mathrm{tan}^2\beta `$ (30)
The $`𝐘_\nu ^+𝐘_\nu `$ matrix, which is therefore the crucial quantity, has been evaluated in sect. 2 in the general and in particular cases. Next, we apply those results to the computation of BR($`l_il_j,\gamma `$),

## 4 Experimental data on neutrinos and predictions for BR($`l_il_j,\gamma `$)

We study here the expectable predictions for $`l_il_j,\gamma `$ from the existing experimental information on neutrinos. Let us mention that the present and forthcoming experimental bounds on $`l_il_j,\gamma `$ make $`\mu e,\gamma `$ the process with higher potential to constrain theories involving flavour violation. The present limit is BR($`\mu e,\gamma `$) $`<`$ $`1.2\times 10^{11}`$  and it is going to improve in the near future (year 2003) until $`10^{14}`$ . Hence, we will pay special attention to $`\mu e,\gamma `$. Nevertheless, as we will see, other processes like $`\tau \mu ,\gamma `$ and $`\tau e,\gamma `$ offer information complementary to that from $`\mu e,\gamma `$, which in some cases can be extremely important. For $`\tau \mu ,\gamma `$, in particular, the present (future) upper limit is $`1.1\times 10^6`$ ($`10^9`$) . Finally, there exist interesting limits on other branching ratios, like BR($`\mu eee`$) $`<`$ $`10^{12}`$ and BR($`\mu TieTi`$) $`<`$ $`6\times 10^{13}`$ . The latter is also to improve in projected experiments . These limits are certainly very small, but the theoretical predictions on the branching ratios are also much smaller due to the presence of extra electromagnetic vertices

From the results of the previous section, namely eq.(30), we know that $`\mathrm{BR}(l_il_j,\gamma )\left|(𝐘_\nu ^+𝐘_\nu )_{ij}\right|^2`$. On the other hand, from the general equation (14), it is clear that the matrix element $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$ will depend on several facts: the low-energy spectrum of neutrinos (contained in $`D_\kappa `$), the neutrino mixing angles (contained in $`U`$), the right-handed neutrino spectrum (contained in $`D_{}`$) and the choice of the $`R`$ matrix. Let us discuss them in order.

The experimental (solar and atmospheric) data  strongly suggest a hierarchy of neutrino mass-splittings, $`\mathrm{\Delta }\kappa _{sol}^2\mathrm{\Delta }\kappa _{atm}^2`$. Numerically, $`\mathrm{\Delta }\kappa _{atm}^2(1.46.1)\times 10^3\mathrm{e}V^2/v_2^4`$  (with central value around $`3\times 10^3\mathrm{e}V^2/v_2^4`$), while the value of $`\mathrm{\Delta }\kappa _{sol}^2`$ depends on the solution considered to explain the solar neutrino problem . The most favoured one from the recent analyses of data  is the large-angle MSW solution (LAMSW), which requires $`\mathrm{\Delta }\kappa _{sol}^23\times 10^5\mathrm{e}V^2/v_2^4`$. Other solutions are (in order of reliability from these analyses) the LOW MSW-solution, the small angle MSW-solution (SAMSW) and the vacuum oscillations solution (VO). They require, in $`\mathrm{e}V^2/v_2^4`$ units, $`\mathrm{\Delta }\kappa _{sol}^210^7`$, $`5\times 10^6`$ and $`8\times 10^{10}`$ respectively. The corresponding uncertainties are discussed e.g. in ref.. In any case, there are basically three types of neutrino spectra consistent with the hierarchy of mass-splittings : hierarchical ($`\kappa _1^2\kappa _2^2\kappa _3^2`$), “intermediate” ($`\kappa _1^2\kappa _2^2\kappa _3^2`$) and “degenerate” ($`\kappa _1^2\kappa _2^2\kappa _3^2`$). We follow the usual notation in which $`\kappa _3`$ corresponds to the most split mass eigenvalue, so that $`\mathrm{\Delta }\kappa _{atm}^2\mathrm{\Delta }\kappa _{32}^2`$, $`\mathrm{\Delta }\kappa _{sol}^2\mathrm{\Delta }\kappa _{21}^2`$.

The hierarchical scenario seems perhaps the most natural one, since it resembles the other fermion spectra. Indeed, it is the natural spectrum in GUT theories, particularly $`SO(10)`$ models. In any case, it is the most extensively analyzed scenario in the literature and we will pay special attention to it. In this scenario one should take $`\kappa _3^2\mathrm{\Delta }\kappa _{atm}^2`$, $`\kappa _2^2\mathrm{\Delta }\kappa _{sol}^2`$, $`\kappa _1^20`$. The other two scenarios (intermediate and degenerate) are also interesting and we will address them. In particular, the degenerate scenario is the only one compatible with a relevant cosmological role of neutrinos (provided $`\kappa 𝒪(1)\mathrm{e}V/v_2^2`$).

Concerning the mixing angles, $`\theta _{23}`$ and $`\theta _{13}`$ are constrained by the atmospheric and CHOOZ data to be near maximal and minimal, respectively. The $`\theta _{12}`$ angle depends on the solution considered to explain the solar neutrino problem. In the LAMSW case, $`\theta _{12}`$ is near maximal. In the other explanations of the solar neutrino problem $`\theta _{12}`$ can be either minimal (SAMSW), or near maximal (VO and LOW). The precise preferred values for the angles depend on the sets of data used in the fits. For example, from the global analysis of ref. these are $`\mathrm{tan}^2\theta _{13}0.005`$, $`\mathrm{tan}^2\theta _{23}1.4`$ and, for the LAMSW, $`\mathrm{tan}^2\theta _{12}0.36`$. The uncertainties depend also on the type of fit, as well as on the values used for the mass splittings. In summary, the two basic forms that $`U`$ can present are either a single-maximal or (more plausibly) a bimaximal mixing matrix. Schematically,

$$U\left(\begin{array}{ccc}1& 0& 0\\ 0& \frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}\\ 0& \frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}\end{array}\right)\mathrm{o}rU\left(\begin{array}{ccc}\frac{1}{\sqrt{2}}& \frac{1}{\sqrt{2}}& 0\\ \frac{1}{2}& \frac{1}{2}& \frac{1}{\sqrt{2}}\\ \frac{1}{2}& \frac{1}{2}& \frac{1}{\sqrt{2}}\end{array}\right),$$
(31)
With regard to the spectrum of right-handed neutrinos, there are two basic possibilities: either they are hierarchical, i.e. $`_1_2_3`$, or they are degenerate, i.e. $`_1_2_3`$. Both of them were considered in sect. 2, where it was also shown that the first possibility counts the case of an “inverse hierarchy” or any other ordering of the $``$ eigenvalues.

Finally, $`R`$ is the generic (complex) orthogonal matrix discussed in sect. 2. In principle, it contains three independent complex parameters, but we saw that this number can be usually reduced in a drastic way.

Let us analyze now the predictions for BR($`l_il_j,\gamma `$) in the different scenarios above exposed.

### 4.1 $`\nu _L`$’s and $`\nu _R`$’s completely hierarchical

#### 4.1.1 The generic case

We consider first the case where both the left-handed and the right-handed neutrinos are completely hierarchical, i.e.

$$D_\kappa \mathrm{d}iag(0,\kappa _2,\kappa _3),D_{}\mathrm{d}iag(0,0,_3).$$
(32)
We are neglecting for the moment contributions of $`𝒪(_2)`$. (Later on, we will discuss when those contributions would be relevant.) The general form of the $`𝐘_\nu `$ matrix in this case was already worked out in subsect. 2.1. If $`R_{32}0`$ or $`R_{33}0`$ (later on we will consider the case $`R_{32},R_{33}0`$), $`𝐘_\nu `$ is given by eq.(16), so only the last row of $`𝐘_\nu `$ has sizeable entries. Thus, the relevant quantity $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$ has the form

$`(𝐘_\nu ^+𝐘_\nu )_{ij}=(𝐘_\nu )_{3i}^{}(𝐘_\nu )_{3j}=_3\left[{\displaystyle \underset{l=2,3}{}}R_{3l}^{}\sqrt{\kappa _l}U_{il}\right]\left[{\displaystyle \underset{l^{}=2,3}{}}R_{3l}\sqrt{\kappa _l^{}}U_{jl^{}}^{}\right]`$ (33)
Focusing on $`\mu e,\gamma `$, and using the parametrization of $`R`$ of eq.(17), eq.(30) reads

$`\mathrm{BR}(\mu e,\gamma )`$ $``$ $`\left|(𝐘_\nu ^+𝐘_\nu )_{21}\right|^2`$ (34)
$`=`$ $`\left|_3|\widehat{c}_2|^2\left[\widehat{s}_1^{}\sqrt{\kappa _2}U_{22}+\widehat{c}_1^{}\sqrt{\kappa _3}U_{23}\right]\left[\widehat{s}_1\sqrt{\kappa _2}U_{12}^{}+\widehat{c}_1\sqrt{\kappa _3}U_{13}^{}\right]\right|^2`$
Let us recall here that, since $`\kappa _2`$ and $`\kappa _3`$ represent the solar and atmospheric splitting, respectively, then $`\kappa _2\kappa _3`$. Also, the matrix element $`U_{13}s_{13}1`$, as $`\theta _{13}`$ is the mixing angle constrained by CHOOZ. In the SAMSW scenario, $`U_{12}`$ and $`U_{23}`$ are also extremely small (see eq.(9)), so $`\mu e,\gamma `$ is suppressed. However, in the other scenarios (LAMSW, VO and LOW) $`U`$ is basically a bimaximal mixing matrix. Then $`U_{12},U_{23}1/\sqrt{2}`$, $`U_{22}1/2`$. Let us recall that this precisely includes the scenarios favoured from present data, i.e. LAMSW and LOW.

The prefactor $`_3|\widehat{c}_2|^2`$ in eq.(30) can be written as

$`_3|\widehat{c}_2|^2={\displaystyle \frac{|Y_0|^2}{|\widehat{s}_1|^2\kappa _2+|\widehat{c}_1|^2\kappa _3}}`$ (35)
where $`|Y_0|^2`$ is the largest eigenvalue of $`(𝐘_\nu ^+𝐘_\nu )`$. So, the right hand side of eq.(34) depends just on two independent parameters, $`|Y_0|`$ and $`\widehat{\theta }_1`$ (the remaining parameters have a precise physical meaning and are measurable at low energy). If we further impose a unification condition, like $`|Y_0(M_X)|=|Y_t(M_X)|`$, BR($`\mu e,\gamma `$) depends just on one parameter, $`\widehat{\theta }_1`$. This conclusion could also be obtained from eq.(18).

Now, as discussed in sect. 2, the results of ref.  correspond to take $`R=\mathrm{𝟏}`$, i.e. $`\widehat{c}_i=1`$. Then, BR($`\mu e,\gamma `$), as written in eq.(34) is suppressed by $`|U_{13}|^21`$. The authors of ref. considered however the contributions of $`𝒪(_2)`$ (which we have not considered yet). Then, working in the LAMSW scenario, they obtained a non-vanishing result for BR($`\mu e,\gamma `$), as a function of $`_2`$. What is apparent from eqs.(34, 35) is that that result corresponds to a very special point in the parameter space, namely $`\widehat{c}_1=1`$. In general, $`\widehat{c}_11,\widehat{s}_10`$, so the prediction for BR($`\mu e,\gamma `$) is much larger. For random values of $`\widehat{\theta }_1`$ we expect

$`(𝐘_\nu ^+𝐘_\nu )_{21}{\displaystyle \frac{|Y_0|^2}{|\widehat{s}_1|^2\kappa _2+|\widehat{c}_1|^2\kappa _3}}\widehat{c}_1^{}\widehat{s}_1\sqrt{\kappa _3\kappa _2}U_{23}U_{12}^{}`$ (36)
This formula allows to compare the foreseeable predictions for BR($`\mu e,\gamma `$) in different scenarios, taking into account the different values for $`\kappa _2`$, $`U_{ij}`$. Clearly, the largest values for BR($`\mu e,\gamma `$) are obtained for the LAMSW scenario (precisely, the most favoured one by the present analyses). For the LOW and SAMSW scenarios one expects predictions $`20`$ and $`1000`$ times smaller respectively. Indeed, for the LAMSW scenario the previous equation generically gives a branching ratio above the present experimental limits, at least for $`Y_0=𝒪(1)`$, as it occurs in the unified scenarios!

This is apparent from the plots of Fig. 3, where it is shown BR($`\mu e,\gamma `$) vs. $`\widehat{\theta }_1`$ in the LAMSW scenario for different typical sets of supersymmetric parameters, using a unification condition for the largest neutrino Yukawa coupling, i.e. $`|Y_0(M_X)|^2=|Y_t(M_X)|^2`$, with $`M_X=2\times 10^{16}`$ GeV. Incidentally, the value of the supersymmetric $`\mu `$-parameter is extracted (here and throughout the paper) from the electroweak breaking condition, choosing $`\mu >0`$. We have not shown the results for $`\mu <0`$ since they are always very similar. From Fig. 3, it is clear that most of the parameter space (i.e. the range of $`\widehat{\theta }_1`$) is already excluded by the present bounds on BR($`\mu e,\gamma `$), or it will become probed by the next generation of $`\mu e,\gamma `$ experiments (PSI), scheduled for 2003. For the sake of clarity, we took $`\widehat{\theta }_1`$ real in Fig. 3. The results for $`\widehat{\theta }_1`$ complex are very similar, as can be seen e.g. from eq.(36). Fig. 4 is completely analogous to Fig. 3, but for $`\mathrm{tan}\beta =10`$ (instead of $`\mathrm{tan}\beta =3`$). It illustrates the strong dependence of BR($`\mu e,\gamma `$) on $`\mathrm{tan}\beta `$ in the expected way ($`\mathrm{tan}^2\beta `$). An alternative representation of the allowed parameter space is given in Fig.5, where $`\widehat{\theta }_1`$ is given in logarithmic units (the set of supersymmetric parameters corresponds to that of the second plot of Fig.3). These can be more meaningful in particular model-building constructions. Of course, the basic results remain: $`\widehat{\theta }_1`$ should be $`<0.12`$ ($`<3\times 10^3`$) with the present (future) upper bounds on BR($`\mu e,\gamma `$). This corresponds to $`𝐘_{31}<0.04𝐘_{33}`$ ($`𝐘_{31}<9\times 10^4𝐘_{33}`$). This gives a measure of the constraints on the textures. (In a logarithmic plot it is not possible to say what percent of the parameter space is allowed since this includes the origin.)

It is worth-noticing from the curves of Figs. 3, 4, the existence of two values of $`\widehat{\theta }_1`$, namely $`\widehat{\theta }_1=0`$, $`\widehat{\theta }_12`$, in whose neighborhood the value of BR($`\mu e,\gamma `$) is drastically suppressed. We will comment shortly on this feature. For large enough values of the universal soft scalar mass, $`m_0`$, BR($`\mu e,\gamma `$) decreases, as expected: for large $`m_0`$, the branching ratio goes as $`m_0^4`$, as it is apparent from Fig. 6. Typically, for $`\mathrm{tan}\beta =3`$ BR($`\mu e,\gamma `$) does not fall below the present or forthcoming experimental upper bounds until $`m_0\stackrel{>}{_{}}1.5`$ TeV, i.e. beyond the reasonable range to avoid fine-tuning problems. For larger $`\mathrm{tan}\beta `$, this extremal value increases, since BR($`\mu e,\gamma `$) $``$ $`\mathrm{tan}^2\beta /m_0^4`$. Notice that there is a very narrow window at around 120 GeV, for $`A_0=0`$ GeV, where the branching ratio is very suppressed. This happens for every value of $`\widehat{\theta }_1`$ and is the result of a cancellation between the chargino and the neutralino amplitudes.

Finally, if one relaxes the unification condition on $`|Y_0(M_X)|`$, BR($`\mu e,\gamma `$) logically decreases, as shown in Fig. 7. Anyway, it is noticeable that even for $`Y_0(M_X)`$ ten times smaller than $`Y_t(M_X)`$, the conclusion remains that most of the parameter space will be probed in the forthcoming experiments.

#### 4.1.2 Special textures

Let us now turn to the two values of $`\widehat{\theta }_1`$ where BR($`\mu e,\gamma `$) is suppressed. They correspond to the values for which $`(𝐘_\nu ^+𝐘_\nu )_{21}`$, as given by eq.(33), is small. Notice that there are two possible ways to fulfill this requirement, namely $`(𝐘_\nu )_{31}0`$ or $`(𝐘_\nu )_{32}0`$, which translate into

$`\mathrm{tan}\widehat{\theta }_1`$ $``$ $`\sqrt{{\displaystyle \frac{\kappa _3}{\kappa _2}}}{\displaystyle \frac{U_{13}^{}}{U_{12}^{}}}=\sqrt{{\displaystyle \frac{\kappa _3}{\kappa _2}}}{\displaystyle \frac{V_{13}^{}}{V_{12}^{}}}0`$ (37)
$`\mathrm{tan}\widehat{\theta }_1`$ $``$ $`\sqrt{{\displaystyle \frac{\kappa _3}{\kappa _2}}}{\displaystyle \frac{U_{23}^{}}{U_{22}^{}}}=\sqrt{{\displaystyle \frac{\kappa _3}{\kappa _2}}}{\displaystyle \frac{V_{23}^{}}{V_{22}^{}}},`$ (38)
where we have used eqs.(18, 8). We do not know any reason why any of the previous conditions should be satisfied in a particular model, although of course it is not excluded. To this respect, the first condition, i.e. $`\widehat{\theta }_10`$, seems to be more “natural” to obtain than the second one: note that $`\widehat{\theta }_i=0`$, i.e. $`R=\mathrm{𝟏}`$, corresponds to the case where all the mixing can be attributed to the sector of charged leptons, as discussed in sect. 2, which is a physical possibility. It is interesting to explicitly show the textures of $`𝐘_\nu `$ that correspond to these ‘privileged’ values of $`\widehat{\theta }_1`$ from eqs.(37, 38). They read, respectively,

$`𝐘_\nu `$ $``$ $`\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 0& V_{31}& V_{21}\end{array}\right)\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 0& 1& 1\end{array}\right),`$ (39)
$`𝐘_\nu `$ $``$ $`\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ V_{31}& 0& V_{11}\end{array}\right)\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 1& 0& \sqrt{2}\end{array}\right).`$ (40)
Here we have used $`\mathrm{I}mV_{11}=0`$ and $`s_{13}1`$, so that $`\mathrm{I}mV_{21},\mathrm{I}mV_{31}1`$. Besides, the numerical values of the entries in (39, 40) correspond to a bimaximal mixing $`V`$ matrix, which includes the preferred LAMSW scenario. Let us stress that the previous textures are to be understood in the basis we have chosen to work, i.e. where both $`𝐘_𝐞`$ and $``$ are diagonal, and the latter has the eigenvalues ordered in the usual hierarchical way shown in eq.(32). It is interesting to note that texture (39) has been advocated in ref. on different grounds.

Some comments are in order here. As was discussed in subsect. 2.1, the first two rows of $`𝐘_\nu `$, though suppressed, are normally important to reconstruct $`\kappa `$ from $`𝐘_\nu `$, $``$. In particular, for texture (40) the precise values of the first two rows are crucial to reproduce the bimaximal structure of $`\kappa `$. On the other hand, this is not necessarily the case for texture (39). Actually, it is remarkable that texture (39) works fine, both to reproduce $`\kappa `$ and to give suppressed $`\mathrm{BR}(\mu e,\gamma )`$, even if it is given in a basis where $``$ is not diagonal.

In general, it is clear from eq.(33) that the $`l_il_j,\gamma `$ process is suppressed at the points in the parameter space where $`(𝐘_\nu )_{3i}0`$ or $`(𝐘_\nu )_{3j}0`$. Therefore, there are three special textures: the two given in the above eqs.(39, 40), which correspond to $`(𝐘_\nu )_{31}0`$ and $`(𝐘_\nu )_{32}0`$, and the one corresponding to $`(𝐘_\nu )_{33}0`$. The latter implies

$`\mathrm{tan}\widehat{\theta }_1\sqrt{{\displaystyle \frac{\kappa _3}{\kappa _2}}}{\displaystyle \frac{U_{33}^{}}{U_{32}^{}}}=\sqrt{{\displaystyle \frac{\kappa _3}{\kappa _2}}}{\displaystyle \frac{V_{33}^{}}{V_{32}^{}}},`$ (41)
and the associated texture is

$`𝐘_\nu \left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ V_{21}& V_{11}& 0\end{array}\right)\left(\begin{array}{ccc}0& 0& 0\\ 0& 0& 0\\ 1& \sqrt{2}& 0\end{array}\right),`$ (42)
where, again, the numerical values values correspond to the bimaximal scenario. If the texture of $`𝐘_\nu `$ is one of the special forms (39, 40, 42), there are two $`l_il_j,\gamma `$ processes suppressed and one unsuppressed. This means that for those privileged cases where $`\mu e,\gamma `$ is unusually suppressed, the $`\tau \mu ,\gamma `$ and $`\tau e,\gamma `$ processes are relevant and may give the most important restrictions on models. In particular, if $`𝐘_\nu `$ is as in eq.(39) or eq.(40) (i.e. $`(𝐘_\nu )_{31}0`$ or $`(𝐘_\nu )_{32}0`$), then $`\mathrm{BR}(\tau \mu ,\gamma )`$ or $`\mathrm{BR}(\tau e,\gamma )`$ respectively are unsuppressed.

We have illustrated this interesting aspect by plotting BR($`\tau \mu ,\gamma `$) vs. $`\widehat{\theta }_1`$ in Fig. 8 for one of the typical supersymmetric scenarios considered in Fig. 3. Notice that for the special value $`\widehat{\theta }_10`$, at which $`\mathrm{BR}(\mu e,\gamma )`$ is suppressed, BR($`\tau \mu ,\gamma `$) is not, lying above the forthcoming experimental upper bound.

#### 4.1.3 The contribution $`_2`$

So far, all the results of this subsection have been obtained neglecting the value of $`_2`$ (see eq.(32)). Let us consider here a value of $`_20`$, though smaller than $`_3`$, as it was done in ref..

Clearly, for generic values $`\widehat{\theta }_i`$, the contributions already computed \[see eqs.(3436)\] will be the dominant ones. However, if $`\widehat{\theta }_1`$ corresponds to one of the particular values of eqs.(37, 38), the contributions $`_2`$ can be relevant or even dominant. Normally, this will happen only for $`\widehat{\theta }_1`$ quite close to eqs.(37, 38). Suppose for example that $`R\mathrm{𝟏}`$, i.e. $`\widehat{s}_i`$ small and real. This corresponds to eq.(37) and it is probably the most interesting case (and the one considered in ref. ). Then $`(𝐘_\nu ^+𝐘_\nu )_{ij}=(𝐘_\nu )_{2i}^{}(𝐘_\nu )_{2j}+(𝐘_\nu )_{3i}^{}(𝐘_\nu )_{3j}`$ with $`𝐘_\nu `$ given by eq.(13). Keeping only terms up to $`𝒪(\widehat{s}_i)`$ we get

$`(𝐘_\nu ^+𝐘_\nu )_{21}`$ $``$ $`_2\kappa _2U_{12}^{}U_{22}+_3\kappa _3U_{13}^{}U_{23}`$ (43)
$`+`$ $`\widehat{s}_1(_3_2)\sqrt{\kappa _2\kappa _3}\left(U_{13}^{}U_{22}+U_{12}^{}U_{23}\right),`$
so that in a bimaximal case, even taking $`U_{13}=0`$, the part proportional to $`_3`$ will be indeed the dominant one except for $`\widehat{s}_1\stackrel{<}{_{}}(_2/_3)(\sqrt{\kappa _2/\kappa _3})`$, i.e. very close to zero.

For the specific case $`R=\mathrm{𝟏}`$, i.e. $`\widehat{s}_1=0`$, we get $`(𝐘_\nu ^+𝐘_\nu )_{21}_2\kappa _2U_{12}^{}U_{22}`$. Taking into account from eq.(35) $`|Y_0|^2=_3\kappa _3`$, the value of $`(𝐘_\nu ^+𝐘_\nu )_{21}`$ is exactly as the one we will obtain in another context in subsect. 4.2, see eq.(48) below, but multiplied by $`_2/_3`$. In consequence, the plots of $`\mathrm{BR}(\mu e,\gamma )`$ are identical to those of Fig. 9 below, but multiplied by $`(_2/_3)^2`$. For $`_2/_3\stackrel{>}{_{}}10^2`$ the branching ratio is still testable in the forthcoming generation of experiments.

#### 4.1.4 $`R_{32}R_{33}0`$

So far we have assumed $`R_{32}0`$ or $`R_{33}0`$ in eq.(13). Recall that $`R`$ is an unknown orthogonal matrix. The experimental data do not restrict its form, although $`R`$ may be much more definite in specific models. We consider now the special case when $`R_{32}R_{33}0`$. It should be kept in mind that this represents a narrow region of the whole parameter space of $`R`$.

However, as we are about to see, this possibility includes the particular case where all the leptonic flavour violation can be attributed to the sector of right-handed neutrinos. This happens when there exists a basis in which $`𝐘_𝐞`$ and $`𝐘_\nu `$ get simultaneously diagonal, while $``$ gets non-diagonal. This is physically interesting and it can occur in particular models. As explained in sect. 2, in the basis we have chosen to work ($`𝐘_𝐞`$ and $``$ diagonal), this implies $`𝐘_\nu =WD_Y`$, where $`D_Y`$ is a diagonal matrix and $`W`$ is a unitary matrix. Hence, $`𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu =D_{|Y^2|}`$. Then, from the general equation (13) or (14), we get

$`R^+D_{}R=D_{\sqrt{\kappa ^1}}U^+D_{|Y^2|}UD_{\sqrt{\kappa ^1}}.`$ (44)
If $`U`$ is of the bimaximal type, as it occurs in the preferred LAMSW scenario, then $`U_{i1}0`$ for any $`i`$ \[see eq.(31)\]. This means that, in first approximation, the only non-vanishing entry of the right hand side of eq.(44) is the (1,1), so the same holds for the left hand side. Since $`D_{}`$ is hierarchical, $`(R^+D_{}R)_{ij}R_{3i}^{}_3R_{3j}`$. Thus, $`R_{32}R_{33}0`$ q.e.d. A special case occurs when $`D_{|Y^2|}=|Y_0^2|\mathrm{𝟏}`$. Then eq.(44) becomes $`R^+D_{}R=|Y_0^2|D_{\kappa ^1}`$ for any form of $`U`$; and $`R`$ turns out to be simply the “permutation matrix” that interchanges the 1– and 3–planes: $`R=R_p=\{(0,0,\pm 1),(0,\pm 1,0),(\pm 1,0,0)\}`$. Of course, in all these cases, since $`𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu `$ is diagonal, there is no RG generation of non-diagonal soft terms, and the predictions for BR($`l_il_j,\gamma `$) are as in the SM, i.e. negligible.

Nevertheless, the possibility $`R_{32}R_{33}0`$ includes other cases different than the previous (almost trivial) one. As mentioned, this represents a small region of the parameter space, which in principle has no special physical significance; but for the sake of completeness we will analyze it. $`R`$ will admit the parametrization

$`R\left(\begin{array}{ccc}0& \pm \widehat{c}& \widehat{s}\\ 0& \widehat{s}& \widehat{c}\\ \pm 1& 0& 0\end{array}\right),`$ (45)
where $`\widehat{\theta }`$ is an arbitrary complex angle and the sign of the $`\pm 1`$ entry is not correlated to the signs appearing in the second column. Thus, from eq.(13), and keeping only the dominant terms,

$`𝐘_\nu \left(\begin{array}{ccc}0& \pm \sqrt{_1}\widehat{c}\sqrt{\kappa _2}& \sqrt{_1}\widehat{s}\sqrt{\kappa _3}\\ 0& \sqrt{_2}\widehat{s}\sqrt{\kappa _2}& \sqrt{_2}\widehat{c}\sqrt{\kappa _3}\\ \pm \sqrt{_3\kappa _1}& 0& 0\end{array}\right)U^+,`$ (46)
Therefore, the values of $`(𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu )_{ij}`$, $`ij`$, which are the relevant quantities for $`l_il_j,\gamma `$, depend on the value of $`\widehat{s}`$. In general, if a unification condition is imposed for the largest eigenvalue of $`𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu `$, then $`(𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu )_{21}`$ is quite large and BR($`\mu e,\gamma `$) is above the experimental limit.

Let us illustrate the last point with an example: suppose that $`\widehat{s}=\pm 1`$ in eq.(46). Then the eigenvalues of $`𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu `$ are $`\{_3\kappa _1,_2\kappa _2,_1\kappa _3\}`$. If they are degenerate, $`𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu `$ is diagonal and we are in the above-considered case: all the mixing can be attributed to the sector of right-handed neutrinos. If not, since $`U_{13}0`$, then $`(𝐘_{\nu }^{}{}_{}{}^{+}𝐘_\nu )_{21}U_{11}^{}U_{21}_3\kappa _1+U_{12}^{}U_{22}_2\kappa _2`$. This is at least as large as the matrix element computed in the generic case, eq.(36), thus leading to unacceptably large BR($`\mu e,\gamma `$). For generic values of $`\widehat{s}`$ the analysis is more involved, but at the end of the day the conclusion is similar.

### 4.2 $`\nu _L`$’s hierarchical and $`\nu _R`$’s degenerate

As was discussed in subsect. 2.2, if the $`R`$ matrix is real, this scenario is extremely predictive. In particular $`𝐘_\nu ^+𝐘_\nu `$ is given by eq.(21), which does not depend on the form of $`R`$. More precisely

$`(𝐘_\nu ^+𝐘_\nu )_{ij}=\kappa _lU_{il}U_{lj}^+\left[\kappa _2U_{i2}U_{j2}^{}+\kappa _3U_{i3}U_{j3}^{}\right],`$ (47)
where $`\kappa _2\kappa _{sol}`$ and $`\kappa _3\kappa _{atm}`$. The free parameter $``$ can be absorbed in the value of the largest eigenvalue of the $`𝐘_\nu ^+𝐘_\nu `$ matrix, $`|Y_0|^2\kappa _3`$. In particular, the matrix element $`(𝐘_\nu ^+𝐘_\nu )_{21}`$ becomes

$`(𝐘_\nu ^+𝐘_\nu )_{21}\kappa _2U_{22}U_{12}^{}=|Y_0|^2{\displaystyle \frac{\kappa _2}{\kappa _3}}U_{22}U_{12}^{},`$ (48)
where we have used $`U_{13}0`$. Comparing (48) with (36), we see that in this scenario we expect predictions for $`\mathrm{BR}(\mu e,\gamma )`$ and the other processes of the same order as for hierarchical neutrinos (more precisely, $`\kappa _2/\kappa _3`$ times smaller), but with no special textures where the branching ratio becomes suppressed. All this is illustrated in Fig. 9. for the LAMSW scenario using $`|Y_0(M_X)|=|Y_t(M_X)|`$.

From the plots, it is clear that in this scenario $`\mathrm{BR}(\mu e,\gamma )`$ is already above the present experimental limits except for a rather small region of $`m_0`$-values which should be probed by the next generation of experiments.

If $`R`$ is complex, the analysis is more involved since it contains more arbitrary parameters; but in general the conclusion is the same: $`\mathrm{BR}(\mu e,\gamma )`$ is at least of the same order as in the real case. To see this note that from the general equation (14)

$`(𝐘_\nu ^+𝐘_\nu )_{21}`$ $`=`$ $`\left[U_{2p}\sqrt{\kappa _p}(R^+R)_{pq}\sqrt{\kappa _q}U_{1q}^{}\right]`$ (49)
$``$ $`\left[U_{23}U_{12}^{}\sqrt{\kappa _3\kappa _2}R_{q3}^{}R_{q2}+U_{22}U_{12}^{}\kappa _2R_{q2}^{}R_{q2}\right].`$
Normally, this is at least of the same order as in the real case, eq.(48), as we have checked numerically. Now there exists, however, the possibility of a (fine-tuned) cancellation between the various contributions of (49). In particular such cancellation will occur if all the mixing can be attributed to the sector of right-handed neutrinos, which implies $`𝐘_\nu ^+𝐘_\nu `$ diagonal.

### 4.3 $`\nu _L`$’s quasi-degenerate

As it was discussed in subsect. 2.3, in this case it is logical to assume that $``$ has degenerate eigenvalues. If, furthermore, the $`R`$ matrix is real, we get from eq.(14)

$`𝐘_\nu ^+𝐘_\nu =UD_\kappa U^+,`$ (50)
Incidentally, for $`R`$ real the requirement that for quasi-degenerate neutrinos the two CP phases, $`\varphi `$ and $`\varphi ^{}`$, are opposite in order to be consistent with $`\nu `$-less double $`\beta `$ decay , implies that $`U`$ and thus $`𝐘_\nu `$ is complex, see eqs.(8, 9, 13).

Since $`\kappa _1\kappa _2\kappa _3\kappa `$ (typically $`\kappa =𝒪(1\mathrm{e}V/v_2^2)`$), the eigenvalues of $`𝐘_\nu ^+𝐘_\nu `$ are quasi-degenerate $`\kappa |Y_0|^2`$. Clearly, eq.(50) presents a GIM-like suppression, so the off-diagonal entries of $`𝐘_\nu ^+𝐘_\nu `$ are small and proportional to the differences of neutrinos masses. More precisely, for $`ij`$

$`(𝐘_\nu ^+𝐘_\nu )_{ij}=U_{ip}U_{jp}^{}\kappa _p=U_{ip}U_{jp}^{}(\kappa _p\kappa _3).`$ (51)
In particular, the matrix element $`(𝐘_\nu ^+𝐘_\nu )_{21}`$, which is the relevant one for $`\mu e,\gamma `$, becomes

$`(𝐘_\nu ^+𝐘_\nu )_{21}=\left[U_{21}U_{11}^{}(\kappa _1\kappa _3)+U_{22}U_{12}^{}(\kappa _2\kappa _3)\right].`$ (52)
Since $`U_{13}0`$, 666Since we are working here with quasi-degenerate neutrinos, one may wonder whether the low-energy condition $`U_{13}0`$ might change at the $``$ scale. However, at is was shown in ref., the $`U_{13}`$ element is stable under the RG running between $``$ and low-energy for any phenomenologically viable scenario. Therefore, the argument that follows is completely general. from unitarity follows that $`U_{21}U_{11}^{}U_{22}U_{12}^{}`$. So

$`(𝐘_\nu ^+𝐘_\nu )_{21}=\left[U_{21}U_{11}^{}(\kappa _1\kappa _2)\right]=|Y_0|^2U_{21}^{}U_{11}{\displaystyle \frac{\mathrm{\Delta }\kappa _{sol}^2}{\kappa ^2}}.`$ (53)
The factor $`\frac{\mathrm{\Delta }\kappa _{sol}^2}{\kappa ^2}`$ is at most $`10^5`$ (this occurs for the LAMSW solution to the solar neutrino problem), which implies a drastic suppression of $`(𝐘_\nu ^+𝐘_\nu )_{21}`$, and thus of BR($`\mu e,\gamma `$) in this kind of scenario. More quantitatively, eq.(53) is identical to eq.(48) (i.e. the scenario where $`\nu _L`$’s are hierarchical and $`\nu _R`$’s degenerate), multiplied by $`\kappa ^2\sqrt{\mathrm{\Delta }\kappa _{sol}^2\mathrm{\Delta }\kappa _{atm}^2}`$. This is a factor $`10^4`$ for the LAMSW, which means that all the plots representing BR($`\mu e,\gamma `$) in that scenario (Fig. 9) are identical here, but with the vertical axis re-scaled eight orders of magnitude smaller.

If the $`R`$ matrix is complex, then $`R^+R`$ is in general non-diagonal, so the $`𝐘_\nu ^+𝐘_\nu `$ matrix becomes

$`𝐘_\nu ^+𝐘_\nu =UD_\sqrt{\kappa }R^+RD_\sqrt{\kappa }U^+,`$ (54)
which may have sizeable off-diagonal entries. Hence, even for quasi-degenerate neutrinos, $`(𝐘_\nu ^+𝐘_\nu )_{21}`$, and thus BR($`\mu e,\gamma `$), could be very large if $`R`$ is complex.

Finally, we would like to comment on the scenario where the neutrinos are partially degenerate, i.e. when their masses present an inverse hierarchy: $`\kappa _3\kappa _1\kappa _2\kappa `$. \[It is convenient to maintain $`\kappa _3`$ as the most split mass eigenvalue, so that the interpretation of $`\theta _{23}`$, $`\theta _{12}`$ as the atmospheric and solar angles respectively remains valid.\] Notice that in this case $`\kappa ^2\mathrm{\Delta }\kappa _{atm}^2`$.

Now, the discussion of the size of $`(𝐘_\nu ^+𝐘_\nu )_{21}`$, and thus of BR($`\mu e,\gamma `$), is similar to that in the complete degenerate case, though the conclusions are different. In particular, suppose that $``$ has degenerate eigenvalues, $`D_{}=\mathrm{d}iag(,,)`$, and that the $`R`$ matrix is real. Then eq.(50) holds, so $`𝐘_\nu ^+𝐘_\nu `$ has two large quasi-degenerate eigenvalues $`\kappa |Y_0|^2`$. Likewise, all the discussion between the eqs. (51)–(53) and the equations themselves remain valid. Therefore, the matrix element $`(𝐘_\nu ^+𝐘_\nu )_{21}`$ becomes

$`(𝐘_\nu ^+𝐘_\nu )_{21}=|Y_0|^2U_{21}U_{11}^{}{\displaystyle \frac{\mathrm{\Delta }\kappa _{sol}^2}{\kappa ^2}}=|Y_0|^2U_{21}U_{11}^{}{\displaystyle \frac{\mathrm{\Delta }\kappa _{sol}^2}{\mathrm{\Delta }\kappa _{atm}^2}}.`$ (55)
Again, this equation is identical to eq.(48), multiplied now by $`\sqrt{\frac{\mathrm{\Delta }\kappa _{sol}^2}{\mathrm{\Delta }\kappa _{atm}^2}}`$. This suppression factor is not as strong as before. Namely, it is $`10^1`$ for the LAMSW, which means that all the plots representing BR($`\mu e,\gamma `$) in this scenario are as in the one where $`\nu _L`$’s are hierarchical and $`\nu _R`$’s degenerate, re-scaled by a factor $`10^2`$. Hence, from Fig. 9 it is clear that BR($`\mu e,\gamma `$) for the partially degenerate case should be testable within the next generation of experiments.

The matrix element $`(𝐘_\nu ^+𝐘_\nu )_{31}`$ gets a similar suppression. However $`(𝐘_\nu ^+𝐘_\nu )_{32}`$ is not suppressed at all:

$`(𝐘_\nu ^+𝐘_\nu )_{32}=\left[U_{31}U_{21}^{}(\kappa _1\kappa _3)+U_{32}U_{22}^{}(\kappa _2\kappa _3)\right]\kappa \left[U_{31}U_{21}^{}+U_{32}U_{22}^{}\right],`$ (56)
where $`\kappa =|Y_0|^2`$. This is a similar size to that from the hierarchical scenarios. So, for the partially degenerate case, BR($`\tau \mu ,\gamma `$) should also be testable in the next generation of experiments.

Again, if the $`R`$ matrix is complex $`𝐘_\nu ^+𝐘_\nu `$ is given by eq.(54) and we expect much larger values for $`(𝐘_\nu ^+𝐘_\nu )_{21}`$, $`(𝐘_\nu ^+𝐘_\nu )_{31}`$.

## 5 Summary of predictions for BR($`l_il_j,\gamma `$)

If the origin of the neutrino masses is a see-saw mechanism, implemented in a supersymmetric theory, the soft breaking terms get off-diagonal contributions through the RG running, even if one starts with a universality condition for the soft terms. E.g. for the slepton doublet mass-matrix, at the leading-log approximation,

$`\left(m_L^2\right)_{ij}`$ $``$ $`{\displaystyle \frac{1}{8\pi ^2}}(3m_0^2+A_0^2)(𝐘_\nu ^+𝐘_\nu )_{ij}\mathrm{log}{\displaystyle \frac{M_X}{M}},`$ (57)
where $`𝐘_\nu `$ is the matrix of neutrino Yukawa couplings, $`M`$ is the typical (Majorana) mass of the right-handed neutrinos and $`M_X`$ is the initial scale at which universality is imposed. Therefore, the predictions on BR($`l_il_j,\gamma `$) are directly linked to the size of $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$. Moreover, the larger (smaller) the initial scale at which universality is imposed the larger the branching ratios. Also, due to the structure of the diagrams, BR($`l_il_j,\gamma `$) has a strong dependence on $`\mathrm{tan}\beta `$ ($`\mathrm{tan}^2\beta `$). The most general form of $`𝐘_\nu `$ and thus of $`𝐘_\nu ^+𝐘_\nu `$, in terms of measurable low-energy parameters has been determined in sect. 2. In particular, working in the flavour basis in which the gauge interactions and the charged-lepton Yukawa matrix, $`𝐘_𝐞`$, are flavour-diagonal,

$`𝐘_\nu ^+𝐘_\nu =UD_\sqrt{\kappa }R^+D_{}RD_\sqrt{\kappa }U^+`$ (58)
Here $`D_\sqrt{\kappa }`$ is the diagonal matrix of positive neutrino mass eigenvalues ($`\kappa =_\nu /H_2^0^2`$), while $`U`$ contains the mixing angles and CP phases \[see eqs.(6, 7)\]. In addition, $`𝐘_\nu ^+𝐘_\nu `$ depends on unknown parameters: the three positive mass eigenvalues of the righthanded neutrinos (contained in $`D_{}`$) and the three complex parameters defining the complex orthogonal matrix $`R`$. In practical cases, however, the number of relevant free parameters becomes drastically reduced.

Next, we enumerate the possible low-energy scenarios and the corresponding predictions for $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$ and BR($`l_il_j,\gamma `$), as has been worked out in detail in sects. 2, 3, 4.

### $`\nu _L`$’s and $`\nu _R`$’s completely hierarchical

* If $`R`$ is a generic matrix, $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$ is given by eq.(33). Using the parametrization of $`R`$ given in eq.(17), which is sufficiently general for this case, one obtains in particular

$`(𝐘_\nu ^+𝐘_\nu )_{21}{\displaystyle \frac{|Y_0|^2}{|\widehat{s}_1|^2\kappa _2+|\widehat{c}_1|^2\kappa _3}}\widehat{c}_1^{}\widehat{s}_1\sqrt{\kappa _3\kappa _2}U_{23}U_{12}^{}`$ (59)
Here $`|Y_0|^2`$ is the largest eigenvalue of $`𝐘_\nu ^+𝐘_\nu `$ \[see eq.(35)\] and $`\widehat{\theta }_1`$ is an arbitrary complex angle.

For the LAMSW scenario the previous equation generically gives BR($`\mu e,\gamma `$) above the present experimental limits, at least for $`Y_0=𝒪(1)`$, as it occurs in the unified scenarios. This is illustrated in Figs. 3, 4, 5, 6 using $`|Y_0(M_X)|=|Y_t(M_X)|`$, with $`M_X=2\times 10^{16}`$ GeV. But even for $`Y_0=𝒪(10^1)`$, most of the parameter space will be probed in the forthcoming generation of experiments experiments (see Fig. 7).
* The only exceptions to the previous result are

+ If $`\widehat{\theta }_1`$ takes one of the two values given in eqs.(37, 38), which correspond to the textures of eqs.(39, 40). Then, BR($`\mu e,\gamma `$) is drastically reduced, but other processes, as $`\mathrm{BR}(\tau \mu ,\gamma )`$, are not, normally lying above the forthcoming experimental upper bound, see Fig. 8.
    + If $`R`$ is such that $`𝐘_\nu ^+𝐘_\nu `$ is diagonal. This means that in a certain basis $`𝐘_𝐞`$ and $`𝐘_\nu `$ are diagonal whereas the right-handed Majorana mass matrix, $``$, is not. This requires a very special form of $`R`$, which in particular has $`R_{32},R_{33}0`$. Of course, in this special case there is no RG generation of non-diagonal soft terms, and the predictions for BR($`l_il_j,\gamma `$) are as in the SM, i.e. negligible.

### $`\nu _L`$’s hierarchical and $`\nu _R`$’s degenerate

* If $`R`$ is real, this scenario is very predictive. $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$ is given by eq.(47). In particular, taking into account $`U_{13}0`$,

$`(𝐘_\nu ^+𝐘_\nu )_{21}\kappa _2U_{22}U_{12}^{}=|Y_0|^2{\displaystyle \frac{\kappa _2}{\kappa _3}}U_{22}U_{12}^{},`$ (60)
where $`|Y_0|^2\kappa _3`$ is the largest eigenvalue of $`𝐘_\nu ^+𝐘_\nu `$.

The corresponding $`\mathrm{BR}(\mu e,\gamma )`$ is shown in Fig. 9. for the LAMSW scenario, and using $`|Y_0(M_X)|=|Y_t(M_X)|`$. 777For other values of $`|Y_0|`$ note that $`\mathrm{BR}(\mu e,\gamma )`$) $``$ $`|Y_0|^4`$. Actually, this dependence is milder, since the larger $`|Y_0(M_X)|`$, the faster it decreases with the scale. The branching ratio turns out to be already above the present experimental limits except for a rather small region of $`m_0`$ values which should be probed by the next generation of experiments. Here are not special textures where the branching ratio becomes suppressed.
* If $`R`$ is complex, the analysis is more involved since it contains more arbitrary parameters; $`(𝐘_\nu ^+𝐘_\nu )_{21}`$ is given by eq.(49). But in general the conclusion is the same: $`\mathrm{BR}(\mu e,\gamma )`$ is at least of the same order as in the real case. Now there exists, however, the possibility of a (fine-tuned) cancellation between the various contributions of (49). In particular such cancellation will occur if all the mixing can be attributed to the sector of right-handed neutrinos, which implies $`𝐘_\nu ^+𝐘_\nu `$ diagonal.

### $`\nu _L`$’s quasi-degenerate

* If $`R`$ is real, $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$ is given by eq.(50). In particular $`(𝐘_\nu ^+𝐘_\nu )_{21}`$ is given by eq.(53)

$`(𝐘_\nu ^+𝐘_\nu )_{21}=\left[U_{21}U_{11}^{}(\kappa _1\kappa _2)\right]=|Y_0|^2U_{21}U_{11}^{}{\displaystyle \frac{\mathrm{\Delta }\kappa _{sol}^2}{\kappa ^2}}.`$ (61)
where, again, $`|Y_0|^2\kappa `$ is the largest eigenvalue of $`𝐘_\nu ^+𝐘_\nu `$. This equation is identical to eq.(60), multiplied by $`\kappa ^2\sqrt{\mathrm{\Delta }\kappa _{sol}^2\mathrm{\Delta }\kappa _{atm}^2}`$. This is a factor $`10^4`$ for the LAMSW. Therefore all the plots representing BR($`\mu e,\gamma `$) in the previous scenario (Fig. 9) are valid here, but with the vertical axis re-scaled eight orders of magnitude smaller. Consequently, BR($`\mu e,\gamma `$) is naturally suppressed below the present (and even forthcoming) limits.
* If $`R`$ is complex, $`𝐘_\nu ^+𝐘_\nu `$ is given by eq.(54), which may have sizeable off-diagonal entries. Hence, BR($`\mu e,\gamma `$), could be very large in this case.
* If the (quasi-) degeneracy is only partial: $`\kappa _3\kappa _1\kappa _2\kappa \sqrt{\mathrm{\Delta }\kappa _{atm}^2}`$, $`(𝐘_\nu ^+𝐘_\nu )_{21}`$ is given (for $`R`$ real) by eq.(55). Again, this equation is identical to eq.(60), multiplied now by $`\sqrt{\frac{\mathrm{\Delta }\kappa _{sol}^2}{\mathrm{\Delta }\kappa _{atm}^2}}`$. This represents a suppression factor $`10^1`$ for the LAMSW, which means that all the plots representing BR($`\mu e,\gamma `$) are as those of Fig. 9, re-scaled by a factor $`10^2`$. As a consequence, BR($`\mu e,\gamma `$) for this partially degenerate scenario should be testable within the next generation of experiments. The conclusion is similar for BR($`\tau \mu ,\gamma `$).

For generic complex $`R`$, the value of BR($`\mu e,\gamma `$) does not get any suppression and falls naturally above the present experimental limits.

## 6 Concluding remarks

If the origin of the neutrino masses is a supersymmetric see-saw, which is probably the most attractive scenario to explain their smallness, then the leptonic soft breaking terms acquire off-diagonal contributions through the RG running, which drive non-vanishing BR($`l_il_j,\gamma `$). These contributions are proportional to $`(𝐘_\nu ^+𝐘_\nu )_{ij}`$, where $`𝐘_\nu `$ is the neutrino Yukawa matrix,

Therefore, in order to make predictions for these branching ratios, one has first to determine the most general form of $`𝐘_\nu `$ and $`𝐘_\nu ^+𝐘_\nu `$, compatible with all the phenomenological requirements. This has been the first task of this paper, and the result is summarized in eqs.(13, 14).

Then, we have shown that the predictions for BR($`\mu e,\gamma `$) are normally above the present experimental limits if the three following conditions occur

1. The solution to the solar neutrino problem is the LAMSW, as favoured by the most recent analyses.
2. $`Y_0(M_X)=𝒪(1)`$, where $`|Y_0|^2`$ is the largest eigenvalue of $`(𝐘_\nu ^+𝐘_\nu )`$. This occurs e.g. in most grand-unified scenarios.
3. The soft-breaking terms are generated at a high-energy scale, e.g. $`M_X`$, above the Majorana mass of the right-handed neutrinos, $`M`$.

These conditions are very plausible. In our opinion, the most natural scenarios fulfill them, but certainly there exists other possibilities. E.g. concerning condition 1, the solution to the solar neutrino problem could be LOW. In this case, BR($`\mu e,\gamma `$) is ($`20`$ times) smaller, but still testable in the forthcoming experiments. Condition 2 is satisfied in $`SO(10)`$ models, but it could be unfulfilled in $`SU(5)`$ models or in string scenarios with no GUT group. Still, we have noted in this work that even for $`Y_0=𝒪(10^1)`$, most of the parameter space will be probed in the forthcoming generation of experiments (year 2003). Finally, condition 3 seems a natural one, as the see-saw mechanism requires the existence of high-energy scales, and it is hardly compatible e.g. with models where the fundamental scale is $`𝒪`$ (TeV). Actually, if the initial scale at which universality is imposed is larger than the ordinary GUT $`M_X`$ scale (as it is common in gravity mediated supersymmetry breaking), the predictions for the branching ratios increase. However, it may happen that supersymmetry is broken at a scale below $`M`$. This is the case of gauge-mediated scenarios, where there would be no generation of off-diagonal leptonic soft terms through the RG running. In this sense, the gauge-mediated models are free from $`\mu e,\gamma `$ constraints and represent a way-out to the previous conclusion.

Even under the previous 1–3 conditions, there are physical scenarios compatible with the present BR($`\mu e,\gamma `$) experimental limits. Namely

* Whenever all the leptonic flavour violation can be attributed to the sector of right-handed neutrinos. This happens when there exists a basis in which the gauge couplings, $`𝐘_𝐞`$ and $`𝐘_\nu `$ get simultaneously flavour-diagonal, while the right-handed Majorana matrix, $``$, gets non-diagonal. In this case there is no RG generation of non-diagonal soft terms, and the predictions for BR($`l_il_i,\gamma `$) are as in the SM, i.e. negligible.

This possibility certainly exists and might be considered an attractive one. However, it should be noticed that it does not occur in the quark sector since the $`𝐘_𝐮`$ and $`𝐘_𝐝`$ Yukawa matrices cannot be simultaneously diagonalized in a basis where the gauge couplings remain diagonal. In our opinion, it would be a big chance that this happens in the leptonic sector.
* In the scenario where both, the left-handed and the right-handed neutrinos have hierarchical masses, there are two special textures, shown in eqs.(39, 40), where BR($`\mu e,\gamma `$) becomes suppressed. The first one is probably the most appealing, since it takes place for the simple case $`R=\mathrm{𝟏}`$ in eq.(13). This is equivalent to assume that there exists a basis of $`L_i`$ and $`\nu _{R}^{}{}_{i}{}^{}`$ in which $`𝐘_\nu `$ and $``$ are simultaneously diagonal, though not $`𝐘_𝐞`$. In other words, all the leptonic flavour violation can be attributed to the sector of charged leptons.

Then, we have shown that for this privileged texture BR($`\tau \mu ,\gamma `$) is not suppressed, lying in general above the forthcoming experimental upper bound.
* If the left-handed neutrinos are quasi-degenerate, then BR($`l_il_j,\gamma `$) are naturally suppressed. Namely, if the $`R`$ matrix in eq.(13) is real, there is a GIM-like mechanism that drives BR($`l_il_j,\gamma `$) below the experimental limits, even the forthcoming ones. (This is not the case, however, if $`R`$ is complex.)

If the degeneracy is partial (case of inverse hierarchy of neutrino masses), BR($`\mu e,\gamma `$) is somewhat suppressed, but still testable in the forthcoming experiments. The latter is also true for BR($`\tau \mu ,\gamma `$)

In our opinion, if the above 1–3 conditions are fulfilled, the scenario of quasi-degenerate neutrinos and the one with gauge mediated supersymmetry breaking represent the most plausible explanations to the absence of $`\mu e,\gamma `$ observations, specially if the absence persists after the next generation of experiments.

As a final conclusion, the discovery of neutrino oscillations makes much more plausible the possibility of observing lepton-flavour-violation processes, specially $`\mu e,\gamma `$, if the theory is supersymmetric and the neutrino masses are generated by a see-saw mechanism. Large regions of the parameter space are already excluded on these grounds, and there exists great chances to observe $`\mu e,\gamma `$ in the near future (PSI, 2003). This means that, hopefully, we will have signals of supersymmetry before LHC.

## Acknowledgements

The work of J.A.C. was supported in part by the CICYT (contract AEN98-0816). A.I. would like to thank JoAnne Hewett and Gudrun Hiller, and specially Sacha Davidson and Graham Ross for very interesting discussions and invaluable suggestions. A.I. would also like to thank the SLAC theory group for hospitality during the first stages of this work.

## Appendix

In the appendix, we summarize the relevant formulas concerning RGEs and the coefficients of the amplitude appearing in eq.(27).

Between $`M_X`$ and $`M`$ the evolution of the Yukawa couplings is given by

$`{\displaystyle \frac{d𝐘_\nu }{dt}}={\displaystyle \frac{1}{16\pi ^2}}𝐘_\nu \left[\left(3g_2^2+{\displaystyle \frac{3}{5}}g_1^2\mathrm{T}_2\right)𝐈_\mathrm{𝟑}\left(3𝐘_\nu ^+𝐘_\nu +𝐘_𝐞^+𝐘_𝐞\right)\right],`$ (62)
$`{\displaystyle \frac{d𝐘_𝐞}{dt}}={\displaystyle \frac{1}{16\pi ^2}}𝐘_𝐞\left[\left(3g_2^2+{\displaystyle \frac{9}{5}}g_1^2\mathrm{T}_1\right)𝐈_\mathrm{𝟑}\left(𝐘_\nu ^+𝐘_\nu +3𝐘_𝐞^+𝐘_𝐞\right)\right],`$ (63)
$`{\displaystyle \frac{d𝐘_𝐔}{dt}}={\displaystyle \frac{1}{16\pi ^2}}𝐘_𝐔\left[\left({\displaystyle \frac{13}{15}}g_1^2+3g_2^2+{\displaystyle \frac{16}{3}}g_3^2T_2\right)𝐈_33𝐘_𝐔^+𝐘_𝐔𝐘_𝐃^+𝐘_𝐃\right],`$ (64)
where

$`\mathrm{T}_1=\mathrm{T}r(3𝐘_𝐃^+𝐘_𝐃+𝐘_𝐞^+𝐘_𝐞),\mathrm{T}_2=\mathrm{T}r(3𝐘_𝐔^+𝐘_𝐔+𝐘_\nu ^+𝐘_\nu ),`$ (65)
and

$`{\displaystyle \frac{d}{dt}}={\displaystyle \frac{1}{8\pi ^2}}\left[(𝐘_\nu 𝐘_{\nu }^{}{}_{}{}^{+})^T+𝐘_\nu 𝐘_{\nu }^{}{}_{}{}^{+}\right];`$ (66)
The renormalization group equations for the soft masses and trilinear terms are:

$`{\displaystyle \frac{d𝐦_𝐋^\mathrm{𝟐}}{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}[(𝐦_𝐋^\mathrm{𝟐}𝐘_𝐞^+𝐘_𝐞+𝐘_𝐞^+𝐘_𝐞𝐦_𝐋^\mathrm{𝟐})+(𝐦_𝐋^\mathrm{𝟐}𝐘_\nu ^+𝐘_\nu +𝐘_\nu ^+𝐘_\nu 𝐦_𝐋^\mathrm{𝟐})`$ (67)
$`+2\left(𝐘_𝐞^+𝐦_𝐞^\mathrm{𝟐}𝐘_𝐞+m_{H_1}^2𝐘_𝐞^+𝐘_𝐞+𝐀_𝐞^+𝐀_𝐞\right)+2\left(𝐘_\nu ^+𝐦_\nu ^\mathrm{𝟐}𝐘_\nu +m_{H_2}^2𝐘_\nu ^+𝐘_\nu +𝐀_\nu ^+𝐀_\nu \right)`$
$`\left({\displaystyle \frac{6}{5}}g_1^2\right|M_1|^2+6g_2^2\left|M_2|^2\right)𝐈_\mathrm{𝟑}{\displaystyle \frac{3}{5}}g_1^2S𝐈_\mathrm{𝟑}],`$
$`{\displaystyle \frac{d𝐦_𝐞^\mathrm{𝟐}}{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}[2(𝐦_𝐞^\mathrm{𝟐}𝐘_𝐞𝐘_𝐞^++𝐘_𝐞𝐘_𝐞^+𝐦_𝐞^\mathrm{𝟐})+4(𝐘_𝐞𝐦_𝐋^\mathrm{𝟐}𝐘_𝐞^++m_{H_1}^2𝐘_𝐞𝐘_𝐞^++𝐀_𝐞𝐀_𝐞^+)`$ (68)
$`{\displaystyle \frac{24}{5}}g_1^2|M_1|^2𝐈_\mathrm{𝟑}+{\displaystyle \frac{6}{5}}g_1^2S𝐈_\mathrm{𝟑}],`$
$`{\displaystyle \frac{d𝐦_\nu ^\mathrm{𝟐}}{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}\left[2\left(𝐦_\nu ^\mathrm{𝟐}𝐘_\nu 𝐘_\nu ^++𝐘_\nu 𝐘_\nu ^+𝐦_\nu ^\mathrm{𝟐}\right)+4\left(𝐘_\nu 𝐦_𝐋^\mathrm{𝟐}𝐘_\nu ^++m_{H_2}^2𝐘_\nu 𝐘_\nu ^++𝐀_\nu 𝐀_\nu ^+\right)\right],`$ (69)
$`{\displaystyle \frac{dm_{H_2}^2}{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}[6\mathrm{T}r(𝐘_𝐔^+(𝐦_𝐐^\mathrm{𝟐}+𝐦_𝐔^\mathrm{𝟐}+m_{H_2}^2𝐈_\mathrm{𝟑})𝐘_𝐔+𝐀_𝐔^+𝐀_𝐔)`$ (70)
$`+2\mathrm{T}r\left(𝐘_\nu ^+(𝐦_𝐋^\mathrm{𝟐}+𝐦_\nu ^\mathrm{𝟐}+m_{H_2}^2𝐈_\mathrm{𝟑})𝐘_\nu +𝐀_\nu ^+𝐀_\nu \right)`$
$`\left({\displaystyle \frac{6}{5}}g_1^2\right|M_1|^2+6g_2^2\left|M_2|^2\right)+{\displaystyle \frac{3}{5}}g_1^2S],`$
$`{\displaystyle \frac{d𝐀_𝐞}{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}[\{{\displaystyle \frac{9}{5}}g_1^23g_2^2+3\mathrm{T}r(𝐘_𝐃^+𝐘_𝐃)+\mathrm{T}r(𝐘_𝐞^+𝐘_𝐞)\}𝐀_𝐞`$ (71)
$`+2\left\{{\displaystyle \frac{9}{5}}g_1^2M_13g_2^2M_2+3\mathrm{T}r(𝐘_𝐃^+𝐀_𝐃)+\mathrm{T}r(𝐘_𝐞^+𝐀_𝐞)\right\}𝐘_𝐞`$
$`+4(𝐘_𝐞𝐘_𝐞^+𝐀_𝐞)+5(𝐀_𝐞𝐘_𝐞^+𝐘_𝐞)+2(𝐘_𝐞𝐘_\nu ^+𝐀_\nu )+(𝐀_𝐞𝐘_\nu ^+𝐘_\nu )],`$
$`{\displaystyle \frac{d𝐀_\nu }{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}[\{{\displaystyle \frac{3}{5}}g_1^23g_2^2+3\mathrm{T}r(𝐘_𝐔^+𝐘_𝐔)+\mathrm{T}r(𝐘_\nu ^+𝐘_\nu )\}𝐀_\nu `$ (72)
$`+2\left\{{\displaystyle \frac{3}{5}}g_1^2M_13g_2^2M_2+3\mathrm{T}r(𝐘_𝐔^+𝐀_𝐔)+\mathrm{T}r(𝐘_\nu ^+𝐀_\nu )\right\}𝐘_\nu `$
$`+4(𝐘_\nu 𝐘_\nu ^+𝐀_\nu )+5(𝐀_\nu 𝐘_\nu ^+𝐘_\nu )+2(𝐘_\nu 𝐘_𝐞^+𝐀_𝐞)+(𝐀_\nu 𝐘_𝐞^+𝐘_𝐞)],`$
$`{\displaystyle \frac{d𝐀_𝐔}{dt}}`$ $`=`$ $`{\displaystyle \frac{1}{16\pi ^2}}[\{{\displaystyle \frac{13}{15}}g_1^23g_2^2{\displaystyle \frac{16}{3}}g_3^2+3\mathrm{T}r(𝐘_𝐔^+𝐘_𝐔)+\mathrm{T}r(𝐘_\nu ^+𝐘_\nu )\}𝐀_𝐔`$ (73)
$`+2\left\{{\displaystyle \frac{13}{15}}g_1^2M_13g_2^2M_2{\displaystyle \frac{16}{3}}g_3^2M_3+3\mathrm{T}r(𝐘_𝐔^+𝐀_𝐔)+\mathrm{T}r(𝐘_\nu ^+𝐀_\nu )\right\}𝐘_𝐔`$
$`+4(𝐘_𝐔𝐘_𝐔^+𝐀_𝐔)+5(𝐀_𝐔𝐘_𝐔^+𝐘_𝐔)+2(𝐘_𝐔𝐘_𝐃^+𝐀_𝐃)+(𝐀_𝐔𝐘_𝐃^+𝐘_𝐃)],`$
where

$`S=\mathrm{T}r(𝐦_𝐐^\mathrm{𝟐}+𝐦_𝐝^\mathrm{𝟐}2𝐦_𝐮^\mathrm{𝟐}𝐦_𝐋^\mathrm{𝟐}+𝐦_𝐞^\mathrm{𝟐})m_{H_1}^2+m_{H_2}^2.`$ (74)
Here $`g_2`$ and $`g_1`$ are the $`SU(2)_L`$ and $`U(1)_Y`$ gauge coupling constants, and $`𝐘_{𝐔,𝐃,𝐞,\nu }`$ are the Yukawa matrices for up quarks, down quarks, charged leptons and neutrinos respectively. The RGEs for $`𝐘_𝐃`$ and the rest of the soft masses or trilinear terms are the same as in the MSSM. Below $`M`$, the RGEs are the same, except that the right neutrinos decouple, and hence $`𝐘_\nu `$ disappears from the equations.

On the other hand, following ref. the amplitude for the process $`l_il_j\gamma `$ can be written as

$`T=ϵ^\alpha \overline{l}_jm_{l_i}i\sigma _{\alpha \beta }q^\beta (A_LP_L+A_RP_R)l_i,`$ (75)
$`A_{L,R}`$ $`=`$ $`A_{L,R}^{(c)}+A_{L,R}^{(n)}`$
$`A_L^{(c)}`$ $`=`$ $`{\displaystyle \frac{e}{16\pi ^2}}{\displaystyle \frac{1}{m_{\stackrel{~}{\nu }_X}^2}}\left[C_{jAX}^LC_{iAX}^LI_1(M_{\stackrel{~}{\chi }_A^{}}^2/m_{\stackrel{~}{\nu }_X}^2)+C_{jAX}^LC_{iAX}^R{\displaystyle \frac{M_{\stackrel{~}{\chi }_A^{}}}{m_{l_i}}}I_1(M_{\stackrel{~}{\chi }_A^{}}^2/m_{\stackrel{~}{\nu }_X}^2)\right],`$ (76)
$`A_L^{(n)}`$ $`=`$ $`{\displaystyle \frac{e}{16\pi ^2}}{\displaystyle \frac{1}{m_{\stackrel{~}{l}_X}^2}}\left[N_{jAX}^LN_{iAX}^LJ_1(M_{\stackrel{~}{\chi }_A^0}^2/m_{\stackrel{~}{l}_X}^2)+N_{jAX}^LN_{iAX}^R{\displaystyle \frac{M_{\stackrel{~}{\chi }_A^0}}{m_{l_i}}}J_2(M_{\stackrel{~}{\chi }_A^0}^2/m_{\stackrel{~}{l}_X}^2)\right],`$ (77)
$`A_R^{(c,n)}`$ $`=`$ $`A_L^{(c,n)}|_{LR},`$ (78)
where $`M_{\stackrel{~}{\chi }_A^{}}(M_{\stackrel{~}{\chi }_A^0})`$ is the chargino (neutralino) mass and $`m_{\stackrel{~}{\nu }_X}^2(m_{\stackrel{~}{l}_X}^2)`$ is the sneutrino (charged slepton) mass squared.

The functions $`I_{1,2}`$ and $`J_{1,2}`$ are defined as follows:

$`I_1(r)`$ $`=`$ $`{\displaystyle \frac{1}{12(1r)^4}}(2+3r6r^2+r^3+6r\mathrm{ln}r),`$ (79)
$`I_2(r)`$ $`=`$ $`{\displaystyle \frac{1}{2(1r)^3}}(3+4rr^22\mathrm{ln}r),`$ (80)
$`J_1(r)`$ $`=`$ $`{\displaystyle \frac{1}{12(1r)^4}}(16r+3r^2+2r^36r^2\mathrm{ln}r),`$ (81)
$`J_2(r)`$ $`=`$ $`{\displaystyle \frac{1}{2(1r)^3}}(1r^2+2r\mathrm{ln}r),`$ (82)
Finally,

$`C_{iAX}^R`$ $`=`$ $`g_2(O_R)_{A1}U_{X,i}^\nu ,`$ (83)
$`C_{iAX}^L`$ $`=`$ $`g_2{\displaystyle \frac{m_{l_i}}{\sqrt{2}m_W\mathrm{cos}\beta }}(O_L)_{A2}U_{X,i}^\nu ,`$ (84)
where $`A=1,2`$ and $`X=1,2,3`$, and

$`N_{iAX}^R`$ $`=`$ $`{\displaystyle \frac{g_2}{\sqrt{2}}}\{[(O_N)_{A2}(O_N)_{A1}\mathrm{tan}\theta _W]U_{X,i}^l+{\displaystyle \frac{m_{l_i}}{m_W\mathrm{cos}\beta }}(O_N)_{A3}U_{X,i+3}^l\},`$ (85)
$`N_{iAX}^L`$ $`=`$ $`{\displaystyle \frac{g_2}{\sqrt{2}}}\{{\displaystyle \frac{m_{l_i}}{m_W\mathrm{cos}\beta }}(O_N)_{A3}U_{X,i}^l+2(O_N)_{A1}\mathrm{tan}\theta _WU_{X,i+3}^l\},`$ (86)
where $`A=1,\mathrm{},4`$ and $`X=1,\mathrm{},6`$. The different matrices in these equations are defined in Appendix B of ref..