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:1803.05599 · The Journal of Chemical Physics

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

Mario Motta, Shiwei Zhang

Published May 8, 2018• 32 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1063/1.5029508

Related papers

Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo

Siyuan Chen, Shiwei Zhang

The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. It also paves the way for many-body computations of the phonon spectrum in solids.

arXiv:2302.07460•Physical review. B./Physical review. B•May 26, 2023•12 citations

Computing accurate forces in quantum Monte Carlo using Pulay’s corrections and energy minimization

Mosè Casalegno, Massimo Mella, Andrew M. Rappe

In order to overcome the difficulty of optimizing molecular geometry using quantum Monte Carlo methods, we introduce various approximations to the exact force expectation value. We follow Pulay’s suggestion [Mol. Phys. 17, 153 (1969)] to correct the Hellmann–Feynman estimator by introducing the contributions due to the changes in the wave function with respect to the nuclear positions. When used in conjunction with energy-optimized explicitly correlated trial wave functions for H2 and LiH, these approximations appear to yield accurate forces using both the variational and diffusion Monte Carlo methods. Also, the accuracy of the second-order estimate of the Hellmann–Feynman force estimator was investigated employing our energy-optimized trial wave functions, and an erratic behavior was uncovered for some of the studied bond lengths. The additional computational cost required to compute the corrections to the Hellmann–Feynman estimator was found to be only a small fraction of the cost for a simple mean energy calculation. The same approach could be exploited also in computing the derivative of other energy-dependent quantum-mechanical observables.

The Journal of Chemical Physics•Apr 11, 2003•51 citations

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

Juntao Lai, Yi Fan, Qiang Fu +2 more

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.

The Journal of Chemical Physics•Sep 19, 2023•14 citations

Nuclear gradients from auxiliary-field quantum Monte Carlo and their application in geometry optimization and transition state search

Jo S. Kurian, Ankit Mahajan, Sandeep Sharma

In this article, we present a method for computing accurate and scalable nuclear forces within the phaseless auxiliary-field quantum Monte Carlo (AFQMC) framework. Our approach leverages automatic differentiation of the energy functional to obtain nuclear gradients at a computational cost comparable to that of energy evaluation. The accuracy of the method is validated against finite difference calculations, showing excellent agreement. We then explore several machine learning (ML) strategies for learning noisy AFQMC data. These ML potentials are subsequently used to perform geometry optimizations and nudged elastic band (NEB) calculations, successfully identifying the transition state of the formamide-formimidic acid tautomerization. The resulting transition state geometry and barrier heights are in close agreement with coupled-cluster reference values. This work paves the way for highly accurate geometry optimization, molecular dynamics, or reaction path calculations.

arXiv:2602.13187•Feb 13, 2026•0 citations

Electronic Quantum Monte Carlo Calculations of Energies and Atomic Forces for Diatomic and Polyatomic Molecules

Myung Won Lee, Massimo Mella, Andrew M. Rappe

We calculated the energies and atomic forces of first-row monohydrides, carbon monoxide, and small organic polyatomic molecules using quantum Monte Carlo (QMC) method. Accurate forces were obtained with the method of Casalegno, Mella, and Rappe, combining the Hellmann- Feynman theorem forces calculated by the Assaraf-Caffarel method with a many-body Pulay correction. Improved algorithms for the minimization of the variational integral were useful in the force calculations.

ACS symposium series•Dec 31, 2006•0 citations

Nuclear Gradients from Auxiliary-Field Quantum Monte Carlo and Their Applications in ML-Driven Geometry Optimization and Transition State Search

Jo S. Kurian, Ankit Mahajan, Sandeep Sharma

