为什么DFT 计算不收敛？

说明：本文华算科技主要介绍 DFT 自洽迭代失稳的数学来源、诱发长波电荷振荡和电子态切换的材料特征，以及残差轨迹包含的诊断信息。

电子步未达阈值，常称为“不收敛”。残差轨迹可能呈现密度在两组状态间往返、长波误差持续放大或磁矩反复换向，或由离子移动后的电子态重排引起。相似迭代步数可以来自不同反馈机制，诊断时应追踪残差的形状、空间分布和电子谱，而非只改某个输入参数。

DFT自洽的对象为什么是非线性固定点

Kohn–Sham DFT 把输入电子密度 ρ_in 写入有效势，求解单电子轨道和占据数，再产生输出密度 ρ_out。两者的差 R(ρ)=ρ_out−ρ_in 是自洽残差；当 R(ρ)=0 对应同一个密度固定点时，势、轨道、占据和密度才彼此相容。程序把新旧密度按一定比例混合，实际操作是沿残差方向更新下一步输入。

迭代能否靠近固定点，取决于密度扰动经过“密度—势—轨道—新密度”循环后是缩小还是放大。固定点附近可用雅可比算符描述这种响应；若某些本征模式的更新幅度过大，残差便会交替变号或持续增长。混合参数控制更新幅度，预条件器改写不同空间模式的更新比例，两者处理的是反馈环路的谱，而非单纯增加计算次数。

图1中的三条响应曲线对应三类电子谱。金属在费米能级 E_F 处保留有限态密度，长波极限 q→0 时库仑核 4π/q² 会把密度响应强烈放大；半导体和绝缘体的带隙抑制低能电子—空穴激发，介电响应保持有限。同一套混合规则面对的反馈增益随材料谱而变，金属大超胞通常比有隙小晶胞更容易出现长波振荡。

电子自洽与离子优化属于内外两重迭代。内层寻找给定原子位置下的电子固定点，外层依据力移动原子。若电子残差尚大，力中会夹带密度误差；外层优化器随即依据带噪声的力产生新构型，下一离子步又从偏离固定点的密度起跑。能量变化偶尔很小，只约束相邻步的一个标量差值，无法排除局域残差、磁矩和力仍在摆动。

自洽方程可能存在多个固定点。不同自旋排布、轨道占据或电荷局域方式都能满足方程，并形成近能驻点。迭代历史既影响速度，还可能决定吸引域；前一步留下的密度、磁矩和混合历史会把后续计算引向某个解，或让迭代在几个吸引域之间切换。

哪些材料特征会放大密度误差？

金属长超胞中的电荷晃动常由最小波矢模式主导。沿长度 L 方向可表示的最小 q 约为 2π/L，L 增大后 q 变小，金属介电响应随 1/q² 急剧增大。误差并非均匀分布在全部傅里叶分量中，少数跨越整个晶胞的长波模式就能控制总残差，局域成键区域却可能已经接近自洽。

图2把 Al 超胞中介电算符的五个主导模式画成实空间分布，波长都与晶胞尺度相当。简单线性混合给所有 q 分量相近的更新权重，长波分量的更新幅度便可能过大；Kerker 型预条件用 q²/(q²+q₀²) 压低小 q 更新。晶胞变长时，响应最大的模式会继续向更小 q 移动，固定阻尼对体系尺寸的适应范围随之收窄。

名义上的半导体表面有时呈现金属式响应。未钝化悬挂键、缺陷带、掺杂能级或部分占据态贴近 E_F 时，带隙材料的局部区域仍有低能电子激发。决定小 q 屏蔽的是费米能级附近的实际谱，不能由体相材料名称替代。表面重构或吸附过程若让该态跨过 E_F，相邻离子步适用的混合权重还会突然改变。

图3中，裸 Si(111) 薄膜在 E_F 附近出现表面态，部分电荷集中于两侧表面；H 钝化消除悬挂键后，E_F 上方形成清晰能隙。两种模型的元素组成很接近，低能响应却不同。钝化改变的是迭代映射的谱条件，这也是表面化学状态会改变收敛表现的直接原因。

金属/真空、金属/氧化物和电极/分子模型还具有空间不均匀屏蔽。图4的 Al 区域在 E_F 附近保留连续 LDOS，真空区没有电子态，SiO₂ 区域呈现能隙。单一 q₀ 无法覆盖三个区域的响应：针对金属设置的长波抑制可能让绝缘区域更新过慢，针对绝缘体设置的权重又可能放大金属侧电荷晃动。

