Neutrino decoherence from generalised uncertainty

Quantum gravity models predict a minimal measurable length which gives rise to a modification in the uncertainty principle. One of the simplest manifestations of these generalised uncertainty principles is the linear quadratic generalised uncertainty principle which leads to a modified Heisenberg algebra. This can alter the usual von-Neumann evolution of density matrix to a Lindblad-type equation. We show how this can give rise to a decoherence in neutrino propagation in vacuum. The decoherence effects due to the linear quadratic generalised uncertainty principle are extremely minimal and is unlikely to be detectable in the existing or upcoming experimental facilities for any of the natural sources of neutrinos. We also show that, in principle, there can be other variants of generalised uncertainty principle which predicts verifiable decoherence effects for the cosmic neutrino background.


Introduction
The discovery and subsequent verification of the neutrino oscillation phenomenon ushered a new age in the neutrino physics [1][2][3][4][5].Current development in the experimental front has already brought neutrino physics into the precision era and thus opened the possibilities to explore physics beyond the standard model (BSM).In fact, the very observation of oscillation phenomenon takes away the masslessness nature of neutrinos, necessitating BSM physics.The mixing between neutrino mass and flavour eigenstates is at the root of the neutrino oscillation.The non-diagonal nature of the chargedcurrent weak interaction, in the mass eigenbasis for both the charged leptons and the neutrinos, is the underlying reason for the aforementioned mixing.More precisely, in the diagonal basis of the charged leptons, the neutrino flavour eigenstates (ν α , α = e, µ, τ ) are a linear superposition of the mass eigenstates (ν i , i = 1, 2, 3), and they are related as ν α = U αi ν i , where U αi are the elements of the unitary lepton mixing matrix popularly known as Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix.Now, due to the feeble interactions with matter neutrinos can travel enormous distances in isolation from its environment.This makes neutrino oscillation effectively coherent, i.e., the wavefunctions of two neutrinos of identical energy travelling along an identical path evolve identically.However, this coherence can be lost over large distances if the neutrinos interact with the environment.Note that this is different from the wave-packet decoherence that occurs since the packets of the individual propagation eigenstates would get separated in space over long distances [6][7][8].Therefore we can dub the former scenario as environmental decoherence which is what we are going to focus in this work.The usual matter effects e.g., Mikheev-Smirnov-Wolfenstein (MSW) effect [9,10] and parametric resonance [11] etc. are the result from the influence of the multitude of matter particles on propagating neutrinos and therefore these effects consistently work on all traversing neutrinos maintaining the coherence.However, the environmental decoherence can, in principle, alter the known matter effect and can have implications in experimental observations [12,13].
Among other things a stochastic environment can be a source of environmental decoherence for neutrinos.In general quantum gravity models predict fluctuating space-time at Planck scales.This space-time "foam" results in a stochastic medium which can affect the propagation of the particles and effectively result in a decoherence in the neutrinos [14][15][16][17][18][19].To correctly demonstrate the effect of decoherence we need to go beyond the standard quantum mechanical description of neutrino oscillation since due to the system-environment interactions an initial pure quantum state evolves to a mixed state and the ideas from open quantum systems have to be employed.Quantum decoherence has been a subject of interest since the 1970s [20].It was proposed as an explanation of the transition from quantum to classical in case of macroscopic behaviour of a system.It explains how a purely coherent quantum system interacts with the environment and loses its characteristic quantum feature, i.e., the wavefuncions' ability to interfere with themselves.In this process the population of a quantum system changes due to the entanglement with the environment which can be a thermal bath or due to some other phenomenon [21].The evolution of the density matrix of a state due to decoherence is governed by Gorini-Kossakowski-Sudarshan-Lindblad master equation instead of the von-Neumann equation [22,23].Using the framework of open quantum systems neutrino decoherence has been extensively studied from phenomenological and experimental point of view [14,15,18,19,[24][25][26][27][28][29][30][31][32][33][34][35].
In this paper we explore the possibility neutrino decoherence that can arise from a different manifestation of quantum gravity ideas.The most promising candidates of quantum gravity, namely, string theory, loop quantum gravity, modified dispersion relations and doubly special relativity etc. all predicts a minimum measurable length approximately of the order of the Planck length l p ∼ 10 −35 m.A constraint of this kind modifies the Heisenberg uncertainty principle to generalised uncertainty principle (GUP) [36][37][38][39][40][41][42][43][44].Some of the interesting phenomenological and theoretical implications of GUP can be found in [45][46][47][48][49][50][51][52][53][54].The manifestation of the generalisation is contained in the deformed canonical commutation relations and they in turn gives rise to modified equations for the evolution of a quantum state which are Lindblad-type master equations [55,56] which will determine the time evolution of the neutrino state and ultimately give rise to the decoherence.The effect of decoherence is in general parametrised by the energy-dependent decoherence parameter Γ. Due to the ambiguity in the theory of quantum gravity sometimes an empirical approach is taken and a power law energy dependence of Γ is assumed [26,29,30], In some studies, experimental observations have been used to determine the validity of different values of n in a model independent way [14,15,18,19].In this study, we take a model-dependent approach to use the Lindblad-type master equation arising from a GUP to investigate the decoherence effect on the neutrino sector.
The paper is organised as follows: in section 2 we take a short recap on the conventional treatment of the neutrino oscillation, in section 3 we discuss the generalised uncertainty principle and the corresponding deformed canonical commutation relation used in this study.In section 4 decoherence effects on neutrino due to this model has been explained.Possible extensions of the current scenario are discussed in section 5. We explain the results in section 6 and finally we summarize and conclude.

