Toward Quantum Monte Carlo Forces On Heavier Ions: Scaling Properties

Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 <= Z(eff) <= 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Z(eff), increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Z(eff). The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Z(eff)(6.5(3)) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Z(eff)(-2), insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Z(eff).