Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Chemical Theory and Computation

Multicomponent Unitary Coupled Cluster and Equation-of-Motion for Quantum Computation

Fabijan Pavošević, Sharon Hammes‐Schiffer

Published May 4, 2021• 37 citations
Publisher

Abstract

The variational quantum eigensolver (VQE) algorithm combined with the unitary coupled cluster (UCC) ansatz has been developed for the quantum computation of molecular energies and wave functions within the Born–Oppenheimer approximation. Herein, this approach is extended to multicomponent systems to enable the quantum mechanical treatment of more than one type of particle, such as electrons and positrons or electrons and nuclei, without invoking the Born–Oppenheimer approximation. Specifically, we introduce the multicomponent unitary coupled cluster (mcUCC) method combined with the VQE algorithm for the calculation of ground-state energies and wave functions as well as the multicomponent equation-of-motion (mcEOM) method for the calculation of excitation energies. These methods are developed within the nuclear-electronic orbital (NEO) framework and are formulated in the qubit basis to enable implementations on quantum computers. Moreover, these methods are used to calculate the ground-state energy and excitation energies of positronium hydride, where both electrons and the positron are treated quantum mechanically, as well as the H 2 molecule, where both electrons and one proton are treated quantum mechanically. These applications validate the implementation and provide benchmark data for future calculations. The errors due to Trotterization of the mcUCC ansatz are also analyzed. This formalism, as well as the accompanying computer code, will serve as the basis for applications to more complex multicomponent systems, such as simulations of photoinduced nonadiabatic molecular processes, on both classical and quantum computers.

DOI
10.1021/acs.jctc.1c00220

Related papers

Error-Mitigation Enabled Multicomponent Quantum Simulations Beyond the Born-Oppenheimer Approximation

Delmar G. A. Cabral, Brandon Allen, Fabijan Pavošević +6 more

We introduce a multicomponent unitary coupled cluster framework for quantum simulations of molecular systems that incorporate both electronic and nuclear quantum effects beyond the Born-Oppenheimer approximation. Using the nuclear-electronic orbital formalism, we construct mcUCC ansätze for positronium hydride and molecular hydrogen with a quantum proton, and analyze hardware requirements for different excitation truncations. To further reduce resource costs effectively, we employ the local unitary cluster Jastrow ansatz and implement it experimentally on IBM Q's Heron superconducting hardware. With the Physics-Inspired Extrapolation error mitigation protocol, the computed ground-state energies remain within chemical accuracy, consistent with the stated uncertainty level. These results provide the first demonstration of error-mitigated multicomponent correlated simulations on quantum hardware and outline a path toward scalable algorithms unifying electronic and nuclear degrees of freedom.

quant-pharXiv:2511.11941•Journal of Chemical Theory and Computation•Nov 18, 2025•2 citations

Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz

Jonathan Romero, Ryan Babbush, Jarrod R. McClean +3 more

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.

arXiv:1701.02691•Quantum Science and Technology•Jul 17, 2018•726 citations

Quantum computing approach to atomic and molecular three-body systems

Mohammad Haidar, Hugo Demattos Nogueira, Jean‐Philippe Karr

We present high-precision quantum computing simulations of three-body atoms (He, H$^-$) and molecules (H$_2^+$, HD$^+$), the latter being studied beyond the Born-Oppenheimer approximation. The Non-Iterative Disentangled Unitary Coupled Cluster Variational Quantum Eigensolver (NI-DUCC-VQE) [M. Haidar et al., Quantum Sci. Technol. 10, 025031 (2025)] is used. By combining a first-quantized Hamiltonian with a Minimal Complete Pool (MCP) of Lie-algebraic excitations, we construct a compact ansatz with a gradient-independent construction, avoiding costly gradient evaluations and yielding efficient computational scaling with both basis size and electron number. It avoids barren plateaus and enables rapid convergence, achieving energy errors as low as 10$^{-11}$ a.u. with state fidelities only limited by arithmetic precision in only a few thousand function evaluations in all four systems. These results make three-body atoms and molecules excellent candidates for benchmarking and testing on current Noisy Intermediate-Scale Quantum (NISQ) devices. Further, our approach can be extended to more complex systems with larger basis sets, taking advantage of the efficient scaling of qubit requirements to study electronic correlations and non-adiabatic effects with high precision. We also demonstrate the applicability of NI-DUCC-VQE for simulating higher-order effects such as relativistic corrections and hyperfine interactions.

arXiv:2510.18005•Physical Review A•Apr 30, 2026•0 citations

Generalized unitary coupled cluster excitations for multireference molecular states optimized by the variational quantum eigensolver

Gabriel Greene‐Diniz, David Muñoz Ramo

