The GW approximation in 90 minutes or so. F. Bruneval Service de Recherches de Métallurgie Physique CEA, DEN
|
|
- Magnus Ward
- 6 years ago
- Views:
Transcription
1 The GW approxmaton n 90 mnutes or so Servce de Recherches de Métallurge Physque CEA, DEN DFT tutoral, Lyon december 2012
2 Outlne I. Standard DFT suffers from the band gap problem II. Introducton of the Green's functon III. The GW approxmaton IV. The GW code n ABINIT and the G0W0 method V. Some applcatons DFT tutoral, Lyon december 2012
3 Standard DFT has unfortunately some shortcomngs band gap Band gap problem! after van Schlfgaarde et al PRL (2008) DFT tutoral, Lyon december 2012
4 A pervasve problem Effectve masses for transport n semconductors Optcal absorpton onset Defect formaton energy, dopant solublty Photoemsson Exp. DFT tutoral, Lyon december 2012
5 How do go beyond wthn the DFT framework? Not easy to fnd mprovement wthn DFT framework There s no such thng as a perturbatve expanson Perdew's Jacob's ladder does not help for the band gap after J. Perdew JCP (2005). Need to change the overall framework! DFT tutoral, Lyon december 2012
6 Outlne I. Standard DFT suffers from the band gap problem II. Introducton of the Green's functon III. The GW approxmaton IV. The GW code n ABINIT and the G0W0 method V. Some applcatons DFT tutoral, Lyon december 2012
7 Many-body perturbaton theory Hstorcally older than the DFT (from the 40-50's)! Bg names: Feynman, Schwnger, Hubbard, Hedn, Lundqvst Green's functons = propagator G r t, r ' t ' = DFT tutoral, Lyon december 2012
8 The Green's functon N,0 Exact ground state wavefuncton: Creaton, annhlaton operator: 1 2 r t N,0 r t, r t s a (N+1) electron wavefuncton not necessarly n the ground state r ' t ' N,0 s another (N+1) electron wavefuncton Let's compare the two of them! DFT tutoral, Lyon december 2012
9 Green's functon defnton N,0 r t r ' t ' N, e = G r t, r ' t ' for t t ' Mesures how an extra electron propagates from (r't') to (rt). DFT tutoral, Lyon december 2012
10 Green's functon defnton N,0 r ' t ' r t N, h for = G r ' t ', r t t ' t Mesures how a mssng electron (= a hole) propagates from (rt) to (r't'). DFT tutoral, Lyon december 2012
11 Fnal expresson for the Green's functon G r t, r ' t ' = N,0 T [ r t r ' t ' ] N, 0 tme-orderng operator e G r t, r ' t ' =G r t, r ' t ' h G r ' t ', r t Compact expresson that descrbes both the propagaton of an extra electron and an extra hole DFT tutoral, Lyon december 2012
12 Lehman representaton G r, r ', t t ' = N,0 T [ r t r ' t ' ] N,0 Closure relaton M, M, M, Lehman representaton: G r, r ', = where { f r f r ' ± E N 1, E N,0 = E N,0 E N 1, Exact exctaton energes! DFT tutoral, Lyon december 2012
13 Related to photoemsson spectroscopy Ekn hν Energy conservaton: before after h E N,0 =E kn E N 1, Quaspartcle energy: =E N,0 E N 1, =E kn h DFT tutoral, Lyon december 2012
14 And nverse photoemsson spectroscopy hν Ekn Energy conservaton: before after E kn E N,0 =h E N 1, Quaspartcle energy: =E N 1, E N,0 =E kn h DFT tutoral, Lyon december 2012
15 Other propertes of the Green's functon Galtsk-Mgdal formula for the total energy: E total = 1 d Tr [ h0 Im G ] Expectaton value of any 1 partcle operator (local or non-local) O =lm Tr [ OG ] t t ' DFT tutoral, Lyon december 2012
16 Outlne I. Standard DFT suffers from the band gap problem II. Introducton of the Green's functon III. The GW approxmaton IV. The GW code n ABINIT and the G0W0 method V. Some applcatons DFT tutoral, Lyon december 2012
17 How to calculate the Green's functon? Feynman dagrams Hedn's functonal approach PRA (1965). DFT tutoral, Lyon december 2012
18 Hedn's coupled equatons 6 coupled equatons: 1= r 1 t 1 1 2= r 2 t 2 2 G 1,2 =G 0 1,2 d34 G0 1,3 3,4 G 4,2 1,2 = d34 G 1,3 W 1,4 4,2,3 1,2,3 = 1,2 1,3 d 4567 Dyson equaton self-energy 1,2 G 4,6 G 5,7 6,7,3 G 4,5 vertex 0 1,2 = d34 G 1,3 G 4,1 3,4,2 polarzablty 1,2 = 1,2 d3 v 1,3 0 3,2 delectrc matrx 1 W 1,2 = d3 1,3 v 3,2 screened Coulomb nteracton DFT tutoral, Lyon december 2012
19 Smplest approxmaton 1,2 = G 1,2 v 1,2 t Fock exchange Dyson equaton: v 1,2 G=G0 G 0 G G=G0 G 0 G0... G 1,2 Not enough: Hartree-Fock s known to perform poorly for solds DFT tutoral, Lyon december 2012
20 Hartree-Fock approxmaton for band gaps DFT tutoral, Lyon december 2012
21 Hedn's coupled equatons 6 coupled equatons: G 1,2 =G 0 1,2 d34 G0 1,3 3,4 G 4,2 1,2 = d34 G 1,3 W 1,4 4,2,3 1,2,3 = 1,2 1,3 d 4567 Dyson equaton self-energy 1,2 G 4,6 G 5,7 6,7,3 G 4,5 0 1,2 = d34 G 1,3 G 4,1 3,4,2 1,2 = 1,2 d3 v 1,3 0 3,2 1 W 1,2 = d3 1,3 v 3,2 screened Coulomb nteracton DFT tutoral, Lyon december 2012
22 Hedn's coupled equatons 6 coupled equatons: G 1,2 =G 0 1,2 d34 G0 1,3 3,4 G 4,2 1,2 = d34 G 1,3 W 1,4 4,2,3 1,2,3 = 1,2 1,3 d 4567 Dyson equaton self-energy 1,2 G 4,6 G 5,7 6,7,3 G 4,5 0 1,2 = d34 G 1,3 G 4,1 3,4,2 1,2 = 1,2 d3 v 1,3 0 3,2 1 W 1,2 = d3 1,3 v 3,2 screened Coulomb nteracton DFT tutoral, Lyon december 2012
23 Hedn's coupled equatons 6 coupled equatons: G 1,2 =G 0 1,2 d34 G0 1,3 3,4 G 4,2 1,2 = d34 G 1,3 W 2 1,4 2 4,2,3 1,2,3 = 1,2 1,3 d 4567 Dyson equaton self-energy 1,2 G 4,6 G 5,7 6,7,3 G 4,5 0 1,2 = d34 G 1,3 G 3,4,2 2 4,1 2 1,2 = 1,2 d3 v 1,3 0 3,2 1 W 1,2 = d3 1,3 v 3,2 screened Coulomb nteracton DFT tutoral, Lyon december 2012
24 Here comes the GW approxmaton 1,2 = G 1,2 W 1,2 GW approxmaton 0 1,2 = G 1,2 G 2,1 RPA approxmaton 1,2 = 1,2 d3 v 1,3 0 3,2 W 1,2 = d3 1 1,3 v 3,2 DFT tutoral, Lyon december 2012
25 What s W? Interacton between electrons n vacuum: 2 1 e v ( r, r ' )= 4 πε0 r r ' Interacton between electrons n a homogeneous polarzable medum: 2 1 e W (r, r ' )= 4 π ε0 ε r r r ' Delectrc constant of the medum Dynamcally screened nteracton between electrons n a general medum: 2 1 e ε ( r, r ' ', ω) W (r, r ', ω)= d r ' ' 4 π ε0 r ' ' r ' DFT tutoral, Lyon december 2012
26 W s frequency dependent W can measured drectly by Inelastc X-ray Scatterng W q=0.80 a.u, Slcon Plasmon frequency ω [ev] Zero below the band gap H-C Wessker et al. PRB (2010) DFT tutoral, Lyon december 2012
27 GW has a super Hartree-Fock GW Approxmaton Hartree-Fock Approxmaton x r 1, r 2 = xc r 1, r 2, = d ' G r 1, r 2, ' v r 1, r 2 2 d ' G r 1, r 2, ' W r 1, r 2, ' 2 = bare exchange x r 1, r 2 Bare exchange c r 1, r 2, + correlaton GW s nothng else but a screened verson of Hartree-Fock. Non Hermtan dynamc DFT tutoral, Lyon december 2012
28 Summary: DFT vs GW Electronc densty r Local and statc exchange-correlaton potental v xc r Approxmatons:, GGA, hybrds Green's functon G r t, r ' t ' Non-local, dynamc Depends onto empty states exchange-correlaton operator = self-energy xc r, r ', GW approxmaton GW r t, r ' t ' =G r t, r ', t ' W r t, r ' t ' DFT tutoral, Lyon december 2012
29 GW approxmaton gets good band gap No more a band gap problem! after van Schlfgaarde et al PRL (2008) DFT tutoral, Lyon december 2012
30 Outlne I. Standard DFT suffers from the band gap problem II. Introducton of the Green's functon III. The GW approxmaton IV. The GW code n ABINIT and the G0W0 method V. Some applcatons DFT tutoral, Lyon december 2012
31 Avalable GW codes DFT tutoral, Lyon december 2012
32 Avalable GW codes has a GW code nsde DFT tutoral, Lyon december 2012
33 How to get G? Remember the Lehman representaton: f r f r ' ± G r, r ', = where the f r and the are complcated quanttes But for ndependent electrons lke Kohn-Sham electrons: G r, r ', = KS KS KS KS r r' ± Ths can be consdered as the best guess for G One can get W and GW DFT tutoral, Lyon december 2012
34 GW as a perturbaton wth respect to GW quaspartcle equaton: [ h0 xc ] = GW GW GW GW KS equaton: [h 0 v Approxmaton : xc ] = GW DFT tutoral, Lyon december 2012
35 GW as a perturbaton wth respect to GW quaspartcle equaton: [ h 0 xc GW ] = GW KS equaton: [ h v ] = GW = xc 0 [ xc GW v xc ] DFT tutoral, Lyon december 2012
36 Lnearzaton of the energy dependance GW = [ xc GW v xc ] Not yet known Taylor expanson: GW xc = xc GW xc... Fnal result: GW = where Z [ xc xc Z =1/ 1 v xc ] DFT tutoral, Lyon december 2012
37 Quaspartcle equaton A typcal ABINIT ouptput for Slcon at Gamma pont k = Band E0 <Vxc> SgX SgC(E0) E^0_gap E^GW_gap GW = Z dsgc/de Sg(E) E-E E Z [ xc v xc ] DFT tutoral, Lyon december 2012
38 Flow chart of a typcal GW calculaton DFT, occuped AND empty states calculate W If self-consstent GW GW, Egenvalues calculate G * W GW k DFT tutoral, Lyon december 2012
39 Outlne I. Standard DFT suffers from the band gap problem II. Introducton of the Green's functon III. The GW approxmaton IV. The GW code n ABINIT and the G0W0 method V. Some applcatons DFT tutoral, Lyon december 2012
40 GW approxmaton gets good band gap No more a band gap problem! after van Schlfgaarde et al PRL (2008) DFT tutoral, Lyon december 2012
41 Clusters de sodum Na 4 e Na 4 4 E 0 Na 4 E 0 Na = { HOMO, Na4 LUMO, Na4 + Na4 /Na4 Bruneval PRL (2009) DFT tutoral, Lyon december 2012
42 Defect calculaton wthn GW approxmaton Up to 215 atoms Cubc slcon carbde DFT tutoral, Lyon december 2012
43 Photolumnescence of VS Bruneval and Roma PRB (2011) DFT tutoral, Lyon december 2012
44 Band Offset at the nterface between two semconductors Very mportant for electroncs! Example: S/SO2 nterface for transstors GW correcton wth respect to R. Shaltaf PRL (2008). DFT tutoral, Lyon december 2012
45 Summary The GW approxmaton solves the band gap problem! The calculatons are extremely heavy, so that we resort to many addtonal techncal approxmatons: method named G0W0 The complexty comes from Dependance upon empty states Non-local operators Dynamc operators that requres freq. convolutons There are stll some other approxmatons lke the Plasmon-Pole model... that I'll dscuss durng the practcal sesson... DFT tutoral, Lyon december 2012
GW approximation in ABINIT
GW approxmaton n ABINIT Servce de Recherches de Métallurge Physque CEA Saclay Outlne I. Introducton: gong beyond DFT II. Introducton of the Green's functon III. Exact Hedn's equatons and the GW approxmaton
More informationThe GW Approximation. Lucia Reining, Fabien Bruneval
, Faben Bruneval Laboratore des Soldes Irradés Ecole Polytechnque, Palaseau - France European Theoretcal Spectroscopy Faclty (ETSF) Belfast, June 2007 Outlne 1 Remnder 2 GW approxmaton 3 GW n practce 4
More informationIntroduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)
Introducton to Densty Functonal Theory Jereme Zaffran nd year-msc. (anochemstry) A- Hartree appromatons Born- Oppenhemer appromaton H H H e The goal of computatonal chemstry H e??? Let s remnd H e T e
More informationProbabilistic method to determine electron correlation energy
Probablstc method to determne electron elaton energy T.R.S. Prasanna Department of Metallurgcal Engneerng and Materals Scence Indan Insttute of Technology, Bombay Mumba 400076 Inda A new method to determne
More information5. Response properties in ab initio schemes
5. Response propertes n ab nto schemes A number of mportant physcal observables s expressed va dervatves of total energy (or free energy) E. Examples are: E R 2 E R a R b forces on the nucle; crtcal ponts
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More information8. Superfluid to Mott-insulator transition
8. Superflud to Mott-nsulator transton Overvew Optcal lattce potentals Soluton of the Schrödnger equaton for perodc potentals Band structure Bloch oscllaton of bosonc and fermonc atoms n optcal lattces
More informationLecture 14: Forces and Stresses
The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal
More informationHomework 4. 1 Electromagnetic surface waves (55 pts.) Nano Optics, Fall Semester 2015 Photonics Laboratory, ETH Zürich
Homework 4 Contact: frmmerm@ethz.ch Due date: December 04, 015 Nano Optcs, Fall Semester 015 Photoncs Laboratory, ETH Zürch www.photoncs.ethz.ch The goal of ths problem set s to understand how surface
More informationNote on the Electron EDM
Note on the Electron EDM W R Johnson October 25, 2002 Abstract Ths s a note on the setup of an electron EDM calculaton and Schff s Theorem 1 Basc Relatons The well-known relatvstc nteracton of the electron
More informationThis chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.
1 Unform Electron Gas Ths chapter llustrates the dea that all propertes of the homogeneous electron gas (HEG) can be calculated from electron densty. Intutve Representaton of Densty Electron densty n s
More informationMolecular Dynamics and Density Functional Theory
Molecular Dynamcs and Densty Functonal Theory What do we need? An account n pemfc cluster: Host name: pemfc.chem.sfu.ca I wll take care of that. Ths can be usually a common account for all of you but please
More informationThree-dimensional Yang-Mills theory as a testbed for truncations of Dyson-Schwinger equations
Introducton Dyson-Schwnger equatons d3 Summary & conclusons Three-dmensonal Yang-Mlls theory as a testbed for truncatons of Dyson-Schwnger equatons Marus Q. Huber Insttute of Physcs, Unversty of Graz ACHT05,
More informationψ ij has the eigenvalue
Moller Plesset Perturbaton Theory In Moller-Plesset (MP) perturbaton theory one taes the unperturbed Hamltonan for an atom or molecule as the sum of the one partcle Foc operators H F() where the egenfunctons
More informationRöntgen s experiment in X-ray Spectroscopy. Röntgen s experiment. Interaction of x-rays x. x-rays. with matter. Wavelength: m
X-ray Spectroscopy Röntgen s experment n 1895 Lecture 1: Introducton & expermental aspects Lecture : Atomc Multplet Theory Crystal Feld Theory CTM4XAS program Lecture 3: Charge Transfer Multplet Theory
More informationLecture 3. Interaction of radiation with surfaces. Upcoming classes
Radaton transfer n envronmental scences Lecture 3. Interacton of radaton wth surfaces Upcomng classes When a ray of lght nteracts wth a surface several nteractons are possble: 1. It s absorbed. 2. It s
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationHW #6, due Oct Toy Dirac Model, Wick s theorem, LSZ reduction formula. Consider the following quantum mechanics Lagrangian,
HW #6, due Oct 5. Toy Drac Model, Wck s theorem, LSZ reducton formula. Consder the followng quantum mechancs Lagrangan, L ψ(σ 3 t m)ψ, () where σ 3 s a Paul matrx, and ψ s defned by ψ ψ σ 3. ψ s a twocomponent
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationarxiv:cond-mat/ v2 [cond-mat.mes-hall] 3 Jan 2006
arxv:cond-mat/0210519v2 [cond-mat.mes-hall] 3 Jan 2006 Non Equlbrum Green s Functons for Dummes: Introducton to the One Partcle NEGF equatons Magnus Paulsson Dept. of mcro- and nano-technology, NanoDTU,
More informationUsing TranSIESTA (II): Integration contour and tbtrans
Usng TranSIESTA (II): Integraton contour and tbtrans Frederco D. Novaes December 15, 2009 Outlne Usng TranSIESTA The ntegraton contour Electron densty from GF Why go complex? The non-equlbrum stuaton Usng
More informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More information1) Silicon oxide has a typical surface potential in an aqueous medium of ϕ,0
1) Slcon oxde has a typcal surface potental n an aqueous medum of ϕ, = 7 mv n 5 mm l at ph 9. Whch concentraton of catons do you roughly expect close to the surface? What s the average dstance between
More informationNeutral-Current Neutrino-Nucleus Inelastic Reactions for Core Collapse Supernovae
Neutral-Current Neutrno-Nucleus Inelastc Reactons for Core Collapse Supernovae A. Juodagalvs Teornės Fzkos r Astronomjos Insttutas, Lthuana E-mal: andrusj@tpa.lt J. M. Sampao Centro de Físca Nuclear da
More informationExactly solvable Richardson-Gaudin models in nuclear structure
Exactly solvable Rchardson-Gaudn models n nuclear structure Jorge Dukelsky In collaboraton wth people n audence: S. Pttel, P. Schuck, P. Van Isacker. And many others Rchardson s Exact Soluton Exact Soluton
More informationRate of Absorption and Stimulated Emission
MIT Department of Chemstry 5.74, Sprng 005: Introductory Quantum Mechancs II Instructor: Professor Andre Tokmakoff p. 81 Rate of Absorpton and Stmulated Emsson The rate of absorpton nduced by the feld
More informationSUPPLEMENTARY INFORMATION
do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of
More informationA how to guide to second quantization method.
Phys. 67 (Graduate Quantum Mechancs Sprng 2009 Prof. Pu K. Lam. Verson 3 (4/3/2009 A how to gude to second quantzaton method. -> Second quantzaton s a mathematcal notaton desgned to handle dentcal partcle
More informationMulti-electron atoms (11) 2010 update Extend the H-atom picture to more than 1 electron: H-atom sol'n use for N-elect., assume product wavefct.
Mult-electron atoms (11) 2010 update Extend the H-atom pcture to more than 1 electron: VII 33 H-atom sol'n use for -elect., assume product wavefct. n ψ = φn l m where: ψ mult electron w/fct φ n l m one
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationMeasurement of Radiation: Exposure. Purpose. Quantitative description of radiation
Measurement of Radaton: Exposure George Starkschall, Ph.D. Department of Radaton Physcs U.T. M.D. Anderson Cancer Center Purpose To ntroduce the concept of radaton exposure and to descrbe and evaluate
More information1. Mean-Field Theory. 2. Bjerrum length
1. Mean-Feld Theory Contnuum models lke the Posson-Nernst-Planck equatons are mean-feld approxmatons whch descrbe how dscrete ons are affected by the mean concentratons c and potental φ. Each on mgrates
More informationA non-perturbative study of the correlation functions of three-dimensional Yang-Mills theory
Introducton Dyson-Schwnger equatons YM n d=3 Summary & conclusons A non-perturbatve study of the correlaton functons of three-dmensonal Yang-Mlls theory arxv:1602.02038 Marus Q. Huber Insttute of Physcs,
More informationTowards three-loop QCD corrections to the time-like splitting functions
Towards three-loop QCD correctons to the tme-lke splttng functons Insttute of Nuclear Physcs (Cracow) Sven-Olaf Moch Unversty of Hamburg Radcor-Loopfest Symposum Los Angeles (USA) 15 Jun 2015 Three-loop
More informationGW Approximation in ωk space
GW Approxmaton n ωk space Valero Olevano March 26, 208 Contents GREEN FUNCTION G 0 2 Green functon G 0 G 0 ζ,ζ 2 ;ω = ζ ζ 2 +ηsgn ǫ 0 G 0 r,r 2 ;ω = r r 2 +ηsgn ǫ 0 G 0 r,r 2 ;ω = ξ ξ 2 G 0 ζ,ζ 2 ;ω =
More information14 The Postulates of Quantum mechanics
14 The Postulates of Quantum mechancs Postulate 1: The state of a system s descrbed completely n terms of a state vector Ψ(r, t), whch s quadratcally ntegrable. Postulate 2: To every physcally observable
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationV.C The Niemeijer van Leeuwen Cumulant Approximation
V.C The Nemejer van Leeuwen Cumulant Approxmaton Unfortunately, the decmaton procedure cannot be performed exactly n hgher dmensons. For example, the square lattce can be dvded nto two sublattces. For
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More information) is the unite step-function, which signifies that the second term of the right-hand side of the
Casmr nteracton of excted meda n electromagnetc felds Yury Sherkunov Introducton The long-range electrc dpole nteracton between an excted atom and a ground-state atom s consdered n ref. [1,] wth the help
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationMultiscale Modeling UCLA Prof. N. Ghoniem
Interdscplnary graduate course Fall 003 Multscale Modelng UCLA Prof. N. Ghonem Introducton n Electronc Structure Calculatons Ncholas Kousss Department of Physcs Calforna State Unversty Northrdge The Fundamental
More informationDepartment of Chemistry Purdue University Garth J. Simpson
Objectves: 1. Develop a smple conceptual 1D model for NLO effects. Extend to 3D and relate to computatonal chemcal calculatons of adabatc NLO polarzabltes. 2. Introduce Sum-Over-States (SOS) approaches
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationProgress & challenges with Luttinger-Ward approaches for going beyond DFT
Progress & challenges with Luttinger-Ward approaches for going beyond DFT Sohrab Ismail-Beigi Yale University Dept. of Applied Physics and Physics & CRISP (NSF MRSEC) Ismail-Beigi, Phys. Rev. B (2010)
More informationSnce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t
8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes
More informationQuantum Field Theory III
Quantum Feld Theory III Prof. Erck Wenberg February 16, 011 1 Lecture 9 Last tme we showed that f we just look at weak nteractons and currents, strong nteracton has very good SU() SU() chral symmetry,
More informationSTATISTICAL MECHANICS
STATISTICAL MECHANICS Thermal Energy Recall that KE can always be separated nto 2 terms: KE system = 1 2 M 2 total v CM KE nternal Rgd-body rotaton and elastc / sound waves Use smplfyng assumptons KE of
More informationBasics of density functional perturbation theory
CECAM Tutoral: Dynamcal, delectrc and magnetc propertes of solds wth ABINIT Lyon, 12-16 may 2014 Bascs of densty functonal perturbaton theory Gan-Marco Rgnanese European Theoretcal Spectroscopy Faclty
More informationRelaxation in water /spin ice models
Relaxaton n water /spn ce models Ivan A. Ryzhkn Insttute of Sold State Physcs of Russan Academy of Scences, Chernogolovka, Moscow Dstrct, 1443 Russa Outlne specfcty of approach quaspartcles Jaccard s theory
More informationThe Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory
The Electronic Structure of Dye- Sensitized TiO 2 Clusters from Many- Body Perturbation Theory Noa Marom Center for Computational Materials Institute for Computational Engineering and Sciences The University
More informationThis model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:
1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors
More informationMany-Body Perturbation Theory. Lucia Reining, Fabien Bruneval
, Fabien Bruneval Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Belfast, 27.6.2007 Outline 1 Reminder 2 Perturbation Theory
More informationIntroduction to circuit analysis. Classification of Materials
Introducton to crcut analyss OUTLINE Electrcal quanttes Charge Current Voltage Power The deal basc crcut element Sgn conventons Current versus voltage (I-V) graph Readng: 1.2, 1.3,1.6 Lecture 2, Slde 1
More informationAnnouncements. Lecture #2
Announcements Lectures wll be n 4 LeConte begnnng Frday 8/29 Addtonal dscusson TA Denns Chang (Sectons 101, 105) Offce hours: Mo 2-3 PM; Th 5-6 PM Lab sectons begn Tuesday 9/2 Read Experment #1 onlne Download
More informationElectronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory
Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.
More informationNon-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields
EJTP 6, No. 0 009) 43 56 Electronc Journal of Theoretcal Physcs Non-nteractng Spn-1/ Partcles n Non-commutng External Magnetc Felds Kunle Adegoke Physcs Department, Obafem Awolowo Unversty, Ile-Ife, Ngera
More information5.76 Lecture #5 2/07/94 Page 1 of 10 pages. Lecture #5: Atoms: 1e and Alkali. centrifugal term ( +1)
5.76 Lecture #5 /07/94 Page 1 of 10 pages 1e Atoms: H, H + e, L +, etc. coupled and uncoupled bass sets Lecture #5: Atoms: 1e and Alkal centrfugal term (+1) r radal Schrödnger Equaton spn-orbt l s r 3
More informationTurbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.
Turbulence Lecture 1 Non-lnear Dynamcs Strong non-lnearty s a key feature of turbulence. 1. Unstable, chaotc behavor.. Strongly vortcal (vortex stretchng) 3 s & 4 s Taylor s work on homogeneous turbulence
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More information5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR
5.0, Prncples of Inorganc Chemstry II MIT Department of Chemstry Lecture 3: Vbratonal Spectroscopy and the IR Vbratonal spectroscopy s confned to the 00-5000 cm - spectral regon. The absorpton of a photon
More informationApplied Nuclear Physics (Fall 2004) Lecture 23 (12/3/04) Nuclear Reactions: Energetics and Compound Nucleus
.101 Appled Nuclear Physcs (Fall 004) Lecture 3 (1/3/04) Nuclear Reactons: Energetcs and Compound Nucleus References: W. E. Meyerhof, Elements of Nuclear Physcs (McGraw-Hll, New York, 1967), Chap 5. Among
More informationEAM series and their history
EM seres and ther hstory Embedded tom Method (EM and Modfed EM (MEM By Kuno TKHSHI Tokyo Insttute of Technology, Tokyo 5-855, JPN Phone/Fax +8-3-5734-395 takahak@mep.ttech.ac.p http://www.mep.ttech.ac.p/ktakahash/
More informationCI/CEPA. Introduction CI Size Consistency Derivation CEPA EPV Results Remarks
CI/CEPA Introducton CI Sze Consstency Dervaton CEPA EPV Results Remarks 1 HY=EY H e Ψ e = E e Ψe Expanson of the Many-Electron Wave Functon Ψ HF Ψ e = Φ o + a a C Φ + ab ab C j Φ j +... Φo = A φ 1 φ 2...
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationMany electrons: Density functional theory Part II. Bedřich Velický VI.
Many electrons: Density functional theory Part II. Bedřich Velický velicky@karlov.mff.cuni.cz VI. NEVF 514 Surface Physics Winter Term 013-014 Troja 1 st November 013 This class is the second devoted to
More information1 Rabi oscillations. Physical Chemistry V Solution II 8 March 2013
Physcal Chemstry V Soluton II 8 March 013 1 Rab oscllatons a The key to ths part of the exercse s correctly substtutng c = b e ωt. You wll need the followng equatons: b = c e ωt 1 db dc = dt dt ωc e ωt.
More informationNeutral Electronic Excitations:
Neutral Electronic Excitations: a Many-body approach to the optical absorption spectra Claudio Attaccalite http://abineel.grenoble.cnrs.fr/ Second Les Houches school in computational physics: ab-initio
More informationŜ z (s 1 ) m (s 1 ), m 1/2, m 1/2. Note that only ion states
JOURAL OF CHEMICAL PHYSICS VOLUME 119, UMBER 4 22 JULY 2003 Physcal nterpretaton and evaluaton of the Kohn Sham and Dyson components of the I relatons between the Kohn Sham orbtal energes and the onzaton
More informationIntermolecular force fields and how they can be determined
Intermolecular force felds and how they can be determned Ad van der Avord Unversty of Njmegen Han-sur-Lesse, December 23 p.1 Equaton of state (Van der Waals) of non-deal gas ( p + a )( ) V 2 V b = kt repulson
More informationIndependent electrons in an effective potential
ABC of DFT Adiabatic approximation Independent electrons in an effective potential Hartree Fock Density Functional Theory MBPT - GW Density Functional Theory in a nutshell Every observable quantity of
More information1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:
68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse
More informationLecture 10. Reading: Notes and Brennan Chapter 5
Lecture tatstcal Mechancs and Densty of tates Concepts Readng: otes and Brennan Chapter 5 Georga Tech C 645 - Dr. Alan Doolttle C 645 - Dr. Alan Doolttle Georga Tech How do electrons and holes populate
More informationOne can coose te bass n te 'bg' space V n te form of symmetrzed products of sngle partcle wavefunctons ' p(x) drawn from an ortonormal complete set of
8.54: Many-body penomena n condensed matter and atomc pyscs Last moded: September 4, 3 Lecture 3. Second Quantzaton, Bosons In ts lecture we dscuss second quantzaton, a formalsm tat s commonly used to
More informationSolutions to Problems Fundamentals of Condensed Matter Physics
Solutons to Problems Fundamentals of Condensed Matter Physcs Marvn L. Cohen Unversty of Calforna, Berkeley Steven G. Loue Unversty of Calforna, Berkeley c Cambrdge Unversty Press 016 1 Acknowledgement
More informationAverage Matrix Relative Sensitivity Factors (AMRSFs) for Auger Electron Spectroscopy (AES)
erage Matrx Relatve Senstvty Factors (AMRSFs) for Auger Electron Spectroscopy (AES) These tables and plots contan AMRSFs for AES calculated for the total peak area for the K, L, M and N shells for both
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationPhotoelectronic properties of chalcopyrites for photovoltaic conversion:
Photoelectronic properties of chalcopyrites for photovoltaic conversion: self-consistent GW calculations Silvana Botti 1 LSI, CNRS-CEA-École Polytechnique, Palaiseau, France 2 LPMCN, CNRS-Université Lyon
More informationOPTICAL POTENTIAL APPROACH TO THE SLOW POSITRON SCATTERING FROM HELIUM ATOM UDC ; A. R. Tančić 1, M. R. Nikolić 2
FACTA UNIVERSITATIS Seres: Physcs, Chemstry and Technology Vol., N o 4, 00, pp. 183-188 OPTICAL POTENTIAL APPROACH TO THE SLOW POSITRON SCATTERING FROM HELIUM ATOM UDC 530.145; 539.1 A. R. Tančć 1, M.
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationAmplification and Relaxation of Electron Spin Polarization in Semiconductor Devices
Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton
More informationLecture 8 Modal Analysis
Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More informationAPPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam
APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development
More informationLecture 4. Macrostates and Microstates (Ch. 2 )
Lecture 4. Macrostates and Mcrostates (Ch. ) The past three lectures: we have learned about thermal energy, how t s stored at the mcroscopc level, and how t can be transferred from one system to another.
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationLinear-response excitations. Silvana Botti
from finite to extended systems 1 LSI, CNRS-CEA-École Polytechnique, Palaiseau, France 2 LPMCN, CNRS-Université Lyon 1, France 3 European Theoretical Spectroscopy Facility September 3, 2008 Benasque, TDDFT
More informationAutomated calculations for MPI
Automated calculatons for MPI Andreas van Hameren Insttute of Nuclear Physcs Polsh Academy of Scences Kraków presented at the 8 th Internatonal Workshop on MPI at the LHC Former Convent of San Agustn,
More information( ) + + REFLECTION FROM A METALLIC SURFACE
REFLECTION FROM A METALLIC SURFACE For a metallc medum the delectrc functon and the ndex of refracton are complex valued functons. Ths s also the case for semconductors and nsulators n certan frequency
More informationElectronic Structure for Excited States (multiconfigurational methods) Spiridoula Matsika
Electronc Structure for Excted States (multconfguratonal methods) Sprdoula Matska Excted Electronc States Theoretcal treatment of excted states s needed for: UV/Vs electronc spectroscopy Photochemstry
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationElectron-Impact Double Ionization of the H 2
I R A P 6(), Dec. 5, pp. 9- Electron-Impact Double Ionzaton of the H olecule Internatonal Scence Press ISSN: 9-59 Electron-Impact Double Ionzaton of the H olecule. S. PINDZOLA AND J. COLGAN Department
More information