Dissipative Effect in Long Baseline Neutrino Experiments

The propagation of neutrinos in long baselines experiments may be influenced by dissipation effects. Using Lindblad Master Equation we evolve neutrinos taking into account these dissipative effects. The MSW and the dissipative effects may change the probabilities behavior. 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 decoherence parameters change the probabilities analyzing the possible modification using numeric and analytic approach.

The propagation of neutrinos in long baselines experiments may be influenced by dissipation effects. Using Lindblad Master Equation we evolve neutrinos taking into account these dissipative effects. The MSW and the dissipative effects may change the probabilities behavior. 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 decoherence parameters change the probabilities analyzing the possible modification using numeric and analytic approach.

I. INTRODUCTION
In the near future there will be new long baselines experiments [1][2][3] to test the standard pattern of neutrino oscillation as never before. These experiments will have great sensitivity to determine the parameters that describe the standard oscillation pattern. Through them will be possible investigate the open questions, like CPviolation, mass hierarchy and octant problem [4][5][6][7][8][9][10][11][12][13], and maybe new phenomena.
As it is well known, the quantum mechanics explains how the neutrinos are able to change your flavors during their propagation. The success of the model opened opportunity for testing other kinds of effects in oscillation experiments [14][15][16][17][18][19][20][21]. So, considering that the precision level that the next generation of experiments will achieve, the limit on each usual neutrino oscillation parameter will be very stringent and the space for new physics may also be lessened.
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 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 to neutrino propagation 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 Lindblad Master Equation is non-unitary and it adds in the evolution the possibility of occurring dissipative effects like decoherence, relaxation and other. The quantum evolu- * Electronic address: robertol@ifi.unicamp.br tion through Lindblad Master Equation is non-unitary and it adds in the evolution the possibility of occurring dissipative effects like decoherence, relaxation and other.
In neutrino oscillation, the decoherence effect has been studied more than the other dissipative effects and the most part 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 that is much less studied in the literature even when the propagation is in vacuum. This effect are not important for terrestrial long baseline and only the solar neutrinos are responsible for putting the most stringent bound on 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 two neutrino approximation and will point out which 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 three neutrino oscillations. Then, we will evolve the neutrinos in an exact and analytical scheme considering the baseline of DUNE [1].
In analytical case, there is a study that can be found in Ref. [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 a long and short analytical expressions for the survival and appearance probabilities will be introduced. Besides, we will include a discussion about the range of valid of these new probabilities.
We will show the behavior of these two versions of 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 short version of the analytical probabilities for DUNE, we are able to understand how the behavior 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 range energy important in DUNE.
Then, we conclude with a simple study about CPviolation where we show how the CP-violation phenomena may be changed depending on the magnitude of the decoherence effects.

II. 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 being the decoherence and relaxation effects 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 energy dimension being defined as where V k are dissipative operators that arises from interaction between the subsystem of interest and the quantum reservoir. So, V k may describe the dissipative effect of the neutrino propagating in quantum gravity spacetime [27-29, 41, 42] or in flutuating matter [30,31]. Besides, V k operators act only on the N -dimensional ρ ν space. In the evolution equation (1), the first term on right side evolves the quantum state like Liouville equation, while the second term, which depends on V k operators, becomes the evolution non-unitary. The second term is also responsible for introducing many kinds of the dissipative effects [22][23][24][25] We assume that the standard oscillation Hamiltonian is the same for neutrinos as subsystem of interest [36]. For next generation of long baseline experiment, 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 right side of the last equality, we are assuming to be possible to write the Hamiltonian in effective mass basis. For our aim, it is useful to expand the 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 dimension, respectively. Thus, the evolution equation can be rewritten as with D µ0 = D 0ν = 0 to keep the probability conservation and we change t → x as usual for ultra-relativistic approximation. Asρ 0 (t) = 0, its solution is trivial and it is just ρ 0 (t) = 1/N , where N is the number of families. We do not include this component in equation above and besides, the component a 0 from the V k operators must be null. The D kl matrix in Eq. (4) has many parameters even in two-neutrino approximation. However, if we impose [H S , V k ] = 0 the D kl is written as for three-neutrino approach and for two-neutrino approximation. We relabeled by Γ ij 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 Γ ij 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 indice [20,32,43].
The constraint [H S , V k ] = 0 imposes energy conservation on subsystem of interest which, in this case, are the neutrinos. However, this constraint is very stringent and physically there is not any guarantee that this must occur, because 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 that is called relaxation effect.
As the decoherence and relaxation effects depend on propagation distance [36,43], solar and astrophysical neutrinos are only 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 about 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 k ] = 0, the dissipator considering 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 k ] = 0. However, we can ignore all offdiagonal terms in (7) due to complete positivity [45]. As D kl must be positive definite, the diagonal parameters must be larger than the off-diagonal 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 main diagonal of the dissipators as well. The dissipators obtained in Eqs. (5) and (7) have the following constraints: and respectively. It is important to note that each Γ ij may have a specific value.
As we do not assume any microscopic environment model [24], it is not possible to know which are the energy dependence on the Γ parameters [20]. In this work, we are going to consider Γ in eq. (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 relation with a particular physical mechanism 1 .
In the literature there are not any analysis from experiments using the dissipative model for three-neutrino oscillations presented in this section. Although for twoneutrino 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 that complete analysis of the decoherence effect must be made assuming that Γ ij may have different values for each parameter.

III. RESONANCE AND DISSIPATION IN TWO FAMILIES
In order to observe the new behavior in the situation where neutrinos propagate in density constant 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 by in Eq. (6) and (7). We also consider a hypothetical neutrino source whose oscillation channel would be ν e − ν τ . In standard matter effect, ν e interacts with ordinary matter via charge-current while ν τ does not. Even considering this 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 a 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 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 that is the possible decoherence effects due to the interaction between the matter with the space-time foam [29,42].
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 Γ 1 and Γ 2 parameters, respectively. One can obtain the survival probability to the case of the quantum dissipator in Eq. (6) straightforward from Eq. (12) setting Γ 2 = 0. The standard oscillation probability is obtained when all the Γ i parameters are null.
We analyze the probability (12) supposing cases allowed for constraints given in Eq. (9). So, we hold the following cases: a-) The Fig. 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 that is the current order of magnitude obtained from accelerator and atmospheric experiments in 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 shows that relaxation effect becomes the survival probability smaller than the standard one for each energy point, while in the appearance probability occurs the opposite.
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 that is 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 from 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 In this form, G(θ, Γ 2 , x) in Eq. (15) and h(∆,θ, Γ 1 , x) in Eq. (17) are resonance functions with opposite sign and they are different only by Γ 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 peak 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 the 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-) tends 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 Γ 1 = 0.
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 decoherence effect in the case c-). So, in this particular case, the decoherence effect suppresses weakly 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 the Eq. (14) is clear that the new peaks are not expected for antineutrinos once this equation does not have a resonance behavior in this case. So, the asymmetry neutrino-antineutrino 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 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-).

IV. 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 it is considered a long baseline neutrino experiment, we can find many approximated solution, as in Refs. [48,49].
We are going to use the same condition used by M. 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 we was reproduced in the 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 matrix and we can write the general probability 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. (5) and the ∆ jk are the approximate eigenvalues from Hamiltonian 2 More details can be found in Ref [48]. (18). Explicitly, ∆ jk are given by With the current values for α and θ 13 the probability in (19) is an acceptable approximation when they compared with the exact solution [48], even when the decoherence effects are taken into account, as it is 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 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 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 proportional to α 2 , α sin θ 13 and sin 2 θ 13 [48,49]. This procedure allows us to write the survival probability as where α cos δ cos 2 θ 23 sin 2θ 12 sin θ 13 sin 2θ 23 and the appearance probability is written as where P 0 = sin 2 2θ 13 sin 2 θ 23 P cos δ = α cos δ sin 2θ 12 sin θ 13 sin 2θ 23 and P sin δ = α sin δ sin 2θ 12 sin θ 13 sin 2θ 23 where ∆ = ∆m 2 31 /4E. Disregarding the decoherence parameters, the approximate probabilities in Eqs. (24) and (28) are continuous function and the apparent divergences for A → 1 and A → 0 are canceled by composition of the terms of these probabilities. This same situation was discussed in Ref. [49] and besides, without the decorehence parameters in Eq. (28) , it is possible to obtain the same expression for the appearance probability found in Ref. [48,49].
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 range energy important for DUNE experiment, where the use of these probabilities will be interesting, the approximate and the exact probabilities have a similar behavior in most part of the range energy even when the decoherence effect is taken into account, as it is possible to see in Fig. 2. The larger difference just occurs for the appearance case at the resonance region depending on the decoherence magnitude.
In concrete case, 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 present many terms. However, the probabilities in Eqs. (24) and (28) are able to show details about 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 are consider the DUNE baseline and use the exact approach to show the behaviors of the probabilities and analytical approach to explain the modifications. It is possible since the Fig. 2 show the agreement between the analytical and exact approach on the DUNE range energy. For simplicity, we have used the following values for oscillation 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.
The Fig. 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 Γ ij for each probability. Besides, these probabilities are limited between 0 and 1 as consequence of the complete positivity which is guaranteed by the inequalities in (9).
In the Fig. 4 we show the changes due to each term of the decoherence effects. In this case, in order to see as each one these effects alone, the inequalities in the Eq. (8) were not respected.

