Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states

Abstract

We propose a discrete transition-based reweighting analysis method (dTRAM) for analyzing configuration-space-discretized simulation trajectories produced at different thermodynamic states (temperatures, Hamiltonians, etc.) dTRAM provides maximum-likelihood estimates of stationary quantities (probabilities, free energies, expectation values) at any thermodynamic state. In contrast to the weighted histogram analysis method (WHAM), dTRAM does not require data to be sampled from global equilibrium, and can thus produce superior estimates for enhanced sampling data such as parallel/simulated tempering, replica exchange, umbrella sampling, or metadynamics. In addition, dTRAM provides optimal estimates of Markov state models (MSMs) from the discretized state-space trajectories at all thermodynamic states. Under suitable conditions, these MSMs can be used to calculate kinetic quantities (e.g. rates, timescales). In the limit of a single thermodynamic state, dTRAM estimates a maximum likelihood reversible MSM, while in the limit of uncorrelated sampling data, dTRAM is identical to WHAM. dTRAM is thus a generalization to both estimators.

References in corpus (7)

• Escaping free-energy minima
• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Statistically optimal analysis of samples from multiple equilibrium states
• Simulated Tempering: A New Monte Carlo Scheme
• Parallel Tempering Algorithm for Conformational Studies of Biological Molecules
• Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages
• xTRAM: Estimating equilibrium expectations from time-correlated simulation data at multiple thermodynamic states

Cited by in corpus (14)

• VAMPnets: Deep learning of molecular kinetics
• Time-lagged autoencoders: Deep learning of slow collective variables for molecular kinetics
• Special Topic: Markov Models of Molecular Kinetics
• Adaptive Markov State Model estimation using short reseeding trajectories
• Girsanov reweighting for path ensembles and Markov state models
• Particle-based membrane model for mesoscopic simulation of cellular dynamics
• Dynamical reweighting methods for Markov models
• MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations
• Efficient estimation of rare-event kinetics
• A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models
• Path probability ratios for Langevin dynamics -- exact and approximate
• An efficient multi-scale Green's Functions Reaction Dynamics scheme
• Optimization of Non-Equilibrium Self-Assembly Protocols Using Markov State Models
• Microscopic reweighting for non-equilibrium steady states dynamics