No Arabic abstract
An {em ab initio} (i.e., from first principles) theoretical framework capable of providing a unified description of the structure and low-energy reaction properties of light nuclei is desirable to further our understanding of the fundamental interactions among nucleons, and provide accurate predictions of crucial reaction rates for nuclear astrophysics, fusion-energy research, and other applications. In this contribution we review {em ab initio} calculations for nucleon and deuterium scattering on light nuclei starting from chiral two- and three-body Hamiltonians, obtained within the framework of the {em ab initio} no-core shell model with continuum. This is a unified approach to nuclear bound and scattering states, in which square-integrable energy eigenstates of the $A$-nucleon system are coupled to $(A-a)+a$ target-plus-projectile wave functions in the spirit of the resonating group method to obtain an efficient description of the many-body nuclear dynamics both at short and medium distances and at long ranges.
In any finite system, the presence of a non-zero permanent electric dipole moment (EDM) would indicate CP violation beyond the small violation predicted in the Standard Model. Here, we use the ab initio no-core shell model (NCSM) framework to theoretically investigate the magnitude of the nuclear EDM. We calculate EDMs of several light nuclei using chiral two- and three-body interactions and a PT-violating Hamiltonian based on a one-meson-exchange model. We present a benchmark calculation for $^3$He, as well as results for the more complex nuclei $^{6,7}$Li, $^9$Be, $^{10,11}$B, $^{13}$C, $^{14,15}$N, and $^{19}$F. Our results suggest that different nuclei can be used to probe different terms of the PT violating interaction. These calculations allow us to suggest which nuclei may be good candidates in the search for a measurable permanent electric dipole moment.
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_{3max}=28$, well beyond the previous limit of 18. We demonstrate that this is sufficient to obtain ground-state energies in $^{132}$Sn converged to within a few MeV with respect to the $E_{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$ 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.
Emergent properties such as nuclear saturation and deformation, and the effects on shell structure due to the proximity of the scattering continuum and particle decay channels are fascinating phenomena in atomic nuclei. In recent years, ab initio approaches to nuclei have taken the first steps towards tackling the computational challenge of describing these phenomena from Hamiltonians with microscopic degrees of freedom. This endeavor is now possible due to ideas from effective field theories, novel optimization strategies for nuclear interactions, ab initio methods exhibiting a soft scaling with mass number, and ever-increasing computational power. This paper reviews some of the recent accomplishments. We also present new results. The recently optimized chiral interaction NNLO$_{rm sat}$ is shown to provide an accurate description of both charge radii and binding energies in selected light- and medium-mass nuclei up to $^{56}$Ni. We derive an efficient scheme for including continuum effects in coupled-cluster computations of nuclei based on chiral nucleon-nucleon and three-nucleon forces, and present new results for unbound states in the neutron-rich isotopes of oxygen and calcium. The coupling to the continuum impacts the energies of the $J^pi = {1/2}^-,{3/2}^-,{7/2}^-,{3/2}^+$ states in $^{17,23,25}$O, and - contrary to naive shell-model expectations - the level ordering of the $J^pi = {3/2}^+,{5/2}^+,{9/2}^+$ states in $^{53,55,61}$Ca.
A quantitative and predictive microscopic theoretical framework that can describe reactions induced by $alpha$ particles ($^4$He nuclei) and heavier projectiles is currently lacking. Such a framework would contribute to reducing uncertainty in the modeling of stellar evolution and nucleosynthesis and provide the basis for achieving a comprehensive understanding of the phenomenon of nuclear clustering (the organization of protons and neutrons into distinct substructures within a nucleus). We have developed an efficient and general configuration-interaction framework for the description of low-energy reactions and clustering in light nuclei. The new formalism takes full advantage of powerful second-quantization techniques, enabling the description of $alpha$-$alpha$ scattering and an exploration of clustering in the exotic $^{12}$Be nucleus. We find that the $^4$He($alpha$, $alpha$)$^4$He differential cross section computed with non-locally regulated chiral interactions is in good agreement with experimental data. Our results for $^{12}$Be indicate the presence of strongly mixed helium-cluster states consistent with a molecular-like picture surviving far above the $^6$He+$^6$He threshold, and reveal the strong influence of neutron decay in both the $^{12}$Be spectrum and in the $^6$He($^6$He,$alpha$)$^8$He cross section. We expect that this approach will enable the description of helium burning cross sections and provide insight on how three-nucleon forces influence the emergence of clustering in nuclei.
We report ab initio benchmark calculations of nuclear matrix elements (NMEs) for neutrinoless double-beta ($0 ubetabeta$) 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 ($Delta T=0$) and isospin-changing ($Delta 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 $Delta T=0$ transitions are in good agreement among all calculations, at the level of 10%. For $Delta 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 ubetabeta$ decays in heavier candidate nuclei, where quasi-exact benchmarks are not available.