arxiv: v1 [hep-ph] 27 Apr 2018

Transcription

1 HUPD181 Revisiting A 4 model for leptons in light of NuFIT 3. arxiv: v1 [hep-ph] 7 Apr 18 Sin Kyu Kang 1,, Yusuke Shimizu,, Kenta Takagi,, Shunya Takahashi,, Morimitsu Tanimoto 3, 1 School of Liberal Arts, Seoul-Tech, Seoul , Korea Graduate School of Science, Hiroshima University, Higashi-Hiroshima , Japan 3 Department of Physics, Niigata University, Niigata , Japan Abstract We revisit the A 4 model for leptons in light of new result of NuFIT 3.. We introduce a new flavon η transforming as A 4 singlet 1 or 1 which couples to both charged leptons and neutrinos in next-leading order operators. The model consists of the five parameters: the lightest neutrino mass m 1, the vacuum expectation value of η and three CP violating phases after inputting the experimental values of m atm and m sol. The model with the 1 singlet flavon gives the prediction of sin θ 1 around the best fit of NuFIT 3. while keeping near the maximal mixing of θ 3. Inputting the experimental mixing angles with the 1 σ error-bar, the Dirac CP violating phase is clearly predicted to be δ CP = 5 1, which will be tested by the precise observed value in the future. In order to get the best fit value sin θ 3 =.538, the sum of three neutrino masses is predicted to be larger than 9 mev. The cosmological observation for the sum of neutrino masses will also provide a crucial test of our predictions. It is remarked that the model is consistent with the experimental data only for the normal hierarchy of neutrino masses. address: skkang@seoultech.ac.kr address: yu-shimizu@hiroshima-u.ac.jp address: takagi-kenta@hiroshima-u.ac.jp address: s-takahashi@hiroshima-u.ac.jp address: tanimoto@muse.sc.niigata-u.ac.jp

2 1 Introduction The origin of the quark/lepton flavor is still unknown in spite of the remarkable success of the standard model (SM). To reveal the underlying physics of flavors is the challenging work. The recent development of the neutrino oscillation experiments provides us with important clues to investigate the flavor physics. Indeed, the neutrino oscillation experiments have determined precisely two neutrino mass squared differences and three neutrino mixing angles. Especially, the recent data of both TK [1, ] and NOνA [3, 4] give us that the atmospheric neutrino mixing angle θ 3 favors near the maximal angle 45. The global analysis by NuFIT 3. presents the best fit θ 3 = 47. for the normal hierarchy (NH) of neutrino masses [5]. The closer the observed θ 3 is to the maximal mixing, the more likely we are to believe in some flavor symmetry behind it. In addition to the precise measurements of the mixing angles, the TK and NOνA strongly indicate the CP violation in the neutrino oscillation [, 4]. Thus, we are in the era to develop the flavor structure of the lepton mass matrices with focus on the leptonic flavor mixing angles and CP violating phase. Before the reactor experiments measured non-zero value of θ 13 in 1 [6, 7], the paradigm of the tri-bimaximal (TBM) mixing [8, 9], a highly symmetric mixing pattern for leptons, has attracted much attention. It is well known that this mixing pattern is derived in the framework of the A 4 flavor symmetry [1]-[13]. Therefore, non-abelian discrete groups have become the center of attention in the flavor symmetry [14]-[17]. In order to obtain non-vanishing θ 13, two of the authors improved the A 4 model by the minimal modification through introducing another flavon which transforms as 1 ( ) of A 4 and couples only to the neutrino sector [18]. Then, the predicted values of θ 13 are consistent with the experimental data. This pattern is essentially the trimaximal mixing TM [19,, 1] which leads to sin θ 1 1/3. However, the predicted sin θ 1 is outside of σ interval of the experimental data in the NuFIT 3. result [5]. Therefore, the A 4 model should be reconsidered in light of the new data of TK and NOνA as the implication of the new data has been changed. In this work, we introduce a new flavon transforming as A 4 singlet, η (1 or 1 ) which couples to both charged leptons and neutrinos in next-leading order operators. The model consists of the five parameters: the lightest neutrino mass m 1, the vacuum expectation value (VEV) of η and three CP violating phases after inputting the observed values of m atm and m sol. The model with a 1 singlet flavon gives the prediction of sin θ 1 around the best fit of NuFIT 3. with keeping near the maximal mixing of θ 3. The non-vanishing θ 13 is derived from both charged lepton and neutrino sectors. Inputting the observed mixing angles with the 1 σ error-bar, the CP violating Dirac phase is clearly predicted to be δ CP = 5 1. Therefore, the observation of the CP violating phase is essential to test the model in the future. It is remarked that the model is consistent with the experimental data only for NH of neutrino masses. The inverted hierarchy (IH) of neutrino masses is not allowed in the recent experimental data. This situation comes from that the singlet 1 or 1 flavon couples to leptons in the nextleading order. It is contrast with the model in Ref.[18] where both NH and IH are allowed. We present our framework of the A 4 model in section where lepton mass matrices and VEVs of scalars are discussed. The numerical results are shown in section 3. The section 4 is devoted to the summary. Appendix A shows the lepton mixing matrix and CP violating measures which are used in this work. The relevant multiplication rules of A 4 are represented in Appendix B. The derivation of the lepton mixing matrix is given in Appendix C. Appendix D presents the 1

