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

Accurate and efficient calculations of Hellmann–Feynman forces for quantum computation

Juntao Lai, Yi Fan, Qiang Fu, Zhenyu Li, Jinlong Yang

Published Sep 19, 2023• 14 citations
Publisher

Abstract

First-order derivatives of energies with respect to atomic coordinates are widely computed and used in quantum chemistry simulations. The rapidly emerging technology of quantum computing offers a new paradigm for solving relevant quantum chemistry equations. In this work, we have achieved analytical calculations of atomic forces based on the Hellmann-Feynman theorem within the framework of the variational quantum eigensolver. The accuracy of the approach is demonstrated by calculating the atomic forces of H2, LiH, H2O, and NH3 molecules, which are in excellent agreement with values obtained from full configuration interaction calculations. In particular, for systems with degenerate molecular orbitals, the analytical approach has a significant accuracy advantage over finite-difference-based methods and will not involve additional computational effort on a quantum computer. The calculated forces are further used to optimize the geometries of NH3 and CH4 molecules and to perform ab initio molecular dynamics simulations for the umbrella inversion of NH3, demonstrating the feasibility of the approach in practical quantum chemistry simulations.

DOI
10.1063/5.0167444

Related papers

Communication: Calculation of interatomic forces and optimization of molecular geometry with auxiliary-field quantum Monte Carlo

Mario Motta, Shiwei Zhang

We propose an algorithm for accurate, systematic, and scalable computation of interatomic forces within the auxiliary-field quantum Monte Carlo (AFQMC) method. The algorithm relies on the Hellmann-Feynman theorem and incorporates Pulay corrections in the presence of atomic orbital basis sets. We benchmark the method for small molecules by comparing the computed forces with the derivatives of the AFQMC potential energy surface and by direct comparison with other quantum chemistry methods. We then perform geometry optimizations using the steepest descent algorithm in larger molecules. With realistic basis sets, we obtain equilibrium geometries in agreement, within statistical error bars, with experimental values. The increase in computational cost for computing forces in this approach is only a small prefactor over that of calculating the total energy. This paves the way for a general and efficient approach for geometry optimization and molecular dynamics within AFQMC.

arXiv:1803.05599•The Journal of Chemical Physics•May 8, 2018•32 citations

Accurate Vibrational Frequency Calculations for Quantum Computing via an Analytic Second-Order Energy Derivative Framework

Juntao Lai, Qiang Fu, Zhenyu Li +1 more

Quantum computing has emerged as a promising paradigm for tackling electronic structure problems, with most efforts to date focused on molecular energies and, more recently, first-order derivatives. However, the extension to second-order energy derivatives with respect to nuclear coordinates─essential for predicting vibrational spectra and identifying transition states─has remained relatively limited. Here, we present an analytic implementation for computing nuclear Hessians within the variational quantum eigensolver framework. Our approach produces harmonic vibrational frequencies and normal modes in excellent agreement with full configuration interaction (FCI) benchmarks, even for systems with challenging cases of orbital degeneracy or involving weak intermolecular interactions. We further assess the quantum measurement cost of these second-order derivative calculations and compare it with that of variational quantum eigensolver (VQE) energy calculations. Additionally, we demonstrate how point group symmetry can be incorporated to reduce measurement cost without loss of accuracy. This work extends quantum computing capabilities toward advanced quantum chemical simulations, such as the characterization of noncovalent interactions via low-frequency vibrational signatures and the identification of transition states.

Journal of Chemical Theory and Computation•Apr 17, 2026•0 citations

Ab initio molecular dynamics on quantum computers

Dmitry A. Fedorov, Matthew Otten, Stephen K. Gray +1 more

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann–Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

arXiv:2008.06562•The Journal of Chemical Physics•Apr 23, 2021•34 citations

An Interface for Variational Quantum Eigensolver based Energy (VQE-E) and Force (VQE-F) Calculator to Atomic Simulation Environment (ASE)

M R Nirmal, Sarkar, Shampa, Manoj Karunakaran Nambiar

