The Plane-wave Pseudopotential Method

Size: px
Start display at page:

Download "The Plane-wave Pseudopotential Method"

Transcription

1 The Plane-wave Pseudopotential Method k(r) = X G c k,g e i(g+k) r Chris J Pickard

2 Electrons in a Solid

3 Nearly Free Electrons

4 Nearly Free Electrons

5 Nearly Free Electrons

6 Electronic Structures Methods

7 Empirical Pseudopotentials

8 Ab Initio Pseudopotentials

9 Ab Initio Pseudopotentials

10 Ab Initio Pseudopotentials

11 Ab Initio Pseudopotentials

12 Ultrasoft Pseudopotentials (Vanderbilt, 1990)

13 Projector Augmented Waves

14 Deriving the Pseudo Hamiltonian

15 PAW vs PAW C. G. van de Walle and P. E. Bloechl 1993 (I) PAW to calculate all electron properties from PSP calculations" P. E. Bloechl 1994 (II) The PAW electronic structure method" " The two are linked, but should not be confused

16 References David Vanderbilt, Soft self-consistent pseudopotentials in a generalised eigenvalue formalism, PRB " Kari Laasonen et al, Car-Parinello molecular dynamics with Vanderbilt ultrasoft pseudopotentials, PRB " P.E. Bloechl, Projector augmented-wave method, PRB " G. Kresse and D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, PRB

17 On-the-fly Pseudopotentials - All quantities for PAW (I) reconstruction available" - Allows automatic consistency of the potentials with functionals" - Even fewer input files required" - The code in CASTEP is a cut down (fast) pseudopotential code" - There are internal databases" - Other generation codes (Vanderbilt s original, OPIUM, LD1) exist

18 It all works!

19 Kohn-Sham Equations

20 Being (Self) Consistent

21 Eigenproblem in a Basis

22 Just a Few Atoms

23 In a Crystal

24 Plane-waves Represent the orbital in Fourier space" Z (r) = (g)e ig r d 3 r For a periodic system we have Bloch s Theorem" k(r) = X G c k,g e i(g+k) r Where the Gs are reciprocal lattice vectors and k is a symmetry label in the 1 st Brillouin zone"

25 Cutoff Energy Limit the number of plane wave components to those such that (G + k) 2 2 apple E cut This defines a length scale = p Ecut

26 Cutoff Energy "" - The minimum length scale depends on the elements in the system through PSPs" "" - Variational principle -> energy monotonically decreases to ground state energy as E cut increases (not strictly for USP/PAW)" "" - Converge required property with respect to cutoff energy "

27 Cutoff Energy "" - The minimum length scale depends on the elements in the system through PSPs" "" - Variational principle -> energy monotonically decreases to ground state energy as E cut increases (not strictly for USP/PAW)" "" - Converge required property with respect to cutoff energy "

28 Cutoff Energy "" - The minimum length scale depends on the elements in the system through PSPs" "" - Variational principle -> energy monotonically decreases to ground state energy as E cut increases (not strictly for USP/PAW)" "" - Converge required property with respect to cutoff energy "

29 Why Plane-waves? - Systematic convergence with respect to single parameter" - Non-local. Cover all space equally" - Cheap forces (no Pulay term)" - No basis-set superposition error " - Numerically efficient " - Use FFTs to transform between real and reciprocal space" - The obvious choice for periodic and works well for aperiodic systems" - Rapid development of new functionality

30 Progress, but

31 The Total Energy

32 Applying H

33 Evaluating the Energy

34 Iterative Diagonalisation

35 Minimisation

36 Plane-wave Codes CASTEP PEtot VASP PARATEC PWscf Da Capo Abinit CPMD Qbox fhi98md PWPAW SFHIngX DOD-PW NWchem Octopus JDFTx

37 CASTEP Original code due to Mike Payne rewritten by CDG, starting 1999, unveiled at 2001 CASTEP Workshop " Modern modular design, using best features of Fortran90 " Focus on robustness, ease of use, and core functionality " Particular strength in spectroscopy (NMR/IR/Raman/EELS) " Marketed commercially by Biovia (formerly Accelrys) " Low cost European Academic licence available

