No Arabic abstract
Coupled cluster (CC) has established itself as a powerful theory to study correlated quantum many-body systems. Finite temperature generalizations of CC theory have attracted considerable interest and have been shown to work as well as the ground-sate theory. However, most of these recent developments address only fermionic or bosonic systems. The distinct structure of the $su(2)$ algebra requires the development of a similar thermal CC theory for spin degrees of freedom. In this paper, we provide a formulation of our thermofield-inspired thermal CC for SU(2) systems. We apply the thermal CC to the Lipkin-Meshkov-Glick system as well as the one-dimensional transverse field Ising model as benchmark applications to highlight the accuracy of thermal CC in the study of finite-temperature phase diagram in SU(2) systems.
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 $J_1-J_2$ 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.
We describe the use of coupled-cluster theory as an impurity solver in dynamical mean-field theory (DMFT) and its cluster extensions. We present numerical results at the level of coupled-cluster theory with single and double excitations (CCSD) for the density of states and self-energies of cluster impurity problems in the one- and two-dimensional Hubbard models. Comparison to exact diagonalization shows that CCSD produces accurate density of states and self-energies at a variety of values of $U/t$ and filling fractions. However, the low cost allows for the use of many bath sites, which we define by a discretization of the hybridization directly on the real frequency axis. We observe convergence of dynamical quantities using approximately 30 bath sites per impurity site, with our largest 4-site cluster DMFT calculation using 120 bath sites. We suggest coupled cluster impurity solvers will be attractive in ab initio formulations of dynamical mean-field theory.
We propose an orbital optimized method for unitary coupled cluster theory (OO-UCC) within the variational quantum eigensolver (VQE) framework for quantum computers. OO-UCC variationally determines the coupled cluster amplitudes and also molecular orbital coefficients. Owing to its fully variational nature, first-order properties are readily available. This feature allows the optimization of molecular structures in VQE without solving any additional equations. Furthermore, the method requires smaller active space and shallower quantum circuit than UCC to achieve the same accuracy. We present numerical examples of OO-UCC using quantum simulators, which include the geometry optimization of the water and ammonia molecules using analytical first derivatives of the VQE.
We investigate the performance of Greens function coupled cluster singles and doubles (CCSD) method as a solver for Greens function embedding methods. To develop an efficient CC solver, we construct the one-particle Greens function from the coupled cluster (CC) wave function based on a non-hermitian Lanczos algorithm. The major advantage of this method is that its scaling does not depend on the number of frequency points. We have tested the applicability of the CC Greens function solver in the weakly to strongly correlated regimes by employing it for a half-filled 1D Hubbard model projected onto a single site impurity problem and a half-filled 2D Hubbard model projected onto a 4-site impurity problem. For the 1D Hubbard model, for all interaction strengths, we observe an excellent agreement with the full configuration interaction (FCI) technique, both for the self-energy and spectral function. For the 2D Hubbard, we have employed an open-shell version of the current implementation and observed some discrepancies from FCI in the strongly correlated regime. Finally, on an example of a small ammonia cluster, we analyze the performance of the Greens function CCSD solver within the self-energy embedding theory (SEET) with Hartee-Fock (HF) and Greens function second order (GF2) for the treatment of the environment.
Two integrable quantum spin ladder systems will be introduced associated with the fundamental su(2|2) solution of the Yang-Baxter equation. The first model is a generalized quantum Ising system with Ising rung interactions. In the second model the addition of extra interactions allows us to impose Heisenberg rung interactions without violating integrability. The existence of a Bethe ansatz solution for both models allows us to investigate the elementary excitations for antiferromagnetic rung couplings. We find that the first model does not show a gap whilst in the second case there is a gap for all positive values of the rung coupling.