Treatments of the exchange energy in density-functional theory

Size: px
Start display at page:

Download "Treatments of the exchange energy in density-functional theory"

Transcription

1 Int. J. Mod. Phys. B (008) 5 arx: Treatments of the echange energy n densty-functonal theory Tamás Gál Secton of Theoretcal Physcs, Insttute of uclear Research of the Hungaran Academy of Scences, 400 ebrecen, Hungary Abstract: Followng a recent work [Gál, Phys. Re. A 64, (00)], a smple deraton of the densty-functonal correcton of the Hartree-Fock equatons, the Hartree- Fock-Kohn-Sham equatons, s presented, completng an ntegrated ew of quantum mechancal theores, n whch the Kohn-Sham equatons, the Hartree-Fock-Kohn-Sham equatons and the ground-state Schrödnger equaton formally stem from a common ground: densty-functonal theory, through ts Euler equaton for the ground-state densty. Along smlar lnes, the Kohn-Sham formulaton of the Hartree-Fock approach s also consdered. Further, t s ponted out that the echange energy of densty-functonal theory bult from the Kohn-Sham orbtals can be gen by degree-two homogeneous -partcle densty functonals (=,,...), formng a sequence of degree-two homogeneous echange-energy densty functonals, the frst element of whch s mnus the classcal Coulomb-repulson energy functonal. Emal address: galt@phys.undeb.hu

2 I. Introducton ensty-functonal theory (FT) [,] forms the eact theoretcal background for partcle-densty based approaches to the quantum mechancal many-body problem, such as theores of the Thomas-Ferm knd [3-7] and of the Slater knd [8-0]. The man dea behnd FT s to replace the quantum mechancal wae-functon wth the partcle densty as basc arable from whch all propertes of a quantum system can be ganed, whch means the use of densty functonals nstead of wae-functon functonals. Ths dea s justfed for ground states formally by the Hohenberg-Kohn theorems [], but n practcal applcatons FT appears as an appromate theory because of the lack of knowledge of an eact epresson for the ground-state energy densty functonal E [ ρ] = F [ ρ] + ρ( r ) ( r ) dr, () the mnmum poston of whch ges the ground-state densty for a gen eternal potental ( r ) and partcle number = ρ dr, () yeldng the Euler equaton or δe[ ρ] = 0, δ ρ( r ) (3a) δf[ ρ] + r ( ) = µ (3b) δρ( r ) [wth µ beng the Lagrange multpler correspondng to the -conserng constrant on ρ( r ), Eq.()], for the determnaton of the ground-state densty. In Eq.(), F[ ρ ], characterstc to the nteracton between the (dentcal) Fermons of the consdered system, s a unersal functonal of the densty ρ( r ) [not dependng on the eternal potental, ( r ) ], whch represents the unknown n E [ ρ ]. To treat a major part of F[ ρ ] eactly, the Kohn-Sham (KS) method [] ntroduces sngle-partcle orbtals [ u ] nto FT, separatng n F the knetc energy

3 Ts = u ( ) u dr (4) = of a ground-state system of nonnteractng fermons wth the densty ρ( r ) of the gen nteractng fermon system, gettng a set of sngle-partcle equatons, KS u + u = ε u =,, (5) wth KS for the determnaton of the (ground-state) densty KS δ ( F[ ρ] Ts [ ρ]) = + r ( ), (6) δρ( r ) ρ = u u, (7) = nstead of usng Eq.(3) drectly. Wth the subtracton of the classcal Coulomb repulson from F[ ρ] [ ρ], T s J[ ρ] = ρ ρ( r ) drdr r r (8) F[ ρ] = T [ ρ] + J[ ρ] + E [ ρ], (9) s then only the relately small echange-correlaton (c) part of the energy remans to be appromated. The echange part of E c s usually defned by u uj uj( r ) u( r ) E = drdr. (0) r r, j= A further reducton of need for appromaton can be acheed by utlzng the defnton Eq.(0) of the echange energy by sngle-partcle orbtals, and settng up the sngle-partcle equatons uj uj( r ) u( r ) u dr + ( r ) u ( r ) = ε u ( r ) =,,, () r r j= wth the local (multplcate) potental c δ( F[ ρ] Ts [ ρ] E [ ρ]) δj[ ρ] δec [ ρ] = + r ( ) = + + r ( ) δρ( r ) δρ( r ) δρ( r ), () whch ges a densty-functonal correcton of the Hartree-Fock (HF) equatons, yeldng the Hartree-Fock-Kohn-Sham method [,3]. The solutons { u } = of the 3

