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 Chemical Physics

Simulating chemical reaction dynamics on quantum computer

Qiankun Gong, Qingmin Man, Jianyu Zhao, Ye Li, Menghan Dou, Qingchun Wang, Yu-Chun Wu, Guo‐Ping Guo

Published Mar 25, 2024• 7 citations
Publisher

Abstract

The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.

DOI
10.1063/5.0192036

Related papers

Simulating Chemical Reaction Dynamics Using Sample-Based Quantum Diagonalization: A Quantum-Centric Supercomputing Approach

Sumit Kumar, Kalpak Ghosh, Deepak Kumar +1 more

Quantum computers offer the potential to simulate chemical systems beyond the capabilities of classical computers. However, current quantum computing approaches to ab initio molecular dynamics (AIMD) and vibrational analysis largely rely on the Variational Quantum Eigensolver (VQE), whose circuit depth, iterative optimization, and measurement overhead hinder the practical simulation of molecular systems. Here, we extend Sample-based Quantum Diagonalization (SQD), a quantum sampling-based electronic-structure method, to AIMD and vibrational analysis on noisy quantum hardware. The resulting framework enables nuclear forces and Hessians to be derived from quantum-computed electronic energies, facilitating molecular dynamics and characterization of potential energy surfaces. We demonstrate the approach using diverse molecular processes, including H2 bond vibration, the proton-exchange reaction (H+ + H2), the bimolecular nucleophilic substitution (SN2) reaction (Cl- + CH3Cl), and proton-transfer dynamics between two water molecules involving transient formation of the Zundel cation (H2O···H+···OH2). Across these systems, SQD-AIMD trajectories and vibrational frequencies obtained with noisy quantum hardware/noise-models show excellent agreement with classical reference calculations and substantially improved accuracy relative to previously reported VQE-based simulations. Importantly, the SQD-based approach eliminates costly variational optimization and significantly reduces measurement requirements relative to VQE, enabling AIMD simulations spanning approximately 500 timesteps (∼200 fs), substantially longer than previously reported AIMD simulations on quantum hardware. These results establish SQD-AIMD and SQD-vibrational analysis as practical and scalable frameworks for simulating chemical reaction dynamics and elucidating reaction mechanisms on near-term quantum hardware, extending quantum computing applications beyond static electronic structure calculations.

Journal of the American Chemical Society•Oct 12, 2026•0 citations

Simulation of chemical reaction dynamics based on quantum computing

Qiankun Gong, Qingmin Man, Ye Li +4 more

The molecular energies of chemical systems have been successfully calculated on quantum computers, however, more attention has been paid to the dynamic process of chemical reactions in practical application, especially in catalyst design, material synthesis. Due to the limited the capabilities of the noisy intermediate scale quantum (NISQ) devices, directly simulating the reaction dynamics and determining reaction pathway still remain a challenge. Here we develop the ab initio molecular dynamics based on quantum computing to simulate reaction dynamics by extending correlated sampling approach. And, we use this approach to calculate Hessian matrix and evaluate computation resources. We test the performance of our approach by simulating hydrogen exchange reaction and bimolecular nucleophilic substitution SN2 reaction. Our results suggest that it is reliable to characterize the molecular structure, property, and reactivity, which is another important expansion of the application of quantum computing

arXiv:2303.08571•Mar 15, 2023•2 citations

Quantum algorithm for a chemical reaction path optimization by using a variational quantum algorithm and a reaction path generation

Shu Kanno

