Pseudopotentials for hybrid density functionals and SCAN

Size: px
Start display at page:

Download "Pseudopotentials for hybrid density functionals and SCAN"

Transcription

1 Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania

2 Why do we need pseudopotentials? N e - n e - Valence electrons - Dominate bonding - Smooth exponential decay Core electrons - Insensitive to chemical environment - Oscillation near nucleus Pseudopotential Replace core region with a finite shallow potential Reduce the computational cost Preserve the accuracy Frozen-core approximation - Combine core electron and nucleus to effective core (pseudo ion) - Explicitly keep valence electrons (N > n) ψ ps ψ AE V ps -Z eff /r 2

3 Pseudopotential developments Hamann-Schlüter-Chiang norm conserving pseudopotential (1979) Generalized norm conserving pseudopotential (1989) RRKJ optimized pseudopotential (1990) Troullier-Martine soft pseudopotential (1990) Ultrasoft pseudopotential (1990) Projector-augmented waves potentials (1994) Links to pseudopotential databases

4 Density functional calculation vs. pseudopotential development Unoccupied {Φ i } ε x " n n n Generalized RPA Hybrid-GGA Meta-GGA GGA LSD Pseudopotentials still at GGA level Hartree World J. Chem. Phys. 123, (2005) 4

5 Hybrid density functional pseudopotentials The exchange-correlation functional must agree between pseudopotential and plane-wave calculation Non self-consistent PBE0 pseudopotential did not show improved accuracy [Phys. Rev. B 80, (2009)] First self-consistent construction of PBE0 pseudopotentials 5

6 Hartree-Fock formalism The nonlocal Hamiltonian ˆT + ˆVion + ˆV H [{ n 0 l 0 }]+ nl ˆV x [{ n 0 l 0 }] nl(r) = nl nl (r) Hartree potential ˆV H [{ n0 l 0 }] nl i = X n 0l 0 Z d 3 r 0 n 0 l 0 (r 0 ) 2 r r 0 nl(r) 6

7 Hartree-Fock formalism The nonlocal Hamiltonian ˆT + ˆVion + ˆV H [{ n 0 l 0 }]+ nl ˆV x [{ n 0 l 0 }] nl(r) = nl nl (r) Exchange potential operator ˆV x nl [{ n0 l 0 }] nl i = X Z n 0 l 0 d 3 r 0 n 0 l 0 (r 0 ) n 0 l 0 (r) r r 0 nl(r 0 ) 7

8 Hartree-Fock formalism The nonlocal Hamiltonian ˆT + ˆVion + ˆV H [{ n 0 l 0 }]+ nl ˆV x [{ n 0 l 0 }] nl(r) = nl nl (r) Hartree and exchange potentials are orbital dependent Treated as spherical average over all orbital configurations via Slater integral X E av = 1 N 2S+1 L J E[{ }2S+1 L J ] E.g., for C atom, valence configurations N=15 with 5 allowed term symbols: 3 P 2 (5), 3 P 1 (3), 3 P 0 (1), 1 D 2 (5), 1 S 0 (1) E av = 1 15 [5E{ }3 P 2 +3E{ }3 P 1 + E{ }3 P 0 +5E{ }1 D 2 + E{ }1 S 0 ] 8

9 Hartree-Fock formalism The nonlocal Hamiltonian ˆT + ˆVion + ˆV H [{ n 0 l 0 }]+ nl ˆV x [{ n 0 l 0 }] nl(r) = nl nl (r) Slater integral is used to obtain the average energy of Hartree and exchange terms E av = I[{ nl (r)}]+f [{ nl (r)}]+g[{ nl (r)}] [Slater integral J. C. Slater (1960)] I[{φ nl (r)}], F[{φ nl (r)}], G[{φ nl (r)}] are slater integrals containing radial wavefunctions with different n,l quantum numbers 9

