Generating CP violation from a modified Fridberg-Lee model

The overall characteristics of the solar and atmospheric neutrino oscillations are approximately consistent with a tribimaximal form of the mixing matrix $U$ of the lepton sector. Exact tribimaximal mixing leads to $\theta_{13}=0$. However, the results from the Daya Bay and RENO experiments have established, such that in comparison to the other neutrino mixing angles, $\theta_{13}$ is small. Moreover, the atmospheric and solar mass splitting differ by two orders of magnitude. These significant differences constitutes the great enthusiasm and main motivation for our research herein reported. Keeping the leading behavior of U as tribimaximal. We would make a response to the following questions: at some level, whether or not the small parameters such as the solar neutrino mass splitting and $U_{e3}$, which vanish in a new framework, can be interpreted as a modified FL neutrino mass model? Subsequently, a minimal single perturbation leads to nonzero values for both of them? Our minimal perturbation matrix is constructed solely from computing the third mass eigenstate, using the rules of perturbation theory. Let us point out that in contrast with~\cite{my}, this matrix is not ad hoc assumed, but is instead built following a series of steps we will outline. Also in compared to the original FL neutrino mass model which generalize it by inserting phase factors, our work is more accurate. Subsequently, we produce the following results that add new contributions to the literature: a) we obtain a realistic neutrino mixing matrix with $\delta\neq0$ and $\theta_{23}=45^\circ$; b) the solar mass splitting term is dominated by an imaginary term, which could induce the existence of Majorana neutrinos, along with explaining a large CP violation in nature.


I. INTRODUCTION
One of the remarkable observational achievements associated with neutrinos has been reported by neutrino oscillation experiments [2][3][4][5][6], which establish the non-zero neutrino masses. Concretely, that data yields information regarding neutrino masses and mixing, which can be summarized as in Table I [6].
In the standard parametrization, the lepton mixing matrix is given by [7][8][9], where c ij ≡ cos θ ij and s ij ≡ sin θ ij [for i, j = (1, 2), (1,3) and (2,3)]. The phase δ is called the Dirac phase, analogous to the CKM phase, and the phases ρ and σ are called the Majorana phases and are relevant for Majorana neutrinos. Experimental results have therefore shown that θ 13 does not vanish but is very small in comparison to the other neutrino mixing angles. This recent observation ushered the possibility of leptonic CP violation, although the CP

