Neutrinos: Majorana or Dirac?

S. M. Bilenky

Joint Institute for Nuclear Research, Dubna, R-141980, Russia

arXiv:2008.02110v1 [hep-ph] 3 Aug 2020

Abstract
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 understand-
ing 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νββ-decay is briefly discussed and recent
experimental data and sensitivity of future experiments are presented.

1 Introduction
The origin of small neutrino masses, discovered in neutrino oscillation ex-
periments, 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 accu-
racies (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 In-

verted 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 Majo-
rana or Dirac particles?

1
The solution of these problems will be extremely important for the under-
standing 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 num-
ber 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 num-
ber (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).1
We 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 SUL (2) × UY (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 sym-
metry breaking this effective Lagrangian generates Majorana neutrinos with
definite masses, three-neutrino mixing and seesaw-type suppression of neu-
1
Notice that for massless two-component neutrino and V − A interaction there is no
difference between Dirac and Majorana cases.

2
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.2
There are many models which (after heavy fields are integrated out) lead
to the Weinberg effective Lagrangian (or its generalization) and the Majo-
rana 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.

2 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
2
“Simplicity is a guide to the theory choice” A. Einstein.

3
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 SUL (2) ×
UY (1) invariant Yukawa interactions. For the charged leptons the Yukawa
interaction has the form
√ X ′
LY (x) = − 2 ψ̄l1 L (x)Yl1 l2 l2R (x)φ(x) + h.c.. (1)
l1 ,l2

Here
 ′
  
νlL (x) φ+ (x)
ψlL (x) = ′ (l = e, µ, τ ), φ(x) = (2)
lL (x) φ0 (x)
′
are lepton and Higgs doublets, lR (x) is a right-handed lepton singlet field
and Y is a 3 × 3 dimensionless complex matrix. The requirements of the
SUL (2) × UY (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)
!
0
φ(x) = v+H(x)
√
. (3)
2
√
Here v = ( 2 GF )−1/2 ≃ 246 GeV (GF is the Fermi constant) is the Higgs
vacuum expectation value (vev). With such a choice the local SUL (2)×UY (1)
gauge symmetry will be spontaneously broken.
From (1) and (3) for the Yukawa Lagrangian we find the following ex-
pression
¯l′ (x)Yl l l′ (x)(v + H(x)) + h.c.
X
LY (x) = − 1L 1 2 2R (4)
l1 ,l2

The 3 × 3 complex matrix Y can be diagonalized by the biunitary transfor-

mation:
Y = VL y VR† , (5)
where VL and VR are unitary matrices and
yl′ l = yl δl′ l , yl ≥ 0 (l, l′ = e, µ, τ ). (6)

4
From (4) and (5) we find that the Yukawa Lagrangian takes the form

ml ¯l(x)l(x) − yl ¯l(x)l(x) H(x)

X X
LY (x) = − (7)
l=e,µ,τ l=e,µ,τ

Here l(x) = lL (x) + lR (x) is the field of the leptons l± (l = e, µ, τ ) with the
mass
ml = yl v. (8)
′
The fields lL,R (x) are connected with primed fields lL,R (x) by the unitary
transformations
X † X †
′ ′
lL (x) = (VL )ll1 l1L (x), lR (x) = (VR )ll1 l1R (x) (9)
l1 l1

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
ml
yl = . (10)
v
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

mq = yq v (q = u, d, c, s, t, b). (11)

The Yukawa constants of the quark-Higgs interaction are determined by

quark masses and are given by
mq
yq = . (12)
v
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 pre-
dicted 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 gen-
eration.

5
 Formally neutrino masses can also be generated by the standard Higgs
mechanism. In fact, let us assume that into the total Lagrangian enters the
SUL (2) × UY (1) invariant Yukawa interaction
√ X
LνY (x) = − 2 ψ̄l1 L (x)Yl1νl2 νl′2 R (x)φ̃(x) + h.c. (13)
l1 ,l2

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
X
LνY (x) = − ν̄l′1 L (x)Ylν1 l2 νl′2 R (x)(v + H(x)) + h.c. (14)
l1 ,l2

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
3
X 3
X
Lν (x) = − mi (ν̄iL (x) νiR (x) + h.c.) = − mi ν̄i (x) νi (x), (15)
i=1 i=1

where νi (x) = νiL (x) + νiR (x) is the field of neutrino with the mass
mi = yiν v. (16)
It is easy to check that the Lagrangian of Standard Model with the Yukawa
interaction (14) is invariant under the global transformation
νiL (x) → eiΛ νiL (x), νiR (x) → eiΛ νiR (x), lL (x) → eiΛ lL (x), lR (x) → eiΛ lR (x),
(17)
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 yiν = mvi .
Absolute values of neutrino masses at present are unknown. However,
from existing neutrino oscillation and cosmological data for the heaviest neu-
trino mass the following conservative bounds can be found
1X
q
5 · 10−2 eV ≃ ( ∆m2A ) ≤ m3 ≤ ( mi ) ≃ 3 · 10−1 eV, (18)
3 i

6
where ∆m2A ≃ 2.5 · 10−3 eV2 is the atmospheric neutrino mass-squared dif-
ference. From (18) for the Yukawa coupling y3ν we have

2 · 10−13 ≤ y3ν ≤ 10−12 . (19)

Yukawa couplings of other particles of the third generation (t, b quarks and
τ -lepton) are equal, respectfully

yt ≃ 7 · 10−1 , yb ≃ 2 · 10−2 , yτ ≃ 7 · 10−3 . (20)

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
SUL (2) × UY (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.

3 The Weinberg Effective Lagrangian Mech-

anism of the Neutrino Mass Generation
The Weinberg effective Lagrangian mechanism of the generation of small neu-
trino 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, neu-
trino 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

jαCC = 2(ν̄eL γα eL + ν̄µL γα µL ). (21)

3
The results of this paper were generalized in [10].

7
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 right-
handed 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.4
Neutrinos 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
c T
νlL = C ν̄lL (22)

is right-handed. Here C is the matrix of the charge conjugation which satisfies

the relations
CγαT C −1 = −γα , C T = −C.
Thus, from the left-handed neutrino fields νlL (l = e, µ, τ ) it was possible to
built the following neutrino mass term
1 X
LM = − ν̄l′ L MlM c
′ l νlL + h.c. (23)
2 ′
l ,l

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 LM 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).
4
At that time some authors claimed that if neutrinos are left-handed, their masses had
to be equal to zero. This assertion is based, however, on the assumption that the total
lepton number is conserved.

8
 2. The mass term LM does not conserve the total lepton number L =
Le + Lµ + Lτ .

The mass term LM is called the Majorana mass term.

The symmetrical matrix M can be diagonalized by the following trans-
formation
MM = U m UT . (24)
Here U U † = 1 and mik = mi δik , mi > 0. From (23) and (24) we find
3
M 1X
L =− mi ν̄i νi , (25)
2 i=1

where X X
νi = Uil† νlL + (Uil† νlL )c . (26)
l l

From (25) and (26) follows that

• νi is the field of neutrino with the mass mi .

• The field νi satisfies the Majorana condition

νi = νic (27)

and is the Majorana field.

• The flavor neutrino fields νlL are mixed fields

3
X
νlL = Uli νiL . (28)
i=1

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 mi 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 La-
grangian approach.

9
 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 SUL (2) × UY (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 SUL (2) × UY (1) invariant

(φ̃† ψlL ) (l = e, µ, τ ), (29)

where φ̃ = iτ2 φ∗ is the conjugated Higgs doublet. After the spontaneous

symmetry breaking for the proportional to v term we have
v ′
(φ̃† ψlL )SSB = √ νlL . (30)
2
From this expression it is obvious that the only possible effective Lagrangian
which generate neutrino mass term has the form [5]
1 X † 1 X
Leff
I = − (φ̃ ψl′ L ) Xl′′ l (φ̃† ψlL )c + h.c. = − c
ψ̄l′ L φ̃ Xl′′ l φ̃† ψlL + h.c.
Λ l′ ,l Λ l′ ,l
(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 Λ1n 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

10
symmetry braking we find the following neutrino mass term

M v2 X ′
L =− ν̄l′ L Xl′′ l νlL
′c
+ h.c. (32)
2Λ l′ ,l

We can present LM in the standard form in which flavor neutrino fields νlL
enter. The SM leptonic charged current has the following form
X 1 X
jαCC = 2 ψ̄lL (τ1 + iτ2 )γα ψlL = 2 ′
ν̄lL γα lL′ . (33)
2
l l

Taking into account (9) from (33) we find

X X
jαCC = 2 ′
ν̄lL γα (VL )ll1 l1L = 2 ν̄lL γα lL . (34)
l,l1 l

Here l is the field of the charged lepton with mass ml and the flavor (current)
neutrino field νlL is given by
X †
νlL = (VL )ll1 νl′1 L . (35)
l1

From (32) and (35) we find the following Majorana mass term
1 X
LM = − ν̄l′ L MlM c
′ l νlL + h.c. (36)
2 l′ ,l

in which the Majorana matrix M M is given by the expression

v2
MM = X, (37)
Λ
where X = VL† X ′ (VL† )T = X T .
The symmetrical, dimensionless matrix X can be presented in the form

X = U x UT , (38)

where U U † = 1 and xik = xi δik , xi > 0.

From (36), (37) and (38) we find
3
1X
LM = − mi ν̄i νi . (39)
2 i=1

11
Here
νi (x) = νic (x) = C ν̄iT (x) (40)
is the field of the Majorana neutrino with the mass
v
mi = (vxi ) (41)
Λ
and
X3
νlL = Uli νiL , (42)
i=1
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 mi → 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
v EW scale
= , (43)
Λ scale of a new physics
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 q
m3 ≃ ( ∆m2A ) ≃ 5 · 10−2 eV. (44)
The parameters xi in (41) are unknown. If we assume (by analogy with
Yukawa couplings of the particles of the third family) that x3 . 1 we obtain
the following estimate
Λ . 1015 GeV. (45)
The effective Lagrangian (31) does not conserve the total lepton number L.
Notice that the global invariance and conservation of L (and the barion num-
ber B) are not proper symmetries of the Quantum Field Theory (constant

12
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 con-
serve 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 non-
renormalizable, beyond the Standard Model effective Lagrangian approach
the residual SUL (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)).

4 On the origin of the Weinberg Effective La-

grangian
Let us consider the lepton number violating, SUL (2) × UY (1) invariant inter-
action [5, 14]
√ X
LI = − 2 (ψ̄lL φ̃)yli NiR + h.c. (46)
li

Here
Ni = Nic = C(N̄i )T , i = 1, 2, ...n (47)
is the field of the Majorana heavy leptons, SUL (2) × UY (1) singlet, ψlL and φ̃
are SM lepton and conjugated Higgs doublets and yli are Yukawa couplings.
If masses of the heavy leptons Mi 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

(−i)2
Z X
(2)
S = 2 T( ψ̄l′ L (x1 )φ̃(x1 )yl′i NiR (x1 )
2! ′ l ,l,i,k
T
× NkR (x2 )ylk φ̃T (x2 )ψlL
T
(x2 ))d4 x1 d4 x2 + ... (48)

In the electroweak region Q2 ≪ Mi2 for the heavy lepton’s propagator we

have
T 1 1 + γ5
h0|T (NiR(x1 )NkR (x2 ))|0i ≃ i δ(x1 − x2 ) Cδik . (49)
Mi 2

13
From (48) and (49) we obtain the effective Lagrangian (31) in which
1 ′ X 1
Xl ′ l = yl′ i yli . (50)
Λ i
Mi

From this relation it follows that a scale of a new physics Λ is determined by

masses of heavy Majorana leptons Ni .5
Thus, the Weinberg effective Lagrangian and (after spontaneous symme-
try breaking) Majorana neutrino mass term can be generated by the exchange
of heavy virtual Majorana leptons, SUL (2) × UY (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 Wein-
berg 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.

5 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 possibili-
ties 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 (Λ . 1015 GeV). This cir-
cumstance 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 Ni ’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]).

14
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. 6
In many radiative neutrino mass models dimension five Weinberg effective
Lagrangian (31) is generated. Some models lead to a dimension 5 + 2n
effective Lagrangian
† n
1 c (φ φ)
Leff
I
′n †
= − (ψ̄L φ̃) X (φ̃ ψL ) 2n + h.c. (51)
Λ Λ
Here n = 1, 2, .... For the Majorana neutrino masses we have in this case
v
mi = v ( )2n+1 xni (52)
Λ
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 Wein-
berg 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 elec-
troweak scale v.

In the electroweak region effects of beyond the SM particles are integrat-

ing out and the effective Weinberg Lagrangian (or its generalizations) and,
6 1
One of the reason is that with every loop enters a suppression factor 16π 2 ≃ 0.006.

15
correspondingly, the Majorana neutrino mass term
1 X
LM = − ν̄l′ L MlM c
′ l νlL + h.c. (53)
2
l′ ,l=e,µ,τ

is generated.
Basic assumptions, listed above, ensure smallness of Majorana neutrino
masses with respect to the SM masses of leptons and quarks. However, val-
ues of neutrino masses depend on unknown parameters and can not be pre-
dicted. 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 νµ → νe MiniBooNE [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 m1,2,3 ≪ m4 (so
called 3 + 1 scheme). From analysis of the data follows that heavy mass m4
is in the range (10−1 . m4 . 10) eV.
Several new short baseline reactor, accelerator and source neutrino ex-
periments 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:

16
1. In new reactor experiments DANSS [29] (with reactor-detector dis-
tances 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 base-
line neutrino oscillations were found. Best-fit point of previous reactor
experiments is excluded in both experiments (at 5σ in the DANSS ex-
periment).

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 ∆m214 <
5 eV2 [31].

3. In the framework of 3 + 1 oscillation scheme with one mass-squared

difference for all oscillation channels there exist a relation between os-
cillation amplitudes which allows to predict from LSND and old reactor
data the oscillation amplitude of νµ → νµ transition [32, 33]. This pre-
diction is in a strong tension with negative results of the search for
(−) (−)
short-baseline νµ → νµ oscillations in MINOS+ and IceCube experi-
ments (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 ar-
gon 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 (Gd-
loaded liquid scintillator) will be used. It is expected that a background
in the JSNS 2 experiment will be much lower than in the LSND ex-
periment. 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.

17
6 On the Phenomenological Theory of the
0νββ-decay
Experiments on the search for the lepton number violating 0νββ-decay of
some even-even nuclei
(A, Z) → (A, Z + 2) + e− + e− (54)
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
1
Fermi constant GF ≃ 1.166 · 10−5 GeV 2 process.

• Because of the V − A nature of the weak interaction the matrix ele-

ment ofPthe 0νββ-decay is proportional to the effective Majorana mass
2
mββ = k Uek mk . 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 Order-
ing or Inverted Ordering (see later). Existing neutrino oscillation data
favor Normal Ordering. In the case of the Normal Ordering the effec-
tive 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 fol-
lowing assumptions (see reviews [37, 38, 39, 40])
1. The CC interaction is the SM interaction
g CC
LCC α
I (x) = − √ jα (x)W (x) + h.c.. (55)
2 2
Here X
jαCC (x) = 2 ν̄lL (x)γα lL (x) + jαCCquark (x). (56)
l=e,µ,τ

18
 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
3
X
νlL (x) = Ulk νkL (x), (57)
k=1

where
νk (x) = νkc (x) = C ν̄kT (x), (k = 1, 2, 3) (58)
is the field of the Majorana neutrino with mass mk and U is the unitary
3 × 3 PMNS mixing matrix.7

Existing weak interaction data are in perfect agreement with (55). All atmo-
spheric, 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
GF
HI (x) = √ 2 ēL (x)γα νeL (x) j α (x) + h.c. (59)
2
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 GF the matrix element
of the 0νββ-decay is given by the following expression
2
(−i)2 GF
 Z
(2)
hf |S |ii = 4 √ Np1 Np2 ūL (p1 )eip1 x1 γα h0|T (νeL (x1 ) νeL
T
(x2 )|0i
2! 2
×γβT ūTL (p2 )eip2 x2 hNf |T (J α (x1 )J β (x2 ))|Ni i d4 x1 d4 x2 − (p1 ⇆ p2 ). (60)

Here p1 and p2 are electron momenta, J α (x) is the weak charged current in
the Heisenberg representation, Ni (Nf ) are states of the initial (final) nuclei
7
We assume that a beyond the SM physics contributes to the matrix element of the
0νββ-decay via Majorana neutrino masses. If neutrinos with definite masses are Majorana
particles their masses definitely contribute to the matrix element of the 0νββ-decay. For
discussion of possible additional contribution of effective Lagrangians to the 0νββ-decay
see, for example, [41].

19
 1√
with 4-momenta Pi = (Ei , p~i ) (Pf = (Ef , ~pf )) and Np = is the
(2π)3/2 2p0
standard normalization factor.
From the Majorana condition (58) follows that

νkT (x) = ν̄k (x)C T = −ν̄k (x)C. (61)

Using this relation, for the neutrino propagator we have

T
X
2 1 − γ5 1 − γ5
h0|T (νeL (x1 ) νeL (x2 )|0i = − Uek h0|T (νk (x1 ) ν̄k (x2 )|0i C
k
2 2
−iq(x1 −x2 )
1 − γ5
Z
i X
2 mk e
=− Uek d4 q C. (62)
(2π)4 k
q 2 − m2k 2

Further, taking into account that

1 − γ5 1 − γ5
ūL (p1 )γα γβ C ūTL (p2 ) = −ūL (p2 )γβ γα C ūTL (p1 ) (63)
2 2
and
T (J α (x1 )J β (x2 )) = T (J β (x2 )J α (x1 )) (64)
from (60) and (62) for the matrix element of 0νββ-decay we obtain the
following expression
2 !

GF
Z
i X Z mk e−iq (x1 −x2 )
hf |S (2) |ii = −4 √ Np1 Np2 ūL (p1 )eip1 x1 γα 4
2
Uek 2 2
d4 q
2 (2π) k
q − mk
1 + γ5
×γβ C ūTL (p2 )eip2 x2 hNf |T (J α (x1 )J β (x2 ))|Ni i d4 x1 d4 x2 . (65)
2
In this expression we can perform integration over x01 , x02 and q 0 . The matrix
element of the 0νββ-decay takes the form
 2 Z
(2) GF
hf |S |ii = 2i √ Np1 Np2 ū(p1 )γα γβ (1 + γ5 )C ū (p2 ) d3 x1 d3 x2 e−i~p1 ~x1 −i~p2~x2
T
2
Z i~q (~x1 −~x2 ) X hNf |J α (~x1 )|Nn ihNn |J β (~x2 ))|Ni i
X
2 1 e 3
× Uek mk d q [
k
(2π)3 qk0 n
En + p02 + qk0 − Ei − iǫ
X hNf |J β (~x2 )|Nn ihNn |J α (~x1 ))|Ni i
+ 0 0
] 2πδ(Ef + p01 + p02 − Ei ). (66)
n
E n + p 1 + qk − E i − iǫ

20
 p
Here qk0 = ~q2 + m2k and |Nn i is the vector of the state of the intermediate
nucleus with 4-momentum Pn = (En , p~n ). In (66) the sum over the total
system of the states |Nn i is assumed. Notice that we used the relation
0 0
J α (x) = eiHx J α (~x)e−iHx , J α (~x) = J α (0, ~x). (67)

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 approx-
imations are standard ones

1. Small neutrino masses mk in the expression for the neutrino energy can
be neglected (qk0 ≃ |~q| = q). In fact, an average neutrino momentum is
given by q̄ ≃ 1r , where r ≃ 10−13 cm is the average distance between
nucleons in a nucleus. We have q̄ ≃ 100 MeV ≫ mk .

2. Long-wave approximation e−i~pi~xi ≃ 1 (i = 1, 2). In fact, we have

|~pi · ~xi | ≤ pi R, where R ≃ 1.2 · 10−13 A1/3 cm is the radius of a nucleus.
Taking into account that pi . 1 MeV we find |~pi · ~xi | ≪ 1.

3. Closure approximation. The energy of neutrino in an intermediate state

q ≃ 100 MeV is much larger than the excitation energy (En − Ei ). This
means that energies of intermediate states En in the denominators of
Eq.(66) can be replaced by an average energy Ē. In this approximation
we can sum over intermediate states |Nn i in (66).

In the laboratory frame for the energy denominators in (66) we have

p01 + p02 p0 − p02 Mi + Mf

En + p01,2 + qk0 − Ei ≃ Ē + ( )±( 1 ) + q − Mi ≃ Ē + q − ,
2 2 2
(68)
where Mi (Mf ) is the mass if initial (final) nucleus. In (68) we neglect nuclear
p0 −p0
recoil and take into account that ( 1 2 2 ) ≪ q. For the matrix element of the
0νββ-decay we find the following expression
 2
(2) GF
hf |S |ii ≃ 2i √ Np1 Np2 ū(p1 )γα γβ (1 + γ5 )C ūT (p2 )mββ
2
1 ei~q (~x1 −~x2 )
Z Z
3 3
× d x1 d x2 Mi +Mf
d3 q [hNf |(J α (x~1 )J β (x~2 )
(2π)3 q(Ē + q − 2 )
+J β (x~2 )J α (x~1 ))|Ni i] 2πδ(Ef + p01 + p02 − Ei ), (69)

21
where
3
X
2
mββ = Uek mk (70)
i=1

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
ln 2
Γ0ν = 0ν = |mββ |2 |M 0ν |2 G0ν (Q, Z). (71)
T1/2
Here M 0ν is the nuclear matrix element (NME) and G0ν (Q, Z) is the known
phase-space factor which includes effects of the Column interaction of elec-
trons and daughter nucleus.
Calculation of the nuclear matrix elements of the 0νββ-decay is a compli-
cated 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]).

7 Effective Majorana Mass

In the case of the three neutrino mixing there are two independent mass-
squared 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 ∆m2S (∆m2A ).
Neutrino masses are labeled in such a way that the solar mass-squared
difference is given by8
∆m2S ≡ ∆m212 . (72)
From the MSW effect [44, 45], observed in solar neutrino experiments, follows
8
We will use the following definition: ∆m2ki = m2i − m2k .

22
that9
∆m212 > 0. (73)
For the mass of the third neutrino m3 there are two possibilities:
1. Normal ordering (NO)
m3 > m2 > m1 (74)

2. Inverted ordering (IO)

m2 > m1 > m3 . (75)

Atmospheric mass-squared difference ∆m2A can be determined as follows10

∆m2A = ∆m223 (NO), ∆m2A = |∆m213 | (IO). (76)

In NO case for the neutrino masses m2 and m3 we obviously have

q q
m2 = m21 + ∆m2S , m3 = m21 + ∆m2S + ∆m2A . (77)

In the case of IO we find

q q
m1 = m23 + ∆m2A , m2 = m23 + ∆m2A + ∆m2S (78)

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
 
c13 c12 c13 s12 s13 e−iδ
U D =  −c23 s12 − s23 c12 s13 eiδ c23 c12 − s23 s12 s13 eiδ c13 s23  . (79)
iδ iδ
s23 s12 − c23 c12 s13 e −s23 c12 − c23 s12 s13 e c13 c23
Here cik = cos θik , sik = sin θik . If in the lepton sector CP is conserved
U D∗ = U D and δ = 0.
If νi are Majorana particles, the mixing matrix is characterized by three
angles and three phases and has the form

U M = U D S M (α), (80)
9
√
In fact, the MSW resonance condition has a form ∆m212 cos 2θ12 = 2 2GF ne E > 0
(ne is the electron number density). From this condition follows that ∆m212 > 0 under
the standard assumption that the principal values of the mixing angle θ12 are in the range
0 ≤ θ12 ≤ π2 .
10
Notice that other definitions of ∆m2A are also used in the literature (see [46]).

23
where
eiα1 0 0
 

S M (α) =  0 eiα2 0  . (81)

0 0 1

Values of ∆m2S and ∆m2A , three neutrino mixing angles and CP phase, ob-
tained from the global fit of the neutrino oscillation data, are presented in
the Table I.
Table I

Values of neutrino oscillation parameters, obtained from the global fit of

the data [47]

Parameter Normal Ordering Inverted Ordering

sin2 θ12 0.310+0.013
−0.012 0.310+0.013
−0.012
sin2 θ23 0.582+0.015
−0.019 0.582+0.015
−0.018
sin2 θ13 0.02240+0.00065
−0.00066 0.02263+0.00065
−0.00066
δ (in ◦ ) (217+40−28 ) (280+25
−28 )
∆m2S (7.39−0.20 ) · 10−5 eV2
+0.21
(7.39−0.20 ) · 10−5 eV2
+0.21

∆m2A (2.525+0.033
−0.031 ) · 10
−3
eV2 (2.512+0.034
−0.031 ) · 10
−3
eV2

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
m1 ≪ m2 ≪ m3 . (82)
In this case we have
q q q
m2 ≃ ∆m2S , m3 ≃ ∆m2A m1 ≪ ∆m2S . (83)

Thus, in the case of the neutrino mass hierarchy neutrino masses m2 and
m3 are determined by the solar and atmospheric mass-squared differences
and lightest mass m1 is very small. Neglecting its contribution and using

24
the standard parametrization of the PMNS mixing matrix for the effective
Majorana mass we have
 q q 
 2 2 2i α 2

2 2
|mββ | ≃  cos θ13 sin θ12 ∆mS + e sin θ13 ∆mA  , (84)

where α is a phase difference.

The first term in (84) is small because of the smallness of ∆m2S . The con-
tribution to |mββ | of the “large” atmospheric mass-squared difference ∆m2A
is suppressed by the smallness of sin2 θ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 exper-
iments on the search for the 0νββ-decay (see later).

Inverted hierarchy of neutrino masses

m3 ≪ m1 < m2 . (86)
In this case from (78) we have

∆m2S
q q q q
2 2
m1 ≃ ∆mA , m2 ≃ ∆mA (1+ ) ≃ ∆mA , m3 ≪ ∆m2A . (87)
2
2 ∆m2A
The lightest mass m3 in the expression for the effective Majorana mass is
multiplied by the small parameter sin2 θ13 . Neglecting the contribution of
this term, we find
q
1
|mββ | ≃ ∆m2A cos2 θ13 (1 − sin2 2 θ12 sin2 α) 2 , (88)

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
q q
2 2
cos θ13 cos 2 θ12 ∆mA ≤ |mββ | ≤ cos θ13 ∆m2A ,
2
(89)

which are realized if there is CP invariance in the lepton sector.

25
 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

1.44 · 10−2 eV ≤ |mββ | ≤ 5.01 · 10−2 eV (90)

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 gA = 1.27 (gA 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 T1/2 (1025 yr) |mββ | (eV)

76
GERDA [49] Ge 9 0.10-0.23
136
KamLAND-Zen [50] Xe 10.7 0.08-0.24
136
EXO-200 [51] Xe 1.8 0.09-0.29
130
CUORE [52] Te 3.2 0.08-0.35
82
CUPID-0 [53] 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, pos-
sibly, 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 T1/2 > 2 · 1027 yr will be reached. In the SNO+ (130
Te), LEGEND (76 Ge), n-EXO (136 Xe ), CUPID (100 Mo ) and NEXT-HD
(136 Xe ) sensitivities T1/2 > 1 · 1027 yr, T1/2 > 2 · 1028 yr, T1/2 > 5.7 · 1027 yr,
T1/2 > 1.1 · 1027 yr and T1/2 > 3 · 1027 yr will be reached after 10 years of
running.

26
 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?

9 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
3
X
νlL = Uli νiL , νi = νic
i=1

and neutrino masses mi , 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 ex-
planation 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,

27
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.

References
[1] L. D. Landau, On the conservation laws for weak interactions, Nucl.
Phys. 3 (1957) 127.

[2] T. D. Lee and C. N. Yang, Parity Nonconservation and a Two Compo-

nent Theory of the Neutrino, Phys. Rev. 105 (1957) 1671.

[3] A. Salam, On parity conservation and neutrino mass, Nuovo Cim. 5

(1957) 299.

[4] M. Goldhaber, L. Grodzins and A. W. Sunyar, Helicity of neutrinos,

Phys. Rev. 109 (1958) 1015.

[5] S. Weinberg, Baryon and Lepton Nonconserving Processes, Phys. Rev.

Lett. 43 (1979) 1566.

[6] Y. Cai, J. Herrero-Garca, M. A. Schmidt, A. Vicente and R. R. Volkas,

From the trees to the forest: a review of radiative neutrino mass models,
Front. in Phys. 5 (2017) 63, arXiv:1706.08524.

[7] A. Aguilar-Arevalo et al. [LSND Collaboration], Evidence for neu-

trino oscillations from the observation of anti-neutrino(electron) appear-
ance in a anti-neutrino(muon) beam, Phys. Rev. D 64 (2001) 112007,
hep-ex/0104049.

[8] M. Tanabashi et al. [Particle Data Group], Review of Particle Physics,

Phys. Rev. D 98 (2018) no.3, 030001.

[9] V. N. Gribov and B. Pontecorvo, Neutrino astronomy and lepton charge,

Phys. Lett. 28B (1969) 493.

[10] S. M. Bilenky and S. T. Petcov, Massive Neutrinos and Neutrino Oscil-

lations, Rev. Mod. Phys. 59 (1987) 671 Erratum: [Rev. Mod. Phys. 61
(1989) 169] Erratum: [Rev. Mod. Phys. 60 (1988) 575].

28
[11] B. Pontecorvo, Neutrino experiments and the question of leptonic-charge
conservation, Zh. Eksp. Teor. Fiz. 53(1967) 1717 ( Sov. Phys. JETP 26
(1968) 984).

[12] B. Pontecorvo, Inverse beta processes and nonconservation of lepton

charge, Sov. Phys. JETP 7 (1958) 172 [Zh. Eksp. Teor. Fiz. 34 (1957)
247].

[13] Z. Maki, M. Nakagawa and S. Sakata, Remarks on the unified model of

elementary particles, Prog. Theor. Phys. 28 (1962) 870.

[14] S. Weinberg, Varieties of Baryon and Lepton Nonconservation, Phys.

Rev. D 22 (1980) 1694.

[15] E. Witten, Symmetry and Emergence, Nature Phys. 14 (2018) 116.

[16] W. Buchmuller, P. Di Bari and M. Plumacher, Some aspects of thermal

leptogenesis, New J. Phys. 6 (2004) 105.

[17] S. Davidson, E. Nardi and Y. Nir, Leptogenesis, Phys. Rept. 466 (2008)
105.

[18] P. Di Bari, An introduction to leptogenesis and neutrino properties, Con-

temp. Phys. 53 (2012) no.4, 315.

[19] P. Minkowski, µ → eγ at a rate of one out of 1-billion muon decays?,

Phys. Lett. B 67 (1977) 421.

[20] M. Gell-Mann, P. Ramond and R. Slansky,Complex spinors and unified

theories in Supergravity, p. 315, edited by F. van Nieuwenhuizen and D.
Freedman, North Holland, Amsterdam, (1979).

[21] T. Yanagida, Horizontal gauge symmetry and masses of neutrinos, Conf.

Proc. C 7902131 (1979) 95.

[22] S.L. Glashow, NATO Adv.Study Inst. Ser. B Phys. 59 (1979) 687.

[23] R.N. Mohapatra and G. Senjanović, Neutrino masses and mixings in

gauge models with spontaneous parity violation, Phys. Rev. D23 (1981)
165.

29
[24] A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration], Signifi-
cant Excess of ElectronLike Events in the MiniBooNE Short-Baseline
Neutrino Experiment, Phys. Rev. Lett. 121 (2018) no.22, 221801,
arXiv:1805.12028.

[25] G. Mention, M. Fechner, T. Lasserre, T. A. Mueller, D. Lhuillier, M. Cri-

bier and A. Letourneau, The Reactor Antineutrino Anomaly, Phys. Rev.
D 83 (2011) 073006, arXiv:1101.2755.

[26] C. Giunti and M. Laveder, Statistical Significance of the Gallium

Anomaly, Phys. Rev. C 83 (2011) 065504, arXiv:1006.3244.

[27] M. Dentler, . Hernndez-Cabezudo, J. Kopp, P. A. N. Machado, M. Mal-

toni, I. Martinez-Soler and T. Schwetz, Updated Global Analysis of Neu-
trino Oscillations in the Presence of eV-Scale Sterile Neutrinos, JHEP
08, 010 (2018), arXiv:1803.10661.

[28] S. Bser, C. Buck, C. Giunti, J. Lesgourgues, L. Ludhova, S. Mertens,

A. Schukraft and M. Wurm, Status of Light Sterile Neutrino Searches,
Prog. Part. Nucl. Phys. 111, 103736 (2020), arXiv:1906.01739.

[29] Y. Shitov,Recent results from the DANSS experiment, talk at the In-
ternational Conference Neutrino 2020, June 28-July 2, 2020, Chicago,
USA, http://nu2020.fnal.gov.

[30] M. Andriamirado et al. [PROSPECT], Improved Short-Baseline Neu-

trino Oscillation Search and Energy Spectrum Measurement with the
PROSPECT Experiment at HFIR, arXiv:2006.11210 [hep-ex].

[31] P. Adamson et al. [MINOS+ and Daya Bay], Improved Constraints on

Sterile Neutrino Mixing from Disappearance Searches in the MINOS,
MINOS+, Daya Bay, and Bugey-3 Experiments, arXiv:2002.00301 [hep-
ex].

[32] S. M. Bilenky, C. Giunti and W. Grimus, Neutrino mass spectrum from

the results of neutrino oscillation experiments, Eur. Phys. J. C 1 (1998)
247, hep-ph/9607372.

[33] N. Okada and O. Yasuda, A Sterile neutrino scenario constrained by

experiments and cosmology, Int. J. Mod. Phys. A 12 (1997) 3669,
hep-ph/9606411.

30
[34] S. Gariazzo, C. Giunti, M. Laveder and Y. F. Li, Updated Global 3+1
Analysis of Short-BaseLine Neutrino Oscillations, JHEP 06, 135 (2017),
arXiv:1703.00860 [hep-ph].

[35] P. A. Machado, O. Palamara and D. W. Schmitz, The Short-Baseline

Neutrino Program at Fermilab, Ann. Rev. Nucl. Part. Sci. 69, 363-387
(2019), arXiv:1903.04608 [hep-ex].

[36] S. Ajimura et al. Technical Design Report (TDR): Searching for a Sterile
Neutrino at J-PARC MLF (E56, JSNS2), arXiv:1705.08629 [physics.ins-
det].

[37] M. Doi, T. Kotani and E. Takasugi, Double Beta Decay and Majorana
Neutrino, Prog. Theor., Phys. Suppl. 83 (1985) 1.

[38] J.D. Vergados, H. Ejiri, F. Simkovic, Neutrinoless double beta decay and
neutrino mass, Int.J.Mod.Phys. E25 (2016) 1630007, arXiv:1612.02924.

[39] S. Dell’Oro, S. Marcocci, M. Viel, F. Vissani, Neutrinoless double beta

decay: 2015 review, Adv.High Energy Phys. 2016 (2016) 2162659,
arXiv:1601.07512.

[40] S.M. Bilenky, C. Giunti, Neutrinoless Double-Beta Decay: a Probe of

Physics Beyond the Standard Model, Int.J.Mod.Phys. A30 (2015) 0001,
arXiv:1411.4791.

[41] F. del Aguila, A. Aparici, S. Bhattacharya, A. Santamaria and J. Wudka,

Effective Lagrangian approach to neutrinoless double beta decay and neu-
trino masses, JHEP 06, 146 (2012), arXiv:1204.5986 [hep-ph].

[42] J. Engel, J. Menendez, Status and Future of Nuclear Matrix Elements for
Neutrinoless Double-Beta Decay: A Review, Rept.Prog.Phys. 80 (2017)
046301, arXiv:1610.06548.

[43] J. Menendez, Double Beta Decay Matrix Elements, talk at the Interna-
tional Conference Neutrino 2020, June 28-July 2, 2020, Chicago, USA,
http://nu2020.fnal.gov.

[44] L. Wolfenstein, Neutrino Oscillations in Matter, Phys. Rev. D 17, 2369-

2374 (1978) doi:10.1103/PhysRevD.17.2369

31
[45] S. P. Mikheev and A. Y. Smirnov, Resonant amplification of neutrino
oscillations in matter and solar neutrino spectroscopy, Nuovo Cim. C 9,
17-26 (1986)

[46] S. Bilenky, On atmospheric neutrino mass-squared difference in the pre-

cision era, arXiv:1512.04172.

[47] I. Esteban, M. C. Gonzalez-Garcia, A. Hernandez-Cabezudo, M. Maltoni

and T. Schwetz, Global analysis of three-flavour neutrino oscillations:
synergies and tensions in the determination of θ23 , δCP , and the mass
ordering, JHEP 01, 106 (2019), arXiv:1811.05487 [hep-ph].

[48] A. Giuliani et al. [APPEC Committee], Double Beta Decay APPEC

Committee Report, arXiv:1910.04688 [hep-ex].

[49] M. Agostini et al. [GERDA], Probing Majorana neutrinos with double-β

decay, Science 365, 1445 (2019), arXiv:1909.02726 [hep-ex].

[50] A. Gando et al. [KamLAND-Zen], Search for Majorana Neutrinos near

the Inverted Mass Hierarchy Region with KamLAND-Zen, Phys. Rev.
Lett. 117, no.8, 082503 (2016), arXiv:1605.02889 [hep-ex].

[51] J. B. Albert et al. [EXO], Search for Neutrinoless Double-Beta De-

cay with the Upgraded EXO-200 Detector, Phys. Rev. Lett. 120, no.7,
072701 (2018), arXiv:1707.08707 [hep-ex].

[52] D. Q. Adams et al. [CUORE], Improved Limit on Neutrinoless Double-

Beta Decay in ?130 Te with CUORE, Phys. Rev. Lett. 124, no.12, 122501
(2020), arXiv:1912.10966 [nucl-ex].

[53] O. Azzolini et al. [CUPID-0], First Result on the Neutrinoless Double-

β Decay of 82 Se with CUPID-0, Phys. Rev. Lett. 120, no.23, 232502
(2018), arXiv:1802.07791 [nucl-ex].

[54] J. Detwiler, Future Neutrinoless Double Beta Decay Experiments, talk

at the International Conference Neutrino 2020, June 28-July 2, 2020,
Chicago, USA, http://nu2020.fnal.gov.

32