Abstract The variational quantum eigensolver (VQE) algorithm, designed to calculate the energy of molecular ground states on near‐term quantum computers, requires specification of symmetries that describe the system, for example, spin state and number of electrons. This opens the possibility of using VQE to obtain excited states as the lowest‐energy solutions of a given set of symmetries. In this paper, the performances of various unitary coupled cluster (UCC) ansätze applied to VQE calculations on excited states are investigated using quantum circuits designed to represent single reference and multireference wavefunctions to calculate energy curves with respect to variations in the molecular geometry. These ansätze include standard Unitary Coupled Cluster Singles and Doubles (UCCSD), as well as modified versions of Unitary Coupled Cluster Generalized Singles and Doubles (UCCGSD) and k‐UpCCGSD, which are engineered to tackle excited states without undesired spin symmetry crossover to lower states during VQE optimization. These studies are carried out on a range of systems, including H2, H3, H4, NH, OH+, CH2, and , covering examples of spin singlet, doublet, and triplet molecular ground states with single and multireference excited states. In most cases, our calculations are in excellent agreement with results from full‐configuration interaction calculations on classical machines, thus showing that the VQE algorithm is capable of calculating the lowest excited state at a certain symmetry, including multireference closed‐ and open‐shell states, by setting appropriate restrictions on the excitations considered in the cluster operator and appropriate constraints in the qubit register encoding the starting mean field state.

arXiv:1910.05168•International Journal of Quantum Chemistry•Nov 27, 2020•42 citations

Quantum-Inspired Algorithm for the Factorized Form of Unitary Coupled Cluster Theory

Jia Chen, Hai‐Ping Cheng, James K. Freericks

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 classical algorithm for UCC based on an exact operator identity for the individual UCC factors. We implement this algorithm for calculations of the H 10 linear chain and the H 2 O molecule with single and double ζ basis sets to provide insights into UCC as a wave-function ansatz. We find that for weakly correlated molecules, the factorized form of the UCC provides similar accuracy to conventional coupled cluster theory (CC); for strongly correlated molecules, where CC often breaks down, UCC significantly outperforms the configuration interaction (CI) ansatz. As a result, the factorized form of the UCC is an accurate, efficient, and reliable electronic structure method in both the weakly and strongly correlated regions. This classical algorithm now allows robust benchmarking of anticipated results from quantum computers and application of coupled-cluster techniques to more strongly correlated molecules.

Journal of Chemical Theory and Computation•Jan 27, 2021•44 citations

Polaritonic Unitary Coupled Cluster for Quantum Computations

Fabijan Pavošević, Johannes Flick

In the field of polaritonic chemistry, strong light-matter interactions are used to alter chemical reactions inside optical cavities. To understand these processes, the development of reliable theoretical models is essential. While traditional methods have to balance accuracy and system size, new developments in quantum computing offer a path for accurate calculations on currently available quantum devices. Here, we introduce the quantum electrodynamics unitary coupled cluster (QED-UCC) method combined with the Variational Quantum Eigensolver algorithm, as well as the quantum electrodynamics equation-of-motion (QED-EOM) method formulated in the qubit basis that allow accurate calculations of ground-state and excited-state properties of strongly coupled light-matter systems suitable for quantum computers. These methods show excellent agreement with the exact reference results and can outperform their traditional counterparts when strong electronic correlations become significant. This work sets the stage for future developments of polaritonic quantum chemistry methods suitable for both classical and quantum computers.

