×
周云凯;詹姆斯·切利科夫斯基（James R.Chelikowsky）。;优素福·萨阿德

求解Kohn-Sham方程的无稀疏对角化的Chebyshev滤波子空间迭代法。 （英语） Zbl 1351.82098号

J.计算。物理学。 274, 770-782 (2014).
总结：第一原则密度泛函理论电子结构问题的DFT计算需要求解Kohn-Sham方程，这需要求解非线性特征值问题。求解特征值问题通常是DFT计算中最昂贵的部分。对于大规模问题，计算显式特征向量的稀疏迭代对角化方法可能很快变得难以使用。Chebyshev滤波子空间迭代（CheFSI）方法避免了大多数特征向量的显式计算，并在DFT自洽场（SCF）计算中大大提高了迭代对角化方法的速度。然而，CheFSI方法的原始公式在第一个SCF步骤中使用稀疏迭代对角化，为后面的SCF步骤的子空间滤波提供初始向量。这种对角化对于大规模问题来说代价高昂。我们开发了一个新的初始滤波步骤，以完全避免这种对角化，从而使CheFSI方法在所有SCF步骤都没有稀疏迭代对角化。我们的新方法节省了内存使用，并且比原来的CheFSI方法快两到三倍。

引用于20文件

MSC公司：

82C80码 时间相关统计力学的数值方法（MSC2010）

关键词：

密度泛函理论;电子结构问题;实空间赝势;自一致性;切比雪夫滤波器;哈密顿量;对角化;非线性特征值问题;子空间滤波

软件：

ARPACK公司;JDQZ公司;TRLan公司;KSSOLV公司;巴黎证券交易所;洛佩克。米;CheFSI公司
PDF格式BibTeX公司 XML格式引用
\textit{Y.Zhou}等人，J.Compute。物理学。274770--782（2014年；Zbl 1351.82098）
全文： 内政部

参考文献：