A. 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 the 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 thorough the neutrino appearance at this region where it occurs a non-usual appearance 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 not any new neutrino appearance (Fig. 4), we can just conclude that there is an interesting transition to P νµντ at this region if Γ 21 has enough magnitude.

B. The appearance probability
In the case of the appearance probability the decoherence parameters may bring exotics 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), that has high order terms, is able to describe this behavior. On the other hand, Γ 21 tends to increase de second peak and the term P 3 in Eq. (30) is responsible for this behavior. The P 3 term do 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 the section III. However, it is clear that the new peak presented before is outside the range energy of the neutrinos in the DUNE experiment, but it will be possible to bound the Γ 31 parameter through this experiment and depending of the magnitude of this effect, this peak may exist and maybe it can investigated experimentally in other energy configuration that DUNE or another future experiment will may have.
So, considering the current prospect for DUNE and its range energy, the Γ 31 parameter tends to decrease the appearance 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 is increased when Γ 31 = 0. Then, the combination between P 0 and P cos gives origin to 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 are 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.

V. CP-VIOLATION
With the next generation of long baselines experiments there will be the possibilty of to study CP-violation with 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 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 thought this definition. So, we are just to use another definition for CP-violation [37] that is given by where the CP-violation effect is more clear when the decoherence effect is taken into account. Using the Eq. (34), we plot the Fig. 5A with δ = 0 and the Fig. 5B with δ = π/2. In both cases appear the behaviors of ∆P with and without the decoherence effects. In Fig. 5A the CP-violation is due to the matter potential, usually called of fake CP-violation because the CP phase is equal to zero (or δ = π). In Fig. 5B the CP-violation is due to CP phase that is π/2 and also due to matter potential. In both Figs. 5A and 5B the peak is due to the MSW and Γ 13 effects as we have seen before and it is weakly influenced by Pν µ→νe that is basically constant at this region.
In order to see how decoherence effects may change the CP-violation behavior, we can use the Eq. (34). However, through this equation we can only see the action of the decoherence effects on the intrinsic and fake CPviolation 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 approximate way we can see how the decoherence effects change the intrinsic CP-violation. For this end, we can use the following definition [50] where we the matter effect is lessened as it is possible to 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 .

VI. 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 supposed a 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 nonstandard models will be probed from 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 in 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 remain the physical evolution.
Due to the new effect introduced by combination between the MSW and the dissipative effects, we made a study in 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 a different behaviors to the probabilities. These behaviors are evidenced trough the new up and down peaks at the resonance region.
For three neutrino oscillations was used the usual approximation 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 analytical way two versions for the dissipative model obtained from the dissipator defined in Eq. (5). The first version defined by Eq. 19, we used only the coefficients in the 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 numerical solution have similar behaviors in most part 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 these case, the addiction of the decoherence parameters become the probabilities divergent at the resonance region, but the Eqs. (24) and (28) are useful to understand the behaviors of the survival and appearance probabilities in the important range energy of the DUNE experiment.
The Fig. 4 presented 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 important in survival case. In particular, the dissipative and MSW effect was investigated in detail in two neutrino oscillation where at the resonance region appears new peaks due to action of the both MSW and dissipative effects. In three neutrino oscillation the peak due to decoherence effect appears too, but in this case, it is described by Γ 31 parameter.
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 it is possible to see in Fig.  5.
So, we have shown how the decoherence effects may be a non-usual standard phenomena interesting to be investigated in 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.