电子结构计算的有限元方法.pdfVIP

中国科学: 化学 2015 年 第45 卷 第8 期: 800 ~ 811 《中国科学》杂志社 SCIENTIA SINICA Chimica SCIENCE CHINA PRESS 评 述 中国科学院学部 科学与技术前沿论坛 理论与计算化学发展战略专题 电子结构计算的有限元方法 戴小英, 周爱辉* 科学与工程计算国家重点实验室; 中国科学院数学与系统科学研究院计算数学与科学工程计算研究所, 北京 100190 *通讯作者, E-mail: azhou@ 收稿日期: 2015-03-26; 接受日期: 2015-06-23; 网络版发表日期: 2015-07-27 doi: 10.1360/N032015-00047 摘要 本文从数学角度认识电子结构计算的困难与挑战, 评述有关计算方法的优缺点以 关键词 及在超级计算机上实现前景, 并扼要地介绍第一原理电子结构计算的有限元方法和理论. 电子结构 有限元方法 自适应 特征值问题 并行计算 1 引言 本文我们将从数学角度认识基于 Kohn-Sham 方 程的电子结构计算的困难与挑战, 评述有关计算方 电子结构可通过包括Hartree-Fock 方程和Kohn- 法的优缺点以及在超级计算机上实现前景. 我们还 Sham 方程在内的非线性特征值模型(或其等价形式) 将特别介绍第一原理电子结构计算的( 自适应)有限 来刻画. 求解这些非线性特征值问题既可通过反复 元方法和理论. 基于能量直接极小化的电子结构计 求解线性化后的线性特征值问题来实现, 也可转化 算研究见[1~3]及所引文献. 有关电子结构模型与计 为最优化问题计算. 电子结构计算方法研究涉及离 算的数学问题见[4]. 电子结构计算的有限元方法实 散方法、非线性迭代法、代数特征值求解器以及相应 现细节可见[5~10]. 的程序实现技术等方面. 现有的非线性特征值问题 (或其等价形式) 的离散格式大体上可分为三类: 倒空 2 Kohn-Sham 方程 间方法、局部基集法以及实空间方法. 这些方法各有 其优势和不足, 其中倒空间方法和局部基集法研究 描述非相对论多电子体系电子结构的基本的数 相对成熟, 而实空间方法尚处于发展阶段. 实空间方 学模型是Schrödinger 方程[11]: 法主要包括有限差分方法、小波方法、有限体方法与 (T V V ) E (in R