Gradient Descent Optimization of Fermion Nodes in Diffusion Monte Carlo

Abstract

We present a method for optimizing the location of the fermion ground-state nodes using a combination of diffusion Monte Carlo (DMC) and projected gradient descent (PGD). A PGD iteration shifts the parameters of an arbitrary node-fixing trial function in the opposite direction of the DMC energy gradient, while maintaining the cusp condition for atomic electrons. The energy gradient is calculated from DMC walker distributions by one of three methods we derive from an exact analytical expression. We combine our energy gradient calculation methods with different gradient descent algorithms and a projection operator that maintains the cusp condition. We apply this stochastic PGD method to trial functions of Be, Li${}_2$, and Ne, all consisting of a single Slater determinant with randomized parameters, and find that the nodes dramatically improve to the same DMC energy as nodes optimized by variational Monte Carlo. Our method, therefore, departs from the standard procedure of optimizing the nodes with a non-DMC scheme such as variational Monte Carlo, Density function theory, or configuration interaction based calculation, which do not directly minimize the DMC energy.

arXiv comment11 double spaced pages 6 figures

References in corpus (12)

• Adam: A Method for Stochastic Optimization
• Quantum Monte Carlo method using phase-free random walks with Slater determinants
• Absence of superconductivity in the pure two-dimensional Hubbard model
• Quantum Monte Carlo calculations of excited states in A = 6--8 nuclei
• A Constrained Path Quantum Monte Carlo Method for Fermion Ground States
• Energy and variance optimization of many body wave functions
• Phase separation at all interaction strengths in the t-J model
• Fixed-Node Diffusion Monte Carlo potential energy curve of the fluorine molecule F2 using selected configuration interaction trial wavefunctions
• Atomic nuclei from quantum Monte Carlo calculations with chiral EFT interactions
• Fixed-node errors in quantum Monte Carlo: interplay of electron density and node nonlinearities
• Self-healing diffusion quantum Monte Carlo algorithms: methods for direct reduction of the fermion sign error in electronic structure calculations
• Optimisation of quantum Monte Carlo wave function: steepest descent method

Cited by in corpus (1)

• Beyond single-reference fixed-node approximation in ab initio Diffusion Monte Carlo using antisymmetrized geminal power applied to systems with hundreds of electrons