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

The Journal of Physical Chemistry A

Two Algorithms for Excited-State Quantum Solvers: Theory and Application to EOM-UCCSD

Yongbin Kim, Anna I. Krylov

Published Jul 28, 2023• 23 citations
Publisher

Abstract

Near-term quantum devices promise to revolutionize quantum chemistry, but simulations using the current noisy intermediate-scale quantum (NISQ) devices are not practical due to their high susceptibility to errors. This motivated the design of NISQ algorithms leveraging classical and quantum resources. While several developments have shown promising results for ground-state simulations, extending the algorithms to excited states remains challenging. This paper presents two cost-efficient excited-state algorithms inspired by the classical Davidson algorithm. We implemented the Davidson method into the quantum self-consistent equation-of-motion unitary coupled-cluster (q-sc-EOM-UCC) excited-state method adapted for quantum hardware. The circuit strategies for generating desired excited states are discussed, implemented, and tested. We demonstrate the performance and accuracy of the proposed algorithms (q-sc-EOM-UCC/Davidson and its variational variant) by simulations of H 2, H 4, LiH, and H 2 O molecules. Similar to the classical Davidson scheme, q-sc-EOM-UCC/Davidson algorithms are capable of targeting a small number of excited states of the desired character.

DOI
10.1021/acs.jpca.3c02480

Related papers

Molecular Excited States using Quantum Subspace Methods: Accuracy, Resource Reduction, and Error-Mitigated Hardware Implementation of q-sc-EOM

Srivathsan Poyyapakkam Sundar, Prince Frederick Kwao, Alexey Galda +1 more

Problems in quantum chemical simulations, especially achieving accurate excited-state potential energy surfaces, are among the primary applications to achieve quantum utility. On near-term quantum hardware, variants of the variational quantum eigensolver (VQE) algorithms are the primary choice for chemistry simulation. In this study, a combination of leading ground and excited state quantum algorithms for general excited states, namely, ADAPT-VQE/LUCJ and q-sc-EOM, are utilized to calculate accurate excited state potential energy surfaces in challenging bond-breaking scenarios and compared with the classical scalable EOM-CCSD method. This work investigates avenues toward quantum utility in excited-state quantum chemistry using the q-sc-EOM approach. We assess its accuracy while mitigating major scaling bottlenecks through the Davidson algorithm and basis rotation grouping, reducing the measurement scaling from O(N$^{12}$) to O(N$^{5}$), and implementing the method on quantum hardware with various error mitigation strategies to reduce gate and measurement errors in excited states. The hardware implementation of the q-sc-EOM algorithm, augmented by mitigation of M3 readout error and symmetry projection, produces reasonably accurate excited-state energies with gate noise identified as the predominant source of error. This paves the way for accurate and scalable, generally applicable quantum excited-state methods with potential for quantum utility while identifying critical problems that require advancements.

quant-pharXiv:2604.05380•Apr 8, 2026•0 citations

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Ashutosh Kumar, Ayush Asthana, Conner Masteran +5 more

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.

Journal of Chemical Theory and Computation•Aug 19, 2022•44 citations

Self-consistent Quantum Linear Response with a Polarizable Embedding Environment

Peter Reinholdt, Erik Rosendahl Kjellgren, Karl Michael Ziems +3 more

Quantum computing presents a promising avenue for solving complex problems, particularly in quantum chemistry, where it could accelerate the computation of molecular properties and excited states. This work focuses on computing excitation energies with hybrid quantum-classical algorithms for near-term quantum devices, combining the quantum linear response (qLR) method with a polarizable embedding (PE) environment. We employ the self-consistent operator manifold of quantum linear response (q-sc-LR) on top of a unitary coupled cluster (UCC) wave function in combination with a Davidson solver. The latter removes the need to construct the entire electronic Hessian, improving computational efficiency when going toward larger molecules. We introduce a new superposition-state-based technique to compute Hessian-vector products and show that this approach is more resilient toward noise than our earlier gradient-based approach. We demonstrate the performance of the PE-UCCSD model on systems such as butadiene and para-nitroaniline in water and find that PE-UCCSD delivers comparable accuracy to classical PE-CCSD methods on such simple closed-shell systems. We also explore the challenges posed by hardware noise and propose simple error mitigation techniques to maintain accurate results on noisy quantum computers.

arXiv:2411.03852•The Journal of Physical Chemistry A•Jan 23, 2025•4 citations

Quantum self-consistent equation-of-motion method for computing molecular excitation energies, ionization potentials, and electron affinities on a quantum computer

Ayush Asthana, Ashutosh Kumar, Vibin Abraham +8 more

