No Arabic abstract
We introduce a unified approach to nuclear bound and continuum states based on the coupling of the no-core shell model (NCSM), a bound-state technique, with the no-core shell model/resonating group method (NCSM/RGM), a nuclear scattering technique. This new ab initio method, no-core shell model with continuum (NCSMC), leads to convergence properties superior to either NCSM or NCSM/RGM while providing a balanced approach to different classes of states. In the NCSMC, the ansatz for the many-nucleon wave function includes: i) a square-integrable A-nucleon component expanded in a complete harmonic oscillator basis; ii) a binary-cluster component with asymptotic boundary conditions that can properly describe weakly-bound states, resonances and scattering; and, in principle, iii) a three-cluster component suitable for the description of, e.g., Borromean halo nuclei and reactions with final three-body states. The Schroedinger equation is transformed into a system of coupled-channel integral-differential equations that we solve using a modified microscopic R-matrix formalism within a Lagrange mesh basis. We demonstrate the usefulness of the approach by investigating the unbound 7He nucleus.
Nuclear structure and reaction theory is undergoing a major renaissance with advances in many-body methods, strong interactions with greatly improved links to Quantum Chromodynamics (QCD), the advent of high performance computing, and improved computational algorithms. Predictive power, with well-quantified uncertainty, is emerging from non-perturbative approaches along with the potential for guiding experiments to new discoveries. We present an overview of some of our recent developments and discuss challenges that lie ahead. Our foci include: (1) strong interactions derived from chiral effective field theory; (2) advances in solving the large sparse matrix eigenvalue problem on leadership-class supercomputers; (3) selected observables in light nuclei with the JISP16 interaction; (4) effective electroweak operators consistent with the Hamiltonian; and, (5) discussion of A=48 system as an opportunity for the no-core approach with the reintroduction of the core.
The neutron rich exotic unbound 7He nucleus has been the subject of many experimental investigations. While the ground-state 3/2- resonance is well established, there is a controversy concerning the excited 1/2- resonance reported in some experiments as low-lying and narrow (E_R ~ 1 MeV, Gamma < 1 MeV) while in others as very broad and located at a higher energy. This issue cannot be addressed by ab initio theoretical calculations based on traditional bound-state methods. We introduce a new unified approach to nuclear bound and continuum states based on the coupling of the no-core shell model, a bound-state technique, with the no-core shell model/resonating group method, a nuclear scattering technique. Our calculations describe the ground-state resonance in agreement with experiment and, at the same time, predict a broad 1/2- resonance above 2 MeV.
Gamow shell model (GSM) is usually performed within the Woods-Saxon (WS) basis in which the WS parameters need to be determined by fitting experimental single-particle energies including their resonance widths. In the multi-shell case, such a fit is difficult due to the lack of experimental data of cross-shell single-particle energies and widths. In this paper, we develop an {it ab-initio} GSM by introducing the Gamow Hartree-Fock (GHF) basis that is obtained using the same interaction as the one used in the construction of the shell-model Hamiltonian. GSM makes use of the complex-momentum Berggren representation, then including resonance and continuum components. Hence, GSM gives a good description of weakly bound and unbound nuclei. Starting from chiral effective field theory and employing many-body perturbation theory (MBPT) (called nondegenerate $hat Q$-box folded-diagram renormalization) in the GHF basis, a multi-shell Hamiltonian ({it sd-pf} shells in this work) can be constructed. The single-particle energies and their resonance widths can also been obtained using MBPT. We investigated $^{23-28}$O and $^{23-31}$F isotopes, for which multi-shell calculations are necessary. Calculations show that continuum effects and the inclusion of the {it pf} shell are important elements to understand the structure of nuclei close to and beyond driplines.
We present an ab initio approach for the description of collective excitations and transition strength distributions of arbitrary nuclei up into the sd-shell that based on the No-Core Shell Model in combination with the Lanczos strength-function method. Starting from two- and three-nucleon interactions from chiral effective field theory, we investigate the electric monopole, dipole, and quadrupole response of the even oxygen isotopes from 16-O to 24-O. The method describes the full energy range from low-lying excitations to the giant resonance region and beyond in a unified and consistent framework, including a complete description of fragmentation and fine-structure. This opens unique opportunities for understanding dynamic properties of nuclei from first principles and to further constrain nuclear interactions. We demonstrate the computational efficiency and the robust model-space convergence of our approach and compare to established approximate methods, such as the Random Phase Approximation, shedding new light on their deficiencies.
A novel machine learning approach is used to provide further insight into atomic nuclei and to detect orderly patterns amidst a vast data of large-scale calculations. The method utilizes a neural network that is trained on ab initio results from the symmetry-adapted no-core shell model (SA-NCSM) for light nuclei. We show that the SA-NCSM, which expands ab initio applications up to medium-mass nuclei by using dominant symmetries of nuclear dynamics, can reach heavier nuclei when coupled with the machine learning approach. In particular, we find that a neural network trained on probability amplitudes for $s$-and $p$-shell nuclear wave functions not only predicts dominant configurations for heavier nuclei but in addition, when tested for the $^{20}$Ne ground state, it accurately reproduces the probability distribution. The nonnegligible configurations predicted by the network provide an important input to the SA-NCSM for reducing ultra-large model spaces to manageable sizes that can be, in turn, utilized in SA-NCSM calculations to obtain accurate observables. The neural network is capable of describing nuclear deformation and is used to track the shape evolution along the $^{20-42}$Mg isotopic chain, suggesting a shape-coexistence that is more pronounced toward the very neutron-rich isotopes. We provide first descriptions of the structure and deformation of $^{24}$Si and $^{40}$Mg of interest to x-ray burst nucleosynthesis, and even of the extremely heavy nuclei such as $^{166,168}$Er and $^{236}$U, that build upon first principles considerations.