4 dfferent sets of sngle-partcle equatons, Eqs.(5) and (), correspondng to the same ρ( r ), are dfferent. ote that though wth the approach a hgher degree of eactness of practcal calculatons can be reached, t means another step of backng away from a pure densty-functonal theory; n addton to the ntroducton of orbtals, now a nonlocal potental appearng. (The net step would be to ntroduce manypartcle wae-functons to treat correlaton eactly as well, reganng the Schrödnger equaton tself.) Followng the sprt of a recent work by the author [4], n ths artcle a smple deraton of the equatons wll be presented, placng them n a unfed scheme of quantum mechancal theores n whch the Kohn-Sham equatons, the Hartree-Fock- Kohn-Sham equatons and the ground-state Schrödnger equaton orgnate from one root: the Hohenberg-Kohn theorems, through the Euler equaton Eq.(3). Besde the scheme, another mture of FT and the Hartree-Fock theory, the Kohn-Sham formulaton of the Hartree-Fock scheme, wll also be consdered. In the second part of the paper, the concept of -partcle echange-energy densty functonals wll be ntroduced, followng also [4] (see also [5]), and t wll be shown that they emerge naturally as second-degree homogeneous functonals, for eery partcle number. For smplcty, spn s not consdered throughout the paper. II. The Hartree-Fock-Kohn-Sham scheme, and the Kohn-Sham formulaton of the Hartree-Fock theory Buldng on the dea of the constraned search defnton of densty functonals [6-8], the eact ncorporaton of correlaton nto the Hartree-Fock method n the form of a multplcate, densty-dependent potental (whch s a densty-functonal derate), that s, the so-called Hartree-Fock-Kohn-Sham method, can be dered drectly from densty-functonal theory n a smple way. The startng pont s the defnton of the functonal [7] F [ ρ] = mn{ T[ ψ ] + V[ ψ ]}, (3) ψ ρ 4

5 whch, for a ρ( r ), ges the mnmum of the sum T[ ψ ] + V[ ψ ] of wae-functon functonals and = T [ ψ] = ψ ψ (4) V [ ψ] = ψ ψ (5) < j r rj oer the set of normalzed Slater determnants (ψ ) that yeld ρ( r ), that s, ψ $( ρ r) ψ = ρ whle F[ ρ ], defned by F[ ρ] = mn{ T[ ψ ] + V[ ψ ]} ψ ρ, (6) searches oer the space of normalzed (antsymmetrc) wae-functons not restrcted to be determnants. Smlar to the case of the nonnteractng knetc-energy densty functonal Eq.(3), by the mappng T[ ρ] = mn T[ ψ ], (7) s ψ ρ ρ ψ : ρ[ ψ ] = ρ, T[ ψ ] + V[ ψ ] = F[ ρ], (8) assocates eery ρ( r ) wth a fcte nonnteractng system of densty ρ( r ), composed of sngle-partcle orbtals u ( r ) of number = ρ dr. Wrtng Eq.(3) eplctly wth the orbtals, F [ ρ] = mn T [ u, u,..., u, u ] + E [ u, u,..., u, u ] + J[ ρ], (9) { s } { u } ρ wth Ts[ u, u,..., u, u] beng defned by Eq.(4) and E[ u, u,..., u, u] by Eq.(0). The orbtals u ( r ) correspondng to a densty ρ( r ) can be determned through an Euler-Lagrange mnmzaton procedure [6,9,0] (see also p. 5 and p. 87 of []), wth the constrants of fed densty, and orthonormalzaton, u uj dr = δ j, (0) yeldng the Euler-Lagrange equatons u ( r ) j= u ( r ) u ( r j j ) u r r ( r ) dr + λ ρ ( r ) u ( r ) = ε u ( r ) =,, () 5

