VASP计算吸附能时气相分子盒子多大才合适？

引言

吸附能和反应自由能里常有一个被忽略的参照：孤立气相分子。O₂、H₂O、CO₂ 和 NH₃ 看起来只是小分子，真正进入 VASP 时却会牵涉盒子尺寸、K 点、自旋态、偶极相互作用和热力学校正。气相参照不是随手跑一个单点能，而是决定整组吸附能能否比较的能量基准。

气相盒子先控制什么问题？

孤立分子计算要模拟的是一个不与周期镜像相互作用的分子。对 O₂、CO₂、H₂O、NH₃ 这类小分子，最基本的做法是把分子放在足够大的立方或近似立方盒子中心，只使用 Γ 点，并让分子在盒子内完成结构优化。盒子太小会让分子和自己的周期镜像发生静电或轨道尾部相互作用，最后进入 E_ads 时表现为系统性偏移。

实际建模时，可以先把 O₂、H₂O、CO₂ 和 NH₃ 分别放入 15 Å 左右的孤立分子盒子，再用 18 Å 或 20 Å 盒子核对总能变化。若继续增大盒子后能量差已经低于吸附能比较所需精度，就可以把该尺寸作为本组计算的统一口径。盒子尺寸收敛的对象不是绝对总能大小，而是气相分子能量随盒子继续增大时的变化量。

图1 N 掺杂聚三嗪模型中 CO₂ 吸附参照的结构示意。图中同时展示了平面周期模型和分子吸附前的几何尺度。DOI: 10.1021/acsomega.3c00395。

图1 的平面模型提醒我们，吸附计算里的分子参照和表面模型处在两套边界条件中。气相分子用大盒子和 Γ 点近似孤立态，表面吸附物则处在 slab 的周期场中。二者要通过同一套 POTCAR、相近 ENCUT 和明确的能量口径联系起来，不能在一组材料里混用不同赝势或不同泛函。

如果分子带有明显偶极，例如 H₂O 和 NH₃，盒子收敛比 CO₂ 更值得单独检查。对中性小分子，通常不需要像 slab 那样引入表面偶极修正，但分子取向改变时，有限盒子内的偶极镜像作用可能让最后几 meV 不易稳定。实际处理中，把分子放在盒子中心、关闭不必要的对称性，并用同一取向规则生成对比组，会让数据更容易追踪。

O₂ 自旋为什么要单独处理？

O₂ 是最容易让气相参照出错的小分子。它的基态是三重态，不能按闭壳层分子直接用非自旋极化计算。VASP 中通常应设置 ISPIN = 2，并给 O 原子合理的初始磁矩；若计算收敛到低磁矩或错误占据态，后续 ORR、OER 或氧空位形成能都会带入偏差。O₂ 的参照能首先要确认自旋态，而不是只确认几何优化是否结束。

O₂ 还会出现在表面吸附、O₂ 活化和 H₂O₂ 路径中。图2 展示了双原子位点上 O₂ 相关轨道和电荷转移的判读方式：吸附态中 Fe–O 键、O₂ 反键轨道和电子转移共同决定活化程度。这个图不能替代气相 O₂ 计算，却能说明为什么气相 O₂ 的自旋参照必须先算对；只有参照态清楚，吸附后轨道变化才有比较意义。

图2 Fe₂N₆ 位点上 O₂ 与 OH 吸附的电荷转移和 Fe-O 轨道分析。图中给出 O₂ 活化与 OH 脱附相关的轨道占据证据。DOI: 10.1039/d5sc04842g。

O₂ 气相参照的 INCAR 参考示例

SYSTEM = O₂_gas
ISPIN = 2
MAGMOM = 1.0 1.0
ISMEAR = 0
SIGMA = 0.05
EDIFF = 1E-6

这组参数只用于说明 O₂ 气相参照的基本口径。实际使用时，MAGMOM 初值应与 POSCAR 中两个 O 原子的顺序对应，KPOINTS 采用 Γ 点，盒子尺寸仍需按前文方法做收敛测试。真正决定参照是否可用的是 O-O 键长、总磁矩和总能随盒长变化同时满足要求。若后续 slab 采用 DFT+U、HSE06 或特定 vdW 修正，气相参照也要用同一计算口径重算。

很多电催化自由能会绕开直接使用 DFT O₂ 总能，改用 H₂O/H₂ 或实验热力学关系校正 O₂ 误差。即便如此，O₂ 吸附和 O₂ 活化本身仍需要自旋态明确的气相或弱相互作用参照。是否直接把 E_O₂ 代入公式，是自由能口径问题；O₂ 是否按正确自旋态计算，是电子结构问题。

CO₂、H₂O 和 NH₃ 的构型怎样影响参照？

