2.1 Kohn-Sham 方程：映射为单粒子有效势问题

2.1 Kohn-Sham 方程：把多体问题映射为单粒子有效势问题

本节摘要：Kohn-Sham 方程是 DFT 从「理论可行」到「工程可算」的临门一脚。它引入一个与真实体系共享相同基态密度的无相互作用参考系，把总能量拆成动能、外势、Hartree、交换关联四项，使原本高维的多体问题退化为单粒子方程。本节拆解 KS 方程的构造、有效势的三部分、密度重建与自洽循环，并讨论 KS 轨道能量的物理意义边界。

先立靶子：读完这一节你能做什么

阅读完本节，你应当能够：

1. 说出 KS 方程的核心思想：无相互作用参考系 + 相同基态密度
2. 写出 KS 总能量四项分解，并解释每一项的含义
3. 写出 KS 有效势的三部分：外势、Hartree 势、交换关联势
4. 理解自洽循环：密度 → 有效势 → 轨道 → 密度
5. 知道 KS 轨道与轨道能量的物理意义边界（Janak 定理、能隙低估）

问题与直觉

第 1 章的困境可以一句话概括：多体波函数算不动，平均场近似（HF）又抓不住动态关联。HK 定理（下一节细讲）证明了密度可以当基本变量，但留下一个没解的问题——普适泛函 F[n] 没人知道。尤其动能泛函，用密度写出来是什么样？这个问题从 1964 年问到今天，没有答案。

1965 年，Kohn 和 Sham 做了一个绕开这个死结的漂亮动作。他们问了一个看起来很奇怪的问题：能不能找一个无相互作用的参考体系，让它的基态密度和真实体系一模一样？

如果能，那这个参考体系的所有性质都好算——因为电子之间没有相互作用，它的波函数就是单电子轨道的乘积，动能可以精确算。真实体系的难处（相互作用）全被甩到「参考系和真实系的差距」里，这个差距被命名为交换关联能 E_{xc}。这样一来，难题没有消失，但它从「我要直接解多体方程」变成了「我要近似 E_{xc}」——前者无解，后者有大量成功的工程手段。

用一个类比：你想知道一家公司的总资产，但公司内部结构盘根错节，直接审计太复杂。KS 的做法是——先建一个「简化版公司」，结构清晰、账目透明，让它对外展示的资产总额和真实公司完全一致。真实公司的所有「隐藏结构」（复杂股权、交叉持股）被打包成一个黑箱，叫 E_{xc}。你不再审计真实公司，而是不断精修这个黑箱。这个「参考系 + 黑箱」的组合，就是 KS 方程的全部哲学。它的成功也反过来证明了一件事：求解电子结构问题，最实用的路径往往不是「正面对抗复杂性」，而是「把复杂性重新组织到可管理的位置」。理解了这个逻辑，你就能预测第 3 章泛函发展的走向——每一次泛函改进，都是在给这个黑箱塞进更多物理。

核心原理

2.1 总能量四项分解

KS 把真实体系的总能量写成：

E[n] = T_s[n]（无相互作用动能）
     + E_ext[n]（外势：电子与原子核）
     + E_H[n]（Hartree 库仑：经典电子-电子排斥）
     + E_xc[n]（交换关联：所有「剩余」效应）
​

注意 E_{xc}[n] 装了三样东西：交换能（泡利排斥）、关联能（动态关联）、以及真实体系动能与参考系动能之差。最后这一项很关键——它意味着 KS 并没有「假装动能好算」，而是把「动能差」也塞进了交换关联项里。这就是为什么 KS 交换关联比人们直觉上以为的「小」：它补的窟窿不止一个。

2.1b 为什么不做无轨道 DFT

看到 E_{xc} 里混着「动能差」，你可能会问：为什么不能直接把动能也写成密度的显式泛函，彻底甩掉轨道？这确实是 DFT 历史上的一条支线——无轨道密度泛函理论（OF-DFT）。它的思路诱人：没有了 KS 轨道，就没有对角化那一步，计算标度可以降到接近线性，能处理超大体系。

