Benchmark studies using quantum monte carlo: Pressure estimators, energy, and entanglement

Jeremy McMinis

Published May 24, 2013• 1 citation

Abstract

Trial wave function based quantum Monte Carlo is a promising family of methods for the solution of quantum mechanical Hamiltonians. These fundamentally non-perturbative methods can be used to treat bosons or fermions in weakly or strong interacting systems and in any phase. Properties computed for particles with bosonic statistics can be converged to their exact value at relatively moderate computational cost using diffusion quantum Monte Carlo. For fermions the problem is slightly more complicated. Exact methods, such as released-node diffusion quantum Monte Carlo, exist for computing unbiased fermion properties, however their computational expense increases rapidly with the Bose-Fermi energy gap and system size. It is possible to use a fixed-node version of diffusion quantum Monte Carlo, but it introduces a bias due to the nodes of the trial wave function. In this thesis we work towards reducing trial wave function bias in fixed-node calculations and then perform several benchmark studies. We begin with a pedagogical overview of the three methods used throughout the thesis, variational, diffusion, and reptation quantum Monte Carlo. Then we discuss the problem, trial wave function bias in the energy and other observables. Next we outline the Hamiltonians we are typically interested in for electronic systems and trial wave function forms used to solve them. The algorithm used to optimize the trial wave functions is presented along with special considerations for some particular cases. Next we present three studies of the pressure of the electron gas using improved estimators in variational, diffusion, and reptation quantum Monte Carlo. Benchmark results for a set of molecules are then presented for the massive multi-determinant expansion and optimization algorithm tailored to suit it. We conclude with two studies of the interacting Fermi liquid, one on the momentum properties of the electron gas in three dimensions and another on the entanglement properties of interacting Fermi liquids in two dimensions.

Related papers

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.3417Physical Review BMay 19, 200930 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.06006The Journal of Chemical PhysicsDec 28, 20249 citations

Computational Studies Of Strongly Correlated Quantum Matter

Hao Shi

The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.

Jan 1, 20170 citations

Monte Carlo studies in weakly interacting quantum systems

Ryan Springall

This thesis is composed of two parts. In part one we report results of two quantum Monte Carlo methods - variational Monte Carlo and diffusion Monte Carlo - on the potential energy curve of the helium dimer, the prototypical van der Waals system. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all. In part two we switch to investigations of many-body systems at finite temperature. We use the path integral representation of statistical mechanics to investigate the symmetry properties of the canonical ensemble partition function and find a representation in terms of irreducible representations of the symmetric group. We use this as a foundation to propose a novel technique for the stochastic sampling of the bosonic partition function for quantum gases. It is shown that in principle we are able to use two operators which enable us to construct a Markov chain through a graph of the irreducible representation of the symmetric group. As an illustration of this method, a test calculation of four particles in a harmonic trap is performed.

Jun 27, 20120 citations

A novel method for optimizing quantum Monte Carlo wave functions

Hongxin Huang, Zexing Cao

We propose an algorithm for optimizing quantum Monte Carlo wave functions. An improved steepest-descent technique is used, with step size automatically adjustable to obtain a procedure that converges superlinearly. We also propose a novel trial function, which has both correct electron–electron cusp conditions and correct electron–nucleus cusp conditions. To test the optimization procedure and the optimized trial function, the ground states for CH4 and H2O molecules were investigated using variational Monte Carlo (VMC) and fixed-node quantum Monte Carlo (FNQMC) calculations. For CH4 and H2O, the VMC recovered 73.3% and 57.9% of the correlation energy, respectively, and the FNQMC recovered 99.3% and 92.8%, respectively. The optimization procedure is three to five times faster than conventional usual steepest-descent procedures. The trial functions optimized are more accurate than prior trial functions of similar complexity.

The Journal of Chemical PhysicsJan 1, 199626 citations

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, 20200 citations