×
本尼迪克特·莱姆库勒；尚、小程

成对自适应恒温器可提高耗散粒子动力学的准确性和稳定性。 （英语） Zbl 1360.82089号

J.计算。物理学。 324, 174-193 (2016).
小结：我们研究了耗散粒子动力学（DPD）和动量守恒分子动力学的公式和数值处理。我们表明，通过使用两两自适应Langevin恒温器可以提高DPD的准确性和稳定性，该恒温器精确匹配DPD模拟的动力学特性（例如，自相关函数），同时使用负反馈回路自动校正热力学平均值。在低摩擦状态下，可以使用基于仅对相互作用进行恒温的Nosé-Hoover-Langevin方法的更简单的动量守恒变体来替代DPD；我们证明了这种方法对于一类重要的可观察性（超收敛结果）具有额外的精度，同时也允许比替代方法更大的时间步长。本文中提到的所有方法都很容易实现。在平衡和非平衡条件下进行了数值实验；使用Lees-Edwards边界条件诱导剪切流。

引用于5文件

MSC公司：

82C80码 时间相关统计力学的数值方法（MSC2010）
82C22型 含时统计力学中的相互作用粒子系统

关键词：

耗散粒子动力学；成对编号-胡佛-朗格文恒温器；成对自适应Langevin恒温器；收敛阶；构型温度；动量守恒；随机微分方程；分子动力学

软件：

CP2K系统
PDF格式BibTeX公司 XML格式引用
\textit{B.Leimkuhler}和\textit{X.Shang}，J.Compute。物理学。324174-193（2016年；Zbl 1360.82089）
全文： 内政部 arXiv公司

参考文献：