O, and LiH molecules to test our quantum self-consistent equation-of-motion (q-sc-EOM) method and compare it to other current state-of-the-art methods. q-sc-EOM makes use of self-consistent operators to satisfy the vacuum annihilation condition, a critical property for accurate calculations. It provides real and size-intensive energy differences corresponding to vertical excitation energies, ionization potentials and electron affinities. We also find that q-sc-EOM is more suitable for implementation on NISQ devices as it is expected to be more resilient to noise compared with the currently available methods.

arXiv:2206.10502•Chemical Science•Jan 1, 2023•59 citations

Resolving molecular biexcited states using classical leading-order triples corrections to quantum equation-of-motion UCCSD

Zachary W. Windom, Prince Frederick Kwao, Luke W. Bertels +3 more

Excited electronic states are difficult to accurately and tractably model, particularly those exhibiting multiexcitonic character. Scalable single-reference excited state methods typically fail to describe such states as they suffer from the defects of the mean-field approximation and/or neglect higher-rank excitation operators. To address this, we derive the [T] perturbative correction to the quantum self-consistent equation-of-motion unitary coupled cluster singles and doubles (q-sc-EOM-UCCSD) method. The method captures triple excitation effects through a classical step of second-order in many-body perturbation theory following a q-sc-EOM-UCCSD calculation on a quantum computer. The perturbative correction developed has a $\mathcal{O}(N^7)$ classical computational time scaling. We incorporate [T] into a hybrid compute strategy wherein eigenvectors of the effective Hamiltonian from the hybrid quantum-classical q-sc-EOM-UCCSD algorithm are relayed to a classical computer which performs the postprocessing required for the [T] correction. The benefits of [T] are quantified in a benchmark covering a variety of electronically excited states of the isoelectronic CH$^+$ and BH molecules. We find that the addition of the [T] correction to both trotterized and full operator variants of q-sc-EOM-CCSD can offer dramatic improvements over baseline q-sc-EOM-UCCSD. Our assessment clearly demonstrates the importance of considering higher-rank excitation operators in q-sc-EOM-UCCSD, quantifies the overall success of the [T] correction, and discusses some of its limitations.

arXiv:2610.09260•Oct 7, 2026•0 citations

Molecular excited state VQE simulations with iteratively constructed qubit-excitation-based ans\"atze

Yordan S. Yordanov, C. H. W. Barnes, David R. M. Arvidsson-Shukur

Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ansatze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. In this paper, we introduce the ``excited iterative qubit-excitation-based variational quantum eigensolver'' (e-IQEB-VQE) protocol to calculate molecular excited state energies. The e-IQEB-VQE constructs efficient problem-tailored ansatze by iteratively appending evolutions of qubit excitation operators. The e-IQEB-VQE is an adaptation of the IQEB-VQE protocol, designed to be less dependent on the choice of an initial reference state. We perform classical numerical simulations for LiH and BeH$_2$ to benchmark the performance of the IQEB-VQE. We demonstrate that the e-IQEB-VQE can construct highly accurate ansatze that require at least an order of magnitude fewer $CNOT$s than standard fixed UCC ansatze, such as the UCCSD and the GUCCSD.

