Basis set incompleteness errors in fixed-node diffusion Monte Carlo calculations on non-covalent interactions

Abstract

Basis set incompleteness error (BSIE) is a common source of error in quantum chemistry (QC) calculations, but it has not been comprehensively studied in fixed-node Diffusion Monte Carlo (FN-DMC) calculations. FN-DMC, being a projection method, is often considered minimally affected by basis set biases. Here, we show that this assumption is not always valid. While the relative error introduced by a small basis set in the total FN-DMC energy is minor, it can become significant in binding energy ($E_{\rm b}$) evaluations of weakly interacting systems. We systematically investigated BSIEs in FN-DMC-based binding energy ($E_{\rm b}$) evaluations using the A24 dataset, a well-known benchmark set of 24 non-covalently bound dimers. Contrary to common expectations, we found that BSIEs in FN-DMC evaluations of $E_{\rm b}$ are indeed significant when small localized basis sets, such as cc-pVDZ, are employed. We observed that BSIEs are larger in dimers with hydrogen-bonding interactions and smaller in dispersion-dominated interactions. We also found that augmenting the basis sets with diffuse orbitals, using counterpoise (CP) correction, or both, effectively mitigates BSIEs.