Ific/1745 QuasiDirac neutrino oscillations
Abstract
Dirac neutrino masses require two distinct neutral Weyl spinors per generation, with a special arrangement of masses and interactions with charged leptons. Once this arrangement is perturbed, lepton number is no longer conserved and neutrinos become Majorana particles. If these lepton number violating perturbations are small compared to the Dirac mass terms, neutrinos are quasiDirac particles. Alternatively, this scenario can be characterized by the existence of pairs of neutrinos with almost degenerate masses, and a lepton mixing matrix which has 12 angles and 12 phases. In this work we discuss the phenomenology of quasiDirac neutrino oscillations and derive limits on the relevant parameter space from various experiments. In one parameter perturbations of the Dirac limit, very stringent bounds can be derived on the mass splittings between the almost degenerate pairs of neutrinos. However, we also demonstrate that with suitable changes to the lepton mixing matrix, limits on such mass splittings are much weaker, or even completely absent. Finally, we consider the possibility that the mass splittings are too small to be measured and discuss bounds on the new, nonstandard lepton mixing angles from current experiments for this case.
Keywords: Neutrinos, quasiDirac, pseudoDirac, Majorana, experimental constraints.
AHEP Group, Instituto de Física Corpuscular, C.S.I.C./Universitat de València
Edifício de Institutos de Paterna, Apartado 22085, E–46071 València, Spain
1 Introduction
Neutrino oscillation experiments cannot distinguish Dirac from Majorana neutrinos, hence it is still unknown whether or not lepton number is conserved. Other processes, such as neutrinoless double beta decay [1, 2], need to be probed in order to answer this question. However, while the nature of neutrinos is often seen as a dichotomy, presenting two sharply distinct scenarios, the Dirac neutrino case can be seen as a limit of the more general Majorana case in which lepton number violating mass terms are zero, and this limit can be approached smoothly.
In practice, one can start with a model of Majorana neutrinos and get a phenomenology arbitrarily close to the one of a model of Dirac neutrinos. This can already be seen with only one generation of active () and sterile neutrinos (). In the basis the most general mass matrix reads:
(1) 
If , lepton number is preserved and neutrinos are Dirac
particles. This limit can alternatively be characterized by two exactly degenerate mass eigenstates composed in equal parts of
and : and
.
Let us rewrite eq. (1) using:
(2)  
(3) 
As long as and are much smaller than one, we obtain:
(4)  
(5)  
(6) 
Departures from the Dirac case therefore can manifest themselves as either new mass splittings or new mixing angles (or, in general, both). Moreover, as this simple example shows, mass splittings and mixing angles are completely independent of each other. Note that for small values of and , lepton number violation is naturally suppressed, as expected. This can be most easily seen in our one generation scenario for the double beta decay observable : for () it is straightforwardly calculated to be ().
We have therefore the following situation. Oscillation experiments cannot distinguish a model with Majorana neutrinos (containing Weyl spinors) from one with Dirac neutrinos (containing Weyl spinors) with matching masses and mixing angles. Nevertheless, once we add to a model with Dirac neutrinos small sources of lepton number violation, oscillation probabilities will change. Some illustrative examples are shown in fig. (1). We plot there the electron neutrino survival probability for lowenergy (reactor) neutrinos at distances up to (and slightly larger than) the typical distances of the KamLAND experiment [3]. In all plots the black lines show the expectation for the current global best fit point [4] for the ordinary neutrino parameters in the standard three generation case, to which we have added either a nonzero mass splitting to a Dirac state (top row) or one particular new quasiDirac angle (bottom row).
In section (2) we will discuss the general parametrization of masses and mixing angles for scenarios with three generations of quasiDirac neutrinos. However, from the examples shown in fig. (1) one can read off already some basic facts about oscillations of quasiDirac neutrinos, which we will work out in greater detail in section (3). First, small nonzero values of ’s are equivalent to introducing new, large oscillation lengths. Thus, the best constraints on will come from oscillation experiments with the largest possible baselines. And secondly, even if mass splittings are negligibly small, the new, nonstandard angles which appear in this setup (called above) may affect oscillation probabilities in a way similar to standard angles, hence creating parameter degeneracies. For example, as fig. (1) shows, from alone one cannot provide limits on a single angle. (In this example variations of can be compensated by varying .) Even by combining more than one oscillation probability, constraints can only be derived for certain combinations of angles and phases of the mixing matrix. We will discuss this in detail in section (3.2). Constraints on mass splittings are discussed in section (3.1).
A word on nomenclature. The terminologies quasiDirac and pseudoDirac neutrinos appear nearly interchangeably in the
literature. We prefer to define quasiDirac (QD) neutrinos as being a
mixture of active and sterile states, in contrast with pseudoDirac
(PD) neutrinos
With this definition, quasiDirac and pseudoDirac neutrinos are then
very different objects, both theoretically and phenomenologically. Let
us briefly mention that various aspects of pseudoDirac neutrinos have
been considered in the literature: Magnetic moments and double beta
decay [7], possible mass textures
[8, 9, 10], and oscillatory
behavior [11, 12, 13, 14].
We note in passing that models of pseudoDirac neutrinos require
neutrino mass matrices which no longer fit the solar and atmospheric
neutrino oscillation data [15, 16, 17].
Many more papers discussed the phenomenology of quasiDirac neutrinos. For example, double beta decay was first discussed in this context in [6], while [19] and [20, 21] consider quasiDirac neutrinos as a possible explanation of the atmospheric and solar neutrino problems, respectively. More ambitiously, explaining atmospheric, solar and LSND neutrino oscillations simultaneously was discussed in [22, 23]. However, all these proposals are by now ruled out experimentally, since they predict too much oscillations into sterile neutrinos. Limits on quasiDirac neutrino parameters, on the other hand, have been derived from solar neutrino data [24] as well as from solar, atmospheric neutrino data and cosmology [25]. Furthermore, in [26, 27, 28, 29] QD neutrinos have been discussed in the context of neutrino telescopes, such as IceCube.
QuasiDirac neutrino oscillations were also discussed in [30], where it was claimed that to leading order in the flavor composition of the mass eigenstates does not change (only mass splittings appear), hence oscillations for pairs of quasiDirac neutrinos are described by the standard mixing matrix. This assertion was taken to be true by others [26, 27, 28, 31], yet we want to stress that this claim is not correct, as can be seen from the eqs. (4)–(6). Already for one generation, these expressions show that the mass splitting and the departure from maximal mixing are both linearly dependent on and, more importantly, they are controlled by orthogonal combinations of these two parameters. As such, it is even possible to have no mass splittings at all and at the same time have arbitrary mixing angles.
There are also a number of more theoretical papers discussing how quasiDirac neutrinos could arise. One possibility is the socalled “singular” seesaw where the mass matrix for the singlet neutrinos () has a determinant equal or close to zero [32]. QuasiDirac neutrinos from such a singular seesaw with additional typeII seesaw contributions have been discussed in [33]. Another possibility [31] involves introducing additional singlets (), as it is done for the inverse seesaw mechanism [34]. A double seesaw is then responsible for producing very light states which, together with the active states, form quasiDirac neutrinos [31]. The authors of [35] use a Dirac seesaw to explain the necessary smallness of the Dirac neutrino mass terms first, and then generate quasiDirac states by the addition of a very small seesaw typeII term. The “mirror world” model of [36] is another way to obtain these particles.
In models with extended gauge groups quasiDirac neutrinos can also appear. An example is the inspired model of [37]. Here, several electroweak triplets of the gauge group are needed to accomodate the Standard Model leptons, and the observed active light neutrinos are automatically quasiDirac states [38]. A very different idea, based on supergravity has been discussed in [39]. There it was pointed out that if neutrino Dirac terms are generated from the Kähler potential (instead of the superpotential), neutrinos would be quasiDirac, since Majorana terms come from higher order Kähler potential terms and thus are expected to be suppressed. This idea [39] is particularly attractive, since it could, at least in principle, explain the observed smallness of the Dirac neutrino mass terms.
In addition to active neutrinos, models of Dirac neutrinos require the introduction of Weyl spinors transforming trivially under the electroweak gauge group. For this reason, the study of quasiDirac neutrinos necessarily has some overlap with the physics of sterile neutrinos. Many experiments have searched for sterile neutrinos. Most famously, the SNO neutral current measurement rules out dominant contributions of sterile neutrinos to the solar neutrino oscillations [40]. SuperKamiokande searched for steriles in atmospheric neutrinos [41]. OPERA [42], MINOS and DayaBay [43], IceCube[44] and NOA [45] published searches for sterile neutrinos. For a more complete list of references see the recent reviews [46, 47]. Note, however, that constraints on sterile are usually derived assuming best fit point values for the standard oscillation parameters, to which two new parameters (one angle and one mass splitting) are added in the fit. This approach does not cover the general quasiDirac neutrino parameter space. In particular, keeping the standard neutrino parameters fixed can lead to misleading conclusions about limits for the new/extra parameters.
There are also some hints for the existence of sterile neutrinos. However, all these hints point to a new and much larger mass scale in oscillations, i.e. (1) eV. Since these indications imply masses and mixings very different from those of the standard oscillations, they can not be explained by quasiDirac neutrinos. We thus do not discuss these hints any further and refer only to the recent review [47].
The rest of this paper is organized as follows. In section (2) we discuss the basics of quasiDirac oscillations, constructing general expressions for the mixing matrix for the three generation case. In section (3) we discuss constraints on the new, nonstandard parameters from various neutrino experiments. Constraints on quasiDirac mass splittings are discussed in section (3.1), while in section (3.2) we discuss the constraints on angles, for the case in which mass splittings are negligible. We then close with a short summary and discussion.
2 Definitions for quasiDirac neutrino oscillations
Dirac neutrinos can be described either in the weak or in the mass
basis. The two pictures are equivalent. We will choose the latter
one. Consider then a leptonnumber preserving model with three active
and three sterile neutrinos ( and ).
(7) 
In order to diagonalize the matrix , both active and sterile neutrinos must be rotated, and , such that :
(8) 
Strictly speaking, the neutrino mass matrix is not yet diagonal since it is still mixing different states (active and sterile neutrinos). This can be solved by rewriting and () as , and :
(9) 
where the masses and the mixing matrix have a special form ( is a square matrix):
(10)  
(11) 
If the pattern of masses and mixing in eqs. (10) and (11) is perturbed, neutrinos are no longer Dirac particles and lepton number is violated. Note that this is equivalent to switching on the lepton number violating masses and in eq. (1). We shall now look into the possible departures from the Dirac limit as seen from the mass basis.
In the case of masses, it is possible to split the three pairs of , hence we may introduce three such that
(12) 
with the understanding that, for quasiDirac neutrinos, the are small in comparison to the atmospheric and solar mass scales. In total there are now five mass parameters relevant for oscillation experiments: the usual and , plus three new mass splittings. (As usual, the overall mass scale of neutrinos does not enter the oscillation probabilities.)
Let us now turn our attention to a generic mixing matrix with dimensions . Such a matrix can be described by real numbers, yet orthonormality of rows () imposes conditions on them, and furthermore it is possible to absorb phases into the charged lepton fields, hence there is a total of real physical degrees of freedom in . For a matrix, this corresponds to 12 angles and 12 phases, but note that 5 of these phases cannot be observed in neutrino oscillation experiments (they correspond to column phases). The matrix can be explicitly parametrized as follows [48] (called below the SV parametrization). First, consider an elementary rotation in the entries given by the complex number such that, in the (1,2) case, it has the form
(13) 
In the SV parametrization, the th row of () is then given by the expression
(14) 
where is a column vector with entries . We do not give here in full because the expression is very lengthy.
For a particular arrangement of the 24 angles and phases in eq. (14), takes the special form (11) which is associated with the Dirac limit. Note that, as usual, one can write the square matrix with three angles and one phase:
(15) 
Unfortunately, it is very complicated to describe the Dirac limit in the SV parametrization. Hence we make a small modification by introducing the following rotation matrix:
(16) 
With this definition, the mixing matrix in eq. (11), with parametrized as in eq. (15), corresponds to , as in (16), with (), and . In other words, with this definition the Dirac limit for simply corresponds to keeping only the standard three generation neutrino mixing angles nonzero.
We can then write the probability of neutrino oscillation from a flavor
to a flavor for an energy and after a length
as:
(17) 
Note that this expression is insensitive to column rephasings . It is easy to show that eq. (17) reduces to the standard oscillation formula in the Dirac limit.
3 Current experimental limits and future prospects
As discussed in the previous section, the full parameter space for a system of 3 pairs of QD neutrinos has 30 free parameters: Two independent plus one overall mass scale, three , twelve angles and twelve phases. Even discounting the five Majorana phases and the overall mass scale, which can not be probed in oscillation experiments, the remaining number of parameters is much too large to fit simultaneously.
Nearly all experimental data, on the other hand, is consistent with the standard picture of only three active neutrino species participating in oscillations [4], i.e. two mass squared differences ( and ), three mixing angles (, and ) plus one phase () are sufficient to describe the data. As mentioned in the introduction, there are also some hints for sterile neutrinos with a mass scale of the order of (eV) [46, 47]. However, all these hints are at most of the order of () , we will thus not take them into account in the following. Instead, since the standard three generation picture seems to describe the data well, we will consider “small” perturbations and derive limits on particular combinations of nonstandard parameters.
In order to deal effectively with the large number of parameters controlling quasiDirac neutrino oscillations, we will consider two simplified scenarios:

First, we take one nonzero at a time. In these oneparameter extensions, very stringent limits on are found, in agreement with earlier analysis, see for example [24, 25]. We then extend this analysis to two new parameters: One mass splitting plus one new angle. This second step allows us to identify “blind spots” in the oscillation experiments, i.e. degenerate minima in particular directions in parameter space, where limits on mass splittings are much worse than in the one parameter fits. We then discuss a particular parametrization of these degenerate directions in parameter space, where the effects of can be decoupled from oscillation experiments nearly completely.

In the second setup, we discuss the limit where mass splittings are too small to be measured in oscillation experiments, hence there are just angles and phases of to deal with. In this situation, it can be shown that from the 24 parameters in only 13 combinations enter the oscillation probabilities of active neutrinos. Moreover, since there is only very limited information on oscillations involving , we can in practice restrict ourselves to experiments involving ’s and ’s. There are then only 7 combinations of the 24 angles and phases which appear in the oscillation probabilities. We discuss the construction of these 7 quantities, the current constraints and possible tests for quasiDirac neutrinos in this limit.
In our analysis we do not take into account the data from every existing oscillation experiment. Given the scarcity of data on neutrinos, we ignored it altogether, concentrating instead on the available charged current data for and neutrinos and antineutrinos. Also, we focus on those experiments, which should provide the most important constraints for quasiDirac neutrinos. First, we consider KamLAND [3] since it fixes most accurately the socalled solar mass splitting . From the solar neutrino experiments we fit SuperK elastic scattering data [49] and from Borexino the measured [50] and fluxes [51]. This choice is motivated by the fact that SuperK [49] has the most accurate data in the high energy range, while [50, 51] fix the low energy part of the solar neutrino spectrum.
We take the atmospheric neutrino data from [52]. We concentrate in our fit on the subGeV sample, since the lowest energetic neutrinos will be most sensitive to small values of , as demonstrated in fig. (2). In addition, we use data from the MINOS collaboration (muon and antimuon neutrino survival) [53] and T2K (muon neutrino survival and muon to electron neutrino transition) [54], since these two experiments determine better than atmospheric data. And, finally, we take into account data from DayaBay [55], since from the three current reactor neutrino experiments DayaBay determines with the smallest error.
In our fits, we use a simple method to determine the allowed ranges of model parameters. We take statistical and systematic errors from the experimental publications, to which we added a further (small) systematic error for the uncertainties in our theoretical calculations. This latter systematic error was chosen such that our simulations reproduce the allowed parameter ranges for the standard oscillation parameters, determined by the respective experiment, within typically (11.5) c.l. ranges. Note that we do not attempt to do a precision global fit for standard neutrino oscillation parameters. Rather, we consider reproducing the experimental results for the standard case as a test for the reliability of our derived limits.
3.1 Limits on mass splittings
One parameter limits
We will first discuss limits derived on assuming one at a time and taking all nonstandard angles to be zero. Table (1) shows limits on and the corresponding experimental data sets used to derive the limits.
Experiment  [eV]  [eV]  [eV] 
KamLAND  –  
Solar + KamLAND  –  
DayaBay + MINOS + T2K  –  
SuperK + DayaBay + MINOS + T2K  –  
JUNO 
For each case listed in table (1), we have calculated the upper limits on the twice: (a) marginalizing over two of the standard neutrino oscillation parameters and (b) for the best fit point value of the standard parameters. Marginalization over standard oscillation parameters leads to less stringent limits. However, the importance of this marginalization procedure differs widely for different experiments. For example, in the case of KamLAND, bounds on of the order of roughly are derived marginalizing over the allowed ranges of and , while for the best fit values of these last two parameters, the limits are more stringent by “only” roughly a factor 2.
As the table shows, the strongest constraints on and
come from solar neutrino data. This is easily
understood from fig. (3), which shows the electron
neutrino survival probability as a function of neutrino energy for
different values of and for the best fit values of
the standard solar oscillation parameters. For low values of neutrino
energies, vacuum oscillations dominate and so very small
can be probed up to a scale essentially determined
by the EarthSun distance (). Note
that a nonzero reduces , but this reduction
could be hidden in the relatively large error bar of the lowenergy
measurements.
Solar neutrino experiments have essentially no sensitivity to . This is simply due to the smallness of ( [4]). Thus, we have to rely on experiments testing the atmospheric scale to derive limits on . Table (1) quotes numbers for two cases.
In the first scenario, we have combined data from DayaBay [55], T2K [54] and MINOS [53]: DayaBay fixes most accurately , while both MINOS and T2K measure with rather small errors. Here, the limit on is (not surprisingly) less stringent than the one derived from KamLAND (or solar). Depending on whether or not and are used at their best fit value or marginalized over, we get very different limits on . This is due to the fact that when scanning over the standard oscillation parameters, the function has two almost degenerate minima: one for small values of and another for of the order eV. However, as the table also shows (second case), this nonstandard solution is excluded, once we add SuperK atmospheric neutrino subGeV data to the fit. With the combination of these four experiments limits on and are again of order eV.
In the last line of table (1) we give our forecast of the sensitivity of the planned experiment JUNO [56]. JUNO will measure and very precisely and thus it will also be able to derive limits on any . However, our results indicate that, despite being a very precise experiment, JUNO will not lead to a major improvement over existing limits on . Here, it is important to stress that limits using the best fit point and limits marginalizing over standard parameter uncertainties are very different. This can be traced back again to a neardegeneracy in the function: For of the order of one has two only slightly different oscillation lengths contributing in the fit, which can give a better description than a single oscillation length.
In summary, strong limits on mass splittings can be derived from atmospheric and solar neutrino data ( eV and eV, respectively) in the case where no other extra parameter is added to the standard neutrino oscillation picture.
Two parameter case
While the discussion in the previous subsection seems to show that constraints on QD mass splittings are very stringent, we will now see that this conclusion is valid only under the assumption that no other nonstandard parameter is different from zero.
As a simple example, consider the electron neutrino survival
probability at distances short enough that the effects of can be neglected.
(18) 
where and are shorthands for and and
(19) 
It is straightforward to see that the above expression for the neutrino survival probability remains (nearly) unchanged if we swap by . This is true up to very small terms proportional to . More specifically, in the limit where have the same value, is only a function of the combination . Thus, there will be a neardegeneracy of the relevant function involving these two angles, and so values (or limits) derived for one of these parameters, without varying the other, will be misleading.
There is, however, another more interesting degeneracy associated to eq. (18). In calculating this expression we have used a certain parametrization for the mass splitting, which we may call the symmetric parametrization: . Choosing , the second term inside the bracket in eq. (18) vanishes (this choice corresponds to ). So, by adjusting and we can keep constant and equal to the best fit point value of , in which case there will be no upper limit on itself coming from the electron neutrino survival experiments.
Note that we could have defined .
We can break this particular degeneracy in parameter space, by adding more experiments. T2K measures two probabilities: (a) The muon neutrino survival probability, , and (b) the electron neutrino appearance probability, , both at values of which give access to the atmospheric neutrino mass scale, . If the only nonstandard angle different from zero is , then will not depend on at all, while will have dependence which is different from the one of . Thus, adding T2K data to the scan is enough to break the degeneracy in hence an upper limit on reappears. The middle column of fig. (4) illustrates this point; it shows the results of a combined scan over , and for DayaBay plus T2K data. By comparison of the right with the middle column of fig. (4), one can clearly see that the addition of SuperK data generates a strong upper limit on , for this particular choice of parameter subspace.
Given these results, one might wonder if there are particular directions in parameter space for which oscillation experiments become completely blind to QD neutrino mass splittings. Recall that the blind (or: degenerate) direction discussed above for corresponds to the particular choice of . In a similar way, for example, would loose any sensitivity to if . Thus, with some special choice of and such that both and are zero at the same time, one can indeed make DayaBay and T2K blind to variations of .
While it is possible, in principle, to calculate the combination of angles and phases (defined in section (2)) associated to these blind directions, in the following we will consider a simpler alternative. Consider a unitary rotation of the columns and of the mixing matrix . Since we are not interested in column phases, such a rotation is governed by just two parameters ( and ):
(20) 
Now, recall that in the Dirac limit (see eq. (11)) the columns and of the mixing matrix are proportional to each other, . This means that applying the transformation to the Dirac neutrino mixing matrix and, without loss of generality setting , we obtain:
(21) 
From the last of these equations it can be seen that the ’th column (the ()’th column) of vanishes, if one chooses ().
Fig. (5) shows a scan over the allowed parameter
space in the plane () using DayaBay, T2K,
MINOS and SuperK atmospheric neutrino data. In agreement with the
above discussion, there is a blind spot where no limit on
exist. This blind direction corresponds to the
choice of .
We now turn to a discussion of and
. For these two parameters, again solar neutrino
physics provides the most important constraints. As above for
, we can define a rotation angle
() between the columns 1 and 4 (2 and 5) of the mixing
angle which will mitigate the effects of a nonzero
(). Fig. (6) shows the
probability for solar neutrinos as a function of neutrino energy for
two different values of and various values of
.
Fig. (7) shows the allowed parameter space in the two planes () and () using solar data and combining solar data with KamLAND. These plots have been calculated using the best fit point values for and from the global fit [4]. The plots show in all cases that there exists a slight preference, between (22.5) in all cases, for nonzero values of . Note that the preferred solution of the solar data in the region of and is ruled out by KamLAND. However, even combining solar and KamLAND data some preference for nonzero of the order of very roughly eV remains.
We have traced back this preference for nonzero mass splittings in solar data to the wellknown difference in the best fit points from in solar and KamLAND data. As can be seen also in the latest global fits [4], solar data prefers a around eV, while KamLAND prefers eV. This tension between the two data sets is roughly of the order of 2 , with the error bar dominated by the larger error on in the solar data set. We have therefore recalculated the constraints on () from solar data for a value of eV. Fig. (8) shows the results of such a scan. As can be seen, in this calculation there is no longer any preference for a nonzero value of .
Note that such a low value of is ruled out by many from the KamLAND data. Thus, a small nonzero mass splitting could provide, in principle, a solution for the observed tension between solar and KamLAND data.
In summary, by introducing one mass splitting at a time, we extracted bounds for the pairs of parameters , . In the limit where is , one column of the mixing matrix vanishes and therefore the mass splitting becomes unobservable. For this reason, one expects that for reasonably large values of there must be tight limits on , meaning that has to be quite far from the Dirac limit (). On the other hand, for a small enough value of the mass splitting , the associated oscillation length eventually become larger than the baseline of the relevant experiments, and in that case becomes unconstrained.
3.2 QuasiDirac neutrinos in the limit
If masses are degenerate as in eq. (10), then the oscillation probability formula will not change under unitary rotations of the columns and of the mixing matrix, see eq. (20). In other words,
(22) 
for unitary matrices () leaves
unchanged. Hence,
there is a redundancy in our description of , and
in turn this means that out of the 24 parameters describing the mixing
matrix, oscillation experiments are only sensitive to 13.
(23)  
(24)  
(25) 
with the oscillating factors and and the 7 parameters defined as follows:
(26)  
(27)  
(28)  
(29) 
As a side remark, we would like to point out here that a similar approach could, in principle, be used in the presence of
one : in this case, instead of 7, there would be 9 combinations of angles and phases to take into account.
Note that the defined above can take any value in our framework, provided that the following constraints are obeyed:

the first six are nonnegative real numbers;

neither nor can be larger than 1;

and ;

the norm of is fixed by and (), so even though is a complex parameter, only is an independent degree of freedom;

cannot be bigger than .
These conditions are a consequence of the definitions of the and the fact that the rows of the mixing matrix are orthonormal (). Taking them into account, we are able to pick out all valid points in the parameter space, without ever referencing back to specific entries of the mixing matrix.
For reference, the values of these parameters in the Dirac limit as a function of the standard , , and parameters (see eq. (15)) are the following:
(34)  
(35)  
(36)  
(37) 
Using eV and , we performed a 7dimensional scan over all . Electron neutrino survival data from at KamLAND, DayaBay, SuperK and Borexino was used, together with muon neutrino survival data at MINOS and T2K and T2K data. The allowed values in the planes , and are shown in fig. (9). The marginalized function for is shown in fig. (10); the marginalized function is not shown as it is essentially flat from 0 to .
Overall, the bounds on these 7 parameters are broadly consistent with the standard three neutrino oscillation picture. In other words, by substituting in expressions (34)–(37) the numbers obtained for , , and from global fits [4], we get values for the roughly in agreement with figs. (9) and (10). In order to see clearly that current data is consistent with the Dirac limit, note that in this latter case there are only 4 independent parameters. Thus, it follows that the standard three neutrino oscillation picture must correspond to three relations among the 7 . These are
(38) 
and from fig. (11) one can see that oscillation data is compatible with each of these equalities within . The three together are disfavored only at so, assuming no mass splittings , there is currently no significant indication for quasiDirac neutrinos.
We would like to point out that even in the absence of new mass scales, quasiDirac neutrinos can, in principle, be distinguished from some other scenarios through oscillation experiments. In particular, consider 3 active neutrinos and a nonunitary mixing matrix . The oscillation probabilities are then given by the expressions
(39)  
(40)  
(41) 
where . The exact form of the functions and which control the 0distance neutrino behavior is not important for the present discussion. The more important point is that with a nonunitary one can have an oscillatory behavior which is impossible to reproduce with quasiDirac neutrinos, and viceversa.
For example, the coefficients in do not need to be related for a nonunitary , while for quasiDirac neutrinos they must be the same (up to a minus sign — see eq. (25)). On the other hand, note that from the and one can extract the modulus of the absolute value of all , hence by measuring and , as well as the coefficients in , the coefficients of the oscillatory factors in are fixed for a nonunitary (up to signs). However, for quasiDirac neutrinos no such constraint exists. So, with this short theoretical argument, one can conclude that in principle these two nonstandard neutrino scenarios can be distinguished through oscillation experiments.
4 Summary
In general, neutrinos can have lepton number violating (Majorana) and lepton number conserving (Dirac) mass terms. If the lepton number violating mass terms are smaller than the lepton number preserving ones, neutrinos are quasiDirac particles. Phenomenologically, this corresponds to the existence of three pairs of neutrinos with slightly different masses, hence oscillation experiments are sensitive not only to the usual solar and atmospheric mass scales, but also to three small mass splittings . Furthermore, for quasiDirac neutrinos there are more than 3+1 angles and phases to be considered. In this work, we have analyzed the constraints on these quasiDirac neutrino parameters imposed by current neutrino oscillation data and also briefly discussed the potential of the future JUNO experiment to improve upon existing constraints.
In section (2) we have discussed a fully general parametrization of the lepton sector for three generations of quasiDirac neutrinos. In addition to the charged lepton masses, there is a total of 6 masses, 12 angles and 12 phases. Oscillation experiments are not sensitive to the overall neutrino mass scale nor to 5 of the phases (which are of the Majorana type). Hence we are left with a 24dimensional model space, compared to the sixdimensional space for an ordinary three generation case (, , , , and ).
It is numerically too costly to handle such a large number of parameters at the same time, hence we analyzed several different special cases. First, we took a single mass splitting . If we split two neutrinos with mass into a quasidegenerate pair of particles with masses and then a new oscillation length appears which is associated to the conversion of active to sterile neutrinos. Very stringent limits on in such one parameter extensions can be derived, of the order of for (from solar neutrino data) and