Quantum coherence and entanglement in neutral-current neutrino oscillation in matter

Although neutrino–antineutrino states originating from neutral-current interactions are blind concerning the flavor state, an oscillation pattern is predicted provided that both neutrino and antineutrino are detected. This issue arises from both the coherence and entanglement of the neutrino–antineutrino states. Based on quantum resource theory, we use the l1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l_1$$\end{document}-norm and concurrence to quantify quantum coherence and entanglement, respectively. Considering the localization properties by the wave packet approach and a matter potential which appears when neutrino and antineutrino propagate in a material medium, we obtain the l1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l_1$$\end{document}-norm and concurrence. We see that the neutrino and antineutrino remain entangled for larger baseline lengths when they propagate in a material medium. In the case of the coherence property, the l1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l_1$$\end{document}-norm decreases in comparison to the corresponding one in vacuum. However, its damping occurs for larger distances.


I. INTRODUCTION
Quantum coherence and entanglement are two important quantum resources that play crucial roles in quantum information and quantum computation science and they are investigated in various branches of condensed matter, atomic physics, and quantum optics.Additionally, several particle physics phenomena, such as meson mixing and neutrino oscillation, could potentially provide suitable opportunities for their explorations [1][2][3][4][5][6][7][8][9][10][11][12][13].For instance, violation of the Leggett-Garg inequality as a manifestation of quantum coherence has been investigated using MINOS data [1].The l 1 -norm value as a measure of quantum coherence has also been studied using data from the Day Bay, KamLAND, MINOS and T2K neutrino experiments [2].For instance, in the case of the Daya Bay experiment, in which detectors are near the source and the measured transition probability is low, the l 1 -norm takes a value much lower than its maximum, while in the case of KamLAND, in which the neutrino baseline is about 180 km, this coherence measure reaches the maximal value.Certainly, the neutrino oscillation process, which is a quantum mechanical phenomenon occurring on the macroscopic scale, is an appropriate medium in which to study the quantum foundation aspects and quantum correlations such as quantum coherence and entanglement.Using the wave packet approach, one can give a meticulous description of neutrino oscillation [14][15][16][17].Specifically, the production, propagation and detection of a neutrino should be considered as localized processes, and this localization is very well fulfilled using the wave packet approach.To be more precise, we must emphasize that the observation of neutrino oscillation depends on the coherence of neutrinos during production, propagation, and detection processes.The production and detection coherence conditions are satisfied provided that the intrinsic quantum mechanical energy uncertainties during these processes are large compared to the energy difference ∆E ij of different neutrino mass eigenstates: where σ E = min{σ prod E , σ det E }.This condition implies that during the production and detection processes, one cannot discriminate the neutrino mass eigenstates.Conservation of coherence during propagation means that wave packets describing the mass eigenstates overlap from the production to the detection regions.The wave packets describing different neutrino mass eigenstates propagate with different group velocities.After propagating L, the separation of different mass wave packets is where v g is the average group velocity of wave packets of different neutrino mass eigenstates, and σ xν is their common effective spatial width.In other words, similar to the double-slit experiment, if one could determine which mass eigenstate is created or detected, the neutrino oscillation pattern would disappear.In particular, the conservation of energy and momentum implies that exact determination of the energy-momentum of charged leptons leads to determination of the mass eigenstate of the corresponding neutrino (in fact, exact momentum conservation causes the neutrino state to be kinematically entangled with the corresponding charged lepton state), and the neutrino oscillation ceases [18,19].Furthermore, the neutrino oscillation process is observed provided that neutrinos have a specific initial flavor state.The kinematic analysis shows that a neutrino state created through charged current (CC) interactions exhibits this situation.For instance, a muon neutrino is created by pion decay while muon decay only gives an electron neutrino.In contrast, the neutral current (NC) processes or the Z 0 decay are blind with respect to the neutrino flavor states.In other words, every flavor eigenstate and every mass eigenstate is created with equal probability.Therefore, one can write, in general, the state of the created neutrino and antineutrino as follows: Here, the second equality is satisfied because the transformation matrix between mass and flavor eigenstates is unitary.Thus, the usual neutrino oscillation in which either neutrino or antineutrino is detected cannot be observed for NC neutrinos.However, Eq. ( 3) shows that there is another property that is noticeable: the neutrino and antineutrino originating from the NC interaction are maximally entangled due to the conservation of energy-momentum and lepton number in electroweak interactions.Hence, if both neutrino and antineutrino are detected, it is possible that oscillation patterns will be observed between detectors [22][23][24].However, if only either neutrino or antineutrino is detected, we will have an entirely classical ensemble and will be unable to observe any oscillation pattern.Indeed, this phenomenon can be interpreted as a CPT transformation of half of the process maps of the propagation of neutrino-antineutrino pairs onto the traditional neutrino oscillation experiment.Of course, its requirement is that the initial state is an entangled state according to Eq. ( 3).Otherwise, if the initial state is separable as follows: using the CPT transformation on the antineutrino process, one can show that the probability of detecting a neutrino with flavor state |ν α at a detector and an antineutrino with flavor state |ν β at another detector is which is not the traditional neutrino oscillation probability.Although NC neutrino oscillation is beyond the reach of any experiment, it can be implemented analogically by the double-double-slit experiment performed in Ref. [25].
In this experiment, the path-entangled photons are passed through opposite screens with a double slit.Due to entanglement, each photon can reveal the which-slit-path information of the other photon.Therefore, the two-photon interference pattern appears if detection locations of photons are correlated without revealing the which-slit-path information.However, it has also been shown experimentally and theoretically that two-photon quantum interference disappears when the which-slit-path of one photon is detected in the double-double-slit.When neutrinos (antineutrinos) propagate in a material medium, they might have a forward coherent scattering off electrons and nucleons of the matter via both weak CC and weak NC interactions [26,27].All three neutrino (antineutrino) flavor states can interact with the matter through NC, but the electron neutrino (antineutrino) has an additional interaction (i.e.interaction through the CC) which causes it to feel an additional potential.In the calculations, we only consider the effective potential due to the CC interaction because the potentials coming from the NC interactions induce common phases for all three neutrino (antineutrino) flavor states and do not alter neutrino oscillating behaviors.Therefore, the effective potential of matter which must be considered is , where G F and N e are the Fermi coupling constant and the density of electrons in the medium, respectively.If the matter density of the medium is nonuniform, the potential depends on the coordinate; otherwise, it is constant.
According to Eq. (3), the neutrino and antineutrino produced by NC interactions are entangled with respect to either their flavor or mass modes.Indeed, this state has a maximal entanglement similar to the Bell states 1 .Moreover, they 1 Bell states are the four states that can be created when two qubits are maximally entangled.The four states are represented as are entangled due to energy-momentum conservation.In fact, the predicted oscillation pattern between two detectors is based on these features.Also, as a result of entanglement, the coherence condition obtained by considering localization properties via the wave packet approach is stronger than that obtained by imposing the non-separation mass wave packet constraint analogous to Eq. (2) [23,24].In fact, the coherent propagation of both neutrino and antineutrino is not sufficient because the oscillation pattern ceases if the distance between the detectors is larger than the coherence length.In this paper, we consider the matter refraction to reanalyze this problem in the context of twoflavor neutrino oscillations.In this case, the NC neutrino oscillation cannot be concluded by a CPT transformation of half the process.However, since the mixing in the source is independent of matter potential, the initial state can be written in mass eigenstate according to Eq. (3).Therefore, as we will show, one can obtain the standard oscillation pattern in a material medium for the NC neutrinos.Moreover, using quantum resource theory, we study the quantum entanglement and quantum coherence under these conditions.For this purpose, concurrence and l 1 -norm are defined as a measure for entanglement and quantum coherence, respectively.Coherence is the quantum resource of neutrino oscillation, and its origin during production and detection processes, as well as its conservation during propagation, has been discussed extensively and accurately in the literature; for instance, see [14][15][16][17][18][19][20][21].Meanwhile, using quantum resource theory, one can quantify quantum coherence, entanglement, or another quantum resource by introducing some measure quantities [31][32][33][34][35][36].In the case of usual neutrino oscillation, the refraction of the neutrino in a dense medium leads to a decrease in the quantum coherence in standard neutrino oscillation [9].However, as we will see, the value of entanglement of neutrino-antineutrino pairs shows different behavior in terms of matter potential.
In the next section, treating the neutrino and antineutrino as a wave packet, we obtain the probability of the NC neutrino oscillation in a material medium with constant density for the two-flavor schema.In section III, we explore the variation in entanglement and quantum coherence for this issue in the quantum resource theory framework.For this purpose, we use concurrence and l 1 -norm measures for entanglement and coherence, respectively.We summarize our results in the last section.