CO₂ 的气相结构应接近线性，吸附后常发生弯曲或局域电荷转移。图3 显示金属掺杂 NPT 模型上 CO₂ 的吸附构型差异，说明 CO₂ 参照态和吸附态不能混为一谈。气相 CO₂ 能量对应线性中性分子，吸附态能量则包含分子弯曲、表面极化和化学键形成的贡献。

图3 不同金属掺杂 NPT 模型上的 CO₂ 吸附构型。图中展示了 CO₂ 在不同位点上的几何取向和局域相互作用差异。DOI: 10.1021/acsomega.3c00395。

H₂O 的参照要兼顾几何和热力学。VASP 优化给出的是 0 K 电子总能，电化学自由能表里常需要再加入 ZPE、熵和溶剂相关修正。图4 中 H₂O 以不同姿态接近氧化物表面，吸附后可以保持分子态，也可以进入解离路径。若只用一个 E_H₂O(g) 参照解释所有水相关步骤，必须说明它对应的是气相水分子，而不是界面水网络。

图4 H₂O 分子接近 PuO₂(111) 表面的不同姿态。图中给出了水分子取向与表面位点之间的几何关系。DOI: 10.1039/d4ra01015a。

NH₃ 和 H₂O 类似，都是闭壳层极性分子，但 NH₃ 的 N 端配位、H 端氢键和逐步脱氢路径会引入不同吸附态。图5 展示 NH₃ 脱氢相关能量路径，说明同一个气相 NH₃ 参照会连接 NH₃*、NH₂*、NH* 等不同中间体。气相参照只提供起点能量，不能替代对吸附取向、解离端点和反应路径的构型筛选。

图5 NH₃ 在氧化物催化剂上的逐步脱氢能量路径。图中给出 NH₃、NH₂、NH 和 N 相关态的结构与能量变化。DOI: 10.3390/nano8110958。

吸附能公式里哪些量要保持同口径？

吸附能常写为：
式(1) E_ads = E_slab+X – E_slab – E_X(gas)。
其中 X 可以是 O₂、H₂O、CO₂、NH₃ 或其他分子。这个公式看似简单，真正需要统一的是三个计算的泛函、POTCAR、ENCUT、自旋设置、能量收敛和色散修正。slab 与 slab+X 还应使用相同晶胞、相同 K 点密度和相同固定层规则。

若进一步转为自由能，可写为：
式(2) ΔG_ads = E_ads + ΔZPE – TΔS + ΔG_solv + ΔG_U,pH。
式(2) 中后几项不是所有体系都必须加入，但只要加入，就应在同一批材料中使用同一套来源。O₂、H₂O、CO₂、NH₃ 的热力学校正不能随意从不同论文摘取后混合代入。

实际任务可以分成三步。第一步，分别建立 O₂、H₂O、CO₂、NH₃ 的孤立分子盒子，确认几何、自旋和盒子收敛。第二步，用表面 slab 放置多个吸附取向，保留同一 slab 参照和相同 K 点密度。第三步，把电子总能、ZPE、熵和电位或 pH 修正分列，不把修正项提前揉进总能。只有参照态、吸附态和修正项各自独立记录，吸附能差异才有可解释性。

频率计算也要和参照态一一对应。气相 CO₂、H₂O 和 NH₃ 的振动频率可以用于 ZPE 项，但平动和转动熵通常不直接由 VASP 普通频率输出完整给出，需要结合热化学后处理或采用可追溯的实验热力学数据。若一个自由能表只使用电子总能，则应明确它是 0 K 电子能比较；若加入熵项，就应说明温度、压力和数据来源。电子总能表和吉布斯自由能表不能在同一张图中混用同一个符号。

在同一材料上同时比较 CO₂ 产 CO 路径和 NH₃ 吸附时，气相参照的差异会直接进入能量差。若 CO₂ 用 20 Å 盒子、NH₃ 用 10 Å 盒子，而 O₂ 又采用非自旋极化参照，不同反应之间的能量差就会混入设置差异。反过来，四类分子都完成盒子收敛，O₂ 自旋和闭壳层分子口径清楚，后续的位点比较才主要反映化学相互作用。

本章要点总结

要点：气相分子盒子应通过总能随盒长变化来收敛，而不是只给一个固定尺寸。

要点：O₂ 必须检查自旋态；CO₂、H₂O 和 NH₃ 重点检查几何、偶极和热力学校正口径。

要点：E_ads 中的 slab、slab+X 和 X(gas) 要使用可比的泛函、赝势、ENCUT 和收敛标准。

要点：自由能修正应独立列出，不能把不同来源的 ZPE、熵和电位项混在一组比较中。