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 |
---|---|
Автор: | |
Формат: | Стаття |
Мова: | 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 Ukraineid |
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
|