Docsity
Docsity

Prepare-se para as provas
Prepare-se para as provas

Estude fácil! Tem muito documento disponível na Docsity


Ganhe pontos para baixar
Ganhe pontos para baixar

Ganhe pontos ajudando outros esrudantes ou compre um plano Premium


Guias e Dicas
Guias e Dicas


Teoria de Densidade Funcional (DFT) em Física Quântica: Aplicações do Hamiltoniano Eletrôn, Notas de estudo de Engenharia de Produção

Um problema relacionado à teoria do hamiltoniano geral de matéria condensada, especificamente sobre a aproximação de born-oppenheimer para a separação de variáveis electrônicas e iônicas, e as limitações desta abordagem. O texto também discute a teoria de kohn-sham e a definição de um funcional universal (f) com base no hamiltoniano de kohn-sham. Além disso, são discutidos os equações de kohn-sham, o princípio variacional e a definição de um estado de densidade de elétrons estacionário.

Tipologia: Notas de estudo

Antes de 2010

Compartilhado em 09/11/2009

igor-donini-9
igor-donini-9 🇧🇷

4.5

(4)

419 documentos

1 / 22

Toggle sidebar

Esta página não é visível na pré-visualização

Não perca as partes importantes!

bg1
A fast guide
to density functional calculations
P. Kratzer
Fritz-Haber-Institut der MPG
D-14195 Berlin-Dahlem, Germany
I. From the many-particle problem to the
Kohn-Sham functional
II. How to perform a total energy calculation
III. From the total energy to materials science
pf3
pf4
pf5
pf8
pf9
pfa
pfd
pfe
pff
pf12
pf13
pf14
pf15
pf16

Pré-visualização parcial do texto

Baixe Teoria de Densidade Funcional (DFT) em Física Quântica: Aplicações do Hamiltoniano Eletrôn e outras Notas de estudo em PDF para Engenharia de Produção, somente na Docsity!

A fast guide

to density functional calculations

P Kratzer

FritzHaberInstitut der MPG

D BerlinDahlem Germany

I From the manyparticle problem to the

KohnSham functional

I I How to p erform a total energy calculation

I I I From the total energy to materials science

General condensedmatter Hamiltonian

T

e 

T

ion 

V

ee 

V

eion 

V

ionion   E 

wavefunction r      r N

 R

I

    R

M

electronic co ordinates r k  k      N

ionic co ordintes R I

 I      M

 T

e 

N X

k 

p



k

m

 T

ion 

M X

I 

P



I

MI

 V

ee 



  

N N X

k  k



e



jr k  r k  j

 V

ionion 



 

M M X

I I 

ZI Z I 

jRI  R I

 j

 V

eion rk  RI  

N X

k

M X

I

v

ion

I

jRI  rk j

Limitations of this approach

 do esnt account for correlated dynamics of ionic and

electronic co ordinates

Example

supra uid He  p olaroninduced sup erconductivity

 breakdown of the restriction to a single groundstate

BornOpp enheimer surface

Example chemoluminescence

0 1 2

distance

0

1

energy

 breakdown of the adiabatic approximation

Example

excitation of surface plasmons during scattering of

an ion from a metal surface

 timedep endent theories eg TDDFT

electronic

manypa

rticle

Hamiltonian

 

N

k X 

r

k

m

v (^) 

^

r

k 

N

(^) N

X k  

k   

(^) W

r

k  r

k   

r      r N   E

r      r N 

W  r  r 

^

e 



(^) 

 j r  r 

j

v

(^) 

r 

M

I X (^) 

v

I (^) ion

j R I  r j 

still

many

fo

r

a

t ypical

solid



degrees

of

freedom

The

manypa

rticle

p

roblem

can

b

(^) e

solved

only

fo

r

rather

small

systems

atoms

molecules

and

clusters

using

established

metho

(^) ds

of

quantum

chemistry

eg

conguration

interaction

KohnSham single particle Hamiltonian

To nd the stationary p oint we do variations at xed

N 

R

dr n r  which leads to

 E

v



 n r

  Lagrange parameter

If we write the density as a sum over singleparticle

functions

n r 

N

X

j

X

kBZ

j jk r j

 

the variational principle  E v

  

 r   leads

to the KohnSham equations

r



m

 V

e

n r

jk r   jk

jk r

with the eective p otential

V

e

n r  v



R I

r 

Z

dr



e

 n r



 jr  r

 j

 E

X C n

 n

r 

The

total

energy

fo

r

static

ions

E

tot

n   T   n   Z d r v

R (^) 

i  r  n  r   

 Z d r Z d r  e 

n^

 r  n  r 

j r



r  j^

E

X

(^) C

n

V

ion



ion

R

I

E

tot

n

N

j X (^) 

X

k 

BZ

j

k

E

(^) e 

e ^ n





E

X

(^) C

n

V

ion



ion

R

I

E

(^) e 

e ^ n





 Z d r Z d r  e 

n^

 r  n  r 

^

j r



r  j^



E

(^) e 

e ^ n



E

X

(^) C

n



E

X

(^) C

n   Z d r V X

(^) C

n



r  n  r 

E

tot

is

stationa

ry

with

resp

to

va

riations

of

n

a round

n

 (^) 

but

the

individual

terms

a re

not



P

erio

dic

systems

crystal

structure

Bravais

lattice

basis

translational

symmetry

symmetries

of

the

basis

compatible

with

the

Bravais

lattice

group

T

p

(^) oint

group

P

Example

t w

o

(^) dimensional

honeycomb

lattice

Bravais

