Towards unification of quark and lepton flavors in $A_4$ modular invariance

We study quark and lepton mass matrices in the $A_4$ modular symmetry towards the unification of the quark and lepton flavors. We adopt modular forms of weights $2$ and $6$ for quarks and charged leptons, while we use modular forms of weight $4$ for the neutrino mass matrix which is generated by the Weinberg operator. We obtain the successful quark mass matrices, in which the down-type quark mass matrix is constructed by modular forms of weight $2$, but the up-type quark mass matrix is constructed by modular forms of weight $6$. Two regions of $\tau$ are consistent with observed CKM matrix elements. The one is close to $\tau=i$ and the other is in the larger ${\rm Im }[\tau]$. On the other hand, lepton mass matrices work well only at nearby $\tau=i$, which overlaps with the one of the quark sector, for the normal hierarchy of neutrino masses. In the common $\tau$ region for quarks and leptons, the predicted sum of neutrino masses is $87$--$120$meV taking account of its cosmological bound. Since both the Dirac CP phase $\delta_{CP}^\ell$ and $\sin^2\theta_{23}$ are correlated with the sum of neutrino masses, improving its cosmological bound provides crucial tests for our scheme as well as the precise measurement of $\sin^2\theta_{23}$ and $\delta_{CP}^\ell$. The effective neutrino mass of the $0\nu\beta\beta$ decay is $\langle m_{ee}\rangle=15$--$31$\,meV. It is remarked that the modulus $\tau$ is fixed at nearby $\tau=i$ in the fundamental domain of SL$(2,Z)$, which suggests the residual symmetry $Z_2$ in the quark and lepton mass matrices. The inverted hierarchy of neutrino masses is excluded by the cosmological bound of the sum of neutrino masses.


