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

Physical Review A

Coupled-cluster method in Fock space. III. On similarity transformation of operators in Fock space

Leszek Z. Stolarczyk, Hendrik J. Monkhorst

Published Mar 1, 1988• 55 citations
Publisher

Abstract

A detailed and systematic presentation of algebraic and diagrammatic techniques introduced in paper I [L. Z. Stolarczyk and H. J. Monkhorst, Phys. Rev. A 32, 725 (1985)] is given. Performing a similarity transformation of an operator in Fock space, e.g., a Hamiltonian in the generalized coupled-cluster method (see paper I), leads to lengthy algebraic formulas for the linear parameters (amplitudes) of the transformed operator. These formulas can be expressed in a compact form by means of the so-called reduced diagrams. Precise rules are given for constructing such diagrams (related to the so-called directed graphs of the graph theory) and the algebraic expressions corresponding to them. In particular, rules for generating a connected-diagram expansion of the amplitudes of the transformed operator are formulated. Several examples illustrating the use of this diagrammatic-algebraic approach are considered.

DOI
10.1103/physreva.37.1908

Related papers

Many-body similarity transformations generated by normal ordered exponential excitation operators

Marcel Nooijen

Normal ordered exponential operators have been used extensively in open-shell formulations of coupled cluster theory. The inverse of such an operator is known to exist, but a closed form explicit expression for the inverse is not available. We will address the evaluation of many-body similarity transformations generated by normal ordered exponential transformation operators without explicit use of the inverse. The similarity transform can be obtained as the solution of a linear system of equations that can be solved trivially using backward substitution. In addition a closed form diagrammatic expression for the similarity transformed operator is presented. Using the many-body similarity transformation strategy a simple and more general formulation of Fock space coupled cluster theory is presented which is akin in spirit to the formulation by Stolarczyk and Monkhorst [Phys. Rev. A 32, 725, 743 (1985); 37, 1908, 1926 (1988)], but which on the other hand is completely equivalent to the conventional wave operator formulation of Fock space coupled cluster theory (under suitable conditions). Other possible applications of the many-body similarity transformation will be briefly discussed.

The Journal of Chemical Physics•Feb 15, 1996•99 citations

Coupled-cluster method in Fock space. I. General formalism

Leszek Z. Stolarczyk, Hendrik J. Monkhorst

The problem of finding the spectrum of the Fock-space Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{} for a system of many fermions is analyzed. The quasiparticle formalism is employed, with ``holes'' and ``particles'' defined with respect to some N-particle determinantal wave function (the model vacuum). The basic idea behind the proposed approach is to perform a similarity transformation of Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{}, such that the resulting effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{} is, unlike \ifmmode \hat{H}\else \^{H}\fi{}, a quasiparticle-number--conserving operator. It is shown that eigenvalues of \ifmmode \hat{G}\else \^{G}\fi{}, corresponding to small numbers of quasiparticles (0,1,2) can be easily calculated. This is equivalent to finding eigenvalues of \ifmmode \hat{H}\else \^{H}\fi{} for certain states of N, N\ifmmode\pm\else\textpm\fi{}1, and N\ifmmode\pm\else\textpm\fi{}2 particles. The construction of the operator transforming \ifmmode \hat{H}\else \^{H}\fi{} into \ifmmode \hat{G}\else \^{G}\fi{} (the wave operator) stems from an analysis of the structure of the algebra of linear operators acting in a (finite-dimensional) Fock space. The exponential Ansatz for the wave operator is used, resulting in a generalization of the coupled-cluster (CC) method of Coester [Nucl. Phys. 7, 421 (1958)]. The generalized CC equations determining the wave operator, and equations determining the effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{}, are presented in a diagrammatic form. An effort has been made to obtain a concise notation for expressing these equations in an algebraic form. Approximation schemes, necessary for practical applications of the proposed method, are also studied.

Physical Review A•Aug 1, 1985•171 citations

Quantum chemistry in Fock space. II. Effective Hamiltonians in Fock space

Werner Kutzelnigg, Sigurd Koch

