Estimating equilibrium ensemble averages using multiple time slices from driven nonequilibrium processes: theory and application to free energies, moments, and thermodynamic length in single-molecule pulling experiments.

Recently discovered identities in statistical mechanics have enabled the calculation of equilibrium ensemble averages from realizations of driven nonequilibrium processes, including single-molecule pulling experiments and analogous computer simulations. Challenges in collecting large data sets motivate the pursuit of efficient statistical estimators that maximize use of available information. Along these lines, Hummer and Szabo developed an estimator that combines data from multiple time slices along a driven nonequilibrium process to compute the potential of mean force. Here, we generalize their approach, pooling information from multiple time slices to estimate arbitrary equilibrium expectations. Our expression may be combined with estimators of path-ensemble averages, including existing optimal estimators that use data collected by unidirectional and bidirectional protocols. We demonstrate the estimator by calculating free energies, moments of the polymer extension, the thermodynamic metric tensor, and the thermodynamic length in a model single-molecule pulling experiment. Compared to estimators that only use individual time slices, our multiple time-slice estimators yield substantially smoother estimates and achieve lower variance for higher-order moments.

[1]  G. Crooks.  Nonequilibrium Measurements of Free Energy Differences for Microscopically Reversible Markovian Systems , 1998 .

[2]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[4]  David D L Minh,et al.  Optimized free energies from bidirectional single-molecule force spectroscopy. , 2008, Physical review letters.

[5]  C. Jarzynski.  Rare events and the convergence of exponentially averaged work values. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  F. Ritort,et al.  Recovery of free energy branches in single molecule experiments. , 2009, Physical review letters.

[7]  Felix Ritort,et al.  Measurement of work in single-molecule pulling experiments. , 2009, The Journal of chemical physics.

[8]  Stefan Schaefer,et al.  Speeding up parallel tempering simulations. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  C. Jarzynski,et al.  Verification of the Crooks fluctuation theorem and recovery of RNA folding free energies , 2005, Nature.

[10]  Peter Salamon,et al.  Thermodynamic length and dissipated availability , 1983 .

[11]  Ned S. Wingreen,et al.  Self-Organization of the Escherichia coli Chemotaxis Network Imaged with Super-Resolution Light Microscopy , 2009, PLoS biology.

[12]  Peter Salamon,et al.  The geometry of separation processes: A horse-carrot theorem for steady flow systems , 1998 .

[13]  Leonidas J. Guibas,et al.  Optimally combining sampling techniques for Monte Carlo rendering , 1995, SIGGRAPH.

[14]  H. Then,et al.  Measuring the convergence of Monte Carlo free-energy calculations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Manuel Athènes,et al.  Free energy reconstruction from steered dynamics without post-processing , 2010, J. Comput. Phys..

[16]  Gavin E Crooks,et al.  Measuring thermodynamic length. , 2007, Physical review letters.

[17]  Leonard M. Sander,et al.  Scaling and river networks: A Landau theory for erosion , 1997 .

[18]  Andrew J Ballard,et al.  Replica exchange with nonequilibrium switches , 2009, Proceedings of the National Academy of Sciences.

[19]  G. Crooks.  Path-ensemble averages in systems driven far from equilibrium , 1999, cond-mat/9908420.

[20]  Radford M. Neal.  Annealed importance sampling , 1998, Stat. Comput..

[21]  Rory A. Fisher,et al.  Theory of Statistical Estimation , 1925, Mathematical Proceedings of the Cambridge Philosophical Society.

[22]  G. Voth,et al.  Extending the fluctuation theorem to describe reaction coordinates. , 2007, The Journal of chemical physics.

[23]  Michael R. Shirts,et al.  Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. , 2003, Physical review letters.

[24]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[25]  Christoph Dellago,et al.  Efficient extraction of free energy profiles from nonequilibrium experiments , 2009, J. Comput. Chem..

[26]  I. Tinoco,et al.  Equilibrium Information from Nonequilibrium Measurements in an Experimental Test of Jarzynski's Equality , 2002, Science.

[27]  Jianhua Lin,et al.  Divergence measures based on the Shannon entropy , 1991, IEEE Trans. Inf. Theory.

[28]  Peter A. Kollman,et al.  The lag between the Hamiltonian and the system configuration in free energy perturbation calculations , 1989 .

[29]  Leonidas J. Guibas,et al.  Robust Monte Carlo methods for light transport simulation , 1997 .

[30]  David D. L. Minh,et al.  Optimal estimators and asymptotic variances for nonequilibrium path-ensemble averages. , 2009, The Journal of chemical physics.