Decay of the Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} gauge boson with lepton flavor violation

The flavor-violating decay of a new neutral massive gauge boson Z′→μe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime \rightarrow \mu e$$\end{document} is analyzed in the context of extended models, in which this particle emerges. By means of the analysis of the μ→eγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \rightarrow e\gamma $$\end{document} decay, μ-e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu -e$$\end{document} conversion process in nuclei and the μ→ee+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \rightarrow e e^{+}e^{-}$$\end{document} decay, the strength of the Z′μe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime \mu e$$\end{document} coupling is estimated and used to calculate the branching ratio of the Z′→μe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime \rightarrow \mu e$$\end{document} decay. This is done for the so-called ZS,ZLR,Zχ,Zψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_S,\,Z_{LR},\,Z_\chi , Z_\psi $$\end{document} and Zη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_\eta $$\end{document} bosons. We found that, through the μ-e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu -e$$\end{document} conversion process, the most restrictive bound for the coupling is provided by the ZLR\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_{LR}$$\end{document} boson. Meanwhile, by means of the μ→ee+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \rightarrow e e^{+}e^{-}$$\end{document} decay, the most restrictive bound for the Z′μe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime \mu e$$\end{document} coupling is provided by the Zχ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_\chi $$\end{document} boson. However, if we concentrate on the less restrictive prediction for the Br(Z′→μe\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime \rightarrow \mu e$$\end{document}), this comes from the Zη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_\eta $$\end{document} boson and the resulting branching ratio is less than 10-4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-4}$$\end{document}. On the other hand, if we consider the most restrictive bound, the branching ratio for the process is below 2×10-7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\times 10^{-7}$$\end{document}, which results from the ZRL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_{RL}$$\end{document} boson, and is obtained through the μ-e\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu -e$$\end{document} conversion process.


I. INTRODUCTION
It is well known that Flavor Changing Neutral Currents (FCNC) produce very suppressed flavor transition processes in the quark sector of the Standard Model (SM).Additionally to a factor of suppression due to the GIM mechanism, these transitions are induced at one-loop level [1].Beyond the standard model (SM), there are many types of new physics (NP) that can lead FCNC, even at the tree level.These types of NP become important since associated transitions could be enhanced at tree level and consequently competitive with respect to the corresponding ones in the SM at the one-loop level.Moreover, within the NP context, it could arise effects that are definitively forbidden in the SM [2][3][4].
It draws the attention that even with minimal extensions to the SM, for example, by including right-handed components of neutrinos fields into it, some non-trivial physical properties can be explained, such as the neutrino oscillations, which preserve lepton number only at the global level [2][3][4].Nevertheless, within the frame of this extension, some predictions result quite suppressed and the theory requires the use of more elaborated extensions that could include lepton flavor violation in its Lagrangian.Concretely, the branching ratio Br(τ → µγ) < 10 −40 , that is calculated with minimal extensions to the SM [5], is very suppressed with respect to the experimental bound Br(τ → µγ) < 10 −8 reported in the Particle Data Group (PDG) [6].
An interesting process that violates the flavor symmetry in the lepton sector that has been studied in different scenarios is the Z → li l j (l i = l j = e, µ, τ ) decay [7].Although some branching ratios of these decays are more suppressed with respect to the corresponding experimental bounds, others are of the same order of magnitude [6].Similar flavor-violating decays involving a new neutral massive gauge boson, denoted as Z ′ [8][9][10][11][12][13][14][15][16][17][18], show that the branching ratio of these processes could be less suppressed than the respective decays of the Z boson.Indeed, in Ref. [19] the flavor-violating decay Z ′ → τ µ has been studied in models where the so-called Z S , Z LR , Z χ , Z ψ and Z η bosons arise.It was found that the most restrictive branching ratio for the process in question is less than 10 −2 .
On the other hand, the CMS and ATLAS [20,21] collaborations have been searching for new heavy particles via lepton flavor violation, in particular, the Z ′ boson decaying into pair of leptons µe, τ e, τ µ.Especially, the CMS [20] collaboration performed an analysis on the results of the cross-section and the branching ratios of the Z ′ boson decaying into µe in proton-proton collisions at the CM energy √ s = 13 TeV with an integrated luminosity of 35.9 fb −1 ; no evidence on physics beyond the SM related to the µe spectrum was found.
The ATLAS [21] collaboration studied the cross-section times the branching ratio of the Z ′ → µe, τ e, τ µ decays as a function of the Z ′ boson mass, in proton-proton collisions with an integrated luminosity of 36.1 fb −1 .As before, none evidence on excess over the SM predictions was found.The CMS Collaboration concludes that a Z ′ boson with a 10% branching ratio to the µe channel is excluded for masses below 4.4 TeV [20].While, the ATLAS Collaboration concludes that, from the eµ, eτ and µτ final states, Bayesian lower limits at 95 % credibility level on the mass of a Z ′ gauge boson with lepton flavor-violating (LFV) couplings are set at 4.5, 3.7, and 3.5 TeV, respectively [21].However, to completely discard these lepton flavor-violating processes, it is necessary to increase the confidence level on the measurements of such decays as well as the mass search range of the Z ′ gauge boson.
Hence, there is still an open window for theoretical analysis.In fact, a calculation of these flavor-violating decays, in a broad mass interval could be helpful to exclude or leave open the possibility of the existence of the Z ′ gauge boson.
In this work, we study the lepton flavor violation through the Z ′ → µe decay process, by using several models that predict a Z ′ gauge boson.This is done by using the Package-X software [22] that facilitates the calculation of the Feynman diagrams involved in the process, since it is possible to obtain analytical expressions for the amplitudes under analysis.The process in question is analyzed and bounded by using the µ → eγ decay, the µ −e conversion process in nuclei and the µ → ee + e − decay [6], which are useful for extracting bounds of the Z ′ µe coupling and then to constrain the decay under study.In fact, our approach introduces the µe decay mode promoted by the presence of the Z ′ boson as a feasible possibility to be considered in direct searches.However, complementarily, an experimental spin analysis would be required in order to establish or discard it.
The organization of this paper is as follows: Section II describes the flavor-violating Lagrangian in the context of extended models.In section III, we analyze the µ → eγ decay mediated by a Z ′ boson along with the µ − e conversion process in nuclei and the µ → ee + e − decay to estimate the strength of the Z ′ µe coupling.In section IV, we present a discussion on the Z ′ → µe decay.Finally, in section V the conclusions are presented.

