Theory of Pseudopotentials. Outline of Talk

Size: px
Start display at page:

Download "Theory of Pseudopotentials. Outline of Talk"

Transcription

1 Theory of Pseudopotentials David Vanderbilt Rutgers University Outline of Talk Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties Norm Conservation Transferability Tests Relativistic Case Computational Considerations: Softness Ultrasoft Pseudopotentials and PAW Resources Reference list Web resources

2 Motivation Bandstructure of Si: Basic idea of pseudopotentials Pseudo wavefunction All-electron wavefunction Pseudopotential Cutoff radius All-electron potential

3 Pseudopotentials: History Pseudopotentials: Terminology

4 Pseudopotentials: Terminology s, p, and d electrons all feel the same potential Pseudopotentials: Terminology (For molybdenum) s, p, and d electrons feel different potentials

5 Pseudopotentials: Terminology These terms vanish outside r c s, p, and d electrons feel different nonlocal operators First-principles pseudopotentials Philosophy:

6 First-principles pseudopotentials Outline of Talk Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties Norm Conservation Transferability Tests Relativistic Case Computational Considerations: Softness Ultrasoft Pseudopotentials and PAW Resources Reference list Web resources

7 First-principles PSP construction Same for r > r c First-principles PSP construction Beyond r c : Also! AE =! PS

8 First-principles PSP construction Scattering properties For AE and PS separately: Choose channel l and energy! Find solution of SE that is regular at the origin at this! Compare beyond r c If match! good scattering properties! =! bound! >! bound (Used for construction) (Used for testing)

9 Scattering properties Quantify: Logarithmic derivatives D R (!) D R (!) R R! bound Correct by construction at! bound Norm conservation at! =! bound (Used for construction) r c Charge contained in this region is the same (AE vs. PS)

10 Norm conservation " Scattering properties Fundamental advance of Hamann, Schlüter and Chang, 1979: If norm conservation is imposed, then pseudo D R (!) matches all-electron D R (!) to second order in (!! bound ) R! bound Slope automatically matches First-principles PSP construction

11 Example: Hamann, Schlüter, and Chang (Semilocal PSP), 1979 (Mo) Charge self-consistency in PSP construction

12 Transferability tests PSP was generated in reference configuration, e.g.: [core]s 2 p 2 for Si Now, pick a couple of excited configurations, e.g.: [core]sp 3 [core]s 2 p (+1 ion) For each excited configuration, compare: All-electron calculation Pseudopotential calculation using previously generated PSP Points of comparison: Total energies Energy eigenvalues Logarithmic derivatives Transferability tests Example: HSC pseudopotential for oxygen

13 Relativisitic pseudopotentials Relativisitic pseudopotentials Bachelet, Hamann & Schluter, PRB 26, 4199 (1982)

14 Issues of computational expense The expense is in the target calculation (PSP construction is extremely cheap) First consideration: Compatibility with FFT approach to H"? H "(G) = KE "(G) + (V")(G) FFT FFT "(r) Multiply on Real space mesh (V")(r) Issues of computational expense Cheap, but poor transferability Good transferability, but expensive ~ All modern calculations done this way with one or two projectors

15 Expense vs. accuracy Compare different functional forms: Computational expense Want Semilocal PSP Local PSP Nonlocal (KB) PSP Accuracy Improved softness Now assume nonlocal (KB): Computational expense: N PW Smaller r c Larger r c Smooth cutoff functions Accuracy

16 Softness and plane-wave convergence Usual HSC Vanderbilt, 1985 q cut! 11 au -1 q cut! 7 au -1 Softness and plane-wave convergence Apply maximal smoothness to V ps construction!vanderbilt, 1985 This was only marginally successful in lowering the cutoff needed for the wavefunction Apply maximal smoothness to " ps construction!rappe, Rabe, Kaxiras, Joannopoulos (RRKJ, 1990)!Troullier and Martins (TM, 1991) Much more successful These (especially TM) are standard kind of potentials in use today

