DFT指导的N4V2-石墨烯负载3d过渡金属单原子催化剂用于可持续电合成氨

《Nanoscale》：DFT-guided design of 3d transition-metal single-atom catalysts on N4V2-graphene for sustainable ammonia electrosynthesisOpen Access

　　

电化学氮还原反应（Nitrogen Reduction Reaction，NRR）是能源密集型的Haber–Bosch工艺在环境条件下绿色合成氨的一种可持续替代方案。然而，由于对底层电子结构–活性关系的理解有限，合理的催化剂设计仍然具有挑战性。在此，研究人员采用第一性原理密度泛函理论（Density Functional Theory，DFT）计算，系统研究了一系列石墨烯负载的3d过渡金属（Transition Metal，TM = Sc， Ti， V， Cr， Mn， Fe， Co， Ni， Cu）。每种金属均嵌入氮掺杂石墨烯模型载体上的双碳空位（V₂）中，形成TM–N₄V₂构型。N₄V₂配位环境不仅锚定了孤立的金属原子，还通过配体场效应和电荷再分布调控了其电子结构。研究人员的分析揭示了火山型活性趋势，其中Cr在所研究系列中表现出最低的极限电位（?0.65 V）。然而，其略微吸热的N₂吸附表明活性位点上的反应物可用性有限，说明催化性能不能仅凭极限电位来评估。相反，最优催化剂必须在有利的反应能量学和充分的N₂吸附之间取得平衡。在该框架内，V–N₄V₂成为一个特别有前景的候选者，它兼具放热的N₂吸附和相对较低的极限电位，从而在反应物可及性与催化活性之间提供了更现实的平衡。虽然观察到d带中心与N₂吸附强度之间存在单调相关性，但仅凭d带中心并不能完全反映反应能量学，后者取决于活性位点由配位诱导的电子结构。这些发现强调了配位诱导的电子效应在决定催化行为中的重要性，并为设计高效的NRR单原子催化剂提供了合理依据。

氨（NH₃）是现代农业肥料和碳中性能源载体的基础化学品，但其工业合成主要依赖Haber–Bosch工艺，需要在150–350 atm和350–550 °C的苛刻条件下进行，消耗全球1%–2%的能源供应，并产生大量CO₂排放。电化学氮还原反应（Nitrogen Reduction Reaction，NRR）能够在常温常压下利用可再生能源将N₂和水转化为NH₃，是传统工艺的可持续替代方案。然而，N≡N三键键能高达945 kJ mol^?1，导致反应动力学缓慢且活化能垒高，同时水溶液中的竞争性析氢反应（Hydrogen Evolution Reaction，HER）会降低法拉第效率和选择性。因此，需要开发能有效活化N₂并抑制HER的高效、稳定、低成本电催化剂。碳基纳米材料如石墨烯具有高导电性和大比表面积，但本征石墨烯对N₂活化呈化学惰性，需要通过杂原子掺杂（如氮掺杂）调控其电子结构。单原子催化剂（Single-Atom Catalysts，SACs）锚定在掺杂碳载体上（如Fe–N₃、Mo–N₁C₂）已展现出潜力，但局部配位环境如何调控不同过渡金属的电子结构，以及完整3d过渡金属系列在固定配位构型下的催化趋势仍缺乏系统理解。

研究人员采用第一性原理密度泛函理论（Density Functional Theory，DFT）计算，系统研究了九种3d过渡金属（TM = Sc、Ti、V、Cr、Mn、Fe、Co、Ni、Cu）嵌入氮掺杂石墨烯双碳空位（N₄V₂）位点（记为TM–N₄V₂）的NRR催化性能。研究发现，活性呈火山型趋势，Cr具有最低极限电位（?0.65 V），但其N₂吸附略吸热，可能限制反应物可用性；而V–N₄V₂兼具放热N₂吸附和相对较低极限电位（?0.87 V），是更平衡的候选物。d带中心与N₂吸附强度呈单调关系，但不能单独预测完整反应能量学。这些发现揭示了配位诱导电子效应对催化行为的决定性作用。

