×

缺陷三维晶格长程势求和的快速张量方法。 (英语) Zbl 1413.65141号

小结:在本文中,我们提出了一种基于秩结构张量表示的具有多缺陷和非矩形几何的三维晶格长程势的快速求和方法。这是我们最近使用正则张量分解对矩形(L×L×L)晶格上基于网格的静电势求和技术的一个重要推广,并通过传统方法产生了计算复杂性(O(L)而不是(O(L^3)。生成的晶格和计算为Tucker或正则表示,其方向向量由移位参考张量的生成向量的1D求和组合而成,一旦在3D边界框中的大(N倍N倍N)表示网格上进行预计算。缺陷的张量数值处理是通过对正则或Tucker格式的张量进行简单求和,以代数方式进行的。为了减少生成的势和张量秩的显著增加,基于广义约化高阶SVD格式应用了(varepsilon)秩约简过程。对于正则/Tucker张量之和的简化高阶奇异值分解近似,我们证明了边矩阵的丢弃奇异值在相对范数中的稳定误差界。所需的存储容量在一维网格大小(O(N))中呈线性扩展,而数值成本则由(O(NL))估算。该方法适用于一类一般的核函数,包括牛顿、斯莱特、Yukawa、Lennard-Jones和偶极相互作用的核函数。数值试验证实了所提出的张量求和方法的有效性;我们证明,在Matlab实现中,可以在几秒钟内计算出3D晶格上含有缺陷/杂质的数百万牛顿核的总和。张量方法在使用1D复杂度的张量算法的3D网格上表示的晶格势和的进一步泛函微积分中是有利的,如积分或微分。

MSC公司:

65楼30 其他矩阵算法(MSC2010)
15A69号 多线性代数,张量演算
2015财年65 矩阵特征值和特征向量的数值计算
82B20型 格系统(伊辛、二聚体、波茨等)和平衡统计力学中出现的图上系统
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Deserno,《如何计算Ewald和》。I.各种粒子网格程序的理论和数值比较,《化学物理杂志》109(18)pp 7678–(1998)·doi:10.1063/1.477414
[2] Pisani,CRYSCOR:周期系统哈特雷阻塞后处理程序,物理化学化学物理14页7615–(2012)·doi:10.1039/c2cp23927b
[3] Savostyanov,使用张量列形式主义对蛋白质大小自旋系统的精确NMP模拟,《物理评论B:凝聚物质和材料物理学》90 pp 085139-1-085139–(2014)·doi:10.1103/PhysRevB.90.085139
[4] Dolgov S Khoromskij BN Litvinenko A Matthies HG 2015年Tensor列车数据格式响应面计算
[5] 坎塞斯,具有局部缺陷和光谱污染的周期薛定谔算子,SIAM数值分析杂志50(6),第3016页–·Zbl 1259.49075号 ·数字对象标识代码:10.1137/10855545
[6] 坎塞斯,《材料科学中点缺陷的数学建模,应用科学中的数学模型和方法》,23页1795–(2013)·兹比尔1282.35035 ·doi:10.1142/S021820513500528
[7] Ehrlacher,晶体缺陷原子模拟的边界条件分析,电子版ArXiv:1306.5334,(2013)
[8] Losilla,周期系统引力和静电方法的直接方法,化学物理杂志132(2)pp 024102-1-024102–(2010)·数字对象标识代码:10.1063/1.3291027
[9] Luskin,基于能量的混合准连续介质方法的制定和优化,《应用力学和工程中的计算机方法》253,第160页–(2013)·Zbl 1297.74123号 ·doi:10.1016/j.cma.2012.09.07
[10] Beck,密度泛函理论中的实空间网格技术,《现代物理评论》72 pp 1041-1-1041–(2000)·doi:10.1103/RevModPhys.72.1041
[11] Zuzovski,用于计算真实空间网格上密度泛函理论中静电项的辅助网格方法,物理化学化学物理(2015)·doi:10.1039/C5CP01090J
[12] Ewald,Die Berechnung optische und elektrostatischer Gitterpotentiale,《物理学年鉴》369第253页–(1921)·JFM 48.0566.02号 ·doi:10.1002/和p.19213690304
[13] Darten,粒子网格Ewald:大系统中Ewald和的O(NlogN)方法,化学物理杂志98 pp 10089–(1993)·数字对象标识代码:10.1063/1.464397
[14] Hünenberger,分子模拟中计算静电相互作用的格-和方法,美国物理研究所(AIP)会议论文集492 pp 17–(1999)
[15] 库丁,用周期快速多极子方法重温无限晶格和,化学物理杂志121页2886–(2004)·doi:10.1063/1.1771634
[16] Pollock,《关于P3M、FMM和大周期库仑系统的Ewald方法的评论》,《计算机物理通信》95第93页–(1996)·兹比尔0921.65081 ·doi:10.1016/0010-4655(96)00043-4
[17] Toukmaji,《透视中的Ewald求和技术:一项调查》,《计算机物理通信95》第73页–(1996)·Zbl 0923.65090号 ·doi:10.1016/0010-4655(96)00016-1
[18] Greengard,粒子模拟的快速算法,《计算物理杂志》73 pp 325–(1987)·Zbl 0629.65005号 ·doi:10.1016/0021-991(87)90140-9
[19] Khoromskaia,通过组装秩结构张量近似的基于网格的静电势晶格求和,《计算机物理通信》185 pp 3162–(2014)·Zbl 1360.78016号 ·doi:10.1016/j.cpc.2014.08.015
[20] Bertoglio,Galerkin投影Newton/Yukawa核的基于低秩求积的张量近似,计算机物理通信183(4)pp 904–(2012)·Zbl 1308.65213号 ·doi:10.1016/j.cpc.2011.12.016
[21] Braess,用指数和逼近波函数的渐近性,逼近理论杂志83页93–(1995)·Zbl 0868.41012号 ·doi:10.1006/jath.1995.1110
[22] Gavrilyuk,一类算子值函数的数据解析近似,《计算数学》74 pp 681–(2005)·Zbl 1066.65060号 ·doi:10.1090/S0025-5718-04-01703-X
[23] Hackbusch,多维非局部算子的低秩Kronecker乘积近似。第一部分多元函数的可分离近似,计算76 pp 177–(2006)·Zbl 1087.65049号 ·doi:10.1007/s00607-005-0144-0
[24] Grasedyck,低阶张量近似技术的文献综述,arXiv:1302.7121v1(2013)·Zbl 1279.65045号
[25] Hackbusch,计算机科学与工程课堂讲稿,第237页–(2014年)
[26] Khoromskaia,量子化学中的张量数值方法:从Hartree-Fock到激发能,物理化学化学物理(2015)·doi:10.1039/C5CP01215E
[27] Kolda,张量分解与应用,SIAM Review 51(3)pp 455–(2009)·Zbl 1173.65029号 ·doi:10.1137/07070111X
[28] Oseledets,张量-应变分解,SIAM科学计算杂志33(5)pp 2295–(2011)·Zbl 1232.15018号 ·doi:10.1137/090752286
[29] Schollwöck,矩阵产品状态时代的密度矩阵重整化群,《物理学年鉴》326(1),第96页–(2011)·Zbl 1213.81178号 ·doi:10.1016/j.aop.2010.09.012
[30] Khoromskaia V Hartree-Fock方程的多层张量结构数值解柏林工业大学2010年学位论文https://opus4.kobv.de/opus4-tuberlin/files/2780/khoromskaia_edera.pdf
[31] Dolgov-SV张量乘积方法在高维动力学问题数值模拟中的应用德国博士论文2014http://nbn-resolution.de/urn:nbn:de:bsz:15-库科萨-151129
[32] Khoromskij,函数相关多维数组的多重网格张量近似,SIAM科学计算杂志31(4)pp 3002–(2009)·Zbl 1197.65215号 ·数字对象标识代码:10.1137/080730408
[33] Khoromskaia,《Fock算子在一般基础上的快速准确三维张量计算》,《计算机物理通信》183 pp 2392–(2012)·Zbl 1302.81009号 ·doi:10.1016/j.cpc.2012.06.007
[34] Khoromskaia,张量数值方法的Black-box Hartree-Fock解算器,应用数学中的计算机方法14(1),第89页–(2014)·Zbl 1285.65091号 ·doi:10.1515/cmam-2013-0023
[35] Litsarev M Oseledets I 2015年离子-原子碰撞建模中多维积分的快速低秩近似
[36] Rakhuba M Oseledets I基于网格的电子结构计算:张量分解方法2015
[37] De Lathauwer,多重线性奇异值分解,《SIAM矩阵分析与应用杂志》21页1253–(2000)·Zbl 0962.15005号 ·doi:10.1137/S0895479896305696
[38] Khoromskij,O(dlogN)-高维数值建模中N-d张量的量化近似,构造近似34 pp 257–(2011)·Zbl 1228.65069号 ·doi:10.1007/s00365-011-9131-1
[39] Khoromskij,经典势的低秩塔克张量近似,中欧数学杂志5(3)第1页–(2007)·Zbl 1130.65060号 ·doi:10.2478/s11533-007-0018-0
[40] 男孩,分子结构的自动基本计算,《自然》178第1207页–(1956)·数字对象标识代码:10.1038/1781207a0
[41] Braess,非线性近似理论(1986)·doi:10.1007/978-3-642-61609-9
[42] Stenger,基于Sinc和解析函数的数值方法(1993)·兹比尔0803.65141 ·doi:10.1007/978-1-4612-2706-9
[43] Khoromskij,关于Kohn-Sham方程格林迭代的张量近似,《科学计算与可视化》11,第259页–(2008)·doi:10.1007/s00791-008-0077-x
[44] Zeidler,牛津数学用户指南(2003)
[45] Khoromskij,结构化秩-(r1,…,rd)rd中函数相关张量的分解,应用数学中的计算机方法6(2)第194–(2006)页·Zbl 1120.65052号 ·doi:10.2478/cmam-2006-0010
[46] Kolda,正交张量分解,SIAM矩阵分析与应用杂志23 pp 243–(2001)·Zbl 1005.15020号 ·doi:10.1137/S0895479800368354
[47] 2014年格型和周期系统线性化Hartree-Fock方程的Khoromskia V Khoromskij BN张量方法
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。