Normal hierarchy neutrino mass model revisited with leptogenesis
Abstract
We have studied the scenario of baryogenesis via leptogenesis in an flavor symmetric framework considering type I seesaw as the origin of neutrino mass. Because of the presence of the fifth generation right handed neutrino the model naturally generates nonzero reactor mixing angle. We have considered two vev alignments for the extra flavon and studied the consequences in detail. As a whole the additional flavon along with the extra right handed neutrinos allow us to study thermal leptogenesis by the decay of the lightest right handed neutrino present in the model. We have computed the matterantimatter asymmetry for both flavor dependent and flavor independent leptogenesis by considering a considerably wider range of right handed neutrino mass. Finally, we correlate the baryon asymmetry of the universe (BAU) with the model parameters and light neutrino masses.
pacs:
14.60.Pq, 11.30.QcI Introduction
The Standard Model(SM)with its enormous success stands at the top after the milestone discovery of a 125 GeV neutral scalar boson, the Higgs Boson. Still, there remained some unanswered phenomena in this ballpark of particles and forces like matterantimatter asymmetry, origin of neutrino mass and the observed dark matter, which keeps the window open for physics beyond Standard Model. It will be compelling if we can bring all of these occurrences into a single frame. With this motivation we put forward a theory which addresses the neutrino phenomenology consistent with recent global fit oscillation data along with a possible origin of matterantimatter asymmetry.
Oscillation parameters  (NO)  (IO) 

