Dark matter and flavor changing in the flipped 3-3-1 model

The flipped 3-3-1 model discriminates lepton families instead of the quark ones in normal sense, where the left-handed leptons are in two triplets plus one sextet while the left-handed quarks are in antitriplets, under $SU(3)_L$. We investigate a minimal setup of this model and determine novel consequences of dark matter stability, neutrino mass generation, and lepton flavor violation. Indeed, the model conserves a noncommutative $B-L$ symmetry, which prevents the unwanted vacua and interactions and provides the matter parity and dark matter candidates that along with normal matter form gauge multiplets. The neutrinos obtain suitable masses via a type I and II seesaw mechanism. The nonuniversal couplings of $Z'$ with leptons govern lepton flavor violating processes such as $\mu \rightarrow 3e$, $\mu\rightarrow e \bar{\nu}_\mu\nu_e$, $\mu$-$e$ conversion in nuclei, semileptonic $\tau\rightarrow \mu(e)$ decays, as well as the nonstandard interactions of neutrinos with matter. This $Z'$ may also set the dark matter observables and give rise to the LHC dilepton and dijet signals.

As supposed in [58], the flipped 3-3-1 model realizes a special content of fermions. In this case, one can verify that the gravitational anomaly [Gravity] 2 U (1) X vanishes. This makes the model valid up to the Planck scale, where the effect of quantum gravity becomes important [67,68]. This condition is necessary in order to ensure a manifest B −L symmetry at high energy, analogous to the standard model. Additionally, B − L neither commutes nor closes algebraically with SU (3) L , similar to the electric charge. The algebraic closure condition results in a complete gauge symmetry SU (3) C ⊗ SU (3) L ⊗ U (1) X ⊗ U (1) N , where N determines B − L, in the same situation that X defines the electric charge [44-48, 51, 52].
The B − L breaking leads to a matter parity which characterizes and stabilizes dark matter candidates. It is noteworthy that dark matter is unified with normal matter in gauge multiplets due to the noncommutativity of B − L symmetry. Additionally, the matter parity cures the unwanted vacua and interactions, which otherwise imply large flavor-changing neutral currents (FCNCs) in the lepton sector. The type I seesaw can be realized due to the B − L dynamics. However, this model contains naturally type II and III seesaws too.
The rest of this work is organized as follows. In Sec. II we examine the flipped 3-3-1 model when imposing the B − L symmetry and matter parity. The dark matter candidates are identified. The masses of lepton sextet including neutrinos are obtained. The gauge sector is diagonalized. In Sec. III we determine the tree-level FCNC coupled to Z . The setup also implies the tree-level FCNC coupled to the standard model Higgs boson, but it is subleading and neglected. In Sec. IV, the lepton flavor violating processes, the nonstandard interactions of neutrinos with matter, the LHC dilepton and dijet searches, and the dark matter observables are obtained. Finally, we conclude this work in Sec. V.

