Exponential parameterization of neutrino mixing matrix with account of CP-violation data

The exponential parameterization of the Pontecorvo–Maki–Nakagawa–Sakata mixing matrix for neutrinos is discussed. The exponential form allows easy factorization and separate analysis of the CP-violating and Majorana terms. Based upon the recent experimental data on the neutrino mixing, the values for the exponential parameterization matrix for neutrinos are determined. The matrix entries for the pure rotational part in charge of the mixing without CP violation are derived. The complementarity hypothesis for quarks and neutrinos is demonstrated. A comparison of the results based on most recent and on old data is presented. The CP-violating parameter value is estimated, based on the so far imprecise experimental indications, regarding CP violation for neutrinos. The unitarity of the exponential parameterization and the CP-violating term transform is confirmed. The transform of the neutrino mass state vector by the exponential matrix with account of CP violation is shown.


Introduction
One of the paramount achievements of physics of the 20th century was certainly the formulation of the Standard Model [1][2][3], which unifies the description of electromagnetic and weak interactions in one theory. An important role in the Standard Model is played by neutrinos. In the framework of the Standard Model neutrinos may have three flavors, matching three charged leptons, with which they interact by means of weak interaction. The proper states form full normalized orthogonal basis. Originally the Standard Model assumed massless neutrinos; later it was adapted to incorporate their mass. The existence of the mass of the neutrinos means the existence of at least three massive neutrino states ν 1 , ν 2 , ν 3 , and, also, it means the existence of the neutrino oscillations [4], i.e., neutrino flavors change while they a e-mail: zhukovsk@physics.msu.ru b e-mail: francisco.melazzini@physics.msu.ru are propagating. Evidence for that was found in experiments for mixing of solar neutrinos [5], atmospheric neutrinos [6], and reactor neutrinos [7]. The phenomenon of neutrino mixing was predicted by Pontecorvo [8][9][10]. The 2015 Nobel Prize in physics recognized the discovery of neutrino oscillations. The transforms from the mass state basis to the flavor state basis and vice versa is performed with the help of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix. This transform can be considered similarly to the quark mixing by the CKM matrix. The familiar neutrino states, ν e , ν μ , ν τ , are linear combinations of neutrinos with different masses, ν 1 , ν 2 , ν 3 : where V P M N S is the unitary PMNS mixing matrix [11]. Note the analogy with the mixing of the low elements of left components of the quark spinors, Lepton mixing presumes that a charged W-boson can couple to any mass state of charged leptons (e, μ, τ ) with any mass state of neutrino. For example, V αi would be the amplitude of the boson W + decay in a pair of lepton type α and neutrino type i and, hence, the production of the lepton α and neutrino state α implies that all neutrino mass states participate in it in a superposition. In the following we do not consider a sterile neutrino, [12][13][14], which does not interact with Wand Z-bosons. Thus, we end up with a unitary 3 × 3 mixing matrix V, factorized by the matrix: c i j = cos θ i j , s i j = sin θ i j , i, j = 1, 2, 3, and P M jr = diag e iα 1 /2 , e iα 2 /2 , 1 describes the possible Majoranas nature of neutrinos by the phases α 1 and α 2 . If α 1,2 = 0, then the neutrinos are the Majorana particles, i.e. they are identical to their antiparticles; the phases α 1 , α 2 play a role in the processes, which do not preserve the lepton number. The role of the matrix V in the parameterization (3) is similar to the one that the CKM matrix plays in quark mixing [15,[17][18][19][20]. The PMNS matrix is fully determined by four parameters: three mixing angles θ 12 , θ 23 , θ 13 , and the phase δ, in charge of the CP violation description [15]. The experimental values of the mixing angles are relatively well determined [15,16,[21][22][23][24]: Contrary to quark mixing angles, these are not small and the expansion in series of the only small parameter is not possible. Thus, there is no small parameter, like λ = sin θ Cabibbo ≈ 0.22 [25], for neutrino mixing. Experimentally determined absolute values for the elements of the PMNS matrix read as follows [15]: Moreover, there are indications that the CP-violating phase may have a non-zero value; moreover, very approximately it is supposed to be as large as δ ≈ 300 • (see [26,27]).

Exponential mixing matrix
The exponential parameterization for the neutrino mixing matrix was outlined in [28] and then in [29]; it is constructed similarly to that for quarks [30,31]. Unitarity of the exponential mixing matrix V = exp A is guaranteed by the anti-Hermitian form of the exponent [33]), which depends on the mixing parameters λ i , and on the CP-violation phase δ CP . Note that for δ CP = 2π n we obtain simply a rotation matrix around the axis in space [31]. The matrix also becomes real for δ CP = π(2n + 1). The most important advantage of the exponential parameterization for the mixing matrix with respect to the commonly known standard parameterization is that the exponential parameterization allows easy separation of the contributions of the rotation part, the CP violation and possible other terms in stand alone factors. This separation can be made in a variety of modes, of which we omit the details here; proper discussion was presented, for example, in [31,32]. The following unitary parameterization was proposed in [29]: where the rotation part is given by the real exponential matrix the CP violation is accounted for by and it contains an imaginary component, and the Majorana part, depends on the Majorana phases in the exponential: The details of the splitting between CP-conserving and CPviolating terms in the above parameterization can be found in [29] (also compare with [31,32]). The values of the Majorana phases α 1 and α 2 are at present undetermined; the value of the CP-violating phase δ CP can be figured from the existing experimental indications and estimations (see, for example, [26,27]). The rotation matrix can be conveniently presented in the form of the rotation in the angle around the axis, given by the vector n = (n x , n y , n z ): In this form it represents a three-dimensional rotation generator. Then (15) links the entries of the rotational matrix (10) in the exponential parameterization (9) with the rotation angle around the axis n = (n x , n y , n z ) with the following coordinates: Thus, Eqs. (17)- (16) relate the elements of the exponential parameterization with the single axis-rotation angle matrix components. Note that the angles of the rotation in the standard parameterization matrix, (3) c i j = cos θ i j and s i j = sin θ i j , are different from those in our parameterization (9).
Omitting the Majorana part, the exponential parameterization (9) reads as follows: The values of the entries of the rotation matrix M i j can be derived from the following tensor identity (see [34]): where we denote δ i j the Kronecker symbol, ε i jk is the Levi-Civita symbol, n i are the components of the vector n = (n x , n y , n z ) and is the rotation angle. The expressions, relating the entries of the standard parameterization c i j and s i j with n and , can be derived from (18), (19), (20), and (3), but they are very cumbersome and we omit them for brevity. From the experimental data [26] by using matrix equations and expansions we obtain for the rotation vector in 3D space the following coordinates: and for the rotation angle around this axis we obtain The precision of the above values is determined by the errors in the experimental data evaluation and is about of 5 %. From (15)-(17) we obtain the following values for the entries of the rotation matrix (15) in the exponential parameterization: This result is in agreement with the experimental data and with the respective standard parameterization in the PMNS matrix (2). Moreover, based upon the values of the mixing angles for quarks θ Q12 = 13.14 • , θ Q23 = 2.43 • , θ Q13 = 0.23 • , we determine the direction of the rotation vector in space (17) for the exponential parameterization for quarks as follows: n quark = (0.1829,0.0206,0.9831).
Now, upon the comparison with the above determined coordinates of the rotation vector for neutrinos we note that n quark and n neutrino constitute the angle of ≈44 • , as demonstrated in Fig. 1. This fact is interesting itself and it is the demonstration of the so-called hypothesis of complementarity for neutrinos and quarks [35,36], according to which the rotation axes for quarks and neutrinos form the 45 • angle; however, this last statement is rather an observation since it does not have solid theoretical foundation and there are no physical reasons. Note that the obtained value of 44 • differs from 45 • by ≈2 %, which is within the margin of errors of the original experimental data sets, which determines the entries of the exponential mixing matrix and the rotation vectors directions.

Exponential parameterization and CP violation
Now let us take advantage of the possibility, given by the exponential parameterization, which allows us to factorize separately the contributions of the rotation, the CP-violation, and the Majorana term. We can write the matrix product P CP P Mjr in the following form, which reminds one of the rotation by the angle 2 with the proper weights for each entry: Omitting the Majorana part, the complexity due to the CP violation vanishes if either = ±π n/2, or if ω ± ∈ Reals. The value of is limited by μ, which is small, and the above expression (26) becomes real in the case of = 0, which, moreover, returns the unity matrix P CP = I for δ CP = ±2π n.
The absolute values of the entries of the CP-violating matrix are at their maximum for δ CP = π ± 2π n, but then the whole P CP matrix becomes real and the CP violation in fact vanishes. The major contribution of the complex term due to the CP violation is achieved for δ CP = π/2 ± π n. All of the above is in complete analogy with the behavior of the CP-violating phase δ in the standard parameterization. Based upon the present knowledge of the CP violation for neutrinos and on experimental indications and estimations [26,27] and accounting for the above discussion, the CPviolating phase in the exponential parameterization has the value δ CP ∼ = −60 • identical to δ ≈ 300 • . Then we obtain Note that the accuracy of the existing indications with regard to the CP violation is rather low (+60 − 120 %) and thus the above values are also approximate. The parameter is not very sensitive to the exact value of the angle δ CP and its absolute value varies from 0 to the maximum of |μ| = 0.342 for δ CP = π ± 2π n. Note that in any case the value of is small and cos ∼ = 1. For δ CP = −60 • we obtain = μ sin δ CP 2 = μ sin(−30 • ) = −μ/2 = 0.171, cos 2 = cos μ ∼ = 0.942, sin 2 = − sin μ ∼ = 0.335. (28) Then with account of the Majorana phases and of δ CP ∼ = −60 • , we obtain the following matrix: The Majorana phases interplay with the CP-phase only in the entry (3,1) in the factor e i α 1 2 ω − . Otherwise, the Majorana phases just bring about more complexity in the result. For non-Majorana, but Dirac particles, the matrix obtained above, (29), represents just a slight deviation from the unitary matrix I: It is now evident that for non-Majorana particles the CPviolating term (29) is the mixing matrix for two lepton generations: electron, tauon, and proper neutrinos. Employing the generating functions for the Bessel functions cos (x sin α) = ∞ n=−∞ J n (x) cos nα and sin (x sin α) = ∞ n=−∞ J n (x) sin nα, we easily derive the following expression for the general form of the matrix V MCP , in which the contributions of the CP-violating phase δ C P and of the rotation matrix parameter μ are separated: The above result (32) is valid for arbitrary values of δ CP , α 1 , α 2 . For α 1 = 0, (32) reduces to the symmetric form.
With account of the obtained values of the entries of the exponential parameterization of the CP-violating term and for non-zero Majorana phases α 1 , α 2 , the vector of the mixed neutrino state is transformed as follows: We underline that the above transform by V MCP , as well as the transform by the purely rotational part P Rot , is unitary. It can be verified directly with the help of the Hermiteconjugated matrix: This ensures the unitarity of the whole exponential parameterization (9) of the PMNS mixing matrix:

Conclusions
The exponential parameterization of the mixing matrix for neutrinos is explored with account of the present experimental data. The proper entries of the exponential mixing matrix are determined; the CP-violating term in the exponential parameterization is estimated. Based upon the accuracy of the experimental data, the range of the values for the parameters of the neutrino mixing matrix is given. Without CP violation the neutrino mixing represents in fact the geometric rotation in three-dimensional space. In this simple case mixing can be viewed as the rotation by the angle around the axis in threedimensional space. This interpretation follows straightforwardly from the structure of the exponential mixing matrix. Evidently, there is no mixing for = 0, when the mixing matrix without CP violation reduces to the unit matrix, I. Based upon the recent data, we have obtained the value for the rotation angle ∼ = 49.8 • and the coordinates of the rotation vector n = (0.702, 0.394, 0.593). This value of the rotation angle is somewhat smaller than that based upon the tribimaximal parameterization: TBM ∼ = 56.6 • . Moreover, the direction of the rotation vector differs from that of the vector for the TBM matrix: n TBM = (0.7858, 0.2235, 0.5777). The difference in their directions in 3D space is 11 • . Interestingly, the angle between the axes of rotation for quarks and neutrinos remains unchanged and equals ≈45 • , despite the change of 11 • in the direction of the neutrinos rotation axis, verified in the last 10 years. This demonstrates the hypothesis of complementarity for quarks and neutrinos [35,36].
The exponential parameterization allows factorization of the CP and the Majorana contributions and is evidence that the CP-term can also be viewed as a sort of rotation with different weights for the matrix entries. We have calculated the entries of the rotational mixing matrix λ ∼ = 0.516, μ ∼ = −0.342, ν ∼ = −0.611 (see (23)) and we have estimated the entries of the CP-violating matrix in the exponential parameterization (see (29), (31)) from the current indications on CP violation: δ CP ∼-60 • . This value is quite approximate due to uncertain experimental data, regarding CP violation for neutrinos. We calculated the CP-violating exponential matrix for the extremities of the range of δ CP from 0 • to 180 • . The rotation angle 2 for the CP-violating matrix is rather small: = μ/2. By means of the exponential parameterization one can easily transform the neutrino state vector, distinguishing the CP-violating terms for each type of neutrino. In the case of Dirac neutrinos and if δ CP = ±2π n, we get pure rotation, since P CP = I. If δ CP = π ± 2π n, then we end with the real P CP matrix for the CP-violating term, which means the absence of CP violation in this case. For Majorana neutrinos the mass state vector ν 1 , ν 2 , ν 3 is transformed with complex weights, as demonstrated in (33).
We addressed an exponential presentation of the mixing matrix and obtained with its help results and interpretations that can be useful for the treatment and analysis of new experimental data, regarding the neutrino oscillations in currently running experiments as well as in planned experimental projects.