问题在于动能泛函。Weizsäcker 提出的动能密度形式 \tau_W = |\nabla n|^2/(8n) 只对单电子精确，实际体系的动能比它大得多，剩下的差是「关联动能」，无法简单地用密度梯度表达。半个多世纪过去了，无轨道动能泛函的精度始终不够，只能处理近自由电子体系（如简单金属的粗略性质）。KS 方程放弃「纯密度化」的执念、保留轨道来算动能，正是因为这个折中在工程上胜出。这个历史细节解释了本书所有计算都依赖 KS 轨道的原因。

2.2 KS 方程与有效势

对上面这个能量泛函做变分（保持电子数守恒），就得到 KS 方程：

[ 动能算符 + V_ext + V_H + V_xc ] φ_i = ε_i φ_i
​

三个势的含义：

• V_ext：原子核对电子的吸引势，只依赖核位置
• V_H（Hartree 势）：由当前密度产生的经典库仑排斥，V_H(\mathbf{r}) = \int \frac{n(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|} d\mathbf{r}'
• V_xc（交换关联势）：E_{xc} 对密度的泛函导数，V_{xc}(\mathbf{r}) = \delta E_{xc}[n]/\delta n(\mathbf{r})

方程的形式和单电子薛定谔方程长得几乎一样，这是 KS 刻意为之——它把多体问题「映射」成单粒子问题。这个「映射」正是本节标题的意思：一个看似有 N 个独立电子的方程，其实在解一个 N 体问题，只是多体性全藏在 V_{xc} 里。

2.2b Hartree 势的一个细节：自相互作用

Hartree 势描述的是「总密度与自己相互作用」，这必然包含一个物理上荒谬的项——每个电子与它自己的密度也发生了库仑排斥。这个「自相互作用」在精确的交换关联泛函里应该被交换项精确抵消，但近似泛函抵消不干净，就留下了自相互作用误差（SIE）。

SIE 的危害很具体：它让电子倾向于「离域」，因为把一个电子的密度摊薄可以减少自排斥能。后果是 HOMO-LUMO 能隙被低估、电荷转移态被错误地过度稳定、某些局域电子体系（如过渡金属氧化物）被错误描述。这个误差会贯穿第 3 章的泛函讨论和第 6 章的 SIC 校正，在这里先埋个钉子：Hartree 项既是 KS 方程里最大的一项，也是误差的重要来源之一。

2.3 自洽循环：为什么必须迭代

现在你该看出一个循环依赖：V_{xc} 和 V_H 依赖密度，密度又由轨道重建，轨道又由方程解出。这是典型的「鸡生蛋」。解法只能是迭代：

1. 猜一个初始密度（通常用原子密度叠加）
2. 用它构建 V_H 和 V_{xc}
3. 解 KS 方程得到轨道
4. 用占据轨道重建新密度
5. 比较新旧密度：足够接近就收敛，否则混合后回到第 2 步

这个流程和 HF 的自洽循环同构，只是算符不同。第 4 章会详细讲混合策略（Pulay/DIIS）、收敛判据（能量/密度/受力）和不收敛的应对，这里先建立「KS 必须自洽求解」的认知。

2.4 KS 轨道与能量的物理意义：边界在哪

KS 轨道 \phi_i 严格说是「辅助构造」，唯一职责是重建密度。它不是真实电子的波函数。但实践中它太有用了：形状与 HF 轨道很像，可以用来分析分子轨道、成键、能带结构。HOMO/LUMO 也常被用来讨论化学活性。

轨道能量 \epsilon_i 的物理意义更微妙：

• 对于精确的交换关联泛函，HOMO 能量负值等于电离能——这是 Janak 定理
• 对于 LDA/GGA 这类近似泛函，KS 能隙（HOMO-LUMO 差）会系统性低估真实光学带隙，因为忽略了电子-空穴相互作用

这条边界在第 5 章凝聚态、第 6 章会反复出现：用 LDA/GGA 算带隙，得到的是「低估了的能隙」，需要杂化泛函或 GW 校正。

工程实践要点

1. KS 方程是整个 DFT 计算的引擎：无论你用 VASP、QE 还是 Gaussian，内核都在解 KS 方程。理解本节的自洽循环，等于理解了所有 DFT 软件的运行逻辑。
2. 「映射」不是「近似」：KS 方程本身是精确的（在给定 E_{xc} 的意义下）。近似只发生在 E_{xc} 的选择。这意味着：计算不准，先怀疑泛函，再怀疑实现。
3. 能隙要打折扣：用 LDA/GGA 算半导体带隙，心里默认结果会偏低。需要定量带隙时，提前计划用杂化泛函或 GW。
4. 别把 KS 轨道当物理波函数：做定性分析（轨道形状、对称性）可以，做定量结论（占据数、跃迁概率）要谨慎。
5. 自相互作用处处要警惕：凡是电子局域化问题（过渡金属、自由基、电荷转移），都要考虑 SIE 是否在干扰你的结论。可以交叉验证：同样的计算换杂化泛函（引入部分精确交换）再看结果变不变。

关于第 4 点再展开一句：KS 轨道对「形状」的预测通常相当可靠——HOMO 在哪里、对称性如何、跟什么轨道成键，这些定性图像与更高精度方法吻合良好。真正的陷阱在「能量差」：KS 能隙、虚轨道能量这些量，会被近似泛函系统性扭曲。所以在实践里记住一个分工——用 KS 轨道看形状，用更可靠的方法（或谨慎校准）看能量。

一个加深理解的小练习

取一个只有两个电子的氦原子，把它想象成 KS 参考系：两个「虚拟电子」在有效势里独立运动，各自占据一个轨道，重建出的密度等于真实氦原子密度。现在问题来了——这两个虚拟电子的动能和真实氦的动能差了多少？这个差值就是 E_{xc} 里的「动能差」部分。带着这个问题去读第 3 章 LDA 的均匀电子气图像，你会更清楚交换关联项到底在补什么。这类「拿单原子做思想实验」的方法，是理解电子结构理论最便宜也最有效的手段。

⚠️ 常见坑：直接用 LDA/GGA 的 KS 能隙当「带隙」写进论文，与实验对比发现差一半，又找不到原因。根因就是 KS 能隙的固有低估，不是程序 bug。

💡 关键直觉：KS 方程的成功可以用一句话概括——把「不知道的」都塞进一个小而局部的项，剩下的全部精确计算。这个「隔离未知」的策略贯穿整个计算科学，DFT 只是最成功的例子之一。

概念性伪代码：KS 自洽循环

def ks_scf(initial_density, xc):
    n = initial_density
    for step in range(max_steps):
        V_eff = V_ext + hartree_potential(n) + xc.potential(n)
        orbitals, eps = solve_ks(V_eff)
        n_out = sum(orbital_density(o) for o in occupied(orbitals))
        if norm(n_out - n) < tol:
            break
        n = mix(n, n_out)   # 简单混合或 Pulay
    return n, total_energy(n)
​

回顾自查清单

• 核心思想：引入无相互作用参考系，共享真实体系的基态密度
• 四项分解：动能、外势、Hartree、交换关联——前三项可算，难题全在 E_{xc}
• 有效势三部分：外势 + Hartree 势 + 交换关联势
• 自洽循环：密度 → 势 → 轨道 → 密度，必须迭代
• KS 轨道：辅助构造，用于重建密度，也可做定性分析
• Janak 定理：精确泛函下 HOMO 负值 = 电离能
• 能隙低估：近似泛函的 KS 能隙系统性偏低，需杂化或 GW 校正
• SIE 隐患：Hartree 项的自相互作用靠交换项抵消，近似泛函抵消不干净
• 无轨道 DFT：动能泛函无法简单密度化，这是 KS 保留轨道的原因
• 隔离未知：把难题压进局部的小项，是 KS 的方法论精髓

下一节回头补上理论地基——HK 定理如何证明「密度决定一切」，以及它和 KS 方程如何衔接。