II. PROBABILITY OF NC NEUTRINO OSCILLATION IN MATTER
According to the quantum wave packet approach, we should describe the neutrino and antineutrino by a localized wave function at the space-time coordinates (x, t) and (x, t), respectively.Moreover, it is assumed that the neutrino and antineutrino detectors are located at distances L and L from the source, respectively.The neutrino and antineutrino coming from the Z 0 boson must be described by a bipartite entangled state.Therefore, this state can be written as follows: where The superscript m refers to the propagation in a material medium, and the subscript S indicates that this state is related to the source.Here, f S (p, p m i ) and fS (p, pm i ) are the momentum distribution functions of the neutrino and antineutrino with mean momenta p m i and pm i .Also, E m i (p) and Ēm i (p) are the energies of neutrino and antineutrino with mass m i , respectively, and t p and x p denote the production coordinates.The Dirac delta factor is added in Eq. ( 7) to ensure that the momentum is conserved (here, we consider those Z 0 bosons as being almost at rest in the frame where the matter effect is constant).Indeed, the time evolution of neutrino states in a matter medium is governed by the following effective Hamiltonian [28]: In the case of antineutrino, we denote the corresponding Hamiltonian by δ Hm , and it is obtained from Eq. ( 8) by replacing V (x) with −V (x).We consider the density of matter to be constant.The unitary mixing matrix U for a two-generation platform is parameterized as follows: Therefore, one can obtain the following eigenvalues for the Hamiltonian given by Eq. ( 8): and and the corresponding mixing angle is given by In the case of the antineutrino, we denote the eigenvalues of δ Hm and corresponding mixing angle by δ Ēm 1,2 and θm , respectively.These parameters can be represented by the corresponding expressions given in Eqs.(10), (11) and (12) with the difference that V must be replaced by −V .It is clear from Eq. ( 12) that there is a resonance in mixing of neutrinos provided that V ∆m 2 21 cos 2θ/2E.However, no resonance exists for antineutrino mixing.For the momentum distribution function of the neutrino and antineutrino in Eq. ( 7), we use the Gaussian momentum wave function as follows: Hereafter, the momentum uncertainties of the neutrino and antineutrino are denoted by σ p and σp , respectively.It is appropriate to assume that σ p and σp are much smaller than the corresponding mean momenta.Therefore, one can expand the energy up to the first order of the departure from the mean momentum as follows: where and v m i is the group velocity of the ith mass eigenstate and is given by v m i = dE m /dp| p=p m i .Likewise, we can write a similar expression for antineutrinos.Consequently, one can write the position wave packet given in Eq. ( 7) as follows: where σ 2 xν and σ2 xν are the position uncertainties of the neutrino and antineutrino in the source, and they are given by σ x = 1/2σ p and σx = 1/2σ p , respectively.
On the other hand, since the detection processes are essentially time-independent, detected states have no time dependence.Therefore, the wave function of the detected neutrino and antineutrino states is described by where Ψ m iD (x − L) and Ψm j D (x − L) are the wave functions of the detected neutrino and antineutrino at positions L and L, respectively.Moreover, U m αi and Ū m βi in Eq. ( 17) are mixing matrices of neutrinos and antineutrinos in a material medium.In fact, U m ( Ū m ) has a form similar to Eq.( 9) with the difference that θ is replaced by θ m ( θm ) given by Eq. (12).Again, we assume a localized Gaussian wave function in momentum space similar to (13) for both the neutrino and antineutrino detected in the corresponding detectors.Under this scenario, the detected state Eq.( 17) can be written as follows where σ 2 xD and σ2 x D are the uncertainties of the detected neutrino and antineutrino processes.In general, the mean momenta of the produced particles, p i and pi , are different from one of the detected particles, p i and p i .However, we assume that they coincide.Then, the transition amplitude of detecting ν α and νβ in the corresponding detectors becomes: Substituting Eq.( 6) with Ψ m iS given in Eqs. ( 16) and ( 18) into the above equation and calculating the integrals over the position coordinates, we obtain where we assumed that the position uncertainties of neutrinos and antineutrinos which are defined as are approximately equal σ 2 x σ2 x .We want to obtain the probability of detecting a neutrino with flavor state |ν α at one detector and the corresponding antineutrino with flavor state |ν β at another.In general, the distance of the two detectors from the production location may be different.In this case, the arrival times of the neutrino and antineutrino at their corresponding detectors are not equal.It is suitable that instead of t and t, we select the following time variables: Meanwhile, the emission time and the corresponding arrival times of the neutrino and antineutrino are not measured.Hence, we integrate over T , which is the arrival time of the neutrino or antineutrino from its own detector to another detector.Therefore, the probability of the oscillation pattern between detectors can be obtained as follows: in which ∆E m ij (∆ Ēm ij ) and ∆v m ij (∆v m ij ) are the energy and group velocity differences in the material medium between the two mass eigenstates of the neutrino (antineutrino), respectively.By using relativistic approximations, one can write the mean energies as follows: and where E is the neutrino and antineutrino energy in the limit of zero mass.Also, ξ denotes a dimensionless quantity whose value can be estimated from energy-momentum conservation in the production process [29,30].In the case of Z 0 boson decay in the rest frame, we have ξ ≈ 0. Similarly, the corresponding momenta can be given by and Furthermore, in this approximation, the group velocity of a mass eigenstate is also written as for neutrinos and for antineutrinos.Therefore, Eq. ( 23) can be simplified as follows: where the oscillation length L mosc ij and the coherence length L mcoh ij , for i = j, are defined by in which ∆δE m ij ≡ δE m i − δE m j and ∆δ Ēm ij ≡ δ Ēm i − δ Ēm j , and respectively.According to Eq. ( 28), the differences in the neutrino group velocities appearing in Eq. ( 32) become For ∆v m ij , we obtain a similar expression, with the difference that ∆δE m ij and V must be replaced by ∆δ Ēm ij and −V .For neutrinos, ∆v ij becomes zero provided that V = ∆m 2  21 /2E cos 2θ.But in the case of antineutrinos, it does not happen.Therefore, according to Eq. ( 32), the greatest value of L mcoh ij occurs when ∆v m ij = 0.The last exponential factor appearing in Eq. ( 30) is due to the non-simultaneous neutrino and antineutrino detection processes.This factor depends on the subtraction of the difference between the energies and group velocities of the neutrino and antineutrino mass eigenstates.This difference is zero in vacuum, and in the material medium with usual densities its effect can be completely ignored compared to the oscillation and decoherence factors.Therefore, with accuracy such that this factor can be ignored, the processes of neutrino and antineutrino detection can be considered to be simultaneous.Furthermore, the first three terms appearing in the first exponential, in Eq. ( 30), are the same as the results of Refs.[23,24] with the difference that L mosc ij and L mcoh ij are modified due to the neutrino and antineutrino scattering off matter.However, the last term in this exponential is coming absolutely from the difference between the potential of neutrinos and antineutrinos.Similar to the terms coming from the non-simultaneity of neutrino and antineutrino detection processes, the effect of this term can also be ignored.
Since the lifetime of Z 0 is so small (about τ Z0 3 × 10 −25 s), the interval between two nearest collisions becomes smaller than the mean distance between particles in the nucleus.Therefore, σ x is determined by the Z 0 decay width Γ Z0 as follows [16]: The energy of the neutrino and antineutrino coming from the decay of a rest Z 0 is about E = 46GeV.If we take ∆m 2 21 = 7.53 × 10 −5 eV 2 , and θ = 33.46• [37], σ x of approximately 10 −16 m is obtained.With these values for relevant parameters, using (32) and (31), coherence and oscillation length of about 10 10 m and 10 15 m are obtained, respectively.Therefore, under this condition, the oscillation pattern ceases due to the separation of mass eigenstate wave packets.This situation is not improved by choosing the oscillation parameters between the second and third generations, and the oscillation length will still be several orders of magnitude larger than the coherence length.In the following, however, we want to investigate quantum coherence and entanglement as a resource of the oscillation pattern in this problem.Therefore, we suppose σ x to be about the atomic distance (5 × 10 −10 m) so that the oscillating behavior does not cease.For greater clarification, we illustrate the probabilities of the survival (a) and transition patterns (b) between detectors versus their distance L + L in Fig. 1.Specifically, we plot the probability in terms of L + L for four cases: vacuum (V = 0, blue curve), for a value of the potential which corresponds to the neutrino resonance mixing (V = 3.208 × 10 −16 eV, pink curve), for a value of the potential in which L mcoh 12 takes the maximum value (V = 2.087 × 10 −15 eV, green curve) and for a matter-dominated potential (V = 5 × 10 −14 eV, brown curve).