A. The different models
There are several extended models that predict the presence of new massive neutral gauge bosons, commonly denoted as Z ′ .Here, we describe briefly the models used in this work.The grand unified group E 6 [8][9][10] breaks directly: E 6 → SO(10) × U(1) ψ and SO(10) → SU(5) × U(1) χ .In general, the E 6 gauge symmetry group predicts at most four neutral gauge bosons, however, we only consider the Z ψ boson and the Z χ gauge boson that arises from symmetry breaking of the group SO (10).On the other hand, Hewett and Rizzo [11] established that the group U(1) ψ × U(1) χ can be reduced to an effective group , where U(1) θ is a linear combination of U(1) ψ and U(1) χ .
The gauge fields Z ψ and Z χ corresponding to U(1) ψ and U(1) χ symmetries respectively, are massive and the mass eigenstates are defined as [11,23,24] By varying θ, it can be identified different bosons Z ′ , for instance: the Z χ boson with θ = 0, the Z ψ boson with θ = π/2.In particular, with θ = tan −1 (− 5/3), the Z η boson results as a linear combination of the Z χ and Z ψ gauge bosons, given explicitly Z η = 3 8 Z χ − 5 8 Z ψ .These three bosons arise in many superstring-inspired models in which the E 6 group breaks directly to a rank-5 group [9,10].At the last stage of the breaking symmetry of grand unification theories, other that SU(5), the breaking pattern generally finishes with at least one neutral gauge boson Z ′ , it is expected that its mass be of order of TeV's.This situation can be described by the theory based on the symmetry group SU(2) L × U(1) × U ′ (1).Following Refs.[15,25], in the gauge eigenstate basis, the neutral current Lagrangian can be written as where J µ EM is the electromagnetic current, Z 1 is the SU(2) × U(1) neutral gauge boson, and Z 2 is the new gauge boson associated with the additional Abelian symmetry U ′ (1).The g 1,2 couplings are of the form where g 1 = g/ cos θ W , being g the weak coupling, and θ W the Weinberg angle and λ g depends of the symmetry breaking pattern, which is usually assumed O(1) [8].For simplicity, it can be assumed that there is no mixing between Z 1 and Z 2 , and consequently they are the mass eigenstates of Z and Z ′ , respectively.The current associated with the group U ′ (1) can be cast as where i runs over all quarks and leptons and P L,R = 1 2 (1 ∓ γ 5 ) are the quiral projectors.The quiral couplings are ǫ f iL,R = Q f i L,R where the various Q ′ s are the quiral charges.For the case of the Z ψ , Z χ and Z η bosons, the chiral charges were determined in [14].
The Z LR gauge boson emerges in the left-right symmetric model in a scheme of breaking symmetry just as SO (10) [25][26][27][28][29]. Let us consider the model in which Z LR is orthogonal to the Z 1 in Eq. (2) that couples to the "left-right" current where the currents J µ 3R and J µ B−L are associated with the third component of SU(2) R and B − L symmetry, respectively; B(L) denotes the barionic (leptonic) number and α = (1 − 2 sin 2 θ W )/ sin 2 θ W ≈ 1.53.The J µ 3R current is constructed in such a way that all the right-handed fermions are doublets and all the left-handed fermions are singlets.In this way, the chiral couplings are given by [9,14]: where T f i 3R is the third component of its right-handed isospin in the SU(2) R group.The sequential Z S boson is defined to have the same couplings to fermions as the SM Z boson.In this model, g 2 = g 1 .It is useful as reference when comparing constraints from different sources.
Finally, the chiral charges Q f i L,R of the different models aforementioned are obtained and discussed in Refs.[9,14,15,26,29]; these are summarized in Table I.

