Decoherence effects in reactor and Gallium neutrino oscillation experiments: a QFT approach

We adopt the quantum field theoretical method to calculate the amplitude and event rate for a neutrino oscillation experiment, considering neutrino production, propagation and detection as a single process. This method allows to take into account decoherence effects in the transition amplitude induced by the quantum mechanical uncertainties of all particles involved in the process. We extend the method to include coherence loss due to interactions with the environment, similar to collisional line broadening. In addition to generic decoherence induced at the amplitude level, the formalism allows to include, in a straightforward way, additional damping effects related to phase-space integrals over momenta of unobserved particles as well as other classical averaging effects. We apply this method to neutrino oscillation searches at reactor and Gallium experiments and confirm that quantum decoherence is many orders of magnitudes smaller than classical averaging effects and therefore unobservable. The method used here can be applied with minimal modifications also to other types of oscillation experiments, e.g., accelerator based beam experiments.

In this paper we contribute to this discussion by adopting the quantum field theoretical (QFT) approach to neutrino oscillations [4,6,12,38,39].In this approach the combined process of neutrino production, propagation, and detection is considered as a single process whose amplitude is calculated by usual S-matrix methods, adapted to the situation of macroscopically separated production and detection regions.The neutrino is treated as an internal line and coherence properties of the flavour transition are completely determined by the localization of the external particles at the source and at the detector.
We will follow largely the methods discussed by Beuthe in the review article Ref. [12], with some modifications relevant to the experimental situations of interest to us.In particular, we generalize the formalism to take into account decaying particles (see also [14,[40][41][42]) as well as coherence loss due to interactions of external particles with the environment.We focus on quantities which are actually observed in experiments, for instance the energy of the positron produced by the inverse beta-decay reaction in the detectors of reactor neutrino experiments.The formalism naturally allows a distinction of "quantum decoherence" (at the amplitude level) compared to "classical averaging" (at the probability level), which are, however, indistinguishable observationally.Our work is complementary to the recent paper by Akhmedov & Smirnov [27], who base their argumentation on the neutrino wave packet approach, reaching very similar conclusions as we do.A somewhat different approach has been pursued by Jones, Marzec & Spitz [43] whose results for the decoherence parameters differ quantitatively from ours.
We provide here some guideline on how to read this paper.In section 2 we introduce the notation by discussing external wave packets for single vertex processes, such as scattering or decay.Our ansatz of how to include decoherence effects due to interactions with the environment as well as decaying particles is introduced in section 2.1.In section 3 we sketch the QFT calculations for the oscillation amplitudes.A reader mainly interested in the final result may, after a look at the central expression in eq.(3.14), skip directly to section 3.2, which contains a brief discussion of the decoherence effects, and section 3.4, where we comment on classical averaging and stress the equivalence of quantum and classical decoherence.Section 4 contains the numerical estimates for reactor and Gallium source experiments, where in section 4.1 we give some details of how we estimate the localizations of all involved particles.Readers interested only in the main numerical results may proceed directly to section 4.2, where we discuss the effective localization and energy spreads and show that they are many orders of magnitude below the observable level.We summarize our findings in section 5. Supplementary material is provided in appendices A to C. In particular, appendix B contains a discussion of how to derive standard expressions for scattering cross sections and decay rates within our formalism, and in appendix B.3 we discuss how the amplitude for the neutrino oscillation process derived in the QFT formalism can be related to differential event rates in a neutrino oscillation experiment.

Single-vertex calculation
In the QFT approach to neutrino oscillations the combined process of neutrino production, propagation, and detection is evaluated.There, we face the unusual situation that two vertices of the process (production and detection) are macroscopically separated in spacetime.Before we calculate the amplitude for the oscillation process in section 3, we discuss first the conceptually simpler case, where there is only one interaction region, as in standard particle physics interaction calculations.This will serve to introduce the formalism, fix the notation, and to discuss the modifications we introduce to model the specific physical situations of decaying particles as well as particles confined by frequent interactions with the environment.We consider wave packets for both, initial and final state particles.In appendix B we show that, in the limit of plane waves and the appropriate normalization, we can recover standard text book expressions for scattering cross sections and decay rates from our approach.
We define states as superpositions of momentum eigenstates |k as Three-vectors are denoted by bold-face letters.Momentum states are normalized as and the wave packets as In the following we will consider the specific case of Gaussian wave packets.While this is not true in general, it serves as a useful approximation for the purpose of describing the relevant physics and for order-of-magnitude estimates.The technical advantage is that many integrals can be performed analytically in the case of Gaussian wave packtes.Taking into account the above normalization condition, we have Here, p is the mean momentum and σ is the momentum spread.We will use the symbol "σ" to denote uncertainties of dimension 1 (i.e., momentum or energy) and the symbol "δ" to denote uncertainties of dimension −1 (i.e., space or time); they are related via σδ = 1/2.
Let us now consider a process with a set of N i initial and N f final state particles, |φ i and |φ f , respectively.They are the product of states, each having the form as defined in eq.(2.1).Then, the total transition probability for the process i → f is given by the S-matrix [44,45]: where iT is the non-trivial part of the S matrix: S = 1 + iT .We use eq.(2.1) and the standard definition of the matrix element for momentum states [44], where in general M depends on all momenta.Then, we obtain the transition amplitude as In the second step we used the Fourier-transform of the δ-function to introduce an integral over configuration space d 4 x.This integral will receive contributions from the "interaction region" in space-time, determined by the overlap of the wave packets.We now proceed by performing the momentum integrals by adopting the following approximations.First, we assume that the matrix element M varies slowly with momenta on the scales σ i,f and we can replace M(k i , k f ) ≈ M(p i , p f ).Then it can be pulled out of the integrals.Second, consider the exponentials e ik i,f x .For a generic k i,f we have kx = E k t−kx.Again we assume that the momentum is not too far from its mean value p and approximate with the group velocity The physical implication of the approximation eq.(2.9) is that the wave packets are not spreading in time. 1 Stationary states correspond to v = 0, whereas for free particles we have v = p/E p .At this point we do not need to specify the dispersion relation and can allow for arbitrary k dependence of E k .With these two approximations all integrals become integrals over Gaussians, which can be performed, see e.g., [12].The momentum integrals are all of the type Using this in eq.(2.8), we obtain . (2.12) In order to simplify the exponential we introduce some notation: such that ∆E = 0 and ∆p = 0 correspond to exact energy and momentum conservation, respectively.We define the total momentum spread by and a weighted velocity and velocity-squared: Using these definitions to rewrite the argument of the exponential, the d 4 x integration can be performed as well: where the effective energy spread is obtained as We notice that the quantities σ e and σ p correspond to the amounts within which energy and momentum conservation can be violated.It is easy to see that they fulfill 0 ≤ σ e ≤ σ p .The momentum uncertainty, eq.(2.14), is just the sum of the squares of the momentum uncertainty of all particles.Hence, it is dominated by the particle with the largest uncertainty.
For hierarchical uncertainties, σ e is generally dominated by the particle with the second largest uncertainty, see e.g., [12] for a discussion.Note that "energy-momentum violation" happens for the mean quantities p i , p f , whereas for the actual 4-momenta picked at the interaction vertex, energy-momentum conservation is exact, as manifest by the δ-function in eq.(2.7).Therefore, it is not surprising that the mean quantities need to fulfill 4-momentum conservation only approximately.

Interactions with the environment and particle decay
We consider now the case that the involved particles propagate through a medium, with frequent collisions.When an interaction happens, the particle wave function picks up a random phase, which leads to loss of coherence.We model this effect by restricting the time integration in eq.(2.16) from t = 0 till the time when the next interaction happens by introducing an exponential e −tΓ col /2 , where Γ col is the collision rate.This is motivated by considering the scattering as a Poissonian process. 2 Hence, the integral in eq.(2.16) becomes with Ω ≡ ∆E − ∆pv, and in the last step we assumed σ e Γ col .Hence, after squaring the amplitude, instead of the energy Gaussian factor in eq.(2.17) we obtain the typical Lorentzian shape ∝ 1/(Ω 2 + Γ 2 col /4).This modelling describes correctly the collisional line broadening [46], see appendix B.2.
Actually, we will take this picture as the physical origin of wave packets for our external particles: the interactions with the environment lead to the localization of the particle within a length δ x which in turn determines the momentum spread σ via σδ x = 1/2 as well as the mean time between two scatterings as δ x /|v|.We adopt now the following ansatz for the effective collision rate Γ col Γ 2 col ≡ i,j i.e., adding the interaction rates in squares. 3Hence, Γ col is dominated by the particle with the fastest interaction rate: as soon as one of the involved particles interacts with the environment, coherence is lost.Comparing eq.(2.21) with the expression for the energy spread in eq.(2.18), we see that σ e ≤ Γ col .In section 4 we will find for our cases of interest that σ e ≈ Γ col because in most cases v 2 Σ.As we have seen, our physical picture leads to a constraint on Ω with Lorentzian shape.In the following we will, however, replace the Lorenzian with a corresponding Gaussian factor with width Γ col .This approximation will still capture the relevant physics, but simplify the calculations significantly, because integrals can be taken analytically.Hence, we recover an expression for the amplitude, similar to eq. (2.17), but with the replacement In neutrino oscillation experiments neutrinos are very often produced in particle decays, e.g., beta-decay of a nucleus in reactor experiments or pion decay in accelerator experiments.We can model the decay of a particle in complete analogy to the interaction with the environment by introducing an additional factor e −tΓ dec /2 in (2.19).Hence, the width of the Lorentzian becomes just Γ col +Γ dec .As we show in appendix B.2, in the limit Γ col Γ dec (frequent collisions) this ansatz describes the collisional line broadening of a decay line, whereas in the limit Γ col Γ dec (fast decay and negligible collisions) it reproduces the standard definition of the decay width in terms of the matrix element.In the Gaussian approximation we simply obtain the effective energy spread by For neutrino production in a nuclear reactor as well as the source in Gallium experiments, the lifetimes of the relevant beta decays are typically much longer than the interaction times with the environment (see section 4), and therefore, it is save to negect Γ dec compared to Γ col .This may not be the case for experiments using pion decay, where pions decay in a decay tunnel essentially without interacting.
Although we can combine the energy spread σ e , the collision rate, and the decay rate in the same effective quantity σ E (thanks to the Gaussian approximation for all of them) their physical origins are different.1/σ e corresponds to the effective interaction time due to the wave packet overlap of all involved particels, Γ col describes the coherence loss due to interactions of each of the involved particles with the environment, and Γ dec takes into account the energy spread due to the finite lifetime of a decaying particle.

Neutrino oscillation amplitude and event rate
We now move to the discussion of the amplitude relevant for neutrino oscillation experiments consisting of neutrino production, propagation and detection.To be specific, we consider neutrino production by the decay of a particle A into two final state particles and an antineutrino, A → 1 + 2 + ν, and anti-neutrino detection via the process B + ν → 3 + 4. We have in mind reactor neutrinos, where the production process corresponds to the beta decay of a nucleus (A), and the detection process is the inverse beta decay reaction on a proton (B), but many of our considerations will apply also in other circumstances with minor modifications.The total process A + B → 1 + 2 + 3 + 4 is illustrated in fig. 1.The neutrino is considered as an internal propagtor and does not appear as external particle [4,6,38].We now follow the common approach [12] and calculate the amplitude for the total process, by assuming wave packets for all external particles, both initial state (A, B) as well as final state particles (1,2,3,4).
We proceed in complete analogy to the discussion in section 2 but generalize it to the case of macroscopically separated production (P ) and detection (D) regions.In analogy to eq. (2.17), we obtain the following expression for the amplitude describing production of an anti-neutrino with flavour α and detection of an anti-neutrino with flavour β4 (see e.g., [12,14] for explicit derivations): Here, U αj are elements of the PMNS mixing matrix, the normalization factors N i are defined in eq.(2.11), MP,D are the reduced matrix elements of the production and detection processes, the sum over j runs over the neutrino mass states with neutrino mass m j , x P,D are space-time points located in the production and detection region.In the second line of eq.(3.1) we obtain two Gaussian factors related to the approximate energy-momentum conservation at production and detection points, with the momentum spreads σ pI , energy spreads σ EI , and velocities v I (with I = P, D) defined as in eq.(2.14), eq.(2.23), eq.(2.15), respectively.Furthermore, we have defined the kinematic 4-momenta of the neutrino at the production and detection vertices: and E P , E D are the time-components of the corresponding 4-vectors.
Next we perform the integral over d 3 p by using the Grimus-Stockinger theorem [6], which allows us to take into account the macroscopic separation of source and detector.In the relevant limit the propagating neutrinos go on-shell, and we obtain where ) is obtained from the exponential in the second line of eq.(3.1),where according to the Grimus-Stockinger theorem p → p j with Following Ref. [12], we rewrite f j (p 0 ) by decomposing the vectors in components parallel and orthogonal to l: where p I and v I denote the components parallel to l and we have redefined In order to keep notation concise we do not introduce a different symbol for this new variable.

