密度泛函理论 (DFT)

玻尔百科

核心要点

密度泛函理论 (DFT) 提出，量子系统的所有基态性质都可以由其电子密度确定，这是对多电子波函数的一种更简单的替代方案。
Kohn-Sham 方法为 DFT 提供了一条实用途径，它使用一个虚构的无相互作用电子体系来计算动能并生成正确的电子密度。
任何 DFT 计算的准确性都取决于近似的交换相关泛函，该泛函包含了所有复杂的多体量子效应。
DFT 及其扩展是功能强大的工具，用于预测各种性质，从分子颜色和反应动力学到电池电压和晶体的振动模式。

引言

在量子力学领域，描述一个包含许多相互作用电子的系统是一项艰巨的挑战，即所谓的“维度灾难”，其中波函数的复杂性随粒子数的增加而呈指数级增长。这使得对大多数原子、分子和材料进行直接求解在实践中变得不可能。密度泛函理论 (DFT) 提供了一种革命性的替代方案。它通过关注一个简单得多的量——电子密度，回避了复杂的多电子波函数。这种开创性的方法假定，基态密度包含确定系统性质所需的所有信息。本文将探讨这一优雅思想是如何付诸实践的。第一部分原理与机制将揭示 Hohenberg-Kohn 定理的理论保证，以及使 DFT 成为实用计算工具的巧妙的 Kohn-Sham 形式体系。随后，应用与跨学科联系部分将展示 DFT 的巨大威力，说明它如何被用于设计新分子、预测先进材料的性质，甚至与其他科学学科建立联系。

原理与机制

假设你面临一项看似不可能的任务：同时描述一个房间里所有舞者精确而复杂的舞蹈。你可以尝试写下每一位舞者的确切路径——他们的每一步、每一次转身和每一次跳跃。信息量将是惊人的、无法承受的。对于少数几个舞者，这或许还可管理。但对于一千个，这就毫无希望了。这正是量子化学家和物理学家在试图求解多电子波函数 $\Psi$ 时所面临的挑战。这个波函数是量子力学中的至高对象，它是系统中*每一个电子*坐标的函数，即 $\Psi(\mathbf{r}_1, \mathbf{r}_2, \dots, \mathbf{r}_N)$ 。这个对象的复杂性随电子数 $N$ 呈指数级增长，这是一个被称为“维度灾难”的灾难性问题。

现在，如果你发现了一条神奇的捷径呢？如果，你不需要追踪每一位舞者，而只需要知道房间里每一点的舞者*密度*呢？想象一下，你知道舞池中央每平方米有五位舞者，而靠近墙壁的地方只有一位。你已经将关于 $3N$ 个变量的海量信息压缩成一个只包含三个空间变量的简单函数——密度 $\rho(\mathbf{r})$ 。这就是密度泛函理论（DFT）的概念核心。它提出，对于一个量子系统的基态，我们想要的所有信息都编码在这个简单得多的电子密度中，而不是在那个极其复杂的波函数中。

本章将讲述这个好到几乎不真实的想法是如何实现的。我们将从保证其有效性的大胆定理出发，一路探索到使其成为当今应用最广泛的电子结构计算工具的巧妙而实用的机制。

理论保证：Hohenberg-Kohn 定理

你可能会理所当然地感到怀疑。一个简单的密度如何能包含完整波函数的所有信息？其证明是整个物理科学中最优雅、最深刻的证明之一，它被阐述在 Hohenberg-Kohn (HK) 定理中。

第一条定理是唯一性证明。它确立了系统的基态电子密度 $\rho(\mathbf{r})$ 唯一地决定了外势 $v_{\text{ext}}(\mathbf{r})$ 。可以这样想：密度就像一块化石。古生物学家看到一个化石足迹，不仅可以推断出留下足迹的脚的形状，还能推断出该生物的体重、步态及其生活环境的细节。同样地，电子密度是外势（原子核）留下的足迹。由于势和电子数完全定义了系统的哈密顿算符，因此密度必然隐含地决定了关于基态的一切，包括总能量和完整的、复杂的波函数本身。