6 (n canoncal form), wth the Lagrange multplers ε securng the normalzaton of u and λ ρ ( r ) securng the fulflment of Eq.(7) (wth ρ( r ) fed). Ths approach of the mnmzaton problem posed by Eq.(9), howeer, tells nothng more about λ ρ ( r ), about ts connecton to ρ( r ), and hence to ( r ) ; ths s because, as has been ponted out n [4], t does not utlze the defnton Eq.(9), the mplct defnton of the u s correspondng to a gen ρ( r ), fully. Wth a dfferent approach, namely, mnmzng the dfference F = s + [ u, u,..., u, u ]: T [ u, u,..., u, u ] E [ u, u,..., u, u ] ( F J)[ ρ[ u, u,..., u, u ]] () under the orthonormalzaton constrant only (nstead of the aboe mnmzaton of Ts[ u, u,..., u, u] + E[ u, u,..., u, u] ), the ρ-dependence of the multpler that forces the orbtals u to yeld a gen ρ( r ) can be dentfed: the Euler-Lagrange equatons resultng from the mnmzaton of [ u, u,..., u, u ] F, under normalzaton constrant on u ( r ) s, are u j= uj uj( r ) δ( F[ ρ] J[ ρ]) u ( r ) dr u = ε u =,,. (3) r r δρ The bass for ths mnmzaton s that, followng from the defnton of F [ ρ ], [,,...,, ] F 0 u u u u (4) s, the equalty holdng for the orbtals { } for any normalzed u u ( r ) assocated = to ρ( r ) by Eq.(8); thus only the normalzaton of u s has to be reckoned wth as eternal constrant. Wth Eq.(3) then the equatons Eq.() of the method arse mmedatelly through the Hohenberg-Kohn Euler equaton, Eq.(3), that s, δ ( F[ ρ] J[ ρ]) + = µ δρ( r ) wth δ{ F[ ρ] ( F[ ρ] J[ ρ])} = + r ( ) δρ( r ) [Eq.()], whch brngs the eternal potental ( r ) that correspondng to ( r ) ). (5) (6) nto the equatons, fng ρ( r ) (as 6

7 It s mportant to underlne that the correlaton energy functonal of the scheme s not the same as that of the Kohn-Sham theory, snce t s defned by whle n the Kohn-Sham scheme E [ ρ] = F[ ρ] F [ ρ], (7) c E [ ρ] = F[ ρ] T [ ρ] J[ ρ] E [ ρ], (8) c s but F [ ρ] T [ ρ] + J[ ρ] + E [ ρ] (9) s [ E [ ρ ] s defned through Eq.(0) wth u ( r ) s beng the (ρ-dependent) orbtals of the Kohn-Sham scheme]. In partcular, snce F [ ρ ] mnmzes T [ u, u,..., u, u ] + E [ u, u,..., u, u ] n total [under Eq.(7)], s from whch F [ ρ] T [ ρ] + J[ ρ] + E [ ρ], (30) s E c [ ρ] E [ ρ], (3) whch means that, besdes the need for appromaton of E [ ρ ], n the Kohn-Sham method, there s also a greater correlaton part n the energy densty functonal to be appromated. (ote that both knds of correlaton energes are always negate or zero [3].) Besde the densty-functonal amplfcaton of Hartree-Fock theory, dscussed so far, there s another mture of FT and Hartree-Fock theory: the Kohn-Sham formulaton of the Hartree-Fock scheme [9,7,,3] (the Kohn-Sham-Hartree-Fock scheme), where the nonlocal echange operator of the Hartree-Fock equatons s replaced by a local potental whch s a densty-functonal derate. Ths theory ges the sound justfcaton of early attempts n that drecton. As the Hartree-Fock method s an appromate theory by defnton, neglectng correlaton completely by narrowng the space of wae-functons to Slater determnants, ts Kohn-Sham formulaton s not eact ether. Here, the densty functonal F[ ρ ] of eact FT s replaced by F [ ρ ], and the energy functonal HF E [ ρ ] = F[ ρ] + ρ( r ) ( r ) dr (3) c 7

8 s mnmzed, under the -conseraton constrant, gng the Hartree-Fock groundstate energy, and densty. Ths yelds the Euler equaton δf [ ρ] δρ( r ) HF + ( r ) = µ (33) for the determnaton of the Hartree-Fock ground-state densty. ote that, wth E HF [ρ], E [ρ ] arses as HF E [ ρ] = E [ ρ] + E [ ρ]. (34) Separatng T s [ ρ ] n E HF [ ρ ], analogously to eact FT, a major part of the unknown F [ ρ ] can be treated eactly, obtanng the sngle-partcle Schrödnger equatons HF KSHF u + u = ε u =,,, (35) wth Eq.(7), for the determnaton of the densty, where HF KS KS c δ ( F[ ρ] Ts[ ρ]) ( r ) = + ( r ). (36) δρ( r ) These equatons can be obtaned n a smple manner, utlzng the dea descrbed aboe n connecton wth the method; namely, mnmzng the dfference T [ u s, u,..., u, u ]: = T s [ u, u,..., u, u ] T s [ ρ [ u, u,..., u, u ]] (37) under orthonormalzaton constrant on u s (just as n the case of the deraton of the Kohn-Sham equatons gen n [4]), then usng the pure densty-functonal Euler δt [ ] equaton Eq.(33) to substtute for s ρ n the obtaned Euler-Lagrange equatons, δρ( r ) hereby fng ρ( r ) by the eternal potental. The tradtonal quantum mechancal equatons of the Hartree-Fock approach, the Hartree-Fock equatons themseles, can be recoered as well from Eq.(33), followng the aboe procedure, but now F [ u, u,..., u, u] = Ts [ u, u,..., u, u ] + E[ u, u,..., u, u ] ( F J)[ ρ [ u, u,..., u, u]] (38) has to be mnmzed (followed by the nserton of Eq.(33)); just as can the (groundstate) Schrödnger equaton of eact quantum mechancs be dered from the Hohenberg-Kohn Euler equaton [Eq.(3)] [4]. ote that the dfference functonals Eq.() and Eq.(38) are the same; the dfference between the sngle-partcle equatons of the and the HF theores arses from the dfferent [eact, or appromate 8

9 (HF)] physcs behnd the equatons, that s, dfferent FT backgrounds, characterzed by the densty-functonal Euler equatons Eq.(3) and Eq.(33), respectely. Regardng the deratons of Schrödnger equatons from FT, some words need to be sad about the densty-functonal derates (local potentals) appearng n these equatons snce the densty functonals defned through wae-functons are defned only for denstes of nteger norm. Ths problem can be elmnated wth the help of the concept of -partcle densty functonals and -conserng functonal dfferentaton [] (for a dscusson of ts mathematcs, see [3]). In [4], t has been ponted out that the Euler-Lagrange mnmzaton procedure proposed there, and used aboe as well, to determne the correspondng orbtals for a densty ρ( r ), can be carred out by usng an -partcle densty functonal A [ ρ ] nstead of the gen densty functonal A[ ρ ], defnng the dfference functonal wth A [ ρ ], A [ u, u,..., u, u ]: = Au [, u,..., u, u ] A [ ρ [ u, u,..., u, u ]], snce the normalzaton of the orbtals and the faton of ther number consere the norm of ρ( r ). An A [ ρ ] defned for eery ρ( r ) can then be chosen; the tral constant shftng [] (or other homogeneous etenson) of A[ ρ ], e.g., s one proper choce, beng (unconstraned) dfferentable f A[ ρ ] s dfferentable oer the space of ρ( r ) s of the gen norm. In the orbtal equatons resultng from the mnmzaton of [ u, u,..., u, u ], addng δa [ ρ] ρ( r ) dr u r ( ) to both sdes of the equatons δρ( r ) A [], the unconstraned derate of A [ ρ ] s replaced by the -conserng functonal derate, δ A ρ [ ], for whch [] δ ρ( r ) δa[ ρ] δa[ ρ ] =. (39) δ ρ( r ) δ ρ( r ) After usng Eq.(39), fnally, the replacement of the unconstraned derate by the correspondng -conserng derate (of A[ ρ ]) s acheed n the orbtal equatons (for ρ (r ) 's), thereby restrctng the dfferentaton to the allowed doman of ρ( r ) s. ote that smlar consderatons hold also f A[ ρ ] comes from a wae-functon 9

10 functonal A[ ψ ] that s not an orbtal functonal. Wth -conserng densty-functonal derates, the deratons of Schrödnger equatons from FT thus look lke: from (e.g.), usng u u j= j= uj uj( r ) δ( F[ ρ] J[ ρ]) u ( r ) dr u = ε u (40) r r δ ρ δ F[ ρ] + r ( ) = µ, (4) δ ρ( r ) u u ( r j j ) δ J[ ρ ] δ Ec [ ρ ] u ( r ) dr ( r) u = ε u. (4) r r δ ρ( r) δ ρ( r) Another way of treatng the problem of A[ ρ ] (= T s [ ρ ], F [ρ ], or F[ ρ ]) beng defned only for nteger 's s to utlze that a chan rule, namely, embracng δa[ ρ[ g]] δa[ ρ] δρ( ) = d, δ g( ) δρ( ) δ g( K ) δa[ ρ[ g]] δa[ ρ] δρ( ) = d, δ g( ) δ ρ( ) δ g( K ) stll holds (as proed n the Append of [3]) for cases A[ ρ [ g]] where the functonal ρ [g] s such that ρ ( )[ g] satsfes the gen K-conseraton constrant here ρ ( r ) dr = for any (allowed) g ( ). In the present cases, g u, u,..., u, u ) [or = ( g = ( ψ, ψ ) ], and the normalzaton of u (r) 's ensures the fulflment of ρ ( r ) dr =. Thus, δ F[ ρ [ u j, u j ]] δ u u ( r ) u ( r ) dr = u ( r ) u ( r ) dr = δ F[ ρ ] δ ρ( r ) = dr, (43a) δ ρ( r ) δ u or δ F[ ρ [ u j, u j ]] δ u δ F[ ρ ] δ ρ( r ) = dr δ ρ( r ) δ u u ( r ) u ( r ) dr = u ( r ) u ( r ) dr =. (43b) ote that ths chan rule s essental also for the aratonal deraton of the Kohn- Sham equatons themseles. 0

11 It s mportant to note here fnally that, hang the -conserng derates (as potentals) n the dered wae-functon equatons, the fractonal-partcle-number generalzaton [4] of densty functonals can be substtuted for the orgnal densty functonals, obtanng the possblty of relang the constrant of -conseraton on the functonal dfferentaton and usng the unconstraned (left or rght) derates of densty functonals. III. Second-degree homogeneous -partcle echange-energy densty functonals To construct appromatons for densty functonals, knowledge about the structure of the functonals s essental. In ths secton, a natural formal constructon of the echange-energy densty functonal E [ ρ ] of Kohn-Sham FT, a second-degree homogeneous -partcle echange-energy densty functonals, wll be proposed. It s known that E [ ρ ], defned through Eq.(0) wth u ( r ) as the (ρdependent) Kohn-Sham orbtals, scales homogeneously of degree one wth coordnate scalng [5], that s, from whch 3 E [ λρλ ] = λe [ ρ], (44) E δe [ ρ] [ ρ] = ρ r dr δρ( r ) emerges [6]. Besde ths property, the queston of homogenety of E [ ρ ] n ρ tself may also arse, consderng that the classcal Coulomb repulson functonal J[ ρ ], the negate of whch s an eact epresson for E for one-partcle denstes, s (45) homogeneous of degree two,.e. J[ λρ] = λ J[ ρ]. (46) The queston of homogenety wth respect to the densty, and related, or affntate, propertes of the arous components of the energy, has been addressed n seeral works recently [7-39], homogenety beng a ery attracte (strong) analytcal property for functonals. For the nonnteractng knetc-energy densty functonal

12 T s [ ρ ], the frst-degree homogenety n ρ [33] has been shown to be an ncorrect proposal [37,38], though the Wezsäcker (one-partcle knetc-energy) functonal, r T [ ] ρ( ) W ρ = dr 8 ρ( r ) (47) (whch s an eact epresson for T s [ ρ ] for one-partcle denstes), s homogeneous of degree one (n ρ). In the case of E [ ρ ], the naldty of homogenety (of second, or any other, degree) n ρ can be proen n a smple way, as for T s [ ρ ] [37]; snce f homogenety of a degree k held then k ρ ρ E[ ρ] = E J = J[ ρ] =, (48) that s, E [ ρ ] would be k J[ ρ]. ote that from the second-degree homogenety of E [ ρ ] (n ρ), Eq.(5) of [39], dered from the result of [3] for the full electronelectron-nteracton energy densty functonal, would also follow. (Chan and Handy [37] hang shown that the homogenety results of [3-34] are ncorrect n the case of the knetc energy, Joubert [39] has also shown that for the correlaton energy functonal after, n [38], t was ponted out that an essental queston was not treated properly n [3-34], therefore ther results cannot be consdered to be ald.) Wth the ntroducton of the concept of -partcle densty functonals, that s, functonals that are ald epressons for a gen densty functonal for -partcle denstes, howeer, t can be shown that the echange energy can be gen by densty functonals homogeneous of degree two, as n the case of T s, whch can be gen by frst-degree homogeneous -partcle densty functonals, for all partcle numbers [4]. By utlzng the result n [4] that the KS orbtals arse most naturally as degree- homogeneous functonals of the densty, that s, ( ),..., ( ) [ ] = ( ),..., ( ) [ ρ], (49) ( u r u r ) kρ k( u r u r ) E arses as a second-degree homogeneous functonal of the densty, through Eq.(0). The resultant densty functonal s an -partcle densty functonal, E [ ρ ], snce ts form depends on the number of orbtals used n ts constructon. Ths E [ ρ ] thus satsfes

13 that s, E [ λρ] = λ E [ ρ], (50) E δe [ ρ] [ ρ] = ρ dr. (5) δρ( r ) ote that an -partcle echange-energy densty functonal (as any other -partcle densty functonal) can be consdered as a two-arable functonal, E[, ρ ], and the general echange-energy densty functonal emerges as E [ ρ] = E [ ρ, ρ], (5) for whch thus the homogenety Eq.(50) yelds the property E[ ρ] = E[ ρ, λρ] λ. (53) Besde Eq.(5), also the analogue of Eq.(45) can be wrtten for E [ ρ ], for - partcle denstes [ ρ ( r ) ], E δe [ ρ ] [ ρ] = ρ r dr, (54) δρ( r ) snce the dfference between the derates δ E ρ [ ] and δ E ρ [ ] s r -ndependent δρ( r ) δρ( r ) [,3]. Eq.(54) s n accordance wth the fact that coordnate scalng conseres the normalzaton of the densty, wth the use of whch E [ λρ 3 ( λr)] E[ λρ 3 = ( λr)] = λe[ ρ] = λe [ ρ ] (55) [gng also Eq.(54)]; see the Append for the proof for a general ρ( r ). It s worth notng that requrng -partcle densty functonals to be homogeneous of some degree elmnates the ambguty present n ther defnton [5], meanng that the only possble -partcle nonnteractng knetc-energy and echangeenergy densty functonals homogeneous of degree one and two, respectely, are those constructed n [4] and here (see also the Append). The smplest echange energy densty functonal arsng from the condtons of homogenety of degree two n densty scalng [Eq.(50)] and homogenety of degree one n coordnate scalng [Eq.(55), see also Eq.(A6)] s mnus the Coulomb functonal J[ ρ ] [whch s just (the degree-two homogeneous) E [ ρ ], [ ρ] = J[ ρ] ]; just as the E 3

14 Wezsäcker functonal s the smplest epresson for a knetc energy densty functonal satsfyng the two homogenety scalng-condtons for T s [ ρ ] [4] (see also [40]). That J[ρ] s the smplest echange energy functonal satsfyng the homogenety condtons s of course stated from a physcal pont of ew: a "real" echange-energy densty functonal, on one hand, has to be a two-partcle functonal (regardng ts structure), and on the other hand, a term, characterzng the r r nteracton, has to appear eplctly besde the densty-functonal part, snce the echange energy s obtaned a Eq.(0), the orbtals u formng the ρ-dependent part n the functonal; thus, t has to hae a form (or a combnaton of forms lke) f [ ρ] f [ ρ( r )] dr dr. (56) r r ote that wth the separaton of a factor ρ( ) ρ( r r ) n the ntegrand, r r ρ ρ( r ) g[ ρ] g[ ρ( r )] dr dr, (57) r r the remanng part g[ ρ( r )] g[ ρ( r )] (or, smply g [ ρ( r )]) has to be homogeneous of degree zero n ρ( r ) for Eq.(57) to satsfy the scalng condtons. Fnally, the queston of separablty [4] for -partcle densty functonals may be worth takng a look at, referrng also to a recent work [4], whch nestgates the separablty problem of (eplctly) -dependent densty functonals. It can be seen easly that -partcle densty functonals A [ ρ ] correspondng to a separable densty functonal A[ ρ ] are separable for -partcle denstes (only for whch they hae physcal releance): A [ ρ + ρ ] = A[ ρ + ρ ] = A[ ρ ] + A[ ρ ] = A [ ρ ] + A [ ρ ], (58) where + =. IV. Summary Varous treatments of echange n the framework of densty-functonal theory hae been nestgated. Frst, the so-called Hartree-Fock-Kohn-Sham method, whch s 4

15 the densty-functonal amplfcaton of the Hartree-Fock approach, has been consdered. In Sec.II, a smple deraton of the Hartree-Fock-Kohn-Sham equatons has been presented, followng the sprt of [4], showng how these equatons, together wth the Kohn-Sham equatons and the ground-state Schrödnger equaton [4], formally orgnate from one common root: densty-functonal theory, through ts Hohenberg-Kohn Euler equaton. Besde the Hartree-Fock-Kohn-Sham scheme, along the same lnes, n Sec.II another mture of densty-functonal theory and Hartree-Fock theory, the Kohn-Sham formulaton of the Hartree-Fock approach (wth the Kohn- Sham-Hartree-Fock equatons), has also been consdered. In Sec.III, t has been ponted out that the echange energy of ground-state systems of dentcal fermons naturally emerges as a second-degree homogeneous -partcle densty functonal E [ ρ ], for any partcle number, yeldng a = ρ constructon E [ ρ] E [ ρ] of the echange-energy densty functonal, for whch Eq.(53) holds. The frst element of the arsng sequence (that s, =,,...) of degreetwo homogeneous -partcle echange-energy densty functonals E [ ρ ] s just mnus the classcal Coulomb-repulson energy functonal. It has further been shown that E [ ρ ] (homogeneous n ρ) retans the coordnate scalng property of the general echange-energy densty functonal, E [ ρ ], and generally, any -partcle densty functonal homogeneous n the densty has the same coordnate scalng behaour as the general densty functonal to whch t corresponds. Thus, the scalng propertes of E [ ρ ] proposed n Sec.III are just the reerse of those of T s [ ρ ] of [4], whch s frst-degree homogeneous n ρ and second-degree homogeneous wth respect to coordnate scalng. Append: Unqueness of homogeneous -partcle densty functonals, and ther coordnate scalng Gen a densty functonal A[ ρ ], -partcle densty functonals A [ ρ ] can be defned, for whch A [ ρ ] = A[ ρ ] (A) 5

16 (wth ρ dr = ), gng the alue of A for -partcle denstes ρ. The defnton Eq.(A) s of course ambguous. E.g., the free choce of generalzaton of normalzaton represents the ambguty n the case of -partcle nonnteractng knetcenergy densty functonals constructed as descrbed n [4]. Requrng an A [ ρ ] to be homogeneous of some degree k, whch means A k [ λρ] = λ A [ ρ], (A) howeer, leads to a unque defnton for ths A [ ρ ]. For, utlzng Eqs.(A) and (A), k ρ ρ ρ A[ ρ] = A A ρ =, (A3) ρ makng use of that the norm of ρ s [5]. From Eq.(A3) then ρ A [ ρ] = ρ k ρ A, (A4) ρ whch shows that the faton of (on the rght) leads to the -partcle nature. Eq.(A4) s the unque epresson for a degree-k homogeneous -partcle densty functonal correspondng to an A[ ρ ]. Ths unqueness, n the case of the densty functonals T s [ ρ ] and E [ ρ ], means that the homogeneous -partcle densty functonals T s [ ρ ] and E [ ρ ] constructed n [4] and Sec.III of the present artcle are the only frst-degree homogeneous -partcle nonnteractng knetc-energy and second-degree homogeneous -partcle echange-energy densty functonals, respectely. Thus, e.g., T W [ ρ ] and J[ρ] are the only possble frst-degree homogeneous and second-degree homogeneous densty functonals that ge the eact nonnteractng knetc energy and the eact echange energy, respectely, for oneelectron systems. Wth Eq.(A4), the aldty of Eq.(54) [and Eq.(55)] for arbtrary ρ( r ) (not just for ρ( r ) s of norm ),.e. E [ ρ] [ ρ] δe = ρ r δρ( r ) dr, (A5) 6

17 3 can be justfed easly. Utlzng λρλ ( r ) dr = ρ( r ) dr and the homogenety of degree one of E [ ρ ] n coordnate scalng, E ρ [ λρλ ] E λ ( r ) E E [ ] ρ ρλ ρ 3 3 λ ρ ρ λ ρ = = =, (A6) whch means that E [ ρ ] tself, too, s homogeneous of degree one n coordnate scalng, yeldng Eq.(A5). A smlar proof apples for T s [ ρ ] of [4] as well, meanng that T s [ ρ ] has the same coordnate scalng property as T s [ ρ ], that s, homogenety of degree two. Acknowledgements: Grants from the Hungaran Hgher Educaton And Research Foundaton and from OTKA (048675) are gratefully acknowledged. References [] Parr, R. G., and Yang, W., ensty Functonal Theory of Atoms and Molecules (Oford Un. Press, ew York, 989). [] agy, Á., Phys. Rep. 98, (998). [3] Thomas, L. H., Proc. Cambrdge Phl. Soc. 3, 54 (97). [4] Ferm, E., Z. Phys. 48, 73 (98). [5] rac, P. A. M., Proc. Cambrdge. Phl. Soc. 6, 376 (930). [6] on Wezsäcker, C. F., Z. Phys. 96, 43 (935). [7] March,. H., Ad. Phys. 6, (957). [8] Slater, J. C., Phys. Re. 8, 385 (95). [9] Gáspár, R., Acta Phys. Hung. 3, 63 (954) [for Englsh translaton, see J. Mol. Struct. Theochem 50, (000)]. [0] Slater, J. C., The Self-Consstent Feld for Molecules and Solds, Quantum Theory of Molecules and Solds, Vol. 4 (McGraw-Hll, ew York, 974). [] Hohenberg, P., and Kohn, W., Phys. Re. 36, B864 (964). [] Kohn, W., and Sham, L. J., Phys. Re. 40, A33 (965). [3] Holas, A., and March,. H., n ensty Functonal Theory, Topcs n Current 7

18 Chemstry, Vol. 80, ed. R. F. alewajsk (Sprnger, Hedelberg, 996), p. 57. [4] Gál, T., Phys. Re. A 64, (00). [5] Gál, T., Int. J. Quantum Chem. 07, 586 (007). [6] Percus, J. K., Int. J. Quantum Chem. 3, 89 (978). [7] Ley, M., Proc. atl. Acad. Sc. USA 76, 606 (979). [8] Leb, E. H., Int. J. Quantum Chem. 4, 43 (983). [9] Payne, P. W., J. Chem. Phys. 7, 490 (979). [0] Freed, K. F., and Ley, M., J. Chem. Phys. 77, 396 (98). [] Holas, A., March,. H., Takahash, Y., Zhang, C., Phys. Re. A 48, 708 (993). [] Gál, T., Phys. Re. A 63, 0506 (00). [3] Gál, T., J. Math. Chem. 4, 66 (007) [eprnt arx:math-ph/060307]. [4] Perdew, J.P., Parr, R.G., Ley, M., Balduz, J.L., Phys. Re. Lett. 49, 69 (98). [5] Ley, M., and Perdew, J. P., Phys. Re. A 3, 00 (985). [6] Append n Ghosh, S. K., and Parr, R. G., J. Chem. Phys. 8, 3307 (985). [7] Zhao, Q., Morrson, R. C., and Parr, R. G., Phys. Re. A 50, 38 (994). [8] Parr, R. G., and Ghosh, S. K., Phys. Re. A 5, 3564 (995). [9] Lu, S., and Parr, R. G., Phys. Re. A 53, (996). [30] Morrson, R. C., and Parr, R. G., Phys. Re. A 53, R98 (996). [3] Lu, S., and Parr, R. G., Phys. Re. A 55, 79 (997). [3] Parr, R. G., and Lu, S., Chem. Phys. Lett. 76, 64 (997). [33] Lu, S., and Parr, R. G., Chem. Phys. Lett. 78, 34 (997). [34] Parr, R. G., and Lu, S., Chem. Phys. Lett. 80, 59 (997). [35] Joubert,. P., Chem. Phys. Lett. 88, 338 (998). [36] Tozer,. J., Phys. Re. A 58, 354 (998). [37] Chan, G. K.-L., and Handy,. C., Phys. Re. A 59, 670 (999). [38] Gál, T., Phys. Re. A 6, (000). [39] Joubert,., Phys. Re. A 64, (00). [40] Gál, T., and agy, Á., J. Mol. Struct. Theochem 50, 67 (000). [4] Perdew, J. P., Ad. Quantum Chem., 3 (990). [4] Pérez-Jménez, A. J., Moscardó, F., Sancho-García, J. C., Aba, L. P., San-Fabán, E., and Pérez-Jordá, J. M., J. Chem. Phys. 4, 0 (00). 8

19 Append B. Summary of the connectons between the densty-functonal Euler equaton and the arous wae-functon equatons, ndcatng the functonal the constranedsearch defnton of whch s utlzed to get to the gen wae-functon equaton(s) Eact quantum mechancs: Kohn-Sham equatons T s [ ρ ] δf[ ρ] + r ( ) = µ F δρ( r ) [ ρ ] Hartree-Fock-Kohn-Sham eqs. F[ ρ ] Schrödnger equaton Hartree-Fock appromaton: Kohn-Sham-Hartree-Fock eqs. T s [ ρ ] δf[ ρ] + ( r ) = µ δρ( r ) HF F [ ρ ] Hartree-Fock equatons 9

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

Introduction 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 information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 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 information

where 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

where 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 information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

Probabilistic method to determine electron correlation energy

Probabilistic 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 information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Advanced Quantum Mechanics

Advanced 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 information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The 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 information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 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 information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A 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 information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert 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 information

5. Response properties in ab initio schemes

5. 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 information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC 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 information

The Feynman path integral

The 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 information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Supporting information.

Supporting information. Response to Comment on the paper "Restrcted Geometry Optmzaton: A Dfferent Way to Estmate Stablzaton Energes for Aromatc Molecules of Varous Types" Zhong-Heng Yu* and Peng Bao Supportng nformaton. Contents:

More information

10. Canonical Transformations Michael Fowler

10. Canonical Transformations Michael Fowler 10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2 Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on

More information

On the correction of the h-index for career length

On the correction of the h-index for career length 1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat

More information

The non-negativity of probabilities and the collapse of state

The non-negativity of probabilities and the collapse of state The non-negatvty of probabltes and the collapse of state Slobodan Prvanovć Insttute of Physcs, P.O. Box 57, 11080 Belgrade, Serba Abstract The dynamcal equaton, beng the combnaton of Schrödnger and Louvlle

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Errors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation

Errors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation Errors n Nobel Prze for Physcs (7) Improper Schrodnger Equaton and Drac Equaton u Yuhua (CNOOC Research Insttute, E-mal:fuyh945@sna.com) Abstract: One of the reasons for 933 Nobel Prze for physcs s for

More information

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.

This 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 information

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

9 Characteristic classes

9 Characteristic classes THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

A how to guide to second quantization method.

A 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 information

Integrals and Invariants of

Integrals and Invariants of Lecture 16 Integrals and Invarants of Euler Lagrange Equatons NPTEL Course Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng, Indan Insttute of Scence, Banagalore

More information

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach A Bayes Algorthm for the Multtask Pattern Recognton Problem Drect Approach Edward Puchala Wroclaw Unversty of Technology, Char of Systems and Computer etworks, Wybrzeze Wyspanskego 7, 50-370 Wroclaw, Poland

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

Circuit Variables. Unit: volt (V = J/C)

Circuit Variables. Unit: volt (V = J/C) Crcut Varables Scentfc nestgaton of statc electrcty was done n late 700 s and Coulomb s credted wth most of the dscoeres. He found that electrc charges hae two attrbutes: amount and polarty. There are

More information

Poisson brackets and canonical transformations

Poisson brackets and canonical transformations rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner 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 information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)

More information

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 86 Electrcal Engneerng 6 Volodymyr KONOVAL* Roman PRYTULA** PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY Ths paper provdes a

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Thermodynamics and Kinetics of Solids 33. III. Statistical Thermodynamics. Â N i = N (5.3) N i. i =0. Â e i = E (5.4) has a maximum.

Thermodynamics and Kinetics of Solids 33. III. Statistical Thermodynamics. Â N i = N (5.3) N i. i =0. Â e i = E (5.4) has a maximum. hermodynamcs and Knetcs of Solds 33 III. Statstcal hermodynamcs 5. Statstcal reatment of hermodynamcs 5.1. Statstcs and Phenomenologcal hermodynamcs. Calculaton of the energetc state of each atomc or molecular

More information

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

Neumann Asymptotic Eigenvalues of Sturm-liouville Problem with Three Turning Points

Neumann Asymptotic Eigenvalues of Sturm-liouville Problem with Three Turning Points Australan Journal of Basc and Appled Scences 5(5): 89-96 ISSN 99-878 Neumann Asymptotc Egenalues of Sturm-loulle Problem wth Three Turnng Ponts E.A.Sazgar Department of Mathematcs Yerean state Unersty

More information

The internal structure of natural numbers and one method for the definition of large prime numbers

The internal structure of natural numbers and one method for the definition of large prime numbers The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Report on Image warping

Report on Image warping Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Snce 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

Snce 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 information

Maximal Margin Classifier

Maximal Margin Classifier CS81B/Stat41B: Advanced Topcs n Learnng & Decson Makng Mamal Margn Classfer Lecturer: Mchael Jordan Scrbes: Jana van Greunen Corrected verson - /1/004 1 References/Recommended Readng 1.1 Webstes www.kernel-machnes.org

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter J. Basc. Appl. Sc. Res., (3541-546, 01 01, TextRoad Publcaton ISSN 090-4304 Journal of Basc and Appled Scentfc Research www.textroad.com The Quadratc Trgonometrc Bézer Curve wth Sngle Shape Parameter Uzma

More information

Linear Feature Engineering 11

Linear Feature Engineering 11 Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19

More information

Chapter 7 Partially Balanced Incomplete Block Design (PBIBD)

Chapter 7 Partially Balanced Incomplete Block Design (PBIBD) Chapter 7 Partally Balanced Incomplete Block Desgn (PBIBD) The balanced ncomplete block desgns hae seeral adantages. They are connected desgns as well as the block szes are also equal. A restrcton on usng

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

THE SUMMATION NOTATION Ʃ

THE SUMMATION NOTATION Ʃ Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the

More information

SUMMARY OF STOICHIOMETRIC RELATIONS AND MEASURE OF REACTIONS' PROGRESS AND COMPOSITION FOR MULTIPLE REACTIONS

SUMMARY OF STOICHIOMETRIC RELATIONS AND MEASURE OF REACTIONS' PROGRESS AND COMPOSITION FOR MULTIPLE REACTIONS UMMAY OF TOICHIOMETIC ELATION AND MEAUE OF EACTION' POGE AND COMPOITION FO MULTIPLE EACTION UPDATED 0/4/03 - AW APPENDIX A. In case of multple reactons t s mportant to fnd the number of ndependent reactons.

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

More information

Classical Field Theory

Classical Field Theory Classcal Feld Theory Before we embark on quantzng an nteractng theory, we wll take a dverson nto classcal feld theory and classcal perturbaton theory and see how far we can get. The reader s expected to

More information

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0 Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced,

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced, FREQUENCY DISTRIBUTIONS Page 1 of 6 I. Introducton 1. The dea of a frequency dstrbuton for sets of observatons wll be ntroduced, together wth some of the mechancs for constructng dstrbutons of data. Then

More information

Physics 181. Particle Systems

Physics 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 information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

Density matrix. c α (t)φ α (q)

Density matrix. c α (t)φ α (q) Densty matrx Note: ths s supplementary materal. I strongly recommend that you read t for your own nterest. I beleve t wll help wth understandng the quantum ensembles, but t s not necessary to know t n

More information

Numerical Heat and Mass Transfer

Numerical 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 information

Entropy generation in a chemical reaction

Entropy generation in a chemical reaction Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

coordinates. Then, the position vectors are described by

coordinates. Then, the position vectors are described by Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013 Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

The exponential map of GL(N)

The exponential map of GL(N) The exponental map of GLN arxv:hep-th/9604049v 9 Apr 996 Alexander Laufer Department of physcs Unversty of Konstanz P.O. 5560 M 678 78434 KONSTANZ Aprl 9, 996 Abstract A fnte expanson of the exponental

More information