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:1409.2203 · Physical Review B

Lie algebraic similarity transformed Hamiltonians for lattice model systems

Jacob M. Wahlen-Strothman, Carlos A. Jiménez-Hoyos, Thomas M. Henderson, Gustavo E. Scuseria

Published Jan 23, 2015• 30 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.91.041114

Related papers

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

Biorthogonal projected energies of a Gutzwiller similarity transformed Hamiltonian

Jacob M. Wahlen-Strothman, Gustavo E. Scuseria

We present a method incorporating biorthogonal orbital-optimization, symmetry projection, and double-occupancy screening with a non-unitary similarity transformation generated by the Gutzwiller factor [Formula: see text], and apply it to the Hubbard model. Energies are calculated with mean-field computational scaling with high-quality results comparable to coupled cluster singles and doubles. This builds on previous work performing similarity transformations with more general, two-body Jastrow-style correlators. The theory is tested on 2D lattices ranging from small systems into the thermodynamic limit and is compared to available reference data.

arXiv:1610.02053•Journal of Physics Condensed Matter•Oct 5, 2016•7 citations

Compact numerical solutions to the two-dimensional repulsive Hubbard model obtained via nonunitary similarity transformations

Werner Dobrautz, Hongjun Luo, Ali Alavi

Similarity transformation of the Hubbard Hamiltonian using a Gutzwiller correlator leads to a non-Hermitian effective Hamiltonian, which can be expressed exactly in momentum-space representation and contains three-body interactions. We apply this methodology to study the two-dimensional Hubbard model with repulsive interactions near half filling in the intermediate interaction strength regime $(U/t=4)$. We show that at optimal or near optimal strength of the Gutzwiller correlator, the similarity-transformed Hamiltonian has extremely compact right eigenvectors, which can be sampled to high accuracy using the full configuration interaction quantum Monte Carlo (FCIQMC) method and its initiator approximation. Near-optimal correlators can be obtained using a simple projective equation, thus obviating the need for a numerical optimization of the correlator. The FCIQMC method, as a projective technique, is well suited for such non-Hermitian problems, and its stochastic nature can handle the three-body interactions exactly without undue increase in computational cost. The highly compact nature of the right eigenvectors means that the initiator approximation in FCIQMC is not severe and that large lattices can be simulated, well beyond the reach of the method applied to the original Hubbard Hamiltonian. Results are provided in lattice sizes up to 50 sites and compared to auxiliary-field QMC. New benchmark results are provided in the off-half-filling regime, with no severe sign problem being encountered. In addition, we show that methodology can be used to calculate excited states of the Hubbard model and lay the groundwork for the calculation of observables other than the energy.

arXiv:1811.03607•Physical review. B./Physical review. B•Feb 8, 2019•67 citations

Scalable Implementation of Mean-Field and Correlation Methods Based on Lie-Algebraic Similarity Transformation of Spin Hamiltonians in the Jordan–Wigner Representation

Shadan Ghassemi Tabrizi, Thomas M. Henderson, Thomas D. Kühne +1 more

High Resolution Image Download MS PowerPoint Slide Recent work has highlighted that the strong correlation inherent in spin Hamiltonians can be effectively reduced by mapping spins to Fermions via the Jordan–Wigner transformation (JW). The Hartree–Fock method is straightforward in the Fermionic domain and may provide a reasonable approximation to the ground state. Correlation with respect to the Fermionic mean field can be recovered based on Lie-algebraic similarity transformation (LAST) with two-body correlators. Specifically, a unitary LAST variant eliminates the dependence on site ordering, while a nonunitary LAST yields size-extensive correlation energies. Whereas the first recent demonstration of such methods was restricted to small spin systems, we present efficient implementations using analytical gradients for the optimization with respect to the mean-field reference and the LAST parameters, thereby enabling the treatment of larger clusters, including systems with local spins s > 1 2 .

Journal of Chemical Theory and Computation•Oct 31, 2025•2 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

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 +2 more

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.

arXiv:1512.06111•Physical review. B./Physical review. B•Mar 14, 2016•61 citations
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 ni\ensuremath↑ni\ensuremath↓{n}_{i\ensuremath{\uparrow}}{n}_{i\ensuremath{\downarrow}}ni\ensuremath↑​ni\ensuremath↓​, and two-site products of density (ni\ensuremath↑+ni\ensuremath↓)({n}_{i\ensuremath{\uparrow}}+{n}_{i\ensuremath{\downarrow}})(ni\ensuremath↑​+ni\ensuremath↓​) and spin (ni\ensuremath↑\ensuremath−ni\ensuremath↓)({n}_{i\ensuremath{\uparrow}}\ensuremath{-}{n}_{i\ensuremath{\downarrow}})(ni\ensuremath↑​\ensuremath−ni\ensuremath↓​) 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.
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.
We present a method incorporating biorthogonal orbital-optimization, symmetry projection, and double-occupancy screening with a non-unitary similarity transformation generated by the Gutzwiller factor [Formula: see text], and apply it to the Hubbard model. Energies are calculated with mean-field computational scaling with high-quality results comparable to coupled cluster singles and doubles. This builds on previous work performing similarity transformations with more general, two-body Jastrow-style correlators. The theory is tested on 2D lattices ranging from small systems into the thermodynamic limit and is compared to available reference data.
Similarity transformation of the Hubbard Hamiltonian using a Gutzwiller correlator leads to a non-Hermitian effective Hamiltonian, which can be expressed exactly in momentum-space representation and contains three-body interactions. We apply this methodology to study the two-dimensional Hubbard model with repulsive interactions near half filling in the intermediate interaction strength regime (U/t=4)(U/t=4)(U/t=4). We show that at optimal or near optimal strength of the Gutzwiller correlator, the similarity-transformed Hamiltonian has extremely compact right eigenvectors, which can be sampled to high accuracy using the full configuration interaction quantum Monte Carlo (FCIQMC) method and its initiator approximation. Near-optimal correlators can be obtained using a simple projective equation, thus obviating the need for a numerical optimization of the correlator. The FCIQMC method, as a projective technique, is well suited for such non-Hermitian problems, and its stochastic nature can handle the three-body interactions exactly without undue increase in computational cost. The highly compact nature of the right eigenvectors means that the initiator approximation in FCIQMC is not severe and that large lattices can be simulated, well beyond the reach of the method applied to the original Hubbard Hamiltonian. Results are provided in lattice sizes up to 50 sites and compared to auxiliary-field QMC. New benchmark results are provided in the off-half-filling regime, with no severe sign problem being encountered. In addition, we show that methodology can be used to calculate excited states of the Hubbard model and lay the groundwork for the calculation of observables other than the energy.
High Resolution Image Download MS PowerPoint Slide Recent work has highlighted that the strong correlation inherent in spin Hamiltonians can be effectively reduced by mapping spins to Fermions via the Jordan–Wigner transformation (JW). The Hartree–Fock method is straightforward in the Fermionic domain and may provide a reasonable approximation to the ground state. Correlation with respect to the Fermionic mean field can be recovered based on Lie-algebraic similarity transformation (LAST) with two-body correlators. Specifically, a unitary LAST variant eliminates the dependence on site ordering, while a nonunitary LAST yields size-extensive correlation energies. Whereas the first recent demonstration of such methods was restricted to small spin systems, we present efficient implementations using analytical gradients for the optimization with respect to the mean-field reference and the LAST parameters, thereby enabling the treatment of larger clusters, including systems with local spins s > 1 2 .
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 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.