×
艾哈迈德·卡杜拉;孙淑玉;阿姆加德萨拉马

通过重加权和重建马尔可夫链加速蒙特卡罗分子模拟：对不同温度和密度条件下的正则系综平均值和二阶导数的外推。 （英语） Zbl 1349.65015号

J.计算。物理学。 270, 70-85 (2014).
小结：准确测定石油储层流体的热力学性质在许多应用中具有重要意义，特别是在石油工程和化学工程中。分子模拟具有许多吸引人的特点，特别是它需要较少的可调参数但却具有更好的预测能力；然而众所周知，与状态方程方法相比，分子模拟的CPU开销很大。最近，我们引入了一种高效的热力学一致性技术，以从使用昂贵的经典模拟预先计算的现有数据点，在不同热力学条件下快速再生蒙特卡罗-马尔可夫链（MCMC）。这种技术可以将模拟速度提高一百多万倍，使再生分子模拟的速度几乎与状态方程接近的速度一样快。本文首先简要回顾了这一技术，然后对其在不同相邻热力学条件下预测原始模拟MCMCs初级量系综平均值的能力进行了数值研究。此外，这种外推技术被扩展到预测二阶导数性质（例如热容和流体压缩性）。该方法通过对Lennard-Jones粒子正则系综中生成的MCMC进行重新加权和重构来实现。本文从原始模拟点外推了系统的势能、压力、等容热容和等温压缩系数沿等容线、等温线以及温度和密度变化的路径。最后，针对甲烷、氮和一氧化碳的单站点模型，提出了一组优化的Lennard-Jones参数（varepsilon，sigma）。

MSC公司：

65二氧化碳 蒙特卡罗方法
60J22型 马尔可夫链中的计算方法
80A10号 经典热力学和相对论热力学

关键词：

分子模拟;蒙特卡洛;正则系综;伦纳德-琼斯模型;马尔可夫链重加权与重构
PDF格式BibTeX公司 XML格式引用
\文本{A.Kadoura}等人，J.Comput。物理学。270、70-85（2014年；Zbl 1349.65015）
全文： 内政部

参考文献：

