Entangled Neutrino States in a Toy Model QFT

It has been claimed that wave packets must be covariant and also that decohered neutrino oscillations are always revived during measurement. These conjectures are supported by general arguments which are not specific to the electroweak theory, and so if they are true for neutrinos they will also be true for simplified models. In this paper we produce such a simplified model in which the neutrino wave function, including its entanglement with the source particle and the environment, can be calculated explicitly in quantum field theory. It exhibits neutrino oscillation, which is reduced at late times by decoherence due to interactions of the source with the environment. One simple lesson from this model is that only the difference between the environmental interactions before and after neutrino emission can reduce the amplitude of neutrino oscillations. The model will be used to test the conjectures in a companion paper.


Introduction
Reactor neutrino experiments report lower values of θ 13 than accelerator experiments. It is customary to reduce this tension by assuming the normal hierarchy and a value of the CP-violating phase δ near 270 • . This increases the expected appearance signal at accelerator experiments, allowing the small θ 13 mixing reported by reactor experiments to produce almost as many electron (anti)neutrinos as are observed at muon (anti)neutrino beams. But there is another logically consistent possibility. The reactor neutrinos have lower energy, and so are expected to be more prone to decoherence than accelerator neutrinos [1]. Indeed no decoherence is expected in the case of accelerator neutrinos [2]. In this case the reactor neutrino measurement of θ 13 , based on an analysis with no decoherence, is underestimated and the evidence for the normal hierarchy and maximal CP-violation is weakened. Furthermore, the degradation of the signal observed by JUNO would be considerable [3]. This possibility has been rejected by the Daya Bay collaboration [4]. However their study relied upon a neutrino wave packet model.

Wave Packet Models of Neutrinos
The traditional view of decoherence in neutrino oscillations comes from the quantum mechanical wave packet model. Here neutrinos are produced as a flavor eigenstate wave packet, localized in space and time. The lighter mass eigenstate travels faster than the others and so the wave packets corresponding to different mass eigenstates spatially separate after travelling a distance called the coherence length. This separation leads to decoherence and therefore a decrease in amplitude of neutrino oscillations. The spatially separated mass eigenstates may nonetheless be coherently summed by the detector if the detector has a sufficiently long coherence time, leading to a restoration of neutrino oscillations [5]. The coherence length clearly depends on the spatial size of the wave packet, which is a parameter in such models. It has long been recognized [6] that this spatial size is determined by interactions of the neutrino source particles with the environment. Usually order of magnitude arguments are used to estimate this parameter [6,7,8,9], and the result is substituted into the model.
In quantum mechanics, neutrino wave packets are created by hand. In quantum field theory (QFT) they are created consistently from electroweak interactions. Consistent QFT treatments necessarily create neutrinos entangled to their source particles, such as unstable nuclei or mesons, and also to charged leptons which are created simultaneously. We will refer to all of the particles involved in the interaction which produced the neutrino as source particles, including the charged leptons. Again in this case the environment plays a role. As noted, for example, in Ref. [10] the interactions of the source particles with the environment disentangle the neutrino from the rest of the state and so allow its treatment as a wave packet. This disentanglement is caused by environmental interactions which effectively measure the source particles [11]. It is customary in QFT treatments to apply this interaction by simply projecting the entangled state onto a subsector of the Hilbert space in which the source particles have some definite position or momentum wave function, as if they were actually measured. With the positions of the source particles specified, one can determine a space time region in which the neutrino is created and so the neutrino is again in a localized, flavor eigenstate wave packet. Now, just like the quantum mechanical case, the different mass eigenstates travel at different speeds and so separate, leading to decoherence.
Quite a different QFT treatment appeared in Ref. [12]. Here the different neutrino mass eigenstates were not forcibly created in the same time window. Of course modern neutrino experiments measure neutrinos in a fixed time window, in flavor eigenstates. Therefore the fact that lighter neutrinos travel faster and the travel distance is fixed implies that the lighter mass eigenstates are emitted after the heavy mass eigenstates. So instead of wave packet separation, here the wave packets coelesce, and no decoherence was reported by the authors.
How could QFT produce two such phenomenologically distinct paradigms? In the first case, environmental interactions were imposed by hand, with a simple projection. In the second case, environmental interactions were not included at all.

