U(1)B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(1)_{B-L}$$\end{document} extension of the standard model with S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_3$$\end{document} symmetry

We propose a renormalizable B-L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B-L$$\end{document} Standard Model (SM) extension based on S3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_3$$\end{document} symmetry which successfully accommodates the observed fermion mass spectra and flavor mixing patterns as well as the CP violating phases. The small masses for the light active neutrinos are generated through a type I seesaw mechanism. The obtained physical parameters in the lepton sector are well consistent with the global fit of neutrino oscillations (Esteban et al. in J High Energy Phys 01:106, 2019) for both normal and inverted neutrino mass orderings. The model also predicts effective neutrino mass parameters of ⟨mee⟩=1.02×10-2eV,mβ=1.25×10-2eV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\langle m_{ee}\rangle }= {1.02\times 10^{-2}}\,{\mathrm {eV}},\, m_{\beta }= {1.25}\times 10^{-2}\,{\mathrm {eV}}$$\end{document} for normal hierarchy (NH) and ⟨mee⟩=5.03×10-2eV,mβ=5.05×10-2eV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\langle m_{ee}\rangle } ={5.03}\times 10^{-2}\, {\mathrm {eV}},\, m_{\beta } ={5.05}\times 10^{-2}\, {\mathrm {eV}}$$\end{document} for inverted hierarchy (IH) which are all well consistent with the future large and ultra-low background liquid scintillator detectors which has been discussed in Ref. (Zhao et al. in Chin Phys C 41(5):053001, 2017) or the limit of the effective neutrino mass can be reached by the planning of future experiments. The model results are consistent with and successfully accommodate the recent experimental values of the physical observables of the quark sector, including the six quark masses, the quark mixing angles and the CP violating phase in the quark sector.


Introduction
Although highly successful in describing most of the elementary particle phenomena at low energy scale, the SM leaves many unresolved issues such as the large hierarchy of charged fermion masses, the fermion mixings, and the a e-mail: vovanvien@tdtu.edu.vn (corresponding author) b e-mail: hnlong@iop.vast.ac.vn c e-mail: antonio.carcamo@usm.cl CP-violating phases, etc. Therefore, these issues can be considered as important evidences of physics beyond the SM.
Among the possible extensions of the SM, the one with an extra U (1) B−L gauge symmetry is one of promising extensions which has been considered in Refs. [1][2][3][4][5][6][7][8]. In this type of model, the presence of three right-handed neutrinos is essential to cancel the gauge and the mixed gauge-gravitational anomalies. However, it does not provide a natural explanation for fermion mixings.
Experimentally, the best-fit values for neutrino squared mass splittings, leptonic mixing angles and the Dirac CP violating phase, for both normal and inverted mass hierarchies with Super-Kamiokande atmospheric (wSK-atm) neutrino data, are given in Ref. [9]: The discrete symmetry has revealed many outstanding advantages in explaining the observed pattern of SM fermionic masses and mixing angles given in Eq. (1). There have been many models based on discrete symmetries, see for instance, Refs. . However, in most of these papers, the fermion masses and mixings are generated from non-renormalizable interactions or at loop levels. The renormalizable B − L model with S 3 symmetry was first proposed by Gömez-Izquierdo and Mondragön [52], however, there are substantial differences between the work in Ref. [52] and our current work as follows 1 : (a) Ref. [52] contains two other sym- metries Z 2 and Z e 2 compared to S 3 symmetry where the couplingsL He R andL H N are prevented by Z e 2 to obtain diagonal charged and neutrinos Dirac mass matrices but it does not modify the Majorana mass matrix form; (b) in Ref. [52], the Yukawa coupling constant y q 2 is assumed equal to zero (a q = 0) which corresponds to an extra discrete symmetry (differs from the existing symmetries Z 2 , Z e 2 and S 3 ) is included to prohibit the Yukawa interaction term [52] the inverted and degenerate neutrino mass hierarchies were considered, however, the recent experimental data [9,53] favors the inverted or normal hierarchy; (d) The Dirac CP-violating phase was also not predicted by [52]; (e) The quark mass matrices M u,d in our model seem simpler in the sense that two entries "23" and "32" are vanished.
The aim of this paper is to construct the renormalizable B − L extension of the SM. In this work, the first generation of quarks and leptons are put in 1 while the two others are put in 2 under S 3 symmetry. This paper is organized as follows. In Sect. 2 we present a simple SM extension by adding U (1) B−L and S 3 symmetries. In Sect. 3 we present the lepton sector of the model and introduce necessary Higgs fields responsible for lepton masses and mixings. Section 4 deals with quark mass and mixing. The implications of our model in the Higgs diphoton decay rate are analyzed in Sect. 5. We conclude in Sect. 6. Appendix A provides a brief description of the Clebsch-Gordan coefficients for the S 3 group.

