Self-Consistent Electron-Nucleus Cusp Correction for Molecular Orbitals

Abstract

We describe a method for imposing the correct electron-nucleus (e-n) cusp in molecular orbitals expanded as a linear combination of (cuspless) Gaussian basis functions. Enforcing the e-n cusp in trial wave functions is an important asset in quantum Monte Carlo calculations as it significantly reduces the variance of the local energy during the Monte Carlo sampling. In the method presented here, the Gaussian basis set is augmented with a small number of Slater basis functions. Note that, unlike other e-n cusp correction schemes, the presence of the Slater function is not limited to the vicinity of the nuclei. Both the coefficients of these cuspless Gaussian and cusp-correcting Slater basis functions may be self-consistently optimized by diagonalization of an orbital-dependent effective Fock operator. Illustrative examples are reported for atoms (\ce{H}, \ce{He} and \ce{Ne}) as well as for a small molecular system (\ce{BeH2}). For the simple case of the \ce{He} atom, we observe that, with respect to the cuspless version, the variance is reduced by one order of magnitude by applying our cusp-corrected scheme.

arXiv comment23 pages, 5 figures

References in corpus (21)

• A Mountaineering Strategy to Excited States: Highly-Accurate Reference Energies and Benchmarks
• Jastrow correlation factor for atoms, molecules, and solids
• Optimization of quantum Monte Carlo wave functions by energy minimization
• Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules
• Reference Energies for Double Excitations
• Hybrid stochastic-deterministic calculation of the second-order perturbative contribution of multireference perturbation theory
• Selected configuration interaction dressed by perturbation
• Scheme for adding electron-nucleus cusps to Gaussian orbitals
• Two electrons on a hypersphere: a quasi-exactly solvable model
• Benchmark all-electron ab initio quantum Monte Carlo calculations for small molecules
• Excitation energies from diffusion Monte Carlo using selected Configuration Interaction nodes
• Deterministic construction of nodal surfaces within quantum Monte Carlo: the case of FeS
• Accuracy of Electronic Wave Functions in Quantum Monte Carlo: the Effect of High-Order Correlations
• Framework for constructing generic Jastrow correlation factors
• Exact wave functions of two-electron quantum rings
• Excited states of spherium
• Accurate nonrelativistic ground-state energies of 3d transition metal atoms
• Three- and four-electron integrals involving Gaussian geminals: fundamental integrals, upper bounds and recurrence relations
• Many-electron integrals over gaussian basis functions. I. Recurrence relations for three-electron integrals
• Three-electron coalescence points in two and three dimensions
• Recurrence relations for four-electron integrals over Gaussian basis functions