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:2212.14768

Misleading convergence of the skeleton diagrammatic technique: when the correct solution can be found

Aaram J. Kim, Evgeny Kozik

Published Dec 30, 2022• 3 citations
arXivHTMLPDF

Abstract

Convergence of the self-consistent skeleton diagrammatic technique (SDT) -- in which the full Green's function is determined through summation of Feynman diagrams in terms of itself -- to the wrong answer has been associated with the existence of non-perturbative branches of the Luttinger-Ward functional. Although it has been possible to detect misleading convergence without the knowledge of the exact result, the SDT has remained inapplicable in the regimes where this happens. We show that misleading convergence does not always preclude recovering the exact solution. In addition to the established mechanism, convergence of the SDT to the wrong answer can stem from divergence of the inherent diagrammatic series, which allows us to recover the exact solution by a modified SDT protocol based on controlled analytic continuation. We illustrate this approach by its application to the analytically solvable (0+0)d Hubbard model, the Hubbard atom, and the 2d Hubbard model in a challenging strong-coupling regime, for which the SDT is solved with controlled accuracy by the diagrammatic Monte Carlo (DiagMC) technique.

DOI
10.48550/arxiv.2212.14768

Related papers

Nonexistence of the Luttinger-Ward Functional and Misleading Convergence of Skeleton Diagrammatic Series for Hubbard-Like Models

Evgeny Kozik, Michel Ferrero, Antoine Georges

The Luttinger-Ward functional $\mathrm{\ensuremath{\Phi}}[\mathbf{G}]$, which expresses the thermodynamic grand potential in terms of the interacting single-particle Green's function $\mathbf{G}$, is found to be ill defined for fermionic models with the Hubbard on-site interaction. In particular, we show that the self-energy $\mathbf{\ensuremath{\Sigma}}[\mathbf{G}]\ensuremath{\propto}\ensuremath{\delta}\mathrm{\ensuremath{\Phi}}[\mathbf{G}]/\ensuremath{\delta}\mathbf{G}$ is not a single-valued functional of $\mathbf{G}$: in addition to the physical solution for $\mathbf{\ensuremath{\Sigma}}[\mathbf{G}]$, there exists at least one qualitatively distinct unphysical branch. This result is demonstrated for several models: the Hubbard atom, the Anderson impurity model, and the full two-dimensional Hubbard model. Despite this pathology, the skeleton Feynman diagrammatic series for $\mathbf{\ensuremath{\Sigma}}$ in terms of $\mathbf{G}$ is found to converge at least for moderately low temperatures. However, at strong interactions, its convergence is to the unphysical branch. This reveals a new scenario of breaking down of diagrammatic expansions. In contrast, the bare series in terms of the noninteracting Green's function ${\mathbf{G}}_{0}$ converges to the correct physical branch of $\mathbf{\ensuremath{\Sigma}}$ in all cases currently accessible by diagrammatic Monte Carlo calculations. In addition to their conceptual importance, these observations have important implications for techniques based on the explicit summation of the diagrammatic series.

arXiv:1407.5687•Physical Review Letters•Apr 15, 2015•148 citations

Comment on `Non-existence of the Luttinger-Ward functional and misleading convergence of skeleton diagrammatic series for Hubbard-like models'

R. Eder

The known analytical properties of the Green's function and self-energy rule out an ambiguity of the self-energy. The noninteracting Green's functions obtained by Kozik et al. in arXiv:1407.5687 likely have pathological properties.

arXiv:1407.6599•Jul 24, 2014•3 citations

On the Luttinger-Ward functional and the convergence of skeleton diagrammatic series expansion of the self-energy for Hubbard-like models

Behnam Farid

We consider a number of questions regarding the Luttinger-Ward functional and the many-body perturbation series expansion of the proper self-energy $Σ(\mathbf{k};z)$ specific to uniform ground states (ensemble of states) of interacting fermion systems in terms of skeleton self-energy diagrams and the interacting Green function $G(\mathbf{k};z)$. Utilising a link between the latter series expansion and the classical moment problem (of the Hamburger type), along with the associated continued-fraction expansion, we reaffirm our earlier observation (2007) that for lattice models of fermions interacting through short-range two-body potentials (i.e. for Hubbard-like models) this series is uniformly convergent for almost all wave vectors $\mathbf{k}$ and complex energies $z$. The limit of this series is unique. We inquire into the reasons underlying the contrary observation by Kozik et al. (2015) regarding skeleton-diagrammatic perturbation series expansion of $Σ$. In doing so, we make a number of observations of general interest. Chief amongst these, we observe that contrary to general belief thermal correlation functions calculated in the energy domain (i.e. over the set of the relevant Matsubara frequencies) are generally unreliable, in contrast to those calculated in the imaginary-time domain. In an appendix we present a symbolic computational formalism and short programs for singling out topologically distinct proper self-energy diagrams that are algebraically equivalent up to determinate multiplicative constants. This complements the work presented in our previous publication (2019). [Abridged abstract]

arXiv:2108.10903•Aug 24, 2021•1 citation

Bold Feynman diagrams and the Luttinger-Ward formalism via Gibbs measures. Part II: Non-perturbative analysis

Lin Lin, Michael Lindsey

Many-body perturbation theory (MBPT) is widely used in quantum physics, chemistry, and materials science. At the heart of MBPT is the Feynman diagrammatic expansion, which is, simply speaking, an elegant way of organizing the combinatorially growing number of terms of a certain Taylor expansion. In particular, the construction of the `bold Feynman diagrammatic expansion' involves the partial resummation to infinite order of possibly divergent series of diagrams. This procedure demands investigation from both the combinatorial (perturbative) and the analytical (non-perturbative) viewpoints. In Part II of this two-part series, we approach the analytical investigation of the bold diagrammatic expansion in the simplified setting of Gibbs measures (known as the Euclidean lattice field theory in the physics literature). Using non-perturbative methods, we rigorously construct the Luttinger-Ward formalism for the first time, and we prove that the bold diagrammatic series can be obtained directly via an asymptotic expansion of the Luttinger-Ward functional, circumventing the partial resummation technique. Moreover we prove that the Dyson equation can be derived as the Euler-Lagrange equation associated with a variational problem involving the Luttinger-Ward functional. We also establish a number of key facts about the Luttinger-Ward functional, such as its transformation rule, its form in the setting of the impurity problem, and its continuous extension to the boundary of the domain of physical Green's functions.

