在上一章《第十章：探索密度泛函理论之平面波基组详解！| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了密度泛函理论基础—赝势与基组。VASP通过自洽场计算方法，通过密度泛函理论求解材料的性质。本章将正式引入这一理论核心，从物理图像出发说明自洽场计算在描述不同材料体系时的关键作用。
Kohn-Sham方程的提出，是基于早期的Hohenberg-Kohn定理。Hohenberg-Kohn定理指出，系统的基态能量是电子密度的泛函（即能量是电子密度的函数），并且基态电子密度决定了外部势能和所有系统性质。虽然该定理奠定了理论基础，但它并未给出具体的能量表达式。
Kohn-Sham方程的核心贡献在于，提出了一个辅助系统：非相互作用电子系统。通过引入一个虚构的无相互作用体系，使其产生的电子密度恰好等于真实体系的基态密度。这一假设极大地简化了多体问题，使其可以通过求解一组单电子方程来近似真实系统的性质。
Kohn-Sham方程本质上是一组特征值方程，形式类似于单电子薛定谔方程：
其中，
是有效势能，由三部分组成：
1.外部势能 (
)：由原子核提供的库伦势能（
）和外加场决定。
2.Hartree势 (
)：描述电子之间的经典静电排斥作用，形式为电子密度的库仑积分。
3.交换-关联势 (
)：最关键的部分，包含了量子力学中的交换（Pauli排斥原理导致）和相关（电子动态关联）效应。该项无法精确求解，通常需要使用近似泛函（如LDA或GGA）来代替。
求解完方程后，系统的总能量E并不等于本征能量之和
。总能量的准确计算公式为：
该公式表明，
本征值并不直接对应真实的电子能级（除了最高占据分子轨道能量HOMO），其物理意义有限。
Kohn-Sham方程的求解过程通常采用自洽场（Self-Consistent Field, SCF）迭代法。具体步骤如下：
1.初始化：给定一个初始电子密度
（通常由原子密度叠加得到）。
2.构建有效势能：根据
计算
。
3.求解特征值方程：解一组单电子方程，得到本征函数
和本征值。
4.计算新密度：根据
更新电子密度
。
5.收敛检查：比较
和
的差异（或能量差）。若满足预设阈值，则迭代结束；否则，将作为新的起始密度继续迭代。
Kohn-Sham方程不仅是理论工具，更是现代材料科学和化学研究的工作马。它广泛用于计算：
分子结构：优化分子几何结构、预测键长键角。
光谱性质：计算红外（IR）、拉曼（Raman）光谱。
电子性质：能带结构、态密度（DOS）、电荷分布。
化学反应：反应路径、过渡态、活化能。
理论背景与意义：介绍了周期性边界条件与平面波
数学形式：介绍了KS方程各项物理意义
求解流程：介绍了自洽场迭代的流程
下一章将正式引入本次教程的核心—自洽场计算。我们将从自洽场理论基础、收敛标准、收敛策略，详细介绍自洽场计算，以及它在处理不同材料体系时的独特优势的应用，敬请期待！
