Reweighting non-equilibrium steady-state dynamics along collective variables

Abstract

Computer simulations generate microscopic trajectories of complex systems at a single thermodynamic state point. We recently introduced a Maximum Caliber (MaxCal) approach for dynamical reweighting. Our approach mapped these trajectories to a Markovian description on the configurational coordinates, and reweighted path probabilities as a function of external forces. Trajectory probabilities can be dynamically reweighted both from and to equilibrium or non-equilibrium steady states. As the system's dimensionality increases, an exhaustive description of the microtrajectories becomes prohibitive--even with a Markovian assumption. Instead we reduce the dimensionality of the configurational space to collective variables (CVs). Going from configurational to CV space, we define local entropy productions derived from configurationally averaged mean forces. The entropy production is shown to be a suitable constraint on MaxCal for non-equilibrium steady states expressed as a function of CVs. We test the reweighting procedure on two systems: a particle subject to a two-dimensional potential and a coarse-grained peptide. Our CV-based MaxCal approach expands dynamical reweighting to larger systems, for both static and dynamical properties, and across a large range of driving forces.

12 pages, 7 figures

References in corpus (5)

• Maximum Caliber: a general variational principle for dynamical systems
• Special Topic: Markov Models of Molecular Kinetics
• xTRAM: Estimating equilibrium expectations from time-correlated simulation data at multiple thermodynamic states
• Concurrent parametrization against static and kinetic information leads to more robust coarse-grained force fields
• Free-energy landscapes in magnetic systems from metadynamics