Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0609694 · Physical Review B

Total energies from variational functionals of the Green function and the renormalized four-point vertex

Robert van Leeuwen, Nils Erik Dahlen, Adrian Stan

Published Nov 6, 2006• 23 citations
arXivHTMLPDF

Abstract

We derive variational expressions for the grand potential or action in terms of the many-body Green function $G$ which describes the propagation of particles and the renormalized four-point vertex $\ensuremath{\Gamma}$ which describes the scattering of two particles in many-body systems. The main ingredient of the variational functionals is a term we denote as the $\ensuremath{\Xi}$-functional which plays a role analogously to the usual $\ensuremath{\Phi}$-functional studied by Baym [G. Baym, Phys. Rev. 127, 1391 (1962)] in connection with the conservation laws in many-body systems. We show that any $\ensuremath{\Xi}$-derivable theory is also $\ensuremath{\Phi}$-derivable and therefore respects the conservation laws. We further set up a computational scheme to obtain accurate total energies from our variational functionals without having to solve computationally expensive sets of self-consistent equations. The input of the functional is an approximate Green function $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{G}$ and an approximate four-point vertex $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Gamma}}$ obtained at a relatively low computational cost. The variational property of the functional guarantees that the error in the total energy is only of second order in deviations of the input Green function and vertex from the self-consistent ones that make the functional stationary. The functionals that we will consider for practical applications correspond to infinite order summations of ladder and exchange diagrams and are therefore particularly suited for applications to highly correlated systems. Their practical evaluation is discussed in detail.

DOI
10.1103/physrevb.74.195105

Related papers

Green functions in the renormalized many-body perturbation theory

V. Janiš

I review the way the many-body Green functions are used to renormalize the perturbation theory of correlated fermions. The Green functions are introduced to implement systematically dynamical corrections to the static mean-field theory. The renormalizations enter the perturbation theory via self-consistent evaluations of one-particle and two-particle Green functions. The one-particle self-consistency is discussed within the Baym-Kadanoff construction. We point to an inherent dichotomy of the Baym-Kadanoff approach in the relation between one and two-particle Green functions. They are the Schwinger-Dyson equation, directly derived from the generating Luttinger-Ward functional, and the Ward identity. The latter imposes a severe restriction on the two-particle irreducible vertex and is not guaranteed in the Baym-Kadanoff approach. No approximate theory is capable to obey both relations. I further discuss the two-particle approach with two-particle renormalizations that offers a way to reconcile the Schwinger-Dyson equation and the Ward identity. The two-particle irreducible vertex from the critical Bethe-Salpeter equation becomes a generating functional in the two-particle renormalization of the perturbation theory. The irreducible vertex is used in the Ward identity to determine the order parameter beyond the critical point and in the Bethe-Salpeter equation representing the vertex of the Schwinger-Dyson equation to determine the spectral self-energy in such a way that the thermodynamic critical behavior from both equations is qualitatively the same. Finally, I discuss a static mean-field-like approximation with a renormalized effective interaction determined self-consistently that leads to a qualitatively correct description of strong electron correlations in impurity and lattice models.

arXiv:2109.04231•Sep 9, 2021•0 citations

Generating functional for the full parquet approximation

V. Janiš

Parquet diagrams sum self-consistently Feynman graphs for the vertex function with all two-particle multiple scatterings. We show how the complete parquet equations for the Hubbard-like models can be integrated to a generating functional from which all thermodynamic quantities are derived via (functional) derivatives. An explicit Luttinger-Ward functional $Φ[G;Λ,{\cal K};U]$ is constructed containing only the renormalized one-particle, $G$, irreducible, $Λ$, and reducible, ${\cal K}$, two-particle propagators as independent variational functions. The parquet approximation is proven to be a thermodynamically consistent, $Φ$-derivable theory obeying the necessary conservation laws.

arXiv:cond-mat/9806118•Jun 9, 1998•1 citation

Approximate energy functionals for one-body reduced density matrix functional theory from many-body perturbation theory

Klaas J. H. Giesbertz, Anna-Maija Uimonen, Robert van Leeuwen

We develop a systematic approach to construct energy functionals of the one-particle reduced density matrix (1RDM) for equilibrium systems at finite temperature. The starting point of our formulation is the grand potential Ω [ G ] regarded as variational functional of the Green’s function G based on diagrammatic many-body perturbation theory and for which we consider either the Klein or Luttinger–Ward form. By restricting the input Green’s function to be one-to-one related to a set on one-particle reduced density matrices (1RDM) this functional becomes a functional of the 1RDM. To establish the one-to-one mapping we use that, at any finite temperature and for a given 1RDM γ in a finite basis, there exists a non-interacting system with a spatially non-local potential v [ γ ] which reproduces the given 1RDM. The corresponding set of non-interacting Green’s functions defines the variational domain of the functional Ω . In the zero temperature limit we obtain an energy functional E [ γ ] which by minimisation yields an approximate ground state 1RDM and energy. As an application of the formalism we use the Klein and Luttinger–Ward functionals in the GW-approximation to compute the binding curve of a model hydrogen molecule using an extended Hubbard Hamiltonian. We compare further to the case in which we evaluate the functionals on a Hartree–Fock and a Kohn–Sham Green’s function. We find that the Luttinger–Ward version of the functionals performs the best and is able to reproduce energies close to the GW energy which corresponds to the stationary point.

arXiv:1804.09921•The European Physical Journal B•Nov 1, 2018•15 citations

Non-Equilibrium Green’s Functions: Variational Relations and Approximations for Particle Interactions

Norman J. M. Horing

Abstract Chapter 09 Nonequilibrium Green’s functions (NEGF), including coupled-correlated (C) single- and multi-particle Green’s functions, are defined as averages weighted with the time-development operator U(t0+τ,t0). Linear conductivity is exhibited as a two-particle equilibrium Green’s function (Kubo-type formulation). Admitting particle sources (S:η,η+) and non-conservation of number, the non-equilibrium multi-particle Green’s functions are constructed with numbers of creation and annihilation operators that may differ, and they may be derived as variational derivatives with respect to sources η,η+ of a generating functional eW=TrU(t0+τ,t0)CS/TrU(t0+τ,t0)C. (In the non-interacting case this yields the n-particle Green’s function as a permanent/determinant of single-particle Green’s functions.) These variational relations yield a symmetric set of multi-particle Green’s function equations. Cumulants and the Linked Cluster Theorem are discussed and the Random Phase Approximation (RPA) is derived variationally. Schwinger’s variational differential formulation of perturbation theories for the Green’s function, self-energy, vertex operator, and also shielded potential perturbation theory, are reviewed. The Langreth Algebra arises from analytic continuation of integration of products of Green’s functions in imaginary time to the real-time axis with time-ordering along the integration contour in the complex time plane. An account of the Generalized Kadanoff-Baym Ansatz is presented.

•Jul 27, 2017•0 citations

Functionals in many-particle physics

Richard M. Martin, Lucia Reining, David M. Ceperley

To avoid confusion, it is necessary to be quite explicit about what is assumed and what is to be proved. A. Klein, Phys. Rev. 121, 950 (1961) Summary The topic of this chapter is functionals that provide a concise formulation for thermodynamic quantities and Green's function methods in interacting many-body systems. Building upon the instructive examples of density functional theory and the Hartree– Fock approximation, functionals of the Green's function are developed to provide a framework for practical methods used extensively in many-body perturbation theory (Chs. 9–15), dynamical mean-field theory (Chs. 16–20), and their combination in Ch. 21. Functional concepts give a firm foundation for those methods, for example, they give the conditions for approximations to obey conservation laws. Additional background and specific aspects are in App. H. Many physical quantities can be expressed as functionals , i.e., they depend on an entire function.1 For example, in the quantum variational method to determine a wavefunction, the total energy of a system of many electrons is a functional of a many-body trial wavefunction ( r 1 , r 2 , …): it is the expectation value of the hamiltonian. At the variational minimum, assuming the trial wavefunction is completely flexible, one finds the many-body Schrödinger equation for the ground state. Many of the properties we are interested in are, formally, simple functionals of the ground- or excited-state wavefunctions. In the search for feasible ways to calculate properties such as excitation spectra, we are led to different types of functional. The idea is to deal only with a physical, measurable quantity Q of interest and a variational functional E [ Q ] of this quantity. The variational character of E allows one to derive equations that determine Q 0 , the value of Q for the system in a desired state, typically the ground state or equilibrium state at T ≠ 0. Other properties should be different functionals F m [ Q ] that can then be evaluated at Q = Q 0 . This strategy is outlined in Fig. 8.1. To appreciate the power of such approaches, one can consider the impact of density functional theory, which is formulated as a functional of the density n ( r ) of electrons; it is defined for a range of densities and the variational minimum for the energy functional E [ n ] leads to equations to determine the actual ground-state density n 0 ( r ) and total energy E 0 .

Cambridge University Press eBooks•May 31, 2016•0 citations

Green functions in the renormalized many-body perturbation theory for correlated and disordered electrons

Janiš