Introduction
The standard model (SM) was well established by the discovery of the Higgs boson. However, the flavor theory of quarks and leptons is still unknown. In order to understand the origin of the flavor structure, many works have appeared by using the discrete groups for flavors. In the early models of quark masses and mixing angles, the S 3 symmetry was used [1,2]. It was also discussed to understand the large mixing angle [3] in the oscillation of atmospheric neutrinos [4]. For the last twenty years, the discrete symmetries of flavors have been developed, that is motivated by the precise observation of flavor mixing angles of leptons [5][6][7][8][9][10][11][12][13].
Many models have been proposed by using the non-Abelian discrete groups S 3 , A 4 , S 4 , A 5 and other groups with larger orders to explain the large neutrino mixing angles. Among them, the A 4 flavor model is an attractive one because the A 4 group is the minimal one including a triplet irreducible representation, which allows for a natural explanation of the existence of three families of leptons [14][15][16][17][18][19][20]. However, the variety of models is so wide that it is difficult to obtain a clear evidence of the A 4 flavor symmetry.
Recently, a new approach to the lepton flavor problem appeared based on the invariance under the modular group [21], where the model of the finite modular group Γ 3 ≃ A 4 has been presented. This work inspired further studies of the modular invariance approach to the lepton flavor problem. It should be emphasized that there is a significant difference between the models based on the A 4 modular symmetry and those based on the usual non-Abelian discrete A 4 flavor symmetry. Yukawa couplings transform non-trivially under the modular group and are written in terms of modular forms which are holomorphic functions of the modulus τ .
The modular group includes the finite groups S 3 , A 4 , S 4 , and A 5 [22]. Therefore, an interesting framework for the construction of flavor models has been put forward based on the Γ 3 ≃ A 4 modular group [21], and further, based on Γ 2 ≃ S 3 [23]. The proposed flavor models with modular symmetries Γ 4 ≃ S 4 [24] and Γ 5 ≃ A 5 [25] have also stimulated studies of flavor structures of quarks and leptons. Phenomenological discussions of the neutrino flavor mixing have been done based on A 4 [26,27], S 4 [28,29], A 5 [30], and T ′ [31] modular groups, respectively. In particular, the comprehensive analysis of the A 4 modular group has provided a distinct prediction of the neutrino mixing angles and the CP violating phase [27]. The A 4 modular symmetry has been also applied to the SU(5) grand unified theory (GUT) of quarks and leptons [32], while the residual symmetry of the A 4 modular symmetry has been investigated phenomenologically [33]. Furthermore, modular forms for ∆(96) and ∆(384) were constructed [34], and the extension of the traditional flavor group is discussed with modular symmetries [35]. Moreover, multiple modular symmetries are proposed as the origin of flavor [36]. The modular invariance has been also studied combining with the generalized CP symmetries for theories of flavors [37]. The quark mass matrix has been discussed in the S 3 and A 4 modular symmetries as well [38][39][40]. Besides mass matrices of quarks and leptons, related topics have been discussed in the baryon number violation [38], the dark matter [41,42] and the modular symmetry anomaly [43].
In this work, we study both quark and lepton mass matrices in the A 4 modular symmetry. If flavors of quarks and leptons are originated from a same two-dimensional compact space, quarks and leptons have the same flavor symmetry and the same value of the modulus τ . Therefore, it is challenging to reproduce observed three Cabibbo-Kobayashi-Maskawa (CKM) mixing angles and the CP violating phase while observed large mixing angles in the lepton sector within the framework of the A 4 modular invariance with the common τ . This work provides a new aspect in order to discover the unification theory of the quark and lepton flavors. We have already discussed the quark mass matrices in terms of A 4 modular forms of weight 2 [40]. It has been found that quark mass matrices of A 4 do not work unless Higgs sector is extended by A 4 triplet representations. In this paper, we propose to adopt modular forms of weight 6 in addition to modular ones of weight 2 for quarks with the Higgs sector of SM. We use modular forms of weight 4 for the neutrino mass matrix, which is generated by the Weinberg operator. The common value of the modulus τ is successfully obtained by using observed four CKM matrix elements and three lepton mixing angles. We also predict the CP violating Dirac phase of leptons, which is expected to be observed at T2K and NOνA experiments [44,45], and the sum of neutrino masses 1 .
The paper is organized as follows. In section 2, we give a brief review on the modular symmetry and modular forms of weights 2, 4 and 6. In section 3, we present the model for quark mass matrices in the A 4 modular symmetry. In section 4, we show numerical results for the CKM matrix. In section 5, we discuss the lepton mass matrices and present numerical results. Section 6 is devoted to a summary. In Appendix A, the tensor product of the A 4 group is presented. In Appendix B, we present how to obtain Dirac CP phase, Majorana phases and effective mass of the 0νββ decay.
2 Modular group and modular forms of weights 2, 4, 6 The modular groupΓ is the group of linear fractional transformations γ acting on the modulus τ , belonging to the upper-half complex plane as: which is isomorphic to P SL(2, Z) = SL(2, Z)/{I, −I} transformation. This modular transformation is generated by S and T , which satisfy the following algebraic relations, We introduce the series of groups Γ(N) (N = 1, 2, 3, . . . ), called principal congruence subgroups, defined by For N = 2, we defineΓ(2) ≡ Γ(2)/{I, −I}. Since the element −I does not belong to Γ(N) for N > 2, we haveΓ(N) = Γ(N). The quotient groups defined as Γ N ≡Γ/Γ(N) are finite modular groups. In this finite groups Γ N , T N = I is imposed. The groups Γ N with N = 2, 3, 4, 5 are isomorphic to S 3 , A 4 , S 4 and A 5 , respectively [22]. Modular forms of level N are holomorphic functions f (τ ) transforming under Γ(N) as: where k is the so-called as the modular weight. Superstring theory on the torus T 2 or orbifold T 2 /Z N has the modular symmetry [47][48][49][50][51][52]. Its low energy effective field theory is described in terms of supergravity theory, and string-derived supergravity theory has also the modular symmetry. Under the modular transformation of Eq.(1), chiral superfields φ (I) transform as [53], where −k I is the modular weight and ρ (I) (γ) denotes a unitary representation matrix of γ ∈Γ.
In the present article we study global supersymmetric models, e.g., minimal supersymmetric extensions of the Standard Model (MSSM). The superpotential which is built from matter fields and modular forms is assumed to be modular invariant, i.e., to have a vanishing modular weight. For given modular forms this can be achieved by assigning appropriate weights to the matter superfields.
The kinetic terms are derived from a Kähler potential. The Kähler potential of chiral matter fields φ (I) with the modular weight −k I is given simply by where the superfield and its scalar component are denoted by the same letter, andτ = τ * after taking the vacuum expectation value (VEV) 2 . Therefore, the canonical form of the kinetic terms is obtained by changing the normalization of parameters [27]. For Γ 3 ≃ A 4 , the dimension of the linear space M k (Γ(3)) of modular forms of weight k is k + 1 [55][56][57], i.e., there are three linearly independent modular forms of the lowest non-trivial weight 2. These forms have been explicitly obtained [21] in terms of the Dedekind eta-function η(τ ): where η(τ ) is a so called modular form of weight 1/2. In what follows we will use the following basis of the A 4 generators S and T in the triplet representation: where ω = exp(i 2 3 π) . The modular forms of weight 2 transforming as a triplet of A 4 can be written in terms of η(τ ) and its derivative [21]: The overall coefficient in Eq. (10) is one possible choice. It cannot be uniquely determined. The triplet modular forms of weight 2 have the following q-expansions: They satisfy also the constraint [21]: The modular forms of the higher weight, k, can be obtained by the A 4 tensor products of the modular forms with weight 2, Y By using these modular forms of weights 2, 4 and 6, we discuss quark and lepton mass matrices.

