Texture zeros of lowenergy Majorana neutrino mass matrix in 3+1 scheme
Abstract
In this work we revisit the zero textures in low energy Majorana neutrino mass matrix when the active neutrino sector is extended by a light sterile neutrino in the eV scale i.e., the 3+1 scheme. In 3+1 scenario, the low energy neutrino mass matrix () has ten independent elements. Thus in principle one can have minimum onezero texture to maximum tenzero texture. We summarize the previous results of one, two, three and fourzero textures which already exist in the literature and present our new results on fivezero textures. In our analysis we find that among six possible fivezero textures, only one is allowed by the present data. We discuss possible theoretical model which can explain the origin of the allowed fivezero texture and discuss other possible implications of such a scenario. Our results also concludes that in 3+1 scheme, one can not have more than fivezeros in .
pacs:
12.60.i,12.60.Cn,14.60.PqI Introduction
The possibility of light sterile neutrinos with mass at the eV scale have gathered serious attention in the last two decades following the neutrino anomalies reported by some experiments which could be explained by incorporating additional light neutrinos to which the active neutrinos can oscillate into. For a review, one may refer to Ref. Abazajian:2012ys (). The first such anomaly was reported by the Liquid Scintillator Neutrino Detector (LSND) experiment in their antineutrino flux measurements Athanassopoulos:1996jb (); Aguilar:2001ty (). The LSND experiment searched for oscillations in the appearance mode and reported an excess of interactions that could be explained by incorporating at least one additional light neutrino with mass in the eV range. This result was supported by the subsequent measurements at the MiniBooNE experiment AguilarArevalo:2013pmq (). Similar anomalies have also been observed at reactor neutrino experiments Mention:2011rk () as well as gallium solar neutrino experiments Acero:2007su (); Giunti:2010zu (). Since the precision measurements at the LEP experiment do not allow additional light neutrinos coupling to the standard model (SM) gauge bosons Agashe:2014kda (), such additional light neutrinos are called sterile neutrinos. These anomalies require the presence of a light sterile neutrino at eV scale having nontrivial mixing with the active neutrinos as presented in the global fit studies Kopp:2013vaa (); Giunti:2013aea (); Gariazzo:2015rra ().
Apart from these reactor and accelerator based experiments, there were initial hints from cosmology as well, suggesting the presence of one additional light neutrino. For example, the nine year Wilkinson Microwave Anisotropy Probe (WMAP) data suggested the total number of relativistic degrees of freedom to be Hinshaw:2012aka (). Since the standard value is , the WMAP data could accommodate one additional light species. Such cosmology experiments can constrain the number of such relativistic degrees of freedom as they affect the big bang nucleosynthesis (BBN) predictions as well as cause changes in the cosmic microwave background (CMB) spectrum, which are very accurately measured. Contrary to the WMAP findings, the more recent Planck experiment puts limit on the effective number of relativistic degrees of freedom is Ade:2015xua ()
(1) 
which is consistent with the standard value . Here the keywords in parenthesis refer to different constraints imposed to obtain the bound, the details of which can be found in Ref. Ade:2015xua (). The Planck bound is clearly inconsistent with one additional light neutrino. Although this latest bound from the Planck experiment can not accommodate one additional light sterile neutrino at eV scale within the standard CDM model of cosmology, one can evade these tight bounds by considering the presence of some new physics beyond the standard model (BSM). For example, additional gauge interactions in order to suppress the production of sterile neutrinos through flavour oscillations were studied recently by the authors of Hannestad:2013ana (); Dasgupta:2013zpn (). Recently, the IceCube experiment at the south pole has excluded the three active and one sterile neutrino (the framework where the sterile state is heavier than the active states GomezCadenas:1995sj (); Goswami:1995yq ()) parameter space mentioned in global fit data Kopp:2013vaa () at approximately confidence level TheIceCube:2016oqi (). However, in the presence of nonstandard interactions, the neutrino global fit data can remain consistent with the IceCube observations Liao:2016reh (). Therefore, there is still room for existence of an eV scale sterile neutrino within some specific BSM frameworks that can provide a consistent interpretation of experimental data. The interesting cosmological implications of such light sterile neutrinos can be found in the recent review article Abazajian:2017tcc () and references therein.
Apart from finding a consistent neutrino framework compatible with short baseline neutrino anomalies as well as cosmology, another challenge in particle physics is to explain the origin of this light sterile neutrino and its non trivial mixing with the active neutrinos. Apart from explaining the eV scale mass of sterile neutrino, it is also desirable that the particle physics model predicts some of the neutrino parameters that can undergo further scrutiny at ongoing neutrino oscillation experiments An:2016luf (); MINOS:2016viw (); Adamson:2016jku (). Typically, a BSM framework for explaining neutrino masses and mixing comes with a large number of free parameters lacking predictability. However, if the theory has a well motivated underlying symmetry that gives rise to a very specific structure of neutrino mass matrix, then number of free parameters can be significantly reduced. Here we consider such a possibility where an underlying symmetry can restrict the mass matrix to have nonzero entries only at certain specific locations. Such scenarios are more popularly known as zero texture models, a nice summary of which within three neutrino framework can be found in the review article Ludl:2014axa (). The light neutrino mass matrix in framework is a complex symmetric matrix, assuming the neutrinos to be Majorana fermions. Such a mass matrix can be parametrised by sixteen parameters: four masses, six angles and six phases. In the presence of zero textures, these parameters get related to each other through the zero texture equations resulting in more constrained set of parameters or more predictability. Recently, the possibilities of such zero textures were explored in the framework in Ref. Ghosh:2012pw (); Ghosh:2013nya (); Zhang:2013mb (); Nath:2015emg (); Borah:2016xkc (). The authors in these works pointed out the allowed zero texture mass matrices containing up to fourzeros in framework from the requirement of satisfying recent data of masssquared differences and mixing angles. In present paper we briefly summarise all previous works and also extend them to study the possibility of having five and six zero texture mass matrices. Since the simultaneous existence of zeros in active and sterile sectors is phenomenologically disallowed Ghosh:2012pw (), six is the maximum number of possible zeros in the light neutrino mass matrix. Therefore, our present study is going to give a complete picture of all possible zero texture mass matrices in framework. It should be noted that we stick to a diagonal chargedlepton basis for simplicity and hence all our conclusions are valid in this basis only.
After summarising the earlier works on zero texture mass matrices upto fourzeros, we show that one possible fivezero texture mass matrix is allowed from the present neutrino data while the possibility of sixzero texture is ruled out. Apart from finding the predictions for different neutrino parameters in this particular fivezero texture mass matrix, we also point out one possible symmetry realisation that can naturally generate such a mass matrix. This is based on an abelian gauge symmetry where the relative difference between second and third generation lepton number is gauged. We also discuss other interesting implications of such a scenario related to the anomalous magnetic moment of muon. We also discuss one interesting discrete symmetry which the fivezero texture mass matrix possesses partially and its possible implications.
The paper is organized in the following way. In Section II, we discuss the low energy neutrino mass matrix in the 3+1 framework. In Section III we give a brief summary of the past results on one, two, three and fourzero textures. In Section IV we present our new results on fivezero texture. Section V will contain the theoretical model which explain the origin of the allowed fivezero texture and Section VI contains some possible phenomenological implications of our results. Finally we will conclude in Section VII.
Ii Neutrino mass matrix in 3+1 framework
In presence of an extra sterile neutrino having mass in the eV scale, there will be two possible mass ordering of the neutrinos: Normal hierarchy (NH) i.e., and inverted hierarchy (IH) i.e., , where , , are the masses of the active neutrinos and is the mass of the sterile neutrino. They can also have quasidegenerate spectra (QD) if . Irrespective of the mass spectrum of the neutrinos, the low energy neutrino mass matrix in the 3+1 scheme can be expressed as
(2)  
(3) 
where and is the unitary PMNS matrix which contains six mixing angles i.e., , , , , , , three Dirac type CP phases i.e., , , and three Majorana type CP phases i.e., , , . is the diagonal Majorana phase matrix given by and we parametrize as
(4) 
where are the rotation matrices and can be expressed as
(5)  
(6) 
and so on, with and are the Dirac CP phases. In this parametrization, the six CP phases vary from to .
Iii Previous results of zero textures in in 3+1 scheme
In this section we will discuss briefly the previous results of one, two, three and fourzero textures. Onezero texture in the neutrino mass matrix is given by the condition
(7) 
where , are the flavour indices. Thus there exist ten possible onezero mass matrices. Onezero mass matrices in 3+1 scheme have been discussed in Refs. Ghosh:2013nya (); Nath:2015emg (). The main results of these works are:

All the onezero textures except and are allowed.

The texture is allowed in both NH and IH. This is in sharp contrast to the standard three flavour case where is not allowed in IH. Note that the texture zero condition of in 3+1 scenario depends on the value of . Specifically in this case is a rising function of the lowest neutrino mass. The condition also strong constrains Majorana phase to be around . Thus any two, three and fourzero texture that involve will predict as a rising function of lowest mass and around .

The textures , , , and are allowed in both NH and IH.

The texture is only allowed if the neutrino masses are quasidegenerate.
Before moving forward let us discuss a bit more about the elements in the fourth row/column of . For the elements , , and , the leading order term looks like , , and respectively. Now as is quite large ( eV), the coefficient of needs to be small to obtain (with , , and ). Thus we see that can never be zero. In the analysis of Refs. Ghosh:2013nya (); Nath:2015emg (), the mixing angles and are bounded from below but can be as small as zero. This is the reason why they have concluded that is not allowed and vanishes only in the quasidegenerate regime. But if one assumes the values of and close to zero are allowed, then the conclusions about and may change.
Twozero textures in neutrino mass matrix are obtained when two of the matrix elements are zero simultaneously. In 3+1 scheme the twozero textures are discussed in Ref. Ghosh:2012pw (). The number of possible twozero textures in is 45. Among the 45 cases, there are 30 cases in which the zero texture includes (where , , and ). It was shown that any texture zero which includes an element corresponding to the fourth row or fourth column of the mass matrix is not allowed ^{1}^{1}1In this analysis, the mixing angles and are considered to be bounded from below. However, numerical analysis reveals that even when and are close to zero, the twozero textures involving the elements corresponding to the fourth row or fourth column are not allowed.. Thus we are left with the 15 twozero textures listed in Table 1. Here it is interesting to note that these 15 textures coincide with the 15 possible twozero textures in the standard three generation in the absence of the sterile neutrino. There are numerous studies in the literature which discuss the viability of these 15 twozero textures in 3 generation Xing:2004ik (); onezero (); onezero1 (); onezero2 (); onezero3 (); alltex (); texturesym9 (); twozero (); twozero1 (); twozero2 (); twozero3 (); twozero4 (); twozero4_1 (); twozero5 (); twozero6 (); twozero7 (); twozero8 (); twozero9 (); twozero10 (); twozero11 (); Gautam:2016qyw (). All these analyses show that among these 15 twozero textures in three generations, only seven textures are allowed. These allowed textures belong to , and class as given in Table 1. But the conclusions in these studies change when one includes a light sterile neutrino in addition to the three active neutrinos. The analysis of Ref. Ghosh:2012pw () shows that in the 3+1 scheme only the textures in class are allowed in NH whereas textures belonging to all the other classes (i.e., , , , and ) are allowed in both NH and IH. This is the most remarkable finding of this work. The textures which were disallowed in the threegeneration case become allowed from the contribution of additional terms from the sterile sector.
The analysis of zero textures when three elements of the neutrino mass matrix are simultaneously zero can be found in Ref. Zhang:2013mb (). Note than in the standard threegeneration case the maximum allowed numbers of zero textures are two while threezero textures are phenomenologically disallowed. Thus the possibility of having more than two zeros in is a special feature of the 3+1 scheme. In the 3+1 scenario, there can be 120 possible threezero textures. But among them 100 threezero textures contain an element belonging to the fourth row/fourth column of the mass matrix and hence they are not allowed. The 20 textures are classified in six sets as follows:
(8)  
Numerical analysis shows that the textures , , , , , , , , , and are allowed in both NH and IH whereas the remaining seven textures i.e., , , , , , and prefer NH over IH. The texture is excluded almost completely, being allowed in a very small part of the parameter space for IH.
An analysis of fourzero textures, when four elements of the low energy Majorana neutrino mass matrix can vanish simultaneously is carried out in Ref. Borah:2016xkc (). In this case the number of possible zero textures are 210. But again out of these 210 textures, 195 are readily ruled out because for these textures we have . The remaining 15 textures are then classified in either group where the texture zeros contain or group in which .
(9)  
Numerical analysis reveals that the textures , and are the disallowed in NH whereas and are the disallowed in IH.
Note that results of the one, two, three and fourzero neutrino mass matrix textures are largely consistent with each other. In Ref. Ghosh:2012pw (), it was shown that the twozero textures and are disallowed in IH. In class the elements and vanish whereas in class we have . Thus this work predicts any simultaneous zero texture involving or with will be disallowed in IH. In the analysis of threezero textures, the classes where we obtain or are , , , , , and . According to the analysis of Ref. Zhang:2013mb () though these cases are preferred in NH over IH, they are not completely ruled out in IH. This difference is mainly due to the choice of different ranges of the sterile mixing parameters, in particular , and . In the fourzero case the textures where we have the condition or are , , , , , , , and . According to the analysis of Ref. Borah:2016xkc (), all other textures except and are allowed in IH. This difference occurs because in Ref. Borah:2016xkc (), there are no lower bounds on and and these can be as low as zero. Thus we understand that the viability of zero textures in 3+1 case is extremely sensitive to the choice of the ranges of activesterile mixing parameters.
Iv Results of fivezero textures
In this section we present our results for the fivezero textures. The fivezero texture condition can be expressed mathematically as
(10) 
where , , and are functions of the mixing angles and phases. Note that Eq. 10 is a set of five complex equations. To obtain the real set of equations we put the real part and the imaginary part individually to zero:
(11) 
Parameters  Allowed ranges  Parameters  Allowed ranges 
3  180 to 180  
3  86 to 86  
45  180 and 180  
3  180 to 30  
30 to 180  
3  35.5 to 35.5  
3.4  11  180 to 137  
137 to 180  
0.018  0.22 eV  0.7  2.5 eV 
Now Eq. 11 is a set of ten real equations relating sixteen independent parameters i.e. six mixing angles, six phases and four masses. To solve this set of equations we supply the input values of three active neutrino mixing angles (i.e., , and ), and three mass square differences (i.e., , and ). The ten equations are solved for the remaining ten parameters using Mathematica (which internally implements the multidimensional NewtonRaphson algorithm). For NH, we have expressed , and as , and respectively, keeping the lowest mass free, whereas for IH we have expressed , and as , and keeping the lowest mass free. We have varied our input parameters for the threegeneration parameters in the allowed range as given the global analysis of the world neutrino data Forero:2014bxa (); Esteban:2016qun (); Capozzi:2013csa () and varied from 0.7 eV to 2.5 eV. If the output of , and falls between to , to and to respectively Kopp:2013vaa (); An:2014bik (); Adamson:2011ku () with the condition , then we say this texture is allowed in NH (IH)^{2}^{2}2Note that according to the global analysis of the shortbaseline data Kopp:2013vaa () we have and at . However the Refs. An:2014bik (); Adamson:2011ku (), give only an upper limit on and as they analyse standalone data. Thus for a conservative approach, in our analysis we have taken the upper limits of and from the global analysis and allowed them to have lower limits as zero.. In 3+1 scenario, the number of possible fivezero textures are . But among them, 246 appears with the with one of the elements belonging to the fourth row/column and thus they are not allowed. The remaining six possible fivezero textures are:
(12)  
Among these six possible structures, our analysis shows that only the texture is allowed in NH, and all the textures are ruled out in IH. In Fig. 1, we have given the correlation plots for the allowed texture in NH. In the texture , we have the condition . As mentioned in the previous section, the property that any zero texture mass matrix having the condition , will have as a rising function of the lowest mass and the Majorana phase will be constrained around , as clearly seen in Fig. 1 (top right panel). This property can be simply understood by looking at the expression of .
Note that in the above equation the sterile term is given by . Now the condition of is simply obtained by the cancellation of active and sterile terms. Therefore it is easy to understand that when is small (large) we need smaller (larger) values of to achieve cancellation. At the same time, for the cancellation of sterile and active terms, the coefficient of the sterile term must acquire a negative sign which is only possible if the phase is around . From Eq. IV we also infer that leads to the standard threeflavor neutrino mixing scenario. In that case the lowest mass cannot be zero in order to produce zero texture at . The vanishing of the activesterile mixing angles results in the reduction of expressions of , and from fourneutrino mixing to threeneutrino mixing scenarios. Earlier studies Ghosh:2013nya (); Nath:2015emg () have shown that the activesterile mixing angles , and and the lowest mass cannot simultaneously vanish in order to produce zeros textures , and . We summarize the allowed ranges of neutrino oscillation parameters in Table 2 for texture . As seen from Fig. 1, very small and are disallowed for the current texture under consideration (top left panel). The same figure also shows the constraint on the lightest neutrino mass from the upper bound on the sum of absolute neutrino masses eV given by the latest data from the Planck mission Ade:2015xua (). The bottom left panel shows that . This can be attributed to the fact that a very large value of negates the existence of zero texture at . A rigorous analysis of zero textures at each element in neutrino mass matrix and the interdependency of neutrino parameters is done in Refs. Ghosh:2013nya (); Nath:2015emg (), and their results apply here. From the bottom right panel of Fig. 1, we also see that the maximal value of is disallowed in this texture. Note that the features discussed above are of great importance to probe this texture in the future generation oscillation experiments. For example, if the future experiments measure or or not in the region to or , then this texture can be readily ruled out and hence the possibility of having 5 zeros in .
From the earlier results it is obvious that fivezero textures are the maximum which is allowed in the 3+1 scheme. A texture containing more than fivezeros in the low energy neutrino mass matrix is not allowed in the 3+1 scenario.
V Flavor symmetry origin of fivezero texture
Fermion Fields  

