Abstract The iterative quantum phase estimation algorithm (IQPE) is theoretically appealing in its wide scope of being able to handle electronic correlation. However, the quality of the initial input state strongly enhances the probability of landing on the desired eigenstate. In this work, we systematically study two different parametrization schemes of the unitary coupled cluster (UCC) ansätz in the variational quantum eigensolver (VQE) framework toward the reference state preparation for IQPE. The efficacy of the UCC variants toward an appropriate state preparation is studied with prototypical H4 molecule on a circle. While the conventional UCC ansätz can lead to high success probability across various degrees of electronic complexity, a resource efficient minimally parametrized UCC ansätz consisting of active space excitations is shown to incorporate the essential static correlation in the reference state description. We demonstrate that such a carefully prepared initial state can significantly reduce the effects of noise due to sampling in the estimation of the desired eigenphase.
Abstract The iterative quantum phase estimation algorithm (IQPE) is theoretically appealing in its wide scope of being able to handle electronic correlation. However, the quality of the initial input state strongly enhances the probability of landing on the desired eigenstate. In this work, we systematically study two different parametrization schemes of the unitary coupled cluster (UCC) ansätz in the variational quantum eigensolver (VQE) framework toward the reference state preparation for IQPE. The efficacy of the UCC variants toward an appropriate state preparation is studied with prototypical H4 molecule on a circle. While the conventional UCC ansätz can lead to high success probability across various degrees of electronic complexity, a resource efficient minimally parametrized UCC ansätz consisting of active space excitations is shown to incorporate the essential static correlation in the reference state description. We demonstrate that such a carefully prepared initial state can significantly reduce the effects of noise due to sampling in the estimation of the desired eigenphase.
Digital quantum simulation of strong correlation effects with iterative quantum phase estimation over the variational quantum eigensolver algorithm: H4 on a circle as a case study
The iterative quantum phase estimation algorithm, applied to calculating the ground state energies of quantum chemical systems, is theoretically appealing in its wide scope of being able to handle both weakly and strongly correlated regimes. However, the goodness of the initial state that is sent as an input to the algorithm could strongly decide the quality of the results obtained. In this work, we generate the initial state by using the classical-quantum hybrid variational quantum eigensolver algorithm with unitary coupled cluster ansatz. We apply the procedure to obtain the ground state energies of the H4 molecule on a circle, as the system exhibits an interplay of dynamic as well as static correlation effects at different geometries. Furthermore, we argue on the importance of static correlation in construction of the reference determinant, and propose a minimally parametrized unitary coupled cluster ansatz, which drastically reduces number of variational parameters while incorporating the static correlation effects in the wavefunction. We demonstrate that a carefully and appropriately prepared initial state can greatly reduce the effects of noise due to sampling in the estimation of the desired eigenphase.
We present a novel approach for efficient preparation of electronic ground states, leveraging the optimizer ExcitationSolve [Jäger et al., Comm. Phys. (2025)] and established variational quantum eigensolver-based operator selection methods, such as Energy Sorting (ES). By combining these tools, we demonstrate a computationally efficient protocol that enables the construction of an approximate ground state from a unitary coupled cluster ansatz via a single sweep over the operator pool. Utilizing efficient classical pre-processing to select the majority of relevant operators, this approach reduces the computational complexity associated with traditional variational quantum eigensolver (VQE) optimization methods. We further show that second-order Epstein-Nesbet (EN2) perturbation theory emerges as the first-order Taylor expansion of our protocol in terms of a correlation measure, clarifying why our approach provides a more robust initial guess for the ground state in strongly correlated regimes. We also find that second-order Møller-Plesset perturbation (MP2) theory, which is widely used for unitary coupled cluster (UCC) initialization, performs worse than both EN2 and our protocol. Furthermore, we show that our method can be seamlessly integrated with one-variational-parameter couple exchange operators, thereby further reducing the number of required CNOT operations. Overall, we empirically observe a quadratic convergence speedup beyond state-of-the-art methods, advancing the preparation of high-fidelity electronic ground states - one of the cornerstones of meaningful electronic structure calculations in the noisy intermediate-scale quantum computing (NISQ) era, and a prerequisite for fault-tolerant quantum computing (FTQC) algorithms such as quantum phase estimation.
Abstract The variational quantum eigensolver (VQE) algorithm combines the ability of quantum computers to efficiently compute expectation values with a classical optimization routine in order to approximate ground state energies of quantum systems. In this paper, we study the application of VQE to the simulation of molecular energies using the unitary coupled cluster (UCC) ansatz. We introduce new strategies to reduce the circuit depth for the implementation of UCC and improve the optimization of the wavefunction based on efficient classical approximations of the cluster amplitudes. Additionally, we propose an analytical method to compute the energy gradient that reduces the sampling cost for gradient estimation by several orders of magnitude compared to numerical gradients. We illustrate our methodology with numerical simulations for a system of four hydrogen atoms that exhibit strong correlation and show that the circuit depth of VQE using a UCC ansatz can be reduced without introducing significant loss of accuracy in the final wavefunctions and energies.
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 unitary variant of CC (UCC) with singles and doubles (q-UCCSD) is a popular wavefunction Ansatz for the variational quantum eigensolver 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 H4, H2O, and N2 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 factorized form of unitary coupled cluster theory (UCC) is a promising wave-function ansatz for the variational quantum eigensolver algorithm. Here, we present a quantum inspired algorithm for UCC based on an exact operator identity for the individual UCC factors. We implement this algorithm for calculations of the H10 linear chain and the H2O molecule with single and double ζ basis sets to provide insights about UCC as a wave-function ansatz. We find that as an electronic structure method, UCC could potentially be valuable for strongly correlated systems, for which conventional coupled cluster theory has difficulties. On quantum computers, the factorized form of UCC can also serve as an initial state preparation method for the quantum phase-estimation algorithm, since it yields higher overlap with the ground state than many other common variational ansatzes. This classical algorithm for UCC also enables widespread testing of the technique, which will be useful for benchmarking quantum computations.
This is a case study of the variational quantum eigensolver (VQE) method using numerical simulations to test the influence of noise on the accuracy of the underlying circuit ansatz. We investigate a computational chemistry application of VQE to calculate the electronic ground state and its energy for Sodium Hydride (NaH), a prototypical two-electron problem. Using a one-parameter ansatz derived from unitary coupled cluster (UCC) theory, we simulate the effects of noise on the energy expectation value and variance with respect to the ansatz parameter. These numerical simulations provide insights into the accuracy of the prepared quantum state and the efficiency of the classical optimizer that iteratively refines the ansatz. We conduct a comparative study between analytical results derived for the UCC ansatz in the absence of noise and the noisy numerical simulation results obtained using an isotropic depolarizing noise model for each gate. We also compare the relative increase in noise on logically equivalent UCC ansatz circuits generated by randomized compiling. Notably, we observe that the intrinsic variance in the energy due to the simplicity of the ansatz itself compares with the noise induced by the bare circuit.