Wave Packets from Entanglement
It is our goal to understand when the wave packet treatment of neutrinos is and is not reliable, and to understand how to calculate the wave packet size. We will do this via a first principles, consistent calculation in QFT. Papers on QFT treatments of neutrinos generally calculate the S matrix for neutrino creation and detection, which is the amplitude for the creation of a given state in the asymptotic future, long after the neutrino has been absorbed. However we are interested in the state of the neutrino itself, and so are interested in intermediate states. Such information can not be directly obtained from the S matrix. It is accessible in the Schrodinger picture of QFT, in which operators are time-independent and states evolve via the action of the Hamiltonian operator. An experiment begins with a source state entangled with the environment and the Hamiltonian evolves this initial state into the future. This evolution creates neutrinos.
As was noted in Ref. [12], it is true that different neutrino mass eigenstates may be created at different times. Indeed, evolving the state of a 235 U nucleus for one year in the Schrodinger picture, neutrinos may be emitted at any time during the year and so the neutrino wave function extends for one light year. It is certainly not a localized wave packet. In the calculation of matrix elements, one must sum over each mass eigenstate and separately integrate the interaction times over the entire year. Now the key question is, whether at a fixed time the different mass eigenstates contribute coherently to matrix elements. If they do, one expects to observe neutrino oscillations, if they do not, these oscillations will be damped. Measurements occur in a flavor basis and, in modern experiments, at a reasonably well-determined time. Therefore contributions to the relevant amplitudes come from states in which the different mass eigenstates are localized in space time at detection, meaning that the lighter neutrino was emitted later, again in agreement with [12].
However for a coherent summation of neutrino mass eigenstates it is not sufficient that they spatially overlap. The entire final states must agree, including the source particles and the environment. In other words, if the state is where |E i are the environment plus source particles part of the state and |ν i are the neutrino mass eigenstates, then the summation is fully coherent only if the |E i are equal up to a phase. This condition is the origin of decoherence. The fact that the lighter neutrino was emitted later means that the source particles interacted differently with the environment, for example the unstable particle had more time to interact while the product particles had less. This necessarily implies that the environment part of the state will be different in the case of each mass eigenstate. The bigger the difference in mass or the further the neutrino has traveled, the bigger the difference in time between the emissions of the different mass eigenstates and so the bigger the decoherence.
The conclusion is that while Ref. [12] is correct that the times of the emissions of the various mass eigenstates need not agree, nonetheless if the difference exceeds some threshold then coherence will be lost. We claim that this threshold should be interpreted as the wave packet size in the wave packet model. In this case, decoherence will correspond to the spatial separation of the wave packets. However it is not obvious that long measurements may now restore coherence as in Ref. [5].

Our Approach
For the questions of interest, concerning neutrino oscillations, wave packets, and decoherence, the details of the electroweak interactions do not play any essential role. Therefore, we will work in the simplest toy model which has the features of interest, a scalar field theory in 1+1 dimensions. Here we can, in the Schrodinger picture of QFT, numerically evolve the full entangled state to any desired moment in time to understand it. Thus our approach is similar to that of Ref. [13] but including environmental interactions. To simplify the situation yet further, we will not consider measurements of the neutrinos. Therefore our final states will be the neutrinos themselves and we will calculate transition amplitudes and transition probabilities from states with no neutrinos to states with a neutrino. We will see that these probabilities already have a rich phenomenology of oscillations and decoherence. Of course it means that we cannot tell whether coherence can be revived through measurement, however we feel that a robust study of coherence revival via measurement requires a characterization of the coherence before measurement, which our method provides. We do not model interactions with the environment by projecting on to a definite state for the environment and the source particles. Instead all particles are consistently evolved in the Schrodinger picture of QFT. In the calculation of probabilities, the distinct environment and source final states are incoherently summed.
The phenomenology of wave packet models includes several potentially interesting effects, such as the revival of oscillations ruined by docoherence via long measurements in Ref. [5]. In [14] it was asserted that, presumably as a result of revival, decoherence is unobservable in neutrino oscillation experiments. Another claim [15,16] is that neutrino wave functions are always "covariant wavepackets." This means that they depend on the momentum only via Lorentz scalars. The covariant wave packet hypothesis was assumed in the experimental analysis of decoherence at Daya Bay [4]. We believe that our QFT approach will allow a robust test of these claims.
Our study has three advantages over most quantum field theory (QFT) approaches to neutrino oscillations and decoherence. First, we calculate the full, entangled state consisting of the source, the neutrinos and the environment 1 at arbitrary times and not just the asymptotic S-matrix. This will allow a robust test of the covariant wave packet proposal. Second, we explicitly consider interactions between the source and the environment 2 . Third, we integrate our transition probability over the possible final states of the source and the environment. It is this integration which leads to decoherence, reducing the amplitude of neutrino oscillations in the transition probability.
Perhaps one of the most serious attempts at the determination of the wave packet size, in the case of solar neutrinos, was Ref. [6]. Unlike later estimates, it includes an estimate of the phase angle variation resulting from each interaction instead of merely assuming that an interaction automatically results in decoherence. However, in the case of reactor neutrinos, unlike solar neutrinos, the source nuclei are large and so the Coulomb interactions in some cases are hardly affected by a beta decay. We will see in our example that the decoherence is not determined by the total phase induced by an interaction, but rather by the difference in the phase that would be acquired before and after the beta decay. This difference, in the case of reactor neutrinos, may be one or two orders of magnitude smaller than the total phase, and thus the wave packet size may be expected to be an order or magnitude or two 1 The key role played by the entanglement of the neutrino and the source particles in a QFT treatment has been stressed in Ref. [13]. In Ref. [17] it is claimed that the full entangled QFT treatment leads to the same amplitudes as a wave packet treatment. However neither study included interactions of the source with the environment.
2 Such interactions were included in Ref. [18] by including a phenomenological smearing of energies. We instead consistently treat the interactions in QFT.
larger than may be expected by simply adapting the argument of Ref. [6] to the case of reactor neutrinos. This is one immediate lesson that may be drawn from our simple model.
We begin in Sec. 2 with a simplified model in which the neutrinos are created from a classical source. This model exhibits oscillations. However the neutrinos are always off-shell and also, because the source is classical, it cannot be entangled with the environment and so there is no decoherence. Next in Sec. 3 we introduce our full model. We include both source fields and also environment states. Our analysis of this model is presented in Sec. 4.

