Search

We present the results of a National Science Foundation Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a roadmap for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The roadmap includes specific projects in theoretical and computational physics as well as the use of Bayesian methods to quantify both intra- and inter-model uncertainties. The goal of this ambitious program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play.

Self-conjugate nuclei in the A≈70–80 region have attracted a great deal of attention due to phenomena such as shape coexistence and increasing collectivity along the N=Z line. We investigate the structure of nuclei in this region through lifetime measurements using the GRETINA array. The first implementation of the Differential Recoil Distance Doppler Shift technique with fast radioactive beams is demonstrated and verified through a measurement of the well-known B(E2;21+→01+) transition strength in 74Kr. The method is then applied to determine the B(E2;21+→01+) transition strength in 74Rb, the heaviest odd–odd N=Z nucleus for which this quantity has been determined. This result and extended systematics along N=Z suggest the dominance of the isoscalar part of the quadrupole transition strengths in self-conjugate nuclei, as well as the possible presence of shape coexistence in 74Rb.

Self-conjugate nuclei in the A≈70–80 region have attracted a great deal of attention due to phenomena such as shape coexistence and increasing collectivity along the N=Z line. We investigate the structure of nuclei in this region through lifetime measurements using the GRETINA array. The first implementation of the Differential Recoil Distance Doppler Shift technique with fast radioactive beams is demonstrated and verified through a measurement of the well-known B(E2;21+→01+) transition strength in 74Kr. The method is then applied to determine the B(E2;21+→01+) transition strength in 74Rb, the heaviest odd–odd N=Z nucleus for which this quantity has been determined. This result and extended systematics along N=Z suggest the dominance of the isoscalar part of the quadrupole transition strengths in self-conjugate nuclei, as well as the possible presence of shape coexistence in 74Rb.

The "island of inversion" at N≈20 for the neon, sodium, and magnesium isotopes has long been an area of interest both experimentally and theoretically due to the subtle competition between 0p-0h and np-nh configurations leading to deformed shapes. However, the presence of rotational band structures, which are fingerprints of deformed shapes, have only recently been observed in this region. In this work, we report on a measurement of the low-lying level structure of Mg33 populated by a two-stage projectile fragmentation reaction and studied with the Gamma Ray Energy Tracking In-Beam Nuclear Array (GRETINA). The experimental level energies, ground-state magnetic moment, intrinsic quadrupole moment, and γ-ray intensities show good agreement with the strong-coupling limit of a rotational model.

The calcium isotopes have emerged as a critical testing ground for new microscopically derived shell-model interactions, and a great deal of experimental and theoretical focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with 1f7/2 neutron hole states in Ca47,49 following one-neutron knockout reactions from Ca48,50. The observed reduction of strength populating the 7/21- state in Ca49, as compared to Ca47, is inconsistent with shell-model calculations using both phenomenological interactions such as GXPF1, and interactions derived from microscopically based two- and three-nucleon forces. The result suggests a fragmentation of the l=3 strength to higher-lying states as suggested by the microscopic calculations, but the observed magnitude of the reduction is not reproduced in any shell-model description.

A segmented Si-telescope and HPGe array, STARS-LIBERACE, was used to study the 156Gd(p, tγ)154Gd direct reaction by particle- γ coincidence spectroscopy. New cross sections with a 25MeV proton beam are reported and compared to previous (p,t) and (t,p) studies. Furthermore, additional evidence for coexisting Kπ=01+,21+ and 02+, 22+ configurations at N = 90 is presented. Direct and indirect population patterns of the low-lying states are also explored. Review of the new and existing evidence favors an interpretation based on a configuration-dependent pairing interaction. The weakening of monopole pairing strength and an increase in quadrupole pairing strength could bring 2p-2h 0 + states below 2 Δ. This may account for a large number of the low-lying 0 + states observed in two-nucleon transfer reactions. A hypothesis for the origin of the 02+ and 03+ states is provided.

