📚 课程目录第4章力与能量收敛标准

4. 力与能量收敛标准：如何设置合理的收敛判据

做结构优化，说白了就是在找那个「最舒服」的原子排布。但你怎么知道优化到头了？这就得靠收敛判据来拍板。我见过不少新手，上来就把收敛标准设得特别严，结果算了一个星期还在跑；也有人设得太松，算出来的结构根本不对。

今天咱们就把这事聊透。力、能量、位移，这三个判据到底怎么用，怎么搭配，我心里有个谱。

4.1 三大收敛判据：各司其职

先看一张总览图，把三个判据的关系理清楚。

4.2 能量判据：最直观，也最容易骗人

能量判据看的是相邻两步之间体系总能量的变化。如果 ΔE 小于你设的阈值，比如 1e-5 eV，优化就算收敛了。

听起来挺合理对吧？但我得提醒你，能量判据有个大坑——它太容易满足了。尤其是体系比较大的时候，总能量动辄几千 eV，相邻两步差个 1e-5 eV 根本不算啥。我曾经算一个 200 个原子的表面体系，能量判据三步就过了，但一看原子受力，还有 0.5 eV/Å 呢。

⚠ 注意： 能量判据不能单独使用。它只能告诉你「体系整体变化不大」，但局部可能还差得远。我建议把它当作辅助判据，配合力和位移一起用。

4.3 力判据：最可靠的主判据

力判据看的是所有原子中受力最大的那个。如果最大原子力 F_max 小于阈值，比如 0.01 eV/Å，说明每个原子都基本处于平衡位置了。

我个人习惯把力判据作为主判据。为什么？因为力直接反映了原子是否「想动」。能量可能骗人，但力不会。一个原子受力大，它就是没到位，没什么好商量的。

不同场景下，力判据的阈值怎么选？我整理了一张表：

应用场景	力阈值 (eV/Å)	说明
粗优化	0.05 ~ 0.1	快速扫描，初步筛选结构
标准优化	0.01 ~ 0.05	大多数体系，平衡精度和速度
精细优化	0.001 ~ 0.01	过渡态、反应路径、高精度计算
极端精度	< 0.001	声子谱、弹性常数、需要二阶导数的场景

💡 我的经验： 做表面吸附的时候，力判据设到 0.02 eV/Å 就够了。但做过渡态搜索，我一般会压到 0.005 eV/Å 以下。你想想看，过渡态对结构太敏感了，差一点能量就偏了。

4.4 位移判据：防止原子「荡秋千」

位移判据看的是原子位置的变化量。如果最大位移 ΔR 小于阈值，比如 0.001 Å，说明原子基本不动了。

这个判据在什么情况下特别重要？我遇到过一种情况：体系在能量和力上都收敛了，但原子还在小幅度振荡。尤其是软材料或者分子晶体，原子间的势能面很平缓，力很小但位移不小。这时候位移判据就能派上用场。

位移判据的典型设置：

刚性体系（金属、陶瓷）：0.001 ~ 0.005 Å
柔性体系（有机分子、聚合物）：0.005 ~ 0.02 Å
非常软（分子晶体、层状材料）：0.01 ~ 0.05 Å

4.5 如何搭配使用：我的推荐方案

三个判据各有侧重，怎么搭配最合理？我给出两套方案，你根据场景选。

方案一：标准模式（力 + 能量）

这是最常用的组合，适合大多数晶体和表面优化。

# VASP 示例
EDIFF = 1E-5        # 能量判据
EDIFFG = -0.02      # 力判据（负号表示力的收敛标准）

# Quantum ESPRESSO 示例
etot_conv_thr = 1.0D-5    # 能量判据
forc_conv_thr = 1.0D-3    # 力判据（单位 Ry/Bohr，约 0.026 eV/Å）

方案二：精细模式（力 + 能量 + 位移）

适合过渡态、表面重构、分子吸附等对结构敏感的场景。

# VASP 示例
EDIFF = 1E-6
EDIFFG = -0.005
# VASP 没有直接的位移判据，但可以通过 NSW 和 POTIM 控制

# CP2K 示例
MOTION % GEO_OPT
  MAX_ITER 200
  OPTIMIZER BFGS
  MAX_DR 0.003        # 位移判据 (Bohr)
  RMS_DR 0.0015       # 均方根位移
  MAX_FORCE 0.0005    # 力判据 (Hartree/Bohr)
  RMS_FORCE 0.0003    # 均方根力
END

  🔑 核心原则： 力判据是主，能量判据是辅，位移判据是保险。三者同时满足才算真正收敛。我建议至少同时使用力和能量两个判据。

4.6 避坑指南：我踩过的几个坑

做结构优化这么多年，有些坑我是一步一步踩过来的。分享给你，省得走弯路。

能量判据设得太松：我曾经为了省时间，把 EDIFF 设成 1E-4，结果优化出来的结构力还有 0.1 eV/Å。后来算性质全不对，白费了一个星期。
力判据设得太严：有一次做分子动力学后的弛豫，我把力判据设到 1E-5 eV/Å，跑了 500 步还没收敛。后来发现体系太大，根本没必要。
只看最大力，不看均方根力：最大力反映最差的原子，均方根力反映整体。如果最大力很小但均方根力还很大，说明大部分原子还没到位。我一般两个都看。
忽略位移判据：做有机分子的时候，力已经收敛了，但分子还在慢慢转。加上位移判据后，问题就解决了。

⚠ 特别提醒： 不同软件对收敛判据的定义可能不同。比如 VASP 的 EDIFFG 是负值表示力的标准，正值表示能量的标准。CP2K 的 MAX_FORCE 和 RMS_FORCE 是分开的。一定要看手册，别想当然。

4.7 如何判断收敛是否合理

判据设好了，优化跑完了，怎么知道结果靠谱不？我一般做三步检查：

看输出文件：确认最后几步的力、能量、位移都小于阈值，而且不是突然跳变的。
看结构：用可视化软件打开优化后的结构，看看原子位置是否合理。有没有原子靠得太近？有没有键长异常？
做单点能验证：在优化后的结构上做一个高精度的单点能计算，和优化最后一步的能量对比。如果差很多，说明优化可能没到位。

我记得有一次，优化后的结构看起来没问题，力也收敛了，但单点能比优化最后一步高了 0.1 eV。后来发现是优化过程中能量震荡，最后一步刚好落在了一个不太好的点上。重新优化后就好了。

💡 小技巧： 如果条件允许，可以用不同的初始结构做几次优化，看最终结构是否一致。这叫「收敛性测试」，是验证优化结果可靠性的好方法。

4.8 总结：收敛判据设置速查

最后，我把今天的内容浓缩成一张速查表，方便你以后用的时候参考。

场景	能量判据	力判据 (eV/Å)	位移判据 (Å)	备注
快速扫描	1E-4	0.05 ~ 0.1	0.01	粗筛，不要求高精度
标准晶体	1E-5	0.01 ~ 0.02	0.005	最常用组合
表面/界面	1E-5 ~ 1E-6	0.005 ~ 0.01	0.003	表面原子弛豫敏感
分子/团簇	1E-6	0.001 ~ 0.005	0.001	分子构型变化小
过渡态	1E-6 ~ 1E-7	< 0.005	0.001	精度要求最高

嗯，收敛判据这事，说复杂也复杂，说简单也简单。记住一句话：力判据是王道，能量判据做辅助，位移判据防振荡。按这个思路去设，基本不会出大问题。