Practical Dirac Majorana confusion theorem: issues and applicability

We inspect the model-independent study of practical Dirac Majorana confusion theorem (pDMCT) – a wide spread belief that the difference between Dirac and Majorana neutrinos via any kinematical observable would be practically impossible to determine because of the difference only being proportional to the square of neutrino mass – in context of processes that have at least a neutrino antineutrino pair in their final state. We scrutinize the domain of applicability of pDMCT and also highlight those aspects that are often misunderstood. We try to clarify some of the frequently used concepts that are used to assert pDMCT as a generic feature irrespective of the process, or observable, such as the existence of any analytic continuity between Dirac and Majorana neutrinos in the limit mν→0.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\nu } \rightarrow 0.$$\end{document} In summary, we illustrate that pDMCT is not any fundamental property of neutrinos, instead, it is a phenomenological feature of neutrino non-observation, depending on models and processes.


Introduction
Are neutrinos distinct from their antiparticles like the rest of the known fermions of the Standard Model (SM), or are the neutrino and antineutrino quantum mechanically identical to one another?An affirmative response to the first (second) question would imply neutrinos are Dirac (Majorana) fermions.We are yet to have a definite answer to this fundamental question regarding Dirac or Majorana nature of neutrinos.However, the literature is replete with attempts made on both theoretical and experimental fronts, without much success.The situation is such that there is an ossified belief in the community that the difference between Dirac and Majorana neutrinos via any kinematical observable would be practically impossible to determine due to fact that the observable difference between Dirac and Majorana neutrinos is proportional to the tiny neutrino mass.This is often cited as the "practical Dirac Majorana confusion theorem" (pDMCT, in short) [1] 1 .While the "theorem" has been verified in some cases, there is no general model-independent and process-independent proof.Also there is a general lack of clarity regarding its domain of validity.It is therefore necessary and important to explore whether there are any SM allowed processes and kinematic observables that can directly probe the Majorana nature of neutrinos avoiding this pDMCT constraint.In this article we discuss the domain of validity of pDMCT as well as its exceptions.This does not a E-mail at: cskim@yonsei.ac.kr 1 A few notable precursors to the formulation of pDMCT were the analysis made by Refs.[2,3].invalidate pDMCT but brings more clarity with regard to its applicability as an useful tool.
There are two issues regarding pDMCT to which we would like to draw readers' attention.1.The pDMCT should not be taken out-of-context of its historical development.Historically, only SM allowed neutral current interaction mediated processes as well as those processes mediated by exchange of massive Majorana neutrinos [4], were analyzed.In processes involving neutral current interaction, there had been no way to gain any information regarding individual neutrino antineutrino energies or 3-momenta.This invariably leads to integration over neutrino antineutrino related kinematic variables while proposing any relevant kinematic observables.If one considers a process which is mediated not through weak neutral current interactions, and if one has access to individual information of neutrino and antineutrino momenta without directly measuring them, then one need not take pDMCT for granted while analyzing the relevant observables.2. The usual approach to validate pDMCT as a general theorem is by alluding to a non-existent correspondence between massive Dirac and Majorana fermions in the massless limit where the neutrinos have specific chirality.This strangely overlooks the well-known mathematical impossibility of having a chiral massless Majorana neutrino.In any case it does not make practical sense since the neutrinos have non-zero mass.
For the first time in Refs.[5] and [6], we respectively implemented both model-independent and specific process-independent studies of pDMCT.The current work borrows some of the features from both the above papers in a more accessible manner and directly addresses some pedagogical aspects to bring clarity.Thus in Sec. 2 we first consider a generic model-independent analysis of processes that contain a neutrino antineutrino pair in the final state.This structure in the final state allows for application of Pauli exclusion principle through anti-symmetrisation in the case of Majorana neutrinos.The application of the exclusion principle is independent of the size of the non-zero mass of the neutrino or any other dimensional parameter for that matter.We refer to this as model-independent in the sense that our analysis includes both the SM and new physics (NP) contributions.However, the process itself is allowed in the SM.This is exclusively discussed in Sec. 3. In Sec. 4 we highlight the domain of applicability of pDMCT.This is followed by some pedagogical explanations on pDMCT in Sec. 5. Finally we conclude highlighting the important features in Sec. 6.

A model-independent analysis of processes containing ν ν in the final state 2.1 Details on the process under consideration
Consider a general process with a neutrino and an antineutrino2 of the same flavor in the final state, say where X, Y can be single or multi-particle states, Y can also be null, the contents of X and Y (if it exists) are visible particle/s and the 4-momenta p X , p Y are assumed to be well measured so that one can unambiguously infer the total missing 4-momentum of ν ν, p miss = p 1 + p 2 .The 4-momentum of X must either be fixed by design of the experiment (e.g.X might be a particle produced at rest in the laboratory or be the constituent of a collimated beam of known energy or it could consist of two colliding particles of known 4-momenta) or the 4momentum of X be inferred from the fully-tagged partner particle with which it is pair-produced.The final state Y should not contain any additional neutrinos or antineutrinos.The process could be a decay or scattering depending on whether X is a single particle state or two particle state.Some actual processes that satisfy such criteria are e + e − → ν ν, A word of caution: the process X → Y ν ν is not necessarily a neutral current process, and could proceed through other means such as by doubly weak charged currents.To keep our analysis model-independent we allow the process X → Y ν ν to proceed even via NP interactions.We do not consider any specific NP possibility, but simply ensure that whenever explicit NP contributions are needed there are no Lorentz-symmetry violation as well as CPT violation in the underlying effective Lagrangian.
It should be noted that in this work we discuss processes where the effect of measurements does not destroy the identical nature of Majorana neutrino and antineutrino.This is akin to putting the constraint that in a double-slit experiment, meant to observe the interference of light, no measurement should identify the slit through which the photon has passed.We again stress that direct observation of the state of neutrinos (via neutrino detector) destroys the effect of quantum statistics since the neutrinos after the detection are projected on to specific helicity states (and thus no longer remain indistinguishable) by observation.As is well known, quantum statistics requires absolutely identical indistinguishable particles.

Origin of observable difference between Dirac and Majorana neutrinos and practical Dirac Majorana Confusion Theorem (pDMCT)
We now recall some features from our earlier work for completeness.The transition amplitude is, in general, dependent on all the 4-momenta.For brevity of expression and without loss of generality, we denote the transition amplitude by only mentioning the p 1 , p 2 dependence.For Dirac neutrinos, the transition amplitude for X → Y ν ν can be written as, while for Majorana case the amplitude is anti-symmetrized 3with respect to the exchange of the Majorana neutrino and antineutrino which are quantum mechanically identical fermions, where 1/ √ 2 takes care of the symmetry factor.Note that the amplitudes of Eqs. ( 1) and (2) do not necessarily assume the SM interactions, they can involve NP effects as well, and hence they include the most general structures of the amplitude that are allowed by Lorentz invariance.
The difference between Dirac and Majorana cases that can possibly be probed is obtained after squaring the amplitudes (including the usual summation over final spins of ν, ν and averaging over initial spins 4 ) and taking their difference, which is given by, From Eq. (3) it is easy to conclude that there are essentially two major sources of any possible difference between Dirac and Majorana cases: 1. Unequal contributions from "Direct term" and "Exchange term" in general, i.e. ( It is interesting to note that in the case of the SM the interference term always depends on the size of the neutrino mass, that is In presence of NP contributions, the full interference term need not follow Eq. ( 6).
The above sources of difference between Dirac and Majorana cases at the level of amplitude square, may or may not survive at the level of observables 5 which requires appropriate phase space considerations.We note that in the case when no individual information about ν ν are either known or deducible, the only difference between Dirac and Majorana cases that can be experimentally accessed is obtained after full integration over p 1 and p 2 which gives, which is directly proportional to m 2 ν if only the SM interactions are considered.Here we have used the fact that although, in general, the "Direct" and "Exchange" terms differ as shown in Eq. ( 4), when we fully integrate over the 4-momenta of neutrino and antineutrino we get as p 1 and p 2 act as dummy variables since the range of integration is identical.Therefore, if we are able to find any smart method to deduce (or infer) energy-momenta of missing neutrinos, in order not to make full phase space integral over missing neutrinos, and if there exists non-trivial difference between the "Direct term" and the "Exchange term" as shown in Eq. ( 4), the pDMCT would be nicely avoided.
In the simple processes with the SM mediated interaction alone (e.g.weak neutral current mediated decay Z ( * ) → νν) one finds that, (1) the "Direct term" and "Exchange term" are equal (see section 2.3), (2) the "Interference term" is proportional to m 2 ν and (3) the observable usually requires full phase space integration over p 1,2 .This leads to the conclusion that all kinematical observable differences between Dirac and Majorana cases would be proportional to m 2 ν .This is essentially the statement of the "practical Dirac-Majorana Confusion Theorem" (pDMCT).
Note that our model-independent and process-independent analysis suggests that if (1) the doubly weak charged current processes 6 that possibly lead to a non-zero difference between "Direct term" and "Exchange term" of Eq. ( 3) and (2) some kinematic configurations could be identified where individual information about ν, ν can be accessed so as to avoid doing the full phase space integration in Eq. ( 7), then one might avoid pDMCT constraint.

Z ( * )
→ νν in the SM and special cases of There are certain special cases when Eq. ( 4) is not satisfied.These are (a) collinear case: p 1 = p 2 , (b) symmetric case: M (p 1 , p 2 ) = M (p 2 , p 1 ) and (c) antisymmetric case: M (p 1 , p 2 ) = −M (p 2 , p 1 ).As an example, within the SM for the neutral current mediated processes Z ( * ) → νν, such as Z → νν, e + e − → νν, B → Kνν and etc. [5] Therefore, for these processes we find that even at the level of amplitude square, which suggests that in such a case pDMCT holds true always without any exception.Since in this case pDMCT holds at the amplitude square level, it naturally holds true for all observables 7 .See section 3 to find out how this conclusion changes in presence of NP. 6 For the sequential weak charged current mediated decays that produce neutrino and antineutrino of different flavors, e.g.ℓ − → ν ℓ ν ℓ ′ ℓ ′− , where ℓ, ℓ ′ ∈ {e, µ, τ} and one can never have ℓ = ℓ ′ , the ν ℓ and ν ℓ ′ can never be considered as identical fermions even if they might indeed be Majorana fermions.We do not consider such processes in our analysis. 7The equality of direct and exchange terms, either at amplitude square level as shown in Eq. ( 9) or at the level of experimentally measurable observable that involves full phase space integration shown in Eq. ( 8), is often generalized as existence of one-to-one correspondence between Dirac and Majorana neutrinos in the massless limit m ν → 0. See section 5.1 for more details.

New Physics scenarios and pDMCT in
There is no reason a priori for the "practical DMCT" to hold, if NP contributions in the neutrino interactions are allowed, as in this case Eq. ( 3) "Direct" and "Exchange" terms in general do not need to cancel each other.To illustrate it more clearly using symmetry properties of the transition amplitude, let us assume that some (yet unknown) NP at high energy modifies the low energy effective neutrino interactions Z ( * ) → ν ℓ ν ℓ .
Considering Lorentz invariance, CP and CPT conservation, applying Gordon identities as well as neglecting any m ν dependent terms at the amplitude level, we find that the most general decay amplitude for Z(p) → ν(p 1 ) ν(p 2 ) is as follows (for Dirac neutrinos) [5], where g Z = e/(sin θ W cos θ W ) with θ W being the weak mixing angle and e being the electric charge of positron, and for different lepton family ℓ = e, µ, τ we have the possibility of having different vector and axial-vector coupling parameters C ℓ V , C ℓ A .Since we are considering NP possibilities here, we can write the vector and axial-vector coupling parameters as follows, where ε ℓ V , ε ℓ A parameterise the NP effects, vanishing in the SM case.The amplitude for Majorana case is given8 by It is clear that we can combine the direct and exchange amplitudes in this case and effectively redefine the vertex structure for Z → ν ℓ ν ℓ when Majorana neutrinos are considered.Keeping neutrino mass dependent terms in the amplitude squares, we get different results for Dirac and Majorana neutrinos: such that for the SM alone where we have kept only the leading order contributions of ε ℓ while considering NP effects.It is clear that the SM result is fully in agreement with "practical Dirac Majorana confusion theorem" even at the amplitude-squared level, i.e. in the limit m ν → 0 there is no observable difference between Dirac and Majorana cases in the SM via the process Z → ν ν.It implies within the SM.However, the difference between Dirac and Majorana neutrinos appears in context of NP contributions even when one neglects m ν dependent terms (unless, of course, A in which case the additional NP contributions effectively rescale the SM allowed V − A coupling).Possible example of NP effects in this Z boson decay could arise from kinetic mixing of Z with the neutral gauge bosons from extra gauge groups like additional U(1) or S U(2) R .In this work we are not concerned with any specific model of NP to keep our results and discussions very general.
Note that for the neutral current mediated processes in the SM, such as Z → νν, e + e − → νν, B → Kνν etc., the deduction of energy and/or momentum of the invisible neutrinos does not help to distinguish between Dirac and Majorana neutrinos because the pDMCT holds true even at the amplitude-squared level.In such a case every conceivable observable gives the same result for Dirac and Majorana neutrino, except showing the tiny difference coming from the interference term which is proportional to the neutrino mass.Therefore, only in the presence of physics beyond the SM [5,10,11] a distinction may be made between Dirac and Majorana nature of neutrinos through such processes.

practical Dirac-Majorana Confusion Theorem and its exceptions 4.1 The general strategy to probe nature of neutrinos and pDMCT
The general formalism discussed in previous sections may be suitably illustrated by the chart shown in Fig. 1.The red (color in online) arrows drawn in the figure present the previous works of the pDMCT, all of which have studied the weak neutral current processes within the SM, such as  etc.As can be easily seen, all of these processes have confirmed the pDMCT because (1) due to the weak neutral current processes within the SM and/or (2) there is no way to observe or deduce the 4-momenta of ν and/or ν as these are simple 2-or 3-body processes.
Although there is no general, model-independent, processindependent and observable-independent proof9 of the "practical Dirac Majorana confusion theorem", it is generally assumed to apply to all the probes of Majorana nature of sub-eV neutrinos.The formalism presented in sections 2 and 3 provides a simple model-independent, process-independent and observable independent view of the pathways by which the confusion theorem can be overcome, such as by using the properties of the NP interactions or analysing the "special kinematical scenarios" utilising the chosen parts of neutrino momentum spectra.
In order to illustrate the general strategy within the SM we would like to point out that the 2-body and 3-body νν final state processes are not suitable for the special kinematic scenarios because those decays are weak neutral current processes.On the other hand, 4-body decays such as B, D, K, H, J/ψ, Υ(1s), .. → µ + µ − ν µ ν µ , which are doubly charged weak decay processes [6], could be more relevant for utilising the dependence of decay distributions on kinematic variables to distinguish between Dirac and Majorana neutrino in the case of the SM-like interactions.The main advantage of 4-body decays over the 2/3-body decays is the multitude of kinematic configurations and related observables that can be explored for the purpose of distinguishing between Dirac and Majorana neutrinos.The difference between Dirac and Majorana neutrinos that exists at the level of amplitude square, requires a proper observable so that we can access that difference.This once again highlights how crucial an observable is in this context.In the next subsection we summarize our findings on how to possibly overcome the pDMCT to distinguish between Dirac and Majorana neutrino using the special kinematics in B 0 → µ + µ − ν µ ν µ decay in the SM.

How to overcome pDMCT within the SM by using the special kinematics in
In most of the experimental scenarios, especially true for processes of the form X → Y ν ν, information about individual neutrino momenta is not available.In such a case the difference between the Dirac and the Majorana neutrinos, that may be realised, is given by the integrated interference term in Eq. ( 7).In such a case the evaluation of the squared Feynman diagram for the "Interference term" in the SM necessarily involves two helicity flips which would make it proportional to m 2 ν .Thus, if only the SM interactions are considered and one fully integrates over the neutrino and antineutrino 4-momenta due to unobservable neutrinos , the difference between Dirac and Majorana cases is proportional to m 2 ν .This may be considered as the most general statement of the "practical Dirac Majorana confusion theorem".
However, there is no reason a priori for the pDMCT to hold, if one can consider special kinematic configurations where the 4-momenta (or some components of the 4-momenta) of the neutrino and antineutrino are known, so that full integration over 4-momenta is not necessary before comparison with the experiment.This is valid even in the SM.In the following subsubsections we discuss this scenario where such situation may occur.

pDMCT and the doubly weak charged current decay
The decay B 0 → µ − µ + ν µ νµ takes place via doubly weak charged currents since flavor changing neutral currents are impossible at tree-level in the SM.The branching ratio of this mode gets substantial contributions from intermediate resonances such as π − and D − .Details on the process have been thoroughly investigated in Ref. [6].Similar 4-body decays such as B, D, K, H, J/ψ, Υ(1s), .. → µ + µ − ν µ ν µ could be studied in an analogous manner.
The Eqs. (31, 32) in Ref. [6] clearly show the unequal contributions from "Direct term" and "Exchange term", which satisfy our Eq.( 4) Exchange term , unlike the weak neutral current processes.As explained in section 2.2 if we can measure or deduce the individual energy or momentum of the missing neutrino, then the pDMCT constraint will not apply.
In Fig. 2, we show the general kinematics of B 0 → µ − µ + ν µ νµ in the rest frame of B. The angles θ n and ϕ are indeed inaccessible in general, as the neutrino pair goes missing.Therefore, for a physically useful differential decay rate we must integrate over both θ n and ϕ, i.e.
where Y is the magnitude of 3-momentum of the di-muon system (muon pair with invariant mass m µµ ) or di-neutrino system (neutrino pair with invariant mass m νν ) in the rest frame of B 0 , Y m is the magnitude of 3-momentum of µ ± in the di-muon rest frame, Y n is the magnitude of 3-momentum of ν µ or ν µ in the di-neutrino rest frame.It is straightforward to show that the difference between Dirac and Majorana cases, as shown in Eqs.
(33, 34) of [6], is given by which agrees with the pDMCT.Therefore, even in the case of a doubly weak charged current mediated decay process if we integrate fully over the available phase space of the invisible neutrino pair, we confirm pDMCT as expected from our discussion in Sec.2.2.The situation changes, if and when we access a special kinematic scenario where the individual energy or 3-momentum of the invisible neutrinos can be inferred.In the following we discuss such a special kinematic situation.

Back-to-back muon special kinematic configuration in
Consider the decay of the parent B 0 in its rest frame in which the µ + , µ − back to back with equal but opposite 3-momenta.Experimentally this is an ideal situation since it is easier to detect muons.The neutrino antineutrino pair must also fly away back-to-back since 3-momentum is conserved.This is a much simpler kinematic configuration than the general kinematics for any 4-body decay.Instead of the usual five independent variables one needs to describe any 4-body decay, we only need two independent variables to describe the back-to-back configuration.In this case, the energies of the two muons are the same and let us denote them by E µ .Similarly, the energies of the back-to-back neutrino and antineutrino are the same and let us denote them by E ν .Either E µ or E ν is independent, because from conservation of energy we get, where m B is the mass of the B 0 meson.Let us choose E µ as one independent variable.The other independent variable would then be the angle, say θ, between the muon direction and the neutrino direction.We would like again to point out that direct observation of neutrinos destroys the effect of quantum statistics, as explained before.Hence, any inference from qauntum statistical probe necessarily requires no direct detection of the final neutrinos [17].This is precisely achieved in the back-toback special kinematic configuration where the energy momentum information is inferred purely by observing the muons.
For back-to-back case, with E 1 = E 2 = E ν (say) and the angle between the two neutrinos Θ = π, we get the following, Moreover, for the back-to-back case we have Here X m and X n denote the muon pair and the neutrino pair.
It can be shown that, in general, Whenever E 1 = E 2 for any value of the angle Θ between the neutrino and antineutrino we get cos θ n = 0.By analytic continuation we extend this feature to the back-to-back kinematics for which the cos θ n has a discontinuity otherwise.Moreover, in the back-to-back case we have both the back-to-back muons and the back-to-back neutrino antineutrino pair, in one single plane.This implies that for the back-to-back case we have ϕ = 0 [18] .These choices put the orientation of the coordinate axes in such a way that the back-to-back neutrino and antineutrino fly away defining the x-axis.The xz-plane in Fig. 2 is the one in which the 3-momenta of muons lie, and now the back-to-back neutrino antineutrino define the x-direction.The direction perpendicular to the neutrino direction is the z-direction.
The back-to-back kinematics specified by ⃗ p + + ⃗ p − = ⃗ 0 = ⃗ p 1 +⃗ p 2 has 3 additional constraint equations, which implies that only 2 out of the initial 5 variables would remain independent in back-to-back case, i.e. 3 variables should either be fixed by the back-to-back conditions or they must be dependent on the 2 variables which are independent (ie.E µ , cos θ m ) [18].If we define the angle between the neutrino and muon directions to be θ, then θ m = π/2 − θ.This implies that cos θ m = sin θ. ( The differential decay rate in the back-to-back case is therefore given by, where M D/M ↔ 2 is same as M D/M 2 with the necessary dot product substitutions in the back-to-back case (this is the meaning of the subscript '↔').Please note that the difference between the integrated widths of Γ D ↔ and Γ M ↔ can be very large as shown in Eq. ( 51) of [6], and computable in the SM or any other framework.It is also independent of the magnitude of the unknown neutrino mass for the leading terms.
For simplicity we neglect the masses of muons and neutrinos in comparison with the mass of B 0 as well as the energies.Note again that this does not mean that we consider muons and neutrinos to be massless.With this condition we find that only the non-resonant contributions survive.We consider only the dominant form factor contribution, and assume it to be a constant form factor.The full differential back-to-back decay rates are then given by, where µ is the magnitude of the 3-momentum of the back-to-back muons.There are no m ν dependent terms here.The muon energy distribution obtained after integrating over sin θ shows that there exists non-zero difference between the Dirac and the Majorana cases.Moreover, the corresponding branching ratio for the back-to-back kinematics for Majorana case is more than 15 times bigger than that for the Dirac case.Thus, these results are not in agreement with pDMCT.Strictly speaking, since we are not integrating over the full phase space of neutrinos, the pDMCT need not apply in this case.The backto-back kinematic configuration provides a way of realising this exception.For more details of the back-to-back kinematics and the related issues, please see Appendix A.
This result confirms the discussion in Sec. 2. However, it should be acknowledged that the reader may find the result contradictory or even counter intuitive to the previous understanding of the pDMCT.In the next section we try to clarify some of the frequently used concepts that are used to assert pDMCT as a generic feature irrespective of the process, or observable.

Discussions on concepts usually accompanying explanations of pDMCT
It is generally believed that all observable difference between Dirac and Majorana neutrinos must always be proportional to some power of neutrino mass m ν [1], which is the content of practical Dirac-Majorana Confusion Theorem (pDMCT).However, all processes where the theorem was shown to hold involved either full integration over the 4-momenta of missing neutrinos and/or only for the weak neutral current process within the SM [5,6].
As we have discussed, pDMCT is not a fundamental theorem of neutrinos: pDMCT actually depends on physics models, processes and observables, e.g. even for Z → νν, pDMCT holds within the SM, but can be violated beyond the SM depending on the model parameters.Even within the SM, pDMCT depends on the processes, e.g.B → Kνν confirms pDMCT, but B → µ + µ − νν can violate pDMCT.Therefore, while the quantum statistics of Majorana neutrinos10 is a fundamental property of neutrinos, pDMCT is not.Instead, pDMCT is an emergent phenomenological feature arising out of non-observation of neutrinos.
In this section we give comments on the existence of any analytic continuity between Dirac and Majorana neutrinos in the limit m ν → 0 and the issue on anti-symmetrization of amplitude while dealing with pair of identical Majorana neutrinos.We also address in the appendix some pedagogical issue on massless Majorana neutrino, which can be a fundamental difference between Dirac and Majorana neutrinos.

Is there any one-to-one correspondence between Dirac and Majorana neutrinos in the massless limit,
The issue of one-to-one correspondence between Dirac and Majorana neutrinos in the massless limit m ν → 0 can be analyzed from the context of specific processes and observables.As mentioned in Sec.2.3 if one considers neutral current mediated processes such as Z → ν ν (and include summation over neutrino spins while evaluating amplitude squares) within the SM, the direct and exchange terms in amplitude square become equal and the difference between Dirac and Majorana neutrinos becomes proportional to m 2 ν , which vanishes if we were to simply apply the limit m ν → 0 at the end. 11However, when one considers processes that are not facilitated by the SM neutral current interactions, the direct and exchange terms can have non-zero difference even when we neglect m ν dependent terms (or equivalently put m ν → 0).In such a case some specific observable might be able to probe these important non-zero differences.In these instances, in the context of a specific process and observable, there is indeed no one-to-one correspondence between the Dirac and Majorana neutrino scenarios.However, in all cases with the SM only interactions if the observable includes full phase space integration over the neutrino and antineutrino, we do find that the direct and exchange terms have equal contribution after integration (see Eq. ( 8)) which amounts to no observable difference between the two scenarios in the limit m ν → 0.
As explained in Appendix B, it is a mathematical impossibility to preserve Majorana nature of a fermion when its mass becomes zero.In fact due to Lorentz invariance and conservation of chirality for massless fermions, such chiral fermions have a distinct nomenclature as being Weyl fermions.To describe Weyl fermions it is sufficient to use 2-component complex spinors instead of 4-component complex spinors.Nevertheless, the 2-component Weyl spinors can be used to construct 4-component complex Dirac spinors as well as 4-component Majorana spinors (that are real in the Majorana basis).Both the Dirac and Majorana spinors have both left-and right-chiral components.When one takes the massless limit or when one considers ultra-relativistic fermions one finds that these constructs of Dirac or Majorana spinors prefer specific chirality states.Nevertheless, as long as the mass of the fermion is nonzero both the chiral states are present.However, once the fermion is massless, it becomes fully chiral and it can not have Majorana nature at all.The Dirac nature (meaning its particle and antiparticle states are distinct and distinguishable) survives the massless limit.Therefore, although both Dirac and Majorana 4-component spinors get reduced to 2-component Weyl spinors in the massless limit, the Weyl spinors only show Dirac nature and the Majorana nature is completely lost.This implies there is really no one-to-one correspondence between Dirac and Majorana nature of neutrinos in the massless limit.

Majorana neutrinos of the same flavor with m ν > 0?
It is well known that when two identical particles are present in the final state, the corresponding transition amplitude needs to be symmetrized (or anti-symmetrized) with respect to their exchange if they are bosons (or fermions).Therefore, if a final state has two massive neutrinos or two massive antineutrinos of the same flavor (i.e.ν ℓ ν ℓ or ν ℓ ν ℓ , with ℓ = e, µ, τ) then the transition amplitude would always be anti-symmetric under their exchange (which involves exchange of 4-momenta and spin) irrespective of whether they are Dirac and Majorana fermions.This is to ensure the Fermi-Dirac statistics.However, if we have a final state that has ν ℓ ν ℓ , then it has distinct particles for Dirac neutrinos, but it has identical particles for massive Majorana neutrinos.Thus, when considering Majorana nature of the massive neutrinos, one needs to anti-symmetrize the transition amplitude in this case.This is one of the main differences between Dirac and Majorana neutrinos, and it has been noted by many authors, (see eg. [7,8,13] and etc.) before us.One simple example of this amplitude anti-symmetrization is, as shown in Sec. 3, the most general amplitude of Z → νν for Dirac neutrino in Eq. ( 11) and for Majorana neutrino in Eq. (13).
If the 4-momenta (and spins) of ν ℓ , ν ℓ be denoted by p 1 , p 2 (and s 1 , s 2 ) respectively, then the transition amplitude for Dirac case can be symbolically written as, while the amplitude for Majorana case would be, In the calculation of the observable for the specific process in both Dirac and Majorana cases, one takes the square of the amplitude, does the usual trace calculations by summing over the final spins and averages over the initial spins, except when one is interested in an observable that depends on neutrino spins which is practically impossible to do experimentally for sub-eV active neutrinos.Thus the spin information gets wiped out via the spin summation.In context of the SM, we know that the V − A nature of weak interaction ensures that we always get a left-chiral neutrino and a right-chiral antineutrino.However, despite being produced in specific chiral states, their chirality is not conserved due to non-zero mass, following Eq.(38).Due to non-zero mass, chirality is not same as helicity which is the projection of spin along the direction of motion.Thus, left and right helical massive neutrinos get produced from the SM weak interaction.Thus we consider all spin possibilities of the massive eigenstates in our calculation.Any issues related to relativistic or non-relativistic kinematics are automatically taken care of by the field theoretic calculations for amplitude square with summation over final spins and average over initial spins.
Another approach to include all spin possibilities in the final calculation leading to correct amplitude square is via splitting the full amplitude into all possible helicity amplitudes, where one specifies the individual helicities, say λ 1 and λ 2 instead of s 1 and s 2 .In such a method, one has to exchange the helicities (equivalent to exchange of spins) for the Majorana case, i.e.
In the Majorana case amplitude square, it thus becomes clear that there will be interference terms that require helicity flip and such terms would be proportional to m 2 ν .If the SM interactions are only taken into account, then all the interference terms always involve helicity flips.

Conclusions
In this work, we have revisited the practical Dirac Majorana confusion theorem and studied its domain of applicability.We find that one should always keep the historical context of neutral currents in mind while applying this theorem rigorously.If the process involves doubly weak charged currents, or some new physics contributions, and if one can infer the energy or 3-momentum of neutrino and antineutrino using some special kinematic configurations, then this theorem need not hold true.It might hold true in specific processes, but this theorem does not have any generic, model-independent, process-independent and observable-independent proof.We have also highlighted and addressed some of the most commonly and easily misunderstood concepts that come to mind while thinking of this theorem.
As a final note, it is rather tempting to confirm the pDMCT and/or find out how to overcome the pDMCT from the fundamental Lagrangian level.The effective interaction Lagrangian always respects quantum statistics even though it might not be evident at the level of fundamental interaction Lagrangian.Since the Lagrangian pertaining to neutrino masses have neither a direct bearing on the effective interaction Lagrangian nor carry any signature of the quantum statistical difference we are interested in, the mass generating Lagrangians do not affect our analysis.As is previously explained, the use of basic weak neutral current interaction, Z → νν, will always lead to pDMCT within the SM.And the use of weak charged current processes, such as W ± → ℓ ± ν ℓ and ℓ − → ν ℓ ν ℓ ′ ℓ ′− , do not introduce any identical Majorana neutrino pair, so no difference between Dirac and Majorana neutrinos can be probed using these.Only doubly weak charged 4-body decay processes, e.g.B, D, K, H, J/ψ, Υ(1s), .. → µ + µ − ν µ ν µ , could be considered to find out whether the process gives equal contributions from "Direct term" and "Exchange term" -the only meaningful way to overcome pDMCT within the SM.
As a final conclusion, I would like to stress that quantum statistical effects can be used to pin down the nature of the neutrino -whether it is Majorana or Dirac type.In the context of processes allowed in the Standard Model, this follows when 1. the direct and exchange terms in amplitude square are nontrivially different (see Section 2.2), and 2. the observable defined in terms of neutrino momenta is accessible even when the neutrinos are not directly observed in the final state (see Section 4.2).
The neutrino-less double beta decay (0νββ) [19,20] has a limitation that it is dependent on the unknown tiny mass of the neutrino.If it is too small, there is no possibility of establishing the nature of the neutrino through 0νββ.Our proposal to probe quantum statistics of Majorana neutrinos seems to be the only viable alternative to 0νββ as far as probing Majorana nature of sub-eV active neutrinos is concerned.

A Clarifications on back-to-back kinematics
In back-to-back kinematic configuration of Fig. 2, we note the following points.
1.All the final particles fly away in a single decay plane in the rest frame of the parent particle, i.e. ⃗ p + , ⃗ p − , ⃗ p 1 and ⃗ p 2 lie on a single plane, and 2. Only by assuming all the final particles in Fig. 2 to be massless, or by neglecting their masses, do we get equal energies for the particles flying back-to-back, i.e.
and from conservation of energy Thus knowing either E ν or E µ is sufficient.3. The back-to-back kinematics for a measured event would specify E µ as well as the angle θ shown in Fig. 2. Since ν µ and νµ are invisible in the detector, the angle θ is experimentally unknown and therefore should be integrated out for the final observable.4. The back-to-back configuration is a special case of the general kinematic configuration, and not arrived at by any integration or summation.The general kinematic configuration involves two decay planes (see Fig. 2), and requires five independent variables for complete specification (see Sec.IV.E of [6]).For full specification of the back-to-back kinematics one instead needs to specify, only the energy E ν or E µ (here we are making the massless assumption mentioned above) and the angle θ in Fig. 2. To come from general kinematics to the back-to-back kinematics one needs to fix certain quantities.

A.1 Important issues to address when coming to back-to-back kinematics from general kinematics
The angle between two decay planes: The usual description of general kinematics for the 4-body final state has two decay planes, with an angle ϕ between them.For back-to-back kinematics the two decay planes coincide to form a single decay plane.So for back-to-back kinematics ϕ = 0 [18].No integration over ϕ is involved to arrive at the final observables in back-to-back kinematics.

Discontinuity from general kinematics to back-to-back kinematics:
The expression has a discontinuity when E ν → E ν and the angle Θ νν between ν µ and νµ approaches π.In fact this expression yields 0/0 form if we simply substitute E ν = E ν and Θ νν = π.Thus, one needs to apply L'Hospital's rule to get the limit.The limit E ν → E ν with Θ νν = π yields ±1, i.e. cos θ ν has a discontinuity which needs to be resolved.Here we note that for E ν = E ν we get cos θ ν = 0, as long as Θ νν π.At Θ νν = π this discontinuity appears and it can be resolved by taking average of the two limits, which yields 0 and makes cos θ ν continuous.The approach is similar to what is done to remove the discontinuity in Heaviside step function at x = 0. Once cos θ ν and ϕ are put to zero, the rest of the results needs to be consistent with expectations from helicity arguments, as is found consistently in Ref. [6].

Inferred neutrino energy distribution:
The neutrino energy distribution in back-to-back configuration requires that the neutrino energies be inferred.In the case of back-to-back events we have In our paper [6] in Eq. ( 35), we have shown that the difference between Dirac and Majorana neutrinos vanishes when full integration over neutrino phase space is done.Moreover as we have noted before, the back-toback kinematics can be obtained from general kinematics not via any integration but by taking specific values of the parameters in the general kinematics.This makes the back-to-back kinematics a special case of the general kinematics.
For realistic experimental observation of back-to-back kinematics: Our treatment of back-to-back kinematics in Ref. [6] is purely mathematical, i.e. we have used the following exact conditions, Out of these, the first and second conditions are primary ones, while the rest arise as a consequence of the first two criteria.However, none of these quantities are physically observable.
The first two conditions also imply that E + = E − = E µ ≡ m B /2 − E ν which is experimentally observable.Any error, on the muon energy measurement would imply that the muon energy distribution as shown in Fig. 5(c) of Ref. [6] would involve energy bins with bin width corresponding to the experimental error within which equality of E + and E − is satisfied.
There would also possibly be extremely slight deviation from the 180 • angle between the two final muons, in an experimental back-to-back realization.This can lead to small deviation of Θ from π, say Θ = π ± ∆Θ.This implies cos Θ = − cos ∆Θ ≈ −1 + ∆Θ 2 /2, so that the error in cos Θ measurement is ∆Θ 2 /2 and it gets multiplied to the back-to-back muon energy distribution of Eq. (50) of Ref. [6] to give the amount of smearing one can expect from the measurement.Thus, a real experimental realization of the back-to-back kinematics would lead to a muon energy distribution curve similar to Fig. 5(c) of Ref. [6] but it would be a histogram plot with bin size determined by the muon energy resolution and there will be some vertical smearing arising from the slight deviation from 180 • angle requirement.However, the difference between Dirac and Majorana neutrinos should not get washed away as a result of such a tiny smearing from experimental measurements.

B Can a massless neutrino with the SM interactions be a Majorana neutrino ?
To properly address this question we need to make a small detour and start from the beginning, the Dirac equation itself, where ψ(x) is the 4-component complex Dirac spinor field that describes a spin-1 2 fermion of mass m, and γ µ (µ = 0, 1, 2, 3) denote the set of four complex 4 × 4 matrices which satisfy the anti-commutation relation {γ µ , γ ν } ≡ γ µ γ ν + γ ν γ µ = 2 g µν , and also ensure the hermiticity of the corresponding Dirac Hamiltonian via γ 0 (γ µ ) † γ 0 = γ µ .The important question to ask here is whether one can have a real solution of the Dirac equation.when viewed from any other basis for the gamma matrices, yields that the Majorana spinor field be identical to its charge conjugate spinor field (more accurately it is the Lorentz-covariant conjugate).This implies that a Majorana fermion is one which is indistinguishable from its antiparticle state.From Eq. (34) it seems clear that one could have a massless Majorana fermion as well.However, in context of the fermion being neutrino, which gets produced only by the weak interaction in the SM, the answer is slightly more involved.Before we address the issue with massless Majorana neutrino in the context of the SM weak interactions we need to take another detour.Using the four gamma matrices, one can always define a fifth gamma matrix, γ 5 = i γ 0 γ 1 γ 2 γ 3 , called the chirality matrix and it commutes with other gamma matrices, it is its own Hermitian adjoint and its own inverse.We note that the matrix γ 5 is also fully imaginary in the Majorana basis.The usefulness of γ 5 is that it allows us to split the complex 4-component spinor field ψ into two distinct parts, where ψ L , ψ R are two distinct eigenfunctions of γ 5 , The two parts of ψ, namely ψ L and ψ R are called the left-chiral and right-chiral spinor fields.In terms of these chiral parts, the Dirac equation gets split into two equations, The space-time evolution of either of the chiral spinor fields is dependent on the mass of the fermion as well as the other chiral spinor field.Thus, when the fermion has non-zero mass (m 0) the chirality is not conserved.However, when m = 0, the chirality is not only conserved, but it also has the same physical meaning as helicity (which is the projection of the fermion spin along its direction of flight).Because a massless particle always travels with speed of light, it has the same chirality or helicity in all frames of reference.Therefore, massless one-half spin fermions of definite chirality are distinct particles.Upon charge conjugation (or Lorentz-covariant conjugation) the chirality of the particle gets reversed, implying that the chirality of the antiparticle is opposite to that of its particle.Thus, for a massless chiral spin-1 2 fermion chirality (or helicity) distinguish between particle and antiparticle which is against the requirement one has for it to qualify as a Majorana fermion.In the SM, the weak interaction always produces left-chiral neutrinos and right-chiral antineutrinos.If neutrino is taken to be a massless fermion (m ν = 0), then its chirality would be conserved and remain Lorentz invariant.Thus one can distinguish a massless SM neutrino from the corresponding antineutrino by its chirality.Another way to realize the impossibility of having a massless chiral fermion with Majorana nature is by asking the mathematical question whether a chiral spinor field can ever be real in the Majorana basis.Since γ 5 in Majorana basis is purely imaginary, its eigenfunctions ψ R/L with eigenvalues ±1 in Eq. (37) can never be real.This is a mathematical impossibility.Therefore, one can never have massless chiral Majorana fermions or neutrinos. 12

Fig. 1 .
Fig.1.Distinct pathways that can be explored to probe the Majorana nature of sub-eV neutrinos and overcome the limitations imposed by the "practical Dirac Majorana confusion theorem".

Fig. 2 .
Fig.2.The general kinematics of B 0 µ − µ + ν µ νµ in the rest frame of B, showing the polar angles θ m and θ n , as well as the azimuthal angle ϕ.Here X m and X n denote the muon pair and the neutrino pair.
It turns out that, if one works in Majorana basis which has only imaginary γ matrices, say γ µ = −i γµ where γµ are real 4 × 4 matrices, then one can have a real spinor field ψ(x) which satisfies the equation, γµ ∂ µ − m ψ(x) = 0. (34) Such a basis of fully imaginary gamma matrices is called the Majorana basis and the real solution to Dirac equation is said to describe the Majorana fermion.The reality condition in Majorana basis, ψ(x) = ψ * (x) (35)