17 Softness and plane-wave convergence Plane-wave convergence of energy for free Cu atom (RRKJ, 1990) Outline of Talk Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties Norm Conservation Transferability Tests Relativistic Case Computational Considerations: Softness Ultrasoft Pseudopotentials and PAW Resources Reference list Web resources

18 Ultrasoft pseudopotentials O 2p wavefunction (nodeless!) AE PS Overshoot required by Norm conservation a.u. Ultrasoft pseudopotentials TM or RRKJ for O, Cu, etc: N PW Smaller r c Larger r c Ultrasoft PP (Vanderbilt, 1990) Accuracy

19 Ultrasoft pseudopotentials O 2p wavefunction (nodeless!) AE Usual PSP Ultrasoft a.u. Ultrasoft pseudopotentials Covergence tests for oxygen USPP, r c = 1.2 au USPP, r c = 1.8 au

20 Ultrasoft pseudopotentials: Formalism Ultrasoft pseudopotentials: Formalism

21 Ultrasoft pseudopotentials: Formalism Ultrasoft pseudopotentials: Formalism Typically, #=(1,2) in each angular momentum channel:

22 Terminology and Comparison In current usage, PSPs are classified as either NCPP = Norm-conserving pseudopotentials USPP = Ultrasoft pseudopotentials However, remember that USPP are norm-conserving in a generalized sense Thus, they retain the good features of NCPP In fact, their accuracy is usually better than NCPP Warning: Extra coding required in solid-state code Not all code packages accept USPP Ultrasoft Pseudopotentials

23 USPP and PAW P.E. Blöchl, Projector Augmented-Wave Method PRB 50, (1994) G. Kresse and D. Joubert, From USPP to PAW PRB 59, 1758 (1999) Outline of Talk Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties Norm Conservation Transferability Tests Relativistic Case Computational Considerations: Softness Ultrasoft Pseudopotentials and PAW Resources Reference list Web resources

24 Resources: References Resources: Web Sites Jose Luis Martins site for Troullier-Martins potentials: Octopus web interface for pseudopotential generation Vanderbilt Ultrasoft Pseudopotential site:

25 Octopus Web Site Octopus Web Site

26 Octopus Web Site Octopus Web Site

27 Octopus Web Site Ultrasoft Pseudopotential Web Site

28 Ultrasoft Pseudopotential Web Site Summary Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties Norm Conservation Transferability Tests Relativistic Case Computational Considerations: Softness Ultrasoft Pseudopotentials and PAW Resources Reference list Web resources Talk will be posted on

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

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

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

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

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

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

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

PSEUDOPOTENTIALS FOR BAND STRUCTURE CALCULATIONS

PSEUDOPOTENTIALS FOR BAND STRUCTURE CALCULATIONS TMCSIII: Jan 2012, Leeds PSEUDOPOTENTIALS FOR BAND STRUCTURE CALCULATIONS Rita Magri Physics Department, University of Modena and Reggio Emilia, Modena, Italy CNR-Nano -S 3, Modena, Italy OUTLINE Evolution

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

Construction pseudo-potentials for the projector augmentedwave. CSCI699 Assignment 2

Construction pseudo-potentials for the projector augmentedwave. CSCI699 Assignment 2 Construction pseudo-potentials for the projector augmentedwave Method CSCI699 Assignment 2 Make Your Own PAW Pseudopotentials I. Briefly describe Pseudo-wavefunction (RRKJ2) Local pseudo-potential Non-local

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

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

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

On-the-fly pseudopotential generation in CASTEP

On-the-fly pseudopotential generation in CASTEP On-the-fly pseudopotential generation in CASTEP Chris J. Pickard School of Physics and Astronomy, University of St Andrews St Andrews, KY16 9SS, United Kingdom September 13, 2006 A quick tutorial Default

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 May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials

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

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

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

Pseudopotentials and Basis Sets. How to generate and test them

Pseudopotentials and Basis Sets. How to generate and test them Pseudopotentials and Basis Sets How to generate and test them Pseudopotential idea Atomic Si Core electrons highly localized very depth energy are chemically inert 1s 2 2s 2 2p 6 3s 2 3p 2 Valence wave

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

