静态自洽电子结构计算

静态自洽电子结构计算是所有 DFT 任务的底层循环——它不优化原子位置，只把电子密度在给定结构上迭代到自洽。听起来简单，却是新人卡得最久的环节：不是在”电子步不收敛”，就是在”电荷密度震荡”。我当年为了一个金属表面体系，硬生生调了一周混合参数才跑通。

一、背景：为什么 SCF 是绕不开的地基

金属表面体系是我当年噩梦，电子步死活不收敛，调了一周才摸到门道：本质是混合参数没配。这类问题在新人手里极常见，因为教程里很少讲混合算法，大家默认能跑，跑不动就懵了。

另一类高频不收敛来自初始电荷密度太差。拿随机初猜去解强关联或大带隙体系，前几十步能量剧烈震荡，经常直接发散。我们遇到难体系，先上一轮粗网格原子电荷初猜，把震荡压下去再细化。

能带、态密度、电荷密度、受力，全建立在自洽电荷密度之上。SCF 没收敛，后面一切性质都是建立在沙子上。我们在接任何新体系时，第一件事就是先单独把静态自洽跑稳，再谈别的。一个 静态自洽电子结构计算 起不来，往往不是机器问题，而是初猜、混频或对称性配错了。

二、核心原理：自洽在迭代什么

SCF 的关键是混合。金属体系电子步极易震荡，减小 AMIX 或改用 Pulay 混合常能救回来；难收敛体系我们提前备好 charge mixing 与 occupation 方案。理解”混合在迭代什么”比盲目调参数有用得多。

占据方案（smearing）和收敛是绑定的。金属用 Gaussian 或 Methfessel-Paxton，绝缘体用 Fermi 展宽或固定占据，类型选错带边会抹平或震荡。我们按体系导电性查表选 smearing，不让默认替我们担风险。

SCF 本质是反复求解 Kohn-Sham 方程：用当前电荷密度算势→解本征方程得新波函数→构造新电荷密度→和旧的混合。当输入和输出电荷密度差小于阈值，就自洽了。关键在”混合”这一步——混合系数和算法决定了收敛快不快、稳不稳。我们内部对金属体系几乎必上 Davidson 配合适的 mixing，对难收敛体系会提前准备 charge mixing 与 occupation 方案。

三、关键技术要点：三个决定收敛的因素

初猜和对称性也常被忽视。用原子电荷或前一步波函数做初猜比随机快得多；对称性设错会让 k 点权重异常，自洽差一口气。我们有一套”卡住就按这个顺序排查”的清单，基本覆盖九成不收敛。

NBANDS 和收敛也有关系。导带留太少，金属或近金属体系电子步容易卡在窄能带附近震荡。我们对近金属体系额外开 NBANDS 再测，往往比硬调 AMIX 更治本。

第一是初猜，用原子电荷或前一步波函数做初猜，比随机起步快得多。第二是混合参数，金属体系电子步极易震荡，减小 AMIX 或改用 Pulay 混合常能救回来。第三是占据数处理，带隙小的体系用 smearing 要谨慎，既要收敛又不能平滑掉真实能带。我们做 SCF 有套”卡住就按这个顺序排查”的清单，基本能覆盖九成不收敛。

四、实操流程：从初猜到自洽密度

从结构到收敛电荷密度，我们固定流程：原子电荷初猜→小 AMIX 试收敛→Pulay 混合接管→smearing 标定→对称性与 k 点复核→输出电荷与波函数。卡住时按清单逐项过，不靠反复瞎试。

交付前我们检查电荷密度的物理合理性：Bader 电荷符号、磁矩大小、总能量单调性。这些量一旦反直觉，往往是混合或初猜没配好，比”能量看起来降了”更可靠。

落地步骤：准备结构→设初猜与混频方案→跑静态自洽→看能量与电荷密度收敛曲线→确认无电荷震荡再进后续。当你把 静态自洽电子结构计算 当成”必须先过的关”而非”顺手的一步”，后面做能带和 Bader 分析时就会省掉大量返工，因为底座已经稳了。

五、常见问题：SCF 阶段的典型翻车

最频繁是金属电子步震荡不收敛，根因在混合参数；其次是对称性设错 k 点权重异常；还有 smearing 太大抹平窄带隙。按清单逐项排查都能定位，不是玄学。

另一类典型坑是 ALGO 选错，近金属体系用默认算法死活不进自洽。我们按体系查表换 ALGO，远比反复调 AMIX 高效，也少走弯路。

最频繁的是金属体系电子步震荡不收敛，根因几乎都在混合参数。其次是对称性设错导致 k 点权重异常，自洽永远差一口气。还有人 smearing 设太大，把窄带隙抹平了，态密度看着顺眼实则失真。这些都不是玄学，按清单逐项排查都能定位。

六、复盘总结：SCF 的三条经验

回过头看，静态自洽电子结构计算最怕”默认能跑”的心态。金属、强关联、大带隙三类体系，默认参数必翻车，提前备方案才是正道。

我们的经验：初猜 + 混合 + smearing 三件套配齐，九成不收敛都能在清单里定位。把它当工程问题，而不是玄学调参。

回过头看，静态自洽电子结构计算最忌”硬等它自己收敛”。初猜要给好、混合要调对、占据要合理。我们交付的计算，SCF 日志必附收敛曲线，客户能看清密度是怎么稳定的。把这三件事做扎实，后续所有电子结构性质都有了可信起点。