Derivation of the decoherence terms
Equation (3.4) is the starting point to derive the terms leading to decoherence effects, which will be the main focus of our considerations in the following.We are not concerned with overall factors and focus on the interference terms.First, we have to square the amplitude and perform an integration of the unobservable propagation time T , see discussion in appendix B.3, eq.(B.17).We obtain a δ-function from the integral dT e −i(p 0 −p 0 )T and obtain where we used that neutrino masses are small, m j p 0 , and expand the square root in the neutrino momenta as p j ≈ p 0 − m 2 j /(2p 0 ).We see that the time integration implies that only neutrinos with the same energy can interfere and different neutrino energies are summed incoherently [12,14].
The interference terms correspond to j = k in eq.(3.8), for which we recognize the familiar oscillation phase depending on ∆m 2 kj ≡ m 2 k − m 2 j .Let us now simplify the discussion and specialize to the case of two neutrino states with masses m 1 , m 2 , and write We are interested only in terms proportional to ∆m 2 .Terms depending only on the absolute neutrino mass m 2 will lead to (tiny) irrelevant global corrections which do not affect the interference term.Furthermore, the quantities m 2 I = E 2 I −p 2 I (for I = P, D) are of order of the neutrino mass-squared m 2 and also independent of the neutrino mass indices j, k.Therefore, we neglect also terms proportional to m 2 I and set E I ≈ p I .Using p j ≈ p 0 − m 2 j /(2p 0 ) also in the function f j (p 0 ) and dropping all terms proportional to m 2 and independent of ∆m 2 , we obtain at leading order in ∆m 2 : Already at this stage we obtain the term with σ m which potentially can lead to decoherence.We will comment on its physical interpretation in section 3.2.
Next we perform the integral over p 0 .We expand the first term in eq.(3.10) around its minimum, which is at p 0 = E 0 with The oscillatory phase i∆m 2 L/(2p 0 ) is expanded to first order in (p 0 −E 0 ), and we set p 0 = E 0 in the term with σ m (ignoring higher order correction to this decoherence term).Then, the p 0 -integral can be performed with the method outlined in appendix A and we obtain This is a central result of the QFT approach to the oscillation process; let us briefly comment on the terms appearing here.In the first line we have the standard oscilation phase and a term depending on the net momentum components orthogonal to the neutrino direction.The latter term will constrain the phase space integrals for the orthogonal components within the mometum spreads σ pI around zero.Because of the relabeling of the energy variables E P,D in eq.(3.7), also the oscillation phase depends on p I⊥ and in principle the integral over these components can lead to additional non-trivial effects (we comment on it in appendix C).Note that in eq.(3.14) we have three quantities corresponding to an effective neutrino energy: the kinematic neutrino energies E P and E D at the production and detection vertices defined in eqs.(3.2) and (3.3), respectively, and E 0 which is a weighted mean value of the former two.The last term in eq.(3.14) ensures, that all "three neutrino energies" are the same within quantum mechanical uncertainties, determined by the sum of the effective energy-momentum uncertainties at source and detector.Equation (3.14) is completely symmetric with respect to source and detector.

