Update on two-zero textures of the Majorana neutrino mass matrix in light of recent T2K, Super-Kamiokande and NO$\nu$A results

The latest results from atmospheric and accelerator neutrino experiments indicate that the normal neutrino mass ordering $m^{}_1<m^{}_2<m^{}_3$, a maximal leptonic CP-violating phase $\delta = 270^\circ$ and the second octant of neutrino mixing angle $\theta^{}_{23}>45^\circ$ are favored. In light of new experimental results, we update previous phenomenological studies on two-zero textures of the Majorana neutrino mass matrix $M^{}_\nu$, in the flavor basis where the charged-lepton mass matrix $M^{}_l$ is diagonal. When the $1\sigma$ ranges of neutrino mixing parameters are taken into account, only four (i.e., ${\bf A}^{}_{1, 2}$ and ${\bf B}^{}_{2,4}$) among seven two-zero patterns of $M^{}_\nu$ show the aforementioned features of neutrino mass spectrum, mixing angle $\theta^{}_{23}$ and CP-violating phase $\delta$, and thus are compatible with the latest neutrino oscillation data. The correlative relations among neutrino masses and mixing parameters have been derived analytically for these four patterns, and the allowed regions of neutrino mixing angles and the CP-violating phase are also given. Possible realizations of four viable two-zero textures via non-Abelian discrete flavor symmetries are discussed.

1 Introduction particular, only the parameter space around δ = 90 • or satisfiying θ 23 < 45 • is concentrated on in Ref. [12], which is now inappropriate in light of recent T2K, Super-Kamiokande and NOνA results. Furthermore, we have derived a novel analytical formula for the correlation among neutrino mixing parameters, which is applicable to a special set of two-zero textures and turns out to be much simpler than that existing in the literature. Among seven two-zero textures, we find that only four (i.e., A 1,2 and B 2,4 ) remain consistent with the latest neutrino oscillation data if the 1σ ranges of mixing parameters are taken into account.
The rest of this paper is structured as follows. In Sec. 2, the important and basic formulas for twozero textures are summarized. Sec. 3 is devoted to analytical and approximate results of theoretical predictions, and a complete numerical analysis of the allowed parameter space. We have found a novel and useful formula that has not yet been noticed in previous works. It turns out that only four out of seven patterns (i.e., A 1,2 and B 2,4 ) are compatible with the latest neutrino oscillation data. The precise determination of δ and θ 23 could further single out a unique pattern of two-zero textures, or exclude all of them eventually. In Sec. 4, we present some examples to illustrate how to realize those two-zero textures through non-Abelian discrete flavor symmetries. Finally, our main conclusions are summarized in Sec. 5.

