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

Chemical Science

SHARC meets TEQUILA: mixed quantum-classical dynamics on a quantum computer using a hybrid quantum-classical algorithm

Eduarda Sangiogo Gil, Markus Oppel, Jakob S. Kottmann, Leticia González

Published Nov 28, 2024• 10 citations
Publisher

Abstract

Recent developments in quantum computing are highly promising, particularly in the realm of quantum chemistry. Due to the noisy nature of currently available quantum hardware, hybrid quantum-classical algorithms have emerged as a reliable option for near-term simulations. Mixed quantum-classical dynamics methods effectively capture nonadiabatic effects by integrating classical nuclear dynamics with quantum chemical computations of the electronic properties. However, these methods face challenges due to the high computational cost of the quantum chemistry part. To mitigate the computational demand, we propose a method where the required electronic properties are computed through a hybrid quantum-classical approach on combination of classical and quantum hardware. This framework employs the variational quantum eigensolver and variational quantum deflation algorithms to obtain ground and excited state energies, gradients, nonadiabatic coupling vectors, and transition dipole moments. These quantities are used to propagate the nonadiabatic molecular dynamics using the Tully’s fewest switches surface hopping method, although the implementation is also compatible with other molecular dynamics approaches. The approach, implemented by integrating the molecular dynamics program package SHARC with the TEQUILA quantum computing framework, is validated by studying the cis-trans photoisomerization of methanimine and the electronic relaxation of ethylene. The results show qualitatively accurate molecular dynamics that align with experimental findings and other computational studies. This work is expected to mark a significant step towards achieving a "quantum advantage" for realistic chemical simulations.

DOI
10.1039/d4sc04987j

Related papers

Nonadiabatic Molecular Dynamics with Quantum Algorithms

Silvia Riera-Villapún, Jaime Scharfhausen-Curiel, J. L. Sánchez Toural +2 more

Investigating ultrafast electron-nuclear dynamics in molecules requires robust theoretical methods that offer predictive capabilities beyond interpretation purposes. However, the vast number of electronic and nuclear degrees of freedom often prevents an accurate description of chemical systems on classical computers. To address these challenges, hybrid approaches combining quantum and classical computing can leverage the strengths of both to enhance the computational efficiency and fully exploit hardware capabilities. In this work, we perform excited-state molecular dynamics simulations of the formaldimine molecule by combining the surface hopping nonadiabatic molecular dynamics technique with quantum algorithms. We replace conventional electronic energy and overlap calculations with the variational quantum deflation (VQD) algorithm and the Hadamard test, respectively, employing a configuration state function (CSF) spin-restricted ansatz. Finally, we evaluate the accuracy and feasibility of three different gradient calculation approaches integrated with the quantum computing components, demonstrating the advantages of an analytical gradient method over finite difference approaches for hybrid quantum-classical integration.

The Journal of Physical Chemistry A•Aug 19, 2026•0 citations

Nonadiabatic molecular dynamics on quantum computers: challenges and opportunities

Eduarda Sangiogo Gil, Leticia González

Abstract In this Perspective, we discuss how quantum computers may advance the simulation of nonadiabatic molecular dynamics, a framework central to describing excited-state processes in photochemistry, biology, and materials science. Classical approaches span from exponentially scaling full quantum dynamics to more approximate mixed quantum–classical techniques such as surface hopping and Ehrenfest dynamics. Hybrid quantum–classical algorithms – particularly those based on the variational quantum eigensolver – offer a transformative alternative by providing access to the key electronic properties needed to drive nonadiabatic molecular dynamics simulations, including energies, gradients, and nonadiabatic couplings. We examine recent proof-of-principle quantum simulations of reduced model systems which, despite being restricted to small molecules and limited active spaces due to constraints of qubit number and device noise, already showcase the potential of quantum devices to capture phenomena such as conical intersections and ultrafast relaxation. Although practical applications are not yet feasible in the present noisy intermediate-scale quantum era, these efforts underline the conceptual and methodological advances of quantum algorithms paving the way for large-scale quantum simulations of nonadiabatic processes. Framed within the 2025 International Year of Quantum Science and Technology, such progress exemplifies how quantum computing may open new horizons for chemistry and beyond.

