Probing neutrino and Higgs sectors in $SU(2)_1 \times SU(2)_2 \times U(1)_Y$ model with lepton-flavor non-universality

The neutrino and Higgs sectors in the $\mbox{SU(2)}_1 \times \mbox{SU(2)}_2 \times \mbox{U(1)}_Y $ model with lepton-flavor non-universality are discussed. We show that active neutrinos can get Majorana masses from radiative corrections, after adding only new singly charged Higgs bosons. The mechanism for generation of neutrino masses is the same as in the Zee models. This also gives a hint to solving the dark matter problem based on similar ways discussed recently in many radiative neutrino mass models with dark matter. Except the active neutrinos, the appearance of singly charged Higgs bosons and dark matter does not affect significantly the physical spectrum of all particles in the original model. We indicate this point by investigating the Higgs sector in both cases before and after singly charged scalars are added into it. Many interesting properties of physical Higgs bosons, which were not shown previously, are explored. In particular, the mass matrices of charged and CP-odd Higgs fields are proportional to the coefficient of triple Higgs coupling $\mu$. The mass eigenstates and eigenvalues in the CP-even Higgs sector are also presented. All couplings of the SM-like Higgs boson to normal fermions and gauge bosons are different from the SM predictions by a factor $c_h$, which must satisfy the recent global fit of experimental data, namely $0.995<|c_h|<1$. We have analyzed a more general diagonalization of gauge boson mass matrices, then we show that the ratio of the tangents of the $W-W'$ and $Z-Z'$ mixing angles is exactly the cosine of the Weinberg angle, implying that number of parameters is reduced by 1. Signals of new physics from decays of new heavy fermions and Higgs bosons at LHC and constraints of their masses are also discussed.


I. INTRODUCTION
One of the most important purposes of the LHC is to search for manifestations of new physics (NP). It seems that some clues have appeared with massive neutrinos and recent observations of lepton-flavor non-universality (LNU). Recall that the lepton family replication is assumed in the Standard Model (SM). Therefore, the lepton-flavor is universal in the latter. For the recent two decades, neutrino and Higgs physics are hot topics in Particle Physics. With increasing luminosity and beam energy, the LHC becomes a powerful tool for searching for NP. With larger masses, the third generation seems to be more interesting, in the sense of the sensitivity to NP. Nowadays, there are two kinds of anomalies in the semileptonic B meson decays which are captivating for the LNU. The first one is the class of the following ratios of branching fractions: which show 3.5 σ deviations from the corresponding SM predictions [1], R D * = 0.252 ± 0.004 , R D = 0.305 ± 0.012.
The above results provide hints for violation of the lepton flavor universality (LFU).
From the physical point of view, the mass of a particle plays a quite important role in its characteristic properties. To justify this, let us mention some well-known examples. The first is that the proton and neutron have a tiny mass difference (940 vs 938) MeV, but the proton is long-lived while the neutron is unstable with its mean lifetime of just under 15 min (881.5 ± 1.5 s). The second example is the situation with the electron and the muon.
Both particles are leptons with just a mass difference (0.511 vs 105.6) MeV, the electron is 1 The SM value for R K has been first obtained in Ref. [3] stable while the muon is unstable with the mean lifetime of 2.2 µs. So one may expect that the third generation of quarks and leptons where particles are heavier, has to be different from the first two ones. Within this context, the above data showing the LNU look quite understandable. In other words it is quite natural to expect that the third fermion generation is more strongly coupled to some New Physics than the first two ones. Recently the R D and R D * were subjects of intensive studies mostly in scalar leptoquark models [4,5].
One of the beyond the SM models satisfying the recent experimental data of LNU is the model based on the SU(2) 1 × SU(2) 2 × U(1) Y (G221) gauge group [6] (more kinds of G221 models can be found in Ref. [7]). In Ref. [6] the authors have mainly concentrated on explanation of LNU in the lepton sector. But at present, any theoretical model in particle physics has to deal with neutrino masses, the baryon asymmetry of the universe (BAU) and the dark matter (DM).
The aim of this work is to study further details in the gauge, Higgs and neutrino sectors of the model presented in Ref. [6]. We will show that the problems of the active neutrino mass and DM in this model can be solved without any changes of results of allowed parameter regions satisfying all constraints of the flavor physics, tau decays, electroweak precision data, and recent anomalies in B decays, which were indicated in Ref. [6]. In particular, the active neutrinos get Majorana masses from radiative corrections, where new lepton-number violating interactions have to be introduced. The simplest way is the Zee method [8], where a pair of singly charged scalars transformed as singlets under both the SU(2) gauge groups is introduced. Like in the Zee models, where a second SU(2) L Higgs doublet is necessary for creating a nonzero triple coupling of two Higgs doublets and a singly charged Higgs singlet, the SU(2) 1 Higgs doublet φ ′ in this G221 model plays the role of the second SU(2) L Higgs doublet. Hence, no new breaking scales need to appear, implying that there are no new mass terms contributing to the fermion and gauge boson sectors. This explains why all results investigated in Ref. [6] are unchanged, therefore we can use them to study the coupling properties of the Higgs and gauge bosons with fermions. In addition, it suggests that the ways of generating active neutrino masses in many recent radiative neutrino mass models can be applied to the G221 model. Many of these models have DM candidates that are neutral fermion singlets and have odd charges under a new Z 2 symmetry. To avoid complicate Higgs sectors, where just new charged Higgs bosons are included, we will not pay much attention to models solving the DM problems in this work. We will discuss in detail the mechanism of generating neutrino masses from the Zee mechanism, and the Higgs potential with the appearance of two singly charged Higgs singlets. In the gauge boson sector, we will apply general method to diagonalize neutral and charged gauge boson sectors, and from this we get a consequence that the tangents of the mixing angles in two sectors are proportional.
This will reduce the number of the model parameters by 1. In the Higgs sector, the physical Higgs spectrum is presented. Then the SM-like Higgs boson and its couplings to other SMlike particles are identified and compared with the SM predictions. A comparison between properties of the Higgs spectrum in the G221 model and the minimal supersymmetric model (MSSM) and two Higgs doublet models (THDM) will also be discussed in this work. Based on these properties and the constraints of parameters given in [6], we will discuss the bounds of new Higgs boson masses as well as promoting decay channels of the Higgs bosons and fermions that can be searched for at modern colliders such as the LHC.
This paper is organized as follows. In section II, after a brief review of the model, we present a more careful consideration of charged lepton masses and the Zee method for generation of neutrino masses. In the subsection II B we suggest two possibilities of appearance of DM candidates in the G221 model. The first is based on a radiative neutrino mass model introduced previously. This way will not change the results of parameter constraints in Ref. [6]. The second way is different, because a new scalar SU (2)  charged Higgs bosons are also discussed. In section V, we review briefly the allowed regions of parameters given in [6], which resulted from a specific numerical illustration in the limit of two vector-like fermion generations and simple textures of Yukawa couplings. Following the searches for new heavy particles at the LHC, we use these allowed regions to investigate lower bounds of masses and promoting decay channels of new fermions and Higgs bosons predicted by this model. Conclusions are given in the last section VI.

