Non-perturbative construction of the Luttinger-Ward functional

For a system of correlated electrons, the Luttinger-Ward functional provides a link between static thermodynamic quantities on the one hand and single-particle excitations on the other. The functional is useful in deriving several general properties of the system as well as in formulating the ther...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2006
Автор: Potthoff, M.
Формат: Стаття
Мова:English
Опубліковано: Інститут фізики конденсованих систем НАН України 2006
Назва видання:Condensed Matter Physics
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/121360
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Non-perturbative construction of the Luttinger-Ward functional / M. Potthoff // Condensed Matter Physics. — 2006. — Т. 9, № 3(47). — С. 557–567. — Бібліогр.: 31 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-121360
record_format dspace
spelling irk-123456789-1213602017-06-15T03:02:55Z Non-perturbative construction of the Luttinger-Ward functional Potthoff, M. For a system of correlated electrons, the Luttinger-Ward functional provides a link between static thermodynamic quantities on the one hand and single-particle excitations on the other. The functional is useful in deriving several general properties of the system as well as in formulating the thermodynamically consistent approximations. Its original construction, however, is perturbative as it is based on the weak-coupling skeleton-diagram expansion. Here, it is shown that the Luttinger-Ward functional can be derived within a general functional-integral approach. This alternative and non-perturbative approach stresses the fact that the Luttinger-Ward functional is universal for a large class of models. Для системи скорельованих електронiв функцiонал Люттiнгера-Варда забезпечує зв’язок мiж статичними термодинамiчними величинами, з одного боку, та одночастинковими збудженнями, з iншого боку. Цей функцiонал є придатним для отримання деяких загальних властивостей системи, а також для формулювання термодинамiчно узгоджених наближень. Однак його побудова є пертурбативною, оскiльки вона базується на дiаграмному розкладi при слабкому зв’язку. Показано, що функцiонал Люттiнгера-Варда можна отримати в рамках загального функцiонально-iнтегрального пiдходу. Цей альтернативний непертурбативний пiдхiд пiдкреслює той факт, що функцiонал Лютiнгера-Варда є унiверсальним для широкого класу моделей. 2006 Article Non-perturbative construction of the Luttinger-Ward functional / M. Potthoff // Condensed Matter Physics. — 2006. — Т. 9, № 3(47). — С. 557–567. — Бібліогр.: 31 назв. — англ. 1607-324X PACS: 71.10.-w, 71.15.-m DOI:10.5488/CMP.9.3.557 http://dspace.nbuv.gov.ua/handle/123456789/121360 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description For a system of correlated electrons, the Luttinger-Ward functional provides a link between static thermodynamic quantities on the one hand and single-particle excitations on the other. The functional is useful in deriving several general properties of the system as well as in formulating the thermodynamically consistent approximations. Its original construction, however, is perturbative as it is based on the weak-coupling skeleton-diagram expansion. Here, it is shown that the Luttinger-Ward functional can be derived within a general functional-integral approach. This alternative and non-perturbative approach stresses the fact that the Luttinger-Ward functional is universal for a large class of models.
format Article
author Potthoff, M.
spellingShingle Potthoff, M.
Non-perturbative construction of the Luttinger-Ward functional
Condensed Matter Physics
author_facet Potthoff, M.
author_sort Potthoff, M.
title Non-perturbative construction of the Luttinger-Ward functional
title_short Non-perturbative construction of the Luttinger-Ward functional
title_full Non-perturbative construction of the Luttinger-Ward functional
title_fullStr Non-perturbative construction of the Luttinger-Ward functional
title_full_unstemmed Non-perturbative construction of the Luttinger-Ward functional
title_sort non-perturbative construction of the luttinger-ward functional
publisher Інститут фізики конденсованих систем НАН України
publishDate 2006
url http://dspace.nbuv.gov.ua/handle/123456789/121360
citation_txt Non-perturbative construction of the Luttinger-Ward functional / M. Potthoff // Condensed Matter Physics. — 2006. — Т. 9, № 3(47). — С. 557–567. — Бібліогр.: 31 назв. — англ.
series Condensed Matter Physics
work_keys_str_mv AT potthoffm nonperturbativeconstructionoftheluttingerwardfunctional
first_indexed 2025-07-08T19:42:59Z
last_indexed 2025-07-08T19:42:59Z
_version_ 1837109116157296640
fulltext Condensed Matter Physics 2006, Vol. 9, No 3(47), pp. 557–567 Non-perturbative construction of the Luttinger-Ward functional M.Potthoff Institut für Theoretische Physik und Astrophysik, Universität Würzburg, Am Hubland, D–97074 Würzburg, Germany Received February 10, 2006 For a system of correlated electrons, the Luttinger-Ward functional provides a link between static thermo- dynamic quantities on the one hand and single-particle excitations on the other. The functional is useful in deriving several general properties of the system as well as in formulating the thermodynamically consis- tent approximations. Its original construction, however, is perturbative as it is based on the weak-coupling skeleton-diagram expansion. Here, it is shown that the Luttinger-Ward functional can be derived within a gen- eral functional-integral approach. This alternative and non-perturbative approach stresses the fact that the Luttinger-Ward functional is universal for a large class of models. Key words: Luttinger-Ward functional, conserving approximations, diagrammatic perturbation theory, functional-integral approach, dynamical mean-field theory, self-energy-functional theory PACS: 71.10.-w, 71.15.-m 1. Introduction For a system of correlated electrons in equilibrium, there are several relations [1–3] between static quantities which describe the thermodynamics of the system and dynamic quantities which describe its one-particle excitations. Static quantities are given by the grand potential Ω and its derivatives with respect to temperature T , chemical potential µ etc. The one-electron Green’s function G = G(iωn) or the self-energy Σ = Σ(iωn), on the other hand, are dynamic quantities which yield (equivalent) information on an idealized (photoemission or inverse photoemission) excitation process. The Luttinger-Ward functional Φ̂[G] provides a special relation between static and dynamic quantities with several important properties [4]: Firstly, the grand potential is obtained from the Luttinger-Ward functional evaluated at the exact Green’s function, Φ = Φ̂[G], via Ω = Φ + Tr lnG − TrΣG . (1) Secondly, the functional derivative of the Luttinger-Ward functional, 1 T δΦ̂[G] δG = Σ̂[G] , (2) defines a functional Σ̂[G] which gives the exact self-energy of the system if evaluated at the exact Green’s function. The relation Σ = Σ̂[G] is independent of the Dyson equation G−1 = G−1 0 − Σ. Thirdly, in the non-interacting limit: Φ̂[G] ≡ 0 for U = 0 . (3) Finally, the functional dependence Φ̂[G] is completely determined by the interaction part of the Hamiltonian and independent of the one-particle part: Φ̂[G] universal . (4) c© M.Potthoff 557 M.Potthoff This universality property can also be expressed as follows: Two systems with the same interaction U but different one-particle parameters t (on-site energies and hopping integrals) in the respective Hamiltonian are described by the same Luttinger-Ward functional. Using equation (2), this implies that the functional Σ̂[G] is universal, too. In [4] it is shown by Luttinger and Ward that Φ̂[G] can be constructed order by order in diagrammatic weak-coupling perturbation theory. Φ is obtained as the limit of the infinite series of closed diagrams without any self-energy insertions and with all free propagators in a diagram replaced by fully interacting ones (see figure 1). Generally, this skeleton-diagram expansion cannot be summed up to get a closed form for Φ̂[G]. So, unfortunately, the explicit functional dependence Φ̂[G] is actually unknown – even for the most simple Hamiltonians such as the Hubbard model [5]. The defining properties, equations (1)–(4), however, are easily verified [4]. = + + +Φ Figure 1. Classical definition of the Luttinger-Ward functional � Φ[G]. Double lines: fully inter- acting propagator G. Dashed lines: interaction U. The Luttinger-Ward functional is useful for several general considerations: With the help of Φ̂[G] and the Dyson equation, the grand potential can be considered as a functional of the Green’s function Ω = Ω̂[G] or as functional of the self-energy Ω = Ω̂[Σ], such that Ω is stationary at the physical G or Σ [4,6]. This represents a remarkable variational principle which connects static physical quantities with dynamic ones. The Luttinger-Ward functional is also used in the micro- scopic derivation of some zero- or low-temperature properties of Fermi liquids as discussed in [4,7]. The derivative of the functional, equation (2), shows the self-energy to be gradient field when considered as a functional of the Green’s function, Σ̂[G]. This fact is related to certain symme- try properties of two-particle Green’s functions as originally noted by Baym and Kadanoff [8]. Furthermore, the Luttinger-Ward functional is of great importance in the construction of thermo- dynamically consistent approximations. The so-called conserving (“Φ-derivable”) approximations virtually start from the Luttinger-Ward functional [6,8]. This is essential in order to prove that these approximations obey a number of macroscopic conservation laws. The Hartree-Fock and the random-phase approximations are well-known examples. These “classical” conserving approxima- tions as well as higher-order approaches such as the fluctuation-exchange approximation (FLEX) [9,10] are essentially limited to the weak-coupling regime. In the strong-coupling regime, the non- crossing approximation (NCA) [11,12] and different generalizations (see e. g. [13]) are Φ-derivable. However, the Luttinger-Ward functional can also be used in constructing non-perturbative approx- imations. This was first realized in the context of the dynamical mean-field theory (DMFT) for lattice models of correlated electrons [14–17]. Here, one exploits the universality of the functional, equation (4), to achieve an (approximate) mapping of the original lattice model onto a simpler impurity model with the same interaction part. The fact that Φ̂[G] is the same for a large class of systems, has recently been shown [18,19] to be the key feature that makes it possible to construct several non-perturbative and thermodynamically consistent approximations [20–22]. This idea has been termed “self-energy-functional approach” (SFA). Such general considerations remain valid as long as the Luttinger-Ward functional is well de- fined. This presupposes that the skeleton-diagram expansion is convergent or at least that for- mal manipulations of diagrammatic quantities are consistent in themselves and eventually lead to physically meaningful results. Provided that one can assure that no singular point is passed when starting from the non-interacting Fermi gas and increasing the interaction strength, this seems to be plausible. A strict proof that the skeleton-diagram expansion is well-behaved, however, will hardly be possible in most concrete situations. On the contrary, it is well known that the expansion is questionable in a number of cases, e. g. in case of a symmetry-broken state or a state that is 558 Non-perturbative construction of the LW functional not “adiabatically connected” to the non-interacting limit, such as a Mott insulator. The skeleton- diagram expansion may break down even in the absence of any spontaneous symmetry breaking in a (strongly correlated) state that gradually evolves from a metallic Fermi liquid. This has explicitly been shown by Hofstetter and Kehrein [23] for the narrow-band limit of the single-impurity Ander- son model (see [24,25] for a discussion of possible physical consequences). Generally speaking there is no strict argument available that ensures the convergence of the skeleton-diagram expansion in the strong-coupling regime. The purpose of the present paper is to show that a construction of the Luttinger-Ward functional is possible without making use of the skeleton-diagram expansion. The proposed construction is based on a standard functional-integral approach and avoids the formal complications mentioned above. Thereby, one achieves an alternative and in particular non-perturbative route to the general properties of correlated electron systems derived from the functional, to the dynamical mean-field theory as well as to the self-energy-functional approach. It should be stressed that the intended construction of the Luttinger-Ward functional requires more than a simple definition of the quantity Φ (which could trivially be achieved by using equation (1): Φ ≡ Ω − Tr lnG + TrΣG). The task is rather to provide a functional Φ̂[G] with the properties (1)–(4). Previous approaches are either perturbative or are incapable of proving equations (1)–(4): A construction of the Luttinger-Ward functional different from the original one [4] has been given by Baym [6]: The existence of Φ̂[G] is deduced from a “vanishing curl condition”, δΣ(1, 1′)/δG(2′, 2) = δΣ(2, 2′)/δG(1′, 1), which is derived from an analysis of the functional dependence of G on an arbitrary (time-dependent) external perturbation J. Moreover, an independent functional relation Σ = Σ̂[G] is required. In [6] the latter is assumed to be given by the (full or by a truncated) skeleton-diagram expansion, and consequently this approach is perturbative again. As shown in [6], the Green’s function in the presence of an external field J can be derived from the grand potential Ω̂′[J] as Ĝ(J) = (1/T )δΩ̂′[J]/δJ. Using the inverse functional, Ĵ[G], Legendre transformation yields Ω̂[G] = Ω̂′[Ĵ[G]] − TrGĴ[G]. This (non-perturbative) functional and the Dyson equation can be used to define Φ̂[G] ≡ Ω̂[G] − Tr lnG + Tr (G−1 0 − G−1)G. This idea is in the spirit of the effective action approach [26,27]. Here, however, the problem is that the universality of Φ̂[G], equation (4), cannot be proven. The Luttinger-Ward functional constructed in this way explicitly depends on G0 and hence on the one-particle parameters t. The paper is organized as follows: The next section briefly introduces the notations and the quantities of interest. The construction of the Luttinger-Ward functional is described in section 3. Section 4 gives a brief discussion of the properties of the functional and its use within the dynamical mean-field theory and the self-energy-functional approach. The results are summed up in section 5. 2. Static and dynamic quantities Consider a system of electrons at temperature T and chemical potential µ in thermal equilibrium and let H = H(t,U) = H0(t) + H1(U) be its Hamiltonian where H0(t) = ∑ αβ tαβ c†αcβ , H1(U) = 1 2 ∑ αβγδ Uαβδγ c†αc†βcγcδ . (5) An index α refers to an arbitrary set of quantum numbers characterizing a one-particle basis state. If N is the total particle-number operator, the grand potential of the system is given by Ωt,U = −T lnZt,U where Zt,U = tr exp(−(H(t,U) − µN)/T ) (6) is the partition function. Here and in the following the dependence of all quantities on the one- particle parameters t and the interaction parameters U is made explicit through the subscripts. Using a matrix notation, the free one-particle Green’s function is denoted by Gt,0. Its elements (for fixed µ) are given by: Gt,0,αβ(iωn) = ( 1 iωn + µ − t ) αβ . (7) 559 M.Potthoff Here iωn = i(2n + 1)πT is the n-th Matsubara frequency. The fully interacting Green’s function is denoted by Gt,U. Using Grassmann variables ξα(iωn) = T 1/2 ∫ 1/T 0 dτ eiωnτξα(τ) and ξ∗α(iωn) = T 1/2 ∫ 1/T 0 dτ e−iωnτξ∗α(τ), its elements can be written as [3] Gt,U,αβ(iωn) = −〈ξα(iωn)ξ∗β(iωn)〉t,U = −1 Zt,U ∫ DξDξ∗ξα(iωn)ξ∗β(iωn) exp (At,U,ξξ∗) , (8) where At,U,ξξ∗ = ∑ n,αβ ξ∗α(iωn)((iωn + µ)δαβ − tαβ)ξβ(iωn)− 1 2 ∑ αβγδ Uαβδγ ∫ 1/T 0 dτ ξ∗α(τ)ξ∗β(τ)ξγ(τ)ξδ(τ) (9) is the action. Finally, the self-energy is defined as Σt,U = G−1 t,0 − G−1 t,U . (10) The goal is to construct a functional Φ̂U[G] (where G is considered to be a free variable) which vanishes in the non-interacting case, Φ̂0[G] = 0 [equation (3)], which is universal, i. e. independent of t [equation (4)], which yields Φ̂U[Gt,U] = Ωt,U − Tr lnGt,U + Tr(Σt,UGt,U) if evaluated at the exact Green’s function G = Gt,U [equation (1)], and the derivative of which is a functional Σ̂[G] with Σ̂[Gt,U] = Σt,U [equation (2)]. (Here the notation TrA ≡ T ∑ n ∑ α eiωn0+ Aαα(iωn) is used. 0+ is a positive infinitesimal. Functionals  = Â[ · · · ] are indicated by a hat and should be distinguished clearly from physical quantities A.) For the classical construction of Φ̂U[G] via the skeleton-diagram expansion (figure 1), these properties are easily verified: The universality of the functional [equation (4)] is obvious as any diagram depends on U and on G only; there is no explicit dependence on the free Green’s function Gt,0, i. e. no explicit dependence on t. Since there is no zeroth-order diagram, Φ̂U[G] trivially vanishes for U = 0 [equation (3)]. The functional derivative of Φ̂[G] with respect to G corresponds to the removal of a propagator from each of the Φ diagrams. Taking care of topological factors [1,4], one ends up with the skeleton-diagram expansion for the self-energy, i. e. one gets equation (2). Using equation (2), the Dyson equation (10), and Φt,U ≡ Φ̂U[Gt,U], the µ derivative of the l.h.s. and of the r.h.s of equation (1) are equal for any fixed interaction strength U and temperature T . Namely, (∂/∂µ)(Φt,U +Tr lnGt,U−TrΣt,UGt,U) = TrG−1 t,U(∂Gt,U/∂µ)−TrGt,U(∂Σt,U/∂µ) = −TrGt,U = −〈N〉t,U = ∂Ωt,U/∂µ. Integration over µ yields equation (1). (Note that equation (1) trivially holds for µ → −∞, i. e. for 〈N〉t,U → 0 since Σt,U = 0 and Φt,U = 0 in this limit). An equivalent derivation of equation (1) can be given by a coupling-constant integration [4]. 3. Luttinger-Ward functional The starting point is the standard functional-integral representation of the partition function as given in [3], for example: Define the functional Ω̂U[G−1 0 ] = −T ln ẐU[G−1 0 ] (11) with ẐU[G−1 0 ] = ∫ DξDξ∗ exp ( ÂU,ξξ∗ [G−1 0 ] ) (12) and ÂU,ξξ∗ [G−1 0 ] = ∑ n,αβ ξ∗α(iωn)G−1 0,αβ(iωn)ξβ(iωn) − 1 2 ∑ αβγδ Uαβδγ ∫ 1/T 0 dτ ξ∗α(τ)ξ∗β(τ)ξγ(τ)ξδ(τ) . (13) 560 Non-perturbative construction of the LW functional Ω̂U[G−1 0 ] parametrically depends on U. G−1 0 is considered to be a free variable. At the (matrix inverse of the) exact free Green’s function, G−1 0 = G−1 t,0 , the functional yields the exact grand potential, Ω̂U[G−1 t,0 ] = Ωt,U , (14) of the system with Hamiltonian H = H0(t) + H1(U). Its derivative defines a functional ĝU[G−1 0 ], 1 T δΩ̂U[G−1 0 ] δG−1 0 = − 1 ẐU[G−1 0 ] δẐU[G−1 0 ] δG−1 0 ≡ −ĝU[G−1 0 ] , (15) with the property ĝU[G−1 t,0 ] = Gt,U (16) which is easily verified using equation (8). The strategy to be pursued is the following: ĝU[G−1 0 ] is a universal (t independent) functional and can be used in order to construct a universal relation G = ĜU[Σ] between the one-particle Green’s function and the self-energy independent of the Dyson equation. Using the universal func- tionals Ω̂U[G−1 0 ] and ĜU[Σ], a universal functional F̂U[Σ] is defined the derivative of which essen- tially yields ĜU[Σ]. The Luttinger-Ward functional can be obtained by Legendre transformation and is universal in construction. To start with, consider the equation ĝU [ G−1 + Σ ] = G . (17) This is a relation between the variables G and Σ which, for a given Σ, may be solved for G. This defines a functional ĜU[Σ], i. e. ĝU [ ĜU[Σ]−1 + Σ ] = ĜU[Σ] . (18) For a given self-energy Σ, the Green’s function G = ĜU[Σ] is defined to be the solution of equation (17). From the Dyson equation (10) and equation (16) it is obvious that the relation (17) is satisfied for G and Σ being the exact Green’s function and the exact self-energy, G = Gt,U and Σ = Σt,U, of a system with the interaction U and some set of one-particle parameters t (H = H0(t) + H1(U)). Hence, ĜU[Σt,U] = Gt,U . (19) A brief discussion of the existence and the uniqueness of possible solutions of the relation (17) is given in Appendix A. Using the functionals Ω̂U[G−1 0 ] and ĜU[Σ], a functional F̂U[Σ] can be defined as: F̂U[Σ] = Ω̂U [ ĜU[Σ]−1 + Σ ] − Tr ln ĜU[Σ] . (20) Using equation (15) one finds: 1 T δF̂U[Σ] δΣ = −ĝU [ ĜU[Σ]−1 + Σ ] · ( δĜU[Σ]−1 δΣ + 1 ) − ĜU[Σ]−1 · δĜU[Σ] δΣ , (21) and, using equation (18), 1 T δF̂U[Σ] δΣ = −ĜU[Σ] . (22) So ĜU[Σ] can be considered to be the gradient of the (scalar) self-energy functional F̂U[Σ]. There- with, the Legendre transform of F̂U[Σ] can be constructed: Φ̂U[G] = F̂U [ Σ̂U[G] ] + Tr(Σ̂U[G] G) . (23) Here Σ̂U[G] is the inverse of the functional ĜU[Σ]. The functional can be assumed to be invertible (locally) provided that the system is not at a critical point for a phase transition (see also [18]). Equation (23) defines the Luttinger-Ward functional. 561 M.Potthoff 4. Discussion 4.1. Properties of the Luttinger-Ward functional The properties of the Luttinger-Ward functional, equations (1)–(4), can be easily verified: Equa- tions (10), (14), (19) and (20) imply F̂U[Σt,U] = Ωt,U − Tr lnGt,U , (24) and with Σ̂U[Gt,U] = Σt,U the evaluation of the Luttinger-Ward functional at G = Gt,U yields Φt,U ≡ Φ̂U[Gt,U] = Ωt,U − Tr lnGt,U + Tr(Σt,UGt,U) , (25) i. e. equation (1). From equations (22) and (23), one immediately has: 1 T δΦ̂U[G] δG = Σ̂U[G] , (26) i. e. equation (2). In the limit U = 0, the functionals ĜU=0[Σ] and F̂U=0[Σ] are ill-defined (the domain of the functionals shrinks to a single point, Σ = 0, see Appendix A). However, from equation (25), one directly has ΦU=0[Gt,0] = 0 for any t [see equation (3)] since Σt,0 = 0 and Ωt,0 = Tr lnGt,0 (a proof for the latter can be found in [4]). Finally, the universality of Φ̂U[G], equation (4) is obvious as the definition (23) of the Luttinger-Ward functional involves the universal (t independent) functionals F̂U[Σ] and Σ̂U[G] only. 4.2. Variational principle Using the Legendre transform of the Luttinger-Ward functional F̂U[Σ], one may define Ω̂t,U[Σ] = Tr ln 1 G−1 t,0 − Σ + F̂U[Σ] . (27) The functional derivative is easily calculated: 1 T δΩ̂t,U[Σ] δΣ = 1 G−1 t,0 − Σ − ĜU[Σ] . (28) The equation ĜU[Σ] = 1 G−1 t,0 − Σ (29) is a (highly non-linear) conditional equation for the self-energy of the system H = H0(t)+H1(U): equations (10) and (19) show that it is satisfied by the exact self-energy Σ = Σt,U. Note that the l.h.s. of (29) is independent of t but depends on U (universality of Ĝ[Σ]), while the r.h.s. is independent of U but depends on t via G−1 t,0 . The obvious problem of finding a solution of equation (29) is that there is no closed form for the functional ĜU[Σ]. Solving equation (29) is equivalent to a search for the stationary point of the grand potential as a functional of the self-energy: δΩ̂t,U[Σ] δΣ = 0 . (30) Similarly, one can also construct a variational principle using the Green’s function as the basic variable, δΩ̂t,U[G]/δG = 0. 562 Non-perturbative construction of the LW functional 4.3. Dynamical mean-field theory The dynamical mean-field theory [14–17] basically applies to lattice models of correlated elec- trons with on-site interactions such as the Hubbard model [5], for example. The DMFT aims at an approximate solution of equation (29) and is based on two ingredients: (i) It is important to note that the Luttinger-Ward functional Φ̂U[G] is the same for the lattice (e. g. Hubbard) model and for an impurity model (single-impurity Anderson model). Actually a (decoupled) set of impurity models has to be considered – one impurity model with the respective local interaction at each site of the original lattice. This ensures that the interaction (U) term is the same as in the lattice model. (In case of translational symmetry the a priori different impurity models can be assumed to be equivalent). As U is the same in the lattice and in the impurity model, the Luttinger-Ward functional, as well as ĜU[Σ], is the same. (ii) Let the lattice model be characterized by one-particle parameters t and the impurity model by parameters t′. The fundamental equation (29) for the lattice model would then be solved by the exact self-energy Σt,U. As an ansatz for an approximate solution Σ of equation (29), the self-energy is assumed to be local within the DMFT and to be representable as the exact self-energy of the impurity model for some parameters t′: Σ = Σt′,U . (31) The universality of the Luttinger-Ward functional (i) and the local approximation for the self- energy (ii) are sufficient to derive the DMFT: The insertion of the ansatz (31) into equation (29) yields a conditional equation for the one-particle parameters of the impurity model t′. The l.h.s. becomes ĜU[Σt′,U] = Gt′,U, i. e. the exact Green’s function of the impurity model, while the r.h.s. reads (G−1 t,0 − Σt′,U)−1. The resulting equation for the parameters t′ can be fulfilled only locally, i. e. by equating the local elements of the respective Green’s functions at the impurity and the original site respectively: (Gt′,U)loc = ( 1 G−1 t,0 − Σt′,U ) loc . (32) This is the so-called self-consistency equation of the DMFT [17]. This consideration can be seen as an independent and, in particular, non-perturbative re- derivation of the DMFT which supplements the known approaches such as the cavity method [17]. 4.4. Self-energy-functional approach The universality of the Luttinger-Ward functional or of its Legendre transform F̂U[Σ] is central to the recently developed self-energy-functional approach [18,19]. The SFA is a general variational scheme which includes the DMFT as a special limit. The idea is to take as an ansatz for the self-energy of a model H = H0(t) + H1(U) the exact self-energy Σt′,U of the so-called reference system H ′ = H0(t ′) + H1(U) that shares with the original model the same interaction part. The parameters t′ of the one-particle part are considered to be variational parameters to search for the stationary point of the grand potential as a functional of the self-energy. This means to insert the ansatz Σ = Σt′,U into the general expression (27) and to solve the Euler equation ∂Ω̂t,U[Σt′,U]/∂t′ = 0, i. e.: ∂ ∂t′ ( Tr ln 1 G−1 t,0 − Σt′,U + F̂U[Σt′,U] ) = 0 (33) for t′. If the search for the optimum set of one-particle parameters t′ were unrestricted, the approach would be exact in principle as the Euler equation would then be equivalent to the Euler equation (29) of the general variational principle equation (30). 563 M.Potthoff A restriction of the space of variational parameters becomes necessary in order to evaluate the quantity Ω̂t,U[Σt′,U] which, in general, is impossible as a closed form for the functional F̂U[Σ] is not known. With a proper restriction, however, the reference system H ′ can be made accessible to an exact (numerical) solution which makes it possible to derive the exact grand potential and the exact Green’s function of the system H ′. Therewith, making use of the universality of F̂U[Σ] and using equations (23) and (25) for the reference system, F̂U[Σt′,U] = Ωt′,U − Tr lnGt′,U . (34) Note that this implies that an exact evaluation of F̂U[Σ] is possible for self-energies of an exactly solvable reference system with the same interaction part as the original one. Using this result in equation (33), one obtains: ∂ ∂t′ ( Ωt′,U + Tr ln 1 G−1 t,0 − Σt′,U − Tr lnGt′,U ) = 0 , (35) which can be evaluated to fix t′ and therewith the optimal self-energy and grand potential (see [18–21] for details and concrete examples). 4.5. Luttinger’s theorem Finally, the role of the Luttinger-Ward functional in the derivation of general properties of correlated electron systems will be discussed. As an important example, the Luttinger theorem [4] is considered. For a translationally invariant system, the theorem states that in the limit T → 0 the average particle number is equal to the volume enclosed by the Fermi surface in k space: 〈N〉 = VFS . (36) The Fermi surface is defined by the set of k points in the first Brillouin zone that satisfies µ−ηk = 0 where ηk are the eigenvalues of the matrix t+Σ(ω) at the excitation energy ω = 0 vanishing. Hence, to formulate the Luttinger theorem, one obviously has to presuppose that there is a Fermi surface at all, i. e. that Σ(ω = 0) is Hermitian 1. The original proof of the theorem [4] is perturbative as it makes use of the skeleton-diagram expansion. A non-perturbative proof, based on topological considerations, was proposed recently [29] and is based on the assumption that the system is a Fermi liquid. To discuss the Luttinger theorem in the present context, consider the following shift transfor- mation of the Green’s function S(z)G(iωn) = G(z)(iωn) = G(iωn + iz) (37) with z = 2πkT and k integer (z is a bosonic Matsubara frequency). S(z) is a linear and unitary transformation. The shift transformation leaves the functional integral equation (11) unchanged: Ω̂U[S(z)G−1 0 ] = Ω̂U[G−1 0 ] . (38) To verify this invariance, one has to note that the shift of the Matsubara frequencies in G−1 0 by z can be transformed into a shift ωn → ωn − z in the Grassmann numbers: ξα(iωn) → ξα(iωn − iz) . (39) In imaginary-time representation this shift is equivalent to the multiplication of a phase: ξα(τ) → e−izτξα(τ) , ξ∗α(τ) → eizτξ∗α(τ) . (40) 1For systems without Fermi surface, there is no Luttinger theorem. A nice example is given by the Falicov-Kimball model in infinite dimensions, see [28]. 564 Non-perturbative construction of the LW functional This, however, leaves the functional integral unchanged as the transformation (39) or equation (40) is linear and the Jacobian is unity. Note that antiperiodic boundary conditions ξα(τ = 1/T ) = −ξα(τ = 0) are respected for a bosonic shift frequency z. Denoting Ωt,U(z) ≡ Ω̂U[S(z)G−1 t,0 ], equation (38) states that Ωt,U(z) = Ωt,U(0). Following the steps in the construction of the Luttinger-Ward functional in section 3, one easily verifies that this implies Φt,U(z) = Φt,U(0) where Φt,U(z) ≡ Φ̂U[S(z)Gt,U]. For the Legendre transform, one has Ft,U(z) = Ft,U(0) where Ft,U(z) ≡ F̂U[S(z)Σt,U]. Now, in the limit T → 0, z becomes a continuous variable. Hence, d dz lim T→0 Ft,U(z) = 0 . (41) If the limit and the derivative can be interchanged, d dz lim T→0 Ft,U(z) = lim T→0 d dz Ft,U(z) . (42) Equations (24) and (41) imply − lim T→0 dΩt,U(z) dz = − lim T→0 dTr lnS(z)Gt,U dz . (43) The z dependence of the grand potential is the same as its µ dependence, and thus −(d/dz)Ωt,U(z = 0) = −∂Ωt,U/∂µ = 〈N〉. The evaluation of the r.h.s. in equation (43) is straightforward and can be found in [4], for example. It turns out that at z = 0 the r.h.s. is just the Fermi-surface volume VFS. Consequently, the non-perturbative construction of the Luttinger-Ward functional enables us to reduce the proof of the Luttinger theorem to the proof of equation (42). This, however, requires certain assumptions on the regularity of the T → 0 limit which are generally non-trivial. 5. Summary To summarize, the present paper has shown that the Luttinger-Ward functional can be con- structed within the framework of functional integrals under fairly general assumptions. In partic- ular, there is no need in an adiabatic connection to the non-interacting limit and no expansion is needed in the interaction strength as was required in the original approach by Luttinger and Ward [4]. The construction merely assumes the very existence of the functional integral over Grassmann fields, i. e. the existence of the Trotter limit, for the representation of the partition function. It is well known that the Luttinger-Ward functional can be employed for different purposes, some of which have been discussed here: The functional is used in order to derive some general properties of correlated electron systems, such as the Luttinger theorem. It enables us to formulate a variational principle involving a thermodynamical potential as a functional of the Green’s function or the self-energy and thereby provides a unique and thermodynamically meaningful link between static and dynamic quantities which is helpful for interpretations and for the construction of approximations. An independent derivation of the dynamical mean-field theory is possible using the special properties of the Luttinger-Ward functional and the universality of the functional in particular. The latter is of central importance in the context of the self-energy-functional approach which is a general framework in constructing thermodynamically consistent approximations. Referring to the standard definition of the Luttinger-Ward functional that is based on the weak- coupling skeleton-diagram expansion, the above-mentioned and any further considerations based on the functional and its unique properties meet with criticism when applied to strongly correlated, non-Fermi liquid or symmetry-broken states. This is exactly the point where the presented non- perturbative construction of the Luttinger-Ward functional is beneficial. As the proposed functional-integral construction is rather general, one can think of different future extensions: It suggests generalizing the main idea to many-particle Green’s functions. For example, the grand potential of an interacting electron system can be expressed as a functional of one- and two-particle Green’s functions. The standard construction is diagrammatic again [30] and 565 M.Potthoff could possibly be replaced by a non-perturbative construction along the lines proposed here. Fur- thermore, an extension to the non-equilibrium case might be considered. Finally, disordered (and interacting) electron systems deserve further attention. Here an analogous construction is possible which permits to formulate a new class of non-perturbative and thermodynamically consistent approximations. Work along this direction is in progress [31]. Acknowledgements The author acknowledges helpful discussions with F.F. Assaad, M. Bechmann, W. Hanke and G. Kotliar. The work is supported by the Deutsche Forschungsgemeinschaft (Forschergruppe 538). A. Solutions of relation (17) As the relation (17) is highly non-linear, the existence and the uniqueness of possible solutions have to be discussed: Take U to be fixed and assume that the self-energy given is the exact self-energy of a system H = H0(t) + H1(U) with some hopping parameters t. So the self-energy Σ is assumed to be given from the space SU of t representable self-energies SU ≡ {Σ |Σ = Σt,U , t arbitrary} (U fixed). Using equation (16) it is obvious that the exact Green’s function of this system, G = Gt,U, solves equation (17) as the Dyson equation (10) shows that G−1 t,U + Σt,U is the exact free Green’s function of this system. Concluding, one has ĜU[Σt,U] = Gt,U, and thus the existence of a solution is guaranteed on the space SU. Note that it is very convenient to consider SU as the domain of the functional ĜU[Σ] since this ensures the correct analytical and causal properties of the variable Σ. Under the functional ĜU[Σ] the space SU is mapped onto the space GU of t representable Green’s functions: GU ≡ {G |G = Gt,U , t arbitrary} (U fixed). In general, the map ĜU : SU → GU is not unique globally. An example is given in [18] (see Appendix A therein). However, the uniqueness of the functional ĜU[Σ] can be enforced by a proper restriction of the range GU, i. e. of the solution set of equation (17). This is explained in detail in [18] (Appendix A). The considerations in sections 3 and 4 are unaffected by such a restriction and hold for any choice of the range. The rest of the discussion in Appendix A of [18] can be applied directly to the present situation. In particular, this shows that there is no general reason to question the local uniqueness of ĜU[Σ]. References 1. Abrikosow A.A., Gorkov L.P., Dzyaloshinski I.E. Methods of Quantum Field Theory in Statistical Physics. Prentice-Hall, New Jersey, 1964. 2. Fetter A.L., Walecka J.D. Quantum Theory of Many-Particle Systems. McGraw-Hill, New York, 1971. 3. Negele J.W., Orland H. Quantum Many-Particle Systems. Addison-Wesley, Redwood City, 1988. 4. Luttinger J.M., Ward J.C., Phys. Rev., 1960, 118, 1417. 5. Hubbard J., Proc. R. Soc. London A, 1963, 276, 238. 6. Baym G., Phys. Rev., 1962, 127, 1391. 7. Luttinger J.M., Phys. Rev., 1960, 119, 1153. 8. Baym G., Kadanoff L.P., Phys. Rev., 1961, 124, 287. 9. Bickers N.E., Scalapino D.J., White S.R., Phys. Rev. Lett., 1989, 62, 961. 10. Bickers N.E., Scalapino D.J., Ann. Phys. (N.Y.), 1989, 193, 206. 11. Keiter H., Kimball J.C., J. Appl. Phys., 1971, 42, 1460. 12. Grewe N., Keiter H., Phys. Rev. B, 1981, 24, 4420. 13. Haule K., Kirchner S., Kroha J., Wölfle P., Phys. Rev. B, 2001, 64, 155111. 14. Metzner W., Vollhardt D., Phys. Rev. Lett., 1989, 62, 324. 15. Georges A., Kotliar G., Phys. Rev. B, 1992, 45, 6479. 16. Jarrell M., Phys. Rev. Lett., 1992, 69, 168. 17. Georges A., Kotliar G., Krauth W., Rozenberg M.J., Rev. Mod. Phys., 1996, 68, 13. 18. Potthoff M., Euro. Phys. J. B, 2003, 32, 429. 566 Non-perturbative construction of the LW functional 19. Potthoff M., Euro. Phys. J. B, 2003, 36, 335. 20. Potthoff M., Aichhorn M., Dahnken C., Phys. Rev. Lett., 2003, 91, 206402. 21. Dahnken C., Aichhorn M., Hanke W., Arrigoni E., Potthoff M., Phys. Rev. B, 2004, 70, 245110. 22. Potthoff M., Adv. Solid State Phys., 2005, 45, 135. 23. Hofstetter W., Kehrein S., Phys. Rev. B, 1999, 59, R12732. 24. Kehrein S., Phys. Rev. Lett., 1998, 81, 3912. 25. Kotliar G., Euro. Phys. J. B, 1999, 11, 27. 26. Savrasov S.Y., Kotliar G., Phys. Rev. B, 2004, 69, 245101. 27. Georges A. Strongly correlated electron materials: Dynamical mean-field theory and electronic struc- ture. Lectures on the Physics of Highly Correlated Electron Systems VIII. American Institute of Physics Conference Proceedings 2004, 715, 3. 28. Freericks J.K., Zlatić V., Rev. Mod. Phys., 2003, 75, 1333. 29. Oshikawa M., Phys. Rev. Lett., 2000, 84, 3370. 30. Almbladh C.O., von Barth U., van Leeuwen R., Int. J. Mod. Phys. B, 1999, 13, 535. 31. Potthoff M., Balzer M. (to be published). Непертурбативна побудова функцiоналу Люттiнгера-Варда М.Потхоф Iнститут теоретичної фiзики та астрофiзики, Унiверситет Вюрцбург, Хубланд, D–97074 Вюрцбург, Нiмеччина Отримано 10 лютого 2006 р. Для системи скорельованих електронiв функцiонал Люттiнгера-Варда забезпечує зв’язок мiж стати- чними термодинамiчними величинами, з одного боку, та одночастинковими збудженнями, з iншого боку. Цей функцiонал є придатним для отримання деяких загальних властивостей системи, а також для формулювання термодинамiчно узгоджених наближень. Однак його побудова є пертурбатив- ною, оскiльки вона базується на дiаграмному розкладi при слабкому зв’язку. Показано, що функцiо- нал Люттiнгера-Варда можна отримати в рамках загального функцiонально-iнтегрального пiдходу. Цей альтернативний непертурбативний пiдхiд пiдкреслює той факт, що функцiонал Лютiнгера-Варда є унiверсальним для широкого класу моделей. Ключовi слова: функцiонал Люттiнгера-Варда, наближення збереження, дiаграмна теорiя збурень, функцiонально-iнтегральний пiдхiд, динамiчна теорiя середнього поля, власноенергетична функцiональна теорiя PACS: 71.10.-w, 71.15.-m 567 568