Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry

Abstract

Diffusion Monte Carlo (DMC) is being recognized as a higher-accuracy, albeit more computationally expensive, alternative to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the use of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that this computational cost can be significantly reduced by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like binding energies and mapping potential energy surfaces. We test the error cancellation and convergence in two well-known adsorbate/slab systems, H2O/LiH(001) and CO/Pt(111). Based on this, we identify strategies for obtaining binding energies in the thermodynamic limit that optimize error cancellation to balance accuracy and computational efficiency. We then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for DFT. Our accurate, inexpensive DMC calculations recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.

References in corpus (13)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Beyond the locality approximation in the standard diffusion Monte Carlo method
• Finite-size correction in many-body electronic structure calculations
• Water on hexagonal boron nitride from diffusion Monte Carlo
• Phase Stability of TiO${}_2$ Polymorphs from Diffusion Quantum Monte Carlo
• Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo
• The Stability, Energetics, and Magnetic States of Cobalt Adatoms on Graphene
• Hydrogen dissociation on the Mg(0001) surface from quantum Monte Carlo calculations
• Ab Initio Many-body Study of Cobalt Adatoms Adsorbed on Graphene
• Water on BN doped benzene: A hard test for exchange-correlation functionals and the impact of exact exchange on weak binding
• Nonlocal pseudopotentials and time-step errors in diffusion Monte Carlo
• A pathway towards high throughput Quantum Monte Carlo simulations for alloys: A case study of two-dimensional (2D) $Ga{S}_{x}S{e}_{1-x}$
• The Nature of Interlayer Binding and Stacking of $sp$-$s{p}^2$ Hybridized Carbon Layers: A Quantum Monte Carlo Study