No Arabic abstract
Background: The nuclear kinetic density is one of many fundamental quantities in density functional theory (DFT) dependent on the nonlocal nuclear density. Often, approximations may be made when computing the density that may result in spurious contributions in other DFT quantities. With the ability to compute the nonlocal nuclear density from ab initio wave functions, it is now possible to estimate effects of such spurious contributions. Purpose: We derive the kinetic density using ab initio nonlocal scalar one-body nuclear densities computed within the no-core shell model (NCSM) approach, utilizing two- and three-nucleon chiral interactions as the sole input. With the ability to compute translationally invariant nonlocal densities, it is possible to directly gauge the impact of the spurious center-of-mass (COM) contributions in DFT quantities such as the kinetic density. Methods: The nonlocal nuclear densities are derived from the NCSM one-body densities calculated in second quantization. We present a review of COM contaminated and translationally invariant nuclear densities. We then derive an analytic expression for the kinetic density using these nonlocal densities, producing an ab initio kinetic density. Results: The ground state nonlocal densities of textsuperscript{4,6,8}He, textsuperscript{12}C, and textsuperscript{16}O are used to compute the kinetic densities of the aforementioned nuclei. The impact of the COM removal technique in the densities is discussed. The results of this work can be extended to other fundamental quantities in DFT. Conclusions: The use of a general nonlocal density allows for the calculation of fundamental quantities taken as input in theories such as DFT. This allows benchmarking of procedures for COM removal in different many-body techniques.
We present the first application of a new approach, proposed in [Journal of Physics G: Nuclear and Particle Physics, 43, 04LT01 (2016)] to derive coupling constants of the Skyrme energy density functional (EDF) from ab initio Hamiltonian. By perturbing the ab initio Hamiltonian with several functional generators defining the Skyrme EDF, we create a set of metadata that is then used to constrain the coupling constants of the functional. We use statistical analysis to obtain such an ab initio-equivalent Skyrme EDF. We find that the resulting functional describes properties of atomic nuclei and infinite nuclear matter quite poorly. This may point out to the necessity of building up the ab initio-equivalent functionals from more sophisticated generators. However, we also indicate that the current precision of the ab initio calculations may be insufficient for deriving meaningful nuclear EDFs.
We discuss the construction of a nuclear Energy Density Functional (EDF) from ab initio calculations, and we advocate the need of a methodical approach that is free from ad hoc assumptions. The equations of state (EoS) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO$_{rm sat}$ and the phenomenological AV4$^prime$+UIX$_{c}$ Hamiltonians as inputs in the Self-consistent Greens Function (SCGF) and Auxiliary Field Diffusion Monte Carlo (AFDMC) methods, respectively. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio-based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO$_{rm sat}$-based EDF produces encouraging results, whereas the AV4$^prime$+UIX$_{c}$-based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our work paves the way for a practical and systematic way to merge ab initio nuclear theory and DFT, while at the same time it sheds light on some of the critical aspects of this procedure.
We propose a new Monte Carlo method called the pinhole trace algorithm for {it ab initio} calculations of the thermodynamics of nuclear systems. For typical simulations of interest, the computational speedup relative to conventional grand-canonical ensemble calculations can be as large as a factor of one thousand. Using a leading-order effective interaction that reproduces the properties of many atomic nuclei and neutron matter to a few percent accuracy, we determine the location of the critical point and the liquid-vapor coexistence line for symmetric nuclear matter with equal numbers of protons and neutrons. We also present the first {it ab initio} study of the density and temperature dependence of nuclear clustering.
Ab initio approaches in nuclear theory, such as the no-core shell model (NCSM), have been developed for approximately solving finite nuclei with realistic strong interactions. The NCSM and other approaches require an extrapolation of the results obtained in a finite basis space to the infinite basis space limit and assessment of the uncertainty of those extrapolations. Each observable requires a separate extrapolation and most observables have no proven extrapolation method. We propose a feed-forward artificial neural network (ANN) method as an extrapolation tool to obtain the ground state energy and the ground state point-proton root-mean-square (rms) radius along with their extrapolation uncertainties. The designed ANNs are sufficient to produce results for these two very different observables in $^6$Li from the ab initio NCSM results in small basis spaces that satisfy the following theoretical physics condition: independence of basis space parameters in the limit of extremely large matrices. Comparisons of the ANN results with other extrapolation methods are also provided.
In this work we present the first steps towards benchmarking isospin symmetry breaking in ab initio nuclear theory for calculations of superallowed Fermi $beta$-decay. Using the valence-space in-medium similarity renormalization group, we calculate b and c coefficients of the isobaric multiplet mass equation, starting from two different Hamiltonians constructed from chiral effective field theory. We compare results to experimental measurements for all T=1 isobaric analogue triplets of relevance to superallowed $beta$-decay for masses A=10 to A=74 and find an overall agreement within approximately 250 keV of experimental data for both b and c coefficients. A greater level of accuracy, however, is obtained by a phenomenological Skyrme interaction or a classical charged-sphere estimate. Finally, we show that evolution of the valence-space operator does not meaningfully improve the quality of the coefficients with respect to experimental data, which indicates that higher-order many-body effects are likely not responsible for the observed discrepancies.