Exact one-loop results for $l_i \to l_j\gamma$ in 3-3-1 models

We investigate the decays $l_i\rightarrow l_j \gamma$, with $l_i=e,\mu,\tau$ in a general class of 3-3-1 models with heavy exotic leptons with arbitrary electric charges. We present full and exact analytical results keeping external lepton masses. As a by product, we perform numerical comparisons between exact results and approximate ones where the external lepton masses are neglected. As expected, we found that branching fractions can reach the current experimental limits if mixings and mass differences of the exotic leptons are large enough. We also found unexpectedly that, depending on the parameter values, there can be huge destructive interference between the gauge and Higgs contributions when the gauge bosons connecting the Standard Model leptons to the exotic leptons are light enough. This mechanism should be taken into account when using experimental constraints on the branching fractions to exclude the parameter space of the model.


I. INTRODUCTION
The discovery of flavor neutrino oscillations (see Ref. [1] and the references therein) proves that neutrinos are massive. This leads to an important consequence that the lepton-flavor number violating decay µ → eγ is non-vanishing, being proportional to the neutrino masses and the mixing matrix. Assuming tiny neutrino masses satisfying current experimental constraints [1], extension of the Standard Model (SM) with right-handed neutrinos predicts that the branching ratio is Br ≈ 10 −55 , which will be called the SM contribution from now on. Meanwhile, the current experimental limits read [1] Br(µ − → e − γ) < 4.2 × 10 −13 , From theoretical side, the processes l i → l j γ are loop induced. Given that the SM contribution is strongly suppressed, they can be good places to look for new physics. In this paper, we consider a simple extension of the SM using the local gauge group of SU (3) C ⊗ SU (3) L ⊗ U (1) X (3-3-1) with new exotic leptons. The word exotic here means that they can have arbitrary electric charges and arbitrary masses. In this model, the electron (and similarly for muon and tauon) together with a neutrino and a new exotic lepton are in a triplet (or anti-triplet) representation of SU (3) L . In this work, we calculate both neutrino and exotic-lepton contributions, with special attention to the latter because the former is numerically suppressed as above mentioned.
We remark that 3-3-1 model is an active field of research and has a long history, see Ref. [2] and references therein. In this work, we choose a general class of 3-3-1 models, which are similar to the models presented in Refs. [2][3][4] where new heavy leptons are introduced.
However, there is an important difference: instead of fixing the electric charges of the new leptons to specific values being 0, +1 or −1, we let them be arbitrary. We will then study the dependence of the l i → l j γ branching fractions on this arbitrary charge. That class of 3-3-1 models has been studied in many works, see e.g. [5,6]. If we replace the new leptons with charge-conjugated partners of the SM leptons, we will have different 3-3-1 models with lepton-number violation, see e.g. Refs. [7][8][9][10][11]. The decays l i → l j γ in these models have been discussed in Refs. [12,13], see also the recent review [14] and references therein. We do not discuss these types of models in this work, but rather focus on the case with exotic leptons.
In the general class of 3-3-1 models here considered, there is one important parameter usually called β, which together with X, the new charge corresponding to the group U (1) X , define the electric-charge operator. The electric charges of new particles therefore depend on β. It has been known and widely accepted that β is one of the most important parameters to classify 3-3-1 models.
Recently, there occurred new efforts using 3-3-1 models to understand tensions between experimental measurements and the SM results in B physics, see e.g. [15,16]. Motivated by this work, we want to use 3-3-1 models to understand the l i → l j γ decays. Since the new leptons are assumed to be heavy, we expect large branching fractions. However, this is not totally obvious, because there are two contributions from gauge and Higgs sectors. Does a destructive interference effect occur?
The aim of this paper is manifold. First, we calculate the full and exact result for l i → l j γ partial decay widths for a general class of 3-3-1 models with arbitrary β. As a by product, we will perform numerical comparisons between the exact results (i.e. external lepton masses are kept) and approximate ones where external lepton masses are neglected. We note that approximate results have been almost exclusively used in the literature for the SM and many other models. We found this uncomfortable because the neutrino masses, which are much smaller than the lepton masses, are kept. We therefore want to know to what accuracy the approximate results valid, using the SM with arbitrary neutrino masses to answer this. As far as we know, this important point has never been addressed in the literature. We will also perform numerical studies for 3-3-1 model to see whether destructive interference effects occur and to see the dependence on β, gauge boson and Higss masses. To the best of our knowledge, this is the first study of l i → l j γ in 3-3-1 models with exotic leptons.
The paper is organized as follows. In the next section, we review the model and calculate the Feynman rules needed for l i → l j γ decays. We then summarize the main calculation steps and present analytical results in Section III. Numerical results are discussed in Section IV.
In Section IV A we perform comparisons between the approximate and exact results for the neutrino contribution. In Section IV B we present results for the exotic-lepton contribution.
Conclusions are in Section V. Finally, we provide two Appendices A and B to complete the results of Section III.

