×

为类朗之万方程的显式积分器提供无条件稳定性的补丁。 (英语) Zbl 1430.65003号

小结:本文扩展了作者在《普通应用数学》第63卷第5期第655-696页(2010;Zbl 1214.60031号)]到分子动力学中出现的随机微分方程(SDE)。它实现了对显式积分器的补丁,该补丁由Metropolis-Hastings步骤组成。“修补积分器”保持了SDE的平衡分布,并且在有限的时间间隔内是准确的。作为推论,本文证明了积分器在估计无限长解的有限时间动力学方面的准确性——这在分子动力学中尚属首次。本文还涵盖了多个时间步长、完整约束和可扩展性。最后,本文提供了支持该理论的数值试验。

MSC公司:

65立方米 随机微分和积分方程的数值解
82立方31 随机方法(福克-普朗克、朗之万等)应用于含时统计力学问题
60华氏35 随机方程的计算方法(随机分析方面)

软件:

GSHMC公司
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿赫马茨卡亚,E。;北卡罗来纳州Bou-Rabee。;Reich,S.,《带动量翻转和不带动量翻转的广义混合蒙特卡罗方法的比较》,J.Compute。物理。,228, 2256-2265 (2009) ·Zbl 1161.65302号
[2] Akhmatskaya,E。;Reich,S.,GSHMC:一种有效的分子模拟方法,J.Compute。物理。,227, 4937-4954 (2008) ·Zbl 1148.82316号
[3] 艾伦,M.P。;Tildesley,D.J.,《液体的计算机模拟》(1987),克拉伦登出版社·Zbl 0703.68099号
[4] Andersen,H.C.,恒定压力和/或温度下的分子动力学模拟,J.Chem。物理。,72, 2384 (1980)
[5] N.Bou-Rabee和M.Hairer,MALA算法的非症状混合,<arXiv:1008.3514v1[math.PR]>;N.Bou-Rabee和M.Hairer,MALA算法的非症状混合,<arXiv:1008.3514v1[math.PR]>
[6] 北卡罗来纳州Bou-Rabee。;Owhadi,H.,随机变分积分器,IMA J.Numer。分析。,29, 421-443 (2009) ·Zbl 1171.37027号
[7] 北卡罗来纳州Bou-Rabee。;Owhadi,H.,随机背景下变分积分器的长运行精度,SIAM J.Numer。分析。,48, 278-297 (2010) ·Zbl 1215.65012号
[8] 北卡罗来纳州Bou-Rabee。;Vanden-Eijnden,E.,SDE的都市积分器的路径准确性和遍历性,CPAM,63,655-696(2010)·Zbl 1214.60031号
[9] Brünger,A。;布鲁克斯,C.L。;Karplus,M.,ST2水分子动力学模拟的随机边界条件,化学。物理学。莱特。,105, 495-500 (1984)
[10] Bussi,G。;多纳迪奥,D。;Parrinello,M.,《通过速度缩放的典型采样》,J.Chem。物理。,126, 014101 (2007)
[11] 坎塞斯,E。;勒戈尔,F。;Stoltz,G.,分子动力学一些采样方法的理论和数值比较,数学。模型。数字。分析。,41, 351-389 (2007) ·Zbl 1138.82341号
[12] Ciccotti,G。;Lelievre,T。;Vanden-Eijnden,E.,《子流形上扩散的投影:应用于平均力计算》,CPAM,61,0001-0039(2008)
[13] 杜安,S。;A.D.肯尼迪。;彭德尔顿,B.J。;Roweth,D.,混合蒙特卡罗,Phys。莱特。B、 195216-222(1987)
[14] E.渭南。;Li,D.,分子动力学中的安徒生恒温器,CPAM,61,96-136(2008)·Zbl 1148.82019年
[15] Feng,K。;Wang,D.-L.,《动力系统与算法的几何构造》,《现代数学丛书》,第163卷(1994年),AMS,第1-32页·Zbl 0805.70005号
[16] W.Fong,《粒子模拟的时空多尺度方法》,斯坦福大学博士论文,2009年。;W.Fong,《粒子模拟的时空多尺度方法》,斯坦福大学博士论文,2009年。
[17] 方,W。;Darve,E。;Lew,A.,异步变分积分器的稳定性,J.Compute。物理。,227, 8367-8394 (2008) ·Zbl 1153.65120号
[18] Frenkel,D。;Smit,B.,《理解分子模拟:从算法到应用》(2002),学术出版社
[19] Giles,M.B.,多层蒙特卡罗路径模拟,Oper。研究,56,607-617(2008)·Zbl 1167.65316号
[20] 海尔,E。;卢比奇,C。;Wanner,G.,《几何-数值积分》,《计算数学中的Springer系列》,第31卷(2006年),Springer
[21] Hartmann,C.,《约束哈密顿系统的遍历采样方案及其在分子动力学中的应用》,J.Stat.Phys。,130, 687-712 (2008) ·Zbl 1214.82063号
[22] Hastings,W.K.,使用马尔可夫链的蒙特卡罗方法及其应用,生物医学,57,97-109(1970)·兹比尔0219.65008
[23] 海姆·D·J。;毛,X。;Stuart,A.M.,非线性随机微分方程欧拉型方法的强收敛性,SIAM J.Numer。分析。,40, 1041-1063 (2002) ·Zbl 1026.65003号
[24] Horowitz,A.M.,广义制导蒙特卡罗算法,Phys。莱特。B、 268247-252(1991)
[25] 伊扎吉雷,J.A。;Catarello博士。;沃兹尼亚克,J.M。;Skeel,R.D.,分子动力学的Langevin稳定性,化学杂志。物理。,114, 2090 (2001)
[26] A.D.肯尼迪。;Pendleton,B.,自由场理论广义混合蒙特卡罗算法的成本,Nucl。物理学。B、 607456-510(2001)·兹伯利0969.81639
[27] 莱姆库勒,B。;Noorizadeh,E。;Theil,F.,《分子动力学的温和随机恒温器》,J.Stat.Phys。,135, 261-277 (2009) ·Zbl 1179.82065号
[28] 莱姆库勒,B。;Reich,S.,《模拟哈密顿动力学》,剑桥应用和计算数学专著(2004),剑桥大学出版社
[29] 莱姆库勒,B。;Skeel,R.,约束哈密顿系统中的辛数值积分,JCP,112117-125(1994)·Zbl 0817.65057号
[30] T.Lelièvre,M.Rousset,G.Stoltz,带约束的Langevin动力学和自由能差计算,<arXiv:1006.4914v1>;T.Lelièvre,M.Rousset,G.Stoltz,带约束的Langevin动力学和自由能差计算,<arXiv:1006.4914v1>
[31] 卢,A。;马斯登,J.E。;奥尔蒂斯,M。;West,M.,异步变分积分器,Arch。定额。机械。分析。,167, 85-145 (2003) ·Zbl 1055.74041号
[32] Li,D.,《分子动力学中安徒生恒温器收敛到平衡的速度》,J.Stat.Phys。,129, 265-287 (2007) ·2014年8月11日
[33] 马奇。;伊扎吉雷,J.A。;Skeel,R.D.,VERLET-I/R-RESPA/IMPULSE受到非线性不稳定性的限制,SIAM J.Sci。计算。,24, 1951-1973 (2003) ·Zbl 1163.65317号
[34] J.C.马丁利。;Stuart,A.M。;Higham,D.J.,SDE和近似的遍历性:局部Lipschitz向量场和退化噪声,Stoch。程序。申请。,101, 2, 185-232 (2002) ·Zbl 1075.60072号
[35] Mengersen,K.L。;Tweedie,R.L.,《黑斯廷斯和大都会算法的收敛速度》,《Ann.Stat.》,24,101-121(1996)·Zbl 0854.60065号
[36] 北卡罗来纳州大都会。;罗森布鲁思,A.W。;出纳员,A.H。;Teller,E.,《快速计算机器的状态方程计算》,J.Chem。物理。,21, 1087-1092 (1953) ·Zbl 1431.65006号
[37] 米尔斯坦,G.N。;Tretyakov,M.V.,非全局Lipschitz系数随机微分方程的数值积分,SIAM J.Numer。分析。,43, 1139-1154 (2005) ·Zbl 1102.60059号
[38] Roberts,G.O。;Tweedie,R.L.,多维Hastings和Metropolis算法的几何收敛和中心极限定理,Biometrika,195-110(1996)·兹标0888.60064
[39] Rossky,P.J。;多尔,J.D。;Friedman,H.L.,《作为智能蒙特卡罗模拟的布朗动力学》,化学杂志。物理。,69, 4628 (1978)
[40] Ryckaert,J。;Ciccotti,G。;Berendsen,H.,带约束系统笛卡尔运动方程的数值积分:正构烷烃的分子动力学,JCP,23,327-341(1977)
[41] 萨莫列托夫,A.A。;牧师,文学硕士。;Dettmann,C.P.,“慢速”配置模式的恒温器,J.Stat.Phys。,128, 1321-1336 (2008) ·Zbl 1128.82006年
[42] Scemama,A。;Lelièvre,T。;斯托尔茨,G。;坎塞斯,E。;Caffarel,M.,《变分蒙特卡罗的有效采样算法》,J.Chem。物理。,125, 114105 (2006)
[43] 施耐德,T。;Stoll,E.,扭曲相变三维单组分模型的分子动力学研究,Phys。B版,第17页,第1302-1322页(1978年)
[44] C.Schütte,《构象动力学:建模、理论、算法和生物分子应用》,《习惯化》,柏林自由大学,1999年。;C.Schütte,《构象动力学:建模、理论、算法和生物分子应用》,《习惯化》,柏林自由大学,1999年。
[45] G.Stoltz,《分子和多尺度模拟的一些数学方法》,博士论文,国立桥与肖塞大学,2007年。;G.Stoltz,《分子和多尺度模拟的一些数学方法》,博士论文,国立桥与乔斯学院,2007年。
[46] Talay,D.,《随机微分系统的模拟和数值分析:综述》,(Krèe,P.;Wedig,W.,《应用物理学中的概率方法》,第451卷(1995),Springer-Verlag:Springer-Verlag Berlin),54-96·Zbl 0837.65150号
[47] Talay,D.,随机哈密顿系统:指数收敛到不变测度,并通过隐式Euler方案离散化,Markov过程。相关领域,8,1-36(2002)·Zbl 1011.60039号
[48] 塔莱,D。;Tubaro,L.,求解随机微分方程数值格式的全局误差展开,Stoch。分析。申请。,8, 94-120 (1990)
[49] 塔克曼,M.E。;Berne,B.J.,具有多个时间尺度和记忆摩擦的系统中的随机分子动力学,J.Chem。物理。,95, 4389-4396 (1991)
[50] Vanden-Eijnden,E。;Ciccotti,G.,具有完整约束的Langevin方程的二阶积分器,化学。物理学。莱特。,429, 310-316 (2006)
[51] Vanden-Eijnden,E。;Ventroli,M.,《通过轨迹平行化和扭曲进行精确速率计算》,J.Chem。物理。,131, 044120 (2009)
[52] Vanden-Eijnden,E。;Ventroli,M.,《用Voronoi镶嵌法进行马尔科夫磨石》,J.Chem。物理。,130194101(2009年)
[53] Verlet,L.,经典流体的计算机“实验”。I.Lennard-Jones分子的热力学性质,Phys。修订版,159,98-103(1967)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。