Meson-Nucleus Bound States with Neural-Network Quantum States

We present the first systematic calculations of $\phi$-, $\eta_c$-, and $J/\psi$-nucleus ground states up to mass number $A{=}12$ based on the HAL QCD meson-nucleon potentials at near-physical point. The $(A{+}1)$-body Schr\"odinger equation is solved with a neural-network variational Monte Carlo framework, generalized to incorporate mesonic degrees of freedom. Benchmarking on light nuclei from $^2$H to $^{12}$C yields ground-state energies consistent with experiment. Meson-nucleus bound states emerge at $A\ge2$ for $\phi$, $A\ge4$ for $J/\psi$, and $A\ge6$ for $\eta_c$. The $\phi$-nucleus systems exhibit the strongest binding, with binding energies reaching tens of MeV. The $J/\psi$-nucleus and $\eta_c$-nucleus systems are weakly bound at the few-MeV and sub-MeV scale, respectively. The binding energy per nucleon deepens nearly linearly with $A$ for charmonium systems, whereas the $\phi$-nucleus system exhibits a non-monotonic behavior peaking at $^4$He -- a distinctive hallmark of the short-range and strongly attractive $\phi N$ interaction. The meson compresses the nucleon distribution relative to the parent nucleus, and evolves from a halo configuration to one embedded inside the nucleus with increasing $A$. Our results provide predictions for future experimental searches, and establish a quantitative bridge between lattice QCD meson-nucleon interactions and the emergent many-body phenomena in meson-nucleus bound states.