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

Anwendung und Weiterentwicklung der Diffusions-Quanten-Monte-Carlo-Methode

Raphaël Berner, Arne Lüchow

Published Jan 1, 2010• 0 citations
OpenAlex

Abstract

The diffusion quantum Monte Carlo method (DMC) provides a numeric solution of the time dependent Schrödinger equation by means of a random walk process. To fulfil the pauli principle the nodes of a so called trial wavefunction, obtained with standard ab initio quantum chemistry methods like HF, DFT, or CASSCF, are considered as additional boundary conditions. In case of such nodes coincide with those of the exact wavefunction DMC yields the exact solution of the Schrödinger equation. Generally, this is not true and one ends up with an systematic bias, referred to as the node location error. Theoretical chemistry aims the prediction of reaction, excitation, or activation energies, i.e. energy differences. Therefore, the node location error has to be seen under the light of efficient error cancellation. One can distinguish dynamic and non-dynamic electron correlation. The first one accounts for the instantaneous, short distance interaction among the electrons. The non-dynamic correlation addresses the interaction between different states of the same symmetry. Both kinds of correlation have a significant impact on the nodal structure and therefore govern the node location error. In this thesis the isomerization of bicylo[1.1.0]butane has been studied by exploiting the DMC method. The underlying reaction mechanism poses a great challenge to quantum chemistry because it comprises species with different amounts of non-dynamical correlation. Indeed, very sophisticated methods like for example CCSD(T) fail in giving a proper description. For the DMC method, the construction of the trial function is of crucial importance. Hence, the latter is generated by means of CASSCF calculations. With regard to efficient error cancellation it is most reasonable to take only those CSFs into account which contribute significantly to the proper description of the non-dynamic correlation. Since the DMC method is variational for the fermionic ground state, the incorporation of such a CSFs leads to a decrease of the energy compared with the output of a computation just based on the reference determinant. In the present case, only up to three CSFs very necessary. As a consequence, the direct optimization of the DMC energy with respect to the CASSCF coefficients is possible. By means of this construction scheme, dynamic and non-dynamic electron correlation effects are accounted for in equal measure for every molecule comprising the reaction path. Hence, the error cancellation is excellent. The obtained reaction, and activation energies are in very good agreement with experimental data and completely renormalized coupled cluster calculations. The presented CSF selection scheme yields very short linear combinations of CSFs and therefore should be applicable even for larger organic molecules. As it has already pointed out, quantum chemistry aims the prediction of energy differences rather than absolute values. A very efficient approach to calculate energy differences within a Monte Carlo framework is correlated sampling. The main idea behind this procedure is to describe the two different systems whose energy difference is to be calculated by means of the same random walk process. The resulting statistical correlation gives rise to a reduction of the statistical uncertainty of the energy difference. As far as we know, until now no implementation of the correlated sampling formalism into the DMC method has been reported in the literature. To implement correlated sampling in the DMC method a mathematical formalism based on the theory of stochastic differential equations has been devised. Based on this formalism an appropriate algorithm has been implemented for test systems like the particle in a one dimensional box, and the harmonic oszillator. Correlated sampling is very well suited for the computation of small energy differences with considerable small statistical uncertainties. A possible application would be the calculation of the weak interaction between DNA base pairs. For this purpose the corresponding formalism has to be implemented in our quantum Monte Carlo program Amolqc. In doing so, a major task is the implementation of routines which account for the propagation with respect to a positive guide function what has already been finished in the course of this work.

Related papers

Calculation of Intermolecular Interactions via Diffusion Monte Carlo

Kevin E. Gasperich

Quantum chemistry is a useful tool that provides insight into the properties and behavior of chemical systems. Modern software packages have made quantum chemistry methods more easily accessible, and the continued increase in available computational resources has allowed them to be applied to larger systems at higher levels of theory. Two significant problems that the field faces are the high computational complexity of high-level methods and the shift toward parallelism in high performance computing architectures. This work examines the treatment of weakly interacting molecular systems with the fixed-node diffusion Monte Carlo (DMC) method. DMC and other quantum Monte Carlo (QMC) methods offer a possible solution to both of the aforementioned problems: they can produce near-exact results with a lower scaling (with respect to problem size) than other similarly-accurate methods, and they are inherently parallel, so there is little additional cost associated with distributing the work of a single QMC calculation across a large number of processing units. The only error in DMC that is not systematically improvable is the constraint of a fixed nodal surface of the many-particle wave function of the system being studied. There are many cases in which a single Slater determinant trial function is sufficient to obtain accurate results, but there are others in which more sophisticated multi-determinant trial functions are necessary. Furthermore, it is non-trivial to generate nodal surfaces of similar quality for isolated and interacting molecules, so cancellation of errors is not guaranteed. We examine the use of different single- and multi-determinant trial functions in DMC calculations on small chemical systems with the goal of further understanding how to construct appropriate trial functions for molecules and clusters.

