2/18/2015 1 2/18/2015 PHY 752 Spring 2015 -- Lecture 15 1 PHY 752 Solid State Physics 11-11:50 AM MWF Olin 107 Plan for Lecture 15: Reading: Chapter 10 in MPM Numerical Realizations of Density functional theory 1. Electronic structure of atoms 2. Integration of the radial equations 3. Frozen core approximation 4. Extension of formalism to multi-center analysis 2/18/2015 PHY 752 Spring 2015 -- Lecture 15 2 Note: Take-home exam scheduled for the week of March 2 nd . 2/18/2015 PHY 752 Spring 2015 -- Lecture 15 3
8
Embed
2/18/2015users.wfu.edu/natalie/s15phy752/lecturenote/Lecture15for...2/18/2015 1 2/18/2015 PHY 752 Spring 2015 --Lecture 15 1 PHY 752 Solid State Physics 11-11:50 AM MWF Olin 107 Plan
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
2/18/2015
1
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 1
PHY 752 Solid State Physics11-11:50 AM MWF Olin 107
Plan for Lecture 15:Reading: Chapter 10 in MPM
Numerical Realizations of Density functional theory
1. Electronic structure of atoms
2. Integration of the radial equations
3. Frozen core approximation
4. Extension of formalism to multi-center analysis
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 2
Note: Take-home exam scheduled for the week of March 2nd.
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 3
2/18/2015
2
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 4
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 5
Kohn-Sham equations for spherical atom
2
2 2
( 1)( ) ( ) ( ) ( ) ( )
Equations in Rydberg units
i i i i i i
ie
ie exc n lnl ln
l ldV r V r v r P r P r
dr r
2
2 1 3
0
/1/3
[ ]( ) ( ) ( )
[ ] 2( ) ( ) (
12 ' '
)
[ ] 2( )
3
( )
eeee
excexc c
extex
r
t
r
E nr dr n r dr n r
n
E nr n r r
n
E n Zr v r
V r rr
V V
Vn r
p
p
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 6
Self-consistent solution
Iteration 0
2
( )
( ) ( )
i
ii
n
r
r r
1 12
2 ( ) ( ) ( ) ( ) ( )2 i iee ex iV V vm
r r r r r
\alpha
21 1
1 1
( ) ( )
( ) ( ) (1 ) ( )
temp
temp
ii
n
nn x x n
r r
r r r
1
Numerical methods for solving the Kohn-Sham equations
2/18/2015
3
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 7
Example for carbon
2
22 2
22 2
2
1 2 2
1 2 2
=4 2 2 2
4 =
( ) ( )
( ) ( )
2
( )
( ) (2 ) ( )2
i i i i
s s p
s s
n l
p
l ni
n r w r
r r r
P r P rr
P r
p
p
ˆ(
For
)
spherically symmetric a
( ) ( )
( )
:
( )
tom
i i i i i i i
i i
i i
n l m n l l m
n l
n l
r Y
P rr
r
r r
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 8
Results for carbon
r (Bohr)
Pnl(r
)
1s
2s
2p
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 9
ncore(r)
nval(r)
r (Bohr)
nr2
/4p
2/18/2015
4
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 10
Example for Cu
1s22s22p63s23p63d104s1
r (Bohr)
ncore(r)
nvale(r)
nr2
/4p
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 11
r (Bohr)
4s
4p3d
3s
3p
Results for Copper
Pnl(r
)
2/18/2015 PHY 752 Spring 2015 -- Lecture 15 12
Frozen core approximation
Example for Cu
1s22s22p63s23p63d104s1
r (Bohr)
ncore(r)
nvale(r)
nr2
/4p
core vale
vale
( ) ( )
Variationally optimize energy wrt
)
(
(
)
r n r
n r
n r n
2/18/2015
5
2/18/2015 PHY 752 Spring 2015 -- Lecture 15b 13
Systematic study of frozen core approximation in DFT