VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则

说明:本文华算科技主要介绍 VASP中 EDIFF 与 EDIFFG 的计算对象、停止条件、正负号含义,以及它们在静态计算、结构优化和声子计算中的阈值匹配。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则

VASP 两类收敛标准如何区分

VASP 做结构优化时,电子和离子并不在同一层循环里求解。固定一组原子坐标后,程序先求这个几何上的电子基态;电子密度、有效势、能带占据和总能量经过多轮自洽迭代,达到设定标准后,才用得到的能量、力和应力推动下一次几何更新。EDIFF 控制的就是这层电子自洽内循环,单位是 eV,它判断当前几何上的电子态是否已经稳定到指定数值尺度。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图1. 固定点电子自洽流程把输入密度迭代到收敛密度。

图1 画的是固定点电子自洽的通用逻辑:输入密度生成 Kohn-Sham 势,求解后得到新密度,再比较前后密度是否一致。VASP 的 EDIFF 不直接判断原子该不该移动,也不检查晶胞是否还受应力驱动;它只回答一件事:在这一个几何构型上,电子步能否停下来力和应力属于电子步之后的几何判断,离子步只能使用已经求出的电子态。

EDIFFG 控制的是离子弛豫外循环。当某个几何上的电子态已经求完,程序根据力、应力和能量变化决定下一步坐标或晶胞是否更新。EDIFFG 写成正数时,比较相邻离子步的总能变化;写成负数时,比较所有原子的力范数是否低于 |EDIFFG|;设为 0 时,离子弛豫按 NSW 给出的步数结束。

两个标签常出现在同一个 INCAR 文件里,却约束不同误差。EDIFF 管固定几何下的电子求解,EDIFFG 管几何搜索何时停止。电子步没有足够收敛,力会带着电子误差进入离子步;离子步阈值太松,电子态可以很稳定,但原子仍可能处在受力较大的位置。一个看电子变量,一个看几何变量,这是区分二者的起点。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则

EDIFF怎样影响电子步和力的质量

电子自洽不是把能量算一次就结束,而是不断缩小输入密度和输出密度之间的差别。SCF 论文中的铝晶胞例子显示,残差范数可随迭代步跨越多个数量级。VASP 使用 EDIFF 作为电子步停止条件时,真正要避免的是在电子密度还明显变化时就把力交给离子优化器。力来自电子基态能量对原子坐标的导数,电子态没有稳住,力的方向和大小都会夹杂数值噪声。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图2. 铝晶胞电子自洽残差随迭代步下降。

图2 的残差下降曲线可以帮助理解 EDIFF 的位置:它对应电子步内部的收敛判断,而不是优化后结构的好坏判断。对普通结构预优化,EDIFF = 1E-5 到 1E-6 eV 常见;对有限位移声子、微小能量差和精细力常数,电子步往往要更紧。力阈值越严,电子步误差越不能接近目标力尺度,否则 EDIFFG 看见的就不是干净的几何梯度。

EDIFF 不是唯一精度来源。平面波截断能、k 点、赝势、展宽方式和磁态都会改变电子能量基准。图3 中,PTO 结构的电子误差和能量差随平面波截断能改变,即使某次电子自洽达到了阈值,不同基组或截断条件仍会给出不同的误差水平。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图3. PTO 结构中电子误差和能量差随平面波截断能变化。

这类图提醒我们,EDIFF 只负责把某一套计算条件下的电子迭代收紧。若比较两个吸附构型、两个缺陷电荷态或两个磁态,能量差只有几 meV,电子步阈值要低于判别尺度,k 点和 ENCUT 也要同步检查。把 EDIFF 调小不会自动修复模型、基组或电子态选择问题;它只是把当前求解器在当前条件下的自洽误差压低。

金属、表面吸附、缺陷和含局域 d/f 电子的体系还会受占据和初始磁矩影响。费米能级附近态密度高时,电子步可能达到数值阈值,却收敛到另一个电子态附近。电子态本身仍要核对,磁矩、占据、DOS、局域电荷和结构响应都不能被 EDIFF 达标替代。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则

EDIFFG正负号为何会改变停止依据

在 INCAR 中,正值 EDIFFG 与负值 EDIFFG 不是同一单位下的松紧关系:正值比较相邻离子步总能变化,负值比较残余力;这两个量的单位和物理对象不同。结构优化更常用负值,例如 EDIFFG = -0.02 表示力阈值 0.02 eV Å-1,并不是能量 0.02 eV。负号不是数学装饰,它把判据从能量差切换到力,读 INCAR 时要先看符号。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图4. PTO 与 PZO 软模幅度改变时的能量路径。