研究人员使用Vienna Ab initio Simulation Package（VASP）进行自旋极化DFT计算，采用投影缀加平面波（Projector Augmented-Wave，PAW）方法处理离子-电子相互作用，Perdew–Burke–Ernzerhof（PBE）泛函描述交换关联能，并采用DFT-D3（Becke–Johnson阻尼）校正范德华力。平面波截断能设为450 eV，模型为晶格参数12.35 ?的六方超胞，垂直方向真空层20 ?。结构优化和电子结构计算分别采用2×2×1和8×8×1的Monkhorst–Pack k点网格，电子自洽收敛至10^?5 eV，原子力收敛至0.01 eV ?^?1。HSE杂化泛函用于计算投影态密度（Projected Density of States，PDOS）和d带中心。通过计算氢电极（Computational Hydrogen Electrode，CHE）模型获得各基元反应的吉布斯自由能变化，进而确定极限电位（U_L）和过电位。本研究为纯理论计算，未涉及实验样本队列。

**TM–N₄V₂模型体系**：通过几何优化和电子结构计算，确认所有TM原子均稳定嵌入N₄V₂缺陷中心，形成平面配位结构。TM原子作为电子供体向吡啶氮转移电荷，产生配体场效应，改变d轨道占据和能量分布，为后续催化行为提供结构基础。

**N₂在TM–N₄V₂表面的吸附与活化**：计算了端接和侧接两种N₂吸附构型。结果表明，Sc–Co为化学吸附，Ni和Cu为物理吸附。Ti是唯一以侧接模式化学吸附N₂的金属。化学吸附的金属导致N–N键伸长至1.13–1.18 ?，表明N₂被活化。吸附能从Sc的?1.83 eV到Mn的?0.11 eV不等。弱吸附的Ni、Cu难以引发NRR。

**自由能分析**：基于CHE模型计算了完整NRR路径（N₂ + 6H⁺ + 6e^? → 2NH₃）的自由能图。不同金属遵循远端或交替机制，但所有催化剂的N–N键均在第五步氢化时断裂。极限电位分析显示Cr最低（?0.65 V），其电位决定步骤为第一步氢化（N₂ → NNH）。V–N₄V₂的极限电位为?0.87 V，但N₂吸附放热，因此更平衡。Fe、Co等弱结合金属受NNH形成能垒限制，Ti等强结合金属则因中间体过度稳定而中毒。

**催化活性的电子起源**：通过PDOS分析，早期Ti的d态与N₂态强烈杂化，导致强吸附但过度稳定；晚期Ni、Cu杂化弱，吸附弱且活化不足；Cr表现出适度杂化，有利于NNH形成。Mn的d带中心在吸附后上移1.13 eV，但吸附自由能仍为正（+0.39 eV），其低极限电位被吸热吸附抵消。

**d带模型作为催化活性描述符**：将清洁表面的d带中心（ε_d）与ΔG_N2作图，显示总体单调趋势：ε_d越高，吸附越强。但该描述符无法区分Cr、V、Ti（ε_d相近而ΔG_N2差异大），且对于Co、Ni、Cu，ε_d相差2.85 eV但ΔG_N2相同（+0.20 eV）。因此，d带中心仅是近似描述符，完整反应能量学需结合PDOS分析。

讨论部分强调催化活性并非由单一吸附描述符决定，而是由配位控制的电子结构调控整个反应路径。Sabatier原理在TM–N₄V₂系列中体现为吸附强度与反应能量学的平衡。V–N₄V₂是最平衡的候选者，Cr–N₄V₂具有最低极限电位。这些发现为设计下一代单原子催化剂提供了合理原则。

结论部分翻译为：在这项工作中，研究人员利用第一性原理DFT计算系统研究了氮掺杂石墨烯负载的单原子TM催化剂（TM–N₄V₂）对电化学NRR的催化活性。分析表明，3d TM系列的催化活性遵循明显的火山型趋势，符合Sabatier原理。虽然Cr表现出最低的极限电位（?0.65 V），但其略微吸热的N₂吸附表明在实际条件下活性位点的反应物可用性可能受限。因此，催化性能不能仅凭极限电位评估。有效的催化剂必须在有利的N₂吸附和后续氢化步骤的能量学之间取得平衡。在此背景下，V–N₄V₂成为特别有前景的候选者，它结合了放热的N₂吸附和相对较低的极限电位，从而在反应物活化和催化周转之间提供了更平衡的折衷。此外，研究人员发现清洁表面的d带中心与N₂吸附自由能之间存在单调但近似的关系。d带中心可以捕捉N₂吸附的大致趋势，但应理解为活性位点电子结构的描述符，而非完整反应能量学的预测器。这些发现表明，TM–N₄V₂体系的催化活性不是由单一吸附描述符决定，而是由配位控制的电子结构调控整个反应路径。该研究确定了V–N₄V₂为最平衡的候选者，Cr–N₄V₂为极限电位最低的体系，为开发更高活性和选择性的下一代单原子催化剂提供了合理设计原则。