DFT in practice : Part II. Ersen Mete
|
|
- Eileen Norman
- 5 years ago
- Views:
Transcription
1 pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, NanoDFT 09, İzmir Institute of Technology, İzmir
2 Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials Nonlocal separable Kleinman-Bylander form Ultrasoft pseudopotentials Projector Augmented Waves (PAW) method
3 A one-electron atom For a one-electron atom, the attractive potential is spherically symmetric, V ( r) = V (r) = Z r then, the solutions are separable to radial and angular parts, ψ nlm ( r) = R nl (r)yl m (θ, ϕ) = φ nl(r) Yl m (θ, ϕ) r The radial equation becomes, 1 d 2 [ l(l + 1) 2 dr 2 φ nl + 2r 2 ] + V eff (r) φ nl = εφ nl Kohn-Sham single particle Schrödinger-like equations will be identical if V KS ( r) is spherically symmetric as one-electron Coulomb potential V (r).
4 V KS ( r) = V KS (r)? Pseudopotentials density : n( r) = occ n,l,m occ ψ nlm ( r) 2 = (2l + 1) R n,l (r) 2 = n(r) external potential : υ ext ( r) = Z/r = υ ext (r) n( r ) Hartree potential : r r d r = V H (r) exchange-corr. pot. :V xc ( r) = ɛ xc [n(r)] + n(r) dɛ xc dn [n(r)] =V xc(r) Therefore, total effective potential n,l V KS = υ ext (r) + V H (r) + V xc (r) spherically symmetric. The independent-particle Kohn-Sham equations are analogous to the Schrödinger equation of the one-electron atom, 1 d 2 [ ] l(l + 1) 2 dr 2 φ nl + 2r 2 + V KS (r) φ nl = εφ nl DFT for an atom
5 Why do we need pseudopotentials? DFT calculation with all-electron υ ext is expensive & Core electrons are essentially inert in bonding environments. Pseudopotentials, replace the effect of the core electrons, are smooth in the core region, reproduces all-electron potential behavior out of the core region. Computationally, Sharp oscillations near the core region will be smoothed : Reduction of the number of plane waves The number of electrons will be decreased : Reduction of the number of bands to solve for
6 Pseudopotentials Schematically
7 Pseudopotential terminology transferability : ability to describe the valence electrons in different environments. softness : the need for the number of plane waves. inclusion of semicore states Example : Ti with large core 1s 2 2s 2 2p 6 3s 2 3p 6 4s }{{} 2 3d 2 core Ti with semicore 1s 2 2s 2 2p 6 3s }{{} 2 3p 6 4s 2 3d 2 core locality : all l-channel (s, p, d) electrons feel the same potential efficiency a compromise between accuracy and computational cost
8 Accuracy vs computational load Local PSP V PS = V PS (r) Semilocal PSP V PS = l V PS l (r) χ l χ l Nonlocal separable PSP V PS = Vloc PS (r) + D l β lm β lm lm
9 Total energy in terms of valence electron density Then, N core E[{ψ i }] = ψ i ψ i + i n( r) = n core ( r) + n val ( r) υ ext ( r)n core ( r)d r + 1 ncore ( r)n core ( r ) 2 r r d rd r N val + ψ i ψ i + υ ext ( r)n val ( r)d r + 1 nval ( r)n val ( r ) 2 r r d rd r i ncore ( r)n val ( r ) + r r d rd r + E xc [n core + n val ] N val E val [{ψ i }] = ψ i - 1 ( 2 2 ψ i + - Z r + ncore ( r ) ) r r n val ( r)d r+e xc [n core + n val ] }{{} i }{{} Non linear XC Vion SCR corrections
10 Generic pseudopotential transformation For an atom, assume that the core states, χ n, satisfy H χ n = E n χ n A single valence state, ψ, can be replaced by a smoother pseudofunction, φ, expanding the remaining part in terms of χ n, core ψ = φ + a n χ n Using the orhtogonality of valence and core states, core χ m ψ = χ m φ + a n χ m χ n = 0 = a m = χ m φ n n ψ = φ n χ n φ χ n
11 Then, the eigenvalue equation, ( H φ core ) χ n φ χ n n ( = E φ core H φ + (E E n ) χ n χ n φ = E φ This implies, n (H + V nl ) φ = E φ where l is due to spherical symmetry. E E n > 0 extra potential V nl is repulsive. Cancels the effect of the attractive Coulomb potential. core ) χ n φ χ n Resulting potential is weaker and pseudo eigenstate is smoother. n
12 Norm-conserving pseudopotentials Hamann, Schlüter and Chiang [Phys.Rev.Lett.43,1494(1979)] criteria : AE & PS wavefunctions correspond to the same energy for the reference level (AE valence level with angular momentum l), H φ AE nl = ε nl φ AE nl (H + V nl ) φ PS nl = ε nl φ PS nl AE & PS wavefunctions match beyond a certain radial cutoff, r c, φ AE nl (r) = φ PS nl (r) r r c
13 AE and PS norm squares integrated upto r r c are equal. r 0 φ AE nl (r ) 2 dr = r 0 φ PS nl (r ) 2 dr Equal amount of charge in the core region. Gauss Law is satisfied for r. Normalization constraint is achieved in the limit r. logarithmic derivatives (their respective potentials) agree for r r c. d dr ln[φps nl (r)] d dr ln[φae nl (r)] V PS must reproduce the same scattering phase shifts as V AE for r. Necessary to improve transferability. PS wavefunction is nodeless. It s twice differentiable and satisfies lim r 0 φ nl (r) r l+1 = continuous.
14 Silicon : wave functions
15 Norm-conserving PSP generation steps 1 Solve the all-electron atomic system. 2 Determine the core and valence states. 3 Apply norm-conservation criteria (e.g. Hamann scheme) and derive a PS wavefunction from the reference AE valence level with angular momentum l. 4 Invert the Schrödinger equation for PS wavefunctions to get screened PSP components. V SCR,PS l = ε PS 2 l(l + 1) 1 d l 2r 2 + 2φ PS nl (r) dr 2 φps l 5 Subtract the Hartree and XC contributions to obtain υ ext (unscreened PSP). V PS l = V SCR,PS l V H (r) V XC (r) Model pseuodopotential replaces the potential of the nucleus. (r)
16 Semilocal pseudopotentials For l-dependent model PSP, treat angular momentum l separately (nonlocal). PSP in the semilocal form : Vnl PS = V loc (r) + V SL where V SL = Y lm V nl (r) Y lm lm It s local in r and nonlocal in θ, ϕ. Drawback : Plane wave representation of the non-local part is expensive.
17 Semilocal PSP matrix elements in plane waves q V SL q = 1 Ω l,m where r = (r, θ, ϕ ) e i q r Y lm(θ, ϕ)v nl (r)y lm (θ, ϕ)e i q r r 2 dωdω dr e i q r = 4π lm i l j l (qr)y lm(ˆq)y lm (ˆr) where j l (qr) are spherical Bessel functions and ˆq and ˆr denote angles associated with the vectors q and r, respectively. m Y lm (ˆq)Y lm(ˆq ) = 2l + 1 4π P l(cos θ q q ) q V SL q = 4π Ω (2l + 1) l j l (qr)j l (q r)p l (cos θ q q )V l (r)r 2 dr N 2 PW such integrations needed!
18 Fully separable Kleinman-Bylander form where V PS nl ( r, r ) = V loc (r)δ( r r ) + V NL V NL = l V NL l = lm V SL l Action of this form on the PS wavefunction, Vl NL φ PS lm = V SL l φ PS lm φps lm V l SL φ PS lm V l φ PS lm φ PS lm φps lm V l SL φ PS lm V l φ PS lm φ SL lm = Vl SL φ PS lm Planewave representation of non-local PSP matrix elements in Kleinman-Bylander form q V NL q = lm q V SL l φ PS lm φ PS lm V SL l V SL l φ PS φ PS lm lm q Number of integrals to be evaluated reduces to N PW. [PRL 48,1425 (1982)]
19 Silicon & Titanium NCPPs
20 Ti log derivatives Pseudopotentials
21 Ultrasoft pseudopotentials (USPP) : formalism Vanderbilt [PRB 41, 7892 (1990)] proposed a new method by relaxing norm-conservation constraint, φ PS i φ PS i ψi AE ψi AE USPPs are norm-conserving in a generalized form φ PS i (1 + ˆN NL ) φ PS i = ψi AE ψi AE where ˆN is nonlocal charge augmentation operator. USPPs require much smaller PW cutoff (less N PW ) ultrasoft. Scattering properties remain to be correct transferable.
22 The aim is to minimize the total energy, E e = φ i ˆV NL φ i + d rv loc ( r)n( r)+ 1 2 i d rd r n( r)n( r ) r r +E xc[n] subject to where φ i (1 + ˆN NL ) φ j = φ i ŜNL φ j = δ ij n( r) = i ( ) φ i r r + ˆK NL ( r) φ i and for consistency ˆN NL = d r ˆK NL ( r) so that n( r)d r = N v Then, the eigenvalue equation, (T + Vloc PS PS + ˆV NL ) φ i = ε n (1 + ˆN NL ) φ i
23 Ultrasoft PSP generation steps Screened V AE is obtained through self-consistent solution of atomic Kohn-Sham system. Cutoff radii are chosen : r cl for the wave functions r loc c for the local PSP R large enough that all PS and AE quantities agree. A smooth local potential, V loc ( r) is generated which approaches V AE ( r) beyond r loc c.
24 Then, for each angular momentum channel, a few reference energy values, ε i, are chosen where i = {τlm} and τ is the number of reference energies. (T + V AE ε i ) ψ i = 0 ψ i not determined self-consistently. So, ψ i ψ j R = R ψ i ( r)ψ j( r)d r New orbitals defined, χ i = (ε i T V loc ) φ i which vanish at and beyond R where V loc = V AE and φ i = ψ i. To define a nonlocal PSP, χ i are used as projectors to define new wave functions, β i = j χ i χ j φ i = j (B 1 ) ji χ j
25 To compensate valence charge deficit, generalized augmentation charges needed, q ij = ψ i ψ j R φ i φ j R The nonlocal overlap operator can be defined as, S = 1 + i,j q ij β i β j Then the nonlocal potential operator is V NL = ij D ij β i β j where D ij = B ij + ε j q ij PS wave functions satisfy generalized orthonormality condition, φ i S φ j R = φ i φ j R + q ij = ψ i ψ j R = δ ij Then, PS wave functions satisfy generalized eigenvalue problem, (H ε i S) φ i = 0 where H = T + V loc + V NL
26 Verification of generalized eigenvalue problem ( T + V loc + ) ( D nm β n β m φ i = ε i 1 + ) q nm β n β m φ i nm nm where D nm = B nm + ε m q nm Since, ( ) B nm β n β m φ i nm = φ n χ m nm k χ k χ k φ n l φ m χ l χ l φ i = km χ k n χ k φ n φ n χ m δ im = km χ k δ km δ im = χ i Substitution yields, H φ i = ε i S φ i
27 Valence electron density The electron density is augmented, n v ( r) = [ φ n ( r) 2 + n i,j ] Q ij ( r) φ n β i β j φ n where Q ij ( r) = ψ i ( r)ψ j ( r) φ i ( r)φ j ( r) so that it must integrate to the correct number of valence electrons, n v ( r)d r = φ n ( r) 2 d r + Q ij ( r) φ n β i β j φ n d r n ij = nm δ nm [ φ n φ m + ij ] q ij φ n β i β j φ m = nm δ nm φ n S φ m = nm δ nm δ nm = n δ nn = N valence
28 Minimization of total energy δn( r ) δφ n( r) = φ n( r)δ( r r ) + ij Q ij ( r )β i ( r) β j φ n Then, the modified Kohn-Sham equations, δe e δφ n = = d r δe e δn( r ) δn( r ) δφ n( r) [ V eff + ij ( D (0) ij + V eff ( r )Q ij ( r )d r ) ] β i β j φ n where V eff = V H + V loc + V xc The coefficients in the non-local part of the PSP gets updated self-consistently.
29 Ionic USPP Pseudopotentials Ionic potential are obtained by unscreening, V ion loc = V loc V H V xc D (0) ij = D ij d r V loc ( r )n( r )
30
31 Projector Augmented Waves (PAW) method : basic idea For a particular reference energy, P.E. Blöchl, PRB (1994) the behavior of an arbitrary PS wavefunction ψ PS at the atomic site can be calculated by projection at that site in terms of partial waves (spherical), c lm =< P lm ψ PS In each sphere, ψ PS = lm c lm φ PS lm and ψ AE = c lm φ AE lm lm where φ lm are partial waves. Projectors must be dual to partial waves, P lm φ PS l m = δ ll δ mm ψps = i φ i p m ψ PS = ψ PS
32 Then, the PS transformation is given by, ψn AE = ψn PS φ PS lmε p lmε ψn PS lmε + lmε φ AE lmε p lmε ψ PS n Transformation involves AE wavefunction One can derive AE results.
33 References Pseudopotentials E. Kaxiras, Atomic and electronic structure of solids, Cambridge University Press, Cambridge, R.M. Martin, Electronic Structure : Basic Theory and Methods, Cambridge University Press, Cambridge, D.R.Hamann, M. Schlüter, C. Chiang, Phys. Rev. Lett. 43, 1494 (1990). G.B. Bachelet and M. Schlüter, Phys. Rev. B 25, 2103 (1982). L. Kleinman and D.M. Bylander, Phys. Rev. Lett. 48, 1425 (1982). G.B. Bachelet, D.R.Hamann, M. Schlüter, Phys. Rev. B 26, 4199 (1982). A.M. Rappe, K.M. Rabe, E. Kaxiras, and J.D. Joannopoulos, Phys. Rev. B 41, 1227 (1990). N. Troullier and J.L. Martins, Phys. Rev. B 43, 1993 (1991). D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). P.E. Blöchl, Phys. Rev. B 50, (1994).
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 informationWhy 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 informationTheory 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 informationThe 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 informationPseudopotential 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 informationPseudopotentials 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 informationDesigned 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 informationPlane 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 informationThe 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 informationNorm-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 informationPSEUDOPOTENTIALS 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 informationDensity 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 informationDensity 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 informationAn 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 informationIntroduction 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 informationProjector 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 informationStrategies 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 informationPseudopotentials: 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 informationOpen-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 informationTwo 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 informationLecture 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 informationElectron 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 information1. 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 informationProjector-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 informationElectronic 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 informationKey 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 informationThe 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 informationSolving Many-Body Schrödinger Equation Using Density Functional Theory and Finite Elements
Solving Many-Body Schrödinger Equation Using Density Functional Theory and Finite Elements Institute of Physics, Academy of Sciences of the Czech Republic June 21, 2008 Introduction Contens Density Functional
More informationPseudopotentials 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 information6 A Tutorial on Density Functional Theory
6 A Tutorial on Density Functional Theory Fernando Nogueira, Alberto Castro, and Miguel A.L. Marques Departamento de Física, Universidade de Coimbra, Rua Larga, 3004 516, Coimbra, Portugal fnog@teor.fis.uc.pt
More informationConstruction 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 informationAll electron optimized effective potential method for solids
All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22
More informationNotes 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 informationGW 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 informationChapter 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 information1 r 2 sin 2 θ. This must be the case as we can see by the following argument + L2
PHYS 4 3. The momentum operator in three dimensions is p = i Therefore the momentum-squared operator is [ p 2 = 2 2 = 2 r 2 ) + r 2 r r r 2 sin θ We notice that this can be written as sin θ ) + θ θ r 2
More informationElectrons 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 informationNotes 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 informationarxiv: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 informationData 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 informationDensity 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 informationAtomic 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 informationBasics of Density Functional Theory (DFT)
Basics of Density Functional Theory (DFT) Ari Paavo SEITSONEN Ari.P.Seitsonen@iki.fi Département de Chimie École Normale Supérieure, Paris École de Sidi-BelAbbès de Nanomateriaux // Octobre 8-12, 2016
More informationCP2K: 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 informationTight-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 informationThe 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 informationThe Hydrogen Atom. Chapter 18. P. J. Grandinetti. Nov 6, Chem P. J. Grandinetti (Chem. 4300) The Hydrogen Atom Nov 6, / 41
The Hydrogen Atom Chapter 18 P. J. Grandinetti Chem. 4300 Nov 6, 2017 P. J. Grandinetti (Chem. 4300) The Hydrogen Atom Nov 6, 2017 1 / 41 The Hydrogen Atom Hydrogen atom is simplest atomic system where
More informationSelf-consistent Field
Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced
More informationPseudopotentials: 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 informationElectrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation
Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project
More informationElectronic Properties of Gold Nanoclusters from GW Calculations
Condensed Matter Theory Sector Electronic Properties of Gold Nanoclusters from GW Calculations Thesis submitted for the degree of Doctor Philosophiæ Academic Year 2011/2012 CANDIDATE Jiawei Xian SUPERVISORS
More informationComparison 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 informationElectronic 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 informationAB 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 informationOn-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 informationHomogeneous 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 informationSupplementary Information: Exact double-counting in combining the Dynamical Mean Field Theory and the Density Functional Theory
Supplementary Information: Exact double-counting in combining the Dynamical Mean Field Theory and the Density Functional Theory PACS numbers: THE CORRELATION ENERGY FIT The correlation energy of the electron
More informationCHAPTER 1 NORMCONSERVING PSEUDOPOTENTIALS FOR THE EXACT EXCHANGE FUNCTIONAL. Institut fur Theoretische Physik, Universitat Frankfurt, D-60054
CHAPTER 1 NORMCONSERVING PSEUDOPOTENTIALS FOR THE EXACT EXCHANGE FUNCTIONAL E. Engel, A. Hock, R. N. Schmid and R. M. Dreizler, Institut fur Theoretische Physik, Universitat Frankfurt, D-60054 Frankfurt/Main,
More informationOne-electron Atom. (in spherical coordinates), where Y lm. are spherical harmonics, we arrive at the following Schrödinger equation:
One-electron Atom The atomic orbitals of hydrogen-like atoms are solutions to the Schrödinger equation in a spherically symmetric potential. In this case, the potential term is the potential given by Coulomb's
More informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionals by Michel Côté Université de Montréal Département de physique Eamples Carbazole molecule Inside of diamant Réf: Jean-François Brière http://www.phys.umontreal.ca/~michel_
More information1.2 Rare earth atoms H Ψ=E Ψ, (1.2.1) where the non-relativistic Hamiltonian operator is. H = h2 2m. v ext (r i ) (1.2.
8 1. ELEMENTS OF RARE EARTH MAGNETISM 1.2 Rare earth atoms The starting point for the understanding of the magnetism of the rare earths is the description of the electronic states, particularly of the
More informationElectronic 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 informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
More informationarxiv: v1 [physics.chem-ph] 12 Jul 2017
Hybrid functional pseudopotentials Jing Yang, 1 Liang Z. Tan, 1 and Andrew M. Rappe 1 1 Department of Chemistry, University of Pennsylvania, arxiv:1707.04501v1 [physics.chem-ph] 12 Jul 2017 Philadelphia,
More informationNonlocal exchange correlation in screened-exchange density functional methods
Nonlocal exchange correlation in screened-exchange density functional methods Byounghak Lee and Lin-Wang Wang Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More informationFIRST-PRINCIPLES MODELING OF FUNCTIONAL PEROVSKITE MATERIALS AND SUPERLATTICES
FIRST-PRINCIPLES MODELING OF FUNCTIONAL PEROVSKITE MATERIALS AND SUPERLATTICES By QIBIN ZHOU A dissertation submitted to the Graduate School New Brunswick Rutgers, The State University of New Jersey in
More informationAn 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 informationPseudopotentials: 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 information1 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 informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationLecture #21: Hydrogen Atom II
561 Fall, 217 Lecture #21 Page 1 Lecture #21: Hydrogen Atom II Last time: TISE For H atom: final exactly solved problem Ĥ in spherical polar coordinates Separation: ψ nlml ( r,θ,φ) = R nl (r)y m l (θ,φ)
More informationCLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122
CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX
More informationAlgorithms 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 informationShort-Ranged Central and Tensor Correlations. Nuclear Many-Body Systems. Reaction Theory for Nuclei far from INT Seattle
Short-Ranged Central and Tensor Correlations in Nuclear Many-Body Systems Reaction Theory for Nuclei far from Stability @ INT Seattle September 6-, Hans Feldmeier, Thomas Neff, Robert Roth Contents Motivation
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter ifor a one-dimenional problem to a specific and very important
More informationDensity Functional Theory : Implementation in ABINIT
Sherbrooke 2008 Density Functional Theory : Implementation in ABINIT X. Gonze Université Catholique de Louvain, Louvain-la-neuve, Belgium Density Functional Theory - Implementation, Sherbrooke May 2008
More informationPseudo 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 informationGW Many-Body Theory for Electronic Structure. Rex Godby
GW Many-Body Theory for Electronic Structure Rex Godby Outline Lecture 1 (Monday) Introduction to MBPT The GW approximation (non-sc and SC) Implementation of GW Spectral properties Lecture 2 (Tuesday)
More informationEfficient 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 informationLecture 3. Solving the Non-Relativistic Schroedinger Equation for a spherically symmetric potential
Lecture 3 Last lecture we were in the middle of deriving the energies of the bound states of the Λ in the nucleus. We will continue with solving the non-relativistic Schroedinger equation for a spherically
More informationBehind the "exciting" curtain: The (L)APW+lo method
Behind the "exciting" curtain: The (L)APW+lo method Aug 7, 2016 Andris Gulans Humboldt-Universität zu Berlin Kohn-Sham equation Potential due to nuclei Exchange-correlation potential Potential due to electron
More informationIntroduction 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 informationThe 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 informationAdvanced Solid State Theory SS Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt
Advanced Solid State Theory SS 2010 Roser Valentí and Harald Jeschke Institut für Theoretische Physik, Goethe-Universität Frankfurt i 0. Literatur R. M. Martin, Electronic Structure: Basic Theory and
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationDept 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 informationLocalized 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 informationCombining 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 information2.1 Introduction: The many-body problem
Chapter 2 Smeagol: Density Functional Theory and NEGF s 2.1 Introduction: The many-body problem In solid state physics one is interested in systems comprising many atoms, and consequently many electrons.
More informationNucleon-nucleon interaction
Nucleon-nucleon interaction Shell structure in nuclei and lots more to be explained on the basis of how nucleons interact with each other in free space QCD Lattice calculations Effective field theory Exchange
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter for a one-dimensional problem to a specific and very important
More informationLecture VIII : The pseudopotential
Lecture VIII : The pseudopotentia I. KOHN-SHAM PROBLEM FOR AN ISOLATED ATOM For a one-eectron atom, the Couombic potentia, V ( r) = V (r) = Z/r is sphericay symmetric. The soutions may then be spit into
More informationSolutions to selected problems from Giuliani, Vignale : Quantum Theory of the Electron Liquid
Solutions to selected problems from Giuliani, Vignale : Quantum Theory of the Electron Liquid These problems have been solved as a part of the Independent study class with prof. Neepa Maitra. Compiled
More informationVolume Dependence of N-Body Bound States
Volume Dependence of N-Body Bound States Sebastian König in collaboration with Dean Lee INT Workshop 18-70W, University of Washington Seattle, WA February 9, 2018 arxiv:1701.00279 [hep-lat], to appear
More informationHow 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 informationAb initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms
Ab initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms K. JANKOWSKI a, R. SŁUPSKI a, and J. R. FLORES b a Nicholas Copernicus University 87-100 Toruń,
More informationSet the initial conditions r i. Update neighborlist. Get new forces F i
v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
More information