Pseudopotentials and Basis Sets. How to generate and test them

Size: px
Start display at page:

Download "Pseudopotentials and Basis Sets. How to generate and test them"

Transcription

1 Pseudopotentials and Basis Sets How to generate and test them

2 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 functions must be core valence orthogonal to the core wave functions

3 Cut-off radii Balance between softness and transferability controlled by R c R c SOFTNESS TRANSFERABILITY Larger R c : softer pseudo Si Shorter R c : harder pseudo

4 The atom program pseudopotential generation label xc flavor # valence shells # core orbitals n l Valence configuration Cutoff radii

5 Run the shell script (pg.sh) Procedure (I) Check contents of new directory (Si.tm2) Plot the pseudo-potentials/orbitals

6 Pseudo-wave function Procedure (II) Logarithmic derivative Radial Fourier-T Pseudopotential

7 Procedure (III)

8 Radial charge distribution: Compare all-electron with pseudo-charge

9 Pseudopotential testing (I) The all-electron (ae.sh) and pseudo-test (pt.sh) scripts: run script input pseudopotential to test

10 Pseudopotential testing (II) A transferable pseudo will reproduce the AE energy levels and wave functions in arbitrary environments

11 SUMMER SCHOOL ON COMPUTATIONAL MATERIALS SCIENCE Pseudopotential testing (III) Compute the energy of two different configurations Compute the difference in energy For the pseudopotential to be transferible: 3s2 3p2 (reference) 3s2 3p1 3d1 3s1 3p3 3s1 3p2 3d1 3s0 3p3 3d1

12 Large core-valence overlap Errors due to non-linearity of XC-potential

13 ~ V Non-linear core corrections Standard pseudopotential unscreening: Valence charge only V ps However xc V xc = V xc = V [ ρ ] + v ps scr [ ρ ]( r) V [ ρ ]( r) V [ ρ ]( r) v { V [ ρ + ρ ] V [ ρ ]} xc v c H [ ρ + ρ ]( r) V [ ρ ]( r) + V [ ρ ]( r) v c xc Core-correction v v xc v xc xc v c Keep core charge in pseudopotential generation V ps ps ( r) = V [ ρ + ρ ]( r) V [ ρ ]( r) V [ ρ + ρ ]( r) scr v c H v xc v c

14 New flag PCC input file

15 Pseudo-core & pseudo-valence charge

16 Smooth Fourier Transform The real-space grid required fineness depends on how you define the pseudopotential. The meshcutoff parameter can be determined from the Fourier Transform.

17 Basis generation

18 Key for f linear-scaling scaling: : LOCALITY Large system x2 Divide and Conquer W. Yang, Phys. Rev. Lett. 66, 1438 (1992)

19 Very efficient Atomic Orbitals Lack of systematic for convergence Main features: φ Ilm ( r )=R Il (r I )Y lm ( ˆ r I ) r = r R r I I Numerical Atomic Orbitals (NAOs): Numerical solution of the KS Hamiltonian for the isolated pseudoatom with the same approximations (xc, pseudos) as for the condensed system Size or number of functions Range of localization of these functions Shape or functional form used.

20 Basis Size (I) Depends on the required accuracy and available computational power Quick and dirty calculations Minimal basis set (single-ζ; SZ) Highly converged calculations Complete multiple-ζ + Polarization + Diffuse orbitals

21 Basis Size (II): improving One single radial function per angular momentum shell occupied in the free-atom. Single-ζ (minimal or SZ) Improving the quality? Radial flexibilization: Add more than one radial function within the same angular momentum shell Multiple-ζ Angular flexibilization: Add shells of different angular momentum Polarization

22 Examples

23 Basis Size (III): Polarization Perturbative polarization Apply a small E field to the orbital we want to polarize E Atomic polarization Solve Schrödinger equation for higher angular momentum unbound in the free atom require short cut offs s s+p Si 3d orbitals E. Artacho et al, Phys. Stat. Sol. (b), 215, 809 (1999)

