Efficient calculation of unbiased atomic forces in ab initio variational Monte Carlo

Ab initio quantum Monte Carlo (QMC) is a state-of-the-art numerical approach for evaluating accurate expectation values of many-body wave functions. However, one of the major drawbacks that still hinders widespread QMC applications is the lack of an affordable scheme to compute unbiased atomic forces. In this study, we propose an efficient method to obtain unbiased atomic forces and pressures in the variational Monte Carlo (VMC) framework with the Jastrow-correlated Slater determinant ansatz or the Jastrow antisymmetrized geminal power ansatz, exploiting the gauge-invariant and locality properties of their geminal representation. We demonstrate the effectiveness of our method for $H_2$ and ${\mathrm{Cl}}_2$ molecules and for the cubic boron nitride crystal. Our framework has a better algorithmic scaling with the system size than the traditional finite-difference method and, in practical applications, is as efficient as single-point VMC calculations. Thus, it paves the way to study dynamical properties of materials, such as phonons, and is beneficial for pursuing more reliable machine-learning interatomic potentials based on unbiased VMC forces.