三原子催化剂(TACs)在二氧化碳还原反应(CO2RR)中展现了巨大潜力,设计新型三原子催化剂具有重要意义。本工作利用第一性原理,建立了单层多孔氮化石墨烯(C2N)负载的过渡金属三聚体(3TM-C2N, TM = Mn, Mo, Ru, Ti)催化剂用于将二氧化碳(CO2)还原为甲烷(CH4)。计算结果表明,3TM-C2N催化剂结构稳定,不但能有效吸附和活化CO2,而且对析氢反应(HER)有良好的抑制性。吉布斯自由能分布图显示,CO2RR在3TM-C2N上会以不同的反应路径生成CH4。极限电势(UL)分析显示,3Mn-C2N表现出了最好的催化性能,对应的UL为−0.44 V。这些发现不仅为实验上调控C2N基催化剂提供了理论依据,还对开发其他高效的CO2RR电催化剂有一定的指导意义。 Triple-Atom Catalysts (TACs) have shown great potential in carbon dioxide reduction reactions (CO2RR). It is of great significance to design new Triple-Atom Catalysts. In this work, a monolayer porous nitrogen-doped graphene (C2N) supported transition metal trimers (3TM-C2N, TM = Mn, Mo, Ru, Ti) catalyst models were established to catalyze the reduction of carbon dioxide (CO2) to me-thane (CH4) using first-principles. The results show that the structure of the 3TM-C2N catalyst is stable, not only can effectively adsorb and activate CO2, but also has good inhibition of Hydrogen Evolution Reaction (HER). The Gibbs free energy profiles show that the CO2RR produces CH4products in different reaction paths on 3TM-C2N. The limiting potential (UL) analysis showed that 3Mn-C2N ex-hibited the best catalytic performance, UL of −0.44 V. These findings not only provide a theoretical basis for the experimental regulation of C2N-based catalysts, but also provide guidance for the development of other efficient CO2RR electrocatalysts.
三原子催化剂(TACs)在二氧化碳还原反应(CO2RR)中展现了巨大潜力,设计新型三原子催化剂具有重要意义。本工作利用第一性原理,建立了单层多孔氮化石墨烯(C2N)负载的过渡金属三聚体(3TM-C2N, TM = Mn, Mo, Ru, Ti)催化剂用于将二氧化碳(CO2)还原为甲烷(CH4)。计算结果表明,3TM-C2N催化剂结构稳定,不但能有效吸附和活化CO2,而且对析氢反应(HER)有良好的抑制性。吉布斯自由能分布图显示,CO2RR在3TM-C2N上会以不同的反应路径生成CH4。极限电势(UL)分析显示,3Mn-C2N表现出了最好的催化性能,对应的UL为−0.44 V。这些发现不仅为实验上调控C2N基催化剂提供了理论依据,还对开发其他高效的CO2RR电催化剂有一定的指导意义。
C2N,三原子催化剂,二氧化碳还原反应,CH4,第一性原理
Likai Tong, Yu Zhang, Bo Zhang*
State Key Laboratory of Information Photonics and Optical Communications, School of Integrated Circuits, Beijing University of Posts and Telecommunications, Beijing
Received: Oct. 14th, 2022; accepted: Nov. 4th, 2022; published: Nov. 14th, 2022
Triple-Atom Catalysts (TACs) have shown great potential in carbon dioxide reduction reactions (CO2RR). It is of great significance to design new Triple-Atom Catalysts. In this work, a monolayer porous nitrogen-doped graphene (C2N) supported transition metal trimers (3TM-C2N, TM = Mn, Mo, Ru, Ti) catalyst models were established to catalyze the reduction of carbon dioxide (CO2) to methane (CH4) using first-principles. The results show that the structure of the 3TM-C2N catalyst is stable, not only can effectively adsorb and activate CO2, but also has good inhibition of Hydrogen Evolution Reaction (HER). The Gibbs free energy profiles show that the CO2RR produces CH4products in different reaction paths on 3TM-C2N. The limiting potential (UL) analysis showed that 3Mn-C2N exhibited the best catalytic performance, ULof −0.44 V. These findings not only provide a theoretical basis for the experimental regulation of C2N-based catalysts, but also provide guidance for the development of other efficient CO2RR electrocatalysts.
Keywords:C2N, Triple-Atom Catalyst, Carbon Dioxide Reduction, CH4, First-Principles
Copyright © 2022 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
化石燃料的消耗逐年增加导致二氧化碳(CO2)在大气中不断累积,造成了严重的温室效应和海洋酸化等环境问题,迫使人们寻找有效途径来降低大气中CO2的浓度 [
近年来,由于纳米技术的发展,人们可以在纳米尺度进行结构微调从而控制电催化剂的结构,使其拥有更好的稳定性、催化活性和产物选择性 [
因此,在本工作中,通过第一性原理,建立了C2N负载的三原子催化剂3TM-C2N (TM = Mn, Mo, Ru, Ti)模型,并系统地研究了3TM-C2N催化CO2RR生成CH4的性能。首先,计算了3TM-C2N的稳定性及电子结构;其次,从d带中心理论及电子结构的角度,深入分析了3TM-C2N对CO2的吸附与活化;最后,阐明了CO2在3TM-C2N上还原为CH4的最有利路径及其极限电势。计算结果表明,3TM-C2N催化剂结构稳定,能有效吸附和活化CO2,并且对HER有良好的抑制性。其中,3Mn-C2N表现出最好的催化活性,其极限电势为−0.44 V。这些发现不仅为实验上调控C2N基催化剂提供了理论依据,还对开发其他高效的CO2RR电催化剂有一定的指导意义。
基于具有周期性边界条件的密度泛函理论(DFT),本工作中的所有任务都是在Vienna Ab-initio Simulation Package (VASP)中进行计算的 [
本文采用了计算氢电极模型(CHE)计算了标准条件下质子–电子对与氢气的电化学势。
H + ( aq ) + e − ↔ 1 / 2 H 2 ( g ) (1)
依据该方法,每个基本反应的吉布斯自由能变化值计算公式如下 [
Δ G = Δ E + Δ E ZPE − T Δ S (2)
其中, Δ E 为DFT计算得到的能量差值, Δ E ZPE 为零点能的变化,T为温度(298.15 K), Δ S 为熵的变化。
吸附能大小可以判断CO2吸附在催化剂表面时的稳定程度。吸附能公式如下 [
Δ E ads = E system − ( E surf + E react ) (3)
其中, E system 为总能量, E surf 和 E react 分别代表了与吸附分子作用的表面模型和真空中独立的的吸附分子的能量。通常来说,吸附能为负值,表明吸附过程是放热反应,吸附的系统是稳定的。
CO2吸附在催化剂表面的机理可以由d带中心理论阐述,计算d带中心的公式如下:
ε d = ∫ − ∞ ∞ n d ( ε ) ε d ε ∫ − ∞ ∞ n d ( ε ) d ε (4)
其中, n d ( ε ) 为对应的d带上电子的密度, ε 为能量。
另外,极限电势( U L )是使每个基本步骤都能放热的最小负电位。公式为:
U L = − Δ G max / e (5)
其中, Δ G max 为整个CO2RR路径中自由能变化的最大值,e为电子电量 [
优化后的二维C2N晶格参数为a = b = 8.32 Å,与文献中8.30 Å非常接近,说明计算参数设置的比较合理 [
图1. (a) 优化后2 × 2 C2N超晶胞的稳定结构图;(b)~(e) 优化后3TM-C2N的稳定结构图
高稳定性是催化剂应用的先决条件,所以我们用内聚能对3TM-C2N催化剂的稳定性进行了评估。内聚能越高,对应的结构越稳定。3TM-C2N催化剂的内聚能为3Mn-C2N (6.67 eV/atom) < 3Mo-C2N (6.71 eV/atom) < 3Ti-C2N (6.76 eV/atom) < 3Ru-C2N (6.79 eV/atom)。它们的内聚能虽然比初始C2N (6.82 eV/atom)略小,但是要高于碳磷化物(4.12~6.45 eV/atom) [
为了深入理解金属三聚体和C2N之间的相互作用,我们计算了3TM-C2N电子态密度和差分电荷密度,如图2所示。图2(a)~图2(d)电子态密度图显示,在费米能级附近,金属三聚体的d轨道(图中黄色部分)与C2N的中N的p轨道(图中青色部分)之间有明显的重叠峰(图中绿色部分),说明3TM-C2N上的载流子浓度较高,且金属三聚体与C2N的相互作用强烈。尤其在3Ru-C2N中,重叠峰更多,这说明Ru原子与N原子之间相互作用比其他金属原子和N原子之间相互作用更强烈,这是3Ru-C2N催化剂稳定性最高的原因。同时,也正是因为较强的相互作用,导致Ru原子与N原子之间键能更强,键长更短,与图1中有一个Ru原子与C2N处于同一平面的构型相符。
图2(e)~图2(h)差分电荷密度图显示,正电荷在金属原子周围积累(图中青色部分),在N原子附近减少(图中黄色部分),说明电子从金属原子转移到了N原子上,从而使N原子与TM原子结合。
图2. (a)~(d) 3TM-C2N的电子态密度图;(e)~(h) 3TM-C2N的差分电荷密度图
CO2吸附是CO2RR的第一步,也是至关重要的一步。稳定的CO2吸附是CO2RR能继续下去的先决条件。CO2的活化程度一定程度上反映了催化剂的催化性能。如图3所示,CO2能稳定吸附在催化剂表面,并且CO2分子构型发生了明显的变化说明CO2得到了初始活化。为了定量分析,将CO2吸附能和CO2键长键角等信息在表1中列出。CO2吸附在3Mn-C2N、3Mo-C2N、3Ru-C2N和3Ti-C2N催化剂表面的吸附能分别为−1.27、−2.97、−1.59和−3.73 eV。负值的吸附能表示CO2能稳定吸附在3TM-C2N催化剂表面。C-O键长从1.20~1.42 Å的变化以及O-C-O键角从108.97~128.80˚的弯曲都说明了CO2得到了有效的活化。
为了深入探究CO2在催化剂表面的吸附反应,我们计算了金属三聚体的d带中心和差分电荷密度图。图4(a)显示金属三聚体的d带中心与吸附能呈线性关系,相关系数R = −0.75,与d带中心理论模型符合较好。为了更深入的了解金属三聚体上CO2的吸附作用,我们从三维轨道中提取了alpha(α)态和beta(β)态的d带中心位置如图4(b),图4(c)所示。α态和β态的d带中心变化趋势和总d带中心变化趋势一致,但α态d带中心与吸附能线性关系的相关系数为R = −0.99,相关性更好,说明CO2吸附在催化剂表面时与金属三聚体的α轨道相互作用更强。图4(d)~图4(g)中的差分电荷密度显示,正电荷在金属原子上积累(图中青色部分)。电子由金属原子向CO2分子转移使CO2得到了活化。
Models | E ads ( eV ) | Δ G ( eV ) | d C − O (Å) | O-C-O angles |
---|---|---|---|---|
3Mn-C2N | −1.27 | −0.63 | 1.42/1.31 | 110.58˚ |
3Mo-C2N | −2.97 | −2.33 | 1.37/1.37 | 115.68˚ |
3Ru-C2N | −1.59 | −0.93 | 1.20/1.41 | 128.80˚ |
3Ti-C2N | −3.73 | −3.10 | 1.41/1.41 | 108.97˚ |
表1. 3TM-C2N上吸附CO2最稳定构型的吸附能( E ads )、CO2吸附吉布斯自由能变化值( Δ G )、CO2中C和O原子的键长( d C − O )和O-C-O角度
图3. (a)~(d) 优化后3TM-C2N吸附CO2的稳定结构图
图4. (a) 3TM-C2N中过渡金属三聚体的总d带中心与CO2吸附能的线性关系图;(b) 3TM-C2N中过渡金属三聚体 态轨道的d带中心与CO2吸附能的线性关系图;(c) 3TM-C2N中过渡金属三聚体β态轨道的d带中心与CO2吸附能的线性关系图;(d)~(g) 3TM-C2N吸附CO2时的差分电荷密度图
众所周知,在动力学上容易发生的HER作为副反应阻碍了CO2RR的进行。所以一个好的催化剂能抑制HER是必要的。表1和表2显示,在3Mn-C2N、3Mo-C2N、3Ru-C2N和3Ti-C2N上的氢吸附吉布斯自由能分别为−0.99、−2.44、−0.43和−0.89 eV。 | Δ G ( * H ) | > 0.4 eV说明HER在3TM-C2N催化剂上较难发生 [
图5. 在3TM-C2N上发生二氧化碳还原反应和析氢反应第一步加氢反应的自由能变化
Models | Δ G ( * H ) (eV) | E ads ( * H ) (eV) |
---|---|---|
3Mn-C2N | −0.99 | −1.24 |
3Mo-C2N | −2.44 | −2.67 |
3Ru-C2N | −0.43 | −0.67 |
3Ti-C2N | −0.89 | −1.13 |
表2. 3TM-C2N上吸附H原子最稳定构型的吉布斯自由能变化值( Δ G )和吸附能( E ads )
图6显示了在3TM-C2N催化剂上,CO2通过8电子路径生成CH4的中间体吸附结构和自由能分布图。CO2在3Mn-C2N和3Mo-C2N上生成CH4的路径为:CO2→ *CO2→ *COOH → *CO → *COH → *HCOH → *CH2OH → *CH2→ *CH3→ CH4。对于3Mn-C2N,速率限制步骤是*COH → *HCOH,吉布斯自由能变化最大值为−0.44 eV。对于3Mo-C2N,速率限制步骤是*CH3→ CH4,吉布斯自由能变化最大值为−0.97 eV。CO2在3Ru-C2N上生成CH4的路径为:CO2→ *CO2→ *COOH → *HCOOH → *HCO → *HCOH → *CH → *CH2→ *CH3→ CH4,速率限制步骤是*CH2→ *CH3,吉布斯自由能变化最大值为−0.73 eV。CO2在3Ti-C2N上生成CH4的路径为:CO2→ *CO2→ *COOH → *CO → *COH → *HCOH → *CH → *CH2→ *CH3→ CH4,速率限制步骤是*CO → *COH,吉布斯自由能变化最大值为−1.63 eV。
图6. (a) 反应中间体沿b反应路径的对应稳定结构图;(b) 在3TM-C2N上CO2还原产生CH4的吉布斯自由能分布图(数字0~8代表电子转移数)
图7给出了CO2在3Mn-C2N、3Mo-C2N、3Ru-C2N和3Ti-C2N上还原为CH4的极限电势(UL)。这里近零的UL表示容易生成该产物。生成CH4的UL按递增排列都为3Mn-C2N < 3Ru-C2N < 3Mo-C2N < 3Ti-C2N。综合来看,3Mo-C2N和3Ti-C2N整体性能较差, | U L | > 0.9 V说明反应需要吸收大量的能量而较难发生。性能最好的催化剂是3Mn-C2N。在3Mn-C2N上生成CH4的极限电势是−0.44 V,表现出较好的催化活性,是四种催化剂中所需能量最低的。分析表明,3TM-C2N催化CO2RR是可以实现的。我们的工作探索了3TM-C2N在较低电极电位的电催化下对CO2深度加氢具有优异的性能。尤其是3Mn-C2N在四种催化剂中性能最好。这一结果表明,C2N负载的过渡金属三聚体催化剂具有良好的CO2RR潜力。
图7. CO2在3TM-C2N上还原生成CH4时的UL(V)计算值
总之,本工作通过密度泛函理论计算,设计了3TM-C2N (TM = Mn, Mo, Ru, Ti)的模型作为催化剂,讨论3TM-C2N电催化还原CO2为CH4的催化性能。计算表明,锚定在C2N上的锰、钼、钌和钛三聚体构成了稳定的催化剂。6.67~6.79 eV/atom的内聚能表现出强的稳定性。过渡金属三聚体和C2N之间明显的电荷密度积累、明显的轨道重叠和大量的电荷转移是3TM-C2N整体结构高稳定性的原因。CO2吸附表明,3TM-C2N通过三角形结构和TM活性位点的协同作用,为CO2吸附和初始活化提供了有利的环境。CO2在3TM-C2N上的吸附能在–0.63~–3.10 eV范围内。O-C-O键角的弯曲以及C-O键长从1.20~1.42 Å都说明了CO2的有效活化。吉布斯自由能的计算表明,3TM-C2N对CO2较强的吸附以及对中间体*COOH/*OCHO较低的吉布斯自由能表现出对HER良好的抑制性。极限电势表明,3Mn-C2N表现出最好的催化活性,其将CO2还原为CH4的速率限制步骤是*COH → *HCOH,UL为−0.44 V。
本研究得到国家自然科学基金资助(批准号:12174035)。
佟立凯,张 钰,张 博. C2N负载的过渡金属三聚体催化剂用于CO2电催化还原制CH4:第一性原理研究Transition Metals Trimers on C2N as Electrochemical Catalysts for CO2Reduction to CH4: A First-Principles Study[J]. 纳米技术, 2022, 12(04): 330-339. https://doi.org/10.12677/NAT.2022.124033
https://doi.org/10.1039/c3ee41973h
https://doi.org/10.1002/anie.202012329
https://doi.org/10.3390/catal9070604
https://doi.org/10.1039/C9DT01691K
https://doi.org/10.1021/acs.jpclett.5b01559
https://doi.org/10.1021/acs.accounts.0c00626
https://doi.org/10.1021/acsami.9b09941
https://doi.org/10.1021/acscatal.9b02443
https://doi.org/10.1039/C4CS00236A
https://doi.org/10.1002/adma.201807166
https://doi.org/10.1021/acscatal.0c01950
https://doi.org/10.1039/C9TA13599E
https://doi.org/10.1039/C9SC05236D
https://doi.org/10.1039/C8NR04961K
https://doi.org/10.1002/anie.201712451
https://doi.org/10.1016/j.chempr.2017.09.014
https://doi.org/10.1002/anie.201805871
https://doi.org/10.1016/j.apsusc.2021.150331
https://doi.org/10.1016/j.apmt.2021.101245
https://doi.org/10.1039/C8TA11416A
https://doi.org/10.1021/acscatal.0c01641
https://doi.org/10.1038/ncomms7486
https://doi.org/10.1016/j.carbon.2016.04.059
https://doi.org/10.1002/jcc.21057
https://doi.org/10.1103/PhysRevB.54.11169
https://doi.org/10.1063/1.2338035
https://doi.org/10.1103/PhysRevLett.77.3865
https://doi.org/10.1039/C7CP06155B
https://doi.org/10.1021/acs.jpclett.0c02047
https://doi.org/10.1039/C5NR03895B
https://doi.org/10.1039/D0TA03262J
https://doi.org/10.1039/c3cp00066d
https://doi.org/10.1016/j.carbon.2020.10.003