Fast evaluation of unbiased atomic forces in ab initio variational Monte Carlo via the Lagrangian technique

Ab initio quantum Monte Carlo (QMC) methods are state-of-the-art electronic structure calculations based on highly parallelizable stochastic frameworks for accurate solutions of the many-body Schrödinger equation, suitable for modern many-core supercomputer architectures. Despite its potential, one of the major drawbacks that still hinders QMC applications, especially when targeting dynamical properties of large systems or extensive datasets, is the lack of an affordable method to compute atomic forces that are consistent with the corresponding potential energy surfaces (PESs), also known as unbiased atomic forces. Recently, one of the authors in the present paper proposed a way to obtain unbiased forces with the Jastrow-correlated Slater determinant Ansatz, where the determinant part is frozen to the values obtained by a mean-field method, such as density functional theory [K. Nakano, M. Casula, and G. Tenti, Phys. Rev. B 109, 205151 (2024)]. However, the proposed method has a significant drawback for its applications: for a system with N nuclei, one requires 6N additional density functional theory (DFT) calculations to get unbiased forces, which is not negligible as the system size increases. This paper presents a way to replace the 6N DFT calculations with a single coupled-perturbed Kohn-Sham calculation, following the so-called Lagrangian technique established in quantum chemistry. This improves the computational cost and scalability of the method. We also demonstrate that the developed unbiased variational Monte Carlo (VMC) force calculation improves not only the consistency with PESs but also its accuracy, by investigating three molecules from the rMD17 benchmark set, and comparing the unbiased VMC forces with those obtained by the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] calculations. We found that the bare VMC forces are biased from the CCSD(T) ones, while the unbiased ones give values closer to those of the CCSD(T) ones. Our benchmark test also reveals that the unbiased VMC forces yield very consistent values with hybrid and meta generalized gradient approximations (e.g., ωB97X-D3BJ and ωB97M-D3BJ), but do not necessarily yield values that are very close to those of CCSD(T). Our finding paves the way to generate machine learning interatomic potentials based on VMC forces more efficiently and accurately.