Parameter
The experimental data 3σ range The best fit (±1σ) ∆m 2 21   violating phase δ is not a well measured quantity. Furthermore, there is not any data about the magnitude of the Majorana phases ρ and σ. Moreover, as it has been shown, the solar mass splitting is about two orders smaller than the atmospheric one. The sign of the atmospheric mass splitting has not been determined yet. Therefore, the query is that whether the neutrino mass spectrum either does respect the normal ordering or does obey the inverted ordering. Moreover, the absolute neutrino mass scale is unknown. Theoretically, an important question is that how we can define this distinguished neutrino mixing pattern such that it would be perfectly feasible to obtain probable values of unknown parameters along with the other measured ones. The tribimaximal neutrino mixing matrix U T BM is one of the well-known neutrino mixing matrices [10][11][12], which is given by (1. 2) The general exact tribimaximal mixing matrix U T BM , regardless of the model, fixes the element (U e3 ) T BM = 0. It is important to note that, with the exception of θ 13 , the values of the mixing angles associated with the matrix U T BM are consistent with the data of Table I. The role of a non-zero θ 13 , or equivalently U e3 , is rather relevant to many concepts in the lepton sector. It is necessary for CP violation in neutrino oscillations and may be necessary to explain leptogenesis. For CP violation, of course, both θ 13 and the complex phase δ should be non-zero. Moreover, θ 13 = 0 mirrors an equivalent feature present in the quark sector, where mixing between all three generations is a confirmed result, although the mixing angles in the two sectors are very different. However, the discovery of the θ 13 , whose smallness (in comparison to the other mixing angles) proposes to modifythe the neutrino mixing matrix by means of a small perturbation about the basic tribimaximal structure. Consequently, it can lead to a realistic neutrino mixing matrix. There are many neutrino mass models [13][14][15][16][17][18] which can assist us to obtain the tribimaximal neutrino mixing matrix. In order to produce θ 13 = 0 starting from an initial tribimaximal structure, different approaches have been investigated: In [19,20], a perturbative analysis has been examined, from U T BM . In [21], an alternative method has been proposed in which a sequential 'integrating out' of heavy neutrino states is involved. The authors of [22][23][24] have employed the approach of parametrizing the deviation from the tribimaximal form. In [25][26][27][28][29], the deviations from tribimaximal mixing due to charged lepton effects and Renormalization Group running have been the direction of study. Alternative explorations, based on the A(4) symmetry, have been carried out in [30][31][32][33][34], while in [35][36][37][38][39], other discrete symmetries have been the basis for investigations. The mixing parameters and the mass ordering in Table I are required inputs for recognizing viable models associated to the neutrino masses. A natural option could be to take the mixing angles as either θ 13 = 0 or θ 23 = π 4 , such as those taken to obtain U T BM in Eq. (1.2), and the solar mass splitting is missing. In this regard, one proposal is provided as: which the atmospheric mass splitting and the maximal mixing (of this sector) arise from a unperturbed mass matrix while the smalle solar mass splitting and realistic U 13 , (θ 13 and δ), are generating by a perturbation. Moreover, it applies minor amendments to θ 12 . By employing different methods in widely contexts, a lot of endeavors have been pursued to generate some of the neutrino parameters in perturbation theory [40].
The purpose of our work is to introduce a framework, that constitutes a modification of the neutrino mass setup proposed by Friedberg and Lee (FL) [41]. The FL setting 1 can be regarded as a successful phenomenological neutrino mass model with flavor symmetry, which can be appropriately and equivalently employed for both Dirac and Majorana neutrinos. 2 However, in our work, we will employ instead a fundamental approach, which is different from that used in [1]. Let us be more precise. In our herein paper, the perturbation mass matrix will not be added by hand, but in contrast, it will be thoroughly computed within a series of steps. To this aim, we will be using the third perturbed mass eigenstate within perturbation theory methods. More concretely, the perturbation mass matrix will be thus constructed. We proceed as follows: (i) By employing perturbation theory in the mass basis with real parameters, we obtain the elements of the perturbation matrix which breaks the µ − τ symmetry. It will be seen that we get θ 13 = 0, but we do not have CP violation yet. (ii) We extend our work to the case with CP violation, and show that a complex perturbation matrix will be generated. In this case we have nonzero values for both θ 13 and δ. We also investigate the solar neutrino mass splitting in which an imaginary term will be dominant and lead to the generation of the Majorana phases. Consequently, we obtain CP violation along with a realistic neutrino mixing matrix. (iii) Finally, by comparing our phenomenological results with the corresponding experimental data, we will set up allowed parameter ranges, along with neutrino masses and CP violation phases.
This paper is organized as follows. In the next section we briefly introduce our (modified FL) model, and then present the results of the real and complex perturbation analysis described above in two subsections, separately.
Moreover, for the complex case we compute the perturbation mass matrix generating CP violation and we get a realistic neutrino mixing matrix. In section III, we map two of the experimental data onto the allowed region of our parameter space. Thereafter, we find the presently allowed ranges for all the parameters (especially perturbation parameters) of the model. Finally, not only do we check the consistency of all of the results with the available experimental data, but we also present our predictions for the actual masses and CP violation parameters. In section IV we summarize and analyze the results. In Appendix A, we briefly introduce the FL model.

