×

分级分数步近似和并行动力学蒙特卡罗算法。 (英语) 兹比尔1259.82110

众所周知,动力学蒙特卡罗(KMC)算法是模拟非平衡、空间分布过程的重要工具。这种模型在材料科学、催化和复杂生物过程中有许多应用。通常,它们用于模拟原子和分子的传输微机械模型、吸附、脱附过程和表面扩散。此外,它们还可以应用于流行病学、生态学甚至交通网络中基于代理的进化博弈问题。
本文提出了一个数学框架,用于构造和分析晶格动力学蒙特卡罗(KMC)模拟的并行算法。所获得的算法能够模拟具有复杂化学和传输微机械的空间分布、非平衡物理化学过程中的广泛时空尺度。可以将这些算法用于特定的分层并行体系结构,例如多核处理器或图形处理单元(GPU)集群。
本文分为八个部分。在初步介绍了一般数学框架之后,我们介绍了分步动力学蒙特卡罗算法。在第3节中,描述了处理器通信调度(PCS)中分步方法的关键特征。在此之后,作者提出了KMC分数阶近似误差分析的形式化论证,这表明了格式的收敛阶。第5节显示了图形处理单元(GPU)的可能应用。第6节介绍了一种概率策略,该策略基于以下大众运输方法的思想,在仿真过程中动态地重新平衡工作负载。接下来,我们将对所提出的方法进行基准测试。

MSC公司:

82C80码 时间相关统计力学的数值方法(MSC2010)
82C20个 含时统计力学中的动态晶格系统(动力学伊辛等)和图上系统
82C70码 含时统计力学中的输运过程
65二氧化碳 蒙特卡罗方法
65日元10 特定类别建筑的数值算法
2007年第47天 马尔可夫半群及其在扩散过程中的应用
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] G.Arampatzis,M.A.Katsoulakis,P.Plechánch,并行动力学蒙特卡罗算法的误差分析:精度和处理器通信,2012年,预印本。;G.Arampatzis,M.A.Katsoulakis,P.Plechánch,并行动力学蒙特卡罗算法的误差分析:精度和处理器通信,2012年,预印本。
[2] 是,S。;Katsoulakis,医学硕士。;Plecháč,P。;Rey-Bellet,L.,《扩展系统粗粒度方案中的多体相互作用》,SIAM J.Sci。计算。,31, 987-1015 (2008) ·Zbl 1191.82007号
[3] Auerbach,S.M.,《纳米孔中跳跃动力学、扩散和相平衡的理论和模拟》,《国际物理评论》。化学。,19 (2000)
[4] S.P.C.Battaile、M.Chandross、L.Holm、A.Thompson、V.Tikare、G.Wagner、E.Webb、X.Zhou、C.G.Cardona、A.Slepoy,《通过平行动力学蒙特卡罗穿越中尺度非人类陆地》,桑迪亚报告,2009年。;S.P.C.Battaile、M.Chandross、L.Holm、A.Thompson、V.Tikare、G.Wagner、E.Webb、X.Zhou、C.G.Cardona、A.Slepoy,《通过平行动力学蒙特卡罗穿越中尺度非人类陆地》,桑迪亚报告,2009年。
[5] Baxter,R.J.,《统计力学中的精确求解模型》(1989),学术出版社·Zbl 0825.58051号
[6] Bortz,A.B。;Kalos,M.H。;Lebowitz,J.L.,伊辛自旋系统蒙特卡罗模拟的新算法,J.Compute。物理。,17, 10-18 (1975)
[7] 查特吉,A。;Vlachos,D.,空间微观和加速动力学蒙特卡罗方法概述,计算机J。助理马特。设计。,14, 253-308 (2007)
[8] Christensen,C.H。;Norskov,J.K.,《多相催化的分子观点》,J.Chem。物理。,128 (2008)
[9] 艾克·S·G。;格林伯格,A.G。;卢巴切夫斯基,B.D。;Weiss,A.,《电路开关网络应用并行模拟的同步松弛》,ACM Trans。模型。计算。模拟。,3, 287-314 (1993) ·Zbl 0839.94015号
[10] Evans,L.C.,偏微分方程和Monge-Kantorovich质量传递,(数学的当前发展,1997(马萨诸塞州剑桥)(1999),国际出版社:马萨诸塞州立波士顿国际出版社),65-126·Zbl 0954.35011号
[11] Gardiner,C.W.,《物理、化学和自然科学随机方法手册》(2004),Springer·Zbl 0862.60050号
[12] Gillespie,D.T.,《数值模拟耦合化学反应随机时间演化的通用方法》,J.Compute。物理。,22, 403-434 (1976)
[13] 海尔,E。;卢比奇,C。;Wanner,G.,几何-数值积分,常微分方程的结构保持算法,计算数学中的Springer级数,第31卷(2006),Springer-Verlag:Springer-Verlag Berlin·兹比尔1094.65125
[14] 海德堡,P。;Nicol,D.M.,使用均匀化的连续时间马尔可夫链的保守并行模拟,IEEE Trans。并行分配系统。,4, 906-921 (1993)
[15] Katsoulakis,M。;Majda,A。;Vlachos,D.,微观晶格系统的粗粒度随机过程,Proc。美国国家科学院。科学,100782-782-782(2003)·Zbl 1063.82033号
[16] Katsoulakis,医学硕士。;Plechánch,P。;Sopasakis,A.,随机晶格动力学粗粒度误差分析,SIAM J.Numer。分析。,44, 2270-2296 (2006) ·Zbl 1140.82035号
[17] Kohn,W.,诺贝尔讲座:物质波函数和密度泛函的电子结构,Rev.Mod。物理。,71, 1253-1266 (1999)
[18] Kornis,G。;诺沃特尼,硕士。;Rikvold,P.A.,《动态蒙特卡罗算法的并行化:部分无拒绝保守方法》,J.Compute。物理。,153, 488-508 (1999) ·Zbl 0979.82043号
[19] Kurtz,T.G.,《随机Trotter产品配方》,Proc。美国数学。《社会学杂志》,35,147-154(1972)·Zbl 0255.47051号
[20] Landau,D.P。;Binder,K.,《统计物理蒙特卡罗模拟指南》(2000),剑桥大学出版社:剑桥大学出版社·Zbl 0998.82504号
[21] Liggett,T.M.,交互粒子系统,Grundlehren der Mathematischen Wissenschaften,第276卷(1985),Springer-Verlag:Springer-Verlag纽约,柏林,海德堡,东京·Zbl 0832.60094号
[22] 刘德杰。;Evans,J.W.,Pd(100)和Rh(100)上CO氧化的原子和多尺度模拟:从纳米尺度波动到中尺度反应前沿,表面科学。,603, 1706-1716 (2009)
[23] Lubachevsky,B.D.,《动态伊辛自旋系统的高效并行模拟》,J.Compute。物理。,75, 103-122 (1988) ·Zbl 0642.65001号
[24] 卢基恩,J.J。;Segers,J.P.L。;Hilbers,P.A.J。;盖尔滕,R.J。;Jansen,A.P.J.,《模拟催化表面反应的高效蒙特卡罗方法》,Phys。修订版E,582598-2610(1998)
[25] 马丁内斯,E。;Marian,J。;Kalos,M.H。;Perlado,J.M.,《连续扩散反应系统的同步平行动力学蒙特卡罗》,J.Compute。物理。,227, 3804-3823 (2008) ·Zbl 1139.65003号
[26] Merrick,M。;Fichthorn,K.A.,薄膜生长并行动力学蒙特卡罗模拟的同步松弛算法,Phys。版本E,75,011606(2007)
[27] Metiu,H.,《专题前言:多相催化的一些新进展》,J.Chem。物理。,128 (2008)
[28] M.长崎。;Kondoh,H。;中井,I。;Ohta,T.,通过动态蒙特卡罗方法研究的Pt(111)上的CO氧化反应,包括吸附质的横向相互作用,J.Chem。物理。,126, 044704-044707 (2007)
[29] Nandipati,G。;垫片Y。;阿马尔,J.G。;卡里姆,A。;卡拉,A。;Rahman,T.S。;Trushin,O.,《使用大型数据库对Ag(111)岛粗化进行并行动力学蒙特卡罗模拟》,J.Phys。康登斯。Matter,21084214(2009)
[30] Onsager,L.,《水晶统计》。I.具有序-序转换的二维模型,Phys。修订版,65117-149(1944)·Zbl 0060.46001号
[31] 佩恩,M.C。;特特,M.P。;阿伦·D.C。;阿里亚斯,T.A。;Joannopoulos,J.D.,《生物总能量计算的迭代最小化技术——分子动力学和共轭梯度》,修订版。物理。,64, 1045-1097 (1992)
[32] 路透社,K。;Frenkel,D。;Schefler,M.,《多相催化的稳态》,第一原理统计力学研究,物理学。修订稿。,93(2004年)
[33] 桑德斯,J。;Kandrot,E.,CUDA举例:通用GPU编程简介(2010),Addison Wesley Professional:Addison Wesley Professional剑桥
[34] 垫片Y。;Amar,J.G.,薄膜生长的并行动力学蒙特卡罗模拟的严格同步松弛算法,Phys。B版,71,115436(2005)
[35] 垫片Y。;Amar,J.G.,薄膜生长并行动力学蒙特卡罗模拟的半严格同步松弛算法,Phys。版本B,71,125432(2005)
[36] Szabo,G。;Fath,G.,《图上的进化游戏》,《物理学》。代表,446,97-216(2007)
[37] Trotter,H.F.,关于算子半群的乘积,Proc。美国数学。《社会学杂志》,10545-551(1959)·Zbl 0099.10401号
[38] Wu,T.T。;McCoy,B.M。;特蕾西,C.A。;Barouch,E.,《二维伊辛模型的自旋相关函数:缩放区域的精确理论》,Phys。B版,13,316-374(1976)
[39] L.Xu、M.Taufer、S.Collins和D.G.Vlachos:GPU上Tau-Leap粗粒度蒙特卡罗模拟的并行化。IEEE/ACM国际并行和分布式处理研讨会(IPDPS)会议记录,2010年4月,美国佐治亚州亚特兰大。;L.Xu、M.Taufer、S.Collins和D.G.Vlachos:GPU上Tau-Leap粗训练蒙特卡罗模拟的并行化。IEEE/ACM国际并行和分布式处理研讨会(IPDPS)会议记录,2010年4月,美国佐治亚州亚特兰大。
[40] Ziff,R.M。;古拉里,E。;Barshad,Y.,《不可逆表面反应模型中的动力学相变》,Phys。修订稿。,56, 2553 (1986)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。