Connection between $GW$ and Extended Coupled Cluster

Coupled-cluster (CC) theory and Green’s function many-body perturbation theory (MBPT) have long evolved as distinct yet complementary frameworks for describing electronic correlation. While CC methods employ exponential wavefunction parametrizations that guarantee size extensivity and systematic improvability, Green’s function approaches such as the $GW$ approximation describe quasiparticle and optical excitations through diagrammatic resummations. Recent analyses have established a formal correspondence between these frameworks: the $GW$ approximation is equivalent to an equation-of-motion (EOM) treatment of the direct-ring coupled-cluster doubles (drCCD) method. Within this context, the extended CC (ECC) ansatz offers a unified framework connecting CC and MBPT. This formulation bridges CC-based and Green’s function-based methods, providing novel avenues for incorporating vertex corrections within a CC framework that keep a positive semi-definite self-energy and lead to potentially systematically improvable Green’s function approaches.

💡 Research Summary

This paper investigates the formal relationship between the GW approximation—a cornerstone of Green’s‑function many‑body perturbation theory—and the extended coupled‑cluster (ECC) formalism, a bi‑variational generalization of traditional coupled‑cluster (CC) theory. The authors first review the historical development of CC and GW, emphasizing that CC provides size‑extensive, systematically improvable wave‑function based energies and excited‑state properties via equation‑of‑motion (EOM) or linear‑response approaches, while GW, derived from Hedin’s equations, yields quasiparticle energies and satellite structures through a dynamically screened Coulomb interaction. Recent work has shown that GW is mathematically equivalent to an EOM treatment of the direct‑ring CCD (drCCD) method, i.e., GW contains only the ring‑diagram subset of many‑body perturbation theory. However, drCCD lacks exchange, ladder, and mixed ring‑ladder diagrams—collectively referred to as vertex corrections—so GW cannot capture all correlation effects present in higher‑order CC models.

To bridge this gap, the authors introduce the ECC ansatz, which optimizes both excitation (T) and de‑excitation (Λ) operators simultaneously, leading to a double similarity‑transformed Hamiltonian. This bi‑variational framework yields a true energy functional, guarantees the Hellmann‑Feynman theorem, and naturally incorporates third‑order terms absent in conventional CCD, thereby providing a pathway to embed vertex corrections within a CC‑compatible formalism.

The paper then recasts GW as an electron‑boson coupling problem. In this picture, the bosonic degrees of freedom correspond to the collective electronic response described by the direct random‑phase approximation (RPA). The electron‑boson Hamiltonian consists of a fermionic part (the usual Fock operator), a bosonic part (quadratic in quasiboson operators built from particle‑hole excitations), and a linear electron‑boson interaction term involving effective two‑electron integrals. By performing a Bogoliubov transformation that diagonalizes the RPA bosonic Hamiltonian, the authors obtain a set of ideal boson modes with energies Ωμ.

Applying the ECC double similarity transformation to this electron‑boson Hamiltonian yields an “ECC electron‑boson Hamiltonian” whose matrix representation is exactly the GW supermatrix previously introduced in the literature. Consequently, the eigenvalue problem for charged excitations (ionization potentials and electron affinities) derived from the ECC Hamiltonian is mathematically identical to the GW quasiparticle equation, establishing a rigorous equivalence between the two formalisms.

Having established equivalence, the authors identify the additional terms generated by the ECC transformation that correspond to vertex corrections beyond the standard GW approximation. These terms arise from higher‑order commutators involving both excitation and de‑excitation operators and can be interpreted as three‑body and four‑body contributions to the self‑energy. By linearizing the one‑body density matrix obtained from the double‑similarity‑transformed Hamiltonian, they derive a low‑order correction to the Fock matrix. This correction mimics the effect of a self‑consistent GW (GW‑SC) calculation while preserving the positive semi‑definite nature of the self‑energy.

Numerical tests are performed on a benchmark set of small to medium‑size molecules. The authors compute ionization potentials (IPs) and electron affinities (EAs) using (i) conventional GW, (ii) ECC‑GW (i.e., GW plus the identified vertex corrections), and (iii) the linearized GW‑SC approximation. Results show that ECC‑GW systematically reduces the mean absolute error relative to reference coupled‑cluster singles‑doubles with perturbative triples