Quark mass matrices in the A 4 modular invariance
Let us consider a A 4 modular invariant flavor model for quarks. There are freedoms for the assignments of irreducible representations and modular weights to quarks and Higgs doublets. The simplest one is to assign the triplet of the A 4 group to three left-handed quarks, but three different singlets (1, 1 ′′ , 1 ′ ) of A 4 to the three right-handed quarks, (u c , c c , t c ) and (d c , s c , b c ), respectively, where the sum of weights of the left-handed and the right-handed quarks is −2. The Higgs fields are supposed to be A 4 singlets with weight 0. Then, three independent couplings appear in the superpotential of the up-type and down-type quark sectors, respectively, as follows: where Q is the left-handed A 4 triplet quarks, and H q is the Higgs doublet. The parameters α q , β q , γ q are constant coefficients. Assign the A 4 triplet Q as ((d L , u L ), (s L , c L ), (b L , t L )). By using the decomposition of the A 4 tensor product in Appendix A, the superpotentials in Eqs. (15) and (16) give the mass matrix of quarks, which is written in terms of modular forms of weight 2 as: where τ of the modular forms Y i (τ ) is omitted. The constant v q (q = u, d) is the VEV of the neutral component of the Higgs field H q . Parameters α q , β q , γ q are taken to be real without loss of generality, and they can be adjusted to the observed quark masses. The remained parameter is only the modulus τ . The numerical study of the quark mass matrix in Eq. (17) is rather easy. However, it is impossible to reproduce observed hierarchical three CKM mixing angles by fixing one complex parameter τ . In order to obtain realistic quark mass matrices, we use modular forms of weight 6 in Eq. (14). As a simple model, we take modular forms of weight 6 only for the up-type quark mass matrix while the down-type quark one is still given in terms of modular forms of weight 2 such as Eq. (17) 3 . Then, we have six independent couplings in the superpotential of the up-quark sector as: where assignments of representations and weights for MSSM fields are given in Table 1. The up-type quark mass matrix is written as: where g u1 = α ′ u /α u , g u2 = β ′ u /β u and g u3 = γ ′ u /γ u are complex parameters while α u , β u and γ u are real. On the other hand, the down-type quark mass matrix is given as: We will fix the modulus τ phenomenologically by using quark mass matrices in Eqs. (19) and (20).

