Saturation of Nuclear Binding from Lattice Hamiltonians

Abstract

There is a conundrum regarding the binding of $\alpha$ particles in nuclei. On one hand, auxiliary-field Monte Carlo simulations of Hamiltonians on discrete spatial lattices proposed that attractive two-nucleon potentials, alone or together with attractive three-nucleon potentials, yield accurate nuclear binding. On the other hand, such Hamiltonians typically overbind all but the lightest nuclei in continuum-space approaches. We address this puzzle by performing Hartree-Fock computations of the light nuclei ${}^4$He, ${}^8$Be, ${}^{12}$C, and ${}^{16}$O, and of nuclear and neutron matter using established lattice Hamiltonians. These variational upper bounds for the ground-state energies show that the Hamiltonians with only two-nucleon potentials do not yield accurate binding, in contrast to the results from auxiliary-field Monte Carlo simulations. The case is different for Hamiltonians with three-nucleon potentials although it is the dense packing on the lattice -- and not repulsive potentials -- that yield a constant binding energy per nucleon.

Comments: 14 pages total, 4 figures