Medium-mass nuclei with neural quantum states

Abstract

We compute ground-state energies and charge radii of light- to medium-mass nuclei with up to $A=58$ nucleons, leveraging a variational Monte Carlo method based on Pfaffian-Jastrow neural quantum states. To further understand which elements of the nuclear Hamiltonian are "essential" to predict binding energies and charge radii across the nuclear chart with few-percent errors, we consider different interactions inspired by pionless effective field theory. Specifically, in addition to model "o" of [Phys. Rev. C 103, 054003 (2021)], we study the impact of charge-symmetry-breaking and charge-dependent terms in the nucleon-nucleon force, as well as $p$-wave contributions, which have been found to be critical for the stability of $p$-shell nuclei. In addition to its intrinsic interest, our work assesses the performance of neural quantum states in the medium-mass regime and examines the impact of these interaction modifications. Using the resulting ground-state simulations, we analyze the computational scaling of variational Monte Carlo with neural quantum states as a function of system size and computational resources, enabling projections for future large-scale calculations.