The ways of introducing and handling renormalizations in the many-body perturbation theory are reviewed.We stress the indispensable role the technique of Green functions plays in extrapolating the weak-coupling perturbative approaches to intermediate and strong couplings.We separately discuss mass and charge renormalizations.The former is incorporated in a self-consistent equation for the self-energy derived explicitly from generating Feynman diagrams within the Baym and Kadanoff approach.The latter amounts to self-consistent equations for two-particle irreducible vertices.We analyze the charge renormalization initiated by De Dominicis and Martin and demonstrate that its realization via the parquet approach may become a powerful and viable way of using the many-body diagrammatic approach reliably in non-perturbative regimes with cooperative phenomena induced by either strong interaction or strong static randomness.

Condensed Matter Physics•Jan 1, 2006•7 citations
We derive variational expressions for the grand potential or action in terms of the many-body Green function GGG which describes the propagation of particles and the renormalized four-point vertex \ensuremathΓ\ensuremath{\Gamma}\ensuremathΓ which describes the scattering of two particles in many-body systems. The main ingredient of the variational functionals is a term we denote as the \ensuremathΞ\ensuremath{\Xi}\ensuremathΞ-functional which plays a role analogously to the usual \ensuremathΦ\ensuremath{\Phi}\ensuremathΦ-functional studied by Baym [G. Baym, Phys. Rev. 127, 1391 (1962)] in connection with the conservation laws in many-body systems. We show that any \ensuremathΞ\ensuremath{\Xi}\ensuremathΞ-derivable theory is also \ensuremathΦ\ensuremath{\Phi}\ensuremathΦ-derivable and therefore respects the conservation laws. We further set up a computational scheme to obtain accurate total energies from our variational functionals without having to solve computationally expensive sets of self-consistent equations. The input of the functional is an approximate Green function G\ifmmode~\else˜\fi\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{G}G\ifmmode~\else˜\fi​ and an approximate four-point vertex \ensuremathΓ\ifmmode~\else˜\fi\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{\ensuremath{\Gamma}}\ensuremathΓ\ifmmode~\else˜\fi​ obtained at a relatively low computational cost. The variational property of the functional guarantees that the error in the total energy is only of second order in deviations of the input Green function and vertex from the self-consistent ones that make the functional stationary. The functionals that we will consider for practical applications correspond to infinite order summations of ladder and exchange diagrams and are therefore particularly suited for applications to highly correlated systems. Their practical evaluation is discussed in detail.
I review the way the many-body Green functions are used to renormalize the perturbation theory of correlated fermions. The Green functions are introduced to implement systematically dynamical corrections to the static mean-field theory. The renormalizations enter the perturbation theory via self-consistent evaluations of one-particle and two-particle Green functions. The one-particle self-consistency is discussed within the Baym-Kadanoff construction. We point to an inherent dichotomy of the Baym-Kadanoff approach in the relation between one and two-particle Green functions. They are the Schwinger-Dyson equation, directly derived from the generating Luttinger-Ward functional, and the Ward identity. The latter imposes a severe restriction on the two-particle irreducible vertex and is not guaranteed in the Baym-Kadanoff approach. No approximate theory is capable to obey both relations. I further discuss the two-particle approach with two-particle renormalizations that offers a way to reconcile the Schwinger-Dyson equation and the Ward identity. The two-particle irreducible vertex from the critical Bethe-Salpeter equation becomes a generating functional in the two-particle renormalization of the perturbation theory. The irreducible vertex is used in the Ward identity to determine the order parameter beyond the critical point and in the Bethe-Salpeter equation representing the vertex of the Schwinger-Dyson equation to determine the spectral self-energy in such a way that the thermodynamic critical behavior from both equations is qualitatively the same. Finally, I discuss a static mean-field-like approximation with a renormalized effective interaction determined self-consistently that leads to a qualitatively correct description of strong electron correlations in impurity and lattice models.
Parquet diagrams sum self-consistently Feynman graphs for the vertex function with all two-particle multiple scatterings. We show how the complete parquet equations for the Hubbard-like models can be integrated to a generating functional from which all thermodynamic quantities are derived via (functional) derivatives. An explicit Luttinger-Ward functional Φ[G;Λ,K;U]Φ[G;Λ,{\cal K};U]Φ[G;Λ,K;U] is constructed containing only the renormalized one-particle, GGG, irreducible, ΛΛΛ, and reducible, K{\cal K}K, two-particle propagators as independent variational functions. The parquet approximation is proven to be a thermodynamically consistent, ΦΦΦ-derivable theory obeying the necessary conservation laws.
We develop a systematic approach to construct energy functionals of the one-particle reduced density matrix (1RDM) for equilibrium systems at finite temperature. The starting point of our formulation is the grand potential Ω [ G ] regarded as variational functional of the Green’s function G based on diagrammatic many-body perturbation theory and for which we consider either the Klein or Luttinger–Ward form. By restricting the input Green’s function to be one-to-one related to a set on one-particle reduced density matrices (1RDM) this functional becomes a functional of the 1RDM. To establish the one-to-one mapping we use that, at any finite temperature and for a given 1RDM γ in a finite basis, there exists a non-interacting system with a spatially non-local potential v [ γ ] which reproduces the given 1RDM. The corresponding set of non-interacting Green’s functions defines the variational domain of the functional Ω . In the zero temperature limit we obtain an energy functional E [ γ ] which by minimisation yields an approximate ground state 1RDM and energy. As an application of the formalism we use the Klein and Luttinger–Ward functionals in the GW-approximation to compute the binding curve of a model hydrogen molecule using an extended Hubbard Hamiltonian. We compare further to the case in which we evaluate the functionals on a Hartree–Fock and a Kohn–Sham Green’s function. We find that the Luttinger–Ward version of the functionals performs the best and is able to reproduce energies close to the GW energy which corresponds to the stationary point.
Abstract Chapter 09 Nonequilibrium Green’s functions (NEGF), including coupled-correlated (C) single- and multi-particle Green’s functions, are defined as averages weighted with the time-development operator U(t0+τ,t0). Linear conductivity is exhibited as a two-particle equilibrium Green’s function (Kubo-type formulation). Admitting particle sources (S:η,η+) and non-conservation of number, the non-equilibrium multi-particle Green’s functions are constructed with numbers of creation and annihilation operators that may differ, and they may be derived as variational derivatives with respect to sources η,η+ of a generating functional eW=TrU(t0+τ,t0)CS/TrU(t0+τ,t0)C. (In the non-interacting case this yields the n-particle Green’s function as a permanent/determinant of single-particle Green’s functions.) These variational relations yield a symmetric set of multi-particle Green’s function equations. Cumulants and the Linked Cluster Theorem are discussed and the Random Phase Approximation (RPA) is derived variationally. Schwinger’s variational differential formulation of perturbation theories for the Green’s function, self-energy, vertex operator, and also shielded potential perturbation theory, are reviewed. The Langreth Algebra arises from analytic continuation of integration of products of Green’s functions in imaginary time to the real-time axis with time-ordering along the integration contour in the complex time plane. An account of the Generalized Kadanoff-Baym Ansatz is presented.
To avoid confusion, it is necessary to be quite explicit about what is assumed and what is to be proved. A. Klein, Phys. Rev. 121, 950 (1961) Summary The topic of this chapter is functionals that provide a concise formulation for thermodynamic quantities and Green's function methods in interacting many-body systems. Building upon the instructive examples of density functional theory and the Hartree– Fock approximation, functionals of the Green's function are developed to provide a framework for practical methods used extensively in many-body perturbation theory (Chs. 9–15), dynamical mean-field theory (Chs. 16–20), and their combination in Ch. 21. Functional concepts give a firm foundation for those methods, for example, they give the conditions for approximations to obey conservation laws. Additional background and specific aspects are in App. H. Many physical quantities can be expressed as functionals , i.e., they depend on an entire function.1 For example, in the quantum variational method to determine a wavefunction, the total energy of a system of many electrons is a functional of a many-body trial wavefunction ( r 1 , r 2 , …): it is the expectation value of the hamiltonian. At the variational minimum, assuming the trial wavefunction is completely flexible, one finds the many-body Schrödinger equation for the ground state. Many of the properties we are interested in are, formally, simple functionals of the ground- or excited-state wavefunctions. In the search for feasible ways to calculate properties such as excitation spectra, we are led to different types of functional. The idea is to deal only with a physical, measurable quantity Q of interest and a variational functional E [ Q ] of this quantity. The variational character of E allows one to derive equations that determine Q 0 , the value of Q for the system in a desired state, typically the ground state or equilibrium state at T ≠ 0. Other properties should be different functionals F m [ Q ] that can then be evaluated at Q = Q 0 . This strategy is outlined in Fig. 8.1. To appreciate the power of such approaches, one can consider the impact of density functional theory, which is formulated as a functional of the density n ( r ) of electrons; it is defined for a range of densities and the variational minimum for the energy functional E [ n ] leads to equations to determine the actual ground-state density n 0 ( r ) and total energy E 0 .
The ways of introducing and handling renormalizations in the many-body perturbation theory are reviewed.We stress the indispensable role the technique of Green functions plays in extrapolating the weak-coupling perturbative approaches to intermediate and strong couplings.We separately discuss mass and charge renormalizations.The former is incorporated in a self-consistent equation for the self-energy derived explicitly from generating Feynman diagrams within the Baym and Kadanoff approach.The latter amounts to self-consistent equations for two-particle irreducible vertices.We analyze the charge renormalization initiated by De Dominicis and Martin and demonstrate that its realization via the parquet approach may become a powerful and viable way of using the many-body diagrammatic approach reliably in non-perturbative regimes with cooperative phenomena induced by either strong interaction or strong static randomness.