arXiv:1809.02901•Sep 9, 2018•1 citation

Diagrammatic Monte Carlo approach for diagrammatic extensions of dynamical mean-field theory: Convergence analysis of the dual fermion technique

Jan Gukelberger, Evgeny Kozik, Hartmut Hafermann

The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.

arXiv:1611.07523•Physical review. B./Physical review. B•Jul 28, 2017•55 citations

Homotopic Action: A Pathway to Convergent Diagrammatic Theories

Aaram J. Kim, Nikolay V. Prokof'ev, Boris V. Svistunov +1 more

The major obstacle preventing Feynman diagrammatic expansions from accurately solving many-fermion systems in strongly correlated regimes is the series slow convergence or divergence problem. Several techniques have been proposed to address this issue: series resummation by conformal mapping, changing the nature of the starting point of the expansion by shifted action tools, and applying the homotopy analysis method to the Dyson-Schwinger equation. They emerge as dissimilar mathematical procedures aimed at different aspects of the problem. The proposed homotopic action offers a universal and systematic framework for unifying the existing-and generating new-methods and ideas to formulate a physical system in terms of a convergent diagrammatic series. It eliminates the need for resummation, allows one to introduce effective interactions, enables a controlled ultraviolet regularization of continuous-space theories, and reduces the intrinsic polynomial complexity of the diagrammatic Monte Carlo method. We illustrate this approach by an application to the Hubbard model.

