Fast Evaluation of Unbiased Atomic Forces in ab initio Variational Monte Carlo via the Lagrangian Technique

Abstract

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. However, the proposed method has a significant drawback for its applications: for a system with $N$ nuclei, one requires 6$N$ additional DFT calculations to get unbiased forces. This paper presents a way to replace the 6$N$ DFT calculations with a single coupled-perturbed Kohn-Sham calculation, following the so-called Lagrangian technique established in quantum chemistry. We also demonstrate that the developed unbiased 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 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 GGAs, but do not necessarily yield values that are very close to those of CCSD(T).

26 pages, 5 figures

References in corpus (22)

• Weak binding between two aromatic rings: feeling the van der Waals attraction by quantum Monte Carlo methods
• Stable liquid Hydrogen at high pressure by a novel ab-initio molecular dynamics
• A New Generation of Effective Core Potentials for Correlated Calculations
• Phase diagram of hydrogen and a hydrogen-helium mixture at planetary conditions by Quantum Monte Carlo simulations
• Accurate, efficient and simple forces with Quantum Monte Carlo methods
• Stable solid molecular hydrogen above 900K from a machine-learned potential trained with diffusion Quantum Monte Carlo
• Heavy-tailed random error in quantum Monte Carlo
• TREXIO: A File Format and Library for Quantum Chemistry
• Total forces in the diffusion Monte Carlo method with nonlocal pseudopotentials
• Atomic forces by quantum Monte Carlo: application to phonon dispersion calculation
• Space-warp coordinate transformation for efficient ionic force calculations in quantum Monte Carlo
• Interatomic force from neural network based variational quantum Monte Carlo
• Accurate quantum Monte Carlo forces for machine-learned force fields: Ethanol as a benchmark
• A tail-regression estimator for heavy-tailed distributions of known tail indices and its application to continuum quantum Monte Carlo data
• Alternative sampling for variational quantum Monte Carlo
• Principal deuterium Hugoniot via Quantum Monte Carlo and $\Delta$-learning
• Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo
• Training models using forces computed by stochastic electronic structure methods
• Efficient calculation of unbiased atomic forces in ab initio Variational Monte Carlo
• Reproducibility of fixed-node diffusion Monte Carlo across diverse community codes: The case of water-methane dimer
• Shifting sands of hardware and software in exascale quantum mechanical simulations
• Self-consistency error correction for accurate machine learning potentials from variational Monte Carlo