Chiral oscillations in the non-relativistic regime

Massive Dirac particles are a superposition of left and right chiral components. Since chirality is not a conserved quantity, the free Dirac Hamiltonian evolution induces chiral quantum oscillations, a phenomenon related to the \textit{Zitterbewegung}, the trembling motion of free propagating particles. While not observable for particles in relativistic dynamical regimes, chiral oscillations become relevant when the particle's rest energy is comparable to its momentum. In this paper, we quantify the effect of chiral oscillations on the non-relativistic evolution of a particle state described as a Dirac bispinor and specialize our results to describe the interplay between chiral and flavor oscillations of non-relativistic neutrinos: we compute the time-averaged survival probability and observe an energy-dependent depletion of the quantity when compared to the standard oscillation formula. In the non-relativistic regime, this depletion due to chiral oscillations can be as large as 40$\%$. Finally, we discuss the relevance of chiral oscillations in upcoming experiments which will probe the cosmic neutrino background.


I. INTRODUCTION
Dirac equation has unique dynamical predictions for fermionic particles, from the Klein paradox [1,2], related to pair production in scattering problems, to the Zitterbewegung, the trembling motion of free relativistic states [3]. Formally, the solutions of the Dirac equation are given in terms of 4-component spinors, or Dirac bispinors. From a group-theoretical perspective, those objects belong to the irreducible representations of the complete Lorentz group and are constructed by combining Weyl spinors of different chiralities [4]. Any Dirac bispinor has left and rightchirality components which are dynamically coupled by the mass term of the Dirac equation. Since a bispinor with definite chirality cannot be an eigenvector of the massive Dirac Hamiltonian, the free evolution of such state will induce left-right chiral oscillations [5]. The degree of left and right chirality superposition in a Dirac bispinor depends on the energy to mass ratio, thus chiral oscillations are usually not relevant for the description of relativistic particles, but are prominent for dynamical regimes in which the momentum of the particle is comparable to (or smaller than) its mass.
The non-relativistic regime is also interesting for exploring the connection between the Zitterbewegung effect and chiral oscillations. As discussed in [6], the trembling motion of free Dirac particles has an intimate relation to chiral oscillations. Both are related to the fact that a state with initial definite chirality must be described as a superposition of positive and negative energy solutions of the Dirac equation. Although the Zitterbewegung effect is associated with a fast frequency, which is averaged out in typical high energy systems, it has already been considered to explain the Darwin correction term in hydrogenionic atoms [2] and it has been probed in Dirac-like systems, such as graphene [7] and trapped ions [8,9]. For example, the graphene analogous to the Zitterbewegung [10] can be measured via laser excitations [11]. While mono-layer graphene is described by the massless Dirac equation, bilayer graphene displays effects associated with the mass term of the Dirac equation [7] plus an effective non-minimal coupling [12], and is a system in which chiral oscillations and their relation to the Zitterbewegung could be probed.
In this paper, we are concerned with the non-relativistic limit of chiral oscillations within the framework of the Dirac equation. Since the chiral oscillation amplitude depends on the mass to energy ratio [13,14], the minimum survival probability of an initial state with definite chirality could be averaged out in the non-relativistic limit. Therefore, in the non-relativistic limit we specialize our discussion to (Dirac) neutrinos which indeed are ideal candidates for the study of chiral oscillations due to the chiral nature of weak interaction processes in which they are produced and/or detected [15]. The Cosmic Neutrino Background (CνB) is an important example of neutrinos in non-relativistic regime.
The CνB is the neutrino counterpart of the Cosmic Microwave Background (CMB) for photons: it is composed by decoupled relic neutrinos evolving in an expanding background [15]. The PTOLEMY project [16] aims at detecting the CνB via neutrino capture on tritium, a goal which should be reached within the next years. In previous studies [17,18], it has been noticed that the expected event rate of CνB capture on tritium will exhibit a depletion. We show that such a depletion for Dirac neutrinos can be understood as a manifestation of chiral oscillations and give a precise quantification of it. Furthermore, the capture rate depends on the nature of the neutrinos: Majorana or Dirac. While chiral oscillations are present in both cases, Majorana neutrinos are subjected to more measurable oscillation channels than the Dirac case, e.g. to right-chiral and positive helicity states.
Although the interplay between chiral and flavor oscillations is very small in the ultra-relativistic regime [13,14,19,20], it is relevant for describing dynamical features of non-relativistic neutrinos. We describe chiral oscillations in the context of two-flavor neutrino propagation, including flavor oscillations. In particular, we compute the effects of chiral conversion on the averaged flavor oscillation, that is, the survival probability of a neutrino of a given flavor averaged over one period of flavor oscillation. The chiral oscillations corrections are relevant when the particle's momentum is comparable with the lightest neutrino mass and, in the non-relativistic regime, the maximum difference is 40%, a prediction consistent with the preliminary discussion pointed in [18] in connection with CνB detection. A more in depth discussion of chiral oscillations in the context of the CνB can be found in [21]. To highlight the main dynamical features associated to this phenomenon, we adopt a simple plane wave description and leave a more realistic treatment involving wave packets [13,14,20] for future work. Although we focus on Dirac neutrinos, we discuss chiral oscillations effects in Majorana neutrinos and provide a detailed calculation of them for the general Majorana-Dirac mass term in the appendix.