Fixing τ by observed CKM
In order to obtain the left-handed flavor mixing, we calculate M † d M d and M † u M u . At first, we take a random point of τ and g ui which are scanned in the complex plane by generating random numbers. The modulus τ is scanned in the fundamental domain of the modular symmetry. In practice, the scanned range of Im[τ ] is [ √ 3/2, 2], in which the lower-cut √ 3/2 is at the cusp of the fundamental domain, and the upper-cut 2 is enough large for estimating Y i . On the other hand, Re[τ ] is scanned in the fundamental domain [−1/2, 1/2] of the modular group. Supposing |g ui | is of order one, we scan them in [0, 2] while these phases are scanned in [−π, π]. Then, parameters α q , β q , γ q (q = u, d) are given in terms of τ and g q after inputting six quark masses.
Finally, we calculate three CKM mixing angles and the CP violating phase in terms of the model parameters τ and g ui . We keep the parameter sets, in which the value of each observable is reproduced within the three times of 1σ interval of error-bars. We continue this procedure to obtain enough points for plotting allowed region.
We input quark masses in order to constrain model parameters. Since the modulus τ obtains the expectation value by the breaking of the modular invariance at the high mass scale, the quark masses are put at the GUT scale. The observed masses and CKM parameters run to the GUT scale by the renormalization group equations (RGEs). In our work, we adopt numerical values of Yukawa couplings of quarks at the GUT scale 2 × 10 16 GeV with tan β = 5 in the framework of the minimal SUSY breaking scenarios [58,59]: which give quark masses as m q = y q v H with v H = 174 GeV.
We also use the following CKM mixing angles to focus on parameter regions consistent with the experimental data at the GUT scale 2 × 10 16 GeV, where tan β = 5 is taken [58,59]: Here θ CKM ij is given in the PDG notation of the CKM matrix V CKM [60]. The CP violating phase is also given as: in the PDG notation. The errors in Eqs. (21), (22) and (23) represent 1σ interval. The CKM elements V ij at the GUT scale 2 × 10 16 GeV are given by using these angles and the phase. In our model, we have four complex parameters, τ , g u1 , g u2 and g u3 after inputting six quark masses. These eight real parameters are scanned to reproduce the observed three CKM mixing angles and the CP violating phase with three times 1σ error interval in Eqs. (22) and (23) 4 .
We have succeeded to reproduce completely four observed CKM elements in the parameter ranges of Table 2. The modulus τ is close to i, which is the fixed point of the modular symmetry.
In order to check the consistency of our quark mass matrices and the observed CKM, we show the calculated distribution on the |V ub |-δ CP plane at the GUT scale in Fig. 1. The calculated δ CP is uniformly distributed below the observed central value of |V ub | while it is almost larger than the observed central value for the upper-range of |V ub |.   We also present the distribution of CKM elements |V cb | and |V ub | at the GUT scale in Fig. 2. The magnitude of |V ub | is predicted to be in the whole region of the three times 1σ interval while the calculated |V cb | is mostly distributed in the lower-range of the three times 1σ interval. In Fig. 3, we show the plot of Re[τ ] and Im[τ ], which will be compared with the case of leptons. The allowed region of the modulus τ is close to i, but is clearly deviated from it. The modulus τ = i is the fixed point of the modular symmetry. Indeed, τ = i is invariant under the S transformation τ = −1/τ , where the subgroup Z S 2 = {I, S} of A 4 is preserved. This region of τ is discussed in connection with τ of the lepton mass matrices in the next secton.
In Table 3, we present one parameter set and calculated CKM elements, which is the best-fit point, that is, its χ 2 is minimum. The magnitudes of g qi are at most of order O(0.5). Ratios of α q /γ q and β q /γ q (q = u, d) correspond to the observed quark mass hierarchy.
We also present the mixing matrices of up-type quarks and down-type quarks for the sample of Table 3. They are given as: where V CKM = V † u V d . It is noted that these are presented in the diagonal base of the generator S, where we can see the hierarchical structure of mixing. The diagonal base of S is realized by the unitary transformation of S in Eq. (9) Then, the mixing matrix V q in the original base of S is transformed to V S V q because the quark mass matrix is transformed as V S M † q M q V † S . In conclusion, our quark mass matrices with the A 4 modular symmetry reproduce the observed CKM mixing matrix very well at nearby τ = i. This successful result encourages us to investigate the lepton flavors in the same framework. We discuss the lepton sector with the A 4 modular symmetry in Section 5.

Lepton mass matrix
The modular A 4 invariance also gives the lepton mass matrix in terms of the modulus τ which is common to both quarks and leptons if flavors of quarks and leptons are originated from a same two-dimensional compact space. We assign the A 4 representation and weight for leptons in Table 4, where the left-handed lepton doublets compose a A 4 triplet and the right-handed charged leptons are A 4 singlets. The weights of the leptons are assigned like the down-type quarks in Table 1. Then, the superpotential of the charged lepton mass term is given in terms of modular forms of weight 2, Y 3 since weights of the left-handed leptons and the right-handed charged leptons are −2 and 0, respectively. It is given as: where L is the left-handed A 4 triplet leptons. Taking (e L , µ L , τ L ) in the flavor base the charged lepton mass matrix M E is simply written as: where coefficients α e , β e and γ e are real parameters. Suppose neutrinos to be Majorana particles. By using the Weinberg operator, the superpotential of the neutrino mass term, w ν is given as: where Λ is a relevant cutoff scale. Since the left-handed lepton doublet has weight −2, the superpotential is given in terms of modular forms of weight 4, Y 1 and Y 1 ′ . By using the tensor products of A 4 , we have where Y 1 and Y (4) 1 ′ are given in Eq. (13), and g ν1 , g ν2 are complex parameters. The neutrino mass matrix is written as follows: Then, the model parameters are α e , β e , γ e , g ν1 and g ν2 . Parameters α e , β e and γ e are adjusted by the observed charged lepton masses. Therefore, the lepton mixing angles, the Dirac phase and Majorana phases are given by g ν1 and g ν2 in addition to the value of τ . Since τ is restricted in the narrow range for the quark sector as seen in Fig. 3, we can give some predictions in the lepton sector. Practically, we scan τ in |Re[τ ]| ≤ 0.5 and Im[τ ] ≤ 2 like in the analysis of quark mass matrices.

