K. Van Houcke,E. Kozik,R. Rossi,Y. Deng,F. Werner
Abstract:In the standard framework of self-consistent many-body perturbation theory, the skeleton series for the self-energy is truncated at a finite order $\mathcal{N}$ and plugged into the Dyson equation, which is then solved for the propagator $G_{\mathcal{N}}\,$. We consider two examples of fermionic models, the Hubbard atom at half filling and its zero space-time dimensional simplified version. First, we show that $G_{\mathcal{N}}\,$ converges when $\mathcal{N}\to\infty$ to a limit $G_\infty\,$, which coincides with the exact physical propagator $G_{\rm exact}\ $ at small enough coupling, while $G_\infty \neq G_{\rm exact}\ $ at strong coupling. This follows from the findings of [Kozik, Ferrero and Georges, PRL 114, 156402 (2015)] and an additional subtle mathematical mechanism elucidated here. Second, we demonstrate that it is possible to discriminate between the $G_\infty=G_{\rm exact}\ $ and $G_\infty\neq G_{\rm exact}\ $ regimes thanks to a criterion which does not require the knowledge of $G_{\rm exact}\ $, as proposed in [Rossi et al., PRB 93, 161102(R) (2016)].
What problem does this paper attempt to address?
This paper mainly discusses the problems in the physical and non-physical regions of self-consistent many-body perturbation theory. In the standard framework, the skeleton order of the self-energy is truncated at a finite order N and is substituted into the Dyson equation to solve the propagator GN. The researchers demonstrated through two fermionic models (half-filled Hubbard atom and its simplified version in zero temporal dimension) that as N approaches infinity, GN converges to a limit G∞, but this limit is not always the same as the exact physical propagator Gexact. In the case of strong coupling, G∞ is different from Gexact.
The article points out that although Σ(≤N)bold[Gexact] (the truncation and perturbation of the skeleton order on Gexact) may misleadingly converge to results different from the exact self-energy, the misleading convergence of GN is a non-trivial fact, stemming from a subtle mathematical mechanism. In a specific model, the GN sequence converges to incorrect results as N approaches infinity, which occurs when the free propagator exceeds a critical coupling strength.
To detect this misleading convergence, the paper proposes a criterion, which is to check whether the sequence of ΣN,ξ (the ξ power series of Σ truncated on G) remains convergent as N approaches infinity. If there exists ξ>1 such that ΣN,ξ converges, then the results can be considered reliable. For toy models and Hubbard atoms, this criterion successfully identifies cases of misleading convergence and verifies the results in correct cases.
In conclusion, the paper reveals that in self-consistent many-body theory, even when the skeleton order is expanded infinitely, the solution may not converge to the exact physical result but to non-physical results, especially in the case of strong coupling. In addition, a numerical diagnostic tool is proposed to detect such cases, which is of great significance in dealing with practical physical problems with unknown Gexact.