密度泛函理论的进展与问题.docxVIP

密度泛函理论的进展与问题摘要：本文综述了密度泛函理论发展的基础及其最新进展，介绍了求解具体物理化学问题时用到的几种常用的数值计算方法，另外对密度泛函理论的发展进行了展望。密度泛函理论的发展以寻找合适的交换相关近似为主线，从最初的局域密度近似、广义梯度近似到现在的非局域泛函、自相互作用修正，多种泛函形式的相继出现使得密度泛函理论可以提供越来越精确的计算结果。另外，在密度泛函理论体系发展的同时，相应的数值计算方法的发展也非常迅速。随着密度泛函理论本身及其数值方法的发展，它的应用也越来越广泛，一些新的应用领域和研究方向不断涌现。关键词：密度泛函数值计算发展应用1 研究背景量子力学作为20世纪最伟大的发现之一，是整个现代物理学的基石。量子力学最流行的表述形式是薛定谔的波动力学形式，核心是波函数及其运动方程薛定谔方程。对一个外势场v(r)中的N电子体系，量子力学的波动力学范式可以表示成：即对给定的外势，将其代入薛定谔方程可以得到电子波函数，可以得到所有可观测量的值。当用量子力学处理真实的物理化学体系时，传统的波动力学方法便显得有点力不从心。因为在大多数情况下，人们只是关心与实验相关的一部分信息，如能量、密度等。所以，人们希望使用一些较简单的物理量来构造新的理论[1]。电子密度泛函理论是上个世纪60年代在Thomas-Fermi理论的基础上发展起来的量子理论的一种表述方式。传统的量子理论将波函数作为体系的基本物理量，而密度泛函理论则通过粒子密度来描述体系基态的物理性质。因为粒子密度只是空间坐标的函数，这使得密度泛函理论将3N 维波函数问题简化为3维粒子密度问题,十分简单直观。另外，粒子密度通常是可以通过实验直接观测的物理量。粒子密度的这些优良特性，使得密度泛函理论具有诱人的应用前景。2 密度泛函理论的基础2.1 Thomas-Fermi模型1927 年 Thomas和Fermi分别提出：体系的动能可以通过体系的电子密度表达出来。他们提出了一种的均匀电子气模型，把空间分割成足够小的立方体，通过在这些立方体中求解无限势阱中粒子的 Schrodinger 方程(假设电子之间无相互作用)，得到相应的能量和密度的表达式。但是这个早期的理论有点粗糙，更严重的是用在化学方面，甚至出现分子不会形成键的计算结果，这个结论严重阻碍了Thomas-Fermi理论的发展。2.2 Hohenberg-Kohn定理1964年Hohenberg与Kohn发表了后来极著名的 Hohenberg-Kohn(HK)定理，标志着现代密度泛函理论的产生[2]。Hohenberg-Kohn定理的表述如下：第一定理：对于一个共同的外部势v(r)，相互作用的多粒子系统的所有基态性质都由(非简并)基态的电子密度分布ρ唯一地决定。第二定理：如果ρ是体系正确的密度分布，则E[ρ]是最低的能量，即体系的基态能量。2.3 Kohn-Sham方程有了上述两个定理，剩下的问题就是能量泛函的具体表述形式。Kohn和他的博士后沈吕九引进了一个与相互作用多电子体系有相同电子密度的假想的非相互作用多电子体系。因为电子密度一般可以表示成轨道形式（N-表示问题），这个假想的非相互作用体系的动能算符期望值可以非常简单的写成各电子动能的和。其中φ(r)是密度函数对应的Kohn-Sham（KS）轨道，将能量泛函对KS轨道进行变分可以得到著名的KS方程。3 密度泛函理论的进展密度泛函理论体系及其数值实现方法经过几十年的发展已被广泛地应用在化学、物理、材料和生物等学科中，Kohn也因为他对密度泛函理论的贡献获得1998年的诺贝尔化学奖[3]。在密度泛函理论中，所有的近似都被集中到称为交换相关能的一项上，所以密度泛函理论的精度直接由交换相关能量泛函的近似形式决定。寻找更好的交换相关近似就成为密度泛函理论体系发展的一条主线。另一方面，在现有基础上对密度泛函理论的扩展也是一个活跃的领域。3.1 交换相关能量泛函交换相关能量泛函的一个最初的简单近似是局域密度近似(LDA)，即用具有相同密度的均匀电子气的交换相关泛函作为对应的非均匀系统的近似值。如果进一步分别考虑不同自旋分量的电子密度，可以得到自旋极化的局域密度近似L(S)DA。尽管L(S)DA获得了巨大的成功，但是也有许多不足之处,比如系统地高估结合能。在L(S)DA基础上的改进有广义梯度近似(GGA)。在GGA近似下，交换相关能是电子密度及其梯度的泛函。构造GGA交换相关泛函的方法分为两个流派。一个是以Becke为首的一派，认为“一切都是合法的”，人们可以任何理由选择任何可能的泛函形式，而这种形式的好坏由实际计算来决定。通常这类泛函的参数由拟合大量的计算数据得到。另外一个流派以Perdew为首，他们认为发展交换相关泛函必须以一定的物理规律为基础,这些规律包括标度关系、渐进行