Plane waves, pseudopotentials and PAW. X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium
|
|
- Cecilia Johnston
- 5 years ago
- Views:
Transcription
1 Plane waves, pseudopotentials and PAW X. Gonze Université catholique de Louvain, Louvain-la-neuve, Belgium 1
2 Basic equations in DFT Solve self-consistently the Kohn-Sham equation H ψ n = ε n ψ n!!! ρ(r ) = ψ *n (r )ψ n (r ) occ n ˆ V +V [ρ] H =T+ Hxc or minimize ψ n (r) H ρ(r) δ mn = ψ m ψ n for m, n occupied set occ ˆ V ψ +E [ρ] Eel {ψ } = ψ n T+ n Hxc n with! Z V (r) = -!!κ a aκ r-r κ How to represent these quantities? Let s try plane waves
3 Prerequisites of plane waves Plane waves e ikr : simple but infinite spatial extent Cannot use a finite set of planewaves for finite systems! Need periodic boundary conditions. Primitive vectors R j, primitive cell volume Ω0 R3 R1 R OK for crystalline solids But : finite systems, surfaces, defects, polymers, nanosystems...? 3
4 The supercell technique Molecule, cluster Surface : treatment of a slab Interface Point defect in a bulk solid The supercell must be sufficiently big : convergence study 4
5 Plane wave basis set : the maths ψ k (r) = ( NΩ0 ) -1/ i(k +G)r u (G) e k G 5
6 Periodic system : wavevectors For a periodic Hamiltonian : wavefunctions characterized by a wavevector k (crystal momentum) in Brillouin Zone ψ m,k (r+r j ) = e Bloch s theorem ψ m,k (r) = ( NΩ0 ) -1/ ψ m,k (r) ik.r j e ik.r u m,k (r) u m,k (r+r j )=u m,k (r) u = periodic part of the Bloch wavefunction Normalization? Born-von Karman supercell ψ m,k (r+n jr j ) = ψ m,k (r) supercell vectors N jr j with N=N1N N 3 1 1= Ωo Ωo u*k (r) u k (r) dr 6
7 Planewave basis set Reciprocal lattice : set of G vectors such that e e igr has the periodicity of the real lattice u k (r) = u k (G) e igr ψ k (r) = ( NΩ0 ) -1/ igr j u k =1 (G) e i(k +G)r G G 1 -igr (Fourier transform) u k (G) = e u k (r) dr Ω Ωo o (k+g) Kinetic energy of a plane wave The coefficients u k (G) for the lowest eigenvectors ( k + G ) decrease fast with the kinetic energy Selection of plane waves determined by a cut-off energy Ecut (k + G) < E Plane wave sphere cut 7
8 Plane waves : the density and potential Fourier transform of a periodic function 1 f(g) = Ωo!r Poisson equation VH (r) = Ωo!r f(r) f(r) = e igr f(g) e-igr f(r) dr! G n (G) and VH (G) n(r') dr' VH r-r' r = -4π n(r) Relation between Fourier coefficients: G VH (G) = 4π n(g) VH (G) = 4π n(g) G For G =0 ( G=0 ) divergence of VH ( G=0 ) 1 n(g=0) = n(r) dr Average Ω Ωor or Compensation with ionic potential 8
9 Representation of the density Density associated with one eigenfunction : n nk (r) = u*nk (r) u nk (r) Computation of u*nk (r) u nk (r) * -igr -ig'r = u nk (G) e u nk (G') e G G' = u*nk (G) u nk (G') e i(g'-g)r GG' Non-zero coefficients for k+g sphere k+g' sphere k+g k+g' G'-G The sphere for n(g) has a doubled radius 9
10 From real space to reciprocal space n(r) = n(g) e igr G sphere() Use of the discrete Fourier transform {r i } {G} 1 -igri n(g) = n(r ) e i N ri {ri } Reciprocal lattice Fast Fourier Transform algorithm Real lattice: original cell FFT grid in real space used to evaluate the XC contribution occ ˆ V ψ +E [ρ] Eel {ψ } = ψ n T+ n Hxc n 10
11 Simplicity of PW requires psps The Fourier transform theory teaches us : -details in real space are described if their characteristic length is larger than the inverse of the largest wavevector norm (roughly speaking) -quality of a plane wave basis set can be systematically increased by increasing the cut-off energy Problem : huge number of PWs is required to describe localized features (core orbitals, oscillations of other orbitals close to the nucleus) Pseudopotentials (or, in general, «pseudization») to eliminate the undesirable small wavelength features (will be the subject of a later section) 11
12 Convergence wrt to kinetic energy cutoff Total energy (Ha) bulk silicon Cut-off energy (Ha) Courtesy of F. Bruneval 1
13 Number of planewaves Number of plane waves = function of the kinetic energy cut-off not continuous d N PW 4 0 Also, a (discontinuous) function of lattice parameter at fixed kinetic energy ( E cut )1 13
14 Discontinuities in energy and pressure => Energy (and pressure) also (discontinuous) functions of lattice parameter at fixed kinetic energy Energy Pressure Lattice parameter Lattice parameter 14
15 Removing discontinuities Kinetic energy Kinetic energy u(g) 0 u(g)=0 u(g) 0 u(g) weak u(g)=0 k+g k+g ecutsm 15
16 Representation : wrap-up Choice of a basis (e.g. Plane waves) Truncating of the basis -> finite basis ( k+g ) < E cut Sphere of plane waves Discontinuous increase of the number of plane waves? Smearing of u(g) => Progressive incorporation of new G vectors Representation of the density Sphere with a double radius in the reciprocal space Going from the real space to reciprocal space Discrete Fourier transform Grid of points + Fast Fourier Transform {r i } {G} 16
17 Pseudopotentials 17
18 Core and valence electrons (I) Core electrons occupy orbitals that are "the same" in the atomic environment or in the bonding environment "the same" depends on the accuracy of the calculation! Separation between core and valence orbitals : the density... N n(r) = ψ i* (r)ψ i (r) i = N core * ψ i (r)ψ i (r) + i core N val * ψ i (r)ψ i (r) = ncore (r) + nval (r) i val atom «Frozen core» for i core : ψ i = ψ i 18
19 Separation between core/valence Depends on the target accuracy of the calculation! (remark also valid for pseudopotentials, with similar cores) For some elements, core/valence partitioning is obvious, for some others, it is not. F atom : IP Ti atom : (1s) + 1keV (s) ( p) ev (1s) (s) ( p) (3s) (3p) (4s) (3d) small core (1s) (s) ( p) (3s) (3p) (4s) (3d) large core 6 6 IP Gd atom : ev 43.3eV small core with n=1,,3 shells, might include 4s, 4p, and 4d in the core. 4f partially filled 19
20 Energy : core and valence 1 1 n(r )n(r ) EKS {ψ i } = ψ i ψ i + Vext (r)n(r) dr + 1 dr1dr + Exc [ n ] r1 - r i EKS {ψ i } = N core i core N val + i val 1 1 n (r )n (r ) ψ i ψ i + Vext (r)ncore (r) dr + core 1 core dr1dr r1 - r 1 1 nval (r1 )nval (r ) ψ i ψ i + Vext (r)nval (r) dr + dr1dr r1 - r nval (r1 )ncore (r ) + dr1dr + Exc [ ncore + nval ] r1 - r 0
21 Valence electrons in screened potential The potential of the nuclei κ is screened by the core electrons ncore,κ (r1 ) Zκ Vion,κ (r) = + dr1 r Rκ r - r1 Z core,κ Z val,κ ncore,κ (r1 ) Vion,κ (r) = + + dr1 r Rκ r Rκ r - r1 The total energy becomes Z val,κ Z val,κ ' 1 E = Eval + Ecore,κ + (κ,κ ') Rκ Rκ ' κ κ κ ' with 1 ψ ψ + V (r) i i ion,κ nval (r) dr i val κ Non-linear XC 1 nval (r1 )nval (r ) core correction + dr1dr + Exc [ ncore + nval ] r1 - r Eval,KS {ψ i } = N val 1
22 Removing core electrons (I) From the previous construction : valence orbitals must still be orthogonal to core orbitals ( => oscillations, slope at the nucleus...) Pseudopotentials try to remove completely the core orbitals from the simulation Problem with the number of nodes This is a strong modification of the system... Pseudopotentials confine the strong changes within a «cut-off radius»
23 Removing core electrons (II) Going from 1 + v 1 + vps To ψ ψ i > = εi ψ i > ps,i > = ε ps,i ψps, i > Possible set of conditions (norm-conserving pseudopotentials) NCPP - Hamann D.R., Schlüter M., Chiang C, Phys.Rev.Lett. 43, 1494 (1979) ε i = εps,i ψi (r) = ψps,i(r) for r > rc r < rc ψi (r) dr = r < rc ψps,i(r) dr For the lowest angular momentum channels (s + p... d...f) Generalisation : ultra-soft pseudopotentials (USPP), projector-augmented plane waves (PAW) 3
24 Example of NC pseudopotential 3s Radial wave function of Si s 0.5 p Vps(r) [Ha] Amplitude All-electron wave function Pseudo-wave function Radial distance [a.u.] Z val r d Radial distance [a.u.] 4
25 Forms of pseudopotentials Must be a linear, hermitian operator General form : (V ψ ) (r) = V ps kernel ps (r, r')ψ (r')dr' Vpskernel (r, r') = Vloc (r)δ (r-r') + Vnloc (r, r') includes a "non-local" part + Spherically symmetric! Vloc (r) = V (r) Vnloc (r, r') = Yℓ*m (θ, ϕ )Vℓ (r, r ')Yℓm (θ ', ϕ ') ℓm Semi-local psp Vℓ (r, r ') = Vℓ (r)δ (r r ') see Bachelet, Hamann and Schlüter, Phys.Rev.B 6, 4199 (198) Separable pspvℓ (r, r ') = ξℓ* (r) fℓξℓ (r ') Kleinman L., Bylander D.M., Phys.Rev.Lett. 48, 145 (198) 5
26 Ultrasoft Pseudopotentials and Projector-Augmented Waves (PAW) 6
27 Ultra-soft pseudopotentials : the idea Problem with NC pseudopotentials : norm conservation limits the softness! When orbitals without nodes (1s, p, 3d, 4f) treated as valence => small characteristic length energy cut-off large. Idea (Vanderbilt, Phys. Rev. B 41, 789 (1990)) Suppress norm-conservation condition : - modify normalization, to keep correct scattering properties - introduce charge density corrections, For selected elements, can decrease number of PW/FFT Grid points by a factor of two or three, with even larger speed up. More difficult to implement than norm-conserving PPs. Can be obtained as a particular case of PAW construction... 7
28 Projector-Augmented Waves : the idea Idea P. Blöchl Phys. Rev. B 50, (1994) The true wave-function and a well-behaving pseudowavefunction are linked by a linear transformation Strong oscillations near the nucleus Ψ = TΨ True wave function (all electrons) No oscillation near the nucleus Auxiliary wave function (soft pseudofunction) More rigorous than USPP 8
29 USPP and PAW : common features Generalized Schrödinger Eq., with overlap operator S. + v ps ψ ps,i = ε i Sψ ps,i Charge density of each state to be corrected for the missing norm. 9
30 Projector-Augmented Waves : the math Ψ Ψ Ψ = TΨ True wave-function Well-behaving pseudo-wavefunction Linked by a linear transformation Physical quantities like computed in the pseudo representation ΨAΨ A Ψ, Ψ A = T + AT with Similarly, variational principle for total energy gives E T Ψ = εt +T Ψ Ψ [ ] equivalent to Kohn-Sham equation, for pseudowavefunctions. Search for ground state can also be done in the pseudo space. 30
31 Transformation operator Operator T has to modify the smooth pseudowavefunction in each atomic region, to give it the correct nodal structure. Identity + sum of atomic contributions Ψ = TΨ with T = 1+ SR R Choose : Partial waves (R=atomic site label) φi = basis set, solutions of the Schrödinger Eq. for the isolated atoms within some cut-off radius rc,r Pseudo partial waves φ i Define S such as : φ i = (1+ SR ) φ i, identical to the partial waves beyond the cut-off radius, but smoother inside 31
32 Representation of the wavefunctions ( T = 1+ φ i φ i Ψ = TΨ i ) p i Explicitly, ( + ΨR1 Ψ R1 Ψ=Ψ R ) + = Plane wave part All electron spherical part ΨR1 = φ i p i Ψ i R Pseudized spherical part R1 = φ i p i Ψ Ψ i R 3
33 Wavefunctions, density, energy represented by plane waves (might use other representations) Ψ represented on a radial grid, centered on R, times spherical harmonics R1 ΨR1 and Ψ Outside of the spheres, Note : Inside one sphere, R1 ΨR1 = Ψ 1 = Ψ R Ψ = φ i p i Ψ i R Density: Energy: n(r ) = n~ (r ) + å n1r (r ) - n~r1 (r ) ~ ~1 1 E = E + å ER - ER ( R ) ( ) R 33
34 Approximations (1) Core electrons : usually treated in the frozen-core approximation, and treated on radial grid (spherical harmonics). () Finite PW basis set (same as PPs) (3) The partial wave expansion is truncated : only one or two partial wave(s), for each atom R, and each l,m channel φ i p i 1 i Inside one sphere, = Ψ R1 Ψ is only approximately true! Ψ = φ i p i Ψ i R contributes inside the atomic spheres, and corrects for the missing terms due to truncation of partial wave expansion 34
35 Advantages of PW+PP or PAW method? (1) The basis set does not depend on atomic coordinates : easy computation of forces as numerically exact derivatives of the total energy with respect to atomic coordinates (no Pulay forces). Easy structural optimisation, or MD. This leads also to Car-Parrinello technique. () Systematic way to complete the basis set For PW-PP : simple implementation For PW-PAW : can be a numerically accurate implementation of DFT, including properties related to cores ; usually faster than PW-PP Disadvantages : cannot lead to Order(N) implementation, does not treat efficiently finite systems (vacuum!). Higher lying states (very high in energy) cannot be trusted... 35
36 Pseudopotential and PAW atomic data availability and reliability 36
37 Generators (non-complete list) Norm-conserving pseudos :!ONCVPSP (Don Hamann, (M. Oliveira (A. Rappe USPP : Vanderbilt USPP generator ( PAW ATOMPAW (N. Holzwarth GPAW generator ( + converter USPPPAW ( 37
38 Implementations Plane Waves need pseudopotentials (or PAW), but pseudopotentials can be used with many basis sets! Norm-conserving pseudos :!with planewave basis set : ABINIT, QE, CASTEP!with localized orbitals : SIESTA...!with wavelets : BigDFT... USPP : with planewave basis set : QE, DaCapo, CASTEP PAW with planewave basis set : ABINIT, QE, GPAW, VASP Same pseudopotentials can sometimes be used by different codes. 38
39 Pseudopotentials/PAW data " Example of PAW atomic dataset table : JTH v1.0 Δ = 0.4 mev Jollet, Torrent, Holzwarth, Computer Physics Comm. 185, 146 (014) Also : GPAW table, GBRV v1.0 table, norm-conserving pseudopotential table (e.g. ONCVPSP pseudo generator), or many other pseudos see 39
40 Testing pseudopotentials 40
41 Comparing code/pseudopotential Recently, large effort to improve uncertainty related to psps : «Delta-factor» collaboration Lejaeghere... Cottenier, Science 351, aad3000 (016) Specification of 71 elemental solids for different volumes. 41
42 More about pseudopotentials / PAW datasets Review 4
43 More about pseudopotentials / PAW datasets Review 43
44 More about pseudopotentials / PAW datasets Review 44
45 Improvement with time Review JTH 1.0 => 0.4 mev 45
46 Other comparison with All-Electron data 46
47 Summary - Plane waves basis set - PW need pseudopotentials - Easy computation of forces relaxation of geometry, or molecular dynamics - PAW (also USPP) : more accurate, faster for cases with 3d and 4f orbitals 47
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 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 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 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 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 informationThe ABINIT code. X. Gonze. Thanks to the > 50 ABINIT contributors, and especially to GM Rignanese for contributions to the slides
The ABINIT code X. Gonze Thanks to the > 50 ABINIT contributors, and especially to GM Rignanese for contributions to the slides 1 ABINIT software project Ideas (1997) : 1) Softwares for first-principles
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 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 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 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 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 informationPseudopotential 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 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 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 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 information6: 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 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 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 informationThe Plane-wave Pseudopotential Method
The Plane-wave Pseudopotential Method k(r) = X G c k,g e i(g+k) r Chris J Pickard Electrons in a Solid Nearly Free Electrons Nearly Free Electrons Nearly Free Electrons Electronic Structures Methods Empirical
More 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 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 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 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 informationDFT / 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 informationFundamentals 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 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 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 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 informationBasis 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 informationReferences. 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 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 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 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 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 informationElectronic 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 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 informationSelf Consistent Cycle
Self Consistent Cycle Step 0 : defining your system namelist SYSTEM How to specify the System All periodic systems can be specified by a Bravais Lattice and and atomic basis How to specify the Bravais
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 informationElectronic 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 informationCHEM6085: 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 informationPractical 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 informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationIntroduction 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 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 informationIntegrated 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 informationIntroduction 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 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 Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
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 informationALMA: 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 informationDFT EXERCISES. FELIPE CERVANTES SODI January 2006
DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a
More 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 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 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 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 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 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 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 informationPractical 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 informationElectronic 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 informationBloch s Theorem: Some Notes
Bloch s Theorem: Some Notes MJ Rutter Michaelmas 005 1 Bloch s Theorem + V (r) ] ψ(r) = Eψ(r) If V has translational symmetry, it does not follow that ψ(r) has translation symmetry. At first glance we
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 informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More information2. 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 informationBand 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 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 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 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 informationVASP: running on HPC resources. University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
VASP: running on HPC resources University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria The Many-Body Schrödinger equation 0 @ 1 2 X i i + X i Ĥ (r 1,...,r
More informationElectron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46
Electron levels in a periodic potential: general remarks Daniele Toffoli January 11, 2017 1 / 46 Outline 1 Mathematical tools 2 The periodic potential 3 Bloch s theorem and Born-von Karman boundary conditions
More informationExperimental Determination of Crystal Structure
Experimental Determination of Crystal Structure Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, U.S.A. PHYS 624: Introduction to Solid State Physics http://www.physics.udel.edu/~bnikolic/teaching/phys624/phys624.html
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 informationComparison 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 informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
More informationMODULE 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 informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
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 informationFirst-Principles Wannier Functions of Silicon and Gallium. Arsenide arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 22 Nov 1996.
First-Principles Wannier Functions of Silicon and Gallium Arsenide arxiv:cond-mat/9611176v1 [cond-mat.mtrl-sci] 22 Nov 1996 Pablo Fernández 1, Andrea Dal Corso 1, Francesco Mauri 2, and Alfonso Baldereschi
More informationCHAPTER 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 informationThree Most Important Topics (MIT) Today
Three Most Important Topics (MIT) Today Electrons in periodic potential Energy gap nearly free electron Bloch Theorem Energy gap tight binding Chapter 1 1 Electrons in Periodic Potential We now know the
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More informationSimulations of Li ion diffusion in the electrolyte material Li 3 PO 4
Simulations of Li ion diffusion in the electrolyte material Li 3 PO 4 a, b N. A. W. Holzwarth Wake Forest University, Winston-Salem, NC, USA Motivation Calculational methods Diffusion in crystalline material
More informationTime-dependent density functional theory : direct computation of excitation energies
CECAM Tutorial Electronic excitations and spectroscopies : Theory and Codes Lyon 007 Time-dependent density functional theory : direct computation of excitation energies X. Gonze Université Catholique
More informationDensity 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 informationAll-electron density functional theory on Intel MIC: Elk
All-electron density functional theory on Intel MIC: Elk W. Scott Thornton, R.J. Harrison Abstract We present the results of the porting of the full potential linear augmented plane-wave solver, Elk [1],
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationAll-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana?
All-Electron Full-Potential Calculations at O(ASA) Speed A Fata Morgana? SFB 484, Teilprojekt D6 October 5, 2007 Outline 1 2 3 Outline 1 2 3 Outline 1 2 3 Outline 1 2 3 Back in the 1930 s... John C. Slater
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 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 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 informationKey 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 informationLecture 10. Transition probabilities and photoelectric cross sections
Lecture 10 Transition probabilities and photoelectric cross sections TRANSITION PROBABILITIES AND PHOTOELECTRIC CROSS SECTIONS Cross section = = Transition probability per unit time of exciting a single
More informationElectronic structure of solids: basic concepts and methods
Electronic structure of solids: basic concepts and methods Ondřej Šipr II. NEVF 514 Surface Physics Winter Term 2016-2017 Troja, 21st October 2016 Outline A bit of formal mathematics for the beginning
More informationAtkins & de Paula: Atkins Physical Chemistry 9e Checklist of key ideas. Chapter 7: Quantum Theory: Introduction and Principles
Atkins & de Paula: Atkins Physical Chemistry 9e Checklist of key ideas Chapter 7: Quantum Theory: Introduction and Principles classical mechanics, the laws of motion introduced in the seventeenth century
More informationSelf-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT
Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,
More informationDFT: Exchange-Correlation
DFT: Local functionals, exact exchange and other post-dft methods Stewart Clark University of Outline Introduction What is exchange and correlation? Quick tour of XC functionals (Semi-)local: LDA, PBE,
More 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 informationIntroduction of XPS Absolute binding energies of core states Applications to silicene
Core level binding energies in solids from first-principles Introduction of XPS Absolute binding energies of core states Applications to silicene arxiv:1607.05544 arxiv:1610.03131 Taisuke Ozaki and Chi-Cheng
More informationarxiv: v1 [cond-mat.str-el] 7 Dec 2008
Order-N implementation of exact exchange in extended systems Xifan Wu, Annabella Selloni, and Roberto Car Chemistry Department, Princeton University, Princeton, NJ 08544-0001,USA (Dated: October 28, 2018)
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 information