第二条定理提供了实用的工具：关于密度的变分原理。它指出，存在一个普适的密度能量泛函 $E[\rho]$ ，而真实的基态密度是使该泛函最小化的那一个。这意味着我们可以通过尝试不同的候选密度并找出能量最低的那一个，来寻找正确的密度。

这两条定理共同提供了一个革命性的保证。它们证明了基于电子密度的理论不仅仅是一种近似，而是在原理上对基态是精确的。这与像经典的 Hartree-Fock (HF) 理论等其他方法形成鲜明对比。HF 方法的核心假设——多电子波函数可以由单个斯莱特行列式描述——是一种内在的近似。它从根本上忽略了电子为躲避彼此而进行的复杂、相关的运动。它所忽略的能量被称为相关能。而 HK 定理则承诺，必定存在一个精确的能量泛函，它能完美地解释所有这些效应。寻找这个泛函是现代 DFT 的核心追求。

巧妙的技巧：Kohn 与 Sham 的辅助体系

HK 定理是一个美好的承诺，但它们没有告诉我们如何构建至关重要的能量泛函。其中动能部分尤其棘手。如何仅用密度来表示相互作用电子的动能呢？

这正是 Walter Kohn 和 Lu Jeu Sham 做出他们获得诺贝尔奖的创举的地方。他们提出的方法只能被描述为一次巧妙的概念“骗局”。他们说：我们不要直接尝试求解真实的、混乱的相互作用电子体系。相反，让我们虚构一个辅助的无相互作用电子体系，并巧妙地操控它，使其具有与我们真实的相互作用体系完全相同的基态密度。

为什么这如此巧妙？因为我们知道如何精确求解一个无相互作用电子的体系！它们的波函数是一个简单的斯莱特行列式，动能也很容易计算。这个虚构的体系由一组简单的单粒子方程控制，这些方程现在以Kohn-Sham 方程而闻名。

这一策略引入了 Kohn-Sham 轨道。理解这些轨道是什么——以及它们不是什么——是至关重要的。在 Hartree-Fock 理论中，斯莱特行列式及其轨道是真实体系波函数的一个直接（尽管是近似的）表示。而在 Kohn-Sham DFT 中，行列式及其轨道没有这样的直接物理意义。它们纯粹是数学构造，是虚构世界中的工具，其唯一目的是生成正确的密度并为我们提供无相互作用体系的动能 $T_s$ 。它们是用来搭建最终大厦的脚手架，一旦施工完成就会被拆除。

我们所有无知的总和：交换相关泛函

当然，物理学中没有免费的午餐。Kohn-Sham 方法的精妙之处在于，它将所有困难、混乱的电子相互作用物理学都归入一个单一的、神秘的项中：交换相关（XC）能量泛函 $E_{xc}[\rho]$ 。在 Kohn-Sham 形式体系中，总能量写为：

$E[\rho] = T_s[\rho] + \int v_{\text{ext}}(\mathbf{r}) \rho(\mathbf{r}) d\mathbf{r} + J[\rho] + E_{xc}[\rho]$

在这里， $T_s[\rho]$ 是我们虚构的无相互作用电子的动能，积分项是来自原子核的经典势能，而 $J[\rho]$ 是 Hartree 能量——电子云与自身的经典静电排斥能。最后一项 $E_{xc}[\rho]$ 是个神奇的黑箱。它必须包含所有其他东西：

交换能，它源于泡利不相容原理，该原理使自旋相同的电子彼此分开。
相关能，它描述了任何自旋的电子如何因电荷排斥而动态地相互躲避。这正是 Hartree-Fock 理论完全忽略的项。
真实的相互作用体系与我们虚构的无相互作用体系之间的动能差。

为了体会这个项的作用，考虑最简单的系统：单个氢原子。实际上，孤立的电子只感受到原子核的吸引力，没有电子-电子相互作用。然而，在我们的 DFT 公式中，Hartree 项 $J[\rho]$ 计算了电子自身的电荷云与自身的静电排斥——这是一种毫无意义的自相互作用。为了让 DFT 给出氢原子的正确答案（它必须如此！），对于这个单电子体系，精确的交换相关泛函 $E_{xc}[\rho]$ 必须恰好等于 Hartree 能量的负值，即 $E_{xc}[\rho] = -J[\rho]$ 。它的任务是完美地抵消虚假的自相互作用。这提供了一个深刻的见解：XC 泛函的一个主要作用是校正由经典 Hartree 项引入的非物理的自排斥。

