Please wait a minute...
图/表 详细信息
Pt3Co高指数晶面氧还原过程的密度泛函理论研究
刘磊, 郭瑞华, 王丽, 王艳, 张国芳, 关丽丽
无机材料学报    2025, 40 (1): 39-46.   DOI:10.15541/jim20240296
摘要   (281 HTML13 PDF(pc) (8288KB)(216)  

Pt3Co催化剂是Pt基合金中氧还原反应(ORR)活性最高的催化剂, 合成Pt3Co高指数晶面(HIFs)是一种提高其催化性能的有效策略, 但拥有最高ORR活性的HIFs尚未明确, 并且目前缺乏对Pt3Co HIFs ORR的系统研究。本研究构建了六种不同Pt3Co HIFs, 通过从头算分子动力学(AIMD)计算证明了其稳定性, 通过密度泛函理论(DFT)计算了六种Pt3Co HIFs的ORR过程中间物*O、*OH、*OOH的结合能(BE), 通过d带中心(εd)、Bader电荷及配位数(CN)解释了其在台阶与边缘位点BE不同的原因。同时分析了吸附原子CN与εd的关系, 通过ORR自由能台阶图分析了ORR过程中的过电位(η), 发现η大小主要与*OH结合能(BE-*OH)有关, 其中η最小的晶面为Pt3Co(211), 其在台阶处的η达到了0.294 eV。本工作为高ORR活性HIFs催化剂研发提供了一定的理论依据。


HIFs Pt3Co(211) Pt3Co(310) Pt3Co(331) Pt3Co(511) Pt3Co(320) Pt3Co(332)
ηORR (edge)/eV 0.630 0.889 0.675 0.706 1.125 0.654
ηORR (terrace)/eV 0.294 0.556 0.556 0.631 0.874 0.507
View table in article
表3 六个Pt3Co HIFs在不同位置的ηORR
正文中引用本图/表的段落
Stamenkovic等[15]通过密度泛函理论(DFT)计算与实验得出在不同的Pt3M(M = Ni、Co、Fe、Ti、V)合金中, 催化活性与Pt表面d带中心(εd)位置之间的关系呈现火山形状, 其中Pt3Co的催化活性最高。并且在已研究过的Pt-M催化剂中, PtCo催化剂因其良好的活性和稳定性, 被广泛认为是PEMFC中最有前途的ORR催化剂[16?-18]。Tetteh等[19]制备了一种高耐久性、高活性的PtCo催化剂, DFT计算表示, 增强的ORR性能源于PtCo(221), 高分辨率透射电子显微镜也证实了催化剂中存在该晶面。Dong等[20]用原位拉曼光谱研究了Pt(211)和Pt(311)的ORR过程, 研究表明, 不同中间产物在不同表面上的特定吸附导致Pt(211)的ORR活性高于Pt(311)。Huang等[21]通过无配体固态反应, 使用微量元素锑、铋、铅或碲来稳定HIFs, 模拟和实验都证实这种方法可以稳定(210)晶面。
首先对Pt3Co原胞进行优化, 优化后, Pt3Co晶胞的晶格常数为0.383 nm, 与文献[44]相吻合(Pt3Co: 0.388 nm)。然后对Pt3Co晶胞进行切面、扩胞, 保持Pt、Co原子比为3:1, 再进行结构优化, 得到六个吸附表面模型, 如图2所示。为了评估其稳定性, 计算了六个表面的表面能。结果表明Pt3Co(211)、Pt3Co(310)、Pt3Co(331)、Pt3Co(511)、Pt3Co(320)、Pt3Co(332)的表面能分别为1.537、1.921、1.696、1.881、1.820、1.569 J/m2, 与Zeng等[45]的计算结果基本吻合。同时还对表面模型进行了AIMD计算, 结果如图S1所示。在T=300 K时, Pt3Co(211)、Pt3Co(310)、Pt3Co(331)、Pt3Co(511)、Pt3Co(320)、Pt3Co(332)的表面能量趋于稳定, 证明所构建的六个表面具有良好的热力学稳定性。
为了进一步研究Pt3Co HIFs催化剂的ORR活性, 计算了4电子反应过程中, 0和1.23 V以及理论反应电压下的反应自由能, 六个HIFs在边缘处和台阶处的ORR自由能台阶图如图S9、图5所示。可以发现在0 V下, 六个HIFs表面的ORR自由能均逐级降低, 表明ORR过程是自发进行的,Pt3Co(211)、Pt3Co(310)、Pt3Co(331)、Pt3Co(511)、Pt3Co(320)、Pt3Co(332) HIFs催化剂均具有良好的ORR催化活性。
Pt3Co(310)、Pt3Co(331)、Pt3Co(511)、Pt3Co(320)、Pt3Co(332)在边缘处和六个HIFs在台阶处的ORR步骤中, 当电压达到平衡电压(1.23 V)时, 反应(4)上升变化最大。这表明, *OH获得电子并与氢离子结合生成H2O, 最后从催化剂表面解吸的过程是影响催化剂活性的速率决定步骤(RDS), Pt3Co(211)在边缘处的RDS则为反应(2)。根据自由能台阶图, 可以得出在六个HIFs台阶和边缘处的ηORR, 如表3所示。同时, ηORR与BE-*OH有线性比例关系(图6), 其中Pt3Co(211)催化剂位于线的末端, 表明其ORR催化活性较好。相比之下, Pt3Co(320)位于线的顶端, 表明其催化活性较差。其中Pt3Co(211)的ηORR最小, 其在台阶处达到了0.294 eV, 这是因为Pt3Co(211) 对*O、*OH、*OOH的吸附较弱, 而较强的吸附会导致ORR过程中吸附的分子不易解离, 能垒上升。这表明Pt3Co(211)在热力学上能提供更充分、更快的ORR过程, 从而提高能量转换效率。Pt3Co(211)还可以显著降低催化剂表面的解吸能垒, 提高解吸速率, 从而进一步提高催化效率[29]。通过图6对比边缘与台阶处的ηORR, 可以发现每一个HIFs台阶处的ηORR都低于边缘处, 这可能是边缘原子对*O、*OH、*OOH的吸附过强[56], 导致边缘原子被吸附的*O、*OH、*OOH毒化。然而, 在HIFs中, 台阶处的原子数量远多于边缘处的原子数量, 可以为ORR提供更多的反应位点[46], 这也是总体上HIFs ORR活性较高的原因。可以注意到, 六种Pt3Co HIFs的ORR活性大小顺序为Pt3Co(211)> Pt3Co(332)>Pt3Co(331)≈Pt3Co(310)>Pt3Co(511)> Pt3Co(320)。
(1)本研究系统分析了不同Pt3Co HIFs表面的ORR过程, 依据微面标记法选取了六个Pt3Co HIFs, 通过表面能及AIMD计算证明了所选取Pt3Co HIFs的稳定性。
本文的其它图/表