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:2504.04459 · The Journal of Physical Chemistry A

Excitation Energies from the Entanglement Coupled Cluster Model for Doublets

Sarai Dery Folkestad, Kristine Lauvstad Kruken, Henrik Koch A

Published Jul 21, 2025• 1 citation
arXivHTMLPDF

Abstract

We present a spin-adapted coupled cluster singles and doubles model for the excitation energies of molecular doublets. The entanglement coupled cluster approach represents an unconventional take on the notorious problem of spin adaptation for open-shell species. In this approach, the high-spin open-shell molecular system is coupled to non-interacting bath orbitals to form a total closed-shell system. In entanglement coupled cluster theory, many of the attractive features of the spin-adapted closed-shell coupled cluster are retained: an unambiguous definition of the cluster operator and a terminating Baker-Campbell-Hausdorff expansion. The model produces excitation energies of a quality comparable to that of the closed-shell counterpart. Additionally, some ionized states that cannot be modeled accurately with the alternative equation-of-motion approach for ionized states can be described with the entanglement coupled cluster singles and doubles model.

DOI
10.1021/acs.jpca.5c01190

Related papers

Entanglement coupled cluster theory: Exact spin-adaptation

Sarai Dery Folkestad, Bendik Støa Sannes, Henrik Koch A

We present a novel framework for spin-adapted coupled cluster theory. The approach exploits the entanglement of an open-shell molecule with electrons in a non-interacting bath. Together, the molecule and the bath form a closed-shell system, and electron correlation can be included using the standard spin-adapted closed-shell coupled cluster formalism. A projection operator, which enforces conditions on the electrons in the bath, is used to obtain the desired state of the molecule. This entanglement coupled cluster theory is outlined, and proof-of-concept calculations for doublet states are reported. The approach is further extendable to open-shell systems with other values of the total spin.

arXiv:2301.12228•The Journal of Chemical Physics•Jun 12, 2023•6 citations

Coupled cluster theory for multireference and open-shell molecules

Sarai Dery Folkestad, Henrik Koch

We introduce a simple coupled cluster approach for multireference and open-shell systems. In the restricted active space coupled cluster framework, the cluster operator only includes excitations from doubly occupied to virtual orbitals, and this ensures the individual operators commute among each other. The ground and excited states are obtained by diagonalisation of the coupled cluster similarity transformed Hamiltonian in a restricted active space. We demonstrate the framework for the RAS-CCSD model, where single and double excitations are included in the cluster operator, and likewise for a complete active space. We consider systems with multireference character and open-shell systems in a spin-adapted manner. The approach captures significantly more dynamical correlation than without the coupled cluster similarity transformation.

physics.chem-phCross-listedarXiv:2610.07440•Oct 7, 2026

Formulation and implementation of a unitary group adapted state universal multi-reference coupled cluster (UGA-SUMRCC) theory: Excited and ionized state energies

Sangita Sen, Avijit Shee, Debashis Mukherjee

The traditional state universal multi-reference coupled cluster (SUMRCC) theory uses the Jeziorski-Monkhorst (JM) based Ansatz of the wave operator: Ω = Σ(μ)Ω(μ)|φ(μ)><φ(μ)|, where Ω(μ) = exp(T(μ)) is the cluster representation of the component of Ω inducing virtual excitations from the model function φ(μ). In the first formulations, φ(μ)s were chosen to be single determinants and T(μ)s were defined in terms of spinorbitals. This leads to spin-contamination for the non-singlet cases. In this paper, we propose and implement an explicitly spin-free realization of the SUMRCC theory. This method uses spin-free unitary generators in defining the cluster operators, {T(μ)}, which even at singles-doubles truncation, generates non-commuting cluster operators. We propose the use of normal-ordered exponential parameterization for Ω:Σ(μ){exp(T(μ))}|φ(μ)><φ(μ)|, where {} denotes the normal ordering with respect to a common closed shell vacuum which makes the "direct term" of the SUMRCC equations terminate at the quartic power. We choose our model functions {φ(μ)} as unitary group adapted (UGA) Gel'fand states which is why we call our theory UGA-SUMRCC. In the spirit of the original SUMRCC, we choose exactly the right number of linearly independent cluster operators in {T(μ)} such that no redundancies in the virtual functions {χ(μ)(l)} are involved. Using example applications for electron detached/attached and h-p excited states relative to a closed shell ground state we discuss how to choose the most compact and non-redundant cluster operators. Although there exists a more elaborate spin-adapted JM-like ansatz of Datta and Mukherjee (known as combinatoric open-shell CC (COS-CC), its working equations are more complex. Results are compared with those from COS-CC, equation of motion coupled cluster methods, restricted open-shell Hartree-Fock coupled cluster, and full configuration interaction. We observe that our results are more accurate with respect to most other theories as a result of the use of the cluster expansion structure for our wave operator. Our results are comparable to those from the more involved COS-CC, indicating that our theory captures the most important aspects of physics with a considerably simpler scheme.

