- Access by Xinjiang University
Extended state-space Monte Carlo methods
Phys. Rev. E 63, 056701 – Published 13 April, 2001
DOI: https://doi.org/10.1103/PhysRevE.63.056701
Abstract
In this paper various extensions of the parallel-tempering algorithm are developed and their properties are analyzed. The algorithms are designed to alleviate quasiergodic sampling in systems which have rough energy landscapes by coupling individual Monte Carlo chains to form a composite chain. As with parallel tempering, the procedures are based upon extending the state space to include parameters to encourage sampling mobility. One of the drawbacks of the parallel-tempering method is the stochastic nature of the Monte Carlo dynamics in the auxiliary variables which extend the state space. In this work, the possibility of improving the sampling rate by designing deterministic methods of moving through the parameter space is investigated. The methods developed in this article, which are based upon a statistical quenching and heating procedure similar in spirit to simulated annealing, are tested on a simple two-dimensional spin system model) and on a model in vacuo polypeptide system. In the coupled Monte Carlo chain algorithms, we find that the net mobility of the composite chain is determined by the competition between the characteristic time of coupling between adjacent chains and the degree of overlap of their distributions. Extensive studies of all methods are carried out to obtain optimal sampling conditions. In particular, the most efficient parallel-tempering procedure is to attempt to swap configurations after very few Monte Carlo updates of the composite chains. Furthermore, it is demonstrated that, contrary to expectations, the deterministic procedure does not improve the sampling rate over that of parallel tempering.
References (20)
- N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, J. Chem. Phys. 21, 1087 (1953).
- See, for example, K. Binder and D. W. Heermann, Monte Carlo Simulation in Statistical Physics (Spring-Verlag, Berlin, 1992).
- J. P. Valleau and S. G. Whittington, in Statistical Mechanics, edited by B. J. Berne (Plenum, New York, 1977), Chap. 4, pp. 114–119.
- B. A. Berg and T. Neuhaus, Phys. Lett. B 267, 249 (1991); Phys. Rev. Lett. 68, 9 (1992); J. P. Valleau, J. Chem. Phys. 99, 4718 (1993).
- W. Kerler and P. Rehberg, Phys. Rev. E 50, 4220 (1994); B. A. Berg and T. Celik, Phys. Rev. Lett. 69, 2292 (1992).
- U. H. E. Hansmann and Y. Okamoto, J. Comput. Chem. 14, 1333 (1993); M.-H. Hao and H. A. Scheraga, J. Phys. Chem. 98, 4940 (1994).
- G. M. Torrie and J. P. Valleau, J. Comput. Phys. 23, 187 (1977); J. P. Valleau, J. Chem. Phys. 99, 4718 (1993).
- B. A. Berg, Int. J. Mod. Phys. C 3, 1083 (1992).
- E. Marinari and G. Parisi, Europhys. Lett. 19, 451 (1992); A. P. Lyubartsev, A. A. Martinovski, S. V. Shevkunov, and P. N. Vorontsov-Velyminov, J. Chem. Phys. 96, 1776 (1992).
- S. Kirkpatrick, C. D. Gelatt, and P. M. Vecchi, Science 220, 671 (1983).
- C. J. Geyer, in Computing Science and Statistics: Proceedings of the 23rd Symposium on the Interface (Interface, Fairfax Station, VA, 1991), p. 156.
- R. M. Neal (unpublished).
- C. J. Geyer and E. A. Thompson, J. Am. Stat. Assoc. 90, 909 (1995).
- M. C. Tesi, E. J. Janse van Rensberg, E. Orlandini, and S. G. Whittington, J. Stat. Phys. 82, 155 (1996).
- U. H. E. Hansmann, Chem. Phys. Lett. 281, 140 (1997).
- K. Hukushima and K. Nemoto, J. Phys. Soc. Jpn. 65, 1604 (1996).
- R. M. Neal, Stat. Comput. 6, 353 (1996).
- QUANTA Parameter Handbook (Polygen Corp., Waltham, MA, 1990).
- D. D. Frantz, D. L. Freeman, and J. D. Doll, J. Chem. Phys. 93, 2769 (1990).
- R. Iftimie, D. Salahub, D. Wei, and J. Schofield, J. Chem. Phys. 113, 4852 (2000).