The γ-ray tracking array GRETINA was coupled to the S800 magnetic spectrometer for spectroscopy with fast beams of rare isotopes at the National Superconducting Cyclotron Laboratory on the campus of Michigan State University. We describe the technical details of this powerful setup and report on GRETINA's performance achieved with source and in-beam measurements. The γ-ray multiplicity encountered in experiments with fast beams is usually low, allowing for a simplified and efficient treatment of the data in the γ-ray analysis in terms of Doppler reconstruction and spectral quality. The results reported in this work were obtained from GRETINA consisting of 8 detector modules hosting four high-purity germanium crystals each. Currently, GRETINA consists of 10 detector modules.

Previous experiments observed a 4^{+} state in the N=28 nucleus ^{44}S and suggested that this state may exhibit a hindered E2-decay rate, inconsistent with being a member of the collective ground state band. We populate this state via two-proton knockout from a beam of exotic ^{46}Ar projectiles and measure its lifetime using the recoil distance method with the GRETINA γ-ray spectrometer. The result, 76(14)_{stat}(20)_{syst}  ps, implies a hindered transition of B(E2;4^{+}→2_{1}^{+})=0.61(19) single-particle or Weisskopf units strength and supports the interpretation of the 4^{+} state as a K=4 isomer, the first example of a high-K isomer in a nucleus of such low mass.

Previous experiments observed a 4^{+} state in the N=28 nucleus ^{44}S and suggested that this state may exhibit a hindered E2-decay rate, inconsistent with being a member of the collective ground state band. We populate this state via two-proton knockout from a beam of exotic ^{46}Ar projectiles and measure its lifetime using the recoil distance method with the GRETINA γ-ray spectrometer. The result, 76(14)_{stat}(20)_{syst}  ps, implies a hindered transition of B(E2;4^{+}→2_{1}^{+})=0.61(19) single-particle or Weisskopf units strength and supports the interpretation of the 4^{+} state as a K=4 isomer, the first example of a high-K isomer in a nucleus of such low mass.

The distribution of single-particle strength in Ni67,69 was characterized with one-neutron knockout reactions from intermediate-energy Ni68,70 secondary beams, selectively populating neutron-hole configurations at N=39 and 41, respectively. The spectroscopic strengths deduced from the measured partial cross sections to the individual final states, as tagged by their γ-ray decays, are used to identify and quantify neutron configurations in the wave functions. While Ni69 compares well with shell-model predictions, the results for Ni67 challenge the validity of current effective shell-model Hamiltonians by revealing discrepancies that cannot be explained so far. These results suggest that our understanding of the low-lying states in the neutron-rich, semimagic Ni isotopes may be incomplete and requires further investigation on both the experimental and theoretical sides.

The distribution of single-particle strength in Ni67,69 was characterized with one-neutron knockout reactions from intermediate-energy Ni68,70 secondary beams, selectively populating neutron-hole configurations at N=39 and 41, respectively. The spectroscopic strengths deduced from the measured partial cross sections to the individual final states, as tagged by their γ-ray decays, are used to identify and quantify neutron configurations in the wave functions. While Ni69 compares well with shell-model predictions, the results for Ni67 challenge the validity of current effective shell-model Hamiltonians by revealing discrepancies that cannot be explained so far. These results suggest that our understanding of the low-lying states in the neutron-rich, semimagic Ni isotopes may be incomplete and requires further investigation on both the experimental and theoretical sides.

The calcium isotopes have emerged as an important testing ground for newmicroscopically derived shell-model interactions, and a great deal of focus hasbeen directed toward this region. We investigate the relative spectroscopicstrengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Cafollowing one-neutron knockout reactions from $^{48,50}$Ca. The observedreduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca,as compared to $^{47}$Ca, is consistent with the description given byshell-model calculations based on two- and three-nucleon forces in the neutron$pf$ model space, implying a fragmentation of the $l$=3 strength tohigher-lying states. The experimental result is inconsistent with both theGXPF1 interaction routinely used in this region of the nuclear chart and withmicroscopic calculations in an extended model space including the $u1g_{9/2}$orbital.

