Pseudopotentials: design, testing, typical errors

Size: px
Start display at page:

Download "Pseudopotentials: design, testing, typical errors"

Transcription

1 Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015

2 Parameter free calculations. Ab initio theory Often excellent agreement with experiment Lukanov, Garrity, PRB (2012)

3 Parameter free calculations. Ab initio theory Often excellent agreement with experiment Experimental STM Sr on Ge (001) surface Lukanov, Garrity, PRB (2012)

4 Parameter free calculations. Ab initio theory Often excellent agreement with experiment Experimental STM Sr on Ge (001) surface DFT Low Energy Structure Lukanov, Garrity, PRB (2012)

5 Parameter free calculations. Ab initio theory Often excellent agreement with experiment Experimental STM Sr on Ge (001) surface Simulated STM DFT Low Energy Structure Lukanov, Garrity, PRB (2012)

6 Parameter free calculations. Ab initio theory Often excellent agreement with experiment Experimental STM Sr on Ge (001) surface Simulated STM DFT Low Energy Structure Lukanov, Garrity, PRB (2012)

7 Ab initio theory - questions How accurate are our equations / physical approximations? Validation with experiment Lots of work on this popular topic to study NOT the topic of this work.

8 Ab initio theory - questions How accurate are our equations / physical approximations? Validation with experiment Lots of work on this popular topic to study NOT the topic of this work. How precisely are we solving equations? Numerical verification Typical numerical errors? Do our codes agree? Who is correct? Often ignored.

9 What is Si lattice constant in GGA? It depends on time, apparently: And Si is the easiest element to treat!! Historical published Si lattice constants with PBE. Graphic Kurt Lejaeghere et al

10 Outline Density Functional Theory Background Solving Kohn-Sham eqns. All-electron methods Pseudopotentials Testing different basis sets GBRV pseudopotential library/tests Designing pseudopotentials Tradeoffs Practical example for lab

11 Electronic Structure Calculations We want to solve Schrodinger s equation: (For electrons)

12 Electronic Structure Calculations We want to solve Schrodinger s equation: (For electrons) Energy Electron coordinates Many-body Hamiltonian Many-body Wavefunction

13 Electronic Structure Calculations We want to solve Schrodinger s equation: (For electrons)

14 Electronic Structure Calculations We want to solve Schrodinger s equation: (For electrons) Kinetic Electron-ion Electron-electron

15 Electronic Structure Calculations We want to solve Schrodinger s equation: (For electrons) Contains essentially all of chemistry

16 Electronic Structure Calculations

17 Electronic Structure Calculations

18 Electronic Structure Calculations

19 Density Functional Theory

20 Kohn-Sham Equations

21 Kohn-Sham Equations Single particle eigenvalues, wavefunctions External Non-interacting kinetic Classical electrostatic (Hartree) Everything else (exchange-correlation)

22 Kohn-Sham Equations Given Atomic Positions Solve for n(r), Energy Forces

23 Exchange-correlation Functional The approximation in DFT Otherwise formally exact Many approximations exist, varying accuracy

24 Exchange-correlation Functional The approximation in DFT Otherwise formally exact Many approximations exist Local Density Approximation (LDA) Treat E xc like you have a locally uniform electron gas

25 Exchange-correlation Functional

26 Accuracy (LDA / GGA) Structural Properties: Lattice constants ~ 1% Bulk modulus ~10% Silicon LDA Energy vs. Vol Yin Cohen PRB (1982)

27 Accuracy (LDA / GGA) Structural Properties: Lattice constants ~ 1% Bulk modulus ~10% Atomization energy ~30 kcal/mol = 1.2 ev/atom LDA ~8 kcal/mol = 0.35 ev/atom GGA Solid formation energy ~0.25 ev/atom GGA Less when materials similar