Pure and Applied Chemistry•Nov 1, 2025•2 citations

Molecular Dynamics with Nuclear Effects on Quantum Computers

Lukas Haßfurth, Juliane Heitkämper, Elias Walter +1 more

Nuclear quantum effects are critical for describing proton transfer and hydrogen bonding, but their incorporation into quantum chemistry calculations is often computationally prohibitive on classical hardware. A promising alternative are quantum computers due to their linear scaling in the space requirements with system size. We introduce a novel hybrid quantum-classical algorithm for ab-initio molecular dynamics that incorporates nuclear quantum effects via the nuclear-electronic orbital method. The proposed algorithm evaluates ground state energies and forces on the quantum computer using a variational quantum eigensolver, while the molecular geometries are updated classically. We validate our approach through simulations of $\text{H}_2$, $\text{H}_2\text{O}$ and the Zundel ion $\text{H}_5\text{O}_2^+$, comparing the simulated vibrational spectra with experimental data. Upon inclusion of nuclear quantum effects, the proton shuttling movement in the Zundel ion becomes effectively barrierless, and errors in the simulated frequencies improve significantly. Employing compact hardware-efficient ansatz circuits we achieve results comparable to the more accurate UCCSD ansatzes, which hints towards the feasibility of executing our algorithm on near-term quantum devices.

quant-pharXiv:2610.01549•Oct 2, 2026•0 citations

Hybrid Quantum Algorithms for Computational Chemistry: Application to the Pyridine-Li ion Complex

Fatemeh Ghasemi, Yousung Kang, Yukio Kawashima +1 more

Accurately capturing electron correlation in large-scale molecular systems remains one of the foremost challenges in quantum chemistry and a primary driver for the development of quantum algorithms. Classical configuration-interaction methods, while rigorous, suffer from exponential scaling, rendering them impractical for large or strongly correlated systems. Overcoming this limitation is central to realizing the promise of quantum computing in chemistry. Here, we investigate the pyridine-Li ion complex using three quantum algorithms: the variational quantum eigensolver (VQE), the subspace quantum diagonalization (SQD) method, and the recently introduced handover iterative VQE (HI-VQE). Our results demonstrate how new generations of hybrid quantum-classical frameworks overcome the scalability and noise sensitivity that constrain conventional VQE approaches. SQD and HI-VQE achieve ground-state energy calculations for problem sizes inaccessible to classical computation, marking a clear advance toward quantum advantage. In particular, HI-VQE enables calculations within active spaces as large as (24e,22o), requiring 44 qubits-well beyond the reach of classical CASCI and VQE. This capability provides a systematic pathway for incorporating increasing numbers of electrons into quantum treatment, thereby approaching exact molecular energies. Importantly, both SQD and HI-VQE exhibit robustness against hardware noise, a critical improvement over earlier approaches. By enabling quantum simulations of molecular systems previously deemed intractable, SQD and HI-VQE offer a realistic route toward practical quantum advantage in computational chemistry. The comparison between HI-VQE and SQD shows that optimizing circuit parameters is crucial for accurate simulation.

physics.chem-phCross-listedarXiv:2601.10002•Jan 16, 2026•0 citations

Considerations for evaluating thermodynamic properties with hybrid quantum-classical computing work flows

Spencer T. Stober, Stuart M. Harwood, Dimitar Trenev +3 more

