All-exchanges parallel tempering.

An alternative exchange strategy for parallel tempering simulations is introduced. Instead of attempting to swap configurations between two randomly chosen but adjacent replicas, the acceptance probabilities of all possible swap moves are calculated a priori. One specific swap move is then selected according to its probability and enforced. The efficiency of the method is illustrated first on the case of two Lennard-Jones (LJ) clusters containing 13 and 31 atoms, respectively. The convergence of the caloric curve is seen to be at least twice as fast as in conventional parallel tempering simulations, especially for the difficult case of LJ31. Further evidence for an improved efficiency is reported on the ergodic measure introduced by Mountain and Thirumalai [J. Phys. Chem. 93, 6975 (1989)], calculated here for LJ13 close to the melting point. Finally, tests on two simple spin systems indicate that the method should be particularly useful when a limited number of replicas are available.

[1]  Doll,et al.  Approach to ergodicity in monte carlo simulations , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[2]  Straub,et al.  Generalized simulated annealing algorithms using Tsallis statistics: Application to conformational optimization of a tetrapeptide. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[4]  Rehberg,et al.  Simulated-tempering procedure for spin-glass simulations. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[5]  C. Predescu,et al.  On the efficiency of exchange in parallel tempering monte carlo simulations. , 2004, The journal of physical chemistry. B.

[6]  D. L. Freeman,et al.  Phase changes in 38-atom Lennard-Jones clusters. I. A parallel tempering study in the canonical ensemble , 2000, physics/0003068.

[7]  David A. Kofke,et al.  ARTICLES On the acceptance probability of replica-exchange Monte Carlo trials , 2002 .

[8]  Alan M. Ferrenberg,et al.  Critical behavior of the three-dimensional Ising model: A high-resolution Monte Carlo study. , 1991, Physical review. B, Condensed matter.

[9]  Sanford Weisberg,et al.  Computing science and statistics : proceedings of the 30th Symposium on the Interface, Minneapolis, Minnesota, May 13-16, 1998 : dimension reduction, computational complexity and information , 1998 .

[10]  J. D. Doll,et al.  A heat capacity estimator for Fourier path integral simulations , 2000 .

[11]  D. L. Freeman,et al.  Phase changes in 38-atom Lennard-Jones clusters. II. A parallel tempering study of equilibrium and dynamic properties in the molecular dynamics and microcanonical ensembles , 2000, physics/0003072.

[12]  J. Pablo,et al.  Hyper-parallel tempering Monte Carlo: Application to the Lennard-Jones fluid and the restricted primitive model , 1999 .

[13]  Thirumalai,et al.  Ergodic convergence properties of supercooled liquids and glasses. , 1990, Physical review. A, Atomic, molecular, and optical physics.

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

[15]  Yamamoto,et al.  Replica-exchange molecular dynamics simulation for supercooled liquids , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[16]  J. Doye,et al.  Evolution of the Potential Energy Surface with Size for Lennard-Jones Clusters , 1999, cond-mat/9903305.

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

[18]  F. Calvo,et al.  On the convergence of global thermal properties of clusters extracted from simulations , 2000 .

[19]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[20]  A. Young,et al.  Spin glasses and random fields , 1997 .

[21]  J. Pablo,et al.  Multicanonical parallel tempering , 2002, cond-mat/0201179.

[22]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[23]  L. Piela,et al.  Molecular Dynamics on Deformed Potential Energy Hypersurfaces , 1995 .

[24]  Wolff,et al.  Collective Monte Carlo updating for spin systems. , 1989, Physical review letters.

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

[26]  Jianpeng Ma,et al.  Simulated annealing using the classical density distribution , 1994 .

[27]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[28]  K. Jordan,et al.  On the convergence of parallel tempering Monte Carlo simulations of LJ38. , 2005, The journal of physical chemistry. A.

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

[30]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[31]  Jonathan P. K. Doye,et al.  Entropic tempering: A method for overcoming quasiergodicity in simulation , 2000 .

[32]  J. D. Doll,et al.  Heat Capacity Estimators for Random Series Path-Integral Methods by Finite-Difference Schemes , 2003, cond-mat/0307769.

[33]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[34]  F. Y. Wu.  The Potts model , 1982 .

[35]  I. Andricioaei,et al.  On Monte Carlo and molecular dynamics methods inspired by Tsallis statistics: Methodology, optimization, and application to atomic clusters , 1997 .

[36]  K. Sanbonmatsu,et al.  Structure of Met‐enkephalin in explicit aqueous solution using replica exchange molecular dynamics , 2002, Proteins.

[37]  R. Whetten,et al.  Statistical thermodynamics of the cluster solid-liquid transition. , 1990, Physical review letters.

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

[39]  M. Deem,et al.  A biased Monte Carlo scheme for zeolite structure solution , 1998, cond-mat/9809085.

[40]  D. L. Freeman,et al.  Reducing Quasi-Ergodic Behavior in Monte Carlo Simulations by J-Walking: Applications to Atomic Clusters , 1990 .

[41]  Wang,et al.  Nonuniversal critical dynamics in Monte Carlo simulations. , 1987, Physical review letters.

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

[43]  D. D. Frantz.  Magic number behavior for heat capacities of medium-sized classical Lennard-Jones clusters , 2001, physics/0107013.

[44]  D. Kofke.  Comment on "The incomplete beta function law for parallel tempering sampling of classical canonical systems" [J. Chem. Phys. 120, 4119 (2004)]. , 2004, The Journal of chemical physics.

[45]  L. Fried,et al.  Improved heat capacity estimator for path integral simulations , 2002 .

[46]  D Frenkel,et al.  Designing refoldable model molecules. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[48]  Raymond D. Mountain,et al.  Measures of effective ergodic convergence in liquids , 1989 .