II. CHIRALITY AND CHIRAL OSCILLATIONS IN BISPINOR DYNAMICS
Through this paper we will describe dynamical features of free massive fermionic particles in the context of relativistic quantum mechanics. We thus consider the temporal evolution as given by the Dirac equation (hereafter we adopt natural units = c = 1)Ĥ wherep is the momentum operator and m is the particle's mass. The operatorsα i (i = x, y, z) andβ are the 4 × 4 Dirac matrices satisfying the anti-commutation relationsα iαj +α iαj = 2δ ijÎ4 ,α iβ +α iβ = 0, andβ 2 =Î 4 . In the language of group theory, Dirac equation is the dynamical equation for the irreducible representations of the complete Lorentz group 1 , which are the well-known Dirac bispinors [4]. These are constructed by means of the irreducible representations of the (proper) Lorentz group, the Weyl spinors. The later belong to the irreps of a group isomorphic to the SU (2) and as such they carry an intrinsic degree of freedom -the spin. The elements of the two disconnected irreps of the proper Lorentz group are labeled as left and right-handed spinors. Since parity connects the left and right representations, the irreps of the complete Lorentz group are obtained by combining these representations, and thus the Dirac bispinors carry not only the spin, but also another intrinsic discrete degree of freedom, the chirality. Any massive Dirac bispinor is thus a combination of left and right-handed spinors.
Turning our attention to the Dirac equation framework, chirality is the average value of the chiral operator γ 5 = −iα xαyαz . While such operator commutes with the momentum term of the Dirac Hamiltonian, it does not commute with the mass term and thus it is not a dynamically conserved quantity for bispinor states describing massive particles. This can be better appreciated by considering the chiral representation of the Dirac matrices, in whicĥ γ 5 = diag{Î 2 , −Î 2 } 2 . Any bispinor |ξ can be written in this representation as where |ξ R,L are, respective, the positive (right-handed) and negative (left-handed) chirality two-component spinors. The Dirac equationĤ D |ξ = i|ξ can then be written aŝ 1 The complete Lorentz group is the proper Lorentz group plus Parity. 2 In the chiral representationα = σ 0 from which it is clear that the mass term mβ connects the left and right-handed components of the bispinor. For a given initial state, evolution under the free Dirac HamiltonianĤ D induces left-right chiral oscillations. In order to describe this dynamical effect we first introduce the positive and negative plane wave solutions of the Dirac equation, |ψ + (x, t) = e ip·x−iEp,mt |u s (p, m) and |ψ − (x, t) = e −ip·x+iEp,mt |v s (p, m) which are given in terms of the bispinors where E p,m = p 2 + m 2 . In the above equations, |η s (p) is a two component spinor that depends on the spin polarization of the particle. We notice that the orthogonality relations read: For now on, we describe the solutions of the Dirac equation with helicity bispinors, that is, we assume that |η s (p) are eigenstates of the Helicity operator p·σ p . This choice is convenient since helicity is a conserved quantity and all the relevant dynamical features will be entirely related to chiral oscillations. Moreover, we simplify our analysis by considering one-dimensional propagation along the e z direction, such that where |± are the eigenstates of the Pauli matrixσ z . For ultra-relativistic particles, m/p → 0 and positive helicity bispinors have only right-handed components while negative helicity bispinors have only left-handed components. On the other hand, for p/m → 0 the left and right-handed components are equal irrespective of the helicity. We now describe chiral oscillations by considering the temporal evolution of the initial state |ψ(0) = [0, 0, 0, 1] T which has negative helicity and negative chirality:γ 5 |ψ(0) = − |ψ(0) . As we are dealing with plane wave states propagating one-dimensionally, we consider the dynamical evolution in momentum space [14]. The time evolved state |ψ m (t) is given by Chiral oscillations are generated by the massive character of the particle and by the fact that the initial state is a superposition of positive and negative energy eigenstates of the Dirac Hamiltonian. The later is also responsible for the relation between chiral oscillations and the Zitterbewegung effect [6]. The survival probability of the initial state P(t) is given by while the (transition) probability of being in a positive chirality state is given by P T (t) = 1 − P(t), and the average value of the chiral operator γ 5 (t) reads According to (7) the minimum survival probability is , and the frequency of the chiral oscillations is (in natural units) E p,m . Therefore, the corresponding period of one chiral oscillation is τ Ch = 2π/E p,m and as the particle propagates freely, the length corresponding to one chiral oscillation can be evaluated as l Ch = pτ /E p,m = 2π p E 2 p,m . In Fig. 1 we show P min as a function of p/m (a) and l Ch as a function of p for several mass values (b). For non-relativistic states, i.e. states for which p ∼ m, chiral oscillations play an important role and affect significantly the probability of the state to be in its initial configuration. In fact, for p m the minimum survival probability is P min ∼ p 2 m 2 which vanishes as p → 0. In this case, since the eigenspinors are a maximal superposition of left and right chiralities, the free evolution induces a complete oscillation from left to right chiral components. This can also be seen in (8) value −1 and 1. For p m, chiral oscillations are less relevant for the state dynamics. In the limit p → ∞ chirality and helicity coincide, thus the eigenspinors (5) have definite chirality. In fact (8) is constant in the later limit. The chiral oscillation length l Ch is shown in Fig. 1 (b) for several values of the mass. For masses in the range of eV/c 2 , the expected chiral oscillation length is of the order of 1 c/eV ∼ 10 −6 m.
Our analysis has focused on Dirac particles with states given in terms of Dirac bispinors. For Majorana particles, one should have in mind that the mass term is proportional to the charge-conjugated spinor. In fact, the Dirac equation with the Majorana mass term, which we call Majorana equation, reads [22][23][24][25] i∂ |ψ =p ·α |ψ + mβ |ψ c , where |ψ c = iβα y |ψ * [31] is the charge conjugated bispinor. If the bispinor |ψ is its own charge conjugated, that is, if |ψ c = |ψ , then its temporal evolution is given by the "usual" Dirac equation (1). This last condition is known as the Majorana condition. It can be shown (see the Appendix) that for any bispinor satisfying the Majorana condition γ 5 = 0. In other words, any bispinor that is its own self conjugate has zero average chirality. Since the Majorana condition is preserved under the time evolution, there is no chiral oscillations in bispinors satisfying the Majorana condition, irrespective if the bispinor is an eigenstate of the Hamiltonian. Notice that this is in stark contrast with the Dirac bispinor case: a Dirac bispinor is an unconstrained object whose time evolution is given by the usual Dirac equation. The average chirality of a Dirac bispinor is only constant if it is an eigenstate of the Dirac Hamiltonian.
Nevertheless, given the intrinsic handedness of the weak interaction, we can consider a bispinor that is initially in a left handed and negative helicity state but whose time evolution is given by the Majorana equation. In this case, we follow the formalism of [26,27] and obtain that the survival probability of such state is the same as the one given in (7) while there is a transition probability to a right-handed component, associated with the charge conjugation of the left-handed initial state. This fact was also briefly quoted in [21]. A more general situation includes both Dirac and Majorana mass terms with two non-degenerate Majorana masses, which we describe in the appendix.