24 Basis Size (IV): Convergence Cohesion curves Bulk Si PW and NAO convergence J. Junquera et. al. PRB 64, (2001).

25 Range (I): How to get sparsity for O(n) Neglecting interactions below a tolerance or beyond some scope of neighbours numerical instablilities for high tolerances. Strictly localized atomic orbitals (zero beyond a given cutoff radius, r c ) Accuracy and computational efficiency depend on the range of the atomic orbitals. Way to define all the cutoff radii in a balanced way.

26 Range (II): Energy Shift Easy approach to define the cutoff radii for the NAOs: 1 2r d dr 2 2 r l( l + 1) + 2 2r + V l ( r) φ l ( r) = ( ε l + δε l ) φl ( r) E. Artacho et al. Phys. Stat. Solidi (b) 215, 809 (1999) A single parameter for all cutoff radii Fireballs O. F. Sankey & D. J. Niklewski, Phys. Rev. B 40, 3979 (1989) BUT, a different cutoff radius for each orbital

27 Range (II): Convergence Bulk Si equal s, p orbitals radii J. Soler et al, J. Phys: Condens. Matter, 14, 2745 (2002)

28 Shape (I) The radial function shape is mainly determined by the pseudopotential. Extra parameters can be introduced to add flexibility: δq : extra charge per atomic specie. Confinement : imposed separately for each angular momentum shell.

29 Shape (II) Soft confinement (J. Junquera et al, Phys. Rev. B 64, (01) ) Shape of the optimal 3s orbital of Mg in MgO for different schemes Corresponding optimal confinement potential Better variational basis sets Removes the discontinuity of the derivative

30 PAO.Basis (I) # shells Add polarization 4s 3d 1 st ζ %block PAO.Basis # Define Basis set Cu 2 # Species label, number of l-shells n=4 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol Rc n=3 2 2 # n, l, Nzeta Rc nd ζ %endblock PAO.Basis

31 PAO.Basis (II): new generation charge %block PAO.Basis # Define Basis set Cu n=4 0 1 E V c r c n=4 1 1 E Polarization orbital n=3 2 1 E %endblock PAO.Basis

32 Procedure 1. Check the difference in energies involved in your problem 2. For semiquantitative results and general trends use SZ 3. Improve the basis: Automatic DZP (Split Valence & Perturbative Polarization): High quality for most systems Good valence between well converged results & computational cost Standard Rule of thumb in Quantum Chemistry: «a basis should always be doubled before being polarized». 4. Functional optimization of the basis

33 Pseudos & Basis repository Pseudopotentials and basis sets available in the SIESTA web page: Uploaded by users input files to generate them Plots of the radial functions Documentation of the tests done Author s contact information The PAO is pseudopotential-dependent. Check also in the user s mailing list.

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

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

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

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

Systematic convergence for realistic projects Fast versus accurate

Systematic convergence for realistic projects Fast versus accurate Systematic convergence for realistic projects Fast versus accurate Daniel Sánchez-Portal Centro de Física de Materiales, Centro Mixto CSIC- UPV/EHU,San Sebastián, Spain Email: sqbsapod@sc.ehu.es Thanks

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

The Siesta program for electronic structure simulations

