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:1609.02703 · Physical review. E

Confidence and efficiency scaling in variational quantum Monte Carlo calculations

Fran�ois Delyon, B. Bernu, Markus Holzmann

Published Feb 17, 2017• 1 citation
arXivHTMLPDF

Abstract

Based on the central limit theorem, we discuss the problem of evaluation of the statistical error of Monte Carlo calculations using a time-discretized diffusion process. We present a robust and practical method to determine the effective variance of general observables and show how to verify the equilibrium hypothesis by the Kolmogorov-Smirnov test. We then derive scaling laws of the efficiency illustrated by variational Monte Carlo calculations on the two-dimensional electron gas.

DOI
10.1103/physreve.95.023307

Related papers

Heavy-tailed random error in quantum Monte Carlo

J. R. Trail

The combination of continuum many-body quantum physics and Monte Carlo methods provide a powerful and well established approach to first principles calculations for large systems. Replacing the exact solution of the problem with a statistical estimate requires a measure of the random error in the estimate for it to be useful. Such a measure of confidence is usually provided by assuming the central limit theorem to hold true. In what follows it is demonstrated that, for the most popular implementation of the variational Monte Carlo method, the central limit theorem has limited validity, or is invalid and must be replaced by a generalized central limit theorem. Estimates of the total energy and the variance of the local energy are examined in detail, and shown to exhibit uncontrolled statistical errors through an explicit derivation of the distribution of the random error. Several examples are given of estimated quantities for which the central limit theorem is not valid. The approach used is generally applicable to characterizing the random error of estimates, and to quantum Monte Carlo methods beyond variational Monte Carlo.

arXiv:0909.5505•Physical Review E•Jan 8, 2008•46 citations

Optimum and efficient sampling for variational quantum Monte Carlo

J. R. Trail, Ryo Maezono

Quantum mechanics for many-body systems may be reduced to the evaluation of integrals in 3N dimensions using Monte Carlo, providing the Quantum Monte Carlo ab initio methods. Here we limit ourselves to expectation values for trial wave functions, that is to variational quantum Monte Carlo. Almost all previous implementations employ samples distributed as the physical probability density of the trial wave function, and assume the central limit theorem to be valid. In this paper we provide an analysis of random error in estimation and optimization that leads naturally to new sampling strategies with improved computational and statistical properties. A rigorous lower limit to the random error is derived, and an efficient sampling strategy presented that significantly increases computational efficiency. In addition the infinite variance heavy tailed random errors of optimum parameters in conventional methods are replaced with a Normal random error, strengthening the theoretical basis of optimization. The method is applied to a number of first row systems and compared with previously published results.

arXiv:1011.4344•The Journal of Chemical Physics•Nov 4, 2010•8 citations

Towards Correlated Sampling for the Fixed-Node Diffusion Quantum Monte Carlo Method

Raphael Berner, René Petz, Arne Lüchow

Most methods of quantum chemistry calculate total energies rather than directly the energy differences that are of interest to chemists. In the case of statistical methods like quantum Monte Carlo the statistical errors in the absolute values need to be considerably smaller than their difference. The calculation of small energy differences is therefore particularly time consuming. Correlated sampling techniques provide the possibility to compute directly energy differences by simulating the underlying systems with the same stochastic process. The smaller the energy difference the smaller its statistical error. Correlated sampling is well established in variational quantum Monte Carlo, but it is much more difficult to implement in diffusion quantum Monte Carlo because of the fixed node approximation. A correlated sampling formalism and a corresponding algorithm based on a transformed Schrödinger equation having the form of a Kolmogorov’s backward equation is derived. The numerical verification of the presented algorithm is given for the harmonic oscillator. The extension of the algorithm to electron structure calculations is discussed.

Zeitschrift für Naturforschung A•Jun 2, 2014•0 citations

Alternative sampling for variational quantum Monte Carlo

J. R. Trail

