Hypernuclei with neural network quantum states

Leveraging complementary machine-learning-based approaches, we compute properties of s - and p -shell Λ hypernuclei—including binding energies, single-particle densities, and radii—starting from the individual interactions among their constituents. These interactions are modeled based on a leading-order finite-cutoff pionless effective field theory expansion. The couplings and the range of the three-nucleon and Λ -nucleon-nucleon potentials are determined via a Gaussian process framework anchored on virtually exact few-body techniques. We solve the many-body Schrödinger equation using a variational Monte Carlo method based on neural network quantum states, extending it for the first time to include Λ particles alongside protons and neutrons. The predicted binding energies show remarkably good agreement with experimental results, given the simplicity of the input Hamiltonian. We also confirm the experimentally observed shrinkage of the proton radius in Li Λ 7 compared to its parent nucleus, Li 6 . This work paves the way for an description of medium-mass and heavy hypernuclei, which is critical for understanding the onset of strange degrees of freedom in the core of neutron stars.