Strongly correlated quantum systems, including frustrated spin lattices, remain challenging for conventional coupled-cluster methods because their ground states often exhibit strong multireference character and symmetry breaking. Although multireference and symmetry-projected coupled-cluster theories offer promising avenues, they are not yet black-box methods and often require substantial system-specific development. Spin Hamiltonians such as the XXZ and J1-J2 models provide compact but demanding benchmarks for these challenges. In particular, spin-BCS-based coupled-cluster methods may require high excitation ranks to recover the exact ground-state physics in strongly correlated regimes. In this work, we introduce a nonlinear spin-fermion-spin (SFS) transformation designed to map a strongly correlated spin Hamiltonian to a new spin Hamiltonian whose ground state is more accurately described at the mean-field level. For a representative SFS transformation, we identify a compact Clifford circuit that maps transformed product states to circuit-dressed, generally entangled states in the original representation. The transformed Hamiltonian can then be treated using broken-symmetry mean-field theory and post-mean-field methods. Benchmarks for the XXZ and J1-J2 Hamiltonians show that spin mean-field theory in the transformed representation reduces the energy error relative to conventional spin-BCS mean field. The spin-BCS-based coupled-cluster calculations for the one-dimensional XXZ model provide substantially improved accuracy relative to the corresponding untransformed methods. In 2D, the accuracy of the transformed mean-field solutions depends on the site ordering and boundary conditions because nonlocal strings appear. These results show that exact transformations of operator representations can serve as many-body preconditioners for mean-field and low-rank coupled-cluster approximations to strongly correlated spin systems.
Strongly correlated quantum systems, including frustrated spin lattices, remain challenging for conventional coupled-cluster methods because their ground states often exhibit strong multireference character and symmetry breaking. Although multireference and symmetry-projected coupled-cluster theories offer promising avenues, they are not yet black-box methods and often require substantial system-specific development. Spin Hamiltonians such as the XXZ and J1-J2 models provide compact but demanding benchmarks for these challenges. In particular, spin-BCS-based coupled-cluster methods may require high excitation ranks to recover the exact ground-state physics in strongly correlated regimes. In this work, we introduce a nonlinear spin-fermion-spin (SFS) transformation designed to map a strongly correlated spin Hamiltonian to a new spin Hamiltonian whose ground state is more accurately described at the mean-field level. For a representative SFS transformation, we identify a compact Clifford circuit that maps transformed product states to circuit-dressed, generally entangled states in the original representation. The transformed Hamiltonian can then be treated using broken-symmetry mean-field theory and post-mean-field methods. Benchmarks for the XXZ and J1-J2 Hamiltonians show that spin mean-field theory in the transformed representation reduces the energy error relative to conventional spin-BCS mean field. The spin-BCS-based coupled-cluster calculations for the one-dimensional XXZ model provide substantially improved accuracy relative to the corresponding untransformed methods. In 2D, the accuracy of the transformed mean-field solutions depends on the site ordering and boundary conditions because nonlocal strings appear. These results show that exact transformations of operator representations can serve as many-body preconditioners for mean-field and low-rank coupled-cluster approximations to strongly correlated spin systems.
Coupled-cluster theory is the method of choice for weakly correlated systems, but in strongly correlated systems where the mean-field reference is qualitatively poor, it can break down badly. These failures can be ameliorated by using symmetry-broken mean-field references. For systems that spontaneously break spin symmetry, we might favor an unrestricted reference. In systems that instead break number symmetry, we turn to Bogoliubov coupled cluster theory. In this work, we apply Bogoliubov coupled-cluster to Hamiltonians with an underlying su(2) algebra to study the pairing Hamiltonian and the spin XXZ and J1-J2 models by exploiting the correspondence between fermionic-pair operators and spin-1/2 operators. For the pairing problem, the reference is the usual Bardeen--Cooper--Schrieffer (BCS) quasiparticle vacuum, while for spin systems we use the analogous spin-BCS reference, which breaks Sz symmetry and provides a flexible starting point for strongly correlated regimes. Correlation is incorporated through a hierarchy of coupled-cluster approximations. Beyond energies, we compute the response density matrices for the evaluation of spin--spin correlation functions for the XXZ and J1-J2 Heisenberg models, and the pairing parameter and particle-number variance for the pairing Hamiltonian. We also derive the corresponding relaxed density matrices, including the effects of orbital relaxation. Benchmark calculations for the half-filled pairing Hamiltonian and spin models demonstrate that Bogoliubov coupled cluster provides a systematically improvable wave-function-based approach for su(2) Hamiltonians, including geometrically frustrated spin systems requiring complex and noncollinear BCS references.
Quantum chemistry aims to understand and predict molecular structure and properties by solving the time-independent Schr\"odinger equation within fixed nuclei (Born-Oppenheimer) approximation. The Hamiltonian may have continuous symmetries (e.g., spin) or discrete symmetries (for example, point group). Mean-field wave functions may break some or all of these symmetries. By taking the symmetry-adapted pieces of these broken symmetry states, one can restore the physical character of the wave function while also treating relevant strong correlations. On the other hand, mean-field methods lack dynamic correlation arising from short-range interactions. Still, single-reference methods such as configuration interaction or coupled cluster theory can handle these correlations efficiently. Ideally, we would like to combine these approaches: performing configuration interaction or coupled cluster theory on a broken symmetry reference while restoring the desired symmetries. This work introduces a new approach to integrating symmetry projection and coupled cluster theory. By borrowing techniques from previous research, the new method can restore continuous symmetries like U(1) and SU(2) and discrete symmetries like point group and time reversal. This method achieves accurate results in challenging systems and has the same computational scaling as traditional coupled cluster theory. Jordan-Wigner transformation is a potent tool in spin models, but its application is challenging due to the inherent strings. The same techniques we use to combine symmetry projection and coupled cluster theory can be used to combine Jordan-Wigner transformation and coupled cluster. The strings are the ' symmetry projection ' operator when constructing the Hartree-Fock and coupled cluster wave functions for a Jordan-Wigner transformed Hamiltonian. Leveraging the method of symmetry projection, we obtain excellent results for both 1-D and 2-D spin models.
We introduce Sz spin-projection based on cluster mean-field theory and apply it to the ground state of strongly correlated spin systems. In cluster mean-fields, the ground state wavefunction is written as a factorized tensor product of optimized cluster states. In previous work, we have focused on unrestricted cluster mean-field, where each cluster is Sz symmetry adapted. We here remove this restriction by introducing a generalized cluster mean-field (GcMF) theory, where each cluster is allowed to access all Sz sectors, breaking Sz symmetry. In addition, a projection scheme is used to restore global Sz, which gives rise to the Sz spin-projected generalized cluster mean-field (SzGcMF). Both of these extensions contribute to accounting for inter-cluster correlations. We benchmark these methods on the 1D, quasi-2D, and 2D J1 - J2 and XXZ Heisenberg models. Our results indicate that the new methods (GcMF and SzGcMF) provide a qualitative and semi-quantitative description of the Heisenberg lattices in the regimes considered, suggesting them as useful references for further inter-cluster correlations, which are discussed in this work.
We introduce cluster-based mean-field, perturbation and coupled-cluster theories to describe the ground state of strongly-correlated spin systems. In cluster mean-field, the ground state wavefunction is written as a simple tensor product of optimized cluster states. The cluster-language and the mean-field nature of the ansatz allows for a straightforward improvement based on perturbation theory and coupled-cluster, to account for inter-cluster correlations. We present benchmark calculations on the 2D square J1−J2 Heisenberg model, using cluster mean-field, second-order perturbation theory and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies very accurately. Our results indicate that, even with relatively small clusters, the correlated methods can provide an accurate description of the Heisenberg model in the regimes considered. Some ways to improve the results presented in this work are discussed.
The Jordan-Wigner transformation permits one to convert spin 1/2 operators into spinless fermion ones, or vice versa. In some cases, it transforms an interacting spin Hamiltonian into a noninteracting fermionic one, which is exactly solved at the mean-field level. Even when the resulting fermionic Hamiltonian is interacting, its mean-field solution can provide surprisingly accurate energies and correlation functions. Jordan-Wigner is, however, only one possible means of interconverting spin and fermionic degrees of freedom. Here, we apply several such techniques to the XXZ and J1-J2 Heisenberg models, as well as to the pairing or reduced Bardeen-Cooper-Schrieffer Hamiltonian, with the aim of discovering which of these mappings is most useful in applying fermionic mean-field theory to the study of spin Hamiltonians.
We introduce perturbation and coupled-cluster theories based on a cluster mean-field reference for describing the ground state of strongly correlated spin systems. In cluster mean-field, the ground state wave function is written as a simple tensor product of optimized cluster states. The cluster language and the mean-field nature of the ansatz allow for a straightforward improvement which uses perturbation theory and coupled-cluster to account for intercluster correlations. We present benchmark calculations on the 1D chain and 2D square J 1 – J 2 Heisenberg model, using cluster mean-field, perturbation theory, and coupled-cluster. We also present an extrapolation scheme that allows us to compute thermodynamic limit energies accurately. Our results indicate that, with sufficiently large clusters, the correlated methods (cPT2, cPT4, and cCCSD) can provide a relatively accurate description of the Heisenberg model in the regimes considered, which suggests that the methods presented can be used for other strongly correlated systems. Some ways to improve upon the methods presented in this work are discussed.