DFT三类收敛判据解析：电子步、能量、力收敛的信号读数与适用场景

说明：本文华算科技主要介绍能量收敛、力收敛和电子步收敛分别对应的计算信号、常用读数和适用范围。

能量/力/电子步收敛的基础对象是什么

一次 DFT结构优化里，程序反复处理两个对象：在固定原子位置上求电子基态，再按能量对原子坐标和晶胞自由度的导数更新结构。电子步收敛描述电子密度、势场和占据是否达到自洽；能量收敛描述相邻电子步或相邻离子步的 E_tot 变化；力收敛描述每个原子沿各自由度受到的残余力是否降到给定阈值附近。三个量共用同一套计算输出，却不回答同一个问题。

电子步的典型信号包括密度残差、能量变化、特征值占据和混合过程是否振荡；离子步能量变化来自一系列已完成电子自洽的构型；力来自 E_tot 对原子位移的导数。能量差很小、最大力很大，常见于势能面平坦但结构仍沿某个软方向有斜率；电子残差尚未稳定、力数值已经变小，则可能把电子噪声带入后续离子步。

VASP 语境中的 EDIFF、EDIFFG 属于阈值名。读结果时应把阈值名还原成物理对象：EDIFF 常限制电子步能量变化，EDIFFG 取负值时常限制最大力，取正值时常限制离子步能量变化。同一个“收敛”词在不同输出行里指向不同停止条件，所以判断结构是否可用于能带、声子、吸附能或过渡态计算，不宜把某一个阈值名当成全部依据。

输出文件里的时间顺序也要分开。电子步是一组固定构型上的自洽迭代，离子步是一组结构更新后的能量、力和应力记录。每走到新的离子步，电子密度又从新的原子势场出发重新自洽；上一轮电子步通过，并不会自动保证下一轮结构更新后的电子态同样平滑。电子内循环和离子外循环的序号应分开读，否则容易把自洽误差误作结构能量变化。

电子步收敛看哪类自洽信号？

电子步收敛的对象是 ρ_in 到 ρ_out 的自洽闭合。固定原子坐标后，程序由当前密度构造有效势，求解 Kohn-Sham 方程，再得到新的电子密度。残差下降、占据稳定、总能变化缩小 共同限定一次 SCF 循环是否到达电子基态附近。金属、低带隙半导体、带电缺陷、磁性过渡金属氧化物和大真空 slab 中，低频电荷摆动会让残差曲线出现平台或振荡。

电子步的迭代次数本身没有精度单位含义。很少的迭代数可能来自初始电荷接近目标解，也可能来自阈值偏松；很多的迭代数可能来自混合参数不合适、占据在费米能级附近反复改变，或局域磁矩在多个电子态之间切换。同一残差阈值下的迭代步数 反映的是电子自洽求解难度；结构弛豫是否完成，还要读后续最大力和离子步能量变化。

残差在阈值附近上下跳动时，每一个离子步拿到的 E_tot 和 F_i 都带有电子求解误差，后面的能量和力会随电子解抖动；在软声子、弱吸附或小能量差比较中，这类误差可能和目标能量差处在同一数量级。静态单点计算、结构优化和声子计算对电子步阈值的容忍度不同，声子和小能量差通常要求电子噪声低于目标信号。

金属体系用较大的展宽后，SCF 残差可能下降得平滑，但自由能、零温外推能和真实占据态并非同一输出量；半导体和绝缘体若出现费米能级附近的异常分数占据，通常要回看带隙、磁矩和初始电荷。电子步收敛还要核对占据是否稳定，不宜只盯住最后一行能量差。

能量收敛看总能还是能量差？

能量收敛常被读成“总能不再变化”，但材料计算里多数结论比较的是 ΔE：缺陷形成前后、吸附前后、不同构型之间、相邻反应中间体之间的能量差。固定结构里的电子步 ΔE 代表电子自洽停止条件；相邻离子步的 ΔE 代表几何优化过程中的势能下降幅度。同样写成 eV 的两个 ΔE，来源可能一个在电子内循环，一个在离子外循环。

能量基准会受截断能、k 点、网格、赝势、展宽和真空层影响。相同方向的系统误差有时在两个相近模型之间部分抵消，ΔE 可先稳定；不同化学组成、不同电荷态、不同超胞大小之间的比较，误差抵消幅度会下降。能量收敛应围绕目标能量差设定；吸附能、相对构型能、缺陷形成能、迁移能垒各自对应不同的误差敏感项。

静态单点的能量收敛偏向电子基态读数，结构优化中的能量变化还受坐标更新步长影响。算法采用线搜索、阻尼动力学或准牛顿更新时，某一步 ΔE 变小可能来自步长缩短；吸附构型、表面重构和缺陷迁移更应记录关键键长、配位环境和下一步 |F_max|。能量停顿不同于几何平衡，这一点在软模和弱相互作用模型里尤其明显。

能量曲线接近平坦时，相邻离子步的 ΔE 很小仍可能伴随某个软模方向的残余斜率。表面吸附物、轻元素氢、二维材料的面外声学模、铁电软模和低频分子转动，常让能量变化比力变化更早变小。|F_max| 仍偏大时，能量收敛只能说明势能下降慢，结构坐标并未到达所有自由度的平衡位置。

势能面很平缓时，较大的位移也可能只带来很小的能量变化；势能面很陡时，微小位移就能给出明显的力。能量变化看高度，力看斜率。这一区别解释了为什么表面扩散、分子转动和软模畸变中，ΔE 与 |F_max| 经常不同步。

力收敛怎样限定结构是否还能移动？

力收敛直接检查原子沿坐标自由度受到的残余驱动力。常用读数包括最大力 |F_max|、RMS force、固定原子是否排除在统计外、晶胞优化时的应力张量 σ。|F_max| 的单位通常写成 eV Å^-1，阈值大小应跟后续任务匹配：普通结构筛选可宽一些，声子、NEB 初末态、吸附能微小差异和弹性常数通常要求更低的残余力。

平均力小，可能只是大部分原子已经平衡；某个吸附氢、表面氧、缺陷邻近原子或晶胞一侧的原子仍可能保留较大的局部力。局部最大力决定结构是否还会沿特定坐标移动。做 slab 或界面模型时，冻结层不应参与最大力阈值判断，非冻结区域的 |F_max| 才对应可弛豫结构的剩余驱动力。平均力和最大力给出的结构信号不同，二者在低对称缺陷和单分子吸附模型中差别尤其明显。

体相晶格常数、二维材料面内应变和外延界面模型中，应力张量会和原子力一起进入弛豫判读；原子坐标已经接近平衡时，晶胞仍可能承受拉伸或剪切应力。ISIF 类自由度改变后，能量、力和 σ 会共同决定晶胞是否继续更新。力阈值约束原子坐标，应力约束晶胞自由度。固定层、二维真空方向和约束晶格常数的模型，要把约束自由度从收敛判读中剥离；比较不同模型时，保留同一约束方案。

收敛阈值的数值必须跟后续任务相称。若后面要比较 0.01 eV 量级的吸附构型差异，电子噪声、离子步能量变化和 |F_max| 都应低于这个能量尺度对应的结构扰动；若只是做粗筛，过紧阈值会把计算资源消耗在不影响排序的微小变化上。

三类信号的用途分开：电子步收敛 看密度残差、占据和电子能量变化；能量收敛 看目标能量差和离子步势能下降；力收敛 看非冻结原子的最大残余力与晶胞应力。过渡态、声子、AIMD 和小能量差比较对初始结构更敏感，|F_max|、低频模式方向和电子占据稳定程度会进入频率、能垒或有限温轨迹。