Document Type


Date of Original Version





A method is introduced that is easy to implement and greatly reduces the systematic error resulting from quasi-ergodicity, or incomplete sampling of configuration space, in Monte Carlo simulations of systems containing large potential energy barriers. The method makes possible the jumping over these barriers by couplingn the usual Metropolis sampling to the Boltzmann distribution generated by another random walker at a higher temperature. the basic techniques are illustrated on some simple classical systems, beginning for heuristic purposes with a simple one-dimensional double well potential based on a quartic polynomial. the method's suitability for typical multidimensional Monte Carlo systems is demonstrated by extending the double well potential to several dimensions, and then by applying the method to a multiparticle cluster system consisting of argon atoms bound by pairwise Lennard-Jones potentials. Remarkable improvements are demonstrated in the convergence rate for the cluster configuration energy, and especially for the heat capacity, at temperatures near the cluster melting transition region. Moreover, these improvements can be obtained even in the worst-case scenario where clusters are initialized from random configurations.

Publisher Statement

© 1990 American Institute of Physics.



To view the content in your browser, please download Adobe Reader or, alternately,
you may Download the file to your hard drive.

NOTE: The latest versions of Adobe Reader do not support viewing PDF files within Firefox on Mac OS and if you are using a modern (Intel) Mac, there is no official plugin for viewing PDF files within the browser window.