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

Physical Chemistry Chemical Physics

Solving complex eigenvalue problems on a quantum annealer with applications to quantum scattering resonances

Alexander Teplukhin, Brian K. Kendrick, Dmitri Babikov

Published Jan 1, 2020• 28 citations
Publisher

Abstract

Quantum computing is a new and rapidly evolving paradigm for solving chemistry problems. In previous work, we developed the Quantum Annealer Eigensolver (QAE) and applied it to the calculation of the vibrational spectrum of a molecule on the D-Wave quantum annealer. However, the original QAE methodology was applicable to real symmetric matrices only. For many physics and chemistry problems, the diagonalization of complex matrices is required. For example, the calculation of quantum scattering resonances can be formulated as a complex eigenvalue problem where the real part of the eigenvalue is the resonance energy and the imaginary part is proportional to the resonance width. In the present work, we generalize the QAE to treat complex matrices: first complex Hermitian matrices and then complex symmetric matrices. These generalizations are then used to compute a quantum scattering resonance state in a 1D model potential for O + O collisions. These calculations are performed using both a software (classical) annealer and hardware annealer (the D-Wave 2000Q). The results of the complex QAE are also benchmarked against a standard linear algebra library (LAPACK). This work presents the first numerical solution of a complex eigenvalue problem of any kind on a quantum annealer, and it is the first treatment of a quantum scattering resonance on any quantum device.

DOI
10.1039/d0cp04272b

Related papers

Electronic structure with direct diagonalization on a D-wave quantum annealer

Alexander Teplukhin, Brian K. Kendrick, Sergei Tretiak +1 more

Quantum chemistry is regarded to be one of the first disciplines that will be revolutionized by quantum computing. Although universal quantum computers of practical scale may be years away, various approaches are currently being pursued to solve quantum chemistry problems on near-term gate-based quantum computers and quantum annealers by developing the appropriate algorithm and software base. This work implements the general Quantum Annealer Eigensolver (QAE) algorithm to solve the molecular electronic Hamiltonian eigenvalue-eigenvector problem on a D-Wave 2000Q quantum annealer. The approach is based on the matrix formulation, efficiently uses qubit resources based on a power-of-two encoding scheme and is hardware-dominant relying on only one classically optimized parameter. We demonstrate the use of D-Wave hardware for obtaining ground and excited electronic states across a variety of small molecular systems. The approach can be adapted for use by a vast majority of electronic structure methods currently implemented in conventional quantum-chemical packages. The results of this work will encourage further development of software such as qbsolv which has promising applications in emerging quantum information processing hardware and has expectation to address large and complex optimization problems intractable for classical computers.

arXiv:2009.01373•Scientific Reports•Nov 27, 2020•36 citations

Calculation of Molecular Vibrational Spectra on a Quantum Annealer

Alexander Teplukhin, Brian K. Kendrick, Dmitri Babikov

Abstract Until recently molecular energy calculations using quantum computing hardware have been limited to gate-based quantum computers. In this paper, a new methodology is presented to calculate the vibrational spectrum of a molecule on a quantum annealer. The key idea of the method is a mapping of the ground state variational problem onto an Ising or quadratic unconstrained binary optimization (QUBO) problem by expressing the expansion coefficients using spins or qubits. The algorithm is general and represents a new revolutionary approach for solving the real symmetric eigenvalue problem on a quantum annealer. The method is applied to two chemically important molecules: O2 (oxygen) and O3 (ozone). The lowest two vibrational states of these molecules are computed using both a hardware quantum annealer and a software based classical annealer. Extension of the algorithm to higher dimensions is explicitly demonstrated for an N-dimensional harmonic oscillator (N ≤ 5). The algorithm scales exponentially with dimensionality if a direct product basis is used but will exhibit polynomial scaling for a nondirect product basis.

arXiv:1812.05211•Journal of Chemical Theory and Computation•Jul 17, 2019•38 citations

Computing molecular excited states on a D-Wave quantum annealer

Alexander Teplukhin, Brian K. Kendrick, Susan M. Mniszewski +6 more

The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree-Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm-Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example [Formula: see text] molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for [Formula: see text] as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

arXiv:2107.00162•Scientific Reports•Sep 22, 2021•34 citations

Solution of the Equation-of-Motion Phonon Method eigenvalue problems on the D-Wave quantum annealer

C. De Lucia, A. Martone, F. A. D'Aniello +7 more