The search for new computational tasks of quantum chemistry that can be performed on current quantum computers is important for the development of quantum computing and quantum chemistry. Although calculations of chemical reactions have a wide range of applications in quantum chemical calculations, a quantum algorithm for obtaining activation energy $E_a$, which determines the rate of chemical reactions, has not been performed. In this study, we propose a quantum algorithm for the chemical reaction path optimization to obtain $E_a$. In our algorithm, quantum circuits can be used not only for the energy evaluation by the variational quantum eigensolver (VQE) but also for chemical reaction path generation. The chemical reaction path is obtained by encoding the initial reaction path to the circuit, operating parameterized gates, and extracting the path information by measurement. The nudged elastic band method was used for optimizing a reaction path, and the ground-state calculation for each state on the path was performed by the VQE or the exact diagonalization (ED). The proposed algorithm was applied to $\mathrm{H}_2 + \mathrm{H} \rightarrow \mathrm{H}_2 + \mathrm{H}$ reaction, and we confirmed that $E_a$ was obtained accurately in the case of both the VQE and the ED. We also obtained numerical results that the entanglement between the images accelerates the path optimization. From these results, we show the feasibility of performing fast and accurate chemical reaction calculations by using quantum algorithms.

arXiv:2009.06803•Sep 15, 2020•1 citation

Assessing Ground State Energy of Molecules and Energy Profile of the NH3 Capturing CO2 System Using the Quantum Computing Algorithms

Yueh‐Lin Lee, Manh Tien Nguyen, Dominic R. Alfonso +5 more

Molecule size correlates with the number of electrons on electronic energies and strength of anharmonicity on vibrational properties, however, it is challenging to address using classical computing. In this study, variational quantum eigensolver (VQE) algorithm was implemented on a quantum simulator to quantify electronic and vibrational energies and reaction pathways of CO2 + NH3 = NH2COOH. The VQE-based Hartree-Fock-Embedding algorithm was adopted to benchmark electronic energies for a series of molecules (doi.org/10.1063/5.0188249) and quantify the reaction energy profile of the CO2 capture reaction (doi.org/10.1116/5.0137750). The generated reaction profile is in good agreement with the classical high-level Coupled-Cluster-Singles-and-Doubles (CCSD) results. The quantum computing algorithm also helps enhance the calculation of vibrational ground-state energies by considering the many-body coupling using the Vibrational Self-Consistent Field method, providing results for CO2 and NH3 molecules with accuracy comparable to the direct diagonalization method. Our approach indicates quantum computing can be applied to solve practical problems.

•Oct 9, 2024•0 citations

Optimization of chemical reaction path via molecular geometry generation using quantum circuits

Shu Kanno

The search for new tasks related to quantum chemical calculations that can be performed on current quantum computers is important for the development of both quantum computing and quantum chemical calculations. Although chemical reaction calculations have a wide range of applications in quantum chemical calculations, there is no quantum algorithm for obtaining the activation energy ($E_a$), which determines the rate of chemical reactions for any given substance. In this study, we propose a quantum algorithm for obtaining $E_a$ for a given substance by optimizing the reaction path via molecular geometry generation using quantum circuits. The nudged elastic band method was used to optimize a reaction path, and the ground state calculation for each state on the reaction path involved either variational quantum eigensolver (VQE) or exact diagonalization (ED). The proposed algorithm was applied to H2 + H -> H + H2 reaction, and the reaction path was correctly optimized. Moreover, the activation energy was obtained with good accuracy regardless of whether VQE or ED is used for the ground state calculation. The results demonstrate the feasibility of performing the chemical reaction calculations using quantum algorithms.

•Sep 15, 2020•0 citations

Theoretical analysis of chemical reactions using a variational quantum eigensolver method without specifying molecular charge

Soichi Shirai, Takahiro Horiba, Hirotoshi Hirai