[1] Hohenberg，P。；Kohn，W.，《非均匀电子气体》，Phys。版本136，B864-871（1964）
[2] 科恩，W。；Sham，L.J.，《包含交换和相关效应的自洽方程》，《物理学》。修订版，140，A1133-1138（1965）
[3] 菲利普斯，J.C。；Kleinman，L.，计算晶体和分子波函数的新方法，物理学。修订版，116287-294（1959年）·Zbl 0095.44601号
[4] 切里科夫斯基，J.R。；科恩，M.L。，从头算半导体赝势，（半导体手册，第1卷（1992），爱思唯尔：爱思唯尔阿姆斯特丹），59
[5] Martin，R.M.，《电子结构：基本理论和实用方法》（2004），剑桥大学出版社·Zbl 1152.74303号
[6] 佩恩，M.C。；特特，M.P。；阿伦·D.C。；阿里亚斯，T.A。；Joannopoulos，J.D.，从头算总能量计算的迭代最小化技术：分子动力学和共轭梯度，修订版。物理。，64, 1045-1097 (1992)
[7] Kresse，G。；Furthmüller，J.，使用平面波基集进行从头算全能量计算的高效迭代方案，Phys。B版，54，11169-11186（1996）
[8] 切里科夫斯基，J.R。；Troullier，N。；Saad，Y.，《有限差分伪势法：无基电子结构计算》，Phys。修订稿。，72, 1240-1243 (1994)
[9] 切里科夫斯基，J.R。；Troullier，N。；Wu，K。；Saad，Y.，《高阶有限差分赝势法：双原子分子的应用》，Phys。B版，5011355-11364（1994）
[10] Seitsonen，A.P。；Puska，M.J。；Nieminen，R.M.，《实空间电子结构计算：有限差分和共轭梯度方法的组合》，《物理学》。B版，51，14057-14061（1995）
[11] Beck，T.L.，密度泛函理论中的实空间网格技术，修订版。物理。，72, 1041-1080 (2000)
[12] Kronik，L。；Makmal，A。；蒂亚戈，M.L。；Alemany，M.M.G。；Jain，M。；黄，X。；萨阿德，Y。；Chelikowsky，J.R.，PARSEC——真实空间电子结构计算的伪势算法：纳米结构的最新进展和新应用，Phys。Solidi B状态，2431063（2006）
[13] 萨阿德，Y。；Stathopoulos，A.公司。；切里科夫斯基，J.R。；Wu，K。；Öğüt，S.，电子结构计算中大特征值问题的解决方案，BIT，36563-578（1996）·Zbl 0862.65059号
[14] Lehoucq，R.B。；索伦森特区。；Yang，C.，ARPACK用户指南：用隐式重启Arnoldi方法解决大规模特征值问题（1998），SIAM·Zbl 0901.65021号
[15] Wu，K。；坎宁，A。；西蒙，H.D。；Wang，L.-W.，电子结构计算的厚重启动Lanczos方法，J.Compute。物理。，154, 156-173 (1999) ·Zbl 0963.82048号
[16] Knyazev，A.V.，朝向最优预处理特征解算器：局部最优块预处理共轭梯度法，SIAM J.Sci。计算。，23, 517-541 (2001) ·Zbl 0992.65028号
[17] 周，Y。；Saad，Y.，大型对称特征值问题的Chebyshev-Davidson算法，SIAM J.矩阵分析。申请。，2954-971年9月29日（2007年）·Zbl 1151.65321号
[18] 周，Y。；萨阿德，Y。；蒂亚戈，M.L。；Chelikowsky，J.R.，使用Chebyshev滤波子空间迭代的自洽场计算，J.Compute。物理。，219, 172-184 (2006) ·Zbl 1105.65111号
[19] 周，Y。；萨阿德，Y。；蒂亚戈，M.L。；Chelikowsky，J.R.，使用Chebyshev滤波子空间加速度的并行自我场计算，Phys。E版，74066704（2006）
[20] 蒂亚戈，M.L。；周，Y。；Alemany，M。；萨阿德，Y。；Chelikowsky，J.R.，《铁从原子到体的磁性演化》，《物理学》。修订稿。，97, 147201 (2006)
[21] Zhou，Y.，一种用于大型特征值问题的带内-外重启的块Chebyshev-Davidson方法，J.Compute。物理。，229, 9188-9200 (2010) ·Zbl 1203.65077号
[22] García-Cervera，C.J。；Lu，J。；Xuan，Y。；Weinan，E.，Kohn-Sham密度泛函理论的线性尺度子空间迭代算法，Phys。B版，79，1-13（2009年）
[23] 渭南，E。；李·T。；Lu，J.，特征子空间的局部化基和算子压缩，Proc。国家。阿卡德。科学。美国，107，4，1273-1278（2010）·Zbl 1205.15017号
[24] Goedecker，S。；Teter，M.P.，紧束缚电子结构计算和局域轨道紧束缚分子动力学，Phys。B版，51，9455-9464（1995）
[25] Ozaki，T.，《利用局部化基函数进行大规模电子结构计算的高效低阶标度方法》，Phys。B版，82，075131（2010）
[26] 林，L。；Lu，J。；Ying，L。；Weinan，E.，Kohn Sham密度泛函理论的优化局部基集，计算机J。物理。，2314515-4529（2012年）·Zbl 1250.82002号
[27] 林，L。；Lu，J。；Ying，L。；Weinan，E.，《非连续Galerkin框架下Kohn-Sham密度泛函理论的自适应局部基集I：总能量计算》，J.Compute。物理。，231, 2140-2154 (2012) ·Zbl 1251.82008年
[28] Ozaki，T.，用于大规模电子结构计算的Fermi-Dirac函数的连续分数表示，Phys。B版，75，035123（2007）
[29] Sidje，R.B。；Saad，Y.，费米-迪拉克函数的有理逼近及其在密度泛函理论中的应用，Numer。算法，56455-479（2011）·Zbl 1211.65026号
[30] 林，L。；Lu，J。；Ying，L。；Weinan，E.，Fermi-Dirac函数的基于极点的近似，Chin。数学安。，序列号。B、 30229-742（2009年）·Zbl 1188.41007号
[31] 佩杜，J.P。；伯克，K。；Ernzerhof，M.，《简化广义梯度近似》，Phys。修订稿。，77, 3865-3868 (1996)
[32] Kümmel，S。；Perdew，J.P.，《优化有效势简化：轨道泛函、轨道位移和精确的Kohn-Sham交换势》，Phys。B版，68，035103（2003）
[33] 阿尼西莫夫，V.I。；Aryasetiawan，F。；Lichtenstein，A.I.，强关联体系电子结构和光谱的第一原理计算：LDA+U方法，J.Phys。康登斯。Matter，9767-808（1997）
[34] 科科奇奥尼，M。；de Gironcoli，S.，LDA+U方法中有效相互作用参数计算的线性响应方法，Phys。B版，71，035105（2005）
[35] Kulik，H.J。；科科奇奥尼，M。；谢利斯，D.A。；Marzari，N.，《过渡金属化学中的密度泛函理论：自洽Hubbard方法》，Phys。修订稿。，97, 103001 (2006)
[36] 科赫，W。；Holthausen，M.C.，《密度泛函理论化学家指南》（2001），Wiley-VCH
[37] 恩格尔，B。；Dreizler，R.M.，《密度泛函理论：高级课程》（理论与数学物理（2011），Springer）·Zbl 1216.81004号
[38] Axelsson，O.，《迭代求解方法》（1994），剑桥大学出版社·Zbl 0795.65014号
[39] Rivlin，T.J.，《函数逼近导论》（2003），多佛·Zbl 0189.06601号
[40] 周，Y。；李，R.-C.，大厄米矩阵谱的边界，线性代数应用。，435480-493（2011年）·兹比尔1221.15022
[41] Saad，Y.，Chebyshev解决非对称特征值问题的加速技术，数学。计算。，42, 567-588 (1984) ·Zbl 0539.65013号
[42] Zhou，Y.，计算HITS ExpertRank向量的实际加速度，J.Compute。申请。数学。，236, 4398-4409 (2012) ·Zbl 1287.65027号
[43] 霍恩，R.A。；Johnson，C.R.，矩阵分析（1985），剑桥大学出版社·Zbl 0576.15001号
[44] 杨，C。；Meza，J.C。；Lee，B。；Wang，L.-W.，KSSOLV-一个用于求解Kohn-Sham方程的MATLAB工具箱，ACM Trans。数学。软质。，36, 1-35 (2009) ·Zbl 1364.65112号
[45] Ceperley，D.M。；Alder，B.J.，用随机方法研究电子气体的基态，物理学。修订稿。，45, 566-569 (1980)
[46] 方，H.-R。；Saad，Y.，非线性加速度的两类多斜面方法，Numer。线性代数应用。，16, 3, 197-221 (2009) ·Zbl 1224.65134号
[47] 斯蒂克斯鲁德，L。；Karki，B.B.，地球下地幔中（MgSiO_3）液体的结构和冻结，科学，310，297-299（2005）
[48] Umemoto，K。；温茨科维奇，R.M。；Allen，P.B.，气态巨行星和类地系外行星核心中（MgSiO_3）的解离，《科学》，311，983-986（2006）
[49] Wu，K。；Simon，H.，《大型对称特征值问题的厚重启动Lanczos方法》，SIAM J.矩阵分析。申请。，22, 602-616 (2000) ·Zbl 0969.65030号
[50] Sleijpen，G.L.G。；van der Vorst，H.A.，线性特征值问题的Jacobi-Davidson迭代法，SIAM J.矩阵分析。申请。，17, 401-425 (1996) ·Zbl 0860.65023号
[51] Stathopoulos，A。；萨阿德，Y。；Wu，K.，戴维森的动态厚重启和隐式重启的阿诺迪方法，SIAM J.Sci。计算。，19, 227-245 (1998) ·Zbl 0924.65028号
[52] Stathopoulos，A。；McCombs，J.R.，有限记忆下厄米特特征问题的近似最优预处理方法。第二部分：寻找许多特征值，SIAM J.Sci。计算。，29, 2162-2188 (2007) ·Zbl 1151.65320号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。