Discussion of the decoherence terms
The first two terms in the second line of eq.(3.14) describe the exponential damping due to decoherence.We provide here a brief review of the two decoherence terms, see e.g., Refs.[12][13][14]47] for more discussions.The two terms correspond to two generic types of damping.We define with We introduced the generic energy spread σ en , where in eq.(3.14) we have σ en = σ eff , and E ν is a relevant neutrino energy, E 0 in the version obtained in eq.(3.14).
Both decoherence terms have a simple physical interpretation.Starting with ξ loc , the interpretation in the energy-representation of this term is that the energy-momentum uncertainty encoded in σ m needs to be large enough, such that individual mass states can neither be resolved at the production nor at the detection process.If energy-momentum was defined with an accuracy better than ∆m 2 /E ν , the individual neutrino mass states would be determined and no interference of different mass states would be possible.
Another interpretation of this term becomes apparent in the spatial representation: defining the oscillation length by the decoherence terms can be written as In this form we see that ξ loc can also be interpreted as the condition δ loc L osc , i.e., that both, source and detection points need to be localized better than the oscillation length, by noting that δ 2 loc = δ 2 P + δ 2 D .Moving now to ξ en , this term says that for an experiment around the oscillation maximum (L/L osc 1), the neutrino energy needs to be sufficiently determined, such that σ en E ν ; for experiments beyond the oscillation maximum the condition becomes correspondingly stronger.Using σ en = σ eff with σ eff defined in eqs.(3.11) and (3.13), this is a condition on the energy and momentum spreads of external particles of the neutrino production and detection processes.For L L osc , decoherence will become relevant only if σ en E ν , i.e., the quantum mechanical uncertainty on the neutrino energy needs to become comparable to the neutrino energy itself, implying that the neutrino would not have a well defined energy.Note that we work under the assumption that momentum spreads are small compared to relevant momenta or energies, and therefore our approximations adopted in section 2 to perform the integrals may not apply if σ en ∼ E ν .
In the above derivation we have integrated first over T and then over p 0 , which is more convenient to derive decoherence terms.If the order of the T and p 0 integrals are exchanged, one can see that an exponential factor makes sure that the T integral is dominated by values of T constrained by |T − L − Lm 2 j /(2E 2 0 )| 1/σ eff , which can be interpreted as relating L T v j , with v j ≈ 1 − m 2 j /(2E 2 0 ) corresponding to the "velocity" of the neutrino with mass m j [12,48].This suggests a wave packet interpretation of the internal neutrino in the QFT approach [14].In that picture the damping due to ξ en can be interpreted as wave packet separation of the propagating neutrinos.
We note that the two terms ξ loc and ξ en have the opposite dependence on the spreads σ m and σ en .This means that for oscillations to be observable, quantum uncertainties have to be big enough that different mass states can interfere (σ m ∆m 2 /E ν ), but small enough that interference is not damped (σ en E ν L osc /L).Assuming that very roughly σ m ∼ σ en , we see that there are many orders of magnitude available to fulfill both requirements, thanks to the smallness of ∆m 2 /E 2 ν or, in other words, due to the macroscopically large oscillation length, L osc E ν 1.