The Journal of Chemical Physics•Aug 15, 2012•26 citations

An Alternative Algebraic Framework for the Simplification of Coupled Cluster Type Expectation Values

Daniel Pape

The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.

•Nov 12, 2013•0 citations

Correlation problems in atomic and molecular systems. V. Spin-adapted coupled cluster many-electron theory

Josef Paldus

The graphical techniques of spin algebras are combined with a diagrammatic approach based on the time independent Wick theorem to yield the spin-adapted form of the coupled cluster theory. The general rules for the implementation of this formalism are formulated and illustrated on the basic coupled-pair many-electron theory pertaining to closed shell ground states, for which case the explicit spin-adapted equations are derived. The advantages of the spin-adapted form of the theory are discussed along with the new insights it affords.

The Journal of Chemical Physics•Jul 1, 1977•242 citations

Unitary group approach to spin‐adapted open‐shell coupled cluster theory

Bogumił Jeziorski, Josef Paldus, Piotr Jankowski

Abstract We show that the irreducible tensor operators of the unitary group provide a natural operator basis for the exponential Ansatz which preserves the spin symmetry of the reference state, requires a minimal number of independent cluster amplitudes for each substitution order, and guarantees the invariance of the correlation energy under unitary transformations of core, open‐shell, and virtual orbitals. When acting on the closed‐shell reference state with nc doubly occupied and nv unoccupied (virtual) orbitals, the irreducible tensor operators of the group U(nc) ⊗ U(nV) generate all Gelfand‐Tsetlin (GT) states corresponding to appropriate irreducible representation of U(nc + nv). The tensor operators generating the M‐tuply excited states are easily constructed by symmetrizing products of M unitary group generators with the Wigner operators of the symmetric group SM. This provides an alternative to the Nagel‐Moshinsky construction of the GT basis. Since the corresponding cluster amplitudes, which are also U(nc) ⊗ U(ns) tensors, can be shown to be connected, the irreducible tensor operators of U(nc) ⊗ U(nv) represent a convenient basis for a spin‐adapted full coupled cluster calculation for closed‐shell systems. For a high‐spin reference determinant with n, singly occupied open‐shell orbitals, the corresponding representation of U(n), n=nc + nv + ns is not simply reducible under the group U(nc) ⊗ U(ns) ⊗ U(nv). The multiplicity problem is resolved using the group chain U(n) ⊃ U(nc + nv) ⊗ U(ns) ⊃ U(nc) ⊗U(ns)⊗ U(nv) ⊗ U(nv). The labeling of the resulting configuration‐state functions (which, in general, are not GT states when nc > 1) by the irreducible representations of the intermediate group U(nc + nv) ⊗U(ns) turns out to be equivalent to the classification based on the order of interaction with the reference state. The irreducible tensor operators defined by the above chain and corresponding to single, double, and triple substitutions from the first‐, second‐, and third‐order interacting spaces are explicitly constructed from the U(n) generators. The connectedness of the corresponding cluster amplitudes and, consequently, the size extensivity of the resulting spin‐adapted open‐shell coupled cluster theory are proved using group theoretical arguments. The perturbation expansion of the resulting coupled cluster equations leads to an explicitly connected form of the spin‐restricted open‐shell many‐body perturbation theory. Approximation schemes leading to manageable computational procedures are proposed and their relation to perturbation theory is discussed. © 1995 John Wiley & Sons, Inc.

