Abstract:Stochastic electronic structure theories, e.g., Quantum Monte Carlo methods, enable highly accurate total energy calculations which in principle can be used to construct highly accurate potential energy surfaces. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged-elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method - previously developed for QMC structural optimization - to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and on saddle points, we show that it is possible to identify MEPs and TSs using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia molecule and an SN2 reaction. We validate our results using Density Functional Theory- and coupled cluster-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities - free energy differences, rate constants, and equilibrium constants - that incorporates stochastically-optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations, at significantly reduced computational expense.

Non-empirical weighted Langevin mechanics for the potential escape problem: parallel algorithm and application to the Argon clusters

Force-free identification of minimum-energy pathways and transition states for stochastic electronic structure theories

Nonlocal stationary probability distributions and escape rates for an active Ornstein-Uhlenbeck particle

Exit Event from a Metastable State and Eyring-Kramers Law for the Overdamped Langevin Dynamics

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker-Planck-Boltzmann Model

A Branching Random Walk Method for Many-Body Wigner Quantum Dynamics

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

Accelerating and parallelizing Lagrangian simulations of mixing-limited reactive transport

An Efficient Convergence Method for Calculating the Angular Distribution of Electron Multiple Elastic Scattering

Action Functional Gradient Descent algorithm for estimating escape paths in Stochastic Chemical Reaction Networks

Markov-chain sampling for long-range systems without evaluating the energy

Effective Langevin Equation Approach to the Molecular Diffusion on Optical Lattices

A single-walker approach for studying quasi-ergodic systems

Effective conservation of energy and momentum algorithm using switching potentials suitable for molecular dynamics simulation of thermodynamical systems

Sampling along reaction coordinates with the Wang-Landau method

A RANDOM BATCH EWALD METHOD FOR PARTICLE SYSTEMS WITH COULOMB INTERACTIONS

Efficient implementation and performance evaluation of the Wigner branching random walk

Escape rate of active particles in the effective equilibrium approach

Escape-from-a-layer approach for simulating the boundary local time in Euclidean domains

Atomistic Simulations Of Rare Events Using Gentlest Ascent Dynamics

Fast and accurate nonadiabatic molecular dynamics enabled through variational interpolation of correlated electron wavefunctions