[1] Crippen，G.M.，《能量嵌入构象分析》，J.Compute。化学。，3, 4, 471-476 (1982)
[2] Duane，S.，混合蒙特卡洛，物理。莱特。B、 195,2216-222（1987年）
[3] Wang，J.S。；Swendsen，R.，《蒙特卡罗模拟中的非通用临界动力学》，Phys。修订稿。，58, 2, 86-88 (1987)
[4] Lyubartsev，A.，《自由能蒙特卡罗计算的新方法：扩展系综方法》，J.Chem。物理。，96, 1776 (1992)
[5] E.马里纳里。；Parisi，G.，模拟回火：一种新的蒙特卡洛方案，Europhys。莱特。，19, 6, 451 (1992)
[6] Geyer，C.J。；汤普森，E.A.，《退火马尔可夫链蒙特卡罗及其在祖先推断中的应用》，美国统计协会，90，431，909-920（1995）·Zbl 0850.62834号
[7] Frenkel，D.，通过拒绝状态采样加速蒙特卡罗模拟，Proc。国家。阿卡德。科学。美国，101，51，17571-17575（2004）
[8] Cortés Morales，A.D.，《模拟协议对吉布斯集合蒙特卡罗模拟效率的影响》，Mol.Simul。，39, 14-15, 1135-1142 (2013)
[9] A.卡杜拉。；萨拉马，A。；Sun，S.，《加速蒙特卡罗分子模拟的保守和混合早期排斥方案方法》，摩尔物理学。(2014)
[10] Frenkel，D。；Smit，B.，《理解分子模拟：从算法到应用》（2001），学术出版社
[11] 费伦伯格，A.M。；Swendsen，R.H.，研究相变的新蒙特卡罗技术，物理学。修订稿。，63, 15, 1635-1638 (1989)
[12] Torrie，G.M。；Valleau，J.P.，《蒙特卡罗自由能估算中的非物理抽样分布：伞式抽样》，J.Compute。物理。，23, 2, 187-199 (1977)
[13] Kiyohara，K。；Gubbins，K.E。；Panagiotopoulos，A.Z.，极化Stockmayer流体的相位共存特性（1996）
[14] 康拉德，P。；De Pablo，J.，《柔性水模型直方图重加权技术的比较》，流体相平衡。，150，51-61（1998年）
[15] 埃林顿，J.R。；Panagiotopoulos，A.Z.，哈密顿标度大正则蒙特卡罗修正白金汉指数-6势的相平衡，J.Chem。物理。，109, 1093 (1998)
[16] 加齐，S.M。；Kanhere，D.，（Ga的热力学性质_｛27｝Si_3\)使用密度泛函分子动力学的簇，J.Phys。化学。A、 116、1、11-17（2011）
[17] Menzl，G。；Köfinger，J。；Dellago，C.，纳米孔膜中水的相变和孔间关联，Phys。修订稿。，109, 2, 20602 (2012)
[18] 托里，G。；Valleau，J.，通过伞形取样对相分离液体混合物的蒙特卡罗研究，J.Chem。物理。，66, 4, 1402-1408 (1977)
[19] Valleau，J.，《密度标度：统计力学中的新蒙特卡罗技术》，J.Compute。物理。，96, 1, 193-216 (1991)
[20] Valleau，J.P.，《库仑相变：密度标度蒙特卡罗》，J.Chem。物理。，95, 1, 584-589 (1991)
[21] Valleau，J.P.，《亚临界Lennard-Jonesium的密度标度蒙特卡罗研究》，J.Chem。物理。，99, 4718 (1993)
[22] Kiyohara，K.，热力学标度吉布斯系综蒙特卡罗：测定流体相共存性质的新方法，摩尔物理学。，89, 4, 965-974 (1996)
[23] Valleau，J.P.，《温度和密度标度蒙特卡罗：Lennard-Jonesium的方法学和正则热力学》，摩尔模拟。，31, 4, 223-253 (2005) ·Zbl 1071.80006号
[24] Valleau，J.P.，《温度和密度标度蒙特卡罗：Lennard-Jonesium的等温等压热力学》，摩尔模拟。，31, 4, 255-275 (2005) ·Zbl 1071.80007号
[25] Sun，S。；A.卡杜拉。；Salama，A.，《重加权和重构蒙特卡罗分子模拟数据以外推到不同温度和密度条件的有效方法》（Procedia Computer Science，Procedia计算机科学，国际计算科学会议，西班牙巴塞罗那（2013））
[26] Möller，D.，通过分子动力学模拟预测流体混合物的热力学性质：甲烷-乙烷，摩尔物理学。，75, 2, 363-378 (1992)
[27] 段，Z。；莫勒，N。；Weare，J.H.，超临界流体混合物的一般状态方程和混合物PVTX特性的分子动力学模拟，Geochim。科斯莫辛。《学报》，60，7，1209-1216（1996）
[28] Mognetti，B.，流体及其混合物的粗粒模型：相行为的蒙特卡罗研究与微扰理论和实验的比较，J.Chem。物理。，130, 4, 044101 (2009)
[29] 用于分子流体模拟的Avendano，C.（SAFT-\gamma）力场。1.二氧化碳的单站点粗粒度模型，J.Phys。化学。B、 115、38、11154-11169（2011）
[30] Li，J.等人。；Calo，V.M.，使用单粒子模型的马尔可夫链蒙特卡罗方法对流体混合物的相共存模拟，J.Compute。物理。，249233-248（2013年）
[31] Sandler，S.I.，《应用统计热力学导论》（2010），Wiley
[32] 塞兹曼，美国。；Wagner，W.，《压力高达1000 MPa时，从熔线到625 K范围内甲烷的新状态方程和热力学性质表》，J.Phys。化学。参考数据，20，6，1061-1155（1991）
[33] Wagner，W.，《温度为63.151至1000 K、压力为2200 MPa的氮气热力学性质的参考状态方程》，J.Phys。化学。参考数据，29，6（2000）
[34] McCarty，R.，《一氧化碳热物理性质的相关性》（1989），美国国家标准与技术研究所：美国国家标准与技术研究所，科罗拉多州博尔德
[35] 跨度，R。；Wagner，W.，《压力高达800 MPa时，覆盖流体区域的二氧化碳从三点温度到1100 K的新状态方程》，J.Phys。化学。参考数据，251509（1996）
[36] 朋友，D.G。；Ingham，H。；Fly，J.F.，乙烷的热物理性质，J.Phys。化学。参考数据，20275（1991）
[37] 宫本茂，H。；Watanabe，K.，《液相丙烷的热力学性质模型》，《国际热力学杂志》。，21, 5, 1045-1072 (2000)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。