Strategies for Solving Kohn- Sham equations
|
|
- Gwenda Smith
- 5 years ago
- Views:
Transcription
1 Strategies for Solving Kohn- Sham equations Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany 1 1
2 Appetizer: high-k oxides McKee for new transitors source gate oxide gate drain n doped n doped Hu et al. p doped p doped channel silicon waver Metal Oxide Field Effect Transistor (MOSFET) Norga et Ojima,Yoshim 2 2
3 Questions Which functional? LDA,GGA, Meta-GGA, Hybrid,LDA+U,... Which bandstructure method? Gaussian, Slater-orbitals, Pseudopotentials, KKR, LMTO, LAPW, PAW Which relativistic treatment? non-relativistic, scalar relativistic, spin-orbit coupling? Which spin description? non-spin-polarized, spinpolarized, non-colllinear Which kpoint set? Which supercells? How to decouple spurious periodic images? 3 3
4 Starting point: DFT Minimize total energy functional: E[ ψ n,f n ] = f n ψ n ˆp 2 2m ψ n n with E kin + E xc [n(r)] E xc n,m n(r) = f n ψ(r) 2 n electron density d 3 r Solve Kohn-Sham equations: 2 2m 2 + v eff (r) ψ n (r) =ψ n (r) n and v eff (r) = d 3 r e2 [n(r)+z(r)][n(r )+Z(r )] 4π 0 r r E Hartree Λ n,m (ψ n ψ m δ n,m ) constraint Z(r) = Z j δ(r R j ) j charge density of nuclei in electron charges d 3 r e2 [n(r )+Z(r )] 4π 0 r r v Hartree (r) + δe xc δn(r) µ xc (r) 4 4
5 Starting point: DFT Minimize total energy functional: E[ ψ n,f n ] = f n ψ n ˆp 2 2m ψ n n with E kin + E xc [n(r)] E xc n,m n(r) = f n ψ(r) 2 n electron density d 3 r Solve Kohn-Sham equations: 2 2m 2 + v eff (r) ψ n (r) =ψ n (r) n and v eff (r) = d 3 r e2 [n(r)+z(r)][n(r )+Z(r )] 4π 0 r r E Hartree Λ n,m (ψ n ψ m δ n,m ) constraint Z(r) = Z j δ(r R j ) j charge density of nuclei in electron charges d 3 r e2 [n(r )+Z(r )] 4π 0 r r v Hartree (r) + δe xc δn(r) µ xc (r) 5 5
6 Starting point: DFT Minimize total energy functional: E[ ψ n,f n ] = f n ψ n ˆp 2 2m ψ n n with E kin + E xc [n(r)] E xc n,m n(r) = f n ψ(r) 2 n electron density d 3 r Solve Kohn-Sham equations: 2 2m 2 + v eff (r) ψ n (r) =ψ n (r) n and v eff (r) = d 3 r e2 [n(r)+z(r)][n(r )+Z(r )] 4π 0 r r E Hartree Λ n,m (ψ n ψ m δ n,m ) constraint Z(r) = Z j δ(r R j ) j charge density of nuclei in electron charges d 3 r e2 [n(r )+Z(r )] 4π 0 r r v Hartree (r) + δe xc δn(r) µ xc (r) 6 6
7 Starting point: DFT Minimize total energy functional: E[ ψ n,f n ] = f n ψ n ˆp 2 2m ψ n n with E kin + E xc [n(r)] E xc n,m n(r) = f n ψ(r) 2 n electron density d 3 r Solve Kohn-Sham equations: 2 2m 2 + v eff (r) ψ n (r) =ψ n (r) n and v eff (r) = d 3 r e2 [n(r)+z(r)][n(r )+Z(r )] 4π 0 r r E Hartree Λ n,m (ψ n ψ m δ n,m ) constraint Z(r) = Z j δ(r R j ) j charge density of nuclei in electron charges d 3 r e2 [n(r )+Z(r )] 4π 0 r r v Hartree (r) + δe xc δn(r) µ xc (r) 7 7
8 What is the problem? 1. wave function oscillates strongly in the atomic region (Coulomb singularity) 2. wave function needs to be very flexible in the bonding region and the tails (Chemistry) K.Schwarz Talk/WWW From APW to LAPW to (L)APW+lo 3. most electrons (core) are irrelevant 4. relativistic effects 5. tiny but finite nucleus r atomic valence wave functions (s and d) of Pt 8 8
9 Brute force: LCAO Use atom-centered basisset Gaussians (Gaussian, Turbomole, etc) Slater orbitals (ADF) Numeric atomic orbitals (DMOL) LCAO Ansatz: ψ n (r) = α wave function solve eigenvalue problem basis orbital χ α (r) coefficient c α,n [H n O]c n = 0 with H α,β = χ α ˆp 2 2m +ˆv χ β and O α,β = χ α χ β 9 9
10 Strategies Pseudopotential Augmented waves Chop off singular potential node-less wave functions no core states no information on inner electrons transferability problems Hamann-Bachelet-Schlueter, Kerker, Kleinman-Bylander, Troullier Martins, Ultrasoft,etc start with envelope function replace incorrect shape with atomic partial waves complex basisset retain full information on wave functions and potential LMTO, ASW, LAPW, APW, PAW, (all-electron) 10 10
11 Energy-dependent partial waves ε-dependent radial Schrödinger equation 2 2m 1 r 2 r r ( + 1) r 2 + v eff (r) φ (,r)=0 Fe-d partial wave ε Linear approximation: r[a 0 ] φ(,r)=φ( ν,r)+( ν ) φ( ν,r) +O(( ν ) 2 ) two functions contain most chemical information φ/ 11 11
12 Logarithmic derivative and Wigner rule Logarithmic derivative function describes the scattering properties of an atom D () = rφ (,r) φ (,r) D[ ε] ϕ () =N nodes () π arctan[d ()] 5 4 s 0 ε B ε A ε B εa ε d p ε [ev] 12 12
13 Augmentation (LAPW) Use envelope functions: e.g. plane waves Define atom-centered augmentation spheres Ω R muffin-tin spheres, atomic spheres expand envelope function into spherical harmonics χ G =e i Gr χ 1 G(r) =e i Gr =4π,m i j ( G r )Y,m( G)Y,m (r) keep value and derivative of radial function for each,m Replace radial part of envelope functions with partial waves from the instantaneous spherically averaged potential 13 13
14 Augmented waves Linear methods (O.K. Andersen75) linear augmented plane wave method (LAPW), linear muffin-tin orbital method (LMTO), augmented spherical waves (ASW) Parents (coming back..) Augmented plane wave (APW) method (Slater53), Korringa-Kohn-Rostocker (KKR) method (Korringa47 and Kohn-Rostocker54) Descendent Projector augmented wave (PAW) method (Blöchl94) 14 14
15 Pseudopotentials Important is the logarithmic derivative and its energy derivative value and derivative result from condition of differentiability energy derivative of the log-derivative is related to the norm dd d = 2m 2 r 2 Ω φ2 (,r Ω ) d 3 r φ(, r) 2 Ω Norm Replace partial waves by nodeless (pseudo-) partial waves with the correct norm Invert pseudo partial waves to obtain potential 2 2m 2 + v ps (r) φ(r)! =0 v ps (r) = 1 φ(r) 2 2m 2 φ(r) 15 15
16 Pseudopotentials Invert pseudo partial waves to obtain potential 2 2m 2 + v ps (r) φ(r)! =0 v ps (r) = 1 φ(r) 2 2m 2 φ(r) one potential for each angular momentum Project onto angular momenta before applying potential semi-local potential ˆv ps =ˆv local +,m ˆP,m (ˆv ps, ˆv local ) ˆP,m 16 16
17 Pseudopotentials pseudo-schrödinger equation 2 2m 2 +ˆv ψ(r) =0 with ˆv = ˆv local + ˆP,m (ˆv ps, ˆv local ) ˆP,m,m + d 3 r e2 [ñ(r)+ Z(r)] + δe xc[ñ] 4π 0 r r δn(r) ñ(r) = n valence Z(r) = (N core Z) v eff (ñ(r)],r)] f n ψ n (r) 2 a π 3 e ar2 pseudo charge density d 3 r e2 [ñ at (r)+ Z(r)] 4π 0 r r v eff (ñ at (r)],r)] pseudo core charge + δe xc[ñ at ] δn(r) 17 17
18 Plane waves ψ k,n (r) = G e i( k+ G)r c k,n ( G) with G = g1 i 1 + g 2 i 2 + g 3 i 3 convergence controlled by a single parameter 1 π r = 2 k + G 2 E PW 2EPW Fast Fourier Transform (FFT) map to real space in N ln(n) operations (<< N 2!) u k,n (r) = e igr c k,n ( G) ; c k,n = 1 e igr u k,n (r) N r G kinetic energy and Coulomb energy are diagonal in G-space E kin = V unit-cell 2 2m ( k + G) 2 c k,n ( G) 2 G potential and exchange correlation is diagonal in r-space r 18 18
19 Self-Consistency Cycle 19 19
20 Self-consistency cycle SCF-Cycle 1.Schrödinger equation 2.density from wavefunction 3.potential from density 4.mixing Mixing is critical potential or density mixing charge-sloshing for large systems Mixing is an art! in v eff =(1 p 2 φ n n v eff (r ) v 2m eff ε n φ n =0 e φ n n(r )=Σ * (r ) (r ) n(r ) φ n out v eff (r ) = v ext(r ) + d 3 e r 2 n(r ) 4 πε 0 r r + µ xc (r ) α)vin out eff +αv eff in in out in out v eff =v eff? v eff (r ) 20 20
21 Paradigm shift: Car-Parrinello method eigenvalue problems minimize energy Access to the molecular dynamics Car, Parrinello, PRL 55, 2471(1985) 21 21
22 Fictitious Lagrangian L = f n ψ n ˆm ψ ψ n M j R 2 j n j fictitious kinetic energy E kin of nuclei Equations of motion Car-Parrinello method ˆm ψ ψ n f n = Ĥ ψ nf n E DF T / ψ n Car, Parrinello, PRL 55, 2471(1985) E DFT {f n, ψ n, R j } Λ n,m (ψ n ψ m δ n,m ) n,m constraints ˆm ψ ψf n m friction ψ m Λ m,n force of constraint M j Rj = Fj M j α R R j Rj E DF T friction Discretize using Verlet algorithm ẍ ẋ x(t + ) 2x(t)+x(t ) 2 x(t + ) x(t ) 2 Constraints using [Ryckaert,Cicotti, Berendsen, J. Comp.Phys. 23, 327 (1977)] 22 22
23 K-points and Brillouin-zone integration 23 23
24 K-points and Band structure Wave functions in a periodic crystal can be written as Bloch states: periodic function modulated by a plane wave periodic ψ k,n (r) = u k,n (r)e i kr ψ k,n (r) = χ α r ( R α + t) α Ω 0 t e i k( R α +t) c α, k,n z L x U S X Z W K y 24 24
25 Monkhorst-Pack Special Points Replace integral by sum over equispaced grid A = n 1 Fermi function d 3 k f( n ( V k)) A n ( k) g V g matrix element Brillouin-zone integral H.J. Monkhorst and J.D.Pack, PRB13, 5188 (1976) 1 N k N k Corresponds to an interpolation by a plane wave expansion j=1 f( n ( k))a n (k j ) k-point sum f( n ( k))a n ( k)= {t} e i kt c n (t) rapid (exponential) convergence poor description of the Fermi surface control grid density with a cutoff t <R x ;(R x 20) 25 25
26 Special points for metals produce instabilities in the self-consistency non-differential total energies Finite temperatures Introduce finite temperatures to smoothen density of states: Mermin functional [Mermin, PR137, A1441 (1965)] F ({ ψ n,f n },T) free energy = E({ ψ n,f n }) + k B T [f n ln(f n )+(1 f n )ln(1 f n )] n internal energy U TS with entropy S produce finite temperature Fermi distribution f n = effectively smears out energies εn by kbt 1 1+e 1 k B T ( n µ) Density of states 26 ε 26
27 Tetrahedron method Jepsen, Andersen, sol. st. comm. 9, 1763 (1971); Lehmann,Taut, phys. stat. sol. 54, 469 (1972) linear interpolation inside tetrahedra analytic integration of interpolated functions. Blöchl correction No misweighting: produces special-point ε result for filled bands corrects for linear over-estimation of partially filled bands Blöchl, Jepsen, Andersen, PRB49, (1994) µ k x k y Writes result as a weighted sum k x from A = tetrahedron method w n ( k i ) A n ( k i ) ki,n 27 27
28 k-point convergence Tetrahedron method + correction (preliminary) 10 0 Al Fe E-E 0 /ev E-E 0 /ev E-E 0 /ev MgO R x E-E 0 /ev R x Si k < 2π R x 28 28
29 Supercells How to deal with non-crystals? Make super-crystals out of them! molecules surfaces interfaces point defects amorphous (disordered) keep 0.6 nm vacuum between hydrogen atoms, - more for lone-pairs, dangling bonds etc, electropositive atoms keep 1 nm between defects in a solid - more for very dispersive bands choose isotropic cells (z.b. face-centered cubic) - maximize distance/volume, i.e. accuracy/cost 29 29
30 Frozen-core approximation 30 30
31 Frozen-core approximation Core electrons do not take part in chemical binding Rounding errors: core energy are 10 7 times larger than chemical accuracy (2 32 =10 9 ) Most electrons are core electrons are savings possible? 31 31
32 Frozen-core approximation Import atomic core wave functions as basisfunctions for the crystal n c ψ fc n = m=1 ψ at mu m,n One-particle expectation values (kinetic energy, density, etc.) are invariant under a unitary transformation of core states  fc = = n c n=1 n c m,o=1 ψ fc n  ψfc n = ψ at m  ψat o n c n c n=1 m,o n c n U o,n U n,m δ m,o U m,nψ at m  ψat o U o,n = n c n=1 ψ at n  ψat n =  at 32 32
33 Frozen-core approximation Frozen-core wave functions are not identical to atomic wave functions! errors in the total energy are second-order in the deviation from the frozen-core wave function quality of the F-C approximation depends on the exact solution of the atomic problem 33 33
34 Spins 34 34
35 Spin orbitals: Spins an electron is represented by two wave functions, one for spin and one for spin spin-density is a three-dimensional vector Obtain non-collinear spin density as 3-d vector field in 3-d s(r) = 2 ψ (r) ψ (r) ψ (r)ψ (r)+ψ (r)ψ (r) iψ (r)ψ (r)+iψ (r)ψ (r) ψ (r)ψ (r) ψ (r)ψ (r) 35 35
36 Spins Non-collinear: use two-component wave functions most general formulation spin-polarized, collinear (local spin density approximation, LSDA) one of the two components are set to zero electrons are in eigenstates of sz spin density points along z-direction non-spin polarized all electrons are paired zero spin density 36 36
37 Done 37 37
The electronic structure of materials 2 - DFT
Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie
More 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 informationDFT 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 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 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 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 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 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 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 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 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 informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationDensity Functional Theory in Practice
A 7 Density Functional Theory in Practice Stefan Blügel Institut für Festkörperforschung Forschungszentrum Jülich GmbH Contents 1 Introduction 2 2 Kohn-Sham Approach in a Nutshell 5 2.1 Total Energy and
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 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 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 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 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 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 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 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 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 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 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 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 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 informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
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 informationFULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY
FULL POTENTIAL LINEARIZED AUGMENTED PLANE WAVE (FP-LAPW) IN THE FRAMEWORK OF DENSITY FUNCTIONAL THEORY C.A. Madu and B.N Onwuagba Department of Physics, Federal University of Technology Owerri, Nigeria
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 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 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 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 informationLecture on First-principles Computations (11): The Linearized Augmented Plane-wave Method
Lecture on First-principles Computations (11): The Linearized Augmented Plane-wave Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2018.10.24 Recall: orthogonalized
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 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 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 informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
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 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 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 informationCHAPTER 6. ELECTRONIC AND MAGNETIC STRUCTURE OF ZINC-BLENDE TYPE CaX (X = P, As and Sb) COMPOUNDS
143 CHAPTER 6 ELECTRONIC AND MAGNETIC STRUCTURE OF ZINC-BLENDE TYPE CaX (X = P, As and Sb) COMPOUNDS 6.1 INTRODUCTION Almost the complete search for possible magnetic materials has been performed utilizing
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? Center for Electronic Correlations and Magnetism Institute for Physics, University of Augsburg February 4, 2008 Outline 1 2 3 Outline
More informationSystematic 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 informationOptimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics
Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Sangeeta Sharma 1,2, J. K. Dewhurst 3, C. Ambrosch-Draxl 4, S. Pittalis 2, S. Kurth 2, N. Helbig 2, S. Shallcross
More informationMODULE 2: QUANTUM MECHANICS. Principles and Theory
MODULE 2: QUANTUM MECHANICS Principles and Theory You are here http://www.lbl.gov/cs/html/exascale4energy/nuclear.html 2 Short Review of Quantum Mechanics Why do we need quantum mechanics? Bonding and
More informationPerformance ofhybrid density functional methods,screened exchange and EXX-OEP methodsin the PAW approach p.1/26
Performance of hybrid density functional methods, screened exchange and EXX-OEP methods in the PAW approach Georg Kresse, J Paier, R Hirschl, M Marsmann Institut für Materialphysik and Centre for Computational
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More 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 informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationLecture 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 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 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 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 informationfile://c:\documents and Settings\lhebditc\Desktop\EAC-116.htm
Page 1 of 40-1- B3.2 Quantum structural methods for the solid state and surfaces Frank Starrost and Emily A Carter B3.2.1 INTRODUCTION We are entering an era when condensed matter chemistry and physics
More informationFirst-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design
First-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design Peter Kratzer 5/2/2018 Peter Kratzer Abeokuta School 5/2/2018 1 / 34 Outline
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 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 informationSloppy Nuclear Energy Density Functionals: effective model optimisation. T. Nikšić and D. Vretenar
Sloppy Nuclear Energy Density Functionals: effective model optimisation T. Nikšić and D. Vretenar Energy Density Functionals the nuclear many-body problem is effectively mapped onto a one-body problem
More informationA Green Function Method for Large Scale Electronic Structure Calculations. Rudolf Zeller
A Green Function Method for Large Scale Electronic Structure Calculations Rudolf Zeller Institute for Advanced Simulation, Forschungszentrum Jülich Electronic structure calculations (density functional
More informationPBS: FROM SOLIDS TO CLUSTERS
PBS: FROM SOLIDS TO CLUSTERS E. HOFFMANN AND P. ENTEL Theoretische Tieftemperaturphysik Gerhard-Mercator-Universität Duisburg, Lotharstraße 1 47048 Duisburg, Germany Semiconducting nanocrystallites like
More 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 informationBasics of DFT applications to solids and surfaces
Basics of DFT applications to solids and surfaces Peter Kratzer Physics Department, University Duisburg-Essen, Duisburg, Germany E-mail: Peter.Kratzer@uni-duisburg-essen.de Periodicity in real space and
More informationAtomic Structure and Atomic Spectra
Atomic Structure and Atomic Spectra Atomic Structure: Hydrogenic Atom Reading: Atkins, Ch. 10 (7 판 Ch. 13) The principles of quantum mechanics internal structure of atoms 1. Hydrogenic atom: one electron
More informationInteraction between a single-molecule
Interaction between a single-molecule magnet Mn 12 monolayer and a gold surface 12 Kyungwha Park Department of Physics, Virginia Tech Salvador Barraza-Lopez (postdoc) Michael C. Avery (undergraduate) Supported
More informationKey concepts in Density Functional Theory
From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationDFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
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 informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More informationBasics of density-functional theory and fast guide to actual calculations Matthias Scheffler
Basics of density-functional theory and fast guide to actual calculations Matthias Scheffler http://www.fhi-berlin.mpg.de/th/th.html I. From the many-particle problem to the Kohn-Sham functional II. From
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 informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More information5 The Band Structure of the Electrons
5 The Band Structure of the Electrons 5.1 Introduction In the Fourier representation, i.e., in the basis of plane waves χ n (k) = 1 Vg e i(k+gn)r, (5.1) the Hamilton operator of the Kohn-Sham equation
More informationOptimized Effective Potential method for non-collinear Spin-DFT: view to spindynamics
Optimized Effective Potential method for non-collinear Spin-DFT: view to spindynamics S. Sharma and E. K. U. Gross nstitut für Theoretische Physik, FU Berlin Fritz-Haber nstitut, Berlin 21 Sept 2006 Synopsis
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 informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationLecture 15 From molecules to solids
Lecture 15 From molecules to solids Background In the last two lectures, we explored quantum mechanics of multi-electron atoms the subject of atomic physics. In this lecture, we will explore how these
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 information2 Theoryand Practice of Density-Functional Theory
2 Theoryand Practice of Density-Functional Theory Peter E. Blöchl Institute for Theoretical Physics Clausthal University of Technology Contents 1 Introduction 2 2 Basicsofdensity-functional theory 3 3
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 information