Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js

Linear seesaw model with a modular S4 flavor symmetry

  • We discuss a linear seesaw model with a field content as minimum as possible, introducing a modular S4 using gauged U(1)BL symmetries. Owing to the rank two neutrino mass matrix, we obtain a vanishing neutrino mass eigenvalue, and only the normal mass hierarchy of neutrinos is favored via the modular S4 symmetry. In our numerical Δχ2 analysis, we especially determine a relatively sharp prediction on the sum of neutrino masses to be approximately 60 meV, in addition to other predictions.
  • The neutrino sector is theoretically unconfirmed in the standard model (SM), because only two mass squared differences and three mixings are experimentally found and the scale of mass is extremely minuscule compared to the other three sectors in the SM. Hence, several scientists expect the neutrino sector to possess new physics. A gauged BL (baryon number minus lepton number) U(1) symmetry; U(1)BL, is a promising prescription for generating such tiny neutrino masses, introducing three right-handed neutrinos with rather heavy masses (MR), and is called a canonical seesaw model [1-4]. Because the mass scale is sometimes expected to be that of a grand unified theory (MGUT1015 GeV) to be small neutrino masses, its scale cannot be verified by our current experiments, and the spontaneous U(1)BL symmetry breaking scale is naturally expected to be the same energy as the cut-off scale; MRMGUT.

    To achieve a successful neutrino mass model within our scale (TeV), heavy neutral fermions (SL) with left-handed chirality are introduced along this line of thought. Inverse seesaw model [5, 6] also requests both of NR and SL and the neutrino mass could be realized within TeV scale. But this model may not require GUT scale. Currently, NR and SL can be embedded into the middle scale with SU(3)CSU(2)LSU(2)RU(1)BL [6-8], which is included in SO(10) group [9]. Then, the neutrino mass matrix would not be suppressed by the middle scale, but by MGUT scale when appropriate charge assignments are assigned for each fields in a supersymmetric theory [9]. Hence, we can test phenomenologies with our current experiments, supposing the middle scale breaking can occur at the TeV scale. This type of model is called the "linear seesaw". Because these models typically require more free parameters than the other three sectors in the SM fermion, flavor symmetries are also frequently introduced in these models, to reduce the parameters and obtain predictions (if possible).

    In 2017, attractive flavor symmetries were proposed in Refs. [10, 11], in which the authors applied modular non-Abelian discrete flavor symmetries to quark and lepton sectors. A remarkable advantage is that any dimensionless coupling can also be transformed as non-trivial representations under these symmetries. Therefore, we do not need so many scalars to determine a predictive mass matrix. Another advantage is that we have a modular weight from the modular origin that can play a role in stabilizing DM when appropriate charge assignments are distributed to each of the fields of the models. Along this line of thought, significantly relevant references have recently appeared in the literature, e.g., A4 [10, 12-48], S3 [49-54], S4 [55-66], A5 [60, 67, 68], double covering of A5 [69-72], larger groups [73], multiple modular symmetries [74], double covering of A4 [75-77], S4 [78, 79], and other types of groups [80-85] in which masses, mixing, and CP phases for the quark and/or lepton have been predicted. For interested readers, we provide some literature reviews, which are useful to understand the non-Abelian group and its applications to flavor structure [86-93]. Moreover, a systematic approach to understanding the origin of CP transformations has been discussed in Ref. [94]; in addition, CP/flavor violation in the models with modular symmetry was discussed in Refs. [95-97], and a possible correction from the Kähler potential was discussed in Ref. [98]. Furthermore, a systematic analysis of the fixed points (stabilizers) was discussed in Ref. [99]. It would be interesting to consider a linear see-saw model with the local U(1)BL(SO(10)) under modular symmetry because these symmetries can be obtained from a string theory. Moreover, the nature of the modular symmetry can be used to realize the linear seesaw mechanism, in addition to constraining the flavor structure.

    In this study, we propose a linear seesaw model under modular S4, using the U(1)BL symmetry, in which we try to construct the predictive model as minimum as possible. Owing to the rank two neutrino mass matrix, we obtain a vanishing neutrino mass eigenvalue. Furthermore, only the normal mass hierarchy of neutrinos is favored via the modular S4 symmetry. In our numerical analysis, we perform Δχ2 analysis in the neutrino sector, considering non-unitarity constraint.

    This remainder of this paper is organized as follows. In Sec. II, we review our model, constructing renormalizable Lagrangian and mass matrices in the lepton sector. Then, we formulate the neutrino mass matrix with rank two, in which we estimate the structure of the neutrino mass matrix in the expansion of modulus. In addition, we derive several observables in the lepton sector. At the end of this section, we discuss the non-unitarity bound. In Sec. III, we perform Δχ2 analysis in the lepton sector, and present some predictions using our model. In Sec. IV, we present a summary and discussion. In the Appendix, we elucidate the modular S4 symmetry.

    In this section, we review our model framework for the linear seesaw mechanism, introducing the BL local Abelian symmetry; U(1)BL and modular A4 symmetry. Regarding the fermion sector, we add two left-handed neutral fermions {SL1,SL2} that belong to the isospin singlet, where they exhibit a zero charge under U(1)BL, {1,1} under S4, and 1 under kI, respectively. In addition, we introduce three right-handed neutral fermions ¯NR that belong to the isospin singlet, where they exhibit 1 charge under U(1)BL, 3 under S4, and 3 under k, respectively. The SM left-handed leptons LL[LLe,LLμ,LLτ]T belong to 1 charge under U(1)BL, 3 under S4, and 1 under k, respectively. However, the SM right-handed leptons {¯eR,[¯μR,¯τR]} belong to the +1 charge under U(1)BL, {1,2} under S4, and {1,3} under kI, respectively.

    Regarding the scalar sector, we adopt two Higgs doublets, H1,H2, and an isospin singlet field φ. The isospin singlet φ has a 1 charge under U(1)BL. Here, H1 is the SM-like Higgs that has a zero charge under U(1)BL and kI while H2 has 1 charge under U(1)BL, and zero under kI. We denote each of the vacuum expectation values (VEVs) to be H1,2[0,v1,2/2]T, and φvφ/2. We summarize our particle content and assignments in Table 1.

    Table 1

    Table 1.  Lepton and boson particle contents and their charge assignments under SU(2)L×U(1)Y×U(1)BL×S4×(kI) where LL[LLe,LLμ,LLτ]T kI represents the number of the modular weight.
    Fermions Scalars
    LL ¯eR,[¯μR,¯τR] ¯NR SL1,SL2 H1 H2 φ
    SU(2)L 2 1 1 1 2 2 1
    U(1)Y 12 1 0 0 12 12 0
    U(1)BL 1 1 1 0 0 1 1
    S4 3' 1, 2 3 1,1 1 1 1
    kI 1 1,3 3 1 0 0 0
    DownLoad: CSV
    Show Table

    Then the valid lepton Yukawa Lagrangian is expressed as

    Llepton=LM+LMD+LMD+LMNS,

    (1)

    where LM denotes the charged lepton Yukawa Lagrangian. LMD belongs to [¯NRLL˜H1] where ˜Hiσ2H. LMD belongs to [¯LCLSLH2]. LMNS belongs to [¯NRSLφ]. [] implies that the concrete flavor structures are manifolded. Each of their structures are presented below.

    The scalar potential of our model is expressed as

    V=m2φφφ+m21H1H1+m22H2H2μ12(H1H2φ+h.c.)+λ1(H1H1)2+λ2(H2H2)2+λφ(φφ)2+λ3(H1H1)(H2H2)+λ4(H1H2)(H2H1)+λφH1(H1H1)(φφ)+λφH2(H2H2)(φφ),

    (2)

    where h.c. stands for Hermitian conjugate. We consider that φ develops a VEV at a significantly larger scale than H1,2. Then, after the φ developes a VEV, the scalar potential matches with that of two Higgs doublet model (THDM) without the (H1H2)2 term owing to the U(1)BL symmetry. In addition, the Yukawa couplings associated with the two Higgs doublet belong to the type-I THDM because only H1 can couple to the SM fermions. Hence, in our analysis below, we do not discuss the THDM part further, but focus on the neutrino mass.

    Before discussing the valid Lagrangians for the lepton sector, we define Yukawa couplings under the modular S4 symmetry as follows: the S4 doublets with kI=2 and 4 are denoted by Y(2)2[y1,y2]T and [y1,y2]T, respectively. The S4 triplet with kI=2 is Y(2)3[y3,y4,y5]T, while the S4 triplets with kI=4 are Y(4)3[y3,y4,y5]T and Y(4)3[y3. Each of these structures are explicitly presented in the Appendix.

    Charged lepton mass matrix:

    The renormalizable Lagrangian for the charged-lepton sector is expressed as

    \begin{aligned}[b]\\[-8pt] -{\cal{L}}_{M_\ell} = &\alpha_\ell \overline{e_{\rm R}} (y_3 L_{{\rm L}_e} +y_4 \bar L_{{\rm L}_\tau} +y_5 \bar L_{{\rm L}_\mu}) \tilde H_1 + \beta_\ell\left[\frac{\sqrt3}{2}\overline{\mu_{\rm R}}(y'_4 L_{{\rm L}_\mu} +y'_5 L_{{\rm L}_\tau}) +\overline{\tau_{\rm R}} \left(-y'_3 L_{{\rm L}_e} + \frac12 (y'_4 L_{{\rm L}_\tau} +y'_5 L_{{\rm L}_\mu}) \right) \right] \tilde H_1 \\ &+ \gamma_\ell\left[\frac{\sqrt3}{2}\overline{\tau_{\rm R}}(y''_4 L_{{\rm L}_\mu} +y''_5 L_{{\rm L}_\tau}) +\overline{\mu_{\rm R}} \left(y''_3 L_{{\rm L}_e} - \frac12 (y''_4 L_{{\rm L}_\tau} +y''_5 L_{{\rm L}_\mu}) \right) \right] \tilde H_1 +{\rm{h.c.}}, \end{aligned}

    (3)

    where \{\alpha_\ell,\beta_\ell,\gamma_\ell \} are real parameters without loss of generality. Subsequently, the charged-lepton mass matrix after the spontaneous symmetry breaking is expressed as

    \begin{aligned} (M_\ell)_{\rm RL} = \frac{v_1}{\sqrt2} \begin{pmatrix} \alpha_\ell y_3 & \alpha_\ell y_5 & \alpha_\ell y_4 \\ \gamma_\ell y''_3 & \dfrac{\sqrt3}{2}\beta_\ell y'_4-\dfrac12 \gamma_\ell y''_5 & \dfrac{\sqrt3}{2}\beta_\ell y'_5-\dfrac12 \gamma_\ell y''_4 \\ -\beta_\ell y'_3 & \dfrac{\sqrt3}{2}\gamma_\ell y''_4+\dfrac12 \beta_\ell y'_5 & \dfrac{\sqrt3}{2}\gamma_\ell y''_5+\dfrac12 \beta_\ell y'_4 \\ \end{pmatrix} . \end{aligned}

    (4)

    The charged-lepton mass eigenvalues are obtained by diagonalizing {\rm{diag}}[m_e,m_\mu,m_\tau] = V_{{\rm R}_\ell}^\dagger M_\ell V_{{\rm L}_\ell} , where V_{{\rm L}_\ell, {\rm R}_\ell} represent unitary matrices. In our numerical analysis, we will determine the free parameters \{\alpha_\ell,\beta_\ell,\gamma_\ell\} , to fit the three charged-lepton mass eigenstates after providing all the numerical values, by applying the following relationship. Here, we fix \alpha_\ell,\beta_\ell,\gamma_\ell , to determine the experimental three charged-lepton masses by applying the following relationship:

    \begin{aligned}[b] &{\rm{Tr}}[{M_\ell}^\dagger M_\ell] = |m_e|^2 + |m_\mu|^2 + |m_\tau|^2,\\& {\rm{Det}}[{M_\ell}^\dagger M_\ell] = |m_e|^2 |m_\mu|^2 |m_\tau|^2, \\ &({\rm{Tr}}[{M_\ell}^\dagger] M_\ell)^2 -{\rm{Tr}}[({M_\ell}^\dagger M_\ell)^2] \\=& 2( |m_e|^2 |m_\mu|^2 + |m_\mu|^2 |m_\tau|^2+ |m_e|^2 |m_\tau|^2 ). \end{aligned}

    (5)

    Neutral fermion mass matrices:

    First, we construct the valid Lagrangian of the Dirac mass matrix {\cal{L}}_{{M_{\rm D}}} . This is given by

    \begin{aligned}[b]\\[-8pt] {\cal{L}}_{M_{\rm D}} =& \alpha_{\rm D} \left[y'_3\left(\overline{N_{{\rm R}_3}} L_{{\rm L}_\mu}-\overline{N_{{\rm R}_2}} L_{{\rm L}_\tau}\right) + y'_4\left(\overline{N_{{\rm R}_1}} L_{{\rm L}_\tau}-\overline{N_{{\rm R}_3}} L_{{\rm L}_e}\right) + y'_5\left(\overline{N_{{\rm 2}}} L_{{\rm L}_e}-\overline{N_{{\rm R}_1}} L_{{\rm L}_\mu}\right)\right] \tilde H_1 \\ &+ \beta_{\rm D} \left[y''_3\left(\overline{N_{{\rm R}_3}} L_{{\rm L}_\tau}-\overline{N_{{\rm R}_2}} L_{{\rm L}_\mu}\right) + y''_4\left(-\overline{N_{{\rm R}_1}} L_{{\rm L}_\mu}-\overline{N_{{\rm R}_2}} L_{{\rm L}_e}\right) + y''_5\left(\overline{N_{{\rm R}_1}} L_{{\rm L}_\tau}+\overline{N_{{\rm R}_3}} L_{{\rm L}_e}\right)\right] \tilde H_1 \\ &+ \gamma_{\rm D} \left[\frac{\sqrt3}{2}y'_1(\overline{N_{{\rm R}_2}} L_{{\rm L}_\mu}+\overline{N_{{\rm R}_3}} L_{{\rm L}_\tau}) + y'_2 \{-\overline{N_{{\rm R}_1}} L_{{\rm L}_e} +\frac12 \left(\overline{N_{{\rm R}_2}} L_{{\rm L}_\tau}+\overline{N_{{\rm R}_3}} L_{{\rm L}_\mu}\right)\}\right] \tilde H_1+{\rm{h.c.}}, \end{aligned}

    (6)

    where we suppose \alpha_{\rm D} to be real while \beta_{\rm D},\ \gamma_{\rm D} are complex after rephasing the fields. Similar to the charged-lepton sector, we determine the Dirac mass matrix as follows:

    \begin{aligned} (M_{\rm D})_{\rm RL} = \frac{v_1}{\sqrt2} \begin{pmatrix} -\gamma_{\rm D} y'_2 & -\alpha_{\rm D} y'_5 - \beta_{\rm D} y''_4& \alpha_{\rm D} y'_4+\beta_{\rm D} y''_5 \\ \alpha_{\rm D} y'_5 - \beta_{\rm D} y''_4 &-\beta_{\rm D} y''_3+ \dfrac{\sqrt3}{2}\gamma_{\rm D} y'_1 &-\alpha_{\rm D} y'_3 + \dfrac12 \gamma_{\rm D} y'_2 \\ - \alpha_{\rm D} y'_4+\beta_{\rm D} y''_5 &\alpha_{\rm D} y'_3 + \dfrac12 \gamma_{\rm D} y'_2 & \beta_{\rm D} y''_3+ \dfrac{\sqrt3}{2}\gamma_{\rm D} y'_1 \\ \end{pmatrix} . \end{aligned}

    (7)

    For convenience, to analyze the neutrino oscillation, we redefine M_{\rm D}\equiv \dfrac{v_1}{\sqrt2}\tilde{M_{\rm D}} .

    Another Dirac Lagrangian is induced via {\cal{L}}_{M'_{\rm D}} , which is given by

    \begin{align} {\cal{L}}_{{M'_{\rm D}}} & = \alpha'_{\rm D} \left[y_3 \overline{L^{\rm C}_{{\rm L}_e}} +y_4 \overline{L^{\rm C}_{{\rm L}_\tau}}+y_5 \overline{L^{\rm C}_{{\rm L}_\mu}}\right]S_{{\rm L}_1} H_2 +{\rm{h.c.}}, \end{align}

    (8)

    Subsequently, we determine another Dirac mass matrix as follows:

    \begin{align} (M'_{\rm D})_{L_{\rm L} S_{\rm L}} = \frac{v_2}{\sqrt2} \alpha'_{\rm D} \begin{pmatrix} y_3 & 0 \\ y_5 & 0 \\ y_4 & 0 \\ \end{pmatrix} . \end{align}

    (9)

    Being the same as the reason for M_{\rm D} , we redefine M'_{\rm D}\equiv\dfrac{v_2}{\sqrt2} \alpha'_{\rm D} \tilde{M'_{\rm D}} .

    The third term of the Lagrangian {\cal{L}}_{M_{NS}} is given by

    \begin{aligned}[b] {\cal{L}}_{M_{NS}} = & \alpha_{NS} \left[y'_3 \overline{N_{{\rm R}_e}} + y'_4 \overline{N_{{\rm R}_\tau}} + y'_5\overline{N_{{\rm R}_\mu}}\right] S_{{\rm L}_1}\varphi \\&+ \beta_{NS} \left[y''_3 \overline{N_{{\rm R}_e}} + y''_4 \overline{N_{{\rm R}_\tau}} + y''_5\overline{N_{{\rm R}_\mu}}\right] S_{{\rm L}_2}\varphi +{\rm{h.c.}}, \end{aligned}

    (10)

    where \alpha_{NS},\; \beta_{NS} are real without loss of generality. Accordingly, we obtain the mass matrix

    \begin{aligned} M_{NS} = \frac{v_\varphi}{\sqrt2} \begin{pmatrix} y'_3 & y''_3 \\ y'_5 & y''_5 \\ y'_4 & y''_4 \\ \end{pmatrix} \begin{pmatrix} \alpha_{NS} & 0 \\ 0 & \beta_{NS} \\ \end{pmatrix} \equiv \frac{v_\varphi}{\sqrt2} \tilde M_{N S} . \end{aligned}

    (11)

    In basis of [\nu_{\rm L},N_{\rm R}^{\rm C},S_{\rm L}]^T , the neutral fermion mass matrix is given by

    \begin{align} M_N = \begin{pmatrix} 0_{3\times3} & M_{\rm D}^T & M'_{\rm D} \\ M_{\rm D} & 0_{3\times3} & M_{NS} \\ m'^T_{\rm D} & M_{NS}^T & 0_{2\times2} \end{pmatrix} . \end{align}

    (12)

    Hence, the active neutrino mass matrix is expressed as

    \begin{aligned}[b] m_\nu =& M'_{\rm D} (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D} \\&+ [M'_{\rm D} (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D}]^T \\ =& \frac{v_1v_2}{\sqrt{2} v_\varphi} \Big( \tilde M'_{\rm D} (\tilde M_{NS}^T \tilde M_{NS})^{-1}\tilde M_{NS}^T\tilde M_{\rm D}\\& + [\tilde M'_{\rm D} (\tilde M_{NS}^T\tilde M_{NS})^{-1}\tilde M_{NS}^T \tilde M_{\rm D}]^T \Big) = \kappa \tilde m_\nu, \end{aligned}

    (13)

    where \kappa\equiv \dfrac{v_1v_2}{\sqrt{2} v_\varphi} and we assume the mass hierarchies among M_{\rm D}, M'_{\rm D} \ll M_{NS} . Mass hierarchies is dynamically achieved in Refs. [100, 101]. The neutrino mass eigenvalues are obtained as follows: D_\nu = \kappa D_\nu = U_\nu^T m_\nu U_\nu = \kappa U_\nu^T \tilde m_\nu U_\nu , where U_\nu is a unitary matrix. Then, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix is given by U_{\rm PMNS}\equiv V^\dagger_{{\rm L}_\ell} U_\nu . Notice that m_\nu is rank two; hence, the lightest neutrino mass is zero. Here κ is described by one experimental value and dimensionless neutrino mass eigenstates as:

    \begin{align} {\rm{(NH)}}:\ \kappa^2 = \frac{|\Delta m_{\rm{atm}}^2|}{\tilde D_{\nu_3}^2}, \quad {\rm{(IH)}}:\ \kappa^2 = \frac{|\Delta m_{\rm{atm}}^2|}{\tilde D_{\nu_2}^2}, \end{align}

    (14)

    where \Delta m_{\rm{atm}}^2 denotes the atmospheric neutrino mass-squared difference and NH (IH) stands for normal (inverted) ordering. Subsequently, the solar mass difference squared can be expressed in terms of κ as follows:

    \begin{align} {\rm{(NH)}}:\ \Delta m_{\rm{sol}}^2 = {\kappa^2} {\tilde D_{\nu_2}^2}, \quad {\rm{(IH)}}:\ \Delta m_{\rm{sol}}^2 = {\kappa^2}({\tilde D_{\nu_2}^2-\tilde D_{\nu_1}^2}), \end{align}

    (15)

    which can be compared to the observed value. In other words, we explicitly express the mass eigenvalues in terms of \Delta m_{\rm{atm}}^2 and \Delta m_{\rm{sol}}^2 as:

    \begin{align} &{\rm{(NH)}}:\ D_{\nu_1}^2 = 0, \quad \ D_{\nu_2}^2 = \Delta m_{\rm{sol}}^2, \quad \ D_{\nu_3}^2 = \Delta m_{\rm{atm}}^2, \end{align}

    (16)

    \begin{align} &{\rm{(IH)}}:\ D_{\nu_1}^2 = \Delta m_{\rm{atm}}^2 - \Delta m_{\rm{sol}}^2, \;\;\; \ D_{\nu_2}^2 = \Delta m_{\rm{sol}}^2,\;\;\; \ D_{\nu_3}^2 = 0, \end{align}

    (17)

    which implies that NH is hierarchical, but IH is degenerate, as \Delta m_{\rm{sol}}^2/ \Delta m_{\rm{atm}}^2<<1 . Here, we expand |\tilde m_\nu|^2 in terms of q\equiv {\rm e}^{2p\pi {\rm i}\tau}<<1 . Then, the mass matrix is given by

    \begin{aligned} |\tilde m_\nu|^2\sim \begin{pmatrix} {\cal{O}}(1) & {\cal{O}}(q) & 0 \\ {\cal{O}}(q) & {\cal{O}}(1) & {\cal{O}}(q^2) \\ 0 & {\cal{O}}(q^2) & 0 \end{pmatrix}. \end{aligned}

    (18)

    The ratio between two nonzero squared eigenvalues R is estimated by

    \begin{align} R = {\cal{O}}(q) <<1 . \end{align}

    (19)

    This suggests that the neutrino mass eigenvalues tend to be hierarchical; hence, NH is favored. In fact, we would not obtain the allowed region within 3\sigma for IH in our numerical analysis. Therefore, we focus on NH hereafter.

    In our model, the PMNS matrix is parametrized by three mixing angles \theta_{ij} (i,j = 1,2,3; i < j) , one CP violating Dirac phase \delta_{\rm CP} , and one Majorana phase \alpha_{21} as follows:

    \begin{equation} U_{\rm PMNS} = \begin{pmatrix} c_{12} c_{13} & s_{12} c_{13} & s_{13} {\rm e}^{-{\rm i} \delta_{\rm CP}} \\ -s_{12} c_{23} - c_{12} s_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & c_{12} c_{23} - s_{12} s_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & s_{23} c_{13} \\ s_{12} s_{23} - c_{12} c_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & -c_{12} s_{23} - s_{12} c_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & c_{23} c_{13} \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 \\ 0 & {\rm e}^{{\rm i} \frac{\alpha_{21}}{2}} & 0 \\ 0 & 0 & 1 \end{pmatrix}, \end{equation}

    (20)

    where c_{ij} and s_{ij} represent \cos \theta_{ij} and \sin \theta_{ij} , respectively. Then, these mixings are given in terms of the components of U_{\rm PMNS} as follows:

    \begin{aligned}[b]& \sin^2\theta_{13} = |(U_{\rm PMNS})_{13}|^2,\quad \sin^2\theta_{23} = \frac{|(U_{\rm PMNS})_{23}|^2}{1-|(U_{\rm PMNS})_{13}|^2},\\& \sin^2\theta_{12} = \frac{|(U_{\rm PMNS})_{12}|^2}{1-|(U_{\rm PMNS})_{13}|^2}. \end{aligned}

    (21)

    Furthermore, we compute the Jarlskog invariant J_{\rm CP} derived from PMNS matrix elements as follows:

    \begin{equation} J_{\rm CP} = {\rm{Im}} [U_{e1} U_{\mu 2} U_{e 2}^* U_{\mu 1}^*] = s_{23} c_{23} s_{12} c_{12} s_{13} c^2_{13} \sin \delta_{\rm CP}. \end{equation}

    (22)

    Majorana phase is estimated in terms of other invariant I_1 as follows:

    \begin{equation} I_1 = {\rm{Im}}[U^*_{e1} U_{e2}] = c_{12} s_{12} c_{13}^2 \sin \left( \frac{\alpha_{21}}{2} \right). \end{equation}

    (23)

    In addition, the effective mass for the neutrinoless double beta decay is expressed as:

    \begin{align} &{\rm{(NH)}}:\ \langle m_{ee}\rangle = \kappa|\tilde D_{\nu_2} s^2_{12} c^2_{13}{\rm e}^{{\rm i}\alpha_{21}}+\tilde D_{\nu_3} s^2_{13}{\rm e}^{-2{\rm i}\delta_{CP}}|, \end{align}

    (24)

    \begin{align} &{\rm{(IH)}}:\ \langle m_{ee}\rangle = \kappa|\tilde D_{\nu_1} c^2_{12} c^2_{13}+\tilde D_{\nu_2} s^2_{12} c^2_{13}{\rm e}^{{\rm i}\alpha_{21}}|, \end{align}

    (25)

    where its value could be measured by KamLAND-Zen in the future [102]. We will adopt the neutrino experimental data at the 3σ interval in Nufit 5.0 [103, 104].

    Non-unitarity:

    Here, let us briefly discuss the non-unitarity matrix U'_{\rm PMNS} . This is typically parametrized by the form

    \begin{align} U'_{\rm PMNS}\equiv \left(1-\frac12 F^\dagger F\right) U_{\rm PMNS}, \end{align}

    (26)

    where F\equiv (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D} is a hermitian matrix and U'_{\rm PMNS} represents the deviation from the unitarity. Applying global constraints [105], one finds [106]

    \begin{align} |FF^\dagger|\le \left[\begin{array}{ccc} 2.5\times 10^{-3} & 2.4\times 10^{-5} & 2.7\times 10^{-3} \\ 2.4\times 10^{-5} & 4.0\times 10^{-4} & 1.2\times 10^{-3} \\ 2.7\times 10^{-3} & 1.2\times 10^{-3} & 5.6\times 10^{-3} \\ \end{array}\right]. \end{align}

    (27)

    In our case, F\equiv (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D} = \dfrac{v_1}{v_\varphi} (\tilde M_{NS}^T\tilde M_{NS})^{-1}\times \tilde M_{NS}^T \tilde M_{\rm D} . Because v_\varphi is freely considered to be large (while v_1 = {\cal{O}}(100) GeV at most), we easily satisfy this bound.

    In this section, we perform a numerical \Delta \chi^2 analysis to determine the parameters that satisfy the neutrino oscillation data and non-unitarity constraint; in addition, we present our predictions, where we adopt best fit values of charged-lepton masses. Here, we focus on NH because IH is disfavored by the analytical estimation, as can be observed in the previous section.

    In our numerical analysis, we randomly scan free parameters in following ranges

    \begin{align} & \{ \alpha_{\rm D}, |\beta_{\rm D}|, |\gamma_{\rm D}|, \alpha_{NS}, \beta_{NS} \} \in [10^{-5},1.0], \quad t_\beta\in[10,100], \end{align}

    (28)

    where τ runs over the fundamental region, t_\beta\equiv v_1/v_2 , and \sqrt{v_1^2+v_2^2} = 246 GeV. We perform numerical analysis under the above regions. Figure 1 illustrates the correlation between the real and imaginary parts of τ, where the blue, green, and red points are allowed within 2, 3, and 5 of the \Delta \chi^2 analysis, respectively, for five accurately known observables \Delta m^2_{\rm{atm}},\Delta m^2_{\rm{sol}},s_{12}^2,s_{23}^2,s_{13}^2 in Nufit 5.0 [103, 104]. The real part runs through the entire range, but the imaginary part is localized at the [1.35-1.5] .

    Figure 1

    Figure 1.  (color online) Allowed region of modulus τ, where the blue, green, and red points are allowed within 2, 3, and 5 of the \Delta \chi^2 analysis, respectively.

    Figure 2 illustrates the correlation between the Majorana phase \alpha_{21} and Dirac CP phase \delta_{\rm CP} . The legend is the same as in the case of Fig. 1. We identify a clear feature in the region where \alpha_{21} is within \{0^\circ - 120^\circ, 240^\circ - 310^\circ \} , while \delta_{\rm CP} is within \{70^\circ - 290^\circ \} .

    Figure 2

    Figure 2.  (color online) Correlation between a Majorana phase \alpha_{21} and \delta_{\rm CP} , where the legend is the same as the case of Fig. 1.

    The upper figures in Fig. 3 present the correlation between the effective mass for the neutrinoless double beta decay \langle m_{ee} \rangle and the sum of neutrino masses \sum m_i in the [eV] unit, where the legend is the same as the case of Fig. 1. We determine that \langle m_{ee} \rangle is allowed within [0.001-0.004] eV. However, \sum m_i is restricted to be approximately 0.06 eV, which would be a sharp prediction of this model.

    Figure 3

    Figure 3.  (color online) Predicted correlation between the effective mass for the neutrinoless double beta decay \langle m_{ee} \rangle and the sum of neutrino masses \sum m_i in [eV], where the legend is the same as the case of Fig. 1.

    Finally, we present a benchmark in Table 2, which is selected, such that \sqrt{\Delta \chi^2} is minimum. The mass matrices for the dimensionless neutrino and charged-lepton are determined as

    Table 2

    Table 2.  Benchmark point of our input parameters and observables, which is selected, such that \sqrt{\Delta \chi^2} is minimum.
    parameter valve
    τ 0.113762 + 1.43906 {\rm i}
    t_\beta 98.6
    [\alpha_\ell, \gamma_{\ell}, \beta_{\ell}] [1.72\times 10^{-5},\ 6.15\times 10^{-4},\ 8.06\times 10^{-4}]
    [\alpha_{\rm D}, \alpha_{NS}, \beta_{NS}] [-0.0112907,\ -0.00078512,\ 0.0203379]
    [\beta_{\rm D}, \gamma_{\rm D}] [-0.01432- 0.00360 {\rm i},-1.90\times 10^{-5} - 2.34\times10^{-6} {\rm i}]
    \Delta m^2_{\rm{atm}} 2.53\times10^{-3} {\rm{eV}}^2
    \Delta m^2_{\rm{sol}} 7.48\times10^{-5} {\rm{eV}}^2
    \sin^2\theta_{12} 0.289
    \sin^2\theta_{23} 0.565
    \sin^2\theta_{13} 0.02207
    [\delta_{\rm CP},\ \alpha_{21}] [248^\circ,\, 80.3^\circ]
    \sum m_i 58.5 meV
    \langle m_{ee} \rangle 1.65 meV
    \sqrt{\Delta\chi^2} 1.40
    DownLoad: CSV
    Show Table

    \begin{align} \tilde m_\nu& = \left[\begin{array}{ccc} 37.1331- 17.9304 {\rm i} & -16.2576 + 55.2032 {\rm i} & 35.2235+ 24.3462 {\rm i} \\ -16.2576 + 55.2032 {\rm i} &-40.0223 -38.1947 {\rm i} & -55.9229 + 44.0753 {\rm i} \\ 35.2235+ 24.3462 {\rm i} & -55.9229 + 44.0753 {\rm i} &-1.23231 + 74.8642 {\rm i} \\ \end{array}\right], \end{align}

    (29)

    \begin{align} M_\ell& = \left[\begin{array}{ccc} 0.00232 + 0.000467 {\rm i} & -0.00111 + 0.00314 {\rm i} & -0.00453 + 0.00274 {\rm i} \\ -0.282 + 0.102 {\rm i} & -0.0524 - 0.469 {\rm i}& 0.214 - 0.559 {\rm i} \\ 0.336 - 0.277 {\rm i} & 0.362 + 0.706 {\rm i}& -0.160 + 1.23 {\rm i} \\ \end{array}\right]. \end{align}

    (30)

    We studied a linear seesaw model with a field content as minimum as possible, introducing a modular S_4 using U(1)_{B-L} symmetries. Owing to the rank two neutrino mass matrix, we obtained a vanishing neutrino mass eigenvalue. Furthermore, only the normal mass hierarchy of neutrinos was favored via the modular S_4 symmetry. In our numerical \Delta \chi^2 analysis, we determined a rather sharp prediction on the sum of neutrino masses to be approximately 60 meV. The imaginary part of τ was restricted at 1.35-1.5, while the real part ran through the entire range in the fundamental region. Other remarks are presented below:

    1. \alpha_{21} was within \{0^\circ - 120^\circ, 240^\circ - 310^\circ \} , while \delta_{\rm CP} was within \{70^\circ - 290^\circ \} .

    2. \langle m_{ee} \rangle was allowed by \{0.001 - 0.004\} eV.

    Therefore our model indicates several predictions in the neutrino sector that was triggered by the minimal structure with the S_4 modular symmetry.

    This research was supported by the Korean Local Governments - Gyeongsangbuk-do Province and Pohang City (H.O.). H. O. is sincerely grateful to the KIAS member.

    Here, we review some properties of the modular S_4 symmetry. In general, the modular group \bar\Gamma is the group of linear fractional transformation γ acting on the modulus τ, which belongs to the upper-half complex plane and transforms as

    \begin{equation} \tau \longrightarrow \gamma\tau = \frac{a\tau + b}{c \tau + d}, \end{equation}\tag{A1}

    {\rm{where}}\; \; a,b,c,d \in \mathbb{Z}\; \; {\rm{and}}\; \; ad-bc = 1, \; \; {\rm{Im}} [\tau]>0 \; .

    This is isomorphic to PSL(2,\mathbb{Z}) = SL(2,\mathbb{Z})/\{I,-I\} transformation. Then modular transformation is generated by two transformations S and T defined as:

    \begin{eqnarray} S:\tau \longrightarrow -\frac{1}{\tau}\ , \qquad T:\tau \longrightarrow \tau + 1\ , \end{eqnarray} \tag{A2}

    and they satisfy the following algebraic relations,

    \begin{equation} S^2 = \mathbb{I}\ , \qquad (ST)^3 = \mathbb{I}\ . \end{equation} \tag{A3}

    Here, we introduce the series of groups \Gamma(N)\; (N = 1,2,3,\dots) which are defined by

    \begin{align} \begin{aligned} \Gamma(N) = \left \{ \begin{pmatrix} a & b \\ c & d \end{pmatrix} \in SL(2,\mathbb{Z})\; , \; \; \begin{pmatrix} a & b \\ c & d \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} \; \; ({\rm{mod}} N) \right \} \end{aligned}, \end{align}\tag{A4}

    and we define \bar\Gamma(2)\equiv \Gamma(2)/\{I,-I\} for N = 2 . Because the element -I does not belong to the \Gamma(N) for N>2 case, we have \bar\Gamma(N) = \Gamma(N) , which are the infinite normal subgroup of \bar \Gamma known as principal congruence subgroups. Hence, we obtain finite modular groups as the quotient groups defined by \Gamma_N\equiv \bar \Gamma/\bar \Gamma(N) . For these finite groups \Gamma_N , T^N = \mathbb{I} is imposed, and the groups \Gamma_N with N = 2,3,4,5 are isomorphic to S_3 , A_4 , S_4 and A_5 , respectively [11].

    Modular forms of level N are holomorphic functions f(\tau) transformed under the action of \Gamma(N) , which is expressed as:

    \begin{equation} f(\gamma\tau) = (c\tau+d)^k f(\tau)\; , \; \; \gamma \in \Gamma(N)\; , \end{equation}\tag{A5}

    where k denotes the reputed modular weight.

    Here, we discuss the modular symmetric theory framework without imposing supersymmetry explicitly, considering the S_4 ( N = 4 ) modular group. Under the modular transformation in Eq. (A1), a field \phi^{(I)} is also transformed as

    \begin{equation} \phi^{(I)} \to (c\tau+d)^{-k_I}\rho^{(I)}(\gamma)\phi^{(I)}, \end{equation} \tag{A6}

    where -k_I represents the modular weight and \rho^{(I)}(\gamma) denotes a unitary representation matrix of \gamma\in\Gamma(4) . Hence, a Lagrangian such as the Yukawa terms can be invariant if the sum of modular weight from fields and the modular form in the corresponding term is zero (also invariant under S_4 and gauge symmetry).

    The kinetic and quadratic terms of the scalar fields can be expressed as:

    \begin{equation} \sum\limits_I\frac{|\partial_\mu\phi^{(I)}|^2}{(-i\tau+i\bar{\tau})^{k_I}} \; , \quad \sum\limits_I\frac{|\phi^{(I)}|^2}{(-i\tau+i\bar{\tau})^{k_I}} \; , \end{equation}\tag{A7}

    which is invariant under the modular transformation, and the overall factor is eventually absorbed by a field redefinition consistently. Therefore, the Lagrangian associated with these terms should be invariant under the modular symmetry.

    The basis of modular forms with weight 2, Y^{(2)}_{\rm{2}} = (y_{1},y_{2}) and Y^{(2)}_{\rm{3'}} = (y_{3},y_{4},y_{5}) , which are transformed as a doublet and triplet under S_4 , are determined in terms of the Dedekind eta-function \eta(\tau) and its derivative [10]:

    \begin{aligned}[b] y_1(\tau) = & \frac{\rm i}{8} \left( 8 \frac{\eta'\left(\tau+\dfrac12\right)}{\eta\left(\tau+\dfrac12\right)} + 32 \frac{\eta'(4\tau)}{\eta(4\tau)} - \frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)}- \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)}\right. \\&\left.- \frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} - \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right), \\ y_2(\tau) = & \frac{ {\rm i} \sqrt3}{8} \left(\frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)} - \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)} + \frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} - \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right), \\ y_3(\tau) = & {\rm i} \left(\frac{\eta'\left(\tau+\dfrac12\right)}{\eta\left(\tau+\dfrac12\right)} - 4 \frac{\eta'(4\tau)}{\eta(4\tau)} \right), \\ y_4(\tau) = & \frac{\rm i}{4\sqrt2} \left( - \frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)}+{\rm i} \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)} +\frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} - {\rm i} \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right), \end{aligned}

    \begin{aligned}[b] y_5(\tau) = \frac{\rm i}{4\sqrt2} \left( - \frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)} - i \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)} + \frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} +{\rm i} \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right) . \end{aligned} \tag{A8}

    y_i 's can be expanded in terms of q as follows:

    \begin{align} y_1 & = - 3 \pi \left(\frac{b_1}{8} + 3 b_5\right),\; \ y_2 = 3 \sqrt{3} \pi b_3,\; \ y_3 = - \pi \left( -\frac{b_1}{4} + 2 b_5\right), \\ y_4 & = - \pi \sqrt{2} b_2, \;\ y_5 = - 4\pi \sqrt{2} b_4, \end{align}\tag{A9}

    where b_i are given by

    \begin{align} b_1 \sim 1, \;\;\ b_2 \sim q,\;\; \ b_3 \sim q^2,\; \;\ b_4 \sim 0,\; \;\ b_5 \sim 0, \end{align} \tag{A10}

    with q = \exp(2 \pi i \tau) and |q| \ll 1 [96].

    Subsequently, Yukawas with higher weights are constructed by the multiplication rules of S_4 , and the following couplings can be obtained:

    \begin{aligned}[b] &Y^{(4)}_{\bf2} = \left[\begin{array}{c} y_2^2-y_1^2 \\ 2 y_1 y_2 \\ \end{array}\right], \;\;\; Y^{(4)}_{\bf3} = \left[\begin{array}{c} -2y_2y_3 \\ \sqrt3 y_1 y_5+ y_2 y_4 \\ \sqrt3 y_1 y_4+ y_2 y_5 \\ \end{array}\right],\\& Y^{(4)}_{\bf3'} = \left[\begin{array}{c} 2y_1y_3 \\ \sqrt3 y_2 y_5 - y_1 y_4 \\ \sqrt3 y_2 y_4- y_1 y_5 \\ \end{array}\right] . \end{aligned} \tag{A11}

    Inverse seesaw model [5, 6] also requests both of \begin{document}$ N_{\rm R} $\end{document} and \begin{document}$ S_{\rm L} $\end{document} and the neutrino mass could be realized within TeV scale. But this model may not require GUT scale.

    Notice here that our model has a different mechanism from the one of original linear seesaw, even though our neutrino model can be realized within TeV scale. Therefore, our neutrino mass matrix directly depends on the scale of \begin{document}$ U(1)_{B-L} $\end{document} breaking and ratio between vacuum expectation values (VEVs) of two Higgs doublet model as can be seen in the main text, while we expect higher energy scale such as GUT in order that the modular field must break and get VEV denoted by \begin{document}$ \tau $\end{document} in our literature.

    For interested readers, we provide some literature reviews, which are useful to understand the non-Abelian group and its applications to flavor structure [86-93].

    Mass hierarchies is dynamically achieved in Refs. [100, 101].

    [1] P. Minkowski, Phys. Lett. B 67, 421 (1977) doi: 10.1016/0370-2693(77)90435-X
    [2] T. Yanagida, in Proceedings of the Workshop on the Unified Theory and the Baryon Number in the Universe (O. Sawada and A. Sugamoto, eds.), KEK, Tsukuba, Japan, 1979, p. 95
    [3] M. Gell-Mann, P. Ramond, and R. Slansky, Supergravity (P. van Nieuwenhuizen et al. eds.), North Holland, Amsterdam, 1979, p. 315; S. L. Glashow, The future of elementary particle physics, in Proceedings of the 1979 Cargèse Summer Institute on Quarks and Leptons (M. Levy et al. eds.), Plenum Press, New York, 1980, p. 687
    [4] R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980) doi: 10.1103/PhysRevLett.44.912
    [5] R. N. Mohapatra and J. W. F. ValleSenjanovic, Phys. Rev. D 34, 1642 (1986)
    [6] D. Wyler and L. Wolfenstein, Nucl. Phys. B 218, 205 (1983) doi: 10.1016/0550-3213(83)90482-0
    [7] E. K. Akhmedov, M. Lindner, E. Schnapka et al., Phys. Lett. B 368, 270 (1996), arXiv:[hep-ph/9507275 doi: 10.1016/0370-2693(95)01504-3
    [8] E. K. Akhmedov, M. Lindner, E. Schnapka et al., Phys. Rev. D 53, 2752 (1996), arXiv:[hep-ph/9509255
    [9] M. Malinsky, J. C. Romao, and J. W. F. Valle, Phys. Rev. Lett. 95, 161801 (2005), arXiv:hep-ph/0506296 [hep-ph doi: 10.1103/PhysRevLett.95.161801
    [10] F. Feruglio, doi: 10.1142/9789813238053_0012 arXiv: 1706.08749[hep-ph]
    [11] R. de Adelhart Toorop, F. Feruglio, and C. Hagedorn, Nucl. Phys. B 858, 437-467 (2012), arXiv:1112.1340[hep-ph doi: 10.1016/j.nuclphysb.2012.01.017
    [12] J. C. Criado and F. Feruglio, SciPost Phys. 5(5), 042 (2018), arXiv:1807.01125[hep-ph doi: 10.21468/SciPostPhys.5.5.042
    [13] T. Kobayashi, N. Omoto, Y. Shimizu et al., JHEP 11, 196 (2018), arXiv:1808.03012[hep-ph doi: 10.1007/JHEP11(2018)196
    [14] H. Okada and M. Tanimoto, Phys. Lett. B 791, 54-61 (2019), arXiv:1812.09677[hep-ph doi: 10.1016/j.physletb.2019.02.028
    [15] T. Nomura and H. Okada, Phys. Lett. B 797, 134799 (2019), arXiv:1904.03937[hep-ph doi: 10.1016/j.physletb.2019.134799
    [16] H. Okada and M. Tanimoto, Eur. Phys. J. C 81(1), 52 (2021), arXiv:1905.13421[hep-ph doi: 10.1140/epjc/s10052-021-08845-y
    [17] F. J. de Anda, S. F. King, and E. Perdomo, Phys. Rev. D 101(1), 015028 (2020), arXiv:1812.05620[hep-ph doi: 10.1103/PhysRevD.101.015028
    [18] P. P. Novichkov, S. T. Petcov, and M. Tanimoto, Phys. Lett. B 793, 247-258 (2019), arXiv:1812.11289[hep-ph doi: 10.1016/j.physletb.2019.04.043
    [19] T. Nomura and H. Okada, Nucl. Phys. B 966, 115372 (2021), arXiv:1906.03927[hep-ph doi: 10.1016/j.nuclphysb.2021.115372
    [20] H. Okada and Y. Orikasa, arXiv: 1907.13520 [hep-ph]
    [21] G. J. Ding, S. F. King, and X. G. Liu, JHEP 09, 074 (2019), arXiv:1907.11714[hep-ph doi: 10.1007/JHEP09(2019)074
    [22] T. Nomura, H. Okada, and O. Popov, Phys. Lett. B 803, 135294 (2020), arXiv:1908.07457[hep-ph doi: 10.1016/j.physletb.2020.135294
    [23] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Rev. D 100(11), 115045 (2019) [Erratum: Phys. Rev. D 101(3), 039904 (2020)] arXiv: 1909.05139[hep-ph]
    [24] T. Asaka, Y. Heo, T. H. Tatsuishi a et al., JHEP 01, 144 (2020), arXiv:1909.06520[hep-ph doi: 10.1007/JHEP01(2020)144
    [25] D. Zhang, Nucl. Phys. B 952, 114935 (2020), arXiv:1910.07869[hep-ph doi: 10.1016/j.nuclphysb.2020.114935
    [26] G. J. Ding, S. F. King, X. G. Liu et al., JHEP 12, 030 (2019), arXiv:1910.03460[hep-ph doi: 10.1007/JHEP12(2019)030
    [27] T. Kobayashi, T. Nomura, and T. Shimomura, Phys. Rev. D 102(3), 035019 (2020), arXiv:1912.00637[hep-ph doi: 10.1103/PhysRevD.102.035019
    [28] T. Nomura, H. Okada, and S. Patra, Nucl. Phys. B 967, 115395 (2021), arXiv:1912.00379[hep-ph doi: 10.1016/j.nuclphysb.2021.115395
    [29] X. Wang, Nucl. Phys. B 957, 115105 (2020), arXiv:1912.13284[hep-ph doi: 10.1016/j.nuclphysb.2020.115105
    [30] H. Okada and Y. Shoji, Nucl. Phys. B 961, 115216 (2020), arXiv:2003.13219[hep-ph doi: 10.1016/j.nuclphysb.2020.115216
    [31] H. Okada and M. Tanimoto, arXiv: 2005.00775 [hep-ph]
    [32] M. K. Behera, S. Singirala, S. Mishra et al., arXiv: 2009.01806 [hep-ph]
    [33] M. K. Behera, S. Mishra, S. Singirala et al., arXiv: 2007.00545 [hep-ph]
    [34] T. Nomura and H. Okada, arXiv: 2007.04801 [hep-ph]
    [35] T. Nomura and H. Okada, arXiv: 2007.15459 [hep-ph]
    [36] T. Asaka, Y. Heo, and T. Yoshida, Phys. Lett. B 811, 135956 (2020), arXiv:2009.12120[hep-ph doi: 10.1016/j.physletb.2020.135956
    [37] H. Okada and M. Tanimoto, Phys. Rev. D 103(1), 015005 (2021), arXiv:2009.14242[hep-ph doi: 10.1103/PhysRevD.103.015005
    [38] K. I. Nagao and H. Okada, arXiv: 2010.03348 [hep-ph]
    [39] H. Okada and M. Tanimoto, JHEP 03, 010 (2021), arXiv:2012.01688[hep-ph doi: 10.1007/JHEP03(2021)010
    [40] C. Y. Yao, J. N. Lu, and G. J. Ding, JHEP 05, 102 (2021), arXiv:2012.13390 [hep-ph doi: 10.1007/JHEP05(2021)102
    [41] P. Chen, G. J. Ding, and S. F. King, JHEP 04, 239 (2021), arXiv:2101.12724 [hep-ph doi: 10.1007/JHEP04(2021)239
    [42] M. Kashav and S. Verma, arXiv: 2103.07207 [hep-ph]
    [43] H. Okada, Y. Shimizu, M. Tanimoto et al., arXiv: 2105.14292 [hep-ph]
    [44] I. de Medeiros Varzielas, and J. Lourenço, arXiv: 2107.04042 [hep-ph]
    [45] T. Nomura, H. Okada, and Y. Orikasa, arXiv: 2106.12375 [hep-ph]
    [46] P. T. P. Hutauruk, D. W. Kang, J. Kim et al., arXiv: 2012.11156 [hep-ph]
    [47] G. J. Ding, S. F. King, and J. N. Lu, arXiv: 2108.09655 [hep-ph]
    [48] K. I. Nagao and H. Okada, arXiv: 2108.09984 [hep-ph]
    [49] T. Kobayashi, K. Tanaka, and T. H. Tatsuishi, Phys. Rev. D 98(1), 016004 (2018), arXiv:1803.10391[hep-ph doi: 10.1103/PhysRevD.98.016004
    [50] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Lett. B 794, 114-121 (2019), arXiv:1812.11072[hep-ph doi: 10.1016/j.physletb.2019.05.034
    [51] T. Kobayashi, Y. Shimizu, K. Takagi et al., PTEP 2020(5), 053B05 (2020), arXiv:1906.10341[hep-ph doi: 10.1093/ptep/ptaa055
    [52] H. Okada and Y. Orikasa, Phys. Rev. D 100(11), 115037 (2019), arXiv:1907.04716[hep-ph doi: 10.1103/PhysRevD.100.115037
    [53] S. Mishra, arXiv: 2008.02095 [hep-ph]
    [54] X. Du and F. Wang, JHEP 02, 221 (2021), arXiv:2012.01397[hep-ph doi: 10.1007/JHEP02(2021)221
    [55] J. T. Penedo and S. T. Petcov, Nucl. Phys. B 939, 292-307 (2019), arXiv:1806.11040[hep-ph doi: 10.1016/j.nuclphysb.2018.12.016
    [56] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 04, 005 (2019), arXiv:1811.04933[hep-ph doi: 10.1007/JHEP04(2019)005
    [57] T. Kobayashi, Y. Shimizu, K. Takagi et al., JHEP 02, 097 (2020), arXiv:1907.09141[hep-ph doi: 10.1007/JHEP02(2020)097
    [58] S. F. King and Y. L. Zhou, Phys. Rev. D 101(1), 015001 (2020), arXiv:1908.02770[hep-ph doi: 10.1103/PhysRevD.101.015001
    [59] H. Okada and Y. Orikasa, arXiv: 1908.08409 [hep-ph]
    [60] J. C. Criado, F. Feruglio, and S. J. D. King, JHEP 02, 001 (2020), arXiv:1908.11867[hep-ph doi: 10.1007/JHEP02(2020)001
    [61] X. Wang and S. Zhou, JHEP 05, 017 (2020), arXiv:1910.09473[hep-ph doi: 10.1007/JHEP05(2020)017
    [62] Y. Zhao and H. H. Zhang, JHEP 03, 002 (2021), arXiv:2101.02266 [hep-ph doi: 10.1007/JHEP03(2021)002
    [63] S. F. King and Y. L. Zhou, JHEP 04, 291 (2021), arXiv:2103.02633 [hep-ph doi: 10.1007/JHEP04(2021)291
    [64] G. J. Ding, S. F. King, and C. Y. Yao, arXiv: 2103.16311 [hep-ph]
    [65] X. Zhang and S. Zhou, arXiv: 2106.03433 [hep-ph]
    [66] Bu-Yao Qu, Xiang-Gan Liu, Ping-Tao Chen et al., arXiv: 2106.11659 [hep-ph]
    [67] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 04, 174 (2019), arXiv:1812.02158[hep-ph doi: 10.1007/JHEP04(2019)174
    [68] G. J. Ding, S. F. King, and X. G. Liu, Phys. Rev. D 100(11), 115005 (2019), arXiv:1903.12588[hep-ph doi: 10.1103/PhysRevD.100.115005
    [69] X. Wang, B. Yu, and S. Zhou, Phys. Rev. D 103(7), 076005 (2021), arXiv:2010.10159[hep-ph doi: 10.1103/PhysRevD.103.076005
    [70] C. Y. Yao, X. G. Liu, and G. J. Ding, Phys. Rev. D 103(9), 095013 (2021), arXiv:2011.03501[hep-ph doi: 10.1103/PhysRevD.103.095013
    [71] X. Wang and S. Zhou, arXiv: 2102.04358 [hep-ph]
    [72] M. K. Behera and R. Mohanta, arXiv: 2108.01059 [hep-ph]
    [73] A. Baur, H. P. Nilles, A. Trautner et al., Phys. Lett. B 795, 7-14 (2019), arXiv:1901.03251[hep-th doi: 10.1016/j.physletb.2019.03.066
    [74] I. de Medeiros Varzielas, S. F. King et al., Phys. Rev. D 101(5), 055033 (2020), arXiv:1906.02208[hep-ph doi: 10.1103/PhysRevD.101.055033
    [75] X. G. Liu and G. J. Ding, JHEP 08, 134 (2019), arXiv:1907.01488[hep-ph doi: 10.1007/JHEP08(2019)134
    [76] P. Chen, G. J. Ding, J. N. Lu et al., Phys. Rev. D 102(9), 095014 (2020), arXiv:2003.02734[hep-ph doi: 10.1103/PhysRevD.102.095014
    [77] C. C. Li, X. G. Liu, and G. J. Ding, arXiv: 2108.02181 [hep-ph]
    [78] P. P. Novichkov, J. T. Penedo, and S. T. Petcov, Nucl. Phys. B 963, 115301 (2021), arXiv:2006.03058[hep-ph doi: 10.1016/j.nuclphysb.2020.115301
    [79] X. G. Liu, C. Y. Yao, and G. J. Ding, Phys. Rev. D 103(5), 056013 (2021), arXiv:2006.10722[hep-ph doi: 10.1103/PhysRevD.103.056013
    [80] S. Kikuchi, T. Kobayashi, H. Otsuka et al., JHEP 11, 101 (2020), arXiv:2007.06188[hep-th doi: 10.1007/JHEP11(2020)101
    [81] Y. Almumin, M. C. Chen, V. Knapp-Pérez et al., JHEP 05, 078 (2021), arXiv:2102.11286 [hep-th doi: 10.1007/JHEP05(2021)078
    [82] G. J. Ding, F. Feruglio, and X. G. Liu, SciPost Phys. 10, 133 (2021), arXiv:2102.06716 [hep-ph doi: 10.21468/SciPostPhys.10.6.133
    [83] F. Feruglio, V. Gherardi, A. Romanino et al., JHEP 05, 242 (2021), arXiv:2101.08718 [hep-ph doi: 10.1007/JHEP05(2021)242
    [84] S. Kikuchi, T. Kobayashi, and H. Uchida, arXiv: 2101.00826 [hep-th]
    [85] P. P. Novichkov, J. T. Penedo, and S. T. Petcov, JHEP 04, 206 (2021), arXiv:2102.07488 [hep-ph doi: 10.1007/JHEP04(2021)206
    [86] G. Altarelli and F. Feruglio, Rev. Mod. Phys. 82, 2701-2729 (2010), arXiv:1002.0211[hep-ph doi: 10.1103/RevModPhys.82.2701
    [87] H. Ishimori, T. Kobayashi, H. Ohki et al., Prog. Theor. Phys. Suppl. 183, 1-163 (2010), arXiv:1003.3552[hep-th doi: 10.1143/PTPS.183.1
    [88] H. Ishimori, T. Kobayashi, H. Ohki et al., Notes Phys. 858, 1-227 (2012) doi: 10.1007/978-3-642-30805-5
    [89] D. Hernandez and A. Y. Smirnov, Phys. Rev. D 86, 053014 (2012), arXiv:1204.0445[hep-ph doi: 10.1103/PhysRevD.86.053014
    [90] S. F. King and C. Luhn, Rept. Prog. Phys. 76, 056201 (2013), arXiv:1301.1340[hep-ph doi: 10.1088/0034-4885/76/5/056201
    [91] S. F. King, A. Merle, S. Morisi et al., New J. Phys. 16, 045018 (2014), arXiv:1402.4271[hep-ph doi: 10.1088/1367-2630/16/4/045018
    [92] S. F. King, Prog. Part. Nucl. Phys. 94, 217-256 (2017), arXiv:1701.04413[hep-ph doi: 10.1016/j.ppnp.2017.01.003
    [93] S. T. Petcov, Eur. Phys. J. C 78(9), 709 (2018), arXiv:1711.10806[hep-ph doi: 10.1140/epjc/s10052-018-6158-5
    [94] A. Baur, H. P. Nilles, A. Trautner et al., Nucl. Phys. B 947, 114737 (2019), arXiv:1908.00805[hep-th doi: 10.1016/j.nuclphysb.2019.114737
    [95] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Rev. D 101(5), 055046 (2020), arXiv:1910.11553[hep-ph doi: 10.1103/PhysRevD.101.055046
    [96] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 07, 165 (2019), arXiv:1905.11970[hep-ph doi: 10.1007/JHEP07(2019)165
    [97] M. Tanimoto and K. Yamamoto, arXiv: 2106.10919 [hep-ph].
    [98] M. C. Chen, S. Ramos-Sánchez, and M. Ratz, Phys. Lett. B 801, 135153 (2020), arXiv:1909.06910[hep-ph doi: 10.1016/j.physletb.2019.135153
    [99] I. de Medeiros Varzielas, M. Levy, and Y. L. Zhou, JHEP 11, 085 (2020), arXiv:2008.05329[hep-ph doi: 10.1007/JHEP11(2020)085
    [100] A. Das, T. Nomura, H. Okada et al., Phys. Rev. D 96(7), 075001 (2017), arXiv:1704.02078[hep-ph doi: 10.1103/PhysRevD.96.075001
    [101] W. Wang and Z. L. Han, Phys. Rev. D 92, 095001 (2015), arXiv:1508.00706 [hep-ph doi: 10.1103/PhysRevD.92.095001
    [102] A. Gando et al. (KamLAND-Zen Collaboration), Phys. Rev. Lett. 117(8), 082503 (2016) [Addendum: Phys. Rev. Lett. 117(10), 109903 (2016)] arXiv: 1605.02889 [hep-ex]
    [103] I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo et al., JHEP 1901, 106 (2019), arXiv:1811.05487[hep-ph
    [104] I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo et al., NuFIT 4.1 (2019), www.nu-fit.org, (2019)
    [105] E. Fernandez-Martinez, J. Hernandez-Garcia, and J. Lopez-Pavon, JHEP 08, 033 (2016), arXiv:1605.08774[hep-ph doi: 10.1007/JHEP08(2016)033
    [106] N. R. Agostinho, G. C. Branco, P. M. F. Pereira et al., Eur. Phys. J. C 78(11), 895 (2018), arXiv:1711.06229[hep-ph doi: 10.1140/epjc/s10052-018-6347-2
  • [1] P. Minkowski, Phys. Lett. B 67, 421 (1977) doi: 10.1016/0370-2693(77)90435-X
    [2] T. Yanagida, in Proceedings of the Workshop on the Unified Theory and the Baryon Number in the Universe (O. Sawada and A. Sugamoto, eds.), KEK, Tsukuba, Japan, 1979, p. 95
    [3] M. Gell-Mann, P. Ramond, and R. Slansky, Supergravity (P. van Nieuwenhuizen et al. eds.), North Holland, Amsterdam, 1979, p. 315; S. L. Glashow, The future of elementary particle physics, in Proceedings of the 1979 Cargèse Summer Institute on Quarks and Leptons (M. Levy et al. eds.), Plenum Press, New York, 1980, p. 687
    [4] R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980) doi: 10.1103/PhysRevLett.44.912
    [5] R. N. Mohapatra and J. W. F. ValleSenjanovic, Phys. Rev. D 34, 1642 (1986)
    [6] D. Wyler and L. Wolfenstein, Nucl. Phys. B 218, 205 (1983) doi: 10.1016/0550-3213(83)90482-0
    [7] E. K. Akhmedov, M. Lindner, E. Schnapka et al., Phys. Lett. B 368, 270 (1996), arXiv:[hep-ph/9507275 doi: 10.1016/0370-2693(95)01504-3
    [8] E. K. Akhmedov, M. Lindner, E. Schnapka et al., Phys. Rev. D 53, 2752 (1996), arXiv:[hep-ph/9509255
    [9] M. Malinsky, J. C. Romao, and J. W. F. Valle, Phys. Rev. Lett. 95, 161801 (2005), arXiv:hep-ph/0506296 [hep-ph doi: 10.1103/PhysRevLett.95.161801
    [10] F. Feruglio, doi: 10.1142/9789813238053_0012 arXiv: 1706.08749[hep-ph]
    [11] R. de Adelhart Toorop, F. Feruglio, and C. Hagedorn, Nucl. Phys. B 858, 437-467 (2012), arXiv:1112.1340[hep-ph doi: 10.1016/j.nuclphysb.2012.01.017
    [12] J. C. Criado and F. Feruglio, SciPost Phys. 5(5), 042 (2018), arXiv:1807.01125[hep-ph doi: 10.21468/SciPostPhys.5.5.042
    [13] T. Kobayashi, N. Omoto, Y. Shimizu et al., JHEP 11, 196 (2018), arXiv:1808.03012[hep-ph doi: 10.1007/JHEP11(2018)196
    [14] H. Okada and M. Tanimoto, Phys. Lett. B 791, 54-61 (2019), arXiv:1812.09677[hep-ph doi: 10.1016/j.physletb.2019.02.028
    [15] T. Nomura and H. Okada, Phys. Lett. B 797, 134799 (2019), arXiv:1904.03937[hep-ph doi: 10.1016/j.physletb.2019.134799
    [16] H. Okada and M. Tanimoto, Eur. Phys. J. C 81(1), 52 (2021), arXiv:1905.13421[hep-ph doi: 10.1140/epjc/s10052-021-08845-y
    [17] F. J. de Anda, S. F. King, and E. Perdomo, Phys. Rev. D 101(1), 015028 (2020), arXiv:1812.05620[hep-ph doi: 10.1103/PhysRevD.101.015028
    [18] P. P. Novichkov, S. T. Petcov, and M. Tanimoto, Phys. Lett. B 793, 247-258 (2019), arXiv:1812.11289[hep-ph doi: 10.1016/j.physletb.2019.04.043
    [19] T. Nomura and H. Okada, Nucl. Phys. B 966, 115372 (2021), arXiv:1906.03927[hep-ph doi: 10.1016/j.nuclphysb.2021.115372
    [20] H. Okada and Y. Orikasa, arXiv: 1907.13520 [hep-ph]
    [21] G. J. Ding, S. F. King, and X. G. Liu, JHEP 09, 074 (2019), arXiv:1907.11714[hep-ph doi: 10.1007/JHEP09(2019)074
    [22] T. Nomura, H. Okada, and O. Popov, Phys. Lett. B 803, 135294 (2020), arXiv:1908.07457[hep-ph doi: 10.1016/j.physletb.2020.135294
    [23] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Rev. D 100(11), 115045 (2019) [Erratum: Phys. Rev. D 101(3), 039904 (2020)] arXiv: 1909.05139[hep-ph]
    [24] T. Asaka, Y. Heo, T. H. Tatsuishi a et al., JHEP 01, 144 (2020), arXiv:1909.06520[hep-ph doi: 10.1007/JHEP01(2020)144
    [25] D. Zhang, Nucl. Phys. B 952, 114935 (2020), arXiv:1910.07869[hep-ph doi: 10.1016/j.nuclphysb.2020.114935
    [26] G. J. Ding, S. F. King, X. G. Liu et al., JHEP 12, 030 (2019), arXiv:1910.03460[hep-ph doi: 10.1007/JHEP12(2019)030
    [27] T. Kobayashi, T. Nomura, and T. Shimomura, Phys. Rev. D 102(3), 035019 (2020), arXiv:1912.00637[hep-ph doi: 10.1103/PhysRevD.102.035019
    [28] T. Nomura, H. Okada, and S. Patra, Nucl. Phys. B 967, 115395 (2021), arXiv:1912.00379[hep-ph doi: 10.1016/j.nuclphysb.2021.115395
    [29] X. Wang, Nucl. Phys. B 957, 115105 (2020), arXiv:1912.13284[hep-ph doi: 10.1016/j.nuclphysb.2020.115105
    [30] H. Okada and Y. Shoji, Nucl. Phys. B 961, 115216 (2020), arXiv:2003.13219[hep-ph doi: 10.1016/j.nuclphysb.2020.115216
    [31] H. Okada and M. Tanimoto, arXiv: 2005.00775 [hep-ph]
    [32] M. K. Behera, S. Singirala, S. Mishra et al., arXiv: 2009.01806 [hep-ph]
    [33] M. K. Behera, S. Mishra, S. Singirala et al., arXiv: 2007.00545 [hep-ph]
    [34] T. Nomura and H. Okada, arXiv: 2007.04801 [hep-ph]
    [35] T. Nomura and H. Okada, arXiv: 2007.15459 [hep-ph]
    [36] T. Asaka, Y. Heo, and T. Yoshida, Phys. Lett. B 811, 135956 (2020), arXiv:2009.12120[hep-ph doi: 10.1016/j.physletb.2020.135956
    [37] H. Okada and M. Tanimoto, Phys. Rev. D 103(1), 015005 (2021), arXiv:2009.14242[hep-ph doi: 10.1103/PhysRevD.103.015005
    [38] K. I. Nagao and H. Okada, arXiv: 2010.03348 [hep-ph]
    [39] H. Okada and M. Tanimoto, JHEP 03, 010 (2021), arXiv:2012.01688[hep-ph doi: 10.1007/JHEP03(2021)010
    [40] C. Y. Yao, J. N. Lu, and G. J. Ding, JHEP 05, 102 (2021), arXiv:2012.13390 [hep-ph doi: 10.1007/JHEP05(2021)102
    [41] P. Chen, G. J. Ding, and S. F. King, JHEP 04, 239 (2021), arXiv:2101.12724 [hep-ph doi: 10.1007/JHEP04(2021)239
    [42] M. Kashav and S. Verma, arXiv: 2103.07207 [hep-ph]
    [43] H. Okada, Y. Shimizu, M. Tanimoto et al., arXiv: 2105.14292 [hep-ph]
    [44] I. de Medeiros Varzielas, and J. Lourenço, arXiv: 2107.04042 [hep-ph]
    [45] T. Nomura, H. Okada, and Y. Orikasa, arXiv: 2106.12375 [hep-ph]
    [46] P. T. P. Hutauruk, D. W. Kang, J. Kim et al., arXiv: 2012.11156 [hep-ph]
    [47] G. J. Ding, S. F. King, and J. N. Lu, arXiv: 2108.09655 [hep-ph]
    [48] K. I. Nagao and H. Okada, arXiv: 2108.09984 [hep-ph]
    [49] T. Kobayashi, K. Tanaka, and T. H. Tatsuishi, Phys. Rev. D 98(1), 016004 (2018), arXiv:1803.10391[hep-ph doi: 10.1103/PhysRevD.98.016004
    [50] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Lett. B 794, 114-121 (2019), arXiv:1812.11072[hep-ph doi: 10.1016/j.physletb.2019.05.034
    [51] T. Kobayashi, Y. Shimizu, K. Takagi et al., PTEP 2020(5), 053B05 (2020), arXiv:1906.10341[hep-ph doi: 10.1093/ptep/ptaa055
    [52] H. Okada and Y. Orikasa, Phys. Rev. D 100(11), 115037 (2019), arXiv:1907.04716[hep-ph doi: 10.1103/PhysRevD.100.115037
    [53] S. Mishra, arXiv: 2008.02095 [hep-ph]
    [54] X. Du and F. Wang, JHEP 02, 221 (2021), arXiv:2012.01397[hep-ph doi: 10.1007/JHEP02(2021)221
    [55] J. T. Penedo and S. T. Petcov, Nucl. Phys. B 939, 292-307 (2019), arXiv:1806.11040[hep-ph doi: 10.1016/j.nuclphysb.2018.12.016
    [56] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 04, 005 (2019), arXiv:1811.04933[hep-ph doi: 10.1007/JHEP04(2019)005
    [57] T. Kobayashi, Y. Shimizu, K. Takagi et al., JHEP 02, 097 (2020), arXiv:1907.09141[hep-ph doi: 10.1007/JHEP02(2020)097
    [58] S. F. King and Y. L. Zhou, Phys. Rev. D 101(1), 015001 (2020), arXiv:1908.02770[hep-ph doi: 10.1103/PhysRevD.101.015001
    [59] H. Okada and Y. Orikasa, arXiv: 1908.08409 [hep-ph]
    [60] J. C. Criado, F. Feruglio, and S. J. D. King, JHEP 02, 001 (2020), arXiv:1908.11867[hep-ph doi: 10.1007/JHEP02(2020)001
    [61] X. Wang and S. Zhou, JHEP 05, 017 (2020), arXiv:1910.09473[hep-ph doi: 10.1007/JHEP05(2020)017
    [62] Y. Zhao and H. H. Zhang, JHEP 03, 002 (2021), arXiv:2101.02266 [hep-ph doi: 10.1007/JHEP03(2021)002
    [63] S. F. King and Y. L. Zhou, JHEP 04, 291 (2021), arXiv:2103.02633 [hep-ph doi: 10.1007/JHEP04(2021)291
    [64] G. J. Ding, S. F. King, and C. Y. Yao, arXiv: 2103.16311 [hep-ph]
    [65] X. Zhang and S. Zhou, arXiv: 2106.03433 [hep-ph]
    [66] Bu-Yao Qu, Xiang-Gan Liu, Ping-Tao Chen et al., arXiv: 2106.11659 [hep-ph]
    [67] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 04, 174 (2019), arXiv:1812.02158[hep-ph doi: 10.1007/JHEP04(2019)174
    [68] G. J. Ding, S. F. King, and X. G. Liu, Phys. Rev. D 100(11), 115005 (2019), arXiv:1903.12588[hep-ph doi: 10.1103/PhysRevD.100.115005
    [69] X. Wang, B. Yu, and S. Zhou, Phys. Rev. D 103(7), 076005 (2021), arXiv:2010.10159[hep-ph doi: 10.1103/PhysRevD.103.076005
    [70] C. Y. Yao, X. G. Liu, and G. J. Ding, Phys. Rev. D 103(9), 095013 (2021), arXiv:2011.03501[hep-ph doi: 10.1103/PhysRevD.103.095013
    [71] X. Wang and S. Zhou, arXiv: 2102.04358 [hep-ph]
    [72] M. K. Behera and R. Mohanta, arXiv: 2108.01059 [hep-ph]
    [73] A. Baur, H. P. Nilles, A. Trautner et al., Phys. Lett. B 795, 7-14 (2019), arXiv:1901.03251[hep-th doi: 10.1016/j.physletb.2019.03.066
    [74] I. de Medeiros Varzielas, S. F. King et al., Phys. Rev. D 101(5), 055033 (2020), arXiv:1906.02208[hep-ph doi: 10.1103/PhysRevD.101.055033
    [75] X. G. Liu and G. J. Ding, JHEP 08, 134 (2019), arXiv:1907.01488[hep-ph doi: 10.1007/JHEP08(2019)134
    [76] P. Chen, G. J. Ding, J. N. Lu et al., Phys. Rev. D 102(9), 095014 (2020), arXiv:2003.02734[hep-ph doi: 10.1103/PhysRevD.102.095014
    [77] C. C. Li, X. G. Liu, and G. J. Ding, arXiv: 2108.02181 [hep-ph]
    [78] P. P. Novichkov, J. T. Penedo, and S. T. Petcov, Nucl. Phys. B 963, 115301 (2021), arXiv:2006.03058[hep-ph doi: 10.1016/j.nuclphysb.2020.115301
    [79] X. G. Liu, C. Y. Yao, and G. J. Ding, Phys. Rev. D 103(5), 056013 (2021), arXiv:2006.10722[hep-ph doi: 10.1103/PhysRevD.103.056013
    [80] S. Kikuchi, T. Kobayashi, H. Otsuka et al., JHEP 11, 101 (2020), arXiv:2007.06188[hep-th doi: 10.1007/JHEP11(2020)101
    [81] Y. Almumin, M. C. Chen, V. Knapp-Pérez et al., JHEP 05, 078 (2021), arXiv:2102.11286 [hep-th doi: 10.1007/JHEP05(2021)078
    [82] G. J. Ding, F. Feruglio, and X. G. Liu, SciPost Phys. 10, 133 (2021), arXiv:2102.06716 [hep-ph doi: 10.21468/SciPostPhys.10.6.133
    [83] F. Feruglio, V. Gherardi, A. Romanino et al., JHEP 05, 242 (2021), arXiv:2101.08718 [hep-ph doi: 10.1007/JHEP05(2021)242
    [84] S. Kikuchi, T. Kobayashi, and H. Uchida, arXiv: 2101.00826 [hep-th]
    [85] P. P. Novichkov, J. T. Penedo, and S. T. Petcov, JHEP 04, 206 (2021), arXiv:2102.07488 [hep-ph doi: 10.1007/JHEP04(2021)206
    [86] G. Altarelli and F. Feruglio, Rev. Mod. Phys. 82, 2701-2729 (2010), arXiv:1002.0211[hep-ph doi: 10.1103/RevModPhys.82.2701
    [87] H. Ishimori, T. Kobayashi, H. Ohki et al., Prog. Theor. Phys. Suppl. 183, 1-163 (2010), arXiv:1003.3552[hep-th doi: 10.1143/PTPS.183.1
    [88] H. Ishimori, T. Kobayashi, H. Ohki et al., Notes Phys. 858, 1-227 (2012) doi: 10.1007/978-3-642-30805-5
    [89] D. Hernandez and A. Y. Smirnov, Phys. Rev. D 86, 053014 (2012), arXiv:1204.0445[hep-ph doi: 10.1103/PhysRevD.86.053014
    [90] S. F. King and C. Luhn, Rept. Prog. Phys. 76, 056201 (2013), arXiv:1301.1340[hep-ph doi: 10.1088/0034-4885/76/5/056201
    [91] S. F. King, A. Merle, S. Morisi et al., New J. Phys. 16, 045018 (2014), arXiv:1402.4271[hep-ph doi: 10.1088/1367-2630/16/4/045018
    [92] S. F. King, Prog. Part. Nucl. Phys. 94, 217-256 (2017), arXiv:1701.04413[hep-ph doi: 10.1016/j.ppnp.2017.01.003
    [93] S. T. Petcov, Eur. Phys. J. C 78(9), 709 (2018), arXiv:1711.10806[hep-ph doi: 10.1140/epjc/s10052-018-6158-5
    [94] A. Baur, H. P. Nilles, A. Trautner et al., Nucl. Phys. B 947, 114737 (2019), arXiv:1908.00805[hep-th doi: 10.1016/j.nuclphysb.2019.114737
    [95] T. Kobayashi, Y. Shimizu, K. Takagi et al., Phys. Rev. D 101(5), 055046 (2020), arXiv:1910.11553[hep-ph doi: 10.1103/PhysRevD.101.055046
    [96] P. P. Novichkov, J. T. Penedo, S. T. Petcov et al., JHEP 07, 165 (2019), arXiv:1905.11970[hep-ph doi: 10.1007/JHEP07(2019)165
    [97] M. Tanimoto and K. Yamamoto, arXiv: 2106.10919 [hep-ph].
    [98] M. C. Chen, S. Ramos-Sánchez, and M. Ratz, Phys. Lett. B 801, 135153 (2020), arXiv:1909.06910[hep-ph doi: 10.1016/j.physletb.2019.135153
    [99] I. de Medeiros Varzielas, M. Levy, and Y. L. Zhou, JHEP 11, 085 (2020), arXiv:2008.05329[hep-ph doi: 10.1007/JHEP11(2020)085
    [100] A. Das, T. Nomura, H. Okada et al., Phys. Rev. D 96(7), 075001 (2017), arXiv:1704.02078[hep-ph doi: 10.1103/PhysRevD.96.075001
    [101] W. Wang and Z. L. Han, Phys. Rev. D 92, 095001 (2015), arXiv:1508.00706 [hep-ph doi: 10.1103/PhysRevD.92.095001
    [102] A. Gando et al. (KamLAND-Zen Collaboration), Phys. Rev. Lett. 117(8), 082503 (2016) [Addendum: Phys. Rev. Lett. 117(10), 109903 (2016)] arXiv: 1605.02889 [hep-ex]
    [103] I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo et al., JHEP 1901, 106 (2019), arXiv:1811.05487[hep-ph
    [104] I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo et al., NuFIT 4.1 (2019), www.nu-fit.org, (2019)
    [105] E. Fernandez-Martinez, J. Hernandez-Garcia, and J. Lopez-Pavon, JHEP 08, 033 (2016), arXiv:1605.08774[hep-ph doi: 10.1007/JHEP08(2016)033
    [106] N. R. Agostinho, G. C. Branco, P. M. F. Pereira et al., Eur. Phys. J. C 78(11), 895 (2018), arXiv:1711.06229[hep-ph doi: 10.1140/epjc/s10052-018-6347-2
  • 加载中

Cited by

1. Ding, G.-J., Lu, J.-N., Petcov, S.T. et al. Non-holomorphic modular S4 lepton flavour models[J]. Journal of High Energy Physics, 2025, 2025(1): 191. doi: 10.1007/JHEP01(2025)191
2. Palavrić, A.. Discrete leptonic flavor symmetries: UV mediators and phenomenology[J]. Physical Review D, 2024, 110(11): 115025. doi: 10.1103/PhysRevD.110.115025
3. Ding, G.-J., King, S.F. Neutrino mass and mixing with modular symmetry[J]. Reports on Progress in Physics, 2024, 87(8): 084201. doi: 10.1088/1361-6633/ad52a3
4. Kobayashi, T., Tanimoto, M. Modular flavor symmetric models[J]. International Journal of Modern Physics A, 2024, 39(9-10): 2441012. doi: 10.1142/S0217751X24410124
5. Kobayashi, T., Tanimoto, M. Modular Flavor Symmetric Models[J]. Harald Fritzsch: Memorial Volume, 2024. doi: 10.1142/9789811292279_0012
6. de Medeiros Varzielas, I., Levy, M., Penedo, J.T. et al. Quarks at the modular S 4 cusp[J]. Journal of High Energy Physics, 2023, 2023(9): 196. doi: 10.1007/JHEP09(2023)196
7. Calò, S., Marinissen, C., Rahn, R. Discrete symmetries and efficient counting of operators[J]. Journal of High Energy Physics, 2023, 2023(5): 215. doi: 10.1007/JHEP05(2023)215
8. Du, X.K., Wang, F. Flavor structures of quarks and leptons from flipped SU(5) GUT with A4 modular flavor symmetry[J]. Journal of High Energy Physics, 2023, 2023(1): 36. doi: 10.1007/JHEP01(2023)036
9. Kikuchi, S., Kobayashi, T., Otsuka, H. et al. 4D modular flavor symmetric models inspired by a higher-dimensional theory[J]. Physical Review D, 2022, 106(3): 034507. doi: 10.1103/PhysRevD.106.035001
10. Kobayashi, T., Otsuka, H., Tanimoto, M. et al. Lepton flavor violation, lepton (g − 2) μ, e and electron EDM in the modular symmetry[J]. Journal of High Energy Physics, 2022, 2022(8): 13. doi: 10.1007/JHEP08(2022)013

Figures(3) / Tables(2)

Get Citation
Takaaki Nomura and Hiroshi Okada. Linear seesaw model with a modular S4 flavor symmetry[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac4975
Takaaki Nomura and Hiroshi Okada. Linear seesaw model with a modular S4 flavor symmetry[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac4975 shu
Milestone
Received: 2021-12-23
Article Metric

Article Views(2152)
PDF Downloads(36)
Cited by(10)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Linear seesaw model with a modular S4 flavor symmetry

    Corresponding author: Takaaki Nomura, nomura@scu.edu.cn
    Corresponding author: Hiroshi Okada, hiroshi.okada@apctp.org
  • 1. College of Physics, Sichuan University, Chengdu 610065, China
  • 2. Asia Pacific Center for Theoretical Physics (APCTP) - Headquarters San 31, Hyoja-dong, Nam-gu, Pohang 790-784, Korea
  • 3. Department of Physics, Pohang University of Science and Technology, Pohang 37673, Korea

Abstract: We discuss a linear seesaw model with a field content as minimum as possible, introducing a modular S_4 using gauged U(1)_{B-L} symmetries. Owing to the rank two neutrino mass matrix, we obtain a vanishing neutrino mass eigenvalue, and only the normal mass hierarchy of neutrinos is favored via the modular S_4 symmetry. In our numerical \Delta \chi^2 analysis, we especially determine a relatively sharp prediction on the sum of neutrino masses to be approximately 60 meV, in addition to other predictions.

    HTML

    I.   INTRODUCTION
    • The neutrino sector is theoretically unconfirmed in the standard model (SM), because only two mass squared differences and three mixings are experimentally found and the scale of mass is extremely minuscule compared to the other three sectors in the SM. Hence, several scientists expect the neutrino sector to possess new physics. A gauged B-L (baryon number minus lepton number) U(1) symmetry; U(1)_{B-L} , is a promising prescription for generating such tiny neutrino masses, introducing three right-handed neutrinos with rather heavy masses ( M_{\rm R} ), and is called a canonical seesaw model [1-4]. Because the mass scale is sometimes expected to be that of a grand unified theory ( M_{\rm GUT}\sim10^{15} GeV) to be small neutrino masses, its scale cannot be verified by our current experiments, and the spontaneous U(1)_{B-L} symmetry breaking scale is naturally expected to be the same energy as the cut-off scale; M_{\rm R}\sim M_{\rm GUT} .

      To achieve a successful neutrino mass model within our scale (\sim TeV), heavy neutral fermions ( S_{\rm L} ) with left-handed chirality are introduced along this line of thought. Inverse seesaw model [5, 6] also requests both of N_{\rm R} and S_{\rm L} and the neutrino mass could be realized within TeV scale. But this model may not require GUT scale. Currently, N_{\rm R} and S_{\rm L} can be embedded into the middle scale with SU(3)_{\rm C}\;\otimes SU(2)_{\rm L}\otimes SU(2)_{\rm R} \otimes U(1)_{B-L} [6-8], which is included in SO(10) group [9]. Then, the neutrino mass matrix would not be suppressed by the middle scale, but by M_{\rm GUT} scale when appropriate charge assignments are assigned for each fields in a supersymmetric theory [9]. Hence, we can test phenomenologies with our current experiments, supposing the middle scale breaking can occur at the TeV scale. This type of model is called the "linear seesaw". Because these models typically require more free parameters than the other three sectors in the SM fermion, flavor symmetries are also frequently introduced in these models, to reduce the parameters and obtain predictions (if possible).

      In 2017, attractive flavor symmetries were proposed in Refs. [10, 11], in which the authors applied modular non-Abelian discrete flavor symmetries to quark and lepton sectors. A remarkable advantage is that any dimensionless coupling can also be transformed as non-trivial representations under these symmetries. Therefore, we do not need so many scalars to determine a predictive mass matrix. Another advantage is that we have a modular weight from the modular origin that can play a role in stabilizing DM when appropriate charge assignments are distributed to each of the fields of the models. Along this line of thought, significantly relevant references have recently appeared in the literature, e.g., A_4 [10, 12-48], S_3 [49-54], S_4 [55-66], A_5 [60, 67, 68], double covering of A_5 [69-72], larger groups [73], multiple modular symmetries [74], double covering of A_4 [75-77], S_4 [78, 79], and other types of groups [80-85] in which masses, mixing, and CP phases for the quark and/or lepton have been predicted. For interested readers, we provide some literature reviews, which are useful to understand the non-Abelian group and its applications to flavor structure [86-93]. Moreover, a systematic approach to understanding the origin of CP transformations has been discussed in Ref. [94]; in addition, CP/flavor violation in the models with modular symmetry was discussed in Refs. [95-97], and a possible correction from the Kähler potential was discussed in Ref. [98]. Furthermore, a systematic analysis of the fixed points (stabilizers) was discussed in Ref. [99]. It would be interesting to consider a linear see-saw model with the local U(1)_{B-L} (\subset SO(10)) under modular symmetry because these symmetries can be obtained from a string theory. Moreover, the nature of the modular symmetry can be used to realize the linear seesaw mechanism, in addition to constraining the flavor structure.

      In this study, we propose a linear seesaw model under modular S_4 , using the U(1)_{B-L} symmetry, in which we try to construct the predictive model as minimum as possible. Owing to the rank two neutrino mass matrix, we obtain a vanishing neutrino mass eigenvalue. Furthermore, only the normal mass hierarchy of neutrinos is favored via the modular S_4 symmetry. In our numerical analysis, we perform \Delta \chi^2 analysis in the neutrino sector, considering non-unitarity constraint.

      This remainder of this paper is organized as follows. In Sec. II, we review our model, constructing renormalizable Lagrangian and mass matrices in the lepton sector. Then, we formulate the neutrino mass matrix with rank two, in which we estimate the structure of the neutrino mass matrix in the expansion of modulus. In addition, we derive several observables in the lepton sector. At the end of this section, we discuss the non-unitarity bound. In Sec. III, we perform \Delta \chi^2 analysis in the lepton sector, and present some predictions using our model. In Sec. IV, we present a summary and discussion. In the Appendix, we elucidate the modular S_4 symmetry.

    II.   MODEL

      A.   Model review

    • In this section, we review our model framework for the linear seesaw mechanism, introducing the B-L local Abelian symmetry; U(1)_{B-L} and modular A_4 symmetry. Regarding the fermion sector, we add two left-handed neutral fermions \{S_{L_1},S_{L_2} \} that belong to the isospin singlet, where they exhibit a zero charge under U(1)_{B-L} , \{{\boldsymbol{1}},{\boldsymbol{1'}}\} under S_4 , and -1 under -k_I , respectively. In addition, we introduce three right-handed neutral fermions \overline{N_{\rm R}} that belong to the isospin singlet, where they exhibit 1 charge under U(1)_{B-L} , {\boldsymbol{3}} under S_4 , and -3 under -k , respectively. The SM left-handed leptons L_{\rm L}\equiv [L_{{\rm L}_e},L_{{\rm L}_\mu}, L_{{\rm L}_\tau}]^T belong to -1 charge under U(1)_{B-L} , {\boldsymbol{3'}} under S_4 , and -1 under -k , respectively. However, the SM right-handed leptons \{ \overline{e_{\rm R}},[\overline{\mu_{\rm R}},\overline{\tau_{\rm R}}] \} belong to the +1 charge under U(1)_{B-L} , \{ {\bf 1}, {\bf 2} \} under S_4 , and \{ -1,-3 \} under -k_I , respectively.

      Regarding the scalar sector, we adopt two Higgs doublets, H_1, \; H_2 , and an isospin singlet field φ. The isospin singlet φ has a -1 charge under U(1)_{B-L} . Here, H_1 is the SM-like Higgs that has a zero charge under U(1)_{B-L} and -k_I while H_2 has 1 charge under U(1)_{B-L} , and zero under -k_I . We denote each of the vacuum expectation values (VEVs) to be \langle H_{1,2} \rangle\equiv [0,v_{1,2}/\sqrt2]^T , and \langle \varphi \rangle\equiv v_{\varphi} /\sqrt2 . We summarize our particle content and assignments in Table 1.

      Fermions Scalars
      L_{\rm L} \overline{e_{\rm R}},[\overline{\mu_{\rm R}},\overline{\tau_{\rm R}}] \overline{N_{\rm R}} S_{{\rm L}_1},S_{{\rm L}_2} H_1 H_2 φ
      SU(2)_{\rm L} {\bf{2}} {\bf{1}} {\bf{1}} {\bf{1}} {\bf{2}} {\bf{2}} {\bf{1}}
      U(1)_{\rm Y} -\dfrac12 1 0 0 \dfrac12 \dfrac12 0
      U(1)_{B-L} -1 1 1 0 0 1 -1
      S_4 3' 1, 2 {\bf{3}} {\bf{1}},{\bf{1'}} {\bf{1}} {\bf{1}} {\bf{1}}
      -k_I -1 -1,-3 -3 -1 0 0 0

      Table 1.  Lepton and boson particle contents and their charge assignments under SU(2)_{\rm L}\times U(1)_{\rm Y}\times U(1)_{B-L}\times S_4\times (-k_I) where L_{\rm L}\equiv [L_{{\rm L}_e},L_{{\rm L}_\mu}, L_{{\rm L}_\tau}]^T~ k_I represents the number of the modular weight.

      Then the valid lepton Yukawa Lagrangian is expressed as

      \begin{align} -{\cal{L}}_{\rm{lepton}} = {\cal{L}}_{M_\ell} + {\cal{L}}_{{M_{\rm D}}}+{\cal{L}}_{{M'_{\rm D}}} +{\cal{L}}_{M_{NS}} , \end{align}

      (1)

      where {\cal{L}}_{M_\ell} denotes the charged lepton Yukawa Lagrangian. {\cal{L}}_{{M_{\rm D}}} belongs to [ \overline{N_{\rm R}} L_{\rm L} \tilde H_1 ] where \tilde H\equiv {\rm i}\sigma_2 H^* . {\cal{L}}_{{M'_{\rm D}}} belongs to [ \overline{L_{\rm L}^{\rm C}} S_{\rm L} H_2 ]. {\cal{L}}_{M_{NS}} belongs to [ \overline{N_{\rm R}} S_{\rm L} \varphi ]. [\cdots] implies that the concrete flavor structures are manifolded. Each of their structures are presented below.

      The scalar potential of our model is expressed as

      \begin{aligned}[b] V =& m_\varphi^2 \varphi^*\varphi + m_1^2 H_1^\dagger H_1 + m_2^2 H_2^\dagger H_2 - \mu_{12} (H_1^\dagger H_2 \varphi + {\rm h.c.}) \\ & + \lambda_1 (H_1^\dagger H_1)^2+ \lambda_2 (H_2^\dagger H_2)^2 + \lambda_\varphi (\varphi^* \varphi)^2 \\&+ \lambda_3 (H_1^\dagger H_1)(H_2^\dagger H_2) + \lambda_4 (H_1^\dagger H_2) (H_2^\dagger H_1) \\ & + \lambda_{\varphi H_1} (H_1^\dagger H_1) (\varphi^* \varphi)+ \lambda_{\varphi H_2} (H_2^\dagger H_2) (\varphi^* \varphi), \end{aligned}

      (2)

      where {\rm h.c.} stands for Hermitian conjugate. We consider that φ develops a VEV at a significantly larger scale than H_{1,2} . Then, after the φ developes a VEV, the scalar potential matches with that of two Higgs doublet model (THDM) without the (H_1^\dagger H_2)^2 term owing to the U(1)_{B-L} symmetry. In addition, the Yukawa couplings associated with the two Higgs doublet belong to the type-I THDM because only H_1 can couple to the SM fermions. Hence, in our analysis below, we do not discuss the THDM part further, but focus on the neutrino mass.

    • B.   Valid Lagrangians

    • Before discussing the valid Lagrangians for the lepton sector, we define Yukawa couplings under the modular S_4 symmetry as follows: the S_4 doublets with -k_I = 2 and 4 are denoted by Y^{(2)}_{\bf 2}\equiv [y_1,y_2]^T and [y'_1,y'_2]^T , respectively. The S_4 triplet with -k_I = 2 is Y^{(2)}_{\bf 3'}\equiv [y_3,y_4,y_5]^T , while the S_4 triplets with -k_I = 4 are Y^{(4)}_{\bf 3}\equiv [y'_3,y'_4,y'_5]^T and Y^{(4)}_{\bf 3'}\equiv [y''_3,y''_4,y''_5]^T . Each of these structures are explicitly presented in the Appendix.

      Charged lepton mass matrix:

      The renormalizable Lagrangian for the charged-lepton sector is expressed as

      \begin{aligned}[b]\\[-8pt] -{\cal{L}}_{M_\ell} = &\alpha_\ell \overline{e_{\rm R}} (y_3 L_{{\rm L}_e} +y_4 \bar L_{{\rm L}_\tau} +y_5 \bar L_{{\rm L}_\mu}) \tilde H_1 + \beta_\ell\left[\frac{\sqrt3}{2}\overline{\mu_{\rm R}}(y'_4 L_{{\rm L}_\mu} +y'_5 L_{{\rm L}_\tau}) +\overline{\tau_{\rm R}} \left(-y'_3 L_{{\rm L}_e} + \frac12 (y'_4 L_{{\rm L}_\tau} +y'_5 L_{{\rm L}_\mu}) \right) \right] \tilde H_1 \\ &+ \gamma_\ell\left[\frac{\sqrt3}{2}\overline{\tau_{\rm R}}(y''_4 L_{{\rm L}_\mu} +y''_5 L_{{\rm L}_\tau}) +\overline{\mu_{\rm R}} \left(y''_3 L_{{\rm L}_e} - \frac12 (y''_4 L_{{\rm L}_\tau} +y''_5 L_{{\rm L}_\mu}) \right) \right] \tilde H_1 +{\rm{h.c.}}, \end{aligned}

      (3)

      where \{\alpha_\ell,\beta_\ell,\gamma_\ell \} are real parameters without loss of generality. Subsequently, the charged-lepton mass matrix after the spontaneous symmetry breaking is expressed as

      \begin{aligned} (M_\ell)_{\rm RL} = \frac{v_1}{\sqrt2} \begin{pmatrix} \alpha_\ell y_3 & \alpha_\ell y_5 & \alpha_\ell y_4 \\ \gamma_\ell y''_3 & \dfrac{\sqrt3}{2}\beta_\ell y'_4-\dfrac12 \gamma_\ell y''_5 & \dfrac{\sqrt3}{2}\beta_\ell y'_5-\dfrac12 \gamma_\ell y''_4 \\ -\beta_\ell y'_3 & \dfrac{\sqrt3}{2}\gamma_\ell y''_4+\dfrac12 \beta_\ell y'_5 & \dfrac{\sqrt3}{2}\gamma_\ell y''_5+\dfrac12 \beta_\ell y'_4 \\ \end{pmatrix} . \end{aligned}

      (4)

      The charged-lepton mass eigenvalues are obtained by diagonalizing {\rm{diag}}[m_e,m_\mu,m_\tau] = V_{{\rm R}_\ell}^\dagger M_\ell V_{{\rm L}_\ell} , where V_{{\rm L}_\ell, {\rm R}_\ell} represent unitary matrices. In our numerical analysis, we will determine the free parameters \{\alpha_\ell,\beta_\ell,\gamma_\ell\} , to fit the three charged-lepton mass eigenstates after providing all the numerical values, by applying the following relationship. Here, we fix \alpha_\ell,\beta_\ell,\gamma_\ell , to determine the experimental three charged-lepton masses by applying the following relationship:

      \begin{aligned}[b] &{\rm{Tr}}[{M_\ell}^\dagger M_\ell] = |m_e|^2 + |m_\mu|^2 + |m_\tau|^2,\\& {\rm{Det}}[{M_\ell}^\dagger M_\ell] = |m_e|^2 |m_\mu|^2 |m_\tau|^2, \\ &({\rm{Tr}}[{M_\ell}^\dagger] M_\ell)^2 -{\rm{Tr}}[({M_\ell}^\dagger M_\ell)^2] \\=& 2( |m_e|^2 |m_\mu|^2 + |m_\mu|^2 |m_\tau|^2+ |m_e|^2 |m_\tau|^2 ). \end{aligned}

      (5)

      Neutral fermion mass matrices:

      First, we construct the valid Lagrangian of the Dirac mass matrix {\cal{L}}_{{M_{\rm D}}} . This is given by

      \begin{aligned}[b]\\[-8pt] {\cal{L}}_{M_{\rm D}} =& \alpha_{\rm D} \left[y'_3\left(\overline{N_{{\rm R}_3}} L_{{\rm L}_\mu}-\overline{N_{{\rm R}_2}} L_{{\rm L}_\tau}\right) + y'_4\left(\overline{N_{{\rm R}_1}} L_{{\rm L}_\tau}-\overline{N_{{\rm R}_3}} L_{{\rm L}_e}\right) + y'_5\left(\overline{N_{{\rm 2}}} L_{{\rm L}_e}-\overline{N_{{\rm R}_1}} L_{{\rm L}_\mu}\right)\right] \tilde H_1 \\ &+ \beta_{\rm D} \left[y''_3\left(\overline{N_{{\rm R}_3}} L_{{\rm L}_\tau}-\overline{N_{{\rm R}_2}} L_{{\rm L}_\mu}\right) + y''_4\left(-\overline{N_{{\rm R}_1}} L_{{\rm L}_\mu}-\overline{N_{{\rm R}_2}} L_{{\rm L}_e}\right) + y''_5\left(\overline{N_{{\rm R}_1}} L_{{\rm L}_\tau}+\overline{N_{{\rm R}_3}} L_{{\rm L}_e}\right)\right] \tilde H_1 \\ &+ \gamma_{\rm D} \left[\frac{\sqrt3}{2}y'_1(\overline{N_{{\rm R}_2}} L_{{\rm L}_\mu}+\overline{N_{{\rm R}_3}} L_{{\rm L}_\tau}) + y'_2 \{-\overline{N_{{\rm R}_1}} L_{{\rm L}_e} +\frac12 \left(\overline{N_{{\rm R}_2}} L_{{\rm L}_\tau}+\overline{N_{{\rm R}_3}} L_{{\rm L}_\mu}\right)\}\right] \tilde H_1+{\rm{h.c.}}, \end{aligned}

      (6)

      where we suppose \alpha_{\rm D} to be real while \beta_{\rm D},\ \gamma_{\rm D} are complex after rephasing the fields. Similar to the charged-lepton sector, we determine the Dirac mass matrix as follows:

      \begin{aligned} (M_{\rm D})_{\rm RL} = \frac{v_1}{\sqrt2} \begin{pmatrix} -\gamma_{\rm D} y'_2 & -\alpha_{\rm D} y'_5 - \beta_{\rm D} y''_4& \alpha_{\rm D} y'_4+\beta_{\rm D} y''_5 \\ \alpha_{\rm D} y'_5 - \beta_{\rm D} y''_4 &-\beta_{\rm D} y''_3+ \dfrac{\sqrt3}{2}\gamma_{\rm D} y'_1 &-\alpha_{\rm D} y'_3 + \dfrac12 \gamma_{\rm D} y'_2 \\ - \alpha_{\rm D} y'_4+\beta_{\rm D} y''_5 &\alpha_{\rm D} y'_3 + \dfrac12 \gamma_{\rm D} y'_2 & \beta_{\rm D} y''_3+ \dfrac{\sqrt3}{2}\gamma_{\rm D} y'_1 \\ \end{pmatrix} . \end{aligned}

      (7)

      For convenience, to analyze the neutrino oscillation, we redefine M_{\rm D}\equiv \dfrac{v_1}{\sqrt2}\tilde{M_{\rm D}} .

      Another Dirac Lagrangian is induced via {\cal{L}}_{M'_{\rm D}} , which is given by

      \begin{align} {\cal{L}}_{{M'_{\rm D}}} & = \alpha'_{\rm D} \left[y_3 \overline{L^{\rm C}_{{\rm L}_e}} +y_4 \overline{L^{\rm C}_{{\rm L}_\tau}}+y_5 \overline{L^{\rm C}_{{\rm L}_\mu}}\right]S_{{\rm L}_1} H_2 +{\rm{h.c.}}, \end{align}

      (8)

      Subsequently, we determine another Dirac mass matrix as follows:

      \begin{align} (M'_{\rm D})_{L_{\rm L} S_{\rm L}} = \frac{v_2}{\sqrt2} \alpha'_{\rm D} \begin{pmatrix} y_3 & 0 \\ y_5 & 0 \\ y_4 & 0 \\ \end{pmatrix} . \end{align}

      (9)

      Being the same as the reason for M_{\rm D} , we redefine M'_{\rm D}\equiv\dfrac{v_2}{\sqrt2} \alpha'_{\rm D} \tilde{M'_{\rm D}} .

      The third term of the Lagrangian {\cal{L}}_{M_{NS}} is given by

      \begin{aligned}[b] {\cal{L}}_{M_{NS}} = & \alpha_{NS} \left[y'_3 \overline{N_{{\rm R}_e}} + y'_4 \overline{N_{{\rm R}_\tau}} + y'_5\overline{N_{{\rm R}_\mu}}\right] S_{{\rm L}_1}\varphi \\&+ \beta_{NS} \left[y''_3 \overline{N_{{\rm R}_e}} + y''_4 \overline{N_{{\rm R}_\tau}} + y''_5\overline{N_{{\rm R}_\mu}}\right] S_{{\rm L}_2}\varphi +{\rm{h.c.}}, \end{aligned}

      (10)

      where \alpha_{NS},\; \beta_{NS} are real without loss of generality. Accordingly, we obtain the mass matrix

      \begin{aligned} M_{NS} = \frac{v_\varphi}{\sqrt2} \begin{pmatrix} y'_3 & y''_3 \\ y'_5 & y''_5 \\ y'_4 & y''_4 \\ \end{pmatrix} \begin{pmatrix} \alpha_{NS} & 0 \\ 0 & \beta_{NS} \\ \end{pmatrix} \equiv \frac{v_\varphi}{\sqrt2} \tilde M_{N S} . \end{aligned}

      (11)

      In basis of [\nu_{\rm L},N_{\rm R}^{\rm C},S_{\rm L}]^T , the neutral fermion mass matrix is given by

      \begin{align} M_N = \begin{pmatrix} 0_{3\times3} & M_{\rm D}^T & M'_{\rm D} \\ M_{\rm D} & 0_{3\times3} & M_{NS} \\ m'^T_{\rm D} & M_{NS}^T & 0_{2\times2} \end{pmatrix} . \end{align}

      (12)

      Hence, the active neutrino mass matrix is expressed as

      \begin{aligned}[b] m_\nu =& M'_{\rm D} (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D} \\&+ [M'_{\rm D} (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D}]^T \\ =& \frac{v_1v_2}{\sqrt{2} v_\varphi} \Big( \tilde M'_{\rm D} (\tilde M_{NS}^T \tilde M_{NS})^{-1}\tilde M_{NS}^T\tilde M_{\rm D}\\& + [\tilde M'_{\rm D} (\tilde M_{NS}^T\tilde M_{NS})^{-1}\tilde M_{NS}^T \tilde M_{\rm D}]^T \Big) = \kappa \tilde m_\nu, \end{aligned}

      (13)

      where \kappa\equiv \dfrac{v_1v_2}{\sqrt{2} v_\varphi} and we assume the mass hierarchies among M_{\rm D}, M'_{\rm D} \ll M_{NS} . Mass hierarchies is dynamically achieved in Refs. [100, 101]. The neutrino mass eigenvalues are obtained as follows: D_\nu = \kappa D_\nu = U_\nu^T m_\nu U_\nu = \kappa U_\nu^T \tilde m_\nu U_\nu , where U_\nu is a unitary matrix. Then, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix is given by U_{\rm PMNS}\equiv V^\dagger_{{\rm L}_\ell} U_\nu . Notice that m_\nu is rank two; hence, the lightest neutrino mass is zero. Here κ is described by one experimental value and dimensionless neutrino mass eigenstates as:

      \begin{align} {\rm{(NH)}}:\ \kappa^2 = \frac{|\Delta m_{\rm{atm}}^2|}{\tilde D_{\nu_3}^2}, \quad {\rm{(IH)}}:\ \kappa^2 = \frac{|\Delta m_{\rm{atm}}^2|}{\tilde D_{\nu_2}^2}, \end{align}

      (14)

      where \Delta m_{\rm{atm}}^2 denotes the atmospheric neutrino mass-squared difference and NH (IH) stands for normal (inverted) ordering. Subsequently, the solar mass difference squared can be expressed in terms of κ as follows:

      \begin{align} {\rm{(NH)}}:\ \Delta m_{\rm{sol}}^2 = {\kappa^2} {\tilde D_{\nu_2}^2}, \quad {\rm{(IH)}}:\ \Delta m_{\rm{sol}}^2 = {\kappa^2}({\tilde D_{\nu_2}^2-\tilde D_{\nu_1}^2}), \end{align}

      (15)

      which can be compared to the observed value. In other words, we explicitly express the mass eigenvalues in terms of \Delta m_{\rm{atm}}^2 and \Delta m_{\rm{sol}}^2 as:

      \begin{align} &{\rm{(NH)}}:\ D_{\nu_1}^2 = 0, \quad \ D_{\nu_2}^2 = \Delta m_{\rm{sol}}^2, \quad \ D_{\nu_3}^2 = \Delta m_{\rm{atm}}^2, \end{align}

      (16)

      \begin{align} &{\rm{(IH)}}:\ D_{\nu_1}^2 = \Delta m_{\rm{atm}}^2 - \Delta m_{\rm{sol}}^2, \;\;\; \ D_{\nu_2}^2 = \Delta m_{\rm{sol}}^2,\;\;\; \ D_{\nu_3}^2 = 0, \end{align}

      (17)

      which implies that NH is hierarchical, but IH is degenerate, as \Delta m_{\rm{sol}}^2/ \Delta m_{\rm{atm}}^2<<1 . Here, we expand |\tilde m_\nu|^2 in terms of q\equiv {\rm e}^{2p\pi {\rm i}\tau}<<1 . Then, the mass matrix is given by

      \begin{aligned} |\tilde m_\nu|^2\sim \begin{pmatrix} {\cal{O}}(1) & {\cal{O}}(q) & 0 \\ {\cal{O}}(q) & {\cal{O}}(1) & {\cal{O}}(q^2) \\ 0 & {\cal{O}}(q^2) & 0 \end{pmatrix}. \end{aligned}

      (18)

      The ratio between two nonzero squared eigenvalues R is estimated by

      \begin{align} R = {\cal{O}}(q) <<1 . \end{align}

      (19)

      This suggests that the neutrino mass eigenvalues tend to be hierarchical; hence, NH is favored. In fact, we would not obtain the allowed region within 3\sigma for IH in our numerical analysis. Therefore, we focus on NH hereafter.

      In our model, the PMNS matrix is parametrized by three mixing angles \theta_{ij} (i,j = 1,2,3; i < j) , one CP violating Dirac phase \delta_{\rm CP} , and one Majorana phase \alpha_{21} as follows:

      \begin{equation} U_{\rm PMNS} = \begin{pmatrix} c_{12} c_{13} & s_{12} c_{13} & s_{13} {\rm e}^{-{\rm i} \delta_{\rm CP}} \\ -s_{12} c_{23} - c_{12} s_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & c_{12} c_{23} - s_{12} s_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & s_{23} c_{13} \\ s_{12} s_{23} - c_{12} c_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & -c_{12} s_{23} - s_{12} c_{23} s_{13} {\rm e}^{{\rm i} \delta_{\rm CP}} & c_{23} c_{13} \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 \\ 0 & {\rm e}^{{\rm i} \frac{\alpha_{21}}{2}} & 0 \\ 0 & 0 & 1 \end{pmatrix}, \end{equation}

      (20)

      where c_{ij} and s_{ij} represent \cos \theta_{ij} and \sin \theta_{ij} , respectively. Then, these mixings are given in terms of the components of U_{\rm PMNS} as follows:

      \begin{aligned}[b]& \sin^2\theta_{13} = |(U_{\rm PMNS})_{13}|^2,\quad \sin^2\theta_{23} = \frac{|(U_{\rm PMNS})_{23}|^2}{1-|(U_{\rm PMNS})_{13}|^2},\\& \sin^2\theta_{12} = \frac{|(U_{\rm PMNS})_{12}|^2}{1-|(U_{\rm PMNS})_{13}|^2}. \end{aligned}

      (21)

      Furthermore, we compute the Jarlskog invariant J_{\rm CP} derived from PMNS matrix elements as follows:

      \begin{equation} J_{\rm CP} = {\rm{Im}} [U_{e1} U_{\mu 2} U_{e 2}^* U_{\mu 1}^*] = s_{23} c_{23} s_{12} c_{12} s_{13} c^2_{13} \sin \delta_{\rm CP}. \end{equation}

      (22)

      Majorana phase is estimated in terms of other invariant I_1 as follows:

      \begin{equation} I_1 = {\rm{Im}}[U^*_{e1} U_{e2}] = c_{12} s_{12} c_{13}^2 \sin \left( \frac{\alpha_{21}}{2} \right). \end{equation}

      (23)

      In addition, the effective mass for the neutrinoless double beta decay is expressed as:

      \begin{align} &{\rm{(NH)}}:\ \langle m_{ee}\rangle = \kappa|\tilde D_{\nu_2} s^2_{12} c^2_{13}{\rm e}^{{\rm i}\alpha_{21}}+\tilde D_{\nu_3} s^2_{13}{\rm e}^{-2{\rm i}\delta_{CP}}|, \end{align}

      (24)

      \begin{align} &{\rm{(IH)}}:\ \langle m_{ee}\rangle = \kappa|\tilde D_{\nu_1} c^2_{12} c^2_{13}+\tilde D_{\nu_2} s^2_{12} c^2_{13}{\rm e}^{{\rm i}\alpha_{21}}|, \end{align}

      (25)

      where its value could be measured by KamLAND-Zen in the future [102]. We will adopt the neutrino experimental data at the 3σ interval in Nufit 5.0 [103, 104].

      Non-unitarity:

      Here, let us briefly discuss the non-unitarity matrix U'_{\rm PMNS} . This is typically parametrized by the form

      \begin{align} U'_{\rm PMNS}\equiv \left(1-\frac12 F^\dagger F\right) U_{\rm PMNS}, \end{align}

      (26)

      where F\equiv (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D} is a hermitian matrix and U'_{\rm PMNS} represents the deviation from the unitarity. Applying global constraints [105], one finds [106]

      \begin{align} |FF^\dagger|\le \left[\begin{array}{ccc} 2.5\times 10^{-3} & 2.4\times 10^{-5} & 2.7\times 10^{-3} \\ 2.4\times 10^{-5} & 4.0\times 10^{-4} & 1.2\times 10^{-3} \\ 2.7\times 10^{-3} & 1.2\times 10^{-3} & 5.6\times 10^{-3} \\ \end{array}\right]. \end{align}

      (27)

      In our case, F\equiv (M_{NS}^T M_{NS})^{-1}M_{NS}^T M_{\rm D} = \dfrac{v_1}{v_\varphi} (\tilde M_{NS}^T\tilde M_{NS})^{-1}\times \tilde M_{NS}^T \tilde M_{\rm D} . Because v_\varphi is freely considered to be large (while v_1 = {\cal{O}}(100) GeV at most), we easily satisfy this bound.

    III.   NUMERICAL ANALYSIS
    • In this section, we perform a numerical \Delta \chi^2 analysis to determine the parameters that satisfy the neutrino oscillation data and non-unitarity constraint; in addition, we present our predictions, where we adopt best fit values of charged-lepton masses. Here, we focus on NH because IH is disfavored by the analytical estimation, as can be observed in the previous section.

      In our numerical analysis, we randomly scan free parameters in following ranges

      \begin{align} & \{ \alpha_{\rm D}, |\beta_{\rm D}|, |\gamma_{\rm D}|, \alpha_{NS}, \beta_{NS} \} \in [10^{-5},1.0], \quad t_\beta\in[10,100], \end{align}

      (28)

      where τ runs over the fundamental region, t_\beta\equiv v_1/v_2 , and \sqrt{v_1^2+v_2^2} = 246 GeV. We perform numerical analysis under the above regions. Figure 1 illustrates the correlation between the real and imaginary parts of τ, where the blue, green, and red points are allowed within 2, 3, and 5 of the \Delta \chi^2 analysis, respectively, for five accurately known observables \Delta m^2_{\rm{atm}},\Delta m^2_{\rm{sol}},s_{12}^2,s_{23}^2,s_{13}^2 in Nufit 5.0 [103, 104]. The real part runs through the entire range, but the imaginary part is localized at the [1.35-1.5] .

      Figure 1.  (color online) Allowed region of modulus τ, where the blue, green, and red points are allowed within 2, 3, and 5 of the \Delta \chi^2 analysis, respectively.

      Figure 2 illustrates the correlation between the Majorana phase \alpha_{21} and Dirac CP phase \delta_{\rm CP} . The legend is the same as in the case of Fig. 1. We identify a clear feature in the region where \alpha_{21} is within \{0^\circ - 120^\circ, 240^\circ - 310^\circ \} , while \delta_{\rm CP} is within \{70^\circ - 290^\circ \} .

      Figure 2.  (color online) Correlation between a Majorana phase \alpha_{21} and \delta_{\rm CP} , where the legend is the same as the case of Fig. 1.

      The upper figures in Fig. 3 present the correlation between the effective mass for the neutrinoless double beta decay \langle m_{ee} \rangle and the sum of neutrino masses \sum m_i in the [eV] unit, where the legend is the same as the case of Fig. 1. We determine that \langle m_{ee} \rangle is allowed within [0.001-0.004] eV. However, \sum m_i is restricted to be approximately 0.06 eV, which would be a sharp prediction of this model.

      Figure 3.  (color online) Predicted correlation between the effective mass for the neutrinoless double beta decay \langle m_{ee} \rangle and the sum of neutrino masses \sum m_i in [eV], where the legend is the same as the case of Fig. 1.

      Finally, we present a benchmark in Table 2, which is selected, such that \sqrt{\Delta \chi^2} is minimum. The mass matrices for the dimensionless neutrino and charged-lepton are determined as

      parameter valve
      τ 0.113762 + 1.43906 {\rm i}
      t_\beta 98.6
      [\alpha_\ell, \gamma_{\ell}, \beta_{\ell}] [1.72\times 10^{-5},\ 6.15\times 10^{-4},\ 8.06\times 10^{-4}]
      [\alpha_{\rm D}, \alpha_{NS}, \beta_{NS}] [-0.0112907,\ -0.00078512,\ 0.0203379]
      [\beta_{\rm D}, \gamma_{\rm D}] [-0.01432- 0.00360 {\rm i},-1.90\times 10^{-5} - 2.34\times10^{-6} {\rm i}]
      \Delta m^2_{\rm{atm}} 2.53\times10^{-3} {\rm{eV}}^2
      \Delta m^2_{\rm{sol}} 7.48\times10^{-5} {\rm{eV}}^2
      \sin^2\theta_{12} 0.289
      \sin^2\theta_{23} 0.565
      \sin^2\theta_{13} 0.02207
      [\delta_{\rm CP},\ \alpha_{21}] [248^\circ,\, 80.3^\circ]
      \sum m_i 58.5 meV
      \langle m_{ee} \rangle 1.65 meV
      \sqrt{\Delta\chi^2} 1.40

      Table 2.  Benchmark point of our input parameters and observables, which is selected, such that \sqrt{\Delta \chi^2} is minimum.

      \begin{align} \tilde m_\nu& = \left[\begin{array}{ccc} 37.1331- 17.9304 {\rm i} & -16.2576 + 55.2032 {\rm i} & 35.2235+ 24.3462 {\rm i} \\ -16.2576 + 55.2032 {\rm i} &-40.0223 -38.1947 {\rm i} & -55.9229 + 44.0753 {\rm i} \\ 35.2235+ 24.3462 {\rm i} & -55.9229 + 44.0753 {\rm i} &-1.23231 + 74.8642 {\rm i} \\ \end{array}\right], \end{align}

      (29)

      \begin{align} M_\ell& = \left[\begin{array}{ccc} 0.00232 + 0.000467 {\rm i} & -0.00111 + 0.00314 {\rm i} & -0.00453 + 0.00274 {\rm i} \\ -0.282 + 0.102 {\rm i} & -0.0524 - 0.469 {\rm i}& 0.214 - 0.559 {\rm i} \\ 0.336 - 0.277 {\rm i} & 0.362 + 0.706 {\rm i}& -0.160 + 1.23 {\rm i} \\ \end{array}\right]. \end{align}

      (30)

    IV.   SUMMARY AND DISCUSSION
    • We studied a linear seesaw model with a field content as minimum as possible, introducing a modular S_4 using U(1)_{B-L} symmetries. Owing to the rank two neutrino mass matrix, we obtained a vanishing neutrino mass eigenvalue. Furthermore, only the normal mass hierarchy of neutrinos was favored via the modular S_4 symmetry. In our numerical \Delta \chi^2 analysis, we determined a rather sharp prediction on the sum of neutrino masses to be approximately 60 meV. The imaginary part of τ was restricted at 1.35-1.5, while the real part ran through the entire range in the fundamental region. Other remarks are presented below:

      1. \alpha_{21} was within \{0^\circ - 120^\circ, 240^\circ - 310^\circ \} , while \delta_{\rm CP} was within \{70^\circ - 290^\circ \} .

      2. \langle m_{ee} \rangle was allowed by \{0.001 - 0.004\} eV.

      Therefore our model indicates several predictions in the neutrino sector that was triggered by the minimal structure with the S_4 modular symmetry.

    ACKNOWLEDGMENTS
    • This research was supported by the Korean Local Governments - Gyeongsangbuk-do Province and Pohang City (H.O.). H. O. is sincerely grateful to the KIAS member.

    APPENDIX A
    • Here, we review some properties of the modular S_4 symmetry. In general, the modular group \bar\Gamma is the group of linear fractional transformation γ acting on the modulus τ, which belongs to the upper-half complex plane and transforms as

      \begin{equation} \tau \longrightarrow \gamma\tau = \frac{a\tau + b}{c \tau + d}, \end{equation}\tag{A1}

      {\rm{where}}\; \; a,b,c,d \in \mathbb{Z}\; \; {\rm{and}}\; \; ad-bc = 1, \; \; {\rm{Im}} [\tau]>0 \; .

      This is isomorphic to PSL(2,\mathbb{Z}) = SL(2,\mathbb{Z})/\{I,-I\} transformation. Then modular transformation is generated by two transformations S and T defined as:

      \begin{eqnarray} S:\tau \longrightarrow -\frac{1}{\tau}\ , \qquad T:\tau \longrightarrow \tau + 1\ , \end{eqnarray} \tag{A2}

      and they satisfy the following algebraic relations,

      \begin{equation} S^2 = \mathbb{I}\ , \qquad (ST)^3 = \mathbb{I}\ . \end{equation} \tag{A3}

      Here, we introduce the series of groups \Gamma(N)\; (N = 1,2,3,\dots) which are defined by

      \begin{align} \begin{aligned} \Gamma(N) = \left \{ \begin{pmatrix} a & b \\ c & d \end{pmatrix} \in SL(2,\mathbb{Z})\; , \; \; \begin{pmatrix} a & b \\ c & d \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} \; \; ({\rm{mod}} N) \right \} \end{aligned}, \end{align}\tag{A4}

      and we define \bar\Gamma(2)\equiv \Gamma(2)/\{I,-I\} for N = 2 . Because the element -I does not belong to the \Gamma(N) for N>2 case, we have \bar\Gamma(N) = \Gamma(N) , which are the infinite normal subgroup of \bar \Gamma known as principal congruence subgroups. Hence, we obtain finite modular groups as the quotient groups defined by \Gamma_N\equiv \bar \Gamma/\bar \Gamma(N) . For these finite groups \Gamma_N , T^N = \mathbb{I} is imposed, and the groups \Gamma_N with N = 2,3,4,5 are isomorphic to S_3 , A_4 , S_4 and A_5 , respectively [11].

      Modular forms of level N are holomorphic functions f(\tau) transformed under the action of \Gamma(N) , which is expressed as:

      \begin{equation} f(\gamma\tau) = (c\tau+d)^k f(\tau)\; , \; \; \gamma \in \Gamma(N)\; , \end{equation}\tag{A5}

      where k denotes the reputed modular weight.

      Here, we discuss the modular symmetric theory framework without imposing supersymmetry explicitly, considering the S_4 ( N = 4 ) modular group. Under the modular transformation in Eq. (A1), a field \phi^{(I)} is also transformed as

      \begin{equation} \phi^{(I)} \to (c\tau+d)^{-k_I}\rho^{(I)}(\gamma)\phi^{(I)}, \end{equation} \tag{A6}

      where -k_I represents the modular weight and \rho^{(I)}(\gamma) denotes a unitary representation matrix of \gamma\in\Gamma(4) . Hence, a Lagrangian such as the Yukawa terms can be invariant if the sum of modular weight from fields and the modular form in the corresponding term is zero (also invariant under S_4 and gauge symmetry).

      The kinetic and quadratic terms of the scalar fields can be expressed as:

      \begin{equation} \sum\limits_I\frac{|\partial_\mu\phi^{(I)}|^2}{(-i\tau+i\bar{\tau})^{k_I}} \; , \quad \sum\limits_I\frac{|\phi^{(I)}|^2}{(-i\tau+i\bar{\tau})^{k_I}} \; , \end{equation}\tag{A7}

      which is invariant under the modular transformation, and the overall factor is eventually absorbed by a field redefinition consistently. Therefore, the Lagrangian associated with these terms should be invariant under the modular symmetry.

      The basis of modular forms with weight 2, Y^{(2)}_{\rm{2}} = (y_{1},y_{2}) and Y^{(2)}_{\rm{3'}} = (y_{3},y_{4},y_{5}) , which are transformed as a doublet and triplet under S_4 , are determined in terms of the Dedekind eta-function \eta(\tau) and its derivative [10]:

      \begin{aligned}[b] y_1(\tau) = & \frac{\rm i}{8} \left( 8 \frac{\eta'\left(\tau+\dfrac12\right)}{\eta\left(\tau+\dfrac12\right)} + 32 \frac{\eta'(4\tau)}{\eta(4\tau)} - \frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)}- \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)}\right. \\&\left.- \frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} - \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right), \\ y_2(\tau) = & \frac{ {\rm i} \sqrt3}{8} \left(\frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)} - \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)} + \frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} - \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right), \\ y_3(\tau) = & {\rm i} \left(\frac{\eta'\left(\tau+\dfrac12\right)}{\eta\left(\tau+\dfrac12\right)} - 4 \frac{\eta'(4\tau)}{\eta(4\tau)} \right), \\ y_4(\tau) = & \frac{\rm i}{4\sqrt2} \left( - \frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)}+{\rm i} \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)} +\frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} - {\rm i} \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right), \end{aligned}

      \begin{aligned}[b] y_5(\tau) = \frac{\rm i}{4\sqrt2} \left( - \frac{\eta'\left(\dfrac\tau4\right)}{\eta\left(\dfrac\tau4\right)} - i \frac{\eta'\left(\dfrac{\tau+1}4\right)}{\eta\left(\dfrac{\tau+1}4\right)} + \frac{\eta'\left(\dfrac{\tau+2}4\right)}{\eta\left(\dfrac{\tau+2}4\right)} +{\rm i} \frac{\eta'\left(\dfrac{\tau+3}4\right)}{\eta\left(\dfrac{\tau+3}4\right)} \right) . \end{aligned} \tag{A8}

      y_i 's can be expanded in terms of q as follows:

      \begin{align} y_1 & = - 3 \pi \left(\frac{b_1}{8} + 3 b_5\right),\; \ y_2 = 3 \sqrt{3} \pi b_3,\; \ y_3 = - \pi \left( -\frac{b_1}{4} + 2 b_5\right), \\ y_4 & = - \pi \sqrt{2} b_2, \;\ y_5 = - 4\pi \sqrt{2} b_4, \end{align}\tag{A9}

      where b_i are given by

      \begin{align} b_1 \sim 1, \;\;\ b_2 \sim q,\;\; \ b_3 \sim q^2,\; \;\ b_4 \sim 0,\; \;\ b_5 \sim 0, \end{align} \tag{A10}

      with q = \exp(2 \pi i \tau) and |q| \ll 1 [96].

      Subsequently, Yukawas with higher weights are constructed by the multiplication rules of S_4 , and the following couplings can be obtained:

      \begin{aligned}[b] &Y^{(4)}_{\bf2} = \left[\begin{array}{c} y_2^2-y_1^2 \\ 2 y_1 y_2 \\ \end{array}\right], \;\;\; Y^{(4)}_{\bf3} = \left[\begin{array}{c} -2y_2y_3 \\ \sqrt3 y_1 y_5+ y_2 y_4 \\ \sqrt3 y_1 y_4+ y_2 y_5 \\ \end{array}\right],\\& Y^{(4)}_{\bf3'} = \left[\begin{array}{c} 2y_1y_3 \\ \sqrt3 y_2 y_5 - y_1 y_4 \\ \sqrt3 y_2 y_4- y_1 y_5 \\ \end{array}\right] . \end{aligned} \tag{A11}

Reference (106)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return