DOCSLIB.ORG

Introduction to Pseudopotentials and Electronic Structure Philip B

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Introduction to Pseudopotentials and Electronic Structure Philip B Introduction to Pseudopotentials and Electronic Structure Philip B. Allen Stony Brook University Nov. 3, 2014 Brookhaven National Laboratory You searched for: TOPIC: (TiO2 110 surface) – Results: 3,551 (from Web of Science Core Collection) Example: rutile (TiO2) Hybrid func,onal studies of the oxygen 6-atom tetragonal unit cell vacancy in TiO2, A. Jano, J. B. Varley, P. Rinke, N. Umezawa, G. Kresse, and C. G. Van de Walle, PRB 81, 085212 (2010) Ti: 4s23d2 (4 electrons, 6 orbitals) O: 2s22p4 (6 electrons, 4 orbitals) 2(TiO2) – 32 valence electrons 28 minimal basis orbitals ~1130 plane waves DFT hybrid 12 empty bands Perdew, Burke, Ernzerhof Heyd, Scuseria, and Ernzerhof 12 filled bands 4 oxygen 2s bands lower in energy Single-particle ideas and philosophies – converging or diverging? 1. Drude, Sommerfeld, Bloch. Classical, Quantum, free to bound in solid. ! 1/ 3 2. Hartree-Fock; Slater ρ ( r ) (average exchange). 3. Interacting many-body theory: Single-particle Green’s function; Landau theory; “GW” approximation 4. Hohenberg-Kohn-Sham density functional theory (DFT) € Single-particle Methods of Solution 1. Local orbitals; LCAO-Hückel 2. Plane waves 3. Augmented plane waves (APW) 4. Othogonalized plane waves (OPW) 5. Green’s function methods (KKR) 6. Pseudopotentials 7. PAW Walter Kohn Pseudopotentials (~1959) were born before DFT (1964) and remain partially independent. DFT “does” ρ(r). What do pseudopotentials do? Bloch’s theorem: can choose eigenstates of a periodic Hamiltonian in the form ! ! ! ! ! ! ! ˆ ! T (R )ψ(r ) ≡ ψ(r + R ) = exp(ik ⋅ R )ψ(r ) (label as ψk ) ˆ ! (eigenvalue of H labeled as εk ) ! ! ! (eigenvalue of translation Tˆ (R ) is exp(ik ⋅ R ) ! ! ! k and k + G give the same translation eigenvalue. € 2 2 1d free electrons ε ! = " k /2m k 1d free electrons with periodicity of k-space (the reciprocal lattice) 1d single-orbital tight-binding band LCAO: ψk=Σexp(ikR)φ(r-R) εk = -2tcos(ka) € volume 51 (1937) Augmented Plane Waves (APW) The wave funcPon is expanded in spherical harmonics and radial soluPons of the wave equaon within the spheres, and in plane waves outside the spheres, joining students: conPnuously at the surface. A single unperturbed funcPon consists of a single plane Leland Allen wave outside the spheres, together with the Don Ellis necessary spherical funcPons within the Art Freeman spheres. ... It is hoped that the method will be Frank Herman useful for comparavely low energy excited George Koster electrons, for which the usual method of Len Mattheiss expansion in plane waves converges too slowly. Dick Watson John Wood Birth of computational electronic structure theory? The pseudopotential emerged ~1959 as a related way of reducing the number of plane waves needed. Richard M. Martin, Electronic Structure: Basic Theory and Practical Methods Cambridge (2004) 3 How many plane waves are needed? ~(Gmax/Gmin) Crude estimate for s-orbitals: λmin ~ 2aB/Z 3 3 (Gmax/Gmin) ~ (Za/2√3aB) ~ 71,400 for Si. With a smooth pseudopotential, λmin ~ aB works. ~ 200 plane waves for Si. But it depends on the kind of pseudopotential. Bloch’s theorem ! “k” is short for k , n , σ (all the quantum numbers of the state). Let there be a “single-particle” Hamiltonian. Then: Need to solve (in a primitive cell): Two€ issues: 1. What is H? 2. How to solve? Most methods of solution expand ψ in a “basis set.” The simplest basis set adapted to Bloch’s theorem is plane-waves. ! ! This is a Bloch wave, since exp( iG ⋅ R ) =1 They are complete, orthogonal, and bias-free. € Dimensionless units “Rydberg atomic units” “Hartree atomic units” In both systems, the unit of length is the Bohr radius: Simplified notation; drop the label “k”. Determine the expansion coefficients cg (secular equation.) This is only one of many matrix versions of the Schrödinger equation ∞ where is a complete, orthonormal set of states. ψ = ∑ α n n { n } n=1 ∞ ∞ ˆ ˆ Hψ = ∑ α n H n = ∑ αn E n = Eψ € € n=1 n=1 ∞ ∞ Left project by m ∑ α n m H n = ∑ α n E m n = Eα m n=1 n=1 € ⎛ 1 H 1 ! 1 H N !⎞⎛ α1 ⎞ ⎛ α1 ⎞ Exact version of ⎜ ⎟⎜ ⎟ ⎜ ⎟ one-electron ⎜ " € # €" !⎟⎜ " ⎟ ⎜ " ⎟ Schroedinger eqn. = E (assuming we know H!) ⎜ N H 1 ! N H N !⎟⎜ α N ⎟ ⎜α N ⎟ ⎜ ⎟⎜ ⎟ ⎜ ⎟ ⎝ " ! ! #⎠⎝ " ⎠ ⎝ " ⎠ € Truncated version of the matrix formulation – variational accuracy Theorem: if Hψ=Eψ and we expand ψ in a (finite, truncated) orthonormal basis, N ψ = ∑ α n n n=1 then the states that are stationary under variation of all coefficients δ * [ ψ H ψ − λ ψ ψ ] = 0 € δαn obey the matrix eigenvalue equation ⎛ ⎞⎛ ⎞ ⎛ ⎞ € 1 H 1 ! 1 H N α1 α1 ⎜ ⎟⎜ ⎟ ⎜ ⎟ ⎜ " # " ⎟⎜ " ⎟ = λ⎜ " ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ N H 1 N H N ⎟ α α ⎝ ! ⎠⎝ n ⎠ ⎝ n ⎠ € ⎛ α1 ⎞ ⎛ α1 ⎞ Basis functions are often not ⎜ ⎟ ⎜ ⎟ orthogonal or normalized. ! ! This creates a “generalized” Hˆ ⎜ ⎟ = ESˆ ⎜ ⎟ Hermitean eigenvalue problem. ⎜α N ⎟ ⎜α N ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ! ⎠ ⎝ ! ⎠ ⎛ 1 H 1 ! 1 H N !⎞ ⎛ 1 1 ! 1 N !⎞ ⎜ ⎟ ⎜ ⎟ " # " ! " # " ! Hˆ = ⎜ € ⎟ Sˆ = ⎜ ⎟ ⎜ N H 1 ! N H N !⎟ ⎜ N 1 ! N N !⎟ ⎜ ⎟ ⎜ ⎟ " ! ! # " ! ! # ⎝ ⎠ ⎝ ⎠ Linear algebra methods solve this with little added difficulty. One approach: € Hˆ ψ = ESˆ ψ and φ€= Sˆ 1/ 2ψ ˆ ˆ ˆ −1/ 2 ˆ ˆ −1/ 2 → H effψ = Eψ where H eff = S HS € More annoyingly, efficient basis functions (or a “pseudopotential”) tend to be energy-dependent. Then it is not a linear problem any more. ⎛ α1 ⎞ ⎛ α1 ⎞ ⎜ ⎟ ⎜ ⎟ ˆ ⎜ ! ⎟ ˆ ⎜ ! ⎟ Solve iteratively? H (E) = ES (E) Lose nice linear algebra theorems. ⎜α N ⎟ ⎜α N ⎟ ⎜ ⎟ ⎜ ⎟ ⎝ ! ⎠ ⎝ ! ⎠ Hˆ (E) ψ = ESˆ (E) ψ € € Hˆ (E) ψ = ESˆ (E) ψ Linearized approximations: LAPW (O. K. Andersen), etc. ˆ ˆ ˆ H (E) ≈ H (E0 ) + (E − E0 )H Eʹ € 0 Sˆ (E) ≈ Sˆ (E ) + (E − E )Sˆ ʹ 0 0 E0 This leads to a generalized linear Hermitean eigenvalue problem: ˆ ˆ ˆ H (E0 ) − E0 ψ = (E − E0 ) 1+ E0S Eʹ − H Eʹ ψ €[ ] [ 0 0 ] and works only in an energy interval around E0. € The soft letter ell (L): “linearized” (LAPW) or “local” (LDA) Janotti, Varley, Rinke, Umezawa, Kresse, and Back to rutile Van de Walle, PRB 81, 085212 (2010) 4 x 3 oxygen 2p levels p(x), p(y), p(z) bonded with 2 x 2 Ti Eg levels d(3z2-r2), d(3x2-r2), d(3y2-r2) Local orbital picture necessary to describe covalent orbital mixing chemistry. covalent bonding/anti-bonding splitting. “Oxygen vacancies in rutile TiO2 were simulated by removing one O atom from a supercell with 72 atoms.” 12 primitive cells ! 336 orbitals in a minimal basis “We used a plane-wave basis set with a cutoff of 400 eV” 2 2 400 eV = kmax = (2π/λmin) corresponds to a minimum wave-length of 0.6 Å 3 ! kmax=5.42 a.u. ! Vk-space (4π/3)kmax = 667 a.u. The super-cell has dimensions V = 2a x 2a x 3c = 748 Å3 = 5,055 a.u. 3 The Brillouin zone volume VBZ = (2π) /V = 0.049 a.u. The matrix dimension is Vk-space/VBZ = 13,600 plane waves Plane waves enlarge the basis set by 40 and the computer time by (40)3 = 64,000 “... integrations over the Brillouin zone were perfomed using a 2 x 2 x 2 mesh of special k points.” Self-consistent potential V(r) depends on charge density ρ(r) 0 (π/a)(1,0,0) Discrete k-mesh of special points used to evaluate ρ(r). Birth of the pseudopotential Phys. Rev. 116, 287 (1959) Conyers Herring ! " real wave - function : ψ(r ) = r ψ ˜ ! " ˜ pseudo wave - function : ψ (r ) = r ψ ψ˜ ≡ ψ − ∑ i i ψ i € i i is the projection operator that orthogonalizes ∑ ψ to the (somewhat arbitrary?) states |i>. i H ψ =€E ψ defines an effective Hamiltonian H˜ ψ ˜ = Eψ˜ . If H = T +V (V = potential), then H˜ = T +V˜ (V˜ = pseudopotential). € DFT is not yet born. People were guessing that V “existed” (Slater showed how to “guess” it quite well). Then V ˜ is a weaker potential, a non-local € operator, that also exists and has pseudo-wavefunctions for its solution. So, brave computational scientist, proceed to guess V ˜ ! € € R. M. Martin, Electronic Structure ... Pseudopotential for a problem where the 1-electron potential is “known”: – real pseudo – Phillips-Kleinman-Antoncik The full self-consistent treatment was not an option until 15 years later. The full apparatus of non-local (and energy-dependent) potentials was over the top. Simpler approaches were encouraged. “SL” = semi-local Phys. Rev. 141, 789–796 (1966) [cited 1336 times] Empirical pseudopotentials “Empirical pseudo-potential method” local version ! 1 1 ! 2 a a 0, , A π (1 11) = ( 2 2 ) = a ! ! fcc structure: b a 1 ,0, 1 B 2π (11 1) = ( 2 2 ) = a ! 1 1 ! 2 c = a , ,0 C = π (111 ) ( 2 2 ) a 5 smallest reciprocal lattice vector (RLV) types ! ! ! ! ! G = 2π (111), G = 2π (200), G = 2π (220), G = 2π (311), G = 2π (222) 3 a 4 a 8 a 11 a 12 a € ! ! ! ! matrix elements k V k + G 3 = V(| G 3 |) ≡ V3 are the same for all 12 vectors € Diamond structure: fcc with a basis of two identical atoms at locations τ1 and τ2.. ! ! ! ! ! ! ! ! ! ! k V k + G = V(G )S(G ) where S(G ) = eiG ⋅τ 1 + eiG ⋅τ 2 /2 € ( ) Note – destructive interference kills structure factor S(G) for certain G’s. ! ! S(G ) S(G ) 0 Three numbers (V3, V8, V11) give a “good” 4 = 12 = band structure for C, Si, Ge, Sn. € € Empirical pseudopotentials are adjusted until theory “agrees” with optical experiment. There are aspects of experiment lying outside “quasi-particle theory:” excitonic effects, for example. Two particle Green’s function, Bethe-Salpeter equation.
Load more

Recommended publications
  • Simplified Hartree-Fock Computations on Second-Row Atoms
    Simplified Hartree-Fock Computations on Second-Row Atoms S. M. Blinder Wolfram Research, Inc. Champaign, IL 61820, USA May 18, 2021 Modern computational quantum chemistry has developed to a large extent beginning with applications of the Hartree-Fock method to atoms and molecules[1]. Density-functional methods have now largely supplanted conventional Hartree-Fock computations in current applications of computational chemistry. Still, from a pedagogical point of view, an un- derstanding of the original methods remains a necessary preliminary. However, since more elaborate Hartree-Fock computations are now no longer needed, it will suffice to consider a computationally simplified application of the method[2]. Accordingly, we will represent the many-electron wavefunction by a single closed-shell Slater determinant and use basis functions which are simple Slater-type orbitals. This is in contrast to the more elaborate Hartree-Fock computations, which might use multi-determinant wavefunctions and double- zeta basis functions. All of the symbolic and numerical operations were carried out using Mathematica. We consider Hartree-Fock computations on the ground states of the second-row atoms He through Ne, Z = 2 to 10, using the simplest set of orthonormalized 1s, 2s and 2p orbital functions: s α3=2 3β5 α + β = p e−αr; = 1 − r e−βr; 1s π 2s π(α2 − αβ + β2) 3 5=2 γ −γr 2pfx;y;zg = p re fsin θ cos φ, sin θ sin φ, cos θg: (1) arXiv:2105.07018v1 [quant-ph] 14 May 2021 π An orbital function multiplied by a spin function α or β is known as a spinorbital, a product of the form ( α φ(x) = (r) : (2) β 1 A simple representation of a many-electron atom is given by a Slater determinant con- structed from N occupied spin-orbitals: φa(1) φb(1) : : : φn(1) 1 φa(2) φb(2) : : : φn(2) Ψ(1;:::;N) = p ; (3) .
    [Show full text]
  • A Self-Interaction-Free Local Hybrid Functional: Accurate Binding
    A self-interaction-free local hybrid functional: accurate binding energies vis-`a-vis accurate ionization potentials from Kohn-Sham eigenvalues Tobias Schmidt,1, ∗ Eli Kraisler,2, ∗ Adi Makmal,2, † Leeor Kronik,2 and Stephan K¨ummel1 1Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany 2Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100,Israel (Dated: September 12, 2018) We present and test a new approximation for the exchange-correlation (xc) energy of Kohn- Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density. It contains one parameter that is not determined ab initio. We investigate whether it is possible to construct a functional that yields accurate binding energies and affords other advantages, specifically Kohn-Sham eigenvalues that reliably reflect ionization potentials. Tests for a set of atoms and small molecules show that within our local-hybrid form accurate binding energies can be achieved by proper optimization of the free parameter in our functional, along with an improvement in dissociation energy curves and in Kohn-Sham eigenvalues. However, the correspondence of the latter to experimental ionization potentials is not yet satisfactory, and if we choose to optimize their prediction, a rather different value of the functional’s parameter is obtained. We put this finding in a larger context by discussing similar observations for other functionals and possible directions for further functional development that our findings suggest.
    [Show full text]
  • Account 0103 - 5053 $6.00+0.00Ac Carbocations on Zeolites
    J. Braz. Chem. Soc., Vol. 22, No. 7, 1197-1205, 2011. Printed in Brazil - ©2011 Sociedade Brasileira de Química Account 0103 - 5053 $6.00+0.00Ac Carbocations on Zeolites. Quo Vadis? Claudio J. A. Mota* and Nilton Rosenbach Jr.# Instituto de Química and INCT de Energia e Ambiente, Universidade Federal do Rio de Janeiro, Cidade Universitária, CT Bloco A, 21941-909 Rio de Janeiro-RJ, Brazil A natureza dos carbocátions adsorvidos na superfície da zeólita é discutida neste artigo, destacando estudos experimentais e teóricos. A adsorção de halogenetos de alquila sobre zeólitas trocadas com metais tem sido usada para estudar o equilíbrio entre alcóxidos, que são espécies covalentes, e carbocátions, que têm natureza iônica. Cálculos teóricos indicam que os carbocátions são mínimos (intermediários) na superfície de energia potencial e estabilizados por ligações hidrogênio com os átomos de oxigênio da estrutura zeolítica. Os resultados indicam que zeólitas se comportam como solventes sólidos, estabilizando a formação de espécies iônicas. The nature of carbocations on the zeolite surface is discussed in this account, highlighting experimental and theoretical studies. The adsorption of alkylhalides over metal-exchanged zeolites has been used to study the equilibrium between covalent alkoxides and ionic carbocations. Theoretical calculations indicated that the carbocations are minima (intermediates) on the potential energy surface and stabilized by hydrogen bonds with the framework oxygen atoms. The results indicate that zeolites behave like solid solvents, stabilizing the formation of ionic species. Keywords: zeolites, carbocations, alkoxides, catalysis, acidity 1. Introduction Thus, NH4Y can be prepared by exchange of the NaY with ammonium salt solutions. An acidic material can Zeolites are crystalline aluminosilicates with pores of be obtained by calcination of the ammonium-exchanged molecular dimensions.
    [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]
  • ACF NATIONALS 2018 ROUND 1 PRELIMS 1 Packet by OXFORD + DUKE
    ACF NATIONALS 2018 ROUND 1 PRELIMS 1 packet by OXFORD + DUKE authors Oxford: Daoud Jackson, George Charlson, Jacob Robertson, Freddy Potts, Chris Stern Duke: Ryan Humphrey, Gabe Guedes, Lucian Li, Annabelle Yang ACF Nationals 2018 | Packet: Oxford + Duke | Page 1 Editors: Jordan Brownstein, Andrew Hart, Stephen Liu, Aaron Rosenberg, Andrew Wang, Ryan Westbrook Tossups 1. ARCIMBOLDO (ar-chim-BOL-do) and the auto-rickshaw platform process data from this technique. The sphere-of- influence algorithm for this technique helps differentiate solvents and the target molecule, and like the free lunch algorithm, is found in the SHELXE (“shell X-E”) package. Data from this technique can be more easily processed by repeating data collection after either using a laser to induce radiation damage in the sample, or soaking the sample in a heavy metal compound. This is the most prominent technique for which Rigaku produces instruments, including the MiniFlex. MAD, SAD, and isomorphous replacement are methods used in this technique to solve the phase problem. The systematic absences from this technique are used to help find the space group of the analyte, which must be determined in order to solve the final structure. For 10 points, name this technique that uses Bragg’s law to determine the structure of a crystal. ANSWER: x-ray diffraction [accept x-ray crystallography or XRD] 2. Five years after this man’s death, his diaries were collected by the manager of the Pomeranian Mortgage Bank, including correspondences produced while he lived with his former lieutenant Wilhelm Junker. This man’s peaceful meeting with Omukama Kabalega (oh-moo-KAH-mah kah-bah-LEH-gah) is commemorated by a cone-shaped structure at the Mparo Royal Tombs.
    [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]
  • Core Electron Binding Energies in Heavy Atoms
    l~m;~, Vol. 21, No. 2, August 1983, pp. 103-110. © Printed in India. Core electron binding energies in heavy atoms M P DAS Department of Physics, Sambalpur University,Jyoti Vihar, Sambalpur 768 017, India MS received 8 October 1982; revised 26 April 1983 Abstract. Inner shell binding of electrons in heavy atoms is studied through the relativistic density functional theory in which many electron interactions are treated in a local density approximation. By using this theory and the A scF procedure binding energies of several core electrons of mercury atom are calculatedin the frozen and relaxed configurations. The results are compared with those carried out by the non-local Dirac-Fock Scheme. K-shellbinding energies of several closed shell atoms are calculated by using the Kolm-Shamand the relativisticexchange potentials. The results are discussed and the discrepancies in our local densityresults, when compared with experimental values, may be attributed to the non-localityand to the many-body effects. Keywords. Atomic structure; binding ener~; relaxed orbitals; density functional; Breit interaction. 1. Introduction Calculations on the electronic structure of heavy atoms based on the relativistic theory, such as multi-oonfigurational Dirac-Fock (MeDF) method (Desolaux 1980) are in reasonably good agreement with experimental findings. This theory employs one- electron Dirao Hamiltonian that contains the kinematics of the electrons and their interactions with the classical nuclear field. The electron-electron interaction is then added in two parts. The first part is treated in ~e Hartree-Fock sense by applying variational principle to a many-electron wavefunction as an antisymmctric product of one-electron orbitals.
    [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]
  • New Method for Measuring Angle-Resolved Phases in Photoemission
    PHYSICAL REVIEW X 10, 031070 (2020) New Method for Measuring Angle-Resolved Phases in Photoemission Daehyun You ,1 Kiyoshi Ueda,1,* Elena V. Gryzlova ,2 Alexei N. Grum-Grzhimailo ,2 Maria M. Popova ,2,3 Ekaterina I. Staroselskaya,3 Oyunbileg Tugs,4 Yuki Orimo ,4 Takeshi Sato ,4,5,6 Kenichi L. Ishikawa ,4,5,6 Paolo Antonio Carpeggiani ,7 Tamás Csizmadia,8 Miklós Füle,8 Giuseppe Sansone ,9 Praveen Kumar Maroju,9 Alessandro D’Elia ,10,11 Tommaso Mazza,12 Michael Meyer ,12 Carlo Callegari ,13 Michele Di Fraia,13 Oksana Plekan ,13 Robert Richter,13 Luca Giannessi,13,14 Enrico Allaria,13 Giovanni De Ninno,13,15 Mauro Trovò,13 ‡ Laura Badano,13 Bruno Diviacco ,13 Giulio Gaio,13 David Gauthier,13, Najmeh Mirian,13 Giuseppe Penco,13 † Primož Rebernik Ribič ,13,15 Simone Spampinati,13 Carlo Spezzani,13 and Kevin C. Prince 13,16, 1Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, Sendai 980-8577, Japan 2Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Moscow 119991, Russia 3Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia 4Department of Nuclear Engineering and Management, Graduate School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan 5Photon Science Center, Graduate School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan 6Research Institute for Photon Science and Laser Technology, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 7Institut für Photonik, Technische
    [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]