[1] 艾伦，M.P。；Tildesley，D.J.，《液体的计算机模拟》（1989），牛津大学出版社·Zbl 0703.68099号
[2] Frenkel，D。；Smit，B.，《理解分子模拟：从算法到应用》（2001），学术出版社
[3] Hoogerbrugge，P。；Koelman，J.，用耗散粒子动力学模拟微观流体动力学现象，Europhys。莱特。，19, 3, 155 (1992)
[4] Dünweg，B.，《分子动力学算法和流体动力学筛选》，J.Chem。物理。，99, 9, 6977-6982 (1993)
[5] 弗里什，美国。；Hasslacher，B。；Pomeau，Y.，Navier-Stokes方程的格-加斯自动机，Phys。修订稿。，56, 14, 1505-1508 (1986)
[6] Evans，D.J。；Morris，G.P.，《简单流体中的剪切增稠和湍流》，Phys。修订稿。，56, 20, 2172 (1986)
[7] Evans，医学博士。；Morris，G.，《非平衡液体的统计力学》（2008），剑桥大学出版社·Zbl 1189.82064号
[8] 索德曼，T。；Dünweg，B。；Kremer，K.，《耗散粒子动力学：平衡和非平衡分子动力学模拟的有用恒温器》，Phys。E版，68，4，第046702条，pp.（2003）
[9] Español，P.，耗散粒子动力学的流体动力学，物理学。修订版E，521734-1742（1995）
[10] Symeonidis，V。；卡尼亚达基斯，G。；Caswell，B.，《多尺度复杂流体模拟的无缝方法》，计算。科学。工程，7，3，39-46（2005）
[11] 潘，W。；卡斯维尔，B。；Karniadakis，G.E.，《布朗胶体悬浮液的流变学、微观结构和迁移》，Langmuir，26，1，133-142（2010）
[12] 潘·D。；Phan Thien，N。；Khoo，B.C.，低毛细数剪切流中液滴悬浮的耗散粒子动力学模拟，J.Non-Newton。流体力学。，212, 63-72 (2014)
[13] 李，X。；弗拉霍夫斯卡，P.M。；Karniadakis，G.E.，《健康和疾病中红细胞形状和动力学的连续和基于粒子的建模》，《软物质》，9，1，28-37（2013）
[14] Pastorino，C。；Krer，T。；米勒，M。；Binder，K.，《耗散粒子动力学和朗之万恒温器在聚合物系统非平衡模拟中的比较》，Phys。修订版E，76，2，第026706页（2007年）
[15] 贝索尔德，G。；瓦图莱宁，I。；卡图恩，M。；Polson，J.M.，《耗散粒子动力学模拟的更好积分器》，物理。E版，62、6，第R7611条，pp.（2000）
[16] 瓦图莱宁，I。；卡图恩，M。；贝索尔德，G。；Polson，J.M.，《耗散粒子动力学模拟的集成方案：从软交互系统到混合模型》，J.Chem。物理。，116, 3967 (2002)
[17] 尼库宁，P。；卡图恩，M。；Vattulainen，I.，你将如何在耗散粒子动力学模拟中整合运动方程？，计算。物理学。社区。，153, 3 (2003)
[18] Chaudhri，A。；Lukes，J.R.，等温耗散粒子动力学中的速度和应力自相关衰减，物理学。版本E，81，2，第026707条pp.（2010）
[19] Lowe，C.，耗散粒子动力学的另一种方法，欧洲物理学。莱特。，47, 2, 145 (1999)
[20] De Fabritis，G。；塞拉诺，M。；埃斯帕诺，P。；Coveney，P.，《随机模型的有效数值积分器》，Phys。A、 统计机械。申请。，361, 2, 429-440 (2006) ·Zbl 1102.82022号
[21] Shardlow，T.，耗散粒子动力学的分裂，SIAM J.Sci。计算。，24, 4, 1267-1282 (2003) ·兹比尔1043.60048
[22] Peters，E.A.J.F.，消除DPD中的时间步长效应，Europhys。莱特。，66, 3, 311 (2004)
[23] 斯托亚诺夫，S.D。；Groot，R.D.，《从分子动力学到流体动力学：新型伽利略恒温器》，J.Chem。物理。，122，第114112条pp.（2005）
[24] 艾伦，M.P。；Schmid，F.，复杂流体分子动力学恒温器，分子模拟。，33, 1-2, 21-26 (2007) ·Zbl 1118.81533号
[25] 莱姆库勒，B。；Shang，X.，《耗散粒子动力学及相关系统的数值处理》，J.Compute。物理。，280, 72-95 (2015) ·Zbl 1349.76705号
[26] Jones，A。；Leimkuhler，B.，采样驱动分子系统的自适应随机方法，J.Chem。物理。，135，8，第084125条pp.（2011）
[27] 莱姆库勒，B。；Shang，X.，《噪声梯度系统的自适应恒温器》，SIAM J.Sci。计算。，38、2、A712-A736（2016）·Zbl 1382.65019号
[28] X.Shang。；朱，Z。；莱姆库勒，B。；Storkey，A.J.，用于大规模贝叶斯采样的协方差控制自适应Langevin恒温器，（神经信息处理系统进展，第28卷（2015）），37-45
[29] 埃斯帕诺，P。；Warren，P.，《耗散粒子动力学的统计力学》，欧罗皮斯出版社。莱特。，1995年4月30日
[30] Groot，R.D。；Warren，P.B.，《耗散粒子动力学：弥合原子模拟和介观模拟之间的差距》，J.Chem。物理。，107, 4423 (1997)
[31] Lyubartsev，A。；Laaksonen，A.，从径向分布函数计算有效相互作用势：反向蒙特卡罗方法，物理学。E版，52、4、3730-3737（1995）
[32] Lei，H。；卡斯维尔，B。；Karniadakis，G.E.，从微观模拟直接构建介观模型，物理学。版本E，81，2，第026704条pp.（2010）
[33] 李，Z。；卞，X。；卡斯维尔，B。；Karniadakis，G.E.，通过Mori-Zwanzig公式构建复杂流体的耗散粒子动力学模型，软物质，10，43，8659-8672（2014）
[34] Suter，J.L。；格伦·D·。；Coveney，P.V.，《粘土-聚合物纳米复合材料的化学特定多尺度建模揭示了插层动力学、塔克体自组装和涌现材料特性》，Adv.Mater。，27, 6, 966-984 (2014)
[35] Moeendarbary，E。；Ng，T。；Zangeneh，M.，《耗散粒子动力学：导论、方法和复杂流体应用——综述》，《国际应用杂志》。机械。，1, 04, 737-763 (2009)
[36] Shardlow，T。；Yan，Y.，耗散粒子动力学的几何遍历性，Stoch。动态。，6, 1, 123-154 (2006) ·兹比尔1127.60066
[37] 莱姆库勒，B。；Noorizadeh，E。；Theil，F.，《分子动力学的温和随机恒温器》，J.Stat.Phys。，135, 2 (2009) ·Zbl 1179.82065号
[38] 莱姆库勒，B。；马修斯，C。；Stoltz，G.，平衡和非平衡Langevin分子动力学平均值的计算，IMA J.Numer。分析。，36，1，13-79（2016）·Zbl 1347.65014号
[39] Sivak，D.A。；乔德拉，J.D。；Crooks，G.E.，《利用非平衡涨落定理理解和纠正离散Langevin动力学平衡和非平衡模拟中的错误》，Phys。第X版，第3、1条，第011007页（2013年）
[40] 李，Z。；科尔莫德，J.R。；De Vita，A.，《分子动力学与量子力学的在线机器学习》，Phys。修订稿。，114，9，第096405条pp.（2015）
[41] 莫内斯，L。；Jones，A。；Götz，A.W。；莱诺，T。；R.C.沃克。；莱姆库勒，B。；Cányi，G。；Bernstein，N.，CP2K和AMBER软件包中的自适应缓冲力QM/MM方法，J.Compute。化学。，36, 9, 633-648 (2015)
[42] Backer，J.A。；Lowe，C.P。；霍夫斯洛特，H.C.J。；Iedema，P.D.，《测量颗粒模型流体粘度的泊松流》，J.Chem。物理。，122，第154503条pp.（2005）
[43] Nosé，S.，《恒温分子动力学方法的统一公式》，J.Chem。物理。，81, 1, 511 (1984)
[44] 胡佛，W.G.，计算统计力学，现代热力学研究（1991），爱思唯尔科学
[45] 莱姆库勒，B。；Matthews，C.，《分子取样随机数值方法的理性构造》，应用。数学。Res.Express，2013，1，34-56（2013）·兹比尔1264.82102
[46] 莱姆库勒，B。；Matthews，C.，《通过朗之万动力学进行稳健高效的构型分子取样》，J.Chem。物理。，138，第174102条pp.（2013）
[47] 莱姆库勒，B。；Matthews，C.，《分子动力学：用确定性和随机数值方法》（2015），Springer·Zbl 1351.82001号
[48] 塞拉诺，M。；De Fabritis，G。；埃斯帕诺，P。；Coveney，P.V.，耗散粒子动力学的随机Trotter积分方案，数学。计算。模拟。，72, 2, 190-194 (2006) ·Zbl 1102.82022号
[49] 克劳登，体育。；Platen，E.，随机微分方程的数值解（1992），施普林格·Zbl 0925.65261号
[50] 北卡罗来纳州Bou-Rabee。；Owhadi，H.，随机背景下变分积分器的长运行精度，SIAM J.Numer。分析。，48, 1, 278-297 (2010) ·Zbl 1215.65012号
[51] 塔莱，D。；Tubaro，L.，求解随机微分方程数值格式的全局误差展开，Stoch。分析。申请。，8，483-509（1990年）·Zbl 0718.60058号
[52] 德布斯切，A。；Faou，E.，SDE的弱向后误差分析，SIAM J.Numer。分析。，50, 3, 1735-1752 (2012) ·Zbl 1256.65002号
[53] 阿卜杜勒。；Vilmart，G。；Zygalakis，K.C.，遍历SDE不变测度的高阶数值逼近，SIAM J.Numer。分析。，52, 4, 1600-1622 (2014) ·Zbl 1310.65007号
[54] 阿卜杜勒。；Vilmart，G。；Zygalakis，K.C.，朗之万动力学的李-罗特分裂方法的长期精度，SIAM J.Numer。分析。，53, 1, 1-16 (2015) ·Zbl 1327.65015号
[55] Givon，D。；Kupferman，R。；Stuart，A.，《提取宏观动力学：模型问题和算法》，非线性，17，6，R55-R127（2004）·Zbl 1073.82038号
[56] Verlet，L.，经典流体的计算机“实验”。I.Lennard-Jones分子的热力学性质，物理学。修订版，159、1、98（1967）
[57] 雅各布森，A.F。；O.G.穆里森。；Besold，G.，《粗颗粒模型脂质双层动力学模拟中的人工制品》，J.Chem。物理。，第122、20条，第204901条，pp.（2005）
[58] Rugh，H.H.，《温度的动力学方法》，Phys。修订稿。，78, 5, 772 (1997)
[59] 布拉加，C。；Travis，K.P.，《结构温度Nosé-Hoover恒温器》，J.Chem。物理。，123，第134101条pp.（2005）
[60] Allen，M.P.，《利用耗散粒子动力学进行膜模拟的结构温度》，J.Phys。化学。B、 110、8、3823-3830（2006）
[61] 特拉维斯，K.P。；Braga，C.，《原子和分子系统的配置温度控制》，J.Chem。物理。，128，第014111条pp.（2008）
[62] 俄克拉荷马州Farago。；Grönbech-Jensen，N.，《耗散粒子动力学与Itó-Stratonovich困境之间的联系》，J.Chem。物理。，144，8，第084102条pp.（2016）
[63] 唐，S。；Evans，G.T。；梅森，C.P。；Allen，M.P.，《坚硬椭球体流体的剪切粘度：动力学理论和分子动力学研究》，J.Chem。物理。，102, 9, 3794-3811 (1995)
[64] Hansen，J.-P。；McDonald，I.，《简单液体理论》（2006），学术出版社
[65] Ashurst，W.T。；胡佛，W.G.，《通过非平衡分子动力学实现的密度流体剪切粘度》，《物理学》。修订版A，11，2，658-678（1975）
[66] 奇科蒂，G。；Tenenbaum，A.，分子动力学标准系综和非平衡态，J.Stat.Phys。，23, 6, 767-772 (1980)
[67] 特罗齐，C。；Ciccotti，G.，分子动力学中的定态非平衡态。二、。牛顿定律，物理学。版本A，29，2，916-925（1984）
[68] Joubaud，R。；Stoltz，G.，用Langevin动力学计算非平衡剪切粘度，多尺度模型。模拟。，10, 1, 191-216 (2012) ·Zbl 1246.82092号
[69] A.Lees。；Edwards，S.，《极端条件下运输过程的计算机研究》，J.Phys。C、 固态物理。，5, 15, 1921 (1972)
[70] Fedosov，D.A。；潘，W。；卡斯维尔，B。；Gompper，G。；Karniadakis，G.E.，预测人类血液在硅胶中的粘度，Proc。国家。阿卡德。科学。，108, 29, 11772-11777 (2011)
[71] Boek，E.S。；Coveney，P.V.公司。；Lekkerkerker，H.N.W.，用耗散粒子动力学模拟稠密胶体悬浮液中的流变现象，J.Phys。康登斯。Matter，8，47，9509-9512（1996）
[72] Boek，E.S。；Coveney，P.V.公司。；Lekkerkerker，H.N.W。；van der Schoot，P.，使用耗散粒子动力学模拟致密胶体悬浮液的流变学，Phys。E版，55、3、3124（1997）
[73] Fedosov，D.A。；Karniadakis，G.E。；Caswell，B.，聚合物流体在反向Poiseuille流中耗散粒子动力学模型的稳态剪切流变仪，J.Chem。物理。，第132、14条，第144103页（2010年）
[74] 勒文加，M。；苏尼加，I。；Español，P.，耗散粒子动力学中的边界条件，计算。物理学。社区。，121-122（1999年）
[75] 皮夫金，I.V。；Karniadakis，G.E.，在耗散粒子动力学中施加无滑移边界条件的新方法，J.Compute。物理。，207, 1, 114-128 (2005) ·Zbl 1177.76325号
[76] 皮夫金，I.V。；卡斯维尔，B。；Karniadakis，G.，耗散粒子动力学，Rev.Compute。化学。，27, 85-110 (2010)
[77] Chatterjee，A.，高耗散率耗散粒子动力学模拟的Lees-Edwards周期边界条件修正，Mol.Simul。，33233-1233-1236（2007年）
[78] 惠特尔，M。；Travis，K.P.，《用核心修正耗散粒子动力学对胶体进行动态模拟》，J.Chem。物理。，第132、12条，第124906页（2010年）
[79] 德霍姆梅尔，J。；Evans，D.J.，密闭流体中的结构温度曲线。I.原子流体，J.Chem。物理。，114, 14, 6229 (2001)
[80] 德霍姆梅尔，J。；Evans，D.J.，密闭流体中的结构温度曲线。二、。分子流体，J.Chem。物理。，114, 14, 6236 (2001)
[81] 卢，L。；杰普斯，O.G。；德霍姆梅尔，J。；Evans，D.J.，分子系统的配置恒温器，分子物理学。，100, 14, 2387-2395 (2002)
[82] 德霍姆梅尔，J。；Evans，D.J.，《剪切流流体的温度恒温器配置：液氯应用》，摩尔物理学。，99, 21, 1825-1829 (2001)
[83] 德霍姆梅尔，J。；Evans，D.J.，构型温度和分子动力学温度恒温器之间的对应，J.Chem。物理。，117, 13, 6016 (2002)
[84] Junghans，C。；普拉普特尼克，M。；Kremer，K.，恒温器控制的传输特性：扩展耗散粒子动力学恒温器，软物质，4，1156-161（2008）
[85] 拉伦佐斯，J.P。；布伦南，J.K。；摩尔，J.D。；Lísal，M。；Mattson，W.D.，使用类Shardlow分裂算法并行实现等温和等能量耗散粒子动力学，计算。物理学。社区。，185, 7, 1987-1998 (2014) ·Zbl 1352.65024号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。