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.
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.
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.
Abstract We develop algebraic geometry for coupled cluster (CC) theory of quantum many-body systems. The high-dimensional eigenvalue problems that encode the electronic Schrödinger equation are approximated by a hierarchy of polynomial systems at various levels of truncation. The exponential parametrization of the eigenstates gives rise to truncation varieties. These generalize Grassmannians in their Plücker embedding. We explain how to derive Hamiltonians, we offer a detailed study of truncation varieties and their CC degrees, and we present the state of the art in solving the CC equations.
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.
We develop and numerically analyze an algebraic-geometric framework for spin-adapted coupled-cluster (CC) theory. Since the electronic Hamiltonian is SU(2)-invariant, physically relevant quantum states lie in the spin singlet sector. We give an explicit description of the SU(2)-invariant (spin singlet) many-body space by identifying it with an Artinian commutative ring, called the excitation ring, whose dimension is governed by a Narayana number. We define spin-adapted truncation varieties via embeddings of graded subspaces of this ring, and we identify the CCS truncation variety with the Veronese square of the Grassmannian. Compared to the spin-generalized formulation, this approach yields a substantial reduction in dimension and degree, with direct computational consequences. In particular, the CC degree of the truncation variety -- governing the number of homotopy paths required to compute all CC solutions -- is reduced by orders of magnitude. We present scaling studies demonstrating asymptotic improvements and we exploit this reduction to compute the full solution landscape of spin-adapted CC equations for water and lithium hydride.
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.
Abstract This article presents an educational overview of the latest mathematical developments in coupled cluster (CC) theory, beginning with Schneider's seminal work from 2009 that introduced the first local analysis of CC theory. We provide a tutorial review of second quantization and the CC ansatz, laying the groundwork for understanding the mathematical basis of the theory. This is followed by a detailed exploration of the most recent mathematical advancements in CC theory. Our review starts with an in‐depth look at the local analysis pioneered by Schneider which has since been applied to various CC methods. We subsequently discuss the graph‐based framework for CC methods developed by Csirik and Laestadius. This framework provides a comprehensive platform for comparing different CC methods, including multireference approaches. Next, we delve into the latest numerical analysis results analyzing the single reference CC method developed by Hassan, Maday, and Wang. This very general approach is based on the invertibility of the CC function's Fréchet derivative. We conclude the article with a discussion on the recent incorporation of algebraic geometry into CC theory, highlighting how this novel and fundamentally different mathematical perspective has furthered our understanding and provides exciting pathways to new computational approaches.