10 Hybrid density functional formalism PBE0 formalism Mix the HF exchange with PBE exchange E PBE0 xc Combine HF formalism and Kohn-Sham equation apple ˆT + ˆV ion + ˆV H [{ n 0 l 0 }]+ 1 4 = aex HF +(1 a)ex PBE nl ˆV x [{ n 0 l 0 }] E PBE c a =0.25 ˆV PBE x Invert the radial equation to obtain a screened pseudopotential + [JCP 105, 9982 (1996)] ˆV PBE c nl (r) = nl nl (r) ˆV nl scr(r) = nl (r) l(l + 1) 2r val (r) d 2 dr 2 [ val(r)] 10

11 Hybrid functional all-electron solver Solve apple ˆT + ˆV ion + ˆV H [{ n 0 l 0}]+1 4 nl ˆV x [{ n 0 l 0}]+3 4 ˆV PBE x + ˆV PBE c nl (r) = nl nl (r) iteratively Start with an initial guess for the wavefunctions ψ (e.g. hydrogenic wavefunctions) Construct Hartree, exact-exchange potential, and DFT exchange-correlation potentials V H, V X, V XC ψ ps Solve the modified HF differential equations to obtain updated wavefunctions ψ, and eigenvalues ε ψ AE V ps -Z eff /r Orthogonalize the wavefunctions ψ 11

12 Pseudopotential construction Constructing the pseudowavefunctions using the RRKJ method Use spherical Bessel functions to construct Pseudowavefunctions [Phys. Rev. B 44, (1990)] ψ ps ψ AE Inverting the Schrödinger equation to obtain V PSP V ps -Z eff /r Design non-local transformation on V PSP [Phys. Rev. B 59, 1247 (1999)] Descreening to obtain V ion 12

13 Pseudopotential construction Constructing the pseudowavefunctions using the RRKJ method Use spherical Bessel functions to construct Pseudowavefunctions [Phys. Rev. B 44, (1990)] ψ ps ψ AE Inverting the Schrödinger equation to obtain V PSP V ps -Z eff /r Design non-local transformation on V PSP [Phys. Rev. B 59, 1247 (1999)] V PS ion = V PS l (r) Y [{ val }](r) 1 4 X[{ val }](r) l(r) 3 4 V x PBE (r) Vc PBE (r) 13

14 Use of OPIUM Only one input file: *.param How to run?./opium f f.log ae ps nl ke rpt upf plot vi [Atom] F [Pseudo] opt [Optinfo] [Configs] 2 # # [XC] hy Atomic reference configuration Pseudopotential radial cut off Pseudowavefunction settings Configurations for transferability test Exchangecorrelation method opium.sourceforge.net/index.html 14

15 PBE0 pseudopotential performance bond lengths Relative Error (%) 2 1 PBE0@PBE PBE0 Relative Error (%) 0 H 2 LiH FH Li 2 LiF CO N 0 2 F 2 BeH CH NH OH CN NO O PBE0@PBE PBE0 err% = l l AE l AE 100% Closed shell Open shell Overall small (but systematic) improvements for bond lengths without further tuning the potential parameters Influence tends to be large in open-shell systems All-electron (AE) calculations are performed using FHI-aims 15

16 PBE0 pseudopotential performance HOMO-LUMO gaps Relative Error (%) H 2 LiH FH Li 2 LiF CO N 2 F 2 Closed shell %err = PBE0@PBE PBE0 AE AE 100 Relative Error (%) PBE0@PBE PBE0 0 BeH CH NH OH NO O 2 Open shell The error in HOMO-LUMO gaps is reduced using PBE0 pseudopotentials by 2% on average 16

17 PBE0 pseudopotential performance bond formation energies Relative error (%) E AE err% = EPW E AE 100% PBE0@PBE PBE0-20 H 2 LiH BeH CH NH OH FH Li 2 LiF CN CO N 2 NO O 2 F 2 Formation energies are improved by up to 4 % using PBE0 pseudopotentials 17

