Deep learning the slow modes for rare events sampling

Abstract

The development of enhanced sampling methods has greatly extended the scope of atomistic simulations, allowing long-time phenomena to be studied with accessible computational resources. Many such methods rely on the identification of an appropriate set of collective variables. These are meant to describe the system's modes that most slowly approach equilibrium. Once identified, the equilibration of these modes is accelerated by the enhanced sampling method of choice. An attractive way of determining the collective variables is to relate them to the eigenfunctions and eigenvalues of the transfer operator. Unfortunately, this requires knowing the long-term dynamics of the system beforehand, which is generally not available. However, we have recently shown that it is indeed possible to determine efficient collective variables starting from biased simulations. In this paper, we bring the power of machine learning and the efficiency of the recently developed on-the-fly probability enhanced sampling method to bear on this approach. The result is a powerful and robust algorithm that, given an initial enhanced sampling simulation performed with trial collective variables or generalized ensembles, extracts transfer operator eigenfunctions using a neural network ansatz and then accelerates them to promote sampling of rare events. To illustrate the generality of this approach we apply it to several systems, ranging from the conformational transition of a small molecule to the folding of a mini-protein and the study of materials crystallization.

References in corpus (29)

• 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
• Identification of slow molecular order parameters for Markov model construction
• VAMPnets: Deep learning of molecular kinetics
• Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
• Rethinking Metadynamics: from bias potentials to probability distributions
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• 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
• Data-Driven Collective Variables for Enhanced Sampling
• Variational Encoding of Complex Dynamics
• Nonlinear Discovery of Slow Molecular Modes using State-Free Reversible VAMPnets
• Silicon liquid structure and crystal nucleation from ab-initio deep Metadynamics
• Variational Koopman models: slow collective variables and molecular kinetics from short off-equilibrium simulations
• A Unified Approach to Enhanced Sampling
• Neural networks-based variationally enhanced sampling
• A variational conformational dynamics approach to the selection of collective variables in metadynamics
• Identification of simple reaction coordinates from complex dynamics
• The role of water in host-guest interaction
• Calculation of phase diagrams in the multithermal-multibaric ensemble
• Folding a Small Protein Using Harmonic Linear Discriminant Analysis
• Gaussian Mixture Based Enhanced Sampling
• Unfolding Hidden Barriers by Active Enhanced Sampling
• Finding Efficient Collective Variables: The Case of Crystallization
• Collective Variables for the Study of Crystallization
• Making the best of a bad situation: a multiscale approach to free energy calculation
• Multiscale reweighted stochastic embedding (MRSE): Deep learning of collective variables for enhanced sampling