No Arabic abstract
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.
The Coupled Cluster (CC) method is used to compute the electronic correlation energy in atoms and molecules and often leads to highly accurate results. However, due to its single-reference nature, standard CC in its projected form fails to describe quantum states characterized by strong electronic correlations and multi-reference projective methods become necessary. On the other hand, quantum algorithms for the solution of many-electron problems have also emerged recently. The quantum UCC with singles and doubles (q-UCCSD) is a popular wavefunction Ansatz for the Variational Quantum Eigensolver (VQE) algorithm. The variational nature of this approach can lead to significant advantages compared to its classical equivalent in the projected form, in particular for the description of strong electronic correlation. However, due to the large number of gate operations required in q-UCCSD, approximations need to be introduced in order to make this approach implementable in a state-of-the-art quantum computer. In this work, we evaluate several variants of the standard q-UCCSD Ansatz in which only a subset of excitations is included. In particular, we investigate the singlet and pair q-UCCD approaches combined with orbital optimization. We show that these approaches can capture the dissociation/distortion profiles of challenging systems such as H$_4$, H$_2$O and N$_2$ molecules, as well as the one-dimensional periodic Fermi-Hubbard chain. These results promote the future use of q-UCC methods for the solution of challenging electronic structure problems in quantum chemistry.
The previously proposed ansatz for density cumulant theory that combines orbital-optimization and a parameterization of the 2-electron reduced density matrix cumulant in terms of unitary coupled cluster amplitudes (OUDCT) is carefully examined. Formally, we elucidate the relationship between OUDCT and orbital-optimized unitary coupled cluster theory and show the existence of near-zero denominators in the stationarity conditions for both the exact and some approximate OUDCT methods. We implement methods of the OUDCT ansatz restricted to double excitations for numerical study, up to the fifth commutator in the Baker-Campbell-Hausdorff expansion. We find that methods derived from the ansatz beyond the previously known ODC-12 method tend to be less accurate for equilibrium properties and less reliable when attempting to describe $H2$ dissociation. New developments are needed to formulate more accurate DCT variants.
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 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 present a quantum algorithm to compute the entanglement spectrum of arbitrary quantum states. The interesting universal part of the entanglement spectrum is typically contained in the largest eigenvalues of the density matrix which can be obtained from the lower Renyi entropies through the Newton-Girard method. Obtaining the $p$ largest eigenvalues ($lambda_1>lambda_2ldots>lambda_p$) requires a parallel circuit depth of $mathcal{O}(p(lambda_1/lambda_p)^p)$ and $mathcal{O}(plog(N))$ qubits where up to $p$ copies of the quantum state defined on a Hilbert space of size $N$ are needed as the input. We validate this procedure for the entanglement spectrum of the topologically-ordered Laughlin wave function corresponding to the quantum Hall state at filling factor $ u=1/3$. Our scaling analysis exposes the tradeoffs between time and number of qubits for obtaining the entanglement spectrum in the thermodynamic limit using finite-size digital quantum computers. We also illustrate the utility of the second Renyi entropy in predicting a topological phase transition and in extracting the localization length in a many-body localized system.