Two-zero Textures
Since neutrino masses cannot be accommodated within the standard model of elementary particles (SM), one can regard the SM as an effective theory at the electroweak scale and introduce higher-dimensional operators, which fully comply with the SM gauge symmetries. As already pointed out by Weinberg [15], the lowest-dimensional operator composed of lepton and Higgs doublets is (ℓ L Φ)(Φ T ℓ C L ) of dimension five, which is unique and leads to lepton number violation by two units. After spontaneous gauge symmetry breaking, neutrinos acquire a tiny Majorana mass term from this five-dimensional operator. Such an operator can actually be realized in renormalizable models, e.g., the seesaw models, in which three singlet right-handed neutrinos [16,17,18,19,20] or one triplet scalar [21,22,23,24,25,26] or three triplet fermions [27] are introduced and the mass scale of these new particles is extremely high.
In the seesaw framework, the mass terms relevant for lepton mass spectra and flavor mixing at lowenergy scales can be written as where M l and M ν stand for the charged-lepton and Majorana neutrino mass matrix, respectively. In the flavor basis, where M l = diag m e , m µ , m τ is diagonal, the lepton flavor mixing matrix V stems from the diagonalization of a complex and symmetric matrix M ν , namely, Here m α (for α = e, µ, τ ) and m i (for i = 1, 2, 3) denote respectively the charged-lepton and neutrino masses. For later convenience, we decompose the mixing matrix as V = U · P , where P ≡ diag{e iρ , e iσ , 1} with {ρ, σ} being two Majorana CP-violating phases, and U can be parametrized in terms of three mixing angles {θ 12 , θ 13 , θ 23 } and one Dirac CP-violating phase δ, namely, with s ij ≡ sin θ ij and c ij ≡ cos θ ij for ij = 12, 13, 23. The three mixing angles {θ 12 , θ 13 , θ 23 }, two neutrino mass-squared differences ∆m 2 21 ≡ m 2 2 − m 2 1 and ∆m 2 31 ≡ m 2 3 − m 2 1 , and the Dirac CP-violating phase δ can be determined from neutrino oscillation experiments. In addition, the constraints on the absolute neutrino masses can be obtained from both beta-decay and neutrinoless double-beta decay experiments, while the latter could also provide us with some useful information on the two Majorana CP-violating phases. The latest observations of cosmic microwave background by the Planck collaboration have placed a restrictive bound on the sum of neutrino masses Σ = m 1 + m 2 + m 3 < 0.23 eV at the 95% confidence level [28]. In this context, any useful experimental information on the mixing matrix V and neutrino masses m i (for i = 1, 2, 3) can be implemented to reconstruct the Majorana neutrino mass matrix via where λ 1 ≡ m 1 e 2iρ , λ 2 ≡ m 2 e 2iσ , and λ 3 ≡ m 3 . On the other hand, any constraints on M ν arising from flavor symmetries or empirical assumptions will result in specific predictions for neutrino mass spectrum and mixing parameters, which are ready to be tested in neutrino oscillation experiments. Now we follow the second way and assume two texture zeros in M ν , which were first considered by Frampton, Glashow and Marfatia in Ref. [29] and thoroughly studied by Xing in a series of papers [30,31,32], and by some other authors [33,34,35]. The two-zero textures have been reexamined after the experimental discovery of a nonzero θ 13 in Refs. [12,13,14,36,37,38,39]. The predictions from two-zero textures and their experimental tests can be found in previous works, however, we summarize the important formulas in this section and establish our notations. Since M ν is symmetric, there are totally six independent matrix elements. If two of them are assumed to be zero, i.e., (M ν ) ab = (M ν ) cd = 0 for ab = cd, we obtain From Eq. (4), one can immediately extract the neutrino mass ratios [30,31,32,12] ξ ≡ m 1 and two Majorana CP-violating phases Now it is instructive to count the number of free parameters in this scenario. The physical parameters include three neutrino masses {m 1 , m 2 , m 3 }, three mixing angles {θ 12 , θ 13 , θ 23 } and three CP-violating phases {δ, ρ, σ}. Given the four constraints or relations in Eqs. (5) and (6), we are left with five degrees of freedom, which can be fixed by two neutrino mass-squared differences ∆m 2 21 and ∆m 2 31 , and three mixing angles {θ 12 , θ 13 , θ 23 }. Therefore, the absolute neutrino mass (e.g., m 3 ) and CP-violating phases {δ, ρ, σ} can be determined by the available neutrino oscillation data on {∆m 2 21 , ∆m 2 31 } and {θ 12 , θ 13 , θ 23 }. To be explicit, neutrino masses can be found by using and the Dirac CP-violating phase δ can be fixed by where the dependence of ζ and ξ on δ is implied in Eq. (5). The approximate analytical results of {δ, ρ, σ} can be found in Refs. [30,31,32,12]. In Ref. [12], the global-fit results of neutrino oscillation data at that time have been used to check the viability of two-zero textures of M ν . It turns out that the following seven patterns survive current neutrino oscillation data: where the nonzero elements are in general complex. However, it is always possible to make three nonzero matrix elements real by redefining the phases of three charged-lepton fields. In this case, we have five free model parameters (e.g., three real elements {a, b, c} and one complex element d), as it should be consistent with our parameter counting. Motivated by the recent results from T2K, Super-Kamiokande and NOνA, we reexamine those seven patterns by focusing on whether NO, δ = 270 • and θ 23 > 45 • are allowed simultaneously. If this is the case for a given pattern, we further completely explore the viable regions of neutrino masses and flavor mixing parameters. Unfortunately, the latest global-fit analysis of neutrino oscillation data has not yet included the preliminary results from T2K, Super-Kamiokande and NOνA, although partial joint analyses of reactor and long-baseline accelerator neutrino data have been performed in Refs. [10,11]. For comparison, we collect the global-fit data from Ref. [40] in Table 1. Due to the lack of a full statistical analysis, we quote the 1σ intervals of ∆m 2 21 and θ 12 from Table 1 in the NO case, while the 68% Bayesian credible regions of ∆m 2 32 , θ 23 , θ 13 and δ are taken from Ref. [7], in which a Bayesian analysis of both disappearance and appearance data from T2K and reactor neutrino data has been performed. The oscillation data used in our numerical calculations are summarized as follows and The best-fit values of three mixing angles (θ 12 = 33 • , θ 13 = 9 • , θ 23 = 47 • ), two mass-squared differences (∆m 2 21 = 7.50 × 10 −5 eV 2 , ∆m 2 31 = 2.58 × 10 −3 eV 2 ), and the CP-violating phase δ = 270 • will also be used for illustration. The best-fit values, together with the 1σ, 2σ and 3σ intervals, for three neutrino mixing angles {θ 12 , θ 13 , θ 23 }, two mass-squared differences {∆m 2 21 , ∆m 2 31 or ∆m 2 32 } and the Dirac CP-violating phase δ from a global analysis of current experimental data [40]. Two independent global-fit analyses can be found in Refs. [41,42], which are in perfect agreement with the results presented here at the 3σ level.

