The recent quantum information boom has effected a resurgence of interest in the unitary coupled cluster (UCC) theory. Our group’s interest in local energy landscapes of unitary ansätze prompted us to investigate the classical approach of truncating the Taylor series expansion (instead of a perturbative expansion) of the UCC with singles and doubles (UCCSD) energy at the second order. This amounts to an approach where the electron correlation energy is estimated by taking a single Newton–Raphson step from Hartree–Fock toward UCCSD. Such an approach has been explored previously, but the accuracy was not extensively studied. In this paper, we investigate the performance and observe similar pathologies to the linearized coupled cluster with singles and doubles. We introduce the use of derivatives of order three or greater to help partially recover the variational lower bound of true UCCSD, restricting these derivatives to those of the “unmixed” category in order to simplify the model. By testing the approach on several potential energy surfaces and reaction energies, we find this “diagonal” approximation to higher order terms to be effective at reducing sensitivity near singularities for strongly correlated regimes, while not significantly diminishing the accuracy of weakly correlated systems.
The recent quantum information boom has effected a resurgence of interest in the unitary coupled cluster (UCC) theory. Our group’s interest in local energy landscapes of unitary ansätze prompted us to investigate the classical approach of truncating the Taylor series expansion (instead of a perturbative expansion) of the UCC with singles and doubles (UCCSD) energy at the second order. This amounts to an approach where the electron correlation energy is estimated by taking a single Newton–Raphson step from Hartree–Fock toward UCCSD. Such an approach has been explored previously, but the accuracy was not extensively studied. In this paper, we investigate the performance and observe similar pathologies to the linearized coupled cluster with singles and doubles. We introduce the use of derivatives of order three or greater to help partially recover the variational lower bound of true UCCSD, restricting these derivatives to those of the “unmixed” category in order to simplify the model. By testing the approach on several potential energy surfaces and reaction energies, we find this “diagonal” approximation to higher order terms to be effective at reducing sensitivity near singularities for strongly correlated regimes, while not significantly diminishing the accuracy of weakly correlated systems.
The unitary coupled cluster (UCC) approach has been one of the most popular wavefunction parametrizations for the variational quantum eigensolver due to the relative ease of optimization compared to hardware-efficient ansätze. In this contribution, we explore the energy landscape of the unitary coupled cluster singles and doubles (UCCSD) wavefunction for two commonly employed benchmark systems, lithium hydride and the nitrogen dimer. We investigate the organization of the solution space in terms of local minima and show how it changes as the number and order of operators of the UCC ansatz are varied. Surprisingly, we find that in all cases, the UCCSD energy has numerous low-lying minima connected by high energy transition states. Additionally, the energy spread of the minima that lie in the same band as the global minimum may exceed chemical accuracy, making the location of the true global minimum especially challenging.
The Trotterized Unitary Coupled Cluster Single and Double (UCCSD) ansatz has recently attracted interest due to its use in Variation Quantum Eigensolver (VQE) molecular simulations on quantum computers. However, when the size of molecules increases, UCCSD becomes less interesting as it cannot achieve sufficient accuracy. Including higher-order excitations is therefore mandatory to recover the UCC's missing correlation effects. In this Letter, we extend the Trotterized UCC approach via the addition of (true) Triple T excitations introducing UCCSDT. We also include both spin and orbital symmetries. Indeed, in practice, these later help to reduce unnecessarily circuit excitations and thus accelerate the optimization process enabling to tackle larger molecules. Our initial numerical tests (12-14 qubits) show that UCCSDT improves the overall accuracy by at least two-orders of magnitudes with respect to standard UCCSD. Overall, the UCCSDT ansatz is shown to reach chemical accuracy and to be competitive with the CCSD(T) gold-standard classical method of quantum chemistry.
The Trotterized Unitary Coupled Cluster Single and Double (UCCSD) ansatz has recently attracted interest due to its use in Variation Quantum Eigensolver (VQE) molecular simulations on quantum computers. However, when the size of molecules increases, UCCSD becomes less interesting as it cannot achieve sufficient accuracy. Including higher-order excitations is therefore mandatory to recover the UCC's missing correlation effects. Here, we extend the Trotterized UCC approach via the addition of (true) Triple T excitations introducing UCCSDT. We also include both spin and orbital symmetries. Indeed, in practice, the latter help to reduce unnecessary circuit excitations and thus accelerate the optimization process enabling researchers to tackle larger molecules. Our initial numerical tests (12-14 qubits) show that UCCSDT improves the overall accuracy by at least two orders of magnitude with respect to standard UCCSD. Overall, the UCCSDT ansatz is shown to reach chemical accuracy and to be competitive with the CCSD(T) gold-standard classical method of quantum chemistry.
Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.
Unitary coupled cluster (UCC), originally developed as a variational alternative to the popular traditional coupled cluster method, has seen a resurgence as a functional form for use on quantum computers. However, the number of excitors present in the Ansatz often presents a barrier to implementation on quantum computers. Given the natural sparsity of wavefunctions obtained from quantum Monte Carlo methods, we consider here a stochastic solution to the UCC problem. Using the coupled cluster Monte Carlo framework, we develop cluster selection schemes that capture the structure of the UCC wavefunction, as well as its Trotterized approximation, and use these to solve the corresponding projected equations. Due to the fast convergence of the equations with order in the cluster expansion, this approach scales polynomially with the size of the system. Unlike traditional UCC implementations, our approach naturally produces a non-variational estimator for the energy in the form of the projected energy. For unitary coupled cluster singles and doubles (UCCSD) in small systems, we find that this agrees well with the expectation value of the energy and, in the case of two electrons, with full configuration interaction results. For the larger N2 system, the two estimators diverge, with the projected energy approaching the coupled cluster result, while the expectation value is close to results from traditional UCCSD.
Unitary coupled cluster (UCC), originally developed as a variational alternative to the popular traditional coupled cluster method, has seen a resurgence as a functional form for use on quantum computers. However, the number of excitors present in the Ansatz often presents a barrier to implementation on quantum computers. Given the natural sparsity of wavefunctions obtained from quantum Monte Carlo methods, we consider here a stochastic solution to the UCC problem. Using the coupled cluster Monte Carlo framework, we develop cluster selection schemes that capture the structure of the UCC wavefunction, as well as its Trotterized approximation, and use these to solve the corresponding projected equations. Due to the fast convergence of the equations with order in the cluster expansion, this approach scales polynomially with the size of the system. Unlike traditional UCC implementations, our approach naturally produces a non-variational estimator for the energy in the form of the projected energy. For unitary coupled cluster singles and doubles (UCCSD) in small systems, we find that this agrees well with the expectation value of the energy and, in the case of two electrons, with full configuration interaction results. For the larger N2 system, the two estimators diverge, with the projected energy approaching the coupled cluster result, while the expectation value is close to results from traditional UCCSD.