The concept of an effective Hamiltonian in Fock space is introduced. It is based on the division of the entire one-particle space into subspaces of ‘‘active’’ and ‘‘inactive’’ orbitals. The effective Fock space Hamiltonian has—for active model states—the same eigenvalues as the full Hamiltonian. The theory outlined in this context differs from that of paper I mainly in a different definition of the ‘‘diagonal part’’ of an operator, and in the fact that the ‘‘quasidegenerate case’’ applies throughout. The separation theorem, and as a consequence the connected diagram theorem, is shown to hold, in a more limited sense though, even for those normalizations where it did not in the context of universal wave and energy operators. Unlike in the theory of the ‘‘universal’’ operators of paper I the Fock space and n-particle Hilbert space approaches with analogous normalizations are no longer equivalent. In particular, the Primas normalization with a fully Lie-algebraic structure does not lead to a connected diagram expansion if it is formulated in n-particle Hilbert space, only so in a Fock space formulation. In n-particle Hilbert space with the present definition of the diagonal part of an operator the normalizations b (‘‘canonical’’) and c (‘‘Primas’’) happen to agree. As an alternative to the construction of the wave and energy operator W and L by perturbation theory the nonperturbative approach is presented as a generalization of the coupled-cluster method, in detail both in the intermediate and in the unitary normalization. In the unitary variant only a linear system for σ (the logarithm of the wave operator) has to be solved in order to get L correct through fifth order in perturbation theory with important contributions of higher orders included. A generalization of the Hartree–Fock method to Fock space theory is outlined, which guarantees stationarity of all (active) eigenstates with respect to one-particle transformations. A generalized electron pair theory is also defined. An analysis of the necessary computational steps shows that the nonperturbative approaches do not require significantly more computational effort than perturbation theory to the corresponding order. As a numerical example the H2 molecule in a small basis is discussed.

The Journal of Chemical Physics•Nov 1, 1983•171 citations

Green’s Functions in Coupled-Cluster Form

Hu Cheng-zheng

The Green’s function and coupled-cluster (CC) methods are two important tools of quantum many-body theory. It seems worthwhile to find their connection. In the present paper we demonstrate this connection and give the expressions of the Green’s functions in the CC form. Diagrammatic techniques are used to write explicit algebraic expressions.

Communications in Theoretical Physics•Apr 1, 1990•3 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

Diagrammatic Coupled Cluster Monte Carlo

Charles J. C. Scott, Roberto Di Remigio, T. Daniel Crawford +1 more

We propose a modified coupled cluster Monte Carlo algorithm that stochastically samples connected terms within the truncated Baker-Campbell-Hausdorff expansion of the similarity-transformed Hamiltonian by construction of coupled cluster diagrams on the fly. Our new approach-diagCCMC-allows propagation to be performed using only the connected components of the similarity-transformed Hamiltonian, greatly reducing the memory cost associated with the stochastic solution of the coupled cluster equations. We show that for perfectly local, noninteracting systems diagCCMC is able to represent the coupled cluster wavefunction with a memory cost that scales linearly with system size. The favorable memory cost is observed with the only assumption of fixed stochastic granularity and is valid for arbitrary levels of coupled cluster theory. Significant reduction in memory cost is also shown to smoothly appear with dissociation of a finite chain of helium atoms. This approach is also shown not to break down in the presence of strong correlation through the example of a stretched nitrogen molecule. Our novel methodology moves the theoretical basis of coupled cluster Monte Carlo closer to deterministic approaches.