II. BRIEF REVIEW OF THE MODEL
The model is based on the gauge group SU(2) 1 ×SU(2) 2 ×U(1) Y with the following gauge couplings, fields and generators [6]: where i = 1, 2, 3 is the SU(2) index. All the chiral fermions transform as where the numbers in brackets refer to SU(3) C , SU(2) 1 , SU(2) 2 , and the hypercharge. The electric charge operator is determined in the form For the subgroup SU(2) 1 there are n V L generations of vector-like fermions which are transformed as its doublets, while they are singlets for the SU(2) 2 , The vector-like fermion generation number is greater than one in order to explain successfully the LNU, and it was fixed by n V L = 2 for simplicity in numerical illustration [6].
The Higgs sector consists of two doublets φ and φ ′ and one self-dual bidoublet Φ (i.e., with components as withΦ 0 = (Φ 0 ) * . The scalar fields develop VEVs The spontaneous symmetry breaking (SSB) of the model follows the pattern The main phenomenology of the model concerned B-decay anomalies and the lepton-flavor non-universality has been presented in [6]. However, the current physical model has to satisfy Higgs and neutrino physics as well as DM candidate.
With the above breaking chain, the VEVs are assumed to satisfy the relation Yukawa Lagrangian, fermion mass matrices, and diagonalization steps to construct physical states and masses of fermions were presented in detail in [6]. Hence, we will summarize here only important results and focus on new features of generating active neutrino masses from loop corrections.

A. Charged fermion masses
The chiral fermions couple to the SM Higgs-like φ doublet whereφ ≡ iσ 2 φ * . The matrices y d , y u , y ℓ are 3 × 3 matrices. The vector-like fermions can have gauge-invariant Dirac mass terms Other contributions are where λ † q, ℓ andỹ u, d, ℓ are n V L × 3 Yukawa matrices. After the SSB, the above couplings will induce mixing between the vector-like and the SM chiral fermions. This is crucial for the phenomenology of the model.
For the sake of simplicity one can assume a softly broken discrete Z 2 symmetry under which only φ ′ is odd, making unnecessary Yukawa couplings vanish, i.e.,ỹ u, d, ℓ ≃ 0 [6]. There is another charge assignment that also forbids Lagrangian in (16), while keeps φ ′ even: only Q L and L L are odd. This is necessary for generating active neutrino masses by the Zee method considered in this work.
We combine the chiral and vector-like fermions as where i = 1, 2, 3, k = 1, · · · , n V L and I = 1, · · · , 3 + n V L . After the SSB, the fermion mass Lagrangian has the form Here, all above mass matrices are (3 + n V L ) × (3 + n V L ) and have the form In the limit ǫ = v/u ≪ 1, these matrices are blocked-diagonalized perturbatively via two steps. After that, the SM parts are separated from the total. The transformations of fermion states are: , and W f are (3 + n V L ) × (3 + n V L ) unitary matrices [6]. At the first step where v = 0, every M F (F = U, D, E) is diagonalized by an exact V F depending on u, M F and λ ℓ,q . At the second step, transformations V f and W f are expanded in terms of power series of ǫ, V f = 1 + iǫ 2 H f V + ... and W f = 1 + iǫH f W + 1/2(iǫH f W ) 2 + .... They were listed precisely in [6]. After the two steps, all original mass matrices in (19) will be transformed into block-diagonal formsM F = V f V F M F W † f . One of the blocks in everyM F is identified with a SM fermion block, which is diagonalized by 3 × 3 Only the CKM matrix, V CKM = S u S † d , appears in the gauge couplings [6]. We can fix S e = U e = I 3 .
For studying Higgs boson phenomenology satisfying the allowed regions of parameters given in [6], which resulted from a specific assumption of two new lepton families and textures of Yukawa couplings λ q,ℓ , we will present more detailed masses and eigenstates of charged leptons. The quark sector can be derived similarly. In the flavor basis E of charged leptons, the mass matrix M E in (19) is 5 × 5. Following Ref. [6], a simple texture of λ ℓ is chosen as where new parameters ∆ µ and ∆ τ will be considered as free parameters; while M L 1 , M L 2 are "reduced" masses of new charged leptons, m E k ≃ u M L k [6], We recall here important properties of charged lepton parameters used in constructing radiative active neutrino masses. According to [6], physical masses (

the mass bases of left-and right-handed leptons E
. Non-diagonal elements of V L may be large because those of M E are at the SU(2) 1 scale. In contrast, those of V e and W e are at least one order of v φ u , because these elements of V L M E are order of the electroweak scale. Hence, V e and W e are nearly identical when u ≫ v φ . They only play the role of generating light charged lepton masses of e, µ, and τ . Hence, in many cases we can use the approximations We can see that the V L is exactly the mixing matrix of neutrinos if they are all considered as the pure Dirac particles. Formula of V L is written in the block form, namely [6] where analytic expression of V ij L , with i, j = 1, 2, corresponding to λ ℓ in Eq. (20) are given in Appendix A. The Yukawa coupling matrix y ℓ (13) is also mentioned, with a requirement that the SM block of the charged leptons is diagonal after the block-diagonalization. It does not affect results obtained in Ref. [6], which depend mainly on the gauge couplings.
Hereafter, many calculations to discuss on phenomenology of Higgs bosons will ignore small mixing between different flavor quarks. We will apply the same results of the charged lepton sector to the quarks. The equivalences between notations are: V ij L , λ ℓ , M L 1,2 , ∆ µ,τ → V ij Q , λ q , M Q 1,2 , ∆ b,s , which were given in [6]. Next, we will discuss another possibility that neutrinos can get Majorana mass terms.

B. Neutral lepton masses
Unlike charged leptons, where the SM-like charged leptons have their own right-handed partners, the SM-like neutrinos do not. In addition, the neutral leptons may inherit Majorana mass terms, for example 1 2 (ν L ) c m ν ν L for active neutrinos. Hence, it is more convenient to write the mass matrix of neutral leptons in the form discussed in the seesaw models [9], which is different from [6]. At the beginning ν L , N R and N L will be considered as independent fields, where the left-and right-handed bases are N ′I . For n VL = 2, the mass matrix of the neutral leptons is a 7 × 7 symmetric matrix having the following form: where m D ≡ 1 2 λ ℓ u and M L are 3 × 2 and 2 × 2 matrices, respectively. Similarly to seesaw models, (N L ) c and N R are additional right-handed neutrinos. The matrix (23)  A pair of two degenerate values corresponds to one Dirac mass of a heavy Dirac neutrino, the same as that mentioned in [6]. The mixing matrix of neutrinos is derived from (22) as follows: where new neutrino masses are pure Dirac. In addition, new lepton masses in each family are nearly degenerate. Equation (24) gives the relations between the original and mass bases R , which are the same as those of the charged leptons. To keep the lepton spectrum being unchanged and looking for a solution of active neutrino mass problem, the mass terms of active neutrinos must come from the effective Majorana Because the active neutrino masses are tiny, their effect on the mixing parameters with heavy neutrinos is negligible. Based on the mechanism of the neutrino mass generation in the Zee model [8], in this model only one pair of new singly charged Higgs bosons, denoted as δ ± ∼ (1, 1, 1) ±1 carrying even Z 2 charges, is introduced.
New couplings for generating one-loop radiative neutrino masses are where i, j = 1, 2, 3; and k, l = 1, 2. In the general case, k, l = 1, 2, .., n V L . We stress that all terms in (25) are simultaneously survival only when both φ ′ and δ ± are even under Z 2 symmetry.
The terms in the first line of (25) violate the lepton numbers, exactly in the same way as in the Zee model, where φ ′ plays a similar role as the second Higgs doublet. Similarly to the Zee model, the trilinear coupling is λ δ u after the first step of the spontaneous breaking. An one-loop diagram generating active neutrino masses is shown in Fig. 1. Following [8,10], the effective mass matrix of light neutrinos is derived in Appendix B, where ϕ ± and δ ± are assumed to be the physical Higgs bosons. But in the model under consideration, ϕ ± and δ ± are not mass eigenstates. As we will discuss later, the physical fields in the Higgs sector are h ± 1,2 , and there are some useful relations: The parameters c ξ , s ξ , c ζ , and s ζ involve with mixing parameters ξ and ζ of the Higgs bosons, defined in Eqs. (60) and (76), as we will present below. The Higgs  (25) can be rewritten as follows: where v φ = vs β , v φ ′ = vc β and t β ≡ tan β = s β /c β , which are defined in [6].
Charged leptons e c in the loop will be considered as mass eigenstates with masses m ec .
Therefore, the Yukawa terms should be written in terms of physical charged lepton states, and light neutrinos are massless states after the rotation V L . For simplicity, we will assume that V 11 L is real and V e , W e ≃ I. In addition, we ignore one-loop contributions to heavy neutrino masses because they are extremely smaller than the tree level masses. Then the one-loop corrections are mainly from light leptons, namely Here a, c, g, h = 1, 2, 3; I, J = 1, 2, ..., 7; f and f ′ are 3 × 3 and 2 × 2 antisymmetric matrices, respectively. Similarly, we have The effective mass matrix m ν of active neutrinos is derived based on (B3), where M e ≡ diag(m e , m µ , m τ ).
Including loop contributions (27), the SM block of (24) will be changed from zero into consisting of three active neutrino masses, and U PMNS is the well-known neutrino mixing matrix. If V 11 L = I 3 and V 12 L = 0, the Eq. (27) is the same (as B3). Like in the Zee models, the parameters arising from the Higgs sector affect the order of the neutrino masses only. But the masses and mixing angles of active neutrinos depend on unknown parameters in f, f ′ , and V 11 L . As a result, the model under consideration is less restrictive in fitting the neutrino data than the Zee models. Because these models are still valid [11], the neutrino sector mentioned here is realistic. In general, fitting recent neutrino data needs at least five free parameters, in agreement with three mixing angles and two squared mass differences. Because two of four parameters, namely m 0 and three f ab , determine the order of the lightest neutrino mass, there are two free parameters left. When n V L ≥ 3, there are at least three additional parameters f ′ kl , enough for fitting neutrino data without constraints on V 11 L . Interestingly, the neutrino fitting results in [11] would be applied to the model under consider ration if L L carries even Z 2 charge which will survive the lepton coupling matrixỹ ℓ in (16).
L may be involved with fitting neutrino data. Our numerical investigation showed that the allowed regions in Ref. [6], controlled by the texture λ ℓ (20), seems much more constrained.
Note that the m ν in (27) keep only main contributions from loops containing light charged leptons, where mixing terms with order O(ǫ) are ignored. With u around 1 TeV and light new charged leptons, contributions from these lepton mediations to m ν will be significant, implying that their masses can be free parameters for fitting neutrinos data without much changes of ∆ µ,τ . Finding exact allowed regions should be done elsewhere.
When the neutrino data is fitted, the results in Ref. [6] for B-decay anomalies are still unchanged because the analysis considered here addresses only effects of tree contri-butions from heavy gauge bosons, where other contributions from the light lepton masses are suppressed. The unique changes may come from the gauge couplings of active neutrinos with charged gauge bosons. Following [6], after the block-diagonalization these gauge couplings are proportional to W µ l ν L γ µ e L and W µ h ν L γ µ ∆ ℓ e L , where W l and W h are light and heavy charged gauge bosons. In the neutrino mass basis they become W µ l ν L U † PMNS γ µ e L and W µ h ν L U † PMNS ∆ ℓ γ µ e L , resulting in the same factor (U † PMNS ) ii for coupling ν i e i with a diagonal ∆ ℓ obtained from the texture of λ ℓ in (20). This factor will not appear in final results presenting the ratios of B-decay anomalies, as given in [6].
In general, active neutrino mass generation from radiative corrections mentioned above affects only the lepton sector. Furthermore, it does not affect mixing parameters controlling the λ ℓ structure at the first breaking step, hence suggests that the orders of numerical values in allowed regions will not change after neutrino data is fitted.
The above discussion just refers to a simple extension that can generate active neutrino masses through radiative corrections. The problems of neutrino masses and DM can be solved by models with more charged Higgs bosons and singlet right-handed neutral leptons, such as [12]. Following the structures of these models, apart from δ ± , at least one pair of singly charged Higgs bosons S ± and a neutral lepton F R ∼ (1, 1, 1) 0 have to be introduced, where S + ∼ (1, 1, 1) 1 . In addition, only S ± and F R are odd under a new Z 2 discrete symme- Therefore, F R can play the role of DM. It has a Majorana mass term of the form 1 2 (F R ) c m F F R . Active neutrinos get mass from loop corrections, which arise from a new Yukawa term, , and a coupling of charged Higgs bosons, 1 4 . This kind of models seems to be less interesting because the origin of neutrino masses is not related to the new leptons.
New ingredients for generating radiative corrections to active neutrinos do not change both results of the gauge sector and LNU discussed in [6], because no new breaking scale contributes to the masses of gauge Higgs bosons. If we add a new SU(2) 1 triplet, denoted as ∆ ∼ (1, 3, 1) 1 , creating an Yukawa term like −Y ∆ (L L ) c iσ 2 ∆L L + H.c., a neutral component of this triplet will develop a non-zero VEV v ∆ , which contributes a new mass term of the to the neutrino mass matrix (23). This matrix has the same form shown in the inverse seesaw models [9,13]. Hence the active neutrino masses will be non-zero. In addition, some new neutrinos may get light masses and play the role of DM [14]. These models seem interesting because they may give connections between the SU(2) 1 leptons with neutrino masses and DM. But the appearance of the new vev v ∆ will contribute to masses and mixing parameters of the Higgs and gauge bosons, consequently it will affect the results shown in [6]. This extension is beyond our scope, and should be thoroughly studied in another work.
Now we turn to one of the most important elements: gauge bosons.

C. Gauge boson masses
Gauge boson masses arise from the piece where the covariant derivative of Φ is determined as With the help of the notation contributions to masses of gauge bosons are From this, masses and eigenstates of gauge bosons can be found in agreement with those presented in Ref. [6]. We will review important aspect then discuss some new properties when masses and mixing angles are calculated up to order of O(ǫ 2 ).

D. Neutral gauge bosons
In the basis (W 1 3 , W 2 3 , B) the squared mass matrix of neutral gauge boson is M 2 nb . At the first step, where v φ , v φ ′ → 0, only two states W 1 3 and W 2 3 are rotated through a rotation C 1 where n 1 = n 2 1 + n 2 2 was used already in [6]. The first breaking step implies the following transformation of the neutral gauge bosons: [6]. Note that v ≃ 246 GeV, g ′ = gs W /c W , and s W is the sine of the Weinberg angle. From now on, n 1 , g 1 , g 2 and n 2 will be written as At the second step, the mixing matrix C 2 is the SM rotation of only B and W 3 , giving new where A and Z l are the photon and SM gauge boson. The respective matrix M ′2 nb is The mass eigenstates (Z, Z ′ ) relates with the Z l − Z h mixing angle defined as where ǫ ≡ v u . The Z l − Z h mixing vanishes when β ′ = β, where tan β = s β /c β . The masses of the physical eigenstates (Z, Z ′ ) are The relation between the two bases (W 1 , W 2 , B) and (A, where (C 2 C 1 ) T is the first matrix in the right hand side of (35).
Using the new notations of (32), the parameter ζ in [6] can be expressed as ζ ≡ s 2 β − 2β ′ u 2 and g n 2 s β ′ , we can deduce an approximate form ξ Z ≃ 1/2 tan 2ξ Z in the limit ǫ ≪ 1, consistent with the expression of tan 2ξ Z shown in (34).

E. Charged gauge bosons
In the basis (W + 1 , W + 2 ) the squared mass matrix of charged gauge bosons was given in Ref [6]. Setting v = 0, we can define a new basis: where the corresponding squared mass matrix is The SM-like boson W ± is identified with W ± ≡ W ± l with mass m W l = gv/2. The mixing W + l − W + h is defined through the mixing angle ξ W satisfying From (37), it follows that the ratio of the tangents of W − W ′ and Z − Z ′ mixing angles is c W . This will reduce the number of parameters in the model by 1.
The physical mass eigenstates (W ± , W ′± ) are given by  with c ξ W ≡ cos ξ W , s ξ W ≡ sin ξ W , and masses Note that Z and W are the SM-like gauge bosons.
We will derive the approximate formulas for the mixing angles and masses of the SMlike gauge boson up to the order of v 2 × O(ǫ 2 ) because the corrections at this order to the masses may contribute significantly to precision tests such as the ρ parameter. Because (34) and (37), we get This means that s 2 For this reason, the masses of the gauge bosons in (35) and (39) can be written as Then we have In addition, at the tree level the ρ parameter satisfies Using e = g ′ c W = gs W , etc., expression in (42) gives the well-known electromagnetic current where Q is the electric charge operator defined in (6).
Neutral currents are defined as where J µ (Z l ) and J µ (Z h ) can be found from (43) and (44). Remind that physical neutral gauge bosons are Z and Z ′ defined from the Z l − Z h mixing angle (34), leading to the respective neutral currents and The second term in (45) is the NP contribution.
Let us write explicitly the neutral current of the Z boson where T 1 3 = 0 and T 2 3 = 1 2 σ 3 for the SM fermion doublets, and T 1 3 = 1 2 σ 3 , T 2 3 = 0 for the extra fermion doublets. Only interactions in the first line of (47) are the SM ones. The remaining provides NP effects. Note that interactions of new vector-like fermions include both P L and P R parts (as vector).
Let us write the couplings of the Z boson with physical fermion states in the form where P L,R = (1 ∓ γ 5 )/2. The couplings g L and g R are listed in Table I, where we denote Here we keep only significant contributions to g R , in which they contain both factors of heavy masses and ǫ, as shown in the two last lines of Table I. We can see that although new fermions are all vector-like in the flavor bases, they are not vector-like in the mass bases because they are mixed with the chiral SU(2) 2 leptons through Yukawa interactions (16).
In contrast to [6], in our work the neutral currents are written in the basis of physical neutral gauge bosons, SM Z and extra Z ′ , from which their decays can easily be studied.

A. Charged currents
The Lagrangian of charged currents is In the physical states of charged gauge bosons, it is If the W boson part of the Lagrangian is written as L = gc ξ W √ 2 W µ f γ µ (g L P L + g R P R ) f ′ + H.c., the couplings of W boson with physical fermions are shown in Table II. New-physics interactions are in (51). Within the experimental data on the W decay width, ones can get constraints on the mixing angles. That was discussed in detail in Ref. [6].

From
the potential is given as Because the µ parameter is proportional to the squared masses of the charged and CP-odd Higgs bosons, it must be positive with the minus sign before it in the potential (54).
The neutral scalars are expanded as At the tree level, the minimum conditions of the Higgs potential are similar to the ones in Ref. [6], except the opposite signs of µ.
Based on the minimum conditions, the parameter µ 2 φ,φ ′ ,Φ can be expressed as a function of the Higgs-self couplings u, v and β. Next, the masses, mass eigenstates, and couplings of Higgs bosons will be calculated by inserting these functions into the Higgs potential (54).

A. Squared mass matrices of the Higgs bosons
In the original bases of singly charged and CP-odd neutral Higgs bosons φ ± = (ϕ ± , ϕ ′± , Φ ± ) T and A = (A φ , A φ ′ , A Φ ) T , the corresponding squared mass matrices are In the basis of CP-even Higgs bosons S = (S φ , S φ ′ , S Φ ) T the squared mass matrix M 2 S corresponding the mass term 1 The above matrices are consistent with those given in [6] after using the relations (56).

B. Physical spectrum of Higgs bosons and their couplings
We will find the Higgs bosons masses in two steps. At the first step, where v → 0, all the three squared mass matrices are diagonalized through the same transformation In the second step, it is easy to determine rotations diagonalizing the squared mass matrices of charged and CP-odd neutral Higgs bosons. By defining the mixing angle ζ as the total mixing matrices used to diagonalize mass matrices in (57) are Mass eigenstates of the charged and CP-odd Higgs bosons, denoted as H ± = (G ± 1 , G ± 2 , h ± ) T and H A = (G Z 1 , G Z 2 , h a ) T , relate with the original states through the following equations: Two linear combinations of G ± 1 and G ± 2 are Goldstone bosons eaten up by the W ′± and W ± gauge bosons. Similarly, linear combinations of G Z 1 and G Z 2 are eaten up by Z and Z ′ . There are two physical charged Higgs bosons h ± and one physical CP-odd neutral Higgs boson h a with masses m 2 h ± and m 2 A , respectively. They satisfy Regarding the CP-even neutral Higgs bosons, after the rotation (59), the squared mass matrix is M ′2 S = C 1 M 2 S C T 1 , which is a 3 × 3 matrix with following elements: In general, M ′2 S is complicated and it cannot be diagonalized exactly. Instead, using the parameter ǫ ≡ v/u ≪ 1, we will find approximate solutions for mass eigenvalues, keep terms up to the order of the electroweak scale. This is reasonable because the SM-like Higgs boson mass was found to be 125 GeV. Approximate solutions was used earlier to find consistent masses of the lightest CP-even neutral Higgs bosons in supersymmetric models [22]. The mixing matrix will also be determined approximately, corresponding to the mass eigenvalues.
We start from finding the eigenvalues of the matrix M ′2 S by solving the equation Det (M ′2 S − λI 3 ) = 0, where λ is expanded as λ = u 2 (λ 0 + λ 1 ǫ 2 ) to keep it up to the order of electroweak scale v. We assume that λ 0 , λ 1 ∼ O(1). Using v = uǫ, we can write where a 0 = a 0 (λ 0 ) and a 1 = a 1 (λ 0 , λ 1 ). We will consider only the two following equations: The first equation in (65) shows that the largest contributions to Higgs masses are the solutions of a 0 (λ 0 ) = 0, giving one zero and two non-zero values, λ 0 = µu 2s β c β and λ 0 = λ 3 u 2 . Hence there are two heavy CP-even neutral Higgs bosons with the corresponding masses , which equal the largest contributions of the two last diagonal entries of M ′2 S shown in (64). A light CP-even neutral Higgs boson corresponds to λ 0 = 0. Its mass comes from the second equation of (65): We stress that m 2 It can be checked that after this rotation the light Higgs boson mass is consistent with (66). Therefore, the mixing matrix relating two original and physical bases S and H 0 = The light Higgs boson h 0 1 is identified with the SM-like Higgs boson found by the LHC. The recent experimental data shows that the SM predictions perfectly agree with the observation within 1 sigma [16]. Hence, couplings of h 0 1 with other SM particles must be consistent with this data. The relevant couplings of h 0 1 are shown in Table III, including couplings with h ± 2 needed to generate active neutrino masses. We can see easily that all couplings with the SM-like particles are different from the SM predictions by a common factor c h . So, |c h | should be close to unity, i.e., |s h | should be small. Its upper bound can be found as follows. Consider the h 0 1 productions at LHC, new heavy quarks can play the roles of the top quark in gluon-gluon fusion mechanism, where their couplings are proportional to s h or ǫ. A significant contribution related to ǫ may come from the quarks U 2 where the couplings contain a factor (∆ b m t ) 2 ǫ/v 2 . But the constraint from [6] gives (∆ b m t ) 2 ǫ/v 2 ∼ 10 −4 ǫ, which is suppressed. Now, the production of h 0 1 through gluon-gluon fusion at lowest order is [23] σ 0 where t q = for SM-like quarks; new quarks U 1 , D 1 ; and new quarks U 2 , D 2 , respectively. The form factor where for t > 1.
Using m Q 1,2 , m Q 1,2 ≃ M Q 1,2 u, as given in [6], Eq. (69) is written as where t 0 = . The condition m t , m Q 1,2 > m h 0 1 = 125.09 GeV gives the limit A 1/2 (t) → 4/3 for all t = 0, 1, 2. The respective signal strength of Higgs production is where we follow the notations of signal strengths defined in [16]. Similarly, the partial decay width of the channel h 0 1 → gg is determined as where µ hgg = µ ggF . Because s h ǫ = O(ǫ 2 ) and the branching ratio of this decay is smaller than 9%, we will use the naive approximation µ ggF ≃ c 2 h to find a lower bound of |c h |. For all remaining decay channels of the SM-like Higgs boson into SM particles, the treelevel couplings are always different from the SM prediction the factor c h , therefore µ f = c 2 h for all main decays f = ff , W W * , ZZ * . The global signal strength defined in [16] can be formulated approximately as follows This gives the constraint 0.995 ≤ |c h | ≤ 1, and |s h | ≤ 0.10.
If we use the constraints of ∆ b,s and c β ′ given in [6], the values of |s h | satisfying (74) are reasonable for the approximation we have just discussed. In addition, the constraint of s h results in small couplings of h 0 1 with the heavy fermions and W ′ gauge boson, giving their suppressed contributions to the decay rate h 0 1 → γγ. Couplings h 0 1 h ± 1,2 h ± 1,2 depend on many unknown Higgs-self couplings, implying that charged Higgs masses are not constrained from the experimental data of the decay h 0 1 → γγ, so we will not consider this decay further.

C. Singly charged Higgs bosons with additional δ ±
In this section we consider the model including new singly charged Higgs bosons δ ± discussed in the neutral lepton sector. Apart from the second term in (25), the Higgs potential has new terms, The appearance of δ ± does not change the allowed regions of parameters discussed in [6].
Also, the results derived for Higgs bosons are unchanged, except the singly charged Higgs sector. In the basis (ϕ ± , ϕ ′± , Φ ± , δ ± ) T , the squared mass matrix is denoted as has only the following non-zero elements: This matrix is diagonalized by a transformation relating with a mixing angle ξ satisfying Then the total transformation can be found to be which changes the original basis into the mass eigenstate basis ( We note that the Goldstone bosons G ± 1,2 defined in (62) are not affected by the presence of δ ± .
The masses m 2 (62) and h ± 2 ≡ δ ± which is used for simple approximations because Eq. (76) means t 2ξ ∼ ǫ ≪ 1. The relevant couplings of the charged Higgs bosons to fermions are collected in Table IV. Couplings of the charged Higgs bosons with gauge bosons are shown in Table V. Only couplings of h ± 1 are shown because the couplings of h ± 2 can be derived by the following replacements: We consider here only the case of ξ → 0. Some important properties of h ± 1 are as follows. Differences between couplings of h ± 1 and the SM-like Higgs bosons to normal fermions are g h ± 1 ℓν ℓ /g h 0 Table V, the couplings of h ± 1 to the SM-like bosons, namely h ± 1 ZW and h ± 1 h 0 1 W , are extremely small because they contain factors s Z ǫ 2 ∼ O(ǫ 4 ), s Z ∼ O(ǫ 2 ), and another mixing smaller than 0.02. Other couplings to light fermions are also small because s ξ → 0. With c ζ , c ξ → 1, the main decays of h ± 1 into light particles are h + 1 → tb. If m h ± 1 > m W ′ , m Z ′ , there will appear two additional large decay modes h + 1 → Z ′ W + , ZW ′+ . In contrast to h ± 1 , the charged Higgs bosons h ± 2 only couple strongly with leptons and Higgs bosons. Therefore, the main decay modes are h + 2 → (ν e i ) c e j with e j = e, µ, τ . The main processes for h ± 2 production at colliders are We would like to compare the above singly charged Higgs bosons with the ones predicted by Zee models, where the charged Higgs sector was investigated thoroughly in Ref. [17]. The  particles. But the predictions for charged Higgs boson production at colliders like the LHC are different, because of the appearance of new particles, such as new heavy quarks and Higgs bosons, and the constraints from allowed regions of parameters indicated in [6]. We will review these regions before discussing the signal of new particles at colliders.

Vertex
Coupling where g = 2m W v ≃ 0.651 is the SM gauge coupling, and ζ ′ satisfies t ξ W = c W t Z ≃ c 3 β ′ s β ′ ζ ′ ǫ 2 [6]. This gives the constraints We can see that the allowed values of ζ give very small values of t Z,ξ W , even with large ǫ < 1.
For simplicity, we will also use the following simple approximations: The simple texture of λ ℓ in (20) gives M L 1 = u M L 1 s β ′ and M L 2 = u M L 2 ρ µτ . For the quark sector, M Q 1 = u M Q 1 s β ′ and M Q 2 = u M Q 2 ρ sb . Now the masses of the heavy particles are

B. Searches for new fermions at colliders
From the above discussion, if the new fermions are lighter than all new bosons, including W ′± , Z ′ , h 0 2,3 , h ± 1,2 , they have only the following three-body decays: • For the first family of new fermions: • For the second family of new fermions: Because of the suppressed ∆ τ,b , the main decay modes are The partial decay widths of decays F → f h 0 1 , f W, f Z are where V = W, Z, g F f V and Y F f h 0 1 are, respectively, couplings of fermions with gauge and Higgs bosons given in Tables I, II and III. Additional factors 3 are included for quark decays.
The decays listed in (84) and (85) always have , with ζ ′ satisfying (77). Hence, every heavy fermion will decay mainly into a light fermion and a SM-like Higgs boson.
Heavy fermions have been being searched for at the LHC recently, for example the heavy lepton decays into pairs of light leptons and the SM-like gauge bosons [18], and the null result is consistent with this investigation. Other heavy quark decays listed in (85) are [20], and U 2 , D 2 → Zt, Zb [21]. But the promoting channels predicted from this discussion are only U 2 → ch 0 1 and D 2 → bh 0 1 .
In conclusion, we have indicated that the allowed regions of parameters given in [6] predict following main fermion decays: According to our knowledge, these decay channels were not treated experimentally. We emphasize that this discussion is valid for heavy fermions lighter than all other new bosons. Any fermions that are heavier than a heavy gauge boson or a Higgs boson will decay mainly into light fermions and this boson.

C. Searches for new Higgs bosons at colliders
At the LHC, the promoting possibility of detecting h 0 2 coupling strongly with heavy fermions was indicated in [6]. These large couplings are shown in Table VI, where only large couplings of neutral CP-even Higgs bosons are shown for investigating Higgs productions.
Vertex coupling We will focus on the remaining new Higgs bosons, including h 0 3 , h a , h ± 1 and h ± 2 . It turns out that they inherit many properties of the new Higgs bosons predicted in THDMs and the MSSM, except h ± 2 . The Higgs sector of the Zee models can be regarded as the one of a THDM plus a pair of singly charged Higgs bosons, as investigated thoroughly in Ref. [17].
And the complete investigation of the Higgs phenomenology of the MSSM was presented in [23] including a brief comparison with Higgs sector in THDMs. The Higgs bosons h 0 3 , h a , and h ± 1 have degenerate masses containing the factor of the trilinear Higgs self-coupling µ. This property is the same as that of the MSSM, but completely different from THDMs. Both Refs. [17] and [23] considered the Yukawa part of the THDM type II, where up and down right-handed singlets of light fermions couple with different Higgs doublets. In contrast in the model under consideration, all right-handed fermions couple with the same Higgs doublet φ. This explains why the couplings of the neutral h 0 3 and h ± 1 with all quarks always contain the same factor 1 t β , and couplings of h a with all SM-like fermions contain the same factor 1/ t β , as shown in Table VII. While, the couplings of up and down quarks in the MSSM and THDM type II have different factors of 1/t β and t β , respectively. The notation β in this

Vertex
Coupling work is equivalent to 1/t β defined in [17,23], where the allowed t β is consistent with the constraint (78).
Now the recent searches for Higgs bosons in THDMs and MSSM will be used for predictions of detecting new Higgs bosons discussed in this work. We consider only Higgs bosons heavier than the top quark. Possible main decays are where . Expressions for couplings Yukawa Y hf f , gauge-Higgs-Higgs g hV V , the Higgs-Higgs-gauge g hhV , and λ hhh were listed in the above Tables. The correlations between the different partial decay widths of a Higgs boson depend only on the last factors of formulas in (87). Hence, they will be used to estimate the largest partial decay widths.
The main decay channels of h ± 1 are h + 1 → tb, Z ′ W, ZW ′ have relative factors as follows: where the allowed values of t β are given in (78). Hence, if m h ± 1 is not too larger than the heavy gauge boson masses, the main decay is h + 1 → tb, where the h ± 1 tb coupling is the same as in the MSSM. The LHC has searched for this decay recently [24,25], through the production channel pp → tbh ± , giving the lower bound of 1 TeV for m h ± 1 .

VI. CONCLUSION
Recently, the G221 model has been introduced in Ref. [6] with the main purpose to explain all experimental data in flavor physics, tau decays, electroweak precision data, and LNU phenomenology from the anomalies in B decays. But there are still to crucial questions to this model, namely, how to generate active neutrino masses and DM? This work indicated that these problems can be solved based on the mechanisms of generating the active neutrino masses by radiative corrections. In particular, the simplest way to generate the active neutrino masses based on the Zee models was shown in detail. The model predicts the existence of a new pair of singly charged Higgs bosons that have large couplings only with light leptons and Higgs bosons. The DM problem can be solved by applying similar mechanisms shown in many radiative neutrino mass models with DM that were widely investigated previously.
In this work we have analyzed a more general diagonalization of gauge boson mass matrices. We have found that the ratio of the tangents of Z − Z ′ and W − W ′ mixing angles is the cosine of the Weinberg angle, cos θ W . This leads to the consequence that the number of the model parameters is reduced by 1. Hence, their behaviors can be predicted based on well-known studies of the THDM as well as of the MSSM.
We combined the above results and the allowed regions of parameters indicated in Ref. [6] to predict some promoting decay channels of new fermions and Higgs bosons. We found that the decays of new heavy particles to SM-like gauge bosons are very suppressed, due to the very small mixing of heavy and SM gauge bosons. The main decays of heavy fermions into two SM-like particles are the decays F 1,2 → h 0 1 f 1,2 . Decays into SM-like fermions in the third family are very suppressed because the allowed regions contain the tiny coefficient ∆ τ,b . The main decay of h ± 1 is h ± → tb. The latest searches for this decay channel give a 1 TeV lower bound for the charged Higgs boson mass.
The LHC have searched for many decay channels of new fermions into SM-like fermions of the third family. So the model will be checked by experiments in coming years. If these decay channels are detected, the model must be extended. For example, the third family of new vector-like fermions should be added to release the allowed regions of parameters. If λ ℓ has the form given in Eq. (20), the precise formula of V 11 L defined in (22) is where ρ µτ = 1 − c 2 β ′ ∆ 2 µ + ∆ 2 τ . Other submatrices contained in V L are After the block-diagonalization, the SM blocks of fermions matrices must satisfy the experimental constraints. In general, the SM block of the charged lepton mass matrix V e V L M E W † e = M ′ E will not be diagonal if the matrix y ℓ in (13) is assumed to be diagonal for simplicity. Instead of, y ℓ is chosen so that only mixing on µ − τ sector is non-zero, There exist values of y µτ,τ µ so that the matrix (A2) is diagonal and the result of [6] is unchanged. The diagonal SM block of charged leptons also guarantees that the lepton flavor violating decay h 0 1 → µτ is suppressed, consistent with experimental constraints. Then y µτ and y τ µ are chosen to satisfy the condition (M ℓ ) 23 = (M ℓ ) 32 = 0. Now the elements of the Yukawa coupling matrix y ℓ can be expressed as