



M-N4单原子位点常被当成一个平面四配位中心。加入轴向配体 L 后,金属不再只是被四个 N 固定在平面里,而是变成 M-N4-L 这类五配位或准五配位环境。这个改变会影响金属 d 轨道分裂、价电子分布和反应中间体的进入方向。轴向配体不是旁边多了一个基团,而是改变了单原子中心的配位场。
对 OER 来说,*OH、*O、*OOH 都要和金属中心或邻近配位环境相互作用。L 可能削弱 *OH,也可能加强 *O,甚至让 *OOH 从金属顶位转向氢键或邻位结合。路径变化的根源不只是一句“电荷转移”,而是配位场、轨道占据、吸附构型和自由能台阶共同改变。




M-N4和 M-N4-L 不是同一个位点
平面 M-N4中,金属中心主要受四个 N 的平面配位限制,反应中间体常从垂直方向靠近金属。加入 L 后,轴向位置被占据,金属中心的上下方向电子云和空间位阻都会改变。若 L 是 OH、Cl、吡啶、O 或其他基团,其给电子和拉电子能力也不同。
因此比较时不能把 M-N4与 M-N4-L 只当成同一位点的轻微修饰。它们应作为两类活性中心分别优化,再分别放入 *OH、*O、*OOH。轴向配体改变的是金属中心本身,后续吸附自由能要从新位点重新计算。

图1 Fe-N4位点引入轴向 OH 后的反应自由能变化,展示轴向配体可改变关键台阶。DOI:10.1002/advs.202000176
L 会改变 dz²方向的相互作用
轴向配体最直接影响的是垂直于 M-N4平面的轨道相互作用,尤其是 dz²与吸附物 p 轨道的重叠。若 L 与金属成键较强,另一侧中间体靠近时会受到电子排斥或轨道竞争;若 L 主要改变金属电性,则可能让 *OH、*O、*OOH 的吸附强度按不同方向移动。
ΔGOER= ΔE + ΔZPE − TΔS + ΔGU+ ΔGsolv
这也是为什么同一个 L 对不同中间体不一定同向。*O 需要较强 M-O 成键,*OOH 同时受 M-O、O-O 和 O-H 取向影响,*OH 还容易被氢键和轴向竞争影响。配体效应要按中间体分开看,不能只给一个总趋势。

图2 不同轴向配体下 *OH、*O、*OOH 标度关系与火山图,说明 L 会改变吸附自由能窗口。DOI:10.1002/advs.202000176




*OH 常受轴向竞争和氢键影响
*OH 与金属形成 M-OH 键,同时带有 O-H 方向。若轴向 L 已经占据金属一侧,*OH 可能只能从另一侧靠近,或转向与载体边缘形成氢键。若 L 本身也是 OH,则要区分预先存在的配体 OH 和反应路径里的 *OH,二者不能混成同一个物种。
在结果里,*OH 的变化要同时看 M-OH 键长、O-H 方向、L-M-O 角度和 ΔG*OH。如果 ΔG*OH改变很大,但 *OH 最终不是吸在金属顶位,而是被 L 或邻位 N 拉住,结论应写成构型调节,而不是单纯金属电性变化。
对 *OOH 也一样。它不仅需要一个 M-O 键,还要保留 O-O 键和末端 O-H 方向。轴向 L 若能提供氢键或空间牵引,就可能让 *OOH 更容易形成;若 L 挤占反应通道,则可能让 *OOH 远离金属中心。

图3 Fe-N4相关位点在不同条件下的反应路径示意,展示轴向配体可能改变中间体去向。DOI:10.1002/advs.202000176
*O 更敏感于金属氧键强度
*O 通常更强地考验金属中心的氧亲和力。若 L 提供电子,使金属更容易与 O 成键,ΔG*O可能下降;若 L 与金属强配位占用了关键轨道,*O 反而可能变弱。对过渡金属单原子来说,自旋态和 d 轨道占据也会参与调节。
ΔG2= G*O− G*OH− (GH2O− 1/2GH2)
因此,判断 *O 不能只看吸附能一个数。要把 M-O 键长、金属磁矩、Bader 电荷、PDOS 中金属 d 与 O p 态重叠一起看。若 L 改变了 *O 的键强,就很可能把 OER 限速步从 *OH→*O 或 *O→*OOH 之间移动。

图4 不同单原子模型和自由能路径对比,显示配体环境会改变反应台阶。DOI:10.1016/j.isci.2022.104557