38 The CDG Matthew Segall Main author of the code specification, and responsible for all the low-level communications and basis set coding. Also the population analysis algorithms. " Matt Probert Responsible for the geometry optimisation and molecular dynamics coding, and keeper of the parameters module. " Stewart Clark Responsible for the band structure and exchange-correlation functional coding and keeper of the cell module. Also co-author of the linear response code. " Chris Pickard Responsible all pseudopotential coding, amongst other things. Also co-author of the NMR code. " Phil Hasnip Responsible for the electronic energy minimisation coding (including density mixing and ensemble DFT), band-parallelism, run-time profiling and parallel efficiency report, and wavefunction subroutines. Also works on optimisation and parallel performance. " Keith Refson Responsible for phonon and E-field calculations and co-author of the linear response code. " Jonathan Yates Co-author of the NMR code, responsible for spectral module (dos, optics, eels). " Mike Payne CASTEP creator and general overseer of the whole project.

6: Plane waves, unit cells, k- points and all that

6: Plane waves, unit cells, k- points and all that The Nuts and Bolts of First-Principles Simulation 6: Plane waves, unit cells, k- points and all that Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview

More information

The CASTEP Story. Mike Payne. Cavendish Laboratory University of Cambridge

The CASTEP Story. Mike Payne. Cavendish Laboratory University of Cambridge The CASTEP Story Mike Payne Cavendish Laboratory University of Cambridge mcp1@cam.ac.uk Outline of Talk 1. The foundations of quantum mechanics. 2. The unfulfilled promise. 3. The new dawn. i. Density

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

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

Universities of Leeds, Sheffield and York

Universities of Leeds, Sheffield and York promoting access to White Rose research papers Universities of Leeds, Sheffield and York http://eprints.whiterose.ac.uk/ Paper published in Zeitschrift für Kristallographie. White Rose Research Online

More information

First-Principles Vibrational spectroscopy and lattice dynamics of materials in the solid state

First-Principles Vibrational spectroscopy and lattice dynamics of materials in the solid state First-Principles Vibrational spectroscopy and lattice dynamics of materials in the solid state Keith Refson Computational Science and Engineering Department STFC Rutherford Appleton Laboratory First principles

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

Electron energy loss spectroscopy (EELS)

Electron energy loss spectroscopy (EELS) Electron energy loss spectroscopy (EELS) Phil Hasnip Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~pjh503 Many slides courtesy of Jonathan

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

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

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

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

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

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

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

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

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

and strong interlayer quantum confinement

and strong interlayer quantum confinement Supporting Information GeP3: A small indirect band gap 2D crystal with high carrier mobility and strong interlayer quantum confinement Yu Jing 1,3, Yandong Ma 1, Yafei Li 2, *, Thomas Heine 1,3 * 1 Wilhelm-Ostwald-Institute

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

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

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

Ab Initio modelling of structural and electronic. Matt Probert University of York

Ab Initio modelling of structural and electronic. Matt Probert University of York Ab Initio modelling of structural and electronic properties of semiconductors Matt Probert University of York http://www-users.york.ac.uk/~mijp1 Overview of Talk What is Ab Initio? What can we model? How

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

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

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

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

Band Structure Calculations; Electronic and Optical Properties

Band Structure Calculations; Electronic and Optical Properties ; Electronic and Optical Properties Stewart Clark University of Outline Introduction to band structures Calculating band structures using Castep Calculating optical properties Examples results Some applications

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

Electronic Structure Methodology 1

Electronic Structure Methodology 1 Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E

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

Algorithms and Computational Aspects of DFT Calculations

Algorithms and Computational Aspects of DFT Calculations Algorithms and Computational Aspects of DFT Calculations Part II Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational

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

Practical calculations using first-principles QM Convergence, convergence, convergence

Practical calculations using first-principles QM Convergence, convergence, convergence Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................

More information