The solution of large-scale eigenvalue problems is crucial in nuclear many-body theory, where Hamiltonian matrices often reach extremely large dimensions. Quantum computing opens new perspectives for addressing such demanding problems. Although the Quantum Phase Estimation algorithm offers, in principle, a systematic route to matrix diagonalization, its practical deployment demands levels of coherence and error correction that current quantum hardware cannot yet support. A viable near-term strategy is instead to exploit quantum annealing, which enables the recasting of eigenvalue problems into quadratic unconstrained binary optimization formulations that can be addressed by existing annealing-based processors. Here, we propose a hybrid quantum-classical algorithm that combines quantum annealing and classical deflation to iteratively extract the full eigenspectrum of both standard and generalized eigenvalue problems. We benchmark this method on eigenvalue problems arising from the Equation of Motion Phonon Method performing calculations on real quantum hardware. Our approach illustrates the capabilities and limitations of near-term quantum devices in addressing nuclear eigenvalue problems.

nucl-thCross-listedarXiv:2606.07035•Jun 8, 2026•0 citations

Quantum chemistry on quantum annealers

Scott N. Genin, Ilya G. Ryabinkin, Artur F. Izmaylov

Quantum chemistry calculations for small molecules on quantum hardware have been demonstrated to date only on universal-gate quantum computers, not quantum annealers. The latter devices are limited to finding the lowest eigenstate of the Ising Hamiltonian whereas the electronic Hamiltonian could not be mapped to the Ising form without exponential growth of the Ising Hamiltonian with the size of the system [J. Phys. Chem. B 122, 3384 (2018)]. Here we propose a novel mixed discrete-continuous optimization algorithm, which finds the lowest eigenstate of the qubit coupled cluster (QCC) method using a quantum annealer for solving a discrete part of the problem. The QCC method is a potentially exact approach for constructing the electronic wave function in the qubit space. Therefore, our methodology allows for systematically improvable quantum chemistry calculations using quantum annealears. We illustrate capabilities of our approach by calculating QCC ground electronic states for the LiH, H$_2$O, and C$_6$H$_6$ molecules. C$_6$H$_6$ calculations involve 36 qubits and are the largest quantum chemistry calculations made on a quantum annealer (the D-Wave 2000Q system) to date. Our findings opens up a new perspective for use quantum annealers in high-throughput material discovery.

arXiv:1901.04715•Jan 15, 2019•20 citations

Quantum computing for atomic and molecular resonances

Teng Bian, Sabre Kais

The complex-scaling method can be used to calculate molecular resonances within the Born–Oppenheimer approximation, assuming that the electronic coordinates are dilated independently of the nuclear coordinates. With this method, one will calculate the complex energy of a non-Hermitian Hamiltonian, whose real part is associated with the resonance position and imaginary part is the inverse of the lifetime. In this study, we propose techniques to simulate resonances on a quantum computer. First, we transformed the scaled molecular Hamiltonian to second quantization and then used the Jordan–Wigner transformation to transform the scaled Hamiltonian to the qubit space. To obtain the complex eigenvalues, we introduce the direct measurement method, which is applied to obtain the resonances of a simple one-dimensional model potential that exhibits pre-dissociating resonances analogous to those found in diatomic molecules. Finally, we applied the method to simulate the resonances of the H2− molecule. The numerical results from the IBM Qiskit simulators and IBM quantum computers verify our techniques.