II. MODIFIED FRIEDBERG-LEE MODEL
In this section we construct our model, within the FL framework, based on the basic tribimaximal neutrino mixing matrix. We compute the minimal neutrino mass perturbation matrix, along with a realistic neutrino mixing matrix. It is important to note that the innovative distinguishing characteristic of the preset work is that a neutrino mass perturbation will be erected from the minimal principle of the perturbation theory. Namely, we will not add a neutrino mass perturbation from any ad hoc assumptions, whilst in [1], the neutrino mass perturbation was added by hand by considering a few symmetries. It is important to note that, according to the experimental data reported in table I, the results of our herein improved model indicate an efficiency fitting increase of about 30 percent in contrast to those presented in [1]. In section III, we will further elaborate with more detail concerning this result. The tribimaximal neutrino mixing is a natural consequence of the M F L mass matrix in the case of µ−τ permutation symmetry (namely, the neutrino mass matrix remains invariant under interchanging indices µ and τ [42]). From Eq. (A2) it is apparent that M F L possesses exact µ − τ symmetry only when b = c. The magic property 3 of M F L obviously remains under exact µ − τ symmetry. Setting b = c and using the hermiticity of M F L , a straightforward diagonalization procedure Also, we should notice that in the pure FL model one of the neutrino masses is exactly zero. Moreover, just as in the general FL setting m 0 must be positive [41]. The reported experimental results have shown that the solar neutrino mass difference is tiny and the inequality ∆m 2 21 > 0 is confirmed [6]. Considering the experimental data, Eq. (2.1) yields b < 0, and |b| m 0 . Therefore, we can set λ = −b in Eq. (A3) and employ the transformations a → α = a − b and b → 0. Consequently, in the flavor basis, the unperturbed neutrino mass matrix is a magic µ − τ symmetry 4 and given by which has only two parameters α and m 0 . Of course, the diagonalized matrix which obtained from the mass matrix M 0 ν , Eq. (2.2), is a special case of the mass matrix Eq. (2.1), which has three parameters a, b and m 0 . It is of interest that they shares the same neutrino mixing matrix given by Eq. (1.2) 5 , U T BM [41].
The mass spectrum of M Here m 3 is unknown. In the next section, by comparing the results of our model with the experimental data, we will see that the sign of m 2 ) will be specified. We should mention that, up to now, the shortcomings are: (i) the absence of the solar neutrino mass splitting, (ii) the ordering of neutrino masses is unknown and (iii) the mixing matrix is still U T BM . Thus, the main objective will be to obtain the solar mass splitting by means of a mass perturbation, which is the cause of θ 13 = 0 and CP violation. Moreover, CP violation conditions necessarily mandate that µ − τ symmetry should be broken. An interesting question is: after the µ − τ symmetry breaking, will θ 23 = 45 • remain valid or not?
In summary, up to now, we have proposed that the modified FL neutrino mass matrix in Eq. (2.2) has a combination in which ∆m 2 12 and θ 13 are vanishing, while θ 23 = π 4 . Moreover, the atmospheric mass splitting ∆m 2 31 does not vanish. Furthermore, the solar mixing angle θ 12 can be selected as chosen by the popular mixing matrix as U T BM . This is a good estimate of the observed data although small characteristic are missing here. Therefore, the neutrino mass matrix in Eq. (2.2) has two mass eigenvalues, m 2 , which are degenerate, hence it is highly distinctive from the original FL model, in which all neutrino masses are different [41]. Moreover, to the best of our knowledge, such kind of FL modification has not been performed in the literature yet.
In the next stage, we will consider the attendance of a small contribution, which can be obtained by employing the perturbation theory, which generates small parameters in the neutrino mixing component, namely, U 13 , (θ 13 and δ), ∆m 2 21 and provides minor amendments to θ 12 (but not to θ 23 ). CP violation will be investigated. As previously mentioned, because of small θ 13 and ∆m 2 21 , we believe that the perturbative treatment is a more precise method than others for getting the correction of the U T BM . In our point of view, our work is special even in the perturbative treatment, because our minimal perturbation matrix is constructed solely from computing the third mass eigenstate, using the rules of perturbation theory, see Subsection A and B. Therefore, this perturbation matrix could induce both U 13 and ∆m 2 21 . It is worth mentioning, in the original FL model, by inserting phase factors in the neutrino mass matrix, the CP violation incorporate [41].
In order to establish the structure of the neutrino mass matrix, as noted before, our strategy is to employ the perturbation theory. Thus, we set M ν = M 0 ν + M ν where M ν << M 0 ν . In general, M 0 ν and M ν are symmetric and complex. However, as seen from Eq. (2.2), in this case M 0 ν is symmetric and real, i.e. it is Hermitian. In the following two subsections we will consider first the case where M ν is real and then the case where it is complex, respectively. In either situation we have θ 13 = 0, but CP is conserved when M ν is real. Furthermore, in the complex case the solar neutrino masses are split and CP is violated.
In the mass basis the eigenstates of M 0 ν (the unperturbed mass eigenstates) are as follows:

4)
4 A magic µ − τ symmetry for the mass matrix is synonymous to a TBM mixing structure [41,43]. 5 M 0 ν mass matrix in Eq. (2.2) with µ − τ symmetry and magic symmetry can be diagonalized by U T BM . Because of U T BM is a effect of µ − τ and magic symmetries in the neutrino mass matrix. Therefore U T BM is in itself a mixture of those symmetries. [44] 6 Naturally, if α m 0 , the neutrino masses would approach the quasi-degenerate regime.
in which the first two mass eigenstates are degenerate. We choose M ν such that ν 1 and ν 2 are its nondegenerate eigenstates, namely, ν 2 . Then we take (M ν ) 33 = 0 and consequently need to consider only (M ν ) 13 and (M ν ) 23 . Therefore, in order to reproduce the correct solar mixing, the basis vectors ν 1 and ν 2 are chosen, while the physical basis is fixed by the perturbation. It is straightforward to show that by expressing the mass eigenstates given by Eq. (2.4) in terms of the flavor basis, we can get the columns of U T BM as given by Eq. (1.2). Consequently, in the flavor basis, the eigenstates are given by In this subsection, our aim is to determine the third perturbed mass eigenstate in the CP conserving case. When this eigenstate is expressed in the flavor basis, we must obtain the third column of the neutrino mixing matrix, given by Eq. (1.1). Thus we can compute the elements of the perturbation matrix which will also be employed in the next subsection. Again, it is necessary to mention that in this method, we do not pick up a perturbation mass matrix by hand; instead, we compute it using of the third perturbed mass eigenstate, which is an unique feature and distinguishable from that used in [1]. As stated previously, here we assume M ν , which is symmetric, to also be real, and therefore Hermitian. Hence, while it may generate a nonzero θ 13 , it necessarily yields δ = 0, and so leads to no CP violation. For the perturbation expansion we keep terms up to linear order in s 13 . To first order we have where, To linear order in s 13 , we obtain C 31 = − 2 3 s 13 and C 32 = 1 3 s 13 , where we have used maximality of the 2 − 3 mixing angle, (θ 23 = 45 • ). Therefore, in the mass basis, by using Eq.
In this subsection, let us proceed our discussion, now addressing CP violation. We now assume M ν to be a complex symmetric matrix, thus not Hermitian; then this is also true for the total mass matrix M ν = M 0 ν + M ν . This is accomplished by considering nonzero values for both sin θ 13 and δ. The columns of the mixing matrix U in Eq.
3 > and |ν 3 is reproduced, to first order, by substituting the expressions associated to C 3j from Eq. (2.10) into Eq. (2.6). Consequently, by using an appropriate variant of Eq. (2.8) for this case, we get C 31 = − 2 3 s 13 e −iδ and C 32 = 1 3 s 13 e −iδ . It is important to note that, in the mass basis, due to the symmetric nature of M ν , it is easy to relate C 31 and C 32 to the elements of M ν as Up to now, by using Eq. (2.11) and C 31 , we have focused on deriving θ 13 = 0 via a perturbation analysis starting from an FL setting and the basic tribimaximal neutrino mixing matrix. Now, we investigate the solar neutrino mass splitting. In our framework of minimal perturbation we take (M ν ) 12 = (M ν ) 21 = 0. The first order corrections to the neutrino masses are obtained from m (1) We consider these first-order mass corrections as 2 , to sin θ 13 . In the next section, we will employ this parameter to obtain the order of sin θ 13 . In general, the solar mass splitting can take complex values, so let us mention that the Majorana mass is given by m (2.17) Therefore, in the Majorana case, φ is the origin of the Majorana phases which arise from the perturbation. In the next section, we obtain interesting results associated to ϕ and φ, namely, that m (1) 2 is dominated by its imaginary part, and so φ can take large values.
In order to relate our perturbation mass matrix to the FL model, let us rewrite M ν , [which is given by Eq. (2.15) and it was calculated in the mass basis, see Eq. (2.4)] in the flavor basis. Therefore, employing relation between mass and flavor basis and rewrite M ν , as We observe that the first and the second terms on the right-hand side are responsible for θ 13 and for ∆m 2 21 , respectively. Notice that M (f ) ν in Eq. (2.18) violates both µ − τ symmetry and the magic feature of the total mass matrix. Using degenerate perturbation theory [52], to linear order in s 13 , and the relation associated to C 3j given by Eq. (2.10), and aware of C ij = 0 if i = 3, which we know from the elements of M The nonzero δ indicates CP violation in the lepton sector.
In [53], with a different motivation in view, this same form for U has been discussed and the consistency with the observed mixing angles noted.
A rephasing-invariant measure of CP violation in neutrino oscillation is the universal parameter J [47] [given in Eq. (A4)], where its form is independent of the choice of the Dirac or Majorana neutrinos.
Using Eqs. (2.12) and (2.19), the expression for J can be simplified as (2.20) We should notice that in order to have CP violation in the lepton sector, both s 13 and η must take nonzero values.