B. The generic Lagrangian
In order to study the Z ′ → µe decay, we use of the most general renormalizable Lagrangian that includes flavor violation mediated by the Z ′ gauge boson.Such a Lagrangian can be written as follows [25,27,30]: where f i is any fermion of the SM and Z ′ α is the neutral massive gauge boson predicted by several extensions of the SM [8,9,[12][13][14][15][16]. Thus, the Ω Ll i l j and Ω Rl i l j parameters represent the strength of the Z ′ l i l j coupling, where l i is any lepton of the SM.We will assume that the couplings are symmetric: Ω R,L l i l j = Ω R,L l j l i .It is convenient to express the diagonal components of the Ω parameter in terms of the quiral charges as follows: Ω , where the various Q f i L,R are the quiral charges, shown in Table I.The non-diagonal the Ω Ll i l j and Ω Rl i l j parameters are to be somehow determined.

III. THE PROCESSES A. The µ → eγ decay
The Feynman diagrams corresponding to this decay are shown in Fig. 1, where the fermionic internal line represents a muon or an electron.In order to keep dominant contributions, we only consider diagrams involving just one vertex that induces lepton flavor violation.Specifically, the associated Feynman rules with the Z ′ ee and Z ′ µe vertices are , respectively.This vertex rule is denoted in the Feynman diagram with a dot.If we consider two vertices inducing lepton flavor violation, then it is necessary to consider loop diagrams with internal fermionic lines that corresponds to a tau lepton.Thus, the amplitudes corresponding to these diagrams are proportional to Ω µτ Ω τ e or to Ω 2 µe , since Ω µτ Ω τ e ≈ Ω 2 µe [31].The resulting amplitudes are suppressed by a factor Ω µe (∼ 10 −3 − 10 −2 in the interval m Z ′ = [2, 7] TeV) with respect to the amplitudes obtained by using only one lepton flavor violation vertex and they can be neglected.Notice that the calculation of the amplitude is carried out in the context of the unitary gauge.Thus, we find that the total amplitude can be expressed as where being m µ the muon lepton mass.As usual, e represents the electric charge of the electron.
The form factors F 1 and F 2 are explicitly given in the Appendix.It should be noted that the structure of these form factors was obtained using Package-X [22].The amplitude in Eq. ( 8) is free of ultraviolet divergences.According to the definition of decay width [6], we e e e e FIG.1: Feynman diagrams contributing to the µ → eγ decay.
obtain the branching ratio for the process in question where Γ µ is the total decay width of the muon lepton.

B. The µ − e conversion process
The µ − e conversion process refers to the capture of a muon by a heavy nucleus, in which the lepton flavor conservation is not preserved [32][33][34].Such processes is represented by where the nucleus N has mass and atomic number A and Z, respectively.In order to obtain a useful branching ratio for this process, we focus on the case of neutral current interactions with flavor violation at the tree level mediated by a Z ′ gauge boson (see Fig. 2).Specifically, we analyze the non-photonic process described the lepton-quark effective Lagrangian [34,35]  with the various parameters defined as After algebraic manipulations and introducing the nucleon spinor ψ N = ( p n ) and the Pauli matrix τ 3 , Eq. ( 13) can be rewritten as where the couplings C 1S and C 1V are denoted as the isoscalar vector and the isovector vector, respectively.These parameters are given explicitly as On the other hand, in Refs.[34,36] it has been shown that the couplings of axial quarks C 2S and C 2V do not contribute to the coherent nuclear charge (Q ′ W ), therefore the width of the µ − econversion process does not depend on these parameters.Thus, the width can be written as The most stringent bounds for capture rate are provided by gold and titanium nuclei [6].
The branching ratio for the µ − e [32] conversion is given by where Γ caputure ≡ Γ(µ − + N → capture) is the muon decay width.For our propose of bounding the Ω R,Lµe couplings, we use the branching ratio for the µ − e conversion, which, according to Eq. ( 18), can be expressed as with We consider the Z ′ gauge boson for the different models aforementioned.However, as it can be appreciated from Eqs. ( 20) and ( 21), the resulting branching ratio is proportional to the quiral charges of the quarks u and d, which vanishes for the Z χ , Z ψ and Z η bosons, since Q q L + Q q R = 0 [8,37].For the remaining bosons, Z S and Z LR , we take the following parameters of the titanium nucleus: Z ef f ≃ 17.6, F (q 2 ≃ −m 2 µ ) ≃ 0.54 and Γ capture = (2.59 ± 0.012) × 10 −6 s −1 [35,36,38,39], while the experimental bound for Br(µ − + N → e − +N) exp < 4.3×10 −12 [40].Eqs.(11) and (20) will be our starting point in next subsection to constrain the Ω R,Lµe couplings.
C. The µ → ee + e − decay The µ → ee + e − process can be directly calculated at the tree level.The corresponding Feynman diagram is shown in Fig. 3.By following Ref.[19], we found that Br(µ → ee + e − ) can be written as where In last equations, our approximations are valid for a =  22) and this simplifies the bounding process of the couplings.In this manner the final expression for the branching ratio is expressed as

