×

存在模型误差时基于线性响应的参数估计。 (英语) Zbl 07506536号

摘要:最近,我们提出了一种基于线性响应统计的随机动力学参数估计方法。该方法基于一个非线性最小二乘问题,该问题考虑了波动损耗理论产生的响应特性。在本文中,我们讨论了模型错误时出现的一个重要问题。特别是,当平衡密度函数是高维和非高斯的,并且在某些情况下是未知的时,线性响应统计是不可访问的。我们表明,这个问题可以通过将不完美的模型拟合到适当的边际线性响应统计这可以使用可用数据和参数或非参数模型进行近似。参数估计方法的有效性在温度分布不均匀的分子动力学模型(朗之万动力学)和时空混沌的PDE(非朗之万力学)中得到了证明,其中模型误差是由严重的谱截断引起的。在这些例子中,我们展示了不完美模型,即使用所提方案估计参数的Langevin方程,如何在以下条件下预测潜在动力学的非线性响应统计容许外部干扰.

MSC公司:

60华氏度 随机分析
82立方厘米 时间相关统计力学(动态和非平衡)
65立方厘米 概率方法,随机微分方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 科恩,W。;Sham,L.J.,《包含交换和相关效应的自洽方程》,《物理学》。修订版,140,4A,A1133-A1138(1965)
[2] 伦格,E。;Gross,E.K.,《含时系统的密度-泛函理论》,Phys。修订版Lett。,52, 12, 997 (1984)
[3] Haydock,R。;海涅,V。;Kelly,M.,基于紧密结合带的局部原子环境的电子结构,J.Phys。C、 固态物理。,5, 20, 2845 (1972)
[4] 哈曼,D。;施吕特,M。;Chiang,C.,规范守恒赝势,Phys。修订版Lett。,43, 20, 1494 (1979)
[5] 托达,M。;库博,R。;Hashitsume,N.,《统计物理学II》。非平衡统计力学(1983),Springer
[6] 哈里姆·J。;李,X。;Zhang,H.,使用线性响应统计的参数估计方法,J.Stat.Phys。,168, 146-170 (2017) ·Zbl 1376.82030年
[7] 张,H。;李,X。;Harlim,J.,《使用线性响应统计的参数估计方法:数值方案》,Chaos,29,3,第033101页,(2019年)·兹比尔1462.62510
[8] Evans,D.J。;Morris,G.,《非平衡液体的统计力学》(2008),剑桥大学出版社·Zbl 1189.82064号
[9] 巴登,M。;Marrink,S.J.,蛋白质相互作用的粗粒模型,Curr。操作。结构。生物学,23878-886(2013)
[10] Leach,A.R.,《分子建模:原理与应用》(2001),培生教育
[11] Voth,G.A.,《凝聚相和生物分子系统的粗训练》(2008),CRC出版社
[12] Schlick,T.,《分子建模与模拟:跨学科指南》(2002),Springer-Verlag·Zbl 1011.92019
[13] 张,Z。;卢,L。;诺伊德,W.G。;克里希纳,V。;Pfaendtner,J。;Voth,G.A.,《定义大生物分子中粗粒位点的系统方法》,Biophys。J.,95,11,5073-5083(2008)
[14] Noid,W.,《透视:生物分子系统的粗粒度模型》,J.Chem。物理。,第139、9条,第09B201_1页(2013年)
[15] 诺伊德·W。;朱,J.-W。;艾顿,G.S。;克里希纳,V。;Izvekov,S。;Voth,G.A。;Das,A。;Andersen,H.C.,《多尺度粗粒度方法》。I.原子模型和粗粒度模型之间的严格桥梁,J.Chem。物理。,第128、24条,第244114页(2008年)
[16] Chipot,C。;Pohorille,A.,自由能计算(2007),Springer
[17] Dequidt,A。;Solano Canchaya,J.G.,粗颗粒耗散动力学模型的贝叶斯参数化,J.Chem。物理。,143,8,第084122条pp.(2015)
[18] Katsoulakis,医学硕士。;Plecháć,P.,非平衡扩展系统参数化粗粒度的信息理论工具,J.Chem。物理。,第139、7条,第074115页(2013年)
[19] Lei,H。;贝克,N.A。;Li,X.,广义Langevin方程的数据驱动参数化,Proc。国家。阿卡德。科学。,113, 50, 14183-14188 (2016) ·Zbl 1407.62320号
[20] 马,L。;李,X。;Liu,C.,使用降阶技术的粗粒度Langevin动力学,J.Compute。物理。,380, 170-190 (2019) ·兹比尔1451.65010
[21] 马,L。;李,X。;Liu,C.,《从朗之万动力学推导和近似粗颗粒动力学》,J.Chem。物理。,第145、20条,第204117页(2016年)
[22] 贝斯特,R.B。;李,B。;斯图尔特,A。;达格特,V。;Clarke,J.,非机械蛋白质能承受力吗?通过原子力显微镜和分子动力学模拟拉伸谷胱甘肽酶,生物物理学。J.,81,4,2344-2356(2001)
[23] 卢,H。;Schulten,K.,力诱导蛋白质结构域去折叠的Steered分子动力学模拟,蛋白质,结构。功能。生物信息。,35, 4, 453-463 (1999)
[24] Leith,C.,气候响应和波动耗散,J.Atmos。科学。,32, 10, 2022-2026 (1975)
[25] Pavliotis,G.A.,《随机过程和应用:扩散过程》,Fokker-Planck和Langevin方程,第60卷(2014),Springer·Zbl 1318.60003号
[26] Zwanzig,R.,《非平衡统计力学》(2001),牛津大学出版社·Zbl 1267.82001年
[27] 海尔,M。;Majda,A.J.,证明线性响应理论的简单框架,非线性,23,4,909(2010)·Zbl 1186.82006年
[28] Baiesi,M。;梅斯,C。;Wynents,B.,非平衡态的涨落和响应,物理学。修订版Lett。,第103、1条,第010602页(2009年)
[29] Majda,A.J。;齐,D.,通过统计响应和信息理论提高不完善湍流模型的预测技巧,J.非线性科学。,26, 1, 233-285 (2016) ·Zbl 1331.76096号
[30] 齐,D。;Majda,A.J.,通过经验信息理论预测具有不完美模型的被动标量湍流中的胖尾间歇性概率分布,Commun。数学。科学。,14, 6, 1687-1722 (2016) ·Zbl 1349.76101号
[31] Majda,A.J。;Qi,D.,预测复杂湍流动力系统统计响应和不确定性量化的降阶模型策略,SIAM Rev.,60,3,491-549(2018)·Zbl 1393.76041号
[32] Gottwald,G.A。;Harlim,J.,《加性和乘性噪声在过滤复杂动力系统中的作用》,Proc。R.Soc.A,数学。物理学。工程科学。,469、2155,文章20130096 pp.(2013)·Zbl 1371.37104号
[33] Pavliotis,G。;Stuart,A.,《多尺度方法:平均化和均匀化》(2008),Springer Science&Business Media·Zbl 1160.35006号
[34] Berry,T。;Harlim,J.,带模型误差的非线性多尺度系统滤波线性理论,Proc。R.Soc.A,数学。物理学。工程科学。,4702167,文章20140168第(2014)页·Zbl 1371.60068号
[35] Ishimori,Y.,一维Lennard-Jones晶格中的孤子,Prog。西奥。物理。,68, 2, 402-410 (1982) ·Zbl 1074.82512号
[36] 莱普里,S。;利维,R。;Politi,A.,经典低维晶格中的热传导,物理学。代表,377,1,1-80(2003)
[37] Dellago,C。;Bolhuis,P.G。;Chandler,D.,《有效过渡路径采样:Lennard-Jones簇重排的应用》,J.Chem。物理。,108, 22, 9236-9245 (1998)
[38] 萨克斯,M。;莱姆库勒,B。;Danos,V.,《变系数和非保守力的朗之万动力学:从稳态到数值方法》,《熵》,19,12,647(2017)
[39] Sweet,C.R。;Petrone,P。;潘德,V.S。;Izaguirre,J.A.,生物分子Langevin动力学的正常模式分配,J.Chem。物理。,128,14,第04B606条pp.(2008)
[40] Jaynes,E.,《信息理论和统计力学》,Phys。修订版,106、4、620(1957年)·Zbl 0084.43701号
[41] 郝伟(Hao,W.)。;Harlim,J.,解多维矩约束最大熵问题的方程组方法,Commun。申请。数学。计算。科学。,13, 2, 189-214 (2018) ·Zbl 1409.65034号
[42] Lyubartsev,美联社。;Laaksonen,A.,从径向分布函数计算有效相互作用势:反向蒙特卡罗方法,物理学。E版,52、4、3730(1995年)
[43] 安德森·D·F。;Mattingly,J.C.,一类随机微分方程的弱梯形方法,Commun。数学。科学。,9, 1, 301-318 (2011) ·Zbl 1288.60085号
[44] Muandet,K。;Fukumizu,K。;Sriperumbudur,B。;Schölkopf,B.,Kernel意指发行版的嵌入:回顾和超越,Found。趋势马赫数。学习。,10, 1-2, 1-141 (2017) ·Zbl 1380.68336号
[45] 江,S.W。;Harlim,J.,用数据驱动的非参数似然函数进行参数估计,熵,21,6,559(2019)
[46] 江南伟。;Harlim,J.,《缺失动力系统的建模:使用非参数框架推导参数模型》,Res.Math。科学。,7,3,1-25(2020),施普林格·Zbl 1447.37067号
[47] 张,H。;哈里姆·J。;Li,X.,使用正交多项式估计线性响应统计:RKHS公式,发现。数据科学。,2, 4, 443-485 (2020)
[48] Majda,A。;阿布拉莫夫,R。;Grote,M.,《多尺度非线性系统的信息理论和随机性》,CRM专题丛书,第25卷(2005),美国数学学会:美国数学学会普罗维登斯,美国罗德岛·Zbl 1082.60002
[49] Gritsun,A。;Branstator,G.,基于波动扩散定理的三维算子气候响应,J.Atmos。科学。,64, 7, 2558-2575 (2007)
[50] 拉奎,R.E。;Mahajan,S。;卢瑟福,P。;Tang,W.,陷阱离子模式的非线性饱和,物理学。修订版Lett。,34, 7, 391 (1975)
[51] Kuramoto,Y。;Tsuzuki,T.,远离热平衡的耗散介质中浓度波的持续传播,Prog。西奥。物理。,55, 2, 356-369 (1976)
[52] Stinis,P.,Kuramoto-Sivashinsky方程的随机最优预测,多尺度模型。模拟。,2, 4, 580-612 (2004) ·Zbl 1064.65006号
[53] 卡萨姆,A.-K。;Trefethen,L.N.,刚性脉冲的四阶时间步进,SIAM J.Sci。计算。,26, 4, 1214-1233 (2005) ·Zbl 1077.65105号
[54] Lu,F。;Lin,K.K。;Chorin,A.J.,Kuramoto-Sivashinsky方程的基于数据的随机模型简化,Physica D,340,46-57(2017)·Zbl 1376.35100号
[55] Toh,S.,描述Kuramoto-Sivashinsky方程时空混沌的局部结构统计模型,J.Phys。Soc.Jpn.公司。,56, 3, 949-962 (1987)
[56] 哈里姆·J。;江南伟。;Liang,S。;Yang,H.,用于具有缺失动力学的预测的机器学习,J.Comput。物理。,428,第109922条pp.(2021)·Zbl 07511409号
[57] Mehler,F.G.,Ueber die entwicklung einer function von beliebig vielen variablen nach laplaceschen functionen höherer ordnung,J.Reine Angew。数学。,66, 161-176 (1866)
[58] Xiu,D.,《随机计算的数值方法:谱方法》(2010),普林斯顿大学出版社·Zbl 1210.65002号
[59] 施穆克,M。;普拉达斯,M。;Pavliotis,G.A。;Kalliadasis,S.,广义Kuramoto-Sivashinsky方程的新模式简化策略,IMA J.Appl。数学。,80, 2, 273-301 (2013) ·Zbl 1328.35238号
[60] J.C.马丁利。;斯图尔特,A.M。;Higham,D.J.,sdes和近似的遍历性:局部Lipschitz向量场和退化噪声,Stoch。过程。申请。,101, 2, 185-232 (2002) ·Zbl 1075.60072号
[61] 吴,X。;Li,X.,关于具有多体原子间势的分子动力学模型的动量和能量通量的一致定义,模型。模拟。马特。科学。工程,23,1,文章015003 pp.(2015)
[62] 莱姆库勒,B。;Matthews,C.,《分子动力学》(2016),斯普林格出版社
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。