截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析

说明:本文华算科技主要介绍平面波截断能的计算含义、赝势来源、目标量误差和收敛判断标准。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析

截断能控制的是平面波基组上限

平面波DFT 把 Kohn-Sham 轨道写成一组倒空间平面波的线性组合。截断能 Ecut 规定了可进入基组的最大动能,满足 |G+k|2/2 小于阈值的平面波会保留下来。截断能升高时会加入更短波长的基函数,电子密度在原子附近的快速变化才有机会被描述。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
图1. C 和 Fe2+ 的 PA 与 Troullier-Martins 赝势径向形状对比。DOI:10.48550/arXiv.1803.08180

截断能没有跨体系通用的“多少 eV”。同样 400 eV,对软赝势可能已经足够,对一套硬的规范守恒赝势、含半芯态 PAW、第一排元素、氧化物或高压结构,误差仍可能很大。赝势的平滑半径、价电子选择和角动量分量会直接改变所需平面波数量,元素表里的推荐值只能给出初始范围。

图1中的 C 和 Fe2+ 赝势径向形状说明,芯区势函数越尖锐,波函数在实空间变化越快,倒空间就要包含更高 G 分量。平面波数量大致随 Ecut3/2 增长,把截断能调高一档会显著增加计算量;收敛测试追求目标误差内的最小有效值,而不是单纯取很大的数。

平面波基组对周期体系很干净,但对原子附近的尖锐波函数不友好。赝势把芯区波函数平滑后,截断能才成为可控的基组上限。若改用全电子方法、局域轨道基组或实空间网格,数值参数会换成 RKmax、轨道半径或网格间距,平面波 Ecut 的经验值不再适用。

赝势文件中的推荐截断能只是起点。生产计算前,推荐值通常要乘以 1.2 到 1.5 做检查;混合元素体系还要按最硬元素设定。同一批数据必须使用同一套赝势和同一截断能口径,这样形成能、吸附能或相对能量才不会混入基组差异。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析

赝势推荐值为什么只是起点

赝势生成时会测试原子散射态、对数导数、固体晶格常数、体模量、声子或能带。图2的 Fe2+ 对数导数把 s、p、d、f 角动量分量分开比较,赝势可靠性既来自总能曲线平滑,也来自 不同角动量分量的散射性质

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
图2. Fe2+ 在不同角动量分量中的对数导数随能量变化。DOI:10.48550/arXiv.1803.08180

在材料计算里,截断能测试应固定赝势、泛函、k 点、展宽、结构和电子收敛阈值,只改变 Ecut。若每个截断能都重新优化晶格常数,能量变化会混入结构差异;若 k 点同步改变,误差来源会合并。单变量扫描才能判断平面波基组本身是否足够。

测试点不必无限密集。常见做法是从推荐值附近开始,例如 400、450、500、550、600 eV,或者按 Ry 取等间隔序列。若相邻两个高截断能点的目标量变化已经小于设定阈值,再把最高点作为近似参考。参考点要高于工作点,否则低截断能下的平滑曲线可能只是尚未暴露误差。

在含 O、F、N、C、第一过渡金属或半芯态元素的混合体系中,硬元素往往决定基组上限。若只按金属主体的推荐值取 Ecut,O-H 键、M-O 键长、应力和振动频率可能残留可见误差。截断能收敛的对象是整个体系的最硬局部电子结构,不是化学式里占比最高的元素。

VASP 常用 eV,Quantum ESPRESSO 常用 Ry;超软赝势和 PAW 还会有波函数截断与电荷密度截断。跨代码比较时,不同单位、不同赝势库和双截断设置要先换成各自目标量的残余误差,单独比较 400 eV、60 Ry 或 80 Ry 没有材料意义。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析

总能收敛不一定覆盖力、压力和声子