Phase-space integrals over unobserved external momenta
Equation (3.14) contains the intrinsic quantum mechanical decoherence terms.So far we have only averaged the amplitude squared over the unobservable time T (which actually does not introduce a decoherence term).All other manipulations are performed at the amplitude level; in particular, the integral over p 0 corresponds to an internal particle, which can be performed already at amplitude level.In real experiments there are of course always effects leading to additional averaging at the probability level, i.e., of the amplitude squared (see also section 3.4).Some of these averages are intrinsically unavoidable and related to the physical configuration of the experiment.Let us consider first the case of reactor neutrino oscillation experiments; a very similar discussion applies also to accelerator or atmospheric neutrino experiments.In these experiments usually the neutrino energy is reconstructed in the detector, by measuring all (or some of) the outgoing particles at the detector, i.e., particles 3 and 4 in our example, which allows to reconstruct E D via eqs.(3.3) and (3.7) with some accuracy.In contrast, inital and final state particles at the production point usually are not observed. 5Therefore, in order to calculate event rates R D in the detector, we need to integrate the squared amplitude over the phase space of final state particles in the production reaction, see also appendix B. 3. By a suitable variable transformation, one of these integrals can be chosen to be over E P .As we sketch in appendix C, the decoherence term emerging from this integral has the same shape as ξ en and can be combined with the original term present already in eq.(3.14) such that we obtain for the event rate in the detector Hence, the energy spread in ξ en is given by σ en = σ D,eff , i.e., it depends only on the uncertainties in the detector, while σ m in ξ loc remains unchanged and contains contributions from both, the production and detection process, see eq. (3.12).The integral over E P corresponds to a classical sum, i.e., summing the squared amplitude.The corresponding decoherence can therefore be considered as emerging from classical averaging.Note, however, here this averaging is in principle unavoidable, given the physical observables in a neutrino oscillation experiment.In this way, eq.(3.21) depends only on the "observable neutrino energy" E D .Also note that σ D,eff ≥ σ eff .Therefore, the classical averaging due to the E P integration increases effectively the decoherence, see section 3.4.The fact that only the detection process determines the decoherence in the last term of eq.(3.21) and all production-related uncertainties drop out follows from expanding the oscillation phase around E D and from the Gaussianity of all involved uncertainties; this result may not hold in the most general case.Note that we keep only leading terms in the expansion parameters σ/E with σ ∈ (σ eff , σ P,eff , σ D,eff ) and E ∈ (E 0 , E P , E D ), and differences like σ D,eff /E D − σ P,eff /E P are of higher order in this expansion.However, our result that the relevant energy spread in the ξ en term of eq.(3.21) is larger than σ eff as in eq.(3.14) is robust.
In typical detectors of modern reactor experiments only the energy of the outgoing positron can be determined, whereas the neutron momentum is not measured.Hence, we need an additional phase-space integral over the neutron momentum.This will provide another contribution to σ en , which however, is at most of the same order of magnitude and typically much smaller (see the discussion of the p ⊥ integral in appendix C).
Moving to Gallium radioactive source experiments, these are pure counting experiments, i.e., the detector does not provide energy information, whereas the neutrino source consists of a couple of quasi-monochromatic neutrino lines from an electron-capture decay N → N + ν.In this case one would naturally integrate the phase space over E D instead of E P .The calculation is completely symmetric to the one outlined above and we obtain the same result as in eq.(3.21) with the replacements E D → E P and σ D,eff → σ P,eff .Hence, in this case the energy spread in ξ en is set by the production process and will be determined by the localization of the nuclei N, N as well as the (natural and/or thermal) linewidth of the decay (see section 4).The same comments as above apply regarding the dropping out of σ D,eff .

Classical averaging
In realistic experiments both production point and detection point are known only within some uncertainty, related to the size of the neutrino source and the vertex resolution of the detector.Similarly, detectors can determine particle energy and momentum only within certain resolutions, and hence the neutrino energy E D can be reconstructed only within a finite accuracy.These effects are taken into account in predicted event rates by convoluting the neutrino oscillation probability with the corresponding resolution functions.Hence, a classical average is performed.
Let us approximate these spatial and energy resolutions by Gaussians with widths δ clas and σ clas , correspondingly: ) with R D (L, E D ) given in eq.(3.21).Assuming that the width of the Gaussians is small compared to the other L or E D dependent factors in R D , we can expand the oscillation phase to linear order either in (L − L) or (E D − E D ) and see that the integral takes again the same form as in appendix A. Applying eq.(A.4), we obtain decoherence terms of the same form as ξ loc for the L smearing and ξ en for the energy resolution.Hence, they can be combined with the corresponding terms present in eq.(3.21), which amounts to the replacement Hence, decoherence due to classical averaging has precisely the same effect as intrinsical quantum mechanical decoherence [5,9,50]. 6These results reflect the following two (rather obvious) statements: (i) quantum mechanical uncertainties provide a fundamental lower bound on classical uncertainties, and (ii) in order to observe effects of quantum mechanical decoherence, classical averaging effects have to be suppressed down to the quantum level.

Numerical estimates
The QFT formalism outlined above allows to calculate the relevant uncertainties δ loc (or equivalently σ m ) and σ en relevant for the localization and energy spread decoherence factors ξ loc and ξ en , respectively, from the properties of the involved external particles.The required input for their definitions in eqs.(3.11) and (3.12) are the effective energy and momentum uncertainties eqs.(2.14) and (2.23), which in turn are derived from the momentum spreads of all the external particles in the production and detection processes, as well as their velocities, as defined in eq.(2.10).We will now evaluate δ loc and σ en for reactor neutrino and Gallium radioactive source experiments.

