×

多时间尺度化学动力学系统的嵌套随机模拟算法。 (英语) Zbl 1162.80003号

摘要:我们提出了一种有效的数值算法来模拟多时间尺度的化学动力学系统。该算法是对传统随机模拟算法(SSA)的改进,也称为Gillespie算法。它以嵌套SSA的形式,使用外部SSA来模拟慢反应,其速率由模拟快速反应的内部SSA的实现计算得出。该算法本身是非常通用和无缝的,它相当于对原始SSA的一个小修改。我们对这种多尺度化学动力学系统的分析使我们能够识别系统中的慢变量,推导慢时间尺度上的有效动力学,并为嵌套SSA提供误差估计。利用这些误差估计讨论了嵌套SSA的效率,并通过几个数值例子进行了说明。

MSC公司:

80A30型 热力学和传热中的化学动力学
80平方米 其他数值方法(热力学)(MSC2010)
65二氧化碳 蒙特卡罗方法
92E20型 化学中的经典流动、反应等

软件:

罗德斯
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Bensoussan,A。;狮子,J.-L。;Papanicolaou,G.C.,周期结构的渐近分析,(《数学及其应用研究》,第5卷(1978年),爱思唯尔,北荷兰:爱思唯尔,北荷兰纽约)·Zbl 0411.60078号
[2] Bortz,A.B。;卡洛斯,M.H。;Lebowitz,J.L.,伊辛自旋系统蒙特卡罗模拟的新算法,J.Comp。物理。,17, 10-18 (1975)
[3] 曹毅。;Gillespie,D。;Petzold,L.,《慢规模随机模拟算法》,J.Chem。物理。,122, 014116 (2005)
[4] 曹毅。;Gillespie,D。;Petzold,L.,化学反应系统随机部分平衡假设的多尺度随机模拟算法,J.Comp。物理。,206, 395-411 (2005) ·Zbl 1088.80004号
[5] W.E,B.Engquist,X.Li,W.Ren,E.Vanden-Eijnden,异质多尺度方法:综述,预印本,2005。;W.E,B.Engquist,X.Li,W.Ren,E.Vanden-Eijnden,《异质多尺度方法:综述》,预印本,2005年。
[6] E、 W。;刘,D。;Vanden-Eijnden,E.,随机微分方程的多尺度方法分析,Comm.Pure Appl。数学。,58, 1544-1585 (2005) ·Zbl 1080.60060号
[7] W.E,D.Liu,E.Vanden-Eijnden,不同速率化学动力学系统的嵌套随机模拟算法,J.Chem。物理学。123 (2005) 194107.; W.E,D.Liu,E.Vanden-Eijnden,不同速率化学动力学系统的嵌套随机模拟算法,J.Chem。物理学。123 (2005) 194107.
[8] Engquist,B。;Tsai,Y.-H.,一类刚性常微分方程的异质多尺度方法,数学。公司。,74, 1707-1742 (2005) ·Zbl 1081.65075号
[9] Fedoroff,N。;Fontana,W.,《大分子的小数目》,《科学》,2971129-1131(2002)
[10] Gillespie,D.T.,《数值模拟耦合化学反应随机时间演化的通用方法》,J.Comp。物理。,22, 403-434 (1976)
[11] Gillespie,D.T.,耦合化学反应的精确随机模拟,物理学杂志。化学。,81, 2340-2361 (1977)
[12] 海尔,E。;Wanner,G.,《求解常微分方程II:刚性和微分代数问题》。求解常微分方程II:刚性和微分代数问题,计算数学中的Springer级数(2004),Springer·Zbl 0729.65051号
[13] Haseltine,E.L。;罗林斯,J.B.,随机动力学快速和缓慢耦合反应的近似模拟,J.化学。物理。,117, 6959-6969 (2002)
[14] Khasminskii,R.Z.,《关于小参数微分方程定义的随机过程》,理论探索。申请。,11, 211-228 (1966) ·Zbl 0168.16002号
[15] Khasminski,R.Z.,具有随机右手边的微分方程解的极限定理,理论问题。申请。,11, 390-406 (1966)
[16] Khasminskii,R.Z。;尹,G。;Zhang,Q.,构造具有弱和强相互作用的马尔可夫链概率分布的渐近级数,Quart。申请。数学。,55, 177-200 (1997) ·兹比尔0870.34058
[17] Kurtz,T.G.,扰动算子半群的极限定理及其在随机演化中的应用,J.Funct。分析。,第12页,55-67页(1973年)·Zbl 0246.47053号
[18] Meyn,S.P。;Tweedie,R.L.,马尔可夫过程的稳定性,I,II,和III,高级应用。探针。,24542-574(1992)和25518-5481993·Zbl 0757.60061号
[19] Papanicolaou,G.C.,《随机微分方程渐近分析导论》(DiPrima,R.C.,应用数学讲座,第16卷(1977年),美国数学学会:美国数学学会普罗维登斯,R.I.)·Zbl 0443.60049号
[20] C.A.Petri,Kommunikation mit Automaten,Institut fur Instru-mentelle Mathematik,Bonn,Schriften des IIM Nr.21962。;C.A.Petri,Kommunikation mit Automaten,Institut fur Instru-mentelle Mathematik,Bonn,Schriften des IIM Nr.21962。
[21] 现时,R.D.,气体动力学理论(1958),麦格劳·希尔(McGraw-Hill),(第8章)·Zbl 0065.19602号
[22] Rao,C.V。;Arkin,A.P.,《随机化学动力学和准静态假设:Gillespie算法的应用》,J.Chem。物理。,118, 4999-5010 (2003)
[23] Samant,A。;Vlachos,D.G.,《克服部分平衡随机模拟中的刚度:多尺度蒙特卡罗算法》,J.Chem。物理。,123, 144114 (2005)
[24] 斯利瓦斯塔瓦,R。;你,L。;萨默斯,J。;Yin,J.,细胞内病毒动力学的随机与确定性建模,J.Theor。《生物学》,218309-321(2002)
[25] 高桥,K。;Kaizu,K。;胡,B。;Tomita,M.,《细胞模拟的多算法、多时间尺度方法》,生物信息学,20538-546(2004)
[26] Vanden-Eijnden,E.,《具有随机效应的多尺度动力系统的数值技术》,Comm.Math。科学。,1, 385-391 (2003) ·Zbl 1088.60060号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。