7.05  8.14  7.05  8.14  
2.43  2.67  2.372.61  
0.273  0.379  0.273  0.379  
0.384  0.635  0.388  0.638  
0.0189  0.0239  0.0189  0.0239  
0  2.00  0.000.17 & 0.79  2.00 
The fact of existence of neutrino mass has been steadfastly established by the dedicated neutrino oscillation experiments Fukuda:1998mi (); Ahmad:2002jz (); arthur2016 (); kajita2016 (). The status of neutrino oscillation data is presented in the table 1. The continuous rumble of various seesaw mechanisms seesaw () in the matter of explaining neutrino mass and their smallness in comparison to other fermions in the Standard Model is magnificent. Because of the absence of right handed neutrinos in the SM, neutrino mass is not explainable within this paradigm. Thus one has to go beyond Standard Model(BSM) by extending the fermion sector of the SM including two or more right handed neutrinos in order to implement see saw mechanism to make the theory viable for explaining neutrino mass and mixing. With robust number of evidences it has now become a proven fact that there exists tiny excess of matter over antimatter which is known as the Baryon asymmetry of the Universe(BAU). The dynamical process of production of baryon asymmetry from baryon symmetric era is familiar as baryogenesis. Although there are huge evidences that suggest the tiny excess of matter over antimatter that was produced in the early universe still its origin remains illusive. With this growing evidence there are several ways through which baryon asymmetry is realized. Among them baryogenesis via leptogenesis is considered as one of the most field theoretically consistent ways of explaining baryogenesis as proposed by Fukugita and YanagidaFukugita:1986 (). While realizing this picture it has come to notice that the L violating out of equilibrium decays of singlet neutrinos with mass larger than the critical temperature create an initial excess of lepton number L(for detail one may refer to leptogenesis ()). This excess in lepton number is then partially converted into the baryon asymmetry of the universe(BAU) via (B+L) violating sphaleron transition sphaleron (); Pilaftsis:1998pd (); Flanz:1998kr ().
In the same context one can traditionally search for the possibility of foreseeing the lepton asymmetry generated in the lepton sector due to the presence of heavy right handed neutrinos(RHN) as the key ingredients for type I seesaw to take place. In this regard lot of works Borah:2017qdu (); Borah:2014bda (); Karmakar:2014dva () have been exercised for a common search, addressing these two problems within a single frame. In order to bring this scenario into picture, SMsinglet heavy RH neutrinos are introduced, which through a dimension five operator eventually gives rise to tiny Majorana neutrino masses. The lepton asymmetry is dynamically generated by the L violating out of equilibrium decay of the lightest RHN satisfying Sakharov’s conditions (for detail please see Sakharov (); Davidson:2008bu ()) required for a nonzero baryon asymmetry of the universe.
If we look at the Lagrangian for the leptons we see that it permits the leptonnumberviolating decays of (for i = 1, 2) via: and . Since each decay mode can take place at both tree and oneloop levels, the interference of two decay amplitudes contributes to a CPviolating asymmetry between and its CPconjugated process . If Sakharov’s third condition is satisfied the outofequilibrium decays of the lightest RHN , may result in a net lepton number asymmetry which later on may convert into the observed predominance of matter over antimatter. Such an elegant baryogenesisvialeptogenesis mechanism provides a viable interpretation of the cosmological baryon number asymmetry, which is a ratio of the difference in number densities of baryons () and anti baryons () to the entropy density of the universe, , which has recently been reported by Planck 2015 Ade:2015xua ().
The novel Yukawacoupling texture of our model leads to the normal neutrino mass hierarchy with and a broken tribimaximal neutrino mixing pattern, as a result of which there is an automatic generation of nonvanishing reactor mixing angle. It is interesting to note that the model under consideration can accommodate non zero reactor mixing angle along with the explanation for non zero BAU without the help of any separate perturbation. From the necessity of seesaw to work a number of right handed neutrinos are introduced which are Majorana by nature. These extra right handed neutrinos can address by modifying the light neutrino mass matrix in such a way and at the same time offer explanation for the tiny excess of matter over antimatter via the process of leptogenesis. The complex Dirac Yukawa couplings of this system give rise to a nonzero lepton asymmetry which in turn yields the observed BAU. We have studied the lepton asymmetry generated over a range of right handed neutrino mass thereafter exploring the possibility of having both flavored and unflavored leptogenesis. Several BSM frameworks (Borah:2013bza, ; Borah:2015vra, ; Kalita:2014vxa, ; Kalita:2014mga, ) are available in the literature where baryon asymmetry is produced by thermal leptogenesis with hierarchical righthanded neutrinos. These heavy Majorana neutrinos decay to SM particles violating lepton number, which later on gets converted into baryon number by nonperturbative sphaleron interaction. The goal of this paper is to study the cosmological baryon number asymmetry produced as a consequence of the presence of the fourth() and fifth generation of RHN (), the presence of which also introduces a nonvanishing reactor mixing angle in the theory. Thus, we aim here at presenting the BAU in terms of the set of light neutrino model parameters that gives rise to correct neutrino data. In this piece of work we have kept our analysis only upto finding the matterantimatter asymmetry produced via the mechanism of leptogenesis. Two different kinds of vev alignments for the extra flavon are chosen in order to thoroughly study the affect of the same on observed neutrino parameters and BAU. Although the two different vev alignments lead to two different kinds of Dirac neutrino mass structures, the final light neutrino mass matrix remains the same. This fact permits us to keep the study for neutrino phenomenology as same as it is there in the original work by Meloni:2010sk (). However these two different Dirac mass matrices bring a little modification in lepton asymmetry calculation, which we discuss in numerical analysis section.
This work is solely dedicated to explore the possibility of foreseeing baryogenesis via leptogenesis through the realization of a broken symmetric mass matrix. Rest of the paper has been planned in the following manner, in Section II we present the model. Section III is discussed with type I seesaw and Leptogenesis. Section IV is kept for numerical analysis and results. Finally in Section V we end up with conclusion.
Ii The Model
Among the variant nonAbelian discrete flavor symmetry groups stands out as the most appealing group in the context of understanding neutrino mass and mixing properties. It has shown a promising role in explaining the origin of tribimaximal type of neutrino mixing since long. Keeping this in mind, this symmetry group has been chosen to explain neutrino mass and mixing. Although merely needs an introduction, still we slightly describe some properties of this group and how have they been utilized to give structures to the Dirac and Majorana mass matrices of the model. The nonAbelian group is the first alternating group and is isomorphic to the tetrahedral group . has four irreducible representations, among them there are three singlets and one triplet 3. The group has 12 elements, which can be written in terms of the generators of the group S and T. Where the generators satisfy the following relation (for detail one may refer to discreteRev ())
We consider the model discussed in Meloni:2010sk () for the purpose of studying baryogenesis via leptogenesis through the CP violating decay of the lightest RHN present in the model. In this model we have a total of five right handed neutrinos, among them three are the components of which transform as an triplet and hence are degenerate. The other two viz., and transform as singlets and respectively. Therefore the SM fermion sector has been extended by the inclusion of three fermion singlets. At the same time there is an extra flavon which is kept as an doublet and triplet. The full particle content of the model has been shown in table 2.
2  2  2  1  1  1  1  1  1  2  2  
1  1  3  1  3 
The Yukawa Lagrangian for the neutrino sector can be written as
Following the product rules as mentioned in Appen.A we can arrive at the following structures for Dirac mass matrix by considering two different vev alignments e.g. when takes vev as and respectively .
(1) 
In the same way we obtain the following structure for Majorana neutrino mass,
(2) 
Irrespective of the two different Dirac mass matrix the complete light neutrino mass matrix structure remains almost the same. Therefore, neutrino parameters continue to be similar for the later case. Since the lepton asymmetry parameter depends on the complex Dirac Yukawa couplings as seen from the Eq. (III.1) we study the variation brought out by these two structures by considering two different vev alignments of the the doublet field .
Iii Type I seesaw mechanism and its consequence Leptogenesis
Having set the stage we can write the effective light neutrino mass matrix from the type I seesaw realization.
(3) 
The Light neutrino mass matrices we obtain with the help of Eq.(1) and Eq.(2) are as follows
(4) 
where, the elements are defined as . One of the consequences of the type I seesaw is the lepton asymmetry produced by the decay of the lightest RHN present in the system.
iii.1 Leptogenesis
Since as a crucial requirement of fulfilling the seesaw mechanism, RHNs are already present in the model, this fact gives us the opportunity to study leptogenesis scenario in the model under consideration. Lepton asymmetry is created by the decay of the lightest RHN present in the model. It is to note that all the Dirac Yukawa couplings coming from the type I seesaw are complex and hence can act as a source of CPviolation, as there are no CPviolating phase associated with RHNs. Therefore it can be said that the lone contribution to CPasymmetry is coming from the complex Dirac Yukawa couplings in the present context.
with GeV, the vev of the SM higgs doublet.
with . When we sum over all the flavors popularly which is known as one flavor or unflavored leptogenesis, the second term in the expression vanishes. The sum over all flavors then can be written as
(5) 
From the expression above for lepton asymmetry, one can write the final BAU as,
(6) 
where determines the fraction of lepton asymmetry being converted into baryon asymmetry, the value of is found to be 0.55. is the dilution factor responsible for the wash out processes which erase the generated asymmetry. One can have following expressions for (for detail you may see sphaleron (); Pilaftsis:1998pd (); Flanz:1998kr ())depending on the scale of the wash out factor .
(7)  
(8)  
(9) 
where, quantifies the deviation of the decay rate of the lightest RHN from the expansion rate of the universe and is parametrized as,
(10) 
where, is the decay width of the lightest decaying RHN and H is the Hubble rate of expansion at temperature . We denote the effective number of relativistic degrees of freedom as and is approximately 110. It is worth mentioning that only the decay of the lightest RHN serves effective for a given hierarchical RHN spectra, as the asymmetry produced by the decay of heavier RHNs gets washed out notably. Therefore, is the only germane quantity which contributes to the total lepton asymmetry. There are three regimes of baryogenesis depending on the scale of the decaying right handed neutrino mass. The expression for lepton asymmetry as given in the Eq.(III.1) implies the asymmetry generated by summing over all the flavors. The expression given in the Eq.(5) works in the regime when the decaying RHN mass goes greater than or equal to GeV, when all the flavors behave similarly and are out of equilibrium. Now if the mass of the decaying RHN falls with the range GeV then, the flavor is in equilibrium and distinguishable and the scenario is familiar as tauflavor or two flavored leptogenesis. In this domain there are two relevant lepton asymmetry parameters and . One can write the final BAU for this regime as given by in Eq.(11). Again if the RHN mass falls less than GeV , it is termed as three flavored or fully flavored leptogenesis, where all the flavors are in equilibrium and hence distinguishable. For this mass scale of the RHN the final BAU is estimated by an expression given in Eq.(12).
(11) 
(12) 
where, . One can write the expression for as
For numerical analysis we have diagonalised the right handed neutrino mass matrix for getting the eigenvalues of that we show in numerical analysis section. The Dirac mass matrix has been chosen in a basis where RH neutrino mass matrix is diagonal. For that one can write , where .
Iv Numerical Analysis and results
The light neutrino mass matrix that we obtain with the help of type I seesaw formula using Eq.(1) and Eq.(2) yields a symmetry broken structure which automatically takse nonzero reactor mixing angle into account. The two different vev alignments gives rise to two different Dirac mass matrix and hence two different light neutrino mass matrices. This fact allows us to study leptogenesis for two Dirac Yukawa coupling matrices. From Eq.(4) it is clear that however we have two seperate to feed into the type I seesaw formula, the final light neutrino mass matrices obtained for each case are almost similar and yields the neutrino phenomenology in relatively same manner. The two Dirac mass matrices(, ) are different from each other in a sense that, has more number of nonzero entries in compaison to , which brings modification in the leptogenesis calculation as the CP asymmetry is dependent on the Dirac Yukawa coupling as seen from the Eq.(III.1). Apart from studying leptogenesis, the first kind of vev alignment of the flavon allow us to declare it as a potential dark matter candidate, as shown in Meloni:2010sk (); Hirsch:2010ru (); Mukherjee:2015axj (). It is interesting to note that the presence of the 5th RHN results into the light neutrino mass matrix structure in a way that the symmetry gets broken yielding with a nonzero value which falls within the range. Authors in Meloni:2010sk () has explained the reason of why only normal hierarchy is supported by the presented model. Thus we have calculated the lepton asymmetry only for the parameter space where neutrino mass follow normal hierarchy. There are five eigenvalues for the Majorana mass matrix , among which three belong to the same triplet and hence exactly degenerate. Therefore diagonalization of gives rise to two distinct hierarchical eigenvalues of the RH neutrino mass matrix. Among them we may chose any one to be slightly lighter than the other, CPviolating decay of which to SM leptons and Higgs is supposed to create lepton asymmetry.
Since both kinds of vev alignment leads to almost similar kind of light neutrino mass matrix, we have chosen the same way for calculating the light neutrino model parameters for each case( and ). It is known that a complex symmetric matrix has got 12 independent real parameters, among which there are three unphysical phases after readsorbing whom the number of parameter comes down to nine. In this work the light neutrino mass matrix Eq. (4) has four complex parameters and hence five real independent parameters eg., and the modulous and phase of the combination as derived in Meloni:2010sk (). With the new phase relationship the light neutrino mass matrix one can write in the following form
(13) 
Now we can make three equations in order to relate three parameters e.g. a, b and c with the light neutrino masses using the following set of equations.
(14)  
where, we denote “Tr” and “det” as trace and determinant of the light neutrino mass matrix. We did a random scan of some of the independent parameters in order to find a, b and c. From the above relations we see that a, b and c can be written in terms of the light neutrino mass eigenvalues. Since this analysis is only restricted to normal hierarchy mass pattern thus we can have as the three light neutrino masses. We choose as the independent parameters and are randomly varied within their range as reported by deSalas:2017kay () for numerical computation. On the other hand among the model parameters and have been chosen for random scan within the interval mentioned below.
(15)  
As already mentioned that for the two vev alignments the light neutrino mass matrices are almost similar thus we determine the model parameters with the help of the above mentioned equations. With the values of a, b, c found from the solutions we construct the Dirac Yukawa coupling matrices for each case separately. Once the two different s are constructed, with the help of them lepton asymmetry parameter is calculated for each type i.e., and . To find the Dirac Yukawa couplings for each case following definitions are used during calculation.
For numerical convenience we take which results into . Using this approximation we evaluate . Using this set of solutions we construct for each case.
(16) 
As already mentioned in Sec.III there are three regimes of leptogenesis, one flavor, two flavor and three flavor depending on the RHN mass scale. Now we may investigate whether the presented framework has some parameter space for each of the three regimes or not. For that, we choose the mass scales of RHN to be of three ranges. We choose the parameter g to be slightly larger than f in the majorana mass matrix by denoting f as and g as . For one flavour regime we have chosen GeV and GeV which is the specified mass domain for unflavored leptogenesis as mentioned in Samanta:2017kce (). For tauflavor leptogenesis the RHN mass has been chosen around GeV. As demanded by the corresponding mass regime for RHN, have been chosen less than GeV. We have calculated the lepton asymmetry by changing the ratio of RHN mass squared for each regime of leptogenesis. We choose the RHN masses as required for different region of leptogenesis as shown in table 3.
One flavor  1.5625  
Two flavor  1.0018  
Three flavor  1.00002 
To start with, we first evaluate the model parameters a, b and c by randomly varying the mass splittings in their allowed ranges along with the lightest mass as mentioned in Eq.(15). At the same time we have also varied and in their above mentioned specified ranges. With the help of the values found for a, b and c we found out the Dirac Yukawa couplings by considering three different ranges for heavy RHNs. With the help of the Dirac Yukawa couplings the Dirac mass matrix is constructed and with the help of that the lepton asymmetry for the leptonic sector is calculated. In order to find BAU using that value for lepton asymmetry. We have restricted our analysis for leptogenesis only upto normal hierarchy pattern of light neutrino mass as the model disfavor the IH mass pattern.
iv.1 VEV1 case
As already mentioned in the previous section we have computed the lepton asymmetry parameter for two type of Dirac mass matrices. Using that value of baryon asymmetry is determined. We categorize the corresponding results for each Dirac mass matrices and presented them in separate subsections. In this subsection we show the results for the which results for the VEV alignment of as . We present here the BAU as a function of the model parameters. In Fig 2 we plot for BAU as a function of the model parameters for the case of unflavored leptogenesis. BAU for the two flavor regime has been shown in Fig 3. In this plot we have shown the BAU for two different values of the Majorana mass splitting , one for and another for . And from the two plots it is clear that, as the mass squared ratio decreases we can have more parameter space for baryon asymmetry of the universe for the chosen range of RHN mass. Fig 4 evinces the BAU as a function of the model parameters for three flavored leptogenesis i.e., when all the flavors are out of equilibrium. Here also we have kept two values of and calculated the corresponding BAU for the chosen range of right handed neutrino mass pertinent to three flavored region of leptogenesis. From two and three flavor leptogenesis calculation one significant point is to be noted, that is as we go far below GeV for the concerned RHN mass, the mass squared ratio shall go on decreasing to produce the required amount of CP asymmetry which accounts for the observed BAU. For unflavored leptogenesis scenario the ratio can be bigger compared to that required for flavored leptogenesis as the RHN mass falls larger than or equal to GeV which is quite high. But as we go towards the lower mass regime for RHN the splitting demands a smaller value in order to give rise to adequate CP asymmetry to yield the observed baryon asymmetry. Fig 5 presents the parameter space for and . If we look for a common parameter space for and , a conclusion can be drawn from Fig 5 that, the model prediction for the same is much better for two and three flavor leptogenesis scenarios, while in the case of unflavored leptogenesis the plank bound for the observed BAU does not meet the recent bound for nonzero reactor mixing angle.
iv.2 VEV2 case
We repeat the same procedure for determining BAU considering the Dirac mass matrix of the second kind. While computing lepton asymmetry for this case it is seen that unlike the previous case (with Dirac mass matrix of first kind), the Yukawa coupling associated with the fifth generation of RHN() takes active part in producing a sufficient amount of lepton asymmetry required to generate the expected BAU. In Fig 6 we present the BAU found with respect to the model parameters for unflavored leptogenesis scenario. For flavored leptogenesis we have chosen two possible mass domains of the decyaing RHN, one for two flavor regime and another for three flavor regime. Despite of this fact we need to chose two RHN masses( and ) for each type of leptogenesis scenario(one, two and three flavor)in order to keep the mass squared ratio ()needed for successful leptogenesis. Fig 7 represents BAU as a function of model parameters and nonzero reactor mixing angle. In case of two/three flavored leptogenesis like the earlier case, for / we get the required lepton asymmetry which later on gets converted into baryon asymmetry by electroweak sphaleron. This mass splitting plays a vital role for having a reasonable amount of lepton asymmetry which accounts for successful baryogenesis. For three flavored leptogenesis we present the variation BAU with respect to the model parameters in Fig 8.
V Conclusion
We have exercised a detailed analysis on baryogenesis via leptogenesis scenario considering a framework augmented with flavor symmetry where type I seesaw has been chosen as the mechanism of neutrino mass generation. Here we have considered two vev alignments for the extra flavon one kind of which allow us to accommodate as a stable dark matter candidate by respecting the symmetry. As already mentioned, here we have only shown the results of BAU with respect to the light neutrino model parameters since neutrino phenomenology has already been explored in the work Meloni:2010sk (). While searching for the allowed model parameter space, we have used the latest global fit neutrino oscillation data for two mass squared splittings in their allowed range and found relations among the various model parameters and the known light neutrino parameters. With the help of those light neutrino model parameters we find baryon asymmetry. In order to find the same we have chosen three ranges for the RHN mass as necessary for the corresponding flavor regimes. We have checked the viability of both flavored and unflavored leptogenesis by considering different possible mass regimes concerned with the right handed neutrinos. It is clear from the results that the model addressed leptogenesis at all the possible scales of RHN mass as mentioned in subsection A of Sec.III. Of special importance is the fact that the RHN mass squared ratio () for a particular mass regime plays a crucial role in achieving an adequate amount of lepton asymmetry which accounts for the observed matterantimatter asymmetry. We presented a table 3 summarizing the various RHN mass scales and the mass squared ratios which are demanded by the model under consideration in order to explain matterantimatter asymmetry. Some important observation can be made from the results and analysis:

For all the flavor regimes e.g. one, two and three the model prediction meets the observed BAU.

To account for a required amount of lepton asymmetry for one flavor regime the mass splitting can be a bit larger as compared to that required for two flavor and/or three flavor.

With the same mass splitting as kept for one flavor regime, it is difficult to reproduce an adequate amount of lepton asymmetry for two flavor regime. As we go below GeV for RHN mass which is a criteria for having two/three flavor leptogenesis the factor needs to be a little less to have an enhanced lepton asymmetry in order to generate the observed BAU.

Then for fully flavored regime or three flavor leptogenesis we choose the RHN mass scale to be around GeV. With this RHN mass, the mass squared ratio() requires much more smaller value to have a sufficient amount of lepton asymmetry to meet the observed BAU as reported by Plank 2015 data.

The common parameter space for nonzero reactor mixing angle as envisaged by the model and baryon asymmetry is difficult to coincide in case of one flavor leptogenesis, whereas in case of two and three flavored leptogenesis, the plot between baryon asymmetry versus says that their parameter space matches only for a very narrow region. This fact allows us to have more predictability of the RHN mass scale relevant to that particular range required for two and/or three flavored leptogenesis.
Notwithstanding the prediction for experimentally observed neutrino parameters, this model beautifully sheds light on one of the long standing puzzle of particle physics and cosmology, the baryon asymmetry of the universe. The model also touches the dark sector accommodating a stable dark matter candidate, the detailed exercise of which we keep for our next draft.
Appendix A product rules
is isomorphic to the symmetry group of a tetrahedron. It is a discrete nonAbelian group of even permutations of four objects. It has four irreducible representations: three onedimensional and one three dimensional which are denoted by and respectively. Their product rules are given as
. Presenting two triplets as and respectively, their direct product can be written as
References
 Y. Fukuda et al. [SuperKamiokande Collaboration], Phys. Rev. Lett. 81, 1562 (1998) [hepex/9807003].
 Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 89, 011301 (2002) [nuclex/0204008].
 Arthur B. McDonald, Rev. Mod. Phys. 88, 030502 (2016).
 Takaaki Kajita, Rev. Mod. Phys. 88, 030501 (2016).
 P. F. de Salas, D. V. Forero, C. A. Ternes, M. Tortola and J. W. F. Valle, [arXiv:1708.01186 [hepph]].
 P. Minkowski, Phys. Lett. B67, 421 (1977); R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett 44, 912 (1980); J. Schechter and J. W. F. Valle, Phys. Rev. D22, 2227 (1980).
 M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986).
 A. Riotto and M. Trodden, Ann. Rev. Nucl. Part. Sci. 49, 35 (1999); S. Davidson, E. Nardi and Y. Nir, Phys. Rept. 466, 105 (2008).
 E. W. Kolb and M. S. Turner, Front. Phys. 69, 1 (1990).
 A. Pilaftsis, Int. J. Mod. Phys. A 14, 1811 (1999) [hepph/9812256].
 M. Flanz and E. A. Paschos, Phys. Rev. D 58, 113009 (1998) [hepph/9805427].
 D. Borah, M. K. Das and A. Mukherjee, arXiv:1711.02445 [hepph].
 M. Borah, D. Borah, M. K. Das and S. Patra, Phys. Rev. D 90, no. 9, 095020 (2014) [arXiv:1408.3191 [hepph]].
 B. Karmakar and A. Sil, Phys. Rev. D 91, 013004 (2015) [arXiv:1407.5826 [hepph]]
 A. D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5, 32 (1967).
 S. Davidson, E. Nardi and Y. Nir, Phys. Rept. 466, 105 (2008) [arXiv:0802.2962 [hepph]].
 P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 594, A13 (2016) [arXiv:1502.01589 [astroph.CO]].
 D. Borah and M. K. Das, Phys. Rev. D 90, no. 1, 015006 (2014) [arXiv:1303.1758 [hepph]].
 M. Borah, D. Borah and M. K. Das, Phys. Rev. D 91, 113008 (2015) [arXiv:1503.03431 [hepph]].
 R. Kalita, D. Borah and M. K. Das, Nucl. Phys. B 894, 307 (2015) [arXiv:1412.8333 [hepph]].
 R. Kalita and D. Borah, Int. J. Mod. Phys. A 30, no. 09, 1550045 (2015) [arXiv:1410.8437 [hepph]].
 D. Meloni, S. Morisi and E. Peinado, Phys. Lett. B 697, 339 (2011) [arXiv:1011.1371 [hepph]].
 H. Ishimori, T. Kobayashi, H. Ohki, Y. Shimizu, H. Okada and M. Tanimoto, Prog. Theor. Phys. Suppl. 183, 1 (2010) [arXiv:1003.3552 [hepth]].
 M. Hirsch, S. Morisi, E. Peinado and J. W. F. Valle, Phys. Rev. D 82, 116003 (2010) [arXiv:1007.0871 [hepph]].
 A. Mukherjee and M. K. Das, Nucl. Phys. B 913, 643 (2016) [arXiv:1512.02384 [hepph]].
 R. Samanta, P. Roy and A. Ghosal, arXiv:1712.06555 [hepph].