Analytical and Numerical Results
Before proceeding with a numerical analysis, we briefly comment on the analytical results presented in Refs. [30,31,32,12], where the standard parametrization of U in Eq.
(2) has been inserted into Eqs. (5) and (6) to derive the neutrino mass ratios and CP-violating phases in terms of mixing angles. Of course, this is a correct and straightforward way to achieve an analytical understanding of any numerical results. However, for the patterns A 1 , A 2 , B 3 and B 4 in Eqs. (9) and (10), we can reach the analytical results in a novel and simple way. The two zeros in these four patterns are located in the same row or column of M ν , namely, (M ν ) aa = (M ν ) ab = 0 for a = b. In this case, Eq. (4) can be recast into or equivalently Then, the orthogonality conditions of the unitary matrix U bring us another array of equations, which can be put in a form similar to Eq. (14), i.e., which is generally satisfied for the allowed values of mixing parameters, the solutions to Eq. (14) and Eq. (15) should be identical. By identifying these two solutions, we obtain a set of new relations from which one can extract two neutrino mass ratios {ξ, ζ} and two Majorana CP-violating phases {ρ, σ} in the same manner as in Eqs. (5) and (6). Obviously, the formulas in Eq. (16) are much simpler than those derived in previous works [30,31,32,12].
As proved in Ref. [12], the predictions from A 1 and A 2 are related by a permutation symmetry corresponding to θ 23 → 90 • − θ 23 , and likewise for B 3 and B 4 . Thus, we apply Eq. (16) only to A 1 and B 4 , and discuss their phenomenological implications.
For this pattern, we can specify the flavor indices in Eq. (16) as a = e, b = µ and c = τ . Consequently, one arrives at s 12 s 23 c 12 c 23 e iδ − s 13 e −2iδ , which are identical to those found in Ref. [12] up to different definitions of two Majorana CP-violating phases. To the leading order of s 13 , it has been found [30] ξ ≈ tan θ 23 tan θ 12 sin θ 13 , ζ ≈ tan θ 23 cot θ 12 sin θ 13 .
Given the best-fit values of ∆m 2 21 = 7.50 × 10 −5 eV 2 and ∆m 2 31 = 2.58 × 10 −3 eV 2 below Eq. (12), one obtains R ν = 0.029. Using the best-fit values of three mixing angles, we immediately find that the value on the right-hand side of Eq. (20) is 0.014, which is smaller than R ν by a factor of two. Therefore, the best-fit values seem not to be consistent with the pattern A 1 .
As we have mentioned before, the results for A 2 can be derived from those for A 1 by a replacement θ 23 → 90 • −θ 23 . In this case, it is easy to verify that A 2 is compatible with NO as well. For δ = 270 • , tan 2 θ 23 on the right-hand side of Eq. (20) will be replaced by cot 2 θ 23 , implying a even smaller R ν . Hence, we conclude that both A 1 and A 2 are consistent with NO, but they may run into problems if δ = 270 • is experimentally confirmed. However, if the current uncertainty on δ is taken into account, both A 1 and A 2 are not yet excluded.
Using the best-fit values of three mixing angles and two neutrino mass-squared differences, one can get δ ≈ 269.6 • or 90.4 • , which is very close to the maximal CP-violating phase. The absolute neutrino mass scale is also fixed by neutrino mixing angles and mass-squared differences as where θ 23 = 47 • and ∆m 2 31 = 2.58×10 −3 eV 2 have been input. Obviously, neutrino masses are nearly degenerate, and their sum is Σ ≈ m 3 (1 + 2 cot 2 θ 23 ) ≈ 0.27 eV, which is already in tension with the current cosmological upper bound. In summary, the allowed region of the Dirac CP-violating phase is severely constrained, namely, only δ ≈ 90 • or 270 • is viable. If θ 23 → 45 • , we have δ → 90 • or 270 • and m i → ∞, so the upper bound on absolute neutrino masses forbids the maximal mixing angle θ 23 = 45 • and maximal CP-violating phase δ = 90 • or 270 • [35].
For other cases, Eq. (16) is no longer valid, and the analytical formulas for λ 1 /λ 3 and λ 2 /λ 3 become very complicated, as one can see from Refs. [30,31,12]. To complete the list of two-zero textures that survive the latest neutrino oscillation data, we quote the analytical results for B 2 below from which one can obtain ξ ≈ ζ ≈ cot 2 θ 23 at the leading order. Hence, B 2 is compatible with NO for θ 23 > 45 • , whereas ξ ≈ ζ ≈ tan 2 θ 23 is obtained for B 1 . The latter pattern is then excluded if the hints on NO and θ 23 > 45 • are further confirmed by future neutrino oscillation data. At the next-to-leading order, the neutrino mass difference turns out to be indicating cos δ > 0 due to m 2 > m 1 . Moreover, the CP-violating phase is given by cos δ ≈ − sin 2θ 12 2 sin θ 13 tan 2θ 23 R ν .
Taking the best-fit values of neutrino mixing parameters, we get δ ≈ 270.3 • or 89.7 • . The absolute neutrino mass is the same as given in Eq. (24), and a nearly-degenerate neutrino mass spectrum is expected. To distinguish between B 2 and B 4 , we have to measure the CP-violating phase as precisely as possible, e.g., an accuracy less than 1 • is required if δ ≈ 270 • is confirmed.
Now we perform a thorough numerical analysis of all two-zero textures of M ν by using the 1σ ranges of neutrino mixing angles, Dirac CP-violating phase and neutrino mass-squared differences in Eqs. (11) and (12). The main strategy of our numerical analysis is as follows. (1) We randomly generate the values of {θ 12 , θ 13 , θ 23 } and δ in their 1σ ranges, and then calculate two neutrino mass ratios {ξ, ζ} and R ν . (2) Several constraints are then chosen as criteria for whether a two-zero pattern is compatible with neutrino oscillation data. First of all, the fact of m 2 > m 1 and NO, or equivalently ξ < ζ < 1, should be reproduced. Then, a pair of values for {∆m 2 21 , ∆m 2 31 } are randomly generated and the ratio ∆m 2 21 /∆m 2 31 is compared to R ν that has been computed from neutrino mixing angles and CP-violating phase in the first step. The difference between those two is required to be vanishingly small (e.g., < 10 −4 ). The final numerical results are summarized in Fig. 1, and some comments are in order: 1. Out of seven two-zero textures, only A 1,2 and B 2,4 are found to be consistent with the neutrino oscillation data in Eqs. (11) and (12). In Fig. 1, the allowed regions of {θ 12 , θ 13 , θ 23 } are plotted with respect to those of δ for A 1 (red stars), A 2 (blue squares), B 2 (pink pluses), and B 4 (green crosses), where the best-fit value δ = 270 • and the 1σ interval 202 • < δ < 338 • are represented by a dashed line and a shaded area, respectively. Different from Ref. [12], the whole range of δ is now considered.
2. The patterns A 1 and A 2 are well separated in the parameter space of δ. In particular, the values of δ for A 1 are confined to a narrow area around δ = 180 • , i.e., the CP-conserving limit. Therefore, the sign of cos δ (i.e., cos δ < 0 for A 1 and cos δ > 0 for A 2 ) is of crucial importance to discriminate between these two patterns. This is also the case for B 2 and B 4 , as cos δ > 0 for B 2 and cos δ < 0 for B 4 . However, in the latter case, a very precise determination of δ (e.g., less than 1 • ) is required.
3. As one can observe from the lower plot of Fig. 1, a joint measurement of θ 23 and δ is very useful to tell A 2 apart from B 2,4 . For A 2 , a nearly maximal CP-violating phase δ ≈ 90 • or 270 • is reachable only if θ 23 takes a value close to the upper bound of its 1σ range. Another way to distinguish between A 1,2 and B 2,4 is to discover neutrinoless double-beta decays, which never take place for the former cases because of a vanishing effective neutrino mass m ββ = | (M ν ) ee | = 0. But for B 2 and B 4 , we have m ββ = | (M ν ) ee | ≈ m 3 ≈ 0.1 eV, which is likely to be probed in future experiments of neutrinoless double-beta decays [43,44,45,46].
The Dirac CP-violating phase δ serves as an important discriminator for the two-zero textures of M ν . Although it seems difficult to further differentiate between B 2 and B 4 due to their similar implications at low energies, their high-energy phenomenologies in a realistic flavor model may be quite different.