Projector-Augmented Wave Method:

Projector-Augmented Wave Method: Projector-Augmented Wave Method: An introduction Peter E. Blöchl Clausthal University of Technology Germany http://www.pt.tu-clausthal.de/atp/ 23. Juli 2003 Why PAW all-electron wave functions (EFG s,

More information

Notes on pseudopotential generation

Notes on pseudopotential generation Notes on pseudopotential generation Paolo Giannozzi Università di Udine URL: http://www.fisica.uniud.it/ giannozz October 23, 2017 Contents 1 Introduction 1 1.1 Who needs to generate a pseudopotential?................

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

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

Data structure of OpenMX and ADPACK

Data structure of OpenMX and ADPACK Data structure of OpenMX and ADPACK Structure of ADPACK Solving the 1D Dirac equation Pseudopotentials Localized basis functions Structure of OpenMX Input_std() Total energy Data structure for parallelization

More information

arxiv: v2 [cond-mat.mtrl-sci] 1 Oct 2013

arxiv: v2 [cond-mat.mtrl-sci] 1 Oct 2013 Generation of Projector Augmented-Wave atomic data: a 71 elements validated table in the XML format arxiv:1309.7274v2 [cond-mat.mtrl-sci] 1 Oct 2013 François Jollet a, Marc Torrent a, Natalie Holzwarth

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

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 17 Jan 1998

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 17 Jan 1998 Electronic polarization in the ultrasoft pseudopotential formalism arxiv:cond-mat/9801177v1 [cond-mat.mtrl-sci] 17 Jan 1998 David anderbilt Department of Physics and Astronomy, Rutgers University, Piscataway,

More information

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

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

Table S2. Pseudopotentials PBE 5.2 applied in the calculations using VASP

Table S2. Pseudopotentials PBE 5.2 applied in the calculations using VASP Supporting Information for Understanding the Adsorption of CuPc and ZnPc on Noble Metal Surfaces by Combining Quantum-Mechanical Modelling and Photoelectron Spectroscopy 1. Used vdw Coefficients PBE-vdW

More information

Homogeneous Electric and Magnetic Fields in Periodic Systems

Homogeneous Electric and Magnetic Fields in Periodic Systems Electric and Magnetic Fields in Periodic Systems Josef W. idepartment of Chemistry and Institute for Research in Materials Dalhousie University Halifax, Nova Scotia June 2012 1/24 Acknowledgments NSERC,

More information

Comparison of Different Methods and Codes: (L)APW, LMTO, PAW, Pseudo Potentials, Gaussians, etc. Jörg Behler

Comparison of Different Methods and Codes: (L)APW, LMTO, PAW, Pseudo Potentials, Gaussians, etc. Jörg Behler Comparison of Different Methods and Codes: (L)APW, LMTO, PAW, Pseudo Potentials, Gaussians, etc. Jörg Behler Fritz-Haber-Institut der Max-Planck-Gesellschaft Faradayweg 4-6, D-14195 Berlin, Germany Density-Functional

More information

Lecture on First-Principles Computation (11): The LAPW Method

Lecture on First-Principles Computation (11): The LAPW Method Lecture on First-Principles Computation (11): The LAPW Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息实验室 (Key Laboratory of Quantum Information, USTC) 2015-10-16 Recall: OPW and Pseudopotential Methods Problems:

More information

The Quantum ESPRESSO Software Distribution

The Quantum ESPRESSO Software Distribution The Quantum ESPRESSO Software Distribution The DEMOCRITOS center of Italian INFM is dedicated to atomistic simulations of materials, with a strong emphasis on the development of high-quality scientific

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

TWO-PION EXCHANGE PARITY-VIOLATING NN INTERACTION

TWO-PION EXCHANGE PARITY-VIOLATING NN INTERACTION TWO-PION EXCHANGE PARITY-VIOLATING NN INTERACTION Bertrand DESPLANQUES (LPSC, Grenoble) Some work with Chang Ho Hyun and Shung-ichi Ando (Suwon) MOTIVATIONS New work in the EFT formalism (Zhu et al.) -

