VASP: Basics (DFT, PW, PAW, … )

Total Page:16

File Type:pdf, Size:1020Kb

VASP: Basics (DFT, PW, PAW, … ) VASP: Basics (DFT, PW, PAW, … ) University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria Outline ● Density functional theory ● Translational invariance and periodic boundary conditions ● Plane wave basis set ● The Projector-Augmented-Wave method ● Electronic minimization The Many-Body Schrödinger equation Hˆ (r1,...,rN )=E (r1,...,rN ) 1 1 ∆i + V (ri)+ (r1,...,rN )=E (r1,...,rN ) 0−2 ri rj 1 i i i=j X X X6 | − | @ A For instance, many-body WF storage demands are prohibitive: N (r1,...,rN ) (#grid points) 5 electrons on a 10×10×10 grid ~ 10 PetaBytes ! A solution: map onto “one-electron” theory: (r ,...,r ) (r), (r),..., (r) 1 N ! { 1 2 N } Hohenberg-Kohn-Sham DFT Map onto “one-electron” theory: N (r ,...,r ) (r), (r),..., (r) (r ,...,r )= (r ) 1 N ! { 1 2 N } 1 N i i i Y Total energy is a functional of the density: E[⇢]=T [ [⇢] ]+E [⇢]+E [⇢] + E [⇢]+U[Z] s { i } H xc Z The density is computed using the one-electron orbitals: N ⇢(r)= (r) 2 | i | i X The one-electron orbitals are the solutions of the Kohn-Sham equation: 1 ∆ + V (r)+V [⇢](r)+V [⇢](r) (r)=✏ (r) −2 Z H xc i i i ⇣ ⌘ BUT: Exc[⇢]=??? Vxc[⇢](r)=??? Exchange-Correlation Exc[⇢]=??? Vxc[⇢](r)=??? • Exchange-Correlation functionals are modeled on the uniform-electron-gas (UEG): The correlation energy (and potential) has been calculated by means of Monte- Carlo methods for a wide range of densities, and has been parametrized to yield a density functional. • LDA: we simply pretend that an inhomogeneous electronic density locally behaves like a homogeneous electron gas. • Many, many, many different functionals available: LDA, GGA, meta-GGA, van-der-Waals functionals, etc etc An N-electron system: N = O(1023) Hohenberg-Kohn-Sham DFT takes us a long way: (r ,...,r ) (r), (r),..., (r) 1 N ! { 1 2 N } (#grid points)N N (#grid points) ⇥ Nice for atoms and molecules, but in a realistic piece of solid state material N= O(1023)! Translational invariance: Periodic Boundary Conditions Translational invariance implies: ikR nk(r + R)= nk(r)e and ikr nk(r)=unk(r)e unk(r + R)=unk(r) All states are labeled by Bloch vector k and the band index n: • The Bloch vector k is usually constrained to lie within the first Brillouin zone of the reciprocal space lattice. • The band indeX n is of the order if the number of electrons per unit cell. Reciprocal space & the first Brillouin zone z L 4pi/L z z b b b1 2 1 b2 a2 y y a1 x x a3 y x b3 b3 A B C 2⇡ 2⇡ 2⇡ b = a a b = a a b = a a 1 ⌦ 2 ⇥ 3 2 ⌦ 3 ⇥ 1 3 ⌦ 1 ⇥ 2 ⌦ = a a a a b =2⇡ 1 · 2 ⇥ 3 i · j ij Sampling the 1st BZ The evaluation of many key quantities involves an integral over the 1st BZ. For instance the charge density: 1 ⇢(r)= f (r) 2dk ⌦ nk| nk | BZ n BZ X Z We eXploit the fact that the orbitals at Bloch vectors k that are close together are almost identical and approXimate the integral over the 1st BZ by a weighted sum over a discrete set of k-points: ⇢(r)= w f (r) 2dk, k nk| nk | Xnk 23 Fazit: the intractable task of determining with N=10 (r1,...,rN ) , has been st reduced to calculating at a discrete set of nk(r) k-points in the 1 BZ, for a number of bands that is of the order if the number of electrons in the unit cell. Idea: equally spaced mesh in the 1st Brillouin zone Example: a quadratic 2D lattice r • q1=q2=4, i.e., 16 points in total b • Only 3 symmetry inequivalent points: 44 k1 1 1 1 k2 4 k1 =( , ) !1 = k ⇥ 8 8 ) 4 3 r 3 3 1 b 4 k =( , ) ! = 0 ⇥ 2 8 8 ) 2 4 3 1 1 IBZ 8 k =( , ) ! = ⇥ 3 8 8 ) 3 2 BZ 1 1 1 1 F (k)dk F (k1)+ F (k2)+ F (k3) ⌦BZ ) 4 4 2 BZZ Algorithm: • Calculate equally spaced mesh. • Shift the mesh if desired. • Apply all symmetry operations of the Bravais lattice to all k-points. • Extract the irreducible k-points (IBZ). • Calculate the proper weighting. Common meshes: Two different choices for the center of the mesh. • Centered on Γ • Centered around Γ (can break the symmetry!) before after shifted to symmetrization r r • In certain cell geometries (e.g. heXagonal cells) even meshes break the symmetry. • Symmetrization results in non-uniform distributions of k-points. • Γ-point centered meshes preserve the symmetry. AB C The total energy E[⇢, R,Z ]=T [ [⇢] ]+E [⇢, R,Z ]+E [⇢]+U( R,Z ) { } s { nk } H { } xc { } The kinetic energy 1 T [ [⇢] ]= w f ∆ s { nk } k nkh nk| − 2 | nki Xnk The Hartree energy 1 ⇢eZ (r)⇢eZ (r0) E [⇢, R,Z ]= dr0dr H { } 2 r r ZZ | − 0| where 2 ⇢ (r)=⇢(r)+ Z δ(r R ) ⇢(r)= wkfnk nk(r) dk eZ i − i | | i nk X X The Kohn-Sham equations 1 ∆ + V [⇢ ](r)+V [⇢](r) (r)=✏ (r) −2 H eZ xc nk nk nk ⇣ ⌘ The Hartree potential ⇢eZ (r0) V [⇢ ](r)= dr0 H eZ r r Z | − 0| A plane wave basis set ikr nk(r)=unk(r)e unk(r + R)=unk(r) All cell-periodic functions are eXpanded in plane waves (Fourier analysis): 1 iGr 1 i(G+k)r unk(r)= CGnke (r)= C e ⌦1/2 nk 1/2 Gnk G ⌦ X XG iGr iGr ⇢(r)= ⇢Ge V (r)= VGe XG XG The basis set includes all plane waves for which 1 G + k 2 <E 2| | cuto↵ Transformation by means of FFT between “real” space and “reciprocal” space: iGr FFT 1 iGr C = C e C = C e− rnk Gnk ! Gnk N rnk FFT r XG X Why use plane waves? • Historical reason: Many elements eXhibit a band-structure that can be interpreted in a free electron picture (metallic s and p elements). Pseudopotential theory was initially developed to cope with these elements (pseudopotential perturbation theory). • Practical reason: The total energy eXpressions and the Hamiltonian are easy to implement. • Computational reason: The action of the Hamiltonian on the orbitals can be efficiently evaluated using FFTs. FFT real space reciprocal space G cut τ 2 b2 01 23N−1 0 5 −4012345−3 −2 −1 τ 1 b1 −N/2+1 N/2 x = n / N τ g = n 2 π / τ 1 1 1 1 1 1 iGr FFT 1 iGr C = C e C = C e− rnk Gnk ! Gnk N rnk FFT r XG X FFT real space reciprocal space G cut τ 2 b2 01 23N−1 0 5 −4012345−3 −2 −1 τ 1 b1 −N/2+1 N/2 x = n / N τ g = n 2 π / τ 1 1 1 1 1 1 iGr FFT 1 iGr C = C e C = C e− rnk Gnk ! Gnk N rnk FFT r XG X iGr FFT 1 iGr C = C e C = C e− rnk Gnk ! Gnk N rnk FFT r XG X FFT real space reciprocal space G cut τ 2 b2 01 23N−1 0 5 −4012345−3 −2 −1 τ 1 b1 −N/2+1 N/2 x = n / N τ g = n 2 π / τ 1 1 1 1 1 1 iGr FFT 1 iGr C = C e C = C e− rnk Gnk ! Gnk N rnk FFT r XG X The charge density Gcut FFT ψ ψ ψr G 2 G r cut FFT ρ G ρ r The action of the Hamiltonian The action H | nki 1 ∆ + V (r) (r) −2 nk ✓ ◆ Using the convention 1 r G + k = ei(G+k)r G + k = C h | i ⌦1/2 ! h | nki Gnk • Kinetic energy: 1 1 2 G + k ∆ = G + k C NNPLW h | − 2 | nki 2| | Gnk • Local potential: V = VH[⇢]+Vxc[⇢]+Vext • Exchange-correlation: easily obtained in real space Vxc,r = Vxc[⇢r] • FFT to reciprocal space V V { xc,r} ! { xc,G} 4⇡ • Hartree potential: solve Poisson eq. in reciprocal space VH,G = ⇢G G 2 • Add all contributions VG = VH,G + Vxc,G + Vext,G | | • FFT to real space V V { G} ! { r} The action 1 iGr G + k V = V C e− N log N h | | nki N r rnk FFT FFT FFT r X The action of the local potential 2Gcut 4π e 2 FFT G 2 V V ρ G G r ψr 3G Gcut cut add FFT R R (residual vector) ψG G r The Projector-Augmented-Wave method The number of plane waves needed to describe • tightly bound (spatially strongly localized) states, • and rapid oscillations (nodal features) of the orbitals near the nucleus exceeds any practical limit, except maybe for Li and H. The common solution: • Introduce the frozen core approXimation: Core electrons are pre-calculated in an atomic environment and kept frozen in the course of the remaining calculations. • Use of pseudopotentials instead of eXact potentials: • Norm-conserving pseudopotentials • Ultra-soft pseudopotentials • The Projector-Augmented-Wave (PAW) method [P. E. Blöchl, Phys. Rev. B 50, 17953 (1994)] Pseudopotentials: the general idea ~ Ψ Ψ rc ~ V V Pseudopotentials: cont. 2.5 2.0 1.5 s : E= -0.576 R c=1.9 1.0 wave-function 0.5 Al Al p : E= -0.205 R c=1.9 0.0 0 1 2 3 4 R (a.u.) Exact potential (interstitial) Pseudopotential Al effective Al atom PAW Al atom 3p 2p 3p 3s 1s 3s 2p 2p and 1s are nodeless nodal structure retained 2s 12 The PAW orbitals = + ( φ φ ) p | ni | ni | ii| ii h i| ni i X e e e e φ p φ p ˜ | ni | ni− | iih i| ni | ni | iih i| ni i i X X+ φ p e e e e e e e | ieih i| ni i X e e PW PS-LCAO AE-LCAO = - + AE pseudopseudo-onsite AE-onsite The PAW orbitals (cont.) = + ( φ φ ) p | ni | ni | ii| ii h i| ni i X e e e • n is the pseudo (PS) orbital, eXpressed in a plane wave basis sete | i • φei , , and are atomφ˜ p˜i -centered localized functions | i | ii | i • The all-electron partial waves are solutions to the radial scalar relativistic non-spinpolarized Schrödinger equation: 1 ( ∆ + v ) φ = ✏ φ −2 e↵ | ii i| ii • A pseudization procedure yields: φ φ v v p φ = δ | ii!| ii e↵ ! e↵ h i| ji ij e e e e The PAW orbitals (cont.) • The pseudo partial waves obey: 1 ∆ + v + p D p φ = ✏ 1+ p Q p φ − 2 e↵ | ii ijh j| | ki k | ii ijh j| | ki ij ij ⇣ X ⌘ ⇣ X ⌘ e e e e e e e • with the so-called PAW strength parameters and augmentation charges: 1 1 D = φ ∆ + v φ φ ∆ + v φ ij h i| − 2 e↵ | jih i| − 2 e↵ | ji e e e Q = φ φ φ φ ij h i| jih i| ji e e The all-electron and pseudo eigenvalue spectrum is identical! All-electron scattering properties are reproduced over a wide energy-range.
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]
  • Basis Sets Running a Calculation • in Performing a Computational
    Session 2: Basis Sets • Two of the major methods (ab initio and DFT) require some understanding of basis sets and basis functions • This session describes the essentials of basis sets: – What they are – How they are constructed – How they are used – Significance in choice 1 Running a Calculation • In performing a computational chemistry calculation, the chemist has to make several decisions of input to the code: – The molecular geometry (and spin state) – The basis set used to determine the wavefunction – The properties to be calculated – The type(s) of calculations and any accompanying assumptions 2 Running a calculation •For ab initio or DFT calculations, many programs require a basis set choice to be made – The basis set is an approx- imate representation of the atomic orbitals (AOs) – The program then calculates molecular orbitals (MOs) using the Linear Combin- ation of Atomic Orbitals (LCAO) approximation 3 Computational Chemistry Map Chemist Decides: Computer calculates: Starting Molecular AOs determine the Geometry AOs O wavefunction (ψ) Basis Set (with H H ab initio and DFT) Type of Calculation LCAO (Method and Assumptions) O Properties to HH be Calculated MOs 4 Critical Choices • Choice of the method (and basis set) used is critical • Which method? • Molecular Mechanics, Ab initio, Semiempirical, or DFT • Which approximation? • MM2, MM3, HF, AM1, PM3, or B3LYP, etc. • Which basis set (if applicable)? • Minimal basis set • Split-valence • Polarized, Diffuse, High Angular Momentum, ...... 5 Why is Basis Set Choice Critical? • The
    [Show full text]
  • An Introduction to Hartree-Fock Molecular Orbital Theory
    An Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2000 1 Introduction Hartree-Fock theory is fundamental to much of electronic structure theory. It is the basis of molecular orbital (MO) theory, which posits that each electron's motion can be described by a single-particle function (orbital) which does not depend explicitly on the instantaneous motions of the other electrons. Many of you have probably learned about (and maybe even solved prob- lems with) HucÄ kel MO theory, which takes Hartree-Fock MO theory as an implicit foundation and throws away most of the terms to make it tractable for simple calculations. The ubiquity of orbital concepts in chemistry is a testimony to the predictive power and intuitive appeal of Hartree-Fock MO theory. However, it is important to remember that these orbitals are mathematical constructs which only approximate reality. Only for the hydrogen atom (or other one-electron systems, like He+) are orbitals exact eigenfunctions of the full electronic Hamiltonian. As long as we are content to consider molecules near their equilibrium geometry, Hartree-Fock theory often provides a good starting point for more elaborate theoretical methods which are better approximations to the elec- tronic SchrÄodinger equation (e.g., many-body perturbation theory, single-reference con¯guration interaction). So...how do we calculate molecular orbitals using Hartree-Fock theory? That is the subject of these notes; we will explain Hartree-Fock theory at an introductory level. 2 What Problem Are We Solving? It is always important to remember the context of a theory.
    [Show full text]
  • User Manual for the Uppsala Quantum Chemistry Package UQUANTCHEM V.35
    User manual for the Uppsala Quantum Chemistry package UQUANTCHEM V.35 by Petros Souvatzis Uppsala 2016 Contents 1 Introduction 6 2 Compiling the code 7 3 What Can be done with UQUANTCHEM 9 3.1 Hartree Fock Calculations . 9 3.2 Configurational Interaction Calculations . 9 3.3 M¨ollerPlesset Calculations (MP2) . 9 3.4 Density Functional Theory Calculations (DFT)) . 9 3.5 Time Dependent Density Functional Theory Calculations (TDDFT)) . 10 3.6 Quantum Montecarlo Calculations . 10 3.7 Born Oppenheimer Molecular Dynamics . 10 4 Setting up a UQANTCHEM calculation 12 4.1 The input files . 12 4.1.1 The INPUTFILE-file . 12 4.1.2 The BASISFILE-file . 13 4.1.3 The BASISFILEAUX-file . 14 4.1.4 The DENSMATSTARTGUESS.dat-file . 14 4.1.5 The MOLDYNRESTART.dat-file . 14 4.1.6 The INITVELO.dat-file . 15 4.1.7 Running Uquantchem . 15 4.2 Input parameters . 15 4.2.1 CORRLEVEL ................................. 15 4.2.2 ADEF ..................................... 15 4.2.3 DOTDFT ................................... 16 4.2.4 NSCCORR ................................... 16 4.2.5 SCERR .................................... 16 4.2.6 MIXTDDFT .................................. 16 4.2.7 EPROFILE .................................. 16 4.2.8 DOABSSPECTRUM ............................... 17 4.2.9 NEPERIOD .................................. 17 4.2.10 EFIELDMAX ................................. 17 4.2.11 EDIR ..................................... 17 4.2.12 FIELDDIR .................................. 18 4.2.13 OMEGA .................................... 18 2 CONTENTS 3 4.2.14 NCHEBGAUSS ................................. 18 4.2.15 RIAPPROX .................................. 18 4.2.16 LIMPRECALC (Only openmp-version) . 19 4.2.17 DIAGDG ................................... 19 4.2.18 NLEBEDEV .................................. 19 4.2.19 MOLDYN ................................... 19 4.2.20 DAMPING ................................... 19 4.2.21 XLBOMD ..................................
    [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]
  • 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]
  • 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]
  • Car-Parrinello Molecular Dynamics
    CPMD Car-Parrinello Molecular Dynamics An ab initio Electronic Structure and Molecular Dynamics Program The CPMD consortium WWW: http://www.cpmd.org/ Mailing list: [email protected] E-mail: [email protected] January 29, 2019 Send comments and bug reports to [email protected] This manual is for CPMD version 4.3.0 CPMD 4.3.0 January 29, 2019 Contents I Overview3 1 About this manual3 2 Citation 3 3 Important constants and conversion factors3 4 Recommendations for further reading4 5 History 5 5.1 CPMD Version 1.....................................5 5.2 CPMD Version 2.....................................5 5.2.1 Version 2.0....................................5 5.2.2 Version 2.5....................................5 5.3 CPMD Version 3.....................................5 5.3.1 Version 3.0....................................5 5.3.2 Version 3.1....................................5 5.3.3 Version 3.2....................................5 5.3.4 Version 3.3....................................6 5.3.5 Version 3.4....................................6 5.3.6 Version 3.5....................................6 5.3.7 Version 3.6....................................6 5.3.8 Version 3.7....................................6 5.3.9 Version 3.8....................................6 5.3.10 Version 3.9....................................6 5.3.11 Version 3.10....................................7 5.3.12 Version 3.11....................................7 5.3.13 Version 3.12....................................7 5.3.14 Version 3.13....................................7 5.3.15 Version 3.14....................................7 5.3.16 Version 3.15....................................8 5.3.17 Version 3.17....................................8 5.4 CPMD Version 4.....................................8 5.4.1 Version 4.0....................................8 5.4.2 Version 4.1....................................8 5.4.3 Version 4.3....................................9 6 Installation 10 7 Running CPMD 11 8 Files 12 II Reference Manual 15 9 Input File Reference 15 9.1 Basic rules........................................
    [Show full text]
  • Qbox (Qb@Ll Branch)
    Qbox (qb@ll branch) Summary Version 1.2 Purpose of Benchmark Qbox is a first-principles molecular dynamics code used to compute the properties of materials directly from the underlying physics equations. Density Functional Theory is used to provide a quantum-mechanical description of chemically-active electrons and nonlocal norm-conserving pseudopotentials are used to represent ions and core electrons. The method has O(N3) computational complexity, where N is the total number of valence electrons in the system. Characteristics of Benchmark The computational workload of a Qbox run is split between parallel dense linear algebra, carried out by the ScaLAPACK library, and a custom 3D Fast Fourier Transform. The primary data structure, the electronic wavefunction, is distributed across a 2D MPI process grid. The shape of the process grid is defined by the nrowmax variable set in the Qbox input file, which sets the number of process rows. Parallel performance can vary significantly with process grid shape, although past experience has shown that highly rectangular process grids (e.g. 2048 x 64) give the best performance. The majority of the run time is spent in matrix multiplication, Gram-Schmidt orthogonalization and parallel FFTs. Several ScaLAPACK routines do not scale well beyond 100k+ MPI tasks, so threading can be used to run on more cores. The parallel FFT and a few other loops are threaded in Qbox with OpenMP, but most of the threading speedup comes from the use of threaded linear algebra kernels (blas, lapack, ESSL). Qbox requires a minimum of 1GB of main memory per MPI task. Mechanics of Building Benchmark To compile Qbox: 1.
    [Show full text]
  • Empirical Pseudopotential Method
    Computational Electronics Empirical Pseudopotential Method Prepared by: Dragica Vasileska Associate Professor Arizona State University The Empirical Pseudopotential Method The concept of pseudopotentials was introduced by Fermi [1] to study high-lying atomic states. Afterwards, Hellman proposed that pseudopotentials be used for calculating the energy levels of the alkali metals [2]. The wide spread usage of pseudopotentials did not occur until the late 1950s, when activity in the area of condensed matter physics began to accelerate. The main advantage of using pseudopotentials is that only valence electrons have to be considered. The core electrons are treated as if they are frozen in an atomic-like configuration. As a result, the valence electrons are thought to move in a weak one-electron potential. The pseudopotential method is based on the orthogonalized plane wave (OPW) method due to Herring. In this method, the crystal wavefuntion ψ k is constructed to be orthogonal to the core states. This is accomplished by expanding ψ k as a smooth part of symmetrized combinations of Bloch functions ϕk , augmented with a linear combination of core states. This is expressed as ψ = ϕ + Φ k k ∑bk,t k,t , (1) t where bk,t are orthogonalization coefficients and Φ k,t are core wave functions. For Si-14, the summation over t in Eq. (1) is a sum over the core states 1s2 2s2 2p6. Since the crystal wave function is constructed to be orthogonal to the core wave functions, the orthogonalization coefficients can be calculated, thus yielding the final expression ψ = ϕ − Φ ϕ Φ k k ∑ k,t k k,t .
    [Show full text]
  • Powerpoint Presentation
    Session 2 Basis Sets CCCE 2008 1 Session 2: Basis Sets • Two of the major methods (ab initio and DFT) require some understanding of basis sets and basis functions • This session describes the essentials of basis sets: – What they are – How they are constructed – How they are used – Significance in choice CCCE 2008 2 Running a Calculation • In performing ab initio and DFT computa- tional chemistry calculations, the chemist has to make several decisions of input to the code: – The molecular geometry (and spin state) – The basis set used to determine the wavefunction – The properties to be calculated – The type(s) of calculations and any accompanying assumptions CCCE 2008 3 Running a calculation • For ab initio or DFT calculations, many programs require a basis set choice to be made – The basis set is an approx- imate representation of the atomic orbitals (AOs) – The program then calculates molecular orbitals (MOs) using the Linear Combin- ation of Atomic Orbitals (LCAO) approximation CCCE 2008 4 Computational Chemistry Map Chemist Decides: Computer calculates: Starting Molecular AOs determine the Geometry AOs O wavefunction (ψ) Basis Set (with H H ab initio and DFT) Type of Calculation LCAO (Method and Assumptions) O Properties to HH be Calculated MOs CCCE 2008 5 Critical Choices • Choice of the method (and basis set) used is critical • Which method? • Molecular Mechanics, Ab initio, Semiempirical, or DFT • Which approximation? • MM2, MM3, HF, AM1, PM3, or B3LYP, etc. • Which basis set (if applicable)? • Minimal basis set • Split-valence
    [Show full text]