Mapping the dynamics of multi-dimensional systems onto a nearest-neighbor coupled discrete set of states conserving the mean first-passage times: a projective dynamics approach

We consider classical and semi-classical dynamical systems that start from a given ensemble of configurations and evolve in time until the systems reach a certain fixed stopping criterion, with the mean first-passage time (MFPT) being the quantity of interest. We present a method, projective dynamics, which maps the dynamics of the system onto an arbitrary discrete set of states {ζk}, subject to the constraint that the states ζk are chosen in such a way that only transitions not further than to the neighboring states ζk ± 1 occur. We show that with this imposed condition there exists a master equation with nearest-neighbor coupling with the same MFPT and residence times as the original dynamical system. We show applications of the method for the diffusion process of particles in one- and two-dimensional potential energy landscapes and the folding process of a small biopolymer. We compare results for the MFPT and the mean folding time obtained with the projective dynamics method with those obtained by a direct measurement, and where possible with a semi-analytical solution.

[1]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[2]  A. Fersht.  On the simulation of protein folding by short time scale molecular dynamics and distributed computing , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Olof Görnerup,et al.  A Method for Inferring Hierarchical Dynamics in Stochastic Processes , 2007, Adv. Complex Syst..

[4]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.

[5]  Slow forcing in the projective dynamics method , 1998, cond-mat/9811039.

[6]  M. Karplus,et al.  Hidden complexity of free energy surfaces for peptide (protein) folding. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[7]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[8]  M. Kolesik,et al.  Projection Method for Statics and Dynamics of Lattice Spin Systems , 1998 .

[9]  Miroslav Kolesik,et al.  EXTREME LONG-TIME DYNAMIC MONTE CARLO SIMULATIONS FOR METASTABLE DECAY IN THE d=3 ISING FERROMAGNET , 2002, cond-mat/0207405.

[10]  Sidney Redner,et al.  A guide to first-passage processes , 2001 .

[11]  J. Klafter,et al.  First-passage times in complex scale-invariant media , 2007, Nature.

[12]  M. Novotny,et al.  A Tutorial on Advanced Dynamic Monte Carlo Methods for Systems with Discrete State Spaces , 2001, cond-mat/0109182.

[13]  Martin Karplus,et al.  Aspects of Protein Reaction Dynamics: Deviations from Simple Behavior , 2000 .

[14]  R. Zwanzig,et al.  Diffusion in a rough potential. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[15]  S Banu Ozkan,et al.  Computing the transition state populations in simple protein models. , 2003, Biopolymers.

[16]  R. Elber,et al.  Computing time scales from reaction coordinates by milestoning. , 2004, The Journal of chemical physics.

[17]  M. Novotny,et al.  Advances in the Projective Dynamics Method: A Procedure of Discretizing the Space applied to Markovian Processes , 2010 .

[18]  K. Dill,et al.  Protein folding in the landscape perspective: Chevron plots and non‐arrhenius kinetics , 1998, Proteins.