II. 3-3-1 MODEL WITH ARBITRARY β
One important condition we require is that the 3-3-1 model has to match the SM at the energy of the EW scale, about 250 GeV. This means that the SU (3) L symmetry is valid at a higher energy scale and is spontaneously broken down to the SU (2) L symmetry using the Brout-Englert-Higgs mechanism. In order to match the fermion representation of the SM, the simplest choice is to assign fermions into triplets and anti-triplets of the SU (3) L group.
However, this requires new fermions. In general, the electric charges of these new fermions are unkown. They, however, cannot be totally arbitrary because of the symmetry and of the matching condition with the SM. In most general terms, the electric charge operator can be written as where we have introduced the SU (3) generators T 3 , T 8 . Thus, the charge operator Q depends on two parameters β and X. With this information, we can write down the lepton representation as follows. Left-handed leptons are assigned to anti-triplets and right-handed leptons to singlets: The model includes three RH neutrinos ν aR and exotic leptons E a L,R which are much heavier than the normal leptons. The prime denotes flavor states to be distinguished with mass eigenstates introduced later. The numbers in the parentheses are to label the representation of SU (3) L ⊗ U (1) X group. For singlets, we have Q = X and hence the electric charges of the new leptons can be read off from the above information. The quark sector is not here specified since it is irrelevant to our present work.
We now discuss gauge and Higgs interactions. There are totally 9 EW gauge bosons, included in the following covariant derivative where T 9 = 1/ √ 6, g and g X are coupling constants corresponding to the two groups SU (3) L and U (1) X , respectively. The matrix W a T a , where T a = λ a /2 corresponding to a triplet representation, can be written as where we have defined the mass eigenstates of the charged gauge bosons as From Eq. (2), the electric charges of the gauge bosons are calculated as We note that B is also the electric charge of the new leptons E a .
To generate masses for gauge bosons and fermions, we need three scalar triplets. They are defined as where A, B denote electric charges as defined in Eq. (7). These Higgses develop vacuum expectation values (VEV) defined as 5 The symmetry breaking happens in two steps: It is therefore reasonable to assume that u > v, v . After the first step, five gauge bosons will be massive and the remaining four massless gauge bosons can be identified with the before-symmetry-breaking SM gauge bosons. This leads to the following matching condition for the couplings where g 2 and g 1 are the two couplings of the SM corresponding to SU (2) L and U (1) Y , respectively. From this we get the following important equation which helps to constrain β: where the weak mixing angle is defined as t W = tan θ W = g 1 /g 2 and we denote s W = sin θ W .
Putting in the value of s W , we get approximately which will be used in the numerical analysis.
The masses of the charged gauge bosons are We now discuss the mixings of leptons. In general, the mixing between a SM lepton and a new lepton is allowed if they have the same electric charge. However, since we consider a general class of models with arbitrary β, this mixing effect will be neglected.
This is justified because we will assume that the new leptons are much heavier than the SM leptons. Therefore, only generation mixings as in the SM are allowed. The Yukawa Lagrangian related to these mixings reads where a, b = e, µ, τ are family indices. The corresponding mass terms are: From now on we will work in the basis where the SM charged leptons are in their mass eigenstates. This can always be done without loss of generality. We can therefore set Y e ab to be diagonal and e = e in Eqs. (14,15). The transformations from the flavor states to mass eigenstates are defined as For the Higgs sector, we assume A, B = 0, ±1, so that only the following mixings of scalar fields with the same electric charge are allowed, namely (χ A , η A ), (χ B , ρ B ), (ρ + , η + ) and (χ 0 , ρ 0 , η 0 ). The neutral components are expanded as: The ratios between VEVs are used to define three mixing angles: We will also use the following notation The scalar potential is With the above notation, the mass eigenstates are  are the Goldstone bosons of W ± , Y ±A and V ±B , respectively. The masses of the charged Higgs bosons are The neutral Higgs bosons are not involved in our calculation; hence they have been ignored.
In total, there are six charged Higgs bosons, one neutral pseudoscalar Higgs and three neutral scalar Higgses. Bosonic particles with electric charges of ±B are not involved in the present calculation. Nevertheless, their masses and mixing angles are above provided for the sake of completeness.
From the above information we can obtain all vertices needed for the calculation of l i → l j γ decays. They are listed in Table I.

