Efficient sampling of thermal averages of interacting quantum particle systems with random batches

An efficient sampling method, the pmmLang + RBM, is proposed to compute the quantum thermal average in the interacting quantum particle system. Benefiting from the random batch method (RBM), the pmmLang + RBM has the potential to reduce the complexity due to interaction forces per time step from O(NP²) to O(NP), where N is the number of beads and P is the number of particles. Although the RBM introduces a random perturbation of the interaction forces at each time step, the long time effects of the random perturbations along the sampling process only result in a small bias in the empirical measure of the pmmLang + RBM from the target distribution, which also implies a small error in the thermal average calculation. We numerically study the convergence of the pmmLang + RBM and quantitatively investigate the dependence of the error in computing the thermal average on the parameters such as batch size, time step, and so on. We also propose an extension of the pmmLang + RBM, which is based on the splitting Monte Carlo method and is applicable when the interacting potential contains a singular part.