×

基于自适应权值学习的无限开关极限模拟回火方法。 (英语) Zbl 07382763号

总结:我们研究了模拟回火(ST)方法的理论基础,并利用我们的研究结果设计了一种高效的加速采样算法。首次用于复制交换分子动力学的大偏差论点[N.普拉特纳等,“稀土采样问题的无限交换方法”,J.Chem。物理学。135,文章ID 134111,24 p.(2011;数字对象标识代码:10.1063/1.3643325)],我们证明了模拟回火最有效的方法是无限快速地改变温度。在这个极限条件下,我们可以用后者的平均方程来代替温度和物理变量的运动方程,根据根据温度权重定义的位置相关函数重新缩放力。平均方程与Gao的积分超温方法中使用的方程类似,只是我们表明使用连续温度比使用离散温度更好。我们给出了选择温度权重作为倒数配分函数的理论依据,从而将模拟回火与Wang-Landau取样联系起来。最后,我们描述了一种在仿真过程中同时对正则系综进行采样和学习权重的自持算法。这种无限开关模拟回火(ISST)算法在三个日益复杂的示例上进行了测试:谐波振荡器系统;居里-维斯模型的连续变型,其中ISST的性能优于标准ST,并能准确捕获该模型中观察到的二阶相变;和丙氨酸-12,ISST在计算300-500 K温度范围内分子的均方根偏差(RMSD)和回转半径的自由能分布方面也优于标准ST。

MSC公司:

82M99型 统计力学的基本方法
60层10 大偏差
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Frenkel博士。;Smit,B.,《理解分子模拟:从算法到应用》,第1卷(2001),阿姆斯特丹:爱思唯尔出版社
[2] 麦凯,D.J C.,《信息理论、推理和学习算法》,IEEE Trans。《信息论》,502544-2545(2003)·Zbl 1055.94001号
[3] 理查德,J-F;张伟,高效高维重要性抽样,《计量经济学》,1411385-1411(2007)·Zbl 1420.65005号 ·doi:10.1016/j.jeconom.2007.02.007
[4] Landau,D.P。;Binder,K.,《统计物理蒙特卡罗模拟指南》(2014),剑桥:剑桥大学出版社,剑桥
[5] Wang,F。;Landau,D.P.,计算状态密度的高效多范围随机行走算法,Phys。修订稿。,86, 2050-2053 (2001) ·doi:10.1103/PhysRevLett.86.2050
[6] Junghans,C。;佩雷斯,D。;Vogel,T.,《多正则系综中的分子动力学:wanglandau采样的等效性、统计温度分子动力学和代谢动力学》,J.Chem。理论计算。,10, 1843-1847 (2014) ·doi:10.1021/ct500077d
[7] 莱奥,A。;Parrinello,M.,《逃离自由能量极小值》,Proc。美国国家科学院。科学。,99, 12562-12566 (2002) ·doi:10.1073/pnas.202427399
[8] Barducci,A。;Bussi,G。;Parrinello,M.,《温控元动力学:一种平滑收敛和可调的自由能方法》,Phys。修订稿。,100 (2008) ·doi:10.1103/PhysRevLett.100.020603
[9] Maragliano,L。;Vanden-Eijnden,E.,《稀有事件模拟中自由能取样和确定反应路径的温度加速方法》,《化学》。物理学。莱特。,426, 168-175 (2006) ·doi:10.1016/j.cplett.2006.05.062
[10] 艾布拉姆斯,C.F。;Vanden-Eijnden,E.,使用温度加速分子动力学对蛋白质进行大规模构象取样,Proc。美国国家科学院。科学。美国,1074961-4966(2010)·doi:10.1073/pnas.0914540107
[11] Darve,E。;罗德里格斯-戈麦斯,D。;Pohorille,A.,标量和矢量自由能计算的自适应偏压法,J.Chem。物理。,128 (2008) ·doi:10.1063/1.2829861
[12] Aluffi-Pentini,F。;Parisi,V.公司。;Zirilli,F.,《全局优化与随机微分方程》,J.Optim。理论应用。,47, 1-16 (1985) ·Zbl 0549.65038号 ·doi:10.1007/BF00941312
[13] Swendsen,R.H。;Wang,J.S.,《旋转玻璃的复制蒙特卡罗模拟》,Phys。修订稿。,57, 2607-2609 (1986) ·doi:10.1103/PhysRevLett.57.2607
[14] Geyer,C.J.,马尔可夫链蒙特卡罗最大似然,156-163(1991)
[15] Neal,R.M.,使用回火转变从多峰分布中取样,统计计算。,6, 353-366 (1996) ·doi:10.1007/BF00143556
[16] Hansmann,U.H E.,生物分子构象研究的并行回火算法,化学。物理学。莱特。,281, 140-150 (1997) ·doi:10.1016/S0009-2614(97)01198-6
[17] Y.Sugita。;Okamoto,Y.,蛋白质折叠的副本交换分子动力学方法,化学。物理学。莱特。,314, 141-151 (1999) ·doi:10.1016/S0009-2614(99)01123-9
[18] E.马里纳里。;Parisi,G.,《模拟回火:一种新的蒙特卡罗方案》,Europhys。莱特。,19, 451-458 (1992) ·doi:10.1209/0295-5075/19/6/002
[19] Hansmann,U.H E。;Okamoto,Y.,《蛋白质折叠问题中最近提出的三种算法的数值比较》,J.Compute。化学。,18, 920-933 (1997) ·doi:10.1002/(SICI)1096-987X(199705)18:7<920::AID-JCC5>3.0.CO;2-T型
[20] 班巴,S。;郭,Z。;Brooks,C.L III,《使用λ动力学方法计算自由能时配体取向和构象的有效采样》,J.Phys。化学。B、 104、6903-6910(2000)·doi:10.1021/jp001177i
[21] 伍兹,C.J。;埃塞克斯,J.W。;King,M.A.,《结合自由能计算中的增强构型取样》,社会,13711-13718(2003)·doi:10.1021/jp036162+
[22] 比特蒂·普策(Bitetti-Putzer,R.)。;杨伟(Yang,W.)。;Karplus,M.,广义系综用于提高自由能模拟的收敛性,化学。物理学。莱特。,377, 633-641 (2003) ·doi:10.1016/S0009-2614(03)01057-1
[23] Okamoto,Y.,《广义系综算法:蒙特卡洛和分子动力学模拟的增强采样技术》,J.Mol.Graph。建模,22,425-439(2004)·doi:10.1016/j.jmgm.2003.12.009
[24] Faraldo-Gómez,J.D。;Roux,B.,《通过哈密顿量和温度复制交换模拟表征构象平衡:评估熵和环境效应》,J.Compute。化学。,28, 1634-1647 (2007) ·doi:10.1002/jcc.20652
[25] 赫里茨,J。;Oostenbrink,C.,使用软核相互作用的哈密顿复制交换分子动力学,J.Chem。物理。,128 (2008) ·数字标识代码:10.1063/1.2888998
[26] 冈崎,K-I;北卡罗来纳州科加。;高田,S。;Onuchic,J.N。;Wolynes,P.G.,《蛋白质大幅度构象运动的多池能量景观:基于结构的分子动力学模拟》,Proc。美国国家科学院。科学。美国,10311844-11849(2006)·doi:10.1073/pnas.0604375103
[27] 苗,Y。;辛科,W。;皮尔斯,L。;Bucher,D。;R.C.沃克。;Mccammon,J.A.,《自由能计算加速分子动力学模拟的改进重加权》,J.Chem。理论计算。,10, 2677-2689 (2014) ·doi:10.1021/ct500090q
[28] 刘,P。;Kim,B。;弗里斯纳,R.A。;Berne,B.J.,《溶质回火的复制品交换:显式水中生物系统取样方法》,Proc。美国国家科学院。科学。美国,10213749-13754(2005)·doi:10.1073/pnas.0506346102
[29] 普拉特纳,N。;多尔,J.D。;Dupuis,P。;Wang,H。;刘,Y。;Gubernatis,J.E.,《稀土采样问题的无限交换方法》,J.Chem。物理。,135 (2011) ·doi:10.1063/1.3643325
[30] Dupuis,P。;刘,Y。;普拉特纳,N。;Doll,J.D.,关于平行回火的无限交换极限,多尺度模型。模拟。,10, 986-1022 (2012) ·Zbl 1261.65009号 ·doi:10.1137/10853145
[31] Gao,Y.Q.,《用于增强采样的集成过温方法》,J.Chem。物理。,128 (2008) ·doi:10.1063/1.2825614
[32] 戈博,G。;Leimkuhler,B.J.,《连续回火的扩展哈密顿方法》,物理。E版,91(2015)·doi:10.1103/PhysRevE.91.061301
[33] 帕克,S。;潘德,V.S.,《为模拟回火选择砝码》,物理。E版,76306(2007年)·doi:10.1103/PhysRevE.76.016703
[34] Roberts,G.O。;Rosenthal,J.S.,通过扩散极限最小化MCMC方差,并应用于模拟回火,Ann.Appl。概率。,24, 131-149 (2014) ·Zbl 1298.60078号 ·doi:10.1214/12-AAP918
[35] Tan,Z.,最优调整混合抽样和局部加权直方图分析,J.Compute。图表。Stat.,26,54-65(2017)·doi:10.1080/10618600.2015.1113975
[36] Yang,Y.I。;牛,H。;Parrinello,M.,《综合力学和综合回火取样》,J.Phys。化学。莱特。,9, 6426-6430 (2018) ·doi:10.1021/acs.jpclett.8b03005
[37] 卡尔森,D.E。;Stinson,P。;Pakman,A。;Paninski,L.,Rao-Blackwellized回火取样的配分函数,第48卷(2016)
[38] 马德拉斯,N。;Piccioni,M.,《分布族重要性抽样》,Ann.Appl。概率。,9, 1202-1225 (1999) ·Zbl 0966.60061号 ·doi:10.1214/aoap/1029962870
[39] 医学博士唐斯克。;Varadhan,S.R.S.,大时间特定马尔可夫过程期望的渐近评估,I,Commun。纯应用程序。数学。,28, 1-47 (1975) ·Zbl 0323.60069号 ·doi:10.1002/cpa3160280102
[40] 德乌舍尔,J.D。;斯特罗克,D.W.,《大偏差》(1989),纽约:学术出版社,纽约·Zbl 0705.60029号
[41] Wu,L.,随机阻尼哈密顿系统的大、中偏差和指数收敛,Stoch。过程。申请。,91, 205-238 (2001) ·Zbl 1047.60059号 ·doi:10.1016/S0304-4149(00)00061-2
[42] 博迪诺,T。;Lefevere,R.,耦合到随机恒温器的晶格哈密顿动力学的大偏差,J.Stat.Phys。,133, 1-27 (2008) ·Zbl 1157.82027号 ·doi:10.1007/s10955-008-9601-4
[43] 莱姆库勒,B。;Matthews,C.,《分子动力学与确定性和随机数值方法》(2015),柏林:施普林格出版社,柏林·Zbl 1351.82001号
[44] 莱姆库勒,B。;Matthews,C.,《分子取样随机数值方法的理性构造》,应用。数学。Res.Express,2013,34-56(2013)·Zbl 1264.82102号 ·doi:10.1093/amrx/abs010
[45] Nguyen,P.H。;冈本,Y。;Derreumaux,P.,《通信:模拟回火与快速飞行重量测定》,J.Chem。物理。,138 (2013) ·doi:10.1063/1.4792046
[46] 白求恩,I。;巴尼什,R。;Breitmoser,E。;科利斯,A.B K。;Gibb,G。;戈博,G。;马修斯,C。;Ackland,G.J。;Leimkuhler,B.J.,MIST:用于积分器开发的简单有效的分子动力学抽象库,Comp。物理学。通信。(2018) ·doi:10.1016/j.cpc.2018.10.006
[47] Barth,E。;库泽拉,K。;莱姆库勒,B。;Skeel,R.,《约束分子动力学算法》,J.Comp。化学。,16, 1192-1209 (1995) ·doi:10.1002/jcc.540161003
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。