The calcium isotopes have emerged as an important testing ground for newmicroscopically derived shell-model interactions, and a great deal of focus hasbeen directed toward this region. We investigate the relative spectroscopicstrengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Cafollowing one-neutron knockout reactions from $^{48,50}$Ca. The observedreduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca,as compared to $^{47}$Ca, is consistent with the description given byshell-model calculations based on two- and three-nucleon forces in the neutron$pf$ model space, implying a fragmentation of the $l$=3 strength tohigher-lying states. The experimental result is inconsistent with both theGXPF1 interaction routinely used in this region of the nuclear chart and withmicroscopic calculations in an extended model space including the $\nu1g_{9/2}$orbital.

The quadrupole transition rate for the 41+→21+ transition of Ni58 was determined from an application of the recoil distance method with the Gamma-Ray Energy Tracking In-beam Nuclear Array (GRETINA). The present result of the B(E2;41+→21+) was found to be 50-6+11e2fm4, which is about three times smaller than the literature value, indicating substantially less collectivity than previously believed. Shell model calculations performed with the GXPF1A effective interaction agree with the present data and the validity of the standard effective charges in understanding collectivity in the nickel isotopes is discussed.

The quadrupole transition rate for the 41+→21+ transition of Ni58 was determined from an application of the recoil distance method with the Gamma-Ray Energy Tracking In-beam Nuclear Array (GRETINA). The present result of the B(E2;41+→21+) was found to be 50-6+11e2fm4, which is about three times smaller than the literature value, indicating substantially less collectivity than previously believed. Shell model calculations performed with the GXPF1A effective interaction agree with the present data and the validity of the standard effective charges in understanding collectivity in the nickel isotopes is discussed.

There is significant evidence supporting the existence of deformed ground states within the neutron-rich N≈20 neon, sodium, and magnesium isotopes that make up what is commonly called the "island of inversion." However, the rotational band structures, which are a characteristic fingerprint of a rigid nonspherical shape, have yet to be observed. In this work, we report on a measurement and analysis of the yrast (lowest lying) rotational band in Mg32 up to spin I=6+ produced in a two-step projectile fragmentation reaction and observed using the state-of-the-art γ-ray tracking detector array, GRETINA (γ-ray energy tracking in-beam nuclear array). Large-scale shell-model calculations using the SDPF-U-MIX effective interaction show excellent agreement with the new data. Moreover, a theoretical analysis of the spectrum of rotational states as a function of the pairing gap, together with cranked-shell-model calculations, provides intriguing evidence for a reduction in pairing correlations with increased angular momentum, also in line with the shell-model results.

We studied transition rates for the lowest 1/2+ and 5/2+ excited states of C17 through lifetime measurements with the GRETINA array using the recoil-distance method. The present measurements provide a model-independent determination of transition strengths giving the values of B(M1;1/2+→3/2g.s.+)=1.04-0.12+0.03×10-2μN2 and B(M1;5/2+→3/2g.s.+)=7.12-0.96+1.27×10-2μN2. The quenched M1 transition strength for the 1/2+→3/2g.s.+ transition, with respect to the 5/2+→3/2g.s.+ transition, has been confirmed with greater precision. The current data are compared to importance-truncated no-core shell model calculations addressing effects due to continuum and three-body forces.

We studied transition rates for the lowest 1/2+ and 5/2+ excited states of C17 through lifetime measurements with the GRETINA array using the recoil-distance method. The present measurements provide a model-independent determination of transition strengths giving the values of B(M1;1/2+→3/2g.s.+)=1.04-0.12+0.03×10-2μN2 and B(M1;5/2+→3/2g.s.+)=7.12-0.96+1.27×10-2μN2. The quenched M1 transition strength for the 1/2+→3/2g.s.+ transition, with respect to the 5/2+→3/2g.s.+ transition, has been confirmed with greater precision. The current data are compared to importance-truncated no-core shell model calculations addressing effects due to continuum and three-body forces.

