Wave Function Optimization in Quantum Monte Carlo

Total Page:16

File Type:pdf, Size:1020Kb

Wave Function Optimization in Quantum Monte Carlo Wave Function Optimization in Quantum Monte Carlo Arne Lüchow Theoretical Chemistry Group, RWTH Aachen University Autumn School on Correlated Electrons, FZ Jülich Sept. 16, 2019 1 / 32 Outline Brief overview of quantum Monte Carlo methods (QMC) Optimization of wave functions Examples Recent developments 2 / 32 Two main variants: 1 solve stationary Schrödinger equation with stochastic process: Diffusion Monte Carlo (DMC) Simulate stochastic process with random numbers 2 evaluate E = hΨjHjΨi with Monte Carlo integration: Variational Monte Carlo (VMC) Sample the integrand with random numbers Quantum Monte Carlo – a brief overview Quantum Monte Carlo methods Numerical methods to solve the Schrödinger equation by means of Monte Carlo techniques (i.e. using random sampling) are called quantum Monte Carlo methods (QMC) electron structure QMC for electronic Schrödinger equation 3 / 32 2 evaluate E = hΨjHjΨi with Monte Carlo integration: Variational Monte Carlo (VMC) Sample the integrand with random numbers Quantum Monte Carlo – a brief overview Quantum Monte Carlo methods Numerical methods to solve the Schrödinger equation by means of Monte Carlo techniques (i.e. using random sampling) are called quantum Monte Carlo methods (QMC) electron structure QMC for electronic Schrödinger equation Two main variants: 1 solve stationary Schrödinger equation with stochastic process: Diffusion Monte Carlo (DMC) Simulate stochastic process with random numbers 3 / 32 Quantum Monte Carlo – a brief overview Quantum Monte Carlo methods Numerical methods to solve the Schrödinger equation by means of Monte Carlo techniques (i.e. using random sampling) are called quantum Monte Carlo methods (QMC) electron structure QMC for electronic Schrödinger equation Two main variants: 1 solve stationary Schrödinger equation with stochastic process: Diffusion Monte Carlo (DMC) Simulate stochastic process with random numbers 2 evaluate E = hΨjHjΨi with Monte Carlo integration: Variational Monte Carlo (VMC) Sample the integrand with random numbers 3 / 32 Mathematically, this is a diffusion equation In diffusion Monte Carlo (DMC), this diffusion process is simulated with a stochastic process1 Diffusion Monte Carlo (DMC) The starting point: Schrödinger equation in imaginary time t = it ¶ 1 Ψ = −HΨ = ∆Ψ − V Ψ ¶t 2 Long-time solution is the ground state wave function Ψ0 −Ht lim e Ψ ∝ Ψ0 t!¥ 1A. Lüchow, WIREs Comput Mol Sci 2011, 1, 388-402 4 / 32 Diffusion Monte Carlo (DMC) The starting point: Schrödinger equation in imaginary time t = it ¶ 1 Ψ = −HΨ = ∆Ψ − V Ψ ¶t 2 Long-time solution is the ground state wave function Ψ0 −Ht lim e Ψ ∝ Ψ0 t!¥ Mathematically, this is a diffusion equation In diffusion Monte Carlo (DMC), this diffusion process is simulated with a stochastic process1 1A. Lüchow, WIREs Comput Mol Sci 2011, 1, 388-402 4 / 32 2 stationary distribution of drift-diffusion process is jΨG j for ∆t ! 0. DMC: weights 0 − 1 (E (x)+E (x 0))−E ∆t W (∆t;x;x ) = e [ 2 L L ref ] fluctuate around 1 if local energies EL = Hˆ ΨG =ΨG are close to Eref (nonsingular!) Importance Sampling: Drift-Diffusion process discretized stochastic process: time step ∆t p Xk+1 = Xk + b(Xk )∆t + ∆tx (drift-diffusion process) drift term b(x) = ∇ΨG =ΨG : non-random moves to large values of jΨG j (guide function ΨG ). 5 / 32 DMC: weights 0 − 1 (E (x)+E (x 0))−E ∆t W (∆t;x;x ) = e [ 2 L L ref ] fluctuate around 1 if local energies EL = Hˆ ΨG =ΨG are close to Eref (nonsingular!) Importance Sampling: Drift-Diffusion process discretized stochastic process: time step ∆t p Xk+1 = Xk + b(Xk )∆t + ∆tx (drift-diffusion process) drift term b(x) = ∇ΨG =ΨG : non-random moves to large values of jΨG j (guide function ΨG ). 2 stationary distribution of drift-diffusion process is jΨG j for ∆t ! 0. 5 / 32 Importance Sampling: Drift-Diffusion process discretized stochastic process: time step ∆t p Xk+1 = Xk + b(Xk )∆t + ∆tx (drift-diffusion process) drift term b(x) = ∇ΨG =ΨG : non-random moves to large values of jΨG j (guide function ΨG ). 2 stationary distribution of drift-diffusion process is jΨG j for ∆t ! 0. DMC: weights 0 − 1 (E (x)+E (x 0))−E ∆t W (∆t;x;x ) = e [ 2 L L ref ] fluctuate around 1 if local energies EL = Hˆ ΨG =ΨG are close to Eref (nonsingular!) 5 / 32 Schrödinger equation is solved with additional boundary condition given by nodes of a guide function ΨG : ΨG (R) = 0 R: all electron coords: R = (x1;y1;z1;:::;zn)! Nodes (nodal set): 3n − 1 hypersurface, unrelated to orbital nodes Fermion sign problem Perfect for vibrational ground state (vibrational DMC). BUT: lowest eigenstate of H is not the electronic ground state (but a bosonic state) for n > 2 The stochastic process is forced to an antisymmetric state by fixed node approximation 6 / 32 Fermion sign problem Perfect for vibrational ground state (vibrational DMC). BUT: lowest eigenstate of H is not the electronic ground state (but a bosonic state) for n > 2 The stochastic process is forced to an antisymmetric state by fixed node approximation Schrödinger equation is solved with additional boundary condition given by nodes of a guide function ΨG : ΨG (R) = 0 R: all electron coords: R = (x1;y1;z1;:::;zn)! Nodes (nodal set): 3n − 1 hypersurface, unrelated to orbital nodes 6 / 32 Monte Carlo integration with sample fRk g ∼ p(R): 1 N EVMC = hELip = lim EL(Rk ) VMC N! ∑ ¥ N k=1 Metropolis-Hastings often: drift-diffusion step Variational Monte Carlo (VMC) Variational energy of a trial function ΨT : hΨT jHjΨT i Z EVMC = = EL(R) · p(R)dR hΨT jΨT i with local energy: HΨT (R) EL(R) = ΨT (R) and probability distribution: jΨ(R)j2 p(R) = R jΨ(R)j2dR 7 / 32 Variational Monte Carlo (VMC) Variational energy of a trial function ΨT : hΨT jHjΨT i Z EVMC = = EL(R) · p(R)dR hΨT jΨT i with local energy: HΨT (R) EL(R) = ΨT (R) and probability distribution: jΨ(R)j2 p(R) = R jΨ(R)j2dR Monte Carlo integration with sample fRk g ∼ p(R): 1 N EVMC = hELip = lim EL(Rk ) VMC N! ∑ ¥ N k=1 Metropolis-Hastings often: drift-diffusion step 7 / 32 Variational quantum Monte Carlo (VMC): VMC energy: average of local energies hΨT jHjΨT i EVMC = hELipVMC = hΨT jΨT i DMC more accurate, VMC more efficient VMC and DMC energies Diffusion Monte Carlo (DMC): (FN) weighted/branched drift-diffusion: pDMC (R) ∝ ΨG Ψ0 DMC energy: average of local energies hΨ(FN)jHjΨ i E = hE i = 0 G DMC L pDMC (FN) hΨ0 jΨG i (after time step extrapolation to ∆t ! 0) 8 / 32 VMC and DMC energies Diffusion Monte Carlo (DMC): (FN) weighted/branched drift-diffusion: pDMC (R) ∝ ΨG Ψ0 DMC energy: average of local energies hΨ(FN)jHjΨ i E = hE i = 0 G DMC L pDMC (FN) hΨ0 jΨG i (after time step extrapolation to ∆t ! 0) Variational quantum Monte Carlo (VMC): VMC energy: average of local energies hΨT jHjΨT i EVMC = hELipVMC = hΨT jΨT i DMC more accurate, VMC more efficient 8 / 32 Computational Considerations VMC and DMC require guide or trial function Ψ for local energy no integrals! Only 2nd derivatives! Guide and Trial Functions Compact and physically motivated approximations to exact wave function. Explicitly correlated! Scaling of local energy EL = HΨ=Ψ with number of electrons 2 Variance h(EL − hELi) i determines sample size Scales well to 10,000 of cpu cores 9 / 32 Small parameter changes p0 to p: weighted mean K 2 ∑k=1 wk EL(Rk ;p) jΨ(Rk ;p)j EVMC = Ewm(p) = K ; wk = 2 : ∑k=1 wk jΨ(Rk ;p0)j Stochastic optimization of wave function: minimization of energy Optimization of wave functions: energy minimization large sample of electron configurations fRk g; k = 1;:::;K with p(R) = jΨ(R)j2=R jΨ(R)j2dt parameter vector p wave function Ψ(R;p) ≡ Ψ(p) local energy EL(R;p) ≡ EL(p) 1 K EVMC = Em(p) = ∑ EL(Rk ;p): K k=1 10 / 32 Optimization of wave functions: energy minimization large sample of electron configurations fRk g; k = 1;:::;K with p(R) = jΨ(R)j2=R jΨ(R)j2dt parameter vector p wave function Ψ(R;p) ≡ Ψ(p) local energy EL(R;p) ≡ EL(p) 1 K EVMC = Em(p) = ∑ EL(Rk ;p): K k=1 Small parameter changes p0 to p: weighted mean K 2 ∑k=1 wk EL(Rk ;p) jΨ(Rk ;p)j EVMC = Ewm(p) = K ; wk = 2 : ∑k=1 wk jΨ(Rk ;p0)j Stochastic optimization of wave function: minimization of energy 10 / 32 Optimization of wave functions: variance minimization variance of the local energy K 1 2 Vm(p) = ∑ (EL(Rk ;p) − Em) K k=1 exact eigenfunctions of H have Vm = 0 Advantage: stable, each term positive Variance minimization: very stable and efficient less accurate energies (VMC and DMC) 11 / 32 BUT: (FN) projected Ψ0 depends on (guide) function parameters (FN) pDMC = Ψ0 ΨG only statistically known derivatives? AND: use optimized VMC trial function as guide functions in DMC! VMC optimization can optimize DMC energy parameters modifying nodes of wave function required! FN-DMC energy: optimal wf with given nodes flexible wf: optimal VMC energy close to DMC energy VMC vs DMC wave function optimization DMC more accurate than VMC: DMC optimization desired! 12 / 32 AND: use optimized VMC trial function as guide functions in DMC! VMC optimization can optimize DMC energy parameters modifying nodes of wave function required! FN-DMC energy: optimal wf with given nodes flexible wf: optimal VMC energy close to DMC energy VMC vs DMC wave function optimization DMC more accurate than VMC: DMC optimization desired! BUT: (FN) projected Ψ0 depends on (guide) function parameters (FN) pDMC = Ψ0 ΨG only statistically known derivatives? 12 / 32 VMC vs DMC wave function optimization DMC more accurate than VMC: DMC optimization desired! BUT: (FN) projected Ψ0 depends on (guide) function parameters (FN) pDMC = Ψ0 ΨG only statistically known derivatives? AND: use optimized VMC trial function as guide functions in DMC! VMC optimization
Recommended publications
  • Chapter 5 Quantum Monte Carlo
    Chapter 5 Quantum Monte Carlo This chapter is devoted to the study of quatum many body systems using Monte Carlo techniques. We analyze two of the methods that belong to the large family of the quantum Monte Carlo techniques, namely the Path-Integral Monte Carlo (PIMC) and the Diffusion Monte Carlo (DMC, also named Green’s function Monte Carlo). In the first section we start by introducing PIMC. 5.1 Path Integrals in Quantum Statistical Mechan- ics In this section we introduce the path-integral description of the properties of quantum many-body systems. We show that path integrals permit to calculate the static prop- erties of systems of Bosons at thermal equilibrium by means of Monte Carlo methods. We consider a many-particle system described by the non-relativistic Hamiltonian Hˆ = Tˆ + Vˆ ; (5.1) in coordinate representation the kinetic operator Tˆ and the potential operator Vˆ are defined as: !2 N Tˆ = ∆ , and (5.2) −2m i i=1 ! Vˆ = V (R) . (5.3) In these equations ! is the Plank’s constant divided by 2π, m the particles mass, N the number of particles and the vector R (r1,...,rN ) describes their positions. We consider here systems in d dimensions, with≡ fixed number of particles, temperature T , contained in a volume V . In most case, the potential V (R) is determined by inter-particle interactions, in which case it can be written as the sum of pair contributions V (R)= v (r r ), where i<j i − j v(r) is the inter-particle potential; it can also be due to an external field, call it vext(r), " R r in which case it is just the sum of single particle contributions V ( )= i vex ( i).
    [Show full text]
  • Few-Body Systems in Condensed Matter Physics
    Few-body systems in condensed matter physics. Roman Ya. Kezerashvili1,2 1Physics Department, New York City College of Technology, The City University of New York, Brooklyn, NY 11201, USA 2The Graduate School and University Center, The City University of New York, New York, NY 10016, USA This review focuses on the studies and computations of few-body systems of electrons and holes in condensed matter physics. We analyze and illustrate the application of a variety of methods for description of two- three- and four-body excitonic complexes such as an exciton, trion and biexciton in three-, two- and one-dimensional configuration spaces in various types of materials. We discuss and analyze the contributions made over the years to understanding how the reduction of dimensionality affects the binding energy of excitons, trions and biexcitons in bulk and low- dimensional semiconductors and address the challenges that still remain. I. INTRODUCTION Over the past 50 years, the physics of few-body systems has received substantial development. Starting from the mid-1950’s, the scientific community has made great progress and success towards the devel- opment of methods of theoretical physics for the solution of few-body problems in physics. In 1957, Skornyakov and Ter-Martirosyan [1] solved the quantum three-body problem and derived equations for the determination of the wave function of a system of three identical particles (fermions) in the limiting case of zero-range forces. The integral equation approach [1] was generalized by Faddeev [2] to include finite and long range interactions. It was shown that the eigenfunctions of the Hamiltonian of a three- particle system with pair interaction can be represented in a natural fashion as the sum of three terms, for each of which there exists a coupled set of equations.
    [Show full text]
  • Introduction to the Diffusion Monte Carlo Method
    Introduction to the Diffusion Monte Carlo Method Ioan Kosztin, Byron Faber and Klaus Schulten Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801 (August 25, 1995) A self–contained and tutorial presentation of the diffusion Monte Carlo method for determining the ground state energy and wave function of quantum systems is provided. First, the theoretical basis of the method is derived and then a numerical algorithm is formulated. The algorithm is applied to determine the ground state of the harmonic oscillator, the Morse oscillator, the hydrogen + atom, and the electronic ground state of the H2 ion and of the H2 molecule. A computer program on which the sample calculations are based is available upon request. E-print: physics/9702023 I. INTRODUCTION diffusion equation. The diffusion–reaction equation aris- ing can be solved by employing stochastic calculus as it 1,2 The Schr¨odinger equation provides the accepted de- was first suggested by Fermi around 1945 . Indeed, the scription for microscopic phenomena at non-relativistic imaginary time Schr¨odinger equation can be solved by energies. Many molecular and solid state systems simulating random walks of particles which are subject to are governed by this equation. Unfortunately, the birth/death processes imposed by the source/sink term. Schr¨odinger equation can be solved analytically only in a The probability distribution of the random walks is iden- few highly idealized cases; for most realistic systems one tical to the wave function. This is possible only for wave needs to resort to numerical descriptions. In this paper functions which are positive everywhere, a feature, which we want to introduce the reader to a relatively recent nu- limits the range of applicability of the DMC method.
    [Show full text]
  • Arxiv:Physics/0609191V1 [Physics.Comp-Ph] 22 Sep 2006 Ot Al Ihteueo Yhn Ial,W Round VIII
    MontePython: Implementing Quantum Monte Carlo using Python J.K. Nilsen1,2 1 USIT, Postboks 1059 Blindern, N-0316 Oslo, Norway 2 Department of Physics, University of Oslo, N-0316 Oslo, Norway We present a cross-language C++/Python program for simulations of quantum mechanical sys- tems with the use of Quantum Monte Carlo (QMC) methods. We describe a system for which to apply QMC, the algorithms of variational Monte Carlo and diffusion Monte Carlo and we describe how to implement theses methods in pure C++ and C++/Python. Furthermore we check the effi- ciency of the implementations in serial and parallel cases to show that the overhead using Python can be negligible. I. INTRODUCTION II. THE SYSTEM In scientific programming there has always been a strug- Quantum Monte Carlo (QMC) has a wide range of appli- gle between computational efficiency and programming cations, for example studies of Bose-Einstein condensates efficiency. On one hand, we want a program to go as of dilute atomic gases (bosonic systems) [4] and studies fast as possible, resorting to low-level programming lan- of so-called quantum dots (fermionic systems) [5], elec- guages like FORTRAN77 and C which can be difficult to trons confined between layers in semi-conductors. In this read and even harder to debug. On the other hand, we paper we will focus on a model which is meant to repro- want the programming process to be as efficient as pos- duce the results from an experiment by Anderson & al. sible, turning to high-level software like Matlab, Octave, [6]. Anderson & al.
    [Show full text]
  • Quantum Monte Carlo for Vibrating Molecules
    LBNL-39574 UC-401 ERNEST DRLANDD LAWRENCE BERKELEY NATIDNAL LABDRATDRY 'ERKELEY LAB I Quantum Monte Carlo for Vibrating Molecules W.R. Brown DEC 1 6 Chemical Sciences Division August 1996 Ph.D. Thesis OF THIS DocuMsvr 8 DISCLAIMER This document was prepared as an account of work sponsored by the United States Government. While this document is believed to contain correct information, neither the United States Government nor any agency thereof, nor The Regents of the University of California, nor any of their employees, makes any warranty, express or implied, or assumes any legal responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by its trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or The Regents of the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof, or The Regents of the University of California. Ernest Orlando Lawrence Berkeley National Laboratory is an equal opportunity employer. LBjL-39574 ^UC-401 Quantum Monte Carlo for Vibrating Molecules Willard Roger Brown Ph.D. Thesis Chemistry Department University of California, Berkeley and Chemical Sciences Division Ernest Orlando Lawrence Berkeley National Laboratory University of California Berkeley, CA 94720 August 1996 This work was supported by the Director, Office of Energy Research, Office of Basic Energy Sciences, Chemical Sciences Division, of the U.S.
    [Show full text]
  • Hellmann and Feynman Theorem Versus Diffusion Monte Carlo Experiment
    Solid State Communications 159 (2013) 106–109 Contents lists available at SciVerse ScienceDirect Solid State Communications journal homepage: www.elsevier.com/locate/ssc Hellmann and Feynman theorem versus diffusion Monte Carlo experiment Riccardo Fantoni n National Institute for Theoretical Physics (NITheP) and Institute of Theoretical Physics, Stellenbosch 7600, South Africa article info abstract Article history: In a computer experiment the choice of suitable estimators to measure a physical quantity plays an Received 13 February 2012 important role. We propose a new direct route to determine estimators for observables which do not Received in revised form commute with the Hamiltonian. Our new route makes use of the Hellmann and Feynman theorem and 27 January 2013 in a diffusion Monte Carlo simulation it introduces a new bias to the measure due to the choice of the Accepted 29 January 2013 auxiliary function. This bias is independent from the usual one due to the choice of the trial wave by H. Akai Available online 8 February 2013 function. We used our route to measure the radial distribution function of a spin one half Fermion fluid. & 2013 Elsevier Ltd. All rights reserved. Keywords: Hellmann and Feynman theorem Diffusion Monte Carlo Radial distribution function Jellium In a computer experiment of a many particles system, a fluid, then use our route to measure the radial distribution function of a the determination of suitable estimators to measure, through a spin one half Fermion fluid. statistical average, a given physical quantity, an observable, plays In ground state Monte Carlo simulations [1,2], unlike classical an important role. Whereas the average from different estimators Monte Carlo simulations [3–5] and path integral Monte Carlo must give the same result, the variance, the square of the simulations [6], one has to resort to the use of a trial wave statistical error can be different for different estimators.
    [Show full text]
  • The Uniform Electron Gas! We Have No Simpler Paradigm for the Study of Large Numbers Of
    The uniform electron gas Pierre-Fran¸cois Loos1, a) and Peter M. W. Gill1, b) Research School of Chemistry, Australian National University, Canberra ACT 2601, Australia The uniform electron gas or UEG (also known as jellium) is one of the most funda- mental models in condensed-matter physics and the cornerstone of the most popular approximation — the local-density approximation — within density-functional the- ory. In this article, we provide a detailed review on the energetics of the UEG at high, intermediate and low densities, and in one, two and three dimensions. We also re- port the best quantum Monte Carlo and symmetry-broken Hartree-Fock calculations available in the literature for the UEG and discuss the phase diagrams of jellium. Keywords: uniform electron gas; Wigner crystal; quantum Monte Carlo; phase dia- gram; density-functional theory arXiv:1601.03544v2 [physics.chem-ph] 19 Feb 2016 a)Electronic mail: [email protected]; Corresponding author b)Electronic mail: [email protected] 1 I. INTRODUCTION The final decades of the twentieth century witnessed a major revolution in solid-state and molecular physics, as the introduction of sophisticated exchange-correlation models1 propelled density-functional theory (DFT) from qualitative to quantitative usefulness. The apotheosis of this development was probably the award of the 1998 Nobel Prize for Chemistry to Walter Kohn2 and John Pople3 but its origins can be traced to the prescient efforts by Thomas, Fermi and Dirac, more than 70 years earlier, to understand the behavior of ensembles
    [Show full text]
  • Introduction to the Diffusion Monte Carlo Method
    Introduction to the Diffusion Monte Carlo Method Ioan Kosztin, Byron Faber and Klaus Schulten Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801 (August 25, 1995) A self–contained and tutorial presentation of the diffusion Monte Carlo method for determining the ground state energy and wave function of quantum systems is provided. First, the theoretical basis of the method is derived and then a numerical algorithm is formulated. The algorithm is applied to determine the ground state + of the harmonic oscillator, the Morse oscillator, the hydrogen atom, and the electronic ground state of the H2 ion and of the H2 molecule. A computer program on which the sample calculations are based is available upon request. I. INTRODUCTION by simulating random walks of particles which are subject to birth/death processes imposed by the source/sink term. The The Schrodinger¨ equation provides the accepted description probability distribution of the random walks is identical to the for microscopic phenomena at non-relativistic energies. Many wave function. This is possible only for wave functions which molecular and solid state systems are governed by this equa- are positive everywhere, a feature, which limits the range of applicability of the DMC method. Such a formulation of the tion. Unfortunately, the Schrodinger¨ equation can be solved 3 analytically only in a few highly idealized cases; for most re- DMC method was given for the first time by Anderson who used this method to calculate the ground state energy of small alistic systems one needs to resort to numerical descriptions.
    [Show full text]
  • Quantum Monte Carlo Method for Boson Ground States: Application to Trapped Bosons with Attractive and Repulsive Interactions
    W&M ScholarWorks Dissertations, Theses, and Masters Projects Theses, Dissertations, & Master Projects 2005 Quantum Monte Carlo method for boson ground states: Application to trapped bosons with attractive and repulsive interactions Wirawan Purwanto College of William & Mary - Arts & Sciences Follow this and additional works at: https://scholarworks.wm.edu/etd Part of the Atomic, Molecular and Optical Physics Commons, and the Condensed Matter Physics Commons Recommended Citation Purwanto, Wirawan, "Quantum Monte Carlo method for boson ground states: Application to trapped bosons with attractive and repulsive interactions" (2005). Dissertations, Theses, and Masters Projects. Paper 1539623468. https://dx.doi.org/doi:10.21220/s2-mynw-ys14 This Dissertation is brought to you for free and open access by the Theses, Dissertations, & Master Projects at W&M ScholarWorks. It has been accepted for inclusion in Dissertations, Theses, and Masters Projects by an authorized administrator of W&M ScholarWorks. For more information, please contact [email protected]. NOTE TO USERS This reproduction is the best copy available. ® UMI Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. Reproduced with permission ofof the the copyright copyright owner. owner. FurtherFurther reproduction reproduction prohibited prohibited without without permission. permission. QUANTUM MONTE CARLO METHOD FOR BOSON GROUND STATES Application to Trapped Bosons with Attractive and Repulsive Interactions A Dissertation Presented to The Faculty of the Department of Physics The College of William and Mary in Virginia In Partial Fulfillment Of the Requirements for the Degree of Doctor of Philosophy by Wirawan Purwanto 2005 Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. UMI Number: 3172894 INFORMATION TO USERS The quality of this reproduction is dependent upon the quality of the copy submitted.
    [Show full text]
  • Quantum Monte Carlo: Accuracy, Generality, and Scalability in the Electron Correlation Problem Scuola Normale, Pisa, April 2014
    Quantum Monte Carlo: accuracy, generality, and scalability in the electron correlation problem Scuola Normale, Pisa, April 2014 Mike Towler TCM Group, Cavendish Laboratory, University of Cambridge and UCL QMC web page: vallico.net/casinoqmc/ Email: [email protected] { Typeset by FoilTEX { 1 { Typeset by FoilTEX { 2 What is quantum Monte Carlo? QMC is, essentially, `accurate calculation of the properties of assemblies of interacting quantum particles using random sampling of the full many-body Schr¨odingerwave function..' For my purposes, it is applied to `continuum fermion' cases i.e. realistic atoms, molecules, and crystalline solids just as in regular quantum chemistry. Why use wave functions? We work directly with the many-body wave function Ψ(x1; x2;:::; xN ) rather than the density ρ(x) since we can write down exactly the equation that Ψ obeys (a very helpful thing to do when requiring accuracy. ): H^ Ψ = EΨ • Variational Monte Carlo (VMC) - a cheaper less accurate technique which directly computes expectation values such as the one below using standard Monte Carlo numerical integration. May also involve optimization of Ψ by direct variation of any parameters upon which it depends. Limited by flexibility of parameterization.. R ΨH^ Ψ dx E = VMC R Ψ2 dx • Diffusion Monte Carlo (DMC) - a more sophisticated, accurate (and expensive) method which - in principle - projects out the exact ground state wave function from a given `starting guess' and then computes expectation values as above. { Typeset by FoilTEX { 3 Good things about DMC • `Chemical accuracy' (1 kcal/mol or 0.04 eV) and beyond readily achievable. Can feed results of standard DFT/HF/etc.
    [Show full text]
  • An Introduction to Quantum Monte Carlo
    An introduction to Quantum Monte Carlo Jose Guilherme Vilhena Laboratoire de Physique de la Matière Condensée et Nanostructures, Université Claude Bernard Lyon 1 Outline 1) Quantum many body problem; 2) Metropolis Algorithm; 3) Statistical foundations of Monte Carlo methods; 4) Quantum Monte Carlo methods ; 4.1) Overview; 4.2) Variational Quantum Monte Carlo; 4.3) Diffusion Quantum Monte Carlo; 4.4) Typical QMC calculation; Outline 1) Quantum many body problem; 2) Metropolis Algorithm; 3) Statistical foundations of Monte Carlo methods; 4) Quantum Monte Carlo methods ; 4.1) Overview; 4.2) Variational Quantum Monte Carlo; 4.3) Diffusion Quantum Monte Carlo; 4.4) Typical QMC calculation; 1 - Quantum many body problem Consider of N electrons. Since m <<M , to a good e N approximation, the electronic dynamics is governed by : Z Z 1 2 1 H =− ∑ ∇ i −∑∑ ∑∑ 2 i i ∣r i−d∣ 2 i j ≠i ∣ri−r j∣ which is the Born Oppenheimer Hamiltonian. We want to find the eigenvalues and eigenvectos of this hamiltonian, i.e. : H r1, r 2,. ..,r N =E r 1, r2,. ..,r N QMC allow us to solve numerically this, thus providing E0 and Ψ . 0 1 - Quantum many body problem Consider of N electrons. Since m <<M , to a good e N approximation, the elotonic dynamics is governed by : Z Z 1 2 1 H =− ∑ ∇ i −∑∑ ∑∑ 2 i i ∣r i−d∣ 2 i j ≠i ∣ri−r j∣ which is the Born Oppenheimer Hamiltonian. We want to find the eigenvalues and eigenvectos of this hamiltonian, i.e. : H r1, r 2,.
    [Show full text]
  • Quantum Monte Carlo (QMC)
    Quantum Monte Carlo (QMC) December 14, 2018 Quantum Monte Carlo (QMC) December 14, 2018 1 / 241 Motivation We simulate many-body quantum systems by solving the Schr¨odinger equation stochastically, because there are too many degrees of freedom in anything but trivial systems Before you begin, you need: • Hamiltonian (fix particle masses and their interactions) • particle statistics: bosons, fermions, or sometimes boltzmannons (forget indistinquishability) { A non-relativistic Hamiltonian has no idea if particles are bosons or fermions • Do we want properties for the ground state, for excited states, or for a finite temperature system ? • Do you really need a first-principles result? DFT is faster than QMC. Quantum Monte Carlo (QMC) December 14, 2018 2 / 241 QMC methods • Variational Monte Carlo (VMC) { Sample particle positions from a guessed wave function and simulate the consequences in a system with a certain Hamiltonian { You may improve your guess (for lower energy or variance), but it's not automatic • Diffusion Monte Carlo (DMC) { Usually for sampling particle positions in the ground state { Numerically improves upon a guessed wave function • Path Integral Monte Carlo (PIMC) { Simulate finite temperature quantum systems { Manifest Feynman's path integral representation of quantum mechanics { Also a ground-state variant exists Quantum Monte Carlo (QMC) December 14, 2018 3 / 241 History 1940's E. Fermi: Schr¨odingerequation can written as a diffusion equation; Physical principles of Diffusion Monte Carlo (DMC) discovered. 1942 A. A. Frost: evaluation of the energy of few simple molecules by picking few representative points in the configuration space; Variational Monte Carlo (VMC) is born. 1947 N.
    [Show full text]