The Siesta program for electronic structure simulations Université de Liège The Siesta program for electronic structure simulations Javier Junquera Départament de Physique Université de Liège Linear Scaling CPU load ~ N 3 Early 90 s ~ N ~ 100 N (# atoms) G.

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

2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density

2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density 1. SIESTA DFT In a Nutshell Introduction to SIESTA Atomic Orbitals Capabilities Resources 2. TranSIESTA Transport in the Nanoscale - motivation Boundary Conditions: Open systems Greens functions and charge

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

LDA+U. How to run a simulation within LDA+U framework in SIESTA. Javier Junquera. Daniel Sánchez-Portal

LDA+U. How to run a simulation within LDA+U framework in SIESTA. Javier Junquera. Daniel Sánchez-Portal LDA+U How to run a simulation within LDA+U framework in SIESTA Javier Junquera Daniel Sánchez-Portal Transition metal monoxides are prototypes of highly correlated materials Paradigmatic example: NiO Above

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

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

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

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

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

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

DFT / SIESTA algorithms

DFT / SIESTA algorithms DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang

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

How to run SIESTA. Introduction to input & output files

How to run SIESTA. Introduction to input & output files How to run SIESTA Introduction to input & output files Linear-scaling DFT based on Numerical Atomic Orbitals (NAOs) Born-Oppenheimer DFT Pseudopotentials Numerical atomic orbitals relaxations, MD, phonons.

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

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

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

References. Documentation Manuals Tutorials Publications

References.   Documentation Manuals Tutorials Publications References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =

More information

June 17, Curtin University of Technology, Perth

June 17, Curtin University of Technology, Perth U s e r s G u i d e S I E S T A 2.0.2-c11 June 17, 2007 Emilio Artacho Julian D. Gale Alberto García Javier Junquera Richard M. Martin Pablo Ordejón Daniel Sánchez-Portal José M. Soler University of Cambridge

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

Is Band better than ADF? The core of the issue

Is Band better than ADF? The core of the issue Is Band better than ADF? The core of the issue Introduction Introduction What is Band? Introduction What is Band? Slightly different basis set for ADF/Band Introduction What is Band? Slightly different

More information

Basic introduction to SIESTA

Basic introduction to SIESTA Basic introduction to SIESTA Emilio Artacho Department of Earth Sciences University of Cambridge 1 The method Linear-scaling DFT based on NAOs (Numerical Atomic Orbitals) P. Ordejon, E. Artacho & J. M.

More information

Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method

Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2016.11.11 Recall:

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

CHAPTER 3 WIEN2k. Chapter 3 : WIEN2k 50

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

More information

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

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 (DFT) modelling of C60 and

Density Functional Theory (DFT) modelling of C60 and ISPUB.COM The Internet Journal of Nanotechnology Volume 3 Number 1 Density Functional Theory (DFT) modelling of C60 and N@C60 N Kuganathan Citation N Kuganathan. Density Functional Theory (DFT) modelling

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

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

The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR)

The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR) The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR) 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) Please

More information

Basic introduction to running Siesta. Eduardo Anglada Siesta foundation-uam-nanotec

Basic introduction to running Siesta. Eduardo Anglada Siesta foundation-uam-nanotec Basic introduction to running Siesta Eduardo Anglada Siesta foundation-uam-nanotec eduardo.anglada@uam.es eduardo.anglada@nanotec.es Our method Linear-scaling DFT based on NAOs (Numerical Atomic Orbitals)

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

The LDA+U method: a primer and implementation within SIESTA

The LDA+U method: a primer and implementation within SIESTA The LDA+U method: a primer and implementation within SIESTA Daniel Sánchez-Portal Thanks to Javier Junquera, Sampsa Riikonen and Eduardo Anglada Source of the failure of LDA to describe Mott insulators

More information

Key word: DensityFunctional Theory, Endohedral Energy gap, Electonic properties.

Key word: DensityFunctional Theory, Endohedral Energy gap, Electonic properties. First Principle Studies of Electronic Properties of Nitrogen-doped Endohedral Fullerene 1 M. R. Benam, 2 N. Shahtahmasbi, 1 H. Arabshahi and 1 Z.Zarei 1 Department of Physics, Payame Noor University, P.

More information

2.4. Quantum Mechanical description of hydrogen atom

2.4. Quantum Mechanical description of hydrogen atom 2.4. Quantum Mechanical description of hydrogen atom Atomic units Quantity Atomic unit SI Conversion Ang. mom. h [J s] h = 1, 05459 10 34 Js Mass m e [kg] m e = 9, 1094 10 31 kg Charge e [C] e = 1, 6022

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

Basis Sets and Basis Set Notation

Basis Sets and Basis Set Notation Chemistry 46 Fall 215 Dr. Jean M. Standard November 29, 217 Basis Sets and Basis Set Notation Using the LCAO-MO approximation, molecular orbitals can be represented as linear combinations of atomic orbitals,

More information

Computational Modeling Module: Electronic Band Structure Calculations MSE 410, 10/03 & 10/5

Computational Modeling Module: Electronic Band Structure Calculations MSE 410, 10/03 & 10/5 Computational Modeling Module: Electronic Band Structure Calculations MSE 410, 10/03 & 10/5 Learning Outcomes and Topics Learning outcomes: Get familiar with density functional theory (DFT) and simulation

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

Density-functional calculations of defect formation energies using the supercell method: Brillouin-zone sampling

Density-functional calculations of defect formation energies using the supercell method: Brillouin-zone sampling Density-functional calculations of defect formation energies using the supercell method: Brillouin-zone sampling Jihye Shim and Eok-Kyun Lee Department of Chemistry and School of Molecular Science (BK21),

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

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 8 Nov 2001

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 8 Nov 2001 The Siesta method for ab initio order-n materials simulation arxiv:cond-mat/111138v1 [cond-mat.mtrl-sci] 8 Nov 1 José M. Soler, 1 Emilio Artacho, Julian D. Gale, 3 Alberto García, 4 Javier Junquera, 1,5

More information

Explaining the apparent arbitrariness of the LDA-1/2 self-energy. correction method applied to purely covalent systems

Explaining the apparent arbitrariness of the LDA-1/2 self-energy. correction method applied to purely covalent systems Explaining the apparent arbitrariness of the LDA-1/2 self-energy correction method applied to purely covalent systems Kan-Hao Xue, 1,2 Leonardo R. C. Fonseca, 3 and Xiang-Shui Miao 1,2 1 School of Optical

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

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

ECE440 Nanoelectronics. Lecture 07 Atomic Orbitals

ECE440 Nanoelectronics. Lecture 07 Atomic Orbitals ECE44 Nanoelectronics Lecture 7 Atomic Orbitals Atoms and atomic orbitals It is instructive to compare the simple model of a spherically symmetrical potential for r R V ( r) for r R and the simplest hydrogen

More information

How to compute the projected density of states (PDOS) Javier Junquera

How to compute the projected density of states (PDOS) Javier Junquera How to compute the projected density of states (PDOS) Javier Junquera Density Of States (DOS) the number of one-electron levels between E and E + de SrTiO 3 bulk Units: (Energy) -1 Projected Density Of

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

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

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

More information

Gaussian Basis Sets for Solid-State Calculations

Gaussian Basis Sets for Solid-State Calculations Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction

More information

Efficient projector expansion for the ab initio LCAO method

Efficient projector expansion for the ab initio LCAO method PHYSICAL REVIEW B 72, 045121 2005 Efficient projector expansion for the ab initio LCAO method T. Ozaki Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science

More information

Ab initio study of Mn doped BN nanosheets Tudor Luca Mitran

Ab initio study of Mn doped BN nanosheets Tudor Luca Mitran Ab initio study of Mn doped BN nanosheets Tudor Luca Mitran MDEO Research Center University of Bucharest, Faculty of Physics, Bucharest-Magurele, Romania Oldenburg 20.04.2012 Table of contents 1. Density

More information

Final Exam Tuesday, May 8, 2012 Starting at 8:30 a.m., Hoyt Hall Duration: 2h 30m

Final Exam Tuesday, May 8, 2012 Starting at 8:30 a.m., Hoyt Hall Duration: 2h 30m Final Exam Tuesday, May 8, 2012 Starting at 8:30 a.m., Hoyt Hall. ------------------- Duration: 2h 30m Chapter 39 Quantum Mechanics of Atoms Units of Chapter 39 39-1 Quantum-Mechanical View of Atoms 39-2

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

S I E S T A 4.0. June 23, CIC-Nanogune and University of Cambridge