Scalar Fields  
Since the zero textures appear only in the active block of the mass matrix, they can be explained by different flavor symmetry frameworks. Some possible models are discussed in Refs. texturesym (); texturesym1 (); texturesym2 (); texturesym3 (); texturesym4 (); texturesym5 (); texturesym6 (); texturesym7 (); texturesym8 (); texturesym9 () in the context of zero textures in the threeneutrino picture. Since the sterile neutrino is a singlet under the standard model gauge symmetry, one can not prevent a bare mass term of Majorana type as well as a Dirac mass term involving the active neutrinos and the Higgs field. For a neutrino mass matrix at eV scale, we should be able to keep both the Majorana and the Dirac mass term involving the sterile neutrino at the eV scale, which is unnatural unless some additional symmetries can ensure the smallness of these mass terms. This gives rise to another challenge in addition to generating the active neutrino mass matrix at subeV from the popular seesaw mechanism. Generating a light neutrino mass matrix within different seesaw frameworks have led to several studies in recent times sterileearlier1 (); sterileearlier2 (); sterileearlier3 (); sterileearlier4 (); sterileearlier5 (); sterileearlier6 (); sterileearlier7 (); sterileearlier8 (); sterileearlier9 (); sterileearlier10 (); Borah:2016xkc (); Borah:2016lrl (); sterileearlier12 (). Here we consider a simple extension of the idea proposed in Refs. sterileearlier2 (); sterileearlier3 () in order to accommodate the texture zero criteria or fivezeros in the active block of the light neutrino mass matrix. Instead of giving an effective model based on higher dimensional operators and discrete symmetries as was done in order to explain the fourzero texture mass matrix in Ref. Borah:2016xkc (), here we give a renormalizable model.
The particle content of the proposed model is shown in Table 3. We are showing only the fields responsible for neutrino mass generation here skipping the details of quarks and the charged lepton sector. The gauge symmetry of the standard model is extended by another gauge symmetry LmLt (); LmLt_1 (); LmLt_2 (). Interestingly, the requirement of anomaly cancellation in a model with gauge symmetry does not require any other fermion content apart from the usual standard model ones. Three additional fermions namely, are added with such choices of charges that do not introduce any anomalies. Out of these three fermions, is the light sterile neutrino of eV scale while the other two are heavy neutrinos. The scalar sector of the model also consists of three additional scalar fields apart from the standard model Higgs field . There also exists an approximate global symmetry required to keep the bare mass term of sterile neutrino absent from the Lagrangian. The Yukawa Lagrangian involving the leptonic fields can be written as
(14) 
The relevant part of the scalar potential can be written as
(15) 
Denoting the vacuum expectation values (vev) of the neutral components of the scalar fields as , we can derive the leptonic mass matrices. The charged lepton mass matrix is diagonal and takes the form . The neutral fermion mass matrix in the basis can be written as
(16) 
Here is the Dirac neutrino mass matrix written in basis as
(17) 
The other two matrices are given as
(18) 
In the case where , the effective light neutrino mass matrix in the basis can be written as sterileearlier2 ()
(19) 
Using the above definitions of , the light neutrino mass matrix is