Expectation values of physical quantities may accurately be obtained by the evaluation of integrals within many-body quantum mechanics, and these multidimensional integrals may be estimated using Monte Carlo methods. In a previous publication it has been shown that for the simplest, most commonly applied strategy in continuum quantum Monte Carlo, the random error in the resulting estimates is not well controlled. At best the central limit theorem is valid in its weakest form, and at worst it is invalid and replaced by an alternative generalized central limit theorem and non-normal random error. In both cases the random error is not controlled. Here we consider a new "residual sampling strategy" that reintroduces the central limit theorem in its strongest form, and provides full control of the random error in estimates. Estimates of the total energy and the variance of the local energy within variational Monte Carlo are considered in detail, and the approach presented may be generalized to expectation values of other operators, and to other variants of the quantum Monte Carlo method.

arXiv:0909.5504•Physical Review E•Jan 8, 2008•18 citations

Strategies for improving the efficiency of quantum Monte Carlo calculations

R. M. Lee, G. J. Conduit, Norbert Nemec +2 more

We describe a number of strategies for minimizing and calculating accurately the statistical uncertainty in quantum Monte Carlo calculations. We investigate the impact of the sampling algorithm on the efficiency of the variational Monte Carlo method. We then propose a technique to maximize the efficiency of the linear extrapolation of diffusion Monte Carlo results to zero time step, finding that a relative time-step ratio of 1:4 is optimal. Finally, we discuss the removal of serial correlation from data sets by reblocking, setting out criteria for the choice of block length and quantifying the effects of the uncertainty in the estimated correlation length.

arXiv:1006.1798•Physical Review E•Jun 27, 2011•53 citations

Quantum Monte Carlo methods

Joseph Marie Thijssen

Introduction In Chapters 1 to 4 we studied methods for solving the Schrödinger equation for many-electron systems. Many of the techniques described there carry over to other quantum many-particle systems, such as liquid helium, and the protons and neutrons in a nucleus. The techniques which we discussed there were, however, all of a mean-field type and therefore correlation effects could not be taken into account without introducing approximations. In this chapter, we consider more accurate techniques, which are similar to those studied in Chapter 10 and are based on using (pseudo-)random numbers – hence the name ‘Monte Carlo’ for these methods. In Chapter 10 we applied Monte Carlo techniques to classical many-particle systems; here we use these techniques for studying quantum problems involving many particles. In the next section we shall see how we can apply Monte Carlo techniques to the problem of calculating the quantum mechanical expectation value of the ground state energy. This is used in order to optimise this expectation value by adjusting a trial wave function in a variational type of approach, hence the name variational Monte Carlo (VMC). In the following section we use the similarity between the Schrödinger equation and the diffusion equation in order to calculate the properties of a collection of interacting quantum mechanical particles by simulating a classical particle diffusion process. The resulting method is called diffusion Monte Carlo (DMC).