III. COMPARISON WITH EXPERIMENTAL DATA
In this section, we compare the results obtained through our (modified FL) model with the experimental data [6? ]. It is important to note that in the original FL model [41] almost there is no numerical prediction for neutrino parameters. Therefore, in this section, we compare our herein results with the corresponding ones obtained in [1]. Let us perform this in two steps.
In the first step, we obtain the allowed ranges for the parameters of the neutrino mass matrix along with the perturbation term. We do this by mapping neutrino mass constraints obtained from the experimental data for ∆m where we have denoted these two solutions with upper (plus) and lower (minus) signs. As we see, we always have α + > 0 and α − < 0. In order to visualize the results of Eq. (3.1), in figure 1 we have plotted α versus m 0 . We see that we can have m 3 > 0 or m 3 < 0, corresponding to α + and α − , respectively. The range of m 0 is very impressible and important because when m 0 > 0.197 then we get α + → 0 which yields m + 3 → m 0 and α − → −m 0 which yields m − 3 → −m 0 . Such that for both of the cases, the results are unacceptable by comparing with the experimental data. The interesting point is that these results show that the physical mass spectrum is identical for both cases, the green and magenta curves, corresponding to m + 3 and m − 3 respectively. These curves are symmetric about the m 0 axis which is implying phase choice for the m 3 s, as it is seen in Table II. Namely, the Majorana phases are different for each case. Hence |m + 3 | = |m − 3 |, and the value of m 3 s in both case is the same. Naturally, the physics in the perturbation matrix elements does not depend on the chosen solution. We can see no reason to prefer either specific solution. For both values of α, our model has normal hierarchy, the same result as in [1].
In figure 2 we have plotted the overlap of ∆m 2 21 , by using Eqs.
2 , the solar neutrino mass splitting term, is dominated by its imaginary part. Therefore, due to the allowed range of ϕ [according to (3.2)], the origin of the Majorana phases, φ [in Eq. (2.17)], can take large values. This seems to suggest that the Majorana nature of neutrinos can be responsible for a large value of CP violation in nature [54].
We expect that the different nonzero components of the perturbation matrix Eq. (2.15) are roughly of similar order.
We may then expect ε ≡ m (1) 2 g(η) |∆m 2 31 |s13 ∼ O(1), and could predict the order of sin θ 13 in ε. Therefore, by using the order of α and |m (1) 2 | from the previous stage, we obtain sin θ 13 ∼ O(10 −1 ). In the second step, we obtain the allowed ranges for δ and J, the Jarlskog parameter, as in Eq. (A4) and for this we must first determine the allowed range of η. For this we use the expression of g(η) in Eq. (2.13). Recalling |m 1 | − |m 3 | < g(η) < |m 1 | + |m 3 |, and so 0.049 < g(η) < 0.400, we obtain the allowed range of η as |η| 92.29 • . In the limit m 0 → 0, we get η → δ and g(η) → |m 3 |. In order to get the allowed values of δ, we substitute the allowed ranges of m 0 , α and η into the expression associated to δ in Eq. Not only we have obtained all the parameters of the model [according to (3.2)], but now we can also make predictions for the masses of the neutrinos as well as the phases, see equation (3.4). We emphasize that we have made predictions As mentioned, φ is the origin of the Majorana phases which is retrieved from the perturbation. We could dispense with the overall phase, exp(iφ), and therefore we would have two phases that appear in the mass eigenvalues shown in Eq. (3.4). As we noted, these phases are not the same for |m 3 | ≡ m 3∓ . Namely, for m 3− we have ρ = (− φ 2 ) and σ = ( π 2 − φ 2 ). Whereas, for m 3+ , ρ = σ = (− φ 2 ). For the Dirac neutrinos, these phases can be removed, and for the Majorana neutrinos, these phases remain as Majorana phases and contribute to CP violation.
As mentioned, in the Appendix A, in the FL setting m 0 = 0, while the magnitude of the lower bound in Eq. (3.4) is about m 0 = 0. These predictions are compatible with the neutrino mass predictions in [1]. However, it is important to note that the results in this work (such as ∆m 2 21 and ∆m 2 31 ), established by means of a different approach, fit the experimental data much better than those reported in [1] especially the best fit column in Table I. In order to compare our methodology and approach to that in [1] and employing the current experimental data, the ratio ≈ (2.99 − 3.02) × 10 −2 , than the corresponding result reported in [1]. Namely, the agreement rate using our improved framework regarding the best fit experimental data is about 30 percent better than in [1]. Moreover, in [41], the phenomenological calculations associated to the CP violation case have not been compared with the corresponding experimental data. In Table II, we have displayed all of the relevant experimental data presented at Table I along with the predictions of our model. As is shown in Table  II, we have predictions for some physical quantities for which there are not any experimental data.
In this model, the magnitude of degeneracy associated with the neutrino masses is defined by m3−m1 m3 . Hence, the limit m 1 → 0 (m 3 ) means 0% (100%) degeneracy among the neutrino masses [55]. We should note that by using the allowed ranges of m 1 and m 3 in Eq. (3.4) the magnitude of degeneracy of the neutrino masses in our model is ≈ (0% − 97%).
For the flavor eigenstates, we can just calculate the expectation values associated to the masses. Therefore, we can    sets the following constraint m ν ββ < (0.061 − 0.165) eV at 90% CL [57]. Our prediction (better than those in [1]) for this quantity is m ν ββ ≈ (0.0028 − 0.12) eV which is consistent with the result of kamLAND-Zen experiment.
One of the main experimental result is the sum of the three light neutrino masses which has just been reported by the Planck measurements of the cosmic microwave background (CMB) at 95% CL [58] as m ν < 0.12eV (Plank+WMAP+CMB+BAO). (3.6) In our model, we obtain m ν ≈ (0.058 − 0.597)eV , which is in agreement with (3.6).

