密度泛函理论（Density Functional Theory, DFT）是现代量子化学中最广泛使用的计算方法之一。它通过电子密度而非波函数来描述多电子体系，将复杂的多体问题简化为单电子在有效势场中的运动。DFT的理论基础由Hohenberg和Kohn在1964年建立的两个定理奠定，而Kohn和Sham在1965年提出的Kohn-Sham方程则为其实际应用铺平了道路。
一、Hohenberg-Kohn定理
1964年，Pierre Hohenberg和Walter Kohn发表了密度泛函理论的奠基性论文，提出了两个基本定理：
定理一：唯一性定理
对于处于外势场 (V_{ext}(mathbf{r})) 中的相互作用电子体系，其基态电子密度 ( ho(mathbf{r})) 唯一地决定了外势场（相差一个常数）。这意味着基态的所有性质都是电子密度的泛函。
这个定理的重要性在于，它将 (3N) 个坐标（N个电子的空间位置）的复杂问题简化为3个空间坐标的问题。对于一个含有100个电子的分子，波函数依赖于300个坐标，而电子密度仅依赖于3个坐标。
定理二：能量变分原理
存在一个普适的能量泛函 (E[ ho])，对于正确的基态密度 ( ho_0)，其值等于基态能量 (E_0)；对于任何其他密度 ( ho)，(E[ ho] geq E_0)。
这两个定理表明，原则上我们可以通过最小化能量泛函来找到基态能量和密度。然而，Hohenberg-Kohn定理只证明了能量泛函的存在性，并没有给出其具体形式。
二、Kohn-Sham方程
1965年，Walter Kohn和Lu Sham提出了一个巧妙的方案来解决能量泛函的具体形式问题。他们引入了一个假想的无相互作用电子体系，其电子密度与真实相互作用体系相同。
能量泛函的分解
总能量泛函被分解为以下几个部分：
[ E[ ho] = T_s[ ho] + int V_{ext}(mathbf{r}) ho(mathbf{r})dmathbf{r} + frac{1}{2}intintfrac{ ho(mathbf{r}) ho(mathbf{r}')}{|mathbf{r}-mathbf{r}'|}dmathbf{r}dmathbf{r}' + E_{xc}[ ho] ]
其中：
- (T_s[ ho]) 是无相互作用电子的动能
- (int V_{ext} ho) 是外势场（核-电子吸引）的能量
- (frac{1}{2}intintfrac{ ho ho'}{|mathbf{r}-mathbf{r}'|}) 是经典的电子-电子库仑排斥能（Hartree能）
- (E_{xc}[ ho]) 是交换-相关能，包含了所有多体效应
单电子方程
通过变分原理，可以得到一组单电子方程，即Kohn-Sham方程：
[ left[-frac{hbar^2}{2m} abla^2 + V_{eff}(mathbf{r}) ight]phi_i(mathbf{r}) = varepsilon_iphi_i(mathbf{r}) ]
其中有效势 (V_{eff}) 为：
[ V_{eff}(mathbf{r}) = V_{ext}(mathbf{r}) + intfrac{ ho(mathbf{r}')}{|mathbf{r}-mathbf{r}'|}dmathbf{r}' + V_{xc}(mathbf{r}) ]
交换-相关势 (V_{xc}(mathbf{r})) 是交换-相关能泛函对密度的泛函导数：
[ V_{xc}(mathbf{r}) = frac{delta E_{xc}[ ho]}{delta ho(mathbf{r})} ]
电子密度由Kohn-Sham轨道构造：
[ ho(mathbf{r}) = sum_{i=1}^{N}|phi_i(mathbf{r})|^2 ]
三、交换-相关泛函
Kohn-Sham DFT的精确性完全取决于交换-相关泛函 (E_{xc}[ ho]) 的近似程度。多年来，人们发展了许多不同级别的交换-相关泛函：
1. 局域密度近似（LDA）
最简单的近似是局域密度近似，其中交换-相关能密度取为均匀电子气的值：
[ E_{xc}^{LDA}[ ho] = int ho(mathbf{r})varepsilon_{xc}^{unif}( ho(mathbf{r}))dmathbf{r} ]
LDA在金属和简单半导体中表现良好，但对分子和化学键的描述往往不够准确。
2. 广义梯度近似（GGA）
GGA泛函不仅依赖于局域密度，还依赖于密度的梯度：
[ E_{xc}^{GGA}[ ho] = int ho(mathbf{r})varepsilon_{xc}^{unif}( ho(mathbf{r}))F_{xc}(s(mathbf{r}))dmathbf{r} ]
其中 (s(mathbf{r})) 是无量纲的密度梯度。常用的GGA泛函包括PBE、BLYP、PW91等。GGA显著改善了分子几何结构、解离能和反应势垒的预测。
3. 杂化泛函（Hybrid Functionals）
杂化泛函将一部分精确的Hartree-Fock交换能混合到DFT交换-相关泛函中：
[ E_{xc}^{hybrid} = aE_x^{HF} + (1-a)E_x^{DFT} + E_c^{DFT} ]
最著名的杂化泛函是B3LYP，它在分子化学中被广泛使用。杂化泛函通常能给出更准确的能量和性质预测。
4. 元广义梯度近似（meta-GGA）
meta-GGA泛函进一步引入了动能密度或Laplacian，如M06-L、TPSS等。
5. 范围分离杂化泛函
这类泛函将电子-电子相互作用分为短程和长程部分，对不同部分使用不同的处理方法，如CAM-B3LYP、ωB97X-D等，特别适用于电荷转移激发和长程相互作用。
四、DFT在化学中的应用
1. 分子几何结构优化
DFT可以准确预测分子的平衡几何结构，包括键长、键角和二面角。通过最小化体系的总能量，可以找到稳定的分子构型。
2. 振动频率和红外光谱
通过计算Hessian矩阵（能量的二阶导数），可以得到分子的振动频率，从而预测红外和拉曼光谱。
3. 反应机理和势能面
DFT可以用于搜索过渡态、计算反应路径和反应势垒，帮助理解化学反应机理。
4. 电子结构分析
通过分析Kohn-Sham轨道，可以得到分子的前线轨道（HOMO和LUMO）、电荷分布、键级等信息。
5. 激发态计算（TDDFT）
含时密度泛函理论（TDDFT）是DFT在激发态问题上的扩展，可以用于计算电子激发能、振子强度和激发态性质。
五、DFT的优势与局限
优势
- 计算效率高：与波函数方法相比，DFT的计算量随体系大小的标度更小（通常为 (O(N^3))），可以处理数百甚至数千个原子的体系
- 精度可接受：对于许多化学问题，GGA和杂化泛函可以提供与MP2相当甚至更好的精度
- 适用范围广：从分子到固体、从基态到激发态、从气相到溶液相
- 物理直观：Kohn-Sham轨道提供了类似分子轨道的物理图像
局限
- 强关联体系：对于具有强电子关联的体系（如高温超导体、某些过渡金属氧化物），标准DFT泛函往往失效
- 范德华相互作用：标准DFT泛函无法正确描述长程范德华相互作用，需要引入色散校正
- 带隙问题：DFT（特别是LDA和GGA）倾向于低估半导体的带隙
- 自相互作用误差：DFT存在自相互作用误差，导致离域化倾向和带隙低估
六、未来发展方向
尽管DFT已经取得了巨大成功，但仍有许多研究方向正在积极发展中：
1. 精确交换-相关泛函
寻找更准确的交换-相关泛函是DFT研究的核心目标。Jacob's Ladder方案提供了一个系统改进的框架，从LDA到GGA、meta-GGA、杂化泛函，再到包含更多信息的高级泛函。
2. 双杂化泛函
双杂化泛函在杂化泛函的基础上进一步加入了微扰相关能，如B2PLYP，可以显著提高能量预测的准确性。
3. 色散校正
DFT-D方法通过添加经验性的色散校正项来修正范德华相互作用，如Grimme的D3校正。
4. 机器学习与DFT
机器学习方法正在被用于构建更准确的交换-相关泛函，通过学习高精度波函数方法的结果来改进DFT。
5. 激发态DFT
开发更准确的激发态方法，如多参考态DFT、含时DFT的改进等。
七、常用DFT软件
- Gaussian：商业软件，广泛用于分子量子化学计算
- VASP：商业软件，专注于周期性体系的DFT计算
- Quantum ESPRESSO：开源软件，用于固体和表面的第一性原理计算
- ORCA：功能强大的量子化学程序，支持各种DFT方法
- CP2K：开源软件，支持DFT和分子动力学
- Psi4：开源量子化学软件包
- PySCF：基于Python的量子化学程序
八、结论
密度泛函理论自1964年诞生以来，已经成为量子化学和材料科学中最重要的计算工具。Hohenberg-Kohn定理奠定了理论基础，Kohn-Sham方程提供了实用的计算方法，而交换-相关泛函的不断发展则推动着DFT精度的持续提高。
尽管DFT在某些方面存在局限，但其独特的平衡——在计算成本和预测精度之间——使其在化学、物理、材料科学等领域发挥着不可替代的作用。随着计算能力的提升和理论方法的进步，DFT必将在未来的科学研究中继续发挥关键作用。
Walter Kohn因其对密度泛函理论的贡献获得了1998年诺贝尔化学奖，这充分说明了DFT在现代化学中的重要地位。从分子设计到材料发现，从药物研发到催化机理，DFT正在帮助我们理解和预测物质世界的奥秘。
发表评论 取消回复