Quantum chemical calculations have attracted much attention as a practical application of quantum computing. Quantum computers can prepare superpositions of electronic states with various numbers of electrons on qubits. This special feature could be used to construct an efficient method for analyzing the structural variations of molecules and chemical reactions involving changes in molecular charge. The present work demonstrates a variational quantum eigensolver (VQE) algorithm based on a cost function ($L_{cost}$) having the same form as the grand potential of the grand canonical ensemble of electrons. The chemical potential of the electrons ($w$) is used as an input to these VQE calculations, whereas the molecular charge is not specified in advance but rather is a physical quantity that results from the calculations. Calculations involving model systems are carried out to show the viability of this new approach. Calculations for typical electron-donating and electron-accepting molecules using this technique yielded cationic, neutral or anionic species depending on the value of $w$. Models representing the adsorption of water or ammonia on copper-based catalysts predicted that oxidation would be associated with such adsorption. The molecular structures in which such reactions occurred were found to be dependent on the catalyst model, the adsorbed molecular species, and the value of $w$. These results arise because the electronic state that gives the lowest $L_{cost}$ value depends on the value of $w$ and the molecular structure. This behaviour was successfully simulated by the present VQE calculations.

arXiv:2403.19296•Mar 28, 2024•0 citations
The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.
Quantum computers offer the potential to simulate chemical systems beyond the capabilities of classical computers. However, current quantum computing approaches to ab initio molecular dynamics (AIMD) and vibrational analysis largely rely on the Variational Quantum Eigensolver (VQE), whose circuit depth, iterative optimization, and measurement overhead hinder the practical simulation of molecular systems. Here, we extend Sample-based Quantum Diagonalization (SQD), a quantum sampling-based electronic-structure method, to AIMD and vibrational analysis on noisy quantum hardware. The resulting framework enables nuclear forces and Hessians to be derived from quantum-computed electronic energies, facilitating molecular dynamics and characterization of potential energy surfaces. We demonstrate the approach using diverse molecular processes, including H2 bond vibration, the proton-exchange reaction (H+ + H2), the bimolecular nucleophilic substitution (SN2) reaction (Cl- + CH3Cl), and proton-transfer dynamics between two water molecules involving transient formation of the Zundel cation (H2O···H+···OH2). Across these systems, SQD-AIMD trajectories and vibrational frequencies obtained with noisy quantum hardware/noise-models show excellent agreement with classical reference calculations and substantially improved accuracy relative to previously reported VQE-based simulations. Importantly, the SQD-based approach eliminates costly variational optimization and significantly reduces measurement requirements relative to VQE, enabling AIMD simulations spanning approximately 500 timesteps (∼200 fs), substantially longer than previously reported AIMD simulations on quantum hardware. These results establish SQD-AIMD and SQD-vibrational analysis as practical and scalable frameworks for simulating chemical reaction dynamics and elucidating reaction mechanisms on near-term quantum hardware, extending quantum computing applications beyond static electronic structure calculations.
The molecular energies of chemical systems have been successfully calculated on quantum computers, however, more attention has been paid to the dynamic process of chemical reactions in practical application, especially in catalyst design, material synthesis. Due to the limited the capabilities of the noisy intermediate scale quantum (NISQ) devices, directly simulating the reaction dynamics and determining reaction pathway still remain a challenge. Here we develop the ab initio molecular dynamics based on quantum computing to simulate reaction dynamics by extending correlated sampling approach. And, we use this approach to calculate Hessian matrix and evaluate computation resources. We test the performance of our approach by simulating hydrogen exchange reaction and bimolecular nucleophilic substitution SN2 reaction. Our results suggest that it is reliable to characterize the molecular structure, property, and reactivity, which is another important expansion of the application of quantum computing
The search for new computational tasks of quantum chemistry that can be performed on current quantum computers is important for the development of quantum computing and quantum chemistry. Although calculations of chemical reactions have a wide range of applications in quantum chemical calculations, a quantum algorithm for obtaining activation energy EaE_aEa​, which determines the rate of chemical reactions, has not been performed. In this study, we propose a quantum algorithm for the chemical reaction path optimization to obtain EaE_aEa​. In our algorithm, quantum circuits can be used not only for the energy evaluation by the variational quantum eigensolver (VQE) but also for chemical reaction path generation. The chemical reaction path is obtained by encoding the initial reaction path to the circuit, operating parameterized gates, and extracting the path information by measurement. The nudged elastic band method was used for optimizing a reaction path, and the ground-state calculation for each state on the path was performed by the VQE or the exact diagonalization (ED). The proposed algorithm was applied to H2+H→H2+H\mathrm{H}_2 + \mathrm{H} \rightarrow \mathrm{H}_2 + \mathrm{H}H2​+H→H2​+H reaction, and we confirmed that EaE_aEa​ was obtained accurately in the case of both the VQE and the ED. We also obtained numerical results that the entanglement between the images accelerates the path optimization. From these results, we show the feasibility of performing fast and accurate chemical reaction calculations by using quantum algorithms.
Molecule size correlates with the number of electrons on electronic energies and strength of anharmonicity on vibrational properties, however, it is challenging to address using classical computing. In this study, variational quantum eigensolver (VQE) algorithm was implemented on a quantum simulator to quantify electronic and vibrational energies and reaction pathways of CO2 + NH3 = NH2COOH. The VQE-based Hartree-Fock-Embedding algorithm was adopted to benchmark electronic energies for a series of molecules (doi.org/10.1063/5.0188249) and quantify the reaction energy profile of the CO2 capture reaction (doi.org/10.1116/5.0137750). The generated reaction profile is in good agreement with the classical high-level Coupled-Cluster-Singles-and-Doubles (CCSD) results. The quantum computing algorithm also helps enhance the calculation of vibrational ground-state energies by considering the many-body coupling using the Vibrational Self-Consistent Field method, providing results for CO2 and NH3 molecules with accuracy comparable to the direct diagonalization method. Our approach indicates quantum computing can be applied to solve practical problems.
The search for new tasks related to quantum chemical calculations that can be performed on current quantum computers is important for the development of both quantum computing and quantum chemical calculations. Although chemical reaction calculations have a wide range of applications in quantum chemical calculations, there is no quantum algorithm for obtaining the activation energy (EaE_aEa​), which determines the rate of chemical reactions for any given substance. In this study, we propose a quantum algorithm for obtaining EaE_aEa​ for a given substance by optimizing the reaction path via molecular geometry generation using quantum circuits. The nudged elastic band method was used to optimize a reaction path, and the ground state calculation for each state on the reaction path involved either variational quantum eigensolver (VQE) or exact diagonalization (ED). The proposed algorithm was applied to H2 + H -> H + H2 reaction, and the reaction path was correctly optimized. Moreover, the activation energy was obtained with good accuracy regardless of whether VQE or ED is used for the ground state calculation. The results demonstrate the feasibility of performing the chemical reaction calculations using quantum algorithms.
Quantum chemical calculations have attracted much attention as a practical application of quantum computing. Quantum computers can prepare superpositions of electronic states with various numbers of electrons on qubits. This special feature could be used to construct an efficient method for analyzing the structural variations of molecules and chemical reactions involving changes in molecular charge. The present work demonstrates a variational quantum eigensolver (VQE) algorithm based on a cost function (LcostL_{cost}Lcost​) having the same form as the grand potential of the grand canonical ensemble of electrons. The chemical potential of the electrons (www) is used as an input to these VQE calculations, whereas the molecular charge is not specified in advance but rather is a physical quantity that results from the calculations. Calculations involving model systems are carried out to show the viability of this new approach. Calculations for typical electron-donating and electron-accepting molecules using this technique yielded cationic, neutral or anionic species depending on the value of www. Models representing the adsorption of water or ammonia on copper-based catalysts predicted that oxidation would be associated with such adsorption. The molecular structures in which such reactions occurred were found to be dependent on the catalyst model, the adsorbed molecular species, and the value of www. These results arise because the electronic state that gives the lowest LcostL_{cost}Lcost​ value depends on the value of www and the molecular structure. This behaviour was successfully simulated by the present VQE calculations.