Published in

American Institute of Physics, The Journal of Chemical Physics, 3(156), 2022

DOI: 10.1063/5.0076302

Links

Tools

Export citation

Search in Google Scholar

Space-warp coordinate transformation for efficient ionic force calculations in quantum Monte Carlo

Journal article published in 2022 by Kousuke Nakano ORCID, Abhishek Raghav ORCID, Sandro Sorella ORCID
This paper was not found in any repository, but could be made available legally by the author.
This paper was not found in any repository, but could be made available legally by the author.

Full text: Unavailable

Green circle
Preprint: archiving allowed
Green circle
Postprint: archiving allowed
Orange circle
Published version: archiving restricted
Data provided by SHERPA/RoMEO

Abstract

Ab initio quantum Monte Carlo (QMC) methods are a state-of-the-art computational approach to obtaining highly accurate many-body wave functions. Although QMC methods are widely used in physics and chemistry to compute ground-state energies, calculation of atomic forces is still under technical/algorithmic development. Very recently, force evaluation has started to become of paramount importance for the generation of machine-learning force-field potentials. Nevertheless, there is no consensus regarding whether an efficient algorithm is available for the QMC force evaluation, namely, one that scales well with the number of electrons and the atomic numbers. In this study, we benchmark the accuracy of all-electron variational Monte Carlo (VMC) and lattice-regularized diffusion Monte Carlo (LRDMC) forces for various mono- and heteronuclear dimers (1 ≤ Z ≤ 35, where Z is the atomic number). The VMC and LRDMC forces were calculated with and without the so-called space-warp coordinate transformation (SWCT) and appropriate regularization techniques to remove the infinite variance problem. The LRDMC forces were computed with the Reynolds (RE) and variational-drift (VD) approximations. The potential energy surfaces obtained from the LRDMC energies give equilibrium bond lengths (req) and harmonic frequencies (ω) very close to the experimental values for all dimers, improving the corresponding VMC results. The LRDMC forces with the RE approximation improve the VMC forces, implying that it is worth computing the DMC forces beyond VMC despite the higher computational cost. The LRDMC forces with the VD approximations also show improvement, which unfortunately comes at a much higher computational cost in all-electron calculations. We find that the ratio of computational costs between QMC energy and forces scales as Z∼2.5 without the SWCT. In contrast, the application of the SWCT makes the ratio independent of Z. As such, the accessible QMC system size is not affected by the evaluation of ionic forces but governed by the same scaling as the total energy one.