近似的艺术：泛函“动物园”

如果我们知道 $E_{xc}[\rho]$ 的确切形式，我们就能精确地（在基态范围内）求解任何原子、分子或材料的电子结构。但我们不知道。确切的泛函是未知的，而且可能极其复杂。因此，我们必须对其进行近似。为 $E_{xc}$ 开发更好的近似是现代 DFT 研究的主要战场，并已催生出数百种泛函的“动物园”。

这些近似通常被组织在一个概念性的“雅各布天梯”上，每一级阶梯都增加一个新的要素以增强物理真实性，但这通常伴随着更高的计算成本。

局域密度近似 (LDA): 最简单的一级阶梯假设空间中任意一点的交换相关能与具有相同密度的均匀电子气中的交换相关能相同。它出人意料地有效，但常常在细节上出错。
广义梯度近似 (GGA): 更高一级阶梯不仅考虑某一点的密度，还考虑了它的变化速度（其梯度 $\nabla\rho$ ）。这使得泛函能够更好地区分不同的化学环境。
杂化泛函: 一个重大的突破来自于人们意识到可以混合不同成分。我们知道 Hartree-Fock 理论缺乏相关能，但它以一种特定的、精确的方式处理交换能（对于单个行列式而言）。杂化泛函将一部分这种“精确”的 HF 交换能与一个 GGA 泛函混合。这往往能纠正纯 GGA 泛函的一些最严重的失败，特别是自相互作用误差。

标准 DFT 近似（如 LDA 和 GGA）的一个关键计算优势在于，其对应的势 $v_{xc}(\mathbf{r}) = \frac{\delta E_{xc}[\rho]}{\delta \rho(\mathbf{r})}$ 是一个简单的局域势。它是一个乘法函数——在 $\mathbf{r}$ 点的势仅取决于 $\mathbf{r}$ 点或其附近的密度性质。这比 Hartree-Fock 理论中的交换算符要简单得多，后者是一个非局域算符。非局域算符在计算上是一个噩梦：它在某一点的作用取决于波函数在空间中其他所有地方的值。这种局域性是 DFT 高效率的一个主要原因。杂化泛函通过包含一部分非局域的 HF 交换能，牺牲了部分简洁性以换取更高的准确性。

由于这些泛函的数学形式非常复杂，它们对总能量的贡献无法对分子进行解析计算（用纸和笔）。取而代之的是，计算机程序在空间中的数值积分网格点上计算这些项，将每个点的贡献相加以获得总的 XC 能量。这是优雅理论与计算实践需求相结合的典型例子。

那么，DFT 是一种*从头算*（ab initio，“第一性原理”）方法，还是仅仅是复杂的曲线拟合？学术界普遍认为它是第一性原理方法。为什么？因为这些泛函虽然是近似的，但它们被设计成普适的。像 PBE（一种流行的 GGA）这样的泛函是基于普适的物理原理和模型系统（如均匀电子气）的性质构建的。它不是使用你正在计算的特定分子的实验数据进行调整或参数化的。它在你的体系上的成功（或失败）是对该理论的诚实检验，而不是一个预先设定的结果。

理解数字：能量意味着什么？

当一次 DFT 计算完成后，它会给你一系列数字：一个总能量，以及一组 Kohn-Sham 轨道能。总能量是最大的奖品——它告诉你分子的稳定性。但轨道能又代表什么呢？

在这里，我们再次看到了与 Hartree-Fock 理论的一个微妙但至关重要的区别。在 HF 理论中，Koopmans 定理提供了一个优美而简单的解释：一个占据轨道的能量近似等于从该轨道移除一个电子所需的能量的负值（即电离势）。这是一个近似，因为它假设剩余的 N-1 个电子的轨道不发生改变——即“冻结轨道”近似。

在 DFT 中，相应的论述是 Janak 定理。它是一个精确的定理，而不是一个近似。它指出，轨道能 $\varepsilon_i$ 是总能量关于该轨道占据数 $n_i$ 的偏导数