arXiv:2106.09842•The Journal of Physical Chemistry Letters•Sep 14, 2021•73 citations
The variational quantum eigensolver (VQE) algorithm combined with the unitary coupled cluster (UCC) ansatz has been developed for the quantum computation of molecular energies and wave functions within the Born–Oppenheimer approximation. Herein, this approach is extended to multicomponent systems to enable the quantum mechanical treatment of more than one type of particle, such as electrons and positrons or electrons and nuclei, without invoking the Born–Oppenheimer approximation. Specifically, we introduce the multicomponent unitary coupled cluster (mcUCC) method combined with the VQE algorithm for the calculation of ground-state energies and wave functions as well as the multicomponent equation-of-motion (mcEOM) method for the calculation of excitation energies. These methods are developed within the nuclear-electronic orbital (NEO) framework and are formulated in the qubit basis to enable implementations on quantum computers. Moreover, these methods are used to calculate the ground-state energy and excitation energies of positronium hydride, where both electrons and the positron are treated quantum mechanically, as well as the H 2 molecule, where both electrons and one proton are treated quantum mechanically. These applications validate the implementation and provide benchmark data for future calculations. The errors due to Trotterization of the mcUCC ansatz are also analyzed. This formalism, as well as the accompanying computer code, will serve as the basis for applications to more complex multicomponent systems, such as simulations of photoinduced nonadiabatic molecular processes, on both classical and quantum computers.
We introduce a multicomponent unitary coupled cluster framework for quantum simulations of molecular systems that incorporate both electronic and nuclear quantum effects beyond the Born-Oppenheimer approximation. Using the nuclear-electronic orbital formalism, we construct mcUCC ansätze for positronium hydride and molecular hydrogen with a quantum proton, and analyze hardware requirements for different excitation truncations. To further reduce resource costs effectively, we employ the local unitary cluster Jastrow ansatz and implement it experimentally on IBM Q's Heron superconducting hardware. With the Physics-Inspired Extrapolation error mitigation protocol, the computed ground-state energies remain within chemical accuracy, consistent with the stated uncertainty level. These results provide the first demonstration of error-mitigated multicomponent correlated simulations on quantum hardware and outline a path toward scalable algorithms unifying electronic and nuclear degrees of freedom.
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.
We present high-precision quantum computing simulations of three-body atoms (He, H−^-−) and molecules (H2+_2^+2+​, HD+^++), the latter being studied beyond the Born-Oppenheimer approximation. The Non-Iterative Disentangled Unitary Coupled Cluster Variational Quantum Eigensolver (NI-DUCC-VQE) [M. Haidar et al., Quantum Sci. Technol. 10, 025031 (2025)] is used. By combining a first-quantized Hamiltonian with a Minimal Complete Pool (MCP) of Lie-algebraic excitations, we construct a compact ansatz with a gradient-independent construction, avoiding costly gradient evaluations and yielding efficient computational scaling with both basis size and electron number. It avoids barren plateaus and enables rapid convergence, achieving energy errors as low as 10−11^{-11}−11 a.u. with state fidelities only limited by arithmetic precision in only a few thousand function evaluations in all four systems. These results make three-body atoms and molecules excellent candidates for benchmarking and testing on current Noisy Intermediate-Scale Quantum (NISQ) devices. Further, our approach can be extended to more complex systems with larger basis sets, taking advantage of the efficient scaling of qubit requirements to study electronic correlations and non-adiabatic effects with high precision. We also demonstrate the applicability of NI-DUCC-VQE for simulating higher-order effects such as relativistic corrections and hyperfine interactions.
Abstract The variational quantum eigensolver (VQE) algorithm, designed to calculate the energy of molecular ground states on near‐term quantum computers, requires specification of symmetries that describe the system, for example, spin state and number of electrons. This opens the possibility of using VQE to obtain excited states as the lowest‐energy solutions of a given set of symmetries. In this paper, the performances of various unitary coupled cluster (UCC) ansätze applied to VQE calculations on excited states are investigated using quantum circuits designed to represent single reference and multireference wavefunctions to calculate energy curves with respect to variations in the molecular geometry. These ansätze include standard Unitary Coupled Cluster Singles and Doubles (UCCSD), as well as modified versions of Unitary Coupled Cluster Generalized Singles and Doubles (UCCGSD) and k‐UpCCGSD, which are engineered to tackle excited states without undesired spin symmetry crossover to lower states during VQE optimization. These studies are carried out on a range of systems, including H2, H3, H4, NH, OH+, CH2, and , covering examples of spin singlet, doublet, and triplet molecular ground states with single and multireference excited states. In most cases, our calculations are in excellent agreement with results from full‐configuration interaction calculations on classical machines, thus showing that the VQE algorithm is capable of calculating the lowest excited state at a certain symmetry, including multireference closed‐ and open‐shell states, by setting appropriate restrictions on the excitations considered in the cluster operator and appropriate constraints in the qubit register encoding the starting mean field state.
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 classical algorithm for UCC based on an exact operator identity for the individual UCC factors. We implement this algorithm for calculations of the H 10 linear chain and the H 2 O molecule with single and double ζ basis sets to provide insights into UCC as a wave-function ansatz. We find that for weakly correlated molecules, the factorized form of the UCC provides similar accuracy to conventional coupled cluster theory (CC); for strongly correlated molecules, where CC often breaks down, UCC significantly outperforms the configuration interaction (CI) ansatz. As a result, the factorized form of the UCC is an accurate, efficient, and reliable electronic structure method in both the weakly and strongly correlated regions. This classical algorithm now allows robust benchmarking of anticipated results from quantum computers and application of coupled-cluster techniques to more strongly correlated molecules.
In the field of polaritonic chemistry, strong light-matter interactions are used to alter chemical reactions inside optical cavities. To understand these processes, the development of reliable theoretical models is essential. While traditional methods have to balance accuracy and system size, new developments in quantum computing offer a path for accurate calculations on currently available quantum devices. Here, we introduce the quantum electrodynamics unitary coupled cluster (QED-UCC) method combined with the Variational Quantum Eigensolver algorithm, as well as the quantum electrodynamics equation-of-motion (QED-EOM) method formulated in the qubit basis that allow accurate calculations of ground-state and excited-state properties of strongly coupled light-matter systems suitable for quantum computers. These methods show excellent agreement with the exact reference results and can outperform their traditional counterparts when strong electronic correlations become significant. This work sets the stage for future developments of polaritonic quantum chemistry methods suitable for both classical and quantum computers.