Variational Koopman models: slow collective variables and molecular kinetics from short off-equilibrium simulations

Abstract

Markov state models (MSMs) and Master equation models are popular approaches to approximate molecular kinetics, equilibria, metastable states, and reaction coordinates in terms of a state space discretization usually obtained by clustering. Recently, a powerful generalization of MSMs has been introduced, the variational approach (VA) of molecular kinetics and its special case the time-lagged independent component analysis (TICA), which allow us to approximate slow collective variables and molecular kinetics by linear combinations of smooth basis functions or order parameters. While it is known how to estimate MSMs from trajectories whose starting points are not sampled from an equilibrium ensemble, this has not yet been the case for TICA and the VA. Previous estimates from short trajectories, have been strongly biased and thus not variationally optimal. Here, we employ Koopman operator theory and ideas from dynamic mode decomposition (DMD) to extend the VA and TICA to non-equilibrium data. The main insight is that the VA and TICA provide a coefficient matrix that we call Koopman model, as it approximates the underlying dynamical (Koopman) operator in conjunction with the basis set used. This Koopman model can be used to compute a stationary vector to reweight the data to equilibrium. From such a Koopman-reweighted sample, equilibrium expectation values and variationally optimal reversible Koopman models can be constructed even with short simulations. The Koopman model can be used to propagate densities, and its eigenvalue decomposition provide estimates of relaxation timescales and slow collective variables for dimension reduction. Koopman models are generalizations of Markov state models, TICA and the linear VA and allow molecular kinetics to be described without a cluster discretization.

References in corpus (1)

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

Cited by in corpus (38)

• Machine learning for molecular simulation
• VAMPnets: Deep learning of molecular kinetics
• Extended dynamic mode decomposition with dictionary learning: a data-driven adaptive spectral decomposition of the Koopman operator
• Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
• On Convergence of Extended Dynamic Mode Decomposition to the Koopman Operator
• Data-driven model reduction and transfer operator approximation
• Machine learning force fields and coarse-grained variables in molecular dynamics: application to materials and biological systems
• Deep Neural Networks with Koopman Operators for Modeling and Control of Autonomous Vehicles
• Variational Encoding of Complex Dynamics
• A variational conformational dynamics approach to the selection of collective variables in metadynamics
• Galerkin Approximation of Dynamical Quantities using Trajectory Data
• Water regulates the residence time of Benzamidine in Trypsin
• Variational Selection of Features for Molecular Kinetics
• GraphVAMPNet, using graph neural networks and variational approach to markov processes for dynamical modeling of biomolecules
• Enhanced Sampling in the Age of Machine Learning: Algorithms and Applications
• Variational approach for learning Markov processes from time series data
• Integrated VAC: A robust strategy for identifying eigenfunctions of dynamical operators
• Identification of kinetic order parameters for non-equilibrium dynamics
• Quantifying the dynamics of protein self-organization using deep learning analysis of atomic force microscopy data
• Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics
• Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology
• Deep learning Markov and Koopman models with physical constraints
• Manifold Learning in Atomistic Simulations: A Conceptual Review
• Deep learning the slow modes for rare events sampling
• Deflation reveals dynamical structure in nondominant reaction coordinates
• Computing transition path theory quantities with trajectory stratification
• Elucidating Interfacial Dynamics of Ti-Al Systems Using Molecular Dynamics Simulation and Markov State Modeling
• Simultaneous Coherent Structure Coloring facilitates interpretable clustering of scientific data by amplifying dissimilarity
• Accurate estimates of dynamical statistics using memory
• Flow Matching for Optimal Reaction Coordinates of Biomolecular System
• Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles
• Error bounds for dynamical spectral estimation
• Learning Markovian Dynamics with Spectral Maps
• Physics-aware, probabilistic model order reduction with guaranteed stability
• Kernel Embedding based Variational Approach for Low-dimensional Approximation of Dynamical Systems
• Dimensional Reduction of Dynamical Systems by Machine Learning: Automatic Generation of the Optimum Extensive Variables and Their Time-Evolution Map
• Discovery of Quasi-Integrable Equations from traveling-wave data using the Physics-Informed Neural Networks
• Improving Estimation of the Koopman Operator with Kolmogorov-Smirnov Indicator Functions