The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures.

The "weighted ensemble" method, introduced by Huber and Kim [Biophys. J. 70, 97 (1996)], is one of a handful of rigorous approaches to path sampling of rare events. Expanding earlier discussions, we show that the technique is statistically exact for a wide class of Markovian and non-Markovian dynamics. The derivation is based on standard path-integral (path probability) ideas, but recasts the weighted-ensemble approach as simple "resampling" in path space. Similar reasoning indicates that arbitrary nonstatic binning procedures, which merely guide the resampling process, are also valid. Numerical examples confirm the claims, including the use of bins which can adaptively find the target state in a simple model.

[1]  Martin J Field,et al.  Transition events in one dimension. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Fernando A Escobedo,et al.  Optimizing the sampling and staging for simulations of rare events via forward flux sampling schemes. , 2008, The Journal of chemical physics.

[3]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[4]  Jun S. Liu,et al.  Monte Carlo strategies in scientific computing , 2001 .

[5]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[6]  N. Madras,et al.  THE SELF-AVOIDING WALK , 2006 .

[7]  L Y Chen,et al.  Efficient transition path sampling for systems with multiple reaction pathways. , 2005, The Journal of chemical physics.

[8]  Bin W. Zhang,et al.  Efficient and verified simulation of a path ensemble for conformational change in a united-residue model of calmodulin , 2007, Proceedings of the National Academy of Sciences.

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

[10]  Franz Aurenhammer,et al.  Voronoi diagrams—a survey of a fundamental geometric data structure , 1991, CSUR.

[11]  A. Dinner,et al.  Nonequilibrium umbrella sampling in spaces of many order parameters. , 2009, The Journal of chemical physics.

[12]  D. Sherrington.  Stochastic Processes in Physics and Chemistry , 1983 .

[13]  P. Grassberger.  Go with the Winners: a General Monte Carlo Strategy , 2002, cond-mat/0201313.

[14]  Lawrence R. Pratt,et al.  A statistical method for identifying transition states in high dimensional problems , 1986 .

[15]  D. Wales.  Discrete path sampling , 2002 .

[16]  F. T. Wall,et al.  New Method for the Statistical Computation of Polymer Dimensions , 1959 .

[17]  P. R. ten Wolde,et al.  Sampling rare switching events in biochemical networks. , 2004, Physical review letters.

[18]  D. Frenkel,et al.  Understanding molecular simulation : from algorithms to applications. 2nd ed. , 2002 .

[19]  Aaron R Dinner,et al.  Umbrella sampling for nonequilibrium processes. , 2007, The Journal of chemical physics.

[20]  Roy,et al.  Fast, accurate algorithm for numerical simulation of exponentially correlated colored noise. , 1988, Physical review. A, General physics.

[21]  Sebastian Doniach,et al.  Simulation of protein folding by reaction path annealing , 2001 .

[22]  Zheng Yang,et al.  Allosteric Transitions of Supramolecular Systems Explored by Network Models: Application to Chaperonin GroEL , 2009, PLoS Comput. Biol..

[23]  A Rojnuckarin,et al.  Brownian dynamics simulations of protein folding: access to milliseconds time scale and beyond. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Dieter Forster,et al.  Hydrodynamic fluctuations, broken symmetry, and correlation functions , 1975 .

[25]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

[26]  Stein,et al.  Optimal paths and the prehistory problem for large fluctuations in noise-driven systems. , 1992, Physical review letters.

[27]  Peter G. Bolhuis,et al.  A novel path sampling method for the calculation of rate constants , 2003 .

[28]  J. Araque,et al.  Transition path sampling and forward flux sampling. Applications to biological systems , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[29]  T. Woolf,et al.  DYNAMIC REACTION PATHS AND RATES THROUGH IMPORTANCE-SAMPLED STOCHASTIC DYNAMICS , 1999 .

[30]  P. Grassberger.  Pruned-enriched Rosenbluth method: Simulations of θ polymers of chain length up to 1 000 000 , 1997 .

[31]  Daniel M Zuckerman,et al.  Transition-event durations in one-dimensional activated processes. , 2007, The Journal of chemical physics.

[32]  K. Schulten,et al.  Molecular dynamics study of unbinding of the avidin-biotin complex. , 1997, Biophysical journal.

[33]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[34]  Ernesto E. Borrero,et al.  Reaction coordinates and transition pathways of rare events via forward flux sampling. , 2007, The Journal of chemical physics.

[35]  Jürgen Schlitter,et al.  Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin , 1993 .

[36]  G. Huber,et al.  Weighted-ensemble Brownian dynamics simulations for protein association reactions. , 1996, Biophysical journal.