•Jun 11, 2021•1 citation
Near-term quantum devices promise to revolutionize quantum chemistry, but simulations using the current noisy intermediate-scale quantum (NISQ) devices are not practical due to their high susceptibility to errors. This motivated the design of NISQ algorithms leveraging classical and quantum resources. While several developments have shown promising results for ground-state simulations, extending the algorithms to excited states remains challenging. This paper presents two cost-efficient excited-state algorithms inspired by the classical Davidson algorithm. We implemented the Davidson method into the quantum self-consistent equation-of-motion unitary coupled-cluster (q-sc-EOM-UCC) excited-state method adapted for quantum hardware. The circuit strategies for generating desired excited states are discussed, implemented, and tested. We demonstrate the performance and accuracy of the proposed algorithms (q-sc-EOM-UCC/Davidson and its variational variant) by simulations of H 2, H 4, LiH, and H 2 O molecules. Similar to the classical Davidson scheme, q-sc-EOM-UCC/Davidson algorithms are capable of targeting a small number of excited states of the desired character.
Problems in quantum chemical simulations, especially achieving accurate excited-state potential energy surfaces, are among the primary applications to achieve quantum utility. On near-term quantum hardware, variants of the variational quantum eigensolver (VQE) algorithms are the primary choice for chemistry simulation. In this study, a combination of leading ground and excited state quantum algorithms for general excited states, namely, ADAPT-VQE/LUCJ and q-sc-EOM, are utilized to calculate accurate excited state potential energy surfaces in challenging bond-breaking scenarios and compared with the classical scalable EOM-CCSD method. This work investigates avenues toward quantum utility in excited-state quantum chemistry using the q-sc-EOM approach. We assess its accuracy while mitigating major scaling bottlenecks through the Davidson algorithm and basis rotation grouping, reducing the measurement scaling from O(N12^{12}12) to O(N5^{5}5), and implementing the method on quantum hardware with various error mitigation strategies to reduce gate and measurement errors in excited states. The hardware implementation of the q-sc-EOM algorithm, augmented by mitigation of M3 readout error and symmetry projection, produces reasonably accurate excited-state energies with gate noise identified as the predominant source of error. This paves the way for accurate and scalable, generally applicable quantum excited-state methods with potential for quantum utility while identifying critical problems that require advancements.
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.
Quantum computing presents a promising avenue for solving complex problems, particularly in quantum chemistry, where it could accelerate the computation of molecular properties and excited states. This work focuses on computing excitation energies with hybrid quantum-classical algorithms for near-term quantum devices, combining the quantum linear response (qLR) method with a polarizable embedding (PE) environment. We employ the self-consistent operator manifold of quantum linear response (q-sc-LR) on top of a unitary coupled cluster (UCC) wave function in combination with a Davidson solver. The latter removes the need to construct the entire electronic Hessian, improving computational efficiency when going toward larger molecules. We introduce a new superposition-state-based technique to compute Hessian-vector products and show that this approach is more resilient toward noise than our earlier gradient-based approach. We demonstrate the performance of the PE-UCCSD model on systems such as butadiene and para-nitroaniline in water and find that PE-UCCSD delivers comparable accuracy to classical PE-CCSD methods on such simple closed-shell systems. We also explore the challenges posed by hardware noise and propose simple error mitigation techniques to maintain accurate results on noisy quantum computers.
O, and LiH molecules to test our quantum self-consistent equation-of-motion (q-sc-EOM) method and compare it to other current state-of-the-art methods. q-sc-EOM makes use of self-consistent operators to satisfy the vacuum annihilation condition, a critical property for accurate calculations. It provides real and size-intensive energy differences corresponding to vertical excitation energies, ionization potentials and electron affinities. We also find that q-sc-EOM is more suitable for implementation on NISQ devices as it is expected to be more resilient to noise compared with the currently available methods.
Excited electronic states are difficult to accurately and tractably model, particularly those exhibiting multiexcitonic character. Scalable single-reference excited state methods typically fail to describe such states as they suffer from the defects of the mean-field approximation and/or neglect higher-rank excitation operators. To address this, we derive the [T] perturbative correction to the quantum self-consistent equation-of-motion unitary coupled cluster singles and doubles (q-sc-EOM-UCCSD) method. The method captures triple excitation effects through a classical step of second-order in many-body perturbation theory following a q-sc-EOM-UCCSD calculation on a quantum computer. The perturbative correction developed has a O(N7)\mathcal{O}(N^7)O(N7) classical computational time scaling. We incorporate [T] into a hybrid compute strategy wherein eigenvectors of the effective Hamiltonian from the hybrid quantum-classical q-sc-EOM-UCCSD algorithm are relayed to a classical computer which performs the postprocessing required for the [T] correction. The benefits of [T] are quantified in a benchmark covering a variety of electronically excited states of the isoelectronic CH+^++ and BH molecules. We find that the addition of the [T] correction to both trotterized and full operator variants of q-sc-EOM-CCSD can offer dramatic improvements over baseline q-sc-EOM-UCCSD. Our assessment clearly demonstrates the importance of considering higher-rank excitation operators in q-sc-EOM-UCCSD, quantifies the overall success of the [T] correction, and discusses some of its limitations.
Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ansatze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. In this paper, we introduce the ``excited iterative qubit-excitation-based variational quantum eigensolver'' (e-IQEB-VQE) protocol to calculate molecular excited state energies. The e-IQEB-VQE constructs efficient problem-tailored ansatze by iteratively appending evolutions of qubit excitation operators. The e-IQEB-VQE is an adaptation of the IQEB-VQE protocol, designed to be less dependent on the choice of an initial reference state. We perform classical numerical simulations for LiH and BeH2_22​ to benchmark the performance of the IQEB-VQE. We demonstrate that the e-IQEB-VQE can construct highly accurate ansatze that require at least an order of magnitude fewer CNOTCNOTCNOTs than standard fixed UCC ansatze, such as the UCCSD and the GUCCSD.