Lepton Flavor Violation and Collider Searches in a Type I + II Seesaw Model

Neutrino are massless in the Standard Model. The most popular mechanism to generate neutrino masses are the type I and type II seesaw, where right-handed neutrinos and a scalar triplet are augmented to the Standard Model, respectively. In this work, we discuss a model where a type I + II seesaw mechanism naturally arises via spontaneous symmetry breaking of an enlarged gauge group. Lepton flavor violation is a common feature in such setup and for this reason, we compute the model contribution to the $\mu \rightarrow e\gamma$ and $\mu \rightarrow 3e$ decays. Moreover, we explore the connection between the neutrino mass ordering and lepton flavor violation in perspective with the LHC, HL-LHC and HE-LHC sensitivities to the doubly charged scalar stemming from the Higgs triplet. Our results explicitly show the importance of searching for signs of lepton flavor violation in collider and muon decays. The conclusion about which probe yields stronger bounds depends strongly on the mass ordering adopted, the absolute neutrino masses and which much decay one considers. In the 1-5 TeV mass region of the doubly charged scalar, lepton flavor violation experiments and colliders offer orthogonal and complementary probes. Thus if a signal is observed in one of the two new physics searches, the other will be able to assess whether it stems from a seesaw framework.


Introduction
The observation of neutrino oscillations implies in non-zero neutrino masses. Their masses are much smaller than any other in the Standard Model (SM) spectrum. Several mechanism have surfaced trying to explain the smallness of the neutrino masses [1][2][3][4][5]. In the type I seesaw, the existence of heavy right-handed neutrinos is evoked, whereas in the type II seesaw a scalar triplet is added to the SM. This scalar triplet couples to the SM lepton doublets and features a neutral scalar that develops a small vacuum expectation value giving rise to tiny neutrino masses [6][7][8][9][10]. This mechanism can elegantly explain the masses of the three active neutrinos in the SM and lead to several phenomenological imprints in collider, low energy observables and leptogenesis [11,12,[12][13][14][15][16][17].
On the other side, we do not know why we have three generations of fermions in the SM. Theoretically speaking, it would nice to have a model where these two problems are simultaneously addressed. Several models have been proposed where the number of generations is addressed by anomaly cancellation and asymptotic freedom arguments, and they are known as 3-3-1 models [18][19][20]. These models have an enlarged gauge sector, SU (3) c × SU (3) L × U (1) N (3-3-1 for short). Notice that SU (2) L × U (1) Y in the SM gives place to SU (3) L × U (1) N . Therefore, the fermions will now be arranged in the fundamental representation of SU (3) L , i.e. triplets. The same happens for scalars fields, also sorted in triplets in order to generate fermion masses. Several versions based on the 3-3-1 symmetry have been proposed in the literature where the replication of fermion generations are explained [21][22][23][24][25].
Not all of these realizations are capable of explaining neutrino masses at the same time, though. In this work, we focus on a model which can account for neutrino masses, known as 3-3-1 model with right-handed neutrinos [20,[26][27][28]. Originally the model has only three scalar triplets. With scalar triplets one can generate two mass degenerate neutrinos and a massless one [29], which is in conflict with existing data [30]. The use of high dimensional effective operators has been put forth in the attempt to break the degeneracy [31,32]. In this setup, the smallness of the neutrino masses and the oscillation pattern is not successfully explained. The most simple way to nicely solve this issue is by adding a scalar sextet [33,34]. The interesting aspect of this scalar sextet is that after spontaneous symmetry breaking it breaks down to a scalar triplet, two doublet scalars and a scalar singlet field. The scalar triplet is exactly the one desired to perform the type II seesaw mechanism. Therefore, in this way, the type II seesaw arises as a result of the spontaneous symmetry breaking. Several phenomenological aspects of the 3-3-1 with right-handed neutrinos have been explored in the past [35][36][37][38][39][40], but our work differs from those because, • We explicitly compute the µ → eγ and µ → 3e decays; • We explore their connection to neutrino mass ordering; • We put our results into perspective with collider and lepton flavor violation bounds.
This work is structured as follows: In section II we briefly describe the model and explain how neutrino masses are incorporated, in section III we present the relevant collider bounds, in section IV we derive the lepton flavor violating muon decays and then draw our conclusions. In Appendix, we provide a details discussion of our reasoning and a general derivation of the µ → eγ decay.

Model with right-handed neutrinos and scalar sextet
There are several models based on the 3-3-1 gauge group. The model discussed here differs from the original proposal [18][19][20] because of the presence of three right-handed neutrinos. This model does not suffer from a Landau pole at the TeV scale as the previous version. For this reason, it has attracted lots of attention in the past decade [41][42][43][44][45][46][47][48][49][50][51][52]. In this section, we will briefly discuss the model to ease our reasoning. We start with the fermion content.

Fermion Content
In this model, the left-handed leptons are arranged in the fundamental representation of SU (3) L , whereas right-handed leptons as a singlet, as follows, In a similar vein, the quarks are arranged as follows, where α = 1, 2. We highlight that the first two generations of quarks are placed in the anti-triplet representation of SU (3) L . This has to be the case in order to cancel gauge anomalies. The quarks D 1,2 are down-type heavy quarks, in the sense, they have the same hypercharge as the SM quark down. A similar logic applies to the heavy up-quark U . These quarks have masses proportional to the scale of symmetry breaking of the 3-3-1 symmetry which should lie at several TeV to be consistent with the non-observation of such exotic quarks [53,54]. The scale of symmetry breaking of the 3-3-1 symmetry will be kept sufficiently high to be consistent with this result. We will address the fermion masses below.

Fermion Masses
To generate fermion masses we need to invoke three scalar triplets of the following form, (2.6) These scalar triplets are sufficient to successfully yield masses to all fermions, except neutrinos. The spontaneous symmetry breaking goes as: first In other words, when χ develops vacuum expectation values the 3-3-1 symmetry is broken reproducing the SM gauge group, which then breaks into quantum electrodynamics after η and ρ acquire a vacuum expectation value. For the purpose of this work, it easier to work on the broken phase where, For simplicity we will adopt v η = v ρ . This assumption is typically made to simplify the diagonalization of the mass matrices involving the scalars [43]. Moreover, we take v 2 η + v 2 ρ ∼ 246 2 GeV 2 , this ought to be enforced to generate the correct masses for the W and Z bosons in the SM, as occurs in models with extra scalars contributing to gauge boson masses [55,56].
We have explained in detail how each mass is generated in the Appendix, for this reason, in what follows we will just briefly consider each sector separately.

Charged Leptons Masses
The masses for the charged leptons are generated via the lagrangian, Notice that this term conserves lepton number. One possible term that one could write down is ijk ψ c L i ψ Lj ρ k , which does not conserve lepton flavor. The reason why we did not include this term in Eq.(2.8) is that there is a set of discrete symmetries that will invoke later on to prevent mixing between the SM quarks and the exotic ones. One of them requires ρ → −ρ, which forbids the term above. However, we need to impose e R → −e R to engender masses for charged leptons. When the field ρ 0 2 develops a vacuum expectation value, v ρ , we find, (2.9)

Quarks Masses
The quark masses could stem from two sources, one where flavor is conserved (LFC) and other where it is violated (LFV) as follows, and Nevertheless, the LFV lagrangian is problematic because it induces mixing between the SM quarks and the exotic ones. This mixing could alter the properties of the SM quarks. Thus we need to eliminate them. To do so, we invoke a set of Z 2 symmetries where some fields are odd under: u aR → −u aR , d aR → −d aR , D aR → −D aR , η → −η, ρ → −ρ. The reaming fields transform trivially. We emphasize that these discrete symmetries do not affect the other sectors of the model that is still capable of reproducing the SM features at low energy scales with no prejudice.
In summary, with these discrete symmetries the SM quarks get masses through Eq.(2.10) yielding, and, It is visible that their masses are proportional to v η and v ρ which are set at the electroweak scale to avoid flavor changing interactions and return the correct quark masses [45,[57][58][59][60][61][62][63].
Moreover, from Eq.(2.10) we find the exotic quarks masses, which are proportional to the scale of the 3-3-1 symmetry breaking, v χ , assumed to be sufficiently high to be in agreement with null results from collider searches for exotic quarks [54,64,65]. We will now turn our attention to neutrino masses.

Neutrino Masses-Type I + II Seesaw
As far as neutrino masses are concerned, the scalar triplets do not suffice to successfully generate neutrino masses. Using Lie algebra one may notice thatψ c L ψ L ∼ (3 * + 6, −2/3). Thus, one may generate neutrino masses via the scalar triplet ρ and a scalar sextet [20,66,67] with, (2.14) However, the discrete symmetries mentioned previously prohibit the Yukawa term proportional to ρ and for this reason, only the scalar sextet contributes leading to, Neutrino masses arise after the neutral scalars develop vacuum expectation value as follows, Notice that v s11 , v s13 , and Λ are the vacuum expectation value of the neutral scalars in the sextet.
With this vacuum structure we generate a mass matrix of the form, with, We remind the reader that M ν is a 6 × 6 matrix and M L , M D and M R as 3 × 3 matrices. After the diagonalization procedure we get for the active neutrinos, (2.22) and for the right-handed ones, It is good timing to highlight a few things: (i) the mass ratio m ν /m ν R is independent of the Yukawa coupling; (ii) a scalar sextet breaks down to a triplet plus doublet plus singlet scalar, i.e. 6 → 3 + 2 + 1. In order words, the scalar sextet generates the triplet scalar used in the type II seesaw, the doublet scalar that induces Dirac masses, M D , and a singlet scalar that yields right-handed Majorana masses, M R . Therefore, the scalar sextet naturally gives rise to a type I +II seesaw mechanism via spontaneous symmetry breaking; (iii) Setting v 2 s13 /Λ ∼ v s11 the smallness of the active neutrino masses are justified by taking v s11 to be around 1 eV.
There are different possibilities to successfully generate neutrino masses. Since the vacuum expectation values of the fields in the scalar sextet are in principle arbitrary one can play with them and find different Yukawa couplings leading to the same neutrino masses. Be that as it may, one can draw important and insightful conclusions by adopting some simplifications. We will assume throughout that v 2 s13 /Λ ∼ v s11 and v s13 = v ρ = v η . Therefore, if we choose v s11 = 1 eV this automatically translates into Λ ∼ 10 13 GeV, and v s11 = 100 eV yields Λ ∼ 10 11 GeV. Later on, we will present several benchmark scenarios taking v s11 = 1 eV and v s11 = 100 eV to investigate the impact of the neutrino mass ordering on the lepton flavor violating muon decays. With this information, one can now have an estimate of the right-handed neutrino masses since they are proportional to Λ.
For simplicity we will adopt that v 2 In this way, we have a dominant type II seesaw mechanism with the neutrino masses governed by the vacuum expectation value of S 0 11 . One can think of it as the vacuum expectation value of the triplet under SU (2) L in the broken phase in the usual type II seesaw study.
Looking at Eq.(2.15) we notice that S 0 33 , S 0 11 , S − 12 , S −− 22 carry two units of lepton number and for this reason are called bileptons. The lepton flavor (number) is violated after the neutral components S 0 11 and S 0 33 acquire a vacuum expectation value. The singly charged scalar in the scalar sextet will not mix with other singly charged scalars that do not carry lepton number. The presence of flavor violation will be explored in this work in the context of muon decays.

Gauge Sector
We turn our attention to the gauge sector. The main point is to show that there are new gauge bosons with masses are proportional to the 3-3-1 scale of symmetry breaking and they are subject to stringent collider bounds. Showing this in a pedagogical manner requires us to start with the covariant derivative of SU (3) L ⊗ U (1) N , where g L is the coupling constant of SU L (3) group, g N is the coupling constant of U (1) N , λ m are the Gell-Mann Matrices with m = 1, ..., 8, W m µ are the gauge bosons in the adjoint representation of SU (3) L , W N µ is the gauge field associated to U (1) N and N ϕ is the hypercharge associated to U (1) N . Writing down the term proportional to λ m one finds, Therefore, the 3-3-1 model with right-handed neutrinos predicts the existence of new charged gauge bosons, W ± , which are subject to intense searches at the LHC [68,69], and exotic neutral gauge bosons U 0 and U 0 † [70]. Moreover, from a combination of the W 8 , W 3 and W N fields, we extract the SM photon and Z bosons, as well as a massive Z field. It is clear the model add five new gauge bosons to the SM, as a direct result of the extended gauge sector which predicted N 2 − 1 bosons, where N = 3 in this case.
The masses of SM gauge bosons are slightly altered by the presence of the scalar sextet. This change is proportional to the vacuum expectation value v s11 which is meant to be small. A similar conclusion is found in the usual type II seesaw mechanism. The bound that rises from the ρ parameter enforces [55] The masses of the new gauge bosons are all proportional to the scale of symmetry breaking of the model. For a complete spectrum and equations of the gauge boson masses, we refer to [43]. Although, some relations are quite useful to ease our reasoning. One can find that, Hence, a lower mass bound on the Z boson represents a direct lower mass bound on the W , similarly to what occurs in the minimal left-right model [68,69]. One may wonder about the existence of important limits related collider searches for the U 0 gauge boson. Although, its mass is identical to the W which is subject to much more intensive searches. At the end of the day, the most important gauge boson as far as collider searches are concerned will be the Z field, as we explain further.

Searches for Z bosons
Collider searches for heavy Z bosons are quite popular because they typically feature a clear signal. If they couple to fermions and have a narrow width, they give rise to pronounced bumps in the dilepton or dijet invariant mass [71]. If the couplings to leptons are not very suppressed, the use of dilepton data is more promising because it is subject to a smaller SM background. In this model, the Z couplings to leptons are not small. Using LHC data at 13 TeV center-of-mass energy with 3.2 f b −1 of integrated luminosity the authors in [72] placed a lower mass bound of 3 TeV. Later, in [73] this limit was improved using 36 f b −1 and 3 ab −1 of integrated luminosity finding, (3.1) We emphasize that this limit of 4 TeV relies on the dielectron plus dimuon data with invariant mass in the 500 − 6000 GeV mass range as recommended, with the cuts in transverse energy and momentum as recommended by ATLAS collaboration in [74,75]. The projected bound of 6.4 TeV assumes a similar detector with the same trigger efficiency running at 14 TeV center-of-mass energy and with 3ab −1 of integrated luminosity. This lower mass bound of 4 TeV is rather robust and important because the Z is tied to the scale of symmetry breaking of the model. Since m Z = 0.4v χ (see Eq.(2.27)), we automatically find v χ > 10 TeV. Furthermore, using the mass relations between the W and Z bosons one finds a lower mass bound on the W mass which is much restrictive than other [76,77] on its mass. Anyway, the Z bound aforementioned is the most relevant to our discussion. We highlight that this limit has nothing to do with lepton flavor violation, it is rather simply based on searches for heavy dilepton resonances. Therefore, they are not relevant to our discussion because the mass of the scalars that enter in the lepton flavor violation discussion are not much sensitive to the energy scale of 3-3-1 symmetry breaking. Nevertheless, we emphasize that our entire discussion of lepton flavor violation is fully consistent with these bounds on the Z mass. Now we will address collider searches for lepton flavor violation in what follows.

Searches for a Doubly Charged Scalar at the LHC
The presence of a doubly charged scalar in the spectrum is typical signature of a type II seesaw mechanism [78][79][80][81][82][83][84][85][86]. Typically, this type II seesaw is realized by the addition of a scalar triplet. This popular extension triggered several phenomenological analyses. As we mentioned previously, the scalar triplet that arises after the spontaneous symmetry breaking is key to the type II seesaw mechanism in our model. We emphasize, however, that in our model we have a type I + II seesaw because we do also have Dirac neutrino masses. After spontaneous symmetry breaking, we can single out the Yukawa term involving the scalar triplet and SM particles which reads [55], where L is the SM lepton doublet with, The doubly charged scalar decay width into charged leptons is found to be [87], Thus is clear that the branching ratio into charged leptons can change depending on the Yukawa couplings. Different choices lead to different branching ratios and consequently different lower mass bounds. The dominant decay mode determines cuts, detector efficiency, and backgrounds which the signal is subject to, and consequently yielding different lower mass bounds. Dielectron and dimuon channels offer a cleaner environment and thus yield stronger bounds. If lepton flavor violation is assumed the SM background is suppressed, which again strengthens the limits. This is reasoning behind the limits derived using LHC data.
We used the code fastlim described in [88], and adopted a parton distribution function at next-to-next leading [89,90] in order to project the LHC sensitivity for a high-luminosity (HL) setup, following the recommendations presented in [91]. We stress that the HL-LHC limit refers to a detector similar to LHC running at 14 TeV with 3 ab −1 of integrated-luminosity. The High-Energy LHC configurations represents a 27 TeV colliding beam with 15 ab −1 of data. We highlight that these bounds are based on the simulated signal qq → Z, γ → φ ++ φ −− as outlined in [87], which features bounds stronger than previous previous studies [92][93][94]. In summary we derived, The latest LHC search for doubly charged scalars was with 12.9 f b −1 of data, but notice that with 36 f b −1 we expect LHC to already rule out doubly charged scalars with masses around 940 GeV, assuming a normal mass ordering for the active neutrinos. If we had considered an inverted mass ordering a weaker bound would have been found, lying around 900 GeV. This difference in the lower mass bound from LHC for normal and inverted ordering will not cause a meaningful impact on our conclusions and with this understanding in mind, we will simply quote those in Eq. (3.5). There are other important limits on this scenario [83,[95][96][97][98] but the ones we quote are the most relevant. For a more complete discussion of the collider bounds on the type II seesaw we refer to [99].
Moreover, it is exciting to see that HE-LHC can potentially probe doubly charged scalar with masses up to ∼ 5 TeV. These limits are quite important and serve as an orthogonal test to the type I + type II seesaw scenario we are investigating because that restricts the region in which the seesaw mechanism is viable. Having in mind that the doubly charged scalar is key to the lepton flavor violation observables we are about to discuss, such collider bounds stand as a complementary and important cross-check to lepton flavor violation signatures.

Lepton Flavor Violation
Lepton flavor violation is one of the most interesting probes of physics beyond the SM. The main lepton flavor violation signatures of the seesaw mechanism stem from muon decay namely, µ → eγ and µ → 3e. There are other sources of lepton flavor violation such as µ − e conversion but they are subdominant [100][101][102][103]. Other lepton flavor violating decays involving the τ lepton are less promising, unless one invokes a mechanism to significantly suppress µ → eγ [104][105][106][107][108][109][110][111][112][113]. Anyway, going back to the relevant muon decays, µ → eγ and µ → 3e, one can check that the current bounds read BR(µ → eγ) < 4.2×10 −13 , BR(µ → 3e) < 10 −12 , and future experiments aim BR(µ → eγ) < 4×10 −14 , BR(µ → 3e) < 10 −16 . Therefore, we expect an important experimental improvement in the near future. Eventually, we will superimpose these limits with the model's contribution. That said, having in mind that we have a dominant type II seesaw setup, the first contribution to µ → eγ in our models stems from 1-loop processes involving the doubly and singly charged scalars which lead to, where α EM is the fine-structure constant, G F the Fermi constant, m S ±± the mass of doubly charged scalar, m S ± the mass of singly charged scalar. In what follows will assume that the doubly charged and singly charged scalar have the same mass. This assumption will allow us to connect µ → eγ directly to µ → 3e and LHC limits on the doubly charged scalar mass.
There is an additional source of lepton flavor violation that rises from charged current which reads, where W , is a charged gauge boson defined in Eq.(2.25) and has a mass equal to 0.3v χ , v χ being the energy scale of 3-3-1 symmetry breaking. This interaction results in the following branching ratio, where, with g v , g a being the couplings constants that encompass the constants in Eq.(4.2) and the neutrino mixing matrices, and I ± ± f, 3 given by, As one can see, it is not obvious the computation of the µ → eγ decay in our model. In the regime which the right-handed neutrinos have identical masses to the W , the prediction for µ → eγ is greatly simplified yielding, Keeping right-handed neutrino masses similar to the W mass and around the weak scale one can have a visible µ → eγ decay, but if right-handed neutrino masses are much larger than the TeV scale, one can plugging the numbers in Eq.(4.3) using Eq.(4.4)-(4.6) to show that BR(µ → eγ) < 10 −15 , which is beyond reach current and projected experiments. On the other hand, even in this scenario where right-handed neutrinos are very heavy, the BR(µ → eγ) can be large due to the presence of a type II seesaw mechanism which features doubly charged and singly charged scalars contributions.
Bearing in mind that the right-handed neutrinos in our model will be heavy, and one can neglect the W contribution to the µ → eγ decay and focus on the seesaw component. That said, another important observable is the µ → 3e decay [73]. Since only the doubly charged Higgs contributes to this decay the calculation is simpler and leads to, (4.8) On one hand we can see that the Yukawa couplings f ν ab dictate the lepton flavor violation observables, but on the other hand these Yukawa couplings enter in the neutrino mass matrix. Therefore, neutrino masses and lepton flavor violation observables are correlated.
Going back to Eq.(2.22), if v 2 s13 /Λ ∼ v s11 as we will assume throughout, we have a dominant type II seesaw setup with the active neutrino masses set by vacuum The picture is not so simple because neutrinos oscillate and we need to reproduce the oscillation pattern. Therefore, these Yukawa couplings are found to be [101], where U is the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) neutrino mixing matrix of dimension 3 × 3, parametrized as follows [30],  Table 1. Table with the best-fit parameters that enter in the neutrino mass mixing according to [115].
with the mixing parameters as shown in Table 4. [114], In summary, one needs to incorporate neutrino oscillations to have more solid predictions for the lepton flavor violation observables. We will explore this fact by investigating several benchmark points which encompass normal and inverted mass ordering and different absolute neutrino masses. We will show that the neutrino mass spectrum is rather relevant to the overall lepton flavor violation signatures, a fact that has not been explored in detail in the context of 3-3-1 models. With this input from neutrino oscillations, our reasoning goes as follows: • We choose a neutrino mass ordering; • Then we pick a neutrino mass m ν1 , which then basically fixes m ν2 and m ν3 , for a given vev, v s11 = 1 − 100 eV. From this, we find the Yukawa couplings that reproduce this spectrum taking into account the oscillation patterns using Eq.(4.9); • With these Yukawa couplings we use Eq.(4.1)-(4.8) to compute the lepton flavor violating muon decays.
In Table 2 we summarize our findings using this logic. The table will allow the reader to easily follow our reasoning as we discuss the results in figure 1 where we display four of these benchmark points for m ν1 = 0.1 eV and m ν1 = 0.01 eV including normal and inverted mass ordering. The blue (red) lines are the model predictions for the muon decays for inverted (normal) neutrino mass hierarchies. The difference between the solid and dashed lines is the absolute mass for the neutrino flavor  horizontal lines are the current and projected bounds on the µ → eγ (left-panel) and µ → 3e (right-panel) decays.
One can notice that the neutrino mass ordering has a great impact on the lepton flavor violating muon decays. Looking at the first benchmark scenario with m ν1 = 0.1 eV and v s11 = 1 eV, which assumes an inverted mass ordering (IO) we get BR(µ → eγ) = 0.024/m 4 S ±± , BR(µ → 3e) = 7.5/m 4 S ±± . Having in mind that the current experimental limits, we conclude that we can probe doubly charged scalars with masses of 600 GeV and 1.3 TeV using from µ → eγ and µ → 3e decays. It is exciting to see that in this setup µ → 3e provides stronger bounds than the LHC.
However, the larger the vacuum v s11 the smaller the Yukawa couplings needed to reproduce the same neutrino masses. Hence, when we set v s11 = 100 eV, the predictions change drastically to BR(µ → eγ) = 2.4 × 10 −10 /m 4 S ±± , BR(µ → 3e) = 7.5 × 10 −10 /m 4 S ±± . For these scenarios where the vacuum of is much larger than 1 eV, LHC constitute the best probe. For concreteness, in this second case described above, taking m S ±± = 1000 GeV, would lead to muon decays much smaller than current and projected sensitivity [73], making HL-LHC and HE-LHC the best laboratories, since HE-LHC will probe masses of about ∼ 5 TeV, for instance. All these conclusions are quite visible in figure 1.
Considering the normal mass ordering (NO) we conclude that the qualitative statements do not change. Taking m ν1 = 0.1 eV, v s11 = 1 eV, quantitatively we notice that the while the inverted hierarchy gives BR(µ → eγ) = 0.024/m S ±± we find BR(µ → eγ) = 0.44/m S ±± , which is a factor of 20 larger. A larger much decay into 3e is also found for the NO compared to the IO (see figure 1). An orthogonal way to look at this is by noticing that the region between the current and projected limits delimit a signal region of lepton flavor violation. Looking at both panels of figure 1 we conclude that doubly charged scalars with masses around 1 − 2 TeV might be spotted at the µ → eγ decay, whereas the µ → 3e will be able to detect such scalars with masses of up to ∼ 30 TeV. The experimental progress on the search for the µ → 3e decay is remarkable and it will surpass even the HE-LHC regardless of the mass ordering for benchmark scenarios where v s11 ∼ 1 eV.
In the left-panel of figure 1 one can clearly see the impact of changing the value of the neutrino masses in the µ → eγ decay. This can be checked by comparing the ratio between the dashed blue and red lines with the solid blue and red lines. This is not true for the µ → 3e decay though (see right-panel of figure 1) . Comparing the NO and IO predictions for m ν1 = 0.01 eV and v s11 = 1 eV, we find no much difference, this is because the Yukawa couplings relevant for this observable are similar regardless of the neutrino mass ordering. This would continue to be true if we had taken even smaller values for m ν1 because when m ν1 is sufficiently small, the measured mass differences of the neutrino flavors govern by the neutrino mixings and thus the µ → 3e decay. When we take m ν1 = 0.1 eV, then the difference in the predictions for NO and IO is noticeable. One can easily use Table 2 and figure 1 to validate our conclusions .
We can conclude that regardless of the mass ordering and absolute value of the active neutrino masses that HE-LHC will solidly probe the type II seesaw model with respective to its contributions to the µ → eγ decay. Concerning the µ → 3e decay the situation changes due to the fantastic experimental sensitivity aimed in the near future. The µ → 3e decay will be able to probe this model for doubly charged scalar masses up to 30 TeV, which is way beyond HE-LHC reach, no matter the neutrino mass ordering. If the absolute neutrino masses are much smaller than 0.01 eV, then µ → 3e decay becomes smaller making HE-LHC still the best probe.
Our conclusions explicitly show the importance of searching for signs of lepton flavor violation in collider and muon decays. The conclusion about which probe yields stronger bounds depends strongly on the mass ordering adopted, the absolute neutrino masses and which much decay one considers. In the 1 − 5 TeV mass region of the doubly charged scalar lepton flavor violation experiments and colliders offer orthogonal and complementary probes. Thus if a signal is observed in one of the two new physics searches, the other will be able to assess whether is stems from a seesaw framework.
In summary, within the 3-3-1 model with right-handed neutrinos a type I+II seesaw naturally emerges. In the scenario where we have a dominant type II seesaw, the model offers a clear prediction for the µ → eγ and µ → 3e decays. One may wonder how one could discriminate our model from other type II seesaw proposal and a plausible answer would go as follows: setting aside the type II seesaw, our model predicts the existence of W and Z gauge bosons, as well as heavy exotic quarks. The detection of multiple signals consistent with all these particles could serve a discriminator and favor our models over others.

Conclusions
We have discussed a model which promotes SU (2) L × U (1) Y to SU (3) L × U (1) N . In this extended gauge sector, all fermions get masses via a spontaneous symmetry breaking mechanism that encompasses three scalar triplets but neutrinos. A scalar sextet is added to incorporate neutrino masses. After spontaneous symmetry breaking this scalar sextet breaks down to a doublet and scalar triplet, which play a role in the type I and type II seesaw mechanism. We focus on a scenario of type II seesaw dominance where the relevant lepton flavor violation observables namely, µ → eγ and µ → 3e, are directly tied to neutrino mass ordering and collider bounds on the doubly charged scalar.
We have explicitly shown how the absolute mass scale and neutrino mass ordering change the model predictions for lepton flavor violation within the type II seesaw framework. Combining the LHC, HL-LHC, HE-LHC sensitivity to doubly charged scalars and the experimental sensitivity to these rare muon decays we concluded that for doubly charged scalar with masses around 1 − 5 TeV, these probes are rather complementary. Moreover, regardless of the mass ordering HE-LHC is expected to solidly test any possible signal seen in these muon decays.
One may wonder how one could discriminate our model from other types II seesaw proposals and a plausible answer would rely on the existence of W and Z gauge bosons, as well as heavy exotic quarks, all predicted in our model. The detection of multiple signals consistent with all these particles could serve a discriminator and favor our models over others.

A Appendix
In this appendix will describe in more detail and pedagogical manner some sections of the model which are relevant to our reasoning.

A.1 Lepton Masses
The leptons masses could in principle be generated via the Yukawa lagrangian, where h ν ab is an antisymmetric constant coupling matrix and f ν ab is a symmetric constant coupling matrix. The first and third terms conserve the lepton flavor, while the second term violate it. Some the scalars inside the scalar sextet S do have lepton number, for this reason the third term conserves lepton number as well.
The second term is not problematic but will be removed because a set of discrete that will be invoked. This set of Z 2 symmetries are needed to avoid mixing between the SM quarks and the exotic ones. One of them requires ρ → ρ, which forbids the second term above. Although, we also need to impose e R → −e R to generate masses for charged leptons.
In particular, the charged lepton masses arise from, In this way, in the flavor space the Dirac mass matrix for the charged leptons, at tree level, is given by where v ρ is the vacuum expectation value of the neutral scalar ρ 0 and h l ab is the coupling constant matrix.
The neutrino masses come from the first second term only when the neutral components of scalar sextet acquire a vacuum expectation value as follows, We may arrange the mass terms as, With this result one can now easily understand Eq.(2.17).

A.2 Yukawa Interactions of Quarks
We mentioned in the paper that we needed to invoke some discrete symmetries to prevent mixing between the SM and exotic quarks. We will explain this statement in more detail now. The lepton number conserving terms in the renormalizable Yukawa Lagrangian for the quarks sector are while the lepton number violating terms of quarks are where h and s are constant couplings. One might notice that the terms in Eq.(A.6) will give rise to mass mixing terms involving the SM and exotic quarks, which can be problematic because they will lead to changes in the properties of the SM quarks. Therefore, one needs to prevent that and to do so invoke some discrete symmetries. The discrete symmetries have to be such that keep all the desired mass terms for the SM quarks and neutrino masses but forbid these ones. The set of discrete symmetries is, where α = 1, 2; a = 1, 2, 3. From Eq. ((A.4)), after spontaneous symmetry breaking we find, Thus, we can write the SM quark mass as follows, which leads to, Similarly, for the down-type quarks we get, with, We can see that the the type-up and type-down quarks are associated with the VEVs v η and v ρ , related to the electroweak scale.
From the Eq. (A.7), we can write the U extra quark mass term From the Eq. (A.7), we can write the D extra quark mass lagrangian The D quark mass matrix in the basis (D 1 , D 2 ) is Notice that the masses of D quarks will depend the VEV v χ , related to the TeV-scale.

A.3 Scalar Sector
Considering the three scalar fields (χ, η, ρ), the Higgs potential more general, renormalizable and invariant on the SU (3) L ⊗ U (1) X symmetry group is where µ i are constants, λ i and f are constant couplings. Furthermore, the additional terms of potential with sextet scalar, as combinations with the others scalar triplets, are (A.14) The permitted terms of the scalar potential by the discrete symmetry are In this model, the scalar triplets develop non-trivial vacuum expectation values (VEV), in order to engender correct spontaneous symmetry breaking (SSB), as written below Observe that the scalar triplet χ develops VEV only on the third neutral component, while η develops VEV only on the first neutral component. The steps of symmetry breaking transition is given by i.e., the triplet χ develops VEV breaking the 3 − 3 − 1 gauge symmetry to SM, while the triplets η and ρ develop VEV breaking the SM gauge symmetry to the QED (Quantum Electrodynamics).
Regarding the scalar sextet, the three neutral components develop VEVs in the following way: We will see that v s 1 , v s 3 and Λ are responsible for the mass for the left-handed neutrinos, while Λ is responsible for the right-handed neutrinos Dirac masses. After the SSB of In this model, the lepton number distribution of the scalars is We can see that S − 12 and S −− 22 carry two lepton numbers. Both are essential in our analysis of charged lepton flavor violating decay of muon.

A.4 Charged Lepton Flavor Interactions
The Yukawa Lagrangian with charged lepton flavor violation (CLFV) is given by where a, b = 1, 2, 3 indicates the lepton generations and m, n = 1, 2, 3 indicates the entries of the sextet, f ab is symmetric. The CLFV interactions in the µ − e sector results from taking a = 1, b = 2 in the above equation, Using the following relations: we have, We see that these terms directly induce µ → eγ lepton flavor violating decay mediated by the doubly charged scalar S ++ . Proceeding in an analogous way, we can easily obtain the relevant CLFV terms mediated by the singly charged scalar S + .

B Derivation of the µ → eγ decay
The on-shell amplitude for the l − j → l − i γ process is written in the form where ε µ is the polarization vector of the on-shell photon, e is the electric charge, u i is the spinor of l − i , u j is the spinor of l − j , m l j is the mass of l j lepton, σ µν is the comutator of γ matrices, q is the photon 4-momentum, p is the l j lepton 4-momentum, P L is the chirality projector of the left-handed lepton, P R is the chirality projection of the right-handed lepton, A L 2 and A R 2 are coupling constants. The squared modulus of the amplitude is Summing over the polarization states of the on-shell photon, we have the squared amplitude written in the form Using the completeness relation λ ε λ µ ε * λ α = −g µα , we get In order to obtain the squared amplitude explicitly, we need to calculate the term where P L = (1 − γ 5 ) /2 and P R = (1 + γ 5 ) /2 are the chirality projectors.
By using the following expression the Eq. (B.7) can be written as Reminding that P † L = P L and P † R = P R , so we get Then, we have the following result Given that P L γ 0 = γ 0 P R and P R γ 0 = γ 0 P L , we have Therefore the squared amplitude is written in the following form λ |M| 2 = −e 2 m 2 l j g µα u j A L * 2 P R + A R * 2 P L q ν σ µν u i u n σ αβ q β A L 2 P L + A R 2 P R u m , (B.13) Summing over the spins of the fermions and writing the multiplication matrix in index form, we get spins λ Using the completeness relation a u a u a = / p + m, (B. 16) and doing the average of the initial spin of fermion, we get We need to put back the equation above into normal matrix multiplication order We can write the equation above in terms of a trace Using {γ 5 , γ µ } = 0 and P 2 L = P L , P 2 R = P R , P L P R = P R P L = 0, we have (B.20) Let's calculate, ignoring m 1 , the following term Now we need to calculate this equation below Let's calculate explicitly the equation above: Inserting the g µα , Let's simplify the term above, using the contraction identities It leads us to the result Now, plugging g µα in the expression, We can written the equation above as It follows that Firstly, in order to simplify the evaluation, let's calculate the term only with / p 1 , as follows below Let's calculate the term in the braces explicitly Putting in evidence p 1α in the equation above, it leads to Taking the trace, we have 2Tr γ α γ 5 γ ν γ β }.
In order to write the squared spin-averaged amplitude in a more simplified way, we can write the product of two contractions as, Having in hands the squared amplitude, we can obtain the decay rate of the charged lepton flavor violating process for the muon decay in electron plus one photon. For this, we can need to write the decay rate for the two-bodies decay, which is given by Observe above that the rate decay of charged lepton flavor violating muon depends on the mass muon, the electric charge and it is proportional to the product of sum of the squares of A L 2 and A R 2 . It's well known that the decay rate for the muon decay (µ → eν e ν µ ), in the SM, is where G F is Fermi constant and m µ is the muon's mass. This rate decay of the muon does not violate the charged lepton flavor. The branching ratio for the (µ → eν e ν µ ) charged lepton flavor violating process is Thus, given that α = e 2 /4π (in natural units), we have finally the branching ratio of the charged lepton flavor violating muon decay for the (µ → eγ) process where α is the electromagnetic fine-structure constant, G F is the Fermi constant. Besides to µ → eγ decay, this branching ratio can be used for others decay processes, as τ → µγ.
It's important to point that the branching ratio, given by the Eq. (B.66), is very general, since we can use it for various models, for instance: supersymmetric and 3 − 3 − 1 models. In fact, we need to obtain only A L,R 2 for each model with its specific contributions related on the particle content.