Prediction of NMR parameters in the solid-state

Prediction of NMR parameters in the solid-state 1J PO Prediction of NMR parameters in the solid-state Jonathan Yates Materials Modelling Laboratory, Oxford Materials 2J PSi 1 Nuclear Magnetic Resonance Applied Field Zero Field B Population differences!

More information

Geometry Optimization

Geometry Optimization Geometry Optimization Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~mijp1 Overview of lecture n Background Theory n Hellman-Feynman

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

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

Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering

Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)

More information

NMR-CASTEP. Jonathan Yates. Cavendish Laboratory, Cambridge University. J-coupling cons. NMR-CASTEP York 2007 Jonathan Yates.

NMR-CASTEP. Jonathan Yates. Cavendish Laboratory, Cambridge University. J-coupling cons. NMR-CASTEP York 2007 Jonathan Yates. Jonathan Yates Cavendish Laboratory, Cambridge University 2 2 1 3 2 N1a-N7b a hemical shift Experiment Calculation [ppm] [ppm] -43.8-47.4 "(31P) 29 "( Si1) -213.3-214.8 "(29Si2) -217.0-218.7 29-119.1-128.6

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

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

Introduction to Parallelism in CASTEP

Introduction to Parallelism in CASTEP to ism in CASTEP Stewart Clark Band University of Durham 21 September 2012 Solve for all the bands/electrons (Band-) Band CASTEP solves the Kohn-Sham equations for electrons in a periodic array of nuclei:

More information

Convergence: How to obtain reliable results

Convergence: How to obtain reliable results Convergence: How to obtain reliable results Keith Refson Royal Holloway, University of London August 7, 2015 Results of First-Principles Simulations................................................................................

More information

Anion-redox nanolithia cathodes for Li-ion batteries

Anion-redox nanolithia cathodes for Li-ion batteries ARTICLE NUMBER: 16111 Anion-redox nanolithia cathodes for Li-ion batteries Zhi Zhu 1,2, Akihiro Kushima 1,2, Zongyou Yin 1,2, Lu Qi 3 *, Khalil Amine 4, Jun Lu 4 * and Ju Li 1,2 * 1 Department of Nuclear

More information

GIPAW: A solid-state theory for NMR

GIPAW: A solid-state theory for NMR GIPAW: A solid-state theory for NMR Jonathan Yates jonathan.yates@materials.ox.ac.uk Materials Modelling Laboratory, Oxford Materials NMR parameters we will focus on non-metallic, diamagnetic materials

More information

Weile Jia 1, Long Wang 1, Zongyan Cao 1, Jiyun Fu 1, Xuebin Chi 1, Weiguo Gao 2, Lin-Wang Wang 3

Weile Jia 1, Long Wang 1, Zongyan Cao 1, Jiyun Fu 1, Xuebin Chi 1, Weiguo Gao 2, Lin-Wang Wang 3 A plane wave pseudopotential density functional theory molecular dynamics code on multi-gpu machine - GPU Technology Conference, San Jose, May 17th, 2012 Weile Jia 1, Long Wang 1, Zongyan Cao 1, Jiyun

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

College of Science, Xi an University of Science and Technology, Xi an *Corresponding author

College of Science, Xi an University of Science and Technology, Xi an *Corresponding author 2016 International Conference on Advanced Manufacture Technology and Industrial Application (AMTIA 2016) ISBN: 978-1-60595-387-8 The Study of Coordination Adsorption Effect that CO Adsorption on 4H-SiC

More information

Geometry Optimisation

Geometry Optimisation Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background

More information

1 Setting the stage: why ab initio molecular dynamics?

1 Setting the stage: why ab initio molecular dynamics? 1 Setting the stage: why ab initio molecular dynamics? Classical molecular dynamics using predefined potentials, force fields, either based on empirical data or on independent electronic structure 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

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

CHEM6085: Density Functional Theory

CHEM6085: Density Functional Theory Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues

More information

Intro to ab initio methods

Intro to ab initio methods Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio

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

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

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

Electronic Supplementary Information

