Status and future of nuclear matrix elements for neutrinoless double-beta decay: a review

Experimental neutrino physics in a nuclear landscape

There are profound connections between neutrino physics and nuclear experiments. Exceptionally precise measurements of single and double beta-decay spectra illuminate the scale and nature of neutrino mass and may finally answer the question of whether neutrinos are their own anti-matter counterparts. Neutrino-nucleus scattering underpins oscillation experiments and probes nuclear structure, neutrinos offer a rare vantage point into collapsing stars and nuclear fission reactors and techniques pioneered in neutrino nuclear physics experiments are advancing quantum sensing technologies. In this article, we review current and planned efforts at the intersection of neutrino and nuclear experiments. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

Hexadecapole Deformation of ^{238}U from Relativistic Heavy-Ion Collisions Using a Nonlinear Response Coefficient

The hexadecapole deformation (β_{4}) of the ^{238}U nucleus has not been determined because its effect is overwhelmed by those from the nucleus' large quadrupole deformation (β_{2}) in nuclear electric transition measurements. In this Letter, we identify the nonlinear response of the hexadecapole anisotropy to ellipticity in relativistic U+U collisions that is solely sensitive to β_{4} and insensitive to β_{2}. We demonstrate this by state-of-the-art hydrodynamic calculations and discuss the prospects of discovering the β_{4} of ^{238}U in heavy-ion data at the Relativistic Heavy Ion Collider.

Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in ^{76}Ge

The observation of neutrinoless double-beta (0νββ) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in 0νββ experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the 0νββ-decay NME, in the key nucleus ^{76}Ge. Our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of 2.60_{-1.36}^{+1.28}, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of 187_{-62}^{+205}  meV. The result is important for designing next-generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.

Measurement of Electron-Neutrino Charged-Current Cross Sections on ^{127}I with the COHERENT NaIνE Detector

Using an 185-kg NaI[Tl] array, COHERENT has measured the inclusive electron-neutrino charged-current cross section on ^{127}I with pion decay-at-rest neutrinos produced by the Spallation Neutron Source at Oak Ridge National Laboratory. Iodine is one the heaviest targets for which low-energy (≤50  MeV) inelastic neutrino-nucleus processes have been measured, and this is the first measurement of its inclusive cross section. After a five-year detector exposure, COHERENT reports a flux-averaged cross section for electron neutrinos of 9.2_{-1.8}^{+2.1}×10^{-40}  cm^{2}. This corresponds to a value that is ∼41% lower than predicted using the MARLEY event generator with a measured Gamow-Teller strength distribution. In addition, the observed visible spectrum from charged-current scattering on ^{127}I has been measured between 10 and 55 MeV, and the exclusive zero-neutron and one-or-more-neutron emission cross sections are measured to be 5.2_{-3.1}^{+3.4}×10^{-40} and 2.2_{-0.5}^{+0.4}×10^{-40}  cm^{2}, respectively.

Measurement of the 2νββ Decay Rate and Spectral Shape of ^{100}Mo from the CUPID-Mo Experiment

Neutrinoless double beta decay (0νββ) is a yet unobserved nuclear process that would demonstrate Lepton number violation, a clear evidence of beyond standard model physics. The process two neutrino double beta decay (2νββ) is allowed by the standard model and has been measured in numerous experiments. In this Letter, we report a measurement of 2νββ decay half-life of ^{100}Mo to the ground state of ^{100}Ru of [7.07±0.02(stat)±0.11(syst)]×10^{18}  yr by the CUPID-Mo experiment. With a relative precision of ±1.6% this is the most precise measurement to date of a 2νββ decay rate in ^{100}Mo. In addition, we constrain higher-order corrections to the spectral shape, which provides complementary nuclear structure information. We report a novel measurement of the shape factor ξ_{3,1}=0.45±0.03(stat)±0.05(syst) based on a constraint on the ratio of higher-order terms from theory, which can be reliably calculated. This is compared to theoretical predictions for different nuclear models. We also extract the first value for the effective axial vector coupling constant obtained from a spectral shape study of 2νββ decay.

Final Results of GERDA on the Two-Neutrino Double-β Decay Half-Life of ^{76}Ge

We present the measurement of the two-neutrino double-β decay rate of ^{76}Ge performed with the GERDA Phase II experiment. With a subset of the entire GERDA exposure, 11.8 kg yr, the half-life of the process has been determined: T_{1/2}^{2ν}=(2.022±0.018_{stat}±0.038_{syst})×10^{21}  yr. This is the most precise determination of the ^{76}Ge two-neutrino double-β decay half-life and one of the most precise measurements of a double-β decay process. The relevant nuclear matrix element can be extracted: M_{eff}^{2ν}=(0.101±0.001).

^{138}Ba(d,α) Study of States in ^{136}Cs: Implications for New Physics Searches with Xenon Detectors

We used the ^{138}Ba(d,α) reaction to carry out an in-depth study of states in ^{136}Cs, up to around 2.5 MeV. In this Letter, we place emphasis on hitherto unobserved states below the first 1^{+} level, which are important in the context of solar neutrino and fermionic dark matter (FDM) detection in large-scale xenon-based experiments. We identify for the first time candidate metastable states in ^{136}Cs, which would allow a real-time detection of solar neutrino and FDM events in xenon detectors, with high background suppression. Our results are also compared with shell-model calculations performed with three Hamiltonians that were previously used to evaluate the nuclear matrix element (NME) for ^{136}Xe neutrinoless double beta decay. We find that one of these Hamiltonians, which also systematically underestimates the NME compared with the others, dramatically fails to describe the observed low-energy ^{136}Cs spectrum, while the other two show reasonably good agreement.

β^{+} Gamow-Teller Strengths from Unstable ^{14}O via the (d,^{2}He) Reaction in Inverse Kinematics

For the first time, the (d,^{2}He) reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the β^{+} direction from an unstable nucleus. The new technique was made possible by the use of an active-target time-projection chamber and a magnetic spectrometer, and opens a path to addressing a range of scientific challenges, including in astrophysics and neutrino physics. In this Letter, the nucleus studied was ^{14}O, and the Gamow-Teller transition strength to ^{14}N was extracted up to an excitation energy of 22 MeV. The data were compared to shell-model and state-of-the-art coupled-cluster calculations. Shell-model calculations reproduce the measured Gamow-Teller strength distribution up to about 15 MeV reasonably well, after the application of a phenomenological quenching factor. In a significant step forward to better understand this quenching, the coupled-cluster calculation reproduces the full strength distribution well without such quenching, owing to the large model space, the inclusion of strong correlations, and the coupling of the weak interaction to two nucleons through two-body currents.

Final Result of the Majorana Demonstrator's Search for Neutrinoless Double-β Decay in ^{76}Ge

The Majorana Demonstrator searched for neutrinoless double-β decay (0νββ) of ^{76}Ge using modular arrays of high-purity Ge detectors operated in vacuum cryostats in a low-background shield. The arrays operated with up to 40.4 kg of detectors (27.2 kg enriched to ∼88% in ^{76}Ge). From these measurements, the Demonstrator has accumulated 64.5 kg yr of enriched active exposure. With a world-leading energy resolution of 2.52 keV FWHM at the 2039 keV Q_{ββ} (0.12%), we set a half-life limit of 0νββ in ^{76}Ge at T_{1/2}>8.3×10^{25}  yr (90% C.L.). This provides a range of upper limits on m_{ββ} of (113-269) meV (90% C.L.), depending on the choice of nuclear matrix elements.

Determining g_{A}/g_{V} with High-Resolution Spectral Measurements Using a LiInSe_{2} Bolometer

Neutrinoless double beta decay (0νββ) processes sample a wide range of intermediate forbidden nuclear transitions, which may be impacted by quenching of the axial vector coupling constant (g_{A}/g_{V}), the uncertainty of which plays a pivotal role in determining the sensitivity reach of 0νββ experiments. In this Letter, we present measurements performed on a high-resolution LiInSe_{2} bolometer in a "source=detector" configuration to measure the spectral shape of the fourfold forbidden β decay of ^{115}In. The value of g_{A}/g_{V} is determined by comparing the spectral shape of theoretical predictions to the experimental β spectrum taking into account various simulated background components as well as a variety of detector effects. We find evidence of quenching of g_{A}/g_{V} at >5σ with a model-dependent quenching factor of 0.655±0.002 as compared to the free-nucleon value for the interacting shell model. We also measured the ^{115}In half-life to be [5.18±0.06(stat)_{-0.015}^{+0.005}(sys)]×10^{14}  yr within the interacting shell model framework. This Letter demonstrates the power of the bolometeric technique to perform precision nuclear physics single-β decay measurements, which along with improved nuclear modeling can help reduce the uncertainties in the calculation of several decay nuclear matrix elements including those used in 0νββ sensitivity calculations.

Determination of Double Beta Decay Half-Life of 136Xe with the PandaX-4T Natural Xenon Detector

Precise measurement of two-neutrino double beta decay (DBD) half-life is an important step for the searches of Majorana neutrinos with neutrinoless double beta decay. We report the measurement of DBD half-life of 136Xe using the PandaX-4T dual-phase Time Projection Chamber (TPC) with 3.7-tonne natural xenon and the first 94.9-day physics data release. The background model in the fiducial volume is well constrained in situ by events in the outer active region. With a 136Xe exposure of 15.5 kg-year, we establish the half-life as 2.27 ± 0.03(stat.) ± 0.10(syst.) × 1021 years. This is the first DBD half-life measurement with natural xenon and demonstrates the physics capability of a large-scale liquid xenon TPC in the field of rare event searches.

Final Result on the Neutrinoless Double Beta Decay of ^{82}Se with CUPID-0

CUPID-0, an array of Zn^{82}Se cryogenic calorimeters, was the first medium-scale demonstrator of the scintillating bolometers' technology. The first project phase (March 2017-December 2018) allowed the most stringent limit on the neutrinoless double beta decay half-life of the isotope of interest, ^{82}Se, to be set. After a six month long detector upgrade, CUPID-0 began its second and last phase (June 2019-February 2020). In this Letter, we describe the search for neutrinoless double beta decay of ^{82}Se with a total exposure (phase I+II) of 8.82  kg yr^{-1} of isotope. We set a limit on the half-life of ^{82}Se to the ground state of ^{82}Kr of T_{1/2}^{0ν}(^{82}Se)>4.6×10^{24}  yr (90% credible interval), corresponding to an effective Majorana neutrino mass m_{ββ}<(263-545)  meV. We also set the most stringent lower limits on the neutrinoless decays of ^{82}Se to the 0_{1}^{+}, 2_{1}^{+}, and 2_{2}^{+} excited states of ^{82}Kr, finding 1.8×10^{23}  yr, 3.0×10^{23}  yr, and 3.2×10^{23}  yr (90% credible interval) respectively.

Effective Majorana Neutrino Mass for ΔL = 2 Neutrino Oscillations

It is well known that the observations of neutrinoless double-beta decay prove the Majorana nature of the neutrino. However, with specific values of Majorana phases, the effective Majorana neutrino mass to be estimated from the observation of neutrinoless double-beta decay experiments is strongly suppressed if the neutrino mass pattern adheres to a normal ordering. In this case, double-beta decay might not be observed even though the neutrino is Majorana in nature. We show if neutrinos oscillate to antineutrinos in their propagation; then, the observation of this oscillation proves that neutrinos are Majorana and will provide a measurement of neutrino masses and Majorana phases.

Muon to Positron Conversion

Lepton-flavor violation (LFV) has been discovered in the neutrino sector by neutrino oscillation experiments. The minimal extension of the Standard Model (SM) to include neutrino masses allows LFV in the charged sector (CLFV) at the loop level, but at rates that are too small to be experimentally observed. Lepton-number violation (LNV) is explicitly forbidden even in the minimally extended SM, so the observation of an LNV process would be unambiguous evidence of physics beyond the SM. The search for the LNV and CLFV process μ−+N(A,Z)→e++N′(A,Z−2) (referred to as μ−→e+) complements 0νββ decay searches, and is sensitive to potential flavor effects in the neutrino mass-generation mechanism. A theoretical motivation for μ−→e+ is presented along with a review of the status of past μ−→e+ experiments and future prospects. Special attention is paid to an uncertain and potentially dominant background for these searches, namely, radiative muon capture (RMC). The RMC high energy photon spectrum is theoretically understudied and existing measurements insufficiently constrain this portion of the spectrum, leading to potentially significant impacts on current and future μ−→e+ work.

Ab Initio Calculation of the Contact Operator Contribution in the Standard Mechanism for Neutrinoless Double Beta Decay

Starting from chiral nuclear interactions, we evaluate the contribution of the leading-order contact transition operator to the nuclear matrix element (NME) of neutrinoless double-beta decay, assuming a light Majorana neutrino-exchange mechanism. The corresponding low-energy constant (LEC) is determined by fitting the transition amplitude of the nn→ppe^{-}e^{-} process to a recently proposed synthetic datum. We examine the dependence of the amplitude on similarity renormalization group scale and chiral expansion order of the nuclear interaction, finding that both dependences can be compensated to a large extent by readjusting the LEC. We evaluate the contribution of both the leading-order contact operator and standard long-range operator to the neutrinoless double-beta decays in the light nuclei ^{6,8}He and the candidate nucleus ^{48}Ca. Our results provide the first clear demonstration that the contact term enhances the NME in calculations with commonly used chiral two- plus three-nucleon interactions. In the case of ^{48}Ca, for example, the NME obtained with the EM(1.8/2.0) interaction is enhanced from 0.61 to 0.87(4), where the uncertainty is propagated from the synthetic datum.

Bolometric Double Beta Decay Experiments: Review and Prospects

This review aims to cover the history and recent developments on cryogenic bolometers for neutrinoless double beta decay (0ν2β) searches. A 0ν2β decay observation would confirm the total lepton charge non-conservation, which is related to a global U(1)LC symmetry. This discovery would also provide essential information on neutrino masses and nature, opening the door to new physics beyond the Standard Model. The bolometric technology shows good prospects for future ton-scale experiments that aim to fully investigate the inverted ordering region of neutrino masses. The big advantage of bolometers is the high energy resolution and the possibility of particle identification, as well as various methods of additional background rejection. The CUORE experiment has proved the feasibility of ton-scale cryogenic experiments, setting the most stringent limit on 130Te 0ν2β decay. Two CUPID demonstrators (CUPID-0 and CUPID-Mo) have set the most stringent limits on 82Se and 100Mo isotopes, respectively, with compatibly low exposures. Several experiments are developing new methods to improve the background in the region of interest with bolometric detectors. CUPID and AMoRE experiments aim to cover the inverted hierarchy region, using scintillating bolometers with hundreds of kg of 100Mo. We review all of these efforts here, with a focus on the different types of radioactive background and the measures put in place to mitigate them.

Scintillation in Low-Temperature Particle Detectors

Inorganic crystal scintillators play a crucial role in particle detection for various applications in fundamental physics and applied science. The use of such materials as scintillating bolometers, which operate at temperatures as low as 10 mK and detect both heat (phonon) and scintillation signals, significantly extends detectors performance compared to the conventional scintillation counters. In particular, such low-temperature devices offer a high energy resolution in a wide energy interval thanks to a phonon signal detection, while a simultaneous registration of scintillation emitted provides an efficient particle identification tool. This feature is of great importance for a background identification and rejection. Combined with a large variety of elements of interest, which can be embedded in crystal scintillators, scintillating bolometers represent powerful particle detectors for rare-event searches (e.g., rare alpha and beta decays, double-beta decay, dark matter particles, neutrino detection). Here, we review the features and results of low-temperature scintillation detection achieved over a 30-year history of developments of scintillating bolometers and their use in rare-event search experiments.

Status and Perspectives on Rare Decay Searches in Tellurium Isotopes

Neutrinoless double beta decay (0νββ) is a posited lepton number violating decay whose search is an increasingly active field in modern astroparticle physics. A discovery would imply neutrinos are Majorana particles and inform neutrino physics, cosmology and beyond-standard-model theories. Among the few nuclei where double beta decay (ββ) is allowed, tellurium isotopes stand for their high natural abundance and are currently employed in multiple experiments. The search for 0νββ will provide large exposure data sets in the coming years, paving the way for unprecedented sensitivities. We review the latest rare decay searches in tellurium isotopes and compare past results with theories and prospects from running experiments.

Coupled-Cluster Calculations of Neutrinoless Double-β Decay in ^{48}Ca

We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-β decay of ^{48}Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For ^{48}Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-β decay of ^{48}Ca with a quenching factor deduced from two-body currents in recent ab initio calculation of the Ikeda sum rule in ^{48}Ca [Gysbers et al., Nat. Phys. 15, 428 (2019)NPAHAX1745-247310.1038/s41567-019-0450-7].

Toward Complete Leading-Order Predictions for Neutrinoless Double β Decay

The amplitude for the neutrinoless double β (0νββ) decay of the two-neutron system nn→ppe^{-}e^{-} constitutes a key building block for nuclear-structure calculations of heavy nuclei employed in large-scale 0νββ searches. Assuming that the 0νββ process is mediated by a light-Majorana-neutrino exchange, a systematic analysis in chiral effective field theory shows that already at leading order a contact operator is required to ensure renormalizability. In this Letter, we develop a method to estimate the numerical value of its coefficient (in analogy to the Cottingham formula for electromagnetic contributions to hadron masses) and validate the result by reproducing the charge-independence-breaking contribution to the nucleon-nucleon scattering lengths. Our central result, while derived in dimensional regularization, is given in terms of the renormalized amplitude A_{ν}(|p|,|p^{'}|), matching to which will allow one to determine the contact-term contribution in regularization schemes employed in nuclear-structure calculations. Our results thus greatly reduce a crucial uncertainty in the interpretation of searches for 0νββ decay.

What Is Matter According to Particle Physics, and Why Try to Observe Its Creation in a Lab?

The standard model of elementary interactions has long qualified as a theory of matter, in which the postulated conservation laws (one baryonic and three leptonic) acquire theoretical meaning. However, recent observations of lepton number violations—neutrino oscillations—demonstrate its incompleteness. We discuss why these considerations suggest the correctness of Ettore Majorana’s ideas on the nature of neutrino mass and add further interest to the search for an ultra-rare nuclear process in which two particles of matter (electrons) are created, commonly called neutrinoless double beta decay. The approach of the discussion is mainly historical, and its character is introductory. Some technical considerations, which highlight the usefulness of Majorana’s representation of gamma matrices, are presented in the appendix.

Ab Initio Neutrinoless Double-Beta Decay Matrix Elements for ^{48}Ca, ^{76}Ge, and ^{82}Se

We calculate basis-space converged neutrinoless ββ-decay nuclear matrix elements for the lightest candidates: ^{48}Ca, ^{76}Ge, and ^{82}Se. Starting from initial two- and three-nucleon forces, we apply the ab initio in-medium similarity renormalization group to construct valence-space Hamiltonians and consistently transformed ββ-decay operators. We find that the tensor component is non-negligible in ^{76}Ge and ^{82}Se, and the resulting nuclear matrix elements are overall 25%-45% smaller than those obtained from the phenomenological shell model. While a final matrix element with uncertainties still requires substantial developments, this work nevertheless opens a path toward a true first-principles calculation of neutrinoless ββ decay in all nuclei relevant for ongoing large-scale searches.

Recent experimental activity on heavy-ion induced reactions within the NUMEN project

The possibility to use a special class of heavy-ion induced direct reactions, such as double charge exchange reactions, is discussed in view of their application to extract information that may be helpful to determinate the nuclear matrix elements entering in the expression of neutrinoless double beta decay half-life. The strategies adopted in the experimental campaigns performed at INFN - Laboratori Nazionali del Sud are briefly described, emphasizing the advantages of the multi-channel approach to nuclear reaction data analysis.

Present Status of Nuclear Shell-Model Calculations of 0νββ Decay Matrix Elements

Neutrinoless double beta (0νββ) decay searches are currently among the major foci of experimental physics. The observation of such a decay will have important implications in our understanding of the intrinsic nature of neutrinos and shed light on the limitations of the Standard Model. The rate of this process depends on both the unknown neutrino effective mass and the nuclear matrix element (M0ν) associated with the given 0νββ transition. The latter can only be provided by theoretical calculations, hence the need of accurate theoretical predictions of M0ν for the success of the experimental programs. This need drives the theoretical nuclear physics community to provide the most reliable calculations of M0ν. Among the various computational models adopted to solve the many-body nuclear problem, the shell model is widely considered as the basic framework of the microscopic description of the nucleus. Here, we review the most recent and advanced shell-model calculations of M0ν considering the light-neutrino-exchange channel for nuclei of experimental interest. We report the sensitivity of the theoretical calculations with respect to variations in the model spaces and the shell-model nuclear Hamiltonians.

Neutrino-Mass Sensitivity and Nuclear Matrix Element for Neutrinoless Double Beta Decay

Neutrinoless double beta decay (DBD) is a useful probe to study neutrino properties such as the Majorana nature, the absolute neutrino mass, the CP phase and the others, which are beyond the standard model. The nuclear matrix element (NME) for DBD is crucial to extract the neutrino properties from the experimental transition rate. The neutrino-mass sensitivity, i.e., the minimum neutrino-mass to be measured by the DBD experiment, is very sensitive to the DBD NME. Actually, the NME is one of the key elements for designing the DBD experiment. Theoretical evaluation for the DBD NME, however, is very hard. Recently experimental studies of charge-exchange nuclear and leptonic reactions have shown to be used to get single-β NMEs associated with the DBD NME. Critical discussions are made on the neutrino-mass sensitivity and the NME for the DBD neutrino-mass study and on the experimental studies of the single-β NMEs and nuclear structures associated with DBD NMEs.

Gamow-Teller Strength in ^{48}Ca and ^{78}Ni with the Charge-Exchange Subtracted Second Random-Phase Approximation

We develop a fully self-consistent subtracted second random-phase approximation for charge-exchange processes with Skyrme energy-density functionals. As a first application, we study Gamow-Teller excitations in the doubly magic nucleus ^{48}Ca, the lightest double-β emitter that could be used in an experiment, and in ^{78}Ni, the single-beta-decay rate of which is known. The amount of Gamow-Teller strength below 20 or 30 MeV is considerably smaller than in other energy-density-functional calculations and agrees better with experiment in ^{48}Ca, as does the beta-decay rate in ^{78}Ni. These important results, obtained without ad hoc quenching factors, are due to the presence of two-particle-two-hole configurations. Their density progressively increases with excitation energy, leading to a long high-energy tail in the spectrum, a fact that may have implications for the computation of nuclear matrix elements for neutrinoless double-β decay in the same framework.

Nuclear Response to Second-Order Isospin Probes in Connection to Double Beta Decay

One of the key ingredients needed to extract quantitative information on neutrino absolute mass scale from the possible measurement of the neutrinoless double-beta (0νββ) decay half-lives is the nuclear matrix element (NME) characterizing such transitions. NMEs are not physical observables and can only be deduced by theoretical calculations. However, since the atomic nuclei involved in the decay are many-body systems, only approximated values are available to date. In addition, the value of the coupling constants to be used for the weak interaction vertices is still an open question, which introduces a further indetermination in the calculations of NMEs. Several experimental approaches were developed in the years with the aim of providing useful information to further constrain the theory. Here we give an overview of the role of charge exchange reactions in this scenario, focusing on second-order processes, namely the double charge exchange (DCE) reactions.

Quenched g_{A} in Nuclei and Emergent Scale Symmetry in Baryonic Matter

The recent RIKEN experiment on the quenched g_{A} in the superallowed Gamow-Teller transition from ^{100}Sn indicates the role of scale anomaly encoded in the anomalous dimension β^{'} of the gluonic stress tensor Tr G_{μν}^{2}. This observation provides support to the notion of hidden scale symmetry emerging by strong nuclear correlations with an infrared (IR) fixed point realized-in the chiral limit-in the Nambu-Goldstone mode. We suggest there is an analogy in the way scale symmetry manifests in a nuclear medium to the continuity from the unitarity limit at low density (in light nuclei) to the dilaton limit at high density (in compact stars). In between the limits, say, at normal nuclear matter density, the symmetry is not visible, hence hidden.

Precise Half-Life Values for Two-Neutrino Double-β Decay: 2020 Review

All existing positive results on two-neutrino double beta decay and two-neutrino double electron capture in different nuclei have been analyzed. Weighted average and recommended half-life values for 48Ca, 76Ge, 82Se, 96Zr, 100Mo, 100Mo - 100Ru (01+), 116Cd, 128Te, 130Te, 136Xe, 150Nd, 150Nd - 150Sm (01+), 238U, 78Kr, 124Xe and 130Ba have been obtained. Given the measured half-life values, effective nuclear matrix elements for all these transitions were calculated.

Neutrinos: Majorana or Dirac?

Are neutrinos with definite masses Majorana or Dirac particles? This is one of the most fundamental problems of modern neutrino physics. The solution to this problem could be crucial for understanding the origin of small neutrino masses. We 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.

The NUMEN Heavy Ion Multidetector for a Complementary Approach to the Neutrinoless Double Beta Decay

Neutrinos are so far the most elusive known particles, and in the last decades many sophisticated experiments have been set up in order to clarify several questions about their intrinsic nature, in particular their masses, mass hierarchy, intrinsic nature of Majorana or Dirac particles. Evidence of the Neutrinoless Double-Beta Decay (NDBD) would prove that neutrinos are Majorana particles, thus improving the understanding of the universe itself. Besides the search for several large underground experiments for the direct experimental detection of NDBD, the NUMEN experiment proposes the investigation of a nuclear mechanism strongly linked to this decay: the Double Charge Exchange reactions (DCE). As such reactions share with the NDBD the same initial and final nuclear states, they could shed light on the determination of the Nuclear Matrix Elements (NMEs), which play a relevant role in the decay. The physics of DCE is described elsewhere in this issue, while the focus of this paper will be on the challenging experimental apparatus currently under construction in order to fulfil the requirements of the NUMEN experiment. The overall structure of the technological improvement to the cyclotron, along with the newly developed detection systems required for tracking and identifying the reaction products and their final excitation level are described.

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of ^{48}Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. We use our scheme to compute the matrix element that governs the neutrinoless double beta decay of ^{48}Ca to ^{48}Ti, and find it to have the value 0.61, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as ^{76}Ge, ^{130}Te, and ^{136}Xe.

Improved Limit on Neutrinoless Double-Beta Decay in ^{130} Te with CUORE

We report new results from the search for neutrinoless double-beta decay in ^{130} Te with the CUORE detector. This search benefits from a fourfold increase in exposure, lower trigger thresholds, and analysis improvements relative to our previous results. We observe a background of (1.38±0.07)×10^{-2}  counts/(keV kg yr)) in the 0νββ decay region of interest and, with a total exposure of 372.5 kg yr, we attain a median exclusion sensitivity of 1.7×10^{25}  yr. We find no evidence for 0νββ decay and set a 90% credibility interval Bayesian lower limit of 3.2×10^{25}  yr on the ^{130} Te half-life for this process. In the hypothesis that 0νββ decay is mediated by light Majorana neutrinos, this results in an upper limit on the effective Majorana mass of 75-350 meV, depending on the nuclear matrix elements used.

Relevance of the Nuclear Structure of the Stable Ge Isotopes to the Neutrino-less Double-Beta Decay of 76Ge

Gamma-ray detection following the inelastic neutron scattering reaction on isotopically enriched material was used to study the nuclear structure of 74Ge. From these measurements, low-lying, low-spin excited states were characterized, new states and their decays were identified, level lifetimes were measured with the Doppler-shift attenuation method (DSAM), multipole mixing ratios were established, and transition probabilities were determined. New structural features in 74Ge were identified, and the reanalysis of older 76Ge data led to the placement of the 2+ member of the intruder band. In addition, a number of previously placed states in 74Ge were shown not to exist. A procedure for future work, which will lead to meaningful data for constraining calculations of the neutrinoless double-beta decay matrix element, is suggested.

Neutrinoless Double-Beta Decay and Realistic Shell Model

We report on the calculation of the neutrinoless double-ß decay nuclear matrix element for 76Ge within the framework of the realistic shell model. The effective shell-model Hamiltonian and the two-body transition operator describing the decay are derived by way of many-body perturbation theory. Particular attention is focused on the role played by the so-called Pauli blocking effect in the derivation of the effective operator.

Time-Reversed Particle-Vibration Loops and Nuclear Gamow-Teller Response

The nuclear response theory for charge-exchange modes in the relativistic particle-vibration coupling approach is extended to include for the first time particle-vibration coupling effects in the ground state of the parent nucleus. In a framework based on the effective meson-nucleon Lagrangian, we investigate the role of such complex ground-state correlations in the description of Gamow-Teller transitions in ^{90}Zr in both (p, n) and (n, p) channels. The particle-vibration coupling effects are calculated without introducing new parameters. We find that this new correlation mechanism is fully responsible for the appearance of the strength in the (n, p) branch. Comparison of our results to the available experimental data shows a very good agreement up to excitation energies beyond the giant resonance region when taking into account a phenomenological admixture of the isovector spin monopole transitions. The parent-daughter binding-energy differences are also greatly improved by the inclusion of the new correlations.

Search for Neutrinoless Double-β Decay with the Complete EXO-200 Dataset

A search for neutrinoless double-β decay (0νββ) in ^{136}Xe is performed with the full EXO-200 dataset using a deep neural network to discriminate between 0νββ and background events. Relative to previous analyses, the signal detection efficiency has been raised from 80.8% to 96.4±3.0%, and the energy resolution of the detector at the Q value of ^{136}Xe 0νββ has been improved from σ/E=1.23% to 1.15±0.02% with the upgraded detector. Accounting for the new data, the median 90% confidence level 0νββ half-life sensitivity for this analysis is 5.0×10^{25}  yr with a total ^{136}Xe exposure of 234.1 kg yr. No statistically significant evidence for 0νββ is observed, leading to a lower limit on the 0νββ half-life of 3.5×10^{25}  yr at the 90% confidence level.

Probing Majorana neutrinos with double-β decay

A discovery that neutrinos are Majorana fermions would have profound implications for particle physics and cosmology. The Majorana character of neutrinos would make possible the neutrinoless double-β (0νββ) decay, a matter-creating process without the balancing emission of antimatter. The GERDA Collaboration searches for the 0νββ decay of ⁷⁶Ge by operating bare germanium detectors in an active liquid argon shield. With a total exposure of 82.4 kg⋅year, we observe no signal and derive a lower half-life limit of T _1/2 > 0.9 × 10²⁶ years (90% C.L.). Our T _1/2 sensitivity, assuming no signal, is 1.1 × 10²⁶ years. Combining the latter with those from other 0νββ decay searches yields a sensitivity to the effective Majorana neutrino mass of 0.07 to 0.16 electron volts.

Evidence for Rigid Triaxial Deformation in ^{76}Ge from a Model-Independent Analysis

An extensive, model-independent analysis of the nature of triaxial deformation in ^{76}Ge, a candidate for neutrinoless double-beta (0νββ) decay, was carried out following multistep Coulomb excitation. Shape parameters deduced on the basis of a rotational-invariant sum-rule analysis provided considerable insight into the underlying collectivity of the ground-state and γ bands. Both sequences were determined to be characterized by the same β and γ deformation parameter values. In addition, compelling evidence for low-spin, rigid triaxial deformation in ^{76}Ge was obtained for the first time from the analysis of the statistical fluctuations of the quadrupole asymmetry deduced from the measured E2 matrix elements. These newly determined shape parameters are important input and constraints for calculations aimed at providing, with suitable accuracy, the nuclear matrix elements relevant to 0νββ.

Final Result of CUPID-0 Phase-I in the Search for the ^{82}Se Neutrinoless Double-β Decay

CUPID-0 is the first pilot experiment of CUPID, a next-generation project for the measurement of neutrinoless double beta decay (0νDBD) with scintillating bolometers. The detector, consisting of 24 enriched and 2 natural ZnSe crystals, has been taking data at Laboratori Nazionali del Gran Sasso from June 2017 to December 2018, collecting a ^{82}Se exposure of 5.29  kg×yr. In this Letter we present the phase-I results in the search for 0νDBD. We demonstrate that the technology implemented by CUPID-0 allows us to reach the lowest background for calorimetric experiments: (3.5_{-0.9}^{+1.0})×10^{-3}  counts/(keV kg yr). Monitoring 3.88×10^{25} ^{82}Se nuclei×yr we reach a 90% credible interval median sensitivity of T_{1/2}^{0ν}>5.0×10^{24}  yr and set the most stringent limit on the half-life of ^{82}Se 0νDBD: T_{1/2}^{0ν}>3.5×10^{24}  yr (90% credible interval), corresponding to m_{ββ}<(311-638)  meV depending on the nuclear matrix element calculations.

Precision Analysis of the ^{136}Xe Two-Neutrino ββ Spectrum in KamLAND-Zen and Its Impact on the Quenching of Nuclear Matrix Elements

We present a precision analysis of the ^{136}Xe two-neutrino ββ electron spectrum above 0.8 MeV, based on high-statistics data obtained with the KamLAND-Zen experiment. An improved formalism for the two-neutrino ββ rate allows us to measure the ratio of the leading and subleading 2νββ nuclear matrix elements (NMEs), ξ_{31}^{2ν}=-0.26_{-0.25}^{+0.31}. Theoretical predictions from the nuclear shell model and the majority of the quasiparticle random-phase approximation (QRPA) calculations are consistent with the experimental limit. However, part of the ξ_{31}^{2ν} range allowed by the QRPA is excluded by the present measurement at the 90% confidence level. Our analysis reveals that predicted ξ_{31}^{2ν} values are sensitive to the quenching of NMEs and the competing contributions from low- and high-energy states in the intermediate nucleus. Because these aspects are also at play in neutrinoless ββ decay, ξ_{31}^{2ν} provides new insights toward reliable neutrinoless ββ NMEs.

Light-Neutrino Exchange and Long-Distance Contributions to 0ν2β Decays: An Exploratory Study on ππ→ee

We present an exploratory lattice QCD calculation of the neutrinoless double beta decay ππ→ee. Under the mechanism of light-neutrino exchange, the decay amplitude involves significant long-distance contributions. The calculation reported here, with pion masses m_{π}=420 and 140 MeV, demonstrates that the decay amplitude can be computed from first principles using lattice methods. At unphysical and physical pion masses, we obtain that amplitudes are 24% and 9% smaller than the predication from leading order chiral perturbation theory. Our findings provide the lattice QCD inputs and constraints for effective field theory. A follow-on calculation with fully controlled systematic errors will be possible with adequate computational resources.

Recent results on heavy-ion induced reactions of interest for neutrinoless double beta decay at INFN-LNS

The possibility to use a special class of heavy-ion induced direct reactions, such as double charge exchange reactions, is discussed in view of their application to extract information that may be helpful to determinate the nuclear matrix elements entering in the expression of neutrinoless double beta decay half-life. The methodology of the experimental campaign presently running at INFN - Laboratori Nazionali del Sud is reported andthe experimental challenges characterizing such activity are described.

Contributed Review: The saga of neutrinoless double beta decay search with TeO2 thermal detectors

Neutrinoless double beta decay (0νββ) is a direct probe of physics beyond the standard model. Its discovery would demonstrate that the lepton number is not a symmetry of nature and would provide us with unique information on the nature and mass of the neutrinos. Among the experimental techniques employed in the investigation of this rare process, thermal detectors fulfill the requirements for a powerful search, showing an excellent energy resolution and the possibility of scaling to very large masses. In this work, we review the long chain of bolometric experiments based on TeO₂ crystals that were and continue to be carried out at the Laboratori Nazionali del Gran Sasso (Italy), searching for 0νββ of ¹³⁰Te. We illustrate the progress and improvements achieved in almost thirty years of measurements and compare the various performance and results. We describe the several steps that led to the CUORE detector, the latest of this series and presently in data taking, and we highlight the challenges that a next bolometric experiment will face in order to further improve the sensitivity, especially concerning the background abatement. Finally, we emphasize the advantages of ¹³⁰Te in the search for 0νββ with a further future experiment.

Two neutrino double-beta decay and effective axial-vector coupling constant

The theoretical and experimental study of the two-neutrino double-beta decay (2νββ-decay) is of crucial importance for reliable calculation of matrix elements governing the neutrinoless double-beta decay (0νββ-decay). That will allow to determine masses of neutrinos once the 0νββ-decay, which occurs if the neutrino is a massive Majorana particle and the total lepton number is not onserved, will be observed. Experiments studying the 2νββ-decay are currently approaching a qualitatively new level, where high-precision measurements are performed not only for half-lives but for all other observables of the process. In this context an improved formula for the 2νββ-decay is presented. Further, a novel approach for determining the effective axial-vector coupling constant is proposed.

The nuclear matrix elements of 0νββ decay and the NUMEN project at INFN-LNS

The goal of NUMEN project is to access experimentally driven information on Nuclear Matrix Elements (NME) involved in the neutrinoless double beta decay (0νββ), by high-accuracy measurements of the cross sections of Heavy Ion (HI) induced Double Charge Exchange (DCE) reactions. The knowledge of the nuclear matrix elements is crucial to infer the neutrino average masses from the possible measurement of the half-life of 00νββ decay and to compare experiments on different isotopes. In particular, the (18O, 18Ne) and (20Ne, 20O) reactions are performed as tools for β+β+ and β-β- decays, respectively. The experiments are performed at INFN - Laboratory Nazionali del Sud (LNS) in Catania using the Superconducting Cyclotron (CS) to accelerate the beams and the MAGNEX magnetic spectrometer to detect the reaction products. The measured cross sections are very low, limiting the present exploration to few selected isotopes of interest in the context of typically low-yield experimental runs. In order to make feasible a systematic study of all the candidate nuclei, a major upgrade of the LNS facility is foreseen to increase the experimental yield of about two orders of magnitude. To this purpose, frontier technologies are going to be developed for both the accelerator and the detection systems. In parallel, advanced theoretical models will be developed to extract the nuclear structure information from the measured cross sections.

Heavy Physics Contributions to Neutrinoless Double Beta Decay from QCD

Observation of neutrinoless double beta decay, a lepton number violating process that has been proposed to clarify the nature of neutrino masses, has spawned an enormous world-wide experimental effort. Relating nuclear decay rates to high-energy, beyond the standard model (BSM) physics requires detailed knowledge of nonperturbative QCD effects. Using lattice QCD, we compute the necessary matrix elements of short-range operators, which arise due to heavy BSM mediators, that contribute to this decay via the leading order π^{-}→π^{+} exchange diagrams. Utilizing our result and taking advantage of effective field theory methods will allow for model-independent calculations of the relevant two-nucleon decay, which may then be used as input for nuclear many-body calculations of the relevant experimental decays. Contributions from short-range operators may prove to be equally important to, or even more important than, those from long-range Majorana neutrino exchange.