2.1. Massive Axions in the Peccei and Quinn Picture and Neutrino–Axion Couplings
The gauge field tensor of the QCD Lagrangian is written [
1]
where
is the gauge field tensor,
and
are structure constant indexes,
and
are Lorentz indexes, and
is a constant matrix. If the CP and T invariances are not assumed, one may add the term
where
are the elements of another constant matrix, and
is the Levi-Civita tensor. The term
induces a neutron electric dipole moment, given by the expression
, where
and
are the masses of the pion and of the neutron, respectively, and
e is the absolute value of the electron charge.
Then, if the neutron electric dipole moment is of the order of (or smaller than)
e cm [
11], it implies the constraint
.
The following are a few words about the CP properties of the electroweak Hamiltonian. The total Hamiltonian may be decomposed into two parts, one even (
) and one odd (
), under CP transformations, like
, where
[
12]. The present evidences indicate that the odd component under CP is much weaker than the even one. Then, as explained in [
12], when the ratio between these components of the electroweak Hamiltonian is fixed at the super weak scale (
, the neutron electric dipole moment should then be, as we have mentioned before, of the order of (or smaller than)
e cm.
In order to solve this problem, R. Peccei and H. Quinn (1977) [
2,
3] have proposed the inclusion of a pseudo-scalar field
, the axion, such that
f being a strength constant. The non-vanishing vacuum expectation value of , with the subsequent breaking of the U(1) symmetry associated with it, i.e., , causes to vanish, .
This assumption was later extended, separately, by S. Weinberg and F. Wilczeck in 1978 [
4,
5], who wrote the Lagrangian
where
is a neutral scalar boson field, e.g., the axion, and the last two terms are the interactions of up (
u) and down (
d) quarks with axions.
The departure from the Peccei–Quinn axion [
2,
3] is just the transformation of this Lagrangian to a pion–axion basis. From it, the axion mass naturally arises from triangular vertices of the axion–two-pions type mediated by quark and antiquark loops. This is possible because of the quark–antiquark pair structure of the pions.
The effective Lagrangian which describes this process is written
where
is a 2 × 2 mass-matrix and
. The eigenvalues of
are
and
. The axion mass then becomes [
1]
From the point of view of a fundamental symmetry, the axion of Peccei and Quinn is the Goldstone boson of the
U(1) symmetry which breaks at the scale
GeV, and in the Weinberg–Wilczeck [
4,
5] formulation, it results from the condition
.
With values of the mass of the axion smaller than fractions of meV, its lifetime is of the order of s, which is more than enough to travel cosmological distances without decaying.
The scale factor
f has lower limits which vary from
GeV, based on symmetry arguments, to
GeV, which is the typical value associated with the Peccei–Quinn scale [
1,
2]. Both the Weinberg [
4] and Wilczeck [
5] approaches determine the relevant role of the scale factor
f, and the corresponding limits may be extracted from the measurement of the decay of axions into two pions. In a similar manner, a term describing the decay of axions into two X-ray photons, that is, a triangular vertex mediated by electron–positron pairs, may be added to the effective Lagrangian of Equation (
4). The coupling of neutrinos and axions may also be a way to determine that scale, as we shall discuss next.
Most of the matter in the Universe is dark [
13,
14], its existence is manifest from astronomical evidence [
14]. Basically, it is non-baryonic and collision-less, but its composition is unknown.
Axions may be a dominant part of the cold dark matter, as postulated by Sikivie and Yang in their original paper [
9] and by other authors [
15,
16,
17,
18,
19,
20].
Among the experiments devoted to the direct detection of dark matter particles, we should mention ADMeX (Axion Dark Matter electron-X), which aims to detect axions by the measurement of X-rays produced by the interaction of axions with electrons, mediated by electron–positron pairs. It is the equivalent of the production of two pions by the interaction between axions and quarks, mediated by quark–antiquark pairs, which is allowed as a second-order process by the effective Lagrangian of Weinberg and Wilczek [
4,
5].
If one assumes that axions are indeed the main component of dark matter, and that they have a non-localized distribution in space, it could be possible to couple them with neutrinos by means of a derivative term in the axion sector coupled to the neutrino current. Then, by taking the non-zero vacuum expectation value of the axion field, it results in a mass term for the neutrinos, as it is shown next.
By adding to the Lagrangian of axions an axion–neutrino coupling term, it could be possible to give mass to the neutrinos, as we shall explain next. From From a more fundamental point of view, the extension of the Lagrangian to include the interactions between neutrinos and axions follows from the work of Weinberg [
4] and Wilczek [
5]. It is expressed like Equation (
4), by taking the four-potential as the covariant derivative of the axion field (
and the four-current
of the neutrino sector. Since the time derivative of the axion field in the proper frame is just the mass of the axion, we define the coupling constant
scaled by the mass of the axion. We start from the Lagrangian
which describes the derivative coupling between neutrinos (
) and axions (
), with
being the strength of their coupling divided by the mass of the axion (remember that we have taken
). The time derivative of the axion field then cancels out that mass dependence (see below) since, in the proper frame,
. As explained later, the mass scale resulting from the time derivative is given in terms of the coupling
and of the expectation value of the axion field
(see Equation (
9)).
The breaking of the
U(1) symmetry of the axion sector is represented by the Higgs-like potential
The variation of this potential leads to the extremes
and
As a consequence of it, from the structure of the Lagrangian
, the time derivative and the spatial derivatives can be written separately, leading to the expression
The second term of Equation (
11) can be interpreted as a mass term, after performing the time derivative of the axion field in the proper frame and giving to the axion the non-zero extreme value of Equation (
10). This mechanism is analogous to the conventional Higgs mechanism. The breaking of this extra symmetry, represented by the non-zero value of the expectation value of the axion field, as determined by the potential of Equation (
8), gives mass to the neutrinos. Going beyond this level of approximation requires the exchange of axions at one loop level. The corrections to the zeroth-order mass of the neutrinos
are represented by the mass propagator
, which depends on the momentum exchanged between the axion and the neutrino. In the next paragraphs, we shall show the final results of the calculation of the mass propagator. The calculation is based on the use of conventional methods of quantum field theory and the results can be expressed in terms of scalar (
) and vector (
) components of the propagator.
The physical mass of the neutrino can be computed as
After evaluating
on shell, that is, by taking
, the one-loop correction to the neutrino mass due to the interaction with axions is contained in the kernel
[
21]
where its vector and scalar components are
and
, respectively.
The one-loop neutrino propagator is then written [
10,
21]
The derivation of the previous equations involves the ordering of higher-order corrections to the propagator, as well as fixing the value of the coupling for each mass scale of the axion.
The following are some final words about neutrino–axion coupling:
- (i)
The breaking of the
U(1) symmetry proposed by Peccei and Quinn,
U(1)
PQ, at the level of the Lagrangian which describes the interaction between the axion and the neutrino, at the zeroth order, gives mass to the neutrino (Equation (
12)). That mass is dependent upon the coupling constant of the Lagrangian (
) and of the constant (
f), which determines the value of the mass of the axion;
- (ii)
The one-loop corrections to the zeroth-order neutrino mass are also dependent upon these constants, but they are non-divergent (Equations (
13)–(
15));
- (iii)
In order to complete the scheme, one has to take into account the squared mass differences between the three light-mass eigenstates (both in the normal and inverse ordering) and the amplitudes relating the mass and flavor states in the light-mass sector, as well as the amplitudes of the heavy-mass sector.
To date, the calculations have been restricted to neutrino light-mass eigenstates, but to be consistent with the claim that the representation should have a right-handed channel, one may also have to include heavy-mass neutrinos in the picture, in order to express left- and right-handed lepton doublets [
7].
2.2. Extensions of the Minimal Standard Model of Electroweak Interactions
The minimal version of the Standard Model of electroweak interactions [
1] is based on purely left-handed currents, both for fermions and bosons, massless neutrinos, left-handed lepton doublets for negatively charged
,
, and
, with the corresponding antineutrinos
,
, and
, and singlets of
,
, and
leptons.
The conservation rules of the minimal representation imposed lepton number and lepton flavor conservation. The existence of neutrino flavor oscillations were confirmed experimentally by the SNO and Kamiokande collaborations, the results for which A. McDonald and T. Kajita received the Nobel Prize for Physics in 2015. The existence of neutrino mass oscillations implies the existence of massive neutrinos.
The coefficients of the linear combination between three light-mass neutrino mass eigenstates
are defined in the expression
for the electron neutrino.
The electron neutrino sector depends on the amplitudes and CP phases of the linear combination of neutrino mass eigenstates
and
, respectively; similar expressions are written for the other two flavors.
Table 1 gives the values of the mixing angles and squared mass differences between neutrino mass eigenstates, for normal (NH) and inverse (IH) mass hierarchies.
Starting from these expressions, the amplitudes and probabilities associated with the flavor conversion from flavor
at
to flavor
at time
t (
,
) are written:
and
respectively, where the symbol
is the squared mass differences between the mass eigenstates
k and
.
The pattern of neutrino oscillations in a vacuum may differ from the pattern of neutrino oscillations in the presence of interactions between neutrinos and, for instance, dark matter particles. The interaction mediated by the exchange of virtual massive particles is depicted in the diagram given in
Figure 1 and the results are shown in the panels of
Figure 2. The details of the calculations are given in [
22].
Table 2 shows the results for the survival and conversion probabilities between the three neutrino flavors, where we have included the results obtained by diagonalizing neutrino–neutrino interactions.
In addition to these oscillations between neutrinos of different flavors, which do not require any interaction between the neutrinos or between neutrinos and hadrons, there exist other processes which, if observed, would demonstrate the need to extend the minimal version of the Standard Model. They can be possible if neutrinos are massive. They are the neutrinoless double beta decay, which violates lepton number conservation, and the neutrino flavor violation. Both processes require the interaction of nucleons and pair of leptons, mediated by charged and/or neutral bosons. Among these exotic processes, the neutrinoless double beta decay (
) has been, and still is, the object of the attention of theoreticians and experimentalists for several decades [
7]. In the following paragraphs, we shall introduce the basic elements of the theory related to the double beta decay processes. This process consists of the decay of a mother nucleus with N neutrons and Z protons into a daughter nucleus with N − 2 neutrons and Z + 2 protons, accompanied by the emission of two electrons with an energy equal to the Q-value of the decay. The decay schemes are shown in
Figure 3.
The experimental signal would be the detection of the two electrons flying in opposite directions. The decay implies the exchange of a massive neutrino between the two vertices. It can also take place if right-handed currents are included in the electroweak Lagrangian, also implying the existence of left–right couplings. The observation of this process indicates the need to extend the Standard Model, because it is forbidden, in the current formalism, to accommodate massive neutrinos and/or right-handed currents with the corresponding triplet of right-handed bosons.
Values of the electron neutrino mass can be extracted by comparing the theoretical rates and the experimental lower limits for the non-observation of the neutrinoless double beta decay (
) half-life, which are of the order of
years or larger [
7].
By using average values of the nuclear matrix elements, evaluated in the context of different nuclear structure models, the resulting average neutrino mass is of the order of 0.5 eV, but this value is very much dependent on the nuclear structure models used to evaluate the nuclear matrix elements of the participant electroweak operators.
The expression of the inverse half-life of the decay is separable into three factors, which represent the contributions coming from the neutrino sector, the nuclear structure sector, and phase space factors.
The non-standard electroweak Hamiltonian, which includes left- and right-handed currents and mediators, is written
The simpler mechanism to explain lepton number violation consists of the second-order treatment of the first term of the Hamiltonian of Equation (
19) acting on two neutrons belonging to the initial nucleus, leading to the transformation of them into two protons, followed by the emission of two electrons [
7], as is shown by the diagram of
Figure 3. The matrix elements of the Hamiltonian of Equation (
19) on nuclear states, for operators of the Fermi (F), Gamow–Teller (GT), and Tensorial (T) types are written
The calculation of these matrix elements implies the definition of a neutrino potential
with matrix elements between nuclear states of the form
with two-particle transition amplitudes given by the expressions
In these equations, the index K stands for Fermi, Gamow–Teller, and Tensorial operators.
With them, one can write transition densities depending on nuclear wave functions. In the framework of the Quasiparticle Random Phase Approximation (QRPA) [
7], as an example, they acquire the form
for all possible terms resulting from the multipole expansion of the neutrino potential.
The resulting expression for the half-life, restricted to the mass sector, is given by Equation (
27):
The nuclear structure sector
is written in terms of the matrix elements of the multipole operators which participate in the transitions between nuclear states. The so-called mass term is expressed in terms of nuclear matrix elements of the Gamow–Teller operator
:
The form factors
are written in terms of radial integrals of electron wave functions [
6]. Unlike the double beta decay with neutrinos (
), which is observed with half-lives some orders of magnitude faster, and which is not dependent on neutrino properties [
7], the values of the nuclear matrix elements which participate in the
decay are much more stable along the nuclear systems, with values which are of the order of 3–5, much larger than the values needed to reproduce the
channel, which are strongly suppressed.
The theory and experiments related to both the and decays have been intensively developed during the last three decades.
Some results of the calculated neutrino mass, extracted from the systematics on calculated half-lives for
decay processes, are shown in
Section 3. The cases are those of the decays of the nuclei
and
, which have been, and still are, the object of leading experimental efforts [
23,
24]. For the case of
decay, we show, as examples, the results corresponding to the allowed decay of the mother nuclei
.