Trimaximal Neutrino Mixing from Modular A4 Invariance with Residual Symmetries

We construct phenomenologically viable models of lepton masses and mixing based on modular $A_4$ invariance broken to residual symmetries $\mathbb{Z}^{T}_3$ or $\mathbb{Z}^{ST}_3$ and $\mathbb{Z}^S_2$ respectively in the charged lepton and neutrino sectors. In these models the neutrino mixing matrix is of trimaximal mixing form. In addition to successfully describing the charged lepton masses, neutrino mass-squared differences and the atmospheric and reactor neutrino mixing angles $\theta_{23}$ and $\theta_{13}$, these models predict the values of the lightest neutrino mass (i.e., the absolute neutrino mass scale), of the Dirac and Majorana CP violation(CPV) phases, as well as the existence of specific correlations between i) the values of the solar neutrino mixing angle $\theta_{12}$ and the angle $\theta_{13}$ (which determines $\theta_{12}$), ii) the values of the Dirac CPV phase $\delta$ and of the angles $\theta_{23}$ and $\theta_{13}$, iii) the sum of the neutrino masses and $\theta_{23}$, and iv) between the two Majorana phases.


Introduction
Understanding the origin of the flavour structure of quarks and leptons remains one of the outstanding problems in particle physics. The pattern of two large and one small neutrino (lepton) mixing angles, revealed by the data obtained in neutrino oscillation experiments (see, e.g., [1]), provides an important clue in the investigations of the lepton flavour problem, suggesting the existence of flavour symmetry in the lepton sector. The results of the recent global analyses of the neutrino oscillation data show also that a neutrino mass spectrum with normal ordering (NO) is favoured over the spectrum with inverted ordering (IO), as well as a preference for a value of the Dirac CP violation (CPV) phase δ close to 3π/2 (see, e.g., [2]).
The observed 3-neutrino mixing pattern can naturally be explained by extending the Standard Theory (ST) with a flavour symmetry associated with a non-Abelian discrete symmetry group. Models based on S 3 , A 4 , S 4 , A 5 and other groups of larger orders have been proposed and extensively studied (see, e.g., [3][4][5][6][7][8][9]). In particular, the A 4 flavour model attracted considerable interest because the A 4 group is the minimal one including a triplet unitary irreducible representation, which allows for a natural explanation of the existence of three families of leptons [10][11][12][13][14][15]. In all models based on non-Abelian discrete flavour symmetry, the flavour symmetry must be broken in order to reproduce the measured values of the neutrino mixing angles. This is achieved by introducing typically a large number of ST gauge singlet scalars -the so-called "flavons" -in the Lagrangian of the theory, which have to develope a set of particularly aligned vacuum expectation values (VEVs). Arranging for such an alignment requires the construction of rather elaborate scalar potentials.
An attractive approach to the lepton flavour problem, based on the invariance under the modular group, has been proposed in Ref. [16], where also models of the finite modular group Γ 3 A 4 have been presented. Although the models constructed in Ref. [16] are not realistic and make use of a minimal set of flavon fields, this work inspired further studies of the modular invariance approach to the lepton flavour problem. The modular group includes S 3 , A 4 , S 4 , and A 5 as its principal congruence subgroups, Γ 2 S 3 , Γ 3 A 4 , Γ 4 S 4 and Γ 5 A 5 [17]. However, there is a significant difference between the models based on the modular S 3 , A 4 , S 4 etc. symmetry and those based on the usual non-Abelian discrete S 3 , A 4 , S 4 etc. flavour symmetry. The constants of a theory based on the finite modular symmetry, such as Yukawa couplings and, e.g., the right-handed neutrino mass matrix in type I seesaw scenario, also transform non-trivially under the modular symmetry and are written in terms of modular forms which are holomorphic functions of a complex scalar field -the modulus τ . At the same time the modular forms transform under the usual non-Abelian discrete flavour symmetries. In the most economical versions of the models with modular symmetry, the VEV of the modulus τ is the only source of symmetry breaking without the need of flavon fields.
In Ref. [18] a realistic model with modular Γ 2 S 3 symmetry was built with the help of a minimal set of flavon fields. A realistic model of the charged lepton and neutrino masses and of neutrino mixing without flavons, in which the modular Γ 4 S 4 symmetry was used, was constructed in [19]. Subsequently, lepton flavour models without flavons, based on the modular symmetry Γ 3 A 4 was proposed in Refs. [20,21]. A comprehensive investigation of the simplest viable models of lepton masses and mixing, based on the modular S 4 symmetry, was performed in Ref. [22]. Necessary ingredients for constructing flavour models based, in particular, on the modular symmetries ∆(96) and ∆(384) have been obtained in [23], while for models based on A 5 symmetry they have been derived in [24].
If one of the subgroups of the considered finite modular group is preserved, this residual symmetry fixes τ to a specific value (see, e.g., [22]). Phenomenologically viable models based on the modular S 4 and A 5 symmetries, broken respectively to residual Z 3 and Z 5 symmetries in the charged lepton sector and to a Z 2 symmetry in the neutrino sector, were presented in Refs. [22,24]. So far, apart from these two studies, the implications of residual symmetries have been investigated only in the framework of the usual non-Abelian discrete symmetry approach to the lepton (and quark) flavour problem. It has been shown that they lead, in particular, to specific experimentally testable correlations between the values of some of the neutrino mixing angles and/or between the values of the neutrino mixing angles and of the Dirac CP violation phase in the neutrino mixing [25][26][27][28][29][30][31].
In the present article we construct phenomenologically viable models of lepton masses and mixing based on residual symmetries resulting from the breaking of the A 4 modular symmetry. It is found that the weight 4 modular forms are required to obtain charged lepton and neutrino mass matrices leading to lepton masses and mixing which are consistent with the experimental data on neutrino oscillations. We also find that in these models the PMNS matrix [32][33][34] is predicted to be of the trimaximal mixing form [35,36].
The paper is organized as follows. In section 2, we give a brief review on the modular symmetry. In section 3, we discuss the residual symmetries of A 4 and their modular forms. In section 4, we present the lepton mass matrices in the residual symmetry. In section 5, we present models and their numerical results. Section 6 is devoted to a summary. Appendix A shows the relevant multiplication rules of the A 4 group.