The development of quantum algorithms to solve quantum chemistry problems has offered a promising new paradigm of performing computer simulations at the scale of atoms and molecules. Although majority of the research so far has focused on designing quantum algorithms to compute ground and excited state energies and forces, it is useful to run different simulation tasks, such as geometry optimization, with these algorithms as subroutines. Towards this end, we have created an interface for the Variational Quantum Eigensolver based molecular Energy (VQE-E) and molecular Force (VQE-F) code to the Atomic Simulation Environment (ASE). We demonstrate the working of this hybrid quantum-classical interface by optimizing the geometry of water molecule using a native optimizer implemented in ASE. Furthermore, this interface enables one to compare, combine and use quantum algorithms in conjunction with related classical methods quite easily with minimal coding effort.

arXiv:2209.13876•Sep 28, 2022•2 citations

Accurate Calculation of Interatomic Forces with Neural Networks Based on a Generative Transformer Architecture

Juntao Lai, Bowen Kan, Yangjun Wu +4 more

Using neural networks to express electronic wave functions represents a new paradigm for solving the Schrödinger equation in quantum chemistry. For practical quantum chemistry simulations, one needs to know not only energies of molecules, but also accurate forces acting on constituent atoms. In this work, we achieve the accurate calculation of interatomic forces on QiankunNet, a platform that combines transformer-based deep neural networks with efficient batched autoregressive sampling. Our approach permits the application of the Hellmann-Feynman theorem to force calculations without introducing corrective Pulay terms. The results show that the calculated interatomic forces are in close agreement with those derived from the full configuration interaction method, irrespective of whether the system is a simple molecule or a strongly correlated electron system like a linear hydrogen chain. Furthermore, the calculated interatomic forces are employed for atomic relaxation in the torsional rotation process of ethylene, and the energy barrier obtained from the scanned potential energy surface is in excellent agreement with the experiment. Our work contributes to the application of artificial intelligence to broader quantum chemistry simulations, such as modeling challenging chemical transformations where electron correlations are difficult to describe.

Journal of Chemical Theory and Computation•Oct 23, 2024•12 citations

Calculating energy derivatives for quantum chemistry on a quantum computer

Thomas E. O’Brien, Bruno Senjean, Ramiro Sagastizabal +5 more

Abstract Modeling chemical reactions and complicated molecular systems has been proposed as the “killer application” of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential toward this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within $$0.014$$ 0.014 Å of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H $${}_{2}$$ 2 molecule, finding agreement at the equilibrium bond length to within $$0.06$$ 0.06 a.u. ( $$2 \%$$ 2% relative error).

