Optimal temperature ladders in replica exchange simulations

In replica exchange simulations, a temperature ladder with N rungs spans a given temperature interval. Considering systems with heat capacities independent of the temperature, here we address the question of how large N should be chosen for an optimally fast diffusion of the replicas through the temperature space. Using a simple example we show that choosing average acceptance probabilities of about 45% and computing N accordingly maximizes the round trip rates r across the given temperature range. This result differs from previous analyses which suggested smaller average acceptance probabilities of about 23%. We show that the latter choice maximizes the ratio r/N instead of r.

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

[2]  Ulrich H E Hansmann,et al.  Folding proteins by first-passage-times-optimized replica exchange. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  W. Nadler,et al.  Optimized explicit-solvent replica exchange molecular dynamics from scratch. , 2008, The journal of physical chemistry. B.

[4]  Y. Okamoto,et al.  Alpha-helix folding by Monte Carlo simulated annealing in isolated C-peptide of ribonuclease A. , 1991, Protein engineering.

[5]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[6]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[7]  H. Nymeyer,et al.  Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water. , 2007, Journal of structural biology.

[8]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[9]  C. Predescu,et al.  The incomplete beta function law for parallel tempering sampling of classical canonical systems. , 2003, The Journal of chemical physics.

[10]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[11]  M. Troyer,et al.  Optimized parallel tempering simulations of proteins. , 2006, The Journal of chemical physics.

[12]  D. Kofke,et al.  Selection of temperature intervals for parallel-tempering simulations. , 2005, The Journal of chemical physics.

[13]  Ulrich H E Hansmann,et al.  Generalized ensemble and tempering simulations: a unified view. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  J. D. de Pablo,et al.  Optimal allocation of replicas in parallel tempering simulations. , 2005, The Journal of chemical physics.

[15]  C. Zhang,et al.  Comparison of sampling efficiency between simulated tempering and replica exchange. , 2008, The Journal of chemical physics.

[16]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.