18 PBE0 pseudopotential performance - solids Crystal PBE-PBE0 (err%) PBE0 (err%) All Electron Calc. Lattice constant (Å) PBE0 Si (0.073) (-0.037) GaN (0.066) (0.022) Band Gap (ev) Si 1.79 (9.82) 1.78 (9.20) 1.63 GaN 3.58 (1.13) 3.56 (0.56) 3.54 Lattice constants and band gaps are improved for solids 18

19 Implementation of SCAN SCAN (Strongly constrained and Appropriately Normed) density functional is a meta-gga functional that obeys all currently known constraints Achieves remarkable accuracy for lattice constants of simple solids and weak interactions The computational time is similar to GGA Implemented in FHI-AIMS, VASP, (and pw.x) Phys. Rev. Lett. 115, (2015) 19

20 Implementation of SCAN Start with an initial guess for the wavefunctions ψ (e.g. hydogenic wavefunctions) Initialize the spin-dependent charge density by ρ i = % " ψ i 2 Construct Hartree, and SCAN exchange-correlation potentials V H, V XC Solve spin-dependent Kohn-Sham equation to obtain updated wavefunctions ψ and eigenvalues ε Implement SCAN exchange-correlation Develop a new solver to take care of the spinpolarization Integrate the wavefunctions ψ to obtain charge density ρ 20

21 Apply to real systems Convergence test G2 dataset with bond length and HOMO-LUMO gaps Solids lattice constants, bulk moduli and band gaps Comparison between all electron and plane-wave calculations 21

22 First stage benchmarking G2 data set bond lengths PBE PBEPBE0 BPE0 AE-PBE AE-PBE0 AE-SCAN H LiH BeH CH NH OH FH Li LiF CN CO N NO O F MARE PBE0 pseudopotential plane-wave calculations still show smaller errors on bond lengths using SCAN all-electron data as a reference. 22

23 Conclusion First implementation of fully sefl-consistent PBE0 pseudopotentials Open source OPIUM v4.0 will be published soon on sourceforge ( Systematic improvement of accuracy for PBE0 calculations on molecules and solids Including relativity for spin-orbit coupling studies To implement SCAN functional in OPIUM 23

24 Acknowledgement Dr. Andrew M. Rappe Dr. Liang Z. Tan Dr. Julian Gebhardt Department of Chemistry University of Pennsylvania 24

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

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

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

Pseudopotentials: design, testing, typical errors

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

More information

An Introduction to OPIUM

An Introduction to OPIUM An Introduction to OPIUM Andrew M Rappe Makineni Theoretical Laboratories Department of Chemistry University of Pennsylvania & Eric J Walter Center for Piezoelectrics by Design Department of Physics College

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

Open-Source Pseudopotential Interface/Unification Module (OPIUM): The Basic Ins and Outs of Operation

Open-Source Pseudopotential Interface/Unification Module (OPIUM): The Basic Ins and Outs of Operation Open-Source Pseudopotential Interface/Unification Module (OPIUM): The Basic Ins and Outs of Operation Irene K. Metz, Joseph W. Bennett and Sara E. Mason (Dated: May 31, 2018) Learning Objectives 1. Determine

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

arxiv: v1 [physics.chem-ph] 12 Jul 2017

arxiv: v1 [physics.chem-ph] 12 Jul 2017 Hybrid functional pseudopotentials Jing Yang, 1 Liang Z. Tan, 1 and Andrew M. Rappe 1 1 Department of Chemistry, University of Pennsylvania, arxiv:1707.04501v1 [physics.chem-ph] 12 Jul 2017 Philadelphia,

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

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

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

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

CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122

CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122 CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX

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

Density Functional Theory: from theory to Applications

Density Functional Theory: from theory to Applications Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation

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

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

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

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

Orbital Density Dependent Functionals

Orbital Density Dependent Functionals Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)

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

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

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

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

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

Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory

Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley

More information

Multi-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005

Multi-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson

More information

Introduction to DFTB. Marcus Elstner. July 28, 2006

Introduction to DFTB. Marcus Elstner. July 28, 2006 Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD

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