28 Accuracy (LDA / GGA) Structural Properties: Lattice constants ~ 1% Bulk modulus ~10% Atomization energy ~30 kcal/mol = 1.2 ev/atom LDA ~8 kcal/mol = 0.35 ev/atom GGA Solid formation energy ~0.25 ev/atom GGA Electronic structure Metals good Insulators Qualitatively good, gaps off ~50% Silicon LDA Band Structure Yin Cohen PRB (1982) Gap is 50% too small

29 DFT Use DFT is the workhorse electronic structure method Used extensively in physics, materials science, chemistry, biology Useful combination of accuracy / speed (scales like N 3 ) Used as atomistic level of multi-scale calculations More accurate methods start from DFT GW, DMFT, QMC

30 How to solve Kohn-Sham? N el coupled differential eqns. Should just be a numerical challenge Choose a basis, Transform to linear algebra problem

31 How to solve Kohn-Sham? N el coupled differential eqns. Should just be a numerical challenge Choose a basis, Transform to linear algebra problem Can we just use plane-waves? Complete, periodic, converges systematically w/ G cut

32 Problem with naive plane-waves Strongly different energy/length scales Core vs valence energies 4 orders of magnitude

33 Problem with naive plane-waves

34 Problem with naive plane-waves

35 All-electron calculations Solve this issue by separating space Core regions vs. interstitial Use different basis sets in both regions Atomic-like Plane-waves Issue need to match boundaries, which requires energies

36 LAPW solution Linearized Augmented Plane waves (WIEN2k, etc) Advantages: Accurate if done carefully Disadvantages: Have to set E to expand around Have to set R MT Expensive, hard for user

37 Pseudopotentials

38 Pseudopotentials Zero outside r c

39 Pseudopotentials

40 Steps for generating PSP s

41 Steps for generating PSP s

42 Steps for generating PSP s

43 Steps for generating PSP s Black magic well, hopefully not anymore

44 A clarification Don t think of psp as statistical fit of a model Think of as improvable basis set.

45 A clarification Don t think of psp as statistical fit of a model Think of as improvable basis set. PSP s should: Be insensitive to minor parameter variations Reproduce many properties at once Work in many chemical environments (transferable) Be computationally inexpensive

46 Some nice properties If we require the following: PS and AE WF's match outside r c PS and AE eigenvalues match PS WF is normalized Hamann, et. al. Phys. Rev. Lett. 43, (1979)

47 Some nice properties If we require the following: Then: PS and AE WF's match outside r c PS and AE eigenvalues match PS WF is normalized Get good scattering properties: Matching even outside r c to second order in ε-ε nl Reproduces atomic system by construction Hopefully transferable to solid state Hamann, et. al. Phys. Rev. Lett. 43, (1979)

48 Scattering properties Modern PSP have multiple projectors: Force matching log derivatives at additional energies

49 PSP pros and cons Advantages: Computationally inexpensive Easy to use Disadvantages: Frozen-core approximation Hard to design/test

Pseudopotentials: design, testing, typical errors

Pseudopotentials: design, testing, typical errors Pseudopotentials: design, testing, typical errors Kevin F. Garrity National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.

More information

Pseudopotentials for hybrid density functionals and SCAN

Pseudopotentials for hybrid density functionals and SCAN Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?

More information

Pseudopotentials: design, testing, typical errors

Pseudopotentials: design, testing, typical errors Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 3 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Testing Basis Sets

More information

All electron optimized effective potential method for solids

All electron optimized effective potential method for solids All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22

More information

Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler

Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler http://www.fhi-berlin.mpg.de/th/th.html I. From the many-particle problem to the Kohn-Sham functional II. From

More information

Density Functional Theory. Martin Lüders Daresbury Laboratory

Density Functional Theory. Martin Lüders Daresbury Laboratory Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei

More information

Comparison of various abinitio codes used in periodic calculations

Comparison of various abinitio codes used in periodic calculations Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,

More information

Multi-Scale Modeling from First Principles

Multi-Scale Modeling from First Principles m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations

More information

The electronic structure of materials 2 - DFT

The electronic structure of materials 2 - DFT Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie

More information

