Optimal reaction coordinates and kinetic rates from the projected dynamics of transition paths
arXiv:2302.12497 · doi:10.1021/acs.jctc.3c00158
Abstract
Finding optimal reaction coordinates and predicting accurate kinetic rates for activated processes are two of the foremost challenges of molecular simulations. We introduce an algorithm that tackles the two problems at once: starting from a limited number of reactive molecular dynamics trajectories (transition paths), we automatically generate with a Monte Carlo approach a sequence of different reaction coordinates that progressively reduce the kinetic rate of their projected effective dynamics. Based on a variational principle, the minimal rate accurately approximates the exact one, and it corresponds to the optimal reaction coordinate. After benchmarking the method on an analytic double-well system, we apply it to complex atomistic systems: the interaction of carbon nanoparticles of different sizes in water.
arXiv comment19 pages, 10 figures
References in corpus (16)
- Canonical sampling through velocity-rescaling
- PLUMED 2: New feathers for an old bird
- Identification of slow molecular order parameters for Markov model construction
- Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
- Molecular enhanced sampling with autoencoders: On-the-fly collective variable discovery and accelerated free energy landscape exploration
- Sparse learning of stochastic dynamic equations
- Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems
- Autonomous artificial intelligence discovers mechanisms of molecular self-organization in virtual experiments
- Enhancing entropy and enthalpy fluctuations to drive crystallization in atomistic simulations
- Effective dynamics using conditional expectations
- How the diffusivity profile reduces the arbitrariness of protein folding free energies
- Learning force fields from stochastic trajectories
- Pathways to self-organization: crystallization via nucleation and growth
- Exact dynamical coarse-graining without time-scale separation
- Building general Langevin models from discrete data sets
- Free Energy Landscapes, Diffusion Coefficients, and Kinetic Rates from Transition Paths