Electronic Supplementary Information Electronic Supplementary Information Xenon as a mediator of chemical reactions? ( ) by Dominik Kurzydłowski and Wojciech Grochala Contents. 1. Unit cell vectors, fractional atomic coordinates and enthalpy

More information

Theoretical studies of EELS in carbon structures

Theoretical studies of EELS in carbon structures Theoretical studies of EELS in carbon structures Nobuyuki Fujita Submitted for the Degree of Master of Philosophy University of York Department of Physics October 2012 Abstract Graphene nano-ribbons (GNRs)

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

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

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

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

Sub-Angstrom Edge Relaxations Probed by Electron. Microscopy in Hexagonal Boron Nitride (h-bn)

Sub-Angstrom Edge Relaxations Probed by Electron. Microscopy in Hexagonal Boron Nitride (h-bn) Supplementary Material for Sub-Angstrom Edge Relaxations Probed by Electron Microscopy in Hexagonal Boron Nitride (h-bn) Nasim Alem 1,2,3, Quentin M. Ramasse 4,*, Che R. Seabourne 5, Oleg V. Yazyev 1,3,6,

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

Supporting Information

Supporting Information Supporting Information The Origin of Active Oxygen in a Ternary CuO x /Co 3 O 4 -CeO Catalyst for CO Oxidation Zhigang Liu, *, Zili Wu, *, Xihong Peng, ++ Andrew Binder, Songhai Chai, Sheng Dai *,, School

More information

Practical calculations using first-principles QM Convergence, convergence, convergence

Practical calculations using first-principles QM Convergence, convergence, convergence Practical calculations using first-principles QM, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 20, 2012 - CASTEP Workshop: Frankfurt 2012 1 / 21 First-principles

More information

Experiment and Modelling in Structural NMR

Experiment and Modelling in Structural NMR EPJ Web of Conferences 30, 04002 (2012) DOI: 10.1051/epjconf/20123004002 Owned by the authors, published by EDP Sciences, 2012 Experiment and Modelling in Structural NMR November 28th December 2 nd 2011

More information

Supporting Information for. Structural and Chemical Dynamics of Pyridinic Nitrogen. Defects in Graphene

Supporting Information for. Structural and Chemical Dynamics of Pyridinic Nitrogen. Defects in Graphene Supporting Information for Structural and Chemical Dynamics of Pyridinic Nitrogen Defects in Graphene Yung-Chang Lin, 1* Po-Yuan Teng, 2 Chao-Hui Yeh, 2 Masanori Koshino, 1 Po-Wen Chiu, 2 Kazu Suenaga

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

Density functional theory in solution: Implementing an implicit solvent model for CASTEP and ONETEP

Density functional theory in solution: Implementing an implicit solvent model for CASTEP and ONETEP Density functional theory in solution: Implementing an implicit solvent model for CASTEP and ONETEP, November 2016 James C. Womack University of Southampton, Southampton, UK Overview Background The model

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

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

Supporting Information

Supporting Information Supporting Information A Porous Two-Dimensional Monolayer Metal-Organic Framework Material and its Use for the Size-Selective Separation of Nanoparticles Yi Jiang, 1 Gyeong Hee Ryu, 1, 3 Se Hun Joo, 4

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

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

Introduction to first-principles modelling and CASTEP

Introduction to first-principles modelling and CASTEP to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.

More information

Electronic Structure of Crystalline Solids

Electronic Structure of Crystalline Solids Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem

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

Introduction to VASP (1)

Introduction to VASP (1) Introduction to VASP (1) Yun-Wen Chen http://cms.mpi.univie.ac.at/vasp/vasp/vasp.html https://www.vasp.at/ 1 What is VASP? The Vienna Ab initio Simulation Package (VASP) is a computer program for atomic

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

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

Generation of Thermal Scattering Laws for YH 2 using Ab Initio Methods

