Self-healing diffusion quantum Monte Carlo algorithms: methods for direct reduction of the fermion sign error in electronic structure calculations

Abstract

We develop a formalism and present an algorithm for optimization of the trial wave-function used in fixed-node diffusion quantum Monte Carlo (DMC) methods. We take advantage of a basic property of the walker configuration distribution generated in a DMC calculation, to (i) project-out a multi-determinant expansion of the fixed-node ground-state wave function and (ii) to define a cost function that relates the fixed-node ground-state and the non-interacting trial wave functions. We show that (a) locally smoothing out the kink of the fixed-node ground-state wave function at the node generates a new trial wave-function with better nodal structure and (b) we argue that the noise in the fixed-node wave-function resulting from finite sampling plays a beneficial role, allowing the nodes to adjust towards the ones of the exact many-body ground state in a simulated annealing-like process. We propose a method to improve both single determinant and multi-determinant expansions of the trial wave-function. We test the method in a model system where benchmark configuration interaction calculations can be performed. Comparing the DMC calculations with the exact solutions, we find that the trial wave-function is systematically improved. The overlap of the optimized trial wave function and the exact ground state converges to 100% even starting from wave-functions orthogonal to the exact ground state. In the optimization process we find an optimal non-interacting nodal potential of density-functional-like form whose existence was predicted earlier[Phys.Rev. B {\bf 77}, 245110 (2008)]. We obtain the exact Kohn-Sham effective potential from the DMC data.

Final version of the paper accepted in Physical Review B. The review reports and replies are included in the source

References in corpus (10)

• Inhomogeneous backflow transformations in quantum Monte Carlo calculations
• Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules
• An efficient localized basis set for quantum Monte Carlo calculations on condensed matter
• Neutral and charged excitations in carbon fullerenes from first-principles many-body theories
• Structure of fermion nodes and nodal cells
• Approximate and exact nodes of fermionic wavefunctions: coordinate transformations and topologies
• Molecular hydrogen adsorbed on benzene: insights from a quantum Monte Carlo study
• Total forces in the diffusion Monte Carlo method with nonlocal pseudopotentials
• Quantum corral wave function engineering
• Density-density functionals and effective potentials in many-body electronic structure calculations

Cited by in corpus (14)

• Continuum variational and diffusion quantum Monte Carlo calculations
• Coupling quantum Monte Carlo and independent-particle calculations: self-consistent constraint for the sign problem based on density or density matrix
• Structural Stability and Defect Energetics of ZnO from Diffusion Quantum Monte Carlo
• A new scheme for fixed node diffusion quantum Monte Carlo with pseudopotentials: improving reproducibility and reducing the trial-wave-function bias
• Diffusion quantum Monte Carlo calculations of SrFeO${}_3$ and LaFeO${}_3$
• Fixed-node errors in quantum Monte Carlo: interplay of electron density and node nonlinearities
• Systematic reduction of sign errors in many-body calculations of atoms and molecules
• A Fast and Efficient Algorithm for Slater Determinant Updates in Quantum Monte Carlo Simulations
• Fixed-Node Diffusion Monte Carlo of Lithium Systems
• Systematic reduction of sign errors in many-body problems: generalization of self-healing diffusion Monte Carlo to excited states
• Beyond single-reference fixed-node approximation in ab initio Diffusion Monte Carlo using antisymmetrized geminal power applied to systems with hundreds of electrons
• Gradient Descent Optimization of Fermion Nodes in Diffusion Monte Carlo
• Generalizing the self-healing diffusion Monte Carlo approach to finite temperature: a path for the optimization of low-energy many-body bases
• Simple Impurity Embedded in a Spherical Jellium: Approximations of Density Functional Theory compared to Quantum Monte Carlo Benchmarks