On the basis of an exact perturbational expression for the interacting one-particle Green function
G corresponding to bosons / fermions in terms of the bare interaction potential
v and permanents / determinants of the non-interacting one-particle Green function
G0, we deduce four recursive perturbation expansions for the self-energy
Σ. With
W denoting the dynamic screened interaction potential, these perturbation expansions are identical to those of
Σ in terms of (i) self-energy diagrams and
(v,G0), (ii)
G-skeleton self-energy diagrams and
(v,G), (iii)
W-skeleton self-energy diagrams and
(W,G0), and (iv)
G- and
W-skeleton self-energy diagrams and
(W,G). For the calculation of
W, we rely on a similar exact perturbational expression for the interacting two-particle Green function
G2 as for
G. From this expression, we deduce four recursive perturbation expansions for the polarization function
P, necessary for the calculation of
W, that are similar to those for
Σ specified above. The correlation functions considered in this paper may be corresponding to ground states, and thermal ensemble of states. For thermal ensemble of states, we consider both the imaginary-time formalism of Matsubara, and the real-time formalism of thermo-field dynamics (TFD). The latter is advantageous for the direct calculation of dynamic correlation functions. In an appendix, we apply the formalisms presented in this paper to the Hubbard Hamiltonian for spin-
21 fermions on a lattice in arbitrary
d spatial dimensions. In two further appendices, we present methods and short programs for determining the
νth-order diagrams corresponding to the perturbation expansions of
G in terms of
(v,G0), and
Σ in terms of
(v,G0) and
(v,G) on the basis of the cycle decompositions of the elements of the symmetric group
S2ν.
We extend the self-consistent Green's functions formalism to take into account three-body interactions. We analyze the perturbative expansion in terms of Feynman diagrams and define effective one- and two-body interactions, which allows for a substantial reduction of the number of diagrams. The procedure can be taken as a generalization of the normal ordering of the Hamiltonian to fully correlated density matrices. We give examples up to third order in perturbation theory. To define nonperturbative approximations, we extend the equation-of-motion method in the presence of three-body interactions. We propose schemes that can provide nonperturbative resummation of three-body interactions. We also discuss two different extensions of the Koltun sum rule to compute the ground state of a many-body system.
In order to investigate the validity of the Brueckner and the simple sum of the ladder diagrams approximations to the energy of an infinite system of fermions, we have calculated as a function of density all the terms in the perturbation expansion through the fourth order in the strength of the potential. We have done this for a repulsive, square-well, two-body potential spin-\textonehalf{} fermions. We are able to construct rigorous upper and lower bounds from the coefficients for the value of the ladder approximation and find that the standard solution procedures give reasonably accurate results (within a few percent). We find that the error in the solution as obtained in practice to the equations of the Brueckner approximation is large compared to the size of its departure from the ladder approximation. We further find for low-to-moderate densities and for low-to-moderate potential strengths that the Brueckner approximation both as a sum of a certain class of diagrams and as computed in practice lies above the ladder approximation while the complete perturbation theory lies below it. This result arises from the neglect of the ring diagrams by the Brueckner approximation.
A thorough analytical and numerical characterization of the whole perturbation series of one-particle many-body Green's function (MBGF) theory is presented in a pedagogical manner. Three distinct but equivalent algebraic (first-quantized) recursive definitions of the perturbation series of the Green's function are derived, which can be combined with the well-known recursion for the self-energy. Six general-order algorithms of MBGF are developed, each implementing one of the three recursions, the ΔMPn method (where n is the perturbation order) [S. Hirata et al., J. Chem. Theory Comput. 11, 1595 (2015)], the automatic generation and interpretation of diagrams, or the numerical differentiation of the exact Green's function with a perturbation-scaled Hamiltonian. They all display the identical, nondivergent perturbation series except ΔMPn, which agrees with MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 but converges at the full-configuration-interaction (FCI) limit at n=∞ (unless it diverges). Numerical data of the perturbation series are presented for Koopmans and non-Koopmans states to quantify the rate of convergence towards the FCI limit and the impact of the diagonal, frequency-independent, or ΔMPn approximation. The diagrammatic linkedness and thus size-consistency of the one-particle Green's function and self-energy are demonstrated at any perturbation order on the basis of the algebraic recursions in an entirely time-independent (frequency-domain) framework. The trimming of external lines in a one-particle Green's function to expose a self-energy diagram and the removal of reducible diagrams are also justified mathematically using the factorization theorem of Frantz and Mills. Equivalence of ΔMPn and MBGF in the diagonal and frequency-independent approximations at 1≤n≤3 is algebraically proven, also ascribing the differences at n = 4 to the so-called semi-reducible and linked-disconnected diagrams.
We develop a perturbation expansion in the intersite hopping around the atomic limit of the Hubbard model. It is valid for arbitrary finite temperatures and interaction strengths. Diagrammatic rules that determine the grand-canonical potential and the Green's functions are derived. They reduce the calculation of any finite-order contribution to simple algebra. This opens the way for series extrapolations from computer-aided high-finite-order evaluations. Discrepancies in earlier expansions around the atomic limit are clarified. The present expansion scheme involves only connected diagrams with unrestricted lattice sums. This allows one to perform a vertex renormalization as for the linked-cluster expansion of the Ising model. The renormalized perturbation expansion can be used to construct self-consistent approximations which are automatically exact in the atomic limit. In the limit of high lattice dimensions, only fully two-particle reducible embeddings of diagrams on the lattice contribute. The single-particle properties of the infinite-dimensional Hubbard model reduce to those of independent tight-binding fermions hopping between dressed sites.
Abstract Starting with the equation of motion for the field operator ψ(x,t) of an interacting many-particle system, the n-particle Green’s function (Gn) equation of motion is developed, with interparticle interactions generating an infinite chain of equations coupling it to (n+1)- and (n−1)-particle Green’s functions (Gn+1 and Gn−1, respectively). Particularly important are the one-particle Green’s function equation with its coupling to the two-particle Green’s function and the two-particle Green’s function equation with its coupling to the three-particle Green’s function. To develop solutions, it is necessary to introduce non-correlation decoupling procedures involving the Hartree and Hartree-Fock approximations for G2 in the G1 equation; and a similar factorization “ansatz” for G3 in the G2 equation, resulting in the Sum of Ladder Diagrams integral equation for G2, with multiple Born iterates and finite collisional lifetimes. Similar treatment of the G11-equation for the joint propagation of one-electron and one-hole subject to mutual Coulomb attraction leads to bound electron-hole exciton states having a discrete hydrogen like spectrum of energy eigenstates. Its role in single-particle propagation is also discussed in terms of one-electron self-energy Σ and the T-matrix