路径变化来自自由能台阶重新排序
OER 路径是否改变,最终要看四个质子电子转移台阶的最大 ΔG 是否换了位置。L 可能让 *OH 更容易生成,却让 *O 过强;也可能削弱 *O,使 *OOH 形成成为主要障碍。单看某个中间体吸附增强,很难判断整体路径。
ηOER= max(ΔG1, ΔG2, ΔG3, ΔG4) / e − 1.23 V
若引入 L 后,理论过电位下降,同时限速步从 *OOH 形成转为 *O 形成,说明配体改变了中间体相对能量排序。若过电位几乎不变,但 *OH 和 *O 同时移动,则可能只是沿标度关系平移,并没有真正突破路径限制。

图5 Cl、OH 和吡啶等轴向配体下的自由能路径,说明不同 L 对反应台阶影响不同。DOI:10.1016/j.isci.2022.104557
标度关系没有消失,只是可能被调弯
单原子位点常受 *OH、*O、*OOH 标度关系限制。轴向 L 可以改变金属中心电性和空间构型,让三个中间体的移动幅度不同。如果 *OH 和 *OOH 同步移动,而 *O 单独变化,就可能让火山图位置改变;若三者几乎等比例移动,则改善空间有限。
所以配体效应最好不要只写“优化电子结构”。更清楚的表达是:L 改变 M-N4中心的轴向配位场,使 ΔG*OH、ΔG*O、ΔG*OOH的相对位置发生变化,并进一步改变限速台阶或理论过电位。

图6 轴向 OH 配体引起的金属中心电荷和 DOS 变化,可用于连接配位调节与吸附自由能。DOI:10.1016/j.isci.2022.104557




配体 L 和反应中间体要分清
M-N4-OH 里的 OH 若在反应开始前就存在,它是轴向配体;OER 路径里的 *OH 是中间体。二者能量参照不同,不能在自由能图中随意互换。若计算中 L 脱落、迁移或与中间体合并,就说明初始 M-N4-L 位点没有保持,应重新定义为新的结构。
比较 L 效应时,建议至少保留三类对照:无 L 的 M-N4,带 L 但无中间体的 M-N4-L,带 L 且吸附 *OH/*O/*OOH 的反应态。这样才能分清配体本身结合能、金属中心改变和中间体吸附能。
同一套参照下比较 ΔG
若 L 是外加配体,M-N4-L 的形成能或结合能要单独说明;OER 中间体自由能则应以该带 L 位点作为参照。不要把无 L 位点的 Eslab拿来扣带 L 的吸附态,否则 ΔG 里会混入配体结合能。
SYSTEM = MN4_L_ads_static
ENCUT = 520
EDIFF = 1E-6
IBRION = -1
NSW = 0
LCHARG = .TRUE.
LORBIT = 11
ISPIN = 2
LASPH = .TRUE.
上面片段对应带 L 吸附态的静态单点:结构弛豫完成后,输出电荷密度和投影态密度,用于比较金属中心电荷、磁矩和 d 态变化。若体系含 Fe、Co、Ni、Mn 等过渡金属,还要检查不同初始磁矩是否导向不同电子态。




把配位、吸附和路径串起来
一个清楚的结论应包含三层。第一层:L 与金属形成轴向配位,改变 M-N4的局域配位场和 d 态分布。第二层:这种变化分别影响 *OH、*O、*OOH 的键长、构型和自由能。第三层:自由能台阶重新排序,使限速步或过电位改变。
若只有 Bader 电荷变化,而 *OH、*O、*OOH 自由能没有改善,就不能说路径被优化。若自由能改善但 L 在优化中脱落,也不能继续把结果归因于轴向配体。真正有说服力的结果,是结构保持、电子态变化和自由能台阶变化三者指向同一套解释。
若多个 L 的趋势不同,也不要只按电负性排序。更合适的是分别比较配体结合能、M-L 键长、金属磁矩和三个中间体 ΔG,看哪一个指标真正跟过电位变化同向。




要点一:M-N4-L 与 M-N4是不同配位环境,不能只当成轻微修饰。
要点二:轴向配体主要影响垂直方向轨道相互作用,进而改变 *OH、*O、*OOH 吸附。
要点三:路径变化要看四个 OER 自由能台阶是否重新排序,而不是只看单个吸附能。
要点四:配体 L 与反应中间体 *OH 必须分清,二者不能共用同一能量参照。
要点五:最有力的解释要同时包含配位场变化、电子态变化和自由能路径变化。