图4 展示的是软模幅度变化时的能量路径。能量曲线靠近极小区会变平,离子步之间的能量差可能已经很小,但这并不必然说明每个自由度上的力都低于阈值。能量差变小不等于力已经达标;需要可靠几何基准时,负值力阈值更贴近结构优化的目的。

几何优化算法论文通常不会把图题写成 EDIFFG,但会在方法段给出梯度或最大力阈值。图5 对比了气相分子在 PM6、DFT 和 MP2 势能面上的优化收益;论文图注中对应的收敛标准为梯度范数阈值,其中 DFT 和 MP2 使用 10-3 eV Å-1。这正是 VASP 中负值 EDIFFG 所属的判断语境:看几何梯度是否足够小。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图5. 气相分子几何优化中预条件方法相对普通 LBFGS 的计算节省。

把这种梯度阈值放回 VASP,IBRION 决定怎样搜索下一步,ISIF 决定是否让晶胞自由度进入优化,选择性固定原子会改变可移动自由度集合。EDIFFG 只检查当前自由度集合上的停止条件;固定层、约束晶胞、表面偶极和应力处理都会改变收敛后结构的读法。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图6. 分子晶体 DFT 几何优化中不同预条件策略的计算节省。

图6 进一步把对象换成分子晶体 DFT 优化。不同预条件策略达到同一梯度阈值所需调用次数不同,说明离子步是否顺利结束不只由阈值数值决定,还受初始结构、势能面形状和优化算法影响。若电子步反复达到 NELM 上限,降低 EDIFFG 只会让几何条件更严,却不会解决电子态求解;若最大力长期高于 |EDIFFG|,单纯收紧 EDIFF 也不能让原子自动停在局部极小点。

AIMD 按时间步推进离子运动,关注温度、轨迹、能量漂移和结构演化,EDIFFG 不再作为离子弛豫停止条件。EDIFFG 属于几何最小化语境,不应拿它解释有限温轨迹什么时候结束。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则

不同计算任务中两个阈值怎样匹配

静态 DOS、能带、结构优化、声子和过渡态频率读取的目标量不同,阈值风险也不同。静态 DOS 或能带任务主要依赖固定结构上的电子态,EDIFF 要匹配带边、费米能级附近占据和总能差尺度;结构优化会把电子误差传给力,EDIFF 和 EDIFFG 要配合;声子、弹性常数和过渡态频率从力和二阶响应中取信息,电子步与几何步都要更谨慎。

VASP收敛参数辨析:EDIFF 与 EDIFFG 物理对象、收敛条件及阈值适配规则
图7. 二维超导材料高通量 DFT 筛选工作流包含能隙、磁矩、态密度、动力学稳定性和后续收敛步骤。

图7 的高通量 DFT 工作流把能隙、磁矩、费米能级附近 DOS、动力学稳定性和后续收敛步骤拆成不同筛选环节。这样的流程图说明,材料计算很少只靠一个阈值控制所有结果。阈值选择要服从目标输出量的判别尺度:粗优化用于把结构从高能初态带到合理区域,精优化用于压低最大力,静态计算接在精优化结构之后,并保持展宽、磁态、电荷态和电子阈值一致。

对普通结构预优化,EDIFF = 1E-5 到 1E-6 eV、EDIFFG = -0.03 到 -0.02 eV Å-1 常用于获得可继续计算的几何;对吸附能、缺陷形成能或相稳定性排序,若能量差只有 meV 级,电子步、k 点、截断能、磁态和电荷态都要低于目标差值尺度。末步最大力和电子步历史要随结果一起保存,否则后续很难区分物理差别与数值误差。

阅读 OUTCAR 或 OSZICAR 时,可以把检查对象分开:电子步看能量变化、占据和迭代历史,离子步看最大力、应力、位移和相邻离子步能量变化。声子或 NEB 任务还要保留前序优化结构的力、应力和电子步历史。先分清停止的是哪一层循环,再判断阈值是否与研究对象相称。

EDIFF 和 EDIFFG 的差别可以压缩成一句话:EDIFF 问“当前几何上的电子态是否收敛”,EDIFFG 问“几何搜索是否该停”。前者影响能量、力和应力的电子求解质量,后者决定结构优化终点;把二者混读,最常见的后果就是电子步看似达标,几何还没有到目标力阈值,或几何停下了,却没有核对电子态是否是研究需要的解。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 2025年7月21日 下午4:17
下一篇 2026年7月20日 下午2:38

相关推荐