Web ensemble averages for retrieving relevant information from rejected Monte Carlo moves

Abstract.We study the relevance of including information about rejected Monte-Carlo moves in path-sampling computations of free energies. For this purpose, we define webs as sets of paths linked by the path-sampling scheme and introduce an associated statistical ensemble. Within this web ensemble, we derive and test several statistical averages enabling to include information about configurational and path quantities belonging to the unselected trial moves. We numerically observe that retrieving this information does not always result in variance reduction, as theoretically predicted by Delmas and Jourdain. To explain the possible detrimental effect of information-retrieving from web sampling, an action for the webs is introduced. The behaviour of the statistical variance is observed to correlate to an overlapping area of a web action histogram. This area represents the probability that a generated web is such that the difference of its action between the targeted and reference ensembles is lower than the corresponding difference of free energy. Variance reductions are numerically observed for increased areas, as it is the case for the residence weight method proposed previously. More generally, web ensembles yield a rigorous framework for rationalizing existing methods and also for deriving potentially new methods enabling to retrieve relevant information from rejected trial moves.

[1]  R. Swendsen,et al.  Comparison of free energy methods for molecular systems. , 2006, The Journal of chemical physics.

[2]  D. Kofke,et al.  Rosenbluth-sampled nonequilibrium work method for calculation of free energies in molecular simulation. , 2005, The Journal of chemical physics.

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

[4]  Daan Frenkel,et al.  Virtual-move parallel tempering. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[5]  D. R. Mason,et al.  Stochastic kinetic Monte Carlo algorithms for long-range Hamiltonians , 2004, Comput. Phys. Commun..

[6]  Daan Frenkel,et al.  Monte Carlo Sampling of a Markov Web. , 2005, Journal of chemical theory and computation.

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

[8]  M Athènes Parallel Monte Carlo simulations using a residence weight algorithm. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[10]  M. Athènes,et al.  Gibbs free-energy estimates from direct path-sampling computations. , 2005, The Journal of chemical physics.

[11]  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.

[12]  A. Hauet,et al.  Ordering and precipitation of γ′ phase in low supersaturated Ni-Cr-Al model alloy: An atomic scale investigation , 1997 .

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

[14]  B. Jourdain,et al.  Does waste-recycling really improve Metropolis-Hastings Monte Carlo algorithm? , 2006, math/0611949.

[15]  J. D. de Pablo,et al.  Improved density of states Monte Carlo method based on recycling of rejected states. , 2006, The Journal of chemical physics.

[16]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

[17]  D. Frenkel Speed-up of Monte Carlo simulations by sampling of rejected states. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Christoph Dellago,et al.  On the efficiency of path sampling methods for the calculation of free energies from non-equilibrium simulations , 2007 .

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

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

[21]  Pascal Bellon,et al.  Identification of novel diffusion cycles in B2 ordered phases by Monte Carlo simulation , 1997 .

[22]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[23]  M. Athènes A path-sampling scheme for computing thermodynamic properties of a many-body system in a generalized ensemble , 2004 .

[24]  M. Athènes Computation of a chemical potential using a residence weight algorithm. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[25]  Jeanne-Marie Lanore,et al.  Simulation de l'evolution des defauts dans un reseau par le methode de monte-carlo , 1974 .

[26]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[27]  G. Adjanor,et al.  Free energy landscape from path-sampling: application to the structural transition in LJ38 , 2006 .

[28]  Giovanni Ciccotti,et al.  Computer Simulations in Condensed Matter Systems: From Materials to Chemical Biology. Volume 2 (Lecture Notes in Physics) , 2006 .

[29]  Thomas B Woolf,et al.  Using overlap and funnel sampling to obtain accurate free energies from nonequilibrium work measurements. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.