Background: Triplet energy differences (TED) can be studied to yield information on isospin-non-conserving interactions in nuclei. Purpose: The systematic behavior of triplet energy differences (TED) of T=1, Jπ=2+ states is examined. The A=62 isobar is identified as having a TED value that deviates significantly from an otherwise very consistent trend. This deviation can be attributed to the tentative assignments of the pertinent states in Ga62 and Ge62. Methods: An in-beam γ-ray spectroscopy experiment was performed to identify excited states in Ga62 using Gamma-Ray Energy Tracking In-Beam Nuclear Array with the S800 spectrometer at NSCL using a two-nucleon knockout approach. Cross-section calculations for the knockout process and shell-model calculations have been performed to interpret the population and decay properties observed. Results: Using the systematics as a guide, a candidate for the transition from the T=1, 2+ state is identified. However, previous work has identified similar states with different Jπ assignments. Cross-section calculations indicate that the relevant T=1, 2+ state should be one of the states directly populated in this reaction. Conclusions: As spins and parities were not measurable, it is concluded that an unambiguous identification of the first T=1, 2+ state is required to reconcile our understanding of TED systematics.

An approach is presented to experimentally constrain previously unreachable (p, γ) reaction rates on nuclei far from stability in the astrophysical rp process. Energies of all critical resonances in the (57)Cu(p,γ)(58)Zn reaction are deduced by populating states in (58)Zn with a (d, n) reaction in inverse kinematics at 75 MeV/u, and detecting γ-ray-recoil coincidences with the state-of-the-art γ-ray tracking array GRETINA and the S800 spectrograph at the National Superconducting Cyclotron Laboratory. The results reduce the uncertainty in the (57)Cu(p,γ) reaction rate by several orders of magnitude. The effective lifetime of (56)Ni, an important waiting point in the rp process in x-ray bursts, can now be determined entirely from experimentally constrained reaction rates.

An approach is presented to experimentally constrain previously unreachable (p, γ) reaction rates on nuclei far from stability in the astrophysical rp process. Energies of all critical resonances in the (57)Cu(p,γ)(58)Zn reaction are deduced by populating states in (58)Zn with a (d, n) reaction in inverse kinematics at 75 MeV/u, and detecting γ-ray-recoil coincidences with the state-of-the-art γ-ray tracking array GRETINA and the S800 spectrograph at the National Superconducting Cyclotron Laboratory. The results reduce the uncertainty in the (57)Cu(p,γ) reaction rate by several orders of magnitude. The effective lifetime of (56)Ni, an important waiting point in the rp process in x-ray bursts, can now be determined entirely from experimentally constrained reaction rates.

Excited states in the neutron-rich N = 38, 36 nuclei (60)Ti and (58)Ti were populated in nucleon-removal reactions from (61)V projectiles at 90 MeV/nucleon. The γ-ray transitions from such states in these Ti isotopes were detected with the advanced γ-ray tracking array GRETINA and were corrected event by event for large Doppler shifts (v/c ∼ 0.4) using the γ-ray interaction points deduced from online signal decomposition. The new data indicate that a steep decrease in quadrupole collectivity occurs when moving from neutron-rich N = 36, 38 Fe and Cr toward the Ti and Ca isotones. In fact, (58,60)Ti provide some of the most neutron-rich benchmarks accessible today for calculations attempting to determine the structure of the potentially doubly magic nucleus (60)Ca.

Excited states in the neutron-rich N = 38, 36 nuclei (60)Ti and (58)Ti were populated in nucleon-removal reactions from (61)V projectiles at 90 MeV/nucleon. The γ-ray transitions from such states in these Ti isotopes were detected with the advanced γ-ray tracking array GRETINA and were corrected event by event for large Doppler shifts (v/c ∼ 0.4) using the γ-ray interaction points deduced from online signal decomposition. The new data indicate that a steep decrease in quadrupole collectivity occurs when moving from neutron-rich N = 36, 38 Fe and Cr toward the Ti and Ca isotones. In fact, (58,60)Ti provide some of the most neutron-rich benchmarks accessible today for calculations attempting to determine the structure of the potentially doubly magic nucleus (60)Ca.

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.

