BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
S. R. Stroberg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This white paper was submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting in preparation for the next NSAC Long Range Plan. We advocate to support current and future theoretical and experimental searches for physics beyond the Standard Model using nuclear $β$ decay.
We use chiral interactions and several {\em ab initio} methods to compute the nuclear matrix elements (NMEs) for ground-state to ground-state double Gamow-Teller transitions in a range of isotopes, and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the {\em ab initio} calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double-beta decay.
We propose a novel storage scheme for three-nucleon (3N) interaction matrix elements relevant for the normal-ordered two-body approximation used extensively in ab initio calculations of atomic nuclei. This scheme reduces the required memory by approximately two orders of magnitude, which allows the generation of 3N interaction matrix elements with the standard truncation of $E_{\rm 3max}=28$, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain the ground-state energy of $^{132}$Sn converged to within a few MeV with respect to the $E_{\rm 3max}$ truncation.In addition, we study the asymptotic convergence behavior and perform extrapolations to the un-truncated limit. Finally, we investigate the impact of truncations made when evolving free-space 3N interactions with the similarity renormalization group. We find that the contribution of blocks with angular momentum $J_{\rm rel}>9/2$ to the ground-state energy is dominated by a basis-truncation artifact which vanishes in the large-space limit, so these computationally expensive components can be neglected. For the two sets of nuclear interactions employed in this work, the resulting binding energy of $^{132}$Sn agrees with the experimental value within theoretical uncertainties. This work enables converged ab initio calculations of heavy nuclei.
Background: Electric-quadrupole ($E2$) strengths relate to the underlying quadrupole deformation of a nucleus and present a challenge for many nuclear theories. Mirror nuclei in the vicinity of the line of $N=Z$ represent a convenient laboratory for testing deficiencies in such models, making use of the isospin-symmetry of the systems. Purpose: Uncertainties associated with literature $E2$ strengths in \textsuperscript{23}Mg are some of the largest in $T_z=\left|\frac{1}{2}\right|$ nuclei in the $sd$-shell. The purpose of the present work is to improve the precision with which these values are known, to enable better comparison with theoretical models. Methods: Coulomb-excitation measurements of $^{23}$Mg and $^{23}$Na were performed at the TRIUMF-ISAC facility using the TIGRESS spectrometer. They were used to determine the $E2$ matrix elements of mixed $E2$/$M1$ transitions. Results: Reduced $E2$ transition strengths, $B(E2)$, were extracted for \textsuperscript{23}Mg and \textsuperscript{23}Na. Their precision was improved by factors of approximately six for both isotopes, while agreeing within uncertainties with previous measurements. Conclusions: A comparison was made with both shell-model and {\it ab initio} valence-space in-medium similarity renormalization group calculations. Valence-space in-medium similarity-renormalization-group calculations were found to underpredict the absolute $E2$ strength - in agreement with previous studies.
Background: The nucleus $^{32}$Mg ($N=20$ and $Z=12$) plays a central role in the so-called "island of inversion" where in the ground states $sd$-shell neutrons are promoted to the $fp$-shell orbitals across the shell gap, resulting in the disappearance of the canonical neutron magic number $N=20$. Purpose: The primary goals of this work are to extend the level scheme of $^{32}$Mg, provide spin-parity assignments to excited states, and discuss the microscopic structure of each state through comparisons with theoretical calculations. Method: In-beam $γ$-ray spectroscopy of $^{32}$Mg was performed using two direct-reaction probes, one-neutron (two-proton) knockout reactions on $^{33}$Mg ($^{34}$Si). Final-state exclusive cross sections and parallel momentum distributions were extracted from the experimental data and compared with eikonal-based reaction model calculations combined with shell-model overlap functions. Results: Owing to the remarkable selectivity of the one-neutron and two-proton knockout reactions, a significantly updated level scheme for $^{32}$Mg, which exhibits negative-parity intruder and positive-parity normal states, was constructed. The experimental results were confronted with four different nuclear structure models. Conclusions: In some of these models, different aspects of $^{32}$Mg and the transition into the island of inversion are well described. However, unexplained discrepancies remain, and even with the help of these state-of-the-art theoretical approaches, the structure of this key nucleus is not yet fully captured.
Background: Recent developments in {\it ab initio} nuclear theory demonstrate promising results in medium- to heavy-mass nuclei. A particular challenge for many of the many-body methodologies, however, is an accurate treatment of the electric-quadrupole, $E2$, strength associated with collectivity. Purpose: The valence-space in-medium similarity renormalization group (VS-IMSRG) is a particularly powerful method for accessing medium- and high-mass nuclei but has been found to underpredict $E2$ strengths. The purpose of this work is to evaluate the isospin dependence of this underprediction. Methods: We perform a systematic comparison of valence-space in-medium similarity renormalization group (VS-IMSRG) calculations with available literature. We make use of isoscalar and isovector contributions to the $E2$ matrix elements to assess isoscalar and isovector contributions to the missing strength. Results: It is found that the $E2$ strength is consistent throughout $T_z=\left|\frac{1}{2}\right|$, $T_z=\left|1\right|$, $T_z=\left|\frac{3}{2}\right|$ and $T_z=2$ pairs within the $sd$-shell. Furthermore, no isovector contribution to the deficiency is identified. Conclusions: A comparison with toy-models and coupled-cluster calculations is used to discuss potential origins of the missing strength, which arises from missing many-particle, many-hole excitations out of the model space. The absence of any significant isovector contribution to the missing $E2$ strength indicates that the $E2$ strength discrepancy, and therefore any correction, is largely independent of the isospin of the nuclei in question.
We report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-beta ($0νββ$) decays in light nuclei with mass number ranging from $A=6$ to $A=22$. We use the transition operator derived from light-Majorana neutrino exchange and evaluate the NME with three different methods: two variants of in-medium similarity renormalization group (IMSRG) and importance-truncated no-core shell model (IT-NCSM). The same two-plus-three-nucleon interaction from chiral effective field theory is employed, and both isospin-conserving ($ΔT=0$) and isospin-changing ($ΔT=2$) transitions are studied. We compare our resulting ground-state energies and NMEs to those of recent ab initio no-core shell model and coupled-cluster calculations, also with the same inputs. We show that the NMEs of $ΔT=0$ transitions are in good agreement among all calculations, at the level of 10%. For $ΔT=2$, relative deviations are more significant in some nuclei. The comparison with the exact IT-NCSM result allows us to analyze these cases in detail, and indicates the next steps towards improving the IMSRG-based approaches. The present study clearly demonstrates the power of consistent cross-checks that are made possible by ab initio methodology. This capability is crucial for providing meaningful many-body uncertainties in the NMEs for the $0νββ$ decays in heavier candidate nuclei, where quasi-exact benchmarks are not available.
We calculate basis-space converged neutrinoless $ββ$ decay nuclear matrix elements for the lightest candidates: 48Ca, 76Ge and 82Se. 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 76Ge and 82Se, and 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 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.
We review the status of ab initio calculations of allowed beta decays (both Fermi and Gamow-Teller), within the framework of the valence-space in-medium similarity renormalization group approach.
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$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).
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.
The structure of $^{33}$Si was studied by a one-neutron knockout reaction from a $^{34}$Si beam at 98.5 MeV/u incident on a $^{9}$Be target. The prompt $γ$-rays following the de-excitation of $^{33}$Si were detected using the GRETINA $γ$-ray tracking array while the reaction residues were identified on an event-by-event basis in the focal plane of the S800 spectrometer at NSCL (National Superconducting Cyclotron Laboratory). The presently derived spectroscopic factor values, $C^2S$, for the 3/2$^+$ and 1/2$^+$ states, corresponding to a neutron removal from the $0d_{3/2}$ and $1s_{1/2}$ orbitals, agree with shell model calculations and point to a strong $N=20$ shell closure. Three states arising from the more bound $0d_{5/2}$ orbital are proposed, one of which is unbound by about 930 keV. The sensitivity of this experiment has also confirmed a weak population of 9/2$^-$ and 11/2$_{1,2}^-$ final states, which originate from a higher-order process. This mechanism may also have populated, to some fraction, the 3/2$^-$ and 7/2$^-$ negative-parity states, which hinders a determination of the $C^2S$ values for knockout from the normally unoccupied $1p_{3/2}$ and $0f_{7/2}$ orbits.
A $^{52}$Cr$(p,t)$$^{50}$Cr two-neutron pickup reaction was performed using the Q3D magnetic spectrograph at the Maier-Leibnitz-Laboratorium in Garching, Germany. Excited states in $^{50}$Cr were observed up to an excitation energy of 5.3 MeV. Despite significantly increased sensitivity and resolution over previous work, no evidence of the previously assigned first excited $0^+$ state was found. As a result, the $0^+_2$ state is reassigned at an excitation energy of $E_x=3895.0(5)$ keV in $^{50}$Cr. This reassignment directly impacts direct searches for a non-analogue Fermi $β^+$ decay branch in $^{50}$Mn. These results also show better systematic agreement with the theoretical predictions for the $0^+$ state spectrum in $^{50}$Cr using the same formalism as the isospin-symmetry-breaking correction calculations for superallowed nuclei. The experimental data are also compared to $ab$-$initio$ shell-model predictions using the IM-SRG formalism based on $NN$ and $3N$ forces from chiral-EFT in the $pf$-shell for the first time.
We present ab initio predictions for ground and excited states of doubly open-shell fluorine and neon isotopes based on chiral two- and three-nucleon interactions. We use the in-medium similarity renormalization group, in both flow-equation and Magnus formulations, to derive mass-dependent sd valence-space Hamiltonians. The experimental ground-state energies are reproduced through neutron number N=14, beyond which a new targeted normal-ordering procedure improves agreement with data and large-scale multi-reference calculations. For spectroscopy, we focus on neutron-rich 23-26F and 24-26Ne isotopes near N=14,16 magic numbers. In all cases we find an agreement with experiment competitive with established phenomenology. Moreover, yrast states are well described in 20Ne and 24Mg, providing an ab initio description of deformation in the medium-mass region.
We present a systematic study of both nuclear radii and binding energies in (even) oxygen isotopes from the valley of stability to the neutron drip line. Both charge and matter radii are compared to state-of-the-art {\it ab initio} calculations along with binding energy systematics. Experimental matter radii are obtained through a complete evaluation of the available elastic proton scattering data of oxygen isotopes. We show that, in spite of a good reproduction of binding energies, {\it ab initio} calculations with conventional nuclear interactions derived within chiral effective field theory fail to provide a realistic description of charge and matter radii. A novel version of two- and three-nucleon forces leads to considerable improvement of the simultaneous description of the three observables for stable isotopes, but shows deficiencies for the most neutron-rich systems. Thus, crucial challenges related to the development of nuclear interactions remain.
The calcium isotopes have emerged as an important testing ground for new microscopically derived shell-model interactions, and a great deal of focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Ca following one-neutron knockout reactions from $^{48,50}$Ca. The observed reduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca, as compared to $^{47}$Ca, is consistent with the description given by shell-model calculations based on two- and three-nucleon forces in the neutron $pf$ model space, implying a fragmentation of the $l$=3 strength to higher-lying states. The experimental result is inconsistent with both the GXPF1 interaction routinely used in this region of the nuclear chart and with microscopic calculations in an extended model space including the $\nu1g_{9/2}$ orbital.
Shell evolution is studied in the neutron-rich silicon isotopes 36,38,40 Si using neutron single-particle strengths deduced from one-neutron knockout reactions. Configurations involving neutron excita- tions across the N = 20 and N = 28 shell gaps are quantified experimentally in these rare isotopes. Comparisons with shell model calculations show that the tensor force, understood to drive the col- lective behavior in 42 Si with N = 28, is already important in determining the structure of 40 Si with N = 26. New data relating to cross-shell excitations provide the first quantitative support for repulsive contributions to the cross-shell T = 1 interaction arising from three-nucleon forces.
Direct reaction techniques are powerful tools to study the single-particle nature of nuclei. Performing direct reactions on short-lived nuclei requires radioactive ion beams produced either via fragmentation or the Isotope Separation OnLine (ISOL) method. Some of the most interesting regions to study with direct reactions are close to the magic numbers where changes in shell structure can be tracked. These changes can impact the final abundances of explosive nucleosynthesis. The structure of the chain of tin isotopes is strongly influenced by the Z=50 proton shell closure, as well as the neutron shell closures lying in the neutron-rich, N=82, and neutron-deficient, N=50, regions. Here we present two examples of direct reactions on exotic tin isotopes. The first uses a one-neutron transfer reaction and a low-energy reaccelerated ISOL beam to study states in 131Sn from across the N=82 shell closure. The second example utilizes a one-neutron knockout reaction on fragmentation beams of neutron-deficient 106,108Sn. In both cases, measurements of gamma rays in coincidence with charged particles proved to be invaluable.
Excited states in $^{28}$Na have been studied using the $β$-decay of implanted $^{28}$Ne ions at GANIL/LISE as well as the in-beam $γ$-ray spectroscopy at the NSCL/S800 facility. New states of positive (J$^π$=3,4$^+$) and negative (J$^π$=1-5$^-$) parity are proposed. The former arise from the coupling between 0d$\_{5/2}$ protons and a 0d$\_{3/2}$ neutron, while the latter are due to couplings with 1p$\_{3/2}$ or 0f$\_{7/2}$ neutrons. While the relative energies between the J$^π$=1-4$^+$ states are well reproduced with the USDA interaction in the N=17 isotones, a progressive shift in the ground state binding energy (by about 500 keV) is observed between $^{26}$F and $^{30}$Al. This points to a possible change in the proton-neutron 0d$\_{5/2}$-0d$\_{3/2}$ effective interaction when moving from stability to the drip line. The presence of J$^π$=1-4$^-$ negative parity states around 1.5 MeV as well as of a candidate for a J$^π$=5$^-$ state around 2.5 MeV give further support to the collapse of the N=20 gap and to the inversion between the 0f$\_{7/2}$ and 1p$\_{3/2}$ levels below Z=12. These features are discussed in the framework of Shell Model and EDF calculations, leading to predicted negative parity states in the low energy spectra of the $^{26}$F and $^{25}$O nuclei.
The 9Be(28Mg,27Na) one-proton removal reaction with a large proton separation energy of Sp(28Mg)=16.79 MeV is studied at intermediate beam energy. Coincidences of the bound 27Na residues with protons and other light charged particles are measured. These data are analyzed to determine the percentage contributions to the proton removal cross section from the elastic and inelastic nucleon removal mechanisms. These deduced contributions are compared with the eikonal reaction model predictions and with the previously measured data for reactions involving the re- moval of more weakly-bound protons from lighter nuclei. The role of transitions of the proton between different bound single-particle configurations upon the elastic breakup cross section is also quantified in this well-bound case. The measured and calculated elastic breakup fractions are found to be in good agreement.
Excited states in the neutron-rich N=38,36 nuclei \nuc{60}{Ti} and \nuc{58}{Ti} were populated in nucleon-removal reactions from \nuc{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 \sim 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, \nuc{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 \nuc{60}{Ca}.
We report on the experimental study of quadrupole collectivity in the neutron-deficient nucleus \nuc{104}{Sn} using intermediate-energy Coulomb excitation. The $B(E2; 0^+_1 \rightarrow 2^+_1)$ value for the excitation of the first $2^+$ state in \nuc{104}{Sn} has been measured to be $0.180(37)~e^2$b$^2$ relative to the well-known $B(E2)$ value of \nuc{102}{Cd}. This result disagrees by more than one sigma with a recently published measurement \cite{Gua13}. Our result indicates that the most modern many-body calculations remain unable to describe the enhanced collectivity below mid-shell in Sn approaching $N=Z=50$. We attribute the enhanced collectivity to proton particle-hole configurations beyond the necessarily limited shell-model spaces and suggest the asymmetry of the $B(E2)$-value trend around mid-shell to originate from enhanced proton excitations across $Z=50$ as $N=Z$ is approached.
We report final-state-exclusive measurements of the light charged fragments in coincidence with 26Ne residual nuclei following the direct two-proton removal from a neutron-rich 28Mg secondary beam. A Dalitz-plot analysis and comparisons with simulations show that a majority of the triple- coincidence events with two protons display phase-space correlations consistent with the (two-body) kinematics of a spatially-correlated pair-removal mechanism. The fraction of such correlated events, 56(12) %, is consistent with the fraction of the calculated cross section, 64 %, arising from spin S = 0 two-proton configurations in the entrance-channel (shell-model) 28Mg ground state wave function. This result promises access to an additional and more specific probe of the spin and spatial correlations of valence nucleon pairs in exotic nuclei produced as fast secondary beams.
The two-proton removal reaction from 28Mg projectiles has been studied at 93 MeV/u at the NSCL. First coincidence measurements of the heavy 26Ne projectile residues, the removed protons and other light charged particles enabled the relative cross sections from each of the three possible elastic and inelastic proton removal mechanisms to be determined. These more final-state-exclusive measurements are key for further interrogation of these reaction mechanisms and use of the reaction channel for quantitative spectroscopy of very neutron-rich nuclei. The relative and absolute yields of the three contributing mechanisms are compared to reaction model expectations - based on the use of eikonal dynamics and sd-shell-model structure amplitudes.
New measurements and reaction model calculations are reported for single neutron pickup reactions onto a fast \nuc{22}{Mg} secondary beam at 84 MeV per nucleon. Measurements were made on both carbon and beryllium targets, having very different structures, allowing a first investigation of the likely nature of the pickup reaction mechanism. The measurements involve thick reaction targets and $γ$-ray spectroscopy of the projectile-like reaction residue for final-state resolution, that permit experiments with low incident beam rates compared to traditional low-energy transfer reactions. From measured longitudinal momentum distributions we show that the $\nuc{12}{C} (\nuc{22}{Mg},\nuc{23}{Mg}+γ)X$ reaction largely proceeds as a direct two-body reaction, the neutron transfer producing bound \nuc{11}{C} target residues. The corresponding reaction on the \nuc{9}{Be} target seems to largely leave the \nuc{8}{Be} residual nucleus unbound at excitation energies high in the continuum. We discuss the possible use of such fast-beam one-neutron pickup reactions to track single-particle strength in exotic nuclei, and also their expected sensitivity to neutron high-$\ell$ (intruder) states which are often direct indicators of shell evolution and the disappearance of magic numbers in the exotic regime.