•Jan 16, 2020•0 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.3417•Physical Review B•May 19, 2009•30 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.06006•The Journal of Chemical Physics•Dec 28, 2024•9 citations

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

The quantum Monte Carlo method—electron correlation from random numbers (abstract only)

Richard J. Needs

The fixed-node diffusion quantum Monte Carlo (DMC) method is the most accurate method known for calculating the energies of large many-particle quantum systems. The key element of the method is the development of accurate trial many-body wavefunctions which control the statistical efficiency of the calculations and the accuracy obtained. Accurate wavefunctions can be obtained by building correlation effects on top of mean field descriptions such as density functional theory. The wavefunctions can be improved by introducing multi-determinants, pairing functions, and backflow transformations. The calculations are expensive, but the method scales well with system size and calculations on 1000 particles are possible. Some recent applications of the DMC method to atoms, molecules and solids will be presented.

Journal of Physics Condensed Matter•Jan 24, 2008•1 citation

A Quantum Monte Carlo Study of Quantum Materials in the Intermediate Correlation Regime

Clio Johnson

Variational and diffusion Monte Carlo (VMC and DMC) are widely used techniques in the Monte Carlo study of many-body quantum systems. Through efficient use of ran- dom sampling to solve many-dimensional integrals, they can be used to reliably solve the Schrödinger equation. Within the wider ecosystems of electronic structure calculation, quantum chemistry, and quantum simulation, quantum Monte Carlo (QMC) methods are used in tandem with density functional theory (DFT) and other computational methods to deliver insight into a vast array of condensed matter systems. Of particular interest to this work are systems in the intermediate correlation regime, where methods such as simple band theory and the Hubbard model break down. This thesis proposes and tests an explicit three-body backflow term to improve the accuracy of VMC and DMC calculations by incorporating higher-order correlation in order to reduce fixed-node er- ror in DMC, finding statistically significant improvements in ground state energies and improved variance in a water molecule. Then, QMC is used to calculate accurate charge densities for weak and intermediately correlated materials, to test commonly used dens- ity functional approximations in a controlled environment, providing insight into where these approximations perform well, and where they fail. Finally, DMC is used to generate ground state energy data for a dipolar Fermi gas of dressed Rydberg atoms interacting in the intermediate correlation regime. This data is parametrised in terms of interaction strength, for use in generating exchange-correlation functionals, to enable the DFT study of inhomogeneous systems of dressed Rydberg atoms.