The Model, Fields and States
We do not believe that spin plays a key role in a qualitative understanding of decoherence in neutrino oscillations. Therefore our model will involve only real scalar fields. Similarly, we will restrict our attention to one space and one time dimension. So long as our fields are massive, this assumption leads to only a modest reduction in computational complexity. Finally, as our most significant assumption, we will consider one-body and two-body decays instead of three-body decays. Therefore the scalar fields which we will call "neutrinos" will carry no conserved lepton charge. Nonetheless we will introduce two flavors of neutrinos, so that there will be oscillations.
The neutrinos in our model are described by the canonical real scalar fields where the index i labels the mass eigenstates ψ 1 and ψ 2 . The conjugate momenta are We always work in the Schrodinger picture, so all operators such as fields and their conjugate momenta are time-independent.
The Hamiltonian will be decomposed into a free and interaction term where H 0 is the free real scalar field Hamiltonian density 3 The interaction Hamiltonian describes neutrino creation by a classical source of size 1/(2 √ α) Observe that neutrinos are created not in a mass eigenstate ψ i , but rather in the superposition ψ 1 + ψ 2 which plays the role of a flavor eigenstate in our model.
Let |Ω and |i, p be respectively the ground state and one neutrino states of the free The states |i, p provide an orthogonal basis for the 1-particle states.
In practice one is interested in the measurement of a neutrino at a particular position x. While it is straightforward to define an orthogonal position basis for the 1-particle states, this does not reflect the basis in which neutrinos are usually measured in modern experiments. Usually one measures both a neutrino's momentum and also position. Clearly the uncertainty principle implies that these are each measured with a finite resolution. Let σ be the momentum resolution of a given detector. For simplicity, we will consider a detector which is only sensitive to neutrinos of momentum p 0 , although this can easily be generalized to a multichannel detector. Then the relevant basis of 1-neutrino states will be Note that while these states do form a basis for the 1-neutrino sector of the Hilbert space, they are not orthogonal