The model
The electroweak sector of the SM is supplemented by the S 3 ⊗ Z 4 flavor symmetry and a gauge symmetry U(1) B−L . To the particle content, we add three right-handed neutrinos (ν i R ), two SU (2) L scalar doublets H , H with B − L = 0 put in 1 and 2 under S 3 , respectively, and one flavon φ with B − L = 2 put in 2 under S 3 . The particle content of the model is then given in Table 1.
Since there are many Yukawa couplings in Higgs potential, it is easy to arrange a suitable Higgs potential. In order to generate the remarkable fermion mixing pattern, from the potential minimization conditions of the Higgs potential as presented in Appendix B, we choose the VEVs of scalar fields as follows In order to show that the alignment in Eq. (2) is a automatical solution from the minimization conditions of V total , in the system of minimization equations, let us put v * The system of Eqs. . It is noted that this aligned is only one solution to have the desirable results. Note for the low energy observables, but lowering that scale, for example at the TeV scale, will imply lower values for the neutrino Yukawa couplings needed to reproduce the tiny values of the light active neutrino masses. The scale of breaking of the B − L symmetry does not have implications in the quark sector given that the only scalar field charged under this symmetry is a gauge singlet that only appear in the Yukawa terms for the right handed Majorana neutrinos. Independently on the magnitude of the B − L breaking scale, the vacuum expectation values of the Higgs doublets should satisfy the relation v 2 GeV, which is needed to correctly reproduce the measured values of the W and Z gauge boson masses. If the B − L/S 3 breaking scale is very high, some collider signatures that can be tested in the model are: the production of neutral scalar via gluon fusion mechanism, the associated production of a neutral scalar with a SM gauge boson (via Vector Boson Fusion or Drell-Yan mechanism), the associated production of a charged scalar with a SM gauge boson (via Vector Boson Fusion or Drell-Yan mechanism). Furthermore, flavor violating decay like t → hc (h being the 126 GeV Higgs), h → μτ , μ → eγ can be use to constrain the model parameter space. Some low energy observables that can be used to constrain the model are the muon and electron anomalous magnetic moment, the oblique T and S parameters. Rather than avoiding, to make such process below the current experimental limits, some restrictions on the model parameter space need to be imposed, which could be constraints on the Yukawa couplings, on the mixing angles in the scalar sector and large masses for non SM scalars. Besides that, notice that our model has many free parameters, which allows us freedom to assume that the remaining scalars are sufficiently heavy to fulfill the current experimental bounds. Moreover, the loop effects of the heavy scalars contributing to precision observables can be suppressed by making an appropriate choice of the free parameters in the scalar potential. These adjustments do not affect the physical observables in the quark and lepton sectors, which are determined mainly by the Yukawa couplings. A detailed study of the phenomenology of the model is beyond the scope of this paper and is deferred for a future work.
The Yukawa interactions in charged lepton sector are: With the help of Eqs. (7) and (2), the Lagrangian mass term of the charged leptons can be written in the form: where which has the diagonal form. Thus, the charged lepton diagonalization matrices are U l L = U l R = 1 and the lepton mixing matrix depends only on that of the neutrino.
By comparing Eq. (9) with the experimental values for masses of the charged leptons given in Ref. [53], m e 0.51099 MeV, m μ 105.65837 MeV, m τ 1776.86 MeV, The neutrino masses arise from the couplings transforms as SU (2) L doublets (singlets) and 1 + 1 + 2, 2 and 1 under S 3 , respectively. For the known SU (2) L scalar doublets, available interactions areψ 1L H ν 1R , ψ αL ν α R 1 H and ψ αL ν α R 1 H but they only generate Dirac mass terms. To generate Majorana mass terms for neutrinos we will therefore introduce one new SU (2) L singlets φ put in 2 under S 3 coupling toψ 1L ν 1L andψ αL ν αL . The Yukawa interactions, which are invariant under all the symmetries of the model, in neutrino sector are: With the VEVs given in Eq.
(2), we get the Dirac neutrino mass matrix (M D ) and the right-handed Majorana neutrino mass matrix (M R ) as follows where The effective neutrino mass matrix, obtained through the type-I seesaw mechanism, reads: where It is noted that A 0 , B 0 and C 0 given in Eq. (14) accommodated in the first matrix of Eq. (13) due to the contribution of H and φ only while the last term in Eq. (13) is deviation from the contribution of H only. If there is no contribution of H , the deviations a 1,2 and b 1 will vanish and the matrix M eff in (13) reduces to its first term which generates a μ − τ mixing form. Thus, second term in (13) will take the role for a small deviation of θ 13 and being responsible for the CP violating phase in the lepton sector. The first matrix in Eq. (13) has three eigenvalues, and the corresponding lepton mixing matrix takes the form: where with A 0 , B 0 and m 0 1,2 are defined in Eqs. (14) and (16), respectively.
The matrix U 0 in Eq. (17) implies θ 13 = 0, θ 23 = π/4 and θ 12 = θ which was rule out by the recent data. 3 However, the contribution of the second term in Eq. (13) will improve this. Indeed, at the first order of perturbation theory, the second matrix in Eq. (13) contributes to both eigenvalues and eigenvectors. In this case, the neutrino masses are given as: with b 1 , m 0 1,2 and θ are given in Eqs. (15), (16) and (18), respectively. The corresponding lepton mixing matrix becomes: where U 0 is defined by (17) and δU has the following entries: In the three-neutrino framework, lepton mixing angles can be defined via the elements of the neutrino mixing matrix: where t 12 = s 12 /c 12 , t 23 = s 23 /c 23 , c i j = cos θ i j , s i j = sin θ i j with θ 12 , θ 23 and θ 13 , respectively, being the solar angle, atmospheric angle and the reactor angle and δ is the Dirac CP violating phase. From Eqs. (17), (20)- (22), we found that, in both normal and inverted hierarchies, U 11 and U 12 depend only on two parameters θ, θ 12 : As pointed out in Ref. [9], at 3σ confidence level, U 11 ∈ (0.797, 0.842). Thus, we get a value cos θ = 0.831 (i.e., t 12 = 0.670, θ = 33.83 • ) and U 12 = 0.557, which is consistent with the best-fit value from the global analysis 4 given in Ref. [9]. Furthermore, Eqs. (17), (20)-(22) provide a solution: , Although the global analysis in Ref.
[58] shows a hint in favor of the NH over the inverted one at more than 3σ and the global analysis in Ref. [9] obtain a preference for NH at about 2σ , but, it is the fact that the neutrino mass spectrum is currently unknown and it can be NH or IH depending on the sign of m 2 31 [53]. In the model under consideration, both NH and IH can be found in which the model parameters are in good agreement with the global analysis in Ref. [9].

