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

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

Opportunities for fundamental physics research with radioactive molecules

Molecules containing short-lived, radioactive nuclei are uniquely positioned to enable a wide range of scientific discoveries in the areas of fundamental symmetries, astrophysics, nuclear structure, and chemistry. Recent advances in the ability to create, cool, and control complex molecules down to the quantum level, along with recent and upcoming advances in radioactive species production at several facilities around the world, create a compelling opportunity to coordinate and combine these efforts to bring precision measurement and control to molecules containing extreme nuclei. In this manuscript, we review the scientific case for studying radioactive molecules, discuss recent atomic, molecular, nuclear, astrophysical, and chemical advances which provide the foundation for their study, describe the facilities where these species are and will be produced, and provide an outlook for the future of this nascent field.

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.

Ab Initio Structure Factors for Spin-Dependent Dark Matter Direct Detection

We present converged ab initio calculations of structure factors for elastic spin-dependent WIMP scattering off all nuclei used in dark matter direct-detection searches: ^{19}F, ^{23}Na, ^{27}Al, ^{29}Si, ^{73}Ge, ^{127}I, and ^{129,131}Xe. From a set of established two- and three-nucleon interactions derived within chiral effective field theory, we construct consistent WIMP-nucleon currents at the one-body level, including effects from axial-vector two-body currents. We then apply the in-medium similarity renormalization group to construct effective valence-space Hamiltonians and consistently transformed operators of nuclear responses. Combining the recent advances of natural orbitals with three-nucleon forces expressed in large spaces, we obtain basis-space converged structure factors even in heavy nuclei. Generally results are consistent with previous calculations but large uncertainties in ^{127}I highlight the need for further study.

Interplay of nuclear physics, effective field theories, phenomenology, and lattice QCD in neutrino physics

Experiments in neutrino physics cover a wide range—from deep inelastic scattering, over long base-line oscillation experiments and low-energy coherent neutrino–nucleus scattering (CEνNS), to searches for neutrinoless double β decay (0νββ)—yet in all cases a key aspect in interpreting the results concerns understanding neutrino–nucleus interactions. If the neutrino energy is sufficiently low, the required matrix elements can be constrained in a systematic way by the interplay of effective field theories, phenomenology, and lattice QCD. In these proceedings, we illustrate this strategy focusing on the CEνNS and 0νββ processes.

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.

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.