International Journal of Quantum Chemistry•Nov 5, 1995•79 citations
We present a spin-adapted coupled cluster singles and doubles model for the excitation energies of molecular doublets. The entanglement coupled cluster approach represents an unconventional take on the notorious problem of spin adaptation for open-shell species. In this approach, the high-spin open-shell molecular system is coupled to non-interacting bath orbitals to form a total closed-shell system. In entanglement coupled cluster theory, many of the attractive features of the spin-adapted closed-shell coupled cluster are retained: an unambiguous definition of the cluster operator and a terminating Baker-Campbell-Hausdorff expansion. The model produces excitation energies of a quality comparable to that of the closed-shell counterpart. Additionally, some ionized states that cannot be modeled accurately with the alternative equation-of-motion approach for ionized states can be described with the entanglement coupled cluster singles and doubles model.
We present a novel framework for spin-adapted coupled cluster theory. The approach exploits the entanglement of an open-shell molecule with electrons in a non-interacting bath. Together, the molecule and the bath form a closed-shell system, and electron correlation can be included using the standard spin-adapted closed-shell coupled cluster formalism. A projection operator, which enforces conditions on the electrons in the bath, is used to obtain the desired state of the molecule. This entanglement coupled cluster theory is outlined, and proof-of-concept calculations for doublet states are reported. The approach is further extendable to open-shell systems with other values of the total spin.
We introduce a simple coupled cluster approach for multireference and open-shell systems. In the restricted active space coupled cluster framework, the cluster operator only includes excitations from doubly occupied to virtual orbitals, and this ensures the individual operators commute among each other. The ground and excited states are obtained by diagonalisation of the coupled cluster similarity transformed Hamiltonian in a restricted active space. We demonstrate the framework for the RAS-CCSD model, where single and double excitations are included in the cluster operator, and likewise for a complete active space. We consider systems with multireference character and open-shell systems in a spin-adapted manner. The approach captures significantly more dynamical correlation than without the coupled cluster similarity transformation.
The traditional state universal multi-reference coupled cluster (SUMRCC) theory uses the Jeziorski-Monkhorst (JM) based Ansatz of the wave operator: Ω = Σ(μ)Ω(μ)|φ(μ)><φ(μ)|, where Ω(μ) = exp(T(μ)) is the cluster representation of the component of Ω inducing virtual excitations from the model function φ(μ). In the first formulations, φ(μ)s were chosen to be single determinants and T(μ)s were defined in terms of spinorbitals. This leads to spin-contamination for the non-singlet cases. In this paper, we propose and implement an explicitly spin-free realization of the SUMRCC theory. This method uses spin-free unitary generators in defining the cluster operators, {T(μ)}, which even at singles-doubles truncation, generates non-commuting cluster operators. We propose the use of normal-ordered exponential parameterization for Ω:Σ(μ){exp(T(μ))}|φ(μ)><φ(μ)|, where {} denotes the normal ordering with respect to a common closed shell vacuum which makes the "direct term" of the SUMRCC equations terminate at the quartic power. We choose our model functions {φ(μ)} as unitary group adapted (UGA) Gel'fand states which is why we call our theory UGA-SUMRCC. In the spirit of the original SUMRCC, we choose exactly the right number of linearly independent cluster operators in {T(μ)} such that no redundancies in the virtual functions {χ(μ)(l)} are involved. Using example applications for electron detached/attached and h-p excited states relative to a closed shell ground state we discuss how to choose the most compact and non-redundant cluster operators. Although there exists a more elaborate spin-adapted JM-like ansatz of Datta and Mukherjee (known as combinatoric open-shell CC (COS-CC), its working equations are more complex. Results are compared with those from COS-CC, equation of motion coupled cluster methods, restricted open-shell Hartree-Fock coupled cluster, and full configuration interaction. We observe that our results are more accurate with respect to most other theories as a result of the use of the cluster expansion structure for our wave operator. Our results are comparable to those from the more involved COS-CC, indicating that our theory captures the most important aspects of physics with a considerably simpler scheme.
The Coupled-Cluster (CC) method is one of the most popular and efficient correlation methods in quantum chemistry. \nEspecially the CCSD(T) approximation, which includes single and double excitations by means of the application of the cluster operator \nto a reference determinant and triple excitations via a perturbative treatment, has become a standard tool \nin quantum chemical applications. However, the method is restricted to relatively small system sizes due to its unfavourable scaling (N^7, where N is the number of basis functions applied). \nFor the full treatment of triple excitations the scaling advances to N^8 and every further excitation level increases the exponent by two. \n \nThe goal of this work is to reduce the calculation time for closed shell systems at least by a factor growing with the excitation level for arbitrary truncation levels. \nThis is done by restricting the spatial parts of the spin orbitals and thus treat pairs of spin orbitals on the same footing. The restrictions can be easily constructed \nfor the CCSD model by spin integration. \nThe derivation of the restrictions arising in higher excited case will be done employing \nthe spatial orbital excitation operators E. \n \nIn the first part of this work an algorithm is derived that is capable of the derivation of the energy and amplitude equations for arbitrary excitation levels. In the second part \nan implementation of this algorithm is presented.
The graphical techniques of spin algebras are combined with a diagrammatic approach based on the time independent Wick theorem to yield the spin-adapted form of the coupled cluster theory. The general rules for the implementation of this formalism are formulated and illustrated on the basic coupled-pair many-electron theory pertaining to closed shell ground states, for which case the explicit spin-adapted equations are derived. The advantages of the spin-adapted form of the theory are discussed along with the new insights it affords.
Abstract We show that the irreducible tensor operators of the unitary group provide a natural operator basis for the exponential Ansatz which preserves the spin symmetry of the reference state, requires a minimal number of independent cluster amplitudes for each substitution order, and guarantees the invariance of the correlation energy under unitary transformations of core, open‐shell, and virtual orbitals. When acting on the closed‐shell reference state with nc doubly occupied and nv unoccupied (virtual) orbitals, the irreducible tensor operators of the group U(nc) ⊗ U(nV) generate all Gelfand‐Tsetlin (GT) states corresponding to appropriate irreducible representation of U(nc + nv). The tensor operators generating the M‐tuply excited states are easily constructed by symmetrizing products of M unitary group generators with the Wigner operators of the symmetric group SM. This provides an alternative to the Nagel‐Moshinsky construction of the GT basis. Since the corresponding cluster amplitudes, which are also U(nc) ⊗ U(ns) tensors, can be shown to be connected, the irreducible tensor operators of U(nc) ⊗ U(nv) represent a convenient basis for a spin‐adapted full coupled cluster calculation for closed‐shell systems. For a high‐spin reference determinant with n, singly occupied open‐shell orbitals, the corresponding representation of U(n), n=nc + nv + ns is not simply reducible under the group U(nc) ⊗ U(ns) ⊗ U(nv). The multiplicity problem is resolved using the group chain U(n) ⊃ U(nc + nv) ⊗ U(ns) ⊃ U(nc) ⊗U(ns)⊗ U(nv) ⊗ U(nv). The labeling of the resulting configuration‐state functions (which, in general, are not GT states when nc > 1) by the irreducible representations of the intermediate group U(nc + nv) ⊗U(ns) turns out to be equivalent to the classification based on the order of interaction with the reference state. The irreducible tensor operators defined by the above chain and corresponding to single, double, and triple substitutions from the first‐, second‐, and third‐order interacting spaces are explicitly constructed from the U(n) generators. The connectedness of the corresponding cluster amplitudes and, consequently, the size extensivity of the resulting spin‐adapted open‐shell coupled cluster theory are proved using group theoretical arguments. The perturbation expansion of the resulting coupled cluster equations leads to an explicitly connected form of the spin‐restricted open‐shell many‐body perturbation theory. Approximation schemes leading to manageable computational procedures are proposed and their relation to perturbation theory is discussed. © 1995 John Wiley & Sons, Inc.