Limits on neutrino masses from neutrinoless double-β decay

Enriched Crystal Scintillators for 2β Experiments

The investigation of 2β decay is an important issue in modern physics, allowing the test of the Standard Model of elementary particles and the study of the nature and properties of neutrinos. The crystal scintillators, especially made of isotopically-enriched materials, are powerful detectors for 2β decay experiments thanks to the high radiopurity level and the possibility to realize the calorimetric “source = detector” approach with a high detection efficiency. For the moment, the 2ν2β processes have been observed at the level of 1019–1024 years with enriched crystals; the sensitivity to the 0ν mode have reached the level of 1024–1026 years in some decay channels for different nuclides allowing one to calculate the upper limits on the effective mass of the Majorana neutrino at the level of 0.1–0.6 eV. The paper is intended to be a review on the latest results to investigate 2β processes with crystal scintillators enriched in 48Ca, 106Cd, and 116Cd.

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.

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.

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.

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νββ.

Implementation of Feldman-Cousins corrections and oscillation calculations in the HPC environment for the NOvA Experiment

Analysis of neutrino oscillation data involves a combination of complex fitting procedures and statistical correction techniques that are used to determine the full three-flavor PMNS parameters and constraint contours. These techniques rely on computationally intensive “multi-universe” stochastic modeling. The process of calculating these contours and corrections can dominate final stages of the data analysis and become a bottleneck for examining the ef-fect of systematic variations on the final results. As part of the DOE SciDAC-4 sponsored research program, we present a new implementation of a neutrino oscillation fitting and framework to carry out calculations of Feldman-Cousins corrections. This implementation is specifically designed and optimized to operate on modern High-Performance Computing (HPC) facilities. We present the performance of the system in calculating allowed regions for the NOvA experiment based on 8.9 x 1020 and 6.9 x 1020 protons-on-target (POT) neutrino and antineutrino datasets, and compare the performance of this new implementation run at the NERSC supercomputing facility with that from methods used previously by the NOvA collaboration running on grid computing facilities.

Lepton Number Violation in Higgs Decay at LHC

We show that within the left-right symmetric model, lepton number violating decays of the Higgs boson can be discovered at the LHC. The process is due to the mixing of the Higgs boson with the triplet that breaks parity. As a result, the Higgs boson can act as a gateway to the origin of the heavy Majorana neutrino mass. To assess the LHC reach, a detailed collider study of the same-sign dileptons plus jets channel is provided. This process is complementary to the existing nuclear and collider searches for lepton number violation and can probe the scale of parity restoration even beyond other direct searches.

Nuclear structure aspects of neutrinoless double-β decay

We decompose the neutrinoless double-β decay matrix elements into sums of products over the intermediate nucleus with two less nucleons. We find that the sum is dominated by the J(π)=0(+) ground state of this intermediate nucleus for both the light and heavy neutrino decay processes. This provides a new theoretical tool for comparing and improving nuclear structure models. It also provides the connection to two-nucleon transfer experiments.

Constraint on 0νββ matrix elements from a novel decay channel of the scissors mode: the case of 154Gd

The nucleus (154)Gd is located in a region of the nuclear chart where rapid changes of nuclear deformation occur as a function of particle number. It was investigated using a combination of γ-ray scattering experiments and a γγ-coincidence study following electron capture decay of (154)Tb(m). A novel decay channel from the scissors mode to the first excited 0(+) state was observed. Its transition strength was determined to B(M1;1(sc)(+)→0(2)(+))=0.031(4)μ(N)(2). The properties of the scissors mode of (154)Gd imply a much larger matrix element than previously thought for the neutrinoless double-β decay to the 0(2)(+) state in such a shape-transitional region. Theory indicates an even larger effect for (150)Nd.

Shell-model analysis of the 136Xe double beta decay nuclear matrix elements

Neutrinoless double beta decay, if observed, could distinguish whether the neutrino is a Dirac or a Majorana particle, and it could be used to determine the absolute scale of the neutrino masses. 136Xe is one of the most promising candidates for observing this rare event. However, until recently there were no positive results for the allowed and less rare two-neutrino double beta decay mode. The small nuclear matrix element associated with the long half-life represents a challenge for nuclear structure models used for its calculation. We report a new shell-model analysis of the two-neutrino double beta decay of 136Xe, which takes into account all relevant nuclear orbitals necessary to fully describe the associated Gamow-Teller strength. We further use the new model to analyze the main contributions to the neutrinoless double beta decay matrix element, and show that they are also diminished.

Limit on neutrinoless ββ decay of 136Xe from the first phase of KamLAND-Zen and comparison with the positive claim in 76Ge

We present results from the first phase of the KamLAND-Zen double-beta decay experiment, corresponding to an exposure of 89.5 kg yr of (136)Xe. We obtain a lower limit for the neutrinoless double-beta decay half-life of T(1/2)(0ν)>1.9×10(25) yr at 90% C.L. The combined results from KamLAND-Zen and EXO-200 give T(1/2)(0ν)>3.4×10(25) yr at 90% C.L., which corresponds to a Majorana neutrino mass limit of <m(ββ)> <(120-250) meV based on a representative range of available matrix element calculations. Using those calculations, this result excludes the Majorana neutrino mass range expected from the neutrinoless double-beta decay detection claim in (76)Ge, reported by a part of the Heidelberg-Moscow Collaboration, at more than 97.5% C.L.