Normal spectrum
By taking the best-fit values for neutrino squared mass splittings, leptonic mixing angles and the Dirac CP violating phase for NH with Super-Kamiokande atmospheric neutrino data as given in Eq. (1), we get a solution 5 : where Now, U 2i and U 3i (i = 1, 2, 3) depend only on m 0 3 . In the NH, the estimated value of m 3 is [53] m 3 0.0506 eV thus we find an allowed region of m 0 3 that can reach the constraint on the absolute values of the entries of the lepton mixing matrix given in Ref. [9] which has been depicted in Fig. 1.
In the NH, m 1 < m 2 < m 3 thus m 1 ≡ m N light is the lightest neutrino mass. The effective neutrino masses governing the beta decay and neutrinoless double beta decay [59][60][61][62][63] ei m i and m N light as functions of m 0 3 has been plotted in Fig. 2. In the case m 0 3 = 0.051 eV the other parameters are found in Table 2.
The absolute neutrino mass is found to be 3 i=1 m i = 7.2 × 10 −2 eV which is well consistent with the strongest bound from cosmology, m ν < 0.078 eV [64]. Furthermore, the effective neutrino masses m β = 1.25 × 10 −2 eV and m ee = 1.02 × 10 −2 eV which are well consistent with the results in Ref. [65].
which is consistent with the constraint on the absolute values of the entries of the lepton mixing matrix given in Ref. [9]. The Jarlskog invariant is J N C P = Im(U 23 U * 13 U 12 U * 22 ) = −0.0245. The resulting effective neutrino mass parameter m ee for normal hierarchy is below the upper bound arising from present 0νββ decay experiments. However, it is very well consistent with the future large and ultra-low background liquid scintillator detectors which has been discussed in Ref. [66] or the meV limit of the effective neutrino mass can be reached by the planning of future experiments [67][68][69][70][71][72][73][74].

Inverted spectrum
Similar to the normal spectrum, taking the best-fit values of neutrino oscillation parameters for IH with Super-Kamiokande atmospheric neutrino data as given in Eq. (1), The elements of the lepton mixing matrix U I 2i and U I 3i (i = 1, 2, 3) depend only on m 0 2 . In the IH, the estimated value [53] of m 2 is 0.0504 eV thus we find an allowed region of m 0 2 that can reach the constraint on the absolute values of the entries of the lepton mixing matrix given in Ref. [9] which has been depicted in Fig. 3.
In the IH, m 3 < m 1 < m 2 thus m 3 ≡ m I light is the lightest neutrino mass. The effective neutrino masses governing the beta decay and neutrinoless double beta decay and m I light as functions of m 0 2 has been plotted in Fig. 4. 6 As before, there are four solutions for m 1,3 however they give the same value for m 2 i j and the same absolute values of m 1,2,3 thus we only consider in detail the solution in Eqs. (29) and (30).
In the case m 0 2 = 0.051 eV the other parameters are found in Table 3. The which is consistent with the constraint on the absolute values of the entries of the lepton mixing matrix given in Ref. [9]. The resulting effective neutrino mass parameter for inverted hierarchy is very well consistent with the future large and ultra-low background liquid scintillator detectors which has been discussed in Ref. [66] or the meV limit of the effective neutrino mass can be reached by the above mentioned future experiments [67][68][69][70][71][72][73][74].

