While self-consistent diagrammatic approaches are widely used to compute the physical properties of correlated quantum materials, their applicability may get severely hindered precisely in the parameter regions where the most exciting physics is observed. One of the major issues, referred to as "misleading convergence", is the tendency of iterative schemes to converge to unphysical fixed points for intermediate-to-strong electronic interactions, regardless of the numerical accuracy of the computation. Here, we explicitly verify that the origin of this problem in several established self-consistent many-electron approaches, defined in the general diagrammatic framework of the boson-exchange formalism, resides exclusively in the stability condition of the respective iteration schemes, and not in an intrinsic breakdown of their self-consistent diagrammatic description. This insight enables a simple and general remedy, as recently proposed in Phys. Rev. Lett. 137, 016502 (2026): The redefinition of the iterative procedure, by inverting the unstable eigendirections of the Jacobian associated with the fixed point of the self-consistent algorithm. We illustrate the successful outcome of this procedure by means of systematic calculations performed on exactly solvable testbed models. Our results demonstrate that the physical fixed point of the diagrammatic schemes we considered can be stabilized, de facto, across the entire parameter range, including the most challenging nonperturbative/strong-coupling regimes.
While self-consistent diagrammatic approaches are widely used to compute the physical properties of correlated quantum materials, their applicability may get severely hindered precisely in the parameter regions where the most exciting physics is observed. One of the major issues, referred to as "misleading convergence", is the tendency of iterative schemes to converge to unphysical fixed points for intermediate-to-strong electronic interactions, regardless of the numerical accuracy of the computation. Here, we explicitly verify that the origin of this problem in several established self-consistent many-electron approaches, defined in the general diagrammatic framework of the boson-exchange formalism, resides exclusively in the stability condition of the respective iteration schemes, and not in an intrinsic breakdown of their self-consistent diagrammatic description. This insight enables a simple and general remedy, as recently proposed in Phys. Rev. Lett. 137, 016502 (2026): The redefinition of the iterative procedure, by inverting the unstable eigendirections of the Jacobian associated with the fixed point of the self-consistent algorithm. We illustrate the successful outcome of this procedure by means of systematic calculations performed on exactly solvable testbed models. Our results demonstrate that the physical fixed point of the diagrammatic schemes we considered can be stabilized, de facto, across the entire parameter range, including the most challenging nonperturbative/strong-coupling regimes.
Self-consistent approaches in many-electron problems typically converge to an unphysical solution in strongly correlated regimes. By deriving the mathematical condition for the stability of the physical solution, we unveil the precise relation between two distinct issues previously considered equivalent: the misleading convergence in self-consistent schemes and the multivaluedness of the Luttinger-Ward functional. Although these problems are fundamentally linked through the divergences of the irreducible vertex function, we show that misleading convergence can occur even in the absence of such divergences. Eventually, a systematic procedure for stabilizing the physical solution is proposed.
Accurately describing many-body effects in multi-orbital systems remains a major challenge in theoretical condensed matter physics. At present, there is a significant methodological gap between the numerical tools used in ab initio computational materials science and those developed to study strong electronic correlations. The former can treat realistic, large-scale systems but typically neglect many-body effects, while the latter focus on simplified models with only a few degrees of freedom, as only such models can be solved accurately in the presence of strong interactions. The purpose of this thesis is to bridge these two approaches and establish a systematic theoretical framework for realistic correlated electronic materials. This involves a full-cycle methodology that begins with constructing ab initio interacting models from density-functional theory, solving them using dynamical mean-field theory to capture local correlations, and extending beyond to incorporate non-local collective electronic fluctuations. To this end, we introduce the "dual" approach to strong correlations, which includes the dual fermion, dual boson, and dual triply irreducible local expansion methods. The central idea of the dual theories is to shift the reference point of the conventional Feynman diagrammatic expansion from a non-interacting electronic system to an interacting but exactly solvable one. Integrating out this reference system recasts the expansion in an effective dual space, where all diagrammatic building blocks are renormalized by the corresponding impurity quantities. This procedure transforms a non-perturbative expansion in the original variables into a perturbative expansion in terms of dual fermionic and bosonic fields, exact in both the weak- and strong-coupling limits. In this thesis, we collect and systematize the major developments of dual techniques achieved to date.
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.
Strongly correlated quantum field theories present a formidable challenge to conventional theoretical methods, particularly those rooted in perturbative expansions. The breakdown of perturbation theory in regimes where interaction strengths are comparable to or exceed kinetic energies necessitates the development of non-perturbative approaches. This paper introduces and elaborates on a novel framework for diagrammatic renormalization specifically designed for strongly correlated quantum fields, moving beyond the limitations inherent in iterative or weak-coupling expansions. We explore how a systematic resummation of infinite classes of diagrams, guided by a functional renormalization group perspective, can capture the complex interplay of quantum fluctuations and emergent phenomena characteristic of these systems. Our approach leverages a self-consistent treatment of Green's functions and self-energies, allowing for the faithful description of spectral functions, phase transitions, and exotic ground states in a manner that traditional perturbation theory cannot. We discuss the theoretical underpinnings, computational strategies, and potential applications across diverse domains, including condensed matter physics, quantum chromodynamics, and ultracold atomic systems. The methodology aims to provide a robust, diagrammatically transparent pathway to understanding the non-trivial fixed points and scaling behaviors that govern strongly correlated phenomena, offering a powerful alternative where standard approximations fail.
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
Quantum impurity models, which describe the coupling between interacting orbitals and a non-interacting bath, play a central role in the physics of strongly correlated electron systems. Solving a quantum impurity model in general requires the use of non-perturbative numerical methods. Hamiltonian-based approaches, which rely on an explicit bath discretization, are typically limited to a small number of bath sites or small entanglement, and diagrammatic methods suffer from sign problems, slow convergence, or diagram truncation approximations. Here we show that these two classes of methods can be combined: augmenting diagrammatic methods with a small auxiliary bath can reduce the residual problem to a regime where low-order perturbation theory is highly accurate and rapidly converging. In a simple benchmark, the precision of the hybrid approach surpasses bold-line calculations by several orders of magnitude; for a strongly interacting two-orbital model with a severe sign problem, convergence is achieved at three orders of magnitude lower computational cost than competing methods; and convergence to the unknown exact result is rapidly accelerated in a difficult realistic problem. Our results establish a practical route to high-precision quantum impurity solutions in correlated quantum systems.