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.
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.
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).
In this book, we encounter the hierarchy problem when we apply the equation-of-motion technique to analyze the quantum dynamics of the coupled light–matter system. To obtain systematic approximations, we use the so-called cluster-expansion method where many-body quantities are systematically grouped into cluster classes based on how important they are to the overall quantum dynamics. With increasing complexity, the clusters contain independent single particles (singlets), interacting pairs (doublets), three-particle terms (triplets), coupled four-particle contributions (quadruplets), as well as higher-order correlations. In this context, the N -particle concept is somewhat formal because it refers to generic N -particle expectation values 〈 N 〉, which may consist of an arbitrary mixture of carrier, photon, and phonon operators, as discussed in Section 14.2. In order to truncate the hierarchy problem at a given level, 〈 N 〉 is approximated through a functional structure that includes all clusters up to the predetermined level while all remaining clusters with a higher rank are omitted. It is natural that the corresponding approximations can be systematically improved by increasing the number of clusters included. To the best of our knowledge, the idea of coupled-clusters approaches was first formulated by Fritz Coester and Hermann Kümmel in the 1950s to describe nuclear many-body phenomena. The approach was then modified for the needs of quantum chemistry by Jiri Cizek 1966 to deal with many-body phenomena in atoms and molecules. Currently, it is one of the most accurate methods to compute molecular eigenstates.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
The study of systems with fractional charges and spins has become an extremely important tool to understand errors in approximate electronic structure methods, particularly in the context of density functional theory. Meanwhile, similar studies with wavefunction (WF)-based methods beyond second-order perturbation theory have been lacking. In this contribution, we study the performance of different coupled cluster (CC) and many-body perturbation theory (MBPT)-based methods for fractional charges. The use of the conventional and renormalized formulations of fractional-charge MBPT is discussed. The fractional spin behavior of the coupled cluster doubles (CCD) method is also investigated. Overall, all tested WF methods show very promising performance for the fractional charge problem. CCD is also quite accurate for the fractional spin problem in He+ across most of the range, although it breaks down to near Hartree-Fock quality in the strongly correlated limit. Beyond the study of fractional charge and spin curves, the implementation of CC methods with fractional occupation numbers offers a promising route to treating problems with multi-reference character in a single-reference framework.
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.
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.