从单轴到双轴：电势驱动下 IrN₄ 活性位点配位构型的动态演变与 C-N 偶联前体锁定(Nano Research 2026)

研究背景

甲烷（CH₄）作为天然气的主要成分，是极其丰富的碳基资源。然而，由于甲烷具有极其牢固的 C-H 键（键能高）以及非极性特征，其在温和条件下的选择性活化与高值化利用一直是化学与催化领域极具挑战性的科学难题。传统的甲烷转化工艺往往依赖于高温高压的蒸汽重整或甲醇中介胺化，不仅能量消耗巨大，且伴随着极低的碳利用效率与显著的温室气体排放。近年来，电化学甲烷活化技术重塑了这一领域，实验研究已证实可在温和条件下实现甲烷的氨化反应生成甲胺，但微观反应机理尚不明确，且传统氧化物催化剂的活性位点结构复杂，深层反应机制与选择性调控规律仍不清晰。石墨烯负载的氮掺杂碳单原子催化剂（M-N-C）凭借其明确且可调控的活性中心，成为阐明反应机理的理想平台。其中，具备 IrN₄ 活性中心的 Ir-N-C 单原子催化剂因 5d 过渡金属 Ir 具备空间延伸且能量灵活的 d 轨道，且对氧亲和力相对较弱，在选择性活化甲烷及抑制析氧副反应（OER）方面展现出独特潜力。基于此，清华大学化学系肖海教授团队在期刊《Nano Research》上发表了题为“Potential-induced bis-axial coordination in Ir-N-C single-atom catalyst delivers selective C-N coupling for efficient methane amination”的最新研究论文。作者采用巨正则系综密度泛函理论（GCE-DFT）结合周期性能量分解分析（PEDA/EDA-NOCV）及活化应变模型（ASM
），系统揭示了阳极电势下 IrN₄ 单原子活性位点的轴向配位动态演变机制及其对 C-N 偶联选择性的调控本质。

研究内容