Introduction to Density Functional Theory

Introduction to Density Functional Theory Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary

More information

Density Functional Theory

Density Functional Theory Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function

More information

Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014

Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde

More information

Outline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary

Outline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!

More information

Plane waves, pseudopotentials and PAW. X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium

Plane waves, pseudopotentials and PAW. X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium Plane waves, pseudopotentials and PAW X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium 1 Basic equations in DFT Solve self-consistently the Kohn-Sham equation H ψ n = ε n ψ n!!! ρ(r

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

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

Intermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA

Intermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA Intermediate DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Intermediate DFT Lausanne12

More information

XYZ of ground-state DFT

XYZ of ground-state DFT XYZ of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) XYZ of ground-state

More information

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National

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

Advanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University

Advanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation

More information

Computational Chemistry I

Computational Chemistry I Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of

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

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012 2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.

More information

Introduction to Density Functional Theory in a Plane Wave + Pseudopotential Framework

Introduction to Density Functional Theory in a Plane Wave + Pseudopotential Framework Introduction to Density Functional Theory in a Plane Wave + Pseudopotential Framework Shobhana Narasimhan Theoretical Sciences Unit JNCASR, Bangalore shobhana@jncasr.ac.in School- From the chemical bond

More information

Strategies for Solving Kohn- Sham equations

Strategies for Solving Kohn- Sham equations Strategies for Solving Kohn- Sham equations Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 1 Appetizer: high-k oxides

More information

The Projector Augmented Wave method

The Projector Augmented Wave method The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An

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

Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set

Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless

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

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

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

Density Functional Theory - II part

Density Functional Theory - II part Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not

More information

Supporting Information for. Predicting the Stability of Fullerene Allotropes Throughout the Periodic Table MA 02139

Supporting Information for. Predicting the Stability of Fullerene Allotropes Throughout the Periodic Table MA 02139 Supporting Information for Predicting the Stability of Fullerene Allotropes Throughout the Periodic Table Qing Zhao 1, 2, Stanley S. H. Ng 1, and Heather J. Kulik 1, * 1 Department of Chemical Engineering,

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

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

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

Orbital dependent correlation potentials in ab initio density functional theory

Orbital dependent correlation potentials in ab initio density functional theory Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland

More information

DFT calculations of NMR indirect spin spin coupling constants

DFT calculations of NMR indirect spin spin coupling constants DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling

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

Recent advances in development of single-point orbital-free kinetic energy functionals

Recent advances in development of single-point orbital-free kinetic energy functionals PacifiChem 2010 p. 1/29 Recent advances in development of single-point orbital-free kinetic energy functionals Valentin V. Karasiev vkarasev@qtp.ufl.edu Quantum Theory Project, Departments of Physics and

More information

Density Functional Theory

Density Functional Theory Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using

More information

Pseudopotential. Meaning and role

Pseudopotential. Meaning and role Pseudopotential. Meaning and role Jean-Pierre Flament jean-pierre.flament@univ-lille.fr Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM) Université de Lille-Sciences et technologies MSSC2018

More information

GW quasiparticle energies

GW quasiparticle energies Chapter 4 GW quasiparticle energies Density functional theory provides a good description of ground state properties by mapping the problem of interacting electrons onto a KS system of independent particles

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

Self-consistent Field

Self-consistent Field Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced

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

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

Ab-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials

Ab-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials ISSN 0974-9373 Vol. 15 No.3 (2011) Journal of International Academy of Physical Sciences pp. 337-344 Ab-initio Electronic Structure Calculations of α, β and γ KNO 3 Energetic Materials Pradeep Jain and

More information

Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals

Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University

More information

1 Construction of norm-conserving semi-local pseudopotentials for Si

1 Construction of norm-conserving semi-local pseudopotentials for Si 1 Construction of norm-conserving semi-local pseudopotentials for Si As discussed in class, it is desirable to replace the effective interaction of the valence electrons with the ionic core, i.e. nucleus

More information

5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry.