Vertex
Coupling Vertex Coupling Other notations are defined in the text.

III. ANALYTICAL RESULTS
Equipped with the above Feynman rules, we can proceed to calculate the partial decay width of l 1 → l 2 γ using standard techniques of one-loop calculation. We have done this in a careful way, with at least two independent calculations, and paid special attention to the relative sign between the gauge and Higgs contributions. This relative sign is very important because, as we will see in the numerical results, the interference term can be positive or negative.
In the literature, the calculation of l 1 → l 2 γ is usually done by neglecting the external lepton masses. As stated in the introduction, we found this uneasy because the neutrino masses, which are much smaller than the lepton masses, are kept. We therefore want to check the validity of this approximation. To achieve this we have to keep the external lepton masses.
We have calculated the partial decay width of l 1 → l 2 γ from scratch without approximation. In the following we summarize the key points and present exact analytical results.
Results for the SM case are obtained as a special case and are discussed in Section IV A.
We consider the process where p 1 = p 2 + q and the helicity indices have been omitted for simplicity. The amplitude where λ is the photon's polarization vector, Γ λ are 4 × 4 matrices depending on the gamma matrices, external momenta and coupling constants. After requiring the general conditions that the spinors obey the Dirac equations, q µ µ = 0, and q λū 1 (p 1 )Γ λ u 2 (p 2 ) = 0, we can prove that the amplitude depends on only two form factors as where C L,R are called form factors, P L = (1 − γ 5 )/2, P R = (1 + γ 5 )/2. The partial decay width is then written as This result is well known and has been given in e.g. Ref. [17].
Since we assume that the exotic leptons are much heavier than the SM leptons, the branching fractions of the dominant decays of l 1 → l 2ν2 ν 1 in the 3-3-1 model are the same as those of the SM. Using the well-known tree-level result of Γ(l 1 → l 2ν2 ν 1 ) = G 2 F m 5 1 /(192π 3 ) (see e.g. [18]), where G F is the Fermi coupling constant, we write the branching fraction as where G F = g 2 /(4 √ 2m 2 W ) and we have defined C L,R = m 1 D L,R and the approximation m 2 m 1 has been used for the first factor, but not for D L,R . For later numerical analysis we will use Br(µ → eν e ν µ ) = 100%, Br(τ → eν e ν τ ) = 17.82% and Br(τ → µν µ ν τ ) = 17.39% as given in Ref. [1]. It is noted that (1), and hence, in the approximation This point is important to understand the approximate results discussed in the next sections.
The next step is to calculate D L,R for the 3-3-1 model with arbitrary beta presented in the previous section. Representative Feynman diagrams are shown in Fig. 1. Using the Feynman rules in Table I and summing over all possible Feynman diagrams, we obtain the following results where Finally, we make an important remark on the dependence on coupling constants. From Eq. (30) we have D L,R ∝ eg 2 . Using Eq. (28) and noticing that G F = g 2 /(4 √ 2m 2 W ), we get Br(l 1 → l 2 γ) ∝ e 2 , being independent of g or s W . Clearly, the coupling constant e = √ 4πα should be calculated in the low-energy limit for the processes at hand. Therefore, we will use α(0) as input parameter in our numerical analyses.