Particle localizations and velocities
First we need to estimate the momentum spreads σ of all involved particles, as well as their velocities v. Similar estimates have been performed recently in [27] in the context of neutrino wave packets.The momentum spread is calculated via the spatial localization δ x , assuming the uncertainty principle δ x σ = 1/2.We list the relevant quantities for all the particles involved in the production and detection processes in reactor and Gallium experiments in table 1.They are estimated as follows.
Reactor experiments.For the initial and final state nuclei N , N in the production process, via beta decay within the nuclear fuel, we assume that the localization is determined by a typical interatomic distance [47].We estimate this by using that the lattice parameter of uranium oxide UO 2 is a = 5.471 × 10 −10 m and has 4 U and 8 O atoms in one unit cell [51], which gives δ x a/12 1/3 0.24 nm.For the initial state nucleus we assume a thermal velocity v = k B T /m, where the temperature in the nuclear fuel ranges from 700 K at the outer egde to 2000 K in the center [52].This is justified, as the fission products termalize on time scales much faster than their beta decay lifetimes [27].For nuclei with mass numbers in the range of 80 to 160, the velocity of fission products lies in the range of 6.3×10 −7 to 1.5×10 −6 .For our estimates a typical value of v N 10 −6 is taken.The recoiling nucleus after beta decay, N , is not in thermal equilibrium and we estimate its velocity using v = |p|/E.Assuming typical neutrino energies E 4 MeV and mass numbers from 80 to 160 we find v N 4 × 10 −5 .The protons in the detector are typically bound in carbon molecules.We assume a localization of 10 −10 m corresponding to the typical size of the C-H bound length.For the velocity we take thermal velocities at room temperature.Concerning the neutron, after being produced it undergoes scattering in the liquid scintilator.The spatial localization is therefore estimated by the mean collision length l = 1/nσ int , where we use the neutron on CH 2 scattering cross section from [53].For a typical kinetic neutron energy of E kin 10 keV a cross section of σ int ≈ 50 b leads to a mean collision length of 4.8 mm.Here, the number density of CH 2 was estimated as n 1/(3d 3 ) for an interatomic distance d 20 nm.The velocity is calculated for a neutron with E kin 10 keV.
For the outgoing electron and positron passing through the medium, either in the nuclear fuel or in the detector, we proceed as follows.We consider the mean rate of energy loss −dE/dx using the "Bethe equation" [54], which we numerically integrate.As localization we take the distance which the particle travels until it deposits one mean excitation energy I.This should provide us with a good estimate for the mean free path of the particle.
For the calculation of −dE/dx we need to know the number density of elementary charges, which can be estimated from the density and composition of the stopping medium.For the reactor fuel we use the above mentioned properties of UO 2 .As a typical detector material we assume linear alkylbenzene (LAB), C 6 H 5 C n H 2n+1 , with n ranging from 10 to 13; for definiteness we will assume n = 12.We calculate the number density of elementary charges as 138N A ρ/m mol , where 138 is the number of elementary charges/electrons in one molecule of LAB, and the density of LAB is taken to be ρ = 0.859 g cm −3 [55].Furthermore, we need the mean excitation energies I.As we consider materials made up of different elements, the calculated mean excitation energies were averaged over the contributions from the different atoms.For atoms with high elementary charge Z a good approximation is I = 10 eV • Z [56].Thus, at the source, this is used for UO 2 , where Z = 36 on average and thus I = 360 eV.At the detector we have LAB which we simplify to CH 2 .Then, the average excitation energy of one molecule can be approximated to be I ≈ 48 eV [56].With these numbers we find that an electron (positron) with an initial kinetic energy of 3 MeV deposits one I energy after traveling 260 nm (320 nm), which we take as the localization at the neutrino production (detection).
Gallium source experiments.To estimate the particle localizations at Gallium experiments we take the BEST experiment as an example [24,25].At the source 51 Cr undergoes electron capture to become 51 V and an electron neutrino.Approximately 90 % of the produced neutrinos have an energy of 750 keV while 10 % have energies of 430 keV.For the sake of defineteness we consider the 750 keV neutrinos.The temperature of the source is approximately the one of the surrounding gallium, namely 300 K [24,25].Therefore, the thermal velocity of 51 Cr (compound of nucleus + electron to be captured) is v Cr = k B T /m = 7.3 × 10 −7 .The final 51 V is not thermal, instead its velocity is calculated as v V = p V /E V = 1.5 × 10 −5 .The localizations of the 51 Cr and 51 V are approximately the same and estimated from the cristal lattice dimensions.They are bound in a bcc lattice with 2 atoms in one unit cell of size a = 2.88 × 10 −10 m [57].This leads to a spatial localization of δ Cr,V ≈ a/ √ 2 = 0.20 nm.At the detection we estimate the thermal velocity of Ga at 300 K to v Ga = k B T /m Ga = 6.3 × 10 −7 .The velocity of Ge is not thermal, instead v Ge = p Ge /E Ge = 1.15 × 10 −5 .For the electron we obtain v e − = p e− /E e− = 0.83.All momenta are estimated to be of the order of E ν = 750 keV.With the density of gallium ρ ≈ 6 g cm −3 , we obtain the number density as n = ρN A /m mol and n −1/3 = 2.67 × 10 −10 m, which we take as the localization for both the Ga and Ge atoms.For the localization of the electron we proceed as described for the reactor case and integrate −dE/dx until the electron deposites one mean excitation energy, which for Ga is I = 11 eV • Z = 341 eV [56].This leads to δ e ≈ 310 nm.
Our estimates for the particles involved in the production processes reported in table 1 are in agreement with the results of [27] within about one order of magnitude.Ref. [27] does not provide detailed estimates for the detector particles.

