2023, Vol. 37
2. 长江师范学院 电子信息工程学院 超常配位键工程与新材料技术重庆市重点实验室, 重庆 408100
2. Key Laboratory of Extraordinary Bond Engineering and Advanced Materials Technology of Chongqing, School of Electronic Information Engineering, Yangtze Normal University, Chongqing 408100, China
近年来, 随着工业及交通运输的迅速发展, 全球二氧化碳(CO2)排放量逐年增加, 并于2021年再创排放量新高, 达到3.39×1010 t, 若各国不进行有力的减排治理措施, 未来9年内将耗尽剩余的“碳排放”预算, 这会对全球气候造成灾难性的影响[1]. “碳达峰碳中和”已成全球学术界、工业界的热点问题, 研究机构和学者正在探索高效的CO2治理方式. 在众多CO2的转化利用技术中, CO2催化加氢是理想的低碳技术之一, 以CO2为碳源, 将其转化为甲醇(CH3OH)、甲烷(CH4)、甲酸(HCOOH)等高附加值产品, 一方面可以治理二氧化碳, 减少温室气体的排放; 另一方面催化加氢后的产物具有高经济效益和环境效益.
负载金属催化剂是常见的加氢催化剂, 并已在工业上实现商业化应用. 在众多催化剂载体材料中, 石墨烯(Gr)因其优异的力学性能、高比表面积、高电导率及较低的成本而被认为是一种优秀的载体材料. 但无缺陷的完美石墨烯是化学惰性的, 实验研究表明, 掺杂金属原子的方法可以有效调整石墨烯薄片的化学性质, 使其对CO2吸附增强, 并提高表面反应活性[2−3], 且石墨烯上的空位与金属原子结合足够稳定, 可用于化学反应[4]. 金属原子与周围配位原子的相互作用, 会在金属原子周围产生电荷积累, 这块区域表现出较高的催化活性[5]. Esrafili等[6]采用密度泛函方法(DFT)研究掺杂Ti的石墨烯纳米片(Ti-Gr)氢化CO2的能力, 研究表明Ti原子掺杂能够增加石墨烯的表面缺陷, 改善石墨烯表面反应活性, 同时加强催化剂稳定性. Sirijaraensre等[7]采用密度泛函方法用石墨烯负载金属铜研究CO2氢化生成甲酸的催化反应路径, 研究表明石墨烯表面的Cu对H2具有强吸附性, 为H2异裂提供反应位点, 生成的Cu—H易于和CO2完成后续的氢化. Esrafili等[8]利用石墨烯为载体, 采用第一性原理研究金属Pt及Ni作为负载氢化CO2的催化活性, 由于Pt原子中的正电荷较Ni原子更多, 更显著地影响了表面的反应性能, 因此Pt掺石墨烯具有更高的催化活性. Kumar等[9]研究氮掺石墨烯负载铜原子催化剂将CO2氢化为甲醇的产量, 实验对比结果表明氮掺石墨烯负载铜原子催化剂的甲醇产量为吡啶N掺杂石墨烯的2倍, Cu原子掺杂对于提升反应活性, 增加反应物产量具有正向作用.
此外, 调整载体的结构和组成对改良催化活性也是一种方法[10]. 由于N原子比C原子具有更高的电负性, 将N引入配位环境中, 作为金属原子的配位原子[11], 可以改变电子态的空间和能量分布, 从而显著提高金属掺杂石墨烯薄片的表面活性[12]. N掺石墨烯主要以下3种构型: 吡啶N, 石墨N, 吡咯N, 研究表明, 吡啶N构型最有利于金属原子的结合[13−14]. Deerattrakul等[14]采用水热还原的方法将N掺杂在石墨烯中, 研究其对CO2加氢制甲醇产量的影响, 研究表明在3种N掺杂石墨烯构型中, 吡啶N掺杂石墨烯更有利于改善金属的分散, 能够促进H2的解离, 增强CO2的吸附, 从而促进甲醇的产生, 通过实验结果对比发现, N掺杂石墨烯的甲醇产量是未掺杂石墨烯的4倍. 在N掺杂石墨烯基础上, 进一步引入金属替代N原子, 可以形成性能优良的单原子催化剂. Liu等[15]采用第一性原理对PtN3-Gr的研究表明, N掺杂可以在石墨烯的费米能级附近引入局域缺陷态, 使锚定的Pt原子的扩散与聚集更加困难, 从而有效提升掺杂Pt原子的稳定性. Liu等[16]发现FeN4-Gr在CO2-CO转化过程中表现出良好的性能, 且在吸附表面能促进CO2的活化和质子化. Cai等[2]对多种过渡性金属(Sc、Ti、Cu、Zn等)掺杂N4-Gr进行了研究, 研究结果表明, N4-Gr结构可以改善金属原子对CO2的吸附能力, 通过对比得出, Ti原子掺杂N4-Gr具有最强的吸附能力和最高的电荷转移能力. Wang等[17]采用DFT研究NiN3-Gr中Ni结合位点不同对结构稳定性的影响, 研究表明Ni被3个吡啶N原子包围的位置产生的吸附能最低, 吸附高度最小, 为最有利的吸附点.
研究发现, Zr掺杂可显著提高石墨烯对气体分子的吸附, 学者在开发吸附材料、探测材料、储氢材料等方面进行了研究[18−19]. 然而, 目前对于Zr、N共掺杂石墨烯的催化性能及其在CO2催化转化方面的研究尚为不足. 我们利用密度泛函理论计算, 对比了CO2和H2在ZrNx掺杂石墨烯(ZrNx-Gr)上的吸附, 探讨ZrNx-Gr催化CO2加氢制甲酸和CO的反应机理.
1 计算方法及模型采用Vienna Ab-initio Simulation Package (VASP)软件包[20−21], 基于广义梯度近似(GGA-PBE)[22]和投影增强波(PAW)[23]方法对ZrNx-Gr体系进行计算. N3-Gr和N4-Gr是氮掺杂石墨烯的常见结构, 因此在石墨烯原胞基础上, 采用p(5×5)超胞建模, 构建了ZrN3-Gr和ZrN4-Gr体系. 使用5×5×1的Monkhorst-Pack网格进行弛豫和态密度计算, 截断动能取520 eV, 表面结构弛豫收敛的判定标准为Hellmann-Feynman力小于−0.2 eV/nm. 利用baber电荷来评价反应中的电子转移情况, 采用CI-NEB方法搜索过渡态[24].
ZrN3-Gr和ZrN4-Gr的形成能计算公式如下:
| $ {E}_{\text{form}\text{，} \rm{ZrN_3-Gr}} = {E}_{\rm{ZrN_3-Gr}} + 4{\mu }_{\text{C}} - \left({E}_{\text{Gr}} + 3{\mu }_{\text{N}} + {E}_{\text{Zr}}\right) $ | (1) |
| $ {E}_{\text{form}\text{，} \rm{ZrN_4-Gr}} = {E}_{\rm{ZrN_4-Gr}} + 6{\mu }_{\text{C}} - \left({E}_{\text{Gr}} + 4{\mu }_{\text{N}} + {E}_{\text{Zr}}\right) $ | (2) |
其中Eform为形成能, EZrN3-Gr和EZrN4-Gr分别为ZrN3-Gr和ZrN4-Gr的总能量, μN为N2分子能量的一半, μC为石墨烯中单个碳原子的能量, EGr和EZr分别为石墨烯的能量和Zr金属原胞中单个原子的能量.
将CO2等吸附物定义为X, X在ZrNx-Gr上的吸附能(Eads, eV)定义为
| $ {E_{{\text{ads}}}} = {E_{{\text{XC}}}} - {E_{\text{C}}} - {E_{\text{X}}} $ | (3) |
其中, EXC代表吸附物X吸附在ZrNx-Gr时复合界面的总能量, EC与EX分别代表掺杂ZrNx的石墨烯的能量和吸附物X的能量.
吸附的焓变和吉布斯自由能的变化(25 ℃)采用如下公式计算:
| $ \Delta H = \sum\limits_{{\text{products}}} {\left( {{E_0} + {H_{298}}} \right) - \sum\limits_{{\text{reactants}}} {\left( {{E_0} + {H_{298}}} \right)} } $ | (4) |
| $ \Delta G = \sum\limits_{{\text{products}}} {\left( {{E_0} + {G_{298}}} \right) - \sum\limits_{{\text{reactants}}} {\left( {{E_0} + {G_{298}}} \right)} } $ | (5) |
其中, E0是静态计算得到的体系能量, H298和G298是25 ℃时体系的焓值和吉布斯自由能.
2 结果和讨论 2.1 H2和CO2在ZrN3-Gr和ZrN4-Gr表面的吸附首先计算了ZrN3-Gr和ZrN4-Gr的形成能, 分别为4.19和2.33 eV, 因此ZrN4-Gr是更容易形成的结构. 由于催化反应发生的第一步是反应物吸附在催化剂表面, 因此进一步比较了H2和CO2在ZrN3-Gr和ZrN4-Gr上的吸附情况. 图1和图2分别为H2和CO2在ZrN3-Gr和ZrN4-Gr表面最稳定的吸附结构和差分电荷密度(EDD)图, 由H2和CO2吸附引起的吸附能、焓变、吉布斯自由能和电荷转移等数据如表1所示. H2和CO2在ZrN3-Gr和ZrN4-Gr表面吸附的构型没有很大的差异. H2均平行于石墨烯表面在Zr位点吸附, 吸附能分别为−0.49 eV (ZrN3-Gr)和−0.47 eV (ZrN4-Gr), 两个H原子与Zr形成几乎等距的化学键, H―H键从自由态的0.074分别增加到吸附态的0.088 nm (ZrN3-Gr)和0.082 nm (ZrN4-Gr). 由EDD图可见, H原子之间的电荷密度略微减小, Zr与两个H原子之间的电荷密度明显增加, 说明H―H之间相互作用减弱, Zr―H之间相互作用增强.
|
图 1 CO2和H2分子在ZrN3-Gr上吸附的结构图、EDD图 Fig.1 Optimized structure and EDD diagram of CO2 and H2 adsorption on ZrN3-Gr |
|
图 2 CO2和H2分子在ZrN4-Gr上吸附的结构图、EDD图 Fig.2 Optimized structure and EDD diagram of CO2 and H2 adsorption on ZrN4-Gr |
|
表 1 CO2和H2在ZrN3-Gr和ZrN4-Gr表面的吸附能、焓变、吉布斯自由能与电荷转移量 Table 1 Adsorption energy, change in the enthalpy, Gibbs free energy and charge transfer of CO2 and H2 on ZrN3-Gr and ZrN4-Gr surface |
CO2通过C和O原子形成Zr―C键与Zr―O键吸附在Zr上(图1(b)和图2(b)), 同时该O原子与C原子之间的C―O键长增加, 由自由态的0.126增加到吸附态的0.142 nm (ZrN3-Gr)和0.144 nm (ZrN4-Gr), 吸附能分别为−2.17 eV (ZrN3-Gr)和−1.66 eV (ZrN4-Gr). 在EDD图中, CO2中C原子及靠近Zr的O原子与Zr原子间电荷密度显著增加, 说明CO2吸附时伴随有从表面发生的相对较大的净电荷转移, 并形成了Zr―C与Zr―O离子键, 同时, CO2中的C与吸附在Zr上的O原子间的区域电荷密度减小, 说明C―O键减弱, 与C―O键长增加相符.
当H2和CO2共吸附时, CO2的吸附特性与其单独吸附时特性相同, H2相比于单独吸附远离Zr位点(图1(c)和图2(c)), 表明CO2在ZrN3-Gr上发生氢化时, 优先吸附CO2. 两种分子的共吸附能为−2.24 eV (ZrN3-Gr)和1.67 eV (ZrN4-Gr), 小于H2和CO2单独吸附的吸附能之和, 因此共吸附将削弱两者与表面之间的相互作用. 表1数据表明, 相较于单独吸附, 共吸附中H2的电荷转移量明显下降, ZrN3-Gr由单独吸附的0.30降低至共吸附的0.01 e, ZrN4-Gr由单独吸附的0.17降低至0.02 e. CO2的电荷转移量只有轻微的下降, ZrN3-Gr中由1.10下降至1.07 e, ZrN4-Gr中由1.10降低至0.99 e, 说明H2和CO2的相互作用对于H2影响更大. EDD分析佐证了以上结果, ZrN3-Gr上发生共吸附时, H2和CO2分子内部的电荷密度分布情况与其各自单独吸附时的电荷密度分布情况类似, 但两个H原子与Zr间的电荷密度增量相比于H2单独吸附时显著减小, 证明ZrN3-Gr表面对H2的吸附被削弱; ZrN4-Gr中H2单独吸附时的电荷密度增量相比于明显小于ZrN3-Gr, 共吸附时ZrN4-Gr中CO2的电荷密度增量也有所减少, 说明H2和CO2与ZrN4-Gr的结合弱于ZrN3-Gr. 与文献中数据对比可知, 在ZrN3-Gr表面, CO2单独吸附和共吸附吸附能显著高于CoN3-Gr、FeN3-Gr和ZrN3-Gr表面[12,20], 说明ZrN3-Gr表面更利于CO2加氢反应的发生.
投影态密度(PDOS)分析如图3所示, 其中(a)、(b)、(c)为ZrN3-Gr, (d)、(e)、(f)为ZrN4-Gr. 对于H2的吸附(图3(a)和图3(d)), H2轨道在主要位于−8.32 eV (ZrN3-Gr)和−7.92 eV (ZrN4-Gr)处, 并在费米能级处存在态密度; 对于CO2的吸附(图3(b)和图3(e)), CO2轨道在主要位于−8.68、−5.97、−3.06、−0.57 eV (ZrN3-Gr)和−7.22、−4.72、−4.31、−1.60、1.31 eV (ZrN4-Gr)处. 共吸附后(图3(c)和图3(f)), ZrN3-Gr中, H2在−8.32 eV处的能量峰值比其单独吸附时显著降低, 在−8.18、−8.79和−5.94 eV能量处产生了新峰; CO2在−8.32 eV处产生新峰, 同时, 在−5.97 eV处的峰下降, 在−3.06 eV处的峰上升. 在ZrN4-Gr中, H2单独吸附时在−7.92 eV处的能量峰移至−8.59 eV且峰值显著降低, 在−7.97和−5.33 eV处产生新峰; CO2的能量峰相比于单独吸附也整体左移, 在−7.92 eV处产生新峰, 在−8.59 eV处的峰值下降, 在−5.74 eV处的峰值上升. 在共吸附时, 两种结构的H2与CO2的轨道存在交叠, 说明共吸附后H2和CO2之间产生相互影响. ZrN3-Gr上发生共吸附后, CO2轨道整体变化不大, 而H2轨道显著受到CO2的吸附影响, 说明CO2与Zr的相互作用更强. H2吸附时Zr的电荷转移量为1.86 e, CO2吸附和共吸附时的电荷转移量为2.21 e (表1), PDOS图中H2吸附时Zr的态密度明显高于CO2吸附和共吸附时的态密度, 两者保持一致. ZrN4-Gr上发生共吸附后, H2和CO2的轨道均发生了偏移, 轨道出现明显的交叠, 说明共吸附时两者存在较强的相互作用.
|
图 3 CO2和H2分子在ZrN3-Gr和ZrN4-Gr上吸附的PDOS图 Fig.3 PDOS diagram of CO2 and H2 adsorption on ZrN3-Gr and ZrN4-Gr |
综上, 在ZrN4-Gr和ZrN3-Gr上吸附时的吸附结构相似, ZrN3-Gr上各吸附质的吸附能均大于ZrN4-Gr, 25 ℃的焓变和吸附自由能也具有相同的趋势(表1), 说明CO2、H2与ZrN3-Gr的相互作用更强, 更利于CO2加氢反应的发生. 因此, 虽然ZrN4-Gr的形成能更低, 但是CO2加氢反应在ZrN3-Gr上发生的几率更高, 因此将ZrN3-Gr作为进一步分析反应路径的表面.
2.2 CO2加氢生成甲酸的反应路径CO2加氢生成甲酸反应在MN3-Gr表面通过HCOO*途径进行, 该途径又可分为顺式或反式HCOOH两种路径, 经研究发现, 通过反式HCOOH生成甲酸的途径占优[25], 因此, 我们仅讨论CO2在ZrNx-Gr表面通过HCOO*途径生成反式HCOOH的反应路径.
CO2在ZrN3-Gr表面加氢生成甲酸的反应路径和能垒如图4(a)所示. 对于CO2在ZrN3-Gr上催化加氢制甲酸, 反应从初始态IS-1开始, 此时CO2和H2共吸附在Zr原子上. 随着CO2中的C原子与H2中的H1原子的相互作用, H1原子向C原子靠近, 形成过渡态TS-1, 这一过程导致了Zr原子和CO2中C原子间Zr―C键的断裂, 此时CO2与Zr原子间仅通过O2原子连接, 同时, H2的氢原子间的H―H键长从IS-1的0.077增加到TS-1的0.100 nm. 由IS-1到TS-1所需的能垒仅为0.19 eV. 随后, TS-1中的H原子进一步靠近C原子形成C―H键, 最终获得一个甲酸盐(HCOO*)中间体, 即中间态IM-1, HCOO通过O2原子吸附在Zr原子上, 随后, 甲基基团围绕Zr―O键旋转并略微远离Zr原子, 使Zr―O键键长略增, 与此同时, H2原子开始脱附, 键长由0.187增加到0.221 nm, 形成TS-2, 该过程所需能垒为1.85 eV. 之后TS-2上的Zr―H键断裂, H2完成脱附, HCOO*基团绕Zr―O键旋转的同时吸附H2原子迁移靠近O1原子, 最终形成O―H键, 形成终态FS-1, 最终甲酸从FS-1上脱附得到最终产物, 总反应为吸热反应.
|
图 4 CO2在(a)ZrN3-Gr和(b)ZrN4-Gr上加氢制甲酸的路径 Fig.4 The pathway of CO2 hydrogenation to formic acid on (a) ZrN3-Gr and (b) ZrN4-Gr |
CO2在ZrN4-Gr表面通过HCOO*途径生成反式HCOOH的反应路径与在ZrN3-Gr表面相似(图4(b)), 但在ZrN4-Gr表面的反应路径中增加了IM-2的中间态, 该构型由IM-1上的HCOO*中间体经过一定的空间旋转所得. CO2在ZrN4-Gr表面加氢生成HCOOH决速步骤的能垒为2.48 eV, 高于在ZrN3-Gr表面反应的能垒, 因此ZrN3-Gr更利于催化CO2加氢生成甲酸.
2.3 生成CO与H2O的反应路径CO2加氢生成CO与H2O的反应(CO2+H2→H2O+CO)是加氢制甲酸的一个重要的竞争反应[26], 该反应在ZrN3-Gr表面的反应路径和能垒如图5(a)所示. 反应从CO2和H2共吸附在Zr原子上的初始态IS-1*开始, 经H―H键, C―O键的断裂与O―H的生成形成过渡态TS-1*, 随即形成含OH*与CO*基团的中间态IM-1*, 其中TS-1*与IM-1*两者间键能等相差很小, 仅发生空间结构变化, 随后, H2原子从Zr原子脱离, 转移到O1原子上形成过渡态TS-2*, 该过程所需能垒高达1.86 eV, 略大于CO2加氢制甲酸所需能垒. 因此, 在ZrN3-Gr表面进行催化加氢反应可分两条路径, 即HCOO*途径和CO途径.
|
图 5 CO2在(a)ZrN3-Gr和(b)ZrN4-Gr上加氢产CO与H2O的路径 Fig.5 The pathway of CO2 hydrogenation to produce CO and H2O on (a) ZrN3-Gr and (b) ZrN4-Gr |
在ZrN4-Gr表面的反应路径如图5(b)所示, 其反应决速步骤的能垒为1.73 eV, 小于在ZrN3-Gr表面的反应决速步骤, 表明在ZrN4-Gr表面更有利于该副反应的发生.
3 结论通过密度泛函理论研究了氮掺杂石墨烯负载单原子Zr催化CO2加氢的机理, 得到了H2和CO2在ZrN3-Gr与ZrN4-Gr上单独吸附及共吸附3种状态下的吸附结构、吸附能及电子转移情况. CO2在ZrNx-Gr上加氢反应路径以H2和CO2共吸附的结构为起始结构, 在ZrN3-Gr表面的吸附能显著高于ZrN4-Gr, 说明ZrN3-Gr表面利于CO2加氢反应的发生. CO2在共吸附后保持了其单独吸附时的特性, 削弱了H2分子的吸附. ZrNx-Gr上的CO2加氢反应分为两条路径, 一条路径生成HCOOH, 另一条生成CO与H2O, 其中ZrN3-Gr和ZrN4-Gr表面反式HCOOH形成的能垒分别为1.85和2.48 eV, 产生H2O与CO的反应能垒分别为1.86和1.73 eV. 因此, 相比于ZrN4-Gr, ZrN3-Gr更有利于CO2加氢生成甲酸反应的发生, 而ZrN4-Gr表面更利于CO的产生.
| [1] |
a. BP plc. BP Statistical Review of World Energy 2022[R]. London: BP plc, 2022.b. Ye Zhen(叶 朕), Luo Hao-lin(罗皓霖), Jiang Zhi(江 治), et al. Recent advances of photocatalytic CO2 overall reduction(光催化还原二氧化碳全反应的研究进展)[J]. J Mol Catal(China)(分子催化), 2023, 37(2): 174−186.c. Tang Wen-bin(唐文彬), Zhang Zhi-xiang(张志翔), Chi Jia-sheng(池家晟), et al. The state of the art review on the photo-thermal reactor of CO2 reduction(光热催化还原二氧化碳反应器研究进展)[J]. J Mol Catal (China) (分子催化), 2022, 36(5): 499−512
|
| [2] |
First-principles investigation of carbon dioxide adsorption on MN4 doped graphene
[J]. AIP Adv, 2020, 10(12): 125013.
DOI:10.1063/5.0029724 |
| [3] |
Guo Hui-xia(郭惠霞), Yan Ren-xiang(闫任翔), Xi Xiao-hua(席晓华), et al. DFT study on catalytic oxygen reduction reaction of Pt-doped graphene oxide in acidic medium(Pt掺杂氧化石墨烯在酸性介质中催化氧还原反应DFT研究)[J]. J Mol Catal (China)(分子催化), 2018, 32(5): 471−480.
|
| [4] |
Catalytic conversion of CHx and CO2 on non-noble metallic impurities in graphene
[J]. Phys Chem Chem Phys, 2016, 18(25): 16998–17009.
DOI:10.1039/C6CP01403H |
| [5] |
Electrocatalytic decomposition of formic acid catalyzed by M-embedded graphene (M=Ni and Cu): A DFT study[J]. Top Catal, 2021, 65(5/6): 643–655.
|
| [6] |
A DFT study on the catalytic hydrogenation of CO2 to formic acid over Ti-doped graphene nanoflake
[J]. Chem Phys Lett, 2017, 682: 49–54.
DOI:10.1016/j.cplett.2017.06.011 |
| [7] |
Hydrogenation of CO2 to formic acid over a Cu-embedded graphene: A DFT study
[J]. Appl Surf Sci, 2016, 364: 241–248.
DOI:10.1016/j.apsusc.2015.12.117 |
| [8] |
Catalytic hydrogenation of CO2 over Pt- and Ni-doped graphene: A comparative DFT study
[J]. J Mol Graph Model, 2017, 77: 143–152.
DOI:10.1016/j.jmgm.2017.08.016 |
| [9] |
Nitrogen-doped graphene-supported copper complex: A novel photocatalyst for CO2 reduction under visible light irradiation
[J]. RSC Adv, 2015, 5(68): 54929–54935.
DOI:10.1039/C5RA05319F |
| [10] |
Modulating the local coordination environment of single-atom catalysts for enhanced catalytic performance[J]. Nano Res, 2020, 13(7): 1842–1855.
DOI:10.1007/s12274-020-2755-3 |
| [11] |
Lin Yan(林 研), Yue Hai-rong(岳海荣), Liu Chang-jun(刘长军), et al. Theoretical calculation of CO2 hydrogenation to formic acid catalyzed by Ag monatom supported N modified graphene(Ag单原子负载N改性石墨烯催化CO2加氢制甲酸的理论计算)[J]. Chem Ind Eng(化学工业与工程), 2023, 40(2): 25−32.
|
| [12] |
Theoretical insights into hydrogenation of CO2 to formic acid over a single Co atom incorporated nitrogen-doped graphene A DFT study
[J]. Appl Surf Sci, 2019, 475: 363–371.
DOI:10.1016/j.apsusc.2018.12.302 |
| [13] |
Hollow N-doped carbon spheres with isolated cobalt dingle atomic sites: Superior electrocatalysts for oxygen reduction[J]. J Am Chem Soc, 2017, 139(48): 17269–17272.
DOI:10.1021/jacs.7b10194 |
| [14] |
The roles of nitrogen species on graphene aerogel supported Cu-Zn as efficient catalysts for CO2 hydrogenation to methanol
[J]. Appl Catal A-Gen, 2019, 580: 46–52.
DOI:10.1016/j.apcata.2019.04.030 |
| [15] |
Monodisperse Pt atoms anchored on N-doped graphene as efficient catalysts for CO oxidation: A first-principles investigation[J]. Catal Sci Technol, 2015, 5(3): 1658–1667.
DOI:10.1039/C4CY01327A |
| [16] |
Constructing FeN4/graphitic nitrogen atomic interface for high-efficiency electrochemical CO2 reduction over a broad potential window
[J]. Chem, 2021, 7(5): 1297–1307.
DOI:10.1016/j.chempr.2021.02.001 |
| [17] |
Single Ni atom incorporated with pyridinic nitrogen graphene as an efficient catalyst for CO oxidation: First-principles investigation[J]. RSC Adv, 2017, 7(77): 48819–48824.
DOI:10.1039/C7RA06591D |
| [18] |
Adsorption of phosgene molecule on the transition metal-doped graphene: First principles calculations[J]. Appl Surf Sci, 2017, 425: 340–350.
DOI:10.1016/j.apsusc.2017.06.229 |
| [19] |
Density functional theory calculations of hydrogen adsorption on Ti-, Zn-, Zr-, Al- and N-doped and intrinsic graphene sheets[J]. Int J Hydrogen Energ, 2013, 38(33): 14269–14275.
DOI:10.1016/j.ijhydene.2013.07.098 |
| [20] |
Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set[J]. Phys Rev B, 1996, 54(16): 11169.
DOI:10.1103/PhysRevB.54.11169 |
| [21] |
Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set[J]. Comp Mater Sci, 1996, 6(1): 15–50.
DOI:10.1016/0927-0256(96)00008-0 |
| [22] |
Generalized gradient approximation made simple[J]. Phys Rev Lett, 1996, 77(18): 3865.
DOI:10.1103/PhysRevLett.77.3865 |
| [23] |
Projector augmented-wave method[J]. Phys Rev B, 1994, 50(24): 17953–17979.
DOI:10.1103/PhysRevB.50.17953 |
| [24] |
A climbing image nudged elastic band method for finding saddle points and minimum energy paths[J]. J Chem Phys, 2000, 113(22): 9901–9904.
DOI:10.1063/1.1329672 |
| [25] |
Theoretical insights into catalytic CO2 hydrogenation over single-atom (Fe or Ni) incorporated nitrogen-doped graphene
[J]. J CO2 Util
, 2021, 48: 101532.
DOI:10.1016/j.jcou.2021.101532 |
| [26] |
Hydrogen-assisted transformation of CO2 on nickel: The role of formate and carbon monoxide
[J]. J Phys Chem Lett, 2010, 1(1): 402–406.
DOI:10.1021/jz900221c |