×
内森·沃恩;维克拉姆·加维尼;罗伯特·克拉斯尼

Kohn-Sham密度泛函理论的树型加速格林迭代。 （英语） Zbl 07506532号

J.计算。物理学。 430，文章ID 110101，22 p.（2021）.
摘要：我们提出了一种全电子Kohn-Sham密度泛函理论的实空间计算方法，称为树码加速格林迭代（TAGI）。TAGI基于Kohn-Sham方程的重新公式，在该公式中，微分形式的特征值问题通过与修改的亥姆霍兹-格林函数卷积，转换为积分形式的定点问题。在每个自洽场（SCF）迭代中，通过格林迭代计算不动点，其中离散卷积和通过GPU加速的重心拉格朗日树码有效计算。TAGI中使用的其他技术包括先验自适应网格细化、Fejér求积、奇异性减除、无梯度特征值更新和Anderson混合，以加速SCF和Green迭代的收敛。几个原子（Li，Be，O）和小分子（\（\mathrm）的基态能量计算{H} _2\)，一氧化碳，\（\mathrm{C_6H}_6\))证明TAGI能够有效实现化学精度。

引用于1文件

MSC公司：

65新元 偏微分方程边值问题的数值方法
65Dxx日 数值近似和计算几何（主要是算法）
78轴 光学和电磁理论的一般主题

关键词：

全电子Kohn-Sham密度泛函理论;实空间积分方程;绿色迭代;先验自适应网格细化;重心拉格朗日树码;GPU加速

软件：

Libxc公司;阿比尼特;量子浓缩咖啡;Matlab公司;高斯人;CheFSI公司;DFT-FE型;西北化工;CASTEP公司
PDF格式BibTeX公司 XML格式引用
\textit{N.Vaughn}等人，《计算杂志》。物理学。430，文章ID 110101，22 p.（2021；Zbl 07506532）
全文： 内政部 arXiv公司

参考文献：