Effective energy-momentum spreads and decoherence parameters
Now we are in the position to calculate the relevant effective uncertainies.For convenience we summarize again the corresponding relations, as derived in sections 2 and 3. Here, I = P, D labels production and detection processes, and the sum is over all external particles a (initial and final state) of the respective processes: The upper part of table 2 shows our results for these quantities based on the input from table 1. Momentum spreads σ pI are dominated in all cases by the localization of hadronic particles and are of order few 100 eV to keV.Some comments are in order concerning the velocity.The velocity vector v I defined in eq.(4.1) depends on the relative orientation of the individual particle velocities v a .The quantity ṽI given in table 2  in natural units they are related by σ en δ en = 1/2 and σ m δ loc = 1/2.
given by the momentum spreads of the final state leptons e ± , except for the electron capture production (which has no final state charged lepton).Furthermore, |v I | ≤ ṽI √ Σ I implies that σ eI ≈ Γ I,col .Again, an exception is the production in Gallium experiments, where the approximation ṽI √ Σ I is not very good; nevertheless, Γ P,col is still of the same order as σ eP also in this case.Finally, for typical lifetimes for fission products of order 1 ms or larger, decay widths are at most of order Γ P,dec ∼ 10 −12 eV, and the decay width of 51 Cr is Γ Cr,dec ≈ 4×10 −22 eV.Hence, the decay widths are always completely negligible compared to σ eI and Γ I,col and can be savely neglected in eq. ( 4.4) and we have σ EI ≈ √ 2Σ I σ pI ≈ √ 2σ eI .Therefore, also σ EP and σ ED are determined by the momentum spread (or by the localization) of the outgoing charged lepton (again with the exception of the electron-capture source).Now we can move to the calculation of the effective uncertainties summarized in the lower part of table 2. We start with σ m , defined as Here, v I is actually the length of the component of v I parallel to the neutrino propagation direction.The small values of v I ≤ ṽI imply that for reactors, σ m is dominated by σ pI , whereas in Gallium the small value of ṽP is partially compensated by the small value of σ EP , and the energy term gives a non-negligible contribution to σ m , although the order of magnitude remains unchanged and σ m ∼ σ pI .We find for both type of experiments values of order For the localization term ξ loc from eq. (3.15) this implies Hence, ξ loc = 1 for all practical purposes and localization decoherence is irrelevant in reactor and Gallium experiments.We emphasize, however, that classical spatial averaging in real experiments can be significant, with δ clas ∼ 1 m (corresponding to typical sizes of reactor cores or Gallium detectors), which, depending on the value of ∆m 2 , does play an important role and has to be included in the analyis of these types of experiments.
Moving to the energy spread, we note that In the last relation in eq. ( 4.9) we have used that v I ≤ ṽI 1 and σ EI σ pI to obtain σ I,eff ≈ σ EI .Hence, for the effective energy spread we obtain dominated by the outgoing e ± momentum spreads for reactors and by the electron capture decay for Gallium experiments, the latter leading to a value of σ eff more than one order of magnitude smaller.As discussed above, σ eff corresponds to the "pure quantum mechanical" energy spread, c.f., eq.(3.14).However, even in idealized experimental configurations there is some unavoidable averaging, e.g., due to phase-space integrals over unobserved momenta.Therefore, as argued in section 3.3, the relevant energy spread in reactor experiments is σ D,eff , due to the averaging of unobserved momenta at the production region.This does not change the qualitative picture, as in our approximation they are of the same order of magnitude.Numerically we find for the spread relevant for ξ en the values: For the Gallium source experiments, it is more natural to integrate first over the phasespace of the detector particles, as no momenta are measured in the detector.This would imply σ P,eff as the relevant energy spread.Here the effect of the phase-space integration is even less important, as we anyway have σ P,eff ≈ σ eff ≈ 0.0092 eV or δ en ≈ 11 µm.However, in this case another fundamental (though classical) averaging effect needs to be taken into account, namely the Doppler broadening due to the thermal motions of the source particles [27,47].This leads to an energy smearing with Gaussian shape, with the width set by We can include this effect in a straight forward way, following the discussion in section 3.4.We see that σ Doppler σ P,eff and therefore it dominates the energy spread.Using eq.(3.24) with σ clas = σ Doppler , we obtain quite similar to the values obtained for reactor experiments in eq.(4.11).
Note that Doppler broadening is in principle relevant for reactor experiments as well.In the neutrino source it would, however, not induce additional decoherence, as we already integrate over the effective neutrino energy in the source and the additional smearing due to the Doppler broadening would have no effect.However, it does contribute to the energy resolution of the detector.Here, we consider it as part of the classical energy resolution and do not include it in the spreads given in eq.(4.11).
For the energy spread decoherence ξ en written as in eq.(3.20), we find which again implies ξ en = 1, both for reactor and Gallium experiments to very good accuracy.
The energy resolution of typical detectors in reactor experiments is in the range (0.03 − 0.06) MeV E/MeV [20,55,58] and hence, σ clas 0.1 MeV, about 6 orders of magnitude larger as the intrinsic energy spread eq.(4.11).Equation (4.14) shows that for these values, − ln ξ en can become of order one, which just means that the (classical) averaging due to the energy resolution of neutrino detectors is an important effect.On the other hand consider the upcoming high precision JUNO reactor experiment [55], which aims to observe oscillations due to the mass squared difference ∆m 2 31 2.5 × 10 −3 eV 2 at a distance of about 53 km.This implies L/L osc ∼ 13 and − ln ξ en ∼ 5 × 10 −11 for E ν 4 MeV.Therefore, quantum decoherence effects will remain completely unobservable also for JUNO.
Comparison with previous results.Our results of the uncertainties in eqs.(4.11) and (4.13) are in rough agreement with [27], but differ by a factor of 10 3 − 10 4 from [43] where δ en 10 − 400 pm is obtained.The authors of Ref. [43] consider localization scales induced by nucleon-nucleon interactions within the nucleus as well as the size of the decaying nucleon.However, we argue that dynamics at length scales much smaller than 1/|q|, with q being the typical momentum transfer of the reaction, are irrelevant to the problem, as localizations smaller than 1/|q| cannot be resolved.For the same reason scattering of neutrinos with MeV energies can be considered with the nucleus as a whole, and not on individual nucleons (unless they are effectivley "free", as the proton in Hydrogen), or even quarks.Therefore, it is the localization of the nuclues, which is relevant for defining the quantum mechanical uncertainties for the process.See also Ref. [29].
The authors of Refs.[35,59] consider a term similar to ξ en and perform a phenomenological analysis using reactor data to set a lower bound on a "coherence length", finding δ > 2.1 × 10 −4 nm, which corresponds to σ = 1/(2δ) < 0.47 MeV.This value has also been adopted by [26] to study the impact on sterile neutrino oscillations in short-baseline experiments.This value of the decoherence parameter is of similar order as the energy resolution of the detectors, 5-6 orders of magnitude larger than the intrinsic quantum uncertainty, and therefore fully dominated by the classical averaging effect [27].

