Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2602.16580

On the Coupled Cluster Doubles Truncation Variety of Four Electrons

Fabian M. Faulstich, Vincenzo Galgano, Elke Neuhaus, Irem Portakal

math.AGphysics.chem-phquant-phSubmitted Feb 18, 2026• Announced Feb 19, 2026• 0 citations
arXivHTMLPDF

Abstract

We extend recent algebro-geometric results for coupled cluster theory of quantum many-body systems to the truncation varieties arising from the doubles approximation (CCD), focusing on the first genuinely nonlinear doubles regime of four electrons. Since this doubles truncation variety does not coincide with previously studied varieties, we initiate a systematic investigation of its basic algebro-geometric invariants. Combining theoretical and numerical results, we show that for $4$ electrons on $n\leq 12$ orbitals, the CCD truncation variety is a complete intersection of degree $2^{\binom{n-4}{4}}$. Using representation-theoretic arguments, we uncover a Pfaffian structure governing the quadratic relations that define the truncation variety for any $n$, and show that an exact tensor product factorization holds in a distinguished limit of disconnected doubles. We connect these structural results to the computation of the beryllium insertion into molecular hydrogen ({Be$\cdots$H$_2$ $\to$ H--Be--H}), a small but challenging bond formation process where multiconfigurational effects become pronounced.

Comments
26 pages, 6 figures, code stored in Zenodo. Comments are welcome!

Related papers

Coupled-cluster singles and doubles for extended systems

So Hirata, Rafał Podeszwa, Motoi Tobita +1 more

Coupled-cluster theory with connected single and double excitation operators (CCSD) and related approximations, such as linearized CCSD, quadratic configuration interaction with single and double excitation operators, coupled-cluster with connected double excitation operator (CCD), linearized CCD, approximate CCD, and second- and third-order many-body perturbation theories, are formulated and implemented for infinitely extended one-dimensional systems (polymers), on the basis of the periodic boundary conditions and distance-based screening of integrals, density matrix elements, and excitation amplitudes. The variation of correlation energies with the truncation radii of short- and long-range lattice sums and with the number of wave vector sampling points in the first Brillouin zone is examined for polyethylene, polyacetylene, and polyyne, and is shown to be a function of the degree of pi-electron conjugation or the fundamental band gaps. The t2 and t1 amplitudes in the atomic orbital (AO) basis are obtained by first computing the t amplitudes in the Bloch-orbital basis and subsequently back-transforming them into the AO basis. The plot of these AO-based t amplitudes as a function of unit cells also indicates that the t2 amplitudes of polyacetylene and polyyne exhibit appreciably slower decay than those of polyethylene, although the asymptotic decay behavior is invariably 1/r3. The AO-based t1 amplitudes appear to correlate strongly with the electronic structure, and they decay seemingly exponentially for polyethylene whereas they stay at a constant magnitude across the seventh nearest neighbors of polyacetylene and polyyne, which attests to far reaching effects of nondynamical electron correlation mediated by orbital rotation. Nonetheless, the unit cell contributions to the correlation energies taper below 10(-6) hartree after 15 A for all three polymers. The basis set dependence of the decay behavior of t2 amplitudes is also examined for linear hydrogen fluoride polymer (HF)infinity and linear beryllium polymer (Be)infinity employing the STO-3G, 6-31G, and 6-31G* basis sets, and proves to be rather small.

The Journal of Chemical Physics•Feb 4, 2004•162 citations

Coupled Cluster Theory: Toward an Algebraic Geometry Formulation

Fabian M. Faulstich, Mathias Oster

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.

arXiv:2211.10389•SIAM Journal on Applied Algebra and Geometry•Feb 22, 2024•12 citations

Algebraic Geometry for Spin-Adapted Coupled Cluster Theory

Fabian M. Faulstich, Svala Sverrisdóttir

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.

physics.chem-phCross-listedarXiv:2601.16646•Jan 26, 2026•0 citations

On Coupled Cluster Theory and Beyond

Mohsen Doust Mohammadi

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.

•Sep 8, 2026•0 citations

Tensor decomposition in post-Hartree–Fock methods. II. CCD implementation

Udo Benedikt, Karl-Heinz Böhm, Alexander A. Auer