III. FLAVOR MIXING AND CHIRAL OSCILLATIONS IN NON-RELATIVISTIC REGIME
We now study the effects of chiral oscillations in non-relativistic neutrino mixing. The state of a neutrino of flavor α at a given t is given by the superposition of mass states: where U is the mixing matrix, and |ψ mi (t) are the bispinors describing the temporal evolution of the mass eigenstate |ν i with mass m i [13,14]. The state at t = 0 reads with |ν α ≡ i U α,i |ν i , and thus |ψ mi (t = 0) = |ψ(0) . Since weak interaction processes only create left handed neutrinos, for now own, we take |ψ(0) as the left handed bispinor with negative helicity of the last section. The temporal evolved flavor state is therefore where the |ψ mi (t) are given by eq. (6) with the substitution m → m i . The survival probability, i.e. the probability of the state |ν α (t) to be a left-handed state of α flavor, reads For two flavors mixing, the time evolution of an initial electron neutrino state is [14] |ν and the survival probability can be decomposed as In this formula P S e→e (t) is the standard flavor oscillation formula and are correction terms due to the bispinor structure. Those corrections include an interplay between chiral and flavor oscillations effects and are in agreement with results presented in the literature [19]. While the standard flavor oscillations have a time scale set by the energy difference E p,m2 −E p,m1 , chiral oscillations depend roughly on the mass-momentum ratio (see the previous section). The terms A e and B e depend non-trivially both on the energy difference and on the individual energies E p,m2 and E p,m1 . In the limit p and if the mass difference is |m 2 − m 1 | m 2,1 or m 2,1 , there will be no interference effects between flavor and the chiral oscillations. Those phenomena can be individually identified in the survival probability. If |m 2 − m 1 | is comparable to one of the masses, there will be an interference between flavor and chiral oscillations. We furthermore notice that the term B e (t) is equivalent to corrections obtained via a quantum field description of flavor oscillations [28,29]. In fact, after some algebra we obtain that P S e→e (t) + B e (t) reproduces eq. (30) of [29] (see also [30]). In the quantum field framework, the corrections have a similar origin to the usual Dirac Zitterbewegung, being originated from the fact that the flavor ladder operator has contributions from both particle and antiparticle operators of the massive fields. Since those corrections do not take into account chiral oscillations, we interpret A e (t) as the corrections due to chiral oscillations.
We show the survival probability (15) in Figure 2 for several values of p/m. As already anticipated in the previous section and in [13,14,20], the chiral oscillations induce just small corrections to the state's dynamics in the ultrarelativistic regime p m. For the non-relativistic regime p m 1,2 , the full oscillation formula differs significantly from the standard result due to high amplitude chiral oscillations. For the parameters of Figure 2, the mass difference is very small, and the flavor oscillations have a time scale much longer than the chiral oscillations.
where τ 12 = 4π Ep,m 2 −Ep,m 1 is the period of one flavor oscillation according with the standard formula (16). The time integration of (15) leads tō Considering that the average of the standard oscillation formula (16) is given byP one concludes that the quantityP S e→e −P e→e contains all the chiral oscillation effects on the averaged survival probability, and therefore properly quantifies the effects of chiral oscillations on flavor oscillations.
In Figure 3 we show the differenceP S e→e −P e→e as a function of the neutrino momentum and the squared mass difference and as a function of p/m 1 for fixed values of the squared mass difference. The difference between the averaged standard probability and the full one becomes 0.1 for p m 1 , indicating the influence of chiral oscillations on flavor oscillations. In the non-relativistic regime, the difference between the probabilities can be as big as ∼ 40% for typical values of the mixing angle. Moreover, we notice that for bigger values of the mass difference, oscillations in the probability (with the momentum) are observed. This is due to the terms ∝ f 21 sin (4π/f 21 ) in (19), which depend on the sum of the energies, and oscillate with the momentum when m 1 and m 2 are well separated. In the ultra-relativistic regime p m 2,1 , we can expand the full averaged formula (19) with respect to m 2,1 /p as to havē The corrections to the standard result are thus proportional to both the squared mass difference ∆ 2 21 = m 2 2 − m 2 1 and to the squared mass sum Σ 2 21 = m 2 2 + m 2 1 . The inset of Fig. 3(b) depicts the correctionsP S e→e −P UR e→e in logarithmic scale.

B. Chiral oscillations and tests for Dirac and Majorana neutrinos
Throughout the main text, we have assumed that the initial neutrino state is a left-handed chiral eigenstate, as a consequence of the properties of weak interactions, in which neutrinos are created. We do not address here questions related to the production of the initial state (such as its localization features).
In proposals for non-relativistic neutrinos detection [16][17][18], it was found that the nature of neutrinos (Dirac or Majorana) influences directly the experiment's measurement rate. At low energies Dirac neutrinos can be detected as left-chiral and negative helicity, while Majorana neutrinos can also be detected as right-chiral and positive helicity [21]. When chiral oscillations are considered, the overall flux of detected neutrinos is halved due to chiral oscillations. As we show, Dirac neutrinos chiral oscillations exhibit a strong effect in the non-relativistic regime: oscillations from leftchiral to right-chiral become more prominent. In the case of Majorana neutrinos, the possible transition probabilities (cf. in the appendix) have the same energy dependence as in the Dirac case, and the survival/oscillation probabilities are equal to those for Dirac particles. Therefore, the measurement of non-relativistic neutrinos will unavoidably be affected by chiral oscillations and are a probe of those effects. Since both Majorana and Dirac neutrinos are subjected to chiral oscillations, in the standard scenario, chiral oscillations and their impact in other relevant quantities, such as absorption rates, can be computed once the nature of the neutrinos is known (see the discussion in [21]). The mass term for Majorana particles can be generalized such that left and right-handed components of the neutrino have different masses [31]. In the case of a vanishing Dirac mass term, but for different left and right Majorana mass terms, we recover the same survival probability and chiral oscillations behavior, as discussed in the appendix. In general case of mixed Dirac-Majorana masses, there will be interferences between the different chiral components, which will lead to further imprints in the survival probability. This most general case is also discussed in the appendix. Furthermore, in the presence of an external magnetic field, the coupled flavor-chiral oscillations are modified due to the non-minimal coupling to the field [32]. Majorana neutrinos exhibit different oscillations properties in the presence of a magnetic field, which can also be used as a way to distinguish Dirac from Majorana particles [23,34]. Works within this framework have focused on relativistic neutrinos, an interesting question is whether such effects are more prominent in the non-relativistic regime.
As for the observability of such an effect we note that, in connection to the aforementioned proposal [16][17][18], neutrinos from the CνB are expected to be extremely non-relativistic and therefore subject to chiral oscillation effects. The current temperature of the CνB T CνB is related to the temperature of the cosmic microwave background (CMB) T CMB via T CνB (4/11) 1/3 T CMB ∼ 0.168 meV. Since neutrinos maintain a Fermi-Dirac like distribution after decoupling [35], one can obtain the root mean square momentum asp 0 ∼ 0.603 meV [17]. Considering the lightest mass as m 1 ∼ 0.1 eV, these values would give the maximum effect persisting for m 1 down to the meV order. A more complete recent study of the effects of chiral oscillations in the CνB can be found in [21].

IV. CONCLUSIONS
In this paper we have described dynamical features of chiral oscillations in the non-relativistic regime, both from a general perspective including their relation with Zitterbewegung, to the study of flavor oscillations for non-relativistic neutrinos. Such framework is especially interesting for upcoming experiments that will be capable to measure the CνB [16], for which chiral oscillations should have a prominent influence.
We considered the free plane wave propagation of an initial left-chiral state, as those generated, for example, via weak interactions. The left-right chiral oscillations are induced by the mass term of the Dirac equation, and we computed the characteristic oscillation length in the non-relativistic limit: in this regime, the survival probability of the initial state could vanish. For masses of the order of eV, the chiral oscillation length is ∼ 10 −6 m.
We then specialized our results to describe chiral oscillations in neutrino propagation by considering a two-flavor problem described within the Dirac bispinor framework. While, as discussed in [13,14,19,20], such corrections are negligible for ultra-relativistic particles, we verified that chiral oscillations are relevant in the non-relativistic dynamical regime. We obtained a modified flavor oscillation formula including the chiral oscillations effects, compatible with known results in the literature [19]. While the standard flavor oscillation term goes with the mass eigenstate energy difference, chiral oscillations depend on the individual energies of the mass eigenstates.
To quantify the effects of chiral oscillations, we computed the averaged survival probability and observed that, for non-relativistic neutrinos, chiral oscillations can affect averaged flavor oscillations up to ∼ 40% of the standard value. This is especially relevant for upcoming CνB tests which will probe non-relativistic neutrinos. Furthermore, Majorana neutrinos exhibit the same chiral oscillations properties of Dirac neutrinos, although the framework in the Majorana case is more enhanced. In particular, we show in the appendix that any bispinor that is its own chargeconjugate (that satisfies the Majorana equation), must have vanishing average chirality. Nevertheless, for a neutrino state initially in a left-handed configuration evolving under the Majorana equation, one gets a survival probability compatible with the Dirac case. A further investigations of those effects in the non-relativistic dynamical regime for a mixed Dirac-Majorana mass deserve a more detailed analysis and can shed a light on other possible distinctions to the usual scenario based on oscillation experiments (see for example [33]). While our results were derived for two flavors, the extension of the formalism to three flavors is straightforward with qualitatively similar conclusions. Moreover, wave packets can be readily considered [14]. In these case, there would be additional effects which also have a non-trivial influence in flavor oscillation. Nevertheless, since chiral oscillations depend on the individual energies of the mass eigenstates, decoherence effects due to a wave-packet treatment, which depend on the difference of neutrino masses, do not influence chiral oscillations. In fact, the latter are intrinsic to the Dirac bispinors, and thus to each mass eigenstate. Finally we would like to point that the formalism adopted here is an effective description: flavor oscillations are correctly described in the framework of quantum field theory [28,29], to which we plan to extend our treatment in the future.
where the mass term L mass depends on the nature of field: Dirac or Majorana. For Dirac fermions, the mass term reads which can be rewritten in terms of the left and right-handed components of the field Ψ L = (1 −γ 5 )Ψ/2, Ψ R = (1 +γ 5 )Ψ/2 as The equation of motion for the total Lagrangian in this case is the Dirac equation whose solutions are given in terms of the Dirac bispinors presented in the main text. The mass term is responsible for coupling the left-handed to the right-handed components of the bispinor, yielding to chiral oscillations. Otherwise, the Majorana mass terms reads where Ψ c = iβα y Ψ * denotes the charge-conjugated field and the left and right handed components have different masses. The above Lagrangian term can be rewritten as [31] where Notice that χ c = χ and ω c = ω. A more general mass term includes both Dirac an Majorana mass terms 3 : While in the main text we have considered the "usual" Dirac equation (the one obtained via the mass term (24)), we now consider briefly quantum oscillations under the Majorana mass term and under the general Majorana-Dirac mass term.

Degenerate Majorana mass term in the bispinor formalism
While Majorana fermions are usually described within the framework of quantum field theory, we will consider from now own the framework of single particle relativistic quantum mechanics and, as in the main text, consider the dynamics of the bispinors. The single particle relativistic quantum mechanics of the Dirac equation with a Majorana mass term has its own peculiarities which we will not discuss here, see for example [23][24][25]36].
In this section we consider m L = m R = m M , for which the equation of motion reads or in the language used through the text where |ψ c = iβα y |ψ * is the charge-conjugated bispinor. For now own, we call the Dirac equation with a Majorana mass term as the Majorana equation.
We notice that (30) decouples into two independent sectors: where The |ψ ± are eigenstates of charge conjugation, and following the literature we call them Majorana bispinors. We have now the following scenario: a Dirac bispinor is an unconstrained bispinor with dynamical evolution given by the Dirac equation; a Majorana bispinor satisfies the Majorana condition (is an eigenstate of charge conjugation) AND has the dynamics given by the Dirac equation (31). It is common to discuss the properties of Majorana bispinors in the so-called Majorana representation of the Dirac matrices, where the Dirac matrices are all imaginary and charge conjugation is just a complex conjugation operation. To better compare Majorana and Dirac bispinors, we keep with the chiral representation of the Dirac matrices. The single particle relativistic quantum mechanics for Majorana particles is obtained by considering the sector related to |ψ + . The Majorana condition imposes a constrain in the form of a bispinor. In particular it relates the upper to the lower (or, in our language the left to the right-handed components) of the bispinor. Let ψ (M ) be a Majorana bispinor satisfying the condition ψ (M ) c = ψ (M ) . We can readily show that, if we write the charge conjugation condition implies that ξ and a general bispinor satisfying the Majorana condition must be of the form which is often called a right-handed Majorana bispinor, or of the form called a left-handed Majorana bispinor. Those bispinors correspond to (27), as can be seen by writing ψ R,L in a two component notation and using the explicitly form of the charge conjugation operator. We conclude that any bispinor satisfying the Majorana condition must have both left and right-handed components, this is a well known and commented through the literature (see e.g. [22]).
The average chirality of a generic bispinor is given by Since this is a real quantity, we have where in the last equality we have explicitly assumed that we are working in the chiral representation (thusγ 5 is real). By recalling then thatγ 2 y = −Î 4 we can write where in the last line we have used that {γ y ,γ y } = 0. Finally, since ψ| c = ( ψ|) * iγ y (becauseγ † y = −γ y ) we have where we have defined γ 5 c = ψ| cγ 5 |ψ c . The above equations brings no new information: the average chirality of the charge conjugated spinor is opposite to the average chirality of the spinor. For example, if |ψ is right-handed, then γ 5 = 1, and thus γ 5 c = −1, in other words, |ψ c is a left-handed spinor. Charge conjugation flips the chirality of a spinor.
For a Majorana bispinor |ψ c = |ψ , thus γ 5 c = γ 5 and from (39) we have Thus any bispinor that is its own charge conjugate must have zero average chirality. Notice that these calculations were performed for generic bispinors, and the only assumption is the Majorana condition. Therefore, if a bispinor satisfy the Majorana condition it has zero chirality at any time, independent of the specific time evolution it follows. This conclusion can also be reached by computing γ 5 explicitly for a bispinor satisfying the Majorana condition, that's it, for a bispinor of the forms (34,35). We now illustrate this point by considering the time evolution of a Majorana bispinor constructed with the Dirac bispinor of the main text. The Majorana bispinor ψ (M ) is thus written as In what follows we use the following relations for the bispinors (see Eq. (5)): thus, for our calculations in the momentum space, the charge conjugation has the effect |u ± (p, m) c = ± |v ∓ (−p, m) .
Adopting the shorthand notation notice that f −p,m = g p,m . We can write the plane wave Dirac bispinor describing the evolution of a Dirac state initially in left-handed and negative helicity, see eq. (6), as The corresponding Majorana bispinor state is given by Here we have used that the charge conjugation includes a complex conjugation, which has to be taken into account in the spatial dependence of the state, which can be grouped with the terms from ψ At t = 0 the state reads while for the Dirac case we had Whether ψ (D) m (0) really represents the initial state of a mass eigenstate depends on the specific weak processes that generates the particle, and will not be discussed here.
The survival probability of (46) is given by Thats it, the survival probability of the considered initial state is the same as in the Dirac case. This is consistent with the conclusions drawn in [21] regarding chiral oscillations: The Majorana neutrino also exhibit quantum oscillations. In this case, the upper right-handed and positive helicity component can oscillate to a left-handed positive helicity, while the lower left-handed and negative helicity oscillates to a right-handed negative helicity. This two oscillation channels yield the survival probability obtained in the above equation. This can be further understood by considering the probability P This probability is half of the survival probability for the Dirac bispinor state (7). Since (45) is a superposition of both positive and negative helicities, the probability that the state is in a right-handed and positive helicity is not null. In fact, such probability is given by We therefore notice that Finally, the average chirality of (45) vanishes: as discussed in the beginning of this section, the Majorana condition implies that the bispinor is a superposition of left and right chirality components, thus its average chirality vanishes.

Quantum oscillations of left-handed chiral eigenstates
The calculations above are an illustrative example of the dynamics of a Majorana bispinor obtained with the corresponding Dirac bispinors. We can instead turn our attention to the formalism of [26,27] to evaluate the time evolution of an initially left-handed state.We generalize the result in the aforementioned references and consider through this section the general Dirac-Majorana mass term (28), which can be write in terms of the bispinors χ = Ψ L + Ψ c L and ω = Ψ R + Ψ c R as In other words, considering we are interested in the dynamics of |ψ L (t) under the general mass term (28) 4 . This Lagrangian can be rewritten as where [31] We then write η i (i = 1, 2) in terms of their right and left-chiral components and considering negative helicity eigenstates, we get from the equations of motion: which gives the eigenenergies E To recover the dynamics of |ψ L (t) , we backtrack the problem by using the relations (56) and then Eqs. (27), such that |ψ L (t) = cos (φ) sin (θ 1 )e −iEp,m 1 t |f + 1 + cos (θ 1 )e iEp,m 1 t |f − 1 + sin (φ) sin (θ 2 )e −iEp,m 2 t |f + 2 + cos (θ 2 )e iEp,m 2 t |f − 2 , and therefore |ψ L (t) = A L (t)|ψ L + A R (t)|ψ R + A c L (t)|ψ c L + A c R (t)|ψ c R .
(62) 4 Notice that the connection with the notation used in eq. (28) is: The time dependent coefficients are given explicitly by A L (t) = cos 2 (φ) [cos (E p,m1 t) + i cos (2θ 1 ) sin (E p,m1 t)] + sin 2 (φ) [cos (E p,m2 t) + i cos (2θ 2 ) sin (E p,m2 t)] , A R (t) = i sin (2φ) 2 [sin (2θ 1 ) sin (E p,m1 t) − sin (2θ 2 ) sin (E p,m2 t)] , A c L (t) = − 1 2 sin(2φ) [cos(E p,m1 t) + i cos(2θ 1 ) sin(E p,m1 t)] + 1 2 sin(2φ) [cos(E p,m2 t) + i cos(2θ 2 ) sin(E p,m2 t)] , A c R (t) = −i cos 2 (φ) sin(2θ 1 ) sin (E p,m1 t) + sin 2 (φ) sin(2θ 2 ) sin (E p,m2 t) . (63) The survival probability of the state is given by while the average chirality is calculated by using that ψ L |γ 5 |ψ L = ψ c L |γ 5 |ψ c L = −1 and ψ R |γ 5 |ψ R = ψ c R |γ 5 |ψ c R = 1, such that For a pure Dirac mass term, m R = m L = 0, the explicit temporal evolution (62) simplifies to ψ (D) where A(t) = cos (E p,mD t) + i p E p,m sin(E p,mD t),         L (t) are equal. Those dynamical features are also the same in the more general case of m D = 0 and m R = m L . This is expected: for a vanishing Dirac mass, the equations of motion are fully decoupled in terms of the bispinor χ and ω, and the dynamics of a state initially in a left-handed chirality is determined entirely by the sector related to χ, which only induces oscillations between |ψ L and |ψ c R . In the case m D = m L = 0, there are no quantum oscillations. When both Dirac and Majorana masses are present, the survival probability and the average chirality display a more intricate dynamics. In this case, there will be interferences due to the different values of the masses m 1,2 which drive different oscillation channels. Further investigation of such effects will be addressed in a future work.