D. Bounding the Ω µe coupling
As it was mentioned, we are interested in bounding the flavor-violating branching ratio of the Z ′ gauge boson decaying into the µ and e leptons.In order to do that, we firstly bound the |Ω R,Lµe | 2 couplings by using the branching ratios in Eqs.(11), (20) and (23).The task is carried out by considering each of the previous results as follows.

By using the µ → eγ process
To numerically constrain the couplings, we use Eq. ( 11) along with the experimental bound: Br(µ → eγ) Exp < 4.2 × 10 −13 with 90% C.L [6].Also, we take the following experimental values of the masses of the particles involved: m µ = 0.105658 GeV and m e = 0.00051099 GeV.In addition, we use the most recent experimental limits established for the mass of the various Z ′ gauge bosons, which, most of them, are larger than 3.9 TeV [17,18].
However, the mass of Z LR boson predicted by the left-right symmetric model is larger than 1.162 TeV [41].In this way, in next section we propose an analysis for the couplings between 2 and 7 TeV for m Z ′ .For the study, we consider three different cases: a) Vector-like coupling.For this case, we have that Ω Lµe = Ω Rµe ≡ Ω µe .Therefore, the bound can be expressed as where we have taken the lower bound of Br(µ → eγ) exp in accordance to Eq.( 11).In Fig. 4 (a) it is shown the maxima values of |Ω µe | 2 as a function of m Z ′ for the different models.As it can be observed, the most restrictive bound corresponds to the Z χ boson, while the less restrictive corresponds to the Z η boson, along the mass interval studied.b) Maximal parity violation.For this case, we take Ω Lµe = 0. Thus the corresponding bounding can be written as In Fig.
The maxima values allowed for the sum |Ω Lµe | 2 + |Ω Rµe | 2 are depicted in Fig. 4 (c).As it can be observed, the most restrictive bound is achieved by the Z χ boson.
Finally, let us mention that the resulting growing behavior of square magnitude of the couplings in the mass interval under analysis is consistent with the perturbative regime employed to carry out the study of the Z ′ boson decay.

