Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT
|
|
- Meredith Bryant
- 5 years ago
- Views:
Transcription
1 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, University of Washington, Seattle (b) Pacific Northwest National Laboratory, Richland, USA (c) Faculty of Science, University of Iceland, Reykjavik, Iceland
2 Self-Interaction in DFT Total energy of electronic system in the DFT: ( r) ρ( r ) 1 1 ρ E = Ψ ( r) + V Ψ ( r) dr+ drdr + E [ ρ, ρ ] DFT σ, i ext σ, i xc 1 σ, i r r Interaction of the density of the state Ψ σ, i () r with itself ( self-interaction ): ( r) ρ ( r ) 1 ρσ, i σ, i ESI ρ, i d d σ = + Exc, i,0 0 rr r r ρσ for any approximate exchange-correlation functional. NON-PHYSICAL! In Hartree-Fock approximation, exchange interaction is evaluated exactly, correlation energy is absent, and self-interaction is cancelled out. E HF Exch * * ( r) ( r) ( r ) ( r ) occ occ 1 Ψ Ψ Ψ Ψ { σ, i ( )} Ψ r = drdr r r n m n m n m
3 In this talk: 1. Why self-interaction corrections?. Solving of DFT-SIC functional. Wannier functions. 3. Direct minimization scheme and extension to periodic systems. 4. Exact exchange (hybrid) functional in plane-wave DFT. 5. Applications.
4 Why Apply Self-Interaction? e - H H + 1. Some problems of DFT known (or suspected) to be caused by the SI: Size consistency (dissociation of H + ) Tendency towards charge and spin de-localization. If treated carefully, self-interaction corrections can help calculate accurate band gaps in semiconductors and insulators. Examples of the phenomena studied with DFT-SIC where DFT fails: Band gaps in semiconductors and insulators: Ge becomes a narrow-gap semiconductor; band gaps are generally improved. Qualitatively (and quantitatively) different description of the systems: a) Localization of a vacancy state at TiO rutile (110) surface. b) Complex spin structure of hematite Fe O 3 : formation of a small polaron. c) Localization of spin and charge densities around the structural defects. d) Self-trapping of excitons and holes in SiO.
5 DFT-SIC functional Remove self-interaction in a very straightforward way (Perdew & Zunger, 8): ( r) ρ ( r ) 1 ρ σ, i σ, i E DFT - SIC { ρ, i} σ = EDFT [ ρ()- r ] d d + Exc ρσ, i,0 σ, i r r - r r E SIC E DFT - SIC is an orbital-dependent functional. It can be minimized in two different ways: (1) Variation wrt orbitals () OEP δe () SIC { Ψ, i ( ) i σ } V r δe SIC, i ( ) SIC ( r) = V ρ σ r SIC ( r) = δψ ( r) δρ ( r) σ, i { } In (1) we get an orbital-dependent potential, i.e different Hamiltonians for different states! Although it is an unusual set of Schrodinger equations, we will try to get around having to solve them directly. Scheme () turns out to be much more difficult if at all possible.
6 Calculating the SIC energy ( r) ρ ( r ) 1 ρσ, i σ, i ESIC ψ, i d d σ = + Exc, i,0 rr r r ρσ SIC term is not invariant with respect to the unitary transformations between occupied states Calculating SIC terms with Bloch functions produces zero: SIC scale as, where V is a system volume Orbitals which minimize DFT-SIC functional are known to be spatially localized; For periodic systems, such localized orbitals are Wannier functions. SIC should be calculated with the Wannier functions. Bloch functions vs. Wannier functions delocalized over entire system localized (hopefully with a small unit cell) ikr ( r R) e ( r) Ψ + = Ψ nk nk ( ) The inverse transformation: i i Ψ r = e w r R n, kr k n i R ( ) ( ) i ikr 3 wn d e 3 BZ ( π ) V α V r R = k Ψ r n, k ( )
7 A general unitary transformation mixing Bloch functions with the same crystal momentum but from different bands is used to localize them inside each unit cell: ( ) ( r) U k ( r) Ψ Ψ n, k nm m, k m Marzari and Vanderbilt ( 97) and Silvestrelli (000) developed a scheme for construction of ( ) the transformation producing maximally localized Wannier functions. U k nm In practical implementations, only few different k-vectors are used for Bloch functions. n kr n, k k BZ i ( r R) = Ψ ( r) w% e w(r-r 0 ) w(r-r Γ-point 6 ) k=0 andk=π/r Identify each WF with the piece-wise defined BF
8 Calculating the SIC energy with Wannier functions Wannier functions constructed from single Γ-point (or just limited number of) Kohn-Sham Bloch states, have infinitely many artificial periodic replicas. To use the plane-wave method to calculate SIC, one has to limit integration in the Coulomb SIC term to a single unit cell. To handle extended systems, we replace the Coulomb interaction by a screened Coulomb interaction. A practical example of such screening: f () r N ( 1 exp{ ( r / R) }) 1 = r N 1 (1-(1-exp(-(x/8)**8))**8)/abs(x) 1/abs(x) This kernel is nearly equal to 1/r for r<r and rapidly decays to 0 for r>r
9 Direct Minimization of DFT-SIC functional DFT K-S orbitals Unitary transformation (Maximally) localized Wannier functions Just a useful guess Calculation of SIC terms Self-consistent direct minimization i ( r) ( r) δedft SIC Ψ i δψ Total Energy, Geometric Structure, Ψ r i ( ) Ψ V Ψ = Ψ V Ψ () i ( j) i SIC j j SIC i H ˆ () i ij i H j = Ψ Ψ { } Hermitian at the minimum Self-consistent calculation of the density and wave functions appears to be very important: DFT-SIC total density is typically different from the DFT density; the self-consistent wave functions are well but not necessarily maximally localized. ε i
10 FAQ: 1. Does SIC always lower the total energy? No. LDA total energy is always higher than LDA-SIC; GGA is typically lower than GGA-SIC.. If DFT-SIC raises the DFT total energy why aren t the Bloch functions (giving vanishing SIC) better solutions? DFT-SIC functional calculated with the Bloch functions is NOT size consistent. ε SIC Expect: E =N ε SIC SIC N Wannier functions: Bloch functions: E SIC /N=ε SIC E /N 0 size inconsistent SIC 3. If Bloch functions are excluded from search from the start, how does this formalism describe the metal-insulator transition? The localization character of the Wannier functions appears as part of the solution: exponentially localized WF s indicate insulator, presence of the power-law decaying WF s speaks of a metal; for the latter, SIC is almost zero. 4. What about the OEP for DFT-SIC? As in any KS scheme, OEP solutions will be Bloch functions => standard OEP for DFT-SIC is at least size inconsistent (and has to give zero SIC for periodic systems).
11 Using full SIC strongly overestimates the corrections. On a set of various systems we saw that, in most cases, one needs to use a damping factor of 0.4: EDFT -SIC = EDFT 0.4 ESIC All qualitative effects such as charge/spin localization and opening of the band gaps appear to be insensitive to the value of this factor. Zero-gap semiconductor in LDA Ge Semiconductor in DFT+SIC SIC-driven metal-insulator transition as transition from algebraically to exponentially localized Wannier functions
12 Band gaps in several systems. System Single-Particle Band Gap (ev) Experiment DFT DFT-0.4*SIC Al O TiO (rutile) SiO HfO HfSiO 4? Fe O Si Ge Why Si and Ge are so much overcorrected?
13 Small polarons in hematite (Fe O 3 ) Antiferromagnetic ground state LDA Additional electron density LDA SIC Localization of the vacancy state on the rutile TiO (110) surface
14 1.3 ev N(ε) Conduction band Vacancy state 3 ev N(ε) Conduction band 0.8 ev Vacancy state 3. ev N(ε) Conduction band 1.0 ev Vacancy state Valence band top Valence band top Valence band top DFT (PBE96) Experiment PBE SIC What we learned: 1. DFT-SIC is capable of solving at least some serious problems of DFT.. Using full SIC strongly overcorrects the DFT results in all systems studied and observed that one should use the factor of 0.4 instead in front of the SIC term. 3. Self-consistent SIC are implemented in the NWChem plane-wave package and can be used to study complex systems. 4. Computationally DFT-SIC is about 3 times more expensive than standard DFT (for GGA; the factor is significantly smaller for LDA). SIC part cost scales as N rather than N 3.
15 Exact exchange functionals in plane-wave DFT E Exch * * ( r) ( r) ( r ) ( r ) occ occ 1 Ψ Ψ Ψ Ψ { σ, i ( )} Ψ r = drdr r r This looks much more difficult than n m n m n m * * ( r) ( r) ( r ) ( r ) occ occ Ψn Ψn Ψm Ψ m ρ() r ρ( r ) ECoul [ ρ() r ] = drdr = drdr r r r r n m Even for finite systems, E Exch is very involved, poorly scales with the size, and is a bottleneck in applications to larger molecules or clusters. For periodic systems, Ψ () r n are the Bloch functions, and the sums include all occupied bands as well as integration over k-vectors. Exact exchange term is invariant with respect to the unitary transformation between occupied states. occ m occ * * ( r) ( r ) w ( r R) w ( r R) Ψ Ψ = m m m m m, R
16 E Exch = occ occ drdr n, R m, R * * ( r R) ( r R) ( r R ) ( r R ) w w w w r r n n m m Why is it any better the with the Bloch functions? It gets better if each Wannier function is localized in one unit cell: = R R * w m * ( r R ) w ( r R ) m w n ( r R) n ( r R ) w R R R E Exch = 1 drdr occ occ n m * * ( r) ( r ) ( r) ( r ) w w w w n n m m r r - no k-vector integration w ( r R ) w ( r R) 3 But we can do even better: w1 ( r R) w4 ( r R)
17 If it is possible to create localized Wannier functions for a given system, occ * wn( r R) wn( r R) Frr (, ) = r r is short-range ( r r < L) occ n, R n, R * ( r R) ( r R) ( ) Frr (, ) = w w f r r n n L Furthemore, because of the invariance, Γ-point Bloch functions can be used instead of WF s (only if Coulomb interaction is replaced with the screened Coulomb kernel). Some details of the implementation of exact exchange functionals (HF and hybrids) in the plane-wave NWChem code: Start with maximally localized WF s (for maximum gain) Calculate exact exchange term using these WF s and screened Coulomb kernel Perform direct minimization with respect to Wannier orbitals. WF s do not delocalize during this minimization, possibly because of the effect of the screened kernel.
18 Band gaps in several systems calculated with the hybrid functional PBE0 System Single-Particle Gap (ev) Experiment PBE96 PBE0 Al O TiO (rutile) SiO Si Ge ev N(ε) Conduction band Vacancy state 3 ev N(ε) Conduction band 0.8 ev Vacancy state 3.1 ev N(ε) Conduction band 0.9 ev Vacancy state Valence band top Valence band top Valence band top DFT (PBE96) Experiment PBE0
19
20
21 PBE0 functional { σ, ( r) } [ 1, ] [ 1, ] hybrid Exc = α E Exch Ψ i + βex ρ ρ + Ec ρ ρ α = 0.5 β = 0.75 XC = PBE96 Similar to B3LYP functional.
22 Localization of spin and charge densities by SIC (Al,Si)O system: Al substitutional defect in silica Experiment: Structural symmetry breaking; spin density strongly localized on a single O atom; corresponding bond is longer than the bonds with three other O atoms. PBE96 calculation Cluster calculation in Hartree-Fock
23 Bulk periodic system calculation with DFT - 0.4*SIC
24 (PBE96-0.4SIC) Calculation of Exciton in Rutile (TiO ) ρ ρ singlet triplet ground state excited state unrelaxed structure relaxed structure ρ triplet (electron) 1 ρ singlet not a selftrapped exciton! ρ triplet (hole) 1 ρ singlet
25 Band Gaps in Semiconductors and Insulators True band gap gap = EN ( + 1) EN ( ) + EN ( 1) Kohn-Sham band gap =ε ε KS c v Triplet-Singlet splitting triplet singlet gap = E ( N) E ( N) Discontinuity in V XC δ = gap KS SIC SIC SIC gap =εc εv V DFT-SIC δ 0 v SIC c SIC () r 0 V () r = 0 LDA or GGA δ = 0 Total energy difference between the ground (singlet) state and the lowest triplet state, with the exception of the systems with strongly bound triplet excitons, is a good approximation to a single-particle band gap: for sufficiently large unit cells excitation of ONE electron changes the potential only insignificantly; the lowest excited single-particle state in the absence of excitons, Is at the bottom of conduction band, and has vanishing SIC (as do conduction all band states).
26 In practical implementations, only few different k-vectors are used for Bloch functions. Discretized transformation to Wannier functions becomes: kr i ( ) kr ( ) i i i w% n r R = e Ψ n, = e kr k r e wn( r Ri) k BZ k BZ R For central cell R=0 and in the extreme case of single (Γ-point) k=0 the Wannier function w% ( r n ) is not at all localized and has the periodic images in every cell. If vectors k=0 and k=π/r are included, the periodicity of w% ( r) becomes R. i n Γ-point k=0 and k=π/r
27 ( ) What determines the choice of? U k nm Depends on the type of the problem and on the goals. Heuristic arguments show that the orbitals maximizing the SIC are maximally localized Wannier functions (MLWF). We will try to obtain MLWF as a starting point of the minimization procedure. Once again: n, kr k n i R i i ( r) e w ( r R ) Ψ = ( ) i If we manage to make wn r Ri localized within a cell only one w r R contributes into each Ψ n, k ( r) ( ) n ( ) i Ψn, k r is also maximally localized determined almost piecewise (although it is still a delocalized function). ( ) U k ( r) Our goal is then to find which produces such nm Marzari-Vanderbilt MLWF Minimize the spread functional: Ω= r r n n n Ψ n, k r = 0n r 0n n, r 0 r 0 n = n n ( ) Rn = w n r R ;
28 (PBE96-0.4SIC) Calculation of Exciton in Rutile, and Electron, and Hole In Anatase (TiO ) Self-interaction creates a Coulomb barrier to charge localization: DFT functionals that are not self-interaction free predict either completely delocalized or only partially localized densities (Stokbro et al, Pacchione et al, Jonsson et al). Calculations performed by M. Gabriel Anatase TiO q = + 1 ρ doublet ρ doublet Anatase TiO q = 1 ρ doublet ρ doublet
29 How much SIC is enough? We found that using full SIC strongly overestimates the corrections. On a set of various systems we saw that, in most cases, one needs to use a damping factor of 0.4: EDFT -SIC = EDFT 0.4 ESIC This observation is consistent with the discussion by Perdew et al ( 96) of the exact exchange contribution which is based on the adiabatic connection formula. Qualitative effects such as charge/spin localization and opening of the band gaps appear to be insensitive to the value of this factor. Single-particle energies in DFT-SIC formalism 1. For the orbitals minimizing DFT-SIC functional Ψ V Ψ = Ψ V Ψ () i ( j) i SIC j j SIC i and Hamiltonian matrix is hermitian (Pederson et al, 84). Eigenvalues of this hermitian matrix are also the solutions of a reformulated eigenvalue problem with Hamiltonian with SIC terms (Pederson et al, 84) One can identify the diagonalized Hamiltonian matrix with the single particle energies
30 Atomization energies of selected dimers (kcal/mol). Experiment Hartree-Fock Plane-Wave PBE96 PBE96-0.4*SIC PBE96-1.0*SIC O N P NO Reaction and transition state energies (ev) in silanes. Plane-Wave SiH4 SiH + H SiH4 SiH SiH6 SiH4 + SiH SiH6 SiH4 + H E E Si H SiH Si H 6 SiH 4 + SiH 4 rxn ts E rxn E rxn Si H 6 Si H 4 + H E rxn E ts1 E ts E LDA LDA-0.4*SIC PBE PBE96-1.0*SIC PBE96-0.4*SIC HF B3LYP *) QCISD(T) *) *) P.Nachtigall, K. Jordan, A.Smith, and H.Jonsson (1996)
Orbital Density Dependent Functionals
Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)
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 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 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 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 informationIntermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Intermediate DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Intermediate DFT Lausanne12
More informationXYZ of ground-state DFT
XYZ of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) XYZ of ground-state
More informationThe Self Interaction Correction revisited
The Self Interaction Correction revisited Explicit dynamics of clusters and molecules under irradiation Spectroscopic accuracy at low energy SIC problem : one electron interacts with its own mean-field!
More informationA FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY
A FRESH LOOK AT THE BAND-GAP PROBLEM IN DENSITY FUNCTIONAL THEORY JOHN P. PERDEW PHYSICS & CHEMISTRY, TEMPLE UNIVERSITY PHILADELPHIA, PENNSYLVANIA, USA SUPPORTED BY THE U.S. DEPARTMENT OF ENERGY, EFRC
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 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 informationIntroduction to Density Functional Theory
1 Introduction to Density Functional Theory 21 February 2011; V172 P.Ravindran, FME-course on Ab initio Modelling of solar cell Materials 21 February 2011 Introduction to DFT 2 3 4 Ab initio Computational
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationDFT: Exchange-Correlation
DFT: Exchange-Correlation Local functionals, exact exchange and other post-dft methods Paul Tulip Centre for Materials Physics Department of Physics University of Durham Outline Introduction What is exchange
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
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 informationSpin gaps and spin flip energies in density functional theory
Spin gaps and spin flip energies in density functional theory Klaus Capelle Universidade Federal do ABC Carsten Ullrich University of Missouri-Columbia G. Vignale University of Missouri-Columbia Charge
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 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 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 informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
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 informationModified Becke-Johnson (mbj) exchange potential
Modified Becke-Johnson (mbj) exchange potential Hideyuki Jippo Fujitsu Laboratories LTD. 2015.12.21-22 OpenMX developer s meeting @ Kobe Overview: mbj potential The semilocal exchange potential adding
More informationAdvanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University
Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation
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 informationImproved Electronic Structure and Optical Properties of sp-hybridized Semiconductors Using LDA+U SIC
286 Brazilian Journal of Physics, vol. 36, no. 2A, June, 2006 Improved Electronic Structure and Optical Properties of sp-hybridized Semiconductors Using LDA+U SIC Clas Persson and Susanne Mirbt Department
More informationElectronic correlation and Hubbard approaches
Electronic correlation and Hubbard approaches Matteo Cococcioni Department of Chemical Engineering and Materials Science University of Minnesota Notable failures of LDA/GGA: transition-metal oxides Introduction
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
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 information3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
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 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 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 informationThe effect of the Perdew-Zunger self-interaction correction to density functionals on the energetics of small molecules
The effect of the Perdew-Zunger self-interaction correction to density functionals on the energetics of small molecules Simon Klüpfel, 1, a) Peter Klüpfel, 2 and Hannes Jónsson 2 1) Science Institute of
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 informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
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 informationHECToR CSE technical meeting, Oxford Parallel Algorithms for the Materials Modelling code CRYSTAL
HECToR CSE technical meeting, Oxford 2009 Parallel Algorithms for the Materials Modelling code CRYSTAL Dr Stanko Tomi Computational Science & Engineering Department, STFC Daresbury Laboratory, UK Acknowledgements
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 informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
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 informationDensity functional theory in the solid state
Density functional theory in the solid state Ari P Seitsonen IMPMC, CNRS & Universités 6 et 7 Paris, IPGP Department of Applied Physics, Helsinki University of Technology Physikalisch-Chemisches Institut
More information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
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 informationThe effect of the Perdew-Zunger self-interaction correction to density functionals on the energetics of small molecules
THE JOURNAL OF CHEMICAL PHYSICS 137,124102(2012) The effect of the Perdew-Zunger self-interaction correction to density functionals on the energetics of small molecules Simon Klüpfel, 1,a) Peter Klüpfel,
More informationTowards ab initio device Design via Quasiparticle self-consistent GW theory
Towards ab initio device Design via Quasiparticle self-consistent GW theory Mark van Schilfgaarde and Takao Kotani Arizona State University Limitations to the local density approximation, the GW approximation
More informationCollege of Chemistry, Peking University, Beijing, China. Fritz-Haber-Institut der MPG, Berlin, Germany
KITP Program Excitations in Condensed Matter Localized and Itinerant States in a Unified Picture beyond Density Functional Theory Hong Jiang 1, Patrick Rinke 2 and Matthias Scheffler 2 1 College of Chemistry,
More information2.1 Experimental and theoretical studies
Chapter 2 NiO As stated before, the first-row transition-metal oxides are among the most interesting series of materials, exhibiting wide variations in physical properties related to electronic structure.
More informationNWChem: Hartree-Fock, Density Functional Theory, Time-Dependent Density Functional Theory
NWChem: Hartree-Fock, Density Functional Theory, Time-Depent Density Functional Theory Hartree-Fock! Functionality! Input! Wavefunctions! Initial MO vectors! Direct and semidirect algorithms! Convergence,
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 informationDefects in Semiconductors
Defects in Semiconductors Mater. Res. Soc. Symp. Proc. Vol. 1370 2011 Materials Research Society DOI: 10.1557/opl.2011. 771 Electronic Structure of O-vacancy in High-k Dielectrics and Oxide Semiconductors
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 informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationStudy of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal
Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
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 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 informationDefects in materials. Manish Jain. July 8, Department of Physics Indian Institute of Science Bangalore 1/46
1/46 Defects in materials Manish Jain Department of Physics Indian Institute of Science Bangalore July 8, 2014 Outline 2/46 Motivation. Computational methods. Defects in oxides. Why are defects challenging?
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 informationSemiconductor Physics and Devices Chapter 3.
Introduction to the Quantum Theory of Solids We applied quantum mechanics and Schrödinger s equation to determine the behavior of electrons in a potential. Important findings Semiconductor Physics and
More informationAb-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials
ISSN 0974-9373 Vol. 15 No.3 (2011) Journal of International Academy of Physical Sciences pp. 337-344 Ab-initio Electronic Structure Calculations of α, β and γ KNO 3 Energetic Materials Pradeep Jain and
More informationElectronic structure of correlated electron systems. Lecture 2
Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No
More informationFoster, Adam; Lopez Gejo, F.; Shluger, A. L.; Nieminen, Risto Vacancy and interstitial defects in hafnia
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Foster, Adam; Lopez Gejo, F.; Shluger,
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 informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationElectronic Structure of Surfaces
Electronic Structure of Surfaces When solids made of an infinite number of atoms are formed, it is a common misconception to consider each atom individually. Rather, we must consider the structure of the
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system of
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 informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationPrinciples of Quantum Mechanics
Principles of Quantum Mechanics - indistinguishability of particles: bosons & fermions bosons: total wavefunction is symmetric upon interchange of particle coordinates (space,spin) fermions: total wavefuncftion
More informationSupplemental Material: Experimental and Theoretical Investigations of the Electronic Band Structure of Metal-Organic Framework of HKUST-1 Type
Supplemental Material: Experimental and Theoretical Investigations of the Electronic Band Structure of Metal-Organic Framework of HKUST-1 Type Zhigang Gu, a Lars Heinke, a,* Christof Wöll a, Tobias Neumann,
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 informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationA new generation of density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem
A new generation of density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem Andreas Görling, Patrick Bleiziffer, Daniel Schmidtel, and Andreas Heßelmann Lehrstuhl für
More informationMetal-insulator transitions
Metal-insulator transitions What can we learn from electronic structure calculations? Mike Towler mdt26@phy.cam.ac.uk www.tcm.phy.cam.ac.uk/ mdt26 Theory of Condensed Matter Group Cavendish Laboratory
More informationGeneralized generalized gradient approximation: An improved density-functional theory for accurate orbital eigenvalues
PHYSICAL REVIEW B VOLUME 55, NUMBER 24 15 JUNE 1997-II Generalized generalized gradient approximation: An improved density-functional theory for accurate orbital eigenvalues Xinlei Hua, Xiaojie Chen, and
More informationMany electrons: Density functional theory
Many electrons: Density functional theory Bedřich Velický V. NEVF 54 Surface Physics Winter Term 0-0 Troja, 8 th November 0 This class is devoted to the many-electron aspects of the solid state (computational
More informationTopological Physics in Band Insulators IV
Topological Physics in Band Insulators IV Gene Mele University of Pennsylvania Wannier representation and band projectors Modern view: Gapped electronic states are equivalent Kohn (1964): insulator is
More informationExact exchange (spin current) density-functional methods for magnetic properties, spin-orbit effects, and excited states
Exact exchange (spin current) density-functional methods for magnetic properties, spin-orbit effects, and excited states DFT with orbital-dependent functionals A. Görling Universität Erlangen-Nürnberg
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
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 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 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 informationFirst principles simulations of materials and processes in photocatalysis
First principles simulations of materials and processes in photocatalysis Work with: Annabella Selloni Department of Chemistry, Princeton University Ulrich Aschauer, Jia Chen, Hongzhi Cheng, Cristiana
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationDFT with Hybrid Functionals
DFT with Hybrid Functionals Sanliang Ling University College London 4th CP2K Tutorial, 31st August 4th September 2015, Zurich What are hybrid functionals? Hybrid functionals: mixing non-local Hartree-Fock
More information6.1 Nondegenerate Perturbation Theory
6.1 Nondegenerate Perturbation Theory Analytic solutions to the Schrödinger equation have not been found for many interesting systems. Fortunately, it is often possible to find expressions which are analytic
More informationRandom-phase approximation and beyond for materials: concepts, practice, and future perspectives. Xinguo Ren
Random-phase approximation and beyond for materials: concepts, practice, and future perspectives Xinguo Ren University of Science and Technology of China, Hefei USTC-FHI workshop on frontiers of Advanced
More informationCalculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
More informationCompressed representation of Kohn-Sham orbitals via selected columns of the density matrix
Lin Lin Compressed Kohn-Sham Orbitals 1 Compressed representation of Kohn-Sham orbitals via selected columns of the density matrix Lin Lin Department of Mathematics, UC Berkeley; Computational Research
More informationOne-Electron Hamiltonians
One-Electron Hamiltonians Hartree-Fock and Density Func7onal Theory Christopher J. Cramer @ChemProfCramer 2017 MSSC, July 10, 2017 REVIEW A One-Slide Summary of Quantum Mechanics Fundamental Postulate:
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More information