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:2009.06803

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

Shu Kanno

Published Sep 15, 2020• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2009.06803

Related papers

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

Simulating chemical reaction dynamics on quantum computer

Qiankun Gong, Qingmin Man, Jianyu Zhao +5 more

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.

The Journal of Chemical Physics•Mar 25, 2024•7 citations

A full circuit-based quantum algorithm for excited-states in quantum chemistry

Jingwei Wen, Zhengan Wang, Chi-Tong Chen +7 more

Utilizing quantum computer to investigate quantum chemistry is an important research field nowadays. In addition to the ground-state problems that have been widely studied, the determination of excited-states plays a crucial role in the prediction and modeling of chemical reactions and other physical processes. Here, we propose a non-variational full circuit-based quantum algorithm for obtaining the excited-state spectrum of a quantum chemistry Hamiltonian. Compared with previous classical-quantum hybrid variational algorithms, our method eliminates the classical optimization process, reduces the resource cost caused by the interaction between different systems, and achieves faster convergence rate and stronger robustness against noise without barren plateau. The parameter updating for determining the next energy-level is naturally dependent on the energy measurement outputs of the previous energy-level and can be realized by only modifying the state preparation process of ancillary system, introducing little additional resource overhead. Numerical simulations of the algorithm with hydrogen, LiH, H2O and NH3 molecules are presented. Furthermore, we offer an experimental demonstration of the algorithm on a superconducting quantum computing platform, and the results show a good agreement with theoretical expectations. The algorithm can be widely applied to various Hamiltonian spectrum determination problems on the fault-tolerant quantum computers.

arXiv:2112.14193•Quantum•Jan 4, 2024•9 citations

Computational Analysis of Chemical Reactions Using a Variational Quantum Eigensolver Algorithm without Specifying Spin Multiplicity

Soichi Shirai, Hokuto Iwakiri, Keita Kanno +4 more

The analysis of a chemical reaction along the ground-state potential energy surface in conjunction with an unknown spin state is challenging because electronic states must be separately computed several times using different spin multiplicities to find the lowest energy state. However, in principle, the ground state could be obtained with just a single calculation using a quantum computer without specifying the spin multiplicity in advance. In the present work, ground-state potential energy curves for PtCO were calculated as a proof-of-concept using a variational quantum eigensolver (VQE) algorithm. This system exhibits a singlet-triplet crossover as a consequence of the interaction between Pt and CO. VQE calculations using a statevector simulator were found to converge to a singlet state in the bonding region, while a triplet state was obtained at the dissociation limit. Calculations performed using an actual quantum device provided potential energies within ±2 kcal/mol of the simulated energies after error mitigation techniques were adopted. The spin multiplicities in the bonding and dissociation regions could be clearly distinguished even in the case of a small number of shots. The results of this study suggest that quantum computing can be a powerful tool for the analysis of the chemical reactions of systems for which the spin multiplicity of the ground state and variations in this parameter are not known in advance.

arXiv:2303.05065•ACS Omega•May 25, 2023•11 citations

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

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 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.
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.
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.
Utilizing quantum computer to investigate quantum chemistry is an important research field nowadays. In addition to the ground-state problems that have been widely studied, the determination of excited-states plays a crucial role in the prediction and modeling of chemical reactions and other physical processes. Here, we propose a non-variational full circuit-based quantum algorithm for obtaining the excited-state spectrum of a quantum chemistry Hamiltonian. Compared with previous classical-quantum hybrid variational algorithms, our method eliminates the classical optimization process, reduces the resource cost caused by the interaction between different systems, and achieves faster convergence rate and stronger robustness against noise without barren plateau. The parameter updating for determining the next energy-level is naturally dependent on the energy measurement outputs of the previous energy-level and can be realized by only modifying the state preparation process of ancillary system, introducing little additional resource overhead. Numerical simulations of the algorithm with hydrogen, LiH, H2O and NH3 molecules are presented. Furthermore, we offer an experimental demonstration of the algorithm on a superconducting quantum computing platform, and the results show a good agreement with theoretical expectations. The algorithm can be widely applied to various Hamiltonian spectrum determination problems on the fault-tolerant quantum computers.
The analysis of a chemical reaction along the ground-state potential energy surface in conjunction with an unknown spin state is challenging because electronic states must be separately computed several times using different spin multiplicities to find the lowest energy state. However, in principle, the ground state could be obtained with just a single calculation using a quantum computer without specifying the spin multiplicity in advance. In the present work, ground-state potential energy curves for PtCO were calculated as a proof-of-concept using a variational quantum eigensolver (VQE) algorithm. This system exhibits a singlet-triplet crossover as a consequence of the interaction between Pt and CO. VQE calculations using a statevector simulator were found to converge to a singlet state in the bonding region, while a triplet state was obtained at the dissociation limit. Calculations performed using an actual quantum device provided potential energies within ±2 kcal/mol of the simulated energies after error mitigation techniques were adopted. The spin multiplicities in the bonding and dissociation regions could be clearly distinguished even in the case of a small number of shots. The results of this study suggest that quantum computing can be a powerful tool for the analysis of the chemical reactions of systems for which the spin multiplicity of the ground state and variations in this parameter are not known in advance.
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.
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.