Modular A 4 Group and Modular Forms of Level 3
The modular group Γ is the group of linear fractional transformations γ acting on the complex variable τ belonging to the upper-half complex plane as follows: The group Γ is generated by two transformations S and T satisfying where I is the identity element. Representing S and T as The modular group Γ is isomorphic to the projective special linear group P SL(2, Z) = SL(2, Z)/Z 2 , where SL(2, Z) is the special linear group of 2 × 2 matrices with integer elements and unit determinant, and Z 2 = {I, −I} is its centre (I being the identity element). The special linear group SL(2, Z) = Γ(1) ≡ Γ contains a series of infinite normal subgroups Γ(N ), N = 1, 2, 3, . . . : The quotient groups Γ N ≡ Γ/Γ(N ) are called finite modular groups. For N ≤ 5, these groups are isomorphic to non-Abelian discrete groups widely used in flavour model building (see, e.g., [17]): Γ 2 S 3 , Γ 3 A 4 , Γ 4 S 4 and Γ 5 A 5 . We will be interested in the finite modular group Γ 3 A 4 .
Modular forms of weight k and level N are holomorphic functions f (τ ) transforming under the action of Γ(N ) in the following way: (2.6) Here k is even and non-negative, and N is natural. Modular forms of weight k and level N span a linear space of finite dimension. The dimension of the linear space of modular forms of weight k and level 3, M k (Γ 3 A 4 ), is k + 1. There exists a basis in this space such that a multiplet of modular forms f i (τ ) transforms according to a unitary representation ρ of the finite group Γ N : In the case of N = 3 of interest, the three linear independent weight 2 modular forms form a triplet of A 4 [16]. These forms have been explicitly obtained [16] in terms of the Dedekind eta-function η(τ ): where q = e 2πiτ . In what follows we will use the following basis of the A 4 generators S and T in the triplet representation: where ω = e i 2 3 π . The modular forms (Y 3 ) transforming as a triplet of A 4 can be written in terms of η(τ ) and its derivative [16]: The overall coefficient in eq. (2.10) is one possible choice; it cannot be uniquely determined. The triplet modular forms Y 1,2,3 have the following q-expansions: (2.11) They satisfy also the constraint [16]:

Residual Symmetries of A 4 and Modular Forms
Residual symmetries arise whenever the VEV of the modulus τ breaks the modular group Γ only partially, i.e., the little group (stabiliser) of τ is non-trivial. Residual symmetries have been investigated in the case of modular S 4 invariance in [22], and of A 5 invariance in [24], where viable models of lepton masses and mixing have also been constructed. In the present work we consider models of lepton flavour based on the residual symmetries of the modular A 4 invariance. There are only 2 inequivalent finite points with non-trivial little groups ofΓ, namely, [22]. The first point is the left cusp in the fundamental domain of the modular group, which is invariant under the ST transformation It is possible to calculate the values of the A 4 triplet modular forms of weight 2 at the symmetry points τ L , τ C and τ T . The results are reported in Table 1 in which the values of the modular forms at τ = τ R are also given, to be compared with those at the other two points.
As we have noted, the dimension of the linear space M k (Γ 3 A 4 ) of modular forms of weight k and level 3 is k + 1. The modular forms of weights higher than 2 can be obtained from the modular forms of weight 2. They transform according to certain irreducible representations of the A 4 group. Indeed, for weight 4 we have 5 independent modular forms, which are constructed by the weight 2 modular forms through the tensor product of 3 × 3 (see Appendix A). We obtain one triplet 3 and two singlets 1, 1 , while the third singlet 1" vanishes: where the vanishing Y (4) 1" is due to the condition in Eq. (2.12). Using Eq. (3.2) we can calculate the values of the modular forms of weight 4, transforming as 3 and {1, 1 }, at the symmetry points τ L , τ C and τ T . We show the results also in Table 1.

Lepton Mass Matrices with Residual Symmetry
We will consider next modular invariant lepton flavour models with the A 4 symmetry, assuming that the massive neutrinos are Majorana particles and that the neutrino masses originate from the Weinberg dimension 5 operator. There is a certain freedom for the assignments of irreducible representations and modular weights to leptons. We suppose that three left-handed (LH) lepton doublets form a triplet of the A 4 group. The Higgs doublets are supposed to be zero weight singlets of A 4 . The generic assignments of representations and modular weights k I to the MSSM fields are presented in Table 2. In order to construct models with minimal number of parameters, we introduce no flavons. For the charged leptons, we assign the three right-handed (RH) charged lepton fields for three different singlet representations of A 4 , (1, 1 , 1 ). Therefore, there are three independent coupling constants in the superpotential of the charged lepton sector. These coupling constants can be adjusted to the observed charged lepton masses. Since there are three singlet irreducible representations in the A 4 group, there are six cases for the assignment of the three RH charged lepton fields. However, this ambiguity does not affect the matrix which acts on the LH charged lepton fields and enters into the expression for the PMNS matrix. Thus, effectively we have the following unique form for the superpotential: where the sums of the modular weights should vanish. The parameters α, β, γ and Λ are constant coefficients.