IV. DISCUSSION AND CONCLUSIONS
In the next stage, we will consider the attendance of a small contribution, which can be obtained by employing the perturbation theory, which generates small parameters in the neutrino mixing component, namely, U 13 , (θ 13 and δ), ∆m 2 21 and provides minor amendments to θ 12 (but not to θ 23 ). CP violation will be investigated. We should emphasize that in our present work, the method for retrieving the minimal perturbation mass matrix is completely different from those present in the literature (see, e.g., [1] and references therein) and it can be considered as a more fundamentally based approach. The distinguishing features of our herein model are: i) solely from using the third perturbed mass eigenstate and by employing the rules of the perturbation theory, we constructed the minimal perturbation matrix of the basic tribimaximal mixing matrix, producing a modified Friedberg-Lee model. Therefor, it was produced from the rules of perturbation theory. 9 ii) the perturbation mass matrix is simultaneously responsible for the solar neutrino mass splitting and CP nonconservation in the lepton sector. Consequently, due to these two initiatives and distinctive consideration, our modified framework regarding the best fit experimental data is 30 percent better than in [1]. The model is based on the tribimaximal mixing matrix in which the experimental data of mixing angles (except θ 13 ) is well approximated. Therefore, by employing the Friedberg-Lee neutrino mass framework, we obtained the tribimaximal structure which led us to produce a mass matrix constrained by the elements of the TBM mixing matrix and the experimental data. The mass matrix thus obtained (unperturbed mass matrix) loses the solar neutrino mass splitting whilst it remains as a magic and symmetric matrix under µ − τ symmetry. At this level, by employing perturbation theory, we generate a perturbation matrix which breaks softly both the µ − τ symmetry and the magic feature, and consequently causes CP violation.
Our investigation proceeded in two stages [of section II]: CP conservation and CP violation. In the first stage, we obtained the elements of the perturbation matrix in a non CP violation case. In this case, the elements of the perturbation matrix are real, and therefore δ = 0 while θ 13 = 0.
In the second stage, we extended our study to the case of CP violation, and we obtained the complex elements of the perturbation matrix in both the flavor and mass bases. Moreover, (a) We retrieved a realistic mixing matrix with δ = 0. However, in this case, the µ − τ symmetry is softly broken, but still we have θ 23 = 45 • . (b) We obtained the solar neutrino mass splitting dominated by an imaginary term. Therefore, the most important corresponding results or claims are: the possibility that Majorana neutrinos exist and are more abundant than Dirac neutrinos, and, that this feature would be the reason for the magnitude of CP violation verified in nature.
In order to get valuable predictions concerning neutrino masses, δ, the origin of the Majorana phases and J, we compared the results of our phenomenological model with experimental data. We have shown that how our phenomenological model whether or not is consistent with experimental data. Mapping two sets of experimental data, namely, the allowed ranges of ∆m 2 21 and ∆m 2 31 onto the allowed region of our parameter space can determine valid values for our parameters. The consistency of experimental data with the allowed ranges of our parameters shows that our model has normal hierarchy. We then predict the perturbation mass matrix and the values of three masses, m 1 ≈ (0 − 0.197)eV , |m 2 | ≈ (0.00862 − 0.19719)eV , and m 3 ≈ ∓(0.0490 − 0.2033)eV . Therefore, the magnitude of degeneracy for neutrino masses at the end of the allowed ranges is about 97%. We have shown that the order of sin θ 13 , which is estimated from the order of the mass parameters, is consistent with the experimental data. We also obtain predictions for the CP violation parameters δ, J, ρ and σ. These are |δ| Our predictions for the neutrino masses and CP violation parameters could be tested in future experiments such as the upcoming long baseline neutrino oscillation ones. Moreover, our predictions are entirely consistent with the constraints reported by Planck, WP and high L measurements and the KamLAND-Zen experiment [56,58].
In our model, the minimal perturbation matrix was obtained by means of a fundamental process, and not merely added by hand. Therefore, we can claim that the model presented here can be regarded as presenting a more comprehensive scenario. Another important point is that although our predictions are in correspondence with those in [1], we should emphasize that our outcomes (such as ∆m 2 21 and ∆m 2 31 ), constitute with the best fit experimental data reported in table I. Last and not least to emphasize, the consistency fitting rate concerning the best fit experimental data is about 30 percent more efficient trough the framework introduced in this paper rather than in [1].
Finally, we plan to subsequently proceed and investigate a 6 × 6 neutrino mass matrix, by using the same methodology, to obtain the corresponding perturbation mass matrix and CP violation. In the FL model, the mass eigenstates of the three charged leptons are the same as their corresponding flavor eigenstates. Therefore, the neutrino mixing matrix is simply the 3 × 3 unitary matrix U , which transforms the