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:cond-mat/0703004 · Physical Review Letters

Variational Wave Functions and Their Overlap with the Ground State

Christophe Mora, Xavier Waintal

Published Jul 20, 2007• 18 citations
arXivHTMLPDF

Abstract

An intrinsic measure of the quality of a variational wave function is given by its overlap with the ground state of the system. We derive a general formula to compute this overlap when quantum dynamics in imaginary time is accessible. The overlap is simply related to the area under the E(tau) curve, i.e., the energy as a function of imaginary time. This has important applications to, for example, quantum Monte Carlo simulations where the overlap becomes as a simple by-product of routine simulations. As a result, we find that the practical definition of a good variational wave function for quantum Monte Carlo simulations, i.e., fast convergence to the ground state, is equivalent to a good overlap with the actual ground state of the system.

DOI
10.1103/physrevlett.99.030403

Related papers

How to recognize a good variational wave function

Christophe Mora, Xavier Waintal

We propose a practical skim to calculate the overlap between the true ground state of the system and a variational ansatz within a quantum Monte-Carlo simulation. This overlap is simply related to the area under the $E(\tau)$ curve (energy as a function of imaginary time), and hence comes as a simple byproduct of usual simulations. As a result, we find that the practical definition of a good variational wave function for quantum Monte-Carlo simulations, ``it converges fast to the ground state'', is equivalent to having a good overlap with the actual ground state of the system.

•Feb 28, 2007•0 citations

Reptation Quantum Monte Carlo: A Method for Unbiased Ground-State Averages and Imaginary-Time Correlations

Stefano Baroni A, Saverio Moroni

We introduce a new stochastic method for calculating ground-state properties of quantum systems. Segments of a Langevin random walk guided by a trial wave function are subject to a Metropolis rejection test performed on the time integral of the local energy. The algorithm---which is as simple as variational Monte Carlo---for bosons provides exact expectation values of local observables, as well as their static and dynamic (in imaginary time) response functions, without mixed-estimate nor population-control biases. Our method is demonstrated with a few case applications to ${}^{4}\mathrm{He}$.

Physical Review Letters•Jun 14, 1999•251 citations

Variational Principles in Quantum Monte Carlo: The Troubled Story of Variance Minimization

Alice Cuzzocrea, Anthony Scemama, Wim J. Briels +2 more

We investigate the use of different variational principles in quantum Monte Carlo, namely, energy and variance minimization, prompted by the interest in the robust and accurate estimation of electronic excited states. For two prototypical, challenging molecules, we readily reach the accuracy of the best available reference excitation energies using energy minimization in a state-specific or state-average fashion for states of different or equal symmetry, respectively. On the other hand, in variance minimization, where the use of suitable functionals is expected to target specific states regardless of the symmetry, we encounter severe problems for a variety of wave functions: as the variance converges, the energy drifts away from that of the selected state. This unexpected behavior is sometimes observed even when the target is the ground state and generally prevents the robust estimation of total and excitation energies. We analyze this problem using a very simple wave function and infer that the optimization finds little or no barrier to escape from a local minimum or local plateau, eventually converging to a lower-variance state instead of the target state. For the increasingly complex systems becoming in reach of quantum Monte Carlo simulations, variance minimization with current functionals appears to be an impractical route.

arXiv:2002.05951•Journal of Chemical Theory and Computation•May 17, 2020•38 citations

Exact ground state Monte Carlo method for Bosons without importance sampling

M. Rossi, Marco Nava, L. Reatto +1 more

Generally "exact" quantum Monte Carlo computations for the ground state of many bosons make use of importance sampling. The importance sampling is based either on a guiding function or on an initial variational wave function. Here we investigate the need of importance sampling in the case of path integral ground state (PIGS) Monte Carlo. PIGS is based on a discrete imaginary time evolution of an initial wave function with a nonzero overlap with the ground state, which gives rise to a discrete path which is sampled via a Metropolis-like algorithm. In principle the exact ground state is reached in the limit of an infinite imaginary time evolution, but actual computations are based on finite time evolutions and the question is whether such computations give unbiased exact results. We have studied bulk liquid and solid (4)He with PIGS by considering as initial wave function a constant, i.e., the ground state of an ideal Bose gas. This implies that the evolution toward the ground state is driven only by the imaginary time propagator, i.e., there is no importance sampling. For both phases we obtain results converging to those obtained by considering the best available variational wave function (the shadow wave function) as initial wave function. Moreover we obtain the same results even by considering wave functions with the wrong correlations, for instance, a wave function of a strongly localized Einstein crystal for the liquid phase. This convergence is true not only for diagonal properties such as the energy, the radial distribution function, and the static structure factor, but also for off-diagonal ones, such as the one-body density matrix. This robustness of PIGS can be traced back to the fact that the chosen initial wave function acts only at the beginning of the path without affecting the imaginary time propagator. From this analysis we conclude that zero temperature PIGS calculations can be as unbiased as those of finite temperature path integral Monte Carlo. On the other hand, a judicious choice of the initial wave function greatly improves the rate of convergence to the exact results.

arXiv:0907.4430•The Journal of Chemical Physics•Oct 19, 2009•62 citations

Lower bounds to the overlap

Thomas Hoffmann-Ostenhof, Maria Hoffmann-Ostenhof, Gottfried Olbrich