By using the µ − e conversion process
Let us now use the µ − e conversion branching ratio in Eq. ( 20) to analyze the behavior of the couplings.For the study, we consider the same three cases above discussed: a) The vector-like coupling.For this case: Ω Lµe = Ω Rµe ≡ Ω µe .Therefore, b) The coupling with maximal parity violation.For this case Ω Lµe = 0. Thus b) Maximal parity violation coupling: c) Complete case: In Fig. 6, it is shown the behavior of bounds in question.As it can be observed, the most restrictive bounds for the couplings come from the Z χ boson in the three cases.On the contrary, the less restrictive bounds are provided by the Z η boson.In general, the case of vector-like coupling offers the most restrictive bounds.Although, as the bounds arising from the Z S and Z LR gauge bosons are concerned, the most restrictive ones result from the µ − e conversion process, as it can be appreciated from Figs. 5 (a) and 6 (a).decay.In (a), vector-like coupling case.In (b), maximal parity violation case.In (c), the complete case.
IV. THE Z ′ → µe DECAY The Feynman rules for this process can be obtained from the Lagrangian in (7).By following the usual procedures, we obtain that the general branching ratio for the lepton flavor violating process Z ′ → l i l j can be written as where Γ Z ′ is the total decay width of the Z ′ boson, and it was considered that m µ /m Z ′ ≪ 1.On account of Γ Z ′ , we consider all possible decay modes that include both flavorconserving and flavor-violating decays [15,19,42], that is to say, the decay of Z ′ to ν e νe , ν µ νµ , ν τ ντ , eē, µμ, τ τ , uū, cc, t t, d d, ss, b b, ūc + uc, tc + tc, τ µ + τ μ, τ e + τ ē, and μe + ēµ.
Notice that the total decay width Γ Z ′ is model dependent, since each decay mode of this particle is also model dependent.Previously, we have analyzed three different cases to bounding the couplings, through the three process above discussed, which will be used in the numerical analysis of the Z ′ → µe decay.For the discussion, we restrict ourselves to the vector-like coupling case, which gives the most restrictive bounds for the decay in question.
As matter of fact, the others cases provide bounds of the same order of magnitude for each m Z ′ , along the mass interval.Due to the values of the chiral charges (which make the branching ratio in (32) a non-trivial function of these charges), and the Ω couplings, the bounds for the branching ratio can change up to three orders of magnitude, as it is discussed below.
In Fig. 7 (a) it is shown the bounds of the branching ratio for the Z ′ → µe decay, calculated by using the µ → eγ decay.As we can see, the Z η boson provides the less restrictive bound for the decay, with values less than 10 −4 almost over the entire interval [2,7] TeV.The most suppressed bound corresponds to the Z S boson, which gives values for Br(Z ′ → µe) < 10 −5 for masses below 5 TeV, reaching the value of 1.32 × 10 −4 at m Z ′ = 7 TeV.For the rest of the models, the corresponding branching ratios vary between 10 −5 and 10 −4 .Let us now discuss the results of bounds for the branching ratio by using the predictions for the Ω µe coupling coming from the µ − e conversion process.In Fig. 7 (b), the behavior of such bounds is shown as a function of m Z ′ .As it can be observed, the resulting bounds for the Z S and Z LR bosons are more restrictive than the corresponding ones obtained through the µ → eγ decay (shown in Fig. 7 (a)), since the less restrictive result for the branching ratio is < 10 −6 , even in the most conservative scenario, which is obtained with the Z S boson.
As far as the bounds for Br(Z ′ → µe) resulting from the µ → ee + e − decay is concerned, the results are shown in Fig. 7 (c).The numerical values for the bounds show that the corresponding branching ratios are more suppressed than the obtained through the µ → eγ decay by two order of magnitude for each Z ′ boson.Nevertheless, the only non-vanishing branching ratios (for Z S and Z LR bosons), obtained by means of the µ − e conversion process, result more suppressed along the mass interval analyzed than the obtained by using the µ → eγ decay or the µ → ee + e − decay.
Complementarily, we also present contour plots for the Br(Z ′ → µe) as function of m Z′ and |Ω| 2 (see Figs. 8, 9 and 10) only for the vector-like case, since it provides the most stringent predictions.In these plots (one for each model and bounding approach) it is shown the excluded regions for specific branching ratios.In specific, the values of the branching ratio go from 10 −8 to 10 −6 .However, due to the property Q q L +Q q R = 0 of the chiral charges of quarks the only non-vanishing branching ratios for this process correspond to the Z S and Z LR bosons.For the case of the tree-body decay: µ → ee + e − , the resulting bounds for Br(Z ′ → µe) are more suppressed than the obtained in the cases above discussed, except for that related to the Z S and Z LR bosons.In fact, in this case the bounds of the branching ratio go from 10 −7 to 10 −6 .If we compare these bounds with the current experimental bound for the SM decay Br(Z → µe) < 7.5 × 10 −7 , the bounds for the Br(Z ′ → µe) result less restrictive or at least comparable., where x µ = mµ m Z ′ .