Behind the "exciting" curtain: The (L)APW+lo method

Behind the exciting curtain: The (L)APW+lo method Behind the "exciting" curtain: The (L)APW+lo method Aug 7, 2016 Andris Gulans Humboldt-Universität zu Berlin Kohn-Sham equation Potential due to nuclei Exchange-correlation potential Potential due to electron

More information

DFT: Exchange-Correlation

DFT: Exchange-Correlation DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,

More information

Teoría del Funcional de la Densidad (Density Functional Theory)

Teoría del Funcional de la Densidad (Density Functional Theory) Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals

More information

Introduction to Density Functional Theory

Introduction to Density Functional Theory 1 Introduction to Density Functional Theory 21 February 2011; V172 P.Ravindran, FME-course on Ab initio Modelling of solar cell Materials 21 February 2011 Introduction to DFT 2 3 4 Ab initio Computational

More information

Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić

Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824

More information

DFT: Exchange-Correlation

DFT: Exchange-Correlation DFT: Exchange-Correlation Local functionals, exact exchange and other post-dft methods Paul Tulip Centre for Materials Physics Department of Physics University of Durham Outline Introduction What is exchange

More information

Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory

Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.

More information

Band calculations: Theory and Applications

Band calculations: Theory and Applications Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation

More information

DFT in practice : Part II. Ersen Mete

DFT in practice : Part II. Ersen Mete pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials

More information

Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT

Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,

More information

Ab initio Electronic Structure

Ab initio Electronic Structure Ab initio Electronic Structure M. Alouani IPCMS, UMR 7504, Université Louis Pasteur, Strasbourg France http://www-ipcms.u-strasbg.fr In coll. with: B. Arnaud, O. Bengone, Y. Dappe, and S. Lebègue 1965

More information

DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY

DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES

More information

FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY

FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY C.A. Madu and B.N Onwuagba Department of Physics, Federal University of Technology Owerri, Nigeria

More information

Kevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2.

Kevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2. Quantum Monte Carlo Simulations of a Single Iron Impurity in MgO Kevin Driver 1 Shuai Zhang 1 Burkhard Militzer 1 R. E. Cohen 2 1 Department of Earth & Planetary Science University of California, Berkeley

More information

The Plane-Wave Pseudopotential Method

The Plane-Wave Pseudopotential Method Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste

More information

Practical Guide to Density Functional Theory (DFT)

Practical Guide to Density Functional Theory (DFT) Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between

More information

Solid State Theory: Band Structure Methods

Solid State Theory: Band Structure Methods Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.

More information

Strained Silicon, Electronic Band Structure and Related Issues.

Strained Silicon, Electronic Band Structure and Related Issues. Strained Silicon, Electronic Band Structure and Related Issues. D. Rideau, F. Gilibert, M. Minondo, C. Tavernier and H. Jaouen STMicroelectronics,, Device Modeling 850 rue Jean Monnet, BP 16, F-38926 Crolles

More information

Key concepts in Density Functional Theory (II)

Key concepts in Density Functional Theory (II) Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université

More information

An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method

An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig

More information

Norm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria

Norm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized

More information

Key concepts in Density Functional Theory (I) Silvana Botti

Key concepts in Density Functional Theory (I) Silvana Botti From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre

More information

Electrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation

Electrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project

More information

Density Functional Theory for Electrons in Materials

Density Functional Theory for Electrons in Materials Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for

More information

3: Density Functional Theory

3: Density Functional Theory The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London

More information

nanohub.org learning module: Prelab lecture on bonding and band structure in Si

nanohub.org learning module: Prelab lecture on bonding and band structure in Si nanohub.org learning module: Prelab lecture on bonding and band structure in Si Ravi Vedula, Janam Javerhi, Alejandro Strachan Center for Predictive Materials Modeling and Simulation, School of Materials

More information

Combining quasiparticle energy calculations with exact-exchange density-functional theory