Numerical results of leptons
We input charged lepton masses in order to constrain the model parameters. We take Yukawa couplings of charged leptons at the GUT scale 2 × 10 16 GeV, where tan β = 5 is taken as well as quark Yukawa couplings [58,59]: y e = (1.97 ± 0.024) × 10 −6 , y µ = (4.16 ± 0.050) × 10 −4 , y τ = (7.07 ± 0.073) × 10 −3 , where lepton masses are given by m ℓ = y ℓ v H with v H = 174 GeV. We also input the lepton mixing angles and neutrino mass parameters which are given by NuFit 4.1 in Table 5 [61]. We investigate two possible cases of neutrino masses m i , which are the normal hierarchy (NH), m 3 > m 2 > m 1 , and the inverted hierarchy (IH), m 2 > m 1 > m 3 . Neutrino masses and the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix U PMNS [62,63] are obtained by diagonalizing M † E M E and M * ν M ν . We also investigate the effective mass for the 0νββ decay, m ee (see Appendix B) and the sum of three neutrino masses m i since it is constrained by the recent cosmological data, which is the upperbound m i ≤ 120 meV obtained at the 95% confidence level [64,66].  Let us discuss numerical results for NH of neutrino masses. Since α e /γ e and β e /γ e are obtained by the observed charged lepton masses when τ is fixed, we input charged lepton masses to reduce free parameters. Parameters g ν1 and g ν2 are constrained by four observed quantities; three mixing angles of leptons and observed mass ratio ∆m 2 sol /∆m 2 atm . In practice, the scanned ranges of Im  At first, we show the allowed region on the Re[τ ]-Im[τ ] plane in Fig. 4. Observed three mixing angles of leptons are reproduced at cyan, blue and red points. The sum of neutrino masses is consistent with the cosmological bound 120 meV at cyan points, but not at blue points. At red points, both CKM and PMNS are reproduced with the sum of neutrino masses below 120 meV. For comparison, we add green points for quark CKM of Fig. 3. Allowed points of leptons are almost in Im[τ ] ≤ 1.12 and |Re[τ ]| ≤ 0.1.
The common τ causes the feedback in the quark sector. However, the tendency of |V ub |, |V cb | and δ CP are not so changed compared with Figs. 1 and 2 in the common region of τ .
We show the allowed region on the m isin 2 θ 23 plane in Fig. 5, where colors (cyan, blue and red) of points correspond to points of τ in Fig. 4. Our prediction of the sum of neutrino masses is constrained by the cosmological bound as seen in Fig. 5. The minimal cosmological model, ΛCDM + m i , provides the upperbound m i < 120 meV [64,66] although it becomes weaker when the data are analysed in the context of extended cosmological models [60].
The red region, that is the common τ region for quarks and leptons, is constrained by the cos-    mological bound m i = 120 meV. Then, the predicted sum of neutrino masses is 87-120 meV. The cyan region is inconsistent with τ of quarks while the blue one is excluded by the cosmological bound m i = 120 meV, although both are consistent with the data of NuFIT 4.1 [61]. The calculated sin 2 θ 23 of the red region is distributed in restricted ranges. Therefore, the precise measurement of sin 2 θ 23 and improving the bound of the sum of neutrino masses provide crucial tests for our scheme.
We show the allowed region on the m i -δ ℓ CP plane in Fig. 6. In the region of red points, δ ℓ CP is predicted to be in the restricted ranges, 0 • -50 • , 80 • -100 • , 110 • -130 • , 230 • -250 • , 260 • -280 • and 310 • -360 • . If the cosmological bound for the sum of neutrino masses will be improved, for example, it is 100 meV, δ ℓ CP is predicted in the distict range. In Fig. 7, we plot δ ℓ CP versus sin 2 θ 23 in order to see their correlation. Since there is a significant correlation between them, the precise measurement of sin 2 θ 23 gives the clear prediction of δ ℓ CP . We can also predict the effective mass m ee for the 0νββ decay versus the lightest mass m 1 as seen in Fig. 8. At the red point region, we predict m ee = 15-31 meV. The predicted m ee larger than 31 meV is excluded by the sum of neutrino masses.
In Table 6, we present best-fit values of parameters and outputs, where we input the common value τ = −0.038 + 1.05 i in Table 3. In our χ 2 fit, δ ℓ CP is not included, but it is only an output because T2K and NOνA experiments have presented the best-fit values of δ ℓ CP with opposite sign each other [44,45]. Our predicted δ ℓ CP = 27.4 • in Table 6 is rather small. The systematic χ 2 fit of both quarks and leptons will be needed by including precise data of δ ℓ CP since δ ℓ CP ≃ ±90 • is also predicted in Fig.6.
We also present the mixing matrices of charged leptons and neutrinos for the best-fit sample of  Table 6: Numerical values of parameters and output of PMNS parameters at the best-fit point. Table 6. Those are given as: where U PMNS = U † ℓ U ν . They are also given in the diagonal base of the generator S in order to see the hierarchical structure of mixing like in the quark mixing matrix in Eq. (24), by using the unitary transformation of Eq. (25). It is noticed that the mixing matrix of charged leptons U ℓ is similar to the quark mixing matrices of Eq. (24). On the other hand, two large mixing angles appear in the neutrino mixing matrix U ν .
In our numerical calculations, we have not included the RGE effects in the lepton mixing angles and neutrino mass ratio ∆m 2 sol /∆m 2 atm . We suppose that those corrections are very small between the electroweak and GUT scales for NH of neutrino masses. This assumption is justified well in the case of tan β ≤ 5 unless neutrino masses are almost degenerate [26].
Finally, we discuss briefly the case of IH of neutrino masses. Indeed, there is a very small region of the common τ for quarks and leptons, which is marginal since the sum of neutrino masses is very close to the cosmological bound, 120 meV. Therefore, we omit discussions of this case.

