Neutrinos: Majorana or Dirac?

Are neutrinos with definite masses Majorana or Dirac particles? This is one of the most fundamental problem of the modern neutrino physics. The solution of this problem could be crucial for understanding of the origin of small neutrino masses. We will review here basic arguments in favor of the Majorana nature of massive neutrinos. The phenomenological theory of $0\nu\beta\beta$-decay is briefly discussed and recent experimental data and sensitivity of future experiments are presented.


Introduction
The origin of small neutrino masses, discovered in neutrino oscillation experiments, is the major problem of the modern neutrino physics.From neutrino oscillation data five neutrino oscillation parameters (two neutrino mass-squared differences and three mixing angles) were inferred with accuracies (5-10) %.The aim of current and future experiments is to improve accuracy of the measurement of these parameters and to answer the following basic questions • What is the character of the neutrino mass spectrum (Normal or Inverted Ordering?) • What is the value of the CP phase δ?
• How many neutrinos with definite masses ν i exist in nature?Is the number of ν i is equal to the number of flavor neutrinos ν l (l = e, µ, τ ) or larger (in this case exist additional sterile neutrinos)?
• What is the nature of neutrinos with definite masses?Are they Majorana or Dirac particles?
The solution of these problems will be extremely important for the understanding of the origin of neutrino masses.We will discuss here the problem of the neutrino nature, which, apparently, is the most fundamental one.
Neutrinos with definite masses are Dirac particles if the total lepton number L is conserved.In this case neutrino ν i and antineutrino νi have the same mass (CP T ) and different lepton numbers (L(ν i ) = −L(ν i ) = 1).Neutrinos with definite masses are Majorana particles if there no conserved lepton number (i.e there is no conserved quantum number which allows to distinguish neutrino or antineutrino).
There is a general believe that neutrinos are Majorana particles.We will start with some general argument in favor of Majorana neutrinos.The famous two-component theory of a massless neutrino was proposed in 1957 by Landau [1], Lee and Yang [2] and Salam [3] and was confirmed in in the classical Goldhaber at al. experiment on the measurement of the neutrino helicity [4].The two-component Weil field ν L (x) is the simplest possibility for massless neutrino: two degrees of freedom (instead of four in the case of four-component Dirac neutrino). 1e know now that neutrinos have small masses.For neutrino with a mass Majorana field is a simplest, most economical possibility: two degrees of freedom (left-handed and right-handed neutrinos).The Standard Model teaches us that the nature chooses the simplest possibilities.It looks very plausible that also in the case of neutrino with mass a simplest Majorana possibility is realized.
Neutrino masses are many orders of magnitude smaller then masses of leptons and quarks.It is very unlikely that neutrino masses are of the same Standard Model (SM) Higgs origin as masses of other fundamental fermions.If we assume that Standard Model neutrinos are two-component massless particles in this case neutrino masses are generated by a new, beyond the SM mechanism.The method of the effective Lagrangian allows us to describe effects of a beyond the SM physics in the electroweak region.There exist only one SU L (2) × U Y (1) invariant effective Lagrangian which generates a neutrino mass term [5].We can built such a Lagrangian only if we assume that the total lepton number L is not conserved.After spontaneous symmetry breaking this effective Lagrangian generates Majorana neutrinos with definite masses, three-neutrino mixing and seesaw-type suppression of neu-trino masses with respect to Standard Model masses of lepton and quarks.This is a simplest and most plausible, beyond the SM possibility to generate neutrino masses and mixing. 2here are many models which (after heavy fields are integrated out) lead to the Weinberg effective Lagrangian (or its generalization) and the Majorana mass term (see review [6]).In spite all such models were proposed to explain smallness of neutrino masses, values of masses can not be predicted (many unknown parameters are involved).There are, however, two general features which are common to all models, based on the assumption that SM neutrinos are massless particles and that beyond SM particles, responsible for generation of neutrino masses, are heavy: 1. Neutrinos with definite masses are Majorana particles.
2. The number of neutrinos with definite masses is equal to the number of flavor neutrinos (there is no sterile neutrinos).
The problem of sterile neutrinos, which started about 25 years with first data of the LSND experiment [7], is still open.We will briefly discuss it later.The study of neutrinoless double β-decay of some even-even nuclei (0νββ-decay) is the most sensitive way which could allow us to discover the non conservation of the total lepton number L and to reveal the Majorana nature of neutrinos with definite masses.We will briefly discuss the phenomenological theory of this process and present the latest data.
In conclusion we would like to stress that in spite of strong arguments in favor of Majorana neutrinos the possibility of Dirac neutrinos (of the Standard Model or beyond the SM origin) is not excluded.The observation of the 0νββ-decay will allow to exclude this, apparently, artificial possibility.

On the Higgs Mechanism of the Generation of Fermion Masses
From the discovery of neutrino oscillations, measurement of neutrino masses in tritium experiments and cosmological data follows that neutrino masses are many orders of magnitudes smaller than masses of leptons and quarks.From this basic experimental fact we can conclude that it is very plausible that neutrino masses and masses of other fundamental fermions are of different origin.We will discuss a possible origin of neutrino masses later.In this introductory section we will consider the standard Higgs mechanism of the generation of fermion masses.
The standard Higgs mechanism of the generation of fermion masses is based on the assumption that in the total Lagrangian there are SU L (2) × U Y (1) invariant Yukawa interactions.For the charged leptons the Yukawa interaction has the form Here are lepton and Higgs doublets, l ′ R (x) is a right-handed lepton singlet field and Y is a 3 × 3 dimensionless complex matrix.The requirements of the SU L (2) × U Y (1) invariance do not put any constraints on the matrix Y .
Charged lepton masses are generated after the spontaneous symmetry breaking.Let us introduce the hermitian field of the neutral Higgs particles H(x) and choose the Higgs doublet in the form (the unitary gauge) Here is the Higgs vacuum expectation value (vev).With such a choice the local SU L (2)×U Y (1) gauge symmetry will be spontaneously broken.
From ( 1) and ( 3) for the Yukawa Lagrangian we find the following expression The 3 × 3 complex matrix Y can be diagonalized by the biunitary transformation: where V L and V R are unitary matrices and From ( 4) and ( 5) we find that the Yukawa Lagrangian takes the form Here l(x) = l L (x) + l R (x) is the field of the leptons l ± (l = e, µ, τ ) with the mass The fields l L,R (x) are connected with primed fields l ′ L,R (x) by the unitary transformations Thus, Yukawa interaction (1) after spontaneous symmetry breaking generates the standard mass term of the charged leptons and the Lagrangian of the interaction of charged leptons and the Higgs field.The Yukawa constants of this interaction are determined by the lepton mass and are given by The SM masses of up and down quarks are generated in the same way as charged lepton's masses.For masses of quarks we find The Yukawa constants of the quark-Higgs interaction are determined by quark masses and are given by Let us notice that Yukawa constants for different channels were determined from the LHC study of the decay of the Higgs boson into fermion-antifermion pairs (see [8]).The obtained data are in a good agreement with the SM predicted values (10) and (12).From our point of view this is an important argument in favor of the Higgs mechanism of quark and lepton masses generation.
Formally neutrino masses can also be generated by the standard Higgs mechanism.In fact, let us assume that into the total Lagrangian enters the where φ = iτ 2 φ * is the conjugated Higgs doublet and right-handed fields ν ′ lR are singlets.Thus, we need to assume that not only left-handed flavor neutrino fields ν ′ lL but also right-handed fields ν ′ lR are SM fields.After the spontaneous symmetry breaking we find The proportional to v term of this expression is the neutrino mass term.
After the standard diagonalization of 3 × 3 matrix Y ν for the neutrino mass term we obtain the following expression where ν i (x) = ν iL (x) + ν iR (x) is the field of neutrino with the mass It is easy to check that the Lagrangian of Standard Model with the Yukawa interaction ( 14) is invariant under the global transformation where Λ is an arbitrary constant.Thus, in the Standard Model with massive neutrinos the total lepton number L is conserved and ν i (x) is the Dirac field of neutrinos (L = 1) and antineutrinos (L = −1).
In spite Dirac neutrino masses can be introduced in the Standard Model, this possibility looks extremely implausible.Main reason is connected with the smallness of neutrino masses and neutrino Yukawa couplings y ν i = m i v .Absolute values of neutrino masses at present are unknown.However, from existing neutrino oscillation and cosmological data for the heaviest neutrino mass the following conservative bounds can be found where ∆m 2 A ≃ 2.5 • 10 −3 eV 2 is the atmospheric neutrino mass-squared difference.From (18) for the Yukawa coupling y ν 3 we have Yukawa couplings of other particles of the third generation (t, b quarks and τ -lepton) are equal, respectfully Thus, Yukawa coupling of the heaviest neutrino is more than ten orders of magnitude smaller than Yukawa couplings of other particles of the third family.
Notice also that in the SM Lagrangian, which do not include Yukawa interactions, enter left-handed and right-handed fields of all charged particles.Generation of their masses do not require additional degrees of freedom.The SU L (2) × U Y (1) invariant SM can be built only with left-handed neutrino fields.The SM generation of neutrino masses requires right-handed neutrino fields, additional degrees of freedom.Right-handed neutrino fields are out of line of economy and simplicity of the Standard Model.
Thus, it is very unnatural and unlikely that neutrino masses and masses of leptons and quarks are of the same Standard Model origin.In the next section we will consider the most plausible (and popular) beyond the Standard Model mechanism of the generation of neutrino masses.

The Weinberg Effective Lagrangian Mechanism of the Neutrino Mass Generation
The Weinberg effective Lagrangian mechanism of the generation of small neutrino masses [5] is, apparently, the most popular beyond the SM mechanism.
Before discussing this mechanism we will make the following remark.
In the framework of the approach, based on a neutrino mass term, neutrino masses and mixing were introduced for the first time by Gribov and Pontecorvo in 1969 [9]. 3 At that time only ν e and ν µ were known and it was established that the lepton charged current had a V − A form In 1967 Pontecorvo's paper [11] for the two lepton flavors all possible neutrino oscillations (between active neutrinos and active and sterile neutrinos) were considered.Gribov and Pontecorvo put the following question: is it possible to introduce neutrino masses and oscillations in the case if we assume that exist only left-handed neutrino fields ν eL and ν µL and there are no righthanded sterile fields?The authors of the paper [9] showed that it is possible to introduce neutrino masses in this case if the total lepton number L is not conserved. 4eutrinos have masses and are mixed if a neutrino mass term enters into the total Lagrangian.The mass term is a sum of Lorenz-invariant products of left-handed and right-handed components.Authors of the paper [9] took into account that the conjugated field is right-handed.Here C is the matrix of the charge conjugation which satisfies the relations Thus, from the left-handed neutrino fields ν lL (l = e, µ, τ ) it was possible to built the following neutrino mass term where M M is a complex, non diagonal 3 × 3 matrix.From requirements of the Fermi-Dirac statistics it follows that M M = (M M ) T .Let us stress that 1.The Majorana mass term L M is the only possible neutrino mass term which can be built from the left-handed neutrino fields ν lL .This also means that the Majorana mass term is the most economical general possibility for neutrino masses and mixing (there are no right-handed neutrino fields in the Lagrangian, the number of neutrino degrees of freedom is minimal).
2. The mass term L M does not conserve the total lepton number L = L e + L µ + L τ .
The mass term L M is called the Majorana mass term.
The symmetrical matrix M can be diagonalized by the following transformation Here U U † = 1 and m ik = m i δ ik , m i > 0. From ( 23) and ( 24) we find where From ( 25) and (26) follows that • ν i is the field of neutrino with the mass m i .
• The field ν i satisfies the Majorana condition and is the Majorana field.
• The flavor neutrino fields ν lL are mixed fields The unitary 3 × 3 mixing matrix U is called the PMNS matrix [12,13].Let us stress that in the phenomenological Gribov-Pontecorvo approach neutrino masses m i are parameters.There are no any explanation of the smallness of these parameters.
Weinberg [5] proposed a beyond the SM mechanism of neutrino mass generation which leads to the Majorana mass term and allows to explain smallness of neutrino masses.This mechanism is based on the effective Lagrangian approach.
If SM particles interact with heavy, beyond SM particles with masses much larger than v ≃ 246 GeV, then in the electroweak region fields of heavy particles can be "integrated out" and this new interaction induce a nonrenormalizable interaction which is described by the effective Lagrangian.The effective Lagrangians are dimension five or more operators, invariant under SU L (2) × U Y (1) transformations and built from the Standard Model fields.
We are interested in the effective Lagrangian which generates a neutrino mass term.We assume that only left-handed neutrino fields ν ′ lL , components of the lepton doublets ψ lL , enter in the SM.
Let us consider where φ = iτ 2 φ * is the conjugated Higgs doublet.After the spontaneous symmetry breaking for the proportional to v term we have From this expression it is obvious that the only possible effective Lagrangian which generate neutrino mass term has the form [5] (31) where X ′ is a symmetrical, non diagonal matrix.
The operator in (31) has a dimension M 5 .Because a Lagrangian has dimension M 4 , the parameter Λ, which has a dimesion M, is introduced in (31).The parameter Λ characterizes a scale of a new physics (at Λ → ∞ effects of a new physics disappear).We could expect that Λ ≫ v.
Let us stress again that the operator in (31) is the only dimension five effective Lagrangian.Other effective Lagrangians have dimension six and higher.In such effective Lagrangians enter coefficients 1  Λ n with n ≥ 2. This means that the investigation of effects of neutrino masses and mixing, ( neutrino oscillations, neutrinoless double β-decay etc) is the most sensitive way to probe a beyond the SM new physics.
Let us return back to the effective Lagrangian (31).After the spontaneous symmetry braking we find the following neutrino mass term We can present L M in the standard form in which flavor neutrino fields ν lL enter.The SM leptonic charged current has the following form Taking into account ( 9) from ( 33) we find Here l is the field of the charged lepton with mass m l and the flavor (current) neutrino field ν lL is given by From ( 32) and ( 35) we find the following Majorana mass term in which the Majorana matrix M M is given by the expression where The symmetrical, dimensionless matrix X can be presented in the form where U U † = 1 and From (36), (37) and (38) we find Here is the field of the Majorana neutrino with the mass and where U is the unitary PNMS mixing matrix.Generated by the standard Higgs mechanism masses of quarks and leptons (and also W ± , Z 0 and Higgs bosons) are proportional to the Higgs vacuum expectation value v.This is obviously connected with the fact that v is the only SM parameter which has dimension of mass.The Weinberg mechanism of the neutrino mass generation is characterized by two parameters with the dimension of mass: v and Λ.We have m i → 0 at Λ → ∞.It is obvious also that at v → 0 neutrino masses disappear.Thus, for neutrino masses, generated by the effective Lagrangian mechanism, we naturally come to the expression (41) from which it follows that generated by this mechanism neutrino masses are suppressed with respect to SM masses of quarks and leptons by the factor which is naturally much smaller than one.
From (41) we can try to estimate the parameter Λ, which characterize the scale of a new physics.In accordance with latest neutrino oscillation and cosmological data let us assume hierarchy of neutrino masses.In this case we have The parameters x i in (41) are unknown.If we assume (by analogy with Yukawa couplings of the particles of the third family) that x 3 1 we obtain the following estimate Λ 10 15 GeV.
The effective Lagrangian (31) does not conserve the total lepton number L.
Notice that the global invariance and conservation of L (and the barion number B) are not proper symmetries of the Quantum Field Theory (constant phases are not dynamical variable etc).In the Standard Model local gauge symmetry and renormalizable Yukawa interactions ensure conservation of L and B [14].We could expect that a beyond the SM theory does not conserve L and B (see recent discussion in [15]).This is an additional general argument in favor of Majorana nature of beyond the SM neutrino masses.
In conclusion, let us stress again that we assumed (and this is our basic assumption) that in the Lagrangian of the minimal, renormalizable Standard Model there is no neutrino mass term.Then in the framework of the nonrenormalizable, beyond the Standard Model effective Lagrangian approach the residual SU L (2) × U(1) R symmetry naturally ensure the smallness of neutrino masses(via the additional factor v Λ in the expression for the neutrino mass ( 41)).

On the origin of the Weinberg Effective Lagrangian
Let us consider the lepton number violating, SU L (2) × U Y (1) invariant interaction [5,14] Here is the field of the Majorana heavy leptons, SU L (2) × U Y (1) singlet, ψ lL and φ are SM lepton and conjugated Higgs doublets and y li are Yukawa couplings.
If masses of the heavy leptons M i are much larger than v the Lagrangian (46) in the second order of the perturbation theory generates the Weinberg effective Lagrangian.In fact, for the S-matrix we find In the electroweak region Q 2 ≪ M 2 i for the heavy lepton's propagator we have 0|T ( From ( 48) and ( 49) we obtain the effective Lagrangian (31) in which From this relation it follows that a scale of a new physics Λ is determined by masses of heavy Majorana leptons N i . 5 Thus, the Weinberg effective Lagrangian and (after spontaneous symmetry breaking) Majorana neutrino mass term can be generated by the exchange of heavy virtual Majorana leptons, SU L (2) × U Y (1) singlets, between lepton-Higgs pairs.This mechanism of the generation of Majorana neutrino masses and mixing is called type-I seesaw mechanism [19,20,21,22,23] The interaction (46) is not, however, the only possible interaction which generate the effective Lagrangian (31) in the tree-approximation.The Weinberg effective Lagrangian can be also generated in the tree-approximation by interaction of heavy triplet scalar boson field with pair of lepton doublet and pair of Higgs doublet fields.This scenario is called the type-II seesaw mechanism.
Finally, the Weinberg effective Lagrangian can be generated in the tree approximation by the exchange of heavy virtual Majorana triplet leptons between lepton-Higgs pairs.This mechanism is called the type-III seesaw mechanism.

General Remarks on the Neutrino Mass Generation
The interaction (46) and other seesaw interactions are minimal (in a sense of new degrees of freedom) beyond the SM, lepton number violating possibilities of generation of the effective Weinberg Lagrangian and small Majorana neutrino masses.However, the scale Λ of a new physics is naturally very large , unreachable in laboratory experiments (Λ 10 15 GeV).This circumstance inspired creations of many non minimal radiative neutrino mass 5 Heavy Majorana leptons could explain the barion asymmetry of the Universe.In fact, decays of N i 's, produced in the early Universe, into Higgs-lepton pairs could create a lepton asymmetry, if the interaction (46) violates CP .The lepton asymmetry via QCD sphaleron processes can generate the barion asymmetry.This mechanism of the generation of the barion asymmetry of the Universe is called leptogenesis (see reviews [16,17,18]).models with new physics at much lower scales.Detailed discussion of such models and their classification can be found in the review [6] where references to numerous original papers can be found.
A wide class of radiative models is based on the assumption that neutrinos in the Standard Model are massless, left-handed particles and that there are more new heavy fields and more interactions, than in the tree-level seesaw cases.The Weinberg effective Lagrangian is generated by interactions of SM fields and beyond the SM heavy fields via one or more loops.Smaller Λ then in the the seesaw cases are usually needed to explain values of neutrino masses. 6n many radiative neutrino mass models dimension five Weinberg effective Lagrangian (31) is generated.Some models lead to a dimension 5 + 2n effective Lagrangian Here n = 1, 2, ....For the Majorana neutrino masses we have in this case From this expression it follows that a scale of a new physics Λ, much smaller then in the classical Weinberg case, is required in order to ensure smallness of neutrino masses.
Let us now try to extract some general conclusions from a wide class of neutrino mass models, tree-level and radiative, which generate the Weinberg effective Lagrangian or its generalization and are based on the following assumptions • The Standard Model neutrinos ν l are massless left-handed particles (there are no right-handed neutrino fields in the SM).
• Small neutrino masses are generated by a L-violating, beyond a SM interactions and a scale of a new physics is much larger than the electroweak scale v.
In the electroweak region effects of beyond the SM particles are integrating out and the effective Weinberg Lagrangian (or its generalizations) and, correspondingly, the Majorana neutrino mass term is generated.Basic assumptions, listed above, ensure smallness of Majorana neutrino masses with respect to the SM masses of leptons and quarks.However, values of neutrino masses depend on unknown parameters and can not be predicted.The common features of all models, in which the Weinberg effective Lagrangian is generated, are the following 1.Neutrinos with definite masses ν i are Majorana particles.
2. The number of neutrinos with definite masses ν i is equal to the number of the lepton flavors (three).
The most sensitive experiments, which allow us to probe the Majorana nature of ν i , are experiments on the search for neutrinoless double β-decay of some even-even nuclei.In the next sections we will briefly consider this process.
If the number of massive neutrinos if equal to three, it will be no sterile neutrinos, neutrinos which do not have the Standard Model weak interaction.As it is well known, indications in favor of sterile neutrinos were obtained in different short baseline neutrino experiments: in the appearance νµ → νe LSND [7] and [24] experiments, in the disappearance νe → νe reactor experiments [25] and in the disappearance ν e → ν e source Gallium experiments [26].
The existing data can be explained by neutrino oscillations if we assume that fields of three flavor neutrinos ν e,µ,τ and one sterile neutrino ν s are mixtures of the fields of four massive neutrinos ν 1,2,3,4 with m 1,2,3 ≪ m 4 (so called 3 + 1 scheme).From analysis of the data follows that heavy mass m 4 is in the range (10 −1 m 4 10) eV.
Several new short baseline reactor, accelerator and source neutrino experiments are going on or in preparations at present.Analysis of the status of the light sterile neutrino was done in [27,28].The latest update of the status of experiments on the search for sterile neutrinos can be found in talks presented at the NEUTRINO 2020 conference (see http://nu2020.fnal.gov).From existing data it is not possible to make a definite conclusion on the existence of sterile neutrinos.I will limit myself by the following remarks: 1.In new reactor experiments DANSS [29] (with reactor-detector distances 10.7, 11.7 and 12.7 m) and PROSPECT [30] (with reactor-detector distances in the range 6.7-9.2 m) no indications in favor of short baseline neutrino oscillations were found.Best-fit point of previous reactor experiments is excluded in both experiments (at 5σ in the DANSS experiment).
2. Combined analysis of the data, obtained in the reactor Daya Bay and Bugey-3 experiments and accelerator MINOS+ experiment, allows to exclude at 90 % CL LSND and MiniBooNE allowed regions for ∆m 2 14 < 5 eV 2 [31].
3. In the framework of 3 + 1 oscillation scheme with one mass-squared difference for all oscillation channels there exist a relation between oscillation amplitudes which allows to predict from LSND and old reactor data the oscillation amplitude of ν µ → ν µ transition [32,33].This prediction is in a strong tension with negative results of the search for short-baseline ν µ oscillations in MINOS+ and IceCube experiments (see [34,27]).Notice that this problem does not disappear if the data are analyzed in the framework of schemes with two or more sterile neutrinos.
4. Sensitive searches for light sterile neutrinos will be performed in the nearest future.In the SBN experiment at Fermilab [35] three liquid argon TPC detectors at the distances 110 m, 470 m and 600 m from the target will be used.In the same experiment the search for oscillations in appearance ν µ → ν e and disappearance ν µ → ν µ channels will be performed.It is planned that the full 99% LSND allowed region will be covered at more than 5 σ.In the JSNS 2 experiment at J-PARC (Japan) [36] the direct test of the LSND result will be performed.Much better than in the LSND experiment beam and better detector (Gdloaded liquid scintillator) will be used.It is expected that a background in the JSNS 2 experiment will be much lower than in the LSND experiment.The experiment with one detector would allow to exclude the larger part of the LSND-allowed region.In future experiment with two detectors all LSND-allowed region will be covered.There are all reasons to believe that the sterile neutrino anomalies will be resolved in the nearest years.

On the Phenomenological Theory of the 0νββ-decay
Experiments on the search for the lepton number violating 0νββ-decay of some even-even nuclei ensure an unique probe of the Majorana nature of neutrinos with definite masses.These experiments have many advantages with respect to other possible L-violating experiments: large targets (in future experiments about 1 ton or more), law backgrounds, high energy resolutions etc.However, 0νββexperiments are extremely difficult and challenging.This is connected with the fact that the expected probabilities of the neutrinoless double β-decay of different nuclei (in the case if ν i are Majorana particles) are extremely small.The main reasons are the following: • The 0νββ-decay is the second order of the perturbation theory in the Fermi constant G F ≃ 1.166 • 10 −5 1 GeV 2 process.• Because of the V − A nature of the weak interaction the matrix element of the 0νββ-decay is proportional to the effective Majorana mass m ββ = k U 2 ek m k .Smallness of the neutrino masses is a reason for the additional severe suppression of the probability of the decay.
• There are two possibilities for neutrino mass spectra: Normal Ordering or Inverted Ordering (see later).Existing neutrino oscillation data favor Normal Ordering.In the case of the Normal Ordering the effective Majorana mass is much smaller than in the case of the Inverted Ordering.
The standard theory of the neutrinoless double β-decay is based on the following assumptions (see reviews [37,38,39,40]) 1.The CC interaction is the SM interaction Here Here j CCquark α (x) is the quark charged current, W α (x) is the field of the charged W ± vector bosons and g is the constant of the electroweak interaction.
2. The flavor neutrino field ν lL (x) is given by the mixing relation where is the field of the Majorana neutrino with mass m k and U is the unitary 3 × 3 PMNS mixing matrix. 7xisting weak interaction data are in perfect agreement with (55).All atmospheric, accelerator, solar and long baseline reactor neutrino oscillation data are in agreement with the three neutrino mixing (57).
The effective Hamiltonian of the β-decay is given by the expression where j α (x) is ∆S = 0 hadronic charged current and the field of electron neutrinos ν eL (x) is given by (57).
In the second order of the perturbation theory in G F the matrix element of the 0νββ-decay is given by the following expression Here p 1 and p 2 are electron momenta, J α (x) is the weak charged current in the Heisenberg representation, N i (N f ) are states of the initial (final) nuclei with 4-momenta From the Majorana condition (58) follows that Using this relation, for the neutrino propagator we have Further, taking into account that and from ( 60) and (62) for the matrix element of 0νββ-decay we obtain the following expression In this expression we can perform integration over x 0 1 , x 0 2 and q 0 .The matrix element of the 0νββ-decay takes the form Here q 0 k = q 2 + m 2 k and |N n is the vector of the state of the intermediate nucleus with 4-momentum P n = (E n , p n ).In (66) the sum over the total system of the states |N n is assumed.Notice that we used the relation The equation ( 66) is an exact expression for the matrix element of the 0νββdecay in the second order of the perturbation theory.The following approximations are standard ones 1.Small neutrino masses m k in the expression for the neutrino energy can be neglected (q 0 k ≃ | q| = q).In fact, an average neutrino momentum is given by q ≃ 1 r , where r ≃ 10 −13 cm is the average distance between nucleons in a nucleus.We have q ≃ 100 MeV ≫ m k .
2. Long-wave approximation e −i p i x i ≃ 1 (i = 1, 2).In fact, we have In the laboratory frame for the energy denominators in (66) we have where M i (M f ) is the mass if initial (final) nucleus.In (68) we neglect nuclear recoil and take into account that ( 2 ) ≪ q.For the matrix element of the 0νββ-decay we find the following expression where is the effective Majorana mass.Thus, due to the smallness of neutrino masses the matrix element of the 0νββ-decay is a product the effective Majorana mass m ββ (which depends on absolute values of neutrino masses with unknown lightest neutrino mass and on (known) θ 13 and θ 12 and on (unknown) Majorana phases) and nuclear part (which include the propagator of virtual neutrino).As a result of this factorization, the total decay rate of the 0νββ-decay has the following general form Here M 0ν is the nuclear matrix element (NME) and G 0ν (Q, Z) is the known phase-space factor which includes effects of the Column interaction of electrons and daughter nucleus.Calculation of the nuclear matrix elements of the 0νββ-decay is a complicated many-body nuclear problem.Different approximate methods are used in such calculations.At the moment the results of different calculations of NME differ by 2-3 times.Discussion of these calculations is out of the scope of this review (see reviews [38,42,43]).

Effective Majorana Mass
In the case of the three neutrino mixing there are two independent masssquared differences.From analysis of experimental data it follows that one mass-squared difference is about 30 times smaller than the other one.The small (large) mass-squared difference is usually called solar (atmospheric) and is denoted by ∆m 2 S (∆m 2 A ). Neutrino masses are labeled in such a way that the solar mass-squared difference is given by8 ∆m 2 S ≡ ∆m 2 12 . (72) From the MSW effect [44,45], observed in solar neutrino experiments, follows that 9 For the mass of the third neutrino m 3 there are two possibilities: 1. Normal ordering (NO) Atmospheric mass-squared difference ∆m 2 A can be determined as follows 10 In NO case for the neutrino masses m 2 and m 3 we obviously have In the case of IO we find If ν i are Majorana particles, the mixing matrix is characterized by three angles and three phases and has the form 9 In fact, the MSW resonance condition has a form ∆m 2 12 cos 2θ 12 = 2 √ 2G F n e E > 0 (n e is the electron number density).From this condition follows that ∆m 2  12 > 0 under the standard assumption that the principal values of the mixing angle θ 12 are in the range 0 ≤ θ 12 ≤ π 2 . 10Notice that other definitions of ∆m 2  A are also used in the literature (see [46]). where Values of ∆m Notice that from data of all the latest neutrino oscillation experiments follows that the preferable neutrino mass spectrum is the NO spectrum (see talks at Neutrino 2020, http://nu2020.fnal.gov).The effective Majorana mass strongly depends on the character of the neutrino mass spectrum.For illustration we will consider the following viable neutrino mass spectra Hierarchy of the neutrino masses In this case we have Thus, in the case of the neutrino mass hierarchy neutrino masses m 2 and m 3 are determined by the solar and atmospheric mass-squared differences and lightest mass m 1 is very small.Neglecting its contribution and using the standard parametrization of the PMNS mixing matrix for the effective Majorana mass we have where α is a phase difference.The first term in ( 84) is small because of the smallness of ∆m Inverted hierarchy of neutrino masses In this case from (78) we have The lightest mass m 3 in the expression for the effective Majorana mass is multiplied by the small parameter sin 2 θ 13 .Neglecting the contribution of this term, we find where α = α 2 − α 1 is the Majorana phase difference, the only unknown parameter in the expression (88).From (88) we find the following upper and lower bounds for the effective Majorana mass which are realized if there is CP invariance in the lepton sector.
From (89) and Table I we find the following 3 σ range for the effective Majorana mass in the case of the inverted hierarchy of neutrino masses As we will see in the next section, future experiments will be sensitive to the inverted hierarchy region.
8 Experiments on the Search for 0νββ-decay Up to now neutrinoless double β-decay was not observed.The results of some recent experiments on the search for 0νββ-decay are presented in the Table II (see [48,49] ).In the third column of the Table II the 90% CL lower bounds for the half-live of the decay of different elements are given.In the fourth column upper bound ranges for the effective Majorana mass |m ββ | are presented.It was assumed that g A = 1.27 (g A is the axial constant) and nuclear matrix elements, calculated in different models, were taken from [42].
Table II Lower limits of half-lives and upper limits for the effective Majorana mass, obtained in recent experiments on the search for the 0νββ-decay experiment nucleus T 1/2 (10 25 yr) |m ββ | (eV) GERDA [49] 76 Ge 9 0.10-0.23KamLAND-Zen [50] 136 Xe 10.7 0.08-0.24EXO-200 [51] 136 Xe 1.8 0.09-0.29 CUORE [52] 130 Te 3.2 0.08-0.35CUPID-0 [53] 82 Se 0.24 0.39-0.81 Many new experiments on the search for the 0νββ-decay are in preparation (see [48,54]).In these experiments the inverted hierarchy region and, possibly, part of the normal hierarchy region will be probed.We will mention only a few of them.In the KamLAND-Zen experiment ( 136 Xe ) after 5 years of running the sensitivity T 1/2 > 2 • 10 27 yr will be reached.In the SNO+ ( 130 Future experiments on the search for the 0νββ-decay are planned to solve the most fundamental problem of the modern neutrino physics: Are neutrinos with definite masses Majorana or Dirac particles?

Conclusion
Neutrinos are the only fundamental fermions which can be Dirac or Majorana particles.If neutrino masses and mixing are generated by the Standard Higgs Mechanism in this case into the SM Lagrangian enter right-handed neutrino fields ν lR , the total lepton number L is conserved and neutrino with definite masses are Dirac particles.We know, however, that neutrino masses are many orders of magnitude smaller then masses of other fundamental fermions.It is very unlikely that neutrinos masses are generated by the same Higgs mechanism as masses of leptons and quarks.It is very plausible that Standard Model neutrinos are massless particles and neutrino masses are generated by a beyond the SM mechanism.
The most economical and general possibility to introduce neutrino masses and mixing is based on the assumption that there exist only left-handed flavor neutrino fields ν lL (l = e, µ, τ ).Neutrino masses can be introduced in this case only if the total lepton number L is not conserved and neutrinos with definite masses ν i (i = 1, 2, 3) are Majorana particles (Majorana mass term).
The Weinberg effective Lagrangian (and its generalizations), which are based on numerous beyond the SM models, lead to the Majorana neutrino mass term, the three-neutrino mixing U li ν iL , ν i = ν c i and neutrino masses m i , which are naturally suppressed with respect to SM masses of leptons and quarks, by a ratio of the electroweak scale v and a scale of a new, beyond the SM physics Λ.
Non conservation of the lepton number, which is the basics of such explanation of the smallness of neutrino masses, is an attractive and natural feature of a beyond the SM physics (see [15]).
In spite of very strong arguments in favor of the Majorana nature of neutrinos, a possibility of the Dirac neutrinos is not excluded experimentally.Notice that even a beyond the SM seesaw mechanism for the Dirac neutrinos, which could explain the smallness of neutrino masses, was proposed (see [6]).Observation of the neutrinoless double β-decay and a proof that neutrinos with definite masses are Majorana particle would be unique and crucial test of our basic understanding of the phenomenon of small neutrino masses.

)
If neutrinos with definite masses ν i are Dirac particles the PNMS mixing matrix is characterized by three mixing angles θ 12 , θ 23 , θ 13 and one phase δ.In the standard parametrization the mixing matrix has the form U D =   c 13 c 12 c 13 s 12 s 13 e −iδ −c 23 s 12 − s 23 c 12 s 13 e iδ c 23 c 12 − s 23 s 12 s 13 e iδ c 13 s 23 s 23 s 12 − c 23 c 12 s 13 e iδ −s 23 c 12 − c 23 s 12 s 13 e iδ c 13 c 23   .(79) Here c ik = cos θ ik , s ik = sin θ ik .If in the lepton sector CP is conserved U D * = U D and δ = 0.
is the radius of a nucleus.Taking into account that p i 1 MeV we find | p i • x i | ≪ 1. 3. Closure approximation.The energy of neutrino in an intermediate state q ≃ 100 MeV is much larger than the excitation energy (E n − E i ).This means that energies of intermediate states E n in the denominators of Eq.(66) can be replaced by an average energy Ē.In this approximation we can sum over intermediate states |N n in (66).
2S and ∆m2A , three neutrino mixing angles and CP phase, obtained from the global fit of the neutrino oscillation data, are presented in the TableI.
2 S .The contribution to |m ββ | of the "large" atmospheric mass-squared difference ∆m 2 A is suppressed by the smallness of sin 2 θ 13 .As a result, absolute values of the first and second terms in (84) are of the same order of magnitude.From (84) and Table I we find the 3σ upper bound |m ββ | ≤ 4.29 • 10 −3 eV, (85) which is much smaller than upper bounds of |m ββ | reached in modern experiments on the search for the 0νββ-decay (see later).