Probing neutrino and Higgs sectors in SU(2)1×SU(2)2×U(1)Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{ SU(2) }_1 \times \text{ SU(2) }_2 \times \text{ U(1) }_Y $$\end{document} model with lepton-flavor non-universality

The neutrino and Higgs sectors in the SU(2)1×SU(2)2×U(1)Y\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{ SU(2) }_1 \times \text{ SU(2) }_2 \times \text{ U(1) }_Y $$\end{document} 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 the 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 μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu $$\end{document}. 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 ch\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c_h$$\end{document}, which must satisfy the recent global fit of experimental data, namely 0.995<|ch|<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.995<|c_h|<1$$\end{document}. 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′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W'$$\end{document} and Z–Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z'$$\end{document} 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.


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: = 0.310 ± 0.015 ± 0.008, 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). The second kind of anomalies 1 is the interesting LNU ratio reported recently by LHCb [2], namely, = 0.745 +0.090 −0.074 ± 0.036, which has 2.6σ deviation from the SM value R K = 1.0003± 0.0001 in the dilepton mass squared bin (1 ≤ q 2 ≤ 6) GeV 2 .
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 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][6][7][8][9][10][11].
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 [12,13] (more kinds of G221 models can be found in Ref. [14]). In Refs. [12,13] 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 Refs. [12,13]. 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 Refs. [12,13]. 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 [15,16], where a pair of singly charged scalars transformed as singlets under both 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 1 The SM value for R K has been first obtained in Ref. [3].
φ 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 Refs. [12,13] 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 SM-like 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 [12,13], 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 Sect. 2, after a brief review of the model, we present a more careful consideration of charged lepton masses and the Zee method for the generation of neutrino masses. In Sect. 2.2 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 Refs. [12,13]. The second way is different, because a new scalar SU(2) 2 triplet is included. It contains a new neutral component with nonzero vacuum expectation value (VEV), leading to a new mass scale contributing to gauge boson masses. But this model may predict some active neutrinos playing the role of DM candidates. A more careful diagonalization of squared mass matrices and mixing parameters of gauge bosons is presented in Sect. 2.3. In this section, the relation between the tangents of the W -W and Z -Z mixing angles, is derived. Then a validation of ρ parameter under recent experimental constraint will be shown at TeV scale of SU(2) 1 breaking scale. Section 3 is devoted to charged and neutral currents in the model. Here we notice their difference from the SM ones. For the further discussion of the NP searches at colliders, the couplings of Z and W gauge bosons with fermions are given. In Sect. 4 a detailed analysis of the Higgs sector is presented. This section covers both versions of the Higgs sector content without and with the mentioned charged scalars. And the SM-like Higgs boson is identified. Interesting properties of singly charged Higgs bosons are also discussed. In Sect. 5, we review briefly the allowed regions of parameters given in [12,13], 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 Sect. 6.

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 [12,13]: where i = 1, 2, 3 is the SU(2) index. All the chiral fermions transform as q L ∼ 3,1,2, 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 vectorlike 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 [12,13]. The Higgs sector consists of two doublets φ and φ and one self-dual bidoublet (i.e., = σ 2 * σ 2 where σ 2 is the usual Pauli matrix) 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 [12,13]. 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 [12,13]. Hence, we will summarize here only important results and focus on new features of generating active neutrino masses from loop corrections.

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 [12,13]. 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 blockeddiagonalized perturbatively via two steps. After that, the SM parts are separated from the total. The transformations of fermion states are: [12,13]. 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 , They were listed precisely in [12,13]. After the two steps, all original mass matrices in (19) will be transformed into block-diagonal formŝ One of the blocks in everyM F is identified with a SM fermion block, which is diagonalized by 3×3 unitary transformations: Only the CKM matrix, V CKM = S u S † d , appears in the gauge couplings [12,13]. We can fix S e = U e = I 3 .
For studying Higgs boson phenomenology satisfying the allowed regions of parameters given in [12,13], 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 Refs. [12,13], 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 [12,13], We recall here important properties of charged lepton parameters used in constructing radiative active neutrino masses. According to [12,13] , 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 [12,13] where analytic expression of V i j 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 Refs. [12,13], 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 i j L , λ , M L 1,2 , μ,τ → V i j Q , λ q , M Q 1,2 , b,s , which were given in [12,13].
Next, we will discuss another possibility that neutrinos can get Majorana mass terms.

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 [17], which is different from [12,13]. At the beginning ν L , N R and N L will be considered as independent fields, where the leftand right-handed bases are c. 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) can be generally diagonalized through the transformation A pair of two degenerate values corresponds to one Dirac mass of a heavy Dirac neutrino, the same as that mentioned in [12,13]. 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 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 terms 1 2 (ν L ) c m ν ν L + h.c. 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 [15,16], 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 [15,18,19], 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 coupling in (25) can be rewritten as follows: where v φ = vs β , v φ = vc β and t β ≡ tan β = s β /c β , which are defined in [12,13].
Charged leptons e c in the loop will be considered as mass eigenstates with masses m e c . 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 m ν : consisting of three active neutrino masses, and U PMNS is the wellknown neutrino mixing matrix. If V 11 L = I 3 and V 12 L = 0, 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 [20,21], 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 [20,21] would be applied to the model under consideration if L L carries even Z 2 charge which will survive the lepton coupling matrixỹ in (16).
Regarding n V L = 2, there is only one parameter f 12 = − f 21 . Therefore two parameters in V 11 L may be involved with fitting neutrino data. Our numerical investigation showed that the allowed regions in Refs. [12,13], 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 Refs. [12,13] for B-decay anomalies are still unchanged because the analysis considered here addresses only effects of tree contributions 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 [12,13], 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 [12,13].
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 [22][23][24]. Following the structures of these models, apart from δ ± , at least one pair of singly charged Higgs bosons S ± and a neutral lepton 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 [12,13], 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 form 1 2 μ X N L (N L ) c + H.c. to the neutrino mass matrix (23). This matrix has the same form shown in the inverse seesaw models [17,[25][26][27]. Hence the active neutrino masses will be nonzero. In addition, some new neutrinos may get light masses and play the role of DM [28][29][30]. 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 [12,13]. 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.

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 Refs. [12,13]. We will review important aspect then discuss some new properties when masses and mixing angles are calculated up to order of O( 2 ).

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 . Elements of C 1 depend on a mixing angle β defined by where n 1 = n 2 1 + n 2 2 was used already in [12,13]. The first breaking step implies the following transformation of the neutral gauge bosons: 2 3 , where g = g 1 g 2 /n 1 and g are identified as the SM gauge couplings; v φ = vs β , v φ = vc β [12,13]. 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 a new basis (A, 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 , 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 [12,13] can be expressed as ζ ≡ s 2 β − 2β u 2 and g n 2 , 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).

Charged gauge bosons
In the basis (W + 1 , W + 2 ) the squared mass matrix of charged gauge bosons was given in Refs. [12,13]. Setting v = 0, we can define a new basis: where the corresponding squared mass matrix is 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 SM-like 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 t 2Z , t 2ξ W ∼ 2 , we have (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

Currents
The Lagrangian L fermion = i ff γ μ D μ f contains interactions of the gauge bosons with the fermions. Let us firstly consider neutral currents. From Eq. (35), one gets Using e = g c W = gs W , etc., the expression in (42) gives the well-known electromagnetic current L em = A μ J μ em = e A μ ff γ μ Q f , 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). Recall 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 1, 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 1. 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 Table 1 Couplings of Z boson with fermions, = e, μ, τ and e = d = 1 mixed with the chiral SU(2) 2 leptons through Yukawa interactions (16).
In contrast to [12,13], 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.

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 the W boson with physical fermions are shown in Table 2.
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 Refs. [12,13].

Higgs sector
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 Refs. [12,13], 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).

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 The above matrices are consistent with those given in [12,13] after using the relations (56).

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 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 the 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 [40,41]. 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 0 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 m 2 , 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 [32]. 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 3, 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 [12,13] gives where we follow the notations of signal strengths defined in [32]. Similarly, the partial decay width of the channel h 0 1 → gg is determined as where μ hgg = μ gg F . Because s h = O( 2 ) and the branching ratio of this decay is smaller than 9%, we will use the naive approximation μ gg F 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 tree-level 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 * , Z Z * . The global signal strength defined in [32] 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 [12,13], the values of |s h | satisfying (74) are reasonable for the approxi-mation 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.

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 [12,13]. 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 M 2 h ± . We can find a matrix C h ± so that M 2 h ± = C h ± M 2 h ± C h ± T 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  Table 4.
Couplings of the charged Higgs bosons with gauge bosons are shown in Table 5. Only couplings of h ± 1 are shown because the couplings of h ± 2 can be derived by the following replacements: c ξ h ± 1 → s ξ h ± 2 and s ξ h ± 1 → −c ξ h ± 2 . 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 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. [33]. The equivalent notations are v 1 , the charged components of are Goldstone bosons of W ± . Meanwhile those of φ and φ create two Goldstone bosons of W ± , and two other freedoms that mix with the singlet ones to generate physical states. Hence, the model under consideration and the Zee models predict very similar properties of the charged Higgs couplings to SM-like 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 [12,13]. We will review these regions before discussing the signal of new particles at colliders.

