从占据到磁矩：SCF 不收敛的四类物理干扰与对应修复策略

说明：本文华算科技主要介绍 SCF 电子步不收敛时，能级占据、长程静电响应、轨道求解误差与磁性多解分别怎样干扰电子密度自洽。

SCF 电子步反复求解的是什么？

SCF 是 self-consistent field，即自洽场。Kohn–Sham DFT 在固定原子坐标和晶胞下，用电子密度构造 Hartree 势与交换相关势，再求轨道和能级。轨道按占据数加权得到新的电子密度，输入密度与输出密度相容时，这一组势、轨道和占据才构成自洽解。电子迭代次数记录数值求解进度，不对应电子真实运动的时间。

电子密度可简写为 n(r)=Σ_if_i|ψ_i(r)|²；周期体系还包含 k 点权重和自旋求和。轨道形状 ψ_i 与占据数 f_i都参与密度构造。即使轨道变化很小，几个近简并能级交换占据，也可能使局部电荷明显改变。下列占据函数与展宽核用于构造这类加权求和，展宽参数决定能级附近占据变化的平滑程度。

自洽残差写成 R=n_out−n_in，衡量同一次电子更新前后的密度差。相邻两步总能差衡量的是另一个量：能量沿某些密度变化方向较平坦时，总能差可以很小，密度残差仍较大。增加电子步上限只延长迭代；若残差幅度长期不减，更多电子步仍会重复相似更新。

多数程序还要在给定有效势下迭代求解本征态。固定势下的轨道残差与更新有效势后的密度残差对应两个求解任务。前者尚未降低到当前密度更新所需的误差范围，后者就会使用带有轨道误差的输出；两类残差应按各自定义核对。

能级占据为什么会反复改变？

在零温整数占据下，费米能级以下的态占据、以上的态空置。金属、小带隙材料和含缺陷能级的体系中，能级间隔可能小于一次有效势更新造成的能量移动。两个态一旦换序，占据随之交换；若它们分别局域在不同原子或轨道上，输出密度便会突变，下一次有效势又可能使能级换回原顺序。

载流子注入的 InGaN 氮 Frenkel 缺陷计算提供了具体实例：沿氮原子位移，氮空位 V_N 与氮间隙 N_i 的缺陷能级接近并交叉，受约束 SCF 在交叉区未能收敛。局域态电荷密度显示，两侧低能态的空间归属发生交换。原研究在该区使用分数占据获得连续结果；这里的约束激发态场景与普通基态计算应分别理解。

Fermi–Dirac 占据 f(E)=1/{1+exp[(E−μ)/σ]} 使近费米能级的占据连续变化，其中 μ 是电子化学势、σ 是能量宽度。展宽减小占据对微小能级移动的突变响应，但会改变参与密度构造的态。最高计算能带若仍有明显占据，遗漏的高能尾部会影响电子数与化学势求解；空带数量应覆盖所用占据函数的尾部。

Fermi–Dirac 展宽通过电子自由能 E−T_eS 中的熵项参与变分求解，其他展宽具有各自的广义熵修正。展宽过大造成的能量与受力偏差可以在 SCF 完成后继续存在。六方 Zr 的例子中，大 Gaussian 展宽让弹性常数随 k 点密度较快趋于平坦，平坦值却偏离小展宽极限。

较稀 k 点网格可能漏掉费米面附近的细小能带结构，极小展宽又会使离散能级的占据变化尖锐。k 点采样与展宽应围绕同一目标能量差、磁矩或力比较，并区分“每个网格都完成电子自洽”与“网格加密后目标量已稳定”。粗网格上的自洽解可以存在，采样误差本身未必造成 SCF 发散。

密度更新和轨道求解为何会相互干扰？

在倒易空间，非零波矢 q 处的电子密度误差通过含 1/q² 因子的库仑核改变 Hartree 势，长波分量对势的影响可显著放大。金属电子又能强烈响应这种势变化，一次更新可能过度修正原误差，使密度在连续电子步中交替偏向相反方向。大超胞允许更小的非零 q，这类长波振荡常更难抑制。

