Density Functional Theory Barry T Pickup Department of Chemistry - - PowerPoint PPT Presentation

density functional theory
SMART_READER_LITE
LIVE PREVIEW

Density Functional Theory Barry T Pickup Department of Chemistry - - PowerPoint PPT Presentation

Density Functional Theory Barry T Pickup Department of Chemistry University of Sheffield Density Functional Theory Over 90% of modern Quantum Chemistry calculations use some form of density functional methodology. Over 95% of band theory


slide-1
SLIDE 1

Density Functional Theory

Barry T Pickup Department of Chemistry University of Sheffield

slide-2
SLIDE 2

Over 90% of modern Quantum Chemistry calculations use some form of density functional methodology. Over 95% of band theory calculations in Solid State Physics are based upon some form of density functional approach.

Density Functional Theory

slide-3
SLIDE 3

Two theorems by Pierre Hohenberg and Walter Kohn (Nobel Prize 1998):

  • P. Hohenberg and W. Kohn, Phys. Rev. 136B, 864

(1964)

Density Functional Theory

slide-4
SLIDE 4

N-electron time-independent Schrödinger equation:

H T V W

H E

The N-electron Hamiltonian

Density Functional Theory

Kinetic energy Interaction External potential

slide-5
SLIDE 5

The kinetic energy operator: The two-electron interaction operator: The ‘external’ confining potential:

2 1

1 2

N i i

T

1

1 1 2

N i j ij

W r

1 N A i i i A Ai

Z V V r

B A i j rij rAi

Density Functional Theory

slide-6
SLIDE 6

A ‘confining’ potential implies a spectrum of ‘bound’ states:

Energy E0 E1 E2 E3 Ground state

Density Functional Theory

slide-7
SLIDE 7

The N-electron wave function is a function of N space- spin co-ordinates

The Born Probability interpretation of the wave function

(1,2, , , ) ( , , , ) ( , )

i i i i i i

i N i x y z r   

2

(1,2, ) (1,2, ) P N N  

Density Functional Theory

slide-8
SLIDE 8

The electron density:

2 1 1

( ) (1,2, , ) 2 3 r N N d d d dN   

The electron density is the probability of finding any electron at spatial position .

1

r 

Density Functional Theory

slide-9
SLIDE 9

The Hohenberg-Kohn theorem

Theorem I For any system of interacting particles in an external potential V, the density , is uniquely determined, i.e. the external potential is a unique functional of the density.

Density Functional Theory

slide-10
SLIDE 10

Proof of Theorem I Assume two different external potentials V(1) and V(2) give rise to the same density, . Define two distinct Hamiltonians, Using the variation principle for case 1: Inequality strict if ground state is non-degenerate.

Density Functional Theory

( ) ( )

, 1,2

i i

H T V W i

(1) (1) (1) (1) (2) (1) (2)

E H H

Ground state wave functions

slide-11
SLIDE 11

Proof of Theorem I Inequality strict if ground state is non-degenerate. We can write the r.h.s as So that, finally

Density Functional Theory

1 (1) (1) (1) (2) (1) (2)

E H H

(2) (1) (2) (2) (2) (2) (2) (1) (2) (2) (2) (2) (2) (1) (2)

[ ( ) ( )] ( ) H H H H E dr V r V r r    

1 (2) (1) (2)

[ ( ) ( )] ( ) E E dr V r V r r    

slide-12
SLIDE 12

Proof of Theorem I So that, finally Repeating this process for case 2, leads to Reductio ad absurdum

Density Functional Theory

(1) (2) (1) (2)

[ ( ) ( )] ( ) E E dr V r V r r    

(2) (1) (2) (1)

[ ( ) ( )] ( ) E E dr V r V r r    

(1) (2) (1) (2) (1)

[ ( ) ( )] ( ) E E dr V r V r r E    

slide-13
SLIDE 13

The Hohenberg-Kohn theorem

Theorem II A universal functional of the density , can be

  • defined. The exact ground state is the global

minimum value of this functional.

Density Functional Theory

slide-14
SLIDE 14

The Hohenberg-Kohn theorem

  • 1. The ground state density, , uniquely determines

the ground state wave-function ( ), as well as the external potential V( ).

  • 2. As a consequence, any observable of a static

many-particle system is a functional of the ground- state density.

Density Functional Theory

slide-15
SLIDE 15

The Hohenberg-Kohn theorem

3.The total energy functional

  • f a system characterised by external potential, V, is

equal to the exact ground state energy if and only if the exact ground state density,

0, is inserted.

[ ] [ ] [ ]

V

E T W V

Density Functional Theory

slide-16
SLIDE 16

The Hohenberg-Kohn theorem

  • 4. For all other densities,

, Consequently, the exact density, , and the exact ground state energy can be obtained by solving the Euler-Lagrange equation

[ ]

V

E E

[ ] ( )

V

E r

Density Functional Theory

slide-17
SLIDE 17

The Hohenberg-Kohn theorem

  • 5. The functional

is universal, i.e. it is independent of the external potential, V, of the particular many-particle system. It has the same form for ALL systems with a fixed inter-particle interaction, W (including W = 0).

[ ] [ ] [ ] F T W

Density Functional Theory

slide-18
SLIDE 18

The Kohn-Sham Equations (Kohn & Sham 1966)

Kohn and Sham introduced a fictitious non- interacting system with a Slater determinant as the exact ground state wave-function. This Slater determinant possesses the same density, , and energy E, as the exact interacting system, and must necessarily have a different potential, VKS.

Density Functional Theory

slide-19
SLIDE 19

The Kohn-Sham Equations

The Kohn-Sham satisfies a set of orbital equations analogous to the HF equations where the density

Density Functional Theory

2

1 [ ] [ ]( , ) [ ] [ ]( , ), 1, , 2

KS i i i

V r r i N   

2 1

( ) [ ]( , )

N i i

r r  

slide-20
SLIDE 20

The Kohn-Sham Equations

The Kohn-Sham potential is defined as where the exchange-correlation potential is

Density Functional Theory

3

( ') [ ]( ) ( ) ' [ ]( ) '

KS xc

r V r V r d r V r r r       [ ] [ ]( ) ( )

xc xc

E V r r  

slide-21
SLIDE 21

Implementations

B3LYP (Becke, Lee, Yang and Perdew) where , a0 = 0.2, ax = 0.72, and ac = 0.81, are three empirical parameters determined by fitting the calculated quantities to a set of experimental atomization energies, ionization potentials, proton affinities, and total atomic energies.

Density Functional Theory

B3LYP LDA HF LDA GGA LDA GGA LDA 0(

) ( ) ( ),

xc xc x x x x x c c c

E E a E E a E E a E E

slide-22
SLIDE 22

Current implementations of density functional theory are semi-empirical. They are very accurate in most situations, but tend to break down for weak interactions. There is an overriding need for advances in understanding of the nature of the exact energy functional.

Density Functional Theory

slide-23
SLIDE 23

The exact formula is:

?

Density Functional Theory