S I E S T A 4.0. June 23, CIC-Nanogune and University of Cambridge U s e r s G u i d e S I E S T A 4.0 June 23, 2016 Emilio Artacho José María Cela Julian D. Gale Alberto García Javier Junquera Richard M. Martin Pablo Ordejón Daniel Sánchez-Portal José M. Soler CIC-Nanogune

More information

AB INITIO STUDY OF NANO STRUCTURED FUNCTIONALIZED GRAPHENE WITH 30C ATOMS

AB INITIO STUDY OF NANO STRUCTURED FUNCTIONALIZED GRAPHENE WITH 30C ATOMS International Journal of Science, Environment and Technology, Vol. 1, No 3, 2012, 108-112 AB INITIO STUDY OF NANO STRUCTURED FUNCTIONALIZED GRAPHENE WITH 30C ATOMS Naveen Kumar* and Jyoti Dhar Sharma Deptt.

More information

Introduction to Heisenberg model. Javier Junquera

Introduction to Heisenberg model. Javier Junquera Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics

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

Atomic Structure, Periodic Table, and Other Effects: Chapter 8 of Rex and T. Modern Physics

Atomic Structure, Periodic Table, and Other Effects: Chapter 8 of Rex and T. Modern Physics Atomic Structure, Periodic Table, and Other Effects: Chapter 8 of Rex and T Modern Physics 11/16 and 11/19/2018 1 Introduction In Chapter 7, we studied the hydrogen atom. What about other elements, e.g.,

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

An Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory

An Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement

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

AP Chemistry - Problem Drill 10: Atomic Structures and Electron Configuration

AP Chemistry - Problem Drill 10: Atomic Structures and Electron Configuration AP Chemistry - Problem Drill 10: Atomic Structures and Electron Configuration No. 1 of 10 Instructions: (1) Read the problem statement and answer choices carefully (2) Work the problems on paper as 1.

More information

Tight-Binding Model of Electronic Structures

Tight-Binding Model of Electronic Structures Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to

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

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components. Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative

More information

Physics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =...

Physics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =... Physics 606, Quantum Mechanics, Final Exam NAME Please show all your work. (You are graded on your work, with partial credit where it is deserved.) All problems are, of course, nonrelativistic. 1. Consider

More information

LS coupling. 2 2 n + H s o + H h f + H B. (1) 2m

LS coupling. 2 2 n + H s o + H h f + H B. (1) 2m LS coupling 1 The big picture We start from the Hamiltonian of an atomic system: H = [ ] 2 2 n Ze2 1 + 1 e 2 1 + H s o + H h f + H B. (1) 2m n e 4πɛ 0 r n 2 4πɛ 0 r nm n,m Here n runs pver the electrons,

More information

Conquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale

Conquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale Fortran Expo: 15 Jun 2012 Conquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale Lianheng Tong Overview Overview of Conquest project Brief Introduction

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

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

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

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

Localized basis methods Theory and implementations. Introduction of OpenMX Implementation of OpenMX

Localized basis methods Theory and implementations. Introduction of OpenMX Implementation of OpenMX Localized basis methods Theory and implementations Introduction of OpenMX Implementation of OpenMX Total energy Pseudopontials Basis functions Self-consistency Δ-gauge Taisuke Ozaki (ISSP, Univ. of Tokyo)

More information

arxiv:cond-mat/ v1 17 May 1995

arxiv:cond-mat/ v1 17 May 1995 Projection of plane-wave calculations into atomic orbitals Daniel Sanchez-Portal, Emilio Artacho, and Jose M. Soler Instituto de Ciencia de Materiales Nicolás Cabrera and Departamento de Física de la Materia

More information

Atomic and molecular physics Revision lecture

Atomic and molecular physics Revision lecture Atomic and molecular physics Revision lecture Answer all questions Angular momentum J`2 ` J z j,m = j j+1 j,m j,m =m j,m Allowed values of mgo from j to +jin integer steps If there is no external field,

More information

The 3 dimensional Schrödinger Equation