5.61 Physical Chemistry Exam III 11/29/12. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.61 Physical Chemistry Exam III (1) PRINT your name on the cover page. (2) It is suggested that you READ THE ENTIRE EXAM before

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

Electron Emission from Diamondoids: a DMC Study. Neil D. Drummond Andrew J. Williamson Richard J. Needs and Giulia Galli

Electron Emission from Diamondoids: a DMC Study. Neil D. Drummond Andrew J. Williamson Richard J. Needs and Giulia Galli Electron Emission from Diamondoids: a DMC Study Neil D. Drummond Andrew J. Williamson Richard J. Needs and Giulia Galli October 18, 2005 1 Semiconductor Nanoparticles for Optoelectronic Devices (I) The

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

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

CHEM3023: Spins, Atoms and Molecules

CHEM3023: Spins, Atoms and Molecules CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators

More information

MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )

MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij ) MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation

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

PBS: FROM SOLIDS TO CLUSTERS

PBS: FROM SOLIDS TO CLUSTERS PBS: FROM SOLIDS TO CLUSTERS E. HOFFMANN AND P. ENTEL Theoretische Tieftemperaturphysik Gerhard-Mercator-Universität Duisburg, Lotharstraße 1 47048 Duisburg, Germany Semiconducting nanocrystallites like

More information

Electronic Structure Calculations, Density Functional Theory and its Modern Implementations

Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Tutoriel Big RENOBLE Electronic Structure Calculations, Density Functional Theory and its Modern Implementations Thierry Deutsch L_Sim - CEA renoble 19 October 2011 Outline 1 of Atomistic calculations

More information

Density functional theory in the solid state

Density functional theory in the solid state Density functional theory in the solid state Ari P Seitsonen IMPMC, CNRS & Universités 6 et 7 Paris, IPGP Department of Applied Physics, Helsinki University of Technology Physikalisch-Chemisches Institut

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

Computational Physics. J. M. Thijssen

Computational Physics. J. M. Thijssen Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical

More information

Electron Correlation - Methods beyond Hartree-Fock

Electron Correlation - Methods beyond Hartree-Fock Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014

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

QMC dissociation energy of the water dimer: Time step errors and backflow calculations

QMC dissociation energy of the water dimer: Time step errors and backflow calculations QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.

More information

The Schrödinger equation for many-electron systems

The Schrödinger equation for many-electron systems The Schrödinger equation for many-electron systems Ĥ!( x,, x ) = E!( x,, x ) 1 N 1 1 Z 1 Ĥ = " $ # " $ + $ 2 r 2 A j j A, j RAj i, j < i a linear differential equation in 4N variables (atomic units) (3

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

Introduction to density-functional theory. Emmanuel Fromager

Introduction to density-functional theory. Emmanuel Fromager Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut

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

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

How to generate a pseudopotential with non-linear core corrections

How to generate a pseudopotential with non-linear core corrections How to generate a pseudopotential with non-linear core corrections 14 12 AE core charge AE valence charge PS core charge PS valence charge 10 8 6 4 2 Objectives 0 0 0.5 1 1.5 2 2.5 3 Check whether the

More information

DFT EXERCISES. FELIPE CERVANTES SODI January 2006

DFT EXERCISES. FELIPE CERVANTES SODI January 2006 DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a

More information

Electronic structure, plane waves and pseudopotentials

Electronic structure, plane waves and pseudopotentials Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know

More information

Electronic Structure: Density Functional Theory

Electronic Structure: Density Functional Theory Electronic Structure: Density Functional Theory S. Kurth, M. A. L. Marques, and E. K. U. Gross Institut für Theoretische Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany (Dated:

More information

Understanding electron correlation energy through density functional theory

Understanding electron correlation energy through density functional theory 1 Understanding electron correlation energy through density functional theory by Teepanis Chachiyo, 1,2,* and Hathaithip Chachiyo 2 1 Department of Physics, Faculty of Science, Naresuan University, Phitsanulok

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