Combining quasiparticle energy calculations with exact-exchange density-functional theory Combining quasiparticle energy calculations with exact-exchange density-functional theory Patrick Rinke 1, Abdallah Qteish 1,2, Jörg Neugebauer 1,3,4, Christoph Freysoldt 1 and Matthias Scheffler 1 1 Fritz-Haber-Institut

More information

Designed nonlocal pseudopotentials for enhanced transferability

Designed nonlocal pseudopotentials for enhanced transferability PHYSICAL REVIEW B VOLUME 59, NUMBER 19 15 MAY 1999-I Designed nonlocal pseudopotentials for enhanced transferability Nicholas J. Ramer and Andrew M. Rappe Department of Chemistry and Laboratory for Research

More information

Introduction to First-Principles Method

Introduction to First-Principles Method Joint ICTP/CAS/IAEA School & Workshop on Plasma-Materials Interaction in Fusion Devices, July 18-22, 2016, Hefei Introduction to First-Principles Method by Guang-Hong LU ( 吕广宏 ) Beihang University Computer

More information

Improved Electronic Structure and Optical Properties of sp-hybridized Semiconductors Using LDA+U SIC

Improved Electronic Structure and Optical Properties of sp-hybridized Semiconductors Using LDA+U SIC 286 Brazilian Journal of Physics, vol. 36, no. 2A, June, 2006 Improved Electronic Structure and Optical Properties of sp-hybridized Semiconductors Using LDA+U SIC Clas Persson and Susanne Mirbt Department

More information

Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn

Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,

More information

ELECTRONIC AND MAGNETIC PROPERTIES OF BERKELIUM MONONITRIDE BKN: A FIRST- PRINCIPLES STUDY

ELECTRONIC AND MAGNETIC PROPERTIES OF BERKELIUM MONONITRIDE BKN: A FIRST- PRINCIPLES STUDY ELECTRONIC AND MAGNETIC PROPERTIES OF BERKELIUM MONONITRIDE BKN: A FIRST- PRINCIPLES STUDY Gitanjali Pagare Department of Physics, Sarojini Naidu Govt. Girls P. G. Auto. College, Bhopal ( India) ABSTRACT

More information

Prerequisites for reliable modeling with first-principles methods. P. Kratzer Fritz-Haber-Institut der MPG D Berlin-Dahlem, Germany

Prerequisites for reliable modeling with first-principles methods. P. Kratzer Fritz-Haber-Institut der MPG D Berlin-Dahlem, Germany Prerequisites for reliable modeling with first-principles methods P. Kratzer Fritz-Haber-Institut der MPG D-14195 Berlin-Dahlem, Germany Prerequisites for modeling (I) Issues to consider when applying

More information

Electronic structure calculations with GPAW. Jussi Enkovaara CSC IT Center for Science, Finland

Electronic structure calculations with GPAW. Jussi Enkovaara CSC IT Center for Science, Finland Electronic structure calculations with GPAW Jussi Enkovaara CSC IT Center for Science, Finland Basics of density-functional theory Density-functional theory Many-body Schrödinger equation Can be solved

More information

3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional

More information

Some surprising results of the Kohn-Sham Density Functional

Some surprising results of the Kohn-Sham Density Functional arxiv:1409.3075v1 [cond-mat.mtrl-sci] 10 Sep 2014 Some surprising results of the Kohn-Sham Density Functional L. G. Ferreira 1, M. Marques 2, L. K. Teles 2, R. R. Pelá 2 1 Instituto de Física, Universidade

More information

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)

More information

Many electrons: Density functional theory Part II. Bedřich Velický VI.

Many electrons: Density functional theory Part II. Bedřich Velický VI. Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to

More information

Modified Becke-Johnson (mbj) exchange potential

Modified Becke-Johnson (mbj) exchange potential Modified Becke-Johnson (mbj) exchange potential Hideyuki Jippo Fujitsu Laboratories LTD. 2015.12.21-22 OpenMX developer s meeting @ Kobe Overview: mbj potential The semilocal exchange potential adding

More information

Introduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique

Introduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_

More information

Pseudopotential generation and test by the ld1.x atomic code: an introduction

Pseudopotential generation and test by the ld1.x atomic code: an introduction and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)

More information

Session 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education

Session 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools

More information

MODULE 2: QUANTUM MECHANICS. Principles and Theory

MODULE 2: QUANTUM MECHANICS. Principles and Theory MODULE 2: QUANTUM MECHANICS Principles and Theory You are here http://www.lbl.gov/cs/html/exascale4energy/nuclear.html 2 Short Review of Quantum Mechanics Why do we need quantum mechanics? Bonding and

More information

Comparison of first-principles methods / codes

Comparison of first-principles methods / codes Comparison of first-principles methods / codes Benchmarks run by Jörg Behler (DMol), Jutta Rogal (APW+lo), Mira Todorova (LAPW), Cathy Stampfl (FHI98md and sfhingx)... and Peter Blaha and Peter Blöchl

More information

Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions

Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions Daniele Stradi www.quantumwise.com daniele.stradi@quantumwise.com Introduction

More information

Key concepts in Density Functional Theory (II) Silvana Botti

Key concepts in Density Functional Theory (II) Silvana Botti Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:

More information

Integrated Computational Materials Engineering Education

Integrated Computational Materials Engineering Education Integrated Computational Materials Engineering Education Lecture on Density Functional Theory An Introduction Mark Asta Dept. of Materials Science and Engineering, University of California, Berkeley &

More information

Theory of Pseudopotentials. Outline of Talk

Theory of Pseudopotentials. Outline of Talk Theory of Pseudopotentials David Vanderbilt Rutgers University Outline of Talk Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties

More information

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model

More information

Exchange-Correlation Functional

Exchange-Correlation Functional Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological

More information

Table of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures

Table of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures Table of Contents Table of Contents Spin-orbit splitting of semiconductor band structures Relavistic effects in Kohn-Sham DFT Silicon band splitting with ATK-DFT LSDA initial guess for the ground state

More information

Electronic structure and transport in silicon nanostructures with non-ideal bonding environments

Electronic structure and transport in silicon nanostructures with non-ideal bonding environments Purdue University Purdue e-pubs Other Nanotechnology Publications Birck Nanotechnology Center 9-15-2008 Electronic structure and transport in silicon nanostructures with non-ideal bonding environments

More information

Why use pseudo potentials?

Why use pseudo potentials? Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10

More information

Pseudopotential methods for DFT calculations

Pseudopotential methods for DFT calculations Pseudopotential methods for DFT calculations Lorenzo Paulatto Scuola Internazionale Superiore di Studi Avanzati and CNR-INFM DEMOCRITOS National Simulation Center Tieste Italy July 9, 2008 Outline pseudopotential

More information

Nonlocal exchange correlation in screened-exchange density functional methods

Nonlocal exchange correlation in screened-exchange density functional methods Nonlocal exchange correlation in screened-exchange density functional methods Byounghak Lee and Lin-Wang Wang Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California

More information

exciting in a nutshell

exciting in a nutshell http://exciting-code.org exciting in a nutshell Pasquale Pavone Humboldt-Universität zu Berlin http://exciting-code.org exciting in a (coco)nutshell Pasquale Pavone Humboldt-Universität zu Berlin Outline

More information

Two implementations of the Projector Augmented Wave (PAW) formalism

Two implementations of the Projector Augmented Wave (PAW) formalism Introduction The tools available for detailed first-principles studies of materials have benefited enormously from the development of several international collaborations engaged in developing open source

More information

The Linearized Augmented Planewave (LAPW) Method

The Linearized Augmented Planewave (LAPW) Method The Linearized Augmented Planewave (LAPW) Method David J. Singh Oak Ridge National Laboratory E T [ ]=T s [ ]+E ei [ ]+E H [ ]+E xc [ ]+E ii {T s +V ks [,r]} I (r)= i i (r) Need tools that are reliable