arXiv:1905.03742•npj Quantum Information•Dec 11, 2019•101 citations
First-order derivatives of energies with respect to atomic coordinates are widely computed and used in quantum chemistry simulations. The rapidly emerging technology of quantum computing offers a new paradigm for solving relevant quantum chemistry equations. In this work, we have achieved analytical calculations of atomic forces based on the Hellmann-Feynman theorem within the framework of the variational quantum eigensolver. The accuracy of the approach is demonstrated by calculating the atomic forces of H2, LiH, H2O, and NH3 molecules, which are in excellent agreement with values obtained from full configuration interaction calculations. In particular, for systems with degenerate molecular orbitals, the analytical approach has a significant accuracy advantage over finite-difference-based methods and will not involve additional computational effort on a quantum computer. The calculated forces are further used to optimize the geometries of NH3 and CH4 molecules and to perform ab initio molecular dynamics simulations for the umbrella inversion of NH3, demonstrating the feasibility of the approach in practical quantum chemistry simulations.
We propose an algorithm for accurate, systematic, and scalable computation of interatomic forces within the auxiliary-field quantum Monte Carlo (AFQMC) method. The algorithm relies on the Hellmann-Feynman theorem and incorporates Pulay corrections in the presence of atomic orbital basis sets. We benchmark the method for small molecules by comparing the computed forces with the derivatives of the AFQMC potential energy surface and by direct comparison with other quantum chemistry methods. We then perform geometry optimizations using the steepest descent algorithm in larger molecules. With realistic basis sets, we obtain equilibrium geometries in agreement, within statistical error bars, with experimental values. The increase in computational cost for computing forces in this approach is only a small prefactor over that of calculating the total energy. This paves the way for a general and efficient approach for geometry optimization and molecular dynamics within AFQMC.
Quantum computing has emerged as a promising paradigm for tackling electronic structure problems, with most efforts to date focused on molecular energies and, more recently, first-order derivatives. However, the extension to second-order energy derivatives with respect to nuclear coordinates─essential for predicting vibrational spectra and identifying transition states─has remained relatively limited. Here, we present an analytic implementation for computing nuclear Hessians within the variational quantum eigensolver framework. Our approach produces harmonic vibrational frequencies and normal modes in excellent agreement with full configuration interaction (FCI) benchmarks, even for systems with challenging cases of orbital degeneracy or involving weak intermolecular interactions. We further assess the quantum measurement cost of these second-order derivative calculations and compare it with that of variational quantum eigensolver (VQE) energy calculations. Additionally, we demonstrate how point group symmetry can be incorporated to reduce measurement cost without loss of accuracy. This work extends quantum computing capabilities toward advanced quantum chemical simulations, such as the characterization of noncovalent interactions via low-frequency vibrational signatures and the identification of transition states.
Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann–Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.
The development of quantum algorithms to solve quantum chemistry problems has offered a promising new paradigm of performing computer simulations at the scale of atoms and molecules. Although majority of the research so far has focused on designing quantum algorithms to compute ground and excited state energies and forces, it is useful to run different simulation tasks, such as geometry optimization, with these algorithms as subroutines. Towards this end, we have created an interface for the Variational Quantum Eigensolver based molecular Energy (VQE-E) and molecular Force (VQE-F) code to the Atomic Simulation Environment (ASE). We demonstrate the working of this hybrid quantum-classical interface by optimizing the geometry of water molecule using a native optimizer implemented in ASE. Furthermore, this interface enables one to compare, combine and use quantum algorithms in conjunction with related classical methods quite easily with minimal coding effort.
Using neural networks to express electronic wave functions represents a new paradigm for solving the Schrödinger equation in quantum chemistry. For practical quantum chemistry simulations, one needs to know not only energies of molecules, but also accurate forces acting on constituent atoms. In this work, we achieve the accurate calculation of interatomic forces on QiankunNet, a platform that combines transformer-based deep neural networks with efficient batched autoregressive sampling. Our approach permits the application of the Hellmann-Feynman theorem to force calculations without introducing corrective Pulay terms. The results show that the calculated interatomic forces are in close agreement with those derived from the full configuration interaction method, irrespective of whether the system is a simple molecule or a strongly correlated electron system like a linear hydrogen chain. Furthermore, the calculated interatomic forces are employed for atomic relaxation in the torsional rotation process of ethylene, and the energy barrier obtained from the scanned potential energy surface is in excellent agreement with the experiment. Our work contributes to the application of artificial intelligence to broader quantum chemistry simulations, such as modeling challenging chemical transformations where electron correlations are difficult to describe.
Abstract Modeling chemical reactions and complicated molecular systems has been proposed as the “killer application” of a future quantum computer. Accurate calculations of derivatives of molecular eigenenergies are essential toward this end, allowing for geometry optimization, transition state searches, predictions of the response to an applied electric or magnetic field, and molecular dynamics simulations. In this work, we survey methods to calculate energy derivatives, and present two new methods: one based on quantum phase estimation, the other on a low-order response approximation. We calculate asymptotic error bounds and approximate computational scalings for the methods presented. Implementing these methods, we perform geometry optimization on an experimental quantum processor, estimating the equilibrium bond length of the dihydrogen molecule to within 0.0140.0140.014 0.014 Å of the full configuration interaction value. Within the same experiment, we estimate the polarizability of the H 2{}_{2}2​ 2 molecule, finding agreement at the equilibrium bond length to within 0.060.060.06 0.06 a.u. ( 2%2 \%2% 2% relative error).