IV. NUMERICAL RESULTS
Input parameters are specified as follows. We use, according to Ref.
The neutrino mixing matrix is assumed to be real and is calculated from the above mixing angles as where c ij = cos θ ij , s ij = sin θ ij with i, j = 1, 2, 3.
A. Neutrino contribution: approximate vs. exact The approximate results calculated by neglecting the external lepton masses have been exclusively used in the literature. However, the justification is not totally obvious to us because the neutrino masses, which are much smaller than the lepton masses, are kept. We therefore present here compact formulas for the exact results (i.e. m 1 and m 2 kept) and perform a numerical comparison with the approximate ones.
The SM result includes only the W contribution and is given by D νW L,R . Using the formulas in Appendix A we write the result in terms of scalar one-loop integrals A 0 , B 0 and C 0 , which are also calculated in Appendix A. We obtain where This result was first obtained in Ref. [21] and has been widely used for any values of neutrino masses. We may wonder whether this is justified for the case of m νa m 1 or m νa ≈ m 1 .
This is the reason we perform a numerical comparison here between the exact result and the approximate one with m 1 = m 2 = 0 for many values of m ν 1 from zero to 10 16 GeV. The motivation is of purely mathematical nature and we ignore the physical constraints on the neutrino masses here. The results are shown in Table II and Fig. 2. We have used Eq. (28) to calculate the branching fractions for both cases. We see that the difference is less than permil level for µ → eγ and τ → eγ and is at the permil level for τ → µγ. This result is independent of neutrino masses.
We now take into account the charged Higgs contribution. There are two additional parameters t v v and m H ± (see the D νH + L,R terms in Eq. (30)). We have calculated the difference between the exact and approximate results for four cases of t v v = 1/50 or 50 (we choose  Table II. is very similar to the SM case: the difference is below permil level for µ → eγ and τ → eγ and is at the permil level for τ → µγ. For the absolute value of Br(µ → eγ) the result is We have a technical remark here. Due to the huge hierarchy among the neutrino, charged leptons and W boson masses, the numerical calculation of the exact result is non-trivial because of numerical cancellation. To obtain the µ → eγ results in Table II we have used Mathematica 9 with at least 62 precision digits for m ν 1 = 10 −13 GeV and about 180 precision digits for m ν 1 = 10 16 GeV.

B. Exotic-lepton contribution
In this numerical study we investigate the exotic-lepton contribution, to see how large the branching fractions can reach, what can be the dominant effects and dependence on the parameter β, m Y and m H A . We will also show the gauge-Higgs interference effects.
In the previous section we have shown that the neutrino contribution is well below the current experimental limit. We will therefore neglect the neutrino contribution including interference effects with exotic leptons in the following. The external lepton masses will be neglected as justified in Section IV A.
In the following we choose a benchmark point, which is a typical scenario where the where This means that the terms proportional to t v u in D EH A R in Eq. (30) can be safely neglected and the branching fractions are almost independent of t v v . We note that terms proportional to ct v u are suppressed because they are also proportional to the external lepton masses. We will therefore set t v v = 1 in the following. As a side note, for the choice of v = v there is another good justification: it makes the parameter ρ = m 2 W /(m 2 Z cos 2 θ W ) with θ W being the weak-mixing angle close to unity, as pointed out in Ref. [22] where the same scalar potential is used.
Other parameters related to the exotic leptons are unknown. We choose, as an example, the following default values for the remaining input parameters: The mixing matrix V L is calculated from three mixing angles θ E 12 , θ E 13 , and θ E 23 as in the case of neutrinos. The values of the exotic-lepton masses are chosen within the unitary bound of m E i < 16m Y A as derived from the partial wave unitarity of the E iĒi → E iĒi scattering [23].
A few remarks on the above default input-parameter choice are appropriate here. Concerning gauge bosons, the best ATLAS/CMS limits for 3-3-1 models with exotic leptons are summarized in Table III. We note that, in almost all cases, the contributions from exotic leptons to the Z total width are neglected, except for the case of Ref. [24] where m F = 1 TeV is assumed for all exotic fermions. When those contributions are properly taken into account, the bound on m Z will get weaker, because the branching fractions of Z → l + l − with l = e, µ will decrease. Therefore, the default choice in Eq. (35) may be acceptable. However, one should keep in mind that, strictly speaking, the ATLAS/CMS bound on m Z is unknown for our present numerical analysis, because it depends on the masses and electric charges of the exotic fermions (i.e. leptons and quarks) which have not been properly taken into account. We will therefore relax the constraint on m Y A , varying it from 0.5 to 3 TeV for some plots. In this context, it is noted that, using LEP II data, the authors of Ref. [15] obtained m Z   In the following tables and plots, if not otherwise stated, the above default values are used. Differently from Section IV A, we will use the true values of Br(l 1 → l 2ν2 ν 1 ) as given in the text below Eq. (28) so that one can compare the results in this section with current experimental limits.    In Table IV Table IV. For the decay of µ → eγ, only the cases of β = 0 or m Y A = 3 TeV are below the experimental limit of 4.2 × 10 −13 .
We now focus on the decay µ → eγ and discuss two density plots to see the dependence on β, m Y and m H A . In Fig. 3 we show the density plot of Br(µ → eγ) as a function of β and m Y (left) and of β and m H A (right). We observe from the left plot, consistently with Table IV, the branching fraction are smallest when β is around zero or when m Y is large.
From the right plot, we see a similar dependence on β, but the dependence on m H A is much weaker than on m Y . From those two plots, we conclude that large branching fraction occurs at large |β|, small m Y and small m H A .
In a series of six plots in Fig. 4 we would like to show again the dependence on β, m Y and m H A , but with two-dimensional plots this time and for all three decays. We see clearly that the case of β = 0 is special and different from the other cases of β = 1/ √ 3, contributions are non-negative. Overall, Fig. 5 shows that destructive interference effect tends to occur when |β| and m Y are small.

V. CONCLUSIONS
In this paper, we have provided full and exact analytical results for the l i → l j γ partial decay widths for a general class of 3-3-1 models with exotic leptons and with arbitrary β.
As a by product, we performed numerical comparisons between exact results (i.e. external lepton masses are kept) and approximate ones where m i = m j = 0. We conclude that, for either extremely light neutrinos or very heavy leptons, the difference between exact and approximate results is less than permil level for µ → eγ and τ → eγ and is at the permil level for τ → µγ. Therefore, unsurprisingly, approximation results widely used in the literature are excellently justified.
Concerning the exotic lepton contribution, we found huge destructive interference between the gauge and Higgs contributions. This can happen when |β| and m Y are small enough.
This has an interesting consequence: the branching fractions can be small even for small m Y . Therefore, this destructive interference mechanism must be taken into account when using experimental constraints on Br(l i → l j γ) to exclude parameter space. This in particular means that if one takes into account only the gauge contribution then the results can be completely off. It is likely that this destructive interference mechanism also occurs in b → sγ and other similar processes.
The masses of particles in the loop are written explicitly in the argument list and there is an one-to-one correspondence between those masses and the upper index of the h i (h stands for Higgs) and g i (g for gauge) functions.
The scalar functions A 0 , B 0 , C 0 can be calculated using the techniques of [29]. We have A 0 (m 2 ) = m 2 C U V − log(m 2 ) + 1 , B 0 (0, m 2 , m 2 ) = C U V − log(m 2 ), where C U V = 2/(4 − D) − γ E + log(4πµ 2 ) with γ E being Euler's constant and x σ and y iσ are the roots of the following equations For the case of p 2 1 > 0 and p 2 2 = 0 we have Results for the case of p 2 1 = p 2 2 = 0 have been provided in Ref. [17]. We finally note that the C functions in Eq. (A2) are independent of the auxiliary parameter C U V , meaning that the final results are UV finite. The function B