arXiv:2011.13999•The Journal of Chemical Physics•May 19, 2021•11 citations
Quantum computing is a new and rapidly evolving paradigm for solving chemistry problems. In previous work, we developed the Quantum Annealer Eigensolver (QAE) and applied it to the calculation of the vibrational spectrum of a molecule on the D-Wave quantum annealer. However, the original QAE methodology was applicable to real symmetric matrices only. For many physics and chemistry problems, the diagonalization of complex matrices is required. For example, the calculation of quantum scattering resonances can be formulated as a complex eigenvalue problem where the real part of the eigenvalue is the resonance energy and the imaginary part is proportional to the resonance width. In the present work, we generalize the QAE to treat complex matrices: first complex Hermitian matrices and then complex symmetric matrices. These generalizations are then used to compute a quantum scattering resonance state in a 1D model potential for O + O collisions. These calculations are performed using both a software (classical) annealer and hardware annealer (the D-Wave 2000Q). The results of the complex QAE are also benchmarked against a standard linear algebra library (LAPACK). This work presents the first numerical solution of a complex eigenvalue problem of any kind on a quantum annealer, and it is the first treatment of a quantum scattering resonance on any quantum device.
Quantum chemistry is regarded to be one of the first disciplines that will be revolutionized by quantum computing. Although universal quantum computers of practical scale may be years away, various approaches are currently being pursued to solve quantum chemistry problems on near-term gate-based quantum computers and quantum annealers by developing the appropriate algorithm and software base. This work implements the general Quantum Annealer Eigensolver (QAE) algorithm to solve the molecular electronic Hamiltonian eigenvalue-eigenvector problem on a D-Wave 2000Q quantum annealer. The approach is based on the matrix formulation, efficiently uses qubit resources based on a power-of-two encoding scheme and is hardware-dominant relying on only one classically optimized parameter. We demonstrate the use of D-Wave hardware for obtaining ground and excited electronic states across a variety of small molecular systems. The approach can be adapted for use by a vast majority of electronic structure methods currently implemented in conventional quantum-chemical packages. The results of this work will encourage further development of software such as qbsolv which has promising applications in emerging quantum information processing hardware and has expectation to address large and complex optimization problems intractable for classical computers.
Abstract Until recently molecular energy calculations using quantum computing hardware have been limited to gate-based quantum computers. In this paper, a new methodology is presented to calculate the vibrational spectrum of a molecule on a quantum annealer. The key idea of the method is a mapping of the ground state variational problem onto an Ising or quadratic unconstrained binary optimization (QUBO) problem by expressing the expansion coefficients using spins or qubits. The algorithm is general and represents a new revolutionary approach for solving the real symmetric eigenvalue problem on a quantum annealer. The method is applied to two chemically important molecules: O2 (oxygen) and O3 (ozone). The lowest two vibrational states of these molecules are computed using both a hardware quantum annealer and a software based classical annealer. Extension of the algorithm to higher dimensions is explicitly demonstrated for an N-dimensional harmonic oscillator (N ≤ 5). The algorithm scales exponentially with dimensionality if a direct product basis is used but will exhibit polynomial scaling for a nondirect product basis.
The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree-Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm-Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example [Formula: see text] molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for [Formula: see text] as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.
The solution of large-scale eigenvalue problems is crucial in nuclear many-body theory, where Hamiltonian matrices often reach extremely large dimensions. Quantum computing opens new perspectives for addressing such demanding problems. Although the Quantum Phase Estimation algorithm offers, in principle, a systematic route to matrix diagonalization, its practical deployment demands levels of coherence and error correction that current quantum hardware cannot yet support. A viable near-term strategy is instead to exploit quantum annealing, which enables the recasting of eigenvalue problems into quadratic unconstrained binary optimization formulations that can be addressed by existing annealing-based processors. Here, we propose a hybrid quantum-classical algorithm that combines quantum annealing and classical deflation to iteratively extract the full eigenspectrum of both standard and generalized eigenvalue problems. We benchmark this method on eigenvalue problems arising from the Equation of Motion Phonon Method performing calculations on real quantum hardware. Our approach illustrates the capabilities and limitations of near-term quantum devices in addressing nuclear eigenvalue problems.
Quantum chemistry calculations for small molecules on quantum hardware have been demonstrated to date only on universal-gate quantum computers, not quantum annealers. The latter devices are limited to finding the lowest eigenstate of the Ising Hamiltonian whereas the electronic Hamiltonian could not be mapped to the Ising form without exponential growth of the Ising Hamiltonian with the size of the system [J. Phys. Chem. B 122, 3384 (2018)]. Here we propose a novel mixed discrete-continuous optimization algorithm, which finds the lowest eigenstate of the qubit coupled cluster (QCC) method using a quantum annealer for solving a discrete part of the problem. The QCC method is a potentially exact approach for constructing the electronic wave function in the qubit space. Therefore, our methodology allows for systematically improvable quantum chemistry calculations using quantum annealears. We illustrate capabilities of our approach by calculating QCC ground electronic states for the LiH, H2_22​O, and C6_66​H6_66​ molecules. C6_66​H6_66​ calculations involve 36 qubits and are the largest quantum chemistry calculations made on a quantum annealer (the D-Wave 2000Q system) to date. Our findings opens up a new perspective for use quantum annealers in high-throughput material discovery.
The complex-scaling method can be used to calculate molecular resonances within the Born–Oppenheimer approximation, assuming that the electronic coordinates are dilated independently of the nuclear coordinates. With this method, one will calculate the complex energy of a non-Hermitian Hamiltonian, whose real part is associated with the resonance position and imaginary part is the inverse of the lifetime. In this study, we propose techniques to simulate resonances on a quantum computer. First, we transformed the scaled molecular Hamiltonian to second quantization and then used the Jordan–Wigner transformation to transform the scaled Hamiltonian to the qubit space. To obtain the complex eigenvalues, we introduce the direct measurement method, which is applied to obtain the resonances of a simple one-dimensional model potential that exhibits pre-dissociating resonances analogous to those found in diatomic molecules. Finally, we applied the method to simulate the resonances of the H2− molecule. The numerical results from the IBM Qiskit simulators and IBM quantum computers verify our techniques.