m 2 Z ′ m 2 µ>> 1 ,
which stand for m Z ′ > 1 TeV.Thus, the function h(a) is suppressed with respect to the function f (a) by at lest 8 orders of magnitude.So, we can can drop the term, h(

FIG. 4 :
FIG.4: Bounds for the Z ′ µe coupling as a function of m Z ′ , obtained through the µ → eγ decay.In (a), vector-like coupling case.In (b), maximal parity violation case.In (c), the complete case.
4 (b) it can be appreciated the behavior of maxima of |Ω Rµe | 2 as function of m Z ′ for the different models under consideration.As before, the most restrictive case corresponds precisely to the Z χ boson, and the less restrictive case is for the Z η boson.c) Complete case.This particular case is only achieved for the three different models that share the same absolute value of the left and right chiral charges |Q e L | = |Q e R |, namely, the Z χ , Z ψ and Z η bosons.For this case, we have ) c) Complete case.Here, the parameter under consideration is |Ω Lµe | 2 + |Ω Rµe | 2 .However, this case is the same as to the previous one, since it has the same expression in the right hand of the inequality in (28).In Figs. 5 (a) and (b), the upper bounds of the couplings in question as a function of m Z ′ are shown.As it can be observed in these figures, the most restrictive bounds are provided by the Z LR boson and they result less than 10 −6 almost along the full mass interval.

FIG. 5 :
FIG. 5: Bounds for the |Ω µe | 2 parameter as function of m Z ′ obtained through the µ − e conversion process.(a) For the case of vector-like coupling.(b) Maximal parity violation coupling and complete cases.

3 .
By using the µ → ee + e − decay Once again, we consider the three cases previously discussed.To accomplish the task of bounding the couplings, we use Eq.(23) along with the experimental bound Br(µ → ee + e − ) exp < 1.0 × 10 −12[6].The corresponding results are listed below.a) Vector-like coupling: The corresponding bound is

FIG. 7 :
FIG. 7: The branching ratios for the LFV decay of Z ′ → µe as a function of m Z ′ .The branching ratio was obtained by means of: (a) The µ → eγ decay.(b) The µ − e conversion process in nuclei.(c) The µ → ee + e − decay.

FIG. 8 :
FIG. 8: Contour plots for the Br(Z ′ → µe).The graphs were obtained by using the µ → eγ decay process.In (a) for the Z S boson, in (b) for the Z LR boson, in (c) for the Z χ boson, in (d) for the Z ψ boson, and in (e) for the Z η boson.In all the cases the thick line in color separates the excluded from the allowed region.

FIG. 9 :
FIG. 9: Contour plots for the Br(Z ′ → µe).The graphs were obtained by using the µ − e conversion process.In (a) for the Z S boson and in (b) for the Z LR boson.In both cases the thick line in color separates the excluded from the allowed region.

FIG. 10 :
FIG. 10: Contour plots for the Br(Z ′ → µe).The graphs were obtained by using the µ → ee + e − decay process.In (a) for the Z S boson, in (b) for the Z LR boson, in (c) for the Z χ boson, in (d) for the Z ψ boson, and in (e) for the Z η boson.In all the cases the thick line in color separates the excluded from the allowed region.
The Z ′ gauge boson is proposed in several extensions of the SM that induces FCNC at the tree level.This particle is under quest by the CMS and ATLAS experimental groups at LHC who have established experimental bounds on its mass.The existence of this particle could explain physical processes, such as the flavor-violating decays in the sector of leptons that are very suppressed in the SM, whose numerical values are not compatible with the respective experimental results.In this work, we have studied the flavor-violating branching ratio of decay of the Z ′ gauge boson into µe leptons, using different extended models and three different ways of bounding the decay.By resorting to the µ → eγ process, we found bounds for the Br(Z ′ → µe) as function of m Z ′ ; depending on the model, the values go from 10 −6 to 10 −4 , along the mass interval analyzed.If instead, the µ − e conversion process is used, we obtain more restrictive bounds for the branching ratio of the Z ′ → µe decay.

ACKNOWLEDGMENTS
This work has been partially supported by CIC-UMSNH and CONAHCYT, México.APPENDIXForm Factors: