Neutrino Mass and Mixing with $A_5$ Modular Symmetry

We present a comprehensive analysis of neutrino mass and lepton mixing in theories with $A_5$ modular symmetry. We construct the weight 2, weight 4 and weight 6 modular forms of level 5 in terms of Dedekind eta-functions and Klein forms, and their decomposition into irreducible representation of $A_5$. We construct all the simplest models based on $A_5$ modular symmetry, including scenarios of models with and without flavons in the charged lepton sectors. For each case, the neutrino masses can be generated through either the Weinberg operator or the type I seesaw mechanism. We perform an exhaustive numerical analysis, organising our results in an extensive set of figures and tables.


Introduction
The flavour puzzle, in particular the origin of neutrino mass and lepton mixing, is a major unresolved problem of the Standard Model (SM). The large mixing angles in the lepton sector can be explained using some discrete non-Abelian family symmetry [1,2]. Such a symmetry could either originate from a continuous non-Abelian gauge symmetry [3][4][5][6][7][8][9], or from extra dimensions [10][11][12][13][14][15][16][17][18][19][20][21]. In the case of extra dimensions, the discrete non-Abelian family symmetry could either arise as an accidental symmetry of the orbifold fixed points, or it could appear as a subgroup of the symmetry of the extra dimensional lattice vectors, commonly referred to as modular symmetry [22]. Some time ago it was suggested that a finite subgroup of the modular symmetry group, when interpreted as a family symmetry, might help to provide a possible explanation for the neutrino mass matrices [23,24]. Recently it has been suggested that neutrino masses might be modular forms [25], with constraints on the Yukawa couplings. This has led to a revival of the idea that modular symmetries are symmetries of the extra dimensional spacetime with Yukawa couplings determined by their modular weights [26]. The modular groups Γ(2) [27,28], Γ(3) [25,26,[28][29][30][31], Γ(4) [32,33] and Γ(5) [34] and the phenomenological predictions for lepton mixing have all been discussed in the literature.
In the present work, we focus on the finite modular group Γ(5) ∼ = A 5 . In particular, we construct the weight 2, weight 4 and weight 6 modular forms of level 5 in terms of the Dedekind eta-function and Klein forms, and their decomposition into irreducible representation of A 5 . This is complementary to [34] where only weight 2 and weight 4 modular forms are constructed using the Jacobi theta function. Since the results overlap, we have checked that the same q−expansions of the modular forms are obtained up to an overall irrelevant constant. We construct all the simplest models based on A 5 modular symmetry, both with and without flavons in the charged lepton sectors. For each case, the neutrino mass is generated through both the Weinberg operator or the type I seesaw mechanism. The model building discussed here goes beyond that in [34] where only two models were constructed with diagonal charged lepton mass matrix, one of which is not viable. For each model, we give an extensive numerical analysis, organising our results in a comprehensive set of figures and tables.
The layout of the remainder of the paper is as follows. In section 2 we review the basics of modular symmetry and modular forms of level N = 5 corresponding to A 5 . We also construct the modular space through an infinite product expansion, construct an eleven dimensional basis for the lowest weight modular forms, and decompose this space into the two irreducible triplet and the quintuplet representations of A 5 , from which the higher weight representations may be obtained. In section 3 we systematically construct all the simplest models based on Γ(5) ∼ = A 5 . The charged lepton sector is analysed both with and without flavons. The neutrino sector is analysed without flavons using both the Weinberg operator and the type I seesaw mechanism. In section 4 we perform a comprehensive numerical analysis for each of the models discussed in the previous section. We give the best fit values of the parameters of each model and the corresponding predictions. We also give the possible regions of the parameters and ranges of the predictions in graphical form. In section 5 we summarize the main results of our paper. In Appendix A we present the group theory of A 5 . In Appendix B we show that the charged lepton sector with flavons leads to diagonal charged lepton Yukawa matrices when the flavons are aligned appropriately. 2 2 Modular symmetry and modular forms of level N = 5 The full modular group SL(2, Z) is the group of 2-by-2 matrices with integer entries and determinant 1 [35,36], The modular group can be regarded as the linear fraction transformations of the upper half complex plane H = {τ ∈ C | Im τ > 0}, and it has the following form It's easy to see that γ and −γ lead to the same linear fractional transformation. Therefore the modular group of transformations is is isomorphic to the projective special linear group P SL(2, Z) = SL(2, Z)/{1, −1}, which is the quotient of SL(2, Z) over its center {1, −1}.
The modular group can be generated by two elements S and T [35] S : τ → − 1 τ , S = 0 1 −1 0 , We see that the generators S and T obey the relations: There are some important subgroups of the SL 2 (Z), i.e.
where N is a positive integer, and Γ(N ) is usually called principal congruence subgroup of level N.  [24]. For example, Γ 2 ∼ = S 3 , Γ 3 ∼ = A 4 , Γ 4 ∼ = S 4 and Γ 5 ∼ = A 5 . The multiplication rules of the group Γ N can be obtain by extending Eq. (4) with the condition T N = 1, such that the generators S and T of Γ N satisfy The modular forms f (τ ) of weight 2k and level N are holomorphic functions of the complex variable τ with well-defined transformation properties under the group Γ(N ): where k ≥ 0 is an integer. We shall only consider modular forms of even weight in this work. The function f (τ ) is required to be holomorphic in H and at all the cusps. In particular, the element T N is in Γ(N ) and we have In the present work, we shall study the Γ(5) modular group and construct some models to explain neutrino masses and mixing. The quotient space H/Γ(5) can be described by a fundamental domain F 5 for Γ (5), that is a connected region of H such that every z ∈ H can be mapped into F 5 by a Γ(5) transformation, but no two distinct points in the interior of F 5 are related under Γ(5). The space H/Γ(5) is simply F 5 with certain boundary points identified. A fundamental domain for Γ(5) is shown in figure 1. It can be constructed by the fundamental domain F = {τ ∈ H| − 1/2 ≤ Re τ ≤ 1/2 and |τ | ≥ 1} of SL(2, Z) [37]. The cusps are i∞, −2, −1, − 3 2 , − 1 2 , 0, 2 5 , 1 2 , 1, 3 2 , 2 and 5 2 . H/Γ(5) can be made compact by adding these cusp points. The compactified space H/Γ(5) has genus zero and can be thought of as an icosahedral whose vertices are the cusps. Indeed the cusps are related by transformations of Γ 5 = Γ(1)/Γ(5), which is icosahedral group, given the isomorphism between Γ 5 and A 5 .
The modular forms of weight 2k and level N = 5 form a linear space M 2k (Γ(5)), and its dimension turns out to be 10k + 1 [25,38]. For the lowest nontrivial weight 2k = 2, the dimension is equal to 11. The modular space M 2k (Γ(5)) can be constructed from the Dedekind eta-function and the Klein form. In the following, we present the definition and basic properties of the eta-function and the Klein form which are the building blocks to construct modular forms of level 5.
For any complex number τ with Im τ > 0, the Dedekind eta-function η(τ ) is defined to be the infinite product [35,36,39], This function plays an important role in constructing various modular forms of integral or 4 half-integral weight. It is well-known that η(τ ) satisfies the following identities [35,36,39] η(τ + 1) = e iπ/12 η(τ ), Then we turn to the Klein form. Let Λ = Zτ + Z be a lattice in C and any pair (r 1 , r 2 ) ∈ Q 2 − Z 2 with z = r 1 τ + r 2 , the Klein form k (r 1 ,r 2 ) (τ ) is defined by [39][40][41] where Λ = Λ \ {0}, and the product over ω = mτ + n are performed over n first and then m. Notice that there is an extra factor i/2π in the original definition of [39][40][41]. From the definition of Eq. (11) we know that the Klein form k (r 1 ,r 2 ) (τ ) is a holomorphic function which has no zeros and poles on H. The Klein form has the following infinite product expansion [39][40][41][42]: where q z = e 2πiz . It is proved that the Klein form satisfies the following transformation formulas [39][40][41] with ( r, s) = (−1) s 1 s 2 +s 1 +s 2 e −πi(s 1 r 2 −s 2 r 1 ) . Using Eq. (13), we find the Klein form k (r 1 ,r 2 ) (τ ) fulfills k (r 1 ,r 2 ) (τ + 1) = k (r 1 ,r 1 +r 2 ) (τ ), The modular space M 2k (Γ(5)) has been explicitly constructed through η function and Klein form as follow [43] M 2k (Γ(5)) = a+b=10k, a,b≥0 Obviously the dimension of M 2k (Γ(5)) is 10k + 1. For the lowest weight 2 modular forms, we could choose the basis vectors to bê √ 5 We see that the basis vectorsê i are closed under S and T up to multiplicative factors, and each element is exactly mapped into itself under the action of S 2 , (ST ) 3 and T 5 . Therefore we conclude thatê i really span the whole modular space M 2 (Γ(5)). As shown in [25], the modular space of weight 2k and level N can always be decomposed into different irreducible representations of Γ N . As a consequence, there should exist modular forms f i (τ ) of weight 2 and level 5 which transform under Γ 5 ∼ = A 5 as where is a representative element of Γ N and ρ(γ) is a unitary representation matrix of Γ N . It is sufficient that Eq. (29) is satisfied for the generators S and T , i.e.
where the representation matrix ρ(S) and ρ(T ) are given in Eq. (A.3). Solving the constraints in Eq. (31), we can construct modular form Y 3 , Y 3 and Y 5 which transform as 3, 3 and 5 of A 5 respectively, with e 1 (τ ) =ê 1 (τ ) − 36ê 6 (τ ) −ê 11 (τ ), The q-expansion of the above linearly independent modular forms e i , e i and e i are given by

5,V I = 2Y
(4) Therefore the singlet modular forms Y 1,I , Y 1,II and Y 1,III are exactly identical, the triplet modular functions Y 5,II . Note the linear space of modular forms of weight 2k and level 5 is of dimension 10k + 1 which is equal to 21 for k = 2. Thus we have constructed an explicit basis for the weight 4 modular forms at level 5.
As regards the weight 6 modular forms for k = 3, the dimension of the modular space is equal to 10 × 3 + 1 = 31. That is to say, there should be 31 independent combination Y i Y j Y k where i, j, k = 3, 3 , 5. In the same fashion as the k = 2 case, all the basis vectors and the corresponding constraints can be straightforwardly found although the involved algebra is lengthy and tedious. Firstly we consider the weight 6 modular forms which are A 5 singlet. We can construct the following four singlet modular forms of weight 6 through the tensor products of weight 2 and weight 4 modular forms listed above, which fulfill the relations Hence there is only one singlet modular form of weight 6 and it can be taken as Analogously we can construct the following twelve weight 6 modular forms transforming as 3 under A 5 , Furthermore, we find only two of them are linearly independent and they could be chosen to be Y 1,I = (e 2 1 + 2e 2 e 3 ) (e 1 , e 2 , e 3 ) T , Y In a similar manner, we can construct the following twelve A 5 triplet 3 modular forms of weight 6, Only two out of the above twelve modular forms are linearly independent and they can be taken as 1,I = (e 2 1 + 2e 2 e 3 ) (e 1 , e 2 , e 3 ) T , Eventually we find that there are two linearly independent four-dimensional and five-dimensional modular forms of weight 6 1,I = (e 2 1 + 2e 2 e 3 ) ( e 1 , e 2 , e 3 , e 4 , e 5 ) T , Y In short, a linearly independent basis of weight 6 modular forms for Γ(5) can be chosen to be Y 1,I , Y 3,I , Y 3,II , Y 3 ,I , Y 3 ,II , Y 4,I , Y 4,II , Y 5,I and Y 3 Models based on Γ 5 In this section,we shall construct some minimal modular-invariant supersymmetric model based on Γ 5 . Adopting the N = 1 global supersymmetry, the most general form of the action is where the Kähler potential K(Φ I ,Φ I ; τ,τ ) is a real gauge invariant function of the set of chiral superfields Φ I . The action is required to be invariant under the finite modular group Γ N . The supermultiplets Φ I are assumed to transform in a representation ρ I of the quotient group Γ N with a weight −k I , where ρ I (γ) is the representation matrix of the element γ and k I are integers. Under the modular transformation of Eq. (53), the Kähler potential should be invariant up to Kähler transformations. Thus one can determine the Kähler potential as where the constant h > 0. After the modulus τ develops a vacuum expectation value (VEV), the above Kähler potential leads to the following kinetic term for the scalar components of the supermultiplets Φ I and τ , For each given value of the VEV of τ , the kinetic term of φ I can be made into canonical form by rescaling the fields φ I . This effect amounts to a redefinition of the superpotential parameters in a concrete model. One the other hand, the superpotential W (Φ I , τ ) should be invariant under the group Γ N , and the total weight of W (Φ I , τ ) should be vanishing. In general, the modular invariant superpotential W (Φ I , τ ) can be expanded in power series of the supermultiplets Φ I , For the transformation rule of the field Φ I in Eq. (53), the functions Y I 1 ...In should be modular forms of weight k Y and transform in the presentation ρ Y of Γ N , The modular invariance requires In the following, we shall construct some simple models of lepton masses and mixing based on the Γ 5 modular symmetry. The neutrinos are assumed to be Majorana particles. We shall formulate our models in the framework of supersymmetry, i.e. the minimal supersymmetric standard model (MSSM) and its extension with right-handed neutrinos. In order to be as simple as possible, we suppose the Higgs doublets H u,d are singlets of A 5 and their modular weights are zero. There are freedoms for the assignments of irreducible representation and modular weight to the matter fields. We assume the three left-handed lepton doublets transform as an irreducible three-dimensional representation 3 or 3 of the Γ 5 flavor group. The generic assignments of representations and modular weights to the MSSM fields are listed in table 1. In the following, we proceed to discuss the possible structures of the models in the charged lepton and neutrino sectors.

Charged lepton sector
The most general superpotential for the charged lepton masses can be written as where ϕ denotes a generic flavon field which is a SM singlet, and f E (Y, ϕ) is the most general function of the modular forms Y (τ ) and the flavon field ϕ. Notice that the total weight of each term of W e should be vanishing. In the present work, for simplicity we shall consider the case that f E (Y, ϕ) only depends on either the flavon field ϕ or the modular forms Y . For the scenario that only ϕ is involved in W e , we can choose the modular weights k I to forbid a dependence of W e on the modulus τ . As shown in Appendix B, the flavon multiplets can develop the VEVs along certain directions such that the charged lepton mass matrix is diagonal. Then the lepton flavor mixing completely arises from the neutrino sector.
In the second scenario that the function f E (Y, ϕ) only depends on the modulus τ , the three generations of left-handed lepton doublets are assigned to the triplet representation 3 or 3 while all the right-handed charged leptons E c i transform as 1 under Γ 5 . In order to avoid a charged lepton mass matrix with rank less than 3, we assume that E c 1 , E c 2 and E c 3 have different modular weights such that modular forms of weights 2, 4, 6 (i.e. Y , Y (4) , Y (6) ) are involved. Thus the superpotential for the charged lepton masses is of the following form, where the condition of weight cancellation requires Notice that Y , Y (4) and Y (6) should transform in the same way as L under A 5 in order to form A 5 singlet. All the linearly independent modular forms of weights 2, 4 and 6 have been found out in section 2. The explicit form of the superpotential W e is different for ρ L ∼ 3 and ρ L ∼ 3 . To be more specific, using the Clebsch-Gordan coefficients given in Appendix A, we can expand W e in Eq. (60) as follow, After electroweak symmetry breaking, we can read out the charged lepton mass matrices which are summarized in the table 2. Table 2: The prediction for the charged lepton mass matrix in the models without flavon fields, where the charged lepton mass matrix is given in the right-left basis

Neutrino sector
We shall assume neutrinos are Majorana particles and the left-handed leptons are assigned to a triplet representation 3 or 3 of the modular group Γ 5 . We shall consider two different scenarios: the neutrinos are described by the effective Weinberg operator or are generated by the type-I seesaw mechanism. If the neutrino masses originate from the Weinberg operator, the most general form of the superpotential responsible for neutrino masses is, where f W (Y ) is a generic function of the modular form Y . For the sake of simplicity, we are concerned with the case that f W (Y ) is the modular form of lowest weight 2. Thus the neutrino superpotential reads The modular invariance requires the weight of lepton doublet should be equal to 1, i.e.
Since the linearly independent modular forms of weight 2 decompose into three we know that only the quintuplet modular form Y 5 contributes to W ν in Eq. (65). The concrete form of W ν depends on the transformation property of L under Γ 5 .
which gives rise to the following light neutrino mass matrix with v u = H u .
Then the neutrino mass matrix reads In the second scenario, the neutrino masses are generated via the type-I seesaw mechanism. Three right-handed neutrinos N c = (N c 1 , N c 2 , N c 3 ) T are introduced and they are assigned to transform as a Γ 5 triplet 3 or 3 . The superpotential of the neutrino sector can be generally written as where f N (Y ) and f M (Y ) are generic functions of the modular form Y . In order to build models with minimal number of parameters, we consider the case that both f N (Y ) and f M (Y ) are of weight 2 yet only one of them is involved. If the modulus parameter τ enters into the neutrino masses through the Yukawa coupling, the neutrino superpotential would be In this case the weights of N c and L should be If the complex modulus τ only appears in the right-handed neutrino mass term, the neutrino superpotential would be with the weights The explicit form of W ν depends on the assignments of ρ L and ρ N c . and From the multiplication rule 3⊗3 = 4⊕5 which doesn't contain singlet, we know that the Yukawa coupling N c LH u is not invariant under the Γ 5 modular symmetry. Consequently the superpotential W II ν is absent in this case, and we have Similar to previous case, W II ν is not allowed by the Γ 5 ∼ = A 5 flavor symmetry, and the neutrino superpotential is Analogous to the case of (ρ L , ρ N c ) = (3 , 3), the neutrino superpotential reads and Subsequently we can straightforwardly read out the modular invariant Dirac neutrino mass matrix M D and the right-handed Majorana neutrino mass matrix M N for each case above, and the results are summarized in table 3. After integrating out the heavy neutrinos N c , we can obtain the effective neutrino mass matrix given by the well-known seesaw formula, Since there are a few constructions in the charged lepton and neutrino sectors listed in above, we can obtain plenty of possible models with Γ 5 modular symmetry, as summarized in table 4. We see sixteen models are constructed, eight of them named as Ai (i = 1, . . . , 8) involve flavon fields in the charged lepton mass terms while the remaining eight models Bi (i = 1, . . . , 8) only depend on the modulus τ .

Numerical Analysis
In this section we shall perform a comprehensive numerical analysis for each possible model listed in table 4. Since some phases can be absorbed through field redefinition, some coupling constants of the models can be taken to be real. We first count the number of independent real free parameters of each model. If flavon fields are used, the resulting charged lepton mass matrix would be diagonal in our models and the observed charged lepton masses can be reproduced for certain values of the coupling constants and the flavon VEVs. If the models make no use of any flavon field other than the complex modulus τ , the charged lepton mass matrix is summarized in table 2. We can rephase the charged lepton superfields E c 1 , E c 2 , E c 3 to make the parameters α, β, γ 1 real while the phase of γ 2 can not be removed. Thus the charged lepton mass matrix depends on four independent real parameters β/α, γ 1 /α, |γ 2 /α|, Arg(γ 2 /α) except the overall scale factor αv d . If the neutrino masses originate from the Weinberg operator, the effective neutrino mass matrix would be expressed in terms of modular forms as functions of the modulus τ besides the overall factor v 2 u /Λ. If the neutrino masses are generated through the seesaw mechanism, the light neutrino mass matrix has two independent real parameters |g 1 /g 2 |, Arg(g 1 /g 2 ) and the overall scale factor is g 2 2 v 2 u /Λ for the models S1 and S5. For the cases of S2, S3, S4 and S6, the neutrino

Models
Neutrino mass matrices cases mass matrix uniquely depends on the modulus parameter τ and g 2 v 2 u /Λ controls the absolute scale of neutrino masses, as can seen from table 3. In this way, we can easily read out the independent real input parameters of our models and the results are collected in table 5.
If the modulus τ is changed to −τ , from the q−expansion of the modular forms we see that each modular form would become its complex conjugate. If we also set all the coupling constants to be their complex conjugate, i.e.
the charged lepton and neutrino mass matrices as well as the lepton mixing matrix will become complex conjugate. As a consequence, this transformation leaves lepton masses and mixing angles unchanged while the signs of both Dirac and Majorana CP phases would be flipped. Hence it is sufficient to limit in the range Reτ > 0 during the numerical analysis. We call two sets of input parameters are conjugate if they are related through the transformation in Eq. (83). Moreover we notice that the number of free parameters of each model is less than the number of low energy observables including three charged lepton masses m e,µ,τ , three light neutrino masses m 1,2,3 and three lepton mixing angles θ 12 , θ 13 , θ 23 , one Dirac CP phase δ CP and two Majorana phases α 21 and α 31 . In order to quantitatively measure how well the models can describe the experimental data, we define the χ 2 function to serve as a Models mass matrices assignment weight Table 4: The summary of models and the corresponding predictions for neutrino and charged lepton mass matrices. For the models with flavons in the charged lepton sector, the weights of the right-handed charged leptons should satisfy the constraint in Eq. (B.4a), i.e. k E1 = 5k E3 + 4k L and k E2 = 4k E3 + 3k L .

Models free input parameters
{Re τ, Im τ, β/α, γ 1 /α, |γ 2 /α|, Arg(γ 2 /α), |g 1 /g 2 |, Arg(g 1 /g 2 )} αv d , g 2 2 v 2 u /Λ test-statistic for the goodness-of-fit, where p i denote the input free parameters of a model and they are listed in table 5, and Q i ∈ {sin 2 θ 12 , sin 2 θ 13 , sin 2 θ 23 , r ≡ ∆m 2 21 /|∆m 2 3 |, m e /m µ , m µ /m τ } are observable quantities derived from the neutrino and charged lepton mass matrices as complex nonlinear functions of the free parameters of the models 1 . The parameters Q j,best-fit and σ j refer to the current central values and 1σ deviations respectively of the corresponding observable quantities listed in the table 6. Since the indication of a preferred value of the Dirac CP violation phase δ CP coming from global data analyses is rather weak [44], we do not include any information on δ CP in the χ 2 function. We consider both the normal ordering (NO) neutrino mass spectrum m 1 < m 2 < m 3 and the inverted ordering (IO) neutrino masses m 3 < m 1 < m 2 , where m 1,2,3 denote three light neutrino masses. We define a quantity N σ ≡ χ 2 min where χ 2 min is the global minimum of function χ 2 . The free parameters of the models are scanned over the following ranges, β/α, γ 1 /α, |γ 2 /α|, The modulus τ is taken from the fundamental domain
We see that the observed quantities sin 2 θ 12 , r ≡ ∆m 2 21 /|∆m 2 3 | and δ CP can be reproduced very well, particularly the Dirac phase is close to 3π/2. However, both reactor mixing angle θ 13 and atmospheric angle θ 23 are slightly larger than the 3σ upper bounds, and their pulls dominate the minimum of χ 2 . Similar to Ref. [26], we expect that the discrepancy could be resolved if there is small deviation from the leading order charged lepton flavon alignment such that the charged lepton mass matrix is not diagonal.

(88)
We can see that both sin 2 θ 12 and r ≡ ∆m 2 21 /|∆m 2 3 | are in the 1σ region [44], while sin 2 θ 23 is a bit larger than its 3σ upper limit. Moreover, from Eq. (88) we know the sum of neutrino masses is Under the assumption of ΛCDM cosmology, the latest Planck result on the neutrino mass sum is [48].
Since the upper bound of the neutrino mass sum depends on the cosmological model and whether other experimental data such as baryon acoustic oscillation, gravitational lensing of galaxies and the high multipole TT, TE and EE polarization spectra [48] are considered, we shall not seriously take the constraint in Eq. (90). The neutrinoless double beta (0νββ) decay is the unique probe for the Majorana nature of neutrinos, and it depends on the values of the Majorana CP violation phases. The 0νββ decay experiments can provide valuable information on the neutrino mass spectrum and constrain the Majorana phases. The 0νββ decay rate is proportional to the effective Majorana mass |m ee | defined as [49], |m ee | = |m 1 cos 2 θ 12 cos 2 θ 13 + m 2 sin 2 θ 12 cos 2 θ 13 e iα 21 + m 3 sin 2 θ 13 e i(α 31 −2δ CP ) | . (91) From Eq. (88) we can extract the following predicted value for |m ee | in the A3 model with NO, which is within the reach of future 0νββ experiments. For the model A3 with IO and model A7 with either NO or IO, there exist parameter spaces to be completely consistent with the experimental data on neutrino masses and mixing angles. We find that the allowed ranges of the coupling constants and theoretical predictions crucially depend on the values of the complex modulus τ . In each case, we find several distinct local minima of χ 2 in the fundamental domain of τ . We denote the phenomenologically viable regions of τ as region1, region2, region3 and so on, as shown in figure 2, where we have focused on Reτ > 0 and the conjugate regions are not shown for simplicity. We perform an extensive numerical scan over the parameter space, and the results of the numerical analysis are summarized in tables 7-12. The allowed regions of the input parameters and observables are determined by requiring all the lepton mixing angles and the squared mass splittings ∆m 2 21 and ∆m 2 3 within the 3σ intervals [44]. Since the models depend on few free parameters, we find that the lepton mixing angles, CP violation phases, the lightest neutrino mass m min and the effective Majorana mass |m ee | are generally correlated with each other. We display the different correlations in figures 3-12. We see that the Dirac CP violation phase δ CP is predicted to lie in narrow regions in most cases. The forthcoming long-baseline neutrino oscillation experiments will considerably improve the sensitivity to δ CP if running in both the neutrino and the anti-neutrino modes, we expect the predictions for δ CP in this paper could be tested by future neutrino facilities. Furthermore, we notice that both m min and |m ee | can only obtain values in very limited ranges. The next generation experiments searching for 0νββ decay will be able to probe almost all the IO region, up to m ee 0.02 eV, thus allowing to testing these modular symmetry models as well.

Numerical results of the models without flavons
Taking into account the freedom of field redefinition, the parameters α, β and γ 1 in the charged lepton mass matrix can be set to be real. However, the coupling constant γ 2 is generally complex. The neutrino mass matrix m ν only depends on the modular parameter τ up to an overall factor in the models B1, B2, B4, B5, B6 and B8, and an additional complex parameter g 1 /g 2 is involved in the models B3 and B7. The vacuum expectation value of the complex modulus τ is the unique source of flavor symmetry breaking in these

Region1
Region2 A3 Best         [44]. The present most stringent upper limits |m ee | < 0.120 eV from EXO-200 [50,51] and KamLAND-ZEN [52] is represented by horizontal grey band. The vertical grey exclusion band denotes the current bound from the cosmological data of i m i < 0.130 eV at 95% confidence level obtained by the Planck collaboration [48]. Notice that the bound on i m i sensitively depends on the cosmological model and whether other experimental data such as BAO and gravitational lensing etc are included [48].         models. In the same fashion as section 4.1, we use the package MINUIT to search for the minimum of the χ 2 function. We find that only the models B3 and B7 can accommodate the experimental results, and the neutrino mass spectrum can be either NO or IO. Since the number of the input parameters is larger than that of the models with flavons, it takes a lot of time to comprehensively scan the whole parameter space of the models B3 and B7. We report the numerical results of the χ 2 analysis in table 13, and the values of the lepton mixing parameters and neutrino masses and charged lepton masses at the best fitting points are displayed. It is remarkable that all the observables are quite close to their central values, and the Dirac CP phase is approximately maximal δ CP 1.599π in the B3 model for NO neutrino masses.

Conclusion
In this paper we have provided a comprehensive analysis of neutrino mass and lepton mixing in theories with Γ 5 ∼ = A 5 modular symmetry. We have constructed the weight 2, weight 4 and weight 6 modular forms of level 5 in terms of Dedekind eta-function and Klein forms, and their decomposition into irreducible representation of A 5 . We have provided the necessary mathematical tools which may be useful for future studies of model building based on A 5 modular symmetry.
We have constructed all the simplest models based on A 5 modular symmetry, including models with and without flavons in the charged lepton sectors. For each case, the neutrino masses are considered both using the Weinberg operator and the type I seesaw mechanism. This extends the scope of the previous analyses for A 5 modular symmetry in the literature. The hope is that the compendium of models presented here may be a useful guide to future model building directions.
We have performed an exhaustive numerical analysis for each model, presenting our results in the form of extensive sets of figures and tables. Surveying the comprehensive results in this paper, we observe that the predictions for lepton mixing strongly depend on the modular parameter τ , in other words the assumed vacuum expectation value of the modulus field. The modular symmetry models are therefore not as predictive as might have been expected. If the value of τ is not fixed by any mechanism, the lepton mixing parameters and neutrino masses can vary widely in different regions of τ parameter space. This observation is even true for benchmark modular symmetry models which only depend on τ and overall couplings. For this reason we divided the modular parameter τ into several regions when performing the numerical analysis. This motivates models in which the modulus τ is fixed, as recently discussed in [53] or in the context of string compactifications and supergravity [54][55][56].

B Charged lepton sector with flavons
In this appendix, we shall show that a diagonal charged lepton mass matrix can be achieved if the VEVs of the flavon fields are aligned in appropriate directions. The lefthanded lepton doublets L are assigned to a triplet 3 of A 5 , and the right-handed charged leptons E c i transform as 1. In order to achieve the desired vacuum alignment and to reproduce the observed charged lepton mass hierarchies, the auxiliary symmetry Z 5 is invoked. We use the supersymmetric driving field mechanism to determine the vacuum alignment of the flavons [23], and a continuous U (1) R symmetry related to the usual R−parity is assumed. The R−charge of the matter superfields is +1, the Higgs and flavon fields are uncharged and the so-called driving fields, indicated with the superscript "0" carry two units of R−charge. The field content and the symmetry assignments are shown in table 14. The superpotential of the charged lepton sector can be written as   with 3) The modular weight of each term should be zero such that the following constraints should be satisfied,