过渡金属 d 态、稀土 f 态和强自旋极化体系中的窄能带让少量势变化即可重新排列占据，数个磁态或轨道序的能量差又很小。残差此时包含占据切换与自旋切换，并非只有长波电荷模式；只压低小 q 分量，仍可能留下局域轨道反复换态。

振荡、停滞和换态来自什么数值过程

残差轨迹可区分反馈方向。近似两步一循环、符号交替而幅值相近，通常对应某个模式更新过量；残差单调增大，更新方向或步长已越出吸引域；下降一段后形成平台，常见于预条件器没有改善最慢模式，或数值噪声接近目标阈值。“达到最大电子步”只是停止原因，振荡、发散与平台才对应可区分的数值过程。

图5的 8 原子 Al 模型只改变混合方式与阻尼，轨迹可缓慢下降、失稳，或在自适应步长下恢复。阻尼过小会让全部模式移动缓慢；阻尼过大则放大高响应模式。预条件器负责压缩不同模式的响应跨度，阻尼再决定整体更新尺度，两项功能不能互相替代。

能量差、密度残差、势残差和本征值变化观察的是不同对象。总能量对密度误差在驻点附近可能呈二阶变化，残差仍有可观幅值时，能量差已经显得很小。结构优化更应核对力与电子残差是否同步下降；若能量平缓而力呈锯齿，离子更新正在读取尚未自洽的电子响应。

Fe₂CrGa 的对照计算呈现另一种轨迹：采用 Anderson 历史加速后，多条轨迹可以快速下降，却到达能量不同的固定点。过渡金属窄 d 带与自旋构型造成多个近邻驻点，历史外推会改变进入哪一个吸引域。残差达标只确认方程在该状态下闭合，基态判断仍需比较总能、磁矩、占据和结构特征。

离子步中的突然反弹也有可辨认来源。几何变化若让能带穿过 E_F、局域磁矩翻转或对称性解除，上一离子步积累的混合历史描述的是旧映射，外推方向可能立即失效。反弹恰好伴随占据数或磁矩突变时，清除旧历史、缩小离子位移并重新建立电子响应，可避免旧映射的外推方向继续作用于新电子态。

怎样按失稳来源恢复可信自洽解？

诊断记录应保留每步残差范数、总能、费米能级、磁矩、占据和最大力，并标出离子构型变化。固定几何下的电子失稳，可从残差傅里叶分量判断是否由小 q 模式控制；局域残差集中在界面、缺陷或过渡金属位点时，残差可能来自空间屏蔽或占据切换。残差位置区分界面与体相来源，残差波长区分长波和局域模式，电子步编号不含这些特征。

小 q 分量支配时，应削弱长波更新或采用能描述金属屏蔽的预条件；金属与真空、绝缘体并存时，空间依赖的介电模型比均匀模型更贴近响应分布。若各 q 分量都发生过冲，降低整体阻尼或用受控线搜索限制步长。轨迹平台对应最慢模式没有获得足够更新，继续统一减小阻尼反而可能延长平台。

图7比较 Na 链与真空构成的非均匀模型。Anderson 方法在 32 原子模型中约 60 步下降，在 64 原子模型中 100 步内仍未达到同一残差范围；Kerker 与 Anderson 组合也随尺寸显著变慢。描述空间介电差异的椭圆预条件在两种尺寸下都约用 30 步。体系尺寸依赖可作为长波屏蔽模型失配的证据，因为增大的真空和更小的 q 改变了最慢响应模式。

占据或自旋频繁切换时，展宽只是在有限电子温度下平滑费米面，应检查展宽后的自由能项、磁矩和目标零温能量是否符合研究对象。对称性约束、初始磁矩和电荷态会选择可访问的电子状态。改变这些条件会改变所求固定点；它们属于电子状态的定义条件，作用不限于调整迭代速度。

得到低残差后，应从独立初始密度或不同合理磁态重复计算，并比较总能、局域磁矩、轨道占据、带隙和力。改变初始条件后回到同一电子态，说明该固定点具有较大的吸引域；得到多个自洽解时，应按研究条件比较其能量和序参量。可信结果要求低残差与物理状态可复核，二者共同区分真正的电子解、偶然停滞和亚稳态。