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

arXiv:2110.12096 · Scientific Reports

Molecular dynamics on quantum annealers

Igor Gaidai, Dmitri Babikov, Alexander Teplukhin, Brian K. Kendrick, Susan M. Mniszewski, Yu Zhang, Sergei Tretiak, Pavel A. Dub

Published Oct 7, 2022• 13 citations
arXivHTMLPDF

Abstract

In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.

DOI
10.1038/s41598-022-21163-x

Related papers

Molecular dynamics with time dependent quantum Monte Carlo

Ivan P. Christov

In this paper we propose an ab initio method to solve quantum many-body problems of molecular dynamics where both electronic and nuclear degrees are represented by ensembles of trajectories and guiding waves in physical space. Both electrons and nuclei can be treated quantum mechanically where the guiding waves obey a set of coupled Schrodinger equations (quantum-quantum description) or, alternatively, coupled Schrodinger-Newtonian equations are solved for the quantum-classical approximation. The method takes into account local and nonlocal quantum correlation effects in a self-consistent manner. The general formalism is applied to one- and two-dimensional hydrogen molecules subjected to a strong ultrashort optical pulse. Comparison is made with the results from the "exact" Ehrenfest molecular dynamics for the molecular ionization and for the evolution of the internuclear distance as the molecule dissociates.

arXiv:2501.16132•The Journal of Chemical Physics•Dec 3, 2008•12 citations

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

Alexander Teplukhin, Brian K. Kendrick, Dmitri Babikov

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.

Physical Chemistry Chemical Physics•Jan 1, 2020•28 citations

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

Estimation Of The Quantum Effects Of Nuclei in Large Molecular Systems

Bing Gu

Chemical dynamics, in principle, should be understood by solving the time-dependent Schrödinger equation for a molecular system, describing motion of the nuclei and electrons. However, the computational efforts to solve this partial second-order differential equation scales exponentially with the system size, which prevents us from getting exact numerical solutions for systems larger than 4-5 atoms. Thus, approximations simplifying the picture are necessary. The so-called Born-Oppenheimer approximation, separating motion of the electrons and nuclei is the central one: solution to the electronic Schrödinger equation defines the potential energy surface on which the nuclear motion unfolds, and there are standard quantum chemistry software packages for solving the electronic Schrödinger equation. For the nuclear Schrödinger equation, however, there are no widely applicable quantum-mechanical approaches, and most simulations are performed using classical Newtonian mechanics which is often adequate due to large nuclear masses. However, the nuclear quantum effects are significant for chemical processes involving light nuclei at low energies, and including these effects into simulation, even approximately, is highly desirable. In this dissertation, an approximate methodology of including quantum-mechanical effects within the quantum trajectory or the de Broglie-Bohm formulation of the Schrödinger equations is developed. Use of the trajectory framework makes the approach scalable to hundreds of degrees of freedom. The methodology is applied to study high-dimensional systems (solid He4 and others) relevant to chemistry.

•Jan 1, 2016•0 citations

A Quantum Annealing Approach to Fluid Dynamics Problems Solving Navier-Stokes Equations

Juan De Dios Rodríguez, Alejandro Pujante Pérez, Eduardo Illueca Fernández +1 more

Partial differential equations (PDEs) present important computational challenges in several scientific disciplines. Classical numerical methods, while effective, often encounter scalability issues and high computational costs, particularly for complex problems. Quantum computing seems to address these challenges, offering a potential exponential speedup over classical methods. The presented work studies the application of quantum computing, specifically Quantum Annealing, to the computational complexities inherent in solving the Navier-Stokes equations within Computational Fluid Dynamics (CFD). Our study examines the quantum annealing methodology, conducting two different experiments and analyzing the error between the simulation and the analytical solution. The results obtained demonstrate that quantum annealing paradigm can solve this problems and the errors obtained are acceptable in comparison with other methods in the state of the art. By combining theoretical insights with practical implementation, this research clarifies the potential of quantum annealing, providing a comprehensive overview of the simulation framework, analyzing the obtained results, and suggesting future research avenues.

International Conference on Smart and Sustainable Technologies (SpliTech)•Jun 25, 2024•1 citation
In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.
In this paper we propose an ab initio method to solve quantum many-body problems of molecular dynamics where both electronic and nuclear degrees are represented by ensembles of trajectories and guiding waves in physical space. Both electrons and nuclei can be treated quantum mechanically where the guiding waves obey a set of coupled Schrodinger equations (quantum-quantum description) or, alternatively, coupled Schrodinger-Newtonian equations are solved for the quantum-classical approximation. The method takes into account local and nonlocal quantum correlation effects in a self-consistent manner. The general formalism is applied to one- and two-dimensional hydrogen molecules subjected to a strong ultrashort optical pulse. Comparison is made with the results from the "exact" Ehrenfest molecular dynamics for the molecular ionization and for the evolution of the internuclear distance as the molecule dissociates.
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.
Chemical dynamics, in principle, should be understood by solving the time-dependent Schrödinger equation for a molecular system, describing motion of the nuclei and electrons. However, the computational efforts to solve this partial second-order differential equation scales exponentially with the system size, which prevents us from getting exact numerical solutions for systems larger than 4-5 atoms. Thus, approximations simplifying the picture are necessary. The so-called Born-Oppenheimer approximation, separating motion of the electrons and nuclei is the central one: solution to the electronic Schrödinger equation defines the potential energy surface on which the nuclear motion unfolds, and there are standard quantum chemistry software packages for solving the electronic Schrödinger equation. For the nuclear Schrödinger equation, however, there are no widely applicable quantum-mechanical approaches, and most simulations are performed using classical Newtonian mechanics which is often adequate due to large nuclear masses. However, the nuclear quantum effects are significant for chemical processes involving light nuclei at low energies, and including these effects into simulation, even approximately, is highly desirable. In this dissertation, an approximate methodology of including quantum-mechanical effects within the quantum trajectory or the de Broglie-Bohm formulation of the Schrödinger equations is developed. Use of the trajectory framework makes the approach scalable to hundreds of degrees of freedom. The methodology is applied to study high-dimensional systems (solid He4 and others) relevant to chemistry.
Partial differential equations (PDEs) present important computational challenges in several scientific disciplines. Classical numerical methods, while effective, often encounter scalability issues and high computational costs, particularly for complex problems. Quantum computing seems to address these challenges, offering a potential exponential speedup over classical methods. The presented work studies the application of quantum computing, specifically Quantum Annealing, to the computational complexities inherent in solving the Navier-Stokes equations within Computational Fluid Dynamics (CFD). Our study examines the quantum annealing methodology, conducting two different experiments and analyzing the error between the simulation and the analytical solution. The results obtained demonstrate that quantum annealing paradigm can solve this problems and the errors obtained are acceptable in comparison with other methods in the state of the art. By combining theoretical insights with practical implementation, this research clarifies the potential of quantum annealing, providing a comprehensive overview of the simulation framework, analyzing the obtained results, and suggesting future research avenues.