析氢反应中氮掺杂石墨烯负载金属单/双原子催化活性起源
1.
2.
Activity origin of single/double-atom catalyst for hydrogen evolution reaction
1.
2.
收稿日期: 2021-07-31 修回日期: 2021-08-31 网络出版日期: 2021-11-03
基金资助: |
|
Received: 2021-07-31 Revised: 2021-08-31 Online: 2021-11-03
作者简介 About authors
张诗诗(1997—),女,硕士研究生,主要研究方向为计算电化学催化,E-mail:
苏亚琼,研究员,主要研究方向为计算能源催化/材料,E-mail:
电催化析氢反应(HER)是极具潜力的氢能转化手段,发展高性能、低成本的析氢催化剂是制氢技术的关键。本文选取了以过渡金属(Fe、Ni、Co)为活性中心、以氮掺杂石墨烯(N-graphene)为基底的单、双原子催化剂(SAC、DAC),首先通过理论计算对其结构稳定性做出了讨论,证明了所选催化材料的可行性。然后,采用H吸附能为描述符,对催化剂的HER活性进行了分析讨论。结果表明CoN4表现出最优异的HER活性,FeN3、FeN4、Fe2N6三种结构的Fe原子活性中心次之。最后,系统讨论了催化剂的电子结构,揭示了HER催化活性的来源。本文从理论计算角度对SAC、DAC的HER活性进行了对比,结果表明DAC具有较差的HER活性,SAC中CoN4、FeN3、FeN4在HER中展现出优异的催化性能,可作为取代常用的贵金属(Pt/C)HER催化剂的理想选择。
关键词:
Electrocatalytic hydrogen evolution reaction (HER) is a promising hydrogen energy conversion method. To develop high-performance and low-cost hydrogen evolution electrocatalysts, single- and double-atom catalysts (SACs, DACs) with transition metals (e.g., Fe, Ni, and Co) as the active center and nitrogen-doped graphene (N-graphene) as the substrate are selected for HER utilizing density functional theory calculations. The selected catalytic materials exhibit exceptional stability against sintering. We then chose H adsorption energy as the descriptor for analyzing the HER activity, and the results demonstrate that the CoN4 site exhibits excellent HER activity over other candidates. In contrast, NiN4 and Ni2N6 sites display inferior HER activity. In addition, the electronic structures of the catalysts are systematically discussed to uncover the origin of catalytic activity. This work reveals that DACs have poor HER activity compared to SACs and DACs, while the SACs (e.g., CoN4, FeN3, and FeN4) show low overpotential in HER. Therefore, the SACs can substitute commercial precious metals catalysts (Pt/C) for HER catalysts.
Keywords:
本文引用格式
张诗诗, 秦棪阳, 苏亚琼.
ZHANG Shishi.
氢气具有含量丰富、可再生等特点,可作为取代化石燃料的新型清洁能源。电催化析氢反应(HER)发生在水分解反应中,可实现高纯度、高产量的产氢效果,是一种可持续的安全可靠的氢气生产技术。然而HER通常需要较大的过电位(η),为降低过电位,减少电能消耗,设计高效低廉的HER电催化剂是当前制氢技术的发展趋势。近年来,基于过渡金属的电催化剂被广泛研究,过渡金属具有独特的d电子结构,在电催化反应中表现出优异的催化活性。
单原子催化剂(SAC)仅具有单个活性位点,可实现100%原子利用率,是电催化中极具潜力的贵金属催化剂替代品。通过实验以及理论研究,SAC的催化活性已在多种电催化反应中得到证明[1-2],Xu等人[3]利用理论计算证明了MN4/C结构的SAC在氧还原反应(ORR)中表现出高催化活性;Zhao等人[4]表明SAC在氮气还原反应(NRR)中也表现出高催化性能;HER中,SAC的高催化性能也得到了证明[5-7]。Cheng等人[8]制备出Pt原子催化剂用于HER中,实验结果表明相对于商业中使用的Pt/C催化剂,SAC具有高催化活性、高稳定性;Lei等人[9]利用Ni纳米粒子对Ni-N-C结构的SAC进行调控,使其在HER中催化性能得到提升。由于SAC的单活性位点,其催化性能具有局限性,具有双位点的双原子催化剂(DAC)被引入,DAC中第二个活性位点对第一个活性位点的电子结构起到了调控作用,在NRR、ORR等电化学反应中表现优异[10-11]。
为探明SAC、DAC在HER中的催化性能,本文主要选取了以三种过渡金属(Fe、Co、Ni)为活性中心的SAC、DAC,通过第一性原理计算对其结构稳定性、HER活性、电子结构做出了分析研究,为高效HER电催化剂的研究提供理论支撑,并对新型HER催化剂的设计进行了展望。
1 计算方法
本文所有工作均使用基于密度泛函理论(DFT)的计算软件VASP。体系中波函数采用平面波基组的展开方式,赝势选择PAW,交换相关能使用广义梯度近似(GGA)的PBE表示。所有体系选取的平面波截断能为400 eV,k点抽样选择Monkhorst-Pack抽样法,k网格点为3×3×1。
为评估催化剂结构稳定性,引入结合能,其表达式为:Eb=Etot-Esupport-ETM,其中,Eb为体系结合能,Etot为整个结构的能量,Esupport、ETM分别为基底以及过渡金属原子的能量。结合能可表明各组成部分结合的紧密程度。
为评估催化活性,中间反应步的吉布斯自由能变化由以下公式得到:∆G=∆E+∆ZPE+
2 结果与分析
2.1 结构稳定性
本文选取Fe、Ni、Co三种过渡金属元素作为活性中心,选取氮掺杂石墨烯(N-graphene)作为基底,设计了FeN3、FeN4、Fe2N6、CoN4、Co2N6、NiN4、Ni2N6七种原子催化剂。其结构为如图(1)所示的六方晶格体系,基矢为:a=b=14.76 Å,c=15 Å,真空层厚度为15 Å(1 Å=0.1 nm)。
图1
图1
(a) FeN3、(b) FeN4、(c) Fe2N6结构的侧视及俯视图,橙色、蓝色、灰色小球分别为Fe、N、C原子
Fig. 1
The side and top view of the geometry optimization structures of (a) FeN3, (b) FeN4, (c) Fe2N6-graphene. Cray, blue, orange balls represent C, N, Fe atoms, respectively
表1 各SAC、DAC结构中TM—N键长、结合能
Table 1
Item | FeN3 | FeN4 | CoN4 | NiN4 | Fe2N6 | Co2N6 | Ni2N6 |
---|---|---|---|---|---|---|---|
TM—N/Å | 1.92 | 1.89 | 1.87 | 1.87 | 1.79/1.95 | 1.78/1.94 | 1.83/1.87 |
Eb/eV | -4.91 | -7.59 | -7.42 | -4.64 | -14.41 | -13.90 | -12.89 |
图2
图2
(a) CoN4、(b) Co2N6电荷密度差分图,C、N、Co原子分别用棕色、蓝色、红色小球表示,蓝色、黄色阴影分别代表电荷消耗和聚集
Fig. 2
The differential charge density maps of (a) CoN4, (b) Co2N6, Brown, blue, red balls represent C, N, Co atoms, respectively, blue and yellow shadows characterize the charge loss and accumulation
2.2 HER催化活性
电催化析氢反应(HER)是在电极/电解质界面处发生的多步电化学反应过程,主要是通过还原质子(H+)生成H2:2(H++e-)→H2。其反应步骤分为两步,第一步反应为质子的吸附(Volmer反应),第二步为H2的脱附(Heyrovsky/Tafel反应)。
为保持较优的反应速率,HER中氢原子的吸附和脱附需要达到良好的平衡[14]。H原子吸附能是普遍采用的HER活性描述符[7,15],为研究在HER中所选取SAC、DAC的活性,我们利用VASP计算得出其吸附H的自由能变化(ΔG)。根据Sabatier原理,ΔG接近0时,HER催化效率达到最大,因此,高效HER电催化剂对H的吸附能力要保持适中水平。图3是HER在7种催化剂上进行时的自由能变化图,过电位已在图中标出。可以看出CoN4具有最小的过电位数值(-0.04 V);FeN3、FeN4、Fe2N6过电位分别为0.10、0.29和0.36 V,HER活性较好;Co2N6,Ni2N6,NiN4过电位较高(0.76、1.35、1.28 V),对H吸附能力太弱,导致不理想的HER活性。Tour等人[16]在实验中合成了分散在N掺杂石墨烯上的单Co原子催化剂(Co-NG),并证明在HER中其过电位仅30 mV,这与本文理论计算所得结论相近。Lu等人[6]通过改变N配位数,对CoNx-gra(x=1~4)体系在HER中催化性能进行了对比,从实验角度证明了CoN4-gra具有最高的催化活性。Fe-Nx作为高效电催化剂活性中心,在NRR、ORR等电化学反应中应用广泛,但少有报道将其应用于HER中,本文DFT计算结果表明FeN3、FeN4、Fe2N6三种以Fe原子为活性中心的催化材料都具有作为HER电催化剂的潜能。另外,本文将SAC与DAC进行对比,利用DFT计算证明:对于HER,SAC催化活性高于DAC。Nørskov等人[15]发现金属表面引入H2O分子层会使H吸附能降低0.02 eV,因此H2O分子质子化的微弱影响可以忽略。
图3
图3
HER中各催化剂上能量变化图
Fig. 3
The free energy diagrams of HER on all catalysts studied. The overpotentials (η) are presented
图4
图4
(a) CoN4中Co的d轨道态密度图,(b) d带中心值和HER过电位的近线性关系示意图
Fig. 4
(a) Projected density of state (PDOS) projected on 3d orbitals of Co atom of CoN4. (b) The calculated linear relation of the ORR overpotential (-ηORR) as a function of d-band center (εd)values of Co atom
图5
图5
CoN4吸附H后(a)电荷密度差分图,棕色、蓝色、红色、绿色小球分别代表C、N、Co、H原子,黄色和蓝色阴影分别代表电荷聚集和电荷消耗;(b) Co的3d轨道以及被吸附*H的1s轨道的态密度图
Fig. 5
(a) The differential charge density maps where the brown, blue, red, green balls denote the C, N, Co, H atoms, respectively; (b) The projected density of state (PDOS) projected on 3d orbitals of Co atom and the 1s orbital of adsorbed H of CoN4-H
3 结论
本文选取3种不同金属中心的SAC、DAC,通过DFT计算分析讨论了HER中各催化剂的活性,并对催化剂电荷结构进行了分析。结果表明CoN4吸附H步骤具有最接近于0的自由能变化,表现出最优异的催化性能,FeN3、FeN4次之(∆G=0.10,0.29 eV);DAC在HER中具有较差的催化活性;Ni原子作为活性中心无法很好地催化HER过程。通过电荷分析发现d带中心和催化活性之间存在着近线性关系;吸附H时金属活性中心和H原子之间具有较强的相互作用。
参考文献
/
〈 |
|
〉 |
