Wave-packet treatment of reactor neutrino oscillation experiments and its implications on determining the neutrino mass hierarchy

We derive the neutrino flavor transition probabilities with the neutrino treated as a wave packet. The decoherence and dispersion effects from the wave-packet treatment show up as damping and phase-shifting of the plane-wave neutrino oscillation patterns. If the energy uncertainty in the initial neutrino wave packet is larger than around 0.01 of the neutrino energy, the decoherence and dispersion effects would degrade the sensitivity of reactor neutrino experiments to mass hierarchy measurement to lower than 3 σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} confidence level.


Introduction
Much information regarding neutrino mixing have been revealed in the past few decades. However, in most oscillation data analyses, neutrinos are described as plane waves with definite energy and momentum. Since neutrino production and detection are spatially localized, a wave-packet description is more general and appropriate for a complete understanding of neutrino oscillations. Even if the plane-wave treatment is a good approximation for neutrino flavor transitions, the wave-packet decoherence and dispersion effects could still give rise to small corrections to oscillation parameters. We investigate the wave-packet treatment in detail, constrain the energy uncertainty σ wp of reactor anti-neutrinos, and calculate corrections to the mixing parameters by the a e-mail: ylchan87@gmail.com b e-mail: mcchu@phy.cuhk.edu.hk c e-mail: kmtsui@icrr.u-tokyo.ac.jp d e-mail: wongchf@mail.sysu.edu.cn e e-mail: jyxu@phy.cuhk.edu.hk Daya Bay [1] and KamLAND [2] reactor neutrino experiments.
According to our analyses, the wave-packet treatment does not produce significant modifications of the mixing parameters measured by current reactor neutrino experiments based on the plane-wave analysis. However, current experimental data allows a large possible range in the initial momentum width of the neutrino wave packet (σ ν ). If the initial momentum/energy uncertainty of the neutrino wave packet is larger than around 0.02 of the neutrino energy, the decoherence and dispersion effects could have significant effects on future measurements of the neutrino mass hierarchy.
The article is organized as follows: in Sect. 2, the wavepacket treatment is briefly reviewed. In Sect. 3, we use the survival probability derived in Sect. 2 to explore the sensitivity of potential measurements of mass hierarchy with medium baseline reactor experiments. We summarize and conclude in Sect. 4.

