No Arabic abstract
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.
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.
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.
[Background:] It is well known that effective nuclear interactions are in general nonlocal. Thus if nuclear densities obtained from {it ab initio} no-core-shell-model (NCSM) calculations are to be used in reaction calculations, translationally invariant nonlocal densities must be available. [Purpose:] Though it is standard to extract translationally invariant one-body local densities from NCSM calculations to calculate local nuclear observables like radii and transition amplitudes, the corresponding nonlocal one-body densities have not been considered so far. A major reason for this is that the procedure for removing the center-of-mass component from NCSM wavefunctions up to now has only been developed for local densities. [Results:] A formulation for removing center-of-mass contributions from nonlocal one-body densities obtained from NCSM and symmetry-adapted NCSM (SA-NCSM) calculations is derived, and applied to the ground state densities of $^4$He, $^6$Li, $^{12}$C, and $^{16}$O. The nonlocality is studied as a function of angular momentum components in momentum as well as coordinate space [Conclusions:] We find that the nonlocality for the ground state densities of the nuclei under consideration increases as a function of the angular momentum. The relative magnitude of those contributions decreases with increasing angular momentum. In general, the nonlocal structure of the one-body density matrices we studied is given by the shell structure of the nucleus, and can not be described with simple functional forms.
The existence of multi-neutron systems has always been a debatable question. Indeed, both inter-nucleon correlations and a large continuum coupling occur in these states. We then employ the ab-initio no-core Gamow shell model to calculate the resonant energies and widths of the trineutron and tetraneutron systems with realistic interactions. Our results indicate that trineutron and tetraneutron are both unbound and bear broad widths. The calculated energy and width of tetraneutron are also comparable with recent experimental data. Moreover, our calculations suggest that the energy of trineutron is lower than that of tetraneutron, while its resonance width is also narrower. This strongly suggests that trineutron is more likely to be experimentally observed than tetraneutron. We thus suggest experimentalists to search for trineutron at low energy.
Constructing microscopic effective interactions (`optical potentials) for nucleon-nucleus (NA) elastic scattering requires in first order off-shell nucleon-nucleon (NN) scattering amplitudes between the projectile and the struck target nucleon and nonlocal one-body density matrices. While the NN amplitudes and the {it ab intio} no-core shell-model (NCSM) calculations always contain the full spin structure of the NN problem, one-body density matrices used in traditional microscopic folding potential neglect spin contributions inherent in the one-body density matrix. Here we derive and show the expectation values of the spin-orbit contribution of the struck nucleon with respect to the rest of the nucleus for $^{4}$He, $^{6}$He, $^{12}$C, and $^{16}$O and compare them with the scalar one-body density matrix.