Coarse-Grained Langevin Approximations and Spatiotemporal Acceleration for Kinetic Monte Carlo Simulations of Diffusion of Interacting Particles

被引:1
作者
Are, Sasanka [1 ]
Katsoulakis, Markos A. [1 ,2 ,3 ]
Szepessy, Anders [4 ]
机构
[1] Univ Massachusetts, Dept Math & Stat, Amherst, MA 01003 USA
[2] Univ Crete, Dept Appl Math, Iraklion 71405, Greece
[3] Fdn Res & Technol Hellas, Iraklion 71405, Greece
[4] Kungliga Tekniska Hogskolan Royal Inst Technol, Inst Matemat, SE-10044 Stockholm, Sweden
基金
瑞典研究理事会; 美国国家科学基金会;
关键词
Kinetic Monte Carlo methods; Diffusion; Fluctuations; STOCHASTIC-PROCESSES; EXCLUSION PROCESS; DYNAMICS; SCHEMES; LIMITS;
D O I
10.1007/s11401-009-0219-x
中图分类号
O1 [数学];
学科分类号
0701 ; 070101 ;
摘要
Kinetic Monte Carlo methods provide a powerful computational tool for the simulation of microscopic processes such as the diffusion of interacting paxticles oil a surface, at a detailed atomistic level. However such algorithms are typically computationally expensive and are restricted to fairly small spatiotemporal scales. One approach towards overcoming this problem was the development of coarse-gained Monte Carlo algorithms. In recent literature, these methods were shown to be capable of efficiently describing much larger length scales while still incorporating information on microscopic interactions and fluctuations. In this paper, a coarse-grained Langevin system of stochastic differential equations as approximations of diffusion of interacting particles is derived, based on these earlier coarse-grained models. The authors demonstrate the asymptotic equivalence of transient and long time behavior of the Langevin approximation and the underlying microscopic process, using asymptotics methods such as large deviations for interacting particles systems, and furthermore, present corresponding numerical simulations, comparing statistical quantities like mean paths, auto correlations and power spectra of the microscopic and the approximating Langevin processes. Finally, it is shown that the Langevin approximations presented here are much more computationally efficient than conventional Kinetic Monte Carlo methods, since in addition to the reduction in the number of spatial degrees of freedom in coarse-grained Monte Carlo methods, the Langevin system of stochastic differential equations allows for multiple particle moves in a single timestep.
引用
收藏
页码:653 / 682
页数:30
相关论文
共 32 条
[1]   The effect of bond length on the structure of dense bead-spring polymer melts [J].
Abrams, CF ;
Kremer, K .
JOURNAL OF CHEMICAL PHYSICS, 2001, 115 (06) :2776-2785
[2]  
Allen M. P., 1987, COMPUTER SIMULATION
[3]   MULTIBODY INTERACTIONS IN COARSE-GRAINING SCHEMES FOR EXTENDED SYSTEMS [J].
Are, Sasanka ;
Katsoulakis, Markos A. ;
Plechac, Petr ;
Rey-Bellet, Luc .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2008, 31 (02) :987-1015
[4]   Metastability for the exclusion process with mean-field interaction [J].
Asselah, A ;
Giacomin, G .
JOURNAL OF STATISTICAL PHYSICS, 1998, 93 (5-6) :1051-1110
[5]   Langevin equations for fluctuating surfaces [J].
Chua, ALS ;
Haselwandter, CA ;
Baggio, C ;
Vvedensky, DD .
PHYSICAL REVIEW E, 2005, 72 (05)
[6]  
COMETS F, 1987, ANN I H POINCARE-PR, V23, P135
[7]  
Dupuis P., 1997, WEAK CONVERGENCE APP
[8]  
Friedlin M.I., 1998, RANDOM PERTURBATIONS
[9]  
Gardiner C., 2009, Handbook of stochastic methods for physics, chemistry and the natural sciences, V4th edn
[10]   Phase segregation dynamics in particle systems with long range interactions .1. Macroscopic limits [J].
Giacomin, G ;
Lebowitz, JL .
JOURNAL OF STATISTICAL PHYSICS, 1997, 87 (1-2) :37-61