Path sampling with stochastic dynamics: Some new algorithms

We propose here some new sampling algorithms for path sampling in the case when stochastic dynamics are used. In particular, we present a new proposal function for equilibrium sampling of paths with a Monte-Carlo dynamics (the so-called ''brownian tube'' proposal). This proposal is based on the continuity of the dynamics with respect to the random forcing, and generalizes all previous approaches when stochastic dynamics are used. The efficiency of this proposal is demonstrated using some measure of decorrelation in path space. We also discuss a switching strategy that allows to transform ensemble of paths at a finite rate while remaining at equilibrium, in contrast with the usual Jarzynski like switching. This switching is very interesting to sample constrained paths starting from unconstrained paths, or to perform simulated annealing in a rigorous way.

[1]  C. Jarzynski.  Equilibrium free-energy differences from nonequilibrium measurements: A master-equation approach , 1997, cond-mat/9707325.

[2]  G. Stoltz,et al.  THEORETICAL AND NUMERICAL COMPARISON OF SOME SAMPLING METHODS FOR MOLECULAR DYNAMICS , 2007 .

[3]  G. Stoltz,et al.  Equilibrium Sampling From Nonequilibrium Dynamics , 2005, cond-mat/0511412.

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

[5]  Efficiency analysis of reaction rate calculation methods using analytical models I: The two-dimensional sharp barrier , 2006, cond-mat/0605260.

[6]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[7]  C. Dellago,et al.  Biased sampling of nonequilibrium trajectories: can fast switching simulations outperform conventional free energy calculation methods? , 2005, The journal of physical chemistry. B.

[8]  John E. Boyer,et al.  Variances for adaptive trimmed means , 1983 .

[9]  Thijs J. H. Vlugt,et al.  On the efficient sampling of pathways in the transition path ensemble , 2001 .

[10]  C. Dellago,et al.  Transition Path Sampling , 2005 .

[11]  B. Berne,et al.  Molecular dynamics study of an isomerizing diatomic in a Lennard‐Jones fluid , 1988 .

[12]  Elaborating transition interface sampling methods , 2005, cond-mat/0405116.

[13]  M. Karplus,et al.  Stochastic boundary conditions for molecular dynamics simulations of ST2 water , 1984 .

[14]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[15]  Yukito Iba.  EXTENDED ENSEMBLE MONTE CARLO , 2001 .

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

[17]  D. Zuckerman,et al.  Single-ensemble nonequilibrium path-sampling estimates of free energy differences. , 2004, The Journal of chemical physics.

[18]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[19]  Sean X. Sun.  Equilibrium free energies from path sampling of nonequilibrium trajectories , 2003 .

[20]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[21]  Christoph Dellago,et al.  Transition path sampling: throwing ropes over mountains in the dark , 2000 .

[22]  Richard L. Tweedie,et al.  Markov Chains and Stochastic Stability , 1993, Communications and Control Engineering Series.

[23]  D. Ceperley.  Path integrals in the theory of condensed helium , 1995 .

[24]  J. Rosenthal,et al.  Optimal scaling of discrete approximations to Langevin diffusions , 1998 .

[25]  G. Crooks,et al.  Efficient transition path sampling for nonequilibrium stochastic dynamics. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Christoph Dellago,et al.  Equilibrium Time Correlation Functions from Irreversible Transformations in Trajectory Space , 2004 .

[27]  Caffarel,et al.  Diffusion monte carlo methods with a fixed number of walkers , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[28]  Christoph Dellago,et al.  On the calculation of reaction rate constants in the transition path ensemble , 1999 .

[29]  Mathias Rousset,et al.  On the Control of an Interacting Particle Estimation of Schrödinger Ground States , 2006, SIAM J. Math. Anal..

[30]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[31]  P. Bolhuis,et al.  Sampling the multiple folding mechanisms of Trp-cage in explicit solvent , 2006, Proceedings of the National Academy of Sciences.

[32]  Christoph Dellago,et al.  Activation Energies from Transition Path Sampling Simulations , 2004 .

[33]  D. Frenkel,et al.  Simulating rare events in equilibrium or nonequilibrium stochastic systems. , 2005, The Journal of chemical physics.

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