Wave-packet treatment for neutrino oscillations
The plane-wave description of neutrino oscillation has been developed for almost 40 years [10]. In the standard calculation of plane-wave neutrino oscillations, the probability of detecting a neutrino flavor state |ν β with energy E, evolved from a pure flavor state |ν α , at a distance z from the production point is given by Fig. 1 An illustration of the decoherence and dispersion effects on neutrino flavor oscillations. The different mass eigenstates (red and blue lines) in a flavor neutrino would gradually separate in space due to speed differences, leading to decoherence and reduction of interference.
At the same time, each mass eigenstate wave packet would disperse, partially compensating for the decoherence effect. However, the overlapping fractions of the wave packets remain smaller than those of plane waves where U αi denote the elements of the PMNS matrix and m 2 k j = m 2 k −m 2 j are the differences of the mass eigenvalues squared [11].
However, as neutrino production and detection are spatially localized, there must be finite intrinsic energy/ momentum uncertainties and a neutrino should be described by a wave packet. The wave-packet character of light has been discussed in detail in Ref. [12]. Based on similar arguments as in that reference, all particles are produced and detected as wave packets. In 1976, the wave-packet nature of propagating neutrinos was proposed [13]. A wave-packet description is expected to be more general and appropriate for a complete understanding of neutrino oscillations [6,8,9,[14][15][16][17]. Nevertheless, there are also arguments against the wavepacket treatment. References [18,19] argue that a wavepacket description is unnecessary as the oscillation system is stationary. However, it has been pointed out in Refs. [20][21][22][23][24] that the authors of Refs. [18,19] have mixed the macroscopic stationarity with microscopic stationarity. The wave-packet description of neutrino oscillations is necessary at least in principle.
Therefore, a neutrino is described by a wave packet as it propagates freely [8,9,15]: where |ν i is an energy eigenstate with energy E i , p ν is the mean momentum, σ ν is the width of the wave packet in momentum space, 1 assumed to be independent of the neutrino energy here, and |ν α is a neutrino flavor state. Figure  1 pictorially describes the wave-packet effects on the propagations of neutrino mass eigenstates. In order to calculate the integral in Eq. (2), we expand the energy E i ( p) around the mean momentum p ν up to second order where is the group velocity of the wave packet. We use Eqs. (2), (3), and (4) to calculate the neutrino flavor transition probabilities at baseline L: Detailed derivation of Eq. (5) is shown in Appendix A. In Eq. (5), the terms in the first bracket correspond to the standard plane-wave oscillation probabilities, and those in the second bracket represent the modifications from the wavepacket treatment. The exp(−λ i j ) term corresponds to the decoherence effect due to the fact that different mass states propagate at different speeds v i ( p ν ) and they gradually separate and stop to interfere with each other, resulting in a damping of oscillations. The terms depending on y i j come from the quadratic correction in Eq. (4); they describe the dispersion effects and are dependent on the dispersion length(s) 2 L dis i j . Furthermore, y i j are proportional to σ 2 wp , while x i j ∝ σ wp only. Therefore, if σ wp 1, the dispersion effect is expected to be more suppressed and negligible. Dispersion has two effects on the oscillations. On the one hand the spreading of the wave packet compensates for the spatial separation of the mass states, hence restoring parts of their interferences. On the other hand, dispersion reduces the overlapping fraction of the wave packets, and thus the interference or oscillation effects cannot be fully restored. Moreover, it also modifies the flavor oscillation phases: 2 The "dispersion length" in this report represents the distance where the dispersion effect becomes important in the neutrino oscillation. A different definition of dispersion length and more detailed discussion of dispersion can be found in [7,17,20].
with deviations from the standard plane-wave oscillation phase written in the parentheses. If y i j = 0, then φ i j just reduce to the standard plane-wave oscillation phases. Additional discussions as regards the phenomenological consequences of the wave-packet treatment and details of the derivation of Eq. (5) can be found in [5,7,8,20].
In this paper, we focus on the analyses of reactor neutrino experiments. According to our wave-packet treatment, in reactor neutrino experiments, the anti-electron neutrino survival probability is

The constraints from current reactor experiments
To date, the value of the parameter σ wp has not yet been determined. If σ wp is not negligible, the wave-packet effects could be significant and have important implications on current and future neutrino oscillation experiments. In this article, we constrain σ wp by analyzing the published Daya Bay and KamLAND data shown in references [25] and [2], considering statistical errors only. Figures 2 and 3 show the data points from Daya Bay experiment [25], along with the oscillation curves corresponding to different values of σ wp . Figure  4 shows the data points from KamLAND [2] and the oscillation curves of plane-wave and wave-packet treatments. In Figs. 2 to 4, L effective are the flux-weighted average reactor baselines, 3 and the error bars just show the statistical uncertainties. We use Eq. (7) to fit the data points in Figs. 2 and 3 to get the constraint of σ wp from the Daya Bay Experiment. In these figures, the black (solid) curves correspond to σ wp = 0, sin 2 2θ 13 = 0.084 and m 2 32 = 2.43 × 10 −3 eV 2 , the same as those of the standard plane-wave treatment. The brown (dashed) curve represents the wave-packet result with σ wp = 0.1, sin 2 2θ 13 = 0.084 and m 2 32 = 2.4 × 10 −3 eV 2 (the bestfit mixing parameters when σ wp = 0.1). The green (dashed) curve shows the oscillation pattern with σ wp = 0.3, sin 2 2θ 13 = 0.096 and m 2 32 = 2.22 × 10 −3 eV 2 (the best-fit mixing parameters when σ wp = 0.3). The blue (dot-dashed) curve corresponds to σ wp = 0.5, sin 2 2θ 13 = 0.108 and m 2 32 = 2.10 × 10 −3 eV 2 (the best-fit mixing parameters when σ wp = 0.5). The data points in Figs. 2 and 3 show clearly the existence of neutrino oscillation, and the data agree with the Pēē derived by the plane-wave approach (Eq. (1)) for a certain set of mixing parameters. This is not surprising since the baselines of the Daya Bay Experiment are short compared to the coherence length, so that the oscillations are not washed out yet.
Moreover, the decoherence and dispersion effects are dependent on the baseline L, and thus the wave-packet impact in the far hall of Daya Bay is expected to be more significant than in the near halls. In the plots in Fig. 3, the black (solid) and brown (dashed) curves overlap almost completely, while the differences between the black, green and blue curves are also small compared to the error bars. However, for the far hall EH3, the wave-packet impact becomes significant for σ wp ≥ 0.3. It is because in Eq. (7), the term λ i j depends on the baseline L. Since the effective baseline of EH3 is longer, the damping of oscillation (decoherence effect) in EH3 is more significant.
The result of our data analysis is shown in Fig. 5. Here, we have just considered the statistical errors. The constraints on the parameters could become worse with the systematic errors taken into account. Figure 5 shows that the wavepacket impact is not significant in Daya Bay experiment and it hardly affects the measurement of θ 13 . The vertical black line in the figure represents the best-fit value of sin 2 2θ 13 in the plane-wave model. A larger value of σ wp , or larger decoherence effect, implies that the true value of θ 13 should be larger. 4 From our analysis, there is no strong evidence to suggest non-zero σ wp . Moreover, our result suggests that the modification of θ 13 is not significant even when the wave-packet framework is considered. The best-fit value of sin 2 2θ 13 from plane-wave analysis (vertical black line) is not ruled out even with 1 σ C.L. We believe that it is because the effective baselines of the Daya Bay Experiment are short compared to the coherence length.
We perform a similar analysis with KamLAND data from Fig. 4. The result of our data analysis is shown in Fig. 6. Again, the systematic errors are not considered.
Equation (5) shows that L coh /L dis = 2σ wp , which implies that the dispersion length is much longer than coherence length if σ wp is negligible. This means that the decoherence  The 1 σ (solid), 2 σ (dashed) and 3 σ (dot-dashed) constraints on "σ wp vs. sin 2 2θ 13 " from fitting the published Daya Bay data [25]. Only statistical errors are considered. We have marginalized over m 2 31 in this plot. The vertical black line represents the best-fit sin 2 2θ 13 based on the plane-wave analysis effect from separation of the wave packets is expected to be more significant than the dispersion effect. However, the dispersion effect is also important because it partly restores the oscillation, especially in the case of large σ wp . If the dispersion effect is not considered, the modifications of the true values of mixing angles in Figs. 5 and 6 would be more significant. The bounds in Figs. 5 and 6 come from the combination of decoherence and dispersion effects, but the contribution from decoherence is expected to be dominant. The 1 σ (solid), 2 σ (dashed) and 3 σ (dot-dashed) allowed regions on "σ wp vs sin 2 2θ 12 " from fitting the published KamLAND data [2]. Only statistical errors are considered. m 2 21 has been marginalized in this plot. The vertical black line represents the best-fit sin 2 2θ 12 from the plane-wave analysis Figures 5 and 6 show that wave-packet effects are not significant in the current reactor neutrino experiments. The 1 σ upper bound of the energy uncertainty σ wp ∼ O(0.1), which is larger than some previous theoretical estimations (σ wp ∼ O(10 −7 )) [6,8,26,27]. Although our analyses have not considered systematic errors, our study of the current reactor experiments suggests that σ wp can be around O(0.1) for reactor experiments. We emphasize that σ wp in this article is an effective parameter which include both the production and the detection neutrino energy uncertainties.
The estimation of the value of this parameter or the size of neutrino wave packet has not come to a strong conclusion yet. There have been some estimations of the sizes of neutrino wave packets produced in different situations. For example, Ref. [28] uses the pion decay length to estimate the width of neutrino wave packet in the MINOS experiment and argues that there could be significant decoherence effect in the active to sterile neutrino oscillation in MINOS. 5 Regarding to reactor neutrino experiments, Ref. [26] provides an estimation based on the mean free path and mean thermal velocity of the production process and suggests that σ x ∼ 10 −6 m, which implies that σ wp ∼ 10 −7 in reactor neutrino experiments. Meanwhile, Refs. [29,30] suggest that the neutrino emission process is expected to be localized at the scale of inter-atomic distance, and so σ x 10 −10 m, implying σ wp ∼ 10 −3 or even larger. If, on the other hand, one takes the uncertainty of a nucleon's position in a nucleus as σ x , then σ wp could be much larger, even of order 1. The estimated sizes of neutrino wave packets from different approaches could be different by a few orders of magnitude. Moreover, as pointed out by Ref. [20], the relation between the decay time of the source and the wave-packet size of the oscillating particle is not direct. The decay time only puts an upper bound on the wave packet length. There is still no experimental support for such an assumption. In this paper we do not calculate or suggest the theoretical value of σ wp .
We point out that the wave-packet impact is not significant for current reactor neutrino experiments. Nevertheless, the 1 σ C.L. allowed range of σ wp is O(0.1), within which the potential wave-packet impact could lead to significant effects and additional challenges in future neutrino oscillation experiments.

Measurement of neutrino mass hierarchy
The signs for m 2 31 and m 2 32 have not yet been determined. Normal hierarchy (NH) corresponds to positive m 2 31 and m 2 32 with ν 1 as the lightest mass state. Inverted hierarchy (IH) corresponds to negative m 2 31 and m 2 32 with ν 3 as the lightest mass state [31].
As indicated by the recent data obtained by the Daya Bay Experiment, sin 2 2θ 13 = 0.084 ± 0.005 [25]. With such a (relatively) large value of θ 13 , it is possible to determine the neutrino mass hierarchy (MH) in reactor neutrino experiments at medium baseline [32,33].
For a detector at baseline L, the observed anti-electron neutrino flux of visible energy E vis is given by φ(E) is the reactor neutrino energy spectrum, and σ (E) is the inverse beta decay cross section. R is the detector response function with energy resolution δ E, which will be discussed in more detail in the next subsection. It is well known that at a baseline of around 50 km, which corresponds to the first minimum of θ 12 oscillation for reactor neutrinos, the sensitivity for measuring MH is maximal. The upper panel in Fig. 7 shows f (E vis ) at L = 53 km for NH and IH, with standard oscillation parameters in the plane-wave treatment (σ wp = 0 in Eq. (7)).

Impact of wave-packet treatment
Our wave-packet treatment shows that the amplitudes of neutrino oscillations will be reduced. In particular, damping of the θ 13 oscillations will be significant for an intermediate baseline reactor neutrino experiment if σ wp is O(0.1). The lower panel in Fig. 7 shows the neutrino visible energy spectrum at a baseline of 53 km, with standard neutrino mixing parameters and σ wp = 0.1. If σ wp is large, the neutrino spectra for NH and IH are indistinguishable from each other.
We modify the GLoBES software [34,35] to perform numerical simulations of a 53 km baseline reactor neutrino experiment, using a similar setup as in [36] and [33]: a 20 kton detector with 3 % energy resolution, reactors with a total thermal power of 40 GW and a nominal running time of six years. In the absence of oscillations, the total number of events is about 10 6 . As this paper focuses on the wave-packet impact, the systematic errors were not taken into account in the following simulations. 6 We took the oscillation parameter values from a global analysis [37] as m 2 21 = 7.54 × 10 −5 eV 2 ( m 2 31 + m 2 32 )/2 = 2.43 × 10 −3 eV 2 , sin 2 θ 12 = 0.307, and sin 2 θ 13 = 0.0241.
To distinguish between NH and IH, we quantify the sensitivity of the MH measurement by employing the least-squares method, based on a χ 2 function: where N M i is the measured neutrino events in the ith energy bin, and N T i is the predicted number of neutrino events with oscillations taken into account (without considering the Fig. 7 Upper panel Expected neutrino visible energy spectrum for plane-wave neutrino oscillation at 53 km for NH (solid) and IH (dashed). Lower panel Same as in upper panel, but with σ wp = 0.1. The two curves for NH and IH completely overlap due to the wave-packet impacts systematic errors). 7 The number of bins used N bins is 164, equally spaced between 1.8 and 10 MeV.
We fit the hypothetical data set with m 2 ee as the free variable, defined as m 2 ee = cos 2 θ 12 m 2 31 + sin 2 θ 12 m 2 32 .
The capability to resolve the mass hierarchy is then given by the difference between the minimum χ 2 value for IH and NH: χ 2 MH is used to explore the wave-packet effects as well as the impact of statistics and systematics in measuring the MH. 8 If σ wp = 0, which corresponds to the plane-wave treat- 7 Without loss of generality, in our simulations, NH is assumed to be the true mass hierarchy. The result is identical if we assume IH to be true. 8 Note that min(χ 2 NH ) and min(χ 2 IH ) can be located at different values of | m 2 ee |. ment, we get χ 2 MH ≈ 19.5, implying that we could distinguish the MH with a confidence level of nearly 4 σ , as shown in the upper panel of Fig. 8. However, if σ wp = 0.02, the sensitivity will be reduced due to the damping of oscillations and χ 2 MH will drop to around 3.35, as shown in the bottom panel of Fig. 8.
The solid (black) curve in Fig. 9 further shows the variation of χ 2 MH as a function of σ wp . It shows that χ 2 MH drops rapidly with σ wp , to become smaller than 3 σ C.L. as σ wp > 0.012. In this case, it is difficult to determine MH.
In the rest of this section we will investigate how to increase the sensitivity to MH by improving the experimental setup.

Consequences of energy resolution and statistics
References [32,38,39] suggest that the sensitivity of MH measurement depends on the detector resolution. As shown in Eq. (8), the observedν e flux depends on the detector response function R and energy resolution δ E , which are defined by Fig. 8 The sensitivities of measuring m 2 ee with and without wave-packet effects. NH is assumed to be the true hierarchy in the simulation. The solid (dashed) curve corresponds to the fitting of NH (IH, false hierarchy). σ wp is assumed to be 0 in the upper panel and 0.02 in the lower panel where the detector energy resolution is parameterized as In the previous subsection, we assumed that a = 0.03 and b = 0 in order to achieve a 3 % detector resolution. As mentioned in the footnote in Sect. 2, the detector resolution is different from the energy uncertainty of detection at the microscopic level. The detector resolution is determined by the properties of the macroscopic detector, which should be taken into account even in the plane-wave assumption. Similar to the decoherence effect due to σ wp , a poor energy resolution can also destroy the measured oscillation effect. At L = 53 km, the θ 13 oscillation could be smeared out with a finite energy resolution, and the MH information could be destroyed, particularly in the case of large σ wp . Figure 10 shows 9 that if σ wp > 0.02, which is allowed by the Daya Bay and KamLAND data, the detector resolution (or the parameter a in Eq. 13), must be better than 3 % in order to determine the MH with a C.L. of more than 2 σ .
On the other hand, we can also improve the MH sensitivity by collecting more data. Figure 11 shows the impact of statistics on the measurement of MH as a function of σ wp , assuming a detector energy resolution of 3 %, as suggested Fig. 11 The number of neutrino events (no oscillation) required to resolve MH at 2 σ (dotted), 3 σ (dashed) and 4 σ (solid) C.L., as a function of σ wp . Other details of the experimental setup are described in the text by Refs. [33,38]. Much longer run time would be required for a 2 σ C.L. measurement if σ wp is larger than 0.01.

The optimal baseline
Without considering the wave-packet impact, 50-60 km is the ideal location to measure mass hierarchy for reactor neutrino experiments [33,36]. However, in the presence of significant wave-packet impact, longer baseline would lead to larger damping of the oscillation amplitude. Although reducing the baseline will lead to a loss of maximal phase difference between the NH and IH oscillation curves, this can save back part of the oscillation. Therefore, the optimal baseline of measuring MH could be shorter than 50 km, depending on the value of σ wp . Figure 12 shows the χ 2 MH as a function of baseline for different values of σ wp . In the case of the plane-wave limit (σ wp = 0), the MH can be distinguished with a confidence level of nearly 4 σ at 53 km. However, as σ wp increases, the maximum χ 2 MH shifts to shorter baseline. Figure 13 further shows the value of optimal baseline as a function of σ wp .

Conclusion
The wave-packet impacts on current reactor neutrino oscillation and neutrino mass hierarchy measurements have been discussed. Our analyses show that the wave-packet treatment would not lead to significant modifications of the oscillation parameters from the Daya Bay and KamLAND results based on plane-wave assumptions. Moreover, our analyses The decoherence and dispersion effects depend not only on the initial neutrino energy uncertainty, but also on the values of m 2 and baseline. Since the measurement of the neutrino MH in medium baseline reactor neutrino experiments relies on the fast m 2 32 oscillations, the decoherence and dispersion effects could be significant and make it more difficult. We found that even if σ wp is just around 0.02, the sensitivity of MH measurement would be largely reduced. The optimal baseline shifts to smaller value as σ wp increases, due to the damping of oscillation amplitudes by wave-packet effects.
We have to emphasize that we are not suggesting that the wave-packet impact would be so large as to make the mass hierarchy measurement in medium baseline reactor neutrino experiments impossible. We point out that the planewave model of neutrino oscillation is only an approximation, and the wave-packet treatment is more general. While the wave-packet effect is insignificant in the current reactor neutrino oscillation experiments, its impact on future oscillation experiments needs to be determined.
which is a moving Gaussian wave packet with dispersion.
Here, we have used the approximation Then the transition probability of ν α → ν β at a distance L and time t is given by which is a function of both time (t) and distance (L).
In an oscillation experiment, the neutrino is detected at a fixed baseline L but the time is not measured. In order to obtain the oscillation probability as a function of the baseline, the time has to be integrated out in Eq. (A.4). Since reactor neutrinos propagate almost at the speed of light, P αβ (L , t) is non-zero only around t ∼ L. The transition probability is non-zero only within a small time window t around a time On the other hand, the size of t is constrained by the spatial width of the wave packet which is typically much smaller than the baseline L, which means that t t L . Moreover, Equation (A.6) shows that the factor (t/τ i ) changes slowly with the variable t. Within a small t, these terms can be treated as constants. Therefore, within the small integration region which is constrained by the width of the wave packet, it is acceptable to approximate t = t L ∼ z for this factor, However, since Eq. (A.7) shows that the derivative of (z − v i ( p ν )t) is not negligible even in a small region. We did not use the same approximation in the factors . Therefore, the integral of Eq. (A.4) can be approximated by In the last step of Eq. (A.9), we have ignored the terms proportional to m 2 i m 2 j E( p ν ) 2 , since they are expected to be negligible for ultrarelativistic neutrinos. In the case of reactor neutrino experiments, theν e survival probability is given by The damping factor exp(−γ i j ) corresponds to an decoherence effect from delocalization, which is neglected in the main article. This is because in most circumstances, γ i j ≈ 0. The details of γ i j and the delocalization decoherence effect will be discussed in Appendix B.

Appendix B: The decoherence effect from delocalization
The decoherence effect discussed in Sect. 2 is due to the separation of different neutrino wave packets. With larger values of σ wp , the decoherence effect would be more significant. On the other hand, there exists another decoherence effect which is due to the delocalization of the production and detection processes. Different from what we have studied above, the decoherence effect from delocalization will become significant only if σ wp is extremely small, or the spatial uncertainty σ x is large. In fact, in neutrino oscillations, one of the coherence conditions is that the intrinsic production (and also detection) energy uncertainties are much larger than the energy differences between different mass eigenstates E i j = E i − E j (E i is the energy of mass eigenstate |ν i ) [40], namely, Equation (B.11) implies that in order to measure the interferences between different mass eigenstates, the spatial uncertainty σ x , has to be much smaller than the oscillation length L osc i j . Equation (5) is just an approximate neutrino oscillation probability formula and it does not describe the decoherence effect from delocalization. More precisely, the flavor transition probability in Eq. (5) should be multiplied by an additional factor exp(−γ i j ), (B.12) where With these delocalization terms also taken into account, a more completeν e survival probability should be written as Eq. (A.10) in the previous section. The damping factor in Eq. (B.12) corresponds to the decoherence effect from delocalization. If γ i j ∼ 0, the modifications from delocalization are negligible. In this case Eq. (A.10) just reduces to Eq. (5).
γ i j is proportional to σ 2 x (L osc i j ) 2 , and so the decoherence effect from delocalization matters only when σ x ∼ L osc i j , about 1 km in the measurement of θ 13 oscillation. In this case σ wp is extremely small.
If we use Eq. (A.10) rather than Eq. (5) to do the analysis, we will find that the γ i j term will offer a lower bound on σ wp and the delocalization effect is significant only when σ wp is extremely small, which implies a large σ x . Figure 14 shows our analysis of the Daya Bay data, which is similar to Fig. 5, but this time the delocalization term is considered.
From Fig. 14 we can see that in a large range of the parameter space, there are no modifications on the true value of sin 2 2θ 13 . This means that in most of the parameter region, wave-packet impact can be safely neglected. Moreover, Fig. 14 also suggests that only if σ wp ∼ O(10 −16 ), which means that σ x ∼ O(1 km), the decoherence effect from delocalization is significant.
At this point, we can discuss the wave-packet impact in two different regimes. If σ wp is large (∼ O(10 −1 )), since λ i j ∝ σ 2 wp , the damping factors exp(−λ i j ) in Eq. (5) become significant and the decoherence effect from the separation of the wave packets cannot be neglected. On the other hand, if σ wp ∼ O(10 −16 ) or even smaller, then the additional wp . In this case the decoherence effect from delocalization becomes important. Nevertheless, in most reactor neutrino experiments, the dimensions of the reactor cores and detectors are just around a few meters. It is unlikely that the spatial width of the initial neutrino wave packet would be larger than 1 km.
However, in most reactor experiments, including current ones such as Daya Bay and KamLAND, and also the proposed measurements of neutrino mass hierarchy at medium baseline, the delocalization terms exp(−γ i j ) can be neglected. This is because π 2 (L osc i j ) 2 E 2 σ 2 wp 1 (or, the spatial uncertainty σ x L osc i j ). In this case, Eq. (B.12) just reduces to Eq. (5). Therefore, we neglected the decoherence effect from delocalization in our study and performed the simulations and analyses with Eq. (5).