Cambridge University Press eBooks•Mar 22, 2007•26 citations
Based on the central limit theorem, we discuss the problem of evaluation of the statistical error of Monte Carlo calculations using a time-discretized diffusion process. We present a robust and practical method to determine the effective variance of general observables and show how to verify the equilibrium hypothesis by the Kolmogorov-Smirnov test. We then derive scaling laws of the efficiency illustrated by variational Monte Carlo calculations on the two-dimensional electron gas.
The combination of continuum many-body quantum physics and Monte Carlo methods provide a powerful and well established approach to first principles calculations for large systems. Replacing the exact solution of the problem with a statistical estimate requires a measure of the random error in the estimate for it to be useful. Such a measure of confidence is usually provided by assuming the central limit theorem to hold true. In what follows it is demonstrated that, for the most popular implementation of the variational Monte Carlo method, the central limit theorem has limited validity, or is invalid and must be replaced by a generalized central limit theorem. Estimates of the total energy and the variance of the local energy are examined in detail, and shown to exhibit uncontrolled statistical errors through an explicit derivation of the distribution of the random error. Several examples are given of estimated quantities for which the central limit theorem is not valid. The approach used is generally applicable to characterizing the random error of estimates, and to quantum Monte Carlo methods beyond variational Monte Carlo.
Quantum mechanics for many-body systems may be reduced to the evaluation of integrals in 3N dimensions using Monte Carlo, providing the Quantum Monte Carlo ab initio methods. Here we limit ourselves to expectation values for trial wave functions, that is to variational quantum Monte Carlo. Almost all previous implementations employ samples distributed as the physical probability density of the trial wave function, and assume the central limit theorem to be valid. In this paper we provide an analysis of random error in estimation and optimization that leads naturally to new sampling strategies with improved computational and statistical properties. A rigorous lower limit to the random error is derived, and an efficient sampling strategy presented that significantly increases computational efficiency. In addition the infinite variance heavy tailed random errors of optimum parameters in conventional methods are replaced with a Normal random error, strengthening the theoretical basis of optimization. The method is applied to a number of first row systems and compared with previously published results.
Most methods of quantum chemistry calculate total energies rather than directly the energy differences that are of interest to chemists. In the case of statistical methods like quantum Monte Carlo the statistical errors in the absolute values need to be considerably smaller than their difference. The calculation of small energy differences is therefore particularly time consuming. Correlated sampling techniques provide the possibility to compute directly energy differences by simulating the underlying systems with the same stochastic process. The smaller the energy difference the smaller its statistical error. Correlated sampling is well established in variational quantum Monte Carlo, but it is much more difficult to implement in diffusion quantum Monte Carlo because of the fixed node approximation. A correlated sampling formalism and a corresponding algorithm based on a transformed Schrödinger equation having the form of a Kolmogorov’s backward equation is derived. The numerical verification of the presented algorithm is given for the harmonic oscillator. The extension of the algorithm to electron structure calculations is discussed.
Expectation values of physical quantities may accurately be obtained by the evaluation of integrals within many-body quantum mechanics, and these multidimensional integrals may be estimated using Monte Carlo methods. In a previous publication it has been shown that for the simplest, most commonly applied strategy in continuum quantum Monte Carlo, the random error in the resulting estimates is not well controlled. At best the central limit theorem is valid in its weakest form, and at worst it is invalid and replaced by an alternative generalized central limit theorem and non-normal random error. In both cases the random error is not controlled. Here we consider a new "residual sampling strategy" that reintroduces the central limit theorem in its strongest form, and provides full control of the random error in estimates. Estimates of the total energy and the variance of the local energy within variational Monte Carlo are considered in detail, and the approach presented may be generalized to expectation values of other operators, and to other variants of the quantum Monte Carlo method.
We describe a number of strategies for minimizing and calculating accurately the statistical uncertainty in quantum Monte Carlo calculations. We investigate the impact of the sampling algorithm on the efficiency of the variational Monte Carlo method. We then propose a technique to maximize the efficiency of the linear extrapolation of diffusion Monte Carlo results to zero time step, finding that a relative time-step ratio of 1:4 is optimal. Finally, we discuss the removal of serial correlation from data sets by reblocking, setting out criteria for the choice of block length and quantifying the effects of the uncertainty in the estimated correlation length.
Introduction In Chapters 1 to 4 we studied methods for solving the Schrödinger equation for many-electron systems. Many of the techniques described there carry over to other quantum many-particle systems, such as liquid helium, and the protons and neutrons in a nucleus. The techniques which we discussed there were, however, all of a mean-field type and therefore correlation effects could not be taken into account without introducing approximations. In this chapter, we consider more accurate techniques, which are similar to those studied in Chapter 10 and are based on using (pseudo-)random numbers – hence the name ‘Monte Carlo’ for these methods. In Chapter 10 we applied Monte Carlo techniques to classical many-particle systems; here we use these techniques for studying quantum problems involving many particles. In the next section we shall see how we can apply Monte Carlo techniques to the problem of calculating the quantum mechanical expectation value of the ground state energy. This is used in order to optimise this expectation value by adjusting a trial wave function in a variational type of approach, hence the name variational Monte Carlo (VMC). In the following section we use the similarity between the Schrödinger equation and the diffusion equation in order to calculate the properties of a collection of interacting quantum mechanical particles by simulating a classical particle diffusion process. The resulting method is called diffusion Monte Carlo (DMC).