In a previous publication, we have discussed the usage of tensor decomposition in the canonical polyadic (CP) tensor format for electronic structure methods. There, we focused on two-electron integrals and second order Møller-Plesset perturbation theory (MP2). In this work, we discuss the CP format for Coupled Cluster (CC) theory and present a pilot implementation for the Coupled Cluster Doubles method. We discuss the iterative solution of the CC amplitude equations using tensors in CP representation and present a tensor contraction scheme that minimizes the effort necessary for the rank reductions during the iterations. Furthermore, several details concerning the reduction of complexity of the algorithm, convergence of the CC iterations, truncation errors, and the choice of threshold for chemical accuracy are discussed.

The Journal of Chemical Physics•Dec 9, 2013•56 citations

Algebraic Varieties in Quantum Chemistry

Fabian M. Faulstich, Bernd Sturmfels, Svala Sverrisdóttir

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.

arXiv:2308.05258•Foundations of Computational Mathematics•Jul 17, 2024•5 citations
We extend recent algebro-geometric results for coupled cluster theory of quantum many-body systems to the truncation varieties arising from the doubles approximation (CCD), focusing on the first genuinely nonlinear doubles regime of four electrons. Since this doubles truncation variety does not coincide with previously studied varieties, we initiate a systematic investigation of its basic algebro-geometric invariants. Combining theoretical and numerical results, we show that for 444 electrons on n≤12n\leq 12n≤12 orbitals, the CCD truncation variety is a complete intersection of degree 2(n−44)2^{\binom{n-4}{4}}2(4n−4​). Using representation-theoretic arguments, we uncover a Pfaffian structure governing the quadratic relations that define the truncation variety for any nnn, and show that an exact tensor product factorization holds in a distinguished limit of disconnected doubles. We connect these structural results to the computation of the beryllium insertion into molecular hydrogen ({Be⋯\cdots⋯H2_22​ →\to→ H--Be--H}), a small but challenging bond formation process where multiconfigurational effects become pronounced.
Coupled-cluster theory with connected single and double excitation operators (CCSD) and related approximations, such as linearized CCSD, quadratic configuration interaction with single and double excitation operators, coupled-cluster with connected double excitation operator (CCD), linearized CCD, approximate CCD, and second- and third-order many-body perturbation theories, are formulated and implemented for infinitely extended one-dimensional systems (polymers), on the basis of the periodic boundary conditions and distance-based screening of integrals, density matrix elements, and excitation amplitudes. The variation of correlation energies with the truncation radii of short- and long-range lattice sums and with the number of wave vector sampling points in the first Brillouin zone is examined for polyethylene, polyacetylene, and polyyne, and is shown to be a function of the degree of pi-electron conjugation or the fundamental band gaps. The t2 and t1 amplitudes in the atomic orbital (AO) basis are obtained by first computing the t amplitudes in the Bloch-orbital basis and subsequently back-transforming them into the AO basis. The plot of these AO-based t amplitudes as a function of unit cells also indicates that the t2 amplitudes of polyacetylene and polyyne exhibit appreciably slower decay than those of polyethylene, although the asymptotic decay behavior is invariably 1/r3. The AO-based t1 amplitudes appear to correlate strongly with the electronic structure, and they decay seemingly exponentially for polyethylene whereas they stay at a constant magnitude across the seventh nearest neighbors of polyacetylene and polyyne, which attests to far reaching effects of nondynamical electron correlation mediated by orbital rotation. Nonetheless, the unit cell contributions to the correlation energies taper below 10(-6) hartree after 15 A for all three polymers. The basis set dependence of the decay behavior of t2 amplitudes is also examined for linear hydrogen fluoride polymer (HF)infinity and linear beryllium polymer (Be)infinity employing the STO-3G, 6-31G, and 6-31G* basis sets, and proves to be rather small.
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.
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.
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.
In a previous publication, we have discussed the usage of tensor decomposition in the canonical polyadic (CP) tensor format for electronic structure methods. There, we focused on two-electron integrals and second order Møller-Plesset perturbation theory (MP2). In this work, we discuss the CP format for Coupled Cluster (CC) theory and present a pilot implementation for the Coupled Cluster Doubles method. We discuss the iterative solution of the CC amplitude equations using tensors in CP representation and present a tensor contraction scheme that minimizes the effort necessary for the rank reductions during the iterations. Furthermore, several details concerning the reduction of complexity of the algorithm, convergence of the CC iterations, truncation errors, and the choice of threshold for chemical accuracy are discussed.
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.