Symmetry Realization
It is well known that the texture zeros in the Majorana neutrino mass matrix can always be realized by imposing Abelian flavor symmetries on the generic Lagrangian of neutrino mass models [47]. For instance, it has been demonstrated in Ref. [12] that cyclic symmetries Z n can be used in the framework of type-II seesaw models [21,22,23,24,25,26] to achieve all the viable two-zero textures of M ν . In general, however, it is a nontrivial task to derive two-zero textures by using non-Abelian discrete flavor symmetries.
A predictive neutrino mass model with an A 4 flavor symmetry has been proposed in Ref. [48], where the patterns B 1 and B 2 can be obtained. In this section, we recapitulate the essential idea to produce the two-zero texture of B 2 in Ref. [48], and then propose new scenarios to realize A 1 , A 2 and B 4 . In order to accommodate tiny neutrino masses, one can extend the SM by introducing three righthanded neutrinos N iR (for i = 1, 2, 3) and a Majorana mass term for them. To further explain lepton flavor mixing, we can implement the non-Abelian discrete flavor symmetry A 4 [49,50] and assign relevant fermion and scalar fields into suitable representations of the symmetry group. The assignments of fermion and scalar fields under the full symmetry SU (2) where the lepton doublets ℓ αL (for α = e, µ, τ ) are A 4 singlets, three Higgs doublets Φ i form an A 4 triplet, and one Higgs triplet ∆ is an A 4 singlet. Therefore, the gauge-and A 4 -invariant Lagrangian relevant for lepton masses and flavor mixing reads [48] − L = +a l ℓ eL ( After the scalar fields acquire their vacuum expectation values Φ i = v and ∆ = u, the charged-lepton mass matrix M l and Dirac neutrino mass matrix M D can be obtained from the first and second lines in Eq. (29), respectively. The effective neutrino mass matrix M ν receives one contribution from the type-I seesaw formula −M D M T D /M , and the other from the triplet Higgs. In total, we have which can be identified as B 2 in Eq. (9) with four nonzero matrix elements a = −(a ν v) 2 /M , b = a ∆ u, c = b ∆ u, and d = −b ν c ν v 2 /M . As shown in Ref. [48], if the Higgs triplet is assigned as 1 ′ , the pattern B 1 can be produced. In this model, however, it is impossible to realize any other viable two-zero textures.
To illustrate how to produce the other two-zero textures, we consider the minimal type-(I+II) seesaw model [51], where only one right-handed neutrino N R and one Higgs triplet ∆ have been introduced to account for both tiny neutrino masses and baryon number asymmetry in our Universe. Now this minimal scenario is further extended, and the particle content and assignments under A 4 symmetry are where two Higgs doublets ϕ and φ have been added as A 4 singlets. Except for a proper rearrangement of ℓ eL (ΦE R ) 1 ′′ , ℓ µL (ΦE R ) 1 and ℓ τ L (ΦE R ) 1 ′ , the charged-lepton sector remains unchanged and the mass matrix M l can be diagonalized just by rotating the right-handed charged-lepton fields in the flavor space, which has no impact on the lepton flavor mixing. The gauge-and A 4 -invariant Lagrangian relevant for neutrino masses becomes After the spontaneous symmetry breaking, the Dirac neutrino mass matrix turns out to be a 3 × 1 column vector and is given by which is exactly the pattern A 1 in Eq. (9) with a = a ∆ u, b = b ∆ u−a 2 D /M , c = −a D b D /M and d = −b 2 D /M . In a similar way, one can verify that the pattern A 2 can be achieved by assigning ∆ as 1 ′′ under A 4 . Finally, to obtain B 4 , we need to make a modification of Eq. (31), i.e., leading to a Dirac neutrino mass matrix M D = (a D , b D , 0) T and thus the Majorana neutrino mass matrix The phenomenological implications of these patterns for neutrino oscillation experiments have been explored in previous sections. As we have seen, it is extremely difficult to distinguish between B 2 and B 4 via neutrino oscillation experiments if the future oscillation data confirms a nearly-maximal CP-violating phase δ ≈ 270 • . For instance, the accuracy of the determination of CP-violating phase is required to be below 1 • for this purpose. In a specific model with flavor symmetries, however, new particles or interactions definitely cause additional observable effects that serve as a discriminator. We have not tried to write down the general scalar potential that is invariant under both gauge and flavor symmetries, and investigate whether the desired vacuum expectation values can be reached. This issue can be addressed in a similar way as in a number of works on applying the A 4 flavor symmetry to lepton flavor mixing [52,53,54]. Moreover, it is also interesting to examine if a successful leptogenesis works in our scenarios to explain the observed baryon number asymmetry in the Universe [55,56]. These important questions deserve further dedicated studies, which are beyond the scope of our paper and left for future work.

Summary
In view of the latest results from atmospheric (Super-Kamiokande) and accelerator (T2K and NOνA) neutrino oscillation experiments, which hint at the normal neutrino mass ordering m 1 < m 2 < m 3 , the maximal CP-violating phase δ = 270 • and the second octant θ 23 > 45 • , we reexamine the twozero textures of the Majorana neutrino mass matrix. By using the 1σ ranges of neutrino oscillation parameters, we have found that only four patterns with two texture zeros (i.e., A 1,2 and B 2,4 ) pass all current experimental constraints, and they can be further tested when the precision measurements of neutrino oscillation parameters are available.
For the patterns A 1,2 and B 3,4 , we have derived concise and useful formulas for neutrino mass ratios and Majorana CP-violating phases in Eq. (16), which have not yet been noticed in the literature. The simplicity of these formulas stems from the usage of unitarity conditions, and provides an excellent explanation for why the expressions in these four cases are much simpler than those in the cases of B 1,2 and C. A detailed numerical analysis has also been carried out to check which two-zero texture is compatible with the latest experimental data. Our numerical results show that A 1 , A 2 , B 2 and B 4 survive the data at the 1σ level, and precise determination of θ 23 and δ is of crucial importance to discriminate one pattern from another. Some distinct features should be mentioned: (1) The sign of cos δ is positive for A 2 and B 2 , while negative for A 1 and B 4 ; (2) Only nearly-maximal CP-violating phase δ ≈ 90 • or 270 • is allowed for B 2 and B 4 , and the deviations are less than a few degrees; (3) The neutrino masses are nearly degenerate in the cases of B 2 and B 4 , and substantial decay rates for neutrinoless double-beta decays are expected, while these decays do not occur in the cases of A 1 and A 2 . Hence, in addition to oscillation experiments, the constraints on absolute neutrino masses and the observations of neutrinoless double-beta decays also help verify or disprove two-zero textures. As for model building, we have shown that the non-Abelian discrete flavor symmetry A 4 can be used to realize these four patterns in the type-(I+II) seesaw models.
The ongoing and forthcoming neutrino oscillation experiments will unambiguously determine neutrino mass ordering and precisely measure the flavor mixing angles and the Dirac CP-violating phase. The tritium beta decays, neutrinoless double-beta decays and cosmological observations will probe the absolute neutrino mass scale and hopefully constrain the Majorana CP-violating phases. Any important experimental progress will give us new insights into the structure of lepton mass matrices, which can further help us explore the origin of neutrino masses and the true dynamics of flavor mixing and CP violation.  Figure 1: The allowed regions of neutrino mixing angles {θ 12 , θ 13 , θ 23 } versus the Dirac CP-violating phase δ in the scenarios of A 1 (red stars), A 2 (blue squares), B 2 (pink pluses), and B 4 (green crosses), where the shaded area corresponds to 202 • ≤ δ ≤ 338 • and the vertical dotted line to δ = 270 • , and the 1σ ranges of neutrino mixing parameters in Eqs. (11) and (12) have been used.