ARTICLES On the acceptance probability of replica-exchange Monte Carlo trials

An analysis is presented of the average probability of accepting an exchange trial in the parallel-tempering Monte Carlo molecular simulation method. Arguments are given that this quantity should be related to the entropy difference between the phases, and results from simulations of a simple Lennard-Jones system are presented to support this argument qualitatively. Another analysis based on the energy distributions of a replica pair is presented, and an exact expression for the trial-move acceptance probability in terms of the overlap of these distributions is derived. A more detailed expression is presented using an approximation of constant heat capacity, and an asymptotic form for this result, good for large system sizes, is reported. The detailed analyses are in quantitative agreement with the simulation data. It is further shown that treatment of the energy distributions as Gaussians is an inappropriate way to analyze the acceptance probability. © 2002 American Institute of Physics. @DOI: 10.1063/1.1507776#

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

[2]  Balint Joo,et al.  Parallel tempering in lattice QCD with O(a)-improved Wilson fermions , 1999 .

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

[4]  D. Bedrov,et al.  Exploration of conformational phase space in polymer melts: A comparison of parallel tempering and conventional molecular dynamics simulations , 2001 .

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

[6]  J. Skolnick,et al.  Parallel-hat tempering: A Monte Carlo search scheme for the identification of low-energy structures , 2001 .

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

[8]  J. Skolnick,et al.  Comparison of three Monte Carlo conformational search strategies for a proteinlike homopolymer model: Folding thermodynamics and identification of low-energy structures , 2000 .

[9]  J. Pablo,et al.  Phase equilibria and clustering in size-asymmetric primitive model electrolytes , 2001 .

[10]  Yuko Okamoto,et al.  Protein folding simulations and structure predictions , 2001 .

[11]  J. Pablo,et al.  Hyperparallel tempering Monte Carlo simulation of polymeric systems , 2000 .

[12]  J. Skolnick,et al.  A new combination of replica exchange Monte Carlo and histogram analysis for protein folding and thermodynamics , 2001 .

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

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

[15]  Y. Sugita,et al.  Ab initio replica-exchange Monte Carlo method for cluster studies , 2001 .

[16]  Erik Sandelin,et al.  Monte Carlo study of the phase structure of compact polymer chains , 1998, cond-mat/9812017.

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

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

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

[20]  Michael W. Deem,et al.  Efficient Monte Carlo methods for cyclic peptides , 1999 .

[21]  B. Dünweg,et al.  Parallel excluded volume tempering for polymer melts. , 2000, Physical review. E, Statistical, nonlinear, and soft matter physics.

[22]  Solid–liquid phase diagram of the water octamer , 2002 .