DOCSLIB.ORG

Establishing Quantum Monte Carlo and Hybrid Density Functional Theory As Benchmarking Tools for Complex Solids

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Establishing Quantum Monte Carlo and Hybrid Density Functional Theory As Benchmarking Tools for Complex Solids Establishing Quantum Monte Carlo and Hybrid Density Functional Theory as benchmarking tools for complex solids DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Kevin P. Driver, B.S., M.S. Graduate Program in Physics The Ohio State University 2011 Dissertation Committee: John W. Wilkins, Advisor Richard J. Furnstahl Ciriyam Jayaprakash Arthur J. Epstein c Copyright by Kevin P. Driver 2011 Abstract Quantum mechanics provides an exact description of microscopic matter, but predictions require a solution of the fundamental many-electron Schr¨odingerequation. Since an ex- act solution of Schr¨odinger'sequation is intractable, several numerical methods have been developed to obtain approximate solutions. Currently, the two most successful methods are density functional theory (DFT) and quantum Monte Carlo (QMC). DFT is an exact theory which, which states that ground-state properties of a material can be obtained based on functionals of charge density alone. QMC is stochastic method which explicitly solves the many-body equation. In practice, the DFT method has drawbacks due to the fact that the exchange-correlation functional is not known. A large number of approximate exchange-correlation functionals have been produced to accommodate for this deficiency. Conceptual systematic improve- ments known as \Jacob's Ladder" of functional approximations have been made to the stan- dard local density approximation (LDA) and generalized gradient approximation (GGA). The traditional functionals have many known failures, such as failing to predict band gaps, silicon defect energies, and silica phase transitions. The newer generation functionals in- cluding meta-GGAs and hybrid functionals, such as the screened hybrid, HSE, have been developed to try to improve the flaws of lower-rung functionals. Overall, approximate functionals have generally had much success, but all functionals unpredictably vary in the quality and consistency of their predictions. Often, a failure of one type of DFT functional can be fixed by simply identifying another DFT functional that best describes the system under study. Identifying the best functional for the job is a challenging task, particularly if there is no experimental measurement to ii compare against. Higher accuracy methods, such as QMC, which are vastly more compu- tationally expensive, can be used to benchmark DFT functionals and identify those which work best for a material when experiment is lacking. If no DFT functional can perform adequately, then it is important to show more rigorous methods are capable of handling the task. QMC is high accuracy alternative to DFT, but QMC is too computationally expensive to replace DFT. Hybrid DFT functionals appear to be a good compromise between QMC and standard DFT. Not many large scale computations have been done to test the feasibility or benchmark capability of either QMC or hybrid DFT for complex materials. This thesis presents three applications expanding the scope of QMC and hybrid DFT to large, scale complex materials. QMC computes accurate formation energies for single-, di-, and tri- silicon-self-interstitials. QMC combined with phonon energies from DFT provide the most accurate equations of state, phase boundaries, and elastic properties available for silica. The HSE DFT functional is shown to reproduce QMC results for both silicon defects and high pressure silica phases, establishing its benchmark accuracy compared to other functionals. Standard DFT is still the most efficient and useful for general computation. However, this thesis shows that QMC and hybrid DFT calculations can aid and evaluate shortcomings associated the exchange-correlation potential in DFT by offering a route to benchmark and improve reliability of standard, more efficient DFT predictions. iii To my family and friends for guidance, help, and love. iv Acknowledgments The research and underlying educational enlightenment represented by this thesis are most notably a product of the continuous support and encouragement of my advisor, John Wilkins. John provided impeccable guidance and maintained a high standard of excellence in developing my scientific career. Other than my advisor, a few others deserve specific mention for their critical guidance and support. Richard Hennig was an excellent mentor and source of scientific inspiration through out my entire graduate career. My office mate and group partner, William Parker, provided invaluable amounts of feedback and support throughout my entire graduate career as well. I am also indebted to the help and guidance of Ronald Cohen, whose training made significant portions of this research (silica) possible. There are many other people have taught me or played some role in my scientific ed- ucation. I would like to thank Cyrus Umrigar, Burkhard Militzer, Hyoungki Park, Amita Wadhera, Mike Fellinger, Jeremy Nicklas, Ken Esler, Neil Drummond, Yaojun Du, Jeong- nim Kim, Thomas Lenosky, Shi-Yu Wu, Chakram Jayanthi, David Brown, P. J. Ouseph, and my high school physics teacher { Robert Rollings for general support and advice. Many departmental staff members offered important assistance me in some manner while carrying out this work. I'd like to thank Trisch Longbrake, Shelly Palmer, Carla Allen, Tim Randles, Brian Dunlap, and John Heimaster. I owe much gratitude to my family: Gerald and Patricia and Silvia Driver, Diane and Gerald Link, Jeremy and Leah Driver, Betty and Tom Wells, Robert and Dorothy Driver, and Irene Muller. Thanks for all of the love and support, and the opportunities provided that made my academic career possible. v I also want to acknowledge many important friends that have helped me personally and/or academically persevere in somewhat of a chronological order: Roseanne Cheng, Chuck and Danna Pearsall, Jeff Stevens, Sheldon Bailey, Julia Young, Grayson Williams, Nick and Barbara Harmon, Jake and Nichole Knepper, Brandon and Ester Parks, Chad and Nikki Morris, James Morris, Charlie Ruggiero, Greg Mack, Yuhfen Lin, Becky Daskalova, Kevin Knobbe, Mark and Sara Murphey, Matt Fisher, Yi Yang, Louis Nemzer, Kaden Haz- zard, Shawn Walsh, Justin Link, Chen Zhao, Qiu Weihong, Jia Chen, David Daughton, Kent Qian, Eric Jurgenson, Daniel Clark, Taeyoung Choi, Kerry Highbarger, Anthony Link, Emily Sistrunk, Mike Boss, Mike Hinton, Fred Kuehn, Iulian Hetel, Jen White, Va- lerie Bossow, Jim Potashnik, Deniz Duman, Greg Sollenberger, Patrick Smith, Anastasios Taliotis, Steven Avery, Aaron Sander, Eric Cochran, Hayes Lara, Adam Hauser, Luke Corwin, Srividya Iyer Biswas, Colin Schisler, Borun Chowdhury, Reni Ayachitula, Neesha Anderson, Rakesh Tivari, Nicole De Brabandere, Jim Davis, Rob Guidry, Lee Mosbacker, Don Burdette, Mehul Dixit, Dave Gohlke, Alex Mooney, Greg Vieira, James and Veronica Stapleton, John Draskovic, Mariko Mizuno, Yuval DaYu, Emily Harkins, Sabine Shaikh, Angie Detrow, The two Ashers, The HCGs, Meghan Ruck, Chiaki Ishikawa, April Brown, Kim Pabilona, Eumie Carter, Liesen Parkus, Cassandra Plummer, Nadia Ahmad, Natalie, Emma Brownlee, Claudia Veltze, Savannah Laurel-Zerr, Tom Steele, Alex Gray, Michelle Oglesbee, Kimberly Rousseau, Carlos Rubio, JC Polanco and all my friends from La Fogata, Patrick Roach, Heather Doughty, and many more whose names I've forgotten. I also have much appreciation for several financial agencies that supported me and my work. I was supported for two years at The Ohio State physics department as a Fowler fellow and further supported mostly by the DOE. I'd also like to thank the NSF for supporting my stay at the Carnegie Institution of Washington at the Geophysical Laboratory during the summers of 2007 and 2008. This work was also made possible by generous computational resources from OSC, NERSC, NCSA, and CCNI. vi Vita April 14, 1980 . Born|New Albany, Indiana, USA 2003 . B.S., University of Louisville, Louisville, Kentucky 2003-2005 . .Fowler Fellow, Department of Physics, Ohio State University, Columbus, Ohio 2006 . M.S., Ohio State University, Columbus, Ohio 2005-Present . Graduate Research Associate, Depart- ment of Physics, Ohio State University, Columbus, Ohio Publications K. P. Driver, R. E. Cohen, Zhigang Wu, B. Militzer, P. L´opez R´ıos,M. D. Towler, R. J. Needs, and J. W. Wilkins, Quantum Monte Carlo computations of phase stability, equations of state, and elasticity of high-pressure silica, Proc. Natl. Acad. Sci. USA, 107, 9519 (2010). R. G. Hennig, A. Wadehra, K. P. Driver, W. D. Parker, C. J. Umrigar, and J. W. Wilkins, Phase transformation in Si from semiconducting diamond to metallic beta-Sn phase in QMC and DFT under hydrostatic and anisotropic stress, Phys. Rev. B, 82, 014101 (2010). M. Floyd, Y. Zhang, K. P. Driver, Jeff Drucker, P.A. Crozier, and D.J. Smith, Nanometer- scale composition variations in Ge/Si(100) islands, Appl. Phys. Lett. 82, 1473 (2003). Y. Zhang, M. Floyd, K. P. Driver, Jeff Drucker, P.A. Crozier, and D.J. Smith, Evolution of Ge/Si(100) island morphology at high temperature, Appl. Phys. Lett. 80, 3623 (2002). P. J. Ouseph, K. P. Driver, J. Conklin, Polarization of Light By Reflection and the Brewster Angle, Am. J. Phys. 69, 1166 (2001). vii Fields of Study Major Field: Physics Studies in quantum Monte Carlo claculations of solids: J. W. Wilkins viii Table of Contents Page Abstract . ii Dedication . iv Acknowledgments . v Vita . vii List of Figures ...................................... xii List of Tables ....................................... xiv Chapters 1 Introduction
Load more

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]
  • Quantum Monte Carlo Analysis of Exchange and Correlation in the Strongly Inhomogeneous Electron Gas
    VOLUME 87, NUMBER 3 PHYSICAL REVIEW LETTERS 16JULY 2001 Quantum Monte Carlo Analysis of Exchange and Correlation in the Strongly Inhomogeneous Electron Gas Maziar Nekovee,1,* W. M. C. Foulkes,1 and R. J. Needs2 1The Blackett Laboratory, Imperial College, Prince Consort Road, London SW7 2BZ, United Kingdom 2Cavendish Laboratory, Madingley Road, Cambridge CB3 0HE, United Kingdom (Received 31 July 2000; published 25 June 2001) We use the variational quantum Monte Carlo method to calculate the density-functional exchange- correlation hole nxc, the exchange-correlation energy density exc, and the total exchange-correlation energy Exc of several strongly inhomogeneous electron gas systems. We compare our results with the local density approximation and the generalized gradient approximation. It is found that the nonlocal contributions to exc contain an energetically significant component, the magnitude, shape, and sign of which are controlled by the Laplacian of the electron density. DOI: 10.1103/PhysRevLett.87.036401 PACS numbers: 71.15.Mb, 71.10.–w, 71.45.Gm The Kohn-Sham density-functional theory (DFT) [1] rs ෇ 2a0 (approximately the same as for Al). The strong shows that it is possible to calculate the ground-state prop- density modulations were one dimensional and periodic, erties of interacting many-electron systems by solving only with a roughly sinusoidal profile. Three different modula- 0 0 one-electron Schrödinger-like equations. The results are tion wave vectors q # 2.17kF were investigated, where kF exact in principle, but in practice it is necessary to ap- is the Fermi wave vector corresponding to n0. The strong proximate the unknown exchange-correlation (XC) energy variation of n͑r͒ on the scale of the inverse local Fermi ͓ ͔ 21 2 21͞3 functional, Exc n , which expresses the many-body effects wave vector kF͑r͒ ෇ ͓3p n͑r͔͒ results in a strik- ͑ ͒ in terms of the electron density n r .
    [Show full text]
  • First Principles Investigation Into the Atom in Jellium Model System
    First Principles Investigation into the Atom in Jellium Model System Andrew Ian Duff H. H. Wills Physics Laboratory University of Bristol A thesis submitted to the University of Bristol in accordance with the requirements of the degree of Ph.D. in the Faculty of Science Department of Physics March 2007 Word Count: 34, 000 Abstract The system of an atom immersed in jellium is solved using density functional theory (DFT), in both the local density (LDA) and self-interaction correction (SIC) approxima- tions, Hartree-Fock (HF) and variational quantum Monte Carlo (VQMC). The main aim of the thesis is to establish the quality of the LDA, SIC and HF approximations by com- paring the results obtained using these methods with the VQMC results, which we regard as a benchmark. The second aim of the thesis is to establish the suitability of an atom in jellium as a building block for constructing a theory of the full periodic solid. A hydrogen atom immersed in a finite jellium sphere is solved using the methods listed above. The immersion energy is plotted against the positive background density of the jellium, and from this curve we see that DFT over-binds the electrons as compared to VQMC. This is consistent with the general over-binding one tends to see in DFT calculations. Also, for low values of the positive background density, the SIC immersion energy gets closer to the VQMC immersion energy than does the LDA immersion energy. This is consistent with the fact that the electrons to which the SIC is applied are becoming more localised at these low background densities and therefore the SIC theory is expected to out-perform the LDA here.
    [Show full text]
  • Ab Initio Molecular Dynamics of Water by Quantum Monte Carlo
    Ab initio molecular dynamics of water by quantum Monte Carlo Author: Ye Luo Supervisor: Prof. Sandro Sorella A thesis submitted for the degree of Doctor of Philosophy October 2014 Acknowledgments I would like to acknowledge first my supervisor Prof. Sandro Sorella. During the years in SISSA, he guided me in exploring the world of QMC with great patience and enthusiasm. He's not only a master in Physics but also an expert in high performance computing. He never exhausts new ideas and is always advancing the research at the light speed. For four years, I had an amazing journey with him. I would like to appreciate Dr. Andrea Zen, Prof. Leonardo Guidoni and my classmate Guglielmo Mazzola. I really learned a lot from the collaboration on various projects we did. I would like to thank Michele Casula and W.M.C. Foulkes who read and correct my thesis and also give many suggestive comments. I would like to express the gratitude to my parents who give me unlimited support even though they are very far from me. Without their unconditional love, I can't pass through the hardest time. I feel very sorry for them because I went home so few times in the previous years. Therefore, I would like to dedicate this thesis to them. I remember the pleasure with my friends in Trieste who are from all over the world. Through them, I have access to so many kinds of food and culture. They help me when I meet difficulties and they wipe out my loneliness by sharing the joys and tears of my life.
    [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]
  • Quantum Monte Carlo Methods on Lattices: the Determinantal Approach
    John von Neumann Institute for Computing Quantum Monte Carlo Methods on Lattices: The Determinantal Approach Fakher F. Assaad published in Quantum Simulations of Complex Many-Body Systems: From Theory to Algorithms, Lecture Notes, J. Grotendorst, D. Marx, A. Muramatsu (Eds.), John von Neumann Institute for Computing, Julich,¨ NIC Series, Vol. 10, ISBN 3-00-009057-6, pp. 99-156, 2002. c 2002 by John von Neumann Institute for Computing Permission to make digital or hard copies of portions of this work for personal or classroom use is granted provided that the copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. To copy otherwise requires prior specific permission by the publisher mentioned above. http://www.fz-juelich.de/nic-series/volume10 Quantum Monte Carlo Methods on Lattices: The Determinantal Approach Fakher F. Assaad 1 Institut fur¨ Theoretische Physik III, Universitat¨ Stuttgart Pfaffenwaldring 57, 70550 Stuttgart, Germany 2 Max Planck Institute for Solid State Research Heisenbergstr. 1, 70569, Stuttgart, Germany E-mail: [email protected] We present a review of the auxiliary field (i.e. determinantal) Quantum Monte Carlo method applied to various problems of correlated electron systems. The ground state projector method, the finite temperature approach as well as the Hirsch-Fye impurity algorithm are described in details. It is shown how to apply those methods to a variety of models: Hubbard Hamiltonians, periodic Anderson model, Kondo lattice and impurity problems, as well as hard core bosons and the Heisenberg model. An introduction to the world-line method with loop upgrades as well as an appendix on the Monte Carlo method is provided.
    [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]
  • 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]
  • Science Petascale
    SCIENCE at the PETASCALE 2009 Pioneering Applications Point the Way ORNL 2009-G00700/jbs CONTENTS 1 Introduction DOMAINS 3 Biology 5 Chemistry 9 Climate 13 Combustion 15 Fusion 21 Geoscience 23 Materials Science 33 Nuclear Energy 37 Nuclear Physics 41 Turbulence APPENDICES 43 Specifications 44 Applications Index 46 Acronyms INTRODUCTION Our early science efforts were a resounding success. As the program draws to a close, nearly 30 leading research teams in climate science, energy assurance, and fundamental science will have used more than 360 million processor-hours over these 6 months. Their applications were able to take advantage of the 300-terabyte, 150,000-processor Cray XT5 partition of the Jaguar system to conduct research at a scale and accuracy that would be literally impossible on any other system in existence. The projects themselves were chosen from virtually every science domain; highlighted briefly here is research in weather and climate, nuclear energy, ow is an especially good time to be engaged geosciences, combustion, bioenergy, fusion, and in computational science research. After years materials science. Almost two hundred scientists Nof planning, developing, and building, we have flocked to Jaguar during this period from dozens of entered the petascale era with Jaguar, a Cray XT institutions around the world. supercomputer located at the Oak Ridge Leadership Computing Facility (OLCF). • Climate and Weather Modeling. One team used century-long simulations of atmospheric and By petascale we mean a computer system able to churn ocean circulations to illuminate their role in climate out more than a thousand trillion calculations each variability. Jaguar’s petascale power enabled second—or a petaflop.
    [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]
  • The Early Years of Quantum Monte Carlo (2): Finite- Temperature Simulations
    The Early Years of Quantum Monte Carlo (2): Finite- Temperature Simulations Michel Mareschal1,2, Physics Department, ULB , Bruxelles, Belgium Introduction This is the second part of the historical survey of the development of the use of random numbers to solve quantum mechanical problems in computational physics and chemistry. In the first part, noted as QMC1, we reported on the development of various methods which allowed to project trial wave functions onto the ground state of a many-body system. The generic name for all those techniques is Projection Monte Carlo: they are also known by more explicit names like Variational Monte Carlo (VMC), Green Function Monte Carlo (GFMC) or Diffusion Monte Carlo (DMC). Their different names obviously refer to the peculiar technique used but they share the common feature of solving a quantum many-body problem at zero temperature, i.e. computing the ground state wave function. In this second part, we will report on the technique known as Path Integral Monte Carlo (PIMC) and which extends Quantum Monte Carlo to problems with finite (i.e. non-zero) temperature. The technique is based on Feynman’s theory of liquid Helium and in particular his qualitative explanation of superfluidity [Feynman,1953]. To elaborate his theory, Feynman used a formalism he had himself developed a few years earlier: the space-time approach to quantum mechanics [Feynman,1948] was based on his thesis work made at Princeton under John Wheeler before joining the Manhattan’s project. Feynman’s Path Integral formalism can be cast into a form well adapted for computer simulation. In particular it transforms propagators in imaginary time into Boltzmann weight factors, allowing the use of Monte Carlo sampling.
    [Show full text]