Evolution
To calculate the evolution of this system, we will need to know how the Hamiltonian acts on the various states. In terrestrial neutrino experiments, multineutrino processes are too suppressed to be relevant. Thus we will be interested only in evolution involving a single power of H I and only in 0-neutrino and 1-neutrino states. The action of the Hamiltonian on such states is easily calculated H I |i, p will not arise in the calculation below at first order in H I .
Projecting onto the 1-neutrino sector of the Hilbert space, we then find the evolution of the ground state to an arbitrary time t where the ⊃ on the first line is the restriction to terms with precisely one power of H I . We have chosen to omit terms with no powers of H I from (2.10) as they would not contribute to the matrix elements calculated below and, perhaps more to the point, they contain no neutrinos as so do not contribute to the neutrino wave packet. According to the general arguments in Refs. [15,16], one may identify the state |t with a neutrino wave packet and expect that it is a covariant function of the four-momentum p. No such covariance is manifest in Eq. (2.10). In a sequel, we will investigate whether the wave packets in our models possess the covariance property demanded in these references and assumed by the Daya Bay collaboration in their analysis [4].
Note that we have not explicitly introduced the time t 0 when the neutrino is created. However we may rewrite |t as an integral over t 0 In this note we will not explicitly consider the measurements of neutrinos in our model, these will be included in future work. Our goal for now is to understand neutrino wave ψ 1 H I ψ 2 H I Figure 1: The amplitudes A 1 and A 2 correspond to a neutrino which is created from the ground state of the free Hamiltonian H 0 .
The corresponding probability would be where λ is a normalization constant. However since x is continuous one expects that proba bility of finding a neutrino at any given x to vanish, implying that λ = 0. For a continuous x one is interested instead in the probability density dP(x, t)/dx.
It is therefore tempting to identify (2.14) with a probability density where to cancel the normalization in (2.8). This is not quite right, due to the fact that the |i, x basis is not orthogonal and so if the neutrino is observed at x it has a nonzero probability to also be observed at y. Therefore one cannot define a normalized probability distribution function (PDF) for x. However, such double-valued position probabilities are exponentially suppressed at distances larger than the de Broglie wavelength corresponding to the momen tum resolution. The position resolution of any neutrino detector is much larger than this distance, and so for all practical purposes (2.14) is a PDF. functions. These are already sufficient for constructing amplitudes and probabilities which will eventually be related to measurements in our companion paper. We will be interested in the following amplitude, which corresponds to a transition to a neutrino at a position x at time t This amplitude, when σ = ∞, is the wave function of a single neutrino at time t.
Neutrinos are created in the flavor basis ψ 1 + ψ 2 . Of course, it reality one also measures them in the flavor basis. While we do not consider the measurement here, this does motivate us to introduce the flavor basis matrix element (shown in Fig. 1) (2.13) One can also define a transition probability from the H 0 ground state to a one neutrino state. This is not the probability of a measurement, since there is no term in our Hamiltonian which measures a neutrino. It is simply the probability that a neutrino exists at time t and position x, given that the system began in the H 0 ground state at time t = 0. Naively the transition probability would be where λ is a normalization constant. However since x is continuous one expects that probability of finding a neutrino at any given x to vanish, implying that λ = 0. For a continuous x one is interested instead in the probability density dP (x, t)/dx.
It is therefore tempting to identify (2.14) with a probability density where to cancel the normalization in (2.8). This is not quite right, due to the fact that the |i, x basis is not orthogonal and so if the neutrino is observed at x it has a nonzero probability to also be observed at y. Therefore one cannot define a normalized probability distribution function (PDF) for x. However, such double-valued position probabilities are exponentially suppressed at distances larger than the de Broglie wavelength corresponding to the momentum resolution. The position resolution of any neutrino detector is much larger than this distance, and so for all practical purposes (2.14) is a PDF.
Summarizing, we have argued that the transition probability density for the creation of a neutrino at (x, t) in the flavor basis is This can again be written as an integral over the interaction time t 0 The largest contribution to the integral over p comes from the stationary point of the phase which yields the usual condition that the group velocity is equal to the average velocity of the neutrino in the time t − t 0 since its creation.

Numerical Results
We will now consider the case corresponding to classical source of width 0.5, a measured neutrino momentum of 1 ± 0.1 and neutrino masses of 0.3 and 0.4.
The amplitudes A i (x, 100) defined in Eq. (2.12) are shown in Fig. 2. Three peculiar features are evident in the left panel. First, the maximum amplitude occurs near x = 0. This is a consequence of the fact that the initial energy of the system is equal to zero, since H 0 annihilates the initial state |Ω . The final energy is therefore also equal to zero, as H is time-independent and so time evolution conserves energy. However the neutrinos are massive, and so they will always be off-shell. This is reflected in the ω in the denominator, which vanishes only if ω = 0, as is never the case. The smallest ω is the least off-shell, and therefore the highest amplitude. As a result the highest amplitude arises for the neutrinos with the smallest momentum, which cannot travel far.
The second and least physical peculiar feature is the peak near x = t corresponding to neutrinos created at t 0 = 0. Recall from Eq. (2.11) that one integrates over t 0 , and so why should most of the neutrinos observed arise from t 0 ∼ 0? This is another consequence of the fact that the neutrinos are off-shell. As ω = 0, the phase e −iωt in Eq. (2.11) always oscillates, damping the integral and so the amplitude. This damping is reduced at t 0 = 0 just because this is a boundary of the domain of integration, and so there is no oscillation at t 0 < 0.
The third peculiar feature is the small tail at x > t. One may attribute this tail to the finite size 1/(2 √ α) of the classical source. However the tail is too large to be created by this alone. It is also a consequence of the fact that A is essentially the Feynman propagator Ω|ψ(t)ψ(t 0 )|Ω , albeit with some additional factors. Recall that in quantum field theory only the retarded propagator is causal. The causality of the retarded propagator results from the presence of a commutator term − Ω|ψ(t 0 )ψ(t)|Ω . However no such term is present in A. The physical explanation for the lack of causality of the Feynman propagator is that a particle of mass m cannot be kept in a box of size beneath 1/m, and so a leaking of order 1/m is inevitable [19]. Despite the small mass of the neutrino, the length scale 1/m is far smaller than the position resolution of any experiment and so this tail is irrelevant in neutrino physics.  Fig. 2, that the more massive neutrino travels more slowly and so its amplitude is smaller than that of the lighter neutrino near the light cone. Such kinematic damping is far too small to observe at present day neutrino experiments.

The Model
We are interested in decoherence resulting from interactions of the source particle with the environment, together with quantum entanglement between the neutrino, the source and the environment. The source above was classical and so could not be entangled. Therefore, to incorporate decoherence in our model we must introduce quantum source fields φ I and environment fields E α .

The Fields and Their Interactions
In oscillation experiments the neutrinos travel macroscopic distances and so are observed on-shell. While we do not assert that our final states are on-shell, it will be clear from our expressions that off-shell final states will generally provide a small contribution. The simplest on-shell decay in a Lorentz-invariant theory is the decay of a heavy source particle φ H into a slightly lighter yet still heavy particle φ L and our so-called neutrino ψ i , which is actually a Unlike real-world β decay, the neutrinos in our model are created in a two-body process in which φ H decays to φ L and ψ i .
Decoherence requires coupling to environment fields E α . While two fields would be sufficient, we will consider four, indexed by α ∈ [0, 3]. These will interact with φ H via interactions of the form α φ 2 H E 2 α . We will consider a nonrelativistic approximation of this interaction, so that it is of the form of that in Ref. [11]. In this approximation, we simply add a perturbation to the Hamiltonian equal to where N H and N α are the usual particle number operators for the fields ψ H and E α .

The States
We will perform the usual decomposition of the canonical fields φ I (x) = dp 2π where I runs over the indices {H, L}. The decomposition of the environment fields will not be needed due to our nonrelativistic approximation.
We will only be interested in states with one environmental particle E α , one source particle ψ I and zero or one neutrinos φ i . We will not keep track of the momentum or the position of the environmental particle, we will only be interested in its flavor α. Thus a basis of the states of interest may be written |α; I, p; i, q for states with a neutrino of flavor i and momentum q and a source particle of flavor I and momentum p, together with the states |α; I, p which contain no neutrino. The free particle ground states, with an environment field, may be written as simply |α . These are annihilated by all operators a and A and are orthonormal. The normalizations of the other states are fixed by |α; I, p = A † I,p |α , |α; I, p; i, q = A † I,p a † i,q |α . Figure 5: The amplitudes A 1α (left) and A 2α (right) correspond to a neutrino ψ 1 and ψ 2 respectively which is created from a single source particle φ H interacting with an environment in the state |α . If ψ 2 is heavier than ψ 1 , it will be slower at fixed momentum and so, given an observation at a fixed time and position, it is emitted earlier [9]. Therefore φ H has less time to interact with the environment in A 2α than in A 1α .
close to the measured quantities. Since any measurement will be in the flavor basis and not the mass basis, we will sum these matrix elements over i corresponding to a disappearance channel experiment. Also, since a measurement will also measure, with some resolution σ the neutrino momentum q and will find a value q 0 , we will calculate the matrix elements (shown in Fig. 5) (3.7) If σ = ∞ then these will not be precisely orthogonal. However, we will ignore this and define an approximate probability density as where λ is the normalization constant Note that, unlike the wave packet approach [15,16] in which one considers only a single fina state for the source particle, here the final state p is integrated over. This is reasonable as the final state of the source particle is never measured, and as a result the neutrino is never in a localized wave packet.
13 Figure 5: The amplitudes A 1α (left) and A 2α (right) correspond to a neutrino ψ 1 and ψ 2 respectively which is created from a single source particle φ H interacting with an environment in the state |α . If ψ 2 is heavier than ψ 1 , it will be slower at fixed momentum and so, given an observation at a fixed time and position, it is emitted earlier [12]. Therefore φ H has less time to interact with the environment in A 2α than in A 1α .
Our initial condition will consist of a heavy source particle in a Gaussian wave packet where β is a parameter which determines the initial width of the wave packet. This state is normalized such that One could fix the c α so that this is equal to unity, but we will instead leave the c α free and correct for this normalization in our formula for the probability.
The initial state |0 will evolve into states |α; L, p; i, q and so we will be interested in matrix elements of the form α; L, p; i, q|e −iHt |0 where H is the total Hamiltonian and t is the time to which the system evolves. This matrix element is the amplitude, calculated in the Schrodinger picture, for the initial state |0 to evolve into the final state |α; L, p; i, q .
In this paper we will not yet introduce neutrino measurements. However, as our interest does nonetheless lie in measurement, we consider matrix elements which are close to the measured quantities. Since any measurement will be in the flavor basis and not the mass basis, we will sum these matrix elements over i corresponding to a disappearance channel experiment. Also, since a measurement will also measure, with some resolution σ, the neutrino momentum q and will find a value q 0 , we will calculate the matrix elements (shown in Fig. 5) (3.7) If σ = ∞ then these will not be precisely orthogonal. However, we will ignore this and define an approximate probability density as where λ is the normalization constant Note that, unlike the wave packet approach [20,10] in which one considers only a single final state for the source particle, here the final state p is integrated over. This is reasonable as the final state of the source particle is never measured, and as a result the neutrino is never in a localized wave packet.
These amplitudes and probability densities correspond to transitions from the heavy particle to the light particle plus a neutrino. These are the usual transition amplitudes and transition probabilities in quantum field theory. These are not equal to the amplitudes or probabilities for neutrino measurement, which would require an additional interaction in which the neutrino is absorbed. Nonetheless, these amplitudes and probability densities are interesting because they already manifest neutrino oscillations and decoherence and therefore provide a simple setting in which these pheneomena may be studied.

Analytical Calculation
As events involving multiple neutrinos are suppressed by the Fermi coupling constant, we will work only to linear order in H I and will consider only 0-neutrino and 1-neutrino states. Therefore it will be convenient to decompose the Hamiltonian into a neutrino-number conserving piece H 0 and the neutrino creating term H I given in Eq. (3.1) where H 0 (x) is the free real scalar Hamiltonian density (4. 2) The neutrino-number conserving Hamiltonian is then where we have defined the eigenvalues The interaction H I interpolates between these two sectors The evolution of a 0-neutrino state is slightly more complicated than in the classical source case because H 0 does not annihilate the initial configuration, which now contains both a source particle and also an environment particle. Again, restricting attention to terms with precisely one H I we find Again we remind the reader that the terms with no H 0 have been dropped as they will not contribute to the matrix elements calculated below and also, as they contain no neutrinos, they will not contribute to our understanding of the neutrino wave packet. The projected state (4.7) can be written in terms of an integral over the time t 0 at which the neutrino was created (4.8) A measurement of a neutrino at a specific (x, t) would allow a determination of t 0 to within some uncertainty. However no measurement is implied here and so all values of t 0 ∈ [0, t] contribute to the amplitudes.
The Hamiltonian is again time-independent and so evolution conserves energy. E 0 and E 1 are not precisely the energies of the initial and final states, but rather the energies that they would have were they on-shell. The phase in (4.8) oscillates rapidly in t 0 unless E 0 = E 1 . Therefore the t 0 integral will be dominated by the stationary phase corresponding to the case in which the particles are on-shell. In this way we naturally recover the fact that particles are on-shell when t is large. This is also apparent in Eq. (4.7), where the (E 1 − E 0 ) in the denominator favors E 0 ∼ E 1 . Note that there is no pole as the numerator vanishes when E 0 = E 1 .
As the evolution operator e −iHt is linear, one can now easily evaluate the state at a time t The 3-momentum of the neutrino is p − q. The covariant wave packet conjecture states that |t only depends on p − q via Lorentz scalars. This is certainly not evident, but we will test this claim in the sequel, beginning with an initial condition which is itself a covariant wave packet.
Again we calculate the matrix elements corresponding to transitions to states with neutrinos in the flavor basis. The momentum space matrix elements arẽ where p is the final momentum of the source and q is the momentum of the neutrino. In neutrino measurements, often both the position and the momentum of the neutrino are determined with some known uncertainty. This motivates us to consider a transition amplitude in which both the momentum and the position of the neutrino are fixed, as in Eq. (3.7) (4.11) Eq. (3.8) then yields the approximate probability density for a transition to a state with a neutrino at position x with momentum q 0 We repeat that this is not the probability that the neutrino is measured, as no neutrinos are measured in our model.
In our model the environment only interacts with φ H , producing a shift in E 0,α (p+q) and so a relative phase between the two terms in the numerator on the right. This is ultimately responsible for decoherence. If, on the other hand, we introduce an additional coupling of the environment to both φ H and φ L with equal coefficients, it would produce an equal shift in both E 1,i (p, q) and E 0,α (p + q). The result would be an overall phase in the amplitude, which of course does not affect P (x, t) as this only depends on the absolute value of the amplitude. Therefore in this simple model we see that it is not the total interaction of the source with the environment which contributes to decoherence, as has been assumed in many calculations of decoherence such as Refs. [6,7], but rather the difference between the interaction with the source state before and after the neutrino production. In the case of a Coulomb interaction with a nucleus that produces a neutrino via β decay, this would correspond to the difference in the Coulomb interaction caused by a shift in the charge Z by one unit and the creation of a positron. We claim that this factorization argument is quite general, and not a specific feature of our model.

Numerical Results: Amplitudes
In this subsection we will fix the neutrino mass m 1 and the source masses M I to be

20
We consider only neutrinos whose momenta are equal to q 0 = 2 (4.14) to within an uncertainty σ. The initial width squared of the source particle will be β = 1. At time t = 50, we plot the amplitudes A iα (p, x, 50) in Fig. 6. Recoil momenta p of the source particles are set to p = −2 and p = −3. As we have assumed that the measured neutrino momentum is equal to 2, the amplitudes are in general supported at x > 0. However the source particle momentum p+q is, within σ, equal to 0 and −1 when p = −2 and p = −3 respectively. Therefore in the later case the φ H moved left and so the measured position of the neutrino tends to lower values of x. The phases oscillate quite rapidly, as can be seen, but it is the beating of the phases which leads to neutrino oscillations. Note that interference is only possible between final states with identical quantum numbers, including the recoil momenta. Therefore it is the beating at fixed p which yields neutrino oscillations. On the other hand, one sees that the large environmental energy shifts α considered here have an appreciable effect on the spectra already at t = 50. As the environmental state is not measured, the corresponding probabilities P α will be incoherently summed, degrading the oscillation signal.
Observe the fairly large fractional difference in the red curves corresponding to the two neutrino flavors in A i0 (−2, x, 50). This difference is due to the different phase space for the two masses. The difference is large because the mass difference is large. The difference in these amplitudes will damp the neutrino oscillations. Below, we will see this purely kinetic damping already in the partial probability distributions P α . Such damping is far too small to be observed at current ultrarelativistic neutrino experiments.
To reduce this purely kinetic source of oscillation damping, we will reduce our mass splitting by setting m 2 = 0.35 and we will worsen our momentum resolution to σ = 0.2 so that the experiment cannot hope to determine the neutrino mass eigenstate from a precise momentum measurement. To keep the similar contributions to the probabilities, we set At the late times at which oscillations occur. This has little effect on the phases, so we show the absolute values of the amplitudes for the smaller splitting in Fig. 7. Notice that the difference between the neutrino mass eigenstates is greatly reduced, as expected. In the case of the environment variable α = 2, one sees that the amplitude is quite small at intermediate x, and in fact vanishingly small at p = −2. This is easy to understand. Recall that the neutrino momentum is q = 2.0 ± 0.2. When p = 2, then p + q = 0.0 ± 0.2 and so E 0,2 (p + q) =∼ M H + 2 = 10.5, E 1,i (p, q) ∼ M 2 L + p 2 + q ∼ 9.8 ± 0.2 (4.19) and so the on-shell condition E 0 = E 1 is only satisfied when the momentum deviates from its measured value at more than the 3σ level. Similarly, when p = 3 one finds E 0,2 (p + q) =∼ M 2 H + 1 2 + 2 = 10.55, E 1,i (p, q) ∼ M 2 L + p 2 + q ∼ 10.1 ± 0.2 (4.20) and so the on-shell condition is excluded at about 2σ. This explains why the amplitude is small when p = 3, and very small when p = 2. The two peaks in the amplitude at low x and near the light cone are artifacts of the boundary conditions, as in the classical source case considered in Sec. 2.
At time t = 50 there are not yet any oscillations and certainly no decoherence. The amplitudes at t = 2000 are shown in Fig. 8. These are qualitatively similar to the t = 50 case. However the off-shell contribution at the boundary has become thinner. Note that while the integral of the off-shell shell region is greatly reduced at later time, as expected, nonetheless in the small region of x-space where it is visible due to boundary effects, the amplitudes at t = 50 and t = 2000 are similar.