•Sep 14, 2026•0 citations
The diffusion quantum Monte Carlo method (DMC) provides a numeric solution of the time dependent Schrödinger equation by means of a random walk process. To fulfil the pauli principle the nodes of a so called trial wavefunction, obtained with standard ab initio quantum chemistry methods like HF, DFT, or CASSCF, are considered as additional boundary conditions. In case of such nodes coincide with those of the exact wavefunction DMC yields the exact solution of the Schrödinger equation. Generally, this is not true and one ends up with an systematic bias, referred to as the node location error. Theoretical chemistry aims the prediction of reaction, excitation, or activation energies, i.e. energy differences. Therefore, the node location error has to be seen under the light of efficient error cancellation. One can distinguish dynamic and non-dynamic electron correlation. The first one accounts for the instantaneous, short distance interaction among the electrons. The non-dynamic correlation addresses the interaction between different states of the same symmetry. Both kinds of correlation have a significant impact on the nodal structure and therefore govern the node location error. In this thesis the isomerization of bicylo[1.1.0]butane has been studied by exploiting the DMC method. The underlying reaction mechanism poses a great challenge to quantum chemistry because it comprises species with different amounts of non-dynamical correlation. Indeed, very sophisticated methods like for example CCSD(T) fail in giving a proper description. For the DMC method, the construction of the trial function is of crucial importance. Hence, the latter is generated by means of CASSCF calculations. With regard to efficient error cancellation it is most reasonable to take only those CSFs into account which contribute significantly to the proper description of the non-dynamic correlation. Since the DMC method is variational for the fermionic ground state, the incorporation of such a CSFs leads to a decrease of the energy compared with the output of a computation just based on the reference determinant. In the present case, only up to three CSFs very necessary. As a consequence, the direct optimization of the DMC energy with respect to the CASSCF coefficients is possible. By means of this construction scheme, dynamic and non-dynamic electron correlation effects are accounted for in equal measure for every molecule comprising the reaction path. Hence, the error cancellation is excellent. The obtained reaction, and activation energies are in very good agreement with experimental data and completely renormalized coupled cluster calculations. The presented CSF selection scheme yields very short linear combinations of CSFs and therefore should be applicable even for larger organic molecules. As it has already pointed out, quantum chemistry aims the prediction of energy differences rather than absolute values. A very efficient approach to calculate energy differences within a Monte Carlo framework is correlated sampling. The main idea behind this procedure is to describe the two different systems whose energy difference is to be calculated by means of the same random walk process. The resulting statistical correlation gives rise to a reduction of the statistical uncertainty of the energy difference. As far as we know, until now no implementation of the correlated sampling formalism into the DMC method has been reported in the literature. To implement correlated sampling in the DMC method a mathematical formalism based on the theory of stochastic differential equations has been devised. Based on this formalism an appropriate algorithm has been implemented for test systems like the particle in a one dimensional box, and the harmonic oszillator. Correlated sampling is very well suited for the computation of small energy differences with considerable small statistical uncertainties. A possible application would be the calculation of the weak interaction between DNA base pairs. For this purpose the corresponding formalism has to be implemented in our quantum Monte Carlo program Amolqc. In doing so, a major task is the implementation of routines which account for the propagation with respect to a positive guide function what has already been finished in the course of this work.
Quantum chemistry is a useful tool that provides insight into the properties and behavior of chemical systems. Modern software packages have made quantum chemistry methods more easily accessible, and the continued increase in available computational resources has allowed them to be applied to larger systems at higher levels of theory. Two significant problems that the field faces are the high computational complexity of high-level methods and the shift toward parallelism in high performance computing architectures. This work examines the treatment of weakly interacting molecular systems with the fixed-node diffusion Monte Carlo (DMC) method. DMC and other quantum Monte Carlo (QMC) methods offer a possible solution to both of the aforementioned problems: they can produce near-exact results with a lower scaling (with respect to problem size) than other similarly-accurate methods, and they are inherently parallel, so there is little additional cost associated with distributing the work of a single QMC calculation across a large number of processing units. The only error in DMC that is not systematically improvable is the constraint of a fixed nodal surface of the many-particle wave function of the system being studied. There are many cases in which a single Slater determinant trial function is sufficient to obtain accurate results, but there are others in which more sophisticated multi-determinant trial functions are necessary. Furthermore, it is non-trivial to generate nodal surfaces of similar quality for isolated and interacting molecules, so cancellation of errors is not guaranteed. We examine the use of different single- and multi-determinant trial functions in DMC calculations on small chemical systems with the goal of further understanding how to construct appropriate trial functions for molecules and clusters.
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.
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.
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.
The fixed-node diffusion quantum Monte Carlo (DMC) method is the most accurate method known for calculating the energies of large many-particle quantum systems. The key element of the method is the development of accurate trial many-body wavefunctions which control the statistical efficiency of the calculations and the accuracy obtained. Accurate wavefunctions can be obtained by building correlation effects on top of mean field descriptions such as density functional theory. The wavefunctions can be improved by introducing multi-determinants, pairing functions, and backflow transformations. The calculations are expensive, but the method scales well with system size and calculations on 1000 particles are possible. Some recent applications of the DMC method to atoms, molecules and solids will be presented.
Variational and diffusion Monte Carlo (VMC and DMC) are widely used techniques in the Monte Carlo study of many-body quantum systems. Through efficient use of ran- dom sampling to solve many-dimensional integrals, they can be used to reliably solve the Schrödinger equation. Within the wider ecosystems of electronic structure calculation, quantum chemistry, and quantum simulation, quantum Monte Carlo (QMC) methods are used in tandem with density functional theory (DFT) and other computational methods to deliver insight into a vast array of condensed matter systems. Of particular interest to this work are systems in the intermediate correlation regime, where methods such as simple band theory and the Hubbard model break down. This thesis proposes and tests an explicit three-body backflow term to improve the accuracy of VMC and DMC calculations by incorporating higher-order correlation in order to reduce fixed-node er- ror in DMC, finding statistically significant improvements in ground state energies and improved variance in a water molecule. Then, QMC is used to calculate accurate charge densities for weak and intermediately correlated materials, to test commonly used dens- ity functional approximations in a controlled environment, providing insight into where these approximations perform well, and where they fail. Finally, DMC is used to generate ground state energy data for a dipolar Fermi gas of dressed Rydberg atoms interacting in the intermediate correlation regime. This data is parametrised in terms of interaction strength, for use in generating exchange-correlation functionals, to enable the DFT study of inhomogeneous systems of dressed Rydberg atoms.