More information

Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics

Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Sangeeta Sharma 1,2, J. K. Dewhurst 3, C. Ambrosch-Draxl 4, S. Pittalis 2, S. Kurth 2, N. Helbig 2, S. Shallcross

More information

Solid State Theory: Band Structure Methods

Solid State Theory: Band Structure Methods Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:00-12:30 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Who am I? Assistant Professor, Institute for Theoretical and Computational Physics,

More information

Pseudo potential exercises

Pseudo potential exercises Pseudo potential exercises Johan M. Carlsson Fritz-Haber-Institut der Max-Planck-Gesellschaft D-14195 Berlin Introduction Castep contains an "on the fly" OTF-pseudo potential generator that can be used

More information

Density Functional Theory (DFT)

Density Functional Theory (DFT) Density Functional Theory (DFT) An Introduction by A.I. Al-Sharif Irbid, Aug, 2 nd, 2009 Density Functional Theory Revolutionized our approach to the electronic structure of atoms, molecules and solid

More information

All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana?

All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana? All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana? SFB 484, Teilprojekt D6 October 5, 2007 Outline 1 2 3 Outline 1 2 3 Outline 1 2 3 Outline 1 2 3 Back in the 1930 s... John C. Slater

More information

1. Hydrogen atom in a box

1. Hydrogen atom in a box 1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom

More information

The Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus

The Nature of the Interlayer Interaction in Bulk. and Few-Layer Phosphorus Supporting Information for: The Nature of the Interlayer Interaction in Bulk and Few-Layer Phosphorus L. Shulenburger, A.D. Baczewski, Z. Zhu, J. Guan, and D. Tománek, Sandia National Laboratories, Albuquerque,

More information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

OVERVIEW OF QUANTUM CHEMISTRY METHODS OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density

More information

Set the initial conditions r i. Update neighborlist. Get new forces F i

Set the initial conditions r i. Update neighborlist. Get new forces F i v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )

More information

CHEM6085: Density Functional Theory

CHEM6085: Density Functional Theory Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron

More information

Quantum Mechanical Simulations

Quantum Mechanical Simulations Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock

More information

Projector augmented wave Implementation

Projector augmented wave Implementation Projector augmented wave Implementation Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 = Projector augmented wave +

More information

Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation

Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Multi-wavelets: High accuracy

More information

Supporting Information for Interfacial Effects on. the Band Edges of Functionalized Si Surfaces in. Liquid Water

Supporting Information for Interfacial Effects on. the Band Edges of Functionalized Si Surfaces in. Liquid Water Supporting Information for Interfacial Effects on the Band Edges of Functionalized Si Surfaces in Liquid Water Tuan Anh Pham,,, Donghwa Lee, Eric Schwegler, and Giulia Galli, Department of Chemistry, University

More information

The Plane-wave Pseudopotential Method

The Plane-wave Pseudopotential Method The Plane-wave Pseudopotential Method k(r) = X G c k,g e i(g+k) r Chris J Pickard Electrons in a Solid Nearly Free Electrons Nearly Free Electrons Nearly Free Electrons Electronic Structures Methods Empirical

More information

Phonon wavefunctions and electron phonon interactions in semiconductors

Phonon wavefunctions and electron phonon interactions in semiconductors Phonon wavefunctions and electron phonon interactions in semiconductors Bartomeu Monserrat bm418@cam.ac.uk University of Cambridge Quantum Monte Carlo in the Apuan Alps VII QMC in the Apuan Alps VII Bartomeu

More information

1 Density functional theory (DFT)

1 Density functional theory (DFT) 1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The

More information

file://c:\documents and Settings\lhebditc\Desktop\EAC-116.htm

file://c:\documents and Settings\lhebditc\Desktop\EAC-116.htm Page 1 of 40-1- B3.2 Quantum structural methods for the solid state and surfaces Frank Starrost and Emily A Carter B3.2.1 INTRODUCTION We are entering an era when condensed matter chemistry and physics

More information

Method development at SUNCAT in general

