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

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 Hellman-Fenyman 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 comment5 pages, 4 figures

References in corpus (23)

• A Constrained Path Monte Carlo Method for Fermion Ground States
• Quantum Monte Carlo method using phase-free random walks with Slater determinants
• Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods
• Ab initio computations of molecular systems by the auxiliary-field quantum Monte Carlo method
• Zero-Variance Zero-Bias Principle for Observables in quantum Monte Carlo: Application to Forces
• Finite-size correction in many-body electronic structure calculations
• Correlated sampling in quantum Monte Carlo: a route to forces
• Algorithmic differentiation and the calculation of forces by quantum Monte Carlo
• Auxiliary-field quantum Monte Carlo calculations of molecular systems with a Gaussian basis
• Phaseless auxiliary-field quantum Monte Carlo calculations with planewaves and pseudopotentials--applications to atoms and molecules
• Quantum Monte Carlo method for the ground state of many-boson systems
• Pressure-induced diamond to beta-tin transition in bulk silicon: a near-exact quantum Monte Carlo study
• Excited state calculations using phaseless auxiliary-field quantum Monte Carlo: potential energy curves of low lying C2 singlet states
• Eliminating spin contamination in auxiliary-field quantum Monte Carlo: realistic potential energy curve of F2
• Assessing weak hydrogen binding on Ca+ centers: An accurate many-body study with large basis sets
• Accurate, efficient and simple forces with Quantum Monte Carlo methods
• Quantum Monte Carlo calculations in solids with downfolded Hamiltonians
• An Auxiliary-Field Quantum Monte Carlo Study of the Chromium Dimer
• Computation of ground-state properties in molecular systems: back-propagation with auxiliary-field quantum Monte Carlo
• Chemical Transformations Approaching Chemical Accuracy via Correlated Sampling in Auxiliary-Field Quantum Monte Carlo
• Electronic Quantum Monte Carlo Calculations of Atomic Forces, Vibrations, and Anharmonicities
• Auxiliary-field quantum Monte Carlo calculations of the molybdenum dimer
• The Path Integral Monte Carlo Calculation of Electronic Forces

Cited by in corpus (16)

• Ab initio computations of molecular systems by the auxiliary-field quantum Monte Carlo method
• On Achieving High Accuracy in Quantum Chemical Calculations of 3d Transition Metal Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules
• Efficient ab initio auxiliary-field quantum Monte Carlo calculations in Gaussian bases via low-rank tensor decomposition
• Overcoming the Memory Bottleneck in Auxiliary Field Quantum Monte Carlo Simulations with Interpolative Separable Density Fitting
• Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C${}_{36}$ Fullerene and an Iron Porphyrin Model Complex
• Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Unit
• Hamiltonian symmetries in auxiliary-field quantum Monte Carlo calculations for electronic structure
• The performance of phaseless auxiliary-field quantum Monte Carlo on the ground state electronic energy of benzene
• Auxiliary-field quantum Monte Carlo calculations of the structural properties of nickel oxide
• Systematic Comparison and Cross-validation of Fixed-Node Diffusion Monte Carlo and Phaseless Auxiliary-Field Quantum Monte Carlo in Solids
• Stochastic Resolution-of-the-Identity Auxiliary-Field Quantum Monte Carlo: Scaling Reduction without Overhead
• Calculating ground state properties of correlated fermionic systems with BCS trial wave functions in Slater determinant path-integral approaches
• Deconvolving the components of the sign problem
• Computation of forces and stresses in solids: Towards accurate structural optimization with auxiliary-field quantum Monte Carlo
• Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids
• Algorithm for branching and population control in correlated sampling