Free Energy Landscapes, Diffusion Coefficients, and Kinetic Rates from Transition Paths
arXiv:2106.05415 · doi:10.1021/acs.jctc.2c00324
Abstract
We address the problem of constructing accurate mathematical models of the dynamics of complex systems projected on a collective variable. To this aim we introduce a conceptually simple yet effective algorithm for estimating the parameters of Langevin and Fokker-Planck equations from a set of short, possibly out-of-equilibrium molecular dynamics trajectories, obtained for instance from transition path sampling or as relaxation from high free-energy configurations. The approach maximizes the model likelihood based on any explicit expression of the short-time propagator, hence it can be applied to different evolution equations. We demonstrate the numerical efficiency and robustness of the algorithm on model systems, and we apply it to reconstruct the projected dynamics of pairs of C60 and C240 fullerene molecules in explicit water. Our methodology allows reconstructing the accurate thermodynamics and kinetics of activated processes, namely free energy landscapes, diffusion coefficients, and kinetic rates. Compared to existing enhanced sampling methods, we directly exploit short unbiased trajectories at a competitive computational cost.
References in corpus (19)
- Canonical sampling through velocity-rescaling
- Escaping free-energy minima
- PLUMED 2: New feathers for an old bird
- Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
- A Novel Path Sampling Method for the Calculation of Rate Constants
- Sampling rare switching events in biochemical networks
- State-dependent diffusion: thermodynamic consistency and its path integral formulation
- Elaborating Transition Interface Sampling Methods
- Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
- Coarse Molecular Dynamics of a Peptide Fragment: Free Energy, Kinetics, and Long-Time Dynamics Computations
- Molecular hydrodynamics from memory kernels
- Parameter estimation in nonlinear stochastic differential equations
- Pathways to self-organization: crystallization via nucleation and growth
- Likelihood-based non-Markovian models from molecular dynamics
- Exact dynamical coarse-graining without time-scale separation
- Optimal Langevin modelling of out-of-equilibrium molecular dynamics simulations
- Coarse-graining Molecular Systems by Spectral Matching
- Effective equations for reaction coordinates in polymer transport
- Properties of low-dimensional collective variables in the molecular dynamics of biopolymers