Finite volume analysis on systematics of the derivative expansion in HAL QCD method

We study the convergence of the derivative expansion in HAL QCD method from the finite volume analysis. Employing the (2+1)-flavor lattice QCD data obtained at nearly physical light quark masses $(m_\pi, m_K) \simeq (146, 525)$ MeV and the physical charm quark mass, we study two representative systems, $\Omega\Omega$ and $\Omega_{ccc}\Omega_{ccc}$ in the $^1S_0$ channel, where both systems were found to have a shallow bound state in our previous studies. The HAL QCD potentials are determined at the leading-order in the derivative expansion, from which finite-volume eigenmodes are obtained. Utilizing the eigenmode projection, we find that the correlation functions are dominated by the ground state (first excited state) in the case of $\Omega\Omega$ ($\Omega_{ccc}\Omega_{ccc}$). In both $\Omega\Omega$ and $\Omega_{ccc}\Omega_{ccc}$, the spectra obtained from eigenmode-projected temporal correlators are found to be consistent with those from the HAL QCD potential for both the ground and first excited state. These results show that the derivative expansion is well converged in these systems, and also provide a first explicit evidence that the HAL QCD method enables us to reliably extract the binding energy of the ground state even from the correlator dominated by excited scattering states.