Conventional Treatment of Neutrino Oscillation
In this section we briefly describe the conventional treatment of neutrino oscillation to set our notations and conventions.
In general, neutrino flavour eigenstates can be expressed in terms of the neutrino mass eigenstates with the PMNS mixing matrix as, where n is the number of neutrino mass eigenstates and U is the PMNS matrix.Clearly, for the two neutrino case, n = 2 and α = e, µ, and for the three neutrino case n = 3 and α = e, µ, τ .In general the neutrino oscillation probabilities can be given by the well-known expression [57], where α, β = e, µ, τ, . .., L is the neutrino path length, E is the neutrino energy, ∆m 2 ab = m 2 a − m 2 b , a, b = 1, 2, 3, . . . is the mass squared difference between the neutrino mass eigenstates |ν a and |ν b with masses m a and m b respectively.The oscillation probability for α = β is termed as transition probability and the α = β case is known as survival probability.In the two-flavour case, the mixing matrix U takes the form However, in the realistic three-flavour case the mixing matrix U is the famous PMNS matrix which is given by [58], where the c ij = cos θ ij and the s ij = sin θ ij are used to denote the cosine and sine of the mixing angles θ ij respectively and δ is the CP-violating phase.Conventionally, they are taken to be in the ranges Here we do not mention the possibility of the presence of Majorana phases in this 3 × 3 mixing matrix, since in the oscillation phenomena they play no roles.For our numerical analysis we will use the experimentally obtained best-fit values from the NuFIT collaboration [59], for the normal ordering and for the inverted ordering.
Moreover, for our purposes we will consider just the central values of these parameter.
In passing we note that the transition probability is sometimes taken as Gaussian average.Actually, due to the lack of sharp energy or well-defined propagation length, a Gaussian average over the L/E dependence and other uncertainties are taken.The averaged transition probability takes the form, where σ is the damping parameter which is responsible for the damping of the neutrino oscillation probabilities.It is to be noted that all the variables in the above equation e.g., P αβ , L, E are now averages of the transitions probability, neutrino propagation length and the neutrino energy respectively.The equivalence between this kind of averaging and the neutrino decoherence was shown for the first time in [25].In the next sections we will show this in the present context too.

Generalised Uncertainty Principle and Lindblad Master Equation
In this section we discuss the generalised uncertainty principle (GUP) and set-up the Lindblad equation relevant for our studies.

GUP and DCCR
The Heisenberg's Uncertainty principle, ∆x∆p ≥ /2 is a direct consequence of the non-commutativity of the position and the momentum operators, i.e. [x i , p j ] = i δ ij .As the various theories of quantum gravity suggest, this commutation relation is modified, i.e. we get a deformed canonical commutation relation (DCCR) in order to accommodate the prediction of a minimum measurable length in such theories.The modification presents itself as higher order terms of the momentum observable in the right-hand-side of the commutator relation.An example of such a DCCR, which is motivated by string theory, can be [55], where β = β 0 l 2 P / 2 with l P being the Planck length and β 0 is some constant.DCCRs of this form can give rise to the GUP, A more general approach, supported by doubly special relativity, has been investigated to construct a GUP which contains both linear and quadratic terms of momentum observable in the right hand side (LQGUP).This can be of the form [56], which is resulted in from the following DCCR, Here, the parameter α is defined as α = α 0 l P / .Here α is the GUP parameter and it is often treated as a Gaussian white noise for simplicity to comply with the existence of a fluctuating minimum length scale as predicted by the theories as mentioned before.The quantity α 0 is often taken to be the order of unity.

Generalisation to Lindblad equation
Using the LQGUP mentioned in the previous section, Eq. ( 15), a Lindbland-type master equation can be constructed which determines the evolution of the quantum system corresponding to the decoherence created by the GUP which in turn arises from the quantum gravity aspect of the theory.The DCCR gives rise to modified momentum which satisfies the Heisenberg's Uncertainty Principle [56].This means that taking cue from Eq. ( 16) we can write the operators x i and p j in the following form, where x 0i and p 0i satisfy the conventional relation [x 0i , p 0j ] = iδ ij .In addition the other commutation relations remain intact, i.e., [x i , x j ] = 0 = [p i , p j ].Note that from here onward we will use the natural unit = 1.The Lindbland-type master equation from the LQGUP depends on the Hamiltonian of the system.For example, in [56], the Hamiltonian was that of a non-relativistic particle.Here, since we will be considering neutrinos, we take the Hamiltonian of an ultra-relativistic free particle.The Schrödinger's equation in this case becomes, Using the modified relations in Eq. ( 17), we get Now we expand this and we take into account terms upto first order in α to get, where In this case, the Lindblad-type equation can be derived to be, where γ = t P l 2 P .Here the second term in the right hand side (RHS) of the equation is the dissipator or the Lindblad-type operator.We can express the above equation in a simpler way using the Lindblad In the next sections, we use this Lindblad-type master equation to investigate the decoherence in the neutrino flavour oscillation in vacuum.

Two flavour case
Neutrino flavour transitions through the decoherence mechanisms have been studied for Lindblad operators that are linearly dependent of the hamiltonian and they produce transition probabilities similar to that of Eq. ( 12) [25].In this case we use the Lindblads which are not linearly dependent on the Hamiltonian as shown in the Eqs.( 21) and (22).For the two flavour case, the effective Hamiltonian in the mass eigenstate basis can be written as, where ∆m 2 ≡ m 2 2 − m 2 1 , E is the energy of the neutrinos and ω 2 = ∆m 2 4E .In this case the Lindblad operators takes the form, where Here the index a represents different mass eigenstates.It is evident that since the D a terms are proportional to the identity matrix, there will be no decoherence effect since the second term in the RHS of Eq. ( 22) is vanishing.Nevertheless, below we show this explicitly and in the process lay down the steps of the oscillation probability calculation that will be used in the three flavour case as well.
The relevant density matrix ρ whose elements are functions of time and which follows the condition Tr(ρ) = 1 can be written in the following form, where a, p, q are real valued functions of time, i.e., a = a(t), p = p(t) and q = q(t).Inserting Eqs. ( 24) and ( 25) into Eq.( 22) we get the Lindblad-type master equation for the evolution of the density matrix in the two-flavour case as, where λ 2 = 2 a=1 λ 2 a .From the above equation, the solutions for the parameters can be obtained as, We can now take mass eigenstates to be, and following Eqs.( 2) and ( 4), the flavour eigenstates can be written as, for electron neutrino and muon neutrino, respectively.Now, if the system is initially in the electron neutrino state, then the initial density matrix is, From Eq. ( 25) the initial conditions for the parameters can be determined as, a(0) = cos2θ, p(0) = sin2θ and q(0) = 0. Therefore, the neutrino transition probabilities with the decoherence effects are, Using Eqs.(27) in Eq. (31b), we get transition probability as, Assuming t ∼ L in the above equation we get, As we can see from the above expression that for the two flavour case, the Lindblad-type equation does not give rise to any decoherence upto the first order of α in the Hamiltonian.This, as expected, is due to the fact that the Lindblad operators are proportional to the identity matrix.Hence, the dissipation term gets cancelled.

Three flavour case
In the three flavour case we can have two scenarios, namely the normal and inverted hierarchy between the eigenstates in the mass basis.In the normal mass hierarchy we assume, The Hamiltonian in the mass eigenbasis can be given by, In normal mass hierarchy, ∆m 2 21 , ∆m 2 31 > 0, Hence, the Lindblad operators take the form where ω 3 = 1/2E.Again, similar to Eq. ( 25) we can define the most general density matrix for the three flavour case as, Here, a, b, c, p, q, f, g, x, y are all real-valued functions of time.In order to maintain the condition Tr(ρ) = 1, the constraint on these variables is , Now, using Eqs.( 35)- (37) in Eq. ( 22), we get, Here, we can use t ∼ L in Eqs. ( 48) and ( 49) and can very easily identify that they follow the general form mentioned in Eq. (12).The decoherence effects manifest themselves in the exponentially decaying terms.The survival and transition probabilities for the other generations can also be obtained following the similar prescription.
For the inverted hierarchy case, the mass ordering is different than that of the normal mass hierarchy, i.e., m 2 2 > m 2 1 m 2 3 .In this case, the Hamiltonian takes the form Therefore, the Lindblad operators take the form, From Eqs. ( 36) and ( 51) it is evident that the Lindblad operators are the same for both normal and inverted mass hierarchy.Therefore, the decoherence effect in the survival and transition probabilities will be similar in normal and inverted hierarchies.

Possible extension of the model
Normally in the study of neutrino decoherence, the decoherence parameters, i.e. the argument of the exponentially decaying terms in the transition and survival probabilities are represented as power law of the energy, i.e., Γ ∝ E n [14,18,19].In case of the LQGUP the decoherence parameters, Γ ∝ E −4 as can be seen from Eqs. ( 48) and (49).Motivated by different quantum gravity scenarios there can be other forms of GUPs as well.In principle, they can give rise to different power law relations for the decoherence parameters.In this section, we show how this can be achieved and what modification will be there in the Heisenberg algebra for such a scenario.
For this section, we try to construct a Lindblad-type equation which generates a neutrino decoherence with decoherence parameters which are ∝ E −2 and then we try to generate the modified momentum relationship corresponding to the Heisenberg algebra.In Eq. ( 35), the Hamiltonian for the neutrino is given in the mass eigenstate.Therefore, in order to generate a n = −2 model for the deocherence, we need the Lindblad operators D a ∝ H. Hence, the form of the modified Lindblad-type equation in this case is, where A depends on the GUP parameter α as given in subsection 3.1.Now, in order to achieve such an equation using the treatments given in [55,56], the perturbation Hamiltonian should be of the form, where p 2 0 = 3 i=1 p 0i p 0i and p i and p 0i are related as, p i = p 0i (1 + α 1/6 + α 1/12 + B(α, p 0 )).Here, the function B(α, p 0 ) contains all other terms in α and p 0 .The particular form of B(α, p 0 ) depends on the GUP that generates this modified Heisenberg algebra.This can lead to the Lindblad-type Eq. ( 52), in which case, Similarly, a dependence of the form n = −3 can arise from a relation p i = p 0i (1+α 1/3 +α 1/6 +B(α, p 0 )).This relation can lead to the Lindblad-type equation, in which case, Here we show different power-law dependence of the decoherence parameters, which are usually studied in a model-independent manner [18], can in principle arise from different variants of GUP.Here, it is to be noted that the exact equations and the exact expression of the decoherence parameters depends on the entire form of the GUP.Here we have just tried to obtain the power law behaviour and the order of the decoherence parameters.

Results
Decoherence effects in neutrino oscillations due to quantum gravity have been studied from modelindependent perspective in [14,18,19].In this work, we have used a specific model where the effect of quantum gravity manifests in the LQGUP.In this section, we try to compare the consequences of decoherence due to the LQGUP, as well as possible extensions, to that of the existing modelindependent and other experimental searches of neutrino decoherence.

LQGUP Case
Decoherence parameter: In general model-independent studies, the energy-dependent decoherence parameter, Γ is expressed as a power law, where E 0 is the reference energy, which we take to be O(GeV), and n is the power.Apart from this we can see that unlike some phenomenological models [18], there are multiple decoherence parameters with different expressions depending on the mass-squared differences of the neutrino masses.In Figure (1a), the dependence of these decoherence parameters with respect to the energy are shown.It is to be noted that Γ 31 (E) and Γ 32 (E) are very close in the given range while Γ 21 (E) differs by a small amount.The underlying reason for this is the numerical values of the mass-squared differences which are ∆m 2  21 ∼ 10 −5 eV 2 , ∆m 2 31 ∼ 10 −3 eV 2 and ∆m 2 32 ∼ 10 −3 eV 2 .
Coherence length: The coherence length l coh is defined as the propagation length after which the amplitude of the survival and the transition probabilities of the neutrinos are diminished by a factor of e, i.e., Thus in our case, the coherence length increases rapidly with energy since l coh ∝ E 4 , as is evident from Eq. ( 58).This suggests that in this model, the effects of the quantum gravity becomes more prominent on lower energy regimes as compared to higher energy regimes of the neutrinos.
The coherence length can also be thought of as the interaction mean free path which is the average distance a neutrino covers before interacting with a virtual black hole (VBH) in the space-time foam [18].Therefore, it is likely that in this model, the high energy neutrinos are much less likely to interact with a VBH as compared to the low energy neutrinos as opposed to the models where n > 0 [18].It is also possible that an interaction with the VBH results in more noticeable consequences for low energy neutrinos in comparison with high energy neutrinos.An effect of this kind can also arise from some other model of quantum gravity where space-time uncertainties at very small length scales results in lightcone fluctuations.[60].It has been shown in [19] that lightcone fluctuations can also cause decoherence effects in neutrino oscillations in vacuum.An effect of that kind can lead to an energy dependence on the decoherence parameter ∝ E 2(n−1) where n denotes the dependence of the distance uncertainty on energy.Therefore, a model with n = −1 can lead to an energy dependence of the decoherence effects similar to the one obtained in this study through the LQGUP model.It is also of interest that the decoherence parameters dependent on different frequencies, i.e. different mass splittings (∆m 2  21 , ∆m 2 31 and ∆m 2  32 ) are a common feature in [19] as well.Although the exact form of the decoherence parameter depends on the other variables in the expression of distance fluctuation used in [19], it is a possibility that the lightcone fluctuations can be the underlying reason of the decoherence effects arising from the LQGUP.In Figure (1b) the dependence of decoherence parameters on the neutrino energies are shown.It is to be noted that l coh31 and l coh32 has similar values for the given energy range whereas l coh21 differs slightly which is again due to the hierarchy in the mass-squared differences.The length corresponding to the size of the observable Universe (∼ 10 28 m [61]) is shown by the horizontal dash-dotted line.This corresponds to the fact that according to this model, neutrinos with energy O(10 −23 ) GeV or lower has the coherence length less than that of the observable universe.This implies that the decoherence effect due to the GUP model is predominant for the lower energy neutrinos.

Possible extensions
In section 5, we have shown how for a different power-law the Lindblad-type equation can be modified through the modification in the Heisenberg algebra.The Eqs. ( 52) and ( 55) can now be used instead of Eq. ( 22) to calculate the survival and transition probabilities of the neutrinos.It is to be noted that for the case where n = −2, the Lindblad operators D a ∝ H and for the case where n = −3, the Lindblad operators D a ∝ H 3/2 .

Decoherence parameter:
The decoherence which will take place for the propagation of neutrinos through vacuum in this case will be governed by Eq. ( 52) and this leads to a decoherence parameters with n = −2 power-law, Similarly, for the n = −3 power law, the decoherence parameters are approximately of the form, In figures (2a) and (3a), we show the dependence of these decoherence parameters with respect to the energy.In comparison to the n = −4 case, the values of the decoherence parameters are much larger for the same energy range of the neutrinos.But similar to the previous scenario, this extension also suggests higher decoherence for lower energy neutrinos.Coherence length: The coherence length, as mentioned previously, l coh = 1/Γ(E).For the cases with n = −2 and n = −3 the coherence lengths are proportional to E 2 and E 3 respectively.In figure (2b) and (3b), we show the dependence of the coherence lengths on the energy of the neutrinos.In this case also, higher energy neutrinos will have longer coherence length than that of lower energy neutrinos, but this model predicts a shorter coherence length than the previous model for the same energy neutrinos.
In the scenario where n = −2, for lower energy neutrinos e.g., for solar neutrinos with energy O(1) MeV [57] the decoherence parameter Γ ∼ 10 −55 GeV and the coherence length l coh ∼ 10 39 m.However, for the cosmic neutrino background (CNB) neutrinos with energy O(10 −15 ) GeV [57] the decoherence parameter Γ ∼ 10 −31 GeV and the coherence length l coh ∼ 10 15 m.Clearly, this is in the paradigm of future experiments where finer details of the CNB can be measured.

Effects on free-streaming neutrinos
Neutrinos decoupled from the primordial plasma at T ∼ 1.3 MeV and they were of the free-streaming nature due to their ultra-relativistic nature [57].After decoupling, the temperature of these neutrinos is inversely proportional to the scale factor of the universe(T ν ∝ 1/a(t)).The decoupling occurred in the radiation dominated era.Since then, the universe went through matter-and cosmological constantdominated era respectively.In all these different eras, the scale factor, a(t) had different dependencies on time, i.e. during the radiation dominated era, a(t) ∝ t 1/2 , during the matter dominated era, a(t) ∝ t 2/3 and during the cosmological constant dominated era a(t) ∝ e Ht where H is the Hubble constant [57].Therefore, in all these different eras, the temperature and therefore the energy of these neutrinos will have different dependencies on the propagation length as the propagation length for neutrinos, L ∼ t.Therefore, for all practical purposes, during the radiation dominated era, the energy of the relic neutrinos, E ∝ L −1/2 , during the matter dominated era, E ∝ L −2/3 , and during the cosmological constant dominated era, E ∝ e −H 0 L where H 0 is the Hubble constant today.
According to the standard picture of neutrino oscillation, the flavour ratio of these relic neutrinos is 1 : 1 : 1 [62].We expect that in the presence of decoherence due to interactions with the vacuum, the relic neutrino flavour ratio might differ.To check this, we have used the dependencies of E on the propagation length L during different eras, mentioned in the previous paragraph, in the expressions of transition and survival probabilities, such as Eqs.( 48) and (49).We find that for the LQGUP (n = −4) case and the possible extension where n = −3, the flavour ratio today (t ∼ 14 Ga) remains 1 : 1 : 1 for the different flavour of neutrinos.For the case where n = −2, the flavour ratio changes by O(10 −16 ).In upcoming experiments, e.g., PTOLEMY, the expected number of detection of ν e is O(10) per year [62].Hence, it is not possible to detect a modification of O(10 −16 ) through an experiment with the working precision of that of PTOLEMY.However, the situation will not be so bleak in case of GUP models predicting n > 0, a detailed study of which we leave for a future work.

Summary and Conclusion
The LQGUP model is shown to give rise to a decoherence in the neutrino oscillation with decoherence parameter having a n = −4 power-law dependence on energy.Therefore, lower energy neutrinos will experience significantly more effect of the decoherence than the high energy neutrinos for a given propagation length.This strong quartic suppression of the decoherence effect with energy occurs because the Lindblad operators arising due to the LQGUP are proportional to the square of the Hamiltonian which in this case is inversely proportional to E. Through this approach, the mean free path of the neutrinos can be estimated which gives us the idea regarding the structure of the vacuum in the quantum gravity theories which predicts the LQGUP.Moreover, the decoherence of the neutrinos depends on the energy and therefore the vacuum interacts differently for different energy ranges of the neutrinos, e.g., for CNB neutrinos of energy O(10 −15 ) GeV the mean free path is O(10 63 ) m which suggests that a neutrino in this energy range will on average interact with the vacuum in way that diminishes its transition and survival probabilities by a factor of e after every O(10 63 ) m of propagation.One of the possible ways this vacuum interaction can occur is via VBH.Another possibility can be the lightcone fluctuations.As the energy goes down, the mean free path also reduces rapidly which suggests that lower energy neutrinos will experience more decoherence than higher energy neutrinos in LQGUP model.Since, the coherence length for even CNB neutrinos are very much larger than the size of the observable universe, the LQGUP predictions are impossible to be detected through existing experimental facilities for any natural source of neutrino.
In principle, one can use similar treatments to get other power-law dependencies from different variants of GUP.As examples we have shown the cases of n = −2 and n = −3.There can be different forms of GUPs which can give rise to modified momentum relations and consequently lead to a decoherence parameters with n = −2 and n = −3 power-law dependence on energy.In this case also, the decoherence effects are dominant in case of low energy neutrinos but unlike the LQGUP model, for the case of n = −2, the CNB neutrinos will have a decoherence length O(10 15 ) m which is very much realistic for the purposes of future detection.As far as cosmological observations are concerned, decoherence can, in principle, affect the free streaming properties of the relic neutrinos in the CNB.We find that though decoherence effects lead to no change in the flavour ratio of the CNB for the case of n = −4 and n = −3, but for n = −2 it will lead to a change of O(10 −16 ) in the flavour ratio of the CNB.But, due to the well known difficulties in direct detection of the CNB, this difference is too insignificant to be detected through the upcoming experiments in the near future.This result also serves as a motivation to study the different variants of GUP which originates from different theories of quantum gravity and perform a rigorous analysis to set the stage for a probe of quantum gravity from the neutrino oscillation experiment.
This study connects a quantum gravitational modification in the uncertainty principle to the energy dependence of the interaction of neutrinos with the vacuum.If a GUP model predicts that the Lindblad operators D a are proportional to H −n then that corresponds to the fact that the decoherence parameter Γ(E) is proportional to E 2n .This kind of correlation can be used to pick a particular form of GUP which supports the neutrino oscillation and decoherence due to the interaction with the vacuum from the available data.Such model-dependent work can be used in future to indirectly test the claims of different models of quantum gravity.Similar methods can be followed to investigated other GUP models as well.
For GUP models which gives rise to neutrino decoherence with n > 0, the flavour ratio of high energy neutrinos from different natural and artificial sources could be investigated in order to test those GUP models.Such model-dependent work could confirm, rule out or put stringent bounds on different GUP models and indirectly test the claims of the different quantum gravity models which lead to these GUPs.Similarly, decoherence effects on different sources of high energy neutrinos will be optimized for different models of GUP with positive values of n and one could indirectly test the quantum gravity models through the similar decoherence studies.A generalised GUP model can also be devised using the general power law energy dependence of the decoherence parameter which will provide key insights into the relationship between energy dependence of the decoherence parameter and the modifications to the Heisenberg uncertainty principle.

Figure 1 :
Figure 1: The dependence on energy for (a) the decoherence parameter and (b) the coherence length in the LQGUP scenario.

Figure 2 :
Figure 2: The dependence on energy for (a) the decoherence parameter and (b) the coherence length in the extended GUP scenario with n = −2.

Figure 3 :
Figure 3: The dependence on energy for (a) the decoherence parameter and (b) the coherence length in the extended GUP scenario with n = −3.
13 s 12 c 13 s 13 e −iδ −s 12 c 23 − c 12 s 23 s 13 e iδ c 12 c 23 − s 12 s 23 s 13 e iδ s 23 c 13 s 12 s 23 − c 12 c 23 s 13 e iδ −c 12 s 23 − s 12 c 23 s 13 e iδ c 23 c 13 + 2c 12 s 12 c 23 s 23 c 2 13 s 13 sin δ sin and l p are Planck time and length, respectively.From this, we can note that in this model, the power n = −4.