arXiv:1811.03000 · Physical review. E

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

Alexander A. Kunitsa, So Hirata

Published Jan 27, 2020• 12 citations

Abstract

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. The method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015E_{h} of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10^{-5}E_{h} without an importance sampling or Jastrow factor.

DOI: 10.1103/physreve.101.013311

Related papers

QUANTUM MONTE CARLO SIMULATIONS OF FERMIONS: A MATHEMATICAL ANALYSIS OF THE FIXED-NODE APPROXIMATION

Eric Cancès B, Benjamin Jourdain A, Tony Lelièvre A

The Diffusion Monte Carlo (DMC) method is a powerful strategy to estimate the ground state energy E0 of an N-body Schrödinger Hamiltonian H = -½Δ + V with high accuracy. It consists of writing E0 as the long-time limit of an expectation value of a drift-diffusion process with a source term, and numerically simulating this process by means of a collection of random walkers. As for a number of stochastic methods, a DMC calculation makes use of an importance sampling function ψI which hopefully approximates some ground state ψ0 of H. In the fermionic case, it has been observed that the DMC method is biased, except in the special case when the nodal surfaces of ψI coincide with those of a ground state of H. The approximation due to the fact that, in practice, the nodal surfaces of ψI differ from those of the ground states of H, is referred to as the Fixed Node Approximation (FNA). Our purpose in this paper is to provide a mathematical analysis of the FNA. We prove that, under convenient hypotheses, a DMC calculation performed with the importance sampling function ψI, provides an estimation of the infimum of the energy 〈ψ, Hψ〉 on the set of the fermionic test functions ψ that exactly vanish on the nodal surfaces of ψI.

Mathematical Models and Methods in Applied SciencesAug 23, 200627 citations

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

Fernando Agustín Reboredo, Randolph Q. Hood, Paul R. C. Kent

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. The formalism is based on the DMC mixed estimator of the ground-state probability density. We take advantage of a basic property of the walker configuration distribution generated in a DMC calculation, to (i) project out a multideterminant 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 noninteracting 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 toward the ones of the exact many-body ground state in a simulated annealing-like process. Based on these principles, we propose a method to improve both single determinant and multideterminant expansions of the trial wave function. The method can be generalized to other wave-function forms such as pfaffians. We test the method in a model system where benchmark configuration-interaction calculations can be performed and most components of the Hamiltonian are evaluated analytically. 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. Similarly, the DMC total energy and density converges to the exact solutions for the model. In the optimization process we find an optimal noninteracting nodal potential of density-functional-like form whose existence was predicted in a previous publication [Phys. Rev. B 77, 245110 (2008)]. Tests of the method are extended to a model system with a conventional Coulomb interaction where we show we can obtain the exact Kohn-Sham effective potential from the DMC data.

arXiv:0808.3417Physical Review BMay 19, 200930 citations

Diffusion Monte Carlo for fermions with replica reduction

Myung-Hoon Chung, D. P. Landau

We propose an improved framework of diffusion Monte Carlo method for quantum many-fermion systems. This approach enables us to find the ground state of the system without any a priori information on the nodal structure of the wave function. The main idea in this approach is to use sign-carrying replicas which are the basis vectors in a restricted Hilbert space. For efficiency, we propose to reduce the number of replicas in a proper manner. We test our method on a spinless fermion system, and compare the results with the exact values. We also discuss this method in relation to exact diagonalization.

Physical Review BMar 15, 20123 citations

Phase-free quantum Monte Carlo method: random walks using general basis sets

Henry Krakauer

Fermion quantum Monte Carlo (QMC) methods that work in a general basis space may be more effective for some problems than the traditional diffusion QMC method. Projection of the ground state energy is achieved by random walks in the space of Slater determinants whose one‐particle orbitals are expressed in terms of the chosen basis set. A complication is that the determinants will in general acquire complex phases. This a consequence of ground state Monte Carlo projection using the Hubbard‐Stratonovich transformation of the two‐body interaction. To control the resulting “sign” decay, we describe a method we have recently introduced for the propagation of phaseless determinants. The approximation relies on importance sampling with a trial wave function. The approximation has features in common with diffusion MC fixed node and lattice‐model constrained path methods. Using a planewave basis and non‐local pseudopotentials, we apply the method to Si atom, dimer, and 2, 16, 54 atom (216 electrons) bulk supercells. Single Slater determinant wave functions from density functional theory calculations were used as |ΨT〉 with no additional optimization. The calculated binding energy of Si2 and cohesive energy of bulk Si are in excellent agreement with experiments and are comparable to the best existing theoretical results.

AIP conference proceedingsJan 1, 20030 citations

The fermion Monte Carlo revisited

Roland Assaraf, Michel Caffarel, Anatole Khélif

In this work we present a detailed study of the fermion Monte Carlo algorithm (FMC), a recently proposed stochastic method for calculating fermionic ground-state energies. A proof that the FMC method is an exact method is given. In this work the stability of the method is related to the difference between the lowest (bosonic-type) eigenvalue of the FMC diffusion operator and the exact Fermi energy. It is shown that within a FMC framework the lowest eigenvalue of the new diffusion operator is no longer the bosonic ground-state eigenvalue as in standard exact diffusion Monte Carlo (DMC) schemes but a modified value which is strictly greater. Accordingly, FMC can be viewed as an exact DMC method built from a correlated diffusion process having a reduced Bose–Fermi gap. As a consequence, the FMC method is more stable than any transient method (or nodal release-type approaches). It is shown that the most recent ingredient of the FMC approach (Kalos M H and Pederiva F 2000 Phys. Rev. Lett. 85 3547), namely the introduction of non-symmetric guiding functions, does not necessarily improve the stability of the algorithm. We argue that the stability observed with such guiding functions is in general a finite-size population effect disappearing for a very large population of walkers. The counterpart of this stability is a control population error which is different in nature from the standard diffusion Monte Carlo algorithm and which is at the origin of an uncontrolled approximation in FMC. We illustrate the various ideas presented in this work with calculations performed on a very simple model having only nine states but a full ‘sign problem’. Already for this toy model it is clearly seen that FMC calculations are inherently uncontrolled.

arXiv:cond-mat/0611264Journal of Physics A Mathematical and TheoreticalJan 23, 200719 citations

A brief introduction to the diffusion Monte Carlo method and the fixed-node approximation

Alfonso Annarelli, Dario Alfè, Andrea Zen

Quantum Monte Carlo (QMC) methods represent a powerful family of computational techniques for tackling complex quantum many-body problems and performing calculations of stationary state properties. QMC is among the most accurate and powerful approaches to the study of electronic structure, but its application is often hindered by a steep learning curve; hence it is rarely addressed in undergraduate and postgraduate classes. This tutorial is a step toward filling this gap. We offer an introduction to the diffusion Monte Carlo (DMC) method, which aims to solve the imaginary time Schrödinger equation through stochastic sampling of the configuration space. Starting from the theoretical foundations, the discussion leads naturally to the formulation of a step-by-step algorithm. To illustrate how the method works in simplified scenarios, examples such as the harmonic oscillator and the hydrogen atom are provided. The discussion extends to the fixed-node approximation, a crucial approach for addressing the fermionic sign problem in multi-electron systems. In particular, we examine the influence of trial wave function nodal surfaces on the accuracy of DMC energy by evaluating results from a non-interacting two-fermion system. Extending the method to excited states is feasible in principle, but some additional considerations are needed, supported by practical insights. By addressing the fundamental concepts from a hands-on perspective, we hope this tutorial will serve as a valuable guide for researchers and students approaching DMC for the first time.

arXiv:2412.06006The Journal of Chemical PhysicsDec 28, 20249 citations