林森ACS Catalysis：原始石墨烯的H2活化起源

二维石墨烯因其独特的性质在各领域得到了越来越多的关注。理想的石墨烯模型由均一的sp²杂化的C原子组成，剩下的p_z电子在面内形成离域π键。因此，无缺陷的原始石墨烯通常被认为是惰性的。赋予石墨烯催化活性的普遍做法是引入外界的活性位点，如负载金属，表面官能化。尽管如此，石墨烯自身的催化性能以及调控却鲜有关注。2014年，Primo等人发现(Nat. Commun. 529, 2014, 5, 5291)，原始石墨烯具有类似于甚至优于商业催化剂（Pd/C、Pt/C和Ni/C）的不饱和键加氢性能。同时，他们观察到，包含N/P/S等杂质的石墨烯的整体催化性能整体下降。该研究一个关键的发现是，石墨烯具有H₂活化的能力，且向反应气中混入CO₂/NH₃可以可逆地影响石墨烯的整体活性。因而，他们认为石墨烯表面的活性位点是具备不同的得失电子能力的，并推测石墨烯表面存在活化H₂的受阻路易斯酸碱对（FLP）活性位点。尽管如此，实验方面仍缺乏足够的证据来阐明该催化机制。

福州大学化学学院林森教授等人基于第一性原理方法，探索了原始石墨烯中的褶皱对H₂解离活性的影响。结果表明，褶皱可以导致石墨烯C位点的电子结构差异化，表现出不同的得失电子能力，并降低H₂解离能垒，H₂在活化过程中表现出相反的电荷，即H₂极化后再发生H–H键断裂。

图文解读

石墨烯中C原子都为sp²杂化，每个C最多还能与一个H成键。因此，H₂解离为双位点的反应。根据双位点之间的相对距离，表面位点可以大致分为邻、间、对三种，分别简写为ortho-C–C、meta-C–C和para-C–C（图1a）。

它们催化H₂解离的能垒如图1b所示，即其中，para-C–C上H₂解离所需克服的能垒最低，即使如此，其数值也高达2.88 eV。因此，平整石墨烯对H₂的活化能力几乎可以忽略。

图1 （a）平整石墨烯表面C-C位点分类：邻（ortho-）、间（meta-）、对（para-）。（b）不同位点上H₂解离能垒。

通过常用的对石墨烯平面的晶矢进行等比例的缩放（s = a/a₀ = b/b₀）的操作，我们得到了一系列具备不同褶皱程度的石墨烯（Gr-s）。

随后，对得到的一系列模型进行2 ps室温下的NVT平衡，以使原子位置充分弛豫。平衡中后1 ps的石墨烯厚度如图2a所示，结果表明褶皱越大，其厚度的波动范围也越大。

随后，对得到的Gr-s进行结构优化，这些Gr-s中C–C对距离分布如图2b所示，结构如图2c所示。我们可以从结果中得知，褶皱程度更高的石墨烯中C–C的分布更宽，而平整的石墨烯的距离分布则比较集中，即差异性比较小。

图2（a）不同Gr-s在NVT平衡中厚度的变化。（b）不同Gr-s中C–C对分布函数。（c）高度着色的Gr-s结构。

随后，我们分析了Gr-0.95的电子结构。与平整石墨烯（Gr-1.00）上所有C都几乎等价的情况不同，Gr-0.95中的C原子价态差异性明显（图3a）。同时，计算的Gr-0.95的福井函数（图3b）在表面的分布不再均一，不同C原子之间存在差异，表明这些C得失电子能力的差异化，该函数此前在周期性体系中的应用以及合理性已多有报道。

整体而言，H₂解离的为整个步骤的决速步（图3c）。当H₂物理吸附在具有相反电荷的两个C位点上后，两个C位点都带负电（图3d），而H₂则表现出相反的电荷，即H₂已经发生极化。随着C、H距离的缩短以及H、H距离的边长，C、C之间差异变大，且表现出相反的电荷，H–H物种也继续保持极化的状态，直至H–H断裂以及C–H键形成。结果说明H₂与两个C位点之间的相互作用并非等价的，使得H₂发生极化、活化。

随后，我们对比了其它Gr-s催化H₂解离所需克服的能垒（图3e）。结果表明，Gr-s的厚度与s呈反比，H₂解离能垒与s呈正比，从s = 1.00到s = 0.95，能垒迅速从2.88降到1.74 eV。当s继续减小，下降趋势变缓，s = 0.92时，H₂解离能垒为1.35 eV。这些结果说明褶皱能显著地促进H₂解离。另外，产生的氢物种可以用于乙炔加氢（图3c）。

图3（a）Gr-0.95的原子价态分布。（b）Gr-0.95的f⁺(r), f⁻(r)和Δf(r)的截面图。（c）Gr-0.95催化H₂解离以及乙炔加氢的路径图。（d）H₂活化过程中C和H原子的价态变化。（e）缩放系数s与厚度以及H₂解离能垒的关系。

随后，我们通过slow-growth方法对H₂活化过程进行采样分析。从图4a-b中Gr-1.00和Gr-0.95上H₂活化过程中C、C和H、H之间的距离变化我们可以得知，H₂活化过程中H–H键的断裂伴随着C–C距离的变短。

而相对而言，较平的Gr-1.00的C–C距离降低的程度更低，而Gr-0.95中C–C距离则更灵活，可以减少到更低。这些结果也体现了褶皱结构的灵活性，为H₂解离提供更近的C–C对。

随后，我们横向对比了Gr-s上活化H₂中石墨烯厚度的变化（图4c）后发现，随着H₂活化进行，厚度较小的Gr-s变化较大，而本身具备一定褶皱的Gr-0.97，Gr-0.96和Gr-0.95的厚度变化则比较大，进一步佐证了此前的结论。

而电子结构分析则说明了H₂活化过程中C位点之间电子价态的差异，以及H₂物种在活化过程中的极化状态。这些结果与静态计算得到的结论一致。

图4 （a）Gr-1.00和（b）Gr-0.95催化H₂解离过程中H–H和C–C距离的变化。（c）Gr-s催化H₂解离过程中石墨烯厚度的变化。（d）Gr-0.95催化H₂解离过程中Q_C1 – Q_C2以及Q_H1×Q_H2。

总体而言，褶皱的存在可以造成石墨烯不同C原子电子结构的差异化，并提供具备不同得失电子能力以及不同距离的C–C FLP位点作为H₂解离位点。

因此，本研究的结果不仅揭示了表面形貌对二维材料催化性能的重要影响，而且丰富了人们对FLP概念/内涵的理解。

Corrugation-Induced Active Sites on Pristine Graphene for H2 Activation, ACS Catal. 2022, 12, 14601–14608, DOI: 10.1021/acscatal.2c04894

https://pubs.acs.org/doi/10.1021/acscatal.2c04894

课题组介绍：

林森，教授，博导。2011年在南京大学获得博士学位后加入福州大学，目前在能源与环境光催化国家重点实验室团队，任基础理论课题组组长，2021年入选福建省“雏鹰计划”青年拔尖人才。从事多相催化理论研究，发表论文100余篇，其中第一作者/通讯作者在Nat Catal、Nat Commun、J Am Chem Soc、Angew Chem Int Ed、Sci Adv、Chem、ACS Catal等重要期刊发表论文80余篇，参与编著英文著作2部，所有论文被他引5000余次，主持项目10余项，其中包括国家自然科学基金项目3项，福建省自然科学基金重点项目1项。

课题组网站：https://www.x-mol.com/groups/Sen_Lin