Self-Healing Diffusion Monte Carlo applied to a simple fermionic model: A critical assessment of the method

Abstract

We investigate the Self-Healing Diffusion Monte Carlo (SHDMC) method using a one-dimensional model with periodic boundary conditions. An inversion symmetry is introduced to mimic the antisymmetry property of fermionic wave functions, with the bosonic and fermionic sectors being modeled by the even and odd eigenstates, respectively. As in realistic fermionic systems, the nodal structure is only partially constrained by symmetry, making this model a non-trivial testbed for nodal optimization algorithms such as SHDMC. We show that the nodal evolution under SHDMC iterations can be cast into a dynamical system exhibiting both attractive and repulsive fixed points. When applying SHDMC to this model in the absence of statistical noise, the fixed-node energy is found to increase during iterations and the node converges to an incorrect value. This occurs for any finite basis, even when the exact node is representable in it; the attractive fixed point approaches the exact node only as the basis becomes complete, while the convergence toward it becomes increasingly slow. We further show that this problem can be largely cured by modifying the nodal update criterion to give more importance to the nodal region. Our results suggest that, in the general case, an efficient nodal optimization may require both an improved nodal update criterion and a localized basis set allowing local distortions of the trial wave function, as is the case for this model.

12 pages, 11 figures