Summary
We have used the QFT approach to neutrino oscillations to estimate quantum decoherence effects for reactor neutrino and Gallium source experiments.In this formalism possible decoherence effects on the oscillation probability are fully determined by specifying the momentum spreads (or equivalently the localizations) of the external particles in neutrino production and detection processes.The neutrino is treated as an internal particle, which is integrated out.It is not necessary to introduce the concept of neutrino wave packets, however, the results can also naturally be interpreted in terms of neutrino wave packets [14].Throughout our calculations we adopt the approximation of a Gaussian shape for all functions describing momentum and energy spread.
In section 2 and section 3 we have reviewed the QFT formalism and derived the coherence factors in terms of external particle localizations and velocities.We introduced the concept of coherence loss due to frequent interactions with the environment, and we find that for the cases of interest this effect is of a similar size as the finite interaction time due to wave packet overlap of the particles involved in the production and detection processes.Furthermore, in our approach we focus on experimentally observable quantities, such as the positron energy after inverse beta decay in the detector, which is used to reconstruct the neutrino energy.We show that phase-space integrals over unobserved particle momenta lead to additional decoherence effects, which correspond to classical averaging of the amplitude squared, but nevertheless are intrinsic to the experimental configuration and in principle unavoidable.Finally, we recover the well-known result that quantum and classical decoherence effects are experimentally indistinguishable, and phenomenologically this distinction is unphysical.Quantum decoherence can only be observed if all classical effects leading to energy and localization averaging can be suppressed down to the quantum level, which of course is not possible in real-life neutrino experiments.
In section 4 we have performed numerical estimates of decoherence effects for reactor and Gallium radioactive source experiments.First we estimated the localizations and velocities of all involved particles in the production and detection processes.From these we calculate the effective decoherence parameters relevant for the two types of decoherence related to localization and energy spread, eqs.(3.19) and (3.20), respectively.The main results are summarized in section 4.2 and table 2. Localization decoherence is controlled by the factor (δ loc /L osc ) 2 .With typical values of δ loc ∼ 0.2 nm, this is many orders of magnitude smaller than any relevant oscillation length L osc and therefore, localization decoherence is completely negligible for practical purposes.The second decoherence effect is related to the energy spread and around the first oscillation maximum it is controled by (σ en /E ν ) 2 .In the wave packet picture this term can be interpreted as decoherence due to wave packet separation.Both, for reactor and Gallium experiments we find values of order σ en ∼ 0.5 eV, about 6 to 7 orders of magnitude smaller than typical neutrino energies E ν , leading again to completely negligible decoherence effects.
Hence, in both cases (localization and energy spread) quantum decoherence effects are completely unobservable and many orders of magnitude smaller than classical averaging effects due to finite neutrino source and detector sizes as well as classical energy resolution effects.This remains true even for the upcoming high precision JUNO reactor experiment.If at some point data might provide evidence for decoherence effects beyond classical averaging, standard quantum mechanical uncertainties as discussed here cannot be responsible and this would definitly require exotic new physics [60][61][62][63].

B Decay rate and scattering cross section
We depart from the amplitude eq.(2.17) derived in section 2 for a single-vertex process, considering the effective energy spread σ E from eq. (2.23).From this result we obtain the total transition probability as: Now we follow Weinberg [45] and outline how to derive decay rates or cross sections from this expression, taking into account the wave packet localization of all external particles.Weinberg considers a cubic box around the interaction point in which all particles are localized.This leads to a discretization of all momenta.Here, we adapt the argumentation to the case when particles are localized as Gaussian wave packets, allowing for different degrees of localization for each particle, see also [65,66].
For each particle a we can consider V a ∼ 1/σ 3 a as the volume in which the particle is confined due to its wave packet.It turns out that a suitable definition is The numerical factor can be motivated as follows: define the configuration space spread δ x by δ x σ = 1/2.Then V = ( √ 2πδ x ) 3 and 1/V would be just the prefactor of a Gaussian in configuration space with width δ x .
Let us now relate the transition probability in eq.(B.1) to more familiar quantities.Consider first the final state particles.P if is the total transition probability for i → f where all particles are defined via the respective wave packets.Usually we are interested in a differential quantity, for a transition into an infinitesimal phase space element for the outgoing particles.Such a quantity can be obtained by multiplying the probability by the differential number dN f of outgoing particles within the phase-space element d 3 p f around momentum p f : with where in the second step we have used eq.(B.2).Hence, each final state factor in eq.(B.1) becomes just the usual phase-space factor: where we have used the definition of N from eq. (2.11).
Next we consider the following factors from eq. (B.1): where we have defined the effective interaction volume and interaction time in analogy to eq. (B.2) by using the effective momentum and energy spreads from eqs. (2.14) and (2.23): In eq.(B.7) we have introduced the notation δ G to denote a Gaussian, which will converge to a delta function in the limit σ → 0 and the limit in eq.(B.8) is obtained for σ p → 0.7 Finally, we use eq.(B.2) also to rewrite the initial state factors N i in terms of volumes and obtain We interpret T I as the effective interaction time.Therefore, we obtain the differential interaction rate dΓ if by dividing the differential transition probability by the interaction time.Equation (B.10) is sometimes called Fermi's Golden Rule.

B.1 Scattering cross section
Consider now the case of the scattering of two initial state particles to a set of final state particles AB → f .We get the factor V I /V A V B = V I n A n B in eq.(B.10), with n A , n B denoting the number density of particles A, B. Let us consider particle A as "target".Then V I n A is the number of target particles inside the interaction volume, and the flux of incoming B particles is v rel n B , with v rel being the relative velocity.The cross section per target particle is given by (total interaction rate) / [(number of target particles A) × (flux of incoming particles B)].Hence, which, in the limit σ p , σ E → 0, where the Gaussian δ G functions become true δ-functions, converges to the standard expression for a cross section, e.g., [44].

B.2 Particle decay
Consider now the case of a single particle in the initial state, i.e., the decay of a particle.In this case we obtain a factor V I /V i in eq.(B.10).We can consider n i = 1/V i as the density of initial state particles and V I as the effective interaction volume, with V I ≤ V i .Hence, V I n i is the number of initial-state particles inside the interaction volume.Usually we are interested in the decay rate of a single particle, denoted just by Γ dec .Therefore, we obtain (B.12) In the limit σ p → 0, where the Gaussian δ G functions become true delta-functions, this expression converges to the standard expression for the decay rate of a particle, see e.g., [44].
Let us now comment on the considerations in section 2.1 related to collisions with the environment.We restore the Lorentzian shape of the energy spread, see eq. (2.20), and instead of eq.(B.10) we obtain where now we identify T I = 1/Γ.We recognise that the energy conservation function δ G (Ω) is now replaced by the Lorentzian, whose width is set by the rate Γ.We consider now two physically different situations: • First we assume that the collision rate Γ col is much larger than the decay rate Γ dec and set Γ = Γ col in eq.(B.13).Let us consider as simple example a massive scalar particle (mass M i ) at rest decaying into two massless scalar particles.The matrix element is a constant, and we can perform the phase space integrals, apart from the energy of one of the final state particles and obtain We observe the collisional line broadening [46], i.e., the width of the decay line is set by the collision rate.Integrating over E f leads to the same total decay rate as for Γ col → 0, namely Γ dec = |M| 2 /(16πM i ).
• Let us now assume that Γ dec Γ col , i.e., we can neglect collisions.Now we can set Γ = Γ dec in eq.(B.13) and use that the probability that the particle decays into any final state is one: dP if = 1.Considering again the example of a massive scalar particle at rest decaying into two massless scalar particles, we obtain from eq. (B.13) in the limit of exact momentum conservation δ In the limit Γ dec M i we can solve for Γ dec and recover the standard result for the decay width Γ dec = |M| 2 /(16πM i ).