Quark mass and mixing
The Yukawa interactions in quark sector are: With the VEV alignments of H, H and H in Eq. (2), the mass Lagrangian of quarks reads where M u and M d take the form: with It is noted that without contribution of H , the matrices M u , M d in Eq. (34) become diagonal ones, i.e, the quark mixing matrices are identity matrices, however, the quark masses can well fit to the experimental data given in Ref. [53].
The quark mass matrices satisfy the following relation: where Here q = u, d, x iq = |a iq |, κ iq = arg(a iq ) (i = 1, 2, . . . , 5) and X q , Y q , Z q , W q , R q and S q are real parameters. In order to fit the measured values of the SM quark masses and CKM parameters given in Refs. [53,75] as shown in Tables 4 and 5, 7 we proceed by solving the eigenvalue problem for the SM quark masses. The following solution has been found: This show that our model is consistent with and successfully accommodate the experimental values of the physical observables of the quark sector: the six quark masses,the quark mixing angles and the CP violating phase in the quark sector. 7 To obtain the best-fit values of (U C K M ) bf i j in Tables 4 and 5 we use the Wolfenstein parameters given by Ref. [75]: λ C = 0.2250,ρ = 0.148,η = 0.348 and A = 0.826. Table 4 The best-fit results for quark masses and mixings taken from [76] and [53,75] Table 5 The best-fit results for CKM parameters taken from Refs. [53,75] Observable Experimental best-fit value

Higgs diphoton decay rate
The decay rate for the h → γ γ process takes the form: and a hW W are the deviation factors from the SM Higgs-top quark coupling and the SM Higgs-W gauge boson coupling, respectively (in the SM these factors are unity). Such deviation factors are close to unity in our model, which is a consequence of the numerical analysis of its scalar, Yukawa and gauge sectors.
Furthermore, F 1/2 (z) and F 1 (z) are the dimensionless loop factors for spin-1/2 and spin-1 particles running in the internal lines of the loops. They are given by: with In order to study the implications of our model in the decay of the 126 GeV Higgs into a photon pair, one introduces the Higgs diphoton signal strength R γ γ , which is defined as: That Higgs diphoton signal strength, normalizes the γ γ signal predicted by our model in relation to the one given by the SM. Here we have used the fact that in our model, single Higgs production is also dominated by gluon fusion as in the Standard Model. The ratio R γ γ has been measured by CMS and ATLAS collaborations with the best fit signals [77,78]: 17 −0.14 and R AT L AS γ γ = 0.96 ± 0.14.
The correlations of the Higgs diphoton signal strength with the a hW W deviation factor from the SM Higgs-W gauge boson coupling and with charged scalar masses are shown in Figs. 5 and 6, respectively. Such correlations indicate that our model successfully accommodates the current Higgs diphoton decay rate constraints.

Conclusions
We have proposed a renormalizable B − L standard model extension based on the S 3 symmetry which successfully accommodates the current data on SM fermion masses, fermionic mixing angles and CP violating phases of both quark and lepton sectors. The tiny values of the light active neutrino masses are generated through a type I seesaw mechanism. The obtained physical observables of the lepton sector are well consistent with the global fit of neutrino oscillation experiments [9] for both normal and inverted neutrino mass hierarchies. The model also predicts effective neutrino mass parameters of m ee = 1.02 × 10 −2 eV, m β = 1.25 × 10 −2 eV for NH and m ee = 5.03 × 10 −2 eV, m β = 5.05 × 10 −2 eV for IH which are all well consistent with the future large and ultra-low background liquid scintillator detectors which has been discussed in Ref. [66] or the limit of the effective neutrino mass can be reached by the planning of future experiments. The model results are consistent with and successfully accommodate the recent experimental values of the physical observables of the quark sector, including the six quark masses, the quark mixing angles and the CP violating phase in the quark sector. Finally, we have also shown that our model successfully accommodates the current Higgs diphoton decay rate constraints.
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: The only experimental data used in our work are the ones reported in Table. Table IV, V, Eq. (45) as well as the experimental values of the charged lepton masses and the best-fit values of neutrino oscillation parameters given in Eq.
(1) that we use to compare the predictions of our model in SM fermion masses and mixings with the experimental results.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .