说明：本文华算科技主要介绍计算不收敛的常见来源，包括 SCF 残差振荡、初猜密度、金属占据、混合参数、结构初态、力和应力、k 点与截断能，以及怎样从输出症状判断问题发生在哪个环节。
计算停在中途时，输出文件里的报错通常只是终点信号。能量不降、电子步反复震荡、残余力长期停在高值、晶胞应力忽大忽小、磁矩在几个状态之间切换，这些现象对应的计算环节并不相同。把它们都写成“参数没调好”，会遮蔽模型问题、电子态问题和数值问题之间的差别。
计算收敛到底指哪几个量达到稳定？
Kohn-Sham DFT 的电子步固定原子坐标，反复更新电子密度、有效势、轨道占据和总能量。SCF 收敛指输入密度和输出密度达到自洽，在程序输出里常表现为 总能变化、密度残差或本征值变化 低于给定阈值。这个判断只属于当前原子构型，不自动保证原子已经处在势能面极小点。
结构优化在电子步之外多了一层几何搜索。每一次离子步结束后，程序根据力和应力更新原子位置或晶胞参数。电子步收敛、力收敛、应力收敛和目标性质收敛不是同一个停止条件：静态能量排序关心 meV 量级能量差，声子计算关心残余力噪声，slab 功函数还会受到真空电势平台和偶极修正影响。
输出症状可以先按计算环节分组。SCF 残差上下摆动 多半要检查混合、展宽、初始电荷和磁态；最大力不降 则要检查初始结构、约束、优化算法和电子步力噪声；总能已经平滑，但目标吸附能或缺陷形成能随 k 点变化，说明数值参数还没有达到目标差值尺度。同一个“不收敛”字样背后可能是三个不同问题。
电子步为什么会在振荡、平台和假收敛之间切换？
铝 fcc 四原子晶胞这类固定构型中，SCF 残差随循环次数下降时，残差 L2 范数持续变小，说明新旧密度正在靠近同一个固定点；总能变化同步缩小，当前构型下的电子态也更接近稳定解。正常下降曲线通常同时具备残差减小、能量变化减小和占据不剧烈跳变，三者缺一项，后续力和应力就可能带有电子求解噪声。
硅 fcc 晶胞的阻尼对比曲线显示，程序把上一轮输出密度的一部分加回输入密度时，混合太强会越过固定点，下一轮又被拉向相反方向，残差曲线就出现先降后升的振荡。阻尼迭代把每轮密度更新幅度压小，牺牲部分速度，换取更平滑的电子态接近过程；大真空 slab、带电缺陷和金属表面比小型绝缘体更容易暴露这类问题。
金属、小带隙半导体、表面态密集的 slab 容易出现占据跳动。费米能级附近有很多能量相近的态时，微小势场变化会改变轨道占据，电子密度随之改变，SCF 循环难以稳定。展宽温度改变的是费米能级附近占据函数斜率；展宽过小会让占据像阶跃一样敏感，展宽过大又会改变自由能修正和能量排序的读法。
初猜密度与目标密度相差很大时，电子步会先消耗大量迭代去修正长波密度误差。低维材料、异质界面、强极化 slab、带电体系和磁性体系中，孤立原子密度叠加得到的初猜往往离收敛密度有明显空间差异。这里应检查 初始电荷、初始磁矩、电荷态、对称性约束和前一轮 WAVECAR/CHGCAR 是否匹配当前结构。
过渡金属氧化物、含缺陷氧化物、自由基吸附物和单原子位点中，磁性和强关联效应会给出多个接近能量的电子解；初始磁矩设置 会把计算带向不同局域态。某个电子态形式上满足阈值，却不是目标研究对象的最低能或可比电子态。比较吸附能、形成能或能垒前，应保存末步磁矩、占据、总电荷和自旋态信息。
结构优化为啥会出现力不降/能量反复跳
几何优化失败不一定来自 SCF。初始原子距离过短、吸附物放在高排斥位点、slab 层数太薄、真空层不足、晶胞形变自由度设置错误，都会让势能面局部斜率很大。最大力长期集中在少数原子上，通常指向局域结构问题；所有原子力都乱跳，则更可能是电子步噪声、约束冲突或优化步长过大。
几何搜索算法根据当前力给出下一次位移时，步长太小会让结构移动很慢，步长太大又可能跨过局部低能区，下一步能量随之升高。线搜索曲线显示，同一个搜索方向上，能量随步长并非线性下降，步长选择和势能面曲率共同决定一次离子步是否有效。软模式、表面重构和分子吸附构型转换 时，这种曲率变化尤其明显。
电子密度残差没有稳定时，Hellmann-Feynman 力会带有噪声；力阈值设得很紧，而 EDIFF、k 点、截断能或展宽 仍然偏粗，结构优化就会在阈值附近反复摆动。做声子、弹性常数、NEB 初末态或弱吸附能差时，力噪声必须低于目标物理量允许的误差，否则末态结构看似完成，后续图谱会把数值噪声当成材料信号。
对 slab、缺陷和界面模型，结构自由度本身也会改变收敛表现。固定底层过多会把表面应力集中到少数可动原子，固定太少又可能引入整体漂移；非对称 slab 如果没有处理面外偶极，平面平均电势和电子密度更新都会受到周期镜像影响。带电缺陷还涉及 背景电荷、超胞尺寸和电荷校正，几何优化中的力和静态能量都要按同一电荷态保存。
数值参数/任务目标怎样改变收敛压力
大胞金属、表面和低维模型中，Kerker 预条件处理密度更新中不同波矢分量的响应，长波电荷振荡会让 SCF 残差迟迟不降；Kerker 预条件通过改变 |G| 较小分量的更新权重，降低 charge sloshing 对迭代的破坏。混合参数和预条件不是普通精度旋钮，而是电子密度迭代路径的一部分。
金属体系的费米面采样不足时，k 点和截断能会影响总能、力、应力和态密度，相邻电子步中的占据和总能可能同步波动；截断能不足会让力和应力误差进入结构优化。对于吸附能、缺陷形成能和相稳定排序，参数误差应低于被比较能量差；若候选结构相差只有几 meV，粗 k 点或低截断能足以改变排序。
含半芯态元素、重元素 SOC、强磁性 d/f 电子、带电缺陷和高压结构中，赝势选择会改变平面波需求、局域磁矩和能量排序。某些体系用默认赝势可完成普通结构优化，但在高氧化态、短键长或磁态比较中出现能量漂移，后续静态计算就应保持 同一赝势、同一泛函、同一 U 值和同一展宽方案。
粗优化、精优化、静态能量比较、声子计算和 NEB 对“够用”的要求不同：粗优化只需要把结构带到合理区域，精优化关注最大力和应力；静态能量比较关注电子步阈值、k 点、截断能和展宽修正；声子计算关注力常数噪声；NEB 关注初末态力、插值路径和图像间距。收敛检查应绑定目标输出量：DOS 关心费米能级附近态密度是否稳定，功函数关心真空电势平台，AIMD 关心每一步电子迭代失败是否积累成轨迹偏差。
排查顺序可以从输出症状出发，但每一步都要对应可核对的量。SCF 振荡看残差、占据、混合和展宽；力不降看最大力位置、约束、步长、初始键长和电子步噪声；能量排序不稳看 k 点、截断能、赝势、磁态和电荷态；后处理图谱异常看静态结构是否继承了优化末态的电子条件。正文到这里停在一个具体检查习惯上：每个最终结构旁边保留 OSZICAR、OUTCAR 末步力、磁矩、应力、k 点和截断能记录。