Charged Lepton Mass Matrix with Residual Symmetry
By using the decomposition of the A 4 tensor products given in Appendix A, the superpotential in Eq. (4.1) leads to a mass matrix of charged leptons, which is written in terms of modular forms of A 4 triplet with weight k: where v d ≡ H 0 d . Without loss of generality the coefficients α, β, and γ can be made real positive by rephasing the RH charged lepton fields.
We will discuss next the charged lepton mass matrix at the specific points of τ = τ L , τ R τ C , τ T in the case of weight k = 2. At τ = τ L , the matrix M † E M E , which is relevant for the left-handed mixing, is given as: It is easily noticed that this matrix commutes with ST , which is guranteed by the residual Both matrices M † E M E and ST are diagonalized by the unitary matrix U E : where U E is independent of parameters α, β, γ.
On the other hand, at τ = τ R , we have: (4.7) The matrix M † E M E in Eq. (4.7) commutes with (4.8) The fact that M † E M E and T S commute is a consequence of the residual symmetry Z T S 3 at τ = τ R . The matrix M † E M E and ST is diagonalized by the unitary matrix: At τ = τ C , the determinant of M E vanishes. Indeed, this mass matrix leads to a massless charged lepton, and thus cannot be used for model building.
Finally, at τ = τ T we obtain the real diagonal matrix: In the case of modular forms of weight 4 we can obtain a charged lepton mass matrix in which the modular forms transforming as 1 and 1 do not contribute. As seen in Table  1, the weight 4 triplet modular forms coincide with weight 2 ones at τ = τ L , τ R . Indeed, M † E M E is obtained by replacing parameters (α, β, γ) of the mass matrices in Eqs. (4.4) and (4.7) with (γ, α, β), respectively. Therefore, the mixing matrices in Eqs. (4.6) and (4.9) are the same.
At τ = τ C , the charged lepton mass matrix is of rank one, i.e., two massless charged leptons appear since the triplet modular forms are proportional to (1, 1, 1). At τ = τ T , the charged lepton mass matrix is equal to the diagonal one given in Eq.(4.10) since the triplet weight 4 modular forms coincide with the weight 2 modular forms.

Neutrino Mass Matrix (Weinberg Operator)
The neutrino mass matrix is written in terms of A 4 triplet modular forms of weight k by using the supertoptential in Eq. (4.2): where v u ≡ H 0 u . In the case of weight 2 modular forms it is easily checked that two lightest neutrino masses are degenerate at τ = τ L , τ R , while the determinant of M ν vanishes at τ = τ C . In the latter case one neutrino is massless and two neutrino masses are degenerate. The two lightest neutrino masses are degenerate also at τ = τ T . It may be helpful to add a comment: these degeneracies of neutrino masses still hold even if we use the seesaw mechanism by introducing the three right-handed neutrino fields as A 4 triplet. Thus, the realistic neutrino mass matrix is not obtained as far as we take weight 2 modular forms at τ = τ L , τ R , τ C , τ T .
In the case of weight 4 modular forms, there is one candidate that can be consistent with the observed neutrino masses. At τ = τ L , τ R , the neutrino mass term 3 L 3 L Y (4) 3 is similar to the case of weight 2, where two neutrino masses are degenerate. In the case of weight 4, the singlet 1 also contributes to the neutrino mass matrix through the coupling 1 . However, this additional term cannot resolve the degeneracy. It is easily noticed that two neutrino masses are degenerated also at τ = τ T since Y At τ = τ C , the triplet modular form, as seen in Table 1, is Y (4) 3 ∼ (1, 1, 1), which allows to get large mixing angles. Moreover, we have 1 and 1 modular functions. Therefore, we expect nearly tri-bimaximal mixing pattern of PMNS matrix with three different masive neutrinos. The LH weak-eigenstate neutrino fields couple to Y (4) 3 . This coupling leads to the following neutrino Majorana mass matrix: Moreover, the LH neutrino fields couple also to Y 1 and Y 1 , which gives the following additional contributions to the neutrino Majorana mass matrix M ν : where each of these two terms enters M ν with its own arbitrary constant.
To summarise, the charged lepton mass matrix could be consistent with observed masses at τ = τ L , τ R , τ = τ T for both cases of weight 2 and 4 modular forms. On the other hand, the neutrino Majorana mass matrix is consistent with observed masses only at τ = τ C for weight 4 modular forms. There is no common symmetry value of τ , which leads to charged lepton and neutrino masses that are consistent with the data.

Models with Residual Symmetry
As seen in the previous section, we could not find models with one modulus τ and with residual symmetry, which are phenomenologically viable. Therefore, we consider the case having two moduli in the theory: one τ , responsible via its VEV for the breaking of the modular A 4 symmetry in the charged lepton sector, and the another one τ ν , breaking the modular symmetry in the neutrino sector.
We present next our setup. For the charged lepton mass matrix, we take weight 2 modular forms at τ = τ T (Case I) or at τ = τ L (Case II) 1 . At the same time we use weight 4 modular forms at τ ν = τ C for constructing the neutrino Majorana mass term. In order for the modular weight in the superpotential to vanish, we assign the following weights to the LH lepton and RH charged lepton fields: where the notations are self-explanatory.
Then, the charged lepton mass matrix is obtained by using as input the expressions for the weight 2 modular forms given in Table 1. At τ T , it is a diagonal matrix: At τ = τ L , the charged lepton mass matrix has the form: The matrix M † E M E , which is relevant for the calculation of the left-handed mixing, is given in Eq. (4.4).
The neutrino mass matrix represents a sum of the contributions of modular forms of 3, 1 and 1 , with the terms involving the two singlet modular forms entering the sum with arbitrary complex coefficients A and B: where the constants of the two terms in Eq. (4.13) are absorbed in the parameters A and B.

The Neutrino Mixing
In case I, only the neutrino mass matrix contributes to the PMNS matrix since the charged lepton mass matrix is diagonal. The neutrino mass matrix in this case leads to the so called TM 2 mixing form of PMNS matrix U PMNS [35,36] where the second column of U PMNS is trimaximal: Here θ and φ are arbitrary mixing angle and phase, respectively, and P is a diagonal phase matrix containing contributions to the Majorana phases of U PMNS . Employing the standard parametrisation of U PMNS (see, e.g., [1]), it is possible to show that the trimaximal mixing pattern leads to the following relation between the reactor angle θ 13 and θ, between the atmospheric neutrino mixing angle θ 23 and θ 13 and θ, and sum rules for the solar neutrino mixing angle θ 12 and for the Dirac phase δ [35,36] (see also [9,29]): To leading order in s 13 we obtain from Eq. (5.8): where the numerical values correspond to the maximal allowed value of sin 2 θ 13 at 3σ C.L. for NO (IO) neutrino mass spectrum [2]. The interval of possible values of sin 2 θ 23 in eq. (5.11) is somewhat wider that the 3σ ranges of experimentally allowed values of sin 2 θ 23 for NO and IO spectra given in [2]. Using the 3σ allowed ranges of sin 2 θ 23 and sin 2 θ 13 for NO (IO) spectra from [2] and Eq. (5.8) we also get: where the phase ξ 31 will be specified later, We also have [9]: For further discussion of phenomenology of the neutrino trimaximal mixing (5.5), see, e.g., [9,14,30,37].
In case II, the contribution of the rotation of the charged lepton sector is added to the trimaximal mixing, which is derived from the neutrino mass matrix in Eq. (5.4). The mixing matrix in the charged lepton sector is the matrix U E in Eq. (4.6). The PMNS matrix is given by: It is straightforward to check that after a substitution θ → θ − π/2, φ → −φ, the PMNS matrix (5.19) can be rewritten as The leftmost phase matrix does not contribute to the mixing, since its effect can be absorbed into the charged lepton field phases. The rightmost phase matrix contributes only to the Majorana phases, therefore the numerical predictions in this case are the same as in Case I, apart possibly from the corresponding shift of the Majorana phases. However, as can be shown analytically, and we have confirmed numerically, also the predictions for the Majorana phases in Case II coincide with the predictions in case I.

The Neutrino Masses and Majorana Phases
It follows from (5.4) that the neutrino mass matrix M ν is a linear combination of three basis matrices: To diagonalize M ν , it is convenient to rewrite it in a different basis: where a and b are arbitrary complex coefficients and c is the overall scale factor which can be rendered real positive. M ν is diagonalized by a unitary matrix U • ν of the following form: 3 , where m i e −i2φ i are complex eigenvalues and m i ≥ 0 are the neutrino masses. 2 Extracting the phases φ i from M diag ν , we find: where the phases (φ 2 − φ 1 ) and (φ 3 − φ 1 ) contribute to the Majorana phases α 21 /2 and α 31 /2 of the standard parametrisation of the PMNS matrix [1]. Thus, the PMNS matrix has the form: where the common phase factor e −i2φ 1 is unphysical. The phase ξ 31 /2 in Eq. (5.14) can be identified now with (φ 3 − φ 1 ): ξ 31 /2 = φ 3 − φ 1 . Thus, the Majorana phases α 21 /2 and α 31 /2 are given by: The complex rotation parameters θ and φ are fixed by a choice of a and b, which we will now show explicitly. We find by direct calculation that Thus, the neutrino mass matrix M ν is diagonalized when the corresponding linear combination of the off-diagonal entries vanishes, which leads to cos 2θ + ae iφ sin 2θ + be −iφ sin 2θ = 0 ⇔ ae iφ + be −iφ = − cot 2θ. (5.28) The above condition is equivalent to: It proves convenient to introduce the complex parameter Using (θ, φ, z) is a reparametrisation of (a, b) determined by (5.29) and (5.30). The inverse parameter transformation is given by The neutrino mass matrix eigenvalues are the corresponding linear combinations of the diagonal entries in (5.27): m 3 e −i(2φ 3 +φ) = c z + 1 sin 2θ .

(5.32)
Fitting the mass-squared differences to experimentally observed values, we find the following constraint on z in terms of θ, φ and r ≡ ∆m 2 21 /∆m 2 31 : |z − z 0 | 2 = R 2 , sign (Re z) = ± sign (sin 2θ) , (5.33) where the plus (minus) sign corresponds to NO (IO) spectrum of neutrino masses, and (5.34) Since θ and r are tightly constrained by the experimental data, the set of phenomenologically viable models is effectively described by two angles φ and ψ, with the latter being the angle parameter on the circle (5.33), i.e. z = z 0 + R e iψ . Scanning through φ and ψ numerically, we find that to each set of the experimentally allowed values of the mixing angles and the mass-squared differences corresponds a range of models (parameterised by ψ) with different values of the sum of neutrino masses and the Majorana phases. We report the numerical results in the case of NO spectrum in Fig. 1. The allowed range of the sum of neutrino masses depends on the value of sin 2 θ 23 . The lower bound slightly decreases from 0.097 eV to 0.074 eV as sin 2 θ 23 runs through its 3σ confidence interval of [0.46, 0.58]. 3 On the other hand, the upper bound is highly dependent on the value of sin 2 θ 23 , and tends to infinity as sin 2 θ 23 approaches 0.5, which corresponds to δ = φ = 3π/2. This means that at this point the sum of neutrino masses is allowed to take any value greater than its lower bound of 0.093 eV. There is also a strong correlation between the Majorana phases. The set of best-fit models corresponds to φ = 1.664π and leads to the following values of observables: Similar analysis can be performed in the case of IO neutrino mass spectrum. However, in that case the minimal value of the sum of the three neutrino masses is 0.63 eV, and we do not analyse this case further.

Summary
We have investigated models of lepton masses and mixing based on modular A 4 flavour symmetry broken to residual symmetries in the charged lepton and neutrino sectors. The standard case of three lepton families was considered. In a theory based on finite modular flavour symmetry not only the matter fields, but also the constants such as the Yukawa couplings transform non-trivially under the modular symmetry. These constants are written in terms of modular forms which are holomorphic functions of a complex scalar field -the modulus τ . The modular forms have specific transformation properties under the modular symmetry transformations, which are characterised by a positive even number k called "weight", and depend on the order of the finite modular group via their "level" N . In the case of modular A 4 symmetry we have N = 3 and for the lowest weight modular forms k = 2. The modular forms transform under the usual non-Abelian discrete flavor symmetries as well. Modular forms of weight k and level N span a linear space of finite dimension. There exists a basis in this space such that the modular forms form multiplets transforming according to unitary irreducible representations of the finite modular group. In the case of modular A 4 symmetry, the dimension of the the linear space of modular forms of weight k = 2 is 3, and one can employ modular forms transforming as the triplet irreducible representation of A 4 . Modular forms of higher weights can be obtained as direct products of the modular forms of weight k = 2.
In lepton flavour models with finite modular symmetry, the modular symmetry must be broken in order to distinguish between the electron, muon and tauon, generate three different neutrino masses and reproduce the measured values of the three neutrino mixing angles. In the most economical versions of the flavour models the only source of breaking of the modular symmetry is the VEV of the modulus τ , τ = 0, and there is no need to introduce flavon fields. In the present article we consider models with modular A 4 symmetry without flavons, in which the A 4 symmetry is broken only by τ .
The modular group A 4 has two generators S and T satisfying the presentaion rules: where E is the unit operator. Residual symmetries arise whenever the VEV of the modulus τ breaks the considered finite modular group Γ N , Γ 3 A 4 , only partially, i.e., the little group (stabiliser) of τ is non-trivial. There are only 2 inequivalent finite points with non-trivial little groups, namely, τ = −1/2 + i √ 3/2 ≡ τ L and τ = i ≡ τ C [22]. The first one is the left cusp in the fundamental domain of the modular group, and corresponds to a residual symmetry associated with the subgroup Z ST There is also infinite point τ = i∞ ≡ τ T , in which the subgroup Z T 3 = {I, T, T 2 } of A 4 is preserved. We have constructed phenomenologically viable models of lepton masses and mixing based on modular A 4 invariance broken to residual symmetries Z T 3 or Z ST 3 and Z S 2 respectively in the charged lepton and neutrino sectors. The neutrino Majorana mass term is assumed to be generated by the dimension 5 Weinberg operator. We found that there is no common symmetry value of τ , which leads to charged lepton and neutrino masses that are consistent with the data. For the construction of the charged lepton mass matrix, we used weight 2 modular forms at τ = τ T (Case I) or at τ = τ L (Case II). At the same time we used weight 4 modular forms at τ ν = τ C for constructing the neutrino Majorana mass term. The so constructed two models involve three real parameters fixed by the values of the three charged lepton masses. The three neutrino masses, three neutrino mixing angles and three CPV phases are functions of altogether 2 real constants and two phases. In these models the neutrino mixing matrix is of trimaximal mixing form. In Case I it is given by the tri-bimaximal mixing matrix multiplied on the right by a unitary rotation in the 1-3 plane, which depends on one angle and one phase. Both models lead to the same phenomenology. In addition to successfully describing the charged lepton masses, neutrino mass-squared differences and the atmospheric and reactor neutrino mixing angles θ 23 and θ 13 , these models predict the values of the lightest neutrino mass (i.e., the absolute neutrino mass scale), of the Dirac and Majorana CP violation (CPV) phases and correspondingly of the effective neutrinoless double beta decay Majorana mass, as well as the existence of specific correlations between i) the values of the solar neutrino mixing angle θ 12 and the angle θ 13 , ii) the values of the Dirac CPV phase δ and of the angle θ 23 , iii) the sum of the neutrino masses and θ 23 , and iv) between the two Majorana phases (Fig. 1). These predictions will be tested with future more precise neutrino oscillation data, with results from direct neutrino mass and neutrinoless double beta decay experiments, as well as with improved cosmological measurements.
where ω = e i 2 3 π for a triplet. In this base, the multiplication rule of the A 4 triplet is More details are shown in the review [4,5].