In this article, we present a method for computing accurate and scalable nuclear forces within the phaseless auxiliary-field quantum Monte Carlo (AFQMC) framework. Our approach leverages automatic differentiation of the energy functional to obtain nuclear gradients at a computational cost comparable to that of energy evaluation. The accuracy of the method is validated against finite difference calculations, showing excellent agreement. We then explore several machine learning (ML) strategies for learning noisy AFQMC data. These ML potentials are subsequently used to perform geometry optimizations and nudged elastic band (NEB) calculations, successfully identifying the transition state of the formamide-formimidic acid tautomerization. The resulting transition-state geometry and barrier heights are in close agreement with coupled-cluster reference values. This work paves the way for highly accurate geometry optimization, molecular dynamics, or reaction path calculations.

Journal of Chemical Theory and Computation•Jul 1, 2026•0 citations
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.
The accurate computation of forces and other energy derivatives has been a long-standing challenge for quantum Monte Carlo methods. A number of technical obstacles contribute to this challenge. We discuss how these obstacles can be removed with the auxiliary-field quantum Monte Carlo (AFQMC) approach. AFQMC is a general, high-accuracy, many-body total-energy method for molecules and solids. The implementation of back-propagation for pure estimators allows direct calculation of gradients of the energy via the Hellmann-Feynman theorem. A planewave basis with norm-conserving pseudopotentials is used for the study of periodic bulk materials. Completeness of the planewave basis minimizes the effect of so-called Pulay terms. The ionic pseudopotentials, which can be incorporated in AFQMC in exactly the same manner as in standard independent-electron methods, regulate the force and stress estimators and eliminate any potential divergence of the Monte Carlo variances. The resulting approach allows applications of full geometry optimizations in bulk materials. It also paves the way for many-body computations of the phonon spectrum in solids.
In order to overcome the difficulty of optimizing molecular geometry using quantum Monte Carlo methods, we introduce various approximations to the exact force expectation value. We follow Pulay’s suggestion [Mol. Phys. 17, 153 (1969)] to correct the Hellmann–Feynman estimator by introducing the contributions due to the changes in the wave function with respect to the nuclear positions. When used in conjunction with energy-optimized explicitly correlated trial wave functions for H2 and LiH, these approximations appear to yield accurate forces using both the variational and diffusion Monte Carlo methods. Also, the accuracy of the second-order estimate of the Hellmann–Feynman force estimator was investigated employing our energy-optimized trial wave functions, and an erratic behavior was uncovered for some of the studied bond lengths. The additional computational cost required to compute the corrections to the Hellmann–Feynman estimator was found to be only a small fraction of the cost for a simple mean energy calculation. The same approach could be exploited also in computing the derivative of other energy-dependent quantum-mechanical observables.
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.
In this article, we present a method for computing accurate and scalable nuclear forces within the phaseless auxiliary-field quantum Monte Carlo (AFQMC) framework. Our approach leverages automatic differentiation of the energy functional to obtain nuclear gradients at a computational cost comparable to that of energy evaluation. The accuracy of the method is validated against finite difference calculations, showing excellent agreement. We then explore several machine learning (ML) strategies for learning noisy AFQMC data. These ML potentials are subsequently used to perform geometry optimizations and nudged elastic band (NEB) calculations, successfully identifying the transition state of the formamide-formimidic acid tautomerization. The resulting transition state geometry and barrier heights are in close agreement with coupled-cluster reference values. This work paves the way for highly accurate geometry optimization, molecular dynamics, or reaction path calculations.
In this article, we present a method for computing accurate and scalable nuclear forces within the phaseless auxiliary-field quantum Monte Carlo (AFQMC) framework. Our approach leverages automatic differentiation of the energy functional to obtain nuclear gradients at a computational cost comparable to that of energy evaluation. The accuracy of the method is validated against finite difference calculations, showing excellent agreement. We then explore several machine learning (ML) strategies for learning noisy AFQMC data. These ML potentials are subsequently used to perform geometry optimizations and nudged elastic band (NEB) calculations, successfully identifying the transition state of the formamide-formimidic acid tautomerization. The resulting transition-state geometry and barrier heights are in close agreement with coupled-cluster reference values. This work paves the way for highly accurate geometry optimization, molecular dynamics, or reaction path calculations.