More information

Time-dependent density functional perturbation theory

Time-dependent density functional perturbation theory Time-dependent density functional perturbation theory Iurii Timrov and Tommaso Gorni SISSA Scuola Internazionale Superiore di Studi Avanzati, Trieste, Italy Advanced Quantum ESPRESSO developers' meeting:

More information

Atomic orbitals of finite range as basis sets. Javier Junquera

Atomic orbitals of finite range as basis sets. Javier Junquera Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles

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

The Quantum ESPRESSO Software Distribution

The Quantum ESPRESSO Software Distribution The Quantum ESPRESSO Software Distribution Updated March 2018 This document introduces to the philosophy of Quantum ESPRESSO and describes its main features and capabilities. Further information can be

More information

Notes on pseudopotential generation

Notes on pseudopotential generation Notes on pseudopotential generation Paolo Giannozzi, Scuola Normale Superiore, Pisa e-mail: giannozz@nest.sns.it URL: http://www.nest.sns.it/~giannozz February 27, 2004 1 Introduction When I started to

More information

Basis sets for SIESTA. Emilio Artacho. Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge

Basis sets for SIESTA. Emilio Artacho. Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge Basis sets for SIESTA Emilio Artacho Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge Solving: Basis set Expand in terms of a finite set of basis functions

More information

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 8 Jul 2004

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 8 Jul 2004 Electronic structure methods: Augmented Waves, Pseudopotentials and the Projector Augmented Wave Method arxiv:cond-mat/0407205v1 [cond-mat.mtrl-sci] 8 Jul 2004 Peter E. Blöchl 1, Johannes Kästner 1, and

More information

quantum ESPRESSO stands for Quantum open-source Package for Research in Electronic Structure, Simulation, and Optimization

quantum ESPRESSO stands for Quantum open-source Package for Research in Electronic Structure, Simulation, and Optimization The quantum ESPRESSO distribution The IOM-DEMOCRITOS center of Italian CNR is dedicated to atomistic simulations of materials, with a strong emphasis on the development of high-quality scientific software

More information

Natural Language Processing. Classification. Features. Some Definitions. Classification. Feature Vectors. Classification I. Dan Klein UC Berkeley

Natural Language Processing. Classification. Features. Some Definitions. Classification. Feature Vectors. Classification I. Dan Klein UC Berkeley Natural Language Processing Classification Classification I Dan Klein UC Berkeley Classification Automatically make a decision about inputs Example: document category Example: image of digit digit Example:

More information

Non-Newtonian Fluids and Finite Elements

Non-Newtonian Fluids and Finite Elements Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations

More information

Dept of Mechanical Engineering MIT Nanoengineering group

Dept of Mechanical Engineering MIT Nanoengineering group 1 Dept of Mechanical Engineering MIT Nanoengineering group » Recap of HK theorems and KS equations» The physical meaning of the XC energy» Solution of a one-particle Schroedinger equation» Pseudo Potentials»

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

MODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO

MODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE

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

DISTRIBUTION STATEMENT A Approved for Public Release Distribution Unlimited. Progress Report to the Office of Naval Research

DISTRIBUTION STATEMENT A Approved for Public Release Distribution Unlimited. Progress Report to the Office of Naval Research Progress Report to the Office of Naval Research Grant No. N00014-02-1-0735 Period: September 15. 2004 - May 31, 2005 Program Manager: Dr. Morely Stone Theoretical Analysis of Conductance of Molecular Junctions

More information

Exercises on basis set generation Full control on the definition of the basis set functions: the PAO.Basis block. Javier Junquera

Exercises on basis set generation Full control on the definition of the basis set functions: the PAO.Basis block. Javier Junquera Exercises on basis set generation Full control on the definition of the basis set functions: the PAO.Basis block Javier Junquera 3.50 0.95 # scalefactor(izeta=1,nzeta) (opt) # scalefactor(izeta=1,nzeta)

More information