The 3 dimensional Schrödinger Equation Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum

More information

COMPUTATIONAL TOOL. Fig. 4.1 Opening screen of w2web

COMPUTATIONAL TOOL. Fig. 4.1 Opening screen of w2web CHAPTER -4 COMPUTATIONAL TOOL Ph.D. Thesis: J. Maibam CHAPTER: 4 4.1 The WIEN2k code In this work, all the calculations presented are performed using the WIEN2k software package (Blaha et al., 2001). The

More information

OPENATOM for GW calculations

OPENATOM for GW calculations OPENATOM for GW calculations by OPENATOM developers 1 Introduction The GW method is one of the most accurate ab initio methods for the prediction of electronic band structures. Despite its power, the GW

More information

arxiv: v1 [cond-mat.mes-hall] 13 Feb 2012

arxiv: v1 [cond-mat.mes-hall] 13 Feb 2012 Controlling Band Gap in Silicene Monolayer Using External Electric Field C. Kamal 1 arxiv:1202.2636v1 [cond-mat.mes-hall] 13 Feb 2012 1 Indus Synchrotrons Utilization Division, Raja Ramanna Centre for

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

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

Chapter 5. Atomic spectra

Chapter 5. Atomic spectra Atomic spectra Sommerfelds relativistic model Sommerfeld succeeded partially in explaining fine structure by extending Bohr Theory i) He allowed the possibility of elliptical orbits for the electrons in

More information

Goal: find Lorentz-violating corrections to the spectrum of hydrogen including nonminimal effects

Goal: find Lorentz-violating corrections to the spectrum of hydrogen including nonminimal effects Goal: find Lorentz-violating corrections to the spectrum of hydrogen including nonminimal effects Method: Rayleigh-Schrödinger Perturbation Theory Step 1: Find the eigenvectors ψ n and eigenvalues ε n

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

This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.

This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. 1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis

More information

Section 3 Electronic Configurations, Term Symbols, and States

Section 3 Electronic Configurations, Term Symbols, and States Section 3 Electronic Configurations, Term Symbols, and States Introductory Remarks- The Orbital, Configuration, and State Pictures of Electronic Structure One of the goals of quantum chemistry is to allow

More information

Phys 622 Problems Chapter 5

Phys 622 Problems Chapter 5 1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit

More information

U s e r s G u i d e S I E S T A 3.0-rc2 July 5, 2010 Emilio Artacho José María Cela Julian D. Gale Alberto García Javier Junquera Richard M. Martin Pablo Ordejón Daniel Sánchez-Portal José M. Soler University

More information

Hund s rule for monopole harmonics, or why the composite fermion picture works

Hund s rule for monopole harmonics, or why the composite fermion picture works PERGAMON Solid State Communications 110 (1999) 45 49 Hund s rule for monopole harmonics, or why the composite fermion picture works Arkadiusz Wójs*, John J. Quinn The University of Tennessee, Knoxville,

More information

Kondo effect in multi-level and multi-valley quantum dots. Mikio Eto Faculty of Science and Technology, Keio University, Japan

Kondo effect in multi-level and multi-valley quantum dots. Mikio Eto Faculty of Science and Technology, Keio University, Japan Kondo effect in multi-level and multi-valley quantum dots Mikio Eto Faculty of Science and Technology, Keio University, Japan Outline 1. Introduction: next three slides for quantum dots 2. Kondo effect

More information

LABELING ELECTRONS IN ATOMS

LABELING ELECTRONS IN ATOMS Date: Name: LABELING ELECTRONS IN ATOMS The location of each electron in an atom is determined by a few different factors. Each factor is represented by a QUANTUM NUMBER. Prediction: What do you think

More information

8.1 The hydrogen atom solutions

8.1 The hydrogen atom solutions 8.1 The hydrogen atom solutions Slides: Video 8.1.1 Separating for the radial equation Text reference: Quantum Mechanics for Scientists and Engineers Section 10.4 (up to Solution of the hydrogen radial

More information