A. Proposal
As stated, the 3-3-1 gauge symmetry is given by where the first factor is the usual QCD group, while the last two are a nontrivial extension of the electroweak group. The electric charge and hypercharge are embedded as where T n (n = 1, 2, 3, ..., 8) and X are SU (3) L and U (1) X generators, respectively. The coefficient β determines the electric charge of new particles, which is arbitrary on the theoretical ground and independent of all the anomalies.
The new observation [58] is that the [SU and quarks have three colors. In general, it is proved that the generation number must be a multiple of three in order to embed left-handed fermion doublets in SU (3) L representations while keeping right-handed fermions as singlets similar to the standard model, which provides a partial solution to the flavor question.
That said, the left-handed fermion representations under SU (3) L are generally given by plus right-handed fermion singlets, ξ aR , e aR , k aR , u aR , d aR , and j aR . The generation indices are α = 2, 3 and a = 1, 2, 3. The new fields ξ 1,2,3 , k a , and j a possess electric charges as superscripted, which depend on a basic electric charge parameter q = −(1 + √ 3β)/2. It is easily checked that all the other anomalies vanish. An alternative version can be proposed, such that a lepton generation is in antisextet, two other lepton generations in antitriplets, and all the quark generations in triplets, under SU (3) L . However, this modification should be equivalently physical to the original proposal, which will be skipped.
In this case, ξ 0 2R is a 3-3-1 singlet which can be omitted as ν aR are. Additionally, ξ − 3R can be discarded while ξ − 3L is replaced by (ξ + 1R ) c = ξ − 1L . This restriction results in the most 1 The 3-3-1 model with β = 1/ √ 3 has a Landau pole larger than the Planck scale, as expected [69,70]. economical model, explicitly written under the 3-3-1 symmetry as where we relabel k as E and j as U , since they have the same electric charge as of e and u, respectively. The subscripts of ξ are omitted without confusion. The authors of [58] pointed out that since the model cooperates a real triplet ξ under SU (2) L , it can be made heavy in order to keep the model phenomenologically viable. Further, we will show that the simplest version provides dark matter candidates as well as neutrino masses naturally. Otherwise, for q = −1, the fermions ξ must have right-handed components and gain light Dirac masses in the weak scale, which are strongly disfavored by the electroweak data.
The scalar content responsible for symmetry breaking and mass generation is given by Note that ρ and χ are identical under the gauge symmetry, but distinct under the B − L charge, as shown below.
where N determines B − L in the same situation that X does so for Q [44-48, 51, 52].
This way leads to the group structure SU Additionally, N and B − L are gauged charges since they are related to the gauged charge T 8 . The 3-3-1-1 gauge theory requires ν aR in addition to the existing fermions in order to cancel the B − L anomalies as well as a scalar field φ that couples to ν R ν R and necessarily breaks U (1) N . The N -charges of multiplets are summarized in Tab. I Assuming that U (1) N is broken at high energy due to φ, the heavy particles such as the U (1) N gauge (C) and Higgs (φ) bosons as well as right-handed neutrinos (ν aR ) are all integrated out. The imprint at low energy is only the 3-3-1 model, conserving the matter parity as residual gauge symmetry, which is defined by the vacuum of the mentioned U (1) N breaking field [46]. Because the matter parity is conserved, the lightest W -particle (LWP) is stabilized, responsible for dark matter. The dark matter candidates include a fermion ξ 0 , a vector Y 0 , and a combination of ρ 0 3 , χ 0 2 , and S 0 23 . Due to the gauge interaction, Y 0 annihilates completely into the standard model particles. The realistic candidates that have correct abundance are only the fermion or scalar, as shown below.

C. Lagrangian
Hereafter, we consider the theory at low energy that includes only the light fields. The contribution of heavy fields (C, φ, ν aR ) is separately mentioned if necessary.
The total Lagrangian consists of where the first part composes the kinetic terms plus gauge interactions.
The second part includes Yukawa interactions, obtained by Note that the unwanted Yukawa interactions are L / Yukawa = s e αaψ αL χe aR + s e 1aψ 1L Se aR + s E αaψ αL ρE aR +s u abQ aL χ * u bR + s U abQ aL ρ * U bR + H.c., which are suppressed due to the U (1) N symmetry or matter parity at low energy. In the ordinary 3-3-1 models, they are present, characterizing an approximate B − L symmetry. In other words, they should be small in comparison to the normal couplings s h, respectively.
The last part is the scalar potential, Here the couplings λ's are dimensionless, while the parameters µ's have mass dimension.
The unwanted interactions include which are suppressed by the U (1) N symmetry or matter parity at low energy, where f 's are mass parameters analogous to µ's.

D. Neutrino mass
First note that due to the matter parity conservation, only the even scalars develop vacuum expectation values (VEVs), such that Substituting the VEVs into the Yukawa Lagrangian, the quarks and exotic leptons gain suitable masses as follows The ordinary leptons obtain masses Note that the constraint from the ρ parameter implies κ < ∼ O(1) GeV. Hence, a small mixing between ξ ± and E ± a proportional to κ can be neglected, since κ w, Λ, where the last two VEVs are in the TeV scale. The fields µ and τ get desirable masses. However, the electron and last two neutrinos have vanishing masses, which are inconsistent.
The electron mass vanishes similarly to the original study, which can be radiatively induced [58]. An extra remark is that the interaction ψ 1L ψ 1L S ⊃ (ξ + ξ 0 ξ − )(ν 1L e 1L )(S + 13 S 0 23 ) provides the need for a type III seesaw, where the mediator is a heavy fermion triplet, ξ.
But the relevant neutrino (ν 1 ) mass generated is zero due to S 0 23 = 0 resulting from the matter parity conservation, which does not change the above value of a type II seesaw.
As mentioned, the heavy fields φ, ν R are present and can imply neutrino masses via We achieve Dirac masses Fitting the data m ν ∼ 0.1 eV, we obtain φ ∼ [(h ν ) 2 /h R ]10 14 GeV, since u is proportional to the weak scale. Given that h ν , h R ∼ 1, one has φ ∼ 10 14 GeV, close to the grand unification scale. The right-handed neutrinos ν aR have masses in this scale.
It is clear that two neutrinos ν 2,3L achieve masses via the type I seesaw with the corresponding mixing angle θ 23 comparable to the data, while the neutrino ν 1L has a mass (which one sets h ξ κ ∼ 0.1 eV) via the type II seesaw and does not mix with ν 2,3L . The mixing angles θ 12 and θ 13 can be induced by an effective interaction, such as where M is the new physics scale which can be fixed at M = φ . The mass matrix of observed neutrinos is corrected by where v is proportional to the weak scale like u. These elements can generate appropriate mixing angles for θ 12,13 . In this case, the neutrino mass matrix is generic and small.

E. Gauge sector
The mass Lagrangian of gauge bosons is given by where (t n , T n , X), and (G n , A n , B) correspond to the coupling constants, generators, and gauge bosons of the 3-3-1 groups, respectively 2 . It acts as for a sextet, where T n = 1 2 λ n are Gell-Mann matrices and t X = g X /g.

Define the non-Hermitian gauge bosons,
They are mass eigenstates by themselves with corresponding masses, Because of the matter parity conservation, there is no mixing between W and X as well as A 6 and neutral gauge bosons. Because the type II seesaw requires an infinitesimal κ, its correction to m W and m Y has been suppressed.
The neutral gauge bosons (A 3 , A 8 , B) mix via a 3×3 mass matrix. This yields the massless photon field, where s W = e/g = √ 3t X / 3 + 4t 2 X is the sine of the Weinberg's angle [71]. As usual, we define the remaining neutral fields orthogonal to A as In the new basis (A, Z, Z ), there is only Z-Z mixing via a 2 × 2 mass matrix, while A is decoupled. This mass matrix yields eigenstates with masses and the mixing angle Since κ is tiny, its contribution to the ρ parameter is neglected. The deviation of the ρ parameter from the standard model prediction is due to the Z-Z mixing, obtained by From the W mass, we derive u 2 + v 2 = (246 GeV) 2 . From the global fit, the PDG Collaboration extracts the ρ deviation as ∆ρ = 0.00039 ± 0.00019, which is 2σ above the standard model prediction [72]. We contour ∆ρ as a function of u = 0-246 GeV, since v is related to u. Generally for the whole u range, the new physics scales are bounded by where both ∆ρ and ϕ vanish, i.e. the new physics is always decoupled when w, Λ tend to zero. Therefore, we can close the 3-3-1 symmetry at the weak scale in this regime [74,75].

III. FCNC
Because the first lepton generation transforms differently from the last two lepton generations under SU (3) L , there exist FCNCs at the tree-level associated with leptons. Note that the FCNCs conserve quark flavors, in contrast to the ordinary 3-3-1 models.
Indeed, the neutral currents of fermions depend on the Cartan (diagonal) charges T 3,8 and where F runs over fermion multiplets. It is easily verified that all the quarks u aL , u aR , d aL , d aR , U aL , and U aR as well as the right-handed leptons e aR and E aR correspondingly do not flavor-change. Also, the terms of T 3 and Q conserve all fermion flavors. There remains where note that T 8 ψ αL = 1 2 λ 8 ψ αL , while T 8 ψ 1L = 1 2 λ 8 ψ 1L + ψ 1L 1 2 λ 8 and the corresponding interaction is traced. We obtain where T ν = T e = diag(−1, 1, 1) and T E = diag(−4, −2, −2). We commonly denote l = ν or e or E and change to the mass basis l L,R = V lL,R l L,R . The relevant Lagrangian is where which takes the same form for both l = ν, e and l = E as well as flavor-changes for i = j.
Let us recall that all quark generations transform universally under the 3-3-1 group, which yields the flavor-conserved current with Z as where we denote physical eigenstates for up-quarks as u = (u c t) T , down-quarks as d = (d s b) T , and commonly q = u or d in further investigation.
Integrating out the heavy field, Z , from (42) and/or (44), we obtain an effective Lagrangian that sums over six-dimensional interactions (operators) relevant to the standard model fermions at the tree-level, such that For convenience, we have relabeled the lepton eigenstates (l i ) to be l α , which are determined by a Greece (letter) index α = 1, 2, 3 and should not be confused with a generation index, such that e α = e, µ, τ and ν α = ν 1 , ν 2 , ν 3 . A coefficient η q applies for the up-quarks as η u = 2 and the down-quarks as η d = −1.
The present constraints for the effective Lagrangian come from several processes. The first two terms (45) and (46) provide charged lepton flavor violating processes like µ → 3e, τ → 3e, τ → 3µ, τ → 2eµ, and τ → 2µe. The next four terms (47), (48), (49), and (50) present wrong muon and tau decays as well as the nonstandard neutrino interactions that concern both constraints from oscillation and non-oscillation experiments. The last four terms (51), (52), (53), and (54) describe semileptonic τ → µ(e) decays and µ − e conversion in nuclei as well as the signals for new physics (dilepton, dijet, etc.) at low energy such as the Tevatron. Therefore, in the following, we consider the phenomenological aspect due to the presence of the interactions given above.

A. Leptonic three-body decays
To perform the analysis in the current section, we suppose that the major sources con- Here, we use the notations c ij ≡ cos θ ij and s ij ≡ sin θ ij , and the angles vary in the corresponding ranges, θ ij = [0, π/2] and δ = [0, 2π].
1. τ + → µ + µ + µ − , τ + → e + e + e − and µ + → e + e + e − The lepton flavor violating processes of a lepton decaying into three identical lighter ones, which are usually called as the type I trilepton decays, have been considered before in various models and scenarios of elementary particle physics [76][77][78][79][80][81]. On the experimental side, the temporary upper bound constraints on the branching ratios of those processes have also been obtained, such as [72] Br(µ → 3e) < 1.0 × 10 −12 , Thus, the upper bound on Br(µ → 3e) is about four orders more stringent than the corresponding channels of τ decays.
The lower bound on the new gauge boson mass m Z comes from the LHC dilepton and dijet searches, which is roughly set to be larger than 4 TeV for the Z couplings analogous to the standard model Z couplings [82,83]. However, in this model the Z mass takes a smaller bound, say 2.8 TeV, as shown below. Because the transferred momentum, whose maximal value is about the τ mass, is much smaller than m Z , the type I trilepton decay branching ratios, for instance τ + → µ + µ + µ − , can be calculated with a high enough precision using the effective Lagrangian written as follows where The branching ratio for this process from eq. (59) can be found in [76], In the same way, the effective Lagrangians and branching ratios for τ + → e + e + e − and µ + → e + e + e − channels can be derived. That said, their expressions could be obtained from (59) and (61) Here the minimal (maximal) values frequently occur at many relevant values of (θ ij , δ ), for example, at (sin θ 12 , sin θ 13 , sin θ 23 , δ ): 0.7074, 0.0000, 0.4336, 1.04π (0.1005, 1.0000, 0.8601, 1.70π), 0.0000, 0.7071, 0.2397, 0.45π (0.8498, 0.0000, 0.0464, 0.99π), 1.0000, 0.7072, 0.5033, 0.39π (0.0000, 0.0000, 0.8082, 0.38π), according to Br(µ → 3e), Br(τ → 3e), and Br(τ → 3µ), respectively. Indeed, the minimum is easily seen to be at the points that Γ lZ µe , Γ lZ τ e , and Γ lZ τ µ vanish, respectively. In Fig. 1, we illustrate the dependence of the branching ratio Br( → 3 ) as a single variable function of the new gauge boson mass m Z for θ 12 = π/3, θ 13 = π/6, θ 23 = π/4, 2. τ + → µ + e + e − and τ + → e + µ + µ − Similar to the decays → 3 considered in the previous section, when the new gauge boson mass m Z is at order TeV or larger, the type II three leptonic decay, e.g. τ + → µ + e + e − , can be well described by an effective Lagrangian, which takes the form where Then the branching ratio for τ + → µ + e + e − decay could be expressed as [76], With the aid of the interchange symmetry µ ↔ e, the expression for the branching ratio of the decay τ + → e + µ + µ − could be easily obtained from the above formula by appropriately replacing (Γ lZ τ µ , Γ lZ ee ) by (Γ lZ τ e , Γ lZ µµ ). Taking into account the same set of values of θ ij and δ as in the previous section, we depict in Fig. 2 the behaviors of the branching ratios, τ → eµµ and τ → µee, according to the variation of the new neutral gauge boson mass m Z (green and red lines), respectively. In this figure, the current upper bounds Br(τ → eµµ) ≤ 2.7 × 10 −8 and Br(τ → µee) ≤ 1.8 × 10 −8 are also shown as blue lines [72]. The lower limits for m Z obtained from Fig. 2 responsible for these two processes are roughly 3 TeV. This agrees with the limit from µ → 3e in the case of Fig. 1 right panel, but is about 20 times less stringent than the limit from µ → 3e in the case of Fig. 1 left panel. Particularly, using the current lower bound from µ → 3e in the latter case, i.e. m Z ≥ 65.3 GeV at θ 12 = π/3, θ 13 = π/6, θ 23 = π/4, and δ = 0, the precision calculation shows that One the other hand, with the same strategy as in the previous case when varying θ ij in [0, π/2] and δ in [0, 2π], the bounds derived are 3. τ + → µ + µ + e − and τ + → e + e + µ − In this part of search, we consider the type III leptonic decay modes in which lepton flavors are violated by three units, such as τ + → µ + µ + e − . In the flipped 3-3-1 model, the process is dominantly contributed by an unique effective interaction, which is expressed as where Then, the branching ratio for τ + → µ + µ + e − could be easily written down by generalizing the results given in [76] as On the other hand, the effective Lagrangian and branching ratio for Br(τ + → e + e + µ − ) can be achieved from (75) and (77) simply by interchanging the fields µ ↔ e. Currently, the experimental constraints on the branching ratios of the decays τ + → µ + µ + e − and τ + → e + e + µ − are, indeed, very weak. Therefore, the corresponding lower limits on the new neutral gauge boson mass m Z given by these two channels are less substantial than those obtained from the channels considered above, especially for the µ → 3e decay.
To be concrete, the behaviors of the τ + → µ + µ + e − and τ + → e + e + µ − branching ratios can be found in Fig. 3. Using the same trick as in the previous case, we obtain the theoretical upper bounds as follows for θ 12 = π/3, θ 13 = π/6, θ 23 = π/4, and δ = 0. Hence, these bounds are about six orders below the sensitivities of the current experiments.
Last, but not least, varying the mixing angles and phase in the allowed regimes, the specified ranges for the τ + → µ + µ + e − and τ + → e + e + µ − branching ratios are found as

Comment on wrong µ and τ decays
It is not hard to point out that the wrong muon and tau decays, e.g. µ → eν eνµ and τ → µν µντ , take the same rate as of those in the previous section, respectively. Hence, such decays are far below the experimental limits Br ∼ 0.1 [72].

B. Semileptonic τ → µ and τ → e decays
The next topic we discuss in this paper is the semileptonic decays of τ , say Br(τ + → + P ) and Br(τ + → + V ), (82) in which = e, µ and P, V stand for neutral pseudoscalar and vector mesons, respectively.
These decay channels have been studied formerly in the other models as well as modelindependent scenarios [76,84,85]. In the present model, these processes happen dominantly through the exchange of the new neutral gauge boson with lepton-flavor violating interactions. On the experimental side, the following upper bounds have been obtained [72]: The near future and planned experiments such as LCHb [86], BESIII [87], Belle II [88], and COMET [89] will improve upon these present limits.
As already supplied in (51) and (52), the effective Lagrangian for the semileptonic decays of the τ → µ(e) types can be rewritten in the form where = µ or e and the couplings are The branching ratio for the τ + → + P decay, where P denotes a neutral pseudoscalar meson as π 0 or η 0 or η 0 , can be derived from (84) as follows [76] Br where τ τ is the tau's lifetime and for P = η 0 or P = η 0 , and for P = π 0 . Here the decay constant f π = 135 MeV, while f q η,η and f s η,η are defined as with f q 1.07f π , f s 1.34f π and φ η 0.2183π.  The theoretical maximal values given in (90) are roundly the same orders as the current experimental limits [cf. (83)]. Thus the constraints on the new neutral gauge boson mass derived from these channels have much appealing, validating the model. This can be more visually seen in Fig. 4, where we show the dependence of the branching ratios Br(τ + → + P ) on the new neutral gauge boson mass m Z ≡ M for = e, µ and P = π, η, η . Since the current upper bounds on the branching ratios are around 10 −7 , the lower limit obtained for the new neutral gauge boson mass m Z is about 3 TeV, which is the same limit set by the searches of the LHC dilepton and dijet signals.
Similar conclusions are also obtained for the case of τ + → + V decay, where V is taken as a neutral vector meson ρ or ω or φ. The branching ratios for these processes can be derived from the Lagrangian (84) as [76] Br where G LV is an effective coupling that takes one of the following forms, The detail behaviors of the branching ratios Br(τ + → + V ) are depicted in Fig. 5 for θ 12 = π/3, θ 13 = π/6, θ 23 = π/4, and δ = 0. Comparing to the experimental bounds yields a Z mass around 3 TeV.

C. µ − e conversion in nuclei
In this sector, we consider a hypothetical process, called µ → e conversion in nuclei, in which negative muons are captured in a target of atomic nuclei, such as Titanium (Ti), Aluminum (Al) or Gold (Au), to form muonic atoms. The muon then converses into an electron in the nuclear field without creating a neutrino. There are a number of experiments that have been built or planned to built to search for the process's signals, for instance TRIUMP [90], SINDRUM-II [91], and COMET [92]. The current experimental limits on the branching ratios are 4.3 × 10 −12 at TRIUMF for Titanium target and 7 × 10 −13 for Gold In the considering model, the µ − e conversing ratio can be calculated from a Lagrangian, which has same form as (84) in which The conversing ratio (CR), which is obtained after normalizing to the total nuclear capture rate ω capt , can be simply expressed as where V (n) and V (p) are the overlap integrals and ω capt is the total capture rate. For the cases of 48 22 Ti, 27 13 Al, 197 79 Au, and 208 82 Pb nuclei, they are given in Table III [93].
And, f is the standard model first generation fermion (e, u, d). The NSIs can affect neutrino oscillations in multiple ways, in the production, dectection and propagation of neutrinos.
In this work, we consider the neutrino propagation in matter with the NSIs, assuming no effect of production and detection with the NSIs. The NSIs effect the neutrino propagation via coherent forward scattering in Earth matter. The Hamiltonian, which governs the propagation of neutrino flavor states in matter including the NSIs, is written as followŝ where E is neutrino energy, U is lepton mixing matrix, and m i are neutrino masses. Furthermore, A = 2 √ 2EG F N e is the effective matter potential that is driven by ordinary charged-current weak interaction with electron. N e is the electron density along the neutrino trajectory. The matrix m has a form as where αβ is called a matter NSI parameter and is defined as where N f is the number density of a fermion of type f . In the considering model, the form of f C αβ can be obtained as follows The effective Hamiltonian in (102) governs the neutrino propagation in matter due to the NSIs, hence varying the neutrino oscillation probabilities in comparison to the normal case [98]. However, at present, there is no evidence for the NSI associated with the experimental data of neutrino oscillations. The latest constraints on the NSIs from the global analysis of oscillation data can be found in [99], which give the most stringent constraints on u αβ and which is at order 10 −2 , 10 −4 , and 10 −6 for m Z = 1, 10, and 100 TeV, respectively. The predicted values lie within the experimental limits but do not give significant constraints on the model parameters. Comparing to the u αβ and d αβ bounds, it is hard to probe the nonstandard neutrino interactions in the model for m Z > 3 TeV.

E. LHC dilepton and dijet searches
Since the new neutral gauge boson Z directly interacts with ordinary quarks (q) and leptons (l), the new physics processes pp → Z → ff for f = q, l exist at the LHC dominantly contributed by Z .
The cross-section for producing a dilepton or diquark final state can be computed with the aid of the narrow width approximation [100], where the parton luminosities dL qq /dm 2 Z at the LHC √ s = 13 TeV can be found in [101].
The partonic cross-section and branching ratio Br where N C is the color number of f and assuming that Z decays only to fermions. Indeed, it is easily verified that the other Z decay channels such as to ordinary Higgs and gauge bosons give small contribution to the total Z width.
Because the current bound on dijet signals is less sensitive than the dilepton [83], the corresponding Z mass limit is quite smaller than that obtained from the dilepton, which is not included. In sort, in the present model, the dijet sinals predicted are negligible, given that the dilepton bound applies for Z .

F. Dark matter
The model contains two kinds of dark matter candidates: (i) the fermion triplet ξ which is unified with the standard model lepton doublet (ν 1L e 1L ) in the SU (3) L sextet and (ii) the scalar that is either ρ 3 or a combination (called D) of χ 2 and S 23 , whereas the remaining combination of χ 2 and S 23 is the Goldstone boson of the Y gauge boson. The candidate D transforms as a standard model doublet, which interacts with Z. This gives rise to a large direct dark matter detection cross-section that is already ruled out [102]. The singlet candidate ρ 3 can fit the relic density and detection experiments, which has been studied extensively [45,103]. The fermion candidate ξ is a new observation of this work 3 .
Note that at the tree-level, the components of ξ triplet have degenerate masses, already obtained as m ξ . However, the loop effects of gauge bosons can make the ξ ± mass bigger than the ξ 0 mass by m ξ ± − m ξ 0 = 166 MeV [104]. Therefore, ξ 0 can be regarded as a LWP responsible for dark matter.
In the early universe, the dark matter candidate ξ 0 can (co)annihilate into the standard model particles that set its abundance. Generalizing the result in [104], we obtain the annihilation cross-section, where (α/150 GeV) 2 1 pb. Comparing to the observation, we have Ω ξ h 2 0.1 pb/ σv 0.11 [72], which implies m ξ 2.86 TeV.
The dark matter ξ 0 can scatter off nuclei causing observed effects in the direct detection experiments. At the tree-level, it does not interact with quarks confined in nucleons, since T 3 (ξ 0 ) = Y (ξ 0 ) = 0. The direct detection cross-section starts from the one-loop level, contributed by W, h and ξ ± , leading to σ SI 1.2 × 10 −45 cm 2 [105], in agreement with the experiment for the heavy dark matter mass m ξ 2.86 TeV [106].
The above evaluation is valid when the new gauge (Z ) and Higgs S 33 portals as well as the 3-3-1 model new particles are heavier than the dark matter, so that the dark matter observables are governed by the standard model particles. Alternatively, since both Z and S 33 couple to ξ 0 , this dark matter can annihilate to the standard model Higgs, weak bosons and top quark as well as the appropriate new particles of the 3-3-1 model. In this case, the Z and S 33 resonances set the dark matter density and direct detection cross-section [43,45,103]. That said, we have two viable regimes responsible for the dark matter mass m ξ 0 = 1 2 m Z and m ξ 0 = 1 2 m S 33 , provided that m Z and m S 33 are separated. By contrast, these regimes are coincided. The Z mass bound tells us that the dark matter mass is at TeV or higher scale. This mass is easily evaded the direct detection [106].

V. CONCLUSION
The discovery of the flipped 3-3-1 model has changed the current research of the 3-3-1 model. Indeed, the flavor nonuniversality that is now associated with leptons due to the anomaly cancellation provides the need for realizing the type I and II seasaw mechanism and the matter parity naturally. We have shown that the neutrino masses are produced, yielding the seesaw scales κ at eV and M at 10 14 GeV. Whereas, the dark matter candidate may be a fermion or scalar that has the mass in the TeV scale.
Because of the lepton generation nonuniversality, the charged lepton flavor violating processes and nonstandard neutrino interactions arise at the tree-level due to the exchange of Z boson. We have made a systematic search for such processes and found that the 3-3-1 breaking scale is coincided with those obtained from the LHC dilepton and dark matter constraints, where the Z portal governs the dark matter observables. This is an advantage over the ordinary 3-3-1 models since their FCNC constraints often do not coincide with the collision bounds [45,53,66].
All the results arise from the gauge symmetry principles. Hence, the flipped 3-3-1 model is very predictive that deserves further studies.