Body Centered Cubic Magnesium Niobium Hydride with Facile Room Temperature Absorption and Four Weight Percent Reversible Capacity

Body Centered Cubic Magnesium Niobium Hydride with Facile Room Temperature Absorption and Four Weight Percent Reversible Capacity Electronic Supplementary Information (ESI) for Energy & Environmental Science This journal is The Royal Society of Chemistry 212 Supporting Information Body Centered Cubic Magnesium Niobium Hydride with

More information

Supporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering

Supporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering Supporting Information Tuning Local Electronic Structure of Single Layer MoS2 through Defect Engineering Yan Chen, 1,2,,$, * Shengxi Huang, 3,6, Xiang Ji, 2 Kiran Adepalli, 2 Kedi Yin, 8 Xi Ling, 3,9 Xinwei

More information

GW-like approaches to quantum transport. Rex Godby

GW-like approaches to quantum transport. Rex Godby GW-like approaches to quantum transport Rex Godby Outline Introduction to the quantum transport problem Ab initio quantum conductance in the presence of e e interaction (TDDFT / MBPT) 2 + Bothersome aspects

More information

AB INITIO MODELLING TECHNIQUES APPLIED TO SILICON. P.R.BRIDDON Department of Physics, University of Newcastle, Newcastle upon Tyne, NE1 7RU, UK

AB INITIO MODELLING TECHNIQUES APPLIED TO SILICON. P.R.BRIDDON Department of Physics, University of Newcastle, Newcastle upon Tyne, NE1 7RU, UK AB INITIO MODELLING TECHNIQUES APPLIED TO SILICON P.R.BRIDDON Department of Physics, University of Newcastle, Newcastle upon Tyne, NE1 7RU, UK August 13, 1999 1 Introduction In this chapter we will consider

More information

Numerical atomic basis orbitals from H to Kr

Numerical atomic basis orbitals from H to Kr PHYSICAL REVIEW B 69, 195113 2004 Numerical atomic basis orbitals from H to Kr T. Ozaki Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science and Technology

More information

Human ages are named after materials - stone, bronze, iron, nuclear, silicon

Human ages are named after materials - stone, bronze, iron, nuclear, silicon MATERIALS ARE KEY TO SOCIETAL WELL BEING Human ages are named after materials - stone, bronze, iron, nuclear, silicon MATERIALS ARE KEY TO SOCIETAL WELL BEING We need novel materials for: Energy harvesting,

More information

Electrons in Crystals. Chris J. Pickard

Electrons in Crystals. Chris J. Pickard Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity

More information

ALMA: All-scale predictive design of heat management material structures

ALMA: All-scale predictive design of heat management material structures ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable

More information

QUANTUM CHEMISTRY PROJECT 1

QUANTUM CHEMISTRY PROJECT 1 Chemistry 460 Fall 2017 Dr. Jean M. Standard September 11, 2017 QUANTUM CHEMISTRY PROJECT 1 OUTLINE This project focuses on applications and solutions of quantum mechanical problems involving one-dimensional

More information

Time part of the equation can be separated by substituting independent equation

Time part of the equation can be separated by substituting independent equation Lecture 9 Schrödinger Equation in 3D and Angular Momentum Operator In this section we will construct 3D Schrödinger equation and we give some simple examples. In this course we will consider problems where

More information

This is an author-deposited version published in: Eprints ID: URL:

This is an author-deposited version published in:   Eprints ID: URL: Open Archive Toulouse Archive Ouverte (OATAO) OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible. This is an author-deposited

More information

ATOM User Manual. Version 3.2.8, 29 April Alberto García ICMAB-CSIC, Barcelona vfill

ATOM User Manual. Version 3.2.8, 29 April Alberto García ICMAB-CSIC, Barcelona vfill ATOM User Manual Version 3.2.8, 29 April 2008 Alberto García ICMAB-CSIC, Barcelona albertog@icmab.es vfill Contents 1 PREFACE 2 2 A PRIMER ON AB-INITIO PSEUDOPOTENTIALS 2 3 COMPILING THE PROGRAM 3 4 USING

