Abstract:It is a challenge to obtain an accurate model of the state-to-state dynamics of a complex biological system from molecular dynamics (MD) simulations. In recent years, Markov State Models have gained immense popularity for computing state-to-state dynamics from a pool of short MD simulations. However, the assumption that the underlying dynamics on the reduced space is Markovian induces a systematic bias in the model, especially in biomolecular systems with complicated energy landscapes. To address this problem, we have devised a new approach we call quasi-stationary distribution kinetic Monte Carlo (QSD-KMC) that gives accurate long time state-to-state evolution while retaining the entire time resolution even when the dynamics is highly non-Markovian. The proposed method is a kinetic Monte Carlo approach that takes advantage of two concepts: (i) the quasi-stationary distribution and (ii) dynamical corrections theory. Implementation of QSD-KMC imposes stricter requirements on the lengths of the trajectories than in a Markov State Model approach, as the trajectories must be long enough to dephase. However, the QSD-KMC model produces state-to-state trajectories that are statistically indistinguishable from an MD trajectory mapped onto the discrete set of states, for an arbitrary choice of state decomposition. Furthermore, the aforementioned concepts can be used to construct a Monte Carlo approach to optimize the state boundaries regardless of the initial choice of states. We demonstrate the QSD-KMC method on two one-dimensional model systems, one of which is a driven nonequilibrium system, and on two well-characterized biomolecular systems.

A Tutorial on Advanced Dynamic Monte Carlo Methods for Systems with Discrete State Spaces

Adaptable Parallel Acceleration Strategy For Dynamic Monte Carlo Simulations Of Polymerization With Microscopic Resolution

Computing Long Timescale Biomolecular Dynamics using Quasi-Stationary Distribution Kinetic Monte Carlo (QSD-KMC)

Effective dynamics for a kinetic Monte-Carlo model with slow and fast time scales

Applications of the Monte Carlo Method in Statistical Physics

Off-Lattice Kinetic Monte Carlo Methods

Extreme Long-time Dynamic Monte Carlo Simulations

Efficient Irreversible Monte Carlo Samplers

A quantitative study on the approximation error and speed-up of the multi-scale MCMC (Monte Carlo Markov chain) method for molecular dynamics

Fast Simulation of Multicomponent Dynamic Systems

Accelerated Molecular Dynamics Methods in a Massively Parallel World

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Dynamic Monte Carlo Method for Nuclear Reactor Kinetics Calculations

Semiclassical Monte-Carlo approach for modelling non-adiabatic dynamics in extended molecules

A Review of Enhanced Sampling Approaches for Accelerated Molecular Dynamics

Noise term modeling of dynamic Monte Carlo using stochastic differential equations

System Level Numerical Analysis of a Monte Carlo Simulation of the E. Coli Chemotaxis

Discrete Molecular Dynamics

Markov Chain Monte Carlo for Continuous-Time Switching Dynamical Systems

Monte Carlo and kinetic Monte Carlo methods