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

ACS symposium series

Recent Progress in Fermion Monte Carlo

Malvin H. Kalos, Francesco Pederiva

Published Dec 31, 2006• 1 citation
Publisher

Abstract

Our approach to the solution of the Schrödinger Equation for many-fermion systems has been extensively revised. We have devised a generalization of "acceptance/rejection" that applies to signed random walkers. We have introduced a new class of importance functions for two walkers that better reflects the structure of the enlarged Euclidean space of the pair. For greater flexibility, we no longer rely on the "local energy" of the importance function to determine the dynamics of the walk. We sketch these technical changes and give new results for the two-dimensional electron gas.

DOI
10.1021/bk-2007-0953.ch008

Related papers

Fermion Monte Carlo

Francesco Pederiva, Malvin H. Kalos

We review the fundamental challenge of fermion Monte Carlo for continuous systems, the “sign problem”. We seek that eigenfunction of the many-body Schrödinger equation that is antisymmetric under interchange of the coordinates of pairs of particles. We describe methods that depend upon the use of correlated dynamics for pairs of correlated walkers that carry opposite signs. There is an algorithmic symmetry between such walkers that must be broken to create a method that is both exact and as effective as for symmetric functions. In our new method, it is broken by using different “guiding” functions for walkers of opposite signs, and a geometric correlation between steps of their walks. With a specific process of cancellation of the walkers, overlaps with antisymmetric test functions are preserved. Finally, we describe the progress in treating free-fermion systems and a fermion fluid with 14 3He atoms.

Computer Physics Communications•Sep 1, 1999•0 citations

Exact Monte Carlo Method for Continuum Fermion Systems

Malvin H. Kalos, Francesco Pederiva

We offer a new proposal for the Monte Carlo treatment of many-fermion systems in continuous space. It is based upon diffusion Monte Carlo with significant modifications: correlated pairs of random walkers that carry opposite signs, different functions "guide" walkers of different signs, the Gaussians used for members of a pair are correlated, and walkers can cancel so as to conserve their expected future contributions. We report results for free-fermion systems and a fermion fluid with 14 3He atoms, where it proves stable and correct. Its computational complexity grows with particle number, but slowly enough to make interesting physics within the reach of contemporary computers.

arXiv:cond-mat/0001193•Physical Review Letters•Oct 23, 2000•39 citations

Projector quantum Monte Carlo

Richard M. Martin, Lucia Reining, David M. Ceperley

… as suggested by Fermi, the time-independent Schrödinger equation … can be interpreted as describing the behavior of a system of particles each of which performs a random walk, i.e., diffuses isotropically and at the same time is subject to multiplication, which is determined by the value of the point function V. N. Metropolis and S. Ulam, 1949 Summary In the projector quantum Monte Carlo method, one uses a function of the hamiltonian to sample a distribution proportional to the exact ground-state wavefunction, and thereby computes exact matrix elements of it. An importance sampling transformation makes the algorithm much more efficient. In this chapter we introduce and develop the diffusion Monte Carlo method, which involves drifting, branching random walks. For any excited state, including any system with more than two electrons, one encounters the sign problem, limiting the direct application of these algorithms for most fermion systems. Instead, by using approximate fixed-node or fixed-phase boundary conditions, one can achieve efficiency similar to variational Monte Carlo. We also discuss the application of the projector method in a basis of Slater determinants. In this chapter, we discuss a different quantum Monte Carlo method, projector Monte Carlo (PMC). This general method was first suggested by Fermi [1049]; see the quote at the start of this chapter by two of the inventors of the Monte Carlo method. An implementation of PMC was tried out in the early days of computing [1050]. Advances in methodology, in particular importance sampling, resulted in a significant large-scale application: the exact calculation of the ground-state properties of 256 hard-sphere bosons by Kalos, Levesque, and Verlet [1051] in 1974. Calculations for electronic systems and the fixednode approximation were introduced by Anderson [1052, 1053]. One of the most important projector MC algorithms, the diffusion Monte Carlo algorithm with importance sampling for fermions, was used to compute the correlation energy of homogeneous electron gas by Ceperley and Alder [109] in 1980; the resulting HEG correlation energy was crucial in the development of density functional calculations. Types and properties of projectors In this method, a many-body projector G ( R , R ) = Ĝ is repeatedly applied to filter out the exact many-body ground state from an initial state; the operation of the projector is carried out with a random walk, hence the name of this class of methods.

Cambridge University Press eBooks•May 31, 2016•0 citations

Correlated pairs in model fermion Monte Carlo calculations

Malvin H. Kalos

The issues that prevent the development of efficient and stable algorithms for fermion Monte Carlo calculations in continuum systems are reexamined with special reference to the implications of the ``plus-minus'' symmetry. This is a property of many algorithms that use signed walkers, namely, that the dynamics are unchanged when the signs of the walkers are interchanged. Algorithms that obey this symmetry cannot exhibit the necessary stability. Specifically, estimates of the overlap with any antisymmetric test function cannot be bounded away from zero in the limit of many iterations. Within the framework of a diffusion Monte Carlo treatment of the Schr\"odinger equation, it is shown that this symmetry is easily broken for pairs of walkers while at the same time preserving the correct marginal dynamics for each member of the pair. The key is to create different classes of correlations between members of pairs and to use (at least) two distinct correlations for a pair and for the same pair with signs exchanged. The ideas are applied successfully for a class of simple model problems in two dimensions. \textcopyright{} 1996 The American Physical Society.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•May 1, 1996•12 citations

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

Alexander A. Kunitsa, So Hirata

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.

arXiv:1811.03000•Physical review. E•Jan 27, 2020•12 citations

