Stochastic potential switching algorithm for Monte Carlo simulations of complex systems.

This paper describes a new Monte Carlo method based on a novel stochastic potential switching algorithm. This algorithm enables the equilibrium properties of a system with potential V to be computed using a Monte Carlo simulation for a system with a possibly less complex stochastically altered potential V. By proper choices of the stochastic switching and transition probabilities, it is shown that detailed balance can be strictly maintained with respect to the original potential V. The validity of the method is illustrated with a simple one-dimensional example. The method is then generalized to multidimensional systems with any additive potential, providing a framework for the design of more efficient algorithms to simulate complex systems. A near-critical Lennard-Jones fluid with more than 20,000 particles is used to illustrate the method. The new algorithm produced a much smaller dynamic scaling exponent compared to the Metropolis method and improved sampling efficiency by over an order of magnitude.

[1]  Bruce J. Berne,et al.  Multicanonical jump walking: A method for efficiently sampling rough energy landscapes , 1999 .

[2]  G Stolovitzky,et al.  Catalytic tempering: A method for sampling rough energy landscapes by Monte Carlo. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[3]  C. H. Mak,et al.  A Multigrid Algorithm for Sampling Imaginary-Time Paths in Quantum Monte Carlo Simulations† , 2004 .

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

[5]  Jimmie D. Doll,et al.  A J-WALKING ALGORITHM FOR MICROCANONICAL SIMULATIONS : APPLICATIONS TO LENNARD-JONES CLUSTERS , 1998 .

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

[7]  Wilding.  Critical-point and coexistence-curve properties of the Lennard-Jones fluid: A finite-size scaling study. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[8]  Erik Luijten,et al.  Rejection-free geometric cluster algorithm for complex fluids. , 2004, Physical review letters.

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

[10]  Shang‐keng Ma.  Modern Theory of Critical Phenomena , 1976 .

[11]  B. Alder,et al.  Phase Transition for a Hard Sphere System , 1957 .

[12]  J. D. Doll,et al.  Extending J walking to quantum systems: Applications to atomic clusters , 1992 .

[13]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[14]  Werner Krauth,et al.  Cluster algorithm for hard spheres and related systems , 1995 .

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

[16]  D Chandler,et al.  Van der Waals Picture of Liquids, Solids, and Phase Transformations , 1983, Science.

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

[18]  J. Heringa,et al.  Cluster Monte Carlo: Extending the range , 2002 .

[19]  Lev D. Gelb,et al.  Monte Carlo simulations using sampling from an approximate potential , 2003 .

[20]  C H Mak,et al.  Superfluidity in CH4-doped H2 nanoclusters. , 2005, The Journal of chemical physics.

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

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

[23]  S. Whittington,et al.  Monte carlo study of the interacting self-avoiding walk model in three dimensions , 1996 .

[24]  Ron,et al.  Simulations without critical slowing down. , 1988, Physical review letters.

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

[26]  K. Binder.  Monte Carlo methods in statistical physics , 1979 .