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:1512.06111 · Physical review. B./Physical review. B

Polynomial similarity transformation theory: A smooth interpolation between coupled cluster doubles and projected BCS applied to the reduced BCS Hamiltonian

Matthias Degroote, Thomas M. Henderson, Jinmo Zhao, J. Dukelsky, Gustavo E. Scuseria

Published Mar 14, 2016• 61 citations
arXivHTMLPDF

Abstract

We present a similarity transformation theory based on a polynomial form of a particle-hole pair excitation operator. In the weakly correlated limit, this polynomial becomes an exponential, leading to coupled cluster doubles. In the opposite strongly correlated limit, the polynomial becomes an extended Bessel expansion and yields the projected BCS wave function. In between, we interpolate using a single parameter. The effective Hamiltonian is non-Hermitian and this polynomial similarity transformation theory follows the philosophy of traditional coupled cluster, left projecting the transformed Hamiltonian onto subspaces of the Hilbert space in which the wave function variance is forced to be zero. Similarly, the interpolation parameter is obtained through minimizing the next residual in the projective hierarchy. We rationalize and demonstrate how and why coupled cluster doubles is ill suited to the strongly correlated limit, whereas the Bessel expansion remains well behaved. The model provides accurate wave functions with energy errors that in its best variant are smaller than 1% across all interaction strengths. The numerical cost is polynomial in system size and the theory can be straightforwardly applied to any realistic Hamiltonian.

DOI
10.1103/physrevb.93.125124

Related papers

Communication: Projected Hartree Fock theory as a polynomial similarity transformation theory of single excitations

Yiheng Qiu, Thomas M. Henderson, Gustavo E. Scuseria

Spin-projected Hartree-Fock is written as a particle-hole excitation ansatz over a symmetry-adapted reference determinant. Remarkably, this expansion has an analytic expression that we were able to decipher. While the form of the polynomial expansion is universal, the excitation amplitudes need to be optimized. This is equivalent to the optimization of orbitals in the conventional projected Hartree-Fock framework of non-orthogonal determinants. Using the inverse of the particle-hole expansion, we similarity transform the Hamiltonian in a coupled-cluster style theory. The left eigenvector of the non-Hermitian Hamiltonian is constructed in a similar particle-hole expansion fashion, and we show that to numerically reproduce variational projected Hartree-Fock results, one needs as many pair excitations in the bra as the number of strongly correlated entangled pairs in the system. This single-excitation polynomial similarity transformation theory is an alternative to our recently presented double excitation theory, but supports projected Hartree-Fock and coupled cluster simultaneously rather than interpolating between them.

arXiv:1608.00029•The Journal of Chemical Physics•Sep 19, 2016•33 citations

Lie algebraic similarity transformations: improving wavefunctions for weak and strong correlations

Jacob M. Wahlen-Strothman

We present a class of correlated wavefunctions generated by exponentials of two-body on-site Hermitian operators that can be evaluated with polynomial computational cost via a Hamiltonian similarity transformation. Wavefunctions of this form have been studied with variational Monte Carlo methods, but we present a formalism to perform non-stochastic calculations. The Hausdorff series generated by these Jastrow factors can be summed exactly without truncation resulting in a set of equations with polynomial computational cost. The correlators include the density-density, collinear spin-spin, spin-density cross terms, and on-site double occupancy operators. The resulting non-Hermitian many-body Hamiltonian can be solved in a biorthogonal mean-field approach with only a small set of correlation terms required for accurate calculations in systems with local interactions. Although the energy of the model is unbound, projective equations in the spirit of coupled cluster theory lead to well-defined solutions. The theory is tested on the one and two-dimensional repulsive Hubbard model where it yields accurate results for large systems with low computational cost. Symmetry projection methods are included to further improve the reference wavefunction and results under strong correlation without sacrificing good quantum numbers resulting in very accurate energies for small systems and producing a better ground state for the calculation of other properties.

•Jun 29, 2017•0 citations

Spin-projected generalized Hartree-Fock method as a polynomial of particle-hole excitations

Thomas M. Henderson, Gustavo E. Scuseria

