Daubechies Wavelets As a Basis Set for Density Functional

Total Page:16

File Type:pdf, Size:1020Kb

Daubechies Wavelets As a Basis Set for Density Functional Electronic Structure Calculation Methods on Accelerators, Oak Ridge BigDFT Formalism Wavelets Algorithms Daubechies wavelets as a basis set for density Poisson Solver Features functional pseudopotential calculation Applications Publications Si clusters Boron clusters Conclusion T. Deutsch, L. Genovese, D. Caliste L_Sim – CEA Grenoble February 7, 2012 Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste A basis for nanosciences: The BigDFT project STREP European project: BigDFT(2005-2008) Four partners, 15 contributors: CEA-INAC Grenoble (T. D.), U. Basel (S. Goedecker), U. Louvain-la-Neuve (X .Gonze), U. Kiel (R. Schneider) BigDFT Formalism Aim: To develop an ab-initio DFT code Wavelets based on Daubechies Wavelets Algorithms Poisson Solver and integrated in ABINIT (Gnu-GPL). Features Applications BigDFT 1.0 −! January 2008 Publications Si clusters Boron clusters Conclusion . Why have we done this? Test the potential advantages of a new formalism * A lot of outcomes and interesting results A lot can be done starting from present know-how Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Outline 1 Formalism Wavelets Algorithms Poisson Solver BigDFT Features Formalism Wavelets Algorithms 2 Applications Poisson Solver Features Publications Applications Si clusters Publications Si clusters Boron clusters Boron clusters Conclusion 3 Conclusion Next talk: Performances (MPI+OpenMP+OpenCL) Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste A DFT code based on Daubechies wavelets BigDFT: a PSP Kohn-Sham code A Daubechies wavelets basis has unique properties for DFT usage Systematic, Orthogonal BigDFT Localised, Adaptive Formalism Wavelets Algorithms Kohn-Sham operators are analytic Poisson Solver Daubechies Wavelets Features Short, Separable convolutions 1.5 Applications φ(x) Publications ~c = a c ` ∑j j `−j 1 ψ(x) Si clusters Boron clusters 0.5 Peculiar numerical properties Conclusion 0 Real space based, highly flexible -0.5 -1 Big & inhomogeneous systems -1.5 -6 -4 -2 0 2 4 6 8 x Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste A brief description of wavelet theory Two kind of basis functions A Multi-Resolution real space basis The functions can be classified following the resolution level they span. BigDFT Formalism Scaling Functions Wavelets Algorithms The functions of low resolution level are a linear combination Poisson Solver Features of high-resolution functions Applications Publications Si clusters = + Boron clusters m Conclusion f(x) = ∑ hj f(2x − j) j=−m Centered on a resolution-dependent grid: fj = f0(x − j). Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste A brief description of wavelet theory Wavelets They contain the DoF needed to complete the information which is lacking due to the coarseness of the resolution. = 1 + 1 BigDFT 2 2 m m Formalism (2x) = h~ (x − j) + g~ (x − j) Wavelets f ∑ j f ∑ j y Algorithms j=−m j=−m Poisson Solver Features Applications Increase the resolution without modifying grid space Publications Si clusters Boron clusters SF + W = same DoF of SF of higher resolution Conclusion m y(x) = ∑ gj f(2x − j) j=−m All functions have compact support, centered on grid points. Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Wavelet properties: Adaptivity Adaptivity Resolution can be refined following the grid point. The grid is divided in BigDFT Low (1 DoF) and High Formalism (8 DoF) resolution points. Wavelets Algorithms Points of different resolution Poisson Solver belong to the same grid. Features Applications Empty regions must not be Publications Si clusters “filled” with basis functions. Boron clusters Conclusion Localization property,real space description Optimal for big & inhomogeneous systems, highly flexible Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Adaptivity of the mesh Atomic positions (H2O), hgrid BigDFT Formalism Wavelets Algorithms Poisson Solver Features Applications Publications Si clusters Boron clusters Conclusion Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Adaptivity of the mesh Fine grid (high resolution, frmult) BigDFT Formalism Wavelets Algorithms Poisson Solver Features Applications Publications Si clusters Boron clusters Conclusion Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Adaptivity of the mesh Coarse grid (low resolution, crmult) BigDFT Formalism Wavelets Algorithms Poisson Solver Features Applications Publications Si clusters Boron clusters Conclusion Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste No integration error Orthogonality, scaling relation Z 1 m dx fk (x)fj (x) = dkj f(x) = p ∑ hj f(2x − j) 2 j=−m BigDFT The hamiltonian-related quantities can be calculated up to Formalism machine precision in the given basis. Wavelets Algorithms The accuracy is only limited by the basis set ( h14 ) Poisson Solver O grid Features Applications Publications Exact evaluation of kinetic energy Si clusters Boron clusters Obtained by convolution with filters: Conclusion 2 f (x) = ∑c`f`(x); ∇ f (x) = ∑~c` f`(x) ; ` ` Z 2 ~c` = ∑cj a`−j ; a` ≡ f0(x)¶x f`(x) ; j Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Separability in 3D The 3-dim scaling basis is a tensor product decomposition of 1-dim Scaling Functions/ Wavelets. e ;e ;e e f x y z (x;y;z) = fex (x)f y (y)fez (z) jx ;jy ;jz jx jy jz BigDFT With (jx ;jy ;jz ) the node coordinates, (0) (1) Formalism fj and fj the SF and the W respectively. Wavelets Algorithms Poisson Solver Features Gaussians and wavelets Applications Publications The separability of the basis allows us to save computational Si clusters Boron clusters time when performing scalar products with separable Conclusion functions (e.g. gaussians): Initial wavefunctions (input guess) 1 Poisson solver (tensor decomposition of r ) Non-local pseudopotentials Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Wavelet families used in BigDFT code Daubechies f (x) = ∑` c`f`(x) Interpolating f (x) = ∑j fj jj (x) R Orthogonal c` = dx f`(x)f (x) Dual to dirac deltas fj = f (j) LEAST ASYMMETRIC DAUBECHIES-16 1.5 INTERPOLATIN !"E#LAURIERS-"$%$&'-( scaling function 1 wavelet scaling function wavelet 1 0.5 0.5 BigDFT 0 0 Formalism -0.5 Wavelets -0.5 Algorithms -1 Poisson Solver Features -1.5 -1 -6 -4 -2 0 2 4 6 -4 -2 0 2 4 Applications Used for wavefunctions, Used for charge density, Publications Si clusters scalar products function products Boron clusters Conclusion “Magic Filter” method (A. Neelov, S. Goedecker) The passage between the two basis sets can be performed without losing accuracy. Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Operations performed The SCF cycle Real Space Daub. Wavelets Orbital scheme: Hamiltonian Preconditioner BigDFT Coefficient Scheme: Formalism Wavelets Overlap matrices Algorithms Poisson Solver Features Orthogonalisation Applications Publications Comput. operations Si clusters Boron clusters Convolutions Conclusion BLAS routines FFT (Poisson Solver) Why not GPUs? Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Kohn-Sham Equations: Computing Energies Calculate different integrals E[r] = K [r] + U[r] BigDFT 2 Z 2 ~ ∗ 2 Formalism K [r] = − ∑ dr yi ∇ yi Wavelets 2 me i V Algorithms Poisson Solver Features Z Z 0 1 0 r(r)r(r ) Applications U[r] = dr Vext (r)r(r)+ dr dr 0 + Exc[r] Publications V 2 V jr − r j | {z } Si clusters | {z } exchange−correlation Boron clusters Hartree Conclusion We minimise with the variables yi (r) and fi Z with the constraint dr r(r) = Nel . V Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste inv FWT + Magic Filter occ 2 r(r) = ∑ jyi (r)j Fine non-adaptive mesh: < 8Nf i Poisson solver V (r) = R G(:)r Vxc [r(r)] H VNL(fyi g) Veffective FWT 1 2 − 2 ∇ FWT Kinetic Term i ¶Etotal j dca = − ∗ + ∑Λij ca ¶ci (a) j Λij =< yi jHjyj > i Stop when dca small preconditioning Steepest Descent, new;i i i ca = ca + hstepdca DIIS Direct Minimisation: Flowchart i yi = ∑cafa Basis: adaptive mesh (0;:::;Nf) a Orthonormalized BigDFT Formalism Wavelets Algorithms Poisson Solver Features Applications Publications Si clusters Boron clusters Conclusion Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste Poisson solver V (r) = R G(:)r Vxc [r(r)] H VNL(fyi g) Veffective FWT 1 2 − 2 ∇ FWT Kinetic Term i ¶Etotal j dca = − ∗ + ∑Λij ca ¶ci (a) j Λij =< yi jHjyj > i Stop when dca small preconditioning Steepest Descent, new;i i i ca = ca + hstepdca DIIS Direct Minimisation: Flowchart i yi = ∑cafa Basis: adaptive mesh (0;:::;Nf) a Orthonormalized inv FWT + Magic Filter occ 2 BigDFT r(r) = ∑ jyi (r)j Fine non-adaptive mesh: < 8Nf i Formalism Wavelets Algorithms Poisson Solver Features Applications Publications Si clusters Boron clusters Conclusion Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T. Deutsch, L. Genovese, D. Caliste FWT 1 2 − 2 ∇ FWT Kinetic Term i ¶Etotal j dca = − ∗ + ∑Λij ca ¶ci (a) j Λij =< yi jHjyj > i Stop when dca small preconditioning Steepest Descent, new;i i i ca = ca + hstepdca DIIS Direct Minimisation: Flowchart i yi = ∑cafa Basis: adaptive mesh (0;:::;Nf) a Orthonormalized inv FWT + Magic Filter occ 2 BigDFT r(r) = ∑ jyi (r)j Fine non-adaptive mesh: < 8Nf i Formalism Wavelets Algorithms Poisson solver Poisson Solver Features V (r) = R G(:)r Applications Vxc [r(r)] H VNL(fyi g) Veffective Publications Si clusters Boron clusters Conclusion Laboratoire de Simulation Atomistique http://inac.cea.fr/L_Sim T.
Recommended publications
  • Free and Open Source Software for Computational Chemistry Education
    Free and Open Source Software for Computational Chemistry Education Susi Lehtola∗,y and Antti J. Karttunenz yMolecular Sciences Software Institute, Blacksburg, Virginia 24061, United States zDepartment of Chemistry and Materials Science, Aalto University, Espoo, Finland E-mail: [email protected].fi Abstract Long in the making, computational chemistry for the masses [J. Chem. Educ. 1996, 73, 104] is finally here. We point out the existence of a variety of free and open source software (FOSS) packages for computational chemistry that offer a wide range of functionality all the way from approximate semiempirical calculations with tight- binding density functional theory to sophisticated ab initio wave function methods such as coupled-cluster theory, both for molecular and for solid-state systems. By their very definition, FOSS packages allow usage for whatever purpose by anyone, meaning they can also be used in industrial applications without limitation. Also, FOSS software has no limitations to redistribution in source or binary form, allowing their easy distribution and installation by third parties. Many FOSS scientific software packages are available as part of popular Linux distributions, and other package managers such as pip and conda. Combined with the remarkable increase in the power of personal devices—which rival that of the fastest supercomputers in the world of the 1990s—a decentralized model for teaching computational chemistry is now possible, enabling students to perform reasonable modeling on their own computing devices, in the bring your own device 1 (BYOD) scheme. In addition to the programs’ use for various applications, open access to the programs’ source code also enables comprehensive teaching strategies, as actual algorithms’ implementations can be used in teaching.
    [Show full text]
  • The Impact of Density Functional Theory on Materials Research
    www.mrs.org/bulletin functional. This functional (i.e., a function whose argument is another function) de- scribes the complex kinetic and energetic interactions of an electron with other elec- Toward Computational trons. Although the form of this functional that would make the reformulation of the many-body Schrödinger equation exact is Materials Design: unknown, approximate functionals have proven highly successful in describing many material properties. Efficient algorithms devised for solving The Impact of the Kohn–Sham equations have been imple- mented in increasingly sophisticated codes, tremendously boosting the application of DFT methods. New doors are opening to in- Density Functional novative research on materials across phys- ics, chemistry, materials science, surface science, and nanotechnology, and extend- ing even to earth sciences and molecular Theory on Materials biology. The impact of this fascinating de- velopment has not been restricted to aca- demia, as DFT techniques also find application in many different areas of in- Research dustrial research. The development is so fast that many current applications could Jürgen Hafner, Christopher Wolverton, and not have been realized three years ago and were hardly dreamed of five years ago. Gerbrand Ceder, Guest Editors The articles collected in this issue of MRS Bulletin present a few of these suc- cess stories. However, even if the compu- Abstract tational tools necessary for performing The development of modern materials science has led to a growing need to complex quantum-mechanical calcula- tions relevant to real materials problems understand the phenomena determining the properties of materials and processes on are now readily available, designing a an atomistic level.
    [Show full text]
  • Practice: Quantum ESPRESSO I
    MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW FREE http://www.quantum-espresso.org PW-DFT, PP, PAW FREE http://www.abinit.org DFT PW, PP, Car-Parrinello FREE http://www.cpmd.org DFT PP, US-PP, PAW $3000 [moderate accuracy, fast] http://www.vasp.at DFT full-potential linearized augmented $500 plane-wave (FLAPW) [accurate, slow] http://www.wien2k.at Hartree-Fock, higher order correlated $3000 electron approaches http://www.gaussian.com 3 Quantum ESPRESSO 4 Quantum ESPRESSO Quantum ESPRESSO is an integrated suite of Open- Source computer codes for electronic-structure calculations and materials modeling at the nanoscale. It is based on density-functional theory, plane waves, and pseudopotentials. Core set of codes, plugins for more advanced tasks and third party packages Open initiative coordinated by the Quantum ESPRESSO Foundation, across Italy. Contributed to by developers across the world Regular hands-on workshops in Trieste, Italy Open-source code: FREE (unlike VASP...) 5 Performance Small jobs (a few atoms) can be run on single node Includes determining convergence parameters, lattice constants Can use OpenMP parallelization on multicore machines Large jobs (~10’s to ~100’s atoms) can run in parallel using MPI to 1000’s of cores Includes molecular dynamics, large geometry relaxation, phonons Parallel performance tied to BLAS/LAPACK (linear algebra routines) and 3D FFT (fast Fourier transform) New GPU-enabled version available 6 Usability Documented online:
    [Show full text]
  • The CECAM Electronic Structure Library and the Modular Software Development Paradigm
    The CECAM electronic structure library and the modular software development paradigm Cite as: J. Chem. Phys. 153, 024117 (2020); https://doi.org/10.1063/5.0012901 Submitted: 06 May 2020 . Accepted: 08 June 2020 . Published Online: 13 July 2020 Micael J. T. Oliveira , Nick Papior , Yann Pouillon , Volker Blum , Emilio Artacho , Damien Caliste , Fabiano Corsetti , Stefano de Gironcoli , Alin M. Elena , Alberto García , Víctor M. García-Suárez , Luigi Genovese , William P. Huhn , Georg Huhs , Sebastian Kokott , Emine Küçükbenli , Ask H. Larsen , Alfio Lazzaro , Irina V. Lebedeva , Yingzhou Li , David López- Durán , Pablo López-Tarifa , Martin Lüders , Miguel A. L. Marques , Jan Minar , Stephan Mohr , Arash A. Mostofi , Alan O’Cais , Mike C. Payne, Thomas Ruh, Daniel G. A. Smith , José M. Soler , David A. Strubbe , Nicolas Tancogne-Dejean , Dominic Tildesley, Marc Torrent , and Victor Wen-zhe Yu COLLECTIONS Paper published as part of the special topic on Electronic Structure Software Note: This article is part of the JCP Special Topic on Electronic Structure Software. This paper was selected as Featured ARTICLES YOU MAY BE INTERESTED IN Recent developments in the PySCF program package The Journal of Chemical Physics 153, 024109 (2020); https://doi.org/10.1063/5.0006074 An open-source coding paradigm for electronic structure calculations Scilight 2020, 291101 (2020); https://doi.org/10.1063/10.0001593 Siesta: Recent developments and applications The Journal of Chemical Physics 152, 204108 (2020); https://doi.org/10.1063/5.0005077 J. Chem. Phys. 153, 024117 (2020); https://doi.org/10.1063/5.0012901 153, 024117 © 2020 Author(s). The Journal ARTICLE of Chemical Physics scitation.org/journal/jcp The CECAM electronic structure library and the modular software development paradigm Cite as: J.
    [Show full text]
  • The Empirical Pseudopotential Method
    The empirical pseudopotential method Richard Tran PID: A1070602 1 Background In 1934, Fermi [1] proposed a method to describe high lying atomic states by replacing the complex effects of the core electrons with an effective potential described by the pseudo-wavefunctions of the valence electrons. Shortly after, this pseudopotential (PSP), was realized by Hellmann [2] who implemented it to describe the energy levels of alkali atoms. Despite his initial success, further exploration of the PSP and its application in efficiently solving the Hamiltonian of periodic systems did not reach its zenith until the middle of the twenti- eth century when Herring [3] proposed the orthogonalized plane wave (OPW) method to calculate the wave functions of metals and semiconductors. Using the OPW, Phillips and Kleinman [4] showed that the orthogonality of the va- lence electrons’ wave functions relative to that of the core electrons leads to a cancellation of the attractive and repulsive potentials of the core electrons (with the repulsive potentials being approximated either analytically or empirically), allowing for a rapid convergence of the OPW calculations. Subsequently in the 60s, Marvin Cohen expanded and improved upon the pseudopotential method by using the crystalline energy levels of several semicon- ductor materials [6, 7, 8] to empirically obtain the potentials needed to compute the atomic wave functions, developing what is now known as the empirical pseu- dopotential method (EPM). Although more advanced PSPs such as the ab-initio pseudopotentials exist, the EPM still provides an incredibly accurate method for calculating optical properties and band structures, especially for metals and semiconductors, with a less computationally taxing method.
    [Show full text]
  • Kepler Gpus and NVIDIA's Life and Material Science
    LIFE AND MATERIAL SCIENCES Mark Berger; [email protected] Founded 1993 Invented GPU 1999 – Computer Graphics Visual Computing, Supercomputing, Cloud & Mobile Computing NVIDIA - Core Technologies and Brands GPU Mobile Cloud ® ® GeForce Tegra GRID Quadro® , Tesla® Accelerated Computing Multi-core plus Many-cores GPU Accelerator CPU Optimized for Many Optimized for Parallel Tasks Serial Tasks 3-10X+ Comp Thruput 7X Memory Bandwidth 5x Energy Efficiency How GPU Acceleration Works Application Code Compute-Intensive Functions Rest of Sequential 5% of Code CPU Code GPU CPU + GPUs : Two Year Heart Beat 32 Volta Stacked DRAM 16 Maxwell Unified Virtual Memory 8 Kepler Dynamic Parallelism 4 Fermi 2 FP64 DP GFLOPS GFLOPS per DP Watt 1 Tesla 0.5 CUDA 2008 2010 2012 2014 Kepler Features Make GPU Coding Easier Hyper-Q Dynamic Parallelism Speedup Legacy MPI Apps Less Back-Forth, Simpler Code FERMI 1 Work Queue CPU Fermi GPU CPU Kepler GPU KEPLER 32 Concurrent Work Queues Developer Momentum Continues to Grow 100M 430M CUDA –Capable GPUs CUDA-Capable GPUs 150K 1.6M CUDA Downloads CUDA Downloads 1 50 Supercomputer Supercomputers 60 640 University Courses University Courses 4,000 37,000 Academic Papers Academic Papers 2008 2013 Explosive Growth of GPU Accelerated Apps # of Apps Top Scientific Apps 200 61% Increase Molecular AMBER LAMMPS CHARMM NAMD Dynamics GROMACS DL_POLY 150 Quantum QMCPACK Gaussian 40% Increase Quantum Espresso NWChem Chemistry GAMESS-US VASP CAM-SE 100 Climate & COSMO NIM GEOS-5 Weather WRF Chroma GTS 50 Physics Denovo ENZO GTC MILC ANSYS Mechanical ANSYS Fluent 0 CAE MSC Nastran OpenFOAM 2010 2011 2012 SIMULIA Abaqus LS-DYNA Accelerated, In Development NVIDIA GPU Life Science Focus Molecular Dynamics: All codes are available AMBER, CHARMM, DESMOND, DL_POLY, GROMACS, LAMMPS, NAMD Great multi-GPU performance GPU codes: ACEMD, HOOMD-Blue Focus: scaling to large numbers of GPUs Quantum Chemistry: key codes ported or optimizing Active GPU acceleration projects: VASP, NWChem, Gaussian, GAMESS, ABINIT, Quantum Espresso, BigDFT, CP2K, GPAW, etc.
    [Show full text]
  • The Plane-Wave Pseudopotential Method
    The Plane-wave Pseudopotential Method i(G+k) r k(r)= ck,Ge · XG Chris J Pickard Electrons in a Solid Nearly Free Electrons Nearly Free Electrons Nearly Free Electrons Electronic Structures Methods Empirical Pseudopotentials Ab Initio Pseudopotentials Ab Initio Pseudopotentials Ab Initio Pseudopotentials Ab Initio Pseudopotentials Ultrasoft Pseudopotentials (Vanderbilt, 1990) Projector Augmented Waves Deriving the Pseudo Hamiltonian PAW vs PAW C. G. van de Walle and P. E. Bloechl 1993 (I) PAW to calculate all electron properties from PSP calculations" P. E. Bloechl 1994 (II)! The PAW electronic structure method" " The two are linked, but should not be confused References David Vanderbilt, Soft self-consistent pseudopotentials in a generalised eigenvalue formalism, PRB 1990 41 7892" Kari Laasonen et al, Car-Parinello molecular dynamics with Vanderbilt ultrasoft pseudopotentials, PRB 1993 47 10142" P.E. Bloechl, Projector augmented-wave method, PRB 1994 50 17953" G. Kresse and D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, PRB 1999 59 1758 On-the-fly Pseudopotentials - All quantities for PAW (I) reconstruction available" - Allows automatic consistency of the potentials with functionals" - Even fewer input files required" - The code in CASTEP is a cut down (fast) pseudopotential code" - There are internal databases" - Other generation codes (Vanderbilt’s original, OPIUM, LD1) exist It all works! Kohn-Sham Equations Being (Self) Consistent Eigenproblem in a Basis Just a Few Atoms In a Crystal Plane-waves
    [Show full text]
  • Quantum Chemistry (QC) on Gpus Feb
    Quantum Chemistry (QC) on GPUs Feb. 2, 2017 Overview of Life & Material Accelerated Apps MD: All key codes are GPU-accelerated QC: All key codes are ported or optimizing Great multi-GPU performance Focus on using GPU-accelerated math libraries, OpenACC directives Focus on dense (up to 16) GPU nodes &/or large # of GPU nodes GPU-accelerated and available today: ACEMD*, AMBER (PMEMD)*, BAND, CHARMM, DESMOND, ESPResso, ABINIT, ACES III, ADF, BigDFT, CP2K, GAMESS, GAMESS- Folding@Home, GPUgrid.net, GROMACS, HALMD, HOOMD-Blue*, UK, GPAW, LATTE, LSDalton, LSMS, MOLCAS, MOPAC2012, LAMMPS, Lattice Microbes*, mdcore, MELD, miniMD, NAMD, NWChem, OCTOPUS*, PEtot, QUICK, Q-Chem, QMCPack, OpenMM, PolyFTS, SOP-GPU* & more Quantum Espresso/PWscf, QUICK, TeraChem* Active GPU acceleration projects: CASTEP, GAMESS, Gaussian, ONETEP, Quantum Supercharger Library*, VASP & more green* = application where >90% of the workload is on GPU 2 MD vs. QC on GPUs “Classical” Molecular Dynamics Quantum Chemistry (MO, PW, DFT, Semi-Emp) Simulates positions of atoms over time; Calculates electronic properties; chemical-biological or ground state, excited states, spectral properties, chemical-material behaviors making/breaking bonds, physical properties Forces calculated from simple empirical formulas Forces derived from electron wave function (bond rearrangement generally forbidden) (bond rearrangement OK, e.g., bond energies) Up to millions of atoms Up to a few thousand atoms Solvent included without difficulty Generally in a vacuum but if needed, solvent treated classically
    [Show full text]
  • Introduction to DFT and the Plane-Wave Pseudopotential Method
    Introduction to DFT and the plane-wave pseudopotential method Keith Refson STFC Rutherford Appleton Laboratory Chilton, Didcot, OXON OX11 0QX 23 Apr 2014 Parallel Materials Modelling Packages @ EPCC 1 / 55 Introduction Synopsis Motivation Some ab initio codes Quantum-mechanical approaches Density Functional Theory Electronic Structure of Condensed Phases Total-energy calculations Introduction Basis sets Plane-waves and Pseudopotentials How to solve the equations Parallel Materials Modelling Packages @ EPCC 2 / 55 Synopsis Introduction A guided tour inside the “black box” of ab-initio simulation. Synopsis • Motivation • The rise of quantum-mechanical simulations. Some ab initio codes Wavefunction-based theory • Density-functional theory (DFT) Quantum-mechanical • approaches Quantum theory in periodic boundaries • Plane-wave and other basis sets Density Functional • Theory SCF solvers • Molecular Dynamics Electronic Structure of Condensed Phases Recommended Reading and Further Study Total-energy calculations • Basis sets Jorge Kohanoff Electronic Structure Calculations for Solids and Molecules, Plane-waves and Theory and Computational Methods, Cambridge, ISBN-13: 9780521815918 Pseudopotentials • Dominik Marx, J¨urg Hutter Ab Initio Molecular Dynamics: Basic Theory and How to solve the Advanced Methods Cambridge University Press, ISBN: 0521898633 equations • Richard M. Martin Electronic Structure: Basic Theory and Practical Methods: Basic Theory and Practical Density Functional Approaches Vol 1 Cambridge University Press, ISBN: 0521782856
    [Show full text]
  • The Plane-Wave Pseudopotential Method
    The Plane-Wave Pseudopotential Method Eckhard Pehlke, Institut für Laser- und Plasmaphysik, Universität Essen, 45117 Essen, Germany. Starting Point: Electronic structure problem from physics, chemistry, materials science, geology, biology, ... which can be solved by total-energy calculations. Topics of this (i) how to get rid of the "core electrons": talk: the pseudopotential concept (ii) the plane-wave basis-set and its advantages (iii) supercells, Bloch theorem and Brillouin zone integrals choices we have... Treatment of electron-electron interaction and decisions to make... Hartree-Fock (HF) ¿ ÒÌ Ò· Ú ´ÖµÒ´Öµ Ö Ú ¼ ¼ ½ Ò´ÖµÒ´Ö µ ¿ ¿ ¼ · Ö Ö · Ò Configuration-Interaction (CI) ¼ ¾ Ö Ö Quantum Monte-Carlo (QMC) ÑÒ Ò ¼ Ú Ê ¿ Ò´Ö µ ÖÆ Density-Functional Theory (DFT) ¼ ½ Ò´Ö µ ¾ ¿ ¼ ´ Ö · Ú ´Öµ· Ö · Ú ´Öµµ ´Öµ ´Öµ ¼ ¾ Ö Ö Æ Ò ¾ Ò´Öµ ´Öµ Ú ´Öµ Æ´Öµ Problem: Approximation to XC functional. decisions to make... Simulation of Atomic Geometries Example: chemisorption site & energy of a particular atom on a surface = ? How to simulate adsorption geometry? single molecule or cluster Use Bloch theorem. periodically repeated supercell, slab-geometry Efficient Brillouin zone integration schemes. true half-space geometry, Green-function methods decisions to make... Basis Set to Expand Wave-Functions ½ ¡Ê ´Øµ ´ Öµ Ù ´Ö Ê µ Ô Æ linear combination of Ê atomic orbitals (LCAO) ´ Öµ ´µ ´ Öµ simple, unbiased, ½ ´·µ¡Ö ´ Öµ ´ · µ plane waves (PW) Ô independent of ª atomic positions augmented plane waves . (APW) etc. etc. I. The Pseudopotential Concept Core-States and Chemical Bonding? Validity of the Frozen-Core Approximation U. von Barth, C.D.
    [Show full text]
  • CPMD Tutorial Car–Parrinello Molecular Dynamics
    CPMD Tutorial Car–Parrinello Molecular Dynamics Ari P Seitsonen CNRS & Universit´ePierre et Marie Curie, Paris CSC, October 2006 Introduction CPMD — concept 1 • Car-Parrinello Molecular Dynamics – Roberto Car & Michele Parrinello, Physical Review Letters 55, 2471 (1985) — 20+1 years of Car-Parrinello method! Roberto Car Michele Parrinello CPMD — concept 2 • CPMD code — http://www.cpmd.org/; based on the original code of Roberto Car and Michele Parrinello CPMD What is it, what is it not? • Computer code for performing static, dynamic simulations and analysis of electronic structure • About 200’000 lines of code — FORTRAN77 with a flavour towards F90 • (Freely) available via http://www.cpmd.org/ with source code • Most suitable for dynamical simulations of condensed systems or large molecules — less for small molecules or bulk properties of small crystals • Computationally highly optimised for vector and scalar supercomputers, shared and distributed memory machines and combinations thereof (SMP) Capabilities • Solution of electronic and ionic structure • XC-functional: LDA, GGA, meta-GGA, hybrid • Molecular dynamics in NVE, NVT, NPT ensembles • General constraints (MD and geometry optimisation) • Metadynamics ‡ • Free energy functional ‡ • Path integrals ‡ • QM/MM ‡ • Wannier functions • Response properties: TDDFT, NMR, Raman, IR, . • Norm conserving and ultra-soft pseudo potentials ‡ : Will not be considered during this course CPMD: Examples Amorphous Si CPMD: Examples Amorphous Si CPMD: Examples Liquid water: Solvation and transport of
    [Show full text]
  • Pyscf: the Python-Based Simulations of Chemistry Framework
    PySCF: The Python-based Simulations of Chemistry Framework Qiming Sun∗1, Timothy C. Berkelbach2, Nick S. Blunt3,4, George H. Booth5, Sheng Guo1,6, Zhendong Li1, Junzi Liu7, James D. McClain1,6, Elvira R. Sayfutyarova1,6, Sandeep Sharma8, Sebastian Wouters9, and Garnet Kin-Lic Chany1 1Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena CA 91125, USA 2Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA 3Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 4Department of Chemistry, University of California, Berkeley, California 94720, USA 5Department of Physics, King's College London, Strand, London WC2R 2LS, United Kingdom 6Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA 7Institute of Chemistry Chinese Academy of Sciences, Beijing 100190, P. R. China 8Department of Chemistry and Biochemistry, University of Colorado Boulder, Boulder, CO 80302, USA 9Brantsandpatents, Pauline van Pottelsberghelaan 24, 9051 Sint-Denijs-Westrem, Belgium arXiv:1701.08223v2 [physics.chem-ph] 2 Mar 2017 Abstract PySCF is a general-purpose electronic structure platform designed from the ground up to emphasize code simplicity, so as to facilitate new method development and enable flexible ∗[email protected] [email protected] 1 computational workflows. The package provides a wide range of tools to support simulations of finite-size systems, extended systems with periodic boundary conditions, low-dimensional periodic systems, and custom Hamiltonians, using mean-field and post-mean-field methods with standard Gaussian basis functions. To ensure ease of extensibility, PySCF uses the Python language to implement almost all of its features, while computationally critical paths are implemented with heavily optimized C routines.
    [Show full text]