Quantum chemistry applications on quantum computers currently rely heavily on the variational quantum eigensolver (VQE) algorithm. This hybrid quantum-classical algorithm aims at finding ground-state solutions of molecular systems based on the variational principle. VQE calculations can be systematically implemented for perturbations to each molecular degree of freedom, generating a Born-Oppenheimer potential-energy surface (PES) for the molecule. The PES can then be used to derive thermodynamic properties, which are often desirable for applications in chemical engineering and materials design. It is clear from this process that quantum chemistry applications contain a substantial classical computing component in addition to steps that can be performed using a quantum computer. In order to design efficient work flows that take full advantage of each hardware type, it is critical to consider the entire process so that the high-accuracy electronic energies possible from quantum computing are not squandered in the process of calculating thermodynamic properties. We present a summary of the hybrid quantum-classical work flow to compute thermodynamic properties. This work flow contains many options that can significantly affect the efficiency and the accuracy of the results, including classical optimizer attributes, number of ansatz repetitions, and how the vibrational Schr\"odinger equation is solved to determine vibrational modes. We also analyze the effects of these options by employing robust statistics along with simulations and experiments on actual quantum hardware. We show that, through careful selection of work flow options, nearly order-of-magnitude increases in accuracy are possible at equivalent computing time.

arXiv:2003.02303•Physical Review A•Jan 25, 2022•15 citations

Adiabatic-Inspired Hybrid Quantum-Classical Methods for Molecular Ground State Preparation

Sean Thrasher, Ioannis Kolotouros, Julien Michel +1 more

Quantum computing promises to efficiently and accurately solve many important problems in quantum chemistry which elude classical solvers, such as the electronic structure problem of highly correlated materials. Two leading methods in solving the ground state problem are the Variational Quantum Eigensolver (VQE) and Adiabatic Quantum Computing (AQC) algorithms. VQE often struggles with convergence due to the energy landscape being highly non-convex and the existence of barren plateaux, and implementing AQC is beyond the capabilities of current quantum devices as it requires deep circuits. Adiabatically-inspired algorithms aim to fill this gap. In this paper, we first present a unifying framework for these algorithms and then benchmark the following methods: the Adiabatically Assisted VQE (AAVQE) (Garcia-Saez and Latorre (2018)), the Variational Adiabatic Quantum Computing (VAQC) (Harwood et al (2022)), and the Adiabatic Quantum Computing with Parametrized Quantum Circuits (AQC-PQC) (Kolotouros et al (2025)) algorithms. Second, we introduce a novel hybrid approach termed G-AQC-PQC, which generalizes the AQC-PQC method, and combines adiabatic-inspired initialization with the low-memory BFGS optimizer, reducing the quantum computational cost of the method. Third, we compare the accuracy of the methods for chemistry applications using the beryllium hydride molecule (BeH$_2$). We compare the approaches across a number of different choices (ansätze types, depth, discretization steps, initial Hamiltonian, adiabatic schedules and method used). Our results show that the G-AQC-PQC outperforms conventional VQE. We further discuss limitations such as the zero-gradient problem and identify regimes where adiabatically-inspired methods offer a tangible advantage for near-term quantum chemistry applications.

