Outline. GW approximation. Electrons in solids. The Green Function. Total energy---well solved Single particle excitation---under developing

Similar documents
Modern Energy Functional for Nuclei and Nuclear Matter. By: Alberto Hinojosa, Texas A&M University REU Cyclotron 2008 Mentor: Dr.

ESS 265 Spring Quarter 2005 Kinetic Simulations

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

MCTDH Approach to Strong Field Dynamics

Time-Dependent Density Functional Theory in Condensed Matter Physics

s = rθ Chapter 10: Rotation 10.1: What is physics?

Chapter 3: Vectors and Two-Dimensional Motion

Field due to a collection of N discrete point charges: r is in the direction from

1 Constant Real Rate C 1

A. Inventory model. Why are we interested in it? What do we really study in such cases.

Determining Modern Energy Functional for Nuclei And The Status of The Equation of State of Nuclear Matter. Shalom Shlomo Texas A&M University

Teoría del Funcional de la Densidad (Density Functional Theory)

5-1. We apply Newton s second law (specifically, Eq. 5-2). F = ma = ma sin 20.0 = 1.0 kg 2.00 m/s sin 20.0 = 0.684N. ( ) ( )

Mechanics Physics 151

Mechanics Physics 151

Bethe-Salpeter Equation Green s Function and the Bethe-Salpeter Equation for Effective Interaction in the Ladder Approximation

H = d d q 1 d d q N d d p 1 d d p N exp

CIRCUITS AND ELECTRONICS. The Impedance Model

( ) () we define the interaction representation by the unitary transformation () = ()

Mechanics Physics 151

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

Maximum Likelihood Estimation

FI 3103 Quantum Physics

Reinforcement learning

Chapter 6: AC Circuits

On The Estimation of Two Missing Values in Randomized Complete Block Designs

Course Outline. 1. MATLAB tutorial 2. Motion of systems that can be idealized as particles

Quasiparticle Band Structures and the GW Approximation

Representing Knowledge. CS 188: Artificial Intelligence Fall Properties of BNs. Independence? Reachability (the Bayes Ball) Example

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

Linear Response Theory: The connection between QFT and experiments

Solution in semi infinite diffusion couples (error function analysis)

CS 188: Artificial Intelligence Fall Probabilistic Models

The Unique Solution of Stochastic Differential Equations. Dietrich Ryter. Midartweg 3 CH-4500 Solothurn Switzerland

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Let s treat the problem of the response of a system to an applied external force. Again,

Relative and Circular Motion

Density Functional Theory I

Chapters 2 Kinematics. Position, Distance, Displacement

Physics 11b Lecture #2. Electric Field Electric Flux Gauss s Law

Name of the Student:

Cooling of a hot metal forging. , dt dt

Lecture 5. Plane Wave Reflection and Transmission

Chapter Finite Difference Method for Ordinary Differential Equations

Muffin Tins, Green s Functions, and Nanoscale Transport [ ] Derek Stewart CNF Fall Workshop Cooking Lesson #1

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

Numerical Study of Large-area Anti-Resonant Reflecting Optical Waveguide (ARROW) Vertical-Cavity Semiconductor Optical Amplifiers (VCSOAs)

calculating electromagnetic

EECE 301 Signals & Systems Prof. Mark Fowler

Exponential and Logarithmic Equations and Properties of Logarithms. Properties. Properties. log. Exponential. Logarithmic.

Physics 201 Lecture 15

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Computer Propagation Analysis Tools

TSS = SST + SSE An orthogonal partition of the total SS

May 29, 2018, 8:45~10:15 IB011 Advanced Lecture on Semiconductor Electronics #7

Time-Dependent Density Functional Theory in Optical Sciences

Detection and Estimation Theory

Molecular dynamics modeling of thermal and mechanical properties

EQUATIONS OF MOTION LUCA GUIDO MOLINARI

P a g e 5 1 of R e p o r t P B 4 / 0 9

Density Matrix Description of NMR BCMB/CHEM 8190

Density Matrix Description of NMR BCMB/CHEM 8190

The Production of Polarization

F-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction

Local Density Functional Theory for Superfluid Fermionic Systems. The Unitary Fermi Gas

Motion of Wavepackets in Non-Hermitian. Quantum Mechanics

Slide 1. Quantum Mechanics: the Practice

Suppose we have observed values t 1, t 2, t n of a random variable T.

DENSITY-FUNCTIONAL THEORY, AND DENSITY-FUNCTIONAL PRACTICE

Ferent equation of the Universe

I-POLYA PROCESS AND APPLICATIONS Leda D. Minkova

Lecture 17: Kinetics of Phase Growth in a Two-component System:

19 The Born-Oppenheimer Approximation

( ) ( ) ( ) ( ) ( ) ( ) j ( ) A. b) Theorem

Scattering at an Interface: Oblique Incidence

3.012 Fund of Mat Sci: Bonding Lecture 1 bis. Photo courtesy of Malene Thyssen,

( ) ( ) Weibull Distribution: k ti. u u. Suppose t 1, t 2, t n are times to failure of a group of n mechanisms. The likelihood function is

Physics 240: Worksheet 16 Name

We are estimating the density of long distant migrant (LDM) birds in wetlands along Lake Michigan.

Solving the Dirac Equation: Using Fourier Transform

POISSON S EQUATION 2 V 0

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER 2

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

Chapter Lagrangian Interpolation

Set of square-integrable function 2 L : function space F

6.8 Laplace Transform: General Formulas

X-Ray Notes, Part III

Notes on the stability of dynamic systems and the use of Eigen Values.

Topic 4a Introduction to Root Finding & Bracketing Methods

NATIONAL UNIVERSITY OF SINGAPORE PC5202 ADVANCED STATISTICAL MECHANICS. (Semester II: AY ) Time Allowed: 2 Hours

Chapter Finite Difference Method for Ordinary Differential Equations

Addition & Subtraction of Polynomials

Goodness-of-fit for composite hypotheses.

Basic molecular dynamics

Lesson 5. Chapter 7. Wiener Filters. Bengt Mandersson. r x k We assume uncorrelated noise v(n). LTH. September 2010

Algebra-based Physics II

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

The Unitary Fermi Gas

DFT functionals in WIEN2k

Flux. Area Vector. Flux of Electric Field. Gauss s Law

Transcription:

Peenaon fo Theoecal Condened Mae Phyc n TU Beln Geen-Funcon and GW appoxmaon Xnzheng L Theoy Depamen FHI May.8h 2005

Elecon n old Oulne Toal enegy---well olved Sngle pacle excaon---unde developng The Geen Funcon GW appoxmaon 2

Pa I: Elecon n old Two majo phyc we wan o know fo elecon n old: Toal Enegy; gound ae----geomey of he cyal he elac popey fo magnec yem: magnec popey In expemen we ue he phooemon and nvee phooemon o pobe he occuped ae and unoccuped ae of old. Excaon Specum. exced ae----sngle pacle excaon n h alk 3

4

Exac oluon fo he wave funcon mpoble. Two way fo olvng h poblem. Appoxmae he wave funcon: Haee and Haee-Fock appoxmaon. Go aound he poblem: ue anohe quany nead of wave funcon o olve he poblem: moe accuae Pa one: Deny Funconal Theoy ---deny Pa wo: Geen Funcon mehod ---popagao Geen Funcon 5

Pa one: 6

DFT whn LDA aleady vey good n he oal enegy calculaon. In pecum calculaon hee ae ome poblem. 7

F-pncple calculaon woe han h. 8

Fom: Mebaek Alouan IPCMS hp://www-pcm.uabg.f/gemm/people/alouan/ The eo n LDA-DFT coeced by GW appoxmaon. Snce all DFT ae LDA-lke baed lack a eaonable decpon of exchange and coelaon. GW appoxmaon moe eaonable on h apec. Pa wo: Geen funcon mehod Pa one: DFT whn LDA In pncple Geen funcon mehod could alo be ued o calculae he oal enegy and he expecaon value of one-pacle opeao. Bu n h alk I wll focu on well behavo on he ngle pacle excaon. 9

Summay fo Pa I Wave funcon of many elecon yem no poble o be ge. DFT and Geen funcon mehod amed a go aound h poblem wh anohe quany. DFT wokng vey well on he oal enegy calculaon. Bu fo he ngle pacle excaon. o ye good enough. Geen funcon mehod wokng bee on ngle pacle excaon now. 0

Pa II: Geen Funcon Repeenaon n Quanum mechanc: Schödnge epeenaon. ˆ H H e 0 Wave funcon---me dependen Opeao lke Hamlona---me ndependen Heenbeg epeenaon. O [ O Hˆ ] O e H ˆ Wave funcon---me ndependen Opeao lke Hamlona---me dependen O 0 e Hˆ