The past several years have seen renewed interest in the use of symmetry-projected Hartree-Fock for the description of strong correlations. Unfortunately, these symmetry-projected mean-field methods do not adequately account for dynamic correlation. Presumably, this shortcoming could be addressed if one could combine symmetry-projected Hartree-Fock with a many-body method such as coupled-cluster theory, but this is by no means straightforward because the two techniques are formulated in very different ways. However, we have recently shown that the singlet ${S}^{2}$-projected unrestricted Hartree-Fock wave function can in fact be written in a coupled-cluster-like wave function. That is, the spin-projected unrestricted Hartree-Fock wave function can be written as a polynomial of a double-excitation operator acting on some closed-shell reference determinant. Here, we extend this result and show that the spin-projected generalized Hartree-Fock wave function (which has both ${S}^{2}$ and ${S}_{z}$ projection) is likewise a polynomial of low-order excitation operators acting on a closed-shell determinant and provide a closed-form expression for the resulting polynomial coefficients. The spin projection of the generalized Hartree-Fock wave function introduces connected triple and quadruple excitations which are absent when spin-projecting an unrestricted Hartree-Fock determinant. We include a few preliminary applications of the combination of this spin-projected Hartree-Fock and coupled-cluster theory to the Hubbard Hamiltonian and comment on generalizations of the methodology. Results here are not for production level, but a similarity-transformed theory that combines the two offers the promise of being accurate for both weak and strong correlation, and may offer significant improvements in the intermediate correlation regime where neither projected Hartree-Fock nor coupled cluster is particularly accurate.

arXiv:1706.01548•Physical Review A•Aug 4, 2017•14 citations

Coupled-cluster representation of Green function employing modified spectral resolutions of similarity transformed Hamiltonians

Karol Kowalski, Kiran Bhaskaran‐Nair, W. A. Shelton

In this paper we discuss a new formalism for producing an analytic coupled-cluster (CC) Green's function for an N-electron system by shifting the poles of similarity transformed Hamiltonians represented in N - 1 and N + 1 electron Hilbert spaces. Simple criteria are derived for the states in N - 1 and N + 1 electron spaces that are then corrected in the spectral resolution of the corresponding matrix representations of the similarity transformed Hamiltonian. The accurate description of excited state processes within a Green's function formalism would be of significant importance to a number of scientific communities ranging from physics and chemistry to engineering and the biological sciences. This is because the Green's function methodology provides a direct path for not only calculating properties whose underlying origins come from coupled many-body interactions but also provides a straightforward path for calculating electron transport, response, and correlation functions that allows for a direct link with experiment. As a special case of this general formulation, we discuss the application of this technique for Green's function defined by the CC with singles and doubles representation of the ground-state wave function.

The Journal of Chemical Physics•Sep 2, 2014•42 citations

Non-Hermitian Green’s function theory with N -body interactions: The coupled-cluster similarity transformation

Christopher J. N. Coveney, David P. Tew

We present the diagrammatic theory of the irreducible self-energy and Bethe-Salpeter kernel that naturally arises within the Green’s function formalism for a general N -body non-Hermitian interaction. In this work, we focus specifically on the coupled-cluster (CC) self-energy generated by the similarity transformation of the electronic structure Hamiltonian. The main result of this work is a Green’s function formalism for non-Hermitian many-body quantum systems. We develop the biorthogonal quantum theory to construct dynamical correlation functions where the time dependence of operators is governed by a non-Hermitian Hamiltonian. We extend the Gell-Mann and Low theorem to include non-Hermitian interactions and to generate perturbative expansions of many-body Green’s functions. We introduce the single-particle coupled-cluster Green’s function and derive the perturbative diagrammatic expansion for the non-Hermitian coupled-cluster self-energy in terms of the “noninteracting” reference Green’s function, Σ ̃ [ G 0 ] . From the exact equation of motion of the single-particle coupled-cluster Green’s function, we derive the self-consistent renormalized coupled-cluster self-energy, Σ ̃ [ G ̃ ] , and demonstrate its relationship to the perturbative expansion of the self-energy, Σ ̃ [ G 0 ] . Subsequently, we show that the usual electronic self-energy can be recovered from the coupled-cluster self-energy by neglecting the effects of the similarity transformation. We show how the coupled-cluster ground-state energy can be obtained from the coupled-cluster self-energy and provide an overview of the relationship between approximations for the coupled-cluster self-energy, ionization potential/electron affinity equation-of-motion coupled-cluster theory, and the G 0 W 0 approximation. As a result, we introduce the CC- G 0 W 0 self-energy by leveraging the connections between Green’s function and coupled-cluster theory. Finally, we derive the diagrammatic expansion of the coupled-cluster Bethe-Salpeter kernel.

arXiv:2503.06586•Physical Review Research•Feb 17, 2026•4 citations

Lie algebraic similarity transformed Hamiltonians for lattice model systems

Jacob M. Wahlen-Strothman, Carlos A. Jiménez-Hoyos, Thomas M. Henderson +1 more