quant-pharXiv:2512.14449•Dec 17, 2025•0 citations
Recent developments in quantum computing are highly promising, particularly in the realm of quantum chemistry. Due to the noisy nature of currently available quantum hardware, hybrid quantum-classical algorithms have emerged as a reliable option for near-term simulations. Mixed quantum-classical dynamics methods effectively capture nonadiabatic effects by integrating classical nuclear dynamics with quantum chemical computations of the electronic properties. However, these methods face challenges due to the high computational cost of the quantum chemistry part. To mitigate the computational demand, we propose a method where the required electronic properties are computed through a hybrid quantum-classical approach on combination of classical and quantum hardware. This framework employs the variational quantum eigensolver and variational quantum deflation algorithms to obtain ground and excited state energies, gradients, nonadiabatic coupling vectors, and transition dipole moments. These quantities are used to propagate the nonadiabatic molecular dynamics using the Tully’s fewest switches surface hopping method, although the implementation is also compatible with other molecular dynamics approaches. The approach, implemented by integrating the molecular dynamics program package SHARC with the TEQUILA quantum computing framework, is validated by studying the cis-trans photoisomerization of methanimine and the electronic relaxation of ethylene. The results show qualitatively accurate molecular dynamics that align with experimental findings and other computational studies. This work is expected to mark a significant step towards achieving a "quantum advantage" for realistic chemical simulations.
Investigating ultrafast electron-nuclear dynamics in molecules requires robust theoretical methods that offer predictive capabilities beyond interpretation purposes. However, the vast number of electronic and nuclear degrees of freedom often prevents an accurate description of chemical systems on classical computers. To address these challenges, hybrid approaches combining quantum and classical computing can leverage the strengths of both to enhance the computational efficiency and fully exploit hardware capabilities. In this work, we perform excited-state molecular dynamics simulations of the formaldimine molecule by combining the surface hopping nonadiabatic molecular dynamics technique with quantum algorithms. We replace conventional electronic energy and overlap calculations with the variational quantum deflation (VQD) algorithm and the Hadamard test, respectively, employing a configuration state function (CSF) spin-restricted ansatz. Finally, we evaluate the accuracy and feasibility of three different gradient calculation approaches integrated with the quantum computing components, demonstrating the advantages of an analytical gradient method over finite difference approaches for hybrid quantum-classical integration.
Abstract In this Perspective, we discuss how quantum computers may advance the simulation of nonadiabatic molecular dynamics, a framework central to describing excited-state processes in photochemistry, biology, and materials science. Classical approaches span from exponentially scaling full quantum dynamics to more approximate mixed quantum–classical techniques such as surface hopping and Ehrenfest dynamics. Hybrid quantum–classical algorithms – particularly those based on the variational quantum eigensolver – offer a transformative alternative by providing access to the key electronic properties needed to drive nonadiabatic molecular dynamics simulations, including energies, gradients, and nonadiabatic couplings. We examine recent proof-of-principle quantum simulations of reduced model systems which, despite being restricted to small molecules and limited active spaces due to constraints of qubit number and device noise, already showcase the potential of quantum devices to capture phenomena such as conical intersections and ultrafast relaxation. Although practical applications are not yet feasible in the present noisy intermediate-scale quantum era, these efforts underline the conceptual and methodological advances of quantum algorithms paving the way for large-scale quantum simulations of nonadiabatic processes. Framed within the 2025 International Year of Quantum Science and Technology, such progress exemplifies how quantum computing may open new horizons for chemistry and beyond.
Nuclear quantum effects are critical for describing proton transfer and hydrogen bonding, but their incorporation into quantum chemistry calculations is often computationally prohibitive on classical hardware. A promising alternative are quantum computers due to their linear scaling in the space requirements with system size. We introduce a novel hybrid quantum-classical algorithm for ab-initio molecular dynamics that incorporates nuclear quantum effects via the nuclear-electronic orbital method. The proposed algorithm evaluates ground state energies and forces on the quantum computer using a variational quantum eigensolver, while the molecular geometries are updated classically. We validate our approach through simulations of H2\text{H}_2H2​, H2O\text{H}_2\text{O}H2​O and the Zundel ion H5O2+\text{H}_5\text{O}_2^+H5​O2+​, comparing the simulated vibrational spectra with experimental data. Upon inclusion of nuclear quantum effects, the proton shuttling movement in the Zundel ion becomes effectively barrierless, and errors in the simulated frequencies improve significantly. Employing compact hardware-efficient ansatz circuits we achieve results comparable to the more accurate UCCSD ansatzes, which hints towards the feasibility of executing our algorithm on near-term quantum devices.
Accurately capturing electron correlation in large-scale molecular systems remains one of the foremost challenges in quantum chemistry and a primary driver for the development of quantum algorithms. Classical configuration-interaction methods, while rigorous, suffer from exponential scaling, rendering them impractical for large or strongly correlated systems. Overcoming this limitation is central to realizing the promise of quantum computing in chemistry. Here, we investigate the pyridine-Li ion complex using three quantum algorithms: the variational quantum eigensolver (VQE), the subspace quantum diagonalization (SQD) method, and the recently introduced handover iterative VQE (HI-VQE). Our results demonstrate how new generations of hybrid quantum-classical frameworks overcome the scalability and noise sensitivity that constrain conventional VQE approaches. SQD and HI-VQE achieve ground-state energy calculations for problem sizes inaccessible to classical computation, marking a clear advance toward quantum advantage. In particular, HI-VQE enables calculations within active spaces as large as (24e,22o), requiring 44 qubits-well beyond the reach of classical CASCI and VQE. This capability provides a systematic pathway for incorporating increasing numbers of electrons into quantum treatment, thereby approaching exact molecular energies. Importantly, both SQD and HI-VQE exhibit robustness against hardware noise, a critical improvement over earlier approaches. By enabling quantum simulations of molecular systems previously deemed intractable, SQD and HI-VQE offer a realistic route toward practical quantum advantage in computational chemistry. The comparison between HI-VQE and SQD shows that optimizing circuit parameters is crucial for accurate simulation.
Quantum chemistry applications on quantum computers currently rely heavily on the variational quantum eigensolver (VQE) algorithm. This hybrid quantum-classical algorithm aims at finding ground-state solutions of molecular systems based on the variational principle. VQE calculations can be systematically implemented for perturbations to each molecular degree of freedom, generating a Born-Oppenheimer potential-energy surface (PES) for the molecule. The PES can then be used to derive thermodynamic properties, which are often desirable for applications in chemical engineering and materials design. It is clear from this process that quantum chemistry applications contain a substantial classical computing component in addition to steps that can be performed using a quantum computer. In order to design efficient work flows that take full advantage of each hardware type, it is critical to consider the entire process so that the high-accuracy electronic energies possible from quantum computing are not squandered in the process of calculating thermodynamic properties. We present a summary of the hybrid quantum-classical work flow to compute thermodynamic properties. This work flow contains many options that can significantly affect the efficiency and the accuracy of the results, including classical optimizer attributes, number of ansatz repetitions, and how the vibrational Schr\"odinger equation is solved to determine vibrational modes. We also analyze the effects of these options by employing robust statistics along with simulations and experiments on actual quantum hardware. We show that, through careful selection of work flow options, nearly order-of-magnitude increases in accuracy are possible at equivalent computing time.
Quantum computing promises to efficiently and accurately solve many important problems in quantum chemistry which elude classical solvers, such as the electronic structure problem of highly correlated materials. Two leading methods in solving the ground state problem are the Variational Quantum Eigensolver (VQE) and Adiabatic Quantum Computing (AQC) algorithms. VQE often struggles with convergence due to the energy landscape being highly non-convex and the existence of barren plateaux, and implementing AQC is beyond the capabilities of current quantum devices as it requires deep circuits. Adiabatically-inspired algorithms aim to fill this gap. In this paper, we first present a unifying framework for these algorithms and then benchmark the following methods: the Adiabatically Assisted VQE (AAVQE) (Garcia-Saez and Latorre (2018)), the Variational Adiabatic Quantum Computing (VAQC) (Harwood et al (2022)), and the Adiabatic Quantum Computing with Parametrized Quantum Circuits (AQC-PQC) (Kolotouros et al (2025)) algorithms. Second, we introduce a novel hybrid approach termed G-AQC-PQC, which generalizes the AQC-PQC method, and combines adiabatic-inspired initialization with the low-memory BFGS optimizer, reducing the quantum computational cost of the method. Third, we compare the accuracy of the methods for chemistry applications using the beryllium hydride molecule (BeH2_22​). We compare the approaches across a number of different choices (ansätze types, depth, discretization steps, initial Hamiltonian, adiabatic schedules and method used). Our results show that the G-AQC-PQC outperforms conventional VQE. We further discuss limitations such as the zero-gradient problem and identify regimes where adiabatically-inspired methods offer a tangible advantage for near-term quantum chemistry applications.