Bridging physical intuition and hardware efficiency for correlated electronic states: the local unitary cluster Jastrow ansatz for electronic structure
Mário Motta, Kevin J. Sung, K. Birgitta Whaley, Martin Head‐Gordon, James Shee
A prominent goal in quantum chemistry is to solve the molecular electronic structure problem for the ground state energy with high accuracy. While classical quantum chemistry is a relatively mature field, the accurate and scalable prediction of strongly correlated states found, e.g., in bond breaking and polynuclear transition metal compounds remains an open problem. Within the context of a variational quantum eigensolver, we propose a new family of ansatzes which provides a more physically appropriate description of strongly correlated electrons than unitary coupled cluster with singles and doubles excitations, with vastly reduced quantum resource requirements. Specifically, we present a set of local approximations to the unitary cluster Jastrow wavefunction motivated by Hubbard physics. The resulting ansatz removes the need for SWAP gates and can be tailored to arbitrary qubit topologies (e.g., square, hex, heavy-hex). Our ansatz is well-suited to take advantage of continuous sets of quantum gates recently realized on superconducting devices with tunable couplers, while retaining a unique level of physical transparency and interpretability. As the capabilities of quantum computing devices continue to progress, we expect that the local cluster Jastrow ansatz to be a natural choice to encode both statically and dynamically correlated electronic wavefunctions.
A prominent goal in quantum chemistry is to solve the molecular electronic structure problem for the ground state energy with high accuracy. While classical quantum chemistry is a relatively mature field, the accurate and scalable prediction of strongly correlated states found, e.g., in bond breaking and polynuclear transition metal compounds remains an open problem. Within the context of a variational quantum eigensolver, we propose a new family of ansatzes which provides a more physically appropriate description of strongly correlated electrons than unitary coupled cluster with singles and doubles excitations, with vastly reduced quantum resource requirements. Specifically, we present a set of local approximations to the unitary cluster Jastrow wavefunction motivated by Hubbard physics. The resulting ansatz removes the need for SWAP gates and can be tailored to arbitrary qubit topologies (e.g., square, hex, heavy-hex). Our ansatz is well-suited to take advantage of continuous sets of quantum gates recently realized on superconducting devices with tunable couplers, while retaining a unique level of physical transparency and interpretability. As the capabilities of quantum computing devices continue to progress, we expect that the local cluster Jastrow ansatz to be a natural choice to encode both statically and dynamically correlated electronic wavefunctions.
High Resolution Image Download MS PowerPoint Slide Quantum computers have the potential to efficiently solve the electronic structure problem but are currently limited by noise and shallow circuits. We present an enhanced Variational Quantum Eigensolver (VQE) ansatz based on the Qubit Coupled Cluster (QCC) approach that requires optimization of only n parameters, where n is the number of Pauli string generators, rather than the typical n + 2 m parameters, where m is the number of qubits. We evaluate the ground state energies and molecular dissociation curves of strongly correlated molecules, namely O 3 and Li 4, using active spaces of varying sizes in conjunction with our enhanced QCC ansatz, Unitary Coupled Cluster Single–Double (UCCSD) ansatz, and the classical Coupled Cluster Singles and Doubles (CCSD) method. Compared to UCCSD, our approach significantly reduces the number of parameters while maintaining high accuracy. Numerical simulations demonstrate the effectiveness of our approach, and experiments on superconducting and trapped-ion quantum computers showcase its practicality on real hardware. By eliminating the need for symmetry-restoring gates and reducing the number of parameters, our enhanced QCC ansatz enables accurate quantum chemistry calculations on near-term quantum devices for strongly correlated systems.
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 unitary cluster Jastrow (UCJ) ansatz and its variant known as local UCJ (LUCJ) are promising choices for variational quantum algorithms for chemistry due to their combination of physical motivation and hardware efficiency. The parameters of these ansatzes can be initialized from the output of a coupled cluster, singles and doubles (CCSD) calculation performed on a classical computer. However, truncating the number of repetitions of the ansatz, as well as discarding interactions to accommodate the connectivity constraints of near-term quantum processors, degrade the approximation to CCSD and the resulting energy accuracy. In this work, we propose two methods to improve the parameter initialization. The first method, which is applicable to both expectation value- and sample-based algorithms, uses compressed double factorization of the CCSD amplitudes to improve or recover the CCSD approximation. The second method, which is applicable to sample-based algorithms, uses approximate tensor network simulation to improve the quality of samples produced by the ansatz circuit. We validate our methods using exact state vector simulation on systems of up to 52 qubits, as well as experiments on superconducting quantum processors using up to 65 qubits. Our results indicate that our methods can significantly improve the output of both expectation value- and sample-based quantum algorithms.
In recent years, the Variational Quantum Eigensolver (VQE) has emerged as one of the most popular algorithms for solving the electronic structure problem on near-term quantum computers. The utility of VQE is often hindered by the limitations of current quantum hardware, including, but not limited to, short qubit coherence times and low gate fidelities. These limitations become particularly pronounced when VQE is used along with deep quantum circuits, such as those required by the ”Unitary Coupled Cluster Singles and Doubles” (UCCSD) ansatz, often resulting in significant errors. To address these issues, we propose a low-depth ansatz based on parallelized Givens rotations, which can recover substantial correlation energy while drastically reducing circuit depth and two-qubit gate counts for an arbitrary active space (AS). Also, considering the current hardware architectures with low qubit counts, we introduce a systematic way to select molecular orbitals to define active spaces (ASs) that retain significant electron correlation. We validate our approach by computing bond dissociation profiles of water and strongly correlated systems, such as molecular nitrogen and oxygen, across various ASs. Noiseless simulations using the new ansatz yield ground-state energies comparable to those from the UCCSD ansatz while reducing circuit depth by 50–70%. Moreover, in noisy simulations, our approach achieves energy error rates an order of magnitude lower than that of UCCSD. Considering the efficiency and practical usage of our ansatz, we hope that it becomes a potential choice for performing quantum chemistry calculations on near-term quantum devices.
Recent experiments on quantum computers have challenged the limits of classical computation in chemistry, simulating ground states of strongly correlated molecules. Many of these experiments have utilized the unitary cluster Jastrow ansatz, a quantum circuit inspired by the unitary coupled cluster ansatz that can be tailored to current quantum hardware. Notably, the largest experiment in Sci. Adv. 11, 25 (2025) executed a quantum circuit with 77 qubits and 10,570 gates on an IBM quantum computer and performed classical post-processing with up to 6400 nodes on Fugaku to compute ground state energies better than Hartree-Fock. In this work, we present a polynomial time classical algorithm to compute the energy of any single-layer unitary cluster Jastrow circuit, independent of locality constraints for quantum hardware. Our algorithm can reproduce the largest experiment from Sci. Adv. 11, 25 (2025) in less than a minute on a laptop, and through circuit optimization enabled by fast simulation we achieve a lower ground state energy than the experiment.
Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron-electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. The transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.