Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1209.1663 · Physical Review E

Population size bias in diffusion Monte Carlo

Massimo Boninsegni, Saverio Moroni

Published Nov 28, 2012• 55 citations
arXivHTMLPDF

Abstract

The size of the population of random walkers required to obtain converged estimates in diffusion Monte Carlo (DMC) increases dramatically with system size. We illustrate this by comparing ground state energies of small clusters of parahydrogen (up to 48 molecules) computed by DMC and path integral ground state (PIGS) techniques. We contend that the bias associated with a finite population of walkers is the most likely cause of quantitative numerical discrepancies between PIGS and DMC energy estimates reported in the literature, for this few-body Bose system. We discuss the viability of DMC as a general-purpose ground state technique, and argue that PIGS, and even finite temperature methods, enjoy more favorable scaling, and are therefore a superior option for systems of large size.

DOI
10.1103/physreve.86.056712

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 Sciences•Aug 23, 2006•27 citations

Quantum Monte Carlo: From Few to Many Species

William G. Dawkins

We study systems of particles with varying numbers of distinguishable species ranging from 4-species to the N-species limit, where N is the number of particles in the system itself. We examine these systems using various powerful simulation techniques that belong to the quantum Monte Carlo (QMC) family. The N-species limit is studied via Path Integral Monte Carlo (PIMC), which is a finite temperature, non-perturbative technique. At this limit, we explore the viability of two different thermal density matrices for systems that interact via hard-sphere, hard-cavity potentials. We then study the impact of finite-size effects on calculations of the energy, pressure and specific heat of the system when reaching the thermodynamic limit under periodic boundary conditions. We then move to studying a 4-species fermionic system using Variational Monte Carlo (VMC) and Diffusion Monte Carlo (DMC), which are ground state techniques. We apply these methods to the clustering problem of 4-particle and 8-particle systems. While the 4-particle, 4-species system is relatively simple to calculate as it is physically identical to a bosonic system of 4-particles, the 8-particle fermi system is more complicated. We detail the process of simulating a bound 8-particle state with respect to decay into two independent 4-particle clusters with DMC for the first time.

•May 1, 2019•0 citations

Finite population-size effects in projection Monte Carlo methods

Nicolas J. Cerf, Olivier Martin