Geneao nead of wave funcon The many-body wave funcon no poble o be ge. : v ˆ ˆ gound ae of a - elecon yem : Annhla on of one elecon a place v : Ceaon of one elecon a place me me v ˆ ˆ : ake one elecon a place me ou fom he gound ae of - v : pu one elecon a place me n o he gound ae of - elecon elecon ˆ [ ˆ Hˆ ] ˆ e H ˆ ˆ 0 e Hˆ 2

G ; ˆ ˆ ˆ ˆ fo > fo elecon < hole Heenbeg pcue: Wave funcon ae ndependen of me and opeao ae me dependen. h Hˆ ˆ d ˆ 2 [ ˆ Hˆ ] h d 0 d ˆ ˆ ˆ v ˆ ˆ 3

[ ] ˆ ˆ ˆ ˆ ; 3 0 T v d G h δ A wo-pacle Geen funcon appea hee o olve he wo-pacle Geen funcon we wll alo ge a hee pacle one n un. To beak h heachy we noduce he ma opeao M Haee poenal elf-enegy: [ ] T v d G M d ˆ ˆ ˆ ˆ 3 Then we ge: ; 3 0 G M d G h δ 4

v M H Σ δ Take: ; 3 0 G d G H Σ δ ; 0 0 G H δ on-neacng yem: The meanng of elf-enegy n analogy o one pacle yem: 0 H H H 0 0 E H Whou peubaon: [ ] 0 E H H Wh peubaon: H Σ 5 Peubaon Many body neacon

So fa we gave he defnon of Geen funcon and equaon fo. Wha he ue fo he Geen funcon? ------ Many many n pncple can be ued o ge: The expecaon value of any ngle-pacle opeao n he gound ae The gound ae enegy The excaon pecum------in h alk only h. Refe o: La Hedn Phycal Revew vol 39 A796965 6

Infomaon fom he Geen funcon: G ; ˆ ˆ ˆ ˆ fo > fo elecon < hole 7

ˆ ˆ Relax ew concep: Quapacle ε Phyc meanng of quapacle nk Fo non-neacng yem: o E E E E ε Fo neacng yem: ˆ o ˆ o E E E E ε 8

9

20

2

So fa we ge he meanng of Popee of ε ε Fo non-neacng yem: nk ε E E o E E Fo neacng yem: ε ˆ o ˆ E E o E E Dffeence: Σ......no Heman ew concep: Quapacle And we need o know how o ge hem 22

> <...... ; * * unocc occ e e G ε ε ; 3 0 G M d G h δ Then we ge: 3 0 M d h ε v M H Σ δ Take: 23 3 0 d H ε Σ

So fa: 3 0 d H ε Σ All he complcy of many body yem ae conaned hee. Poblem: Unknown Can be exended n em of V bu he convegence exemely bad. ; Σ GW appoxmaon 24

Pa III: GW appoxmaon o Quapacle neacon Elecon n old Beak when collon happen and fom agan afe ha. When he nex collon happen follow he ame poce agan and agan Elecon polazaon cloud 25

o o When he collon neval longe han he me needed o fom he cloud we can aume he elecon ae quapacle neacng wh each ohe hough a ceened Coulomb neacon. Th mean he elecon ae can no be oo localzed. F ode exenon of elf-enegy accodng o he ceened Coulomb neacon GW appoxmaon Refe o: La Hedn Phycal Revew vol 39 A796965 26

o GW appoxmaon Hedn GW elf - enegy: Σ ω 2π dω e ω δ G ; ω ω W ; ω Sceenedneacon : W ; ω d" ε "; ω v " Deleccfunon: ε ; ω δ d" v " P0 " ; ω Polazably : P0 " ; ω dω G ; ω ω G ; ω 2π Quapacle enegy: ε qp ε H qp Σ ; ε 27

Haee appoxmaon: Σ 0 Haee-Fock appoxmaon: Σ G ; v " δ " ac only exchange Kohn-Sham whn LDAGGA: XC Σ ; " " V δ " δ " coelaon paly decbed bu local and ac We need h whch phycally meanngful dynamc and non-local o expe he elf-enegy GW appoxmaon 28

Elecon n old V XC Kohn-Sham whn LDAGGA: Haee-Fock GW appoxmaon 29

Fom: Mebaek Alouan IPCMS hp://www-pcm.uabg.f/gemm/people/alouan/ So fa GW appoxmaon manly ued o coec he band ucue of old. Some ae gong beyond o udy he fuhe ode exenon of he elf-enegy. Peonally I hnk vey pey phyc whch decbe he complex many-body naue wh a eaonable effo. 30

Thank fo aenon!!! Thank o Rcado Iaac Gomez Abal Pack Rnke and Choph Feyold fo helpful dcuon! 3