Can One have Significant Deviations from Leptonic $3\times 3$ Unitarity in the Framework of Type I Seesaw Mechanism?

We address the question of deviations from $3\times 3$ unitarity of the leptonic mixing matrix showing that, contrary to conventional wisdom, one may have significant deviations from unitarity in the framework of type I seesaw mechanism. In order for this scenario to be feasible, at least one of the heavy neutrinos must have a mass at the TeV scale, while the other two may have much larger masses. We present specific examples where deviations from $3\times 3$ unitarity are sufficiently small to conform to all the present stringent experimental bounds but are sufficiently large to have the potential for being detectable at the next round of experiments.


Introduction
The discovery of neutrino oscillations and at least two non-vanishing neutrino masses, provides clear evidence for Physics Beyond the Standard Model (SM). The simplest extension of the SM accommodating two non-vanishing neutrino masses involves the addition of at least two right-handed neutrinos. The most general gauge invariant Lagrangian includes a right-handed bare Majorana mass matrix M. As a result, the scale of M can be much larger than the electroweak scale, which leads to an elegant explanation for the smallness of neutrino masses, through the seesaw mechanism [1], [2], [3], [4], [5]. The seesaw mechanism necessarily implies violations from 3 × 3 unitarity of the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix, as well as Z-mediated lepton flavour violating couplings. The introduction of these heavy right-handed neutrinos can also have profound cosmological implications since they are a crucial component of the Leptogenesis mechanism to create the observed Baryon Asymmetry of the Universe (BAU) [6]. Leptogenesis is a very appealing scenario [7] but for heavy neutrinos with masses many orders of magnitude higher than the electroweak scale it is difficult or impossible to test it at low energies [8], [9]. Furthermore, without a flavour model, CP violation at high energies, relevant for Leptogenesis, cannot be related to CP violation at low energies [10], [11]. This may be possible in the context of a flavour model, involving symmetries which allow to establish a connection between high and low energies [12], [13], [14], [15], [16], [17], [18], [19].
At this stage, it should be emphasised that within the seesaw type I framework, the observed pattern of neutrino masses and mixing does not require that all the heavy neutrino masses be much larger than the electroweak scale. In this paper, we carefully examine the question of whether it is possible, within the seesaw type I mechanism, to have experimentally detectable violations of 3 × 3 unitarity, taking into account the present experimental constraints. In particular, we address the following questions: i) In the seesaw type I mechanism, is it possible to have significant deviations from 3 × 3 unitarity of the leptonic mixing matrix? By significant, we mean deviations which are sufficiently small to conform to all present stringent experimental constraints on these deviations, but are sufficiently large to be detectable in the next round of experiments. These experimental constraints arise from bounds on rare processes.
ii) In the case the scenario described in (i) can indeed be realised within the framework of seesaw type I, what are the requirements on the pattern of heavy neutrino masses?
For definiteness, we will work in a framework where three right-handed neutrinos are added to the spectrum of the SM. Our analysis starts with the introduction of the unitary 6 × 6 mixing matrix V , characterising all the leptonic mixing. We write this 6 × 6 mixing matrix in terms of four blocks of 3 × 3 matrices. Using unitarity of V , we show that the full matrix V can be expressed in terms of only three blocks of 3 × 3 matrices. Then we apply these results to the diagonalisation of the 6 × 6 neutrino mass matrix, including both Dirac and Majorana mass terms. Through the use of a specially convenient exact parametrisation of the 6 × 6 leptonic unitary mixing matrix, we evaluate deviations from 3 × 3 unitarity and derive the maximum value of the lightest heavy neutrino mass which is required in order to generate significant deviations from unitarity in the framework of the seesaw type I mechanism.
The paper is organised as follows. In the next section, we review the seesaw mechanism, define our notation and introduce a specially convenient exact parametrisation of the 6 × 6 leptonic mixing matrix V . We evaluate the size of the deviations of 3 × 3 unitarity in the present framework and derive a constraint on the magnitude of the mass of the heavy Majorana neutrinos, in order to have significant deviations of unitarity. Numerical examples are given in section 3 and some of the derivations of our results are put in an Appendix. Finally we present our conclusions in the last section.
2 Deviations from Unitarity in the Leptonic Sector

Type I Seesaw mechanism
In the context of the Type I seesaw mechanism, with only three right-handed neutrinos added to the Lagrangian of the SM, the leptonic mass terms are given by: There is no loss of generality in choosing a weak basis where m l is already real and diagonal. The analysis that follows is performed in this basis. The neutrino mass matrix M is a 6 × 6 matrix and has the form: This matrix is diagonalised by the unitary transformation where K, R, S and Z are 3 × 3 matrices. For K and Z non singular, we may write From the unitary relation V V † = 1I (6×6) , we promptly conclude that The matrix V can thus be written: We have thus made clear that the unitary 6 × 6 matrix V can be expressed in terms of three independent 3 × 3 matrices. From the unitarity of V , we obtain: showing that the matrix X parametrizes the deviations from unitary of the matrices K and Z. More explicitly: From Eq. (3) we derive: replacing Z † m T from Eq. (13) into Eq. (11) we get which implies that: where O c is a complex orthogonal matrix, i.e., O T c O c = 1I, or explicitly: It should be stressed that the parametrisation of the 6 × 6 unitary matrix V given by Eq. (8) has the especial property of allowing to connect in a straightforward and simple way the masses of the light and the heavy neutrinos through an orthogonal complex matrix O c , as can be seen from Eqs. (15) and (16). This is an important new result which plays a crucial rôle in our analysis.
Since O c is an orthogonal complex matrix, not all of its elements need to be small; furthermore, not all the M i need to be much larger than the electroweak scale, in order for the seesaw mechanism to lead to naturally suppressed neutrino masses. These observations about the size of the elements of X are specially relevant in view of the fact that some of the important physical implications of the seesaw model depend crucially on X. In particular, the deviations from 3 × 3 unitarity are controlled by X, as shown in Eq. (10).
Given the importance of the matrix X, one may ask whether it is possible to write the 6 × 6 unitary matrix V in terms of 3 × 3 blocks, where only 3 × 3 unitary matrices enter, together with the matrix X. In the Appendix, we show that this is indeed possible, and that the matrix V can be written: where Ω and Σ are 3 × 3 unitary matrices given by: and U, W are the unitary matrices that diagonalise respectively X † X and XX † : It is also shown in the Appendix, that U K and W Z defined by: are in fact unitary matrices.
As will be explained in the next section, this will allow us, in our analysis, to trade the matrix K by the combination U K U † which we identify as the best fit for U PMNS derived under the assumption of unitarity, multiplied by the remaining factor that parametrises the deviations from unitarity.

On the size of deviations from unitarity
In the framework of the type I seesaw, it is the block K of the matrix V that takes the rôle played by U PMNS matrix at low energies in models with only Dirac-type neutrino masses. Clearly, in this framework, K is no longer a unitary matrix. However, present neutrino experiments are putting stringent constraints on the deviations from unitarity. In our search for significant deviations from unitarity of K, we must make appropriate choices for the matrix X in order to comply with the experimental bounds, while at the same time obtain deviations that are sizeable enough to be detected experimentally in the near future. It is our aim to show that, contrary to common wisdom, we can achieve this result with at least one of the heavy neutrinos with a mass at the TeV scale, without requiring unnaturally small Yukawa couplings and still have light neutrino masses not exceeding one eV.
Deviations from unitarity [20], [21], [22], [23], [24] of K have been parametrised as the product of an Hermitian matrix by a unitary matrix [23]: where η is an Hermitian matrix with small entries. In order to identify the different components of our matrix K, given in Eq. (18), with the parametrisation of Eq. (22) we rewrite K as: inside the square brackets we wrote the Hermitian matrix that we identify with (1I − η), and which will parametrise the deviations from unitarity. The matrix V ≡ U K U † is a unitarity matrix which is identified with U PMNS obtained from the standard parametrisation [25] for a unitary matrix. One can also write: (20). Identifying the second expression of Eq. (24) to (1I − η) we derive: for small d 2 X . The matrix U PMNS is then fixed making use of the present best fit values obtained from a global analysis based on the assumption of unitarity. As pointed out in [23], from the phenomenological point of view it is very useful to parametrise K with the unitary matrix on the right, due to the fact that experimentally it is not possible to determine which physical light neutrino is produced, and therefore, one must sum over the neutrino indices. As a result, most observables depend on KK † which depends on the following combination: The standard parametrisation for U PMNS is given by [25]: with P given by where c i j = cos θ i j , s i j = sin θ i j and δ is a Dirac-type CP violating phase, while α 21 , α 31 denote Majorana phases. Neutrino oscillation experiments are not sensitive to these factorisable phases.
The specific bounds vary slightly from group to group. For definiteness we present in Table  1 the present bounds on neutrino masses and leptonic mixing from [26]. The quantities ∆m 2 i j are defined by (m 2 i − m 2 j ).
Ref. [23] also compares these bounds with those of previous studies [22], [29], pointing out that in general there is good agreement. Variations in the scale of the masses of the heavy neutrinos lead to small effects and therefore, do not significantly change our analysis.

The elements of the neutrino Dirac mass matrix m and deviations from unitarity
In this subsection, we show that there is a correlation among: • The size of deviations from unitarity of the 3 × 3 leptonic mixing matrix.
• The mass of the lightest heavy neutrino.
From Eq. (13) we get the experimental fact that K is almost unitary implies that Z is also almost unitary. Therefore the Dirac mass matrix m is of the same order as X times D. Notice that the scale of D may be of the order of the top quark mass, so that indeed the Yukawa couplings need not be extremely small.
The elements of the neutrino Dirac mass matrix m are connected to the deviations from unitarity of the 3 × 3 leptonic mixing matrix. From Eq. (30) together with Eqs. (21), (23) and (24), we obtain: where we have used d X = W † X U from Eq. (20). Thus, we find (32) As previously emphasised, deviations from 3 × 3 unitarity in U PMNS are controlled by the matrix X, as it is clear from Eq. (10). For X = 0, there are no deviations from unitarity. Small deviations from unitarity correspond to d X small and, in that case, one has to a very good approximation, which can be written as: It can be shown that, using the properties of orthogonal complex matrices, only one of the d X i , corresponding to d X 3 , can have a significant value (e.g. d X 3 ≈ 10 −2 ), while the other two are negligible. Thus, we find in good approximation or using the unitary of W From Eq. (37), it is clear that for significant values of d X 3 , M 1 cannot be too large in order to avoid a too large value of Tr mm † , which in turn would imply that at least one of the m i j 2 is too large. This can be seen in both Fig. 1 and Fig. 3 where we plot 1 2 d 2 X 3 versus M 1 . Significant values of d 2 X 3 can only be obtained for M 1 ≤ 1 − 2 TeV. In Figs. 2, 4, we also plot the |η 11 | deviations from unitarity.

Numerical Examples
In this section, we present some illustrative results of our numerical analysis showing that it is possible to obtain the observed pattern of neutrino masses and mixing without requiring that all masses of the heavy Majorana neutrinos, M i , be much larger than the electroweak scale. We have realistic examples even when all the three heavy neutrinos have masses below 2 TeV. One might expect that lowering the scale of the heavy neutrino masses would result in the need for extremely small Yukawa couplings for the Dirac mass terms, thus defeating the rationale for the seesaw mechanism. However, this is not the case and, to illustrate, we include for each example the corresponding moduli of the entries of the neutrino Dirac mass matrices m and the trace of the product mm † .
In Table 2, we include examples with normal ordering of light neutrino masses, for a In all examples given in Table 2 the orthogonal complex matrix is of the form: with different choices of x and τ for the different cases, respectively: for this choice, we made the simplest one, by fixing W Z to be equal to the identity. With a different choice of W Z we could in principle homogenise the orders of magnitude of the entries of m so that all of the Yukawa couplings would be of the same order of magnitude or close.
In Table 3 we include examples with inverted ordering of light neutrino masses, for a particular fixed value of these three masses, common to all examples. We consider the same three different hierarchies for the heavy Majorana neutrino masses M i , as in the cases of The choices of x and τ are respectively: In all our examples the matrix X will have several entries of order at most 10 −2 .

Conclusions
We have studied the possibility of having significant deviations from 3 × 3 unitarity of the leptonic mixing matrix in the framework of type-I seesaw mechanism. The analysis was done in the framework of an extension of the Standard Model where three righthanded neutrinos are added to the spectrum. We have shown that the 6 × 6 unitary leptonic mixing matrix V can be written in terms of two unitary 3 × 3 matrices and a matrix denoted X, which controls the deviations from unitarity of the 3 × 3 PMNS matrix.
This parametrisation of the matrix V , played a crucial role in showing that one may have significant deviations from 3 × 3 unitarity while conforming to all present data on neutrino masses and mixing, as well as respecting all stringent bounds on deviations from 3 × 3 unitarity of the PMNS matrix.
We have presented specific examples where the above deviations from unitarity are sufficiently large to have the potential for being observed at the next round of experiments. An important feature of our analysis is the fact that the mass of the lightest heavy sterile neutrino is, in principle, within experimental reach. This can be achieved without unnaturally small neutrino Yukawa couplings, thus defying the conventional wisdom that    We therefore conclude that K U 1I + d 2 X = U K and Z W 1I + d 2 X = W Z are unitary matrices, and thus also that or using the diagonalisation of the Hermitian matrices 1I + X † X and 1I + X X † in Eq.(42), we write this, as which lead to Eq. (18).