III. QUANTUM CORRELATIONS
As mentioned, the NC neutrino oscillation originates from the neutrino and antineutrino entanglement and quantum coherence due to the overlap of mass eigenstate wave functions.Both of these correlations are affected by the disentangling and decoherence due to the separation of wave packets as well as the scattering of neutrinos and antineutrinos off material medium.Therefore, in this section, we use concurrence and l 1 -norm, respectively, as a measure for entanglement and quantum coherence in the framework of quantum resource theory, in order to investigate the mentioned effects.
Entanglement as a quantum correlation is studied extensively in the framework of the quantum resource theory [31][32][33][34].Concurrence is one of the measures suggested to quantify entanglement.First, let us define ρ as follows: The concurrence measure is generally defined by [33,34]: where λ i 's are the square roots of the four eigenvalues of the non-Hermitian matrix ρρ in decreasing order.Accordingly, considering a general two-qubit state, we have The concurrence for this state is obtained as follows Clearly, for Bell states, we have C(ρ) = 1.Otherwise, if the coefficients in Eq. ( 37) are such that |φ is separable, C(ρ) will be zero.For other cases, one obtains a value between 0 and 1.The more entangled the state, the closer the value of C(ρ) is to 1.
For an entangled neutrino and antineutrino state originating from the Z 0 decay, the corresponding bipartite state is given in Eq. ( 6) which is written based on the wave packet approach.The density matrix operator corresponding to this state, after integrating over the related momenta and time2 , is obtained as Taking into account the issues related to energies, momenta and group velocities, given in Eqs. ( 24)-( 29), and using the definitions given in Eqs. ( 31) and (32), one can simplify this relation as follows: The behavior of this density matrix in terms of x + x is analogous to that obtained in Ref. [21] with some differences coming from the distinct propagation of neutrinos and antineutrinos in a material medium.Now we can calculate the concurrence whose definition is given in Eq. ( 36) as a measure of the neutrino and antineutrino entanglement.By direct calculation, one can obtain the following expression for concurrence: In fact, the first term in the exponential has a decisive role in the behavior of C(ρ) when the propagation length is of the order of the coherence length.Furthermore, we should emphasis that if one uses a plane wave description for neutrinos, the value of concurrence remains equal to 1 for every value of x + x and V .Meanwhile, in the wave packet approach, the concurrence depends on V as well as x + x; it becomes zero when x + x exceeds the coherence length, and its falling is postponed for higher matter density (see Fig. 2).Specifically, according to Eq. ( 41), when the localization properties are considered by the wave packet approach, the overlap of mass eigenstate wave packets is diminished for a propagation length larger than the coherence length because the corresponding velocities are different.However, it can be compensated by propagation of the neutrino in the material medium.In other words, the refraction of neutrinos in a dense medium causes the difference in wave packet velocities to decrease.While the concurrence vanishes for smaller matter potentials for a specific propagation length, it becomes nonzero for larger potentials; for example, see Fig. 3. Now let us study the coherence of neutrinos and antineutrinos originating from the Z 0 decay process.As was noted, in order for oscillation to occur, neutrinos must be coherent during the production, propagation, and detection processes.In the framework of the quantum resource theory, quantum coherence is a quantum correlations that can be quantified by appropriate measure quantities [35,36].Among the several measures for quantum coherence, a very intuitive one is related to the off-diagonal elements of the considered quantum state.Hence, for a given density matrix, ρ, the l 1 -norm as a measure of quantum coherence is defined by In general, the maximum possible value for c(ρ) is d − 1, where d is the dimension of the corresponding density matrix [36].Therefore, the neutrino oscillation shows that the off-diagonal elements of the corresponding density matrix in the basis of the flavor state are nonzero, and consequently we have no vanishing l 1 − norm [9].
In the case of NC neutrinos, the dimension of the density matrix given in Eq. ( 40) is four in two-flavor schema.In this case, therefore, the maximum value of the l 1 -norm is 3.According to the definition given in Eq. (42), one can write this parameter in terms of the transition amplitudes as follows: Here, the first and second terms and the fifth and sixth terms are identical.An explicit form of c(ρ m ) in terms of neutrino oscillation parameters is given in the Appendix.We illustrate the behavior of c(ρ m ) versus the distance of detectors x + x through Fig. 4. In this plot, we consider three values for matter potential: vacuum (V = 0), the value corresponding to the infinite coherence length for neutrinos, and a matter-dominated value.We see that although for baseline less than the coherence length the matter effect causes c(ρ m ) to drop relative to vacuum, the coherence length in a material medium becomes larger than one in the vacuum.Furthermore, in Fig. 5, the behavior of c(ρ m ) in terms of the matter potential is depicted for the results obtained based on both plane wave (a) and wave packet (b) approaches.We consider two values for the baseline: a value smaller than the coherence length, 10 16 m (blue curve), and a value greater than that, 10 18 m (pink curve).In the case of the former (blue curve), c(ρ m ) obtained in both approaches behaves similarly.Although for some values of the potential a slight increase in the value of c(ρ m ) is seen, as the value of the potential increases, c(ρ m ) eventually decreases to a value less than the corresponding one in the vacuum.However, in the case of 10 18 m baseline (pink curve), we see that in the plane wave approach, c(ρ m ) achieves the maximum value for a region of the potential value, while in the wave packet approach, c(ρ m ) is certainly smaller than the maximum value.This is due to the decoherence effects coming from the separation of neutrino wave packets.