MoS₂ 薄膜和富磷 GaP(100) 表面把材料区与真空区放进同一周期晶胞；材料内部的屏蔽与真空中的响应不同，半导体的长波响应也有别于均匀金属。单一混合系数会同等缩放多种密度误差，可能出现某些分量已很小、另一些分量仍持续振荡的情况。

混合更新可写成 n_new=n_in+αPR，其中 α 控制更新幅度，P 是预条件算符。阻尼缩小整体更新，预条件按误差模式调整响应。针对金属长波响应的 Kerker 形式与针对有限介电屏蔽的 Resta 类形式具有不同假设；表面、绝缘体和金属混合体系的屏蔽差异，会影响预条件后的残差下降速度。

给定有效势后，本征态求解要求 Hψ_i−ε_iψ_i 足够小。不同能带的残差下降速度可以相差很大，少数尚未收敛的态若参与占据，就会扰动输出密度。hBN 的固定哈密顿量计算中，改变本征求解器的移位策略即可改变这些轨道残差的下降轨迹；这里尚未涉及有效势的外部更新。

SCF 初期允许较粗的轨道求解以节省计算，接近自洽时则要降低内层误差。内层容差、舍入误差与外层目标误差不匹配时，密度残差可能停在平台，外推历史也可能持续累积受噪声影响的更新方向。固定粗网格或有限截断能仍可产生自洽解，但该解含离散误差；只有当数值噪声或求解容差限制残差继续下降时，才构成这里讨论的停滞。

初始模型中异常短的原子接触、错误的总电子数或不合适的自旋约束，会显著改变局域势和低能谱。沿用旧电荷密度时，晶胞、原子组成和电子态约束应与当前模型相容。几何突变后，旧的密度混合历史未必仍描述当前响应；重新建立电子初猜和更新历史，处理的是变化后的电子问题。

磁态与局域占据为什么让初猜影响结果？

自旋极化计算分别更新 n_↑ 与 n_↓，局域磁矩由两者的差决定。过渡金属或稀土体系中，磁矩大小、磁序和局域轨道占据会反过来改变有效势。若迭代在几种电子态之间切换，总能和磁矩可能反复跳变；若进入某个局部极小值，SCF 又可能很快结束。两种轨迹对应不同的数值状态。

DFT+U 对局域 d、f 轨道的占据矩阵施加能量修正，同一 U、同一磁序下仍可能存在多个自洽占据解。UO₂ 的占据矩阵控制计算显示，不同初始占据可得到能量不同的电子态。应在相同泛函和相同 U 下比较这些解；跨 U 的绝对总能变化同时包含能量模型改变，不能用来直接给不同 U 排列物理优劣。

更改 U 或关闭自旋极化会改变所求电子态，即使随后顺利收敛，也应重新核对原研究对象。多初猜比较要使用一致的数值条件，并记录局域磁矩、占据矩阵与相对能量。人为固定的占据或对称性会限制允许搜索的电子态；释放约束后，密度和磁矩是否继续变化，应在同一模型下重新求解。

单层 CeI₂ 的计算显示，相同铁磁排列可以对应不同 4f 电子态。在所用 DFT+U 模型和未含 SOC 的比较中，S₀ 与亚稳态 S₂₀ 的自旋密度形状、轨道能级和能带均有差别。总磁矩接近时，局域轨道占据仍可能不同，能量已稳定的结果也应保留其电子态身份。

S₀ 的 Ce 4f 态在相对费米能级约 −4.1 eV 处形成窄带，S₂₀ 的 4f 能带分布不同；局域 4f 占据改变还会影响与其杂化的 Ce 5d 态。这组模型中两态的带隙分别约为 386 meV 和 336 meV，差别来自已经自洽的不同电子态。加入 SOC 后还需重新确定轨道态与磁化方向，能量比较应限定在一致的相对论处理下。