arXiv:1901.02753•The Journal of Physical Chemistry Letters•Feb 6, 2019•19 citations
A detailed and systematic presentation of algebraic and diagrammatic techniques introduced in paper I [L. Z. Stolarczyk and H. J. Monkhorst, Phys. Rev. A 32, 725 (1985)] is given. Performing a similarity transformation of an operator in Fock space, e.g., a Hamiltonian in the generalized coupled-cluster method (see paper I), leads to lengthy algebraic formulas for the linear parameters (amplitudes) of the transformed operator. These formulas can be expressed in a compact form by means of the so-called reduced diagrams. Precise rules are given for constructing such diagrams (related to the so-called directed graphs of the graph theory) and the algebraic expressions corresponding to them. In particular, rules for generating a connected-diagram expansion of the amplitudes of the transformed operator are formulated. Several examples illustrating the use of this diagrammatic-algebraic approach are considered.
Normal ordered exponential operators have been used extensively in open-shell formulations of coupled cluster theory. The inverse of such an operator is known to exist, but a closed form explicit expression for the inverse is not available. We will address the evaluation of many-body similarity transformations generated by normal ordered exponential transformation operators without explicit use of the inverse. The similarity transform can be obtained as the solution of a linear system of equations that can be solved trivially using backward substitution. In addition a closed form diagrammatic expression for the similarity transformed operator is presented. Using the many-body similarity transformation strategy a simple and more general formulation of Fock space coupled cluster theory is presented which is akin in spirit to the formulation by Stolarczyk and Monkhorst [Phys. Rev. A 32, 725, 743 (1985); 37, 1908, 1926 (1988)], but which on the other hand is completely equivalent to the conventional wave operator formulation of Fock space coupled cluster theory (under suitable conditions). Other possible applications of the many-body similarity transformation will be briefly discussed.
The problem of finding the spectrum of the Fock-space Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{} for a system of many fermions is analyzed. The quasiparticle formalism is employed, with ``holes'' and ``particles'' defined with respect to some N-particle determinantal wave function (the model vacuum). The basic idea behind the proposed approach is to perform a similarity transformation of Hamiltonian \ifmmode \hat{H}\else \^{H}\fi{}, such that the resulting effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{} is, unlike \ifmmode \hat{H}\else \^{H}\fi{}, a quasiparticle-number--conserving operator. It is shown that eigenvalues of \ifmmode \hat{G}\else \^{G}\fi{}, corresponding to small numbers of quasiparticles (0,1,2) can be easily calculated. This is equivalent to finding eigenvalues of \ifmmode \hat{H}\else \^{H}\fi{} for certain states of N, N\ifmmode\pm\else\textpm\fi{}1, and N\ifmmode\pm\else\textpm\fi{}2 particles. The construction of the operator transforming \ifmmode \hat{H}\else \^{H}\fi{} into \ifmmode \hat{G}\else \^{G}\fi{} (the wave operator) stems from an analysis of the structure of the algebra of linear operators acting in a (finite-dimensional) Fock space. The exponential Ansatz for the wave operator is used, resulting in a generalization of the coupled-cluster (CC) method of Coester [Nucl. Phys. 7, 421 (1958)]. The generalized CC equations determining the wave operator, and equations determining the effective Hamiltonian \ifmmode \hat{G}\else \^{G}\fi{}, are presented in a diagrammatic form. An effort has been made to obtain a concise notation for expressing these equations in an algebraic form. Approximation schemes, necessary for practical applications of the proposed method, are also studied.
The concept of an effective Hamiltonian in Fock space is introduced. It is based on the division of the entire one-particle space into subspaces of ‘‘active’’ and ‘‘inactive’’ orbitals. The effective Fock space Hamiltonian has—for active model states—the same eigenvalues as the full Hamiltonian. The theory outlined in this context differs from that of paper I mainly in a different definition of the ‘‘diagonal part’’ of an operator, and in the fact that the ‘‘quasidegenerate case’’ applies throughout. The separation theorem, and as a consequence the connected diagram theorem, is shown to hold, in a more limited sense though, even for those normalizations where it did not in the context of universal wave and energy operators. Unlike in the theory of the ‘‘universal’’ operators of paper I the Fock space and n-particle Hilbert space approaches with analogous normalizations are no longer equivalent. In particular, the Primas normalization with a fully Lie-algebraic structure does not lead to a connected diagram expansion if it is formulated in n-particle Hilbert space, only so in a Fock space formulation. In n-particle Hilbert space with the present definition of the diagonal part of an operator the normalizations b (‘‘canonical’’) and c (‘‘Primas’’) happen to agree. As an alternative to the construction of the wave and energy operator W and L by perturbation theory the nonperturbative approach is presented as a generalization of the coupled-cluster method, in detail both in the intermediate and in the unitary normalization. In the unitary variant only a linear system for σ (the logarithm of the wave operator) has to be solved in order to get L correct through fifth order in perturbation theory with important contributions of higher orders included. A generalization of the Hartree–Fock method to Fock space theory is outlined, which guarantees stationarity of all (active) eigenstates with respect to one-particle transformations. A generalized electron pair theory is also defined. An analysis of the necessary computational steps shows that the nonperturbative approaches do not require significantly more computational effort than perturbation theory to the corresponding order. As a numerical example the H2 molecule in a small basis is 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.
We propose a modified coupled cluster Monte Carlo algorithm that stochastically samples connected terms within the truncated Baker-Campbell-Hausdorff expansion of the similarity-transformed Hamiltonian by construction of coupled cluster diagrams on the fly. Our new approach-diagCCMC-allows propagation to be performed using only the connected components of the similarity-transformed Hamiltonian, greatly reducing the memory cost associated with the stochastic solution of the coupled cluster equations. We show that for perfectly local, noninteracting systems diagCCMC is able to represent the coupled cluster wavefunction with a memory cost that scales linearly with system size. The favorable memory cost is observed with the only assumption of fixed stochastic granularity and is valid for arbitrary levels of coupled cluster theory. Significant reduction in memory cost is also shown to smoothly appear with dissociation of a finite chain of helium atoms. This approach is also shown not to break down in the presence of strong correlation through the example of a stretched nitrogen molecule. Our novel methodology moves the theoretical basis of coupled cluster Monte Carlo closer to deterministic approaches.