Pair theory revisited ii : time-independent diagrammatic perturbation approach

Josef Paldus, Jakub Čı́žek

Michael A. Robb

M.K. Pal

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.

Vincent Hurtubise, Karl F. Freed

Effective Hamiltonians and effective operators produce, respectively, exact energies and matrix elements of a time-independent operator A for a finite number of eigenstates of a time-independent Hamiltonian H. We obtain degenerate and quasidegenerate perturbative expressions for the particularly useful canonical effective operator ÂC through second order in perturbation theory. The corresponding ÂC diagrammatic expressions are derived for the case where ÂC acts in a complete finite space. Our first order results have been used previously for ab initio computations of dipole and transition dipole moments in diatomic hydrides and for testing the assumptions in semiempirical methods for dipole properties. A brief discussion is also provided on the computational labors required by first and second order ÂC many-body calculations, the derivation of ÂC diagrams when ÂC acts in an incomplete finite space, and on the derivation of diagrammatic rules for ÂC in arbitrary perturbation order.

Yuqi Wang, Wei‐Hai Fang, Zhendong Li

Many-body perturbation theory (MBPT) based on Green’s functions and Feynman diagrams provides a fundamental theoretical framework for various ab initio computational approaches in molecular and materials science, including the random phase approximation (RPA) and GW approximation. Unfortunately, this perturbation expansion often fails in systems with strong multireference characters. Extending diagrammatic MBPT to the multireference case is highly nontrivial and remains largely unexplored, primarily due to the breakdown of Wick’s theorem. In this work, we develop a diagrammatic multireference generalization of MBPT for computing correlation energies of strongly correlated systems, by using the cumulant expansion of many-body Green’s functions in place of Wick’s theorem. This theoretical framework bridges the gap between MBPT in condensed matter physics and multireference perturbation theories (MRPT) in quantum chemistry, which had been almost exclusively formulated within time-independent wave function frameworks prior to this work. Our formulation enables the explicit incorporation of strong correlation effects from the outset as in MRPT, while treating residual weak interactions through a generalized diagrammatic perturbation expansion as in MBPT. As a concrete demonstration, we formulate a multireference (MR) extension of the standard single-reference (SR) RPA by systematically resumming generalized ring diagrams, which naturally leads to a unified set of equations applicable to both SR and MR cases. Benchmark calculations on prototypical molecular systems reveal that MR-RPA successfully resolves the well-known failure of SR-RPA in strongly correlated systems. This theoretical advancement paves the way for advancing ab initio computational methods through diagrammatic resummation techniques in future.