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

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

摘要：本文综述了密度泛函理论发展的基础及其最新进展，介绍了求解具体物理化学问题时用到的几种常用的数值计算方法，另外对密度泛函理论的发展进行了展望。密度泛函理论的发展以寻找合适的交换相关近似为主线，从最初的局域密度近似、广义梯度近似到现在的非局域泛函、自相互作用修正，多种泛函形式的相继出现使得密度泛函理论可以提供越来越精确的计算结果。另外，在密度泛函理论体系发展的同时，相应的数值计算方法的发展也非常迅速。随着密度泛函理论本身及其数值方法的发展，它的应用也越来越广泛，一些新的应用领域和研究方向不断涌现。

关键词：密度泛函数值计算发展应用

1 研究背景

量子力学作为20世纪最伟大的发现之一，是整个现代物理学的基石。量子力学最流行的表述形式是薛定谔的波动力学形式，核心是波函数及其运动方程薛定谔方程。对一个外势场v(r)中的N电子体系，量子力学的波动力学范式可以表示成：

即对给定的外势，将其代入薛定谔方程可以得到电子波函数，可以得到所有可观测量的值。

当用量子力学处理真实的物理化学体系时，传统的波动力学方法便显得有点力不从心。因为在大多数情况下，人们只是关心与实验相关的一部分信息，如能量、密度等。所以，人们希望使用一些较简单的物理量来构造新的理论[1]。

电子密度泛函理论是上个世纪60年代在Thomas-Fermi理论的基础上发展起来的量子理论的一种表述方式。传统的量子理论将波函数作为体系的基本物理量，而密度泛函理论则通过粒子密度来描述体系基态的物理性质。因为粒子密度只是空间坐标的函数，这使得密度泛函理论将3N 维波函数问题简化为3维粒子密度问题,十分简单直观。另外，粒子密度通常是可以通过实验直接观测的物理量。粒子密度的这些优良特性，使得密度泛函理论具有诱人的应用前景。

2 密度泛函理论的基础

Thomas-Fermi模型

1927 年Thomas和Fermi分别提出：体系的动能可以通过体系的电子密度表达出来。他们提出了一种的均匀电子气模型，把空间分割成足够小的立方体，通过在这些立方体中求

解无限势阱中粒子的Schrodinger 方程(假设电子之间无相互作用)，得到相应的能量和密度的表达式。但是这个早期的理论有点粗糙，更严重的是用在化学方面，甚至出现分子不会形成键的计算结果，这个结论严重阻碍了Thomas-Fermi理论的发展。

Hohenberg-Kohn定理

1964年Hohenberg与Kohn发表了后来极著名的Hohenberg-Kohn(HK)定理，标志着现代密度泛函理论的产生[2]。Hohenberg-Kohn定理的表述如下：

第一定理：对于一个共同的外部势v(r)，相互作用的多粒子系统的所有基态性质都由(非简并)基态的电子密度分布ρ唯一地决定。

第二定理：如果ρ是体系正确的密度分布，则E[ρ]是最低的能量，即体系的基态能量。

Kohn-Sham方程

有了上述两个定理，剩下的问题就是能量泛函的具体表述形式。Kohn和他的博士后沈吕九引进了一个与相互作用多电子体系有相同电子密度的假想的非相互作用多电子体系。因为电子密度一般可以表示成轨道形式（N-表示问题），这个假想的非相互作用体系的动能算符期望值可以非常简单的写成各电子动能的和。

其中φ(r)是密度函数对应的Kohn-Sham（KS）轨道，将能量泛函对KS轨道进行变分可以得到著名的KS方程。

3 密度泛函理论的进展

密度泛函理论体系及其数值实现方法经过几十年的发展已被广泛地应用在化学、物理、材料和生物等学科中，Kohn也因为他对密度泛函理论的贡献获得1998年的诺贝尔化学奖[3]。

在密度泛函理论中，所有的近似都被集中到称为交换相关能的一项上，所以密度泛函理论的精度直接由交换相关能量泛函的近似形式决定。寻找更好的交换相关近似就成为密度泛函理论体系发展的一条主线。另一方面，在现有基础上对密度泛函理论的扩展也是一个活跃的领域。

交换相关能量泛函

交换相关能量泛函的一个最初的简单近似是局域密度近似(LDA)，即用具有相同密度的

均匀电子气的交换相关泛函作为对应的非均匀系统的近似值。如果进一步分别考虑不同自旋分量的电子密度，可以得到自旋极化的局域密度近似L(S)DA。尽管L(S)DA获得了巨大的成功，但是也有许多不足之处,比如系统地高估结合能。在L(S)DA基础上的改进有广义梯度近似(GGA)。在GGA近似下，交换相关能是电子密度及其梯度的泛函。构造GGA交换相关泛函的方法分为两个流派。一个是以Becke为首的一派，认为“一切都是合法的”，人们可以任何理由选择任何可能的泛函形式，而这种形式的好坏由实际计算来决定。通常这类泛函的参数由拟合大量的计算数据得到。另外一个流派以Perdew为首，他们认为发展交换相关泛函必须以一定的物理规律为基础,这些规律包括标度关系、渐进行为等。基于这种理念构造的一个著名的GGA 泛函是PBE泛函[4]，也是现在用得最广泛的GGA泛函之一。

LDA与GGA近似的密度泛函理论计算构成了当今最流行的电子结构计算方案，但是这种方案仍有待改进的地方。下面将介绍几种重要的提高精度的交换相关泛函方向。

精确交换和屏蔽交换

精确交换是杂化密度泛函理论中的一个概念，杂化密度泛函正是通过包含精确交换来提高计算精度的。屏蔽交换的LDA(sX-LDA)方法[5]在非局域框架内描述交换相关穴，它通过包含一个屏蔽交换项来提高交换相关近似的质量。扩展系统中的电子屏蔽可以用一个Thomas-Fermi指数衰减因子来考虑。Engel等[6]通过研究均匀电子气的线性响应函数考察了这种方法的理论基础。对一些半导体的自洽计算表明sX-LDA可以相对LDA显著提高计算结果,给出更接近实验值的能隙和光学性质。

自相互作用修正

自相互作用修正(SIC) 最早由Stoll等人[7]提出，用来消除LSDA带结构理论中的不真实的电子与自身的相互作用。现在已经有好几种基于现代密度泛函理论的消除自相互作用的方案，一个典型的代表是Perdew等[8]提出的泛函形式。通过自相互作用修正可以显著改善计算结果，但是SIC中泛函是与轨道有关的，这使得我们需要重新得到自洽方程。

其他交换相关泛函

为了设计一种能处理包含长程关联相互作用体系的交换相关能量泛函，只包括密度和密度梯度信息也许是不够的。这种“不够”是不能够恰当地使用这些信息，所以必须求助于其他信息。这种比GGA 近似包含更多半局域信息的泛函就叫做meta-GGA泛函。这些信息可以是密度的更高阶梯度、KS轨道梯度或者其他一些系统特征变量。比如，PKZB泛函[9]就在PBE泛函的基础上包括了占据轨道的动能密度信息。而最近的TPSS泛函[10]在PKZB泛函的基础上，首次提出了完全不依赖经验参数拟合的meta-GGA泛函，爬到Jacob楼梯[11]的第三阶。