More information

Implementation and performance of the frequency-dependent GW method within the PAW framework

Implementation and performance of the frequency-dependent GW method within the PAW framework PHYSICAL REVIEW B 74, 035101 2006 Implementation and performance of the frequency-dependent method within the PAW framework M. Shishkin and G. Kresse Institut für Materialphysik and Centre for Computational

More information

Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase

Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase Supplementary Information Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase Jia Chen, Ye-Fei Li, Patrick Sit, and Annabella Selloni Department of Chemistry,

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 Theory for Periodic Systems: The Concepts. Christian Ratsch

Electronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3

More information

Rustam Z. Khaliullin University of Zürich

Rustam Z. Khaliullin University of Zürich Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation

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

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

CP2K: the gaussian plane wave (GPW) method

CP2K: the gaussian plane wave (GPW) method CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org

More information

[100] directed Cu-doped h-coo Nanorods: Elucidation of. Growth Mechanism and Application to Lithium-Ion Batteries

[100] directed Cu-doped h-coo Nanorods: Elucidation of. Growth Mechanism and Application to Lithium-Ion Batteries Supplementary Information [100] directed Cu-doped h-coo Nanorods: Elucidation of Growth Mechanism and Application to Lithium-Ion Batteries Ki Min Nam, Young Cheol Choi, Sung Chul Jung, Yong-Il Kim, Mi

More information

arxiv: v1 [cond-mat.mtrl-sci] 12 Aug 2008

arxiv: v1 [cond-mat.mtrl-sci] 12 Aug 2008 A fusion of the LAPW and the LMTO methods: the augmented plane wave plus muffin-tin orbital (PMT) method Takao Kotani and Mark van Schilfgaarde School of Materials, Arizona State University, Tempe, AZ,

More information

Nuclear Structure and Reactions using Lattice Effective Field Theory

Nuclear Structure and Reactions using Lattice Effective Field Theory Nuclear Structure and Reactions using Lattice Effective Field Theory Dean Lee North Carolina State University Nuclear Lattice EFT Collaboration Frontiers of Nuclear Physics Kavli Institute for Theoretical

More information

Projector augmented-wave method: Application to relativistic spin-density functional theory

Projector augmented-wave method: Application to relativistic spin-density functional theory PHYSCAL REVEW B 8, 0756 00 Proector augmented-wave method: Application to relativistic spin-density functional theory Andrea Dal Corso nternational School for Advanced Studies (SSSA), Via Bonomea 65, 3436

More information

6.730 Physics for Solid State Applications

6.730 Physics for Solid State Applications 6.730 Physics for Solid State Applications Lecture 19: Motion of Electronic Wavepackets Outline Review of Last Time Detailed Look at the Translation Operator Electronic Wavepackets Effective Mass Theorem

More information

1 IMEM-CNR, U.O.S. Genova, Via Dodecaneso 33, Genova, IT. 2 Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, Genova, IT

1 IMEM-CNR, U.O.S. Genova, Via Dodecaneso 33, Genova, IT. 2 Dipartimento di Fisica, Università di Genova, Via Dodecaneso 33, Genova, IT Spontaneous Oxidation of Ni Nanoclusters on MgO Monolayers Induced by Segregation of Interfacial Oxygen. M. Smerieri 1, J. Pal 1,2, L. Savio 1*, L. Vattuone 1,2, R. Ferrando 1,3, S. Tosoni 4, L. Giordano

More information

The PWcond code: Complex bands, transmission, and ballistic conductance

The PWcond code: Complex bands, transmission, and ballistic conductance The PWcond code: Complex bands, transmission, and ballistic conductance SISSA and IOM-DEMOCRITOS Trieste (Italy) Outline 1 Ballistic transport: a few concepts 2 Complex band structure 3 Current of a Bloch

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

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

Atomic Models for Anionic Ligand Passivation of Cation- Rich Surfaces of IV-VI, II-VI, and III-V Colloidal Quantum Dots

