1. Introduction
The visible universe (e.g., stars, galaxies) appears to be composed of the ordinary matter just like everything one encounters on Earth or in the solar system, i.e., atomic nuclei made from protons and neutrons (both belonging to the class of particles called baryons) orbited by electrons. There is no evidence for the existence of localized pockets in the universe containing only antimatter, and asymmetry between matter and antimatter appears to be close to 100%. In the framework of the Big Bang genesis, the universe expands from a superdense high-temperature initial state. In these conditions, any existence of localized distinct clusters of matter and antimatter is unlikely. Rather, particles and antiparticles co-exist in an equilibrium where their numbers are replenished continuously via the creation processes and reduced when they annihilate each other in pairs.
Following the Big Bang, the temperature of the rapidly expanding universe was quickly decreasing. Eventually, the energy needed for the creation processes became too low and annihilation took over. Annihilation continued until the particles and antiparticles became sufficiently separated from each other to interact. At this point, the number densities of baryon
and antibaryon
were essentially fixed. However, in the absence of any asymmetry in the laws of nature governing the interactions of particles and antiparticles, these numbers would be catastrophically low. Normalized to the photon number density
, they come to [
1]
This is about
times smaller than the value of baryonic density obtained from the measurements of the abundances of the light elements or inferred from cosmic microwave background measurements [
2]:
Such annihilation catastrophe, where all the matter practically disappears, was evidently avoided. Starting from the symmetric state, without any excess of baryon over antibaryons, a tiny number of baryons survived and seeded the formation of planets, stars, galaxies, etc.—the baryon-dominated Universe. The three necessary conditions that any theory of the baryon creation, or baryogenesis, must fulfill in order that the Baryon Asymmetry of Universe (BAU)—the excess of baryons over antibaryons—could develop was formulated by A. Sakharov in 1967 [
3]. They are
existence of baryon number violation,
existence of Charge (C) and Charge-Parity (CP) symmetry violations,
departure from thermal equilibrium.
The baryon number violation is evidently required to arrive to net baryon excess. Violation of C and CP is also needed to assure that reactions for particles and antiparticles have different rates. Formally, CP transformation converts a particle to an antiparticle, and a theory that is symmetric under CP does not distinguish between the two. If CP is conserved during baryogenesis, every reaction that produces a particle has a counterpart that produces its antiparticle happening at exactly the same rate, thus leading to no net baryon number in the end. As the masses of particles and antiparticles are equal, their concentrations remain forever the same in the thermal equilibrium since energetically there is no preference for the production of one or the other, hence the need for third conditions.
The remainder of this review focuses on CP violation and in particular how it can be observed in long-baseline (LBL) neutrino oscillation experiments with conventionally produced neutrino (or antineutrino) beams at particle accelerators. Typically, such experiments use two detectors, one located near the beam origin (near detector) and the other placed far away (far detector), in order to assess the neutrino beam original characteristics (e.g., energy spectrum and flavor content) and measure its composition after significant flavor oscillations, respectively. A precise analysis of the oscillations of neutrinos in comparison to antineutrinos could then allow determining whether CP symmetry is preserved. The existence of CP violation in the lepton sector is one of the cornerstones of leptogenesis, a framework within which one tries to generate an asymmetry between lepton and antilepton numbers and to convert it into BAU.
2. Quark Mixing and CP Violation
The existence of a CP violation in the quark sector was first observed in 1964 in the decay properties of neutral kaons [
4]. Within the Standard Model (SM) of particle physics, CP violation is embedded in the Cabibbo–Kobayashi–Maskawa (CKM) mixing matrix,
, which relates that quark mass states participate in the flavor-changing weak interactions mediated by
boson exchange (charge current interactions). This
unitary matrix (
),
physically describes the linear combination of left-handed components of down-type quarks
d,
s,
b that participate in the charged current reactions with up-type quarks
u,
c,
t.
The CKM matrix can be parameterized in terms of three mixing angles
,
,
, and a complex phase
[
2]:
where
and
are the short-hands for
and
, respectively.
While non-vanishing
implies the violation of CP in the standard model, this phase in itself, however, is not a good measure. It is dependent on a particular choice of parametrization and could be shifted to different elements of the matrix by rephasing the fields. A parametrization-invariant measure of the CP violation can, however, be constructed from four products of CKM elements,
(
and
). For CP violation to occur,
must be complex and the imaginary part of these products non-vanishing. Moreover,
are all equal up to a sign to a constant called Jarlskog invariant [
5]. Without the loss of generality, one can define it to be given by
and with parametrization of Equation (
2), it can be expressed as
with
representing the maximum possible magnitude of the invariant when
.
The components of the CKM matrix have been extensively measured to a high precision from a variety of meson states. The values of the sines of the mixing angles currently are [
2]
The value of a complex CKM CP phase is [
2]
and that of the CKM Jarlskog invariant is [
2]
However, this (extremely well-measured) CP violation contained in the CKM matrix alone appears to be by far insufficient for BAU generation (e.g., [
1,
6]). It leads to a conjecture that there must be other sources of CP violation in nature, and significant experimental efforts are being dedicated to searching for any deviations from SM description within the quark sector. On the other hand, the evidence of the neutrinos having a non-zero mass, as a consequence to the observations of the neutrino oscillations, opened up an additional avenue to explore CP violation in the lepton sector.
3. Neutrino Mixing and CP Violation
Neutrino mixing arises from the fact that flavor neutrino states
,
, and
which take part in the weak charged current reactions with their charged counterparts
e,
,
are not the proper mass states, but a linear combination of such. Similar to quarks, the mixing is described by a unitary
matrix
U named after Pontecorvo, Maki, Nakagawa, and Sakahata (PMNS):
The elements of
, as for
, can be parameterized in terms of three mixing angles
,
,
and complex phase
yielding the same form as Equation (
2). However, since neutrinos are electrically neutral, they could potentially also be their own antiparticles: so-called Majorana fermions, as opposed to Dirac fermions, such as quarks and charged leptons. In that case, the mixing matrix contains two additional CP violating phases
and
:
where matrix
U has the same elements as Equation (
2). The LBL neutrino oscillation experiments have no sensitivity to the two
Majorana phases as they do not enter into the oscillation probabilities that one can measure. However, some possibility exists to probe their values [
7] if searches for the neutrinoless double beta decay process are successful and neutrinos are indeed Majorana particles.
Neutrino mixing angles show a strikingly different pattern from that of quarks. Based on the values from a global fit to experimental data [
8,
9], the sines of the PMNS mixing angles are
Comparing them to those of CKM (Equation (
6)), one can see that, in the case of neutrino mixing,
and in particular
are rather large. Consequently, if leptonic CP violation is at or near maximal,
, the value of the Jarlskog invariant is close to
and, with the latter currently at [
2]
the CP violation in the lepton sector could potentially be larger than that in the quark sector (Equation (
8)) by up to three orders of magnitude.
Apart from possibly hiding the source of a large CP violation, what also sets neutrinos apart from other SM elementary particles (quarks and charged leptons) is the extreme smallness of their masses. The scale of the neutrino mass could be probed via precise measurements of the endpoint of the energy spectrum of electrons emitted in
decays of tritium. The most stringent limit on the neutrino mass using this approach was set by the KATRIN experiment in 2022 [
10]:
The mass scale can also be constraints from cosmological data. For example, the Planck measurements of the cosmic microwave background set the upper limit on the neutrino masses at 0.12
[
11]. At this level, neutrinos appear to be lighter than electron—itself the lightest of the charged elementary particles—by a factor of almost
.
In the standard model, the particle masses are generated via the interaction with the Higgs boson that couples chirally left-handed fermion fields to right-handed ones. As the right-handed neutrinos have never been observed, within the standard model, neutrinos were naturally considered to be massless. However, the observations of the neutrino oscillations proved otherwise.
Given the smallness of the neutrino mass scale, it appears to be unlikely their masses can be generated by the same mechanism as other fermions. The most popular scenario for the explanation of light neutrino masses is a See-Saw mechanism [
12] that relies on the fact that neutrinos could be Majorana particles and associates an exceptionally heavy right-handed neutral (Majorana) partner
N to the light neutrino. In this model, as the name suggests, the lighter the
, the heavier the
N must be. With its mass somewhere in the range of 10
12 GeV to 10
15 GeV,
N cannot be created in particle accelerators. On the other hand, such particles could be produced after the Big Bang when the temperature of the early universe was still very high and their numbers were effectively frozen once the temperature decreased.
N could then decay both to charged leptons and antileptons. However, if leptonic CP violation exists, these decays would not have the same rates, leading eventually to an asymmetry between the number of leptons and antileptons that, in turn, could be converted into baryon–antibaryon asymmetry.
This procedure for generating BAU from an imbalance in lepton–antilepton numbers as a result of the decays of a heavy neutral particle is leptogenesis [
13]. While one cannot probe directly the existence of
N, we could hope to discover some supporting clues such as whether neutrinos are Majorana fermions and whether the CP symmetry is broken in the leptonic sector.
4. Neutrino Oscillations
Neutrino mixing gives rise to the phenomenon of neutrino oscillations whereby a beam of neutrinos of a given flavor () periodically changes its composition by exhibiting presence of other flavors as it propagates through space. In this context, the neutrino flavor is defined based on the type of the charged lepton that is produced (and detected) in charged current reactions.
The probability for neutrinos remaining in initial state
(survival) or oscillating to another flavor
(appearance) is determined by the elements of the PMNS mixing matrix and the differences between neutrino mass squares
(
). For the case of three active neutrinos (antineutrinos), it can be written compactly:
where
and
is the Kronercker delta,
E is the neutrino energy, and
L is the distance from a neutrino source to a detector (baseline of an experiment). The plus (minus) in front of the last term in Equation (
13) is for neutrinos (antineutrinos). For three neutrino mass states with masses
,
, and
, there are only two independent
. The smaller of the two
has the value of about
and historically has been associated with the measurements of solar neutrino fluxes. The other could be taken to be
or
and has a much larger value of around
. It was initially measured in the oscillation of atmospheric neutrinos propagating through Earth. While the sign of
is known, i.e., one can say that the state with mass
is lighter than the one with
, this is not the case for
for which only the absolute value is currently determined. Thus, two configurations for the hierarchy or ordering of neutrino masses are possible: one where the state with
is the heaviest and the other where it is the lightest. Traditionally, the former is referred to as Normal Ordering (NO), while the latter is called Inverted Ordering (IO).
The four products,
, in Equation (
13) for the survival probability (
) are real, and therefore CP violation can be detected only through appearance channels (
). The imaginary part of any non-zero four-product is equal to the invariant
from Equation (
4) up to a sign. Introducing asymmetry factor
one can trivially see that CP symmetry is preserved when
.
The LBL experiments so far, as well as the ones in the immediate future, rely on conventional neutrino beams [
14] where the neutrinos are produced mainly in the decays of
mesons generated in high-energy collisions of a primary accelerated proton beam with a stationary target. The charge selection of the pions is performed by a system of focusing elements (magnetic horns) allowing production of beams of mostly pure
or
. For CP measurement with such beams, the principal channel is
appearance. The asymmetry factor of interest is then
, which can be written as
For CP violation to be observed, not only the invariant
must be non-zero, but also the three neutrino masses must be distinct. The LBL experiments so far have been tuned to sample the oscillation probability at first maximum driven by
, that is to say when
. From Equation (
14), this condition is fulfilled when
where
is the first oscillation maximum,
is the second, etc. Taking
at the first maximum, the value that the asymmetry in Equation (
16) can attain there is
which is about
when
.
The asymmetry expressed by Equation (
16) is only valid for the neutrinos propagating in a vacuum. In reality, the neutrino beams travel through Earth along the chords that penetrate the crust at varying depths. The neutrinos can then undergo coherent scattering from in the atomic electrons in the surrounding medium, which introduces additional phases in the oscillation amplitudes that are different for neutrinos and antineutrinos. The effect of the matter is included in oscillation probability calculations via introduction of effective matter potential
V, which is proportional to the number density of electrons
in the encountered medium:
where
is the Fermi constant.
Assuming a uniform average electron density, that is to say, constant along the beam path, the oscillation probability for the
appearance channel can be approximated by [
15]
where
and
are defined in Equations (
5) and (
14), respectively, while
For the equivalent antineutrino channel, one must flip the signs of the CP phase and the effective matter potential:
and
.
The first term in Equation (
20) is the leading contribution to the oscillation probability for this channel that allowed the measuring of the
mixing angle in LBL experiments [
16,
17]. Given the sizes of
and
, the second term has a negligible contribution to the oscillation probability as
. The last term introduces sub-leading effects dependent on the value of the
phase, which one then must measure. Additionally, since the matter effects introduce a phase to the oscillations driven by
, the ordering of the neutrino masses could also be determined.
The matter effect generates an asymmetry for neutrino and antineutrino oscillations and its impact grows with the distance the neutrino beam travels (i.e., amount of material it encounters) and neutrino beam energy. The electron number density in Equation (
19) is proportional to the matter density as
, where the factor of
is an approximate value for the ratio of electron to nucleon number densities for ordinary matter and
is the Avogadro number. Dimensionally, the numerator of Equation (
21) must have the same units as
and, with the latter measured in
, one obtains
Equations (
23) and (
21) show the proportionality of matter-induced phase
A to neutrino energy. Alternatively, Equation (
20) can be recast to make the dependence on baseline length
L explicit by expanding
products. Following [
18], one can define distance scale parameter
a,
with a positive (negative) sign for neutrinos (antineutrinos). Then, with
the oscillation probability of Equation (
20) can be rewritten as
With
L typically being in
, the units of
a must be cast to
. Switching from natural units by re-introducing
ℏ and
c factors one finds
As an example, taking a matter density of
(an approximate value if the beam travels through the upper crust layer of the Earth), one obtains
1/3700
as a measure of a distance scale.
The current and the next-generation LBL neutrino oscillation experiments have baseline lengths ranging from about 300 km to 1300 km.
Figure 1 attempts to illustrate qualitatively the evolution of the impact of the matter effect on the
appearance probability as a function of
L. The vertical axis in the figure shows a mean asymmetry factor,
, computed from the averaged oscillation probabilities for neutrinos and antineutrinos over an energy window defined by a full width at half maximum of the first oscillation peak:
This is rather a simplification, as typically one does not perform a single measurement but attempts to fit the shape of the probability by sampling around its maximum value in an energy window defined by the beam energy spectrum, event reconstruction efficiencies and selection cuts, etc. However, the plot still serves as a convenient illustration for a number of features.
Each band in
Figure 1 represents the different possible values of
depending on the value of the
phase for a given mass ordering as a function of baseline length. Some representative cases for the choices of
like the CP-conserving values
or those with
corresponding to the maximal CP violation are also superimposed. Without the knowledge of mass ordering, the experiments with shorter baselines (where the two bands overlap) lose sensitivity to parts of the
parameter space due to a degeneracy introduced by the two unknowns,
and sign of
. In particular, one could even obtain similar asymmetry generated by CP-conserving phase values,
, in one mass ordering scenario as that produced by phases
(i.e., the phase values corresponding to the maximal CP violation) for the other mass ordering. In
Figure 1, such cases correspond to the points where edges of a given band intersect the
lines in the other band. Without the knowledge of mass ordering, one could incorrectly infer the existence or non-existence of CP violation.
As the effect of matter grows with baseline length, acquired phase
in Equation (
26) causes increasing enhancement (suppression) or suppression (enhancement) of appearance probability for neutrinos (antineutrinos) depending on the underlying ordering of neutrino masses. The two bands in
Figure 1 begin to separate at about 800 km until eventually they are so far apart that one can readily distinguish the two possible mass ordering scenarios and measure
at the same time. The matter effect could also lead to breaking the degeneracy between the two CP-conserving
phase values as illustrated by an increasing separation of the two corresponding lines for each choice of mass ordering.
For further illustration,
and
appearance probabilities at the first oscillation maximum for a 1300 km baseline are also plotted in
Figure 2 for the CP conserving case and the maximal CP violation as well as both mass orderings. The enhancement or suppression, depending on the sign of
, of
or
appearance can be clearly seen.
One could notice from Equation (
26) that the second and third (
and CP phase-sensitive) terms disappear when product
is a multiple of
. The first such “magic” baseline is
and, depending on assumed density profile, its value can be
to 8000 km. Experiments searching for neutrino electron appearance on magic baselines were suggested early on to help disentangle contributions from multiple unknowns (
,
, sign of
) to the oscillation probability for this channel.
5. Experimental Landscape
Spurred by the discovery that the value of
, the last mixing measured angle, is rather large, the neutrino oscillation experiments are now entering the precision era focusing on the remaining unknowns such as the existence of CP violation, the ordering of neutrino masses, and determination where
is with respect to
(octant). However, to offer a little historical perspective of this development, one could recall that in the first decade of this century, the knowledge of
was limited to the exclusion region set in 2003 by the Chooz reactor experiment [
19] looking at the survival of electron antineutrinos produced in nuclear reactors. The experiment was configured (the source–detector distance) to probe the oscillations driven by large
that was also measured in a Super-Kamiokande detector from the observations of the disappearance of
component of atmospheric neutrinos [
20].
When
sets dominant frequency, the reactor neutrino survival probability is approximately given by
and this channel provides the most direct access for measuring
. Unfortunately, with
, Chooz was more sensitive to larger
than what it turned out to be. It saw no evidence for the disappearance of reactor antineutrinos over the baseline distance of about a kilometer excluding the values of
greater than about 0.1.
The non-zero value of
is required for the possibility of the existence of leptonic CP violation, and a next generation of experiments, both accelerator and reactor-based ([
21,
22,
23]), was then needed to look for this mixing angle. In Japan, a T2K experiment was proposed [
24] to use a new high-power J-PARC accelerator for the
(
) neutrino source with Super-Kamiokande, a 50
water detector, as the far detector 295 km away. Its goal was to probe
survival and
appearance channels at the first oscillation maximum with a neutrino beam tuned in a narrow energy window (narrow-band beam) around
using the off-axis technique [
25], where the far detector baseline is offset relative to the beam axis by a small angle to decrease the neutrino beam energy spread. In particular for the appearance channel, which to the first order depends on the product of
(Equation (
20)), the aim was to push the limit on
by more than one order of magnitude below
.
In the USA, an accelerator-based neutrino experiment NOvA was proposed [
26] with an off-axis neutrino beam originating from Fermilab and travelling 810 km to a far site equipped with a 14
liquid scintillator detector. Similar to T2K, one of the principal aims was to discover
also with a possibility to probe the mass ordering exploiting the larger matter effect due to the longer baseline.
As an alternative to conventional neutrino beams, the concept of muon storage rings that could be used to produce intense and easily characterized neutrino beams (neutrino factory) was pursued through a rich R&D program (e.g., [
27,
28]). If
turned out to be quite small,
, the neutrino factory as a source appeared as a viable solution for exploring CP violation in the lepton sector. As a substitute for stored muons, ion storage rings for
decaying isotopes producing electron (anti)neutrino beams (beta beams) were also proposed [
29].
In 2011, T2K announced a first hint of a non-zero
[
30]. The value itself was reported a year later by the Daya Bay reactor neutrino experiment [
31] at more than
significance, and it was
which happened to be just below the Chooz exclusion region. Since then, it was measured with increasing precision, particularly by the reactor experiments, with the Daya Bay final result determining
to be better than 3% [
32]:
With such a large and tightly constrained value, the measurement of leptonic CP violation with conventionally produced neutrino beams are a distinct possibility.
In 2020, T2K reported [
33] first indication of a possible existence of CP violation in the lepton sector. Their data disfavored the CP conserving values of phase
at a 95% confidence level showing a preference for
for close maximal CP violation instead. Further analysis of the data collected by the experiment so far [
34] put the constraint on
at
for T2K favored solution of normal mass ordering and including the constraint on
from reactor experiments.
Figure 3 from [
34] illustrates graphically the preference of the T2K data collected with the neutrino (abscissa) and antineutrino (ordinate) beam modes for the solutions that lead to the largest asymmetry in favor of
(
) rather than
appearance (
). The preference for
, which tends to simultaneously increase the expected number of events in both sets, is also visible.
While T2K appears to favor the existence of close to maximal CP violation with
, the results reported by the NOvA experiment do not share the same trend [
35]. The analysis of the data collected up to 2020 does not exhibit strong asymmetry between electron neutrino and antineutrino appearance probabilities.
Figure 4 from [
34] shows the NOvA 90% confidence level contours in the
-
parameter plane for each mass ordering scenario compared to T2K as well as the results from the Super-Kamiokande analysis of atmospheric neutrino sample [
36]. The NOvA best fit point for
is close to the CP-conserving value
:
The T2K preferred solution for and NO appears to be partially disfavored by the current NOvA data, although an overlap between two results remains. In the case of IO, NOvA rules out maximal CP violation at at more than level. On the other hand, its contours for this mass ordering scenario contain entirely the corresponding region from T2K. Although IO is disfavored both by T2K and Super-Kamiokande, in the absence of any further strong constraint on mass ordering, the solution with inverted mass ordering and CP violation close to maximal with is not ruled out. Both NOvA and T2K experiments are scheduled to continue collecting data through 2026 at least doubling statistics of their datasets. This potentially could allow mass ordering reaching sensitivity and exclude CP conservation confidence level within the next few years.
The search for leptonic CP violation is one of the key elements of the physics programs for the next generation of the LBL neutrino oscillation experiments, Hyper-Kamiokande [
37] and DUNE [
38], both currently in a construction stage.
The Hyper-Kamiokande detector is a water detector about five times larger than Super-Kamiokande with a total detector mass of 258 . The detector is located at the same off-axis angle with respect to the T2K beamline as the Super-Kamiokande as well as the same 295 km distance from J-PARC where the planned upgrades to the accelerator complex are expected to raise the beam power to .
Figure 5 shows projected sensitivity of the Hyper-Kamiokande experiment as a function of a number of years of running to exclude CP conservation (
). The constraint, whether from external measurements or analysis of atmospheric neutrino data sample collected Hyper-Kamiokande, on mass ordering is required in order to maximize the sensitivity of the experiment. With the knowledge of mass ordering and if CP violation is close to maximal at
, the value currently favored by T2K, the experiment can discover (
level significance) leptonic CP violation within three years of running with the start of the operation currently projected to be in 2027.
The US-based DUNE experiment utilizes large liquid argon time-projected chambers as the neutrino detectors for the analysis of the neutrino beam produced at the Fermilab accelerator complex located about 1300 km away. The design foresees four modules each containing 17
of liquid argon. The construction is organized in phases. Within the first phase, the first two modules are scheduled to be in operation by mid-2030 [
40] and the neutrino beam running in 2031. The beam power is progressively raised to 1.2
and further upgraded to above 2
as part of the second phase of the experiment also featuring an addition of two more detector modules.
Given the length of the DUNE baseline, the large matter-induced asymmetry allows for the experiment to resolve the mass ordering within few years of running in the first phase of the program independent of the
value. If CP violation is nearly maximal, the experiment could also potentially see it at a level of
in the first phase. However, to increase the precision and enlarge the coverage of the
parameter space, implementation of the second phase of the project is required. This is illustrated in
Figure 6, which shows the evolution of the sensitivity of DUNE to CP violation as a function of time.
Both DUNE and Hyper-Kamiokande are poised to collect a large sample of neutrino events. However, the precision measurement of CP violation also requires a tight control of systematic uncertainties for both experiments. The principal sources of the latter are related to the knowledge of the neutrino fluxes produced at the accelerators and the cross-sections for neutrino interactions. Both of these impact the rates as well as the (energy) spectra of events observable in far detectors from which one tries to infer the parameters (e.g., ) determining the underlying neutrino oscillation probability. Dedicated measurements of hadron production yields from different target materials help put constrains on flux uncertainties. However, a judicious design of near detectors is also required to amply characterize the properties of the neutrino beam at the production point and constraint the interaction models in the neutrino beam energy window. From this information, one then formulates the prediction for the expected observable at the far detector against which the collected data could be compared in order to extract the values of the underlying neutrino mixing parameters.
With largely different baselines of DUNE and Hyper-Kamiokande, there is great complimentarity between the two experiments with each measuring neutrino oscillation in distinct energy regions and being subject to different magnitudes of the matter effects. Ultimately, the combination of results from the two experiments can help shed light on a number of remaining unknowns in the mixing of massive neutrinos in addition to answering the question of whether CP violation exists in the lepton sector.