自洽场

有了总能表达式，下一步就是在程序中高效、稳定地求解 Kohn-Sham 方程。本文讨论自洽场计算的流程与电荷密度更新，利用共轭梯度法求解广义本征值问题，以 RMM-DIIS 为代表的迭代对角化方法及其效率优化，最后介绍用于结构优化和分子动力学的 Hellmann-Feynman 力。

在密度泛函理论（DFT）计算中，**自洽场（Self-Consistent Field，简称 SCF）**的收敛问题是每一个计算化学和凝聚态物理研究者都会遇到的"经典难题"。本质上，DFT 里的 SCF 收敛就是在解一个非线性不动点方程，通过"电荷/势的迭代 + 混合 + 预条件"把电子密度和哈密顿算符做到自洽。本文将深入剖析 SCF 迭代的收敛机制、不收敛的根本原因，以及针对不同体系的具体解决策略。