prospects of lda+dmft - upv/ehunano-bio.ehu.es/files/prospects_of_lda+dmft-biermann.pdf ·...
TRANSCRIPT
![Page 1: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/1.jpg)
Prospects of LDA+DMFT
Silke BiermannCentre de Physique Theorique
Ecole Polytechnique, Palaiseau, France
(*) LDA = the local density approximation (LDA) of density
functional theory
LDA+DMFT = the combination of dynamical mean field theory
(DMFT) with the LDA. – p.1/69
![Page 2: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/2.jpg)
The Mott transition within DMFTSpectral function ofthe Hubbard model
2
0
2
0
2
0
2
0
2
0
−Im
G
ω−4 −2 0 2 4
U/D=1
U/D=2
U/D=2.5
U/D=3
U/D=4
/D
Fig.30
. – p.2/69
![Page 3: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/3.jpg)
Spectral function – survival kitAdd/remove an electron – at which energy?Two “easy” limiting cases:
1. Non-interacting limit:state of N electrons = Slater determinant(N+1)th electron can jump into any (unoccupied) band
probe unoccupied density of states
0
0.2
0.4
0.6
−2 −1 0 1 2E−E
Fermi
DOS
. – p.3/69
![Page 4: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/4.jpg)
Spectral function – survival kit2. “Atomic limit” (complete localization):probe local Coulomb interaction!
0
0.1
0.2
0.3
0.4
0.5
−10 −5 0 5 10E−E
Fermi
U
In the general, interacting case:Spectral function
describes the possibility ofadding an electron with energy (includingrelaxation effects)
. – p.4/69
![Page 5: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/5.jpg)
The Mott transition within DMFTSpectral function ofthe Hubbard model
2
0
2
0
2
0
2
0
2
0
−Im
G
ω−4 −2 0 2 4
U/D=1
U/D=2
U/D=2.5
U/D=3
U/D=4
/D
Fig.30
. – p.5/69
![Page 6: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/6.jpg)
Mott insulator and correl. metal:YTiO and SrVO
Inte
nsity
(ar
b. u
nits
)
2.0 1.0 0.0Binding Energy (eV)
Sr0.5Ca0.5VO3
SrVO3 (x = 0)
CaVO3 (x = 1)
hν = 900 eV hν = 275 eV hν = 40.8 eV hν = 21.2 eV
2
0
2
0
2
0
2
0
2
0−
ImG
ω−4 −2 0 2 4
U/D=1
U/D=2
U/D=2.5
U/D=3
U/D=4
/D
Fig.30
. – p.6/69
![Page 7: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/7.jpg)
Outline
Reminder: Dynamical Mean Field Theory(DMFT)
The “LDA+DMFT” method
:Reminder: an “effective atom” point of view
A functional point of view
Examples
What can we calculate ?The current status
Beyond LDA+DMFT: the GW+DMFT
scheme
LDA = The local density approximation to Density Functional Theory
LDA+DMFT = The combination of LDA and dynamical mean field theory (DMFT)
(**) GW+DMFT = The combination of Hedin’s GW approximation with DMFT. – p.7/69
![Page 8: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/8.jpg)
The local Green’s function ...... is the central object of DMFT
Definition of Green’s function:
Relation to local spectral function:
. – p.8/69
![Page 9: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/9.jpg)
Green’s function – survival kitQuasi-particles are poles of
All correlations are hidden in the self-energy:
. – p.9/69
![Page 10: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/10.jpg)
Dynamical mean field theory ...... maps a lattice problem
onto a single-site (Anderson impurity) problem
with a self-consistency condition
. – p.10/69
![Page 11: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/11.jpg)
Effective dynamics ...... for single-site problem
with the dynamical mean field
. – p.11/69
![Page 12: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/12.jpg)
DMFT (contd.)Green’s function:
Self-energy (k-independent):
DMFT assumption :
Self-consistency condition for
. – p.12/69
![Page 13: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/13.jpg)
The DMFT self-consistency cycleAnderson impurity model solver
Self-consistency condition:
. – p.13/69
![Page 14: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/14.jpg)
Realistic Approach to CorrelationsCombine DMFT with band structure calculations
(Anisimov et al. 1997, Lichtenstein et al. 1998)
effective one-particle Hamiltonian within LDArepresent in localized basis
add Hubbard interaction term for correlated orbitals
solve within Dynamical Mean Field Theory
. – p.14/69
![Page 15: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/15.jpg)
Hamiltonian formulation
(correl. orb.)
(correl. orb.)
in localized basis set( lecture by F. Lechermann)
all valence electrons
Hubbard terms for “correlated orbitals”. – p.15/69
![Page 16: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/16.jpg)
Dynamical mean field theory ...... maps a solid
onto an “effective atom” problem
with a self-consistency condition
. – p.16/69
![Page 17: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/17.jpg)
The DMFT self-consistency cycleAnderson impurity model solver
Self-consistency condition:
. – p.17/69
![Page 18: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/18.jpg)
What do we mean by this?Represent the Green’s function in localized basis,e.g. LMTO’s:
,
where
double counting correctionsis a matrix in orbital space
is a matrix in the space of the correlated orbitals
lecture by F. Lechermann
. – p.18/69
![Page 19: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/19.jpg)
Hubbard interaction termsScreened Coulomb interaction, onsite terms only:
with
Simplification: keep “diagonal density” terms only:
. – p.19/69
![Page 20: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/20.jpg)
Parametrization of interaction:
with the Slater integrals
, e.g. for d-electrons:
. – p.20/69
![Page 21: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/21.jpg)
Double counting corrections
of Hartree type (cf. “LDA+U”, Anisimov et al.)
Lichtenstein et al., 2001:
Tr
( no zero frequency contribution,cf. LDA Fermi surface Luttinger’s theorem)
. – p.21/69
![Page 22: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/22.jpg)
A functional point of view ?Reminder about Density Functional Theory:Energy is a functional of the density:E= E[n(r)]
Kohn-Sham potential can be viewed as a Lagrangemultiplier to fix the density of the non-interactingreference system
. – p.22/69
![Page 23: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/23.jpg)
A functional point of view ...... of LDA+DMFT
“Spectral Density Functional Theory”:Free energy a functional of(1) the density n(r)(2) the local Green’s function of the correlated orbitals
Construction by introduction of source termslecture by M. Katsnelson
Kotliar, Savrasov, PRB 2004
. – p.23/69
![Page 24: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/24.jpg)
LDA+DMFT as a ...... spectral density functional theoryFree energy functional:
tr
tr
Kotliar, Savrasov, 2004 . – p.24/69
![Page 25: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/25.jpg)
LDA+DMFT – the full scheme
DMFT loop
DMFT preludeDFT part
update
PSfrag replacements
VKS = Vext + VH + Vxc
[
−∇2
2 + VKS
]
|ψkν〉 = εkν |ψkν〉
from charge density ρ(r) constructupdate
|χRm
〉
GKS =[
iωn + µ+ ∇2
2 − VKS
]−1
G0
build GKS =[
iωn + µ+ ∇2
2 − VKS
]−1
construct initial G0
impurity solver
Gimpmm′(τ − τ ′) = −〈T dmσ(τ)d†
m′σ′(τ ′)〉Simp
self-consistency condition: construct Gloc
G−10 = G−1
loc + Σimp
Gloc = P(C)R
[
G−1KS −
(
Σimp − Σdc
)]−1
P(C)R
Σimp = G−10 − G−1
imp
ρ
compute new chemical potential µ
ρ(r) = ρKS(r) + ∆ρ(r)
(Appendix A)
(from F. Lechermann, A. Georges, A. Poteryaev, S. B., M. Posternak, A. Yamasaki, O. K. Ander-
sen, Phys. Rev. B 74 125120 (2006))
. – p.25/69
![Page 26: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/26.jpg)
Some more examples
. – p.26/69
![Page 27: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/27.jpg)
Cerium sesquioxide Ce O Mott insulator,
paramagnetic above 10 K
. – p.27/69
![Page 28: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/28.jpg)
Cerium sesquioxide Ce O
Metallic in LDA:
-6
-4
-2
0
2
4
6
8
10
G K M G A L H A
Ene
rgy
(eV
)Ce2O3 bands (Ce 6p in basis)
Note the narrow f-bands!. – p.28/69
![Page 29: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/29.jpg)
Cerium sesquioxide Ce O
LDA+DMFT splits f-states into upper and lowerHubbard bands:
(L. Pourovskii et al., arXiv:0705.2161). – p.29/69
![Page 30: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/30.jpg)
Spectral functionCe O
25
50
0
25
50
Den
sity
of
stat
es (
1/eV
)
LDA
DMFT-nonSC
DMFT-SC
-6 -4 -2 0 2 4Energy (eV)
0
25
50
(L. Pourovskii et al., arXiv:0705.2161)
. – p.30/69
![Page 31: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/31.jpg)
Total energyCe O
3.6 3.7 3.8 3.9 4
a (Ao
)
0
25
50
75
100
E-E
min
(m
Ry)
LDA DMFT (without charge SC)DMFT (with charge SC)
aexp
=3.89 Ao
(L. Pourovskii et al., arXiv:0705.2161)
. – p.31/69
![Page 32: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/32.jpg)
Vanadium dioxide: VO
High temperature rutile phase: metallic
Low temperature monoclinic phase: insulating
. – p.32/69
![Page 33: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/33.jpg)
Metal-insulator Transition
. – p.33/69
![Page 34: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/34.jpg)
VO in LDA and LDA+DMFT
0
0.5
1.0
−2 0 2 4
Spe
ctra
lFun
ctio
n(a
rb.u
nits
)
ω (eV)
VO2
rutile phase
VO2
monoclinic phase
LDA: dashed lineLDA+DMFT: solid line
S.B., A.Poteryaev, A. Georges, A.Lichtenstein, PRL 2005 . – p.34/69
![Page 35: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/35.jpg)
LDA+DMFT simulations ...for the monoclinic phase have to be done with care:
Inter-site (intra-pair)fluctuations ?
"Cluster-DMFT" :embed V pair in bath non-local self-energy !
. – p.35/69
![Page 36: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/36.jpg)
VO recent photoemission
Koethe et al.,PRL 2006
lecture byL.-H. Tjeng
. – p.36/69
![Page 37: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/37.jpg)
What can we calculate ?
. – p.37/69
![Page 38: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/38.jpg)
What can we calculate ?
total energy
Examples: Ce O Pu
3.6 3.7 3.8 3.9 4
a (Ao
)
0
25
50
75
100
E-E
min
(m
Ry)
LDA DMFT (without charge SC)DMFT (with charge SC)
aexp
=3.89 Ao
(Pourovskii et al., arXiv:0705.2161) (Savrasov et al., Nature 2001)
(see also Held et al., PRL 2001, and Amadon et al., PRL 2006, for total energy calculations for
Cerium)
. – p.38/69
![Page 39: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/39.jpg)
Total energy functional
with:
= crystal energy
= Hartree energy
= exchange correlation energysecond line = energy terms stemming from manybody Hamiltonian
(from B. Amadon et al., PRL 2006)
. – p.39/69
![Page 40: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/40.jpg)
What can we calculate ?
local spectral functions
Examples:VO2-M1 and Ce O
25
50
0
25
50
Den
sity
of
stat
es (
1/eV
)
LDA
DMFT-nonSC
DMFT-SC
-6 -4 -2 0 2 4Energy (eV)
0
25
50
0
0.5
1.0
−2 0 2 4
Spe
ctra
lFun
ctio
n(a
rb.u
nits
)
ω (eV)
VO2
rutile phase
VO2
monoclinic phase
. – p.40/69
![Page 41: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/41.jpg)
What can we calculate ?
total energy
local spectral functions
self-energies
-2
-1
0
1
-3 -2 -1 0 1 2 3
Σ(ω
)
ω [eV]
lecture by D. Vollhardt
. – p.41/69
![Page 42: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/42.jpg)
Self-energy and spectral function
0
0.1
0.2
0.3
0.4
0.5
0.6
-3 -2 -1 0 1 2 3
A(ω
)
ω [eV]
-2
-1
0
1
-3 -2 -1 0 1 2 3
Σ(ω
)
ω [eV]
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
-4 -2 0 2 4
A(ω
)
ω [eV]
-5
-4
-3
-2
-1
0
1
2
3
-4 -2 0 2 4
Σ(ω
)ω [eV]
(from J.Tomczak, PhD thesis, 2007) . – p.42/69
![Page 43: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/43.jpg)
What can we calculate ?
self-energies
Example: VO insulating phase
-3
-2
-1
0
1
2
3
4
-3 -2 -1 0 1 2 3 4
ω [eV]
Re
Σ (ω
) [e
V]
a1g
a1g-a1g
egπ1
egπ2
-3
-2
-1
0
1
2
3
4
-3 -2 -1 0 1 2 3 4
ω [eV]Im
Σ (
ω)
[eV
]
. – p.43/69
![Page 44: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/44.jpg)
What can we calculate ?
... and k-resolved spectral functions
Example: VO insulating and metallic phases:
0.1
1
ω [e
V]
ΓZCYΓ-2
-1
0
1
2
3
0.1
1
ω [e
V]
ΓZCYΓ-2
-1
0
1
2
3
effective bandstructures (and absence thereof!)
J.Tomczak, F. Aryasetiawan, SB, arXiv:0704.0902 . – p.44/69
![Page 45: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/45.jpg)
What can we calculate ?
phonons
Example:
-Pu
Dai et al., Science 2003
. – p.45/69
![Page 46: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/46.jpg)
What can we calculate ?
local susceptibilities
0.0 0.5 1.0 1.5 2.00.0
0.5
1.0
0.0
0.5
1.0
Ni
Fe
χ-1M
eff2 /3
Tc
M(T
)/M
(0)
T/Tc
A. Lichtenstein, M. Katsnelson, G. Kotliar,
PRL 2001
. – p.46/69
![Page 47: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/47.jpg)
What can we calculate ?
total energy
local spectral functions
k-resolved spectral functions
phonons
local susceptibilities
optical properties
. – p.47/69
![Page 48: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/48.jpg)
Optical conductivityIn the Hubbard model within DMFT
from Rozenberg et al., 1995
d
f f
tr
k
k
k
k
see also: Pruschke et al., Blümer PhD thesis, Oudovenko et al., Haule et al.
. – p.48/69
![Page 49: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/49.jpg)
Optical conductivityThe example of insulating VO :
from J.M. Tomczak, PhD thesis, 2007(for full formalism Poster by J. Tomczak)
0
1
2
3
4
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5
Re
σ(ω
) [1
03 (Ω
cm)-1
]
ω [eV]
Experiment LDA+CDMFT ’upfolded’ E || [100]E || [010]E || [x0z]
. – p.49/69
![Page 50: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/50.jpg)
LDA+DMFT – current statusMany successful applications:
Correlation effects in transition metals
Mott transitions in oxides
Volume collapse transitions in rare earth/actinidesystems
etc ...
Still (at least partly) open questions/challenges:
influence of full charge self-consistency
basis sets lecture by F. Lechermann
what’s U ? (Is there a U? On what energy scale?)
... and nuisances: double counting. – p.50/69
![Page 51: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/51.jpg)
What’s U in a solid?... an answer from RPA:
Divide where = polarization of thecorrelated orbitals (e.g. 3d orbitals)Then:
where that does not include 3d-3d screening:
Identify = !F. Aryasetiawan, M. Imada, A. Georges, G. Kotliar, S.B., A. I. Lichtenstein PRB 70 195104 (2004)
. – p.51/69
![Page 52: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/52.jpg)
U and W in a solidScreened Coulomb interaction from RPA
and
for Nickel:
0 10 20 30 40 50 60
ω (eV)
0
10
20
30
40
Re
<dd|
W|d
d> (e
V)
W Eg
Wr E
g
Paramagnetic Ni
(Aryasetiawanet al., 2004)
. – p.52/69
![Page 53: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/53.jpg)
What’s U in a solid?
... what about an answer beyond RPA ??
. – p.53/69
![Page 54: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/54.jpg)
Can we calculate ...... from a (dynamical) impurity model?
Question of representability !
DMFT: calculated from impurity model
What about ?
Self-consistency requirement:
= of the solid
= of the solid
“GW+DMFT”
S.B., Aryasetiawan, Georges, PRL 2003
condmat
. – p.54/69
![Page 55: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/55.jpg)
The GW approximation(Hedin, 1965)
dynamically screened Coulomb interaction
GW successful for sp-metals, semiconductors ...
(Reviews:Onida et al., Rev. Mod. Phys. 2002;
Aryasetiawan et al., Rep. Prog. Phys. 1998)
. – p.55/69
![Page 56: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/56.jpg)
A functional point of view[Almbladh et al. 1999]
Tr
Tr
Tr
Tr
Free energy
is a functional of
one-electron Green’s function
the screened Coulomb interaction
= bare (Hartree) Green’s function
. – p.56/69
![Page 57: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/57.jpg)
Approximations to
?GW:
Extended DMFT (“E-DMFT”):
E-DMFT for local part + GW for nonlocal part:
NB: “local” = “onsite” is a basis-set dependent notion!
. – p.57/69
![Page 58: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/58.jpg)
Extended DMFT ...... maps a lattice problem
onto a single-site (Anderson impurity) problem
with a dynamical interaction
. – p.58/69
![Page 59: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/59.jpg)
Self-consistency loop
. – p.59/69
![Page 60: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/60.jpg)
Challenges and questions
Global self-consistency?
Choice of orbitals? Hamiltonian?
Treat all orbitals – localized and delocalized – onequal footing? Downfolding?
How to solve the dynamical impurity model?here: static approximation
. – p.60/69
![Page 61: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/61.jpg)
A simplified implementationNon-selfconsistent GW + local from static impuritymodel
Tr
Nonlocal part: correct Hartree by GW
Local part: correct LDA by DMFT
. – p.61/69
![Page 62: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/62.jpg)
Simplified GW+DMFTNi band structure
−8
−6
−4
−2
0
Γ
Ene
rgy
[eV
]
XX Minority Majority
Circles: GW+DMFTDashed: LDATriangles: photoemission data
(Bünemann et al. 2002, Mårtensson et al., 1984)
. – p.62/69
![Page 63: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/63.jpg)
Simplified GW+DMFT:Spectral function of Ni
0
2
4
6
−10 −5 0 50
2
4
6
−10 −5 0 5E−EF
[eV]
ρ(E)
6eV
Majority and minority spins
Satellite at 6eV correct!
. – p.63/69
![Page 64: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/64.jpg)
Conclusion and perspectivesCombination of LDA and DMFT ...
a useful tool for electronic structure calculationsof strongly correlated materials, such as
transition metals, their oxides, f-electron materials
VO (Correlation-induced Peierls transition)
Ce O , a Mott insulator
Prospects ?
technical and physical questions concerningimplementations
calculate more quantities, more materials ... !
Beyond LDA+DMFT ? GW+DMFT
Plenty of work left! . – p.64/69
![Page 65: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/65.jpg)
References of recent work
VO : Correlation-induced Peierls transitionS.B., A. Poteryaev, A. Georges, A. Lichtenstein, Phys. Rev. Lett. 94, 026404 (2005)
J. M. Tomczak, SB, Journal of Phys. Cond. Mat., in press
J.M. Tomczak, F. Aryasetiawan, SB, condmat07040902
Ce O : A Mott insulator – questions ofcharge-selfconsistencyL. Pourovskii, B. Amadon, S. Bierman, A. Georges arXiv:0705.2161
U and GW+DMFT:F. Aryasetiawan, M. Imada, A. Georges, G. Kotliar, S.B., A. I. Lichtenstein PRB 70195104 (2004),
S. Biermann, F. Aryasetiawan, A. Georges, PRL 2003, condmat 2004
LDA+DMFT - recent Reviews:D. Vollhardt, G. Kotliar, Physics Today 2004
S. B., in Encyclop. of Mat. Science. and Technol., Elsevier 2005
. – p.65/69
![Page 66: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/66.jpg)
GW+DMFT: local part
calculated from local impuritymodel:
. – p.66/69
![Page 67: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/67.jpg)
GW+DMFT (contd)Combine local self-energy and polarization
with non-local self-energy and polarization:
. – p.67/69
![Page 68: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/68.jpg)
Self-consistency condition
Update Weiss field and impurity interaction:
Iterate until self-consistency ...
. – p.68/69
![Page 69: Prospects of LDA+DMFT - UPV/EHUnano-bio.ehu.es/files/Prospects_of_LDA+DMFT-Biermann.pdf · Prospects of LDA+DMFT Silke Biermann Centre de Physique Theorique´ Ecole Polytechnique,](https://reader033.vdocument.in/reader033/viewer/2022042918/5f5d4dd4ef6c3a5a002bfa21/html5/thumbnails/69.jpg)
NickelEvidence for many-body effects:
Exp. LDA GW
LDA+DMFT
bandwidth 3.3 eV 4.5 eV 3.5 eV OK!x-splitting 0.3 eV 0.6 eV 0.6 eV OK!satellite at 6eV NO! NO! YES!
large quasi-particle widthsOpen question: Fermi surface pocket ?
(1) Aryasetiawan, 1992(2) Lichtenstein et al., 2001
. – p.69/69