Summary
We have studied both quark and lepton mass matrices in the A 4 modular symmetry towards the unification of quark and lepton flavors. If flavors of quarks and leptons are originated from a same two-dimensional compact space, quarks and leptons have the same flavor symmetry and the same value of the modulus τ .
For the quark sector, we have adopted modular forms of weights 2 and 6. We have presented the viable model for quark mass matrices, in which the down-type quark mass matrix is constructed by modular forms of weight 2 while the up-type quark mass matrix is constructed by modular forms of weight 6. In the lepton sector, the charged lepton mass matrix is constructed by modular forms of weight 2 while modular forms of weight 4 is used for the neutrino mass matrix, which is generated by the Weinberg operator.
The viable region close to τ = i is obtained in our quark mass matrices. Lepton mass matrices also work well at nearby τ = i, which overlaps with the one of the quark sector, for NH of neutrino masses. In the common τ region for quarks and leptons, the predicted sum of neutrino masses is 87-120meV taking account of its cosmological bound. Since both the Dirac CP phase δ ℓ CP and sin 2 θ 23 are correlated significantly with the sum of neutrino masses, improving its cosmological bound provides crucial tests for our scheme as well as the precise measurement of sin 2 θ 23 and δ ℓ CP . The effective neutrino mass of the 0νββ decay is predicted to be m ee = 15-31meV. The IH of neutrino masses is almost excluded by the cosmological bound of the sum of neutrino masses.
It is remarked that the common τ is fixed at nearby τ = i in the fundamental domain of SL(2, Z), which suggests the residual symmetry Z 2 in the quark and lepton mass matrices. Some corrections could violate the exact symmetry. It is also emphasized that the spontaneous CP violation in Type IIB string theory is possibly realized at nearby τ = i, where the moduli stabilization as well as the calculation of Yukawa couplings is performed in a controlled way [67]. Thus, our phenomenological result of the modulus τ is favored in the theoretical investigation.
There are also other invariants I 1 and I 2 associated with Majorana phases In terms of this parametrization, the effective mass for the 0νββ decay is given as follows: m ee = m 1 c 2 12 c 2 13 + m 2 s 2 12 c 2 13 e iα 21 + m 3 s 2 13 e i(α 31 −2δ ℓ CP ) .