Atomic Models for Anionic Ligand Passivation of Cation- Rich Surfaces of IV-VI, II-VI, and III-V Colloidal Quantum Dots Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2016 Electronic Supplementary Information Atomic Models for Anionic Ligand Passivation of Cation- Rich

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

Supplementary Information: Construction of Hypothetical MOFs using a Graph Theoretical Approach. Peter G. Boyd and Tom K. Woo*

Supplementary Information: Construction of Hypothetical MOFs using a Graph Theoretical Approach. Peter G. Boyd and Tom K. Woo* Electronic Supplementary Material ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2016 Supplementary Information: Construction of Hypothetical MOFs using a Graph Theoretical Approach

More information

Chromium Cluster on Defected Graphene

Chromium Cluster on Defected Graphene Chromium Cluster on Defected Graphene Yuhang Liu June 29, 2017 Abstract In this work, diffusion process of Cr atoms on two types of defected graphene and structure and magnetic properties of Cr cluster

More information

Accurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids

Accurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids Accurate kinetic energy evaluation in electronic structure calculations with localised functions on real space grids Published as Computer Physics Communications 140, 315 (2001) Chris-Kriton Skylaris,

More information

Berry Phases and Curvatures in Electronic-Structure Theory. David Vanderbilt Rutgers University

Berry Phases and Curvatures in Electronic-Structure Theory. David Vanderbilt Rutgers University Berry Phases and Curvatures in Electronic-Structure Theory David Vanderbilt Rutgers University Rahman prize for: Theory of polarization (King-Smith & Vanderbilt) Ultrasoft pseudopotentials Three quick

More information

Plane-wave DFT-LDA calculation of the electronic structure and absorption spectrum of copper

Plane-wave DFT-LDA calculation of the electronic structure and absorption spectrum of copper PHYSICAL REVIEW B, VOLUME 64, 9525 Plane-wave DFT-LDA calculation of the electronic structure and absorption spectrum of copper Andrea Marini, Giovanni Onida, and Rodolfo Del Sole Istituto Nazionale per

More information

Chiral effective field theory on the lattice: Ab initio calculations of nuclei

Chiral effective field theory on the lattice: Ab initio calculations of nuclei Chiral effective field theory on the lattice: Ab initio calculations of nuclei Nuclear Lattice EFT Collaboration Evgeny Epelbaum (Bochum) Hermann Krebs (Bochum) Timo Lähde (Jülich) Dean Lee (NC State)

More information

Supporting information for

Supporting information for Supporting information for Highly Active, Selective and Stable Direct H 2 O 2 Generation by Monodispersive PdAg Nanoalloy Jin Zhang, Bolong Huang, *# Qi Shao, and Xiaoqing Huang * College of Chemistry,

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

STATE-Senri. STATE-Senri

STATE-Senri. STATE-Senri Yoshitada MORIKAWA Department of Precision Science and Technology, Osaka University morikawa@prec.eng.osaka-u.ac.jp (Simulation Tool for Atom TEchnology) Density Functional Theory LDADGGADLDA+U Ultrasoft

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

Electronic Supplementary Information

Electronic Supplementary Information Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Electronic Supplementary Information Rational modifications on champion porphyrin

More information

Modern Physics. Unit 3: Operators, Tunneling and Wave Packets Lecture 3.3: The Momentum Operator

Modern Physics. Unit 3: Operators, Tunneling and Wave Packets Lecture 3.3: The Momentum Operator Modern Physics Unit 3: Operators, Tunneling and Wave Packets Lecture 3.3: The Momentum Operator Ron Reifenberger Professor of Physics Purdue University 1 There are many operators in QM H Ψ= EΨ, or ˆop

More information

Atomic orbitals of finite range as basis sets. Javier Junquera

Atomic orbitals of finite range as basis sets. Javier Junquera Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles

More information

Supplementary Materials for

Supplementary Materials for advances.sciencemag.org/cgi/content/full/3/2/e1601832/dc1 Supplementary Materials for Determination of band offsets, hybridization, and exciton binding in 2D semiconductor heterostructures Neil R. Wilson,

More information