本研究的主要作者为Asnake Sahele Haile、Weldegebriel Yohannes和Yedilfana Setarge Mekonnen，分别隶属于Addis Ababa University的环境科学中心与自然科学学院化学系。该研究于2020年7月21日发表在RSC Advances期刊上，论文标题为《Oxygen reduction reaction on Pt-skin Pt₃V(111) fuel cell cathode: a density functional theory study》。
一、学术背景与研究目的
质子交换膜燃料电池（proton exchange membrane fuel cells, PEMFCs）因其高效、清洁的能量转换特性而被视为极具前景的能源技术。然而，阴极氧还原反应（oxygen reduction reaction, ORR）动力学迟缓，且严重依赖昂贵的铂（Pt）基催化剂，这成为制约PEMFCs商业化应用的主要瓶颈。为降低Pt负载量并提升催化活性，研究者广泛探索将Pt与非贵金属过渡金属（non-precious transition metals, NPTMs）合金化的策略。在众多Pt基合金体系中，Pt₃V因其优异的电子效应和催化性能而受到关注，但此前对Pt-skin Pt₃V(111)表面ORR机理的理论计算研究尚属空白。基于此，作者采用密度泛函理论（density functional theory, DFT）系统研究了Pt-skin Pt₃V(111)电催化剂上ORR的催化活性、反应机理及稳定性，旨在为开发低成本、高效率的PEMFC阴极催化剂提供理论依据。
二、研究流程与方法
本研究采用基于Vienna Ab initio Simulation Package（VASP）的自旋极化DFT计算，并结合原子模拟环境（Atomic Simulation Environment, ASE）进行结构优化与能量计算。计算采用投影缀加波（projector augmented wave, PAW）方法描述核电子，平面波截断能设为400 eV，使用Fermi展宽（宽度0.1 eV）处理电子占据，布里渊区采样采用5×5×1 Monkhorst-Pack k点网格，交换关联泛函选用revised Perdew–Burke–Ernzerhof（RPBE）。所有结构优化至每个离子上的Hellmann-Feynman力小于0.01 eV/Å。
Pt-skin Pt₃V(111)表面模型采用最外层为纯Pt原子、内层为L1₂有序Pt₃V晶体结构构建。计算得到Pt的晶格常数为3.9942 Å，Pt₃V为3.93 Å，均与已有文献吻合。表面模型为2×2超胞，含四层原子，底部两层固定，上部两层弛豫。为模拟溶剂环境，采用VASPsol隐式溶剂化模型，介电常数设为78.54。文章定义了O、OH、OOH等ORR中间体的吸附能（adsorption energy, ΔE_ads）计算公式，并利用CI-NEB方法定位各基元反应的过渡态，进而计算反应能（reaction energy, E_r）与活化能（activation energy, E_a）。此外，通过形成能（formation energy, ΔE_alloy）和电极电位偏移（dissolution potential shift, ΔU）评估Pt₃V合金的热力学稳定性及表面Pt原子的抗溶解能力。
三、主要研究结果
在表面吸附位点分析中，由于次表层V原子的存在，Pt-skin Pt₃V(111)表面形成了f1、f2、t1、t2、b1、b2、b3等多种非等价吸附位点。计算表明，O的最稳定吸附位为f1，OOH和OH的最稳定吸附位均为t1，其吸附能分别为2.06 eV、3.93 eV和0.75 eV。与纯Pt(111)表面相比，Pt-skin Pt₃V(111)上三种ORR中间体的吸附能均较低，表明中间体与表面之间的相互作用更弱。d带中心分析显示，Pt-skin Pt₃V(111)的d带中心相对于纯Pt(111)发生下移，远离Fermi能级，这是导致吸附减弱的主要原因。Bader电荷分析进一步表明，每个Pt原子平均从V原子获得0.41 e⁻的电荷，使Pt-skin表面带更多负电，增强了与含氧中间体之间的静电排斥，从而削弱了吸附。为区分应变效应（strain effect）与配体效应（ligand effect），作者构建了与Pt₃V具有相同晶格常数的应变Pt(111)表面，结果表明对于OOH和O，其弱吸附由应变与配体效应共同贡献；对于OH，应变效应占主导。
在ORR基元反应研究中，作者考察了O₂解离、OOH形成、OOH解离、OH形成、H₂O₂形成、H₂O₂解离及H₂O形成等步骤。计算结果显示，O₂解离的活化能高达1.00 eV，而OOH形成的活化能仅为0.36 eV，因此O₂优先通过加氢生成OOH而非直接解离。OOH随后解离为O和OH的活化能为0.33 eV，低于H₂O₂解离路径中H₂O₂解离为两个OH的活化能（0.20 eV），但OOH解离后的O加氢生成OH的活化能为0.64 eV，成为整个ORR过程的速率决定步骤（rate-determining step）。相比之下，纯Pt(111)表面相同步骤的活化能为0.72 eV，表明Pt-skin Pt₃V(111)能显著降低ORR的能垒，从而提升催化活性。综合活化能分析，Pt-skin Pt₃V(111)上的ORR遵循OOH解离机理，即O₂吸附→OOH形成→OOH解离→O加氢生成OH→OH加氢生成H₂O。
自由能图进一步表明，在零电极电位下，Pt-skin Pt₃V(111)上ORR各基元步骤均为放热反应；当电位升至平衡电位1.23 V时，O加氢步骤变为吸热。计算得到的极限电位为0.59 V，过电位为0.64 V，略高于纯Pt(111)表面，但结合活化能结果，Pt-skin Pt₃V(111)仍表现出更优的催化活性趋势。
在稳定性评估中，Pt₃V合金的形成能为−0.36 eV/atom，负值表明该合金具有热力学稳定性，与已有Pt₃V和Pd₃V的计算结果一致。通过构建表面Pt原子移除模型，计算得到清洁Pt-skin Pt₃V(111)表面相对于纯Pt(111)的电极电位偏移ΔU为0.055 V，正值说明Pt原子从合金表面溶解需要更高的电极电位，即抗溶解能力增强。在吸附0.25单层原子氧的条件下，Pt-skin Pt₃V(111)与纯Pt(111)相比的ΔU为0.065 V，进一步证实即使在氧化性电化学环境下，Pt-skin Pt₃V(111)仍具有优于纯Pt的抗Pt溶解能力。
四、结论与意义
本研究通过DFT计算系统揭示了Pt-skin Pt₃V(111)电催化剂上ORR的反应机理与稳定性特征。结果表明，Pt-skin Pt₃V(111)对ORR中间体具有适中的吸附强度，ORR通过OOH解离路径进行，速率决定步骤为O加氢生成OH，其活化能为0.64 eV，低于纯Pt(111)的0.72 eV，表明该合金催化剂具有更高的ORR催化活性。同时，负的合金形成能和正的电极电位偏移证实了Pt₃V合金及其Pt-skin表面在热力学和电化学条件下均具有良好的稳定性。该研究从理论上证明了Pt₃V合金在降低Pt用量、提升ORR活性与稳定性方面的潜力，为开发低成本高效PEMFC阴极催化剂提供了重要的理论依据和设计思路。
五、研究亮点
本研究的亮点在于：首次通过DFT计算系统阐明了Pt-skin Pt₃V(111)表面上ORR的完整反应路径，明确了OOH解离机理及速率决定步骤；通过分离应变效应与配体效应，揭示了Pt₃V合金弱化中间体吸附的物理机制；同时结合合金形成能与电极电位偏移双重指标，全面评估了催化剂的热力学与电化学稳定性，为Pt基合金催化剂的设计提供了方法学参考。