Projection (Green's function and diffusion) MOnte Carlo techniques sample a wave function by a stochastic iterative procedure. These method converge to a stationary distribution which is biased, i.e., differs from the exact ground state wave function. This bias occurs because of the use of population control procedures. We demonstrate that these biased Monte Carlo algorithms lead to a modified effective mass which is equal to the desired mass only in the limit of an infinite population of walkers. In general, the bias for the energy scales as 1/N for a population of walkers of size N and is proportional to the expectation value of the kinetic energy. Finally, we consider various strategies to reduce this bias.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Apr 1, 1995•22 citations

Path-integral diffusion Monte Carlo: Calculation of observables of many-body systems in the ground state

Balázs Hetényi, Eran Rabani, B. J. Berne

We propose a new method to calculate ground state position dependent observables in quantum many-body systems. The method, which we call the path-integral diffusion Monte Carlo (PI-DMC) method, is essentially a combination of path-integral Monte Carlo (PIMC) and diffusion Monte Carlo (DMC) methods. The distribution resulting from a DMC simulation is further propagated in imaginary time by PIMC sampling. Tests of the new method for simple cases such as the harmonic oscillator, a double well, and a set of ten coupled harmonic oscillators show that the results for several observables converge rapidly to the exact result.

The Journal of Chemical Physics•Apr 1, 1999•21 citations

Persistent Population Biases in Branching Random Walk Algorithms

David J. Geroski

We explore the existence of a residual population bias in Diffusion Monte Carlo algorithms. We develop a model problem, which mimics the single-particle, one dimensional simple harmonic oscillator of unitary ground state energy, to study these biases. We then apply an importance sampling algorithms with a population of a single random walker to calculate the ground state energy and wavefunction of our modeled system in cases of both perfect and imperfect importance sampling. In the former case, we perfectly calculate both the ground state energy and corresponding wavefunction of our system, thus validating our calculations and assumed answer. In the latter case, we show the existence of a residual population bias for populations of small numbers of random walkers. We demonstrate the convergence of our small population calculations and show that we cannot remove this bias through a number of single walker methods. Finally, we demonstrate the overall decay of this residual population bias as the population of random walkers expands, allowing our calculations to converge to the true answer. Key Words: Diffusion Monte Carlo, Population Dynamics, Random Walk

•Jan 1, 2015•0 citations

Diffusion Monte Carlo methods with a fixed number of walkers

Roland Assaraf, Michel Caffarel, Anatole Khélif

In this paper we discuss various aspects of diffusion Monte Carlo methods using a fixed number of walkers. First, a rigorous proof of the divergence of pure diffusion Monte Carlo (PDMC) methods (DMC without branching in which the weights are carried along trajectories) is given. Second, a bias-free Monte Carlo method combining DMC and PDMC approaches, and based on a minimal stochastic reconfiguration of the population, is discussed. Finally, some illustrative calculations for a system of coupled quantum rotators are presented.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Apr 1, 2000•91 citations
The size of the population of random walkers required to obtain converged estimates in diffusion Monte Carlo (DMC) increases dramatically with system size. We illustrate this by comparing ground state energies of small clusters of parahydrogen (up to 48 molecules) computed by DMC and path integral ground state (PIGS) techniques. We contend that the bias associated with a finite population of walkers is the most likely cause of quantitative numerical discrepancies between PIGS and DMC energy estimates reported in the literature, for this few-body Bose system. We discuss the viability of DMC as a general-purpose ground state technique, and argue that PIGS, and even finite temperature methods, enjoy more favorable scaling, and are therefore a superior option for systems of large size.
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.
We study systems of particles with varying numbers of distinguishable species ranging from 4-species to the N-species limit, where N is the number of particles in the system itself. We examine these systems using various powerful simulation techniques that belong to the quantum Monte Carlo (QMC) family. The N-species limit is studied via Path Integral Monte Carlo (PIMC), which is a finite temperature, non-perturbative technique. At this limit, we explore the viability of two different thermal density matrices for systems that interact via hard-sphere, hard-cavity potentials. We then study the impact of finite-size effects on calculations of the energy, pressure and specific heat of the system when reaching the thermodynamic limit under periodic boundary conditions. We then move to studying a 4-species fermionic system using Variational Monte Carlo (VMC) and Diffusion Monte Carlo (DMC), which are ground state techniques. We apply these methods to the clustering problem of 4-particle and 8-particle systems. While the 4-particle, 4-species system is relatively simple to calculate as it is physically identical to a bosonic system of 4-particles, the 8-particle fermi system is more complicated. We detail the process of simulating a bound 8-particle state with respect to decay into two independent 4-particle clusters with DMC for the first time.
Projection (Green's function and diffusion) MOnte Carlo techniques sample a wave function by a stochastic iterative procedure. These method converge to a stationary distribution which is biased, i.e., differs from the exact ground state wave function. This bias occurs because of the use of population control procedures. We demonstrate that these biased Monte Carlo algorithms lead to a modified effective mass which is equal to the desired mass only in the limit of an infinite population of walkers. In general, the bias for the energy scales as 1/N for a population of walkers of size N and is proportional to the expectation value of the kinetic energy. Finally, we consider various strategies to reduce this bias.
We propose a new method to calculate ground state position dependent observables in quantum many-body systems. The method, which we call the path-integral diffusion Monte Carlo (PI-DMC) method, is essentially a combination of path-integral Monte Carlo (PIMC) and diffusion Monte Carlo (DMC) methods. The distribution resulting from a DMC simulation is further propagated in imaginary time by PIMC sampling. Tests of the new method for simple cases such as the harmonic oscillator, a double well, and a set of ten coupled harmonic oscillators show that the results for several observables converge rapidly to the exact result.
We explore the existence of a residual population bias in Diffusion Monte Carlo algorithms. We develop a model problem, which mimics the single-particle, one dimensional simple harmonic oscillator of unitary ground state energy, to study these biases. We then apply an importance sampling algorithms with a population of a single random walker to calculate the ground state energy and wavefunction of our modeled system in cases of both perfect and imperfect importance sampling. In the former case, we perfectly calculate both the ground state energy and corresponding wavefunction of our system, thus validating our calculations and assumed answer. In the latter case, we show the existence of a residual population bias for populations of small numbers of random walkers. We demonstrate the convergence of our small population calculations and show that we cannot remove this bias through a number of single walker methods. Finally, we demonstrate the overall decay of this residual population bias as the population of random walkers expands, allowing our calculations to converge to the true answer. Key Words: Diffusion Monte Carlo, Population Dynamics, Random Walk
In this paper we discuss various aspects of diffusion Monte Carlo methods using a fixed number of walkers. First, a rigorous proof of the divergence of pure diffusion Monte Carlo (PDMC) methods (DMC without branching in which the weights are carried along trajectories) is given. Second, a bias-free Monte Carlo method combining DMC and PDMC approaches, and based on a minimal stochastic reconfiguration of the population, is discussed. Finally, some illustrative calculations for a system of coupled quantum rotators are presented.