lattice

if

a

basis

of

t w

o

atoms

is

used

a 

a 

b 

b^  

Brillouin

zone

BZ

lattice

vecto

rs

a

 (^) 

a

 

a



recip

ro

(^) cal

lattice

vecto

rs

b  

b  

b 

b^  i



a

j



a

k (^) 

f

ijk

g



f



g

cyclic

p

(^) ermutations

of

indices

Blo ch functions

Since H commutes with the elements of T  the wave

functions must b e of the form

 jk r  R  e

ikR  jk r Blo chs theorem

for any R  n  a   n  a   n  a 

The wavevector k is an index sp ecied by a p oint in the

rst Brillouin zone elementary cell of the recipro cal

lattice and the numb er of such p oints is equal to the

numb er of lattice sites in the crystal

Due to Blo chs theorem it suces to calculate  in

just one elementary cell

 jk r  e

ikr u jk r

where u is a latticep erio dic function

In practice sums over k are evaluated by summing

over a discrete set of sp ecial kp oints

Selfconsistency

cycle

Ha

rris

functional

E

tot

n  i    N

j X (^) 

X

k 

BZ

j

k

E

(^) e 

e ^ n

 i 



E

X

(^) C

n

 i 



V

ion



ion

R

I

up

date

cha

rge

densit

y

n

 i  ^ r



N

k

X

k

n

B

X

j

f

j

k (^) j 

 i 

j

k  r  j 

with

f

j

k

exp

j

k

F

k

T



 

In

p

rinciple

w

e

a

re

interested

in

T



but

nite

T

helps

to

stabilize

the

SC

lo

(^) op

E

tot

n

 i 

  T   n  i 

^

Z

d

r

v

R (^) 

I  r  n  i   r 

E

(^) e 

e ^ n

 i   r



E

X

(^) C

n

 i  ^



V

ion



ion

R

I

KohnSham

functional

E

tot

 n    E

tot

n



(^) 

with

selfconsistent

densit

y

n



V

a

riational

p

rop

erties

A

t

the

stationa

ry

p

(^) oint

the

KohnSham

equations

lead

us

to

T

 (^)  n





N

j X (^) 

X

k 

BZ

j

k  Z d r n  r  V

e

n



r



F

o

r

the

nonconsistent

case

w

e

ma

y

intro

duce

the

generalization

T

 (^)  n

V

e

N

j X (^) 

X

k 

BZ

j

k  V

e

  Z d r n  r  V

e

r



This

can

b

(^) e

used

to

dene

the

double

functional

E

tot (^) D

n

V

e

T

 (^)  n

V

e

E

(^) e  e  n   E X

(^) C

n

V

ion



ion

R

I

with

the

va

riational

p

rop

erties

E

tot (^) D

n



(^) n

V

e

n

 (^) 

E

tot (^) D

n

 (^) 

V

e

n

 (^) 

c  (^)   (^) n



E

tot (^) D

n

  V

e

n

 (^) 



(^) v

(^) 



E

tot (^) D

n

 (^) 

V

e

n

 (^) 

c  (^)   (^) v

  

No

rmally



one

has

c    c  

The

Ha

rris

functional

in

this

notation

is

E

tot (^) D

n

 i 



V

e

n

 i 



^



Example Thermal decay of GaAs

GaAs  Ga metal

As 

G

GaAs bulk

Ga

p T   G GaAs

T   G

As p T 

G

As p T  from ideal twoatomic gas S

ion

GaAs

and S

ion

Ga

from Debye mo del of the solid

500 600 700 800 900 1000

temperature [K]

−3.

−3.

−2.

g

Ga

[eV/particle]

μGa(GaAs, p(As 2 )=

− Pa)

μ Ga (GaAs, p(As 2 )=

− Pa)

μ Ga (bulk Ga)

Free enthalpy p er Ga atom in bulk GaAs in

thermo dynamic equilibrium with As  vap or at two

pressures compared to the free enthalpy p er Ga atom

in elemental Ga

P Kratzer et al Phys Rev B   

Finite electronic temp erature

o ccupation numb ers from Fermi distribution

f jk  exp jk

F k T   



 The program actually computes the free energy F

entropy of the electronic system

S

e  k B

n B

X

j 

X

kBZ

f jk lnf jk    f jk ln  f jk

For a freeelectronlike metal we have

F T   E T   
T

  O T

 

E T   E T   
T

  O T

 

extrap olation to zero electronic temp erature

E T    F T   E T   E T   S

e T 

see M J Gillan J Phys Cond Mat   

C second derivatives

Examples

 force constant matrix

 E

tot

R

i

R

j

 calculation of phonon sp ectrum vibrational

entropy   

 particle numb er uctuations

 chemical softness and hardness

sr 

 nr

v T

Z ^

 E tot

nr nr

 



dr



hr 

 nr

v T

Z

 E tot

nr nr

 

nr

 

N

dr



Note When calculating second derivatives the

resp onse of the density must b e taken into account

 Density Functional Perturbation Theory

What materials prop erties are accessable

to calculation 

 structural prop erties

Examples structural phase transitions surface

reconstructions yes

 elastic prop erties

Examples bulk mo dulus C 

 C



 C

     yes

 chemical prop erties

Examples thermo chemical stability of comp ounds

reactivity of surfaces yes

 transp ort prop erties

Examples transp ort eective mass magneto

resistence developing eld

 opticalsp ectroscopic prop erties

Examples photo emission sp ectra cross sections for

light absorption

topic b eyond KohnSham theory encouraging

progress recently

many other applications