Numerical Results: Probabilities
Let us return to the large splitting case m 2 = 0.4, σ = 0.1, c α = 2 3α/2 . The (partial) PDFs are shown in Fig. 9. Note that these PDFs are not localized in x as one would expect from wave packets. This is because all values of t 0 ∈ [0, t] are considered. If the source particles were measured, this would fix t 0 to within some precision and the resulting PDFs would be localized in x. Also a measurement of the neutrino would allow an approximate determination of t 0 .
The fractional amplitude of the oscillations does appreciably decrease with time, as expected. However this decrease is mostly present already in this partial probabilities. It therefore does not result from the environmental interaction, which is not present at all in P 0 (x, t). Rather this is the kinematic decoherence resulting from the fact that the higher mass neutrino has less phase space and so a lower amplitude, as was seen in Fig. 6.
To observe a clear signature of decoherence resulting from environmental interactions, we return to the small splitting case m 2 = 0.35, σ = 0.2, c α = 2 3α/4 . The corresponding (partial) PDFs are shown in Fig. 10. Now the difference in the amplitudes of the two neutrino mass eigenstates is smaller, as was seen in Fig. 8. Thus while the amplitude of the partial PDF oscillation does clearly shrink with time, this effect is less pronounced than it was in the In both cases one may observe that at lower values of x the oscillation phases differ for the various partial probabilities P α . By x ∼ 0 this difference is about 60 • . Therefore the total probability P , which is an incoherent sum of these partial probabilities, has a smaller oscillation amplitude at small x than the partial probabilities. This is the decoherence arising from destructive interference between the various environmental interaction eigenstates. One may observe in Fig. 10 that by x ∼ 0, at t = 3000, it nearly removes the oscillation minimum.
As one might expect, if the environmental interaction is weakened then so is the interference. In Fig. 11 we reduce the environmental interaction to One can see that the various partial probabilities P α oscillate with little phase difference and so constructively interfere. In this note we will not systematically study the necessary environmental interaction for decoherence to set in at a fixed time t. However in this example our results appear to be consistent with the thesis that for the first few oscillations should be of the same order as the neutrino momentum. It is also clear that decoherence has a large effect on the positions where the neutrinos have oscillated more times. In our figures this corresponds to the low values of x, but at JUNO it would correspond to the lower energy part of the spectrum.

Conclusions
In this note we have introduced a simple model of neutrino production, oscillation and decoherence due to environmental interactions of the source particle. This model was treated consistently in quantum field theory and is sufficiently simple that the various wave functions have been calculated explicitly, albeit numerically. Interactions between the source particle(s) and the environment yield a characteristic coherence time. The usual approach is to consider a Gaussian neutrino wave packet with width equal to this coherence time but then to neglect the entanglement with the environment, and often also the entanglement with the source. Following in the suggestion of [13], our approach is different. We have kept the full entangled state consisting of the neutrino, source particle and also the environment. Our first principles calculation of the neutrino wave function can be used to test various conjectures in literature, such as the covariant wave packet conjecture of Refs. [15,16]. We have not yet included a model of measurement, but to do so in the future will be straightforward. A consistent treatment of entanglement and measurement will allow us to test the revival mechanism of Refs. [5,14]. We have worked in a basis in which the environmental interactions H are diagonal. As the Hamiltonian is Hermitian, it may always be diagonalized in principle. While in the case of accelerator neutrinos, the interactions may be relatively simple [2] and so such a diagonalization is straightforward, in the case of reactor neutrinos there are a number of distinct interactions contributing to H and an explicit diagonalization would be difficult. However, our analysis suggests that the environmental interaction is appreciable only if the eigenvalues of α are not too far beneath the neutrino energy, or perhaps the neutrino energy divided by the number of oscillations. In the case of reactor neutrinos, interactions within the nucleus itself after a β decay may be expected to have characteristic energies of hundreds of keV, which would be sufficient. The inner electrons have binding energies of 10s of keV, and so interactions with these electrons may also cause noticeable coherence, at least in experiments such as JUNO that are sensitive to many oscillations. On the other hand interatomic interactions, which are commonly used to set the coherence scale [8,9], have energy scales of eV, and so are unlikely to have noticeable decoherence effects in any proposed reactor neutrino experiment. We have seen that only the difference between the interaction strength before and after the neutrino emission contributes to decoherence, further reducing the impact of interatomic interactions.