The fundamentality of the exponential representation of a second-quantised correlated wave function is emphasised with an accent on the physical sense of cluster amplitudes as cumulants of the correlated ansatz. Three main wave function formalisms, namely, the configuration-interaction theory, the coupled-cluster approach, and the many-body perturbation theory (as well as their extensions, e.g. the equation-of-motion coupled-cluster method, multireference schemes, etc.), are represented in an exponential form, leading to a formulation of the working equations in terms of cluster amplitudes. By expressing the corresponding many-body tensor equations in terms of cluster amplitudes, we could unambiguously check connectivity types and the asymptotic behaviour of all tensors/scalars involved (in the formal limit of an infinite number of correlated particles). In particular, the appearance of disconnected cluster amplitudes corresponds to unphysical correlations. Besides, we demonstrate that the equation-of-motion coupled-cluster approach, as well as certain excited-state configuration-interaction methods, can be recast in a fully connected (exponential) form, thus breaking the common belief that all truncated configuration-interaction methods violate connectivity. Our work is based on the recently developed algebraic framework which can be viewed as a complement to the classical diagrammatic analysis.
The fundamentality of the exponential representation of a second-quantised correlated wave function is emphasised with an accent on the physical sense of cluster amplitudes as cumulants of the correlated ansatz. Three main wave function formalisms, namely, the configuration-interaction theory, the coupled-cluster approach, and the many-body perturbation theory (as well as their extensions, e.g. the equation-of-motion coupled-cluster method, multireference schemes, etc.), are represented in an exponential form, leading to a formulation of the working equations in terms of cluster amplitudes. By expressing the corresponding many-body tensor equations in terms of cluster amplitudes, we could unambiguously check connectivity types and the asymptotic behaviour of all tensors/scalars involved (in the formal limit of an infinite number of correlated particles). In particular, the appearance of disconnected cluster amplitudes corresponds to unphysical correlations. Besides, we demonstrate that the equation-of-motion coupled-cluster approach, as well as certain excited-state configuration-interaction methods, can be recast in a fully connected (exponential) form, thus breaking the common belief that all truncated configuration-interaction methods violate connectivity. Our work is based on the recently developed algebraic framework which can be viewed as a complement to the classical diagrammatic analysis.
Scope The book focuses primarily on many-body (or better, many-electron) methods for electron correlation. These include Rayleigh–Schrödinger perturbation theory (RSPT) , particularly in its diagrammatic representation (referred to as many-body perturbation theory , or MBPT ), and coupled-cluster (CC) theory; their relationship to configuration interaction (CI) is included. Further extensions address properties other than the energy, and also excited states and multireference CC and MBPT methods. The many-body algebraic and diagrammatic methods used in electronic structure theory have their origin in quantum field theory and in the study of nuclear matter and nuclear structure. The second-quantization formalism was first introduced in a treatment of quantized fields by Dirac (1927) and was extended to fermion systems by Jordan and Klein (1927) and by Jordan and Wigner (1928). This formalism is particularly useful in field theory, in scattering problems and in the study of infinite systems because it easily handles problems involving infinite, indefinite or variable numbers of particles. The diagrammatic approach was introduced into field theory by Feynman (1949a,b) and applied to many-body systems by Hugenholtz (1957) and by Goldstone (1957). Many-body perturbation theory and its linked-diagram formalism were first introduced by Brueckner and Levinson (1955) and by Brueckner (1955), and were formalized by Goldstone (1957). Other important contributions to the methodology, first in field theory and then in the theory of nuclear structure, are due to Dyson (1949a,b), Wick (1950), Hubbard (1957, 1958a,b) and Frantz and Mills (1960). Applications to the electronic structure of atoms and molecules began with the work of Kelly (1963, 1964a,b, 1968), and molecular applications using finite analytical basis sets appeared in the work of Bartlett and Silver (1974a, b).
The viability of treating higher excitations in coupled-cluster theory is discussed. An algorithm is presented for solving coupled-cluster (CC) equations which can handle any excitation. Our method combines the formalism of diagrammatic many-body perturbation theory and string-based configuration interaction (CI). CC equations are explicitly put down in terms of antisymmetrized diagrams and a general method is proposed for the factorization of the corresponding algebraic expressions. Contractions between cluster amplitudes and intermediates are evaluated by a string-based algorithm. In contrast to our previous developments [J. Chem. Phys. 113, 1359 (2000)] the operation count of this new method scales roughly as the (2n+2)nd power of the basis set size where n is the highest excitation in the cluster operator. As a by-product we get a completely new CI formalism which is effective for solving both truncated and full CI problems. Generalization for approximate CC models as well as multireference cases is also discussed.
Coupled cluster (CC) theory is the “gold standard” of molecular electronic structure theory, yet the full arc of the subject—from its variational origins in Hartree–Fock theory, through the exponential ansatz and its diagrammatic machinery, to its modern reformulation as a system of polynomial equations on an algebraic variety, its multireference generalizations, and its emerging reincarnation as a quantum circuit—is rarely presented as a single connected narrative. This review builds that narrative completely from first principles and carries it through to today’s research frontier. Starting from the electronic Schrödinger equation, we derive Hartree–Fock theory, second quantization, and Wick’s theorem; construct the coupled cluster singles-and-doubles (CCSD) equations via the similarity-transformed Hamiltonian; and show how these same equations reappear, verbatim, as a truncated eigenvalue problem on a projective variety generalizing the Grassmannian, with a computable CC degree controlling the number of solutions. We then develop the perturbative triples correction CCSD(T), the linear-response equationof-motion (EOM-CC) framework for excited, ionized, and electron-attached states, the effective-Hamiltonian formalism underlying multireference coupled cluster theory and its notorious intruder-state problem, and the unitary coupled cluster (UCC) ansatz that restores the variational principle at the cost of an infinite operator expansion—tamed, on a quantum computer, by Trotterization and the variational quantum eigensolver. We close with a synthesis of what is rigorously known about the resulting nonlinear equations, a map of directions deliberately left out of scope, and an assessment of the open question of quantum advantage in this domain. Every claim is stated formally, with a full statement and proof or proof sketch, and cited to its origin; the computational cost of every method is compared explicitly throughout.
Many-body perturbation theories (MBPT) (both single and multireference), by design, are capable of representing electron correlation very well. However, the practical implementation of these methods for many-electron systems beyond a certain order is extremely difficult, if not impossible. For instance, there are hardly any calculations beyond third order at the multireference and fifth order at the single reference of perturbation theories. Thus, these methods are typically limited to systems where the higher-order electron correlation effect is not significant. A practical solution to this problem is the coupled cluster theory, which was first introduced by Coester [20] and Coester and Kümmel [21] in nuclear physics and was later adopted in the electron structure theory by Čížek [23, 106] and Paldus and Čížek [107]. The coupled-cluster (CC) theory is a non-perturbative all-order theory. This stateof-the-art theory has now been recognized as the basis of modern quantitative quantum chemistry and atomic physics for high-precision calculations. Its superior accuracy is rooted in the exponential ansatz for the wave function, which provides an energy expression that scales correctly with the number of electrons and under certain provisions leads to a faster convergence of the many-body expansion.
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 wave function 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 preserves a sum-over-state representation of the self-energy and leads to potentially systematically improvable Green’s function approaches.
Abstract. Coupled cluster theory produced arguably the most widely used high-accuracy computational quantum chemistry methods. Despite the approach’s overall great computational success, its mathematical understanding is so far limited to results within the realm of functional analysis. The coupled cluster amplitudes, which are the targeted objects in coupled cluster theory, correspond to solutions to the coupled cluster equations, which is a system of polynomial equations of at most degree four. The high dimensionality of the electronic Schrödinger equation and the nonlinearity of the coupled cluster ansatz have so far stalled a formal analysis of this polynomial system. In this article, we present algebraic investigations that shed light on the coupled cluster equations and the root structure of this ansatz. This is of importance for the a posteriori evaluation of coupled cluster calculations. To that end, we investigate the root structure by means of Newton polytopes. We derive a general v-description, which is subsequently turned into an h-description for explicit examples. This perspective reveals an apparent connection between Pauli’s exclusion principle and the geometrical structure of the Newton polytopes. We also propose an alternative characterization of the coupled cluster equations projected onto singles and doubles as cubic polynomials on an algebraic variety with certain sparsity patterns. Moreover, we provide numerical simulations of two computationally tractable systems, namely, the two electrons in four spin-orbitals system and the three electrons in six spin-orbitals system. These simulations provide novel insight into the root structure of the coupled cluster solutions when the coupled cluster ansatz is truncated.