Random Walk Beyond Hartree-Fock

James E. Gubernatis, M. Guerrero

We give a brief discussion of the recently developed Constrained-Path Monte Carlo Method. This method is a quantum Monte Carlo technique that eliminates the fermion sign problem plaguing simulations of systems of interacting electrons. The elimination is accomplished by trading an exact procedure for an approximate one that has been demonstrated to give very accurate estimates of energies and manybody correlation functions. We also give a short review of its applications, a discussion of several strategies for parallelizing it, and some speculation of its future extensions. 1

•Jan 1, 2000•0 citations
Our approach to the solution of the Schrödinger Equation for many-fermion systems has been extensively revised. We have devised a generalization of "acceptance/rejection" that applies to signed random walkers. We have introduced a new class of importance functions for two walkers that better reflects the structure of the enlarged Euclidean space of the pair. For greater flexibility, we no longer rely on the "local energy" of the importance function to determine the dynamics of the walk. We sketch these technical changes and give new results for the two-dimensional electron gas.
We review the fundamental challenge of fermion Monte Carlo for continuous systems, the “sign problem”. We seek that eigenfunction of the many-body Schrödinger equation that is antisymmetric under interchange of the coordinates of pairs of particles. We describe methods that depend upon the use of correlated dynamics for pairs of correlated walkers that carry opposite signs. There is an algorithmic symmetry between such walkers that must be broken to create a method that is both exact and as effective as for symmetric functions. In our new method, it is broken by using different “guiding” functions for walkers of opposite signs, and a geometric correlation between steps of their walks. With a specific process of cancellation of the walkers, overlaps with antisymmetric test functions are preserved. Finally, we describe the progress in treating free-fermion systems and a fermion fluid with 14 3He atoms.
We offer a new proposal for the Monte Carlo treatment of many-fermion systems in continuous space. It is based upon diffusion Monte Carlo with significant modifications: correlated pairs of random walkers that carry opposite signs, different functions "guide" walkers of different signs, the Gaussians used for members of a pair are correlated, and walkers can cancel so as to conserve their expected future contributions. We report results for free-fermion systems and a fermion fluid with 14 3He atoms, where it proves stable and correct. Its computational complexity grows with particle number, but slowly enough to make interesting physics within the reach of contemporary computers.
… as suggested by Fermi, the time-independent Schrödinger equation … can be interpreted as describing the behavior of a system of particles each of which performs a random walk, i.e., diffuses isotropically and at the same time is subject to multiplication, which is determined by the value of the point function V. N. Metropolis and S. Ulam, 1949 Summary In the projector quantum Monte Carlo method, one uses a function of the hamiltonian to sample a distribution proportional to the exact ground-state wavefunction, and thereby computes exact matrix elements of it. An importance sampling transformation makes the algorithm much more efficient. In this chapter we introduce and develop the diffusion Monte Carlo method, which involves drifting, branching random walks. For any excited state, including any system with more than two electrons, one encounters the sign problem, limiting the direct application of these algorithms for most fermion systems. Instead, by using approximate fixed-node or fixed-phase boundary conditions, one can achieve efficiency similar to variational Monte Carlo. We also discuss the application of the projector method in a basis of Slater determinants. In this chapter, we discuss a different quantum Monte Carlo method, projector Monte Carlo (PMC). This general method was first suggested by Fermi [1049]; see the quote at the start of this chapter by two of the inventors of the Monte Carlo method. An implementation of PMC was tried out in the early days of computing [1050]. Advances in methodology, in particular importance sampling, resulted in a significant large-scale application: the exact calculation of the ground-state properties of 256 hard-sphere bosons by Kalos, Levesque, and Verlet [1051] in 1974. Calculations for electronic systems and the fixednode approximation were introduced by Anderson [1052, 1053]. One of the most important projector MC algorithms, the diffusion Monte Carlo algorithm with importance sampling for fermions, was used to compute the correlation energy of homogeneous electron gas by Ceperley and Alder [109] in 1980; the resulting HEG correlation energy was crucial in the development of density functional calculations. Types and properties of projectors In this method, a many-body projector G ( R , R ) = Ĝ is repeatedly applied to filter out the exact many-body ground state from an initial state; the operation of the projector is carried out with a random walk, hence the name of this class of methods.
The issues that prevent the development of efficient and stable algorithms for fermion Monte Carlo calculations in continuum systems are reexamined with special reference to the implications of the ``plus-minus'' symmetry. This is a property of many algorithms that use signed walkers, namely, that the dynamics are unchanged when the signs of the walkers are interchanged. Algorithms that obey this symmetry cannot exhibit the necessary stability. Specifically, estimates of the overlap with any antisymmetric test function cannot be bounded away from zero in the limit of many iterations. Within the framework of a diffusion Monte Carlo treatment of the Schr\"odinger equation, it is shown that this symmetry is easily broken for pairs of walkers while at the same time preserving the correct marginal dynamics for each member of the pair. The key is to create different classes of correlations between members of pairs and to use (at least) two distinct correlations for a pair and for the same pair with signs exchanged. The ideas are applied successfully for a class of simple model problems in two dimensions. \textcopyright{} 1996 The American Physical Society.
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.
We give a brief discussion of the recently developed Constrained-Path Monte Carlo Method. This method is a quantum Monte Carlo technique that eliminates the fermion sign problem plaguing simulations of systems of interacting electrons. The elimination is accomplished by trading an exact procedure for an approximate one that has been demonstrated to give very accurate estimates of energies and manybody correlation functions. We also give a short review of its applications, a discussion of several strategies for parallelizing it, and some speculation of its future extensions. 1