Generation of Thermal Scattering Laws for YH 2 using Ab Initio Methods Generation of Thermal Scattering Laws for YH 2 using Ab Initio Methods Michael L. Zerkle Bettis Atomic Power Laboratory WPEC SG42 Meeting May 18, 2015 May 18-19, 2015 WPEC SG42 Slide 1 Outline Motivation

More information

arxiv: v1 [cond-mat.mtrl-sci] 4 Mar 2015

arxiv: v1 [cond-mat.mtrl-sci] 4 Mar 2015 Auxiliary density functionals: a new class of methods for efficient, stable density functional theory calculations P.J. Hasnip and M.I.J. Probert Department of Physics, University of York, York YO10 5DD,

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

Vibrational Spectroscopy Methods

Vibrational Spectroscopy Methods Vibrational Spectroscopy Methods Keith Refson STFC Rutherford Appleton Laboratory September 10, 2012 Phonons and Spectroscopy: Workshop: Frankfurt 2012 1 / 39 Motivations from experimental spectroscopy:

More information

Algorithms and Computational Aspects of DFT Calculations

Algorithms and Computational Aspects of DFT Calculations Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational

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

Electronic energy optimisation in ONETEP

Electronic energy optimisation in ONETEP Electronic energy optimisation in ONETEP Chris-Kriton Skylaris cks@soton.ac.uk 1 Outline 1. Kohn-Sham calculations Direct energy minimisation versus density mixing 2. ONETEP scheme: optimise both the density

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

Phonon Transport in Nanostructures

Phonon Transport in Nanostructures Phonon Transport in Nanostructures David-Alexander Robinson & Pádraig Ó Conbhuí 08332461 13th December 2011 Contents 1 Introduction & Theory 1 1.1 TiO 2 Nanostructure Production...................... 1

More information

Supplementary material. From cellulose to kerogen: molecular simulation. of a geological process

Supplementary material. From cellulose to kerogen: molecular simulation. of a geological process Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2017 Supplementary material From cellulose to kerogen: molecular simulation of a geological

More information

Estimation of linear thermal expansion coefficient from cohesive energy obtained by ab-initio calculation of metals and ceramics

Estimation of linear thermal expansion coefficient from cohesive energy obtained by ab-initio calculation of metals and ceramics aaaaa Journal of the Ceramic Society of Japan 118 [3] 241-245 2010 Paper Estimation of linear thermal expansion coefficient from cohesive energy obtained by ab-initio calculation of metals and ceramics

More information

Structure and Formation Mechanism of Black TiO 2 Nanoparticles

Structure and Formation Mechanism of Black TiO 2 Nanoparticles Structure and Formation Mechanism of Black TiO 2 Nanoparticles Mengkun Tian 1, Masoud Mahjouri-Samani 2, Gyula Eres 3*, Ritesh Sachan 3, Mina Yoon 2, Matthew F. Chisholm 3, Kai Wang 2, Alexander A. Puretzky

More information

First-principles studies of the structural, electronic, and optical properties of a novel thorium compound Rb 2 Th 7 Se 15

First-principles studies of the structural, electronic, and optical properties of a novel thorium compound Rb 2 Th 7 Se 15 First-principles studies of the structural, electronic, and optical properties of a novel thorium compound Rb 2 Th 7 Se 15 M.G. Brik 1 Institute of Physics, University of Tartu, Riia 142, Tartu 5114, Estonia

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

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

Modern Theory of Solids

Modern Theory of Solids Quantum Mechanical Approach to the Energy Bandgap Knowlton 1 Quantum Mechanical Approach to the Energy Bandgap a+ b = a o = d-spacing of 1D lattice (or plane in 3D) Knowlton 2 Symmetric vs- Asymmetric

More information

Partition-of-unity finite-element method for large scale quantum molecular dynamics on massively parallel computational platforms

Partition-of-unity finite-element method for large scale quantum molecular dynamics on massively parallel computational platforms LLNL-TR-470692 Partition-of-unity finite-element method for large scale quantum molecular dynamics on massively parallel computational platforms J. E. Pask, N. Sukumar, M. Guney, W. Hu March 1, 2011 Disclaimer

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