B.3 Rate for the neutrino oscillation process
In the case of neutrino oscillations, the observable of interest is the differential event rate at the detector dR D , depending on the momenta of the final states in the detector, i.e., p 3 , p 4 in our notation.In the case of reactor neutrino experiments these are the positron and neutron, whose momenta are in principle observable in the detector.We can obtain this quantity from the transition probability P if = |iA αβ | 2 by the following steps: 1.As in eq.(B.3), we multiply by the phase-space elements of all final state particles in order to obtain a differential transition probability, 2. We are interested in the event rate per single detector particle, as well as single decaying particle at the source.Hence, we have to divide by the number of particles in the production and detection regions, V P /V A and V D /V B , respectively, with the volumina defined as in eqs.(B.2) and (B.9).
3. The outgoing particles in the production region are not observed.Hence we have to integrate over their phase space dp 1 dp 2 .
4. The event rate at the detector is obtained by dividing by the effective detection time interval T D defined as in eq.(B.9), in analogy to the first relation in eq.(B.10).
5. The amplitude squared still depends on production and detection times via T = t D −t P .While the detection time in principle is observable in real-time neutrino detectors, the production time typically is not observable (see comments below).Hence, for a given t D , we have to sum the contributions of all possible production times contributing to the amplitude.This is obtained by first dividing by the effective production time interval T P defined as in eq.(B.9), which gives the neutrino production rate, which then has to be integrated over all possible values of t P .By a simple coordinate shift, this integral is transformed into an integral over T and we obtain a time averaged amplitude-squared: Putting everything together, we obtain for the differential event rate at the detector Following comments are in order: • In reactor experiments typically only the energy of the outgoing positron is used to reconstruct the neutrino energy.In this case we need also to integrate over the phase space of the outgoing neutron in the detector.
• In counting experiments, such as Gallium radioactive source experiments, the detection reaction is not observed and we have to integrate in addition over the phase-space of outgoing detector particles dN 3 dN 4 .
• Equation (B.18) corresponds to the hypothetical "event rate" assuming a single decaying particle at production and a single detector particle.For an actual experiment we have to sum over the number of target particles in the detector, as well as the total number of all decays contributing to the neutrino production at the source.The latter is a complicated calculation in case of reactor neutrinos, amounting to an ab initio calculation of the total reactor flux.
• Above we assumed that the neutrino production time is not observable.This is true for a large number of experiments, such as reactor, solar, or atmospheric neutrinos.At many accelerator experiments, neutrinos are produced by a pulsed particle beam, which in principle provides some time information on the production.Typically the duration of the beam spills are of order µs, which however, in practice is much longer compared to other relevant time scales ∼ 1/σ EP , 1/σ ED .Therefore, for practical purposes it is a good approximation to assume the integration in eq.(B.17) over an infinite time interval [41].
• In the case of a so-called monitored neutrino beam, such as proposed by ENUBET [49], in principle time information on neutrino production can be obtained.In such a situation the amplitude needs to be averaged over the accuracies with which both t D and t P can be determined, instead of an infinite time interval.In order to affect the discussion above, the accuracy of both time measurments, t D and t P , need to be comparable to the microscopic time scales 1/σ EP and 1/σ ED .
• Using the formalism discussed in this paper, one can show [48] that in the (realistic) case, when decoherence effects play no role, eq.(B.18) factorizes into the decay rate Γ P of the production process, the standard oscillation probability and the detection cross section σ D as see also [11,42,67] for similar calculations.

C Phase-space integrals
As discussed in section 3.3, depending on the physical configuration, integrals over the momenta of unobserved external particles either in the source or in the detector are unavoidable, which leads to integrals over the effective neutrino energy at production E P or at detection E D (or both).Here we outline the integral of eq.(3.14) over E P ; the integration over E D instead of E P proceeds in complete analogy.We expand the oscillation phase around E D and keep terms up to linear order in (E P − E D ), using E 0 −E D = (σ 2 eff /σ 2 P,eff )(E P −E D ).But we replace E 0 → E D in the two decoherence terms in the second line of eq.(3.14), ignoring higher order corrections to these terms.Then the integral is of the form discussed in appendix A and we obtain We have dropped the term depending on p P ⊥ in eq.(C.2), as the oscillation phase no longer depends on it, and therefore the integration over p P ⊥ becomes trivial.In contrast, E D -and therefore the oscillation phase -does depend on p D⊥ .In principle, p D⊥ is observable, by reconstructing all particles involved in the detection process.However, this is often not the case in realistic situations.For instance, in reactor neutrino experiments, the momentum of the outgoing neutron cannot be observed.Therefore, we have to integrate also over the neutron momentum phase-space, which effectively means integrating over p D⊥ .Once again we can involve appendix A to derive a corresponding decoherence term.There is no simple closed form for this term, however, its size can be estimated to be of order v 2 D σ 2 pD /σ 2 D,eff relative to the last term in eq.(C.2).One can show that v 2 D σ 2 pD /σ 2 D,eff < 1 and in typical cases it is actually 1.Using the numbers from section 4 we obtain for reactor neutrinos v 2 D σ 2 pD /σ 2 D,eff ∼ 10 −4 .Hence this term can at most induce corrections of order one to the last term in eq.(C.2).

4 νFigure 1 :
Figure 1: Feynman diagram for the total process in an oscillation experiment.

Table 1 :
Spatial localization δ x , momentum spread σ = 1/(2δ x ), and velocity v of the external particles involved in the production (P ) and detection (D) processes (first column) of reactor and Gallium source experiments.

Table 2 :
× 10 −5 5.3 × 10 −6 8.2 × 10 −6 6.4 × 10 −6 √ Σ I 6.4 × 10 −4 3.1 × 10 −4 1.1 × 10 −5 5.1 × 10 −4 Effective velocities and energy-momentum spreads at neutrino production (P ) and detection (D)for reactor and Gallium source experiments.In the lower part we give the spreads relevant for the energy and localization decoherence terms ξ en and ξ loc .For convenience we give values for σ's in eV and δ's in nm; Hence, we obtain an upper bound on |v I |.From the table we see that effective velocities and velocity spreads are all 1, which implies σ eI σ pI , and the σ eI are approximately The emerging decoherence term can be combined with the σ eff -term in eq.(3.14).Using the definition of σ eff in eq.(3.13) we obtain:R D (L, E D ) ∝ dE P |A αβ | 2 ∝ exp i