很多人测试截断能只看总能随 Ecut 的变化是否小于 1 meV/atom。这个标准适合筛查相对能量,但结构优化还依赖原子力,表面能和弹性常数还依赖应力,声子计算还依赖能量的二阶导数。导数性质通常比总能更敏感,它们要求的截断能可能更高。

图3给出金刚石 C 的四类量:总能 E、压力 p、原子力 F 和 Γ 点光学声子频率 ω。总能曲线较早趋于平滑,压力、力和声子频率的残余差异仍然存在。同一个截断能对不同目标量对应不同误差,用于结构、声子或应力的计算,不能只沿用能量差的测试结果。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
图3. 金刚石 C 中总能、压力、原子力和光学声子频率随截断能收敛。DOI:10.48550/arXiv.1803.08180

判断收敛时可以把目标量写成 ΔX(Ecut) = X(Ecut) – X(Ecutmax)。若目标是吸附能,比较的是 ΔEads;若目标是相变压力,比较的是压力或焓差;若目标是反应能垒,初态、过渡态和终态要在同一截断能下计算。判断标准跟研究量绑定,固定 eV 数字只负责生成测试点。

Pulay stress 来自平面波基组随晶胞体积变化而改变。低截断能下能量-体积曲线会带入基组误差,体模量、表面应力和弹性常数读数都可能偏移。含应力的计算通常要比单点能量使用更高截断能,尤其是变胞优化、弹性常数和高压相变计算。

绝对总能随截断能变化很大并不奇怪,真正进入材料结论的是能量差。两个结构若化学组成相同、赝势相同、截断能相同,许多基组误差会部分抵消;缺陷形成能、吸附能和反应能垒仍会因局部键型变化留下残余项。相对能量收敛要按同一反应式或同一比较组检查,孤立晶胞的总能曲线只提供基础参照。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析

截断能多大才算收敛取决于目标误差

在 bcc Fe 中,d 电子、自旋极化和赝势硬度都会改变截断能需求。图4中 Fe 的总能、压力、力和光学声子频率随 Ecut 变化,PA 与 Troullier-Martins 赝势的收敛速度并不相同。含 d 电子、自旋极化和较硬赝势时,推荐截断能附近仍可能留下力和声子误差,这类误差会进入磁相能量差和结构稳定判断。

截断能多大才算收敛?平面波截断能物理意义、赝势影响、计算误差与收敛标准解析
图4. bcc Fe 中总能、压力、原子力和光学声子频率随截断能收敛。DOI:10.48550/arXiv.1803.08180

实际判断可把目标误差先定出来:相对能量常用 meV/atom 或 meV/吸附物,力可用 eV/Å,压力可用 GPa,声子可用 cm-1。若 Ecut 从 500 eV 到 550 eV 时,吸附能变化只有 2 meV,而研究结论只比较 0.1 eV 量级差异,500 eV 就可作为该任务的工作值。收敛是误差相对结论尺度足够小

若两个候选结构的能量差只有 5 meV/atom,截断能测试残余误差达到 8 meV/atom,排序就没有数值支撑;若能量差是 200 meV/atom,残余误差 5 meV/atom 通常不会改变主结论。阈值来自待比较能量差的尺度,同一套参数在粗筛和精排中可能对应不同要求。

截断能也要服务于数据一致性。高通量筛选、相图、缺陷形成能、吸附能矩阵和 NEB 路径中,多个结构的能量要放在同一张比较表里,截断能最好统一到最硬体系通过测试的值。同一数据集的可比性会影响排序、筛选阈值和后续复算优先级。

论文或报告里适合写目标量在相邻更高截断能下的变化幅度,而不是只写“ENCUT=520 eV 已收敛”。例如“Ecut 从 520 eV 增至 580 eV 时,关键吸附能变化小于 5 meV,最大残余力变化小于 0.01 eV/Å”。这样的表述把 截断能、目标量和误差阈值 固定在同一句话里,也给后续复算留下清楚参照。

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

相关推荐