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:1606.07486 · Physical review. B./Physical review. B

Partial self-consistency and analyticity in many-body perturbation theory: Particle number conservation and a generalized sum rule

Daniel Karlsson, Robert van Leeuwen

Published Sep 13, 2016• 14 citations
arXivHTMLPDF

Abstract

We consider a general class of approximations which guarantees the conservation of particle number in many-body perturbation theory. To do this we extend the concept of $\mathrm{\ensuremath{\Phi}}$ derivability for the self-energy $\mathrm{\ensuremath{\Sigma}}$ to a larger class of diagrammatic terms in which only some of the Green's function lines contain the fully dressed Green's function $G$. We call the corresponding approximations for $\mathrm{\ensuremath{\Sigma}}$ partially $\mathrm{\ensuremath{\Phi}}$ derivable. A special subclass of such approximations, which are gauge invariant, is obtained by dressing loops in the diagrammatic expansion of $\mathrm{\ensuremath{\Phi}}$ consistently with $G$. These approximations are number conserving but do not have to fulfill other conservation laws, such as the conservation of energy and momentum. From our formalism we can easily deduce whether commonly used approximations will fulfill the continuity equation, which implies particle number conservation. We further show how the concept of partial $\mathrm{\ensuremath{\Phi}}$ derivability plays an important role in the derivation of a generalized sum rule for the particle number, which reduces to the Luttinger-Ward theorem in the case of a homogeneous electron gas, and the Friedel sum rule in the case of the Anderson model. To do this we need to ensure that the Green's function has certain complex analytic properties, which can be guaranteed if the spectral function is positive-semidefinite. The latter property can be ensured for a subset of partially $\mathrm{\ensuremath{\Phi}}$-derivable approximations for the self-energy, namely those that can be constructed from squares of so-called half diagrams. For the case in which the analytic requirements are not fulfilled we highlight a number of subtle issues related to branch cuts, pole structure, and multivaluedness. We also show that various schemes of computing the particle number are consistent for particle number conserving approximations.

DOI
10.1103/physrevb.94.125124

Related papers

Diagrammatic self-energy approximations and the total particle number

Arno Schindlmayr, Pablo García‐González, Rex W. Godby

There is increasing interest in many-body perturbation theory as a practical tool for the calculation of ground-state properties. As a consequence, unambiguous sum rules such as the conservation of particle number under the influence of the Coulomb interaction have acquired an importance that did not exist for calculations of excited-state properties. In this paper we obtain a rigorous, simple relation whose fulfilment guarantees particle-number conservation in a given diagrammatic self-energy approximation. Hedin's ${G}_{0}{W}_{0}$ approximation does not satisfy this relation and hence violates the particle-number sum rule. Very precise calculations for the homogeneous electron gas and a model inhomogeneous electron system allow the extent of the nonconservation to be estimated.

arXiv:cond-mat/0110435•Physical review. B, Condensed matter•Nov 19, 2001•38 citations

Self-Consistent Approximations in Many-Body Systems

Gordon Baym

This paper investigates the criteria for maintenance of the macroscopic conservation laws of number, momentum, and energy by approximate two-particle correlation functions in many-body systems. The methods of generating such approximations are the same as in a previous paper. However, the derivations of the conservation laws given here clarify both why the approximation method works and the connection between the macroscopic conservation laws and those at the vertices.Conserving nonequilibrium approximations are based on self-consistent approximations to the one-particle Green's function. The same condition that ensures that the nonequilibrium theory be conserving also ensures that the equilibrium approximation has the following properties. The several common methods for determining the partition function from the one-particle Green's function all lead to the same result. When applied to a zero-temperature normal fermion system, the approximation procedure maintains the Hugenholtz-Van Hove theorem. Consequently, the self-consistent version of Brueckner's nuclear matter theory obeys this theorem.

Physical Review•Aug 15, 1962•1,367 citations

Study of linear response in Hubbard chains using Many-body Perturbation Theory

Irina Bostan

In this work the basic formalism of non-equilibrium Green’s functions is presented and then applied to study a Ward identity in linear response theory, namely the frequency sum-rule. It can be proven that the frequency sum-rule is satisfied when the quantities involved are calculated using perturbation theory within a conserving approximation for the self-energy. To illustrate this equality along with other properties of the response function, a numerical application that solves the Kadanoff-Baym equations for systems of Hubbard chains was used. The results showed that the frequency sum-rule was satisfied to the same extent by all the conserving approximations used as by the exact diagonalization numerical results. The density response function was analyzed diagrammatically for a series of conserving approximations for the self-energy and this demonstrated that even for a first order in perturbation theory approximation for the self-energy, the response function has a corresponding complex, third order in the perturbation diagrammatic structure.

•Jan 1, 2010•0 citations

Potential energy surfaces from many-body functionals: Analytical benchmarks and conserving many-body approximations

Giovanna Lani, Nicola Marzari

We investigate analytically the performance of many-body energy functionals, derived by Klein, Luttinger, and Ward, at different levels of diagrammatic approximations, ranging from second Born, to , to the so-called T matrix, for the calculation of total energies and potential energy surfaces. We benchmark our theoretical results on the extended two-site Hubbard model, which is analytically solvable and for which several exact properties can be calculated. Despite its simplicity, this model displays the physics of strongly correlated electrons: it is prototypical of the H 2 dissociation, a notoriously difficult problem to solve accurately for the majority of mean-field-based approaches. We show that both functionals exhibit good to excellent variational properties, particularly in the case of the Luttinger-Ward one, which is in close agreement with fully self-consistent calculations, and we elucidate the relation between the accuracy of the results and the different input one-body Green's functions. Provided that these are wisely chosen, we show how the Luttinger-Ward functional can be used as a computationally inexpensive alternative to fully self-consistent many-body calculations, without sacrificing the precision of the results obtained. Furthermore, in virtue of this accuracy, we argue that this functional can also be used to rank different many-body approximations at different regimes of electronic correlation, once again bypassing the need for self-consistency. Published by the American Physical Society 2024

arXiv:2404.16453•Physical Review Research•Dec 23, 2024•2 citations

New approach to the one-particle Green's function for finite Fermi systems

J. Schirmer, Lorenz S. Cederbaum, Olaf Walter

A new approach to the one-particle Green's functions $\mathit{G}$ for finite electronic systems is presented. This approach is based on the diagrammatic perturbation expansions of the Green's function and of the dynamic self-energy part $\mathit{M}$ related to $\mathit{G}$ via the Dyson equation. The exact summation of the latter expansion is reformulated in terms of a simple algebraic form referred to as algebraic diagrammatic construction (ADC). The ADC defines in a systematical way a set of approximation schemes ($n\mathrm{th}$-order ADC schemes) that represent infinite partial summations for $\mathit{M}$ and (via the Dyson equation) for $\mathit{G}$ being complete through $n\mathrm{th}$ order of perturbation theory. The corresponding mathematical procedures are essentially Hermitian eigenvalue problems in restricted configuration spaces of unperturbed ionic configurations. Explicit equations for the second-, third-, and fourth-order ADC schemes are derived and analyzed. While the second- and third-order schemes can be viewed as systematic rederivations of previous approximation schemes, the fourth-order ADC scheme represents a complete fourth-order approximation for the self-energy and the one-particle Green's function which was hitherto not available.

Physical Review A•Sep 1, 1983•725 citations

Nonconvergence of the Feynman-Dyson diagrammatic perturbation expansion of propagators

So Hirata, Ireneusz Grabowski, J. V. Ortiz +1 more

Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (ii) An odd-perturbation-order self-energy has a qualitatively wrong shape and, as a result, many roots of the corresponding Dyson equation are nonphysical in that the poles may be complex or residues can exceed unity or be negative. (iii) A higher even-order self-energy consists of vertical lines at many frequencies, predicting numerous phantom poles with zero residues. (iv) Infinite partial resummations of diagrams by vertex or edge renormalization tend to exacerbate these pathologies. (v) The nonconvergence is caused by the nonanalyticity of the rational-function form of the exact Green's function at many frequencies, where the radius of convergence of its Taylor expansion is zero. This is consistent with the fact that (vi) Pad\'e approximants (power-series expansions of a rational function) can largely restore the correct shape and poles of the Green's function. Nevertheless, not only does the nonconvergence render higher-order Feynman-Dyson diagrammatic perturbation theory useless for many lower-lying ionization or higher-lying electron-attachment states but it also calls into question the validity of its combined use with the Ans\"atze requiring the knowledge of all poles and residues. Such Ans\"atze include the Galitskii-Migdal identity, the self-consistent Green's function methods, and some models of the algebraic diagrammatic construction.