Method development at SUNCAT in general Bayesian error estimation functionals and further method development at SUNCAT GPAW 2016 University of Jyväskylä June 8th, 2016 Johannes Voss vossj@stanford.edu Method development at SUNCAT in general

More information

Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach

Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach Stability, Composition and Function of Palladium Surfaces in Oxidizing Environments: A First-Principles Statistical Mechanics Approach von Diplom-Chemikerin Jutta Rogal im Fachbereich Physik der Freien

More information

CHAPTER 3 WIEN2k. Chapter 3 : WIEN2k 50

CHAPTER 3 WIEN2k. Chapter 3 : WIEN2k 50 CHAPTER 3 WIEN2k WIEN2k is one of the fastest and reliable simulation codes among computational methods. All the computational work presented on lanthanide intermetallic compounds has been performed by

More information

GW Many-Body Theory for Electronic Structure. Rex Godby

GW Many-Body Theory for Electronic Structure. Rex Godby GW Many-Body Theory for Electronic Structure Rex Godby Outline Lecture 1 (Monday) Introduction to MBPT The GW approximation (non-sc and SC) Implementation of GW Spectral properties Lecture 2 (Tuesday)

More information

Photoelectronic properties of chalcopyrites for photovoltaic conversion:

Photoelectronic properties of chalcopyrites for photovoltaic conversion: Photoelectronic properties of chalcopyrites for photovoltaic conversion: self-consistent GW calculations Silvana Botti 1 LSI, CNRS-CEA-École Polytechnique, Palaiseau, France 2 LPMCN, CNRS-Université Lyon

More information

Key concepts in Density Functional Theory

Key concepts in Density Functional Theory From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem

More information

Performance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26

Performance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26 Performance of hybrid density functional methods, screened exchange and EXX-OEP methods in the PAW approach Georg Kresse, J Paier, R Hirschl, M Marsmann Institut für Materialphysik and Centre for Computational

More information

Theoretical and experimental factors affecting measurements of semiconductor mean inner potentials

Theoretical and experimental factors affecting measurements of semiconductor mean inner potentials Journal of Physics: Conference Series Theoretical and experimental factors affecting measurements of semiconductor mean inner potentials Recent citations - Mapping Charge Distribution in Single PbS Core

More information

MBPT and TDDFT Theory and Tools for Electronic-Optical Properties Calculations in Material Science

MBPT and TDDFT Theory and Tools for Electronic-Optical Properties Calculations in Material Science MBPT and TDDFT Theory and Tools for Electronic-Optical Properties Calculations in Material Science Dott.ssa Letizia Chiodo Nano-bio Spectroscopy Group & ETSF - European Theoretical Spectroscopy Facility,

More information

Introduction to DFT and its Application to Defects in Semiconductors

Introduction to DFT and its Application to Defects in Semiconductors Introduction to DFT and its Application to Defects in Semiconductors Noa Marom Physics and Engineering Physics Tulane University New Orleans The Future: Computer-Aided Materials Design Can access the space

More information

André Schleife Department of Materials Science and Engineering

André Schleife Department of Materials Science and Engineering André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic

More information

Magnetism in transition metal oxides by post-dft methods

Magnetism in transition metal oxides by post-dft methods Magnetism in transition metal oxides by post-dft methods Cesare Franchini Faculty of Physics & Center for Computational Materials Science University of Vienna, Austria Workshop on Magnetism in Complex

More information

Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.

Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm

More information

Introduction of XPS Absolute binding energies of core states Applications to silicone Outlook

Introduction of XPS Absolute binding energies of core states Applications to silicone Outlook Core level binding energies in solids from first-principles Introduction of XPS Absolute binding energies of core states Applications to silicone Outlook TO and C.-C. Lee, Phys. Rev. Lett. 118, 026401

More information

All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana?

All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana? All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana? Center for Electronic Correlations and Magnetism Institute for Physics, University of Augsburg February 4, 2008 Outline 1 2 3 Outline

More information