Dissipative effect in long baseline neutrino experiments

The propagation of neutrinos in long baselines experiments may be influenced by dissipation effects. Using the Lindblad master equation we evolve neutrinos taking into account these dissipative effects. The MSW and the dissipative effects may change the behavior of the probabilities. In this work, we show and explain how the behavior of the probabilities can change due to the decoherence and relaxation effects acting individually with the MSW effect. A new exotic peak appears in this case and we show the difference between the decoherence and relaxation effects in the appearance of this peak. We also adapt the usual approximate expression for survival and appearance probabilities with all possible decoherence effects. We suppose the baseline of DUNE and show how each of the decoherence parameters changes the probabilities analyzing the possible modification using a numeric and an analytic approach.


Introduction
In the near future there will be new long baseline experiments [1][2][3] to test the standard pattern of neutrino oscillation as never before. These experiments will have great sensitivity in order to determine the parameters that describe the standard oscillation pattern. Through them it will be possible investigate the open questions, like CP-violation, the mass hierarchy, and the octant problem [4][5][6][7][8][9][10][11][12][13], and maybe new phenomena.
As is well known, quantum mechanics explains how the neutrinos are able to change the flavors during their propagation. The success of the model opened the opportunity for testing other kinds of effects in oscillation experiments [14][15][16][17][18][19][20][21]. So, considering the precision level that the next generation of experiments will achieve, the limits on each usual neutrino oscillation parameter will be very stringent and the space for new physics may also be lessened. a e-mail: robertol@ifi.unicamp.br In this work, we will use the Lindblad master equation to evolve the neutrinos [22,23]. This equation is used when a physical system is considered open to interact with a quantum environment that is treated as a reservoir [24,25]. Due to the interaction between the subsystem of interest and the environment, the quantum behavior of the subsystem of interest may change considerably during its quantum evolution. In the present case the neutrinos are our subsystem of interest, while the current hypotheses of the dissipative sources are quantum foam or quantum gravity when the neutrino propagates in matter or vacuum [20,[26][27][28][29], and there is a phenomenological theory to matter fluctuation as dissipative source [30,31] including a microscopic model in this case, as we can find in Ref. [31]. The quantum evolution through the Lindblad master equation is non-unitary and it adds in the evolution the possibility of dissipative effects occurring like decoherence, relaxation, and others.
In neutrino oscillation, the decoherence effect has been studied more than the other dissipative effects, and in most of these investigations, the neutrino propagation was in vacuum [19,20,[32][33][34][35]. However, considering the next generation of long baseline experiments, the matter interaction will be important and we will study some aspects of this case.
We will also probe the relaxation effect, which is much less studied in the literature even when the propagation is in vacuum. This effect is not important for the terrestrial long baseline case and only the solar neutrinos are responsible for putting the most stringent bound on the relaxation effect [36]. However, in the case of the propagation in constant matter, we are able to verify in another way the difference between the decoherence and relaxation effects. So, we will start with the two-neutrino approximation and will point out the difference between the relaxation and decoherence effects when the MSW effects is present. In this case, the survival and appearance probabilities present new exotic peaks due to the presence of dissipative effect and MSW effect. This behavior will be analyzed and we will show how the differences of these effects become evident in neutrino propagation.
As the θ 13 value is not null, the three neutrino families will have to be taken into account in the new generation of long baseline experiments [1][2][3] as well as the MSW effect [37,38]. In this study, we will present a phenomenological approach where the most effective quantum dissipative operator will be defined to the case of the three-neutrino oscillations. Then we will evolve the neutrinos in an exact and analytical scheme considering the baseline of DUNE [1].
In the analytical case, there is a study [39] where the author considered the approximation m 2 12 << m 2 31 to obtain analytical expressions for the oscillation probabilities in terms of the effective mixing angles. We will adapt these results to the case where the decoherence effects are taken into account and long and short analytical expressions for the survival and appearance probabilities will be introduced. Besides, we will include a discussion as regards the range of validity of these new probabilities.
We will show the behavior of these two versions of the analytical probabilities in comparison with the exact approach. The exact and long analytical expression for the probabilities present the same new peak introduced due to the coupling of the decoherence and MSW effects, while the short analytical expression for the probability due to the decoherence effect is divergent. Besides, we will point out using the exact approach how each decoherence parameter changes the survival and appearance probabilities considering the baseline of DUNE.
Furthermore, considering a short version of the analytical probabilities for DUNE, we are able to understand how the behaviors of the survival and appearance probabilities are influenced by the decoherence effects, because the short expression is an acceptable approximation of the exact probabilities in the relevant energy range in DUNE.
Then we conclude with a simple study of CP-violation where we show how the CP-violation phenomena may be changed depending on the magnitude of the decoherence effects.

Quantum dissipators
We introduce a dissipative formalism that may be used in long baseline experiments of the next generation like DUNE [1]. In this context, we will consider the neutrino as an open quantum system and its propagation will be made using the Lindblad master equation [23][24][25]40].
This evolution equation adds many dissipative effects, the decoherence and relaxation effects being the most effective of them. The physical meaning of these two effects can be found in Ref. [36].
The Lindblad master equation is usually written as where ρ ν is the interesting subsystem state, D[•] comes from the partial trace over the environment states in the evolution equation of the global system [24] and it must have the dimension of energy, being defined as where V b are dissipative operators that arise from interaction between the subsystem of interest and the quantum reservoir. So, V b may describe the dissipative effect of the neutrino propagating in quantum gravity space-time [27][28][29]41,42] or in fluctuating matter [30,31]. Besides, V b operators act only on the N -dimensional ρ ν space. In the evolution equation (1), the first term on the right side evolves the quantum state like the Liouville equation, while the second term, which depends on the V b operators, renders the evolution non-unitary. The second term is also responsible for introducing many kinds of dissipative effects [22][23][24][25].
We assume that the standard oscillation Hamiltonian is the same for the neutrinos as the subsystem of interest [36]. For the next generation of long baseline experiments, like DUNE [1], the matter density is important and the Hamiltonian can be written as where the V cc is usually the matter potential and on the right side of the last equality, we are assuming it to be possible to write the Hamiltonian in an effective mass basis. For our aims, it is useful to expand Eq. (1) in the SU (3) and SU (2) basis matrices. Then each operator in Eq. (1) can be expanded like O μ = a μ F μ , where F μ is composed by an identity matrix and the Gell-Mann and Pauli matrices for three and two dimensions, respectively. Thus, the evolution equation can be rewritten as with D μ0 = D 0ν = 0 to keep probability conservation, and we change t → x as usual for the ultra-relativistic approximation. As the componentρ 0 (t) = 0 we did not include it in the equation above and then the indices in Eq. (4) now run from 1 to k 2 − 1, where k is the number of families. The solution ofρ 0 (t) = 0 is trivial and it is just ρ 0 (t) = 1/k. Besides, the component a 0 from the V b operators must be null.
The D kl matrix in Eq. (4) has many parameters, even in the two-neutrino approximation. However, if we impose [H S , V b ] = 0 D kl is written as for the three-neutrino approach and for the two-neutrino approximation. We relabeled by i j all the compositions of the coefficients a k obtained from expansion of the dissipative terms [19]. The dissipators in Eqs. (5) and (6) add decoherence effects in the propagation. In the case of the dissipator in Eq. (5), each i j parameter may have different values for each quantum interference term and the indices i and j are associated with the quantum decoherence between the families i and j. In two families this association is not necessary and the double indices can be replaced by only one index [20,32,43].
The constraint [H S , V b ] = 0 imposes energy conservation on the subsystem of interest, which, in this case, are the neutrinos. However, this constraint is very stringent and physically there is no guarantee at all that this must occur, because the neutrino energy can fluctuate once it is free to interact with the quantum reservoir. It is important to remember that the energy conservation is always satisfied by the global system [43].
We can define another dissipator to violate the constraint before. It is made replacing the null entries in the main diagonal by new parameters in the dissipators (5) and (6). These parameters describe another dissipation effect, which is called the relaxation effect.
As the decoherence and relaxation effects depend on the propagation distance [36,43], only solar and astrophysical neutrinos are able to put stringent bounds on relaxation effects. In general, the coherence terms are averaged out for neutrinos that come from these sources and hence, it is not possible to have any information as regards the decoherence effect using our model-independent formalism [36,37,44]. On the other hand, the decoherence effect is only one dissipative effect that terrestrial experiments can limit because the oscillation effect is still important in these cases [36,37,44].
So, we rewrite only the dissipator in Eq. (6), adding the relaxation effect. In this last case, we will show how the dissipation and MSW effects together change the usual behavior of the probability in a particular oscillation channel.
Under the condition [H, V b ] = 0, the dissipator considering the two-neutrino approach is written as In the particular case of the dissipators in Eqs. (5) and (6), the off-diagonal elements must be null due to the constraint [H, V b ] = 0. However, we can ignore all off-diagonal terms in (7) due to complete positivity [45]. As D kl must be positive definite, the diagonal parameters must be larger than the offdiagonal ones. Thus, the dissipator in Eq. (7) is enough to study dissipation effects since the off-diagonal parameters only exist if the main diagonal is filled with the decoherence and the relaxation parameters.
Complete positivity constrains the elements in the main diagonal of the dissipators as well. The dissipators obtained in Eqs. (5) and (7) have the following constraints: respectively. It is important to note that each i j may have a specific value. As we do not assume any microscopic environment model [24], it is not possible to know which the energy dependences on the parameters are [20]. In this work, we are going to consider in Eqs. (6) and (7) as being a constant parameter with energy dimension. In general, a power-law dependence for the parameter can be defined as an ansatz [20,46] where each power could have a relation with a particular physical mechanism. 1 In the literature there are no any analyses at all from experiments using the dissipative model for three-neutrino oscillations presented in this section. Yet for two-neutrino oscillations there are bounds on decoherence effects [19,20,47]. However, considering an experiment like DUNE [1], where the oscillation pattern must consider the three-neutrino families, the constraint in Eq. (8) corroborates the idea that a complete analysis of the decoherence effect must be made assuming that i j may have different values for each parameter.

Resonance and dissipation in two families
In order to observe the new behavior in the situation where neutrinos propagate in constant density matter, we are going to use the formalism for two-neutrino oscillations introduced in the previous section.
Let us consider the two models that use the dissipators given in Eqs. (6) and (7). We also consider a hypothetical neutrino source whose oscillation channel would be ν e − ν τ . In the standard matter effect, ν e interacts with ordinary matter via a charge current, while ν τ does not. Even considering this 1 For example, in Ref. [20] the power-law was defined as = γ 0 (E/GeV ) n and with n = 2, the agrees with the typical dimension E 2 /M P , which is associated with the possible decoherence effects due to the interaction between the matter with the space-time foam [29,42]. hypothetical source, this oscillation channel is an important component of the oscillation among the three families since the θ 13 mixing angle is not null. The exotic behavior that we are going to show also depends on this fact.
We can write the Hamiltonian given in Eq. (3) in its diagonal form, H S = diag.{Ẽ 1 ,Ẽ 2 }, using the effective mass basis.
The relation between the flavor and effective mass basis is made through the transformation where U is the usual unitary matrix, which depends on the effective mixing angle. The oscillation probabilities can be obtained from and taking the quantum dissipator in Eq. (7), the survival probability is written as where the transition probability can be obtained from where E is the neutrino energy and the survival probability in Eq. (12) depends on sin 2θ, which is usually expressed as where A = 2 √ 2G F n e E. The probability in Eq. (12) shows the decoherence and relaxation effects that are described by the 1 and 2 parameters, respectively. One can obtain the survival probability in the case of the quantum dissipator in Eq. (6) straightforwardly from Eq. (12) setting 2 = 0. The standard oscillation probability is obtained when all the i parameters are null.
We analyze the probability (12) assuming the cases allowed for the constraints given in Eq. (9). So, we address the following cases: (a) 2 = 2 1 , (b) 2 = 1 , and (c) Figure 1 shows the behavior of the oscillation probabilities in all cases analyzed. We consider m 2 = 2.3 × 10 −23 GeV 2 , θ = 9 o , L = 1300 km, and 1 = 10 −23 GeV, which is the current order of magnitude obtained from accelerator and atmospheric experiments in the two-neutrino approximation [19,20]. In this figure, the survival and appearance probabilities for the cases (a), (b), and (c) are represented by upper and lower curves, respectively.
The cases (a) and (b) in Fig. 1 show that the relaxation effect becomes the survival probability, it being smaller than the standard one for each energy point, while in the appearance probability the opposite occurs.
The decoherence effect can be seen through the difference of the oscillation amplitudes between the standard probability and the probability obtained for all dissipative cases. Furthermore, new peaks arise in the cases (a) and (c). These peaks occur around E ≈ 10.52 GeV, the resonance region considering the matter density of the Earth crust [37].
From the mathematical point view, these new peaks in cases (a) and (c) can be explained by the survival probability in Eq. (12). Let us write this probability in terms of three functions given by such that the survival probability is written as Due to the relaxation effect described by the damping term in the function F( 2 , x), the survival probability assumes smaller values than the standard one. This is a feature of the relaxation effect, where it can change the oscillation probability independently of the oscillation parameters.
At the resonance point (E r ∼ 10.52 GeV), we have cos˜ x ≈ 1. Then the H (˜ ,θ, 1 , x) becomes In this form, G(θ, 2 , x) in Eq. (15) and h(˜ ,θ, 1 , x) in Eq. (17) are resonance functions with opposite signs and they are different only by the 1 and 2 parameters. When 1 > 2 the resonance peak of the G(θ, 2 , x) is smaller than the h(θ, 1 , x) and then we have a peak down (up) in the survival (transition) probability. Taking 1 = 2 , the resonance peaks of the G(θ, 2 , x) and h(θ, 1 , x) have the same amplitude and the sum of them is equal to zero. So, there is not a new peak in this case. On the other hand, when 1 < 2 the resonance peak of G(θ, 2 , x) is larger than the h(θ, 1 , x) and then we have the appearance of a peak up (down) in the survival (transition) probability.
From the phenomenological analysis, the new peaks at the resonance region show that the relation between the MSW and decoherence effects is different from the relation between the MSW and relaxation effects.
In the case (a) the relaxation effect is suppressed by the resonance. As we can see, at the resonance region the term cos 2 (2θ) ∼ 0 in the probability given by Eq. (12). Thus, the new peaks in the probabilities in the case (a) tend to recover the standard behavior at the resonance point. In fact, in Fig.  1 the maximum and minimum values of the peaks do not have the same values of the standard probabilities because The case (c) shows how the decoherence effect still eliminates the quantum interference effect, but in this case, in a very subtle way. The MSW effect changes the amplitude of ν e → ν τ with the energy and at the resonance point, the standard oscillation amplitude is maximal, sin 2 (2θ) = 1, but due to the decoherence effect the oscillation amplitude is smaller. At the resonance region the standard survival probability is close to 1, but the same probability is smaller due to the decoherence effect in the case (c). So, in this particular case, the decoherence effect weakly suppresses the MSW effect, increasing the appearance probability at this region.
So, we can conclude that the relaxation and decoherence act as opposite effects when they are combined with the MSW effect.
Through Eq. (14) is clear that the new peaks are not expected for antineutrinos; this equation does not have a resonance behavior in this case. So, the neutrino-antineutrino asymmetry may be changed with the magnitude of the decoherence and relaxation effects.
As mentioned before, solar and astrophysical neutrinos are responsible by stringent bonds on relaxation effect [36]. Therefore, the terrestrial experiments can be disregarded in the cases (a) and (b) and these experiments are able to put bounds only on the decoherence effect and depending on the magnitude of this effect a new peak may exist such as occurred with the case (c).

Analytic approach for three-neutrino oscillation
The three-neutrino propagation in constant matter does not have a complete analytic solution and some kind of approximation is necessary in order to obtain the solutions. When a long baseline neutrino experiment is considered, we can find many approximate solutions, as in Refs. [48,49].
We are going to use the same condition as used by Freund in Ref. [48] where the solution is obtained from the series expansions up to first order in α = m 2 12 / m 2 31 , and with this approximation, the author mapped the effective mixing angles that was reproduced in Appendix A. Under this condition, we are going to include the decoherence effects in the survival and appearance probabilities.
The complete Hamiltonian in the flavor basis is expressed as The U are the usual unitary matrices and we can write the general proba-bility including the decoherence effects as whereŨ αi is an element of the approximate effective mixing matrix 2 and all the quartic products in the probability in Eq. (19) can be obtained from the expressions that we reproduce in Appendix A. The jk parameters describe the decoherence effects given in Eq.
With the current values for α and θ 13 the probability in (19) is an acceptable approximation when compared with the exact solution [48], even when the decoherence effects are taken into account, as shown in Fig. 2. The standard survival probability obtained from (19) agrees with the exact survival probability, but when the decoherence effects are included, the curves may separate from each other from the resonance region depending only on the decoherence magnitude. The standard appearance probability obtained from Eq. (19) and the exact solution have a small discrepancy at the low energy region and when the decoherence effects are included a new difference appears at the resonance region.
It is not easy to see how the decoherence effects act on the behavior of the probabilities in Eq. (19). It is possible to obtain short expressions for the probabilities keeping terms 2 More details can be found in Ref. [48].
When the decoherence parameters are not null the canceling of the divergences at the resonance region fails even when all decoherence parameters have the same magnitude. Although, if we consider the energy range important for the DUNE experiment, where the use of these probabilities will be interesting, the approximate and the exact probabilities have a similar behavior in most of the energy range even when the decoherence effect is taken into account, as in Fig.  2. The larger difference just occurs for the appearance case at the resonance region depending on the decoherence magnitude.
In concrete cases, all calculations for experimental analysis using the probability in Eq. (19) may not have any advantage over the exact approach; even the shorter approximate probability presents many terms. However, the probabilities in Eqs. (24) and (28) are able to show details as regards the behaviors of the probability in Eq. (19) and the numerical probability. So, we are going to use them to investigate how each decoherence parameter changes the oscillation probabilities.
To this end, we consider the DUNE baseline and use the exact approach to show the behaviors of the probabilities and analytical approach to explain the modifications. This is possible since Fig. 2 shows the agreement between the analytical and exact approach on the DUNE energy range. For simplicity, we have used the following values for oscillation Fig. 3 Numeric behavior obtained for survival and appearance probabilities considering the DUNE baseline. In these plots it was supposed that all the decoherence parameters have the same value parameters: m 2 12 = 8 × 10 −23 GeV, m 2 31 = 2.5 × 10 −23 GeV, θ 23 = 0.74, θ 12 = 0.58, θ 13 = 0.14, δ = 0, and the usual Earth density [37]. However, an important point is that the mixing angles values are close to the current best fit of the literature, but the decoherence effects tend to change the best fit values of the mixing angles. Figure 3 shows the curves for the survival and appearance probabilities using the exact approach with different values for the decoherence parameters. For simplicity, we assumed the same value of i j for each probability. Besides, these probabilities are limited between 0 and 1 as a consequence of the complete positivity which is guaranteed by the inequalities in (9).
In Fig. 4 we show the changes due to each term of the decoherence effects. In this case, in order to see each of these effects on its own, the inequalities in Eq. (8) were not respected.

The survival probability
Considering the probability in Eq. (24), we can see that the part expressed in Eq. (27) has the dominant oscillation term, being proportional only to sin 2 2θ 23 . Then 32 would only be responsible for this decrease in the oscillation amplitude. The decoherence parameters 21 and 31 appear in terms that depend on α 2 or α sin θ 13 , so they are subdominant in this oscillation regime.
As we mentioned before the probability in Eq. (24) does not work outside the DUNE energy. From Fig. 4, we can see the consequences that 21 and 31 bring to the survival probability. In this case, 31 just changes the behavior at the resonance region where the survival probability decreases its value. This is explained through the neutrino appearance at this region where a non-usual appearance occurs due to the MSW effect and the decoherence effect described by 31 . In the same way, 21 also decreases the values for survival probability to high energy. In this case, as P ν μ ν μ + P ν μ ν e + P ν μ ν τ = 1 and there is no new neutrino appearance at all (Fig. 4), we can just conclude that there is an interesting transition to P ν μ ν τ at this region if 21 has high enough magnitude.

The appearance probability
In the case of the appearance probability the decoherence parameters may bring about exotic peaks at the resonance region. We can use the appearance probability in Eq. (28) to explain how the decoherence effects act on the neutrino behavior.
The 21 parameter tends to decrease weakly the first oscillation peak. The approximate probability does not describe this effect. Only with probability in Eq. (19), which has high order terms, one is able to describe this behavior. On the other hand, 21 tends to increase the second peak and the term P 3 in Eq. (30) is responsible for this behavior. The P 3 term does not oscillate at this energy range and its individual value is increased to low energy values; this would contribute to the neutrino appearance.
The 31 parameter is responsible for the phenomenon that was discussed in Sect. I I I . However, it is clear that the new peak presented before is outside the energy range of the neutrinos in the DUNE experiment, but it will be possible to bound the 31 parameter through this experiment; depending of the magnitude of this effect, this peak may exist and maybe it can be investigated experimentally in another energy configuration that DUNE or another future experiment may have.
So, considering the current prospect for DUNE and its energy range, the 31 parameter tends to decrease the appear- new peak due to MSW and the decoherence effect 31 appears at the resonance region ance probability at the first peak. The important term for this phenomenon is P 0 in Eq. (29) because the oscillation amplitude at this energy region is larger than the term P cos in Eq. (31). In addition, when the P cos term has 31 = 0, it tends to generate a phase difference in relation to P cos when 31 = 0, besides the oscillation amplitude of this term being increased when 31 = 0. Then the combination between P 0 and P cos is the origin of a phase difference in relation to the standard case and the maximum value is slightly larger than the standard case.
The phenomenon associated with the 32 parameter is very subtle and high order terms in α are necessary to describe the slight decrease of the amplitude oscillation that is peculiar to the decoherence parameter.

CP-violation
With the next generation of long baseline experiments there will be the possibility to study CP-violation with a sensibility as never before. As each decoherence parameter changes the appearance probability in a specific way, we can expect some change in the CP-violation behavior. In general, the CP-violation can be probed through the quantity called the CP-asymmetry, which is defined usually as and in the case of the dissipative approach, this quantity is also changed due to decoherence effects. However, it is not simple to see the changes through this definition. So, we are just to use another definition for CP-violation [37], given by where the CP-violation effect is clearer when the decoherence effect is taken into account. Using Eq. (34), we plot Fig. 5a with δ = 0 and Fig. 5b with δ = π/2. In both cases the behaviors of P with and without the decoherence effects appear. In Fig. 5a the CP-violation is due to the matter potential, usually called fake CP-violation because the CP phase is equal to zero (or δ = π ). In Fig. 5b the CP-violation is due to the CP phase, that is, π/2 and also due to the matter potential. In both Fig. 5a and b the peak is due to the MSW and 13 effects as we have seen before and it is weakly influenced by Pν μ →ν e , which is basically constant at this region.
In order to see how decoherence effects may change the CP-violation behavior, we can use Eq. (34). However, through this equation we can only see the action of the decoherence effects on the intrinsic and fake CP-violation as a whole. The matter effect can be lessened through the choice of the baseline, but it is not completely eliminated. In any approximation for neutrino propagation in matter, the matter potential is kept and these two types of CP-violation will be involved. Then only in an approximate way can we see how the decoherence effects change the intrinsic CP-violation. To this end, we can use the following definition [50]: where the matter effect is lessened as we see through the short approximate appearance probability in Eq. (28) and so, we can see approximately how the intrinsic CP-violation is changed in Fig. 5c. In this figure, the peak at the resonance region disappears and the decoherence effects decrease the amplitude oscillation of the δ P.