arXiv:2312.03157•Physical Review A•May 16, 2024•11 citations
We consider a general class of approximations which guarantees the conservation of particle number in many-body perturbation theory. To do this we extend the concept of \ensuremathΦ\mathrm{\ensuremath{\Phi}}\ensuremathΦ derivability for the self-energy \ensuremathΣ\mathrm{\ensuremath{\Sigma}}\ensuremathΣ to a larger class of diagrammatic terms in which only some of the Green's function lines contain the fully dressed Green's function GGG. We call the corresponding approximations for \ensuremathΣ\mathrm{\ensuremath{\Sigma}}\ensuremathΣ partially \ensuremathΦ\mathrm{\ensuremath{\Phi}}\ensuremathΦ derivable. A special subclass of such approximations, which are gauge invariant, is obtained by dressing loops in the diagrammatic expansion of \ensuremathΦ\mathrm{\ensuremath{\Phi}}\ensuremathΦ consistently with GGG. These approximations are number conserving but do not have to fulfill other conservation laws, such as the conservation of energy and momentum. From our formalism we can easily deduce whether commonly used approximations will fulfill the continuity equation, which implies particle number conservation. We further show how the concept of partial \ensuremathΦ\mathrm{\ensuremath{\Phi}}\ensuremathΦ derivability plays an important role in the derivation of a generalized sum rule for the particle number, which reduces to the Luttinger-Ward theorem in the case of a homogeneous electron gas, and the Friedel sum rule in the case of the Anderson model. To do this we need to ensure that the Green's function has certain complex analytic properties, which can be guaranteed if the spectral function is positive-semidefinite. The latter property can be ensured for a subset of partially \ensuremathΦ\mathrm{\ensuremath{\Phi}}\ensuremathΦ-derivable approximations for the self-energy, namely those that can be constructed from squares of so-called half diagrams. For the case in which the analytic requirements are not fulfilled we highlight a number of subtle issues related to branch cuts, pole structure, and multivaluedness. We also show that various schemes of computing the particle number are consistent for particle number conserving approximations.
There is increasing interest in many-body perturbation theory as a practical tool for the calculation of ground-state properties. As a consequence, unambiguous sum rules such as the conservation of particle number under the influence of the Coulomb interaction have acquired an importance that did not exist for calculations of excited-state properties. In this paper we obtain a rigorous, simple relation whose fulfilment guarantees particle-number conservation in a given diagrammatic self-energy approximation. Hedin's G0W0{G}_{0}{W}_{0}G0​W0​ approximation does not satisfy this relation and hence violates the particle-number sum rule. Very precise calculations for the homogeneous electron gas and a model inhomogeneous electron system allow the extent of the nonconservation to be estimated.
This paper investigates the criteria for maintenance of the macroscopic conservation laws of number, momentum, and energy by approximate two-particle correlation functions in many-body systems. The methods of generating such approximations are the same as in a previous paper. However, the derivations of the conservation laws given here clarify both why the approximation method works and the connection between the macroscopic conservation laws and those at the vertices.Conserving nonequilibrium approximations are based on self-consistent approximations to the one-particle Green's function. The same condition that ensures that the nonequilibrium theory be conserving also ensures that the equilibrium approximation has the following properties. The several common methods for determining the partition function from the one-particle Green's function all lead to the same result. When applied to a zero-temperature normal fermion system, the approximation procedure maintains the Hugenholtz-Van Hove theorem. Consequently, the self-consistent version of Brueckner's nuclear matter theory obeys this theorem.
In this work the basic formalism of non-equilibrium Green’s functions is presented and then applied to study a Ward identity in linear response theory, namely the frequency sum-rule. It can be proven that the frequency sum-rule is satisfied when the quantities involved are calculated using perturbation theory within a conserving approximation for the self-energy. To illustrate this equality along with other properties of the response function, a numerical application that solves the Kadanoff-Baym equations for systems of Hubbard chains was used. The results showed that the frequency sum-rule was satisfied to the same extent by all the conserving approximations used as by the exact diagonalization numerical results. The density response function was analyzed diagrammatically for a series of conserving approximations for the self-energy and this demonstrated that even for a first order in perturbation theory approximation for the self-energy, the response function has a corresponding complex, third order in the perturbation diagrammatic structure.
We investigate analytically the performance of many-body energy functionals, derived by Klein, Luttinger, and Ward, at different levels of diagrammatic approximations, ranging from second Born, to , to the so-called T matrix, for the calculation of total energies and potential energy surfaces. We benchmark our theoretical results on the extended two-site Hubbard model, which is analytically solvable and for which several exact properties can be calculated. Despite its simplicity, this model displays the physics of strongly correlated electrons: it is prototypical of the H 2 dissociation, a notoriously difficult problem to solve accurately for the majority of mean-field-based approaches. We show that both functionals exhibit good to excellent variational properties, particularly in the case of the Luttinger-Ward one, which is in close agreement with fully self-consistent calculations, and we elucidate the relation between the accuracy of the results and the different input one-body Green's functions. Provided that these are wisely chosen, we show how the Luttinger-Ward functional can be used as a computationally inexpensive alternative to fully self-consistent many-body calculations, without sacrificing the precision of the results obtained. Furthermore, in virtue of this accuracy, we argue that this functional can also be used to rank different many-body approximations at different regimes of electronic correlation, once again bypassing the need for self-consistency. Published by the American Physical Society 2024
A new approach to the one-particle Green's functions G\mathit{G}G for finite electronic systems is presented. This approach is based on the diagrammatic perturbation expansions of the Green's function and of the dynamic self-energy part M\mathit{M}M related to G\mathit{G}G via the Dyson equation. The exact summation of the latter expansion is reformulated in terms of a simple algebraic form referred to as algebraic diagrammatic construction (ADC). The ADC defines in a systematical way a set of approximation schemes (nthn\mathrm{th}nth-order ADC schemes) that represent infinite partial summations for M\mathit{M}M and (via the Dyson equation) for G\mathit{G}G being complete through nthn\mathrm{th}nth order of perturbation theory. The corresponding mathematical procedures are essentially Hermitian eigenvalue problems in restricted configuration spaces of unperturbed ionic configurations. Explicit equations for the second-, third-, and fourth-order ADC schemes are derived and analyzed. While the second- and third-order schemes can be viewed as systematic rederivations of previous approximation schemes, the fourth-order ADC scheme represents a complete fourth-order approximation for the self-energy and the one-particle Green's function which was hitherto not available.
Using a general-order ab initio many-body Green's function method, we numerically illustrate several pathological behaviors of the Feynman-Dyson diagrammatic perturbation expansion of one-particle many-body Green's functions as electron Feynman propagators. (i) The perturbation expansion of the frequency-dependent self-energy is not convergent at the exact self-energy in many frequency domains. (ii) An odd-perturbation-order self-energy has a qualitatively wrong shape and, as a result, many roots of the corresponding Dyson equation are nonphysical in that the poles may be complex or residues can exceed unity or be negative. (iii) A higher even-order self-energy consists of vertical lines at many frequencies, predicting numerous phantom poles with zero residues. (iv) Infinite partial resummations of diagrams by vertex or edge renormalization tend to exacerbate these pathologies. (v) The nonconvergence is caused by the nonanalyticity of the rational-function form of the exact Green's function at many frequencies, where the radius of convergence of its Taylor expansion is zero. This is consistent with the fact that (vi) Pad\'e approximants (power-series expansions of a rational function) can largely restore the correct shape and poles of the Green's function. Nevertheless, not only does the nonconvergence render higher-order Feynman-Dyson diagrammatic perturbation theory useless for many lower-lying ionization or higher-lying electron-attachment states but it also calls into question the validity of its combined use with the Ans\"atze requiring the knowledge of all poles and residues. Such Ans\"atze include the Galitskii-Migdal identity, the self-consistent Green's function methods, and some models of the algebraic diagrammatic construction.