A variational lower bound to the overlap of the ground state psi 0 of a quantum mechanical system with an approximate wavefunction phi is given. A simplified non-variational formula is deduced from this lower bound to mod ( phi mod psi 0 ) mod . Both versions are illustrated with an application to the hydrogen atom and are compared with the lower bounds of Eckart and Wang. The non-variational version turns out to be a useful alternative to Eckart's bound.

Journal of Physics A Mathematical and General•Jan 1, 1976•6 citations

Nonlinear density response from imaginary-time correlation functions: Ab initio path integral Monte Carlo simulations of the warm dense electron gas

Tobias Dornheim, Zhandos A. Moldabekov, Jan Vorberger

The ab initio path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCFs). For example, the well-known density–density ITCF F(q, τ) allows one to estimate the linear response of a given system for all wave vectors q from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor S(q, ω)—a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in the quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor S(q1, ω1; q2, ω2), which evokes the enticing possibility to study dynamic three-body effects on an ab initio level.

arXiv:2106.01619•The Journal of Chemical Physics•Aug 4, 2021•55 citations
An intrinsic measure of the quality of a variational wave function is given by its overlap with the ground state of the system. We derive a general formula to compute this overlap when quantum dynamics in imaginary time is accessible. The overlap is simply related to the area under the E(tau) curve, i.e., the energy as a function of imaginary time. This has important applications to, for example, quantum Monte Carlo simulations where the overlap becomes as a simple by-product of routine simulations. As a result, we find that the practical definition of a good variational wave function for quantum Monte Carlo simulations, i.e., fast convergence to the ground state, is equivalent to a good overlap with the actual ground state of the system.
We propose a practical skim to calculate the overlap between the true ground state of the system and a variational ansatz within a quantum Monte-Carlo simulation. This overlap is simply related to the area under the E(τ)E(\tau)E(τ) curve (energy as a function of imaginary time), and hence comes as a simple byproduct of usual simulations. As a result, we find that the practical definition of a good variational wave function for quantum Monte-Carlo simulations, ``it converges fast to the ground state'', is equivalent to having a good overlap with the actual ground state of the system.
We introduce a new stochastic method for calculating ground-state properties of quantum systems. Segments of a Langevin random walk guided by a trial wave function are subject to a Metropolis rejection test performed on the time integral of the local energy. The algorithm---which is as simple as variational Monte Carlo---for bosons provides exact expectation values of local observables, as well as their static and dynamic (in imaginary time) response functions, without mixed-estimate nor population-control biases. Our method is demonstrated with a few case applications to 4He{}^{4}\mathrm{He}4He.
We investigate the use of different variational principles in quantum Monte Carlo, namely, energy and variance minimization, prompted by the interest in the robust and accurate estimation of electronic excited states. For two prototypical, challenging molecules, we readily reach the accuracy of the best available reference excitation energies using energy minimization in a state-specific or state-average fashion for states of different or equal symmetry, respectively. On the other hand, in variance minimization, where the use of suitable functionals is expected to target specific states regardless of the symmetry, we encounter severe problems for a variety of wave functions: as the variance converges, the energy drifts away from that of the selected state. This unexpected behavior is sometimes observed even when the target is the ground state and generally prevents the robust estimation of total and excitation energies. We analyze this problem using a very simple wave function and infer that the optimization finds little or no barrier to escape from a local minimum or local plateau, eventually converging to a lower-variance state instead of the target state. For the increasingly complex systems becoming in reach of quantum Monte Carlo simulations, variance minimization with current functionals appears to be an impractical route.
Generally "exact" quantum Monte Carlo computations for the ground state of many bosons make use of importance sampling. The importance sampling is based either on a guiding function or on an initial variational wave function. Here we investigate the need of importance sampling in the case of path integral ground state (PIGS) Monte Carlo. PIGS is based on a discrete imaginary time evolution of an initial wave function with a nonzero overlap with the ground state, which gives rise to a discrete path which is sampled via a Metropolis-like algorithm. In principle the exact ground state is reached in the limit of an infinite imaginary time evolution, but actual computations are based on finite time evolutions and the question is whether such computations give unbiased exact results. We have studied bulk liquid and solid (4)He with PIGS by considering as initial wave function a constant, i.e., the ground state of an ideal Bose gas. This implies that the evolution toward the ground state is driven only by the imaginary time propagator, i.e., there is no importance sampling. For both phases we obtain results converging to those obtained by considering the best available variational wave function (the shadow wave function) as initial wave function. Moreover we obtain the same results even by considering wave functions with the wrong correlations, for instance, a wave function of a strongly localized Einstein crystal for the liquid phase. This convergence is true not only for diagonal properties such as the energy, the radial distribution function, and the static structure factor, but also for off-diagonal ones, such as the one-body density matrix. This robustness of PIGS can be traced back to the fact that the chosen initial wave function acts only at the beginning of the path without affecting the imaginary time propagator. From this analysis we conclude that zero temperature PIGS calculations can be as unbiased as those of finite temperature path integral Monte Carlo. On the other hand, a judicious choice of the initial wave function greatly improves the rate of convergence to the exact results.
The ab initio path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCFs). For example, the well-known density–density ITCF F(q, τ) allows one to estimate the linear response of a given system for all wave vectors q from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor S(q, ω)—a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in the quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor S(q1, ω1; q2, ω2), which evokes the enticing possibility to study dynamic three-body effects on an ab initio level.