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, Techncal Unversty of Denmark February 1, 2008 Abstract Non equlbrum Green s functon methods are regularly used to calculate current and charge denstes n nanoscale (both molecular and semconductor conductors under bas. Ths method s manly used for ballstc conducton but may be extended to nclude nelastc scatterng. In ths tutoral paper the NEGF equatons for the current and charge densty matrx are derved and explaned n a hopefully clear way. 1 Introducton Non equlbrum Green s functon methods are regularly used to calculate current and charge denstes n nanoscale (both molecular and semconductor conductors under bas. An overvew of the theory of molecular electroncs can be found n Ref. [1] and for semconductor nanoscale devces see Ref. [2]. The am of ths text s to provde some ntutve explanatons of one partcle Green s functons n a compact form together wth dervatons of the expressons for the current and the densty matrx. It s not ntended as a complete stand-alone tutoral, but rather as a complement to Ref. [1, 2, 3, 4]. 2 Green s functons Dscrete Schrödnger equaton: H n = E n (1 We dvde the Hamltonan and wavefuncton of the system nto contact (H 1,2, ψ 1,2 and devce (H d, ψ d subspaces: H 1 τ 1 0 τ 1 H d τ ψ 1 2 ψ d = E ψ 1 ψ d (2 0 τ 2 H 2 ψ 2 ψ 2 mpn@mc.dtu.dk 1
where τ 1,2 descrbes the nteracton between devce and contacts. In general we have N contacts (H 1,,N connectng (τ 1,,N the devce H d to the reservors. Here we wll assume that the contacts are ndependent,.e., there are no cross terms (τ between the dfferent contacts. We defne the Green s functon 1 : (E HG(E = I (3 2.1 Why do we want to calculate the Green s functon? The Green s functon gves the response of a system to a constant perturbaton v n the Schrödnger equaton: H ψ = E ψ + v (4 The response to ths perturbaton s: (E H ψ = v (5 ψ = G(E v (6 Why do we need the response to ths type of perturbaton? Well, t turns out that t s usually easer (see next secton to calculate the Green s functon than solve the whole egenvalue problem 2 and most (all for the one-partcle system propertes of the system can be calculated from the Green s functon. E.g., the wavefuncton of the contact ( ψ 2 can be calculated f we know the wavefuncton on the devce ( ψ d. From thrd row of Eq. 2: H 2 ψ 2 + τ 2 ψ d = E ψ 2 (7 (E H 2 ψ 2 = τ 2 ψ d (8 ψ 2 = g 2 (Eτ 2 ψ d (9 where g 2 s the Green s functon of the solated contact 2 ((E H 2 g 2 = I. It s mportant to note that snce we have an nfnte system, we obtan two types of solutons for the Green s functons 3, the retarded and the advanced 4 solutons correspondng to outgong and ncomng waves n the contacts. Notaton: We wll denote the retarded Green s functon wth G and the advanced wth G (and maybe G R and G A occasonally. Here, CAPITAL G denotes the full Green s functon and ts sub-matrces G 1, G d, G 1d etc. Lowercase s used for the Green s functons of the solated subsystems, e.g., (E H 2 g 2 = I. Note that by usng the retarded Green s functon of the solated contact (g 2 n Eq. 9 we obtan the soluton correspondng to a outgong wave n the contact. Usng the advanced Green s functon (g 2 would gve the soluton correspondng to an ncomng wave. 1 Others may (and do use the opposte sgn. 2 Especally for nfnte systems. 3 When the energy concdes wth energy band of the contacts there are two solutons correspondng to outgong or ncomng waves n the contacts. 4 In practce these two solutons are usually obtaned by addng an magnary part to the energy. By takng the lmt to zero of the magnary part one of the two solutons s obtaned. If the lmt 0 + s taken the retarded soluton s found, 0 gves the advanced. Ths can be seen from the Fourer transform of the tme dependent Green s functon. 2
2.2 Self-Energy The reason for calculatng the Green s functon s that t s easer that solvng the Schrödnger equaton. Also, the Green s functon of the devce (G d can be calculated separately wthout calculatng the whole Green s functon (G. From the defnton of the Green s functon we obtan: E H 1 τ 1 0 τ 1 E H d τ 2 0 τ 2 E H 2 Selectng the three equatons n the second column: We can solve Eqs. 11 and 13 for G 1d and G 2d : substtuton nto Eq. 12 gves: from whch G d s smple to fnd: G 1 G 1d G 12 G d1 G d G d2 G 21 G 2d G 2 = I 0 0 0 I 0 0 0 I (10 (E H 1 G 1d τ 1 G d = 0 (11 τ 1 G 1d + (E H d G d τ 2 G 2d = I (12 (E H 2 G 2d τ 2 G d = 0 (13 G 1d = g 1 τ 1 G d (14 G 2d = g 2 τ 2 G d (15 τ 1 g 1τ 1 G d + (E H d G d τ 2 g 2τ 2 G d = I (16 G d = (E H d Σ 1 Σ 2 1 (17 where Σ 1 = τ 1 g 1τ 1 and Σ 2 = τ 2 g 2τ 2 are the so called self-energes. Loosely one can say that the effect of the contacts on the devce s to add the self-energes to the devce Hamltonan snce when we calculate the Green s functon on the devce we just calculate the Green s functon for the effectve Hamltonan H effectve = H d +Σ 1 +Σ 2. However, we should keep n mnd that we can only do ths when we calculate the Green s functon. The egen-values and -vectors of ths effectve Hamltonan are not quanttes we can nterpret easly. For normal contacts, the surface Green s functons g 1 and g 2 used to calculate the self-energes are usually calculated usng the perodcty of the contacts, ths method s descrbed n detal n appendx B of Ref. [3] and n secton 3 of Ref. [2]. 2.3 The spectral functon Another mportant use of the Green s functon s the spectral functon: A = ( G G (18 whch gves the DOS and all solutons to the Schrödnger equaton! To see ths we frst note that for any perturbaton v we get two solutons ( ψ R and ψ A to the perturbed Schrödnger equaton: (E H ψ = v (19 3
from the advanced and retarded Green s functons: ψ R = G v (20 ψ A = G v (21 The dfference of these solutons ( ψ R ψ A = (G G v = A v s a soluton to the Schrödnger equaton: (E H( ψ R ψ A = (E H(G G v = (I I v = 0 (22 whch means that ψ = A v s a soluton to the Schrödnger equaton for any vector v! To show that the spectral functon actually gves all solutons to the Schrödnger equaton s a lttle bt more complcated and we need the expanson of the Green s functon n the egenbass: G = 1 E + δ H = k k k E + δ ǫ k (23 where the δ s the small magnary part (see footnote 3, k s are all egenvectors 5 to H wth the correspondng egenvalues ǫ k. Expandng the spectral functon n the egenbass gves: ( 1 A = E + δ H 1 (24 E δ H = ( 1 1 k k (25 E + δ ǫ k E δ ǫ k k = 2δ k k (E ǫ k k 2 (26 + δ 2 where δ s our nfntesmal magnary part of the energy. Lettng δ go to zero gves: A = 2π k δ (E ǫ k k k (27 2δ (here δ(e ǫ k s the delta functon whch can be seen snce goes to zero everywhere but at (E ǫ k 2 +δ 2 E = ǫ k, ntegratng over E (wth a test functon gves the 2πδ (E ǫ k factor. Eq. 27 shows that the spectral functon gves us all solutons to the Schrödnger equaton. 3 Response to an ncomng wave In the non-equlbrum case, reservors wth dfferent chemcal potentals wll nject electrons and occupy the states correspondng to ncomng waves n the contacts. Therefore, we want to fnd the solutons correspondng to these ncomng waves. Consder contact 1 solated from the other contacts and the devce. At a certan energy we have solutons correspondng to an ncomng wave that s totally reflected at the end of the contact. We wll denote these solutons wth ψ 1,n where 1 s the contact number and n s a quantum number (we may 5 Normalzed! 4
have several modes n the contacts. We can fnd all these solutons from the spectral functon a 1 of the solated contact (as descrbed above. Connectng the contacts to the devce we can calculate the wavefuncton on the whole system caused by the ncomng wave n contact 1. To do ths we note that a wavefuncton should be of the form ψ 1,n + ψ R where ψ 1,n s the totally reflected wave and ψ R s the retarded response of the whole system. Puttng n the anzats ψ 1,n + ψ R nto the Schrödnger equaton gves: E ψ 1,n + τ 1 ψ 1,n + 0 H 1 + τ 1 + H d + τ 1 + τ 2 + H 2 + τ 2 + ( ψ 1,n + ψ R = E ( ψ 1,n + ψ R (28 H 1 + τ 1 + H d + τ 1 + τ 2 + H 2 + τ 2 ψ R = E ( ψ 1,n + ψ R (29 H ψ R = E ψ R τ 1 ψ 1,n (30 (note the slght change n notaton and we see that ψ R s notng else the response of the whole system to a perturbaton of τ 1 ψ 1,n, c.f., Eq. 9: ψ R = Gτ 1 ψ 1,n (31 It s mportant to realze that the scatterng states generated from Eq. 31, usng all possble ncomng waves from each contact, form a complete 6 ON set of solutons to the full Schrödnger equaton[5]. Note that we have chosen the retarded response whch means that the only part of the wave that s travelng towards the devce s the ncomng wave (part of ψ 1,n. We wll make full use of ths fact below. It wll be useful to have the expressons for the devce wavefuncton ψ d and contact wavefuncton ( ψ 1,2. The devce part s straghtforward: and from Eq. 9 or Eq. 15: ψ d = G d τ 1 ψ 1,n (32 ψ 2 = g 2 τ 2 ψ d = g 2 τ 2 G d τ 1 ψ 1,n (33 Note, to calculate the wavefuncton n the contact contanng the ncomng wave (contact 1 we need to add the ncomng wave, gvng a slghtly more complcated expresson: ( ψ 1 = 1 + g 1 τ 1 G d τ 1 ψ 1,n (34 Knowng the wavefunctons correspondng to ncomng waves n dfferent contacts enables us to fll up the dfferent solutons accordng to the electron reservors fllng the contacts. 4 Charge densty matrx In the non equlbrum case we are often nterested n two quanttes: the current and the charge densty matrx. Lets start wth the charge densty (whch allows us to use a self-consstent scheme to descrbe chargng. 6 Except for localzed states n the devce regon. 5
The charge densty matrx s defned as: ρ = k f(k, µ ψ k ψ k (35 where the sum runs over all states wth the occupaton number f(e k, µ (pure densty matrx (note the smlarty wth the spectral functon A, n equlbrum you fnd the densty matrx from A and not as descrbed below. In our case, the occupaton number s determned by the reservors fllng the ncomng waves n the contacts such that: 1 f(e k, µ 1 = (36 1 + e (E k µ 1/k BT s the Ferm-Drac functon wth the chemcal potental (µ 1 and temperature (T of the reservor responsble for njectng the electrons nto the contacts. The wavefuncton on the devce gven by an ncomng wave n contact 1 (see Eq. 32 s: Addng up all states from contact 1 gves: ρ d [contact 1] = = = de k = [Eq. 27] = ntroducng the new quantty Γ 1 = τ 1 a 1τ 1 = ψ d,k = G d τ 1 ψ 1,k (37 f(e, µ 1 δ(e E k ψ d,k ψ d,k (38 de f(e, µ 1 δ(e E k G d τ 1 ψ 1,k ψ 1,k τ 1 G d (39 k ( de f(e, µ 1 G d τ 1 δ(e E k ψ 1,k ψ 1,k τ 1 G d (40 ρ[from contact 1] = 1 2π k de f(e, µ 1 G d τ 1 ( Σ 1 Σ 1 The total charge densty thus becomes a sum over all contacts: a 1 2π τ 1G d (41 we obtan the smple formula: de f(e, µ 1 G d Γ 1 G d (42 ρ = 2 (for spn 2π de f(e, µ G d Γ G d (43 6
5 Probablty Current Havng dfferent chemcal potentals n the reservors fllng the contacts gves rse to a current. In the next secton we wll calculate ths current n a smlar way as the charge densty was calculated. But to do ths we need an expresson for the current from the wavefuncton. In the contnuum case we can calculate the current from the velocty operator. However, for a dscrete Hamltonan t s not so clear what the velocty operator s. Therefore, we derve an expresson for the current from the contnuty equaton (usng two contacts. In steady-state, the probablty to fnd an electron on the devce ( ψ 2 where the sum runs over the devce subspace s conserved: ψ 2 0 = = ψ ψ = ( ψ ψ + ψ ψ (44 t t t t = ī ( ψ H ψ ψ H ψ = ī h h ( ψ H ψ d ψ d H ψ (45 = ī h = ī h ( ψ H d + τ 1 + τ 2 ψ d ψ d H d + τ 1 + τ 2 ψ ([ ] [ ] ψ 1 τ 1 ψ d ψ d τ 1 ψ 1 + ψ 2 τ 2 ψ d ψ d τ 2 ψ 2 We nterpret the term n the frst (square bracket as the ncomng probablty current nto the devce from contact 1 and the second bracket from contact 2. Generalzng to an arbtrary contact j gves us the electrc current (at one energy as the charge (-e tmes the probablty current: j = ē h ( ψ j τ j ψ d ψ d τ j ψ j where j s defned as postve for a current from the contacts nto the devce. We can now put n the expressons for the wavefunctons n the same way as for the densty matrx. 6 Electrcal Current To calculate the total current through the devce we only need to put n the wavefuncton of the devce and the contacts ( ψ d, ψ 1, ψ 2 from Eqs. 32, 34 and 9 and add all the contrbutons together. Thus the current nto the devce from a ncomng wave of one energy (E n contact 1 ( ψ 1,n through the couplng defned by τ 2 s: 2 from 1 = ē h ( ψ 2 τ 2 ψ d ψ d τ 2 ψ 2 = ē h ( ψ 1,n τ 1 G d τ 2 g 2 τ 2G d τ 1 ψ 1,n ψ 1,n τ 1 G d τ 2 g 2τ 2 G d τ 1 ψ 1,n (50 = ē ( h ψ 1,n τ 1 G d τ 2 g 2 g 2 τ 2 G d τ 1 ψ 1,n (51 = ē h ψ 1,n τ 1 G d Γ 2G d τ 1 ψ 1,n (52 (46 (47 (48 (49 7
Addng over the modes n and notng that the levels are flled from the reservor connected to contact 1 gves (2 for spn: I 2 from 1 = 2 ē h = 2e h = 2e h = 2e h = e π h de f(e, µ 1 n δ(e E n ψ 1,n τ 1 G d Γ 2G d τ 1 ψ 1,n (53 de f(e, µ 1 m,n δ(e E n ψ 1,n τ 1 m m G d Γ 2G d τ 1 ψ 1,n (54 f(e, µ 1 m de f(e, µ 1 m ( m G d Γ 2G d τ 1 δ(e E n ψ 1,n ψ 1,n τ 1 m (55 m G d Γ 2G d τ 1 de f(e, µ 1 Tr (G d Γ 2G d Γ 1 n a 1 2π τ 1 m (56 To get the total current through the devce the current from contact two have to be subtracted away: I = e π h whch s exactly the Landauer formula for the current. de (f(e, µ 1 f(e, µ 2 Tr (G d Γ 2G d Γ 1 (57 (58 Acknowledgments The author s grateful for the nspraton from Supryo Datta and Mads Brandbyge. The revsed verson was ntated from dscussons wth Casper Krag and Carsten Rostgård. References [1] M. Paulsson, F. Zahd, and S. Datta. Nanoscence, Engneerng and Technology Handbook, chapter Resstance of a Molecule. CRC Press, 2002. Edtors, D. Brenner, S. Lyshevsk and G. Iafrate, condmat/0208183. [2] P. S. Damle, A. W. Ghosh, and S. Datta. Molecular nanoelectroncs, chapter Theory of nanoscale devce modelng. 2002. Edtor M. Reed. [3] F. Zahd, M. Paulsson, and S. Datta. Advanced Semconductors and Organc Nano-Technques, chapter Electrcal Conducton through Molecules. Academc press, 2002. Edtor H. Markoc. [4] S. Datta. Electronc transport n mesoscopc systems. Cambrdge Unversty Press, Cambrdge, UK, 1997. 8
[5] L. E. Ballentne. Quantum mechancs. World Scentfc, 1998. see chapter 16, Scatterng. 9