Efficient kinetic Monte Carlo simulation

This paper concerns kinetic Monte Carlo (KMC) algorithms that have a single-event execution time independent of the system size. Two methods are presented-one that combines the use of inverted-list data structures with rejection Monte Carlo and a second that combines inverted lists with the Marsaglia-Norman-Cannon algorithm. The resulting algorithms apply to models with rates that are determined by the local environment but are otherwise arbitrary, time-dependent and spatially heterogeneous. While especially useful for crystal growth simulation, the algorithms are presented from the point of view that KMC is the numerical task of simulating a single realization of a Markov process, allowing application to a broad range of areas where heterogeneous random walks are the dominate simulation cost.

[1]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[2]  Farid F. Abraham,et al.  Computer Simulation of Vapor Deposition on Two-Dimensional Lattices , 1970 .

[3]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[4]  J. L. Blue,et al.  Faster Monte Carlo simulations. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[5]  Muruhan Rathinam,et al.  Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method , 2003 .

[6]  L. Sander,et al.  Multiscale kinetic Monte Carlo algorithm for simulating epitaxial growth , 2005, cond-mat/0504272.

[7]  G. H. Gilmer,et al.  Simulation of Crystal Growth with Surface Diffusion , 1972 .

[8]  P. Maksym.  Fast Monte Carlo simulation of MBE growth , 1988 .

[9]  T. Schulze.  A hybrid scheme for simulating epitaxial growth , 2004 .

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

[11]  J. E. Norman,et al.  A computer program for the generation of random variables from any discrete distribution , 1972 .

[12]  J. Gentle.  Random number generation and Monte Carlo methods , 1998 .

[13]  Miroslav Kotrla,et al.  Numerical simulations in the theory of crystal growth , 1996 .

[14]  T. Schulze.  Kinetic Monte Carlo simulations with minimal searching. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  George Marsaglia.  Generating discrete random variables in a computer , 1963, CACM.

[16]  Eric Vanden-Eijnden,et al.  Nested stochastic simulation algorithm for chemical kinetic systems with disparate rates. , 2005, The Journal of chemical physics.