作者利用恒定电位下的巨正则系综密度泛函理论(grand-canonical ensemble density functional theory, GCE-DFT)计算，系统解析了工况条件下 IrN₄ 活性中心随电势和 pH 变化的结构演化规律(图1)。在 U = 1.20 V v.s. RHE 和 pH = 10 的典型阳极环境中，IrN₄ 位点展现出独特的化学选择性：CH₄ 优先吸附活化脱氢生成单轴配位的 *CH₃ 化合物，随后捕获氨分子并经进一步脱氢过程，最终演化为双轴配位的 *CH₂-*NH₂ 构型。表面 Pourbaix 图（图1（c））进一步证实，该构型在较宽的电势窗口内为热力学优势相。这一电势诱导的双轴结构具有双重关键优势：一方面通过饱和占据轴向位点有效阻断了 *OH/*O 的吸附，从而抑制了竞争性的析氧反应（oxygen evolution reaction, OER）；另一方面将含碳中间体锁定在具有最佳反应性的
*CH₂状态，为后续高效的选择性 C–N 偶联奠定了结构基础。

在阳极电势条件（U = 1.20 V vs RHE、pH = 10）下，作者结合态密度（Projected density of states, PDOS
）、积分晶体轨道哈密顿布居（integrated crystal orbital Hamilton population, ICOHP）、Bader 电荷以及周期性能量分解分析（Periodic energy decomposition analysis, pEDA），在电子结构层面深入揭示了 IrN₄ 位点轴向配体选择性的起源。如图2，分析表明，Ir-A 键（A = C/N/O）主要由 Ir 的 5d_z² 轨道与吸附物 2p_z 轨道之间的 σ 型相互作用驱动。其中，*OH 的吸附能最弱，归因于氧原子高电负性导致 O 的 2p_z 与 Ir 的 5d_z² 轨道能量严重失配，造成 Ir-O 键共价性降低、离子极化显著（Bader 电荷显示 Ir 正电荷升至 +1.38 |e|，图3），使其在热力学上竞争不过结合力更强、共价性更高的 Ir-C 和 Ir-N 键。在 *NH₂ 与
*CH₃ 配位的竞争中，尽管 PDOS 与 ICOHP 显示 N 2p_z 与 Ir 5d_z² 轨道能量匹配更好且键长更短、结合更强，但热力学上首个轴向配位却优先选择 *CH₃。

为了解释这一现象，作者使用 Amsterdam Modeling Suite（AMS2022.103）软件中的 BAND 模块的
pEDA 方法，将总相互作用能（E_int）精细拆解为泡利排斥（E_Pauli）、静电作用（E_elstat）及轨道相互作用（E_orb），并将轨道相互作用进一步细分为动能项 T_orb、极化静电项 Elst_orb 与交换关联项 XC_orb。pEDA 的详细分析揭示：*NH₂ 上的孤对电子在靠近 Ir 中心时产生了极其强烈的近程泡利排斥（E_Pauli = 291.15 kcal/mol）与极化惩罚（Elst_orb = 1435.78 kcal/mol），这完全抵消了其强轨道相互作用（ T_orb = -1380.77 kcal/mol）所带来的热力学收益；反观 *CH₃，凭借其更弥散的 σ 键电子分布，不仅最小化了泡利排斥（E_Pauli = 230.30 kcal/mol），还在轨道相互作用与电荷转移之间达到了最佳平衡，从而证明首步 PCET 的轴向选择性是由“内在成键”与“泡利排斥最小化”共同调控的。

进一步对双轴配位构型的分析表明，异核 *CH₃-*NH₂ 构型通过 Ir (5d_z²) – N (2p_z) 的 σ 相互作用与 Ir (5d_z²) – N (2p_y) 的 π 键贡献构建了独特的 σ-π 协同作用，辅以适度的反馈键，使其成为全局热力学极小值并稳定了反应中间体(图4)。

针对双轴向 *CH₂-*NH₂ 构型的催化动力学活性，研究证实 C-N 偶联路径具有最低的活化能垒（0.58 eV），确立了其相对于 C-C 和 C-O 偶联的动力学主导地位（图5（a））。反应遵循协同机理，通过活化应变模型（activation strain model, ASM）分析揭示（图5（b）），该选择性优势源于 NH₃ 的高能孤对电子与亲电性卡宾碳（*CH₂）空轨道之间理想的路易斯酸碱匹配，这种相互作用对过渡态的稳定作用远强于其他竞争性排斥；反之，H₂O 因氧的高电负性导致其亲核性不足，而 CH₄ 则因缺乏孤对电子且需克服较大的泡利排斥而难以发生偶联反应，使得催化剂能够精准屏蔽弱亲核试剂，将反应流高效定向至 C–N 偶联生成甲基胺，脱附后经连续质子-电子转移（PCET）步骤重新闭环催化循环，完美实现了甲烷胺化的高活性与高选择性。

总结与展望

本研究利用恒电势巨正则系综密度泛函理论（GCE-DFT）计算揭示了 IrN₄ 单原子催化剂的电势诱导轴向配位演变机制，并通过周期性能量分解分析（pEDA）深入剖析了反应路径的活性根源。利用 pEDA 方法将相互作用能精细拆解，发现 *NH₂ 因强泡利排斥与极化惩罚抵消了轨道优势，而 *CH₃ 凭借弥散电子分布有效最小化排斥项并实现最优成键模式。该分析从电子层面确证了首步 PCET 的轴向选择性由“泡利排斥最小化”与“内在成键”共同调控。进一步分析表明，异核 *CH₃–*NH₂ 构型凭借独特的 σ-π 协同作用成为全局热力学极小值。动力学研究证实，双轴向 *CH₂–*NH₂ 构型介导的 C–N 偶联具有最低反应能垒，ASM 分析揭示其优势源于 NH₃ 与亲电卡宾碳间理想的路易斯酸碱匹配，有效克服了 H₂O 亲核性弱与 CH₄ 泡利排斥强的问题。这些发现不仅确立了“电势诱导轴向配位”作为实现烃类精准转化的设计新策略，更突显了 EDA 方法在解析复杂电催化机理、指导高选择性催化剂理性设计中的关键价值。

参考文献

• Li W, Wang Z, Li J, et al. Potential-induced bis-axial coordination in Ir-NC single-atom catalyst delivers selective CN coupling for efficient methane amination, Nano Research, 2026, 19, 94908484, DOI: 10.26599/NR.2025.94906990

感谢胡钧员老师供稿。