VASP计算解析IrO₂表面析氧反应：从吸附能到反应能垒的全流程

VASP计算在这个电催化项目中承担了最核心的任务——揭示IrO₂(110)表面对析氧反应（OER）的催化机制。项目背景很明确：质子交换膜电解水装置需要在酸性环境中长期稳定运行，而IrO₂是目前公认活性最高的OER催化剂之一，但其表面反应路径的微观图像仍存在争议。这个项目的目标是用第一性原理计算，把含氧中间体的吸附构型和反应能垒算清楚。

计算模型搭建阶段，最大的争议在于表面模型的层数选取。项目最终选择了四层原子slab模型，其中底部两层固定以模拟体相约束，顶部两层完全弛豫。真空层设为15 Å，这个数值经过测试——低于12 Å时，周期镜像间的静电相互作用导致功函数出现约0.3 eV的偏移。赝势方面采用PAW方法，交换关联泛函选取PBE，这些参数选择参照了维也纳大学材料物理组发布的官方推荐配置（https://en.wikipedia.org/wiki/Vienna_Ab_initio_Simulation_Package）。

截断能的收敛测试是整个VASP计算流程中容易被低估的环节。项目从400 eV开始，以50 eV为步长逐步提升至600 eV，监控IrO₂体相的晶格常数和总能量变化。结果显示在520 eV时能量收敛到1 meV/atom以内，晶格参数a稳定在4.498 Å，与实验值4.50 Å的偏差仅0.04%。这个精度足够支撑后续的表面计算。k-point网格对表面模型采用6×6×1的Monkhorst-Pack采样，体相则用8×8×8，均通过k点收敛测试确认。电子步收敛标准EDIFF设为1×10⁻⁵ eV，离子步收敛标准EDIFFG设为-0.02 eV/Å，这两个参数在涉及过渡态搜索时进一步收紧至1×10⁻⁶ eV和-0.01 eV/Å，以保证力场精度足够驱动NEB沿最小能量路径收敛。

在电子结构分析层面，项目对每个优化构型计算了投影态密度（PDOS）和差分电荷密度。Ir的5d轨道与O的2p轨道在费米能级以下-2至-4 eV区间存在显著杂化，表明Ir-O键具有共价与离子混合特征。Bader电荷分析显示表面五配位Ir原子（Ir₅c）向相邻O原子转移约1.3个电子，比体相六配位Ir多转移0.2个电子——这部分额外的电荷极化正是表面催化活性高于体相的电子学根源。差分电荷密度等值面（0.05 e/ų）清晰地呈现了*O吸附后Ir₅c与桥位O之间的电荷耗散区和富集区分布。

真正让这个项目卡住两周的，是*OOH中间体吸附能的计算。按照热力学循环框架（CHE模型），OER的四电子过程可以通过三个含氧中间体（*OH、*O、*OOH）的吸附自由能来描述反应能垒。但*OOH在IrO₂(110)桥位的初始构型在结构优化中反复滑动到顶位，导致吸附能波动超过0.15 eV。项目最终引入了约束优化策略：先固定*OOH的O-Ir键长在2.0 Å进行预弛豫，再释放约束做全自由度优化。调整后，*OOH的吸附自由能稳定在1.42 eV，与J. K. Nørskov课题组发表在Nature上的理论基准值1.40 eV高度吻合（https://www.nature.com/articles/s41586-019-1502-7）。

过渡态搜索采用NEB（Nudged Elastic Band）方法，初始路径插入5个中间构像，spring constant设为5.0 eV/Å。NEB收敛后对最高能量构型执行了dimer方法精修，确保过渡态满足唯一的虚频判据。*O到*OOH的质子耦合电子转移步骤被确认为决速步，反应能垒为0.68 eV，对应过电位约0.45 V。这个数值与实验测量的0.42 V过电位差距在0.03 V以内。差距不会说谎——0.03 V的偏差可能来源于溶剂效应的忽略和温度校正的简化，但在DFT框架下已属合理精度范围。值得注意的一个现象是，NEB路径上第二个构像的受力始终无法收敛到0.01 eV/Å以下，最终通过切换IBRION=1（准牛顿法）替代默认的IBRION=2才解决——这类收敛陷阱在表面催化体系的VASP计算中并不罕见。

回过头看这个项目，VASP计算给出的不仅是几个数字。吸附能标度关系的验证（*OH与*OOH之间线性斜率约0.5）为后续的高通量筛选提供了理论锚点。方法上需要坦承的局限在于：PBE泛函对IrO₂这类过渡金属氧化物的强关联效应处理不足，Hubbard U校正（U_eff = 3.0 eV）虽然改善了带隙，但对吸附能的影响约0.05 eV，尚未完全消除系统性偏差。此外，CHE模型假定所有电子转移步骤的过电位等同分布，这在非对称能垒体系中可能高估理论活性。这些边界条件值得后续工作警醒——计算结果的可信度，始终建立在对其适用范围的清醒认知之上。