A method for accelerating the rate of convergence of the long-time (or small-temperatur) limit of quantum Monte Carlo approaches is presented. To do that, a variation of the Lanczós algorithm suitable for QMC data is introduced. This algorithm allows one to extract more information from correlation functions at small times, thus avoiding large statistical fluctuations associated with large times. It is first applied to an exactly soluble system and then to the LiH molecule. Calculations using both the fixed-node and nodal-release approaches are discussed.
A method for accelerating the rate of convergence of the long-time (or small-temperatur) limit of quantum Monte Carlo approaches is presented. To do that, a variation of the Lanczós algorithm suitable for QMC data is introduced. This algorithm allows one to extract more information from correlation functions at small times, thus avoiding large statistical fluctuations associated with large times. It is first applied to an exactly soluble system and then to the LiH molecule. Calculations using both the fixed-node and nodal-release approaches are discussed.
We show that the standard Lanczos algorithm can be efficiently implemented statistically and self-consistently improved, using the stochastic reconfiguration method, which has been recently introduced to stabilize the Monte Carlo sign problem instability. With this scheme a few Lanczos steps over a given variational wave function are possible even for large size as a particular case of a more general and more accurate technique that allows to obtain lower variational energies. This method has been tested extensively for a strongly correlated model like the $t\ensuremath{-}J$ model. With the standard Lanczos technique it is possible to compute any kind of correlation functions, with no particular computational effort. By using the fact that the variance $〈{H}^{2}〉\ensuremath{-}〈H{〉}^{2}$ is zero for an exact eigenstate, we show that the approach to the exact solution with few Lanczos iterations is indeed possible even for $\ensuremath{\sim}100$ electrons for reasonably good initial wave functions. The variational stochastic reconfiguration technique presented here allows in general a many-parameter energy optimization of any computable many-body wave function, including for instance generic long-range Jastrow factors and arbitrary site-dependent orbital determinants. This scheme improves further the accuracy of the calculation, especially for long-distance correlation functions.
Electronic structure quantum Monte Carlo Quantum Monte Carlo (QMC) is an advanced simulation methodology for studies of manybody quantum systems. The QMC approaches combine analytical insights with stochastic computational techniques for efficient solution of several classes of important many-body problems such as the stationary Schrödinger equation. QMC methods of various flavors have been applied to a great variety of systems spanning continuous and lattice quantum models, molecular and condensed systems, BEC-BCS ultracold condensates, nuclei, etc. In this review, we focus on the electronic structure QMC, i.e., methods relevant for systems described by the electron-ion Hamiltonians. Some of the key QMC achievements include direct treatment of electron correlation, accuracy in predicting energy differences and favorable scaling in the system size. Calculations of atoms, molecules, clusters and solids have demonstrated QMC applicability to real systems with hundreds of electrons while providing 90-95% of the correlation energy and energy differences typically within a few percent of experiments. Advances in accuracy beyond these limits are hampered by the so-called fixed-node approximation which is used to circumvent the notorious fermion sign problem. Many-body nodes of fermion states and their properties have therefore become one of the important topics for further progress in predictive power and efficiency of QMC calculations. Some of our recent results on the wave function nodes and related nodal domain topologies will be briefly reviewed. This includes analysis of few-electron systems and descriptions of exact and approximate nodes using transformations and projections of the highly-dimensional nodal hypersurfaces into the 3D space. Studies of fermion nodes offer new insights into topological properties of eigenstates such as explicit demonstrations that generic fermionic ground states exhibit the minimal number of two nodal domains. Recently proposed trial wave functions based on Pfaffians with pairing orbitals are presented and their nodal properties are tested in calculations of first row atoms and molecules. Finally, backflow "dressed" coordinates are introduced as another possibility for capturing correlation effects and for decreasing the fixed-node bias.
The Z6 scaling of computational effort in quantum Monte Carlo (QMC) results in part from a type of critical slowing down. This insight leads to an attack based on acceleration algorithms which modify the dynamics in Langevin-based simulations without changing the steady-state. Configuration space is more rapidly sampled and steps are inherently more decorrelated. A second approach uses a hybrid Metropolis–Green function Monte Carlo method which factors the core and valence wave function, and treats the former variationally with accurate trial functions, and the latter exactly in the presence of the variationally evoling core. This allows separation of the relevant core and valence time scales, and has sped up quantum simulations by over three orders of magnitude.
Quantum Monte Carlo (QMC) methods have received considerable attention over past decades due to their great promise for providing a direct solution to the many-body Schrodinger equation in electronic systems. Thanks to their low scaling with the number of particles, QMC methods present a compelling competitive alternative for the accurate study of large molecular systems and solid state calculations. In spite of such promise, the method has not permeated the quantum chemistry community broadly, mainly because of the fixed-node error, which can be large and whose control is difficult. In this Perspective, we present a systematic application of large scale multideterminant expansions in QMC and report on its impressive performance with first row dimers and the 55 molecules of the G1 test set. We demonstrate the potential of this strategy for systematically reducing the fixed-node error in the wave function and for achieving chemical accuracy in energy predictions. When compared to traditional quantum chemistry methods like MP2, CCSD(T), and various DFT approximations, the QMC results show a marked improvement over all of them. In fact, only the explicitly correlated CCSD(T) method with a large basis set produces more accurate results. Further developments in trial wave functions and algorithmic improvements appear promising for rendering QMC as the benchmark standard in large electronic systems.
Quantum Monte Carlo (QMC) methods have received considerable attention over the last decades due to their great promise for providing a direct solution to the many-body Schrodinger equation in electronic systems.Thanks to their low scaling with number of particles, QMC methods present a compelling competitive alternative for the accurate study of large molecular systems and solid state calculations.In spite of such promise, the method has not permeated the quantum chemistry community broadly, mainly because of the fixed-node error, which can be large and whose control is difficult.In this Perspective, we present a systematic application of large scale multi-determinant expansions in QMC, and report on its impressive performance with first row dimers and the 55 molecules of the G1 test set.We demonstrate the potential of this strategy for systematically reducing the fixed-node error in the wave function and for achieving chemical accuracy in energy predictions.When compared to traditional quantum chemistry methods like MP2, CCSD(T), and various DFT approximations, the QMC results show a marked improvement over all of them.In fact, only the explicitly-correlated CCSD(T) method with a large basis set produces more accurate results.Further developments in trial wave functions and algorithmic improvements appear promising for rendering QMC as the benchmark standard in large electronic systems.
An alternative approach to conventional quantum chemistry techniques for molecular studies is described.This approach uses the quantum Monte Carlo method, which was developed and used primarily in nuclear and condensed-matter physics.In this approach the many-body Schrodinger equation is re-interpreted as a diffusion equation.Simulation of the appropriate random walk process allows one to calculate expectation values of molecular properties.In principle these expectation values can be calculated exactly, subject only to statistical errors (which may be made arbitrarily small).In preliminary work on small molecules (H 2 , LiH, Li 2 , and H 2 0), we use a simple but accurate approximation to ease the treatment of Fermi statistics.In that approximation, the calculated total energy remains an upper bound to the true energy.How good the bound is, as well as the magnitude of the statistical error, depends on the quality of an "importance function" which guides the diffusion into more probable regions of phase space.With relatively simple importance functions, we obtain from 75-100% of the correlation energy of the above-mentioned molecules.