Properties of masses and mixing parameters of new particles
In this work, the results of parameter constraints reported in Refs. [12,13] where g = 2m W v 0.651 is the SM gauge coupling, and ζ satisfies t ξ W = c W t Z c 3 β s β ζ 2 [12,13]. This gives the constraints s β ∈ [0.186, 0.542], t β ∈ [0.24, 0.654], 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 and the mixing parameters The important property is that the model predicts several groups of new heavy particles having same spins and degenerate masses, therefore forbid many decay modes. The new Yukawa couplings generating heavy mass terms are always the same for both up and down components of the SU(2) 1 fermion doublets in the same families. Mass differences come only from the Yukawa terms of the electroweak sector (13), because of the large V Q,L . As shown in (80), a difference between a pair of new up and down fermions is O(m f ). Hence, the top quark may give the largest difference if b m t < 1.4 is not considered. Hence, the mass differences are always smaller than m W . Because of the kinetic condition and the fermion number conservation, a three-body decay of a new fermion must decay to at least a light fermion, namely a SM fermion plus a boson. If a fermion is the lightest among the new particles, it will decay only to a light fermion and a SM boson like W , Z or h 0 1 , leading to large branching ratios, which can be searched by recent colliders.

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: N 1 → ν e h 0 1 , ν e Z , e ± W ∓ ; E ± 1 → e ± h 0 1 , e ± Z , ν e W ∓ ; U 1 → uh 0 1 , u Z, dW + ; • For the second family of new fermions: Because of the suppressed τ,b , the main decay modes 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 1, 2  , 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 [34], and the null result is consistent with this investigation. Other heavy quark decays listed in (85) are U 2 → h 0 1 t [35], U 2 → W b [36][37][38], and U 2 , D 2 → Zt, Zb [39]. 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 [12,13] predict following main fermion decays: E 1 → h 0 1 e, U 1 , D 1 → h 0 1 u, h 0 1 d, E 2 → h 0 1 μ, and U 2 , D 2 → h 0 1 c, h 0 1 b. 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.

. Expressions for couplings
Yukawa Y h f 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, Z W have relative factors as follows: 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.