We present a class of Lie algebraic similarity transformations generated by exponentials of two-body on-site Hermitian operators whose Hausdorff series can be summed exactly without truncation. The correlators are defined over the entire lattice and include the Gutzwiller factor ${n}_{i\ensuremath{\uparrow}}{n}_{i\ensuremath{\downarrow}}$, and two-site products of density $({n}_{i\ensuremath{\uparrow}}+{n}_{i\ensuremath{\downarrow}})$ and spin $({n}_{i\ensuremath{\uparrow}}\ensuremath{-}{n}_{i\ensuremath{\downarrow}})$ operators. The resulting non-Hermitian many-body Hamiltonian can be solved in a biorthogonal mean-field approach with polynomial computational cost. The proposed similarity transformation generates locally weighted orbital transformations of the reference determinant. Although the energy of the model is unbound, projective equations in the spirit of coupled cluster theory lead to well-defined solutions. The theory is tested on the one- and two-dimensional repulsive Hubbard model where it yields accurate results for small and medium sized interaction strengths.

arXiv:1409.2203•Physical Review B•Jan 23, 2015•30 citations
We present a similarity transformation theory based on a polynomial form of a particle-hole pair excitation operator. In the weakly correlated limit, this polynomial becomes an exponential, leading to coupled cluster doubles. In the opposite strongly correlated limit, the polynomial becomes an extended Bessel expansion and yields the projected BCS wave function. In between, we interpolate using a single parameter. The effective Hamiltonian is non-Hermitian and this polynomial similarity transformation theory follows the philosophy of traditional coupled cluster, left projecting the transformed Hamiltonian onto subspaces of the Hilbert space in which the wave function variance is forced to be zero. Similarly, the interpolation parameter is obtained through minimizing the next residual in the projective hierarchy. We rationalize and demonstrate how and why coupled cluster doubles is ill suited to the strongly correlated limit, whereas the Bessel expansion remains well behaved. The model provides accurate wave functions with energy errors that in its best variant are smaller than 1% across all interaction strengths. The numerical cost is polynomial in system size and the theory can be straightforwardly applied to any realistic Hamiltonian.
Spin-projected Hartree-Fock is written as a particle-hole excitation ansatz over a symmetry-adapted reference determinant. Remarkably, this expansion has an analytic expression that we were able to decipher. While the form of the polynomial expansion is universal, the excitation amplitudes need to be optimized. This is equivalent to the optimization of orbitals in the conventional projected Hartree-Fock framework of non-orthogonal determinants. Using the inverse of the particle-hole expansion, we similarity transform the Hamiltonian in a coupled-cluster style theory. The left eigenvector of the non-Hermitian Hamiltonian is constructed in a similar particle-hole expansion fashion, and we show that to numerically reproduce variational projected Hartree-Fock results, one needs as many pair excitations in the bra as the number of strongly correlated entangled pairs in the system. This single-excitation polynomial similarity transformation theory is an alternative to our recently presented double excitation theory, but supports projected Hartree-Fock and coupled cluster simultaneously rather than interpolating between them.
We present a class of correlated wavefunctions generated by exponentials of two-body on-site Hermitian operators that can be evaluated with polynomial computational cost via a Hamiltonian similarity transformation. Wavefunctions of this form have been studied with variational Monte Carlo methods, but we present a formalism to perform non-stochastic calculations. The Hausdorff series generated by these Jastrow factors can be summed exactly without truncation resulting in a set of equations with polynomial computational cost. The correlators include the density-density, collinear spin-spin, spin-density cross terms, and on-site double occupancy operators. The resulting non-Hermitian many-body Hamiltonian can be solved in a biorthogonal mean-field approach with only a small set of correlation terms required for accurate calculations in systems with local interactions. Although the energy of the model is unbound, projective equations in the spirit of coupled cluster theory lead to well-defined solutions. The theory is tested on the one and two-dimensional repulsive Hubbard model where it yields accurate results for large systems with low computational cost. Symmetry projection methods are included to further improve the reference wavefunction and results under strong correlation without sacrificing good quantum numbers resulting in very accurate energies for small systems and producing a better ground state for the calculation of other properties.
The past several years have seen renewed interest in the use of symmetry-projected Hartree-Fock for the description of strong correlations. Unfortunately, these symmetry-projected mean-field methods do not adequately account for dynamic correlation. Presumably, this shortcoming could be addressed if one could combine symmetry-projected Hartree-Fock with a many-body method such as coupled-cluster theory, but this is by no means straightforward because the two techniques are formulated in very different ways. However, we have recently shown that the singlet S2{S}^{2}S2-projected unrestricted Hartree-Fock wave function can in fact be written in a coupled-cluster-like wave function. That is, the spin-projected unrestricted Hartree-Fock wave function can be written as a polynomial of a double-excitation operator acting on some closed-shell reference determinant. Here, we extend this result and show that the spin-projected generalized Hartree-Fock wave function (which has both S2{S}^{2}S2 and Sz{S}_{z}Sz​ projection) is likewise a polynomial of low-order excitation operators acting on a closed-shell determinant and provide a closed-form expression for the resulting polynomial coefficients. The spin projection of the generalized Hartree-Fock wave function introduces connected triple and quadruple excitations which are absent when spin-projecting an unrestricted Hartree-Fock determinant. We include a few preliminary applications of the combination of this spin-projected Hartree-Fock and coupled-cluster theory to the Hubbard Hamiltonian and comment on generalizations of the methodology. Results here are not for production level, but a similarity-transformed theory that combines the two offers the promise of being accurate for both weak and strong correlation, and may offer significant improvements in the intermediate correlation regime where neither projected Hartree-Fock nor coupled cluster is particularly accurate.
In this paper we discuss a new formalism for producing an analytic coupled-cluster (CC) Green's function for an N-electron system by shifting the poles of similarity transformed Hamiltonians represented in N - 1 and N + 1 electron Hilbert spaces. Simple criteria are derived for the states in N - 1 and N + 1 electron spaces that are then corrected in the spectral resolution of the corresponding matrix representations of the similarity transformed Hamiltonian. The accurate description of excited state processes within a Green's function formalism would be of significant importance to a number of scientific communities ranging from physics and chemistry to engineering and the biological sciences. This is because the Green's function methodology provides a direct path for not only calculating properties whose underlying origins come from coupled many-body interactions but also provides a straightforward path for calculating electron transport, response, and correlation functions that allows for a direct link with experiment. As a special case of this general formulation, we discuss the application of this technique for Green's function defined by the CC with singles and doubles representation of the ground-state wave function.
We present the diagrammatic theory of the irreducible self-energy and Bethe-Salpeter kernel that naturally arises within the Green’s function formalism for a general N -body non-Hermitian interaction. In this work, we focus specifically on the coupled-cluster (CC) self-energy generated by the similarity transformation of the electronic structure Hamiltonian. The main result of this work is a Green’s function formalism for non-Hermitian many-body quantum systems. We develop the biorthogonal quantum theory to construct dynamical correlation functions where the time dependence of operators is governed by a non-Hermitian Hamiltonian. We extend the Gell-Mann and Low theorem to include non-Hermitian interactions and to generate perturbative expansions of many-body Green’s functions. We introduce the single-particle coupled-cluster Green’s function and derive the perturbative diagrammatic expansion for the non-Hermitian coupled-cluster self-energy in terms of the “noninteracting” reference Green’s function, Σ ̃ [ G 0 ] . From the exact equation of motion of the single-particle coupled-cluster Green’s function, we derive the self-consistent renormalized coupled-cluster self-energy, Σ ̃ [ G ̃ ] , and demonstrate its relationship to the perturbative expansion of the self-energy, Σ ̃ [ G 0 ] . Subsequently, we show that the usual electronic self-energy can be recovered from the coupled-cluster self-energy by neglecting the effects of the similarity transformation. We show how the coupled-cluster ground-state energy can be obtained from the coupled-cluster self-energy and provide an overview of the relationship between approximations for the coupled-cluster self-energy, ionization potential/electron affinity equation-of-motion coupled-cluster theory, and the G 0 W 0 approximation. As a result, we introduce the CC- G 0 W 0 self-energy by leveraging the connections between Green’s function and coupled-cluster theory. Finally, we derive the diagrammatic expansion of the coupled-cluster Bethe-Salpeter kernel.
We present a class of Lie algebraic similarity transformations generated by exponentials of two-body on-site Hermitian operators whose Hausdorff series can be summed exactly without truncation. The correlators are defined over the entire lattice and include the Gutzwiller factor ${n}_{i\ensuremath{\uparrow}}{n}_{i\ensuremath{\downarrow}}$, and two-site products of density $({n}_{i\ensuremath{\uparrow}}+{n}_{i\ensuremath{\downarrow}})$ and spin $({n}_{i\ensuremath{\uparrow}}\ensuremath{-}{n}_{i\ensuremath{\downarrow}})$ operators. The resulting non-Hermitian many-body Hamiltonian can be solved in a biorthogonal mean-field approach with polynomial computational cost. The proposed similarity transformation generates locally weighted orbital transformations of the reference determinant. Although the energy of the model is unbound, projective equations in the spirit of coupled cluster theory lead to well-defined solutions. The theory is tested on the one- and two-dimensional repulsive Hubbard model where it yields accurate results for small and medium sized interaction strengths.