Comments and conclusion
We have investigated from a phenomenological point view the dissipative effects in neutrino propagation considering long baseline experiments. In our discussion, we assumed the baseline of the DUNE experiment [1]. This experiment will have greatest sensibility to many oscillation parameters. Besides, it will be able to investigate the open questions in neutrinos physics and many non-standard models will be probed by this experiment [1]. Our approach included dissipative effects with a treatment where all quantum interpretations remain as usual. All dissipative effects were obtained from phenomenology arguments based on a complete positivity constraint. From this, we introduce the quantum dissipator in Eq. (5) as being the most effective among all that we can obtain for terrestrial experiments. It contains three possible decoherence parameters, since the other two possible relaxation parameters are stringently bounded by solar neutrinos. These parameters are correlated by inequalities in Eq. (8) where at least two decoherence parameters have to be non-null in order to keep the physical evolution.
Due to the new effect introduced by the combination between the MSW and the dissipative effects, we made a study in the two-neutrino approximation to explain how the decoherence and relaxation effects react when the MSW is present. In this case, we have shown that the decoherence and relaxation effects brought about different behaviors to the probabilities. These behaviors are evidenced through the new up and down peaks at the resonance region.
For three-neutrino oscillations the usual approximation was used where one assumes m 2 12 << m 2 31 and we rewrite the important mixing coefficients in Appendix A that was mapped in Ref. [48]. So, we presented in an analytical way two versions for the dissipative model obtained from the dissipator defined in Eq. (5). As to the first version, defined by Eq. 19, we used only the coefficients in Appendix A to calculate all quartic product of theŨ αk in Eq. (19). So, we can know the approximate behaviors of the survival and appearance probabilities in this case. These approximate probabilities and the numerical solution have similar behaviors in most of the energy spectrum as we can see in Fig. 2.
In the second version, only terms proportional to α 2 , α sin θ 13 , and sin 2 θ 13 are kept in the probabilities. In this case, on the addition of the decoherence parameters the probabilities become divergent at the resonance region, but Eqs. (24) and (28) are useful to understand the behaviors of the survival and appearance probabilities in the relevant energy range of the DUNE experiment. Figure 4 presents all decoherence parameters individually showing how each decoherence parameter may change the behavior of the probabilities. The parameters 12 and 31 are important in the appearance case, while the 32 is impor-tant in the survival case. In particular, the dissipative and MSW effect was investigated in detail in two-neutrino oscillation where at the resonance region new peaks appear due to the action of both the MSW and the dissipative effects. In three-neutrino oscillation the peak due to decoherence effect appears too, but in this case, it is described by the parameter 31 . We finished this work showing the consequences of the decoherence effects on the CP-violation for specific angles where, in an approximate way, we can distinguish the fake and intrinsic CP-violation case. As we have seen the decoherence effect tends to decrease the amplitude behavior of the CP-violation effects as in Fig. 5.
So, we have shown how the decoherence effects may be non-usual standard phenomena that interestingly may be investigated in the next generation of neutrino experiments. New effects and behavior may be expected and their absences will be due to stringent limits on the decoherence effect. sin 2 2θ 13 = 2αÂ(−Â + cos 2θ 13 ) sin 2 θ 12 sin 2 2θ 13 C 4 + sin 2 2θ 13 sin 2θ 12 = 2 sinθ 12 , (A6) sin 2θ 23 = 2αÂ cos δ sin 2θ 12 sin θ 13 cos θ 23 −1 −Ĉ +Â cos 2θ 13 + sin 2θ 23 , where all the expressions considered are only of first order in α. More simplification can be done if we disregard the subleading terms that come from terms that depend on α × θ 13 .