3 distributions of our parameter which are used in our numerical calculations. Our framework of A 4 model We discuss our A 4 model in the framework of the supersymmetry (SUSY). In the non-abelian finite group A 4, there are four irreducible representations: 1, 1, 1 and 3. The left-handed leptons l and right-handed charged leptons e c, µ c, τ c are assigned to the triplet and singlets, respectively, as seen in Table 1. The two Higgs doublets (h u, h d ) are assigned to the A 4 singlets, and their VEVs are denoted as (v u, v d ) as usual. We introduce several flavons as listed in Table 1. The flavons φ T and φ S are A 4 triplets while ξ and ξ are the same A 4 singlet 1. In addition, η and η are the same non-trivial singlet 1 or 1. The A 4 flavor symmetry is spontaneously broken by VEVs of gauge singlet flavons, φ T, φ S, ξ and η, whereas ξ (1) and η (1, 1 ) are defined to have vanishing VEVs through the linear combinations of ξ and ξ and η and η, respectively, as discussed in Ref. [13]. In the original model proposed by Altarelli and Feruglio [1, 13], φ T, φ S and ξ were introduced, and then the specific vacuum alignments of the triplet flavons lead to the tri-bimaximal mixing where the lepton mixing angle θ 13 vanishes. In 11, two of the authors minimally modified the model by introducing an extra flavon η (1 ) on top of those flavons to generate non-vanishing θ 13 [18]. This modification of the model leads to the trimaximal mixing of neutrino flavors, so called TM which predicts sin θ 1 1/3 [19,, 1]. Unfortunately, this prediction for θ 1 is inconsistent with the data at σ C.L. given in the NuFIT 3. result [5]. In this work, we force the flavon η (1 or 1 ) to couple to both charged lepton and neutrino sectors in next-leading operators by assigning a Z 3 charge to η appropriately. We impose the Z 3 symmetry to control Yukawa couplings in both neutrino sector and charged lepton sector. The third row of Table 1 shows how each chiral multiplet transforms under Z 3 with its charge ω = exp(πi/3). In order to obtain the natural hierarchy among lepton masses m e, m µ and m τ, we resort to the Froggatt-Nielsen mechanism [] with an additional U(1) FN symmetry under which only the right-handed lepton sector is charged. The field Θ denotes the Froggatt-Nielsen flavon in Table 1. The U(1) FN charges are taken as (4,, ) for (e c, µ c, τ c ), respectively. By assuming that Θ, carrying a negative unit charge of U(1) FN, acquires a VEV, the relevant mass ratio is reproduced through the Froggatt-Nielsen charges. We also introduce a U(1) R symmetry in Table 1 to distinguish the flavons and driving fields φ T, φ S, ξ and η, which are required to build a non-trivial scalar potential so as to realize the relevant symmetry breaking. l e c µ c τ c h u,d φ T η η φ S ξ ξ Θ φ T η φ S ξ SU() A (1 ) 1 (1 ) (1 ) 3 1 Z 3 ω ω ω ω ω ω ω 1 1 ω ω ω U(1) FN 4 1 U(1) R Table 1: Assignments of leptons, Higgs, flavons and driving fields, where ω = exp (πi/3).

4 In these setup, the superpotential for respecting A 4 Z 3 U(1) FN U(1) R symmetry is written by introducing the cutoff scale Λ as w = w Y + w d, w Y = w l + w ν, w l = y e (φ T l) 1 e c h d Θ 4 /Λ 5 + y µ (φ T l) 1 µ c h d Θ /Λ 3 + y τ (φ T l) 1 τ c h d /Λ + y e (φ T l) 1 (1 ) ec h d ηθ 4 /Λ 6 + y µ (φ T l) 1 (1) µc h d ηθ /Λ 4 + y τ (φ T l) 1(1 ) τ c h d η/λ, w ν = y S (ll) 3 h u h u φ S /Λ + y ξ (ll) 1 h u h u ξ/λ + y 1(ll) 1 h u h u (φ S φ T ) 1 /Λ 3 + y (ll) 1 h u h u (φ S φ T ) 1 /Λ 3 + y 3(ll) 1 h u h u (φ S φ T ) 1 /Λ 3 + y 4(ll) 3 h u h u (φ S φ T ) 3 /Λ 3 + y 5(ll) 3 h u h u φ S η/λ 3 + y 6(ll) 3 h u h u ξφ T /Λ 3 + y 7(ll) 1 (1 )h u h u ξη/λ 3, w d = w T d + w S d, w T d = Mφ T φ T + gφ T φ T φ T + λφ T φ T η λ 1 η φ T φ S + λ η ηξ + λ 3 η η ξ + λ 4 η ηξ + λ 5 η η ξ, w S d = g 1 φ S φ S φ S + g φ S φ S ξ g3 ξ φ S φ S + g 4 ξ ξξ + g 5 ξ ξ ξ + g 6 ξ ξ ξ, (1) where the subscripts 1 (1 ) in (φ T l) 1 (1 ), etc. correspond to the case of η for 1 (1 ). The Yukawa couplings y s and y s are complex number of order one, M is a complex mass parameter while g s and λ s are trilinear couplings which are also complex number of order one. Both leading operators and next-leading ones are included in w Y, which leads to the flavor structure of lepton mass matrices including next-leading corrections. On the other hand, w d only contains leading operators, where we can force ξ( η) to couple with φ S φ S (φ T φ T ), but not ξ(η) with it since ξ and ξ( η and η) have the same quantum numbers [13]. We can study the vacuum structure and lepton mass matrices with these superpotential..1 Vacuum alignments of flavons Let us investigate the vacuum alignments of flavons. The superpotentials wd T and ws d are written in terms of the components of triplet flavons: in Eq. (1) wd T = M ( ) ( ) φ T 1φ T 1 + φ T φ T 3 + φ T 3φ T + λ φ T 1 φ T + φ T φ T 1 + φ T 3φ T 3 η + g [ ( ) ( ) ( )] φ T 3 1 φ T 1 φ T φ T 3 + φ T φ T φ T 1 φ T 3 + φ T 3 φ T 3 φ T 1 φ T λ 1 η (φ T φ S + φ T 1 φ S3 + φ T 3 φ S1 ) + λ η ηξ + λ 3 η η ξ + λ 4 η ηξ + λ 5 η η ξ, wd S = g 1 [ ( ) ( ) ( )] φ S 3 1 φ S1 φ S φ S3 + φ S φ S φ S1 φ S3 + φ S 3 φ S3 φ S1 φ S ( ) + g φ S 1 φ S1 + φ S φ S3 + φ S 3φ S ξ g 3 ξ ( φ S1 + φ S φ S3 ) + g4 ξ ξ + g 5 ξ ξ ξ + g 6 ξ ξ, () where wd S is the same superpotential given in Ref. [13]. Note that new terms including η and η are added in wd T. 3