IV. CONCLUSION
According to the theory of neutrino oscillation, it is possible to observe an oscillation pattern between the detectors of neutrinos and antineutrinos coming from Z 0 boson decay when the neutrino production, propagation, and detection coherence conditions are satisfied.We see that under a realistic condition for Z 0 boson decay, σ x is estimated to be about 10 −16 m.Thus the coherence length is smaller than the oscillation length, and the propagation coherence conditions are not satisfied.Nevertheless, we have assumed σ x = 5 × 10 −10 so that the coherence conditions are satisfied.However, if only one of them is detected, neutrino oscillation cannot be observed.This phenomenon results from the fact that the neutrino and antineutrino arising from the Z 0 decay are both entangled and coherent.These quantum properties can be quantified according to the quantum resource theory.In this paper, we use the concurrence and l 1 -norm to quantify entanglement and quantum coherence, respectively.On the other hand, considering localization properties by using the wave packet approach, one can see that the wave packet separation of neutrino mass eigenstates due to the difference in the corresponding velocities suppresses the amount of entanglement and coherence measures.Another factor that can affect these quantum correlations is the propagation of neutrinos and antineutrinos in material media.Hence, in this paper, we have reanalyzed the Z 0 decay neutrino oscillation by considering the localization properties and matter potential.In particular, we have obtained the concurrence and l 1 -norm in two-flavor schema.If the localization effects are ignored and in other words, the plane wave approach is used to obtain the transition probability, the concurrence always remains 1.Otherwise, the concurrence ceases when the propagation length exceeds the coherence length.However, the matter potential causes the coherence length to increase, and as a result, this damping occurs in larger propagation lengths (see Fig. 2).Indeed, the propagation in material medium causes wave packet separation, with the result that the damping of entanglement is compensated to some extent (see Fig. 3).Furthermore, we have shown that although the matter potential causes the l 1 -norm value to be small compared to the corresponding one in the vacuum, the damping due to the wave packet separation occurs in larger propagation lengths as well (see Fig. 4).As another point, we have compared the behavior of the l 1 -norm in terms of matter potential for the case when using the plane wave approach and when using the wave packet approach.Obviously, when the propagation length is smaller than the coherence length, both approaches give the same result.Otherwise, in the plane wave approach, the coherence value is maintained, but it is ultimately suppressed due to the interaction with the material medium.In the wave packet approach, its value is suppressed before the significant impact of the potential (see Fig. 5).

FIG. 2 :
FIG.2: Concurrence versus propagation length for various matter potentials.The second and third values of potential have been chosen to correspond to the resonance and the infinite coherence length of a neutrino, respectively.