Remark on neutrino oscillations

The oscillations of ultra-relativistic neutrinos are realized by the propagation of assumed zero-mass on-shell neutrinos with the speed of light in vacuum combined with the phase modulation by the small mass term $\exp[-i(m^{2}_{\nu_{k}}/2|\vec{p}|)\tau]$ with a time parameter $\tau$. This picture is realized in the first quantization by the mass expansion and in field theory by the use of $\delta(x^{0}-y^{0}-\tau) \langle 0|T^{\star}\nu_{L k}(x)\overline{\nu_{L k}(y)}|0\rangle$ with the neutrino mass eigenstates $\nu_{L k}$ and a finite positive $\tau$ after the contour integral of the propagating neutrino energies. By noting that the conventional detectors are insensitive to neutrino masses, the measured energy-momenta of the initial and final states with assumed zero-mass neutrinos are conserved. The propagating neutrinos preserve the three-momentum in this sense but the energies of the massive neutrinos are conserved up to uncertainty relations and thus leading to oscillations. Conceptual complications in the case of Majorana neutrinos due to the charge conjugation in $d=4$ are also discussed.


Neutrino oscillations and mass expansion
The phenomenon of neutrino oscillations [1,2,3] is fundamental to measure the small neutrino masses, and it would be disastrous if the different formulations should lead to different neutrino masses.If one writes the neutrino mixing with the PMNS matrices U αk where |ν k are the mass eigenstates which diagonalize the neutrino mass matrix, and the flavor eigenstates |ν α , (α = e, µ, τ ), are related to each other by the above mixing formula.We define the charged lepton flavor eigenstates by the mass eigenstates.One may start with the production of the flavor eigenstate neutrino ν µ in the energetic pion decay for example, by measuring π + and µ + , and the neutrinos thus produced propagate toward the detector in the oscillation experiment; the direction of each mixed neutrinos may not necessarily be in the exact specific direction considering the accuracy of the measurements of π + and µ + .The oscillation is observed in each direction of the mixed neutrinos.We analyze the Dirac neutrinos in the main part of the present paper for simplicity, and the case of the Majorana neutrinos, which are constrained by the complications of the Majorana fermions in d = 4, shall be discussed in Appendix.It is known that the relation (1), if interpreted as a superposition of on-shell mass eigenstates with identical three-momentum, leads to the standard oscillation formula [4] | ν β (0)|ν α (t) where t = L, the neutrino propagation distance, is assumed together with | p| 2 ≫ m 2 ν k .It is known also that the identical energy assumption of neutrinos, instead of the identical three-momentum assumption in the above derivation, gives essentially the same formula [5].In this paper we want to understand how the oscillation formula is robust against various ways to derive it.
We first mention the idea of the wave packet of neutrinos in the first quantization formalism [6].If one of the mass eigenstates in the neutrino ν µ , for example ν 1 , should be identified immediately after the pion decay (2) such a mass eigenstate due to the reduction of quantum states would propagate without oscillations, although the flavor change β → α would be induced by the (inverse) mixing in (1).The actual values of neutrino masses are, however, very small and thus the specification of the mass of one of neutrinos is practically impossible; in fact, as is explained later, the neutrino masses need to be not measured by the conventional detectors to measure the neutrino oscillations.To treat the un-identified neutrino mass eigenstates consistently, Kayser [6] suggested the idea of the wave packet of particles involved, such as ν µ in (2).The wave packet is more generally understood as a means to incorporate the semi-classical aspects of neutrino oscillations into quantum mechanics in a consistent manner, and it has been successfully incorporated in the field theoretical formulations [5,8,9,10,11,12] and a related quantum mechanical formulation [13].It has been shown also that the kinematics of the neutrino production generally implies the mass dependence of the neutrino momentum such as p = p(m 2 ν k ) and that the mass-dependence of the momentum depends how they are produced [14,15]; for example, the two-body decay π + → µ + + ν or other neutrino production processes.
One may consider the propagating phase of a flavor eigenstate in vacuum determined by the phase of the mass eigenstates initially located at (t, x) = (0, 0) by assuming that the neutrino masses are very small and thus the measured neutrinos are ultra-relativistic, in accord with experimental facts which imply the mass differences ∆m 2 on the order of (10 −2 eV ) 2 .One then obtains the Lorentz invariant where with We assume that the momentum p(0) is common to all the mass eigenstates of neutrinos.The velocity v g is the group velocity of the propagating (now regarded as massless) neutrinos.The basic assumption of the wave packet (although we do not write an explicit form of the wave packet following the analysis by Giunti and Kim [15]) is that the neutrinos are concentrated at the center of the wave packet (see also [5]) The neutrino wave packets are essentially the spreading of the initial and final state weak vertices, since the neutrinos rarely interact with surrounding materials.We assume that the geometrical spreads of the weak vertices, which are large in the microscopic sense so that the spread of the neutrino momentum is negligibly small by the uncertainty principle, but still the geometrical spreads are very small compared with the macroscopic distance between the two weak vertices.In other words, all the neutrino mass eigenstates are assumed to be measured at x − v g t ≃ 0 or at a finite distance away from 0, then the term (5) for large t; this is also ensured by the fact that |m 2 | p(0)| , which is confirmed to be the case.The small quantities with O(m 4 ν k ) or higher powers in the neutrino mass are neglected.One thus measures the oscillations caused by next to the last term of (5) with the common momentum factor When the exponential factor is written in the form is universal, i.e., depends on the intrinsic properties of the neutrinos independently of how the neutrinos were produced.The common factor p(0) x − E(0)t for each mass eigenstate in (5) does not contribute to the oscillation.We thus recover the standard oscillation formula (3) using t ≃ L which arises from x − v g t ≃ 0 in (8) with | v g | = 1.Physically, the semi-classical relation t ≃ L with a small error, which is determined by the conventional detectors, is not influenced by the neutrino masses.In the above analysis, we chose the vanishing masses of propagating neutrinos as the fiducial values and the observed neutrinos are assumed to be essentially massless.
Alternatively, if one assumes the time-to-distance conversion t = L to be a valid ansatz for massless on-shell neutrinos1 , one would obtain the formula (9) directly since the term m 2 ν k p ′ (0) • ( x − v g t) then vanishes.Also this picture is consistent with the exact three-momentum conservation induced by the three dimensional integration at both of the initial and final weak vertices, which are implicitly assumed.We derived the formula (5) by the mass expansion, which is an expansion in terms of a Lorentz scalar quantity of the Lorentz invariant phase p(m 2 ν k ) x − E(m 2 ν k )t; this may imply the Lorentz invariance of the oscillation formula ( 9) [5].The present picture may agree with our intuitive understanding of neutrino oscillations; the ultra-relativistic neutrinos propagate with the speed of light x − v g t = 0 in vacuum for a measured momentum p(0), and the effects of the small mass differences | p(0)| provides an important quantity in this analysis, namely, it needs to be very small and not measured by the conventional detectors; this constraint, namely, not measurable by conventional detectors, generally arises because of the energy non-conservation in neutrino oscillations (or by an analysis of energy-time uncertainty relations).In the next section on the Feynman amplitude approach to neutrino oscillations we discuss how the same criterion arises.

Feynman amplitude approach
To understand the oscillation phenomena in a field theoretical formulation, one may start with an extension of the Standard Model.The leptonic sector is given by with a 3 × 3 unitary mixing matrix U in (1).We ignore the neutral current and electromagnetic interactions.All the particles belong to respective mass eigenstates and the lowest order Feynman amplitudes (using the Fermi approximation) are well-defined without infrared singularities.One can confirm that the conventional tree-level Feynman amplitude, which consists of the production and detection weak vertices connected by the Feynman propagator of massive neutrinos, does not give rise to the neutrino oscillations when one integrates over all the values of the weak interaction points x µ and y µ with x µ > y µ as well as y µ > x µ in the Fermi approximation of weak interactions we work.This integration over interaction points, which incorporates backward moving off-shell anti-neutrinos as well as forward moving off-shell neutrinos, preserves the energy-momentum precisely at each interaction point.We thus have no time scale to measure oscillations which are related to the time translation non-invariance.We shall demonstrate below that the Feynman amplitude with only a part of the forward on-shell neutrino propagator reproduces the neutrino oscillation amplitude (3) by preserving the measured overall energymomentum up to uncertainty relations.Feynman rules are used to specify the quantum mechanically allowed couplings.The Feynman amplitude approach to neutrino oscillations has been discussed in the field theoretical formulation by Kobzarev et al, [7], Grimus and Stockinger [9], Giunti, Kim and Lee [10] who emphasized the wave packets, and using plane waves by Egorov and Volobuev [16,17,18], among others.In the latter approach [16,17,18], the (effective) limit x 0 − y 0 → ∞ was considered in the propagating neutrinos of the form δ(x 0 − y 0 − τ ) T ⋆ ν l L (x)ν k L (y) using a generalization of the Grimus and Stockinger theorem [9] and thus achieving the on-shell condition of all the propagating neutrinos.They emphasized that the momentum space Feynmanlike amplitude thus defined produces the oscillation amplitude [16,17,18] and the probability interpretation of the oscillation amplitude is justified based on the probability interpretation of the conventional Feynman amplitude.They note the simplicity of their formulation compared to those of the past formulations such as [9].We follow the basic ideas of [9,10,16,17,18], and we shall simplify the derivation of oscillation amplitudes and add several remarks on the robustness of the amplitude thus derived.
We write the effective Lagrangian of neutrino processes as by incorporating the neutrino production and detection processes in the sources J α R (y) and J β R (x), respectively, which are chosen generally not to be Hermitian conjugate to each other.We have a generalization of Schwinger's source functions3 For example, } describes the electron production of ν e +n → p+e.The source functions (in a generalized context as above) are expressed in terms of plane waves as in the conventional Feynman amplitudes.In the neutrino oscillation experiments, the macroscopic distance between the production vertex and the detection vertex is one of the main observables.Following [16,17,18] (see also [15]), we work in the framework of plane waves and consider the configurations in a 4-dimensional sense where the neutrino production and absorption points, which are denoted by y µ and x µ , respectively, are correlated by a fixed time difference with a very large macroscopic τ > 0.
A suitable choice of J β R (P f ) and J α R (P i ) specifies the initial and final systems of the neutrino oscillation experiments as described above.We then have the oscillation amplitude for the initial (production) vertex such as and the final (absorption) vertex such as with the neutrino oscillations communicating ν µ (α = µ) to ν e (β = e).Our proposal is to analyze the conventional Feynman amplitude for the weak process, where the neutrino is exchanged, which is written by the prescriptions of source functions described above with an extra δ-function as where P i = q π − p µ is the entering four-momentum in the case of the pion decay, and P f = p p + p e − p n is the four-momentum of the outgoing final system; m 2 ν k stands for the neutrino mass eigenvalue squared.Up to this point, the formula is faithful to what defined by the first line in (18).
We now make an approximation.Since the energy resolution of conventional detectors cannot detect neutrino masses, we neglect the possible neutrino mass dependence in P i and P f ; those four-momenta are written as if all the propagating neutrinos are massless.The summation k over the neutrino masses then operates only on the neutrino propagators, and the formula ( 18) is written after the contour integral over the neutrino energy as The effect of the energy non-conservation induced by δ(x 0 − y 0 − τ ) is still seen by the presence of the time parameter τ in the formula.
We would like to add several comments on the above derivation of the formula (19).It implies that the Feynman amplitude with the modified Feynman propagator [16,17,18,19] describing only a part of the forward propagating neutrinos gives rise to the neutrino oscillation probability for large fixed τ , as is seen in ( 24) later; the formula (19) gives the conventional result (3) if one assumes the timeto-distance conversion τ = L.In the present paper, we adopted the propagator of neutrinos (20) which is the same as in [16,17,18] with a fixed large τ , but we obtain the on-shell condition of neutrinos by performing the contour integral with respect to the neutrino energy [19], instead of taking the (effective) limit x 0 − y 0 → ∞ in [16,17,18] with the help of the Grimus and Stockinger theorem.By this way we performed the above calculations of the amplitude (19) with an approximation stated above, in the lowest order of perturbation.This simplified evaluation was possible since we assumed that the energy-momentum of the initial and final systems, represented by P i and P f , respectively, are independent of the neutrino masses because of the limited accuracy of conventional detectors of weak interactions.The neutrino momentum p = P i also becomes mass independent.By this assumption we were able to take the summation over the massive neutrinos k outside the δfunction as in the final formula (19) 4 .We are assuming that the neutrinos and other do not appear in the above formulation (19), although we used τ = L at several places already.This absence of the analysis of the time-to-distance conversion is analogous to the case of ( 3) and thus we have to remedy the shortcomings, although the direction from y to x in the present case is specified by the given common momentum p = P i .The simplest idea may be to assume that the time-to-distance conversion is a valid ansatz in the analysis of neutrino oscillations for effectively massless on-shell neutrinos, since the conventional detectors cannot recognize the neutrino masses.One would then obtain the desired result from (24) (or ( 5)) directly.This abstract picture is consistent also with the integration over y 0 , namely, the arbitrariness of the origin of initial time in obtaining (19), and also with the exact conservation of the three-momentum at the two weak vertices following from the integration over x and y.Further discussions on this matter shall be given later.
For the specific two-flavor case and the non-diagonal process µ → e, for example, the formula (24) gives the well-known oscillation probability The interval of L(= τ ) to measure the oscillations is then specified by the standard depending on the mass difference of neutrinos, and the momentum p carried by the (massless) neutrinos which is determined by the measured p = P f ; this value is assumed to be independent of neutrino masses.The precise energy-conservation (i.e., time-independence) in ( 25) is given by namely, the vanishing oscillation 1 − cos[ In general, the absolute normalization of the oscillation probability is not wellspecified [5], in particular, in the present case multiplying the Feynman amplitude by δ(x 0 − y 0 − τ ) and thus using only a part of the Feynman amplitude, but the specific oscillation probabilities are well-normalized as above.

Discussion and conclusion
The observed oscillation of ultra-relativistic neutrinos is based on the two basic conditions.The first is that the massless on-shell neutrinos (with various originally small mass eigenstates) propagate for a given momentum p(= P i ) from the position y to another x with the speed of light in vacuum when measured with the conventional detectors which do not recognize the neutrino masses; the present consideration is thus limited to length scales of the neutrino propagation approximately within those of the atmospheric neutrino oscillations.This may be called the semi-classical (particle) aspect.The second is that the oscillations are caused by the small phase modulation with the common momentum factor k such as in ( 24) and ( 5).This phase modulation is the quantum mechanical (wave) effect and supplemented by the subsidiary condition τ = L = | x − y| following from the semi-classical consideration (29).Any formulation satisfying these two conditions (29) and (30) gives the oscillation formula.We assumed in Section 2 that the oscillation formula is, in principle, applicable to any experiments without referring to the specific positions y and x of weak vertices, as the abstract formulations of ( 5) and (24) suggest.The assumption of the time-todistance conversion was accepted as a valid ansatz for semi-classical on-shell massless neutrinos with the (common) given momentum p.In this understanding, the oscillation formulas are valid for any (4-dimensional) configurations parallel transported from each other with fixed p, independently of the specification of the origin of time y 0 .The neutrino oscillations are universal phenomena and applicable to any chosen x and y in the direction of p with τ = L = | x − y|.The simple Feynman diagram approach in [16,17,18] and the present derivation of the oscillation formula (19) may be counted among the schemes based on these assumptions.
Alternatively, one may follow the elaborate wave packet analyses in the past [5,8,9,10,11,12,13].These analyses may be regarded as clarifying the mechanisms how to satisfy these conditions (29) and (30) including the quantum coherence of neutrinos and the specification of positions y and x.The weak vertices at x and y are treated naturally in these wave packet pictures, while the treatment of weak vertices are less transparent in the present field theoretical treatment with Lorentz invariant plane waves.In the wave packet picture, one may assume that the semi-classical relation ( 29) is valid for the points y and x that are by themselves spreading over the three-dimensional domains so that the three-momentum conservation is ensured within the constraints of uncertainty relations.The spreads of the points y and x of two weak vertices are still assumed to be much smaller than the semi-classical distance L = | x− y| between them; as for an explicit wave packet realization of these conditions see, for example, [10].In fact, the wave packet formalism is regarded as an attempt to incorporate semi-classical constraints such as (29) in the framework of quantum mechanics consistently.In this sense, the justification of the present Feynman amplitude approach is also given by the idea of wave packets.
The on-shell neutrinos with the same three-momentum and different masses mean that the energy conservation is not satisfied in the intermediate states of the oscillation in (19) (also in ( 5)), as is well known [5].The time-dependent neutrino oscillation in the case of two flavors in (25), for example, may be regarded as measuring the effective energy-nonconservation in the intermediate states (with the notation with explicit ) where standing for the propagating neutrino energy splitting.The energy non-conservation is manifested as time-dependent oscillations, i.e., the breaking of time translation invariance.But one may regard that the energy non-conservation is superficial since the energy-time uncertainty relation ∆Eτ ≥ /2 is satisfied naturally although rather marginally by (31) with macroscopic τ .This shows that ∆E is actually very small and thus measurable only with oscillations, of which effects vanish on average, but not measurable with the conventional detectors.If ∆E should be measured by conventional detectors, the energy-nonconservation would be confirmed by oscillations; this would be a contradiction 5 .Without the delta-functional constraint δ(x 0 − y 0 − τ ), one would have the conventional Feynman amplitude with the propagating neutrinos where p = P i (now exact), and because of Feynman's iǫ both the forward propagating neutrinos and the backward propagating anti-neutrinos contribute.In contrast, in our formula (19) only a part of the neutrinos propagating forward in time appear and those anti-neutrinos propagating backward in time with negative energy are neglected 6 .With the presence of δ(x 0 − y 0 − τ ) in the measurement of oscillations, the neutrino propagation may be regarded as a large scale quantum effect.One may regard that the microscopic CP symmetry of the oscillation amplitude is an indication of the CP symmetry of the original Feynman amplitude; the microscopic CP violation is described by the phase of the PMNS matrix in the case of three or more leptonic flavors (except for an extra CP violating U(1) phase even for two flavors in the case of Majorana neutrinos, which is however not measured by oscillation experiments [20]).
In conclusion, we first re-formulated the ultra-relativistic neutrino oscillations using the neutrino mass expansion, which simplifies the formulation substantially in the framework of the first quantization.The characteristic property of neutrino oscillations is that the conventional detectors cannot recognize the finite neutrino masses, and thus the absence of energy-momentum conservation in a precise sense.These ideas were then applied to the Feynman amplitude approach to neutrino oscillations with plane waves by constraining the amplitude to a part of forward propagating massive neutrinos.The time-to-distance conversion, which may be justified for effectively massless on-shell neutrinos, was assumed to simplify the analyses of ultra-relativistic neutrino oscillations.A unified picture of ultra-relativistic neutrino oscillations was thus presented.[22].The Weinberg's model is defined by an effective Hermitian Lagrangian [21] L where M ν stands for the 3 × 3 diagonalized neutrino mass matrix and we defined 2 )ν L (x) vanishes under the pseudo-C symmetry.Apparently an idea of the pseudo-C symmetry needs to be understood better [23].
A way to deal with Majorana neutrinos in a consistent manner may be to use a general class of seesaw models [22], which contain the equal number of left-handed and right-handed fermions.One may use a Bogoliubov-type transformation to change the definition of the vacuum of the Weyl fermion to the vacuum of the Majorana fermion [24].The conventional approach gives rise to in the seesaw model, which are Majorana fermions with masses M + = M − if one uses the pseudo-C symmetry.But the conventional parity is not well-defined.
Those fields ψ ± (x) in (39) are the Majorana fermions with well-defined C and P in the conventional sense.

with C = iγ 2 γ 0 .
The field ψ(x) satisfies the classical Majorana condition ψ(x) = Cψ(x) T (36) identically regardless of the choice of ν L .One may define the Majorana fermion by (35) together with the Dirac equation [iγ µ ∂ µ −M ν ]ψ(x) = 0.This is the conventional procedure.In general one cannot define simultaneously a Majorana fermion with well-defined C and P and a Weyl fermion for which C nor P are defined in d = 4.The conceptual complications are how to define an isolated free Majorana neutrino with well-defined C and P, while weak interactions are described by a chiral fermion.Starting with a Dirac fermion, one may obtain under the charge conjugationν L (x) → Cν R (x)T (37) instead of the pseudo-charge conjugation symmetry ν L (x) → Cν L (x) T [23] implicit in the fermion (35).If one adopts the pseudo-C symmetry ν L (x) → Cν L (x) T together with the representation of a Majorana fermion (35), one would encounter various puzzling aspects.For example, the first expression of the effective Lagrangian (34) is not invariant under the P operation; also the mass term ν T L (x)CM ν ν L (x) = ν T L (x)CM ν ( 1−γ 5 If