Synergy between deep neural networks and the variational Monte Carlo method for small ${}^{4}H{e}_N$ clusters

We introduce a neural network-based approach for modeling wave functions that satisfy Bose-Einstein statistics. Applying this model to small ${}^{4}H{e}_N$ clusters (with N ranging from 2 to 14 atoms), we accurately predict ground state energies, pair density functions, and two-body contact parameters $C_2^{(N)}$ related to weak unitarity. The results obtained via the variational Monte Carlo method exhibit remarkable agreement with previous studies using the diffusion Monte Carlo method, which is considered exact within its statistical uncertainties. This indicates the effectiveness of our neural network approach for investigating many-body systems governed by Bose-Einstein statistics.