In spite of the high-density and strongly correlated nature of the atomic nucleus, experimental and theoretical evidence suggests that around particular 'magic' numbers of nucleons, nuclear properties are governed by a single unpaired nucleon 1,2 . A microscopic understanding of the extent of this behaviour and its evolution in neutron-rich nuclei remains an open question in nuclear physics 3-5 . The indium isotopes are considered a textbook example of this phenomenon 6 , in which the constancy of their electromagnetic properties indicated that a single unpaired proton hole can provide the identity of a complex many-nucleon system 6,7 . Here we present precision laser spectroscopy measurements performed to investigate the validity of this simple single-particle picture. Observation of an abrupt change in the dipole moment at N = 82 indicates that, whereas the single-particle picture indeed dominates at neutron magic number N = 82 (refs.  2,8 ), it does not for previously studied isotopes. To investigate the microscopic origin of these observations, our work provides a combined effort with developments in two complementary nuclear many-body methods: ab initio valence-space in-medium similarity renormalization group and density functional theory (DFT). We find that the inclusion of time-symmetry-breaking mean fields is essential for a correct description of nuclear magnetic properties, which were previously poorly constrained. These experimental and theoretical findings are key to understanding how seemingly simple single-particle phenomena naturally emerge from complex interactions among protons and neutrons., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

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.

Collinear laser spectroscopy is performed on the nickel isotopes ^{58-68,70}Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii R_{c} are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO_{sat}, which allows an assessment of their accuracy. We find agreement with experiment in differential radii δ⟨r_{c}^{2}⟩ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO_{sat}. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy(Δr,HFB).

Heavy atomic nuclei have an excess of neutrons over protons, which leads to the formation of a neutron skin whose thickness is sensitive to details of the nuclear force. This links atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by orders of magnitude. The nucleus 208 Pb is of particular interest because it exhibits a simple structure and is experimentally accessible. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools and emulator technology, we make quantitative predictions for the properties of 208 Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore 10 9 different nuclear force parameterizations via history matching, confront them with data in select light nuclei and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of 208 Pb and determine the neutron skin thickness, which is smaller and more precise than a recent extraction from parity-violating electron scattering but in agreement with other experimental probes. This work demonstrates how realistic two- and three-nucleon forces act in a heavy nucleus and allows us to make quantitative predictions across the nuclear landscape., Competing Interests: Competing interestsThe authors declare no competing interests., (© The Author(s) 2022.)

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.

We predict the limits of existence of atomic nuclei, the proton and neutron drip lines, from the light through medium-mass regions. Starting from a chiral two- and three-nucleon interaction with good saturation properties, we use the valence-space in-medium similarity renormalization group to calculate ground-state and separation energies from helium to iron, nearly 700 isotopes in total. We use the available experimental data to quantify the theoretical uncertainties for our ab initio calculations towards the drip lines. Where the drip lines are known experimentally, our predictions are consistent within the estimated uncertainty. For the neutron-rich sodium to chromium isotopes, we provide predictions to be tested at rare-isotope beam facilities.

A recently modified method to enable low-energy nuclear scattering results to be extracted from the discrete energy levels of the target-projectile clusters confined by harmonic potential traps is tested. We report encouraging results for neutron-α and neutron-^{24}O elastic scattering from analyzing the trapped levels computed using two different ab initio nuclear structure methods. The n-α results have also been checked against a direct ab initio reaction calculation. The n-^{24}O results demonstrate the approach's applicability for a large range of systems provided their spectra in traps can be computed by ab initio methods. A key ingredient is a rigorous understanding of the errors in the calculated energy levels caused by inevitable Hilbert-space truncations in the ab initio methods.

Nuclear magic numbers correspond to fully occupied energy shells of protons or neutrons inside atomic nuclei. Doubly magic nuclei, with magic numbers for both protons and neutrons, are spherical and extremely rare across the nuclear landscape. Although the sequence of magic numbers is well established for stable nuclei, experimental evidence has revealed modifications for nuclei with a large asymmetry between proton and neutron numbers. Here we provide a spectroscopic study of the doubly magic nucleus 78 Ni, which contains fourteen neutrons more than the heaviest stable nickel isotope. We provide direct evidence of its doubly magic nature, which is also predicted by ab initio calculations based on chiral effective-field theory interactions and the quasi-particle random-phase approximation. Our results also indicate the breakdown of the neutron magic number 50 and proton magic number 28 beyond this stronghold, caused by a competing deformed structure. State-of-the-art phenomenological shell-model calculations reproduce this shape coexistence, predicting a rapid transition from spherical to deformed ground states, with 78 Ni as the turning point.