arXiv:2010.05301•Physical Review Letters•Jun 23, 2021•18 citations
Convergence of the self-consistent skeleton diagrammatic technique (SDT) -- in which the full Green's function is determined through summation of Feynman diagrams in terms of itself -- to the wrong answer has been associated with the existence of non-perturbative branches of the Luttinger-Ward functional. Although it has been possible to detect misleading convergence without the knowledge of the exact result, the SDT has remained inapplicable in the regimes where this happens. We show that misleading convergence does not always preclude recovering the exact solution. In addition to the established mechanism, convergence of the SDT to the wrong answer can stem from divergence of the inherent diagrammatic series, which allows us to recover the exact solution by a modified SDT protocol based on controlled analytic continuation. We illustrate this approach by its application to the analytically solvable (0+0)d Hubbard model, the Hubbard atom, and the 2d Hubbard model in a challenging strong-coupling regime, for which the SDT is solved with controlled accuracy by the diagrammatic Monte Carlo (DiagMC) technique.
The Luttinger-Ward functional $\mathrm{\ensuremath{\Phi}}[\mathbf{G}]$, which expresses the thermodynamic grand potential in terms of the interacting single-particle Green's function G\mathbf{G}G, is found to be ill defined for fermionic models with the Hubbard on-site interaction. In particular, we show that the self-energy $\mathbf{\ensuremath{\Sigma}}[\mathbf{G}]\ensuremath{\propto}\ensuremath{\delta}\mathrm{\ensuremath{\Phi}}[\mathbf{G}]/\ensuremath{\delta}\mathbf{G}$ is not a single-valued functional of G\mathbf{G}G: in addition to the physical solution for $\mathbf{\ensuremath{\Sigma}}[\mathbf{G}]$, there exists at least one qualitatively distinct unphysical branch. This result is demonstrated for several models: the Hubbard atom, the Anderson impurity model, and the full two-dimensional Hubbard model. Despite this pathology, the skeleton Feynman diagrammatic series for $\mathbf{\ensuremath{\Sigma}}$ in terms of G\mathbf{G}G is found to converge at least for moderately low temperatures. However, at strong interactions, its convergence is to the unphysical branch. This reveals a new scenario of breaking down of diagrammatic expansions. In contrast, the bare series in terms of the noninteracting Green's function G0{\mathbf{G}}_{0}G0​ converges to the correct physical branch of $\mathbf{\ensuremath{\Sigma}}$ in all cases currently accessible by diagrammatic Monte Carlo calculations. In addition to their conceptual importance, these observations have important implications for techniques based on the explicit summation of the diagrammatic series.
We consider a number of questions regarding the Luttinger-Ward functional and the many-body perturbation series expansion of the proper self-energy Σ(k;z)Σ(\mathbf{k};z)Σ(k;z) specific to uniform ground states (ensemble of states) of interacting fermion systems in terms of skeleton self-energy diagrams and the interacting Green function G(k;z)G(\mathbf{k};z)G(k;z). Utilising a link between the latter series expansion and the classical moment problem (of the Hamburger type), along with the associated continued-fraction expansion, we reaffirm our earlier observation (2007) that for lattice models of fermions interacting through short-range two-body potentials (i.e. for Hubbard-like models) this series is uniformly convergent for almost all wave vectors k\mathbf{k}k and complex energies zzz. The limit of this series is unique. We inquire into the reasons underlying the contrary observation by Kozik et al. (2015) regarding skeleton-diagrammatic perturbation series expansion of ΣΣΣ. In doing so, we make a number of observations of general interest. Chief amongst these, we observe that contrary to general belief thermal correlation functions calculated in the energy domain (i.e. over the set of the relevant Matsubara frequencies) are generally unreliable, in contrast to those calculated in the imaginary-time domain. In an appendix we present a symbolic computational formalism and short programs for singling out topologically distinct proper self-energy diagrams that are algebraically equivalent up to determinate multiplicative constants. This complements the work presented in our previous publication (2019). [Abridged abstract]
Many-body perturbation theory (MBPT) is widely used in quantum physics, chemistry, and materials science. At the heart of MBPT is the Feynman diagrammatic expansion, which is, simply speaking, an elegant way of organizing the combinatorially growing number of terms of a certain Taylor expansion. In particular, the construction of the `bold Feynman diagrammatic expansion' involves the partial resummation to infinite order of possibly divergent series of diagrams. This procedure demands investigation from both the combinatorial (perturbative) and the analytical (non-perturbative) viewpoints. In Part II of this two-part series, we approach the analytical investigation of the bold diagrammatic expansion in the simplified setting of Gibbs measures (known as the Euclidean lattice field theory in the physics literature). Using non-perturbative methods, we rigorously construct the Luttinger-Ward formalism for the first time, and we prove that the bold diagrammatic series can be obtained directly via an asymptotic expansion of the Luttinger-Ward functional, circumventing the partial resummation technique. Moreover we prove that the Dyson equation can be derived as the Euler-Lagrange equation associated with a variational problem involving the Luttinger-Ward functional. We also establish a number of key facts about the Luttinger-Ward functional, such as its transformation rule, its form in the setting of the impurity problem, and its continuous extension to the boundary of the domain of physical Green's functions.
The dual fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations, however, neglect higher-order interaction vertices beyond two-particle scattering in the dual effective action and further truncate the diagrammatic expansion in the two-particle scattering vertex to a leading-order or ladder-type approximation. In this work, we compute the dual fermion expansion for the two-dimensional Hubbard model including all diagram topologies with two-particle interactions to high orders by means of a stochastic diagrammatic Monte Carlo algorithm. We benchmark the obtained self-energy against numerically exact diagrammatic determinant Monte Carlo simulations to systematically assess convergence of the dual fermion series and the validity of these approximations. We observe that, from high temperatures down to the vicinity of the DMFT N\'eel transition, the dual fermion series converges very quickly to the exact solution in the whole range of Hubbard interactions considered ($4\ensuremath{\le}U/t\ensuremath{\le}12$), implying that contributions from higher-order vertices are small. As the temperature is lowered further, we observe slower series convergence, convergence to incorrect solutions, and ultimately divergence. This happens in a regime where magnetic correlations become significant. We find, however, that the self-consistent particle-hole ladder approximation yields reasonable and often even highly accurate results in this regime.
The major obstacle preventing Feynman diagrammatic expansions from accurately solving many-fermion systems in strongly correlated regimes is the series slow convergence or divergence problem. Several techniques have been proposed to address this issue: series resummation by conformal mapping, changing the nature of the starting point of the expansion by shifted action tools, and applying the homotopy analysis method to the Dyson-Schwinger equation. They emerge as dissimilar mathematical procedures aimed at different aspects of the problem. The proposed homotopic action offers a universal and systematic framework for unifying the existing-and generating new-methods and ideas to formulate a physical system in terms of a convergent diagrammatic series. It eliminates the need for resummation, allows one to introduce effective interactions, enables a controlled ultraviolet regularization of continuous-space theories, and reduces the intrinsic polynomial complexity of the diagrammatic Monte Carlo method. We illustrate this approach by an application to the Hubbard model.