×

带电荷和偶极子的静电系统的基于NFFT的Ewald求和。 (英语) Zbl 1376.78002号

摘要:带电粒子系统中库仑相互作用的有效计算在分子动力学模拟领域具有重要意义。众所周知,基于Ewald求和方法和快速傅里叶变换(FFT)可以实现近似。在本文中,我们考虑了含有点电荷和点偶极子混合物的粒子系统。推导了电荷-偶极相互作用Ewald求和公式中的新截止误差,并通过数值例子进行了验证。此外,我们首次提出了一种基于非等间距数据FFT(NFFT)的计算混合电荷偶极相互作用的(mathcal{O}(N\log N))粒子网格算法。我们给出了电荷偶极子系统的第一个数值结果,表明所引入的方法可以调整到高精度,并验证了(mathcal{O}(N\log N))标度。为了有效地计算与偶极子的相互作用,推导并实现了NFFT的两个新变体,即Hessian NFFT和伴随梯度NFFT。在NFFT的背景下,这些新变体本身就非常重要。本文提出的粒子网格方法是粒子-粒子NFFT(P^{2})框架的扩展。因此,所有先前推导的NFFT特征,包括周期和非周期边界条件的任意组合的处理,混合电荷偶极子和纯偶极子系统现在也支持处理三斜盒子形状和大规模并行实现。这些算法作为ScaFaCoS库的一部分公开提供。

MSC公司:

78A30型 静电和磁力静力学
78A35型 带电粒子的运动
65T50型 离散和快速傅里叶变换的数值方法
78A50型 光学和电磁理论中的天线、波导
81V70型 多体理论;量子霍尔效应
78M25型 光学中的数值方法(MSC2010)
65兰特 积分变换的数值方法
81V55型 分子物理学
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] af Klinteberg,L。;Shamshirgar,D.S。;Tornberg,A.-K.,自由空间Stokes势的快速Ewald求和,arXiv电子版·Zbl 1360.65295号
[2] 阿诺德,A。;博尔顿,M。;Dachsel,H。;Fahrenberger,F。;Gähler,F。;哈尔弗,R。;Heber,F。;霍夫曼,M。;伊瑟林豪森,J。;卡巴德肖,I。;O.Lenz。;Pippig,M.,ScaFaCoS-可扩展快速Coloumb解算器
[3] 阿诺德,A。;Fahrenberger,F。;霍尔姆,C。;O.Lenz。;波尔滕,M。;Dachsel,H。;哈尔弗,R。;卡巴德肖,I。;Gähler,F。;Heber,F。;Iseringhausen,J。;霍夫曼,M。;Pippig,M。;波茨,D。;Sutmann,G.,《远程交互可扩展快速方法的比较》,Phys。E版,88,第063308条pp.(2013)
[4] Beylkin,G.,关于奇异函数的快速傅里叶变换,应用。计算。哈蒙。分析。,2, 363-381 (1995) ·Zbl 0838.65142号
[5] Bolten,M.,《结构化网格的多重网格方法及其在粒子模拟中的应用》
[6] Cerdá,J.J。;巴伦格,V。;Holm,C.,《偶极相互作用的粒子-粒子-网格方法:关于分析微分和网格交错方案的误差估计和效率》,J.Chem。物理。,135,第184110条pp.(2011)·Zbl 1250.78049号
[7] Cerdá,J.J。;巴伦格,V。;O.Lenz。;Holm,C.,《偶极相互作用的P3M算法》,J.Chem。物理。,129,第234104条pp.(2008)
[8] 达顿,T。;约克·D·。;Pedersen,L.,《粒子网格Ewald:大系统中Ewald和的(n(n))方法》,J.Chem。物理。,98, 10089-10092 (1993)
[9] 德列乌,S.W。;佩拉姆,J.W。;Smith,E.R.,《周期边界条件下静电系统的模拟》。I.晶格和和介电常数,Proc。英国皇家学会。,序列号。A、 数学。物理学。工程科学。,373, 27-56 (1980)
[10] de Leeuw,S.W。;佩拉姆,J.W。;Smith,E.R.,《周期边界条件下静电系统的模拟》。I.晶格和和介电常数,Proc。R.Soc.伦敦。序列号。A、 373、27-56(1980)
[11] Deserno,M。;Holm,C.,《如何对Ewald和进行网格划分》。I.各种粒子网格程序的理论和数值比较,J.Chem。物理。,109, 7678-7693 (1998)
[12] 德塞雷诺,M。;Holm,C.,《如何对Ewald和进行网格划分》。二、。颗粒-颗粒-颗粒-mesh算法的准确误差估计,J.Chem。物理。,109, 7694-7701 (1998)
[13] Duijndam,A.J.W。;Schoneville,M.A.,非均匀快速傅里叶变换,地球物理学,64,539-551(1999)
[14] Dutt,A。;Rokhlin,V.,非等间距数据的快速傅里叶变换,SIAM J.Sci。统计计算。,14, 1368-1393 (1993) ·Zbl 0791.65108号
[15] 伊斯特伍德,J。;霍克尼,R。;Lawrence,D.,P3M3DP-三维周期颗粒-颗粒/颗粒-颗粒程序,计算。物理学。社区。,19, 2, 215-261 (1980)
[16] 埃斯曼,美国。;佩雷拉,L。;Berkowitz,M.L。;达顿,T。;Lee,H。;Pedersen,L.G.,《光滑颗粒网格Ewald方法》,J.Chem。物理。,103, 8577-8593 (1995)
[17] Ewald,P.P.、Die Berechnung optischer und elektrostatischer Gitterpotentiale、Ann.Phys.、。,369, 253-287 (1921)
[18] Fourmont,K.,非等间距快速傅立叶变换及其在断层扫描中的应用,J.傅立叶分析。申请。,9, 431-450 (2003) ·Zbl 1073.65151号
[19] Greengard,L。;Lee,J.-Y.,《加速非均匀快速傅里叶变换》,SIAM Rev.,46,443-454(2004)·Zbl 1064.65156号
[20] Greengard,L。;Rokhlin,V.,《粒子模拟的快速算法》,J.Compute。物理。,73, 325-348 (1987) ·Zbl 0629.65005号
[21] Hedman,F。;Laaksonen,A.,基于非均匀快速傅立叶变换的Ewald求和,化学。物理学。莱特。,425, 142-147 (2006)
[22] 霍克尼,R.W。;Eastwood,J.W.,《使用粒子的计算机模拟》(1988),Taylor&Francis公司:Taylor and Francis,Inc.,美国宾夕法尼亚州布里斯托尔·Zbl 0662.76002号
[23] 霍夫曼,M。;Rünger,G.,耦合并行粒子码的有效数据再分配方法,(第42届国际并行处理会议(ICPP-2013)(2013),IEEE),40-49
[24] J.I.杰克逊。;Meyer,C.H。;西村,D.G。;Macovski,A.,使用网格选择傅里叶反演的卷积函数,IEEE Trans。医学成像,10473-478(1991)
[25] Kabadshow,I.,《周期边界条件和误差控制快速多极方法》(2012年),Bergische Universityät Wuppertal:Bergiche Universität Wuppertal Jülich博士论文
[26] Kaiser,J.F.,《数字滤波器》(Kuo,F.F.;Kaise,J.F,《数字计算机系统分析》(1966),威利出版社:威利纽约)·Zbl 0171.37803号
[27] Keiner,J。;Kunis,S。;Potts,D.,使用NFFT3-一个用于各种非等间距快速傅里叶变换的软件库,ACM Trans。数学。软质。,36,第19条pp.(2009)·Zbl 1364.65303号
[28] Khoromskaia,V。;Khoromskij,B.N.,缺陷三维晶格上长程势求和的快速张量方法,Numer。线性代数应用。,23, 249-271 (2016) ·Zbl 1413.65141号
[29] Kolafa,J。;Perram,J.W.,点电荷系统Ewald求和公式中的截止误差,摩尔模拟。,9, 5, 351-368 (1992)
[30] Nestler,F.,一种基于NFFT的方法,用于在各种周期边界条件下高效计算偶极-偶极相互作用,Appl。数字。数学。,105, 25-46 (2016) ·兹比尔1416.65568
[31] Nestler,F.,基于非等间距FFT的RMS误差的自动参数调整,高级计算。数学。,42, 4, 889-919 (2016) ·Zbl 1348.65184号
[32] Nestler,F.,基于NFFT的快速Ewald求和的参数调整,Front。物理。,4,第28条pp.(2016)
[33] 雀巢,F。;Pippig,M。;Potts,D.,基于混合周期NFFT的快速Ewald求和,J.Compute。物理。,285280-315(2015年)·Zbl 1351.78066号
[34] 雀巢,F。;Pippig,M。;Potts,D.,基于NFFT的各种类型周期边界条件的快速Ewald求和,(Sutmann,G.;Grotendorst,J.;Gompper,G.,Marx,D.,《液体溶化和传输的计算趋势》,IAS系列,第28卷(2015年),Forschungszentrum Jülich股份有限公司:Forschungszentrun Jülch股份有限责任公司Jünich),575-598
[35] Pippig,M.,PFFT——FFTW到大规模并行体系结构的扩展,SIAM J.Sci。计算。,35、C213-C236(2013)·Zbl 1275.65098号
[36] Pippig,M.,《大规模并行、快速傅里叶变换和粒子网格方法》(2015),Chemnitz大学,论文
[37] Pippig,M.,PNFFT-并行非等间距FFT软件库
[38] Pippig,M.,PFFT-并行FFT软件库
[39] Pippig,M。;Potts,D.,基于非等间距快速傅里叶变换的粒子模拟,(Sutmann,G.;Gibbon,P.;Lippert,T.,《复杂系统中长程相互作用的快速方法》,IAS系列,第6卷(2011),Forschungszentrum Jülich GmbH:Forschungszentrum-Jülch GmbHJünich),131-158
[40] Pippig,M。;Potts,D.,《并行三维非等间距快速傅里叶变换及其在粒子模拟中的应用》,SIAM J.Sci。计算。,35、C411-C437(2013)·Zbl 1362.65141号
[41] 丹尼尔·波茨(Daniel Potts);Steidl,Gabriele,NFFTs非等间距结的快速求和,SIAM J.Sci。计算。,24, 2013-2037 (2003) ·Zbl 1040.65110号
[42] 波茨,D。;斯泰德尔,G。;Tasche,M.,《非等间距数据的快速傅里叶变换:教程》(Benedetto,J.J.;Ferreira,P.J.S.G.,《现代采样理论:数学与应用》(2001),Birkhäuser:Birkháuser Boston,MA,USA),247-270
[43] Steidl,G.,《非等间距网格的快速傅里叶变换注释》,《高级计算》。数学。,9, 337-353 (1998) ·Zbl 0917.65123号
[44] 托恩伯格,A.-K。;Greengard,L.,三维Stokes方程的快速多极方法,J.Compute。物理。,2271613-1619(2008年)·Zbl 1290.76116号
[45] 美国特罗滕贝格。;Oosterlee,C.W。;Schuller,A.,Multigrid(2000),学术出版社:学术出版社,美国佛罗里达州奥兰多
[46] 王振伟。;Holm,C.,偶极系统Ewald求和中截止误差的估计,J.Chem。物理。,115, 6277-6798 (2001)
[47] Ware,A.F.,不规则间隔数据的快速近似傅里叶变换,SIAM Rev.,40838-856(1998)·兹比尔0917.65122
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。