Variational Monte Carlo Method for Fermions

Total Page:16

File Type:pdf, Size:1020Kb

Variational Monte Carlo Method for Fermions Variational Monte Carlo Method for Fermions Vijay B. Shenoy [email protected] Indian Institute of Science, Bangalore December 26, 2009 1/91 Acknowledgements Sandeep Pathak Department of Science and Technology 2/91 Plan of the Lectures on Variational Monte Carlo Method Motivation, Overview and Key Ideas Details of VMC Application to Superconductivity 3/91 The Strongly Correlated Panorama... Fig. 1. Phase diagrams ofrepresentativemateri- als of the strongly cor- related electron family (notations are standard and details can be found in the original refer- ences). (A) Temperature versus hole density phase diagram of bilayer man- ganites (74), including several types of antiferro- magnetic (AF) phases, a ferromagnetic (FM) phase, and even a glob- ally disordered region at x 0 0.75. (B) Generic phase diagram for HTSC. SG stands for spin glass. (C) Phase diagram of single layered ruthenates (75, 76), evolving from a superconducting (SC) state at x 0 2 to an AF insulator at x 0 0(x controls the bandwidth rather than the carrier density). Ruthenates are believed to be clean metals at least at large x, thus providing a family of oxides where competition and com- plexity can be studied with less quenched dis- order than in Mn ox- ides. (D) Phase diagram of Co oxides (77), with SC, charge-ordered (CO), and magnetic regimes. (E) Phase diagram of the organic k-(BEDT- TTF)2Cu[N(CN)2]Cl salt (57). The hatched re- gion denotes the co- existence of metal and insulator phases. (F) Schematic phase dia- gram of the Ce-based heavy fermion materials (51). (Dagotto, Science, (2005)) 4/91 ...and now, Cold Atom Quantum Simulators (Bloch et al., RMP, (2009)) Hamiltonian...for here, or to go? 5/91 The Strategy Identify “key degrees of freedom” of the system of interest Identify the key interactions Write a Hamiltonian (eg. tJ-model)...reduce the problem to a model “Solve” the Hamiltonian .... Nature(Science)/Nature Physics/PRL/PRB... 6/91 “Solvers” Mean Field Theories (MFT) Slave Particle Theories (and associated MFTs) Dynamical Mean Field Theories (DMFT) Exact Diagonalization Quantum Monte Carlo Methods (eg. Determinental) Variational Monte Carlo (VMC) ... AdS/CFT 7/91 Great Expectations Prerequisites Graduate quantum mechanics including “second quantization” Ideas of “Mean Field Theory” Some familiarity with “classical” superconductivity At the end of the lectures, you can/will have Understood concepts of VMC Know-how to write your own code A free code (pedestrianVMC) Seen examples of VMC in strongly correlated systems 8/91 References NATO Series Lecture Notes by Nightingale, Umrigar, see http://www.phys.uri.edu/˜nigh/QMC-NATO/webpage/ abstracts/lectures.html Morningstar, hep-lat/070202 Ceperley, Chester and Kalos, PRB 17, 3081 (1977) 9/91 VMC in a Nutshell Strategy to study ground (and some special excited) state(s) “Guess” the ground state Use Schr¨odinger’s variational principle to optimize the ground state Calculate observables of interest with the optimal state ... 10/91 Hall of Fame Bardeen, Cooper, Schrieffer Laughlin Quote/Question from A2Z, Plain Vanilla JPCM (2004) The philosophy of this method (VMC )is analogous to that used by BCS for superconductivity, and by Laughlin for the fractional quantum Hall effect: simply guess a wavefunction. Is there any better way to solve a non-perturbative many-body problem? 11/91 Variational Principle Eigenstates of the Hamiltonian satisfy H Ψ = E Ψ | i | i where E is the energy eigenvalue Alternate way of stating this ◮ Define an “energy functional” for any given state Φ | i E[Φ] = Φ H Φ h | | i ◮ Among all states Φ that satisfy Φ Φ = 1, those that extremize E[Φ] | i h | i are eigenstates of H Exercise: Prove this! ◮ The state(s) that minimises(minimise) E[Φ] is(are) the ground state(s) 12/91 Variational Principle in Action...A Simple Example + Ground state H2 molecule with a given nuclear separation “Guess” the wave function (linear combination of atomic orbitals) ψ = a 1 + b 2 | i | i | i where 1, 2 are atomic orbitals centered around protons 1 and 2 | i Find energy function (taking everything real) 2 2 E(a, b)= H11a + 2H12ab + H22b (symbols have obvious meaning) 1 Minimize (with appropriate constraints) to find ψ g = ( 1 + 2 ) | i √2 | i | i (for H12 < 0) Question: Is this the actual ground state? Example shows key ideas (and limitations) of the variational approach! 13/91 Variational Approach in Many Body Physics Guess the ground state wavefunction Ψ | i Ψ depends on some variational parameters α1, α2... which “contain the| i physics” of the wavefunction ◮ Example, † † BCS = uk + vk c c 0 | i k↑ −k↓ | i " k # Y uk , vk are the “αi s” In general, αi s are chosen ◮ To realize some sort of “order” or “broken symmetry”, (for example, magnetic order)...there could be more than one order present (for example, superconductivity and magnetism) ◮ To enforce some sort of “correlation” that we “simply” know exists and can include all of the “effects” above ...your guess is not as good as mine! 14/91 Variational Principle vs. MFT Mean Field theory ◮ Replace H by HMF which is usually quadratic ◮ Use the Feynman identity F FMF + H HMF MF (F is free energy) ≤ h − i ◮ Kills “fluctuations” Variational Principle ◮ Works with the full Hamiltonian ◮ Non-perturbative ◮ Well designed wavefunctions can include fluctuations above MFT Usually MFT provides us with intuitive starting point...we can then “well design” the variational wavefunction to include the fluctuations that we miss 15/91 Overall Strategy of Variational Approach Summary of Variational Approach 1 Design a ground state Ψ(αi ) | i 2 Energy function Ψ(αi ) H Ψ(αi ) E(αi )= h | | i Ψ(αi ) Ψ(αi ) h | i m 3 Find α that minimize (extremize) E(αi ) and obtain ΨG i | i 4 Study ΨG | i This explains the “Variational” part of VMC ... But...why Monte Carlo? 16/91 Variational Approach: The Difficulty In many cases the variational programme can be carried out analytically where we have a well controlled analytical expression for the state, e. g., BCS state In most problems of interest we do not have a simple enough analytical expression for the wavefunction (e. g., Gutzwiller D wavefunction Ψgutz = g Ψfree ) | i | i In many of these cases we have expressions for the wavefunction which “look simple in real space”, say, Ψα(r1 , r2 ,..., rN, , r1 , r2 ,..., rN, ), where riσ are the coordinates of our↑ electrons↑ on↓ a lattice↓ ↓ (we will↓ work exclusively with lattices, most ideas can be taken to continuum in a straightforward way); Ψα has necessary permutation symmetries, and α subscripts reminds us of the underlying parameters Call C (r1 , r2 ,..., rN, , r1 , r2 ,..., rN, ) a configuration of the ≡ ↑ ↑ ↓ ↓ ↓ ↓ electrons 17/91 Variational Approach: The Difficulty Now, Ψ (C)HΨ (C) E(α )= C α∗ α i Ψ (C)Ψ (C) P C α∗ α i.e., we need to do a sum overP all configurations C of the our 2N electron system Note that in most practical cases we need sum the denominator as well, since we use unnormalized wavefunctions! Sums...how difficult can they be? Consider 50 electrons, and 50 electrons on a lattice with 100 ↑ ↓ sites....How many configurations are there? Num. Configurations 1060!!!! ∼ And this is the difficulty! ... Solution...use Monte Carlo method to perform the sum! 18/91 “Monte Carlo” of VMC How do we use Monte Carlo to evaluate E(αi )? Use the following trick Poor Notation Warning: HΨ(C) C H Ψ ≡ h | | i HΨα(C) Ψ (C)HΨ (C) C Ψα∗ (C)Ψα(C) Ψ (C) E(α ) = C α∗ α = α i Ψ (C)Ψ (C) Ψ (C)Ψ (C) P C α∗ α P C α∗ α HΨ (C) = PP(C) α P Ψ (C) C α X where Ψ (C) 2 P(C)= | α | C Ψα∗ (C)Ψα(C) is the probability of the configurationP C; obviously C P(C) = 1 Note that we do not know the denominator of P(C) P For any operator O OΨ (C) O = P(C) α h i Ψ (C) C α X 19/91 Idea of Monte Carlo Sums Key idea is not to sum over all configurations C, but to visit the “most important” configurations and add up their contribution The most important configurations are those that have “high” probability..., but we do not know the denominator of P(C) Need to develop a scheme where we start with a configuration C1 and “step” to another configuration C2 comparing P(C1) and P(C2)...and continue on Thus we do a “walk” in the space of configurations depending on P(C)...spending “most of our time” in high probability configurations...this is importance sampling Goal: construct a set of rules (algorithm) to “step” in the configuration space that performs the necessary importance sampling Such a stepping produces a sequence configurations a function of steps or “time”...but how do we create this “dynamics”? 20/91 Markov Chains Let us set notation: Label all the configurations (yes, all 1060) using natural numbers C , C ...Cm... etc; we will use P(1) P(C ), and in 1 2 ≡ 1 general P(m) P(Cm) and ≡ C ≡ m We construct transition (Markov) matrix W (n, m) W (n m) P P ≡ ← which defines the probability of moving from the configuration m to configuration n in a step, we consider only W s that are independent of time The probability of taking a step from m n is independent of how → the system got to the state m...such stochastic processes are called as Markov chain Clearly W (n, m) = 1, and W (n, m) 0 Exercise: If λ is a right n ≥ eigenvalue of the matrix W, show that |λ| ≤ 1.
Recommended publications
  • An Overview of Quantum Monte Carlo Methods David M. Ceperley
    An Overview of Quantum Monte Carlo Methods David M. Ceperley Department of Physics and National Center for Supercomputing Applications University of Illinois Urbana-Champaign Urbana, Illinois 61801 In this brief article, I introduce the various types of quantum Monte Carlo (QMC) methods, in particular, those that are applicable to systems in extreme regimes of temperature and pressure. We give references to longer articles where detailed discussion of applications and algorithms appear. MOTIVATION One does QMC for the same reason as one does classical simulations; there is no other method able to treat exactly the quantum many-body problem aside from the direct simulation method where electrons and ions are directly represented as particles, instead of as a “fluid” as is done in mean-field based methods. However, quantum systems are more difficult than classical systems because one does not know the distribution to be sampled, it must be solved for. In fact, it is not known today which quantum problems can be “solved” with simulation on a classical computer in a reasonable amount of computer time. One knows that certain systems, such as most quantum many-body systems in 1D and most bosonic systems are amenable to solution with Monte Carlo methods, however, the “sign problem” prevents making the same statement for systems with electrons in 3D. Some limitations or approximations are needed in practice. On the other hand, in contrast to simulation of classical systems, one does know the Hamiltonian exactly: namely charged particles interacting with a Coulomb potential. Given sufficient computer resources, the results can be of quite high quality and for systems where there is little reliable experimental data.
    [Show full text]
  • Variational Monte Carlo Technique Ground State Energies of Quantum Mechanical Systems
    GENERAL ⎜ ARTICLE Variational Monte Carlo Technique Ground State Energies of Quantum Mechanical Systems Sukanta Deb In this article, I discuss the use of the Variational Monte Carlo technique in the determination of ground state energies of quantum harmonic os- cillators in one and three dimensions. In order to provide a better understanding of this technique, a pre-requisite concept of Monte Carlo integra- Sukanta Deb is an tion is also discussed along with numerical ex- Assistant Professor in the amples. The technique can be easily extended Department of Physics at Acharya Narendra Dev to study the ground state energies of hydrogen College (University of atom, particle in a box, helium atom and other Delhi). His research complex microscopic systems involving N parti- interests include astro- cles in d-dimensions. nomical photometry and spectroscopy, light curve 1. Introduction modeling in astronomy and astrophysics, Monte Carlo There exist many problems in science and engineering methods and simulations whose exact solution either does not exist or is difficult and automated data to find. For the solutions of those problems, one has to analysis. resort to approximate methods. In the context of quan- tum mechanics, the exact solutions of the Schr¨odinger’s equation exist only for a few idealized systems. There are varieties of systems for which the Schr¨odinger’s equa- tion either cannot be solved exactly or it involves very lengthy and cumbersome calculations [1]. For example, solving the Schr¨odinger equation to study the proper- ties of a system with hundreds or thousands of particles often turns out to be impractical [2].
    [Show full text]
  • Unifying Cps Nqs
    Citation for published version: Clark, SR 2018, 'Unifying Neural-network Quantum States and Correlator Product States via Tensor Networks', Journal of Physics A: Mathematical and General, vol. 51, no. 13, 135301. https://doi.org/10.1088/1751- 8121/aaaaf2 DOI: 10.1088/1751-8121/aaaaf2 Publication date: 2018 Document Version Peer reviewed version Link to publication This is an author-created, in-copyedited version of an article published in Journal of Physics A: Mathematical and Theoretical. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://iopscience.iop.org/article/10.1088/1751-8121/aaaaf2/meta University of Bath Alternative formats If you require this document in an alternative format, please contact: [email protected] General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Download date: 04. Oct. 2021 Unifying Neural-network Quantum States and Correlator Product States via Tensor Networks Stephen R Clarkyz yDepartment of Physics, University of Bath, Claverton Down, Bath BA2 7AY, U.K. zMax Planck Institute for the Structure and Dynamics of Matter, University of Hamburg CFEL, Hamburg, Germany E-mail: [email protected] Abstract.
    [Show full text]
  • Arxiv:1902.04057V3 [Cond-Mat.Dis-Nn] 19 Jan 2020 Statistical Mechanics Applications [16], As Well As Den- Spin Systems [13]
    Deep autoregressive models for the efficient variational simulation of many-body quantum systems 1, 1, 1, 2, 1, Or Sharir, ∗ Yoav Levine, † Noam Wies, ‡ Giuseppe Carleo, § and Amnon Shashua ¶ 1The Hebrew University of Jerusalem, Jerusalem, 9190401, Israel 2Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, NY 10010, USA Artificial Neural Networks were recently shown to be an efficient representation of highly-entangled many-body quantum states. In practical applications, neural-network states inherit numerical schemes used in Variational Monte Carlo, most notably the use of Markov-Chain Monte-Carlo (MCMC) sampling to estimate quantum expectations. The local stochastic sampling in MCMC caps the potential advantages of neural networks in two ways: (i) Its intrinsic computational cost sets stringent practical limits on the width and depth of the networks, and therefore limits their expressive capacity; (ii) Its difficulty in generating precise and uncorrelated samples can result in estimations of observables that are very far from their true value. Inspired by the state-of-the-art generative models used in machine learning, we propose a specialized Neural Network architecture that supports efficient and exact sampling, completely circumventing the need for Markov Chain sampling. We demonstrate our approach for two-dimensional interacting spin models, showcas- ing the ability to obtain accurate results on larger system sizes than those currently accessible to neural-network quantum states. Introduction.– The theoretical understanding and been limited to relatively shallow architectures, far from modeling of interacting many-body quantum matter rep- the very deep networks used in modern machine learn- resents an outstanding challenge since the early days of ing applications.
    [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]
  • University of Bath
    Citation for published version: Clark, SR 2018, 'Unifying Neural-network Quantum States and Correlator Product States via Tensor Networks', Journal of Physics A: Mathematical and General, vol. 51, no. 13, 135301. https://doi.org/10.1088/1751- 8121/aaaaf2 DOI: 10.1088/1751-8121/aaaaf2 Publication date: 2018 Document Version Peer reviewed version Link to publication This is an author-created, in-copyedited version of an article published in Journal of Physics A: Mathematical and Theoretical. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://iopscience.iop.org/article/10.1088/1751-8121/aaaaf2/meta University of Bath General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Download date: 23. May. 2019 Unifying Neural-network Quantum States and Correlator Product States via Tensor Networks Stephen R Clarkyz yDepartment of Physics, University of Bath, Claverton Down, Bath BA2 7AY, U.K. zMax Planck Institute for the Structure and Dynamics of Matter, University of Hamburg CFEL, Hamburg, Germany E-mail: [email protected] Abstract. Correlator product states (CPS) are a powerful and very broad class of states for quantum lattice systems whose (unnormalised) amplitudes in a fixed basis can be sampled exactly and efficiently.
    [Show full text]
  • Quantum Monte Carlo Simulations: from Few to Many Species
    Quantum Monte Carlo Simulations: From Few to Many Species by William Gary Dawkins A Thesis Presented to The University of Guelph In partial fulfilment of requirements for the degree of Master of Science in Physics Guelph, Ontario, Canada c William Gary Dawkins, May, 2019 ABSTRACT QUANTUM MONTE CARLO SIMULATIONS: FROM FEW TO MANY SPECIES William Gary Dawkins Advisor: University of Guelph, 2019 Professor Alexandros Gezerlis We study systems of particles with varying numbers of distinguishable species ranging from 4-species to the N-species limit, where N is the number of particles in the system itself. We examine these systems using various powerful simulation techniques that belong to the quantum Monte Carlo (QMC) family. The N-species limit is studied via Path Integral Monte Carlo (PIMC), which is a finite temperature, non- perturbative technique. At this limit, we explore the viability of two different thermal density matrices for systems that interact via hard-sphere, hard-cavity potentials. We then study the impact of finite-size effects on calculations of the energy, pressure and specific heat of the system when reaching the thermodynamic limit under periodic boundary conditions. We then move to studying a 4-species fermionic system using Variational Monte Carlo (VMC) and Diffusion Monte Carlo (DMC), which are ground state techniques. We apply these methods to the clustering problem of 4-particle and 8-particle systems. While the 4-particle, 4-species system is relatively simple to calculate as it is physically identical to a bosonic system of 4-particles, the 8-particle fermi system is more complicated. We detail the process of simulating a bound 8- particle state with respect to decay into two independent 4-particle clusters with DMC for the first time.
    [Show full text]
  • Deep Learning-Enhanced Variational Monte Carlo Method for Quantum Many-Body Physics
    PHYSICAL REVIEW RESEARCH 2, 012039(R) (2020) Rapid Communications Deep learning-enhanced variational Monte Carlo method for quantum many-body physics Li Yang ,1,* Zhaoqi Leng,2 Guangyuan Yu,3 Ankit Patel,4,5 Wen-Jun Hu,1,† and Han Pu1,‡ 1Department of Physics and Astronomy & Rice Center for Quantum Materials, Rice University, Houston, Texas 77005, USA 2Department of Physics, Princeton University, Princeton, New Jersey 08544, USA 3Department of Physics and Astronomy & The Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, USA 4Department of Electrical and Computer Engineering, Rice University, Houston, Texas 77005, USA 5Department of Neuroscience, Baylor College of Medicine, Houston, Texas 77030, USA (Received 26 May 2019; accepted 12 December 2019; published 14 February 2020) Artificial neural networks have been successfully incorporated into the variational Monte Carlo method (VMC) to study quantum many-body systems. However, there have been few systematic studies exploring quantum many-body physics using deep neural networks (DNNs), despite the tremendous success enjoyed by DNNs in many other areas in recent years. One main challenge of implementing DNNs in VMC is the inefficiency of optimizing such networks with a large number of parameters. We introduce an importance sampling gradient optimization (ISGO) algorithm, which significantly improves the computational speed of training DNNs by VMC. We design an efficient convolutional DNN architecture to compute the ground state of a one-dimensional SU(N) spin chain. Our numerical results of the ground-state energies with up to 16 layers of DNNs show excellent agreement with the Bethe ansatz exact solution. Furthermore, we also calculate the loop correlation function using the wave function obtained.
    [Show full text]
  • Variational Quantum Monte Carlo Abstract Introduction
    Variational Quantum Monte Carlo Josh Jiang, The University of Chicago > (1.1) Abstract The many electron system is one that quantum mechanics have no exact solution for because of the correlation between each electron with every other electron. Many methods have been developed to account for this correlation. One family of methods known as quantum Monte Carlo methods takes a numerical approach to solving the many-body problem. By employing a stochastic process to estimate the energy of a system, this method can be used to calculate the ground state energy of a system using the variational principle. The core of this method involves choosing an ansatz with a set of variational parameters that can be optimized to minimize the energy of the system. Recent research in this field is focused on developing methods to faster and more efficiently optimize this trial ansatz and also to make the ansatz more flexible by allowing more parameters. These topics along with the theory behind variational Monte Carlo (VMC) will be the subject of this Maple worksheet. Introduction Quantum mechanics is the study of particles at a very small scale. At this scale, particles can be characterized by wave functions and behave much differently than classical systems. By understanding how molecules behave at the quantum level, researchers can apply that knowledge to gain insights into how macroscopic systems work. While the applications are promising, quantum mechanics still struggles to understand the behavior of large atomic systems. For small, simple systems like particle in a box, harmonic oscillator, and the hydrogen atom, the Schrodinger equation: can be solved exactly, revealing the energy levels and exact wave functions for those systems.
    [Show full text]
  • Variational Monte Carlo Method for Fermions
    Variational Monte Carlo Method for Fermions Vijay B. Shenoy [email protected] Indian Institute of Science, Bangalore Mahabaleswar, December 2009 1/36 Acknowledgements Sandeep Pathak 2/36 Acknowledgements Sandeep Pathak Department of Science and Technology 2/36 Plan of the Lectures on Variational Monte Carlo Method Motivation, Overview and Key Ideas (Lecture 1, VBS) 3/36 Plan of the Lectures on Variational Monte Carlo Method Motivation, Overview and Key Ideas (Lecture 1, VBS) Details of VMC (Lecture 2, VBS) 3/36 Plan of the Lectures on Variational Monte Carlo Method Motivation, Overview and Key Ideas (Lecture 1, VBS) Details of VMC (Lecture 2, VBS) Application to Superconductivity (Lecture 3, VBS) 3/36 Plan of the Lectures on Variational Monte Carlo Method Motivation, Overview and Key Ideas (Lecture 1, VBS) Details of VMC (Lecture 2, VBS) Application to Superconductivity (Lecture 3, VBS) Application to Bosons (Lecture 4, AP) 3/36 Plan of the Lectures on Variational Monte Carlo Method Motivation, Overview and Key Ideas (Lecture 1, VBS) Details of VMC (Lecture 2, VBS) Application to Superconductivity (Lecture 3, VBS) Application to Bosons (Lecture 4, AP) Application to Spin Liquids (Lecture 5, AP) 3/36 The Strongly Correlated Panorama... Fig. 1. Phase diagrams ofrepresentativemateri- als of the strongly cor- related electron family (notations are standard and details can be found in the original refer- ences). (A) Temperature versus hole density phase diagram of bilayer man- ganites (74), including several types of antiferro- magnetic (AF) phases, a ferromagnetic (FM) phase, and even a glob- ally disordered region at x 0 0.75.
    [Show full text]
  • Arxiv:1807.10633V1 [Physics.Chem-Ph] 27 Jul 2018
    Improved speed and scaling in orbital space variational monte carlo Iliya Sabzevari1 and Sandeep Sharma1, ∗ 1Department of Chemistry and Biochemistry, University of Colorado Boulder, Boulder, CO 80302, USA In this work, we introduce three algorithmic improvements to reduce the cost and improve the scal- ing of orbital space variational Monte Carlo (VMC). First, we show that by appropriately screening the one- and two-electron integrals of the Hamiltonian one can improve the efficiency of the al- gorithm by several orders of magnitude. This improved efficiency comes with the added benefit that the resulting algorithm scales as the second power of the system size O(N 2), down from the fourth power O(N 4). Using numerical results, we demonstrate that the practical scaling obtained is in fact O(N 1:5) for a chain of Hydrogen atoms, and O(N 1:2) for the Hubbard model. Second, we introduce the use of the rejection-free continuous time Monte Carlo (CTMC) to sample the determinants. CTMC is usually prohibitively expensive because of the need to calculate a large number of intermediates. Here, we take advantage of the fact that these intermediates are already calculated during the evaluation of the local energy and consequently, just by storing them one can use the CTCM algorithm with virtually no overhead. Third, we show that by using the adaptive stochastic gradient descent algorithm called AMSGrad one can optimize the wavefunction energies robustly and efficiently. The combination of these three improvements allows us to calculate the ground state energy of a chain of 160 hydrogen atoms using a wavefunction containing ∼ 2 × 105 variational parameters with an accuracy of 1 mEh/particle at a cost of just 25 CPU hours, which when split over 2 nodes of 24 processors each amounts to only about half hour of wall time.
    [Show full text]
  • Variational Monte Carlo Algorithm for Solving One Dimensional Harmonic Oscillator Problem the Journal of Mathematics and Comput
    Farshid Mehrdoust, Hossein Aminikhah, Mohammad Ghamgosar / TJMCS Vol .3 No.4 (2011) 376 - 381 The Journal of Mathematics and Computer Science Available online at http://www.TJMCS.com The Journal of Mathematics and Computer Science Vol .3 No.4 (2011) 376 - 381 Variational Monte Carlo algorithm for solving one dimensional harmonic oscillator problem Farshid Mehrdoust Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Guilan, Rasht, Iran, [email protected] Hossein Aminikhah Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Guilan, Rasht, Iran, [email protected] Mohammad Ghamgosar [email protected] Received: September 2011, Revised: November 2011 Online Publication: December 2011 Abstract The goal of this paper is to present an application of variational Monte Carlo method for solving one dimensional harmonic oscillator problem. Keywords: Variational Monte Carlo (VMC); Harmonic oscillator; quantum mechanical systems; Statistical physics. AMS Subject Classification: 11k45 1. Introduction A basic problem in quantum mechanics is to find the energy eigenstates of quantum mechanical systems. The simplest such system is the quantum harmonic oscillator in one dimension. The Monte Carlo method is commonly used in physics to simulate complex systems those have 376 Farshid Mehrdoust, Hossein Aminikhah, Mohammad Ghamgosar / TJMCS Vol .3 No.4 (2011) 376 - 381 random nature in statistical physics [1]. 2. Variational Monte Carlo The variational method is a powerful tool to estimate the ground state energy of a quantum system, and some of its low lying excited states, when an exact solution cannot be found.When combined with Monte Carlo methods, it can be used on a wide range of atomic and solid state systems, as discussed in the review article by Foulkes et al [2].
    [Show full text]