Identification of simple reaction coordinates from complex dynamics

Abstract

Reaction coordinates are widely used throughout chemical physics to model and understand complex chemical transformations. We introduce a definition of the natural reaction coordinate, suitable for condensed phase and biomolecular systems, as a maximally predictive one-dimensional projection. We then show this criterion is uniquely satisfied by a dominant eigenfunction of an integral operator associated with the ensemble dynamics. We present a new sparse estimator for these eigenfunctions which can search through a large candidate pool of structural order parameters and build simple, interpretable approximations that employ only a small number of these order parameters. Example applications with a small molecule's rotational dynamics and simulations of protein conformational change and folding show that this approach can filter through statistical noise to identify simple reaction coordinates from complex dynamics.

18 pages, 10 figures

References in corpus (1)

• Variational cross-validation of slow dynamical modes in molecular kinetics

Cited by in corpus (27)

• Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems
• Variational Encoding of Complex Dynamics
• Variational Koopman models: slow collective variables and molecular kinetics from short off-equilibrium simulations
• A unified framework for machine learning collective variables for enhanced sampling simulations: $\text{mlcolvar}$
• Transition manifolds of complex metastable systems: Theory and data-driven computation of effective dynamics
• Correlation-based feature selection to identify functional dynamics in proteins
• Finding Efficient Collective Variables: The Case of Crystallization
• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Human learning for molecular simulations: the Collective Variables Dashboard in VMD
• Predictive Collective Variable Discovery with Deep Bayesian Models
• Identification and Analysis of Transition and Metastable Markov States
• Riemann metric approach to optimal sampling of multidimensional free-energy landscapes
• Descriptors-free Collective Variables From Geometric Graph Neural Networks
• Characterizing metastable states with the help of machine learning
• Deep learning the slow modes for rare events sampling
• Deflation reveals dynamical structure in nondominant reaction coordinates
• Inexact iterative numerical linear algebra for neural network-based spectral estimation and rare-event prediction
• Integrated Reaction Path Processing from Sampled Structure Sequences
• Efficiently Sampling Conformations and Pathways Using the Concurrent Adaptive Sampling (CAS) Algorithm
• Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles
• Flow Matching for Optimal Reaction Coordinates of Biomolecular System
• Learning Markovian Dynamics with Spectral Maps
• Transient probability currents provide upper and lower bounds on non-equilibrium steady-state currents in the Smoluchowski picture
• Discovery of slow variables in a class of multiscale stochastic systems via neural networks
• Exploring the free energy gain of phase separation via Markov State Modeling
• Recovering Hidden Degrees of Freedom Using Gaussian Processes
• Deep learning of committor and explainable artificial intelligence analysis for identifying reaction coordinates