[1] Challacombe，M。；施韦格勒，E。；Almlöf，J.，《库仑矩阵的快速组装：量子化学树代码》，J.Chem。物理。，104, 4685-4698 (1996)
[2] Challacombe，M。；Schwegler，E.，福克矩阵的线性标度计算，J.Chem。物理。，106, 5526-5536 (1997)
[3] Chew，W.C。；Michielssen，E。；宋，J.M。；Jin，J.M.，《计算电磁学中的快速高效算法》（2001），Artech House
[4] 科恩，O。；Kronik，L.公司。；Brandt，A.，全电子Kohn-Sham方程的局部精细多重网格解，J.Chem。理论计算。，9, 4744-4760 (2013)
[5] Natan，A。；A.本杰米尼。；纳维赫，D。；Kronik，L。；蒂亚戈，M.L。；贝克曼，S.P。；Chelikowsky，J.R.，一般周期和部分周期系统第一性原理计算的实空间赝势方法，Phys。B版，78，第075109条，pp.（2008）
[6] Hohenberg，P。；Kohn，W.，《非均匀电子气体》，Phys。修订版，136、3B、B864-B871（1964年）
[7] 科恩，W。；Sham，L.J.，《包含交换和相关效应的自洽方程》，《物理学》。修订版，140，4A，A1133-A1138（1965）
[8] Burke，K.，《密度泛函理论的观点》，J.Chem。物理。，第136、15条，第150901页（2012年）
[9] Jones，R.O.，《密度泛函理论：起源、崛起和未来》，Rev.Mod。物理。，87, 897-923 (2015)
[10] 北马尔迪罗西安。；Head-Gordon，M.，《计算化学中密度泛函理论的三十年：200个密度泛函的概述和广泛评估》，《分子物理学》。，115, 19, 2315-2372 (2017)
[11] 林，L。；卢，J。；Ying，L.，Kohn-Sham密度泛函理论的数值方法，数值学报。，28, 405-539 (2019) ·Zbl 07099162号
[12] Kresse，G。；Furthmüller，J.，使用平面波基集进行从头算总能量计算的高效迭代方案，Phys。B版，54、16、11169-11186（1996）
[13] Gonze，X。；J.-M.布肯。；加拉加斯，R。；Detraux，F。；Fuchs，M。；Rignanese，G.-M。；辛迪奇，L。；Verstraete，M。；谢拉，G。；Jollet，F。；托伦特，M。；罗伊，A。；Mikami，M。；Ghosez，P。；拉蒂，J.-Y。；Allan，D.，材料特性的第一原理计算：ABINIT软件项目，Compute。马特。科学。，25, 3, 478-492 (2002)
[14] 医学博士Segall。；Lindan，P.J.D。；Probert，M.J。；Pickard，C.J。；Hasnip，P.J。；克拉克·S·J。；Payne，M.C.，《第一原理模拟：想法、插图和CASTEP代码》，J.Phys。康登斯。Matter，14，11，2717-2744（2002）
[15] Giannozzi，P。；巴罗尼，S。；北卡罗来纳州博尼尼。；卡兰德拉，M。；R车。；卡瓦佐尼，C。；Ceresoli，D。；Chiarotti，G.L。；科科奇奥尼，M。；达博，I。；Dal Corso，A。；de Gironcoli，S。；Fabris，S。；弗雷特斯，G。；Gebauer，R。；格斯特曼，美国。；古古西斯，C。；Kokalj，A。；拉泽里，M。；马丁·萨莫斯。；Marzari，N。；Mauri，F。；马扎雷洛（Mazzarello，R.）。；Paolini，S。；Pasquarello，A。；Paulatto，L。；斯布雷西亚，C。；Scandolo，S。；Sclauzero，G。；Seitsonen，A.P。；Smogunov，A。；Umari，P。；Wentzcovitch，R.M.，Quantum ESPRESSO：材料量子模拟的模块化开源软件项目，J.Phys。康登斯。Matter，第21、39条，第395502页（2009年）
[16] Giannozzi，P。；安德烈西，O。；Brumme，T。；布诺，O。；Buongiorno Nardelli，M。；卡兰德拉，M。；R车。；卡瓦佐尼，C。；Ceresoli，D。；科科奇奥尼，M。；科隆纳，N。；卡尼梅奥，I。；Dal Corso，A。；de Gironcoli，S。；Delugas，P。；迪斯塔西奥，R.A。；费雷蒂，A。；Floris，A。；弗雷特斯，G。；Fugallo，G。；Gebauer，R。；格斯特曼，美国。；Giustino，F。；Gorni，T。；贾，J。；川村，M。；Ko，H.-Y。；Kokalj，A。；库苏克本利，E。；拉泽里，M。；马西里，M。；Marzari，N。；Mauri，F。；Nguyen，N.L。；Nguyen，H.-V。；Otero-de-la Roza，A。；Paulatto，L。；南卡罗来纳州蓬塞。；罗卡，D。；Sabatini，R。；桑特拉，B。；施利普夫，M。；Seitsonen，A.P。；Smogunov，A。；蒂姆罗夫，I。；Thonhauser，T。；Umari，P。；Vast，N。；吴，X。；Baroni，S.，《使用量子ESPRESSO进行材料建模的高级功能》，J.Phys。康登斯。Matter，第29、46条，第465901页（2017年）
[17] 劳克斯，T。；Slater，J.C.，《增强平面波法：电子结构计算指南》（Phys.Today，vol.20（1967）），92-93
[18] Andersen，O.K.，《能带理论中的线性方法》，物理学。B版，第12、8、3060-3083页（1975年）
[19] Wimmer，E。；Krakauer，H。；韦纳特，M。；Freeman，A.J.，计算分子和表面电子结构的全势自洽线性化增强平面波方法：{O} _2\)分子，物理。B版，24,2864-875（1981）
[20] 韦纳特，M。；Wimmer，E。；Freeman，A.J.，整体固体和表面的总能量全电子密度泛函方法，Phys。版本B，26，8，4571-4578（1982）
[21] Sjöstedt，E。；Nordström，L。；Singh，D.，增广平面波方法线性化的另一种方法，固态通信。，114, 1, 15-20 (2000)
[22] Madsen，G.K.H。；布拉哈，P。；施瓦兹，K。；Sjöstedt，E。；Nordström，L.，增广平面波方法的有效线性化，物理学。B版，64、19，第195134条，pp.（2001）
[23] 弗里希，M.J。；卡车，通用。；施莱格尔，H.B。；Scuseria，G.E。；Robb，医学硕士。；奇斯曼，J.R。；Scalmani，G。；巴龙，V。；彼得森，G.A。；Nakatsuji，H。；李，X。；Caricato，M。；Marenich，A.V。；布鲁诺，J。；Janesko，B.G。；Gomperts，R。；门努奇，B。；Hratchian，H.P。；奥尔蒂斯，J.V。；伊兹马约洛夫，A.F。；Sonnenberg，J.L。；Williams-Young，D。；丁·F。；Lipparini，F。；Egidi，F。；Goings，J。；彭，B。；Petrone，A。；亨德森，T。；拉纳辛格，D。；扎克泽夫斯基，V.G。；高杰。；雷加，N。；郑庚。；Liang，W。；哈达，M。；埃哈拉，M。；丰田，K。；福田，R。；长谷川，J。；石田，M。；Nakajima，T。；Y.本田。；O.北岛。；Nakai，H。；Vreven，T。；Throssell，K。；蒙哥马利，J.A。；Peralta，J.E。；奥利亚罗，F。；比尔帕克，M.J。；海德·J·J。；兄弟，E.N。；库丁，K.N。；Staroverov，V.N。；Keith，T.A。；小林，R。；诺曼德，J。；Raghavachari，K。；伦德尔，A.P。；伯兰特，J.C。；艾扬格，S.S。；托马西，J。；Cossi，M。；Millam，J.M。；克莱恩，M。；阿达莫，C。；卡米（Cammi），R。；Ochterski，J.W。；马丁·R·L。；莫罗库马，K。；O.Farkas。；Foresman，J.B。；Fox，D.J.（2016），Gaussian 16修订版B.01，Gaussian Inc.Wallingford CT
[24] 瓦利耶夫，M。；Bylaska，E。；戈文德，N。；科瓦尔斯基，K。；斯特拉茨马，T。；Van Dam，H。；王，D。；尼普洛查，J。；阿普拉，E。；温德斯，T。；de Jong，W.，NWChem：用于大规模分子模拟的全面且可扩展的开源解决方案，Compute。物理学。社区。，181, 9, 1477-1489 (2010) ·Zbl 1216.81179号
[25] 津田，E。；Tsukada，M.，基于自适应有限元方法的大尺度电子结构计算，J.Phys。Soc.Jpn.公司。，67, 11, 3844-3858 (1998)
[26] 帕斯克，J.E。；克莱因，B.M。；Fong，C.Y。；Sterne，P.A.，《固态电子结构计算的实空间局部多项式基础：有限元方法》，Phys。B版，59、19、12352-12358（1999）
[27] Lehtovaara，L。；哈夫，V。；Puska，M.，高阶有限元的全电子密度泛函理论和含时密度泛函，J.Chem。物理。，131, 5, 1-10 (2009)
[28] 莫塔马里，P。；Nowak，M。；Leiter，K。；克纳普，J。；Gavini，V.，Kohn-Sham密度泛函理论的高阶自适应有限元方法，J.Compute。物理。，253, 308-343 (2013) ·Zbl 1349.74331号
[29] Motamarri，P。；Gavini，V.，《使用谱有限元离散化进行大规模Kohn-Sham密度泛函理论计算的亚二次尺度子空间投影方法》，Phys。B版，90，11，第115127条pp.（2014）
[30] Kanungo，B。；Gavini，V.，使用丰富的有限元基础的大尺度全电子密度泛函理论计算，Phys。B版，95，3，第035112条pp.（2017）
[31] 莫塔马里，P。；达斯，S。；Rudraraju，S。；Ghosh，K。；达维多夫，D。；Gavini，V.，DFT-FE-用于大规模密度泛函理论计算的大规模并行自适应有限元代码，Compute。物理学。社区。，246，第106853条pp.（2020）·Zbl 07678420号
[32] 达斯，S。；Motamarri，P。；加维尼，V。；Turcksin，B。；李永伟。；Leback，B.，使用混合精度计算进行快速、可扩展和准确的基于有限元的从头计算：金属位错系统的46 PFLOPS模拟，（Int.Conf.High Performance.Compute.Networking，Storage Anal.SC（2019），IEEE Computer Society:IEEE Compute Socility New York，NY，USA），1-11
[33] Rokhlin，V.，《二维散射理论积分方程的快速求解》，J.Compute。物理。，86, 2, 414-439 (1990) ·Zbl 0686.65079号
[34] Jorgenson，R。；Mittra，R.，自由空间周期格林函数的有效计算，IEEE Trans。天线传播。，38, 5, 633-642 (1990) ·Zbl 0947.78526号
[35] 科伊夫曼，R。；Rokhlin，V。；Wandzura，S.，波动方程的快速多极方法：行人处方，IEEE天线提案。Mag.，35，3，7-12（1993）
[36] Medgyesi-Mitschang，L.N。；Putnam，J.M。；Gedera，M.B.，三维穿透散射体的广义矩方法，J.Opt。《美国社会》，第11、4、1383页（1994年）
[37] Bleszynski，E。；Bleszynski，M。；Jaroszewicz，T.，AIM:解决大规模电磁散射和辐射问题的自适应积分方法，无线电科学。，31, 5, 1225-1251 (1996)
[38] Volakis，J.L。；Sertel，K.，《电磁学积分方程方法》（2011），科学技术
[39] Botha，M.M.，《求解电磁散射的体积积分方程》，J.Compute。物理。，218, 1, 141-158 (2006) ·Zbl 1111.78012号
[40] Faddeyev，L.D。；Seckler，B.，《量子散射理论中的反问题》，J.Math。物理。，4, 1, 72-104 (1963) ·Zbl 0112.45101号
[41] 约翰逊，B.R。；Secrest，D.，《无交换非相对论性量子散射问题的解》，J.Math。物理。，2187-2195年7月12日（1966年）
[42] Alt、E。；格拉斯伯格，P。；Sandhas，W.，将三粒子碰撞问题简化为多通道双粒子Lippmann-Schwinger方程，Nucl。物理学。B、 2167-180（1967）
[43] 马塞尔，R.I。；Merrill，R.P。；Miller，W.H.，《正弦硬壁的量子散射：固体表面的原子衍射》，物理学。B版，第12、12、5545-5551页（1975年）
[44] Adhikari，S.K.，《二维量子散射》，美国物理杂志。，54, 4, 362-367 (1986)
[45] Hecht，K.T.，散射格林函数的算子形式和散射问题的积分方程，（散射（2000），纽约斯普林格出版社：纽约斯普林格出版社），477-480
[46] Kalos，M.H。；卡罗，蒙特，三体核和四体核基态的计算，物理学。修订版，128、4、1791-1795（1962）
[47] 赵，Z。；科夫瓦利，N。；Lin，W。；安，C.-H。；库奇曼，L。；Carin，L.，Schrödinger方程建模的快速多极体方法，J.Compute。物理。，224, 2, 941-955 (2007) ·Zbl 1114.65361号
[48] 哈里森·R·J。；Fann，G.I。；Yanai，T。；甘，Z。；Beylkin，G.，《多分辨率量子化学：基本理论和初步应用》，J.Chem。物理。，121, 23, 11587-11598 (2004)
[49] 哈里森·R·J。；Beylkin，G。；比肖夫，F.A。；Calvin，J.A。；Fann，G.I。；福索·坦德，J。；加林多，D。；哈蒙德，J.R。；Hartman-Baker，R。；希尔，J.C。；贾，J。；Kottmann，J.S。；欧，M.J。；裴，J。；Ratcliff，L.E。；路透社，M.G。；Richie-Halford，A.C.公司。；北卡罗来纳州罗梅罗。；Sekino，H。；西澳州谢尔顿。；Sundahl，B.E。；桑顿，W.S。；瓦列夫，E.F。；瓦兹奎兹·马亚戈提亚州。；Vence，N。；Yanai，T。；Yokoi，Y.，《疯狂：科学模拟的多分辨率自适应数值环境》，SIAM J.Sci。计算。，第38卷，S123-S142（2016），工业和应用数学学会出版物·Zbl 1365.65327号
[50] 莫伦坎普，M.J。；Young，T.，薛定谔方程的格林迭代收敛性，最近的高级计算。科学。，201-208 (2008) ·Zbl 1157.65507号
[51] Mohlenkamp，M.J.，多粒子薛定谔方程中单电子函数的函数空间要求，J.Math。物理。，第54、6条，第062105页（2013年）
[52] Khoromskij，B.N.，关于Kohn-Sham方程格林迭代的张量近似，计算。视觉。科学。，11, 4-6, 259-271 (2008) ·Zbl 1522.65249号
[53] 拉胡巴，M.V。；Oseledets，I.V.，《通过交叉近似实现低秩张量格式的快速多维卷积》，SIAM J.Sci。计算。，37、2、A565-A582（2015）·Zbl 1320.65197号
[54] 拉库巴，M.V。；Oseledets，I.V.，《基于网格的电子结构计算：张量分解方法》，J.Compute。物理。，312, 19-30 (2016) ·Zbl 1351.82040号
[55] Fejér，L.，Mechanische Quadratren mit positiven Cotesschen Zahlen，数学。Z.，37，1，287-309（1933）
[56] Trefethen，L.N.，高斯求积比克伦肖-库蒂斯好吗？，SIAM版本，50，1，67-87（2008）·Zbl 1141.65018号
[57] Rabinowitz，P.，高等分析的近似方法。L.V.Kantorovich和V.I.Krylov，《科学》，134、3487、1358（1961年），由柯蒂斯·D·本斯特（Curtis D.Benster）翻译自第三版俄语。Interscience，纽约，1959年
[58] Anselone，P.M.，积分方程数值解中的奇异减法，J.Aust。数学。Soc.序列号。B、 申请。数学。，22, 4, 408-418 (1981) ·Zbl 0477.65095号
[59] Wang，L。；Krasny，R。；Tlupova，S.，基于重心拉格朗日插值的核相关树码，Commun。计算。物理。，28, 1415-1436 (2020) ·Zbl 1473.65025号
[60] 沃恩，N。；Wilson，L。；Krasny，R.，GPU加速重心拉格朗日树码，（2020 IEEE国际并行和分布式处理研讨会研讨会，2020 IEEE International Parallel and Distributed Processing Symposium Workshops，IPDPSW，New Orleans，LA，USA（2020）），701-710
[61] Goedecker，S.，线性标度电子结构方法，修订版。物理。，71, 4, 1085-1123 (1999)
[62] Anderson，D.G.，非线性积分方程的迭代程序，J.ACM，12，4，547-560（1965）·Zbl 0149.11503号
[63] Ceperley，D.M。；Alder，B.J.，用随机方法研究电子气体的基态，物理学。修订稿。，45565-569（1980年）
[64] 佩杜，J.P。；Zunger，A.，多电子系统密度泛函近似的自相互作用校正，Phys。B版，23、10、5048-5079（1981年）
[65] 马奎斯，M.A.L。；Oliveira，M.J.T。；Burnus，T.，Libxc：密度泛函理论的交换和相关泛函库，计算。物理学。社区。，183, 10, 2272-2281 (2012)
[66] 莱托拉，S。；斯泰格曼，C。；奥利维拉，M.J。；Marques，M.A.，《libxc的最新发展——密度泛函理论泛函的综合库》，SoftwareX，7，1-5（2018）
[67] 切里科夫斯基，J.R。；Troullier，N。；Saad，Y.，《有限差分伪势法：无基电子结构计算》，Phys。修订稿。，72, 8, 1240-1243 (1994)
[68] Bernholc，J。；霍达克，M。；Lu，W.，《真实空间多重网格方法的最新发展和应用》，J.Phys。康登斯。Matter，第20、29条，第294205页（2008年）
[69] 萨阿德，Y。；切里科夫斯基，J.R。；Shontz，S.M.，材料电子结构计算的数值方法，SIAM Rev.，52，1，3-54（2010）·邮编：1185.82004
[70] 陈，H。；戴，X。；龚，X。；He，L。；Zhou，A.，Kohn Sham模型的自适应有限元近似，多尺度模型。模拟。，12, 4, 1828-1869 (2014) ·Zbl 1316.35260号
[71] Trefethen，L.N.，MATLAB中的光谱方法（2008），SIAM
[72] Trefethen，L.N。；Bau，D.I.，《数值线性代数》（2000），SIAM工业和应用数学学会·Zbl 0984.65019号
[73] 巴恩斯，J。；Hut，P.，《分层O（N log N）力计算算法》，《自然》，324，6096，446-449（1986）
[74] Greengard，L。；Rokhlin，V.，《粒子模拟的快速算法》，J.Compute。物理。，73, 2, 325-348 (1987) ·Zbl 0629.65005号
[75] 贝鲁特，J.-P。；Trefethen，L.N.，重心拉格朗日插值，SIAM Rev.，46，3，501-517（2004）·兹比尔1061.65006
[76] 李，P。；约翰斯顿，H。；Krasny，R.，《屏蔽库仑相互作用的笛卡尔树码》，J.Compute。物理。，228, 10, 3858-3868 (2009) ·Zbl 1165.78304号
[77] Salzer，H.E.，切比雪夫点的拉格朗日插值（x_{n，nu}\equiv\cos（nu\pi/n），nu=0（1）n\）；一些不可忽视的优势，计算。J.，15，156-159（1972）·Zbl 0242.65007号
[78] Buck，I.，登顶：加速计算为世界上最快的超级计算机提供动力（2018）
[79] 津田，E。；Tsukada，M.，电子结构计算的自适应有限元方法，物理学。B版，54，7602-7605（1996）
[80] 津田，E。；Tsukada，M.，基于自适应有限元方法的大尺度电子结构计算，J.Phys。Soc.Jpn.公司。，67, 11, 3844-3858 (1998)
[81] 周，Y。；萨阿德，Y。；蒂亚戈，M.L。；Chelikowsky，J.R.，使用Chebyshev滤波子空间迭代的自洽场计算，J.Comput。物理。，219, 1, 172-184 (2006) ·Zbl 1105.65111号
[82] 周，Y。；萨阿德，Y。；蒂亚戈，M.L。；Chelikowsky，J.R.，通过Chebyshev滤波子空间加速度进行并行自我场计算，Phys。E版，74，第066704条，pp.（2006）
[83] 周，Y。；切里科夫斯基，J.R。；Saad，Y.，求解Kohn-Sham方程的无稀疏对角化的Chebyshev滤波子空间迭代法，J.Compute。物理。，274, 770-782 (2014) ·Zbl 1351.82098号
[84] Krasny，R。；Wang，L.，基于重心Hermite插值的静电粒子相互作用树码，计算。数学。生物物理学。，7, 73-84 (2019) ·兹比尔1439.78003
[85] 段，Z.-H。；Krasny，R.，计算偶极系统Ewald求和的树码算法，（2003年美国计算机学会应用计算研讨会论文集，2003年美国应用计算学会应用计算会议论文集，SAC’03（2003），计算机械协会：美国纽约州纽约市计算机械协会），172-177
[86] Boateng，H.A.，《周期库仑树方法：平行粒子网格的替代品》，Ewald，J.Chem。理论计算。，16, 1, 7-17 (2020)
[87] Vaughn，N.，GPU加速重心树码及其在Kohn-Sham密度泛函理论中的应用（2020），密歇根大学博士论文
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。