Velocity Jumps for Molecular Dynamics
arXiv:2412.15073 · doi:10.1021/acs.jctc.5c00023
Abstract
We introduce the Velocity Jumps approach, denoted as JUMP, a new class of Molecular dynamics integrators, replacing the Langevin dynamics by a hybrid model combining a classical Langevin diffusion and a piecewise deterministic Markov process, where the expensive computation of long-range pairwise interactions is replaced by a resampling of the velocities at random times. This framework allows for an acceleration in the simulation speed while preserving sampling and dynamical properties such as the diffusion constant. It can also be integrated in classical multi-timestep methods, pushing further the computational speedup, while avoiding some of the resonance issues of the latter thanks to the random nature of jumps. The JUMP, JUMP-RESPA and JUMP-RESPA1 integrators have been implemented in the GPU-accelerated version of the Tinker-HP package and are shown to provide significantly enhanced performances compared to their BAOAB, BAOAB-RESPA and BAOAB-RESPA1 counterparts respectively.
References in corpus (16)
- Robust and efficient configurational molecular sampling via Langevin Dynamics
- Rational Construction of Stochastic Numerical Methods for Molecular Sampling
- The Zig-Zag Process and Super-Efficient Sampling for Bayesian Analysis of Big Data
- Rejection-free Monte-Carlo sampling for general potentials
- The computation of averages from equilibrium and nonequilibrium Langevin molecular dynamics
- Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances
- Tinker-HP : Accelerating Molecular Dynamics Simulations of Large Complex Systems with Advanced Point Dipole Polarizable Force Fields using GPUs and Multi-GPUs systems
- Accurate biomolecular simulations account for electronic polarization
- Stochastic resonance-free multiple time-step algorithm for molecular dynamics with very large time steps
- Piecewise-Deterministic Markov Chain Monte Carlo
- Development of the Quantum Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed Phase Molecular Dynamics Simulations
- Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: from Force Fields to Machine Learning Potentials
- Exact simulation of the jump times of a class of Piecewise Deterministic Markov Processes
- Lambda-ABF: Simplified, Portable, Accurate and Cost-effective Alchemical Free Energy Computations
- ANKH: A Generalized O(N) Interpolated Ewald Strategy for Molecular Dynamics Simulations
- The velocity jump Langevin process and its splitting scheme: long time convergence and numerical accuracy