BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
S. R. Stroberg contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph 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.