Hybridizing pseudo-Hamiltonians and non-local pseudopotentials in diffusion Monte Carlo

Abstract

An accurate treatment of effective core potentials (ECPs) requires care in continuum quantum Monte Carlo (QMC) methods. While most QMC studies have settled on the use of familiar non-local (NL) pseudopotentials with additional localization approximations, these approaches have been shown to result in moderate residual errors for some classes of molecular and solid state applications. Here, we revisit an idea proposed early in the history of QMC ECPs that does not require localization approximations, namely, a differential class of potentials referred to as pseudo-Hamiltonians. We propose to hybridize NL potentials and pseudo-Hamiltonians to reduce residual non-locality of existing potentials. We derive an approach to recast pseudopotentials for 3d elements as hybrid pseudo-Hamiltonians with optimally reduced NL energy. We demonstrate the fidelity of the hybrid potentials by studying atomic ionization potentials of Ti and Fe and the binding properties of TiO and FeO molecules with diffusion Monte Carlo (DMC). We show that localization errors have been reduced relative to potentials with the same NL channels for Sc–Zn by considering the DMC energy change with respect to the choice of approximate localization. While localization error decreases proportionate to the reduced NL energy without a Jastrow, with a Jastrow, the degree of reduction decreases at higher filling of the d-shell. Our results suggest that a subset of existing ECPs may be recast in this hybrid form to reduce the DMC localization error. They also point to the prospect of further reducing this error by generating ECPs within this hybrid form from the start.