5 Then, the scalar potential of the F -term is given as V V T + V S, V T = w T d φ T + h.c. i i = Mφ T 1 + λφ T η + g ( ) φ 3 T 1 φ T φ T 3 + Mφ T 3 + λφ T 1 η + g ( ) φ 3 T φ T 1 φ T 3 + Mφ T + λφ T 3 η + g ( ) φ 3 T 3 φ T 1 φ T + λ 1 (φ T φ S + φ T 1 φ S3 + φ T 3 φ S1 ) + λ ηξ + λ 3 η ξ + λ 4 ηξ + λ 5 η ξ, V S = wd S X + h.c. = g 1 ( ) φ 3 S1 φ S φ S3 + g φ S1 ξ + g 1 ( ) φ 3 S φ S1 φ S3 + g φ S3 ξ + g 1 ( ) φ 3 S3 φ S1 φ S + g φ S ξ ( ) + g 3 φ S1 + φ S φ S3 + g4 ξ + g 5 ξ ξ + g 6 ξ. (3) The vacuum alignments of φ T, φ S and VEVs of η, η, ξ and ξ are derived from the condition of the potential minimum, that is V T = and V S = in Eq.(3) as φ T =v T (1,, ), φ S = v S (1, 1, 1), η = q, η =, ξ = u, ξ =, v T = 3M g, v S = g 4 3g 3 u, q = λ 1v T v S λ u = λ 1 g4 v T, (4) λ 3g 3 where the VEVs of ξ and η are taken to be zero by the linear transformation of ξ and ξ (η and η) without loss of generality. The coefficients λ i and g i are of order one since these flavons have no FN charges. Therefore, the VEVs of η and ξ are of same order as v T and v S, respectively. In our numerical analyses, q/λ is scanned around v T /Λ which is fixed by the tau-lepton mass. On the other hand, the FN flavon Θ is not contained in w d due to the U(1) FN invariance. The VEV of Θ can be derived from the scalar potential of D-term by assuming gauged U(1) FN. The Fayet-Iliopolos term leads to the non-vanishing VEV of Θ as discussed in Ref. [3]. Thus, its VEV is determined independently of v T, v S, u and q.. Lepton Mass Matrices The explicit lepton mass matrices are derived from the superpotentials w l and w ν in Eq. (1) by use of the multiplication rule of A 4 given in Appendix B. Let us begin with writing down the 4

6 charged lepton mass matrices by imposing the vacuum alignments in Eq.(4) as: y e λ 4 y τα η y e λ 4 y µα η λ M l = v d α l y eα η λ 4 y µ λ for η(1 ), v d α l y µ λ y τα η for η(1 ), y µα η λ y τ y eα η λ 4 y τ (5) where α l, α η and λ are defined in terms of the VEVs of φ T, η and Θ, respectively: α l φ T Λ = v T Λ, α η η Λ = q Λ, λ Θ Λ. (6) We note that the off-diagonal elements arise from the next-leading operators. The left-handed mixing matrix of the charged lepton is derived by diagonalizing M l M l. We obtain the mixing matrix U l approximately for the cases of η being 1 or 1 of A 4 as (more explicitly presented in Appendix C): 1 U 1 O(αη) αηe τ iϕ l O(αη) 1 + α τ 1 + α τ η O(α η λ 4 ) for η(1 ), η αηe τ iϕ O(αη) 3 1 U l 1 + α 1 1 τ η 1 + α τ η α ηe µ iϕ O(αηλ 4 ) α 1 + α η µ 1 + α τ ηe µ iϕ α η µ αηe τ iϕ for η(1 ), (7) η O(αη) αηe τ iϕ 1 + α η µ where α τ ηe iϕ y τ y τ α η, α µ ηe iϕ y µ y µ α η. (8) The mass eigenvalues m e, m µ and m τ are obtained by U l M l M l U l as shown in Appendix C. In the leading order approximation, U l depends on one real parameter αη τ and one phase ϕ for the case of η(1 ), whereas it depends on αη, τ α η, µ ϕ and ϕ for the case of η(1 ). The parameter α η is expected to be much less than 1 as discussed in the next section. As seen in Eq.(7), the off-diagonal (1,3) and (3,1) entries in U l are dominant for the case of η(1 ) while the off-diagonal (1,) and (,3) (also (.1) and (3,)) entries in U l are dominant for the case of η(1 ). Thus, it is expected that the assignments of η(1 ) and η(1 ) give rise to different predictions of the mixing and the CP violation. It is found that the effects of the next-leading terms of O(αη), O(αη) 3 and O(α η λ 4 ) in the mixing matrix U l are negligibly small by our numerical estimation. The neutrino mass matrix is derived from the superpotential w ν in Eq. (1) by imposing the vacuum alignments given in Eq.(4). The next-leading operator y 5llh u h u φ S η can be absorbed in the leading one y S llh u h u φ S due to the alignment of φ S (1, 1, 1). Although the next-leading operators llh u h u φ S φ T and llh u h u φ T ξ cannot be absorbed in the leading one, their effects are expected to be suppressed because φ T /Λ is fixed to be small. We have confirmed that the effect of those next-leading operators is negligibly small in our numerical calculations. 5

7 On the other hand, the operator y 7llh u h u ξη leads to the significant contribution to the neutrino mass matrix because η /Λ could be significantly larger than φ T /Λ as discussed in Appendix D. For η(1 ), we have M ν = a 1 + b c 1 + d 1, (9) where the coefficients a, b, c and d are given in terms of the Yukawa couplings and VEVs of flavons as follows: with a = y Sα ν Λ v u, b = y Sα ν 3Λ v u, c = y ξα ξ Λ v u, d = y 7α ξ α η Λ v u, (1) α ν φ S Λ = v S Λ, α ξ ξ Λ = u Λ. (11) Since the parameter d is induced from the next-leading operator llξηh u h u, the magnitude of d is expected to be much smaller than a, b and c. For η(1 ), we get M ν = a 1 + b c 1 + d 1, (1) where the last matrix of the right-hand side is a different one compared with the case of η(1 ). There are three complex parameters in the model since the coefficient b is given in terms of a. We take a to be real without loss of generality and reparametrize them as follows: a a, c c e iφc, d d e iφ d, (13) where a, c and d are real parameters and φ c, φ d are CP violating phases. For the lepton mixing matrix, Harrison-Perkins-Scott proposed a simple form of the mixing matrix, so-called TBM mixing [8, 9], V TBM = , (14) by which M ν is diagonalized in the case of d =. We obtain the neutrino mass matrix in the TBM basis by rotating it with V TBM as: a + ce iφc d ˆM ν = VT T eiφ d 3 deiφ d BMM ν V T BM = ce iφc + de iφ d, (15) 3 deiφ d a ce iφc + d eiφ d where upper (lower) sign in front of (1,3) and (3,1) components correspond to η transforming as 1 (1 ). The neutrino mass eigenvalues are explicitly given in Appendix C. 6

8 The mixing matrix U ν is derived from the diagonalization of ˆMν ˆM ν apart from the Majorana phases such as U ν ( ˆM m ν ˆM ν ) U ν 1 = m. m 3 (16) As shown in Appendix C, we get cos θ sin θe iσ U ν = 1, (17) sin θe iσ cos θ where θ and σ are given in terms of parameters in the neutrino mass matrix. As seen in Eq.(1), the parameter d is related with c as d c = y 7 α η αη ν, (18) y ξ where y 7 and y ξ are coefficients of order one. On the other hand, a and c are given in terms of m 1, αη ν and the experimental data m sol and m atm as shown in Appendix C. Therefore, m 1 and αη ν are free parameters as well as φ c and φ d in our model. It is remarkable that neutrino mass eigenvalues do not satisfy m sol > for the case of IH of neutrino masses as discussed in Appendix C because of the relation, a c and c d, in our model. It is understandable by considering the case of d = limit which corresponds to the exact TBM mixing. It is allowed only for NH of neutrino mass spectrum. Finally, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [4, 5] is given as U PMNS = U l V TBM U ν P, (19) where P is the diagonal matrix responsible for the Majorana phases obtained from P U ν ˆMν U T ν P = diag{m 1, m, m 3 }, () where m 1, m and m 3 are real positive neutrino masses. The effective mass for the neutrinoless double beta (νββ) decay is given as follows: m ee = m1 Ue1 + m Ue + m 3 Ue3, (1) where U ei denotes each component of the PMNS matrix U PMNS, which includes the Majorana phases. From Eq.(19), we can write down the three neutrino mixing angles of Appendix A in terms of our model parameters for the case of 1 singlet η, which shows how experimental results can be accommodated in our model: 1 1 sin θ 1 1 αηe τ iϕ, 1 + α τ 3 η 1 sin θ 13 sin θe iσ 1 α 1 + α τ η τ cos θe iϕ, 6 η sin θ 3 1 cos θ 1 sin θe iσ, () 6 7

9 where the next-leading terms are omitted. It is remarkable that sin θ 13 is composed of contributions from both the charged leptons and neutrinos. On the other hand, the deviation from the trimaximal mixing of θ 1 comes from the charged lepton sector, whereas the deviation from the maximal mixing of θ 3 comes from the neutrino sector. Since these are given in terms of four independent parameters, we cannot obtain the sum rules in the PMNS matrix elements. However, the tau-lepton mass helps us to predict the allowed region of the CP violating Dirac phase δ CP and Majorana phases α 1 and α 31 as discussed in the next section. 3 Numerical results observable best fit and 1 σ 3 σ range m atm ( ) 1 3 ev ( ) 1 3 ev m sol ( ) 1 5 ev (6.8 8.) 1 5 ev sin θ sin θ sin θ Table : The best fit, 1 σ and 3 σ ranges of neutrino oscillation parameters from NuFIT 3. for NH [5]. At first, we present the framework of our calculations to predict the CP violating Dirac phase δ CP and Majorana phases α 1 and α 31. We explain how to get our predictions in terms of three real parameters α τ η, α ν η and m 1 on top of three phases ϕ, φ c and φ d for NH of neutrino masses. We can put for simplicity α η = α τ η = α ν η, (3) that is y 7/y ξ = y τ/y τ = 1 since all Yukawa couplings are of order one. The result of NuFIT 3. [5] is used as input data to constrain the unknown parameters. By taking m 3 m 1 = m atm and m m 1 = m sol with 3 σ and 1 σ data given in Table, a, c and d are fixed in terms of m 1, α η, φ c and φ d. There is also the CP violating phase ϕ in the charged lepton mixing matrix. In our numerical analysis, we perform a parameter scan over those three phases and m 1 by generating random numbers. The scan ranges of the parameters are π (ϕ, φ c, φ d ) π and m 1 5 mev. Note that the range of m 1 is restricted by the lower bound of cosmological data for the sum of neutrino masses, 16 mev [6]. The parameter α η is constrained by the tau-lepton mass: m τ = y τ α l v d, (4) which gives α l =.316 and.1 for the minimal supersymmetric standard model (MSSM) with tan β = 3 and SM, respectively. Here we put y τ = 1. Since α η is of same order as α l as seen in Eq.(4), we vary the parameter α η around α l =.316 (.1) by using the Γ distribution (χ distribution), which is presented in Appendix D. 8

10 We calculate three neutrino mixing angles in terms of the model parameters while keeping the parameter sets leading to values allowed by the experimental data at 1 σ and 3 σ C.L. as given in Table. Then, we calculate the CP violating phases and mee with those selected parameter sets. Accumulating enough parameter sets survived the above procedure, we make various scatter plots to show how observables depend on the model parameters. In subsection 3.1, we show our numerical results for η(1 ). The numerical results for η(1 ) are briefly shown in subsection 3.. Case of a 1 singlet η D D sin Θ sin Θ1 Figure : The allowed region on sin θ1 δcp plane. The meaning of colors is the same as in Fig.1. The red curve represents the model without the rotation to the neutrino mass matrix in the TBM basis. Figure 1: The allowed region on sin θ3 δcp plane, where the blue and green dots correspond to the input of 3 σ and 1 σ data in Table, respectively. The red curve represents the prediction of TM. Let us show numerical results for the case of a 1 singlet η. We analyze only the case of NH of neutrino masses since the case of IH of neutrino masses is inconsistent with experimental data as discussed in Appendix C. At first, we show the prediction of δcp versus sin θ3 in Fig. 1 where the blue and green dots correspond to the input of 3 σ and 1 σ data in Table, respectively. This result is similar to the prediction of TM since the deviation from the maximal mixing of θ3 is due to the extra (1-3) family rotation of the neutrino mass matrix in Eq.(15). In order to compare our prediction with the TM result [7, 8], we draw its prediction by a red curve which is obtained by taking the best fit data in Table. We see that our predicted region is inside of the TM boundary. For the maximal mixing θ3 = π/4, the absolute value of δcp is expected to be 6 9. It is also predicted to be 9. δcp. 11 at the best fit of sin θ3 =.538. All values between 18 and 18 are allowed for δcp in the case of the input data at 3 σ as seen in Fig. 1. However, for the input data at 1 σ, δcp is restricted to be 5 1, which is completely consistent with the present data at 1 σ, 157. δcp. 83 apart from its sign. Thus, the precise data of θ3 and δcp would provide us with a crucial test of our prediction. Next, we show the prediction of δcp versus sin θ1 in Fig.. The deviation from the trimaximal mixing of θ1 is due to the (1-3) family rotation of the charged lepton sector as seen in Eq.(). The model without the additional rotation to the neutrino mass matrix in the TBM 9

11 basis presented a clear correlation between sin θ1 and δcp [7, 8]. We also draw its prediction by a red curve which is obtained by taking the best fit data in Table. Predicted points are scattered around the red curve. Our predicted region is broad for the 3 σ data of mixing angles. However, 1 σ data forces the predicted region to be rather narrow. Then, δcp = 6 1 is predicted at the best fit of sin θ1 =.37, where the maximal CP violation δcp = 9 is still allowed. On the other hand, we cannot find any correlation between δcp and sin θ13 since both phases σ in the neutrino mass matrix and ϕ in the charged lepton mass matrix contribute to sin θ13 as seen in Eq.(). We omit to present the result in a figure sin Θ3 sinθ ΑΗ ΑΗ ΑΗ Figure 4: The allowed region on αη sin θ3 plane. The meaning of colors is the same as in Fig. 1. D sin Θ1 Figure 3: The allowed region on αη sin θ13 plane, where the 3 σ data is taken except for sin θ13. The red lines represent the upper and lower bounds of the experimental data Figure 5: The allowed region on αη sin θ1 plane. The meaning of colors is the same as in Fig ΑΗ.3.4 Figure 6: The allowed region on αη δcp plane. The meaning of colors is the same as in Fig. 1. In order to understand the role of the key parameter αη, we show how the three neutrino mixing angles and the CP violating Dirac phase depend on αη in Figs At first, in Fig. 3, we show the prediction of sin θ13 versus αη where the 3 σ data is taken as input except for sin θ13. The red lines denote the upper and lower bounds of the 3 σ experimental data for sin θ13. Note 1

12 D D..3 Figure 8: The prediction of mee vs. m1. The meaning of colors is the same as in Fig. 1. Figure 7: The predicted Majorana phases on α1 α31 plane. The meaning of colors is the same as in Fig D sin Θ Figure 1: The Σmi dependence of the predicted δcp. The meaning of colors is the same as in Fig. 1. Figure 9: The Σmi dependence of the predicted sin θ3. The meaning of colors is the same as in Fig. 1. that sin θ13 depends on αη crucially as seen in Eq.(). As shown in Fig. 3, the observed value sin θ13 is not reproduced unless αη is larger than.7. The clear dependence between αη and the predicted sin θ3 can be seen in Fig.4. In order to reproduce the maximal mixing of θ3, αη should be larger than.1. The highly probable prediction of sin θ3 is near.47.5 for.1 αη.. The deviation from the trimaximal mixing of sin θ1 explicitly depends on αη as seen in Eq.(). We show the prediction of sin θ1 versus αη in Fig. 5. The predicted sin θ1 is almost independent of αη as far as αη.1. The αη dependence on δcp gives the characteristic prediction as shown in Fig. 6. The CP conservation δcp = is excluded in the smaller region αη.1 for the experimental data with 3 σ. By inputting the 1 σ data in Table, we obtain the prediction of δcp to be ±(5 1 ), which is almost independent of αη for αη =.1.. We show the prediction of the Majorana phases α1 and α31 in Fig. 7. While both Majorana phases are allowed in all region of 18 18, there is a clear correlation between both phases. In Fig. 8, we present the predicted mee, the effective mass for the νββ decay, versus m1 which is another key parameter in our model. The parameter m1 should be larger than 1 mev 11

13 δcp [ ] δcp [ ] sin θ1 sin θ3 Figure 1: The allowed region on sin θ1 δcp plane for η(1 ). The meaning of colors is the same as in Fig. 1. Figure 11: The allowed region on sin θ3 δcp plane for η(1 ). The meaning of colors is the same as in Fig. 1. in order to reproduce the observed mass squared differences, and it is smaller than 46 mev due to the cosmological constraint on the sum of neutrino masses [6]. In the hierarchical case of neutrino masses m1 < m m3, the predicted value mee is at most 1 mev but close to 45 mev for the degenerated neutrino masses. Next, we discuss the sum of three neutrino masses Σmi because the cosmological observation gives us a upper bound for it. We show the predicted region of Σmi sin θ3 plane in Fig. 9. The minimum of the sum of three neutrino masses Σmi is 75 mev in our model. In order to get sin θ3.5, Σmi should be larger than 85 mev. For the best fit of sin θ3 =.538, Σmi is expected to be larger than 9 mev. We show the predicted region of Σmi δcp plane in Fig. 1. The predicted δcp is smaller than 9 if Σmi is smaller than 85 mev. Thus, the cosmological observation for the sum of neutrino masses will be a crucial test of these predictions. We have neglected the next-leading terms llφs φt hu hu and llφt ηhu hu in the neutrino mass matrix of Eq.(9) because α` =.316 (.1) is small compared with αη.1. We have confirmed that those effects are small with our numerical calculation by inputting 1 σ data. Indeed, the prediction of sin θ3 δcp almost remains inside of the red curve in Fig. 1. It is also deserved to comment on the αη distribution in our numerical results. In order to remove the predictions for αη >.3 smoothly, which is about ten times larger than α` =.316, we have used the Gamma distribution for αη given in Eq. (56) of Appendix D. We have confirmed that our results are not changed even if we adopt another Gamma distribution presented in Eq.(57) of Appendix D although the number density of dots gets lower. We have also used α` =.1 which corresponds to SM in our calculations. In this case, the number density of dots significantly gets lower, but the allowed region is almost unchanged. Moreover, we have found that the allowed region is also unchanged even if we use the flat-distribution of αη in the region αη.3. Thus, our results are robust for any distribution of αη. 3. Case of a 1 singlet η We show the numerical results for a 1 singlet η briefly because the correlations of the observables appear to be weak. We show the predicted δcp versus sin θ3 in Fig. 11. The region of δcp 5 is almost excluded while the regions near ±18 are allowed. There is no correlation between 1

14 sin θ 3 and δ CP. We also show the predicted δ CP versus sin θ 1 in Fig. 1. The predicted δ CP increases as sin θ 1 decreases from the trimaximal mixing 1/3, but its correlation is rather weak. Both results in Figs. 11 and 1 are due to both mixing of (1-) and (-3) families in the charged lepton sector. Thus, the model with the 1 singlet η is less attractive than that with the 1 singlet η in light of NuFIT 3. data. 4 Summary The flavor symmetry of leptons can be examined precisely in light of the new data and the upcoming experiments [9]. We study the A 4 model with minimal parameters by using the results of NuFIT 3.. We introduce the A 4 singlet 1 or 1 flavon η which couples to both the charged lepton and neutrino sectors in the next-leading order due to the relevant Z 3 charge for η. The model with the 1 (1 ) flavon is consistent with the experimental data of m sol only for NH of neutrino masses. The key parameter is α η which is derived from the VEV of the flavon η. The parameter α η is distributed around α l =.316 (.1) in the Gamma distribution of the statistic. Our results are robust for different distribution of α η. In the case of the singlet η(1 ), α η should be larger than.7 in order to reproduce the observed value of sin θ 13. The numerical prediction of δ CP versus sin θ 3 is similar to the prediction of TM. However, our predicted region is inside of the TM boundary. The absolute value of the predicted δ CP is 6 9 for the maximal mixing θ 3 = π/4. For the best fit of sin θ 3 =.538, δ CP is in the region of The predicted sin θ 1 is also allowed around the best fit of NuFIT 3. while keeping near the maximal mixing of θ 3. Inputting the data with the 1 σ errorbar, we obtain the clear prediction of the CP violating Dirac phase to be δ CP = 5 1. The lightest neutrino mass m 1 is expected to be 1 mev 46 mev, which leads to the m ee < 45 mev. In order to get the best fit of sin θ 3 =.538, the sum of the three neutrino masses is expected to be larger than 9 mev. The cosmological observation for the sum of the neutrino masses will also provide a crucial test of these predictions. The model with η(1 ) is not attractive in light of NuFIT 3. result because the input data given in Table does not give a severe constraint for the predicted region of δ CP. We expect the precise measurement of the CP violating phase to test the model in the future. Acknowledgement This work is supported by JSPS Grants-in-Aid for Scientific Research 16J533 (YS) and 15K545, 16H86 (MT). S.K.Kang was supported by the National Research Foundation of Korea (NRF) grants (9-8356, 17K1A3A7A91643, 17R1AB46338). Appendix A Lepton mixing matrix Supposing neutrinos to be Majorana particles, the PMNS matrix U PMNS [4, 5] is parametrized in terms of the three mixing angles θ ij (i, j = 1,, 3; i < j), one CP violating Dirac phase δ CP 13

15 and two Majorana phases α 1, α 31 as follows: c 1 c 13 s 1 c 13 s 13 e iδ CP 1 U PMNS = s 1 c 3 c 1 s 3 s 13 e iδ CP c 1 c 3 s 1 s 3 s 13 e iδ CP s 3 c 13 e i α 1, s 1 s 3 c 1 c 3 s 13 e iδ CP c 1 s 3 s 1 c 3 s 13 e iδ CP c 3 c 13 e i α 31 (5) where c ij and s ij denote cos θ ij and sin θ ij, respectively. The rephasing invariant CP violating measure, Jarlskog invariant [3], is defined by the PMNS matrix elements U αi. It is written in terms of the mixing angles and the CP violating phase as: J CP = Im [ U e1 U µ U eu µ1] = s3 c 3 s 1 c 1 s 13 c 13 sin δ CP, (6) where U αi denotes the each component of the PMNS matrix. There are also other invariants I 1 and I associated with Majorana phases [31]-[34], ( I 1 = Im [Ue1U e ] = c 1 c 1 c α1 ) ( 13 sin, I = Im [U α31 ) e1u e3 ] = c 1 s 13 c 13 sin δ CP We calculate δ CP, α 1 and α 31 with these relations.. (7) B Multiplication rule of A 4 group We use the multiplication rule of the A 4 triplet as follow: a 1 a a 3 3 b 1 b b 3 3 = (a 1 b 1 + a b 3 + a 3 b ) 1 (a 3 b 3 + a 1 b + a b 1 ) 1 (a b + a 1 b 3 + a 3 b 1 ) 1 1 a 1 b 1 a b 3 a 3 b a 3 b 3 a 1 b a b 1 3 a b a 1 b 3 a 3 b a b 3 a 3 b a 1 b a b 1 a 3 b 1 a 1 b 3 3, 1 1 = 1, 1 1 = 1, 1 1 = 1, 1 1 = 1. (8) More details are shown in the review [15, 16]. C Charged lepton and neutrino mass matrices The charged lepton mass matrix is given by the multiplication rule of A 4 in Appendix B as follows: y e λ 4 y τα η y e λ 4 y µα η λ M l = v d α l y eα η λ 4 y µ λ for η (1 ), v d α l y µ λ y τα η for η (1 ), y µα η λ y τ y eα η λ 4 y τ (9) 14

16 where α l, α η and λ are written in terms of the VEVs of φ T, η and Θ: α l φ T Λ = v T Λ, α η η Λ = q Λ, λ Θ Λ. (3) The coefficients y i and y i are order one parameters. The left-handed mixing matrix of the charged lepton is derived from the diagonalization of U l M l M l U l. The diagonalizing matrix U l for the charged lepton is given as follows: 1 y µ y τ y U l µ y τ α y τ η y τ α η ( y µ y τ y µ y τ ) αη y 1 µy µ y τ α η λ 4 for η (1 ), ( y τ y τ ) y µ α η y µ y τ y τ αη U l ( y µ y µ ) α η 1 ( y µ y τ y µ y τ ) αη y µ y µ α η y τ y τ y µ y τ ( yµ y τ ) αηλ 4 for η (1 ). (31) y τ y τ α η ( y τ y τ ) α η 1 The mass eigenvalues of the charged leptons are given in a good approximation: m e = y e α l λ 4 v d, m µ = y µ α l λ v d, m τ = y τ α l v d, (3) where the Yukawa couplings are of order one. For η (1 ), the neutrino mass matrix is given as: M ν = a 1 + b c 1 + d 1, (33) where a + 3b = is satisfied. The coefficients a, b, c and d are given in terms of the Yukawa couplings and the VEVs of flavons: with a = y Sα ν Λ v u, b = y Sα ν 3Λ v u, c = y ξα ξ Λ v u, d = y 7α ξ α η Λ v u, (34) α ν φ S Λ = v S Λ, α ξ ξ Λ = u Λ. (35) Taking a to be real without loss of generality, we reparametrize them as follows: a a, c c e iφc, d d e iφ d, (36) where a, c and d are real parameters and φ c, φ d are CP violating phases. For η (1 ), we get M ν = a 1 + b c 1 + d 1. (37)

17 We obtain the neutrino mass matrix in the TBM basis by rotating it with V TBM, V TBM = , (38) as follows: a + ce iφc d ˆM ν = VT T eiφ d 3 deiφ d BMM ν V T BM = ce iφc + de iφ d, (39) 3 deiφ d a ce iφc + d eiφ d where upper (lower) sign in front of (1,3) and (3,1) components correspond to the assignment of 1 and 1 for η, respectively. Next, we consider (1, 1) (1, 3) ˆM ν ˆM ν = ce iφc + de iφ d, (4) (1, 3) (3, 3) where (1, 1) = a + c + d + ac cos φ c cd cos(φ c φ d ) ad cos φ d, (3, 3) = a + c + d ac cos φ c cd cos(φ c φ d ) + ad cos φ d, (1, 3) = 3 [ad cos φ d + icd sin(φ c φ d )]. (41) We obtain the neutrino mass eigenvalues for NH as follows: m 1 = a + c + d cd cos(φ c φ d ) 3c d sin (φ c φ d ) + 4a (c cos φ c + d cos φ d cd cos φ c cos φ d ), m 3 = a + c + d cd cos(φ c φ d ) + 3c d sin (φ c φ d ) + 4a (c cos φ c + d cos φ d cd cos φ c cos φ d ), m = c + d + cd cos(φ c φ d ). (4) The ˆM ν ˆM ν is diagonalized by the (1-3) family rotation as: U ν ( ˆM m ν ˆM ν ) U ν 1 = m, (43) m 3 where cos θ sin θe iσ U ν = 1. (44) sin θe iσ cos θ The θ and σ are given in terms of parameters in the neutrino mass matrix: tan θ = 3 d a cos φ d + c sin (φ c φ d ) a(d cos φ d c cos φ c ), σ = c sin(φ c φ d ) a cos φ d. (45) 16

18 The parameters a, c and d are written in terms of m 1 and α η. As seen in Eq.(34), the parameter d is related with c as d c = y 7 α η αη ν, (46) y ξ where y 7 and y ξ are of order one coefficients. On the other hand, a and c are given in terms of m 1, α ν η, m 31 m 3 m 1 and m 1 m 3 m 1 since we have the following relations in Eq.(4): 1 4 ( m 31) = 3c d sin (φ c φ d ) + 4a (c cos φ c + d cos φ d cd cos φ c cos φ d ), m 1 = c + d + cd cos(φ c φ d ) m 1. (47) Then, putting m atm = m 31 and m sol = m 1, c = m sol + m (α ν η) + α ν η cos(φ c φ d ), a = 1 16c m atm 1c 4 (α ν η) sin (φ c φ d ) cos φ c + (α ν η) cos φ d α ν η cos φ c cos φ d, (48) where m 1 and α ν η are free parameters as well as φ c and φ d. We comment on the case of IH of neutrino masses. In this case, the neutrino mass eigenvalues are given as m 1 = a + c + d cd cos(φ c φ d ) + 3c d sin (φ c φ d ) + 4a (c cos φ c + d cos φ d cd cos φ c cos φ d ), m 3 = a + c + d cd cos(φ c φ d ) 3c d sin (φ c φ d ) + 4a (c cos φ c + d cos φ d cd cos φ c cos φ d ), m = c + d + cd cos(φ c φ d ), (49) where m 1 and m 3 are exchanged each other in Eq.(4). Then, m sol is given as m sol = m m 1 = 3cd cos(φ c φ d ) a 3c d sin (φ c φ d ) + 4a (c cos φ c + d cos φ d cd cos φ c cos φ d ). (5) It is impossible to reproduce the observed value of m sol since a c and c d in our model as seen in Eq.(34). Indeed, d/c is expected to be.1. in our numerical analysis. Finally, the PMNS matrix is given as U PMNS = U l V TBM U ν P, (51) where P is the diagonal phase matrix originated from the Majorana phases. The P is obtained from P U ν ˆMν U T ν P = diag{m 1, m, m 3 }, (5) where m 1, m and m 3 are real positive. The effective mass for the νββ decay is given as follows: m ee = m1 Ue1 + m Ue + m 3 Ue3. (53) 17

19 D Distribution of α η The magnitude of the parameter α l is determined by the tau-lepton mass as seen in Eq.(4). The key parameter α η is related with α l through the vacuum structure as discussed in Eq.(4): α η = λ 1 g4 α l. (54) λ 3g 3 The coefficients λ 1() and g 3(4) are of order one. Then, the factor in front of α l in Eq.(54) could be O(1). We scan α η by using Eq.(54) after fixing α l in the statistical approach. For this purpose, we use the Gamma distribution which is available to find the distribution of the order one parameter: Taking γ = 1 with α = 3/, µ = and β =, we obtain f = (x µ) (αγ 1) e ( x µ β ) γ. (55) f = x e 1 x, (56) which is equivalent to the χ distribution. When we take γ = with α = 1, µ = and β =, we obtain f = x e 1 x, (57) which damps as like the Gaussian distribution at the large x. It is easy to check that f is maximal at x = 1 and f = at x = for both the two types of Gamma distribution. We obtain the distribution of α η by multiplying α l by f, which is used in our numerical calculations. We show the distribution of α η in Figs. 13 and 14 for α l =.316 (MSSM tan β = 3) and α l =.1 (SM) in the case of the distributions of Eqs.(56) and (57), respectively. 35 frequency frequency Α Η Figure 13: α η distribution for α l =.316 (blue) and α l =.1 (red) in Eq.(56) (α = 3/, β =, γ = 1, µ = ). Figure 14: α η distribution for α l =.316 (blue) and α l =.1 (red) in Eq.(57) (α = 1, β =, γ =, µ = ). 18

20 References [1] K. Abe et al. [TK Collaboration], arxiv: [hep-ex]. [] TK report, August 4, 17. [3] P. Adamson et al. [NOvA Collaboration], Phys. Rev. Lett. 118 (17) no.3, 3181 [arxiv: [hep-ex]]. [4] A. Radovic, Latest oscillation results from NOvA. Joint Experimental-Theoretical Physics Seminar, Fermilab, USA, January 1, 18. [5] NuFIT 3. (18), [6] F. P. An et al. [DAYA-BAY Collaboration], Phys. Rev. Lett. 18 (1) [arxiv: [hep-ex]]. [7] J. K. Ahn et al. [RENO Collaboration], Phys. Rev. Lett. 18 (1) 1918 [arxiv:14.66 [hep-ex]]. [8] P. F. Harrison, D. H. Perkins, W. G. Scott, Phys. Lett. B 53 () 167 [hep-ph/74]. [9] P. F. Harrison, W. G. Scott, Phys. Lett. B 535 () [hep-ph/39]. [1] E. Ma and G. Rajasekaran, Phys. Rev. D 64, 1131 (1) [arxiv:hep-ph/1691]. [11] K. S. Babu, E. Ma and J. W. F. Valle, Phys. Lett. B 55, 7 (3) [arxiv:hep-ph/69]. [1] G. Altarelli and F. Feruglio, Nucl. Phys. B 7 (5) 64 [hep-ph/54165]. [13] G. Altarelli and F. Feruglio, Nucl. Phys. B 741 (6) 15 [hep-ph/5113]. [14] G. Altarelli and F. Feruglio, arxiv:1.11 [hep-ph]. [15] H. Ishimori, T. Kobayashi, H. Ohki, Y. Shimizu, H. Okada and M. Tanimoto, Prog. Theor. Phys. Suppl. 183 (1) 1 [arxiv: [hep-th]]. [16] H. Ishimori, T. Kobayashi, H. Ohki, H. Okada, Y. Shimizu and M. Tanimoto, Lect. Notes Phys. 858 (1) 1, Springer. [17] S. F. King, A. Merle, S. Morisi, Y. Shimizu and M. Tanimoto, arxiv: [hep-ph]. [18] Y. Shimizu, M. Tanimoto and A. Watanabe, Prog. Theor. Phys. 16 (11) 81 [arxiv: [hep-ph]]. [19] W. Grimus and L. Lavoura, JHEP 89 (8) 16 [arxiv:89.6 [hep-ph]]. [] C. H. Albright, A. Dueck and W. Rodejohann, Eur. Phys. J. C 7 (1) 199 [arxiv: [hep-ph]]. [1] W. Rodejohann and H. Zhang, Phys. Rev. D 86 (1) 938 [arxiv:17.15]. 19

21 [] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B 147 (1979) 77. [3] G. Altarelli, F. Feruglio and C. Hagedorn, JHEP 83 (8) 5 [arxiv:8.9 [hepph]]. [4] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 8 (196) 87. [5] B. Pontecorvo, Sov. Phys. JETP 6 (1968) 984 [Zh. Eksp. Teor. Fiz. 53 (1967) 1717]. [6] E. Giusarma, M. Gerbino, O. Mena, S. Vagnozzi, S. Ho and K. Freese, Phys. Rev. D 94 (16) no.8, 835 [arxiv: [astro-ph.co]]. [7] Y. Shimizu, M. Tanimoto and K. Yamamoto, Mod. Phys. Lett. A 3 (15) 155 [arxiv: [hep-ph]]. [8] S. K. Kang and M. Tanimoto, Phys. Rev. D 91 (15) no.7, 731 [arxiv: [hep-ph]]. [9] S. T. Petcov and A. V. Titov, arxiv: [hep-ph]. [3] C. Jarlskog, Phys. Rev. Lett. 55 (1985) 139. [31] S. M. Bilenky, S. Pascoli and S. T. Petcov, Phys. Rev. D 64 (1) 531 [hep-ph/165]. [3] J. F. Nieves and P. B. Pal, Phys. Rev. D 36 (1987) 315; Phys. Rev. D 64 (1) 765 [hep-ph/1535]. [33] J. A. Aguilar-Saavedra and G. C. Branco, Phys. Rev. D 6 () 969 [hep-ph/75]. [34] I. Girardi, S. T. Petcov and A. V. Titov, Nucl. Phys. B 911 (16) 754 [arxiv: [hep-ph]].