NiFe 合金做 NRR 时，很多人会先问 N₂ 吸附能和 *NNH 形成自由能。但如果要解释为什么某个 Ni-rich、Fe-rich 或 Ni-Fe 邻位更容易活化 N₂，只看一个吸附能往往不够。功函数和静电势可以提供另一层信息：电子从表面逸出有多难，界面附近是否存在有利于电荷重排的电势梯度。功函数和静电势只能辅助解释 N₂ 活化，不能替代吸附自由能和反应路径自由能。
图1 展示了功函数对齐、电荷转移示意和静电势分布怎样放在同一组 DFT 结果里。虽然案例体系不是 NiFe NRR，但图中的后处理逻辑适合迁移到 NiFe 合金：先计算干净表面的功函数，再比较 N₂ 吸附前后的电势与 Bader 电荷，最后回到 *NNH 形成自由能判断是否真的有利。
图1 DOI: 10.1021/acsomega.6c00736
它描述电子逸出边界，不直接给活性排名
功函数 Φ 描述电子从费米能级移到真空能级需要跨过的能量。对 NiFe 合金 slab 来说，它更像一个表面边界指标：表层 Ni/Fe 比例、终止面、吸附物偶极、真空层厚度和偶极修正都会影响它。功函数较低通常提示表面更容易供电子，但不等于 N₂ 一定更容易还原，因为 N₂ 活化还需要合适的吸附构型和后续加氢路径。
图2 给出了用 UPS 和功函数读数比较表面电子性质的例子。迁移到 VASP 时，它提醒我们不要把一个 Φ 数值直接写成活性结论，而要先说明费米能级、真空平台和表面偶极是怎样取得的。
图2 DOI: 10.1038/s41467-025-65273-2
只有在模型边界统一以后，下面这个定义才有可比性。也就是说，Ni-rich 和 Fe-rich 表面要用相同层数、相同真空层、相同固定层和相同电势读取方式；非对称 slab 还要说明 LDIPOL、IDIPOL 与真空平台取值位置。
式(1) Φ = Evac – EF
这里 Evac 是真空平台能级，EF 是费米能级。实操时要从 LOCPOT 的平面平均静电势中取稳定真空平台，再从 OUTCAR 或 DOS 结果中读取费米能级。若 slab 两侧不对称，需要检查偶极修正是否一致。
平面平均看边界，局域分布看位点
NiFe slab 的静电势分析至少分两层。第一层是沿真空方向的平面平均电势，用来判断真空平台、表面偶极和两侧不对称程度。第二层是局域电势或电荷密度差，用来判断 N₂ 附近是否出现电子累积和耗尽。平面平均电势解决边界条件，局域电势才接近位点解释，二者不要混用。
图3 把电荷、自由能和差分电荷密度放在一组结果里。NiFe NRR 可以借这个结构组织证据：先列 Ni、Fe 及 N₂ 的 Bader 电荷变化，再展示 N₂ 端位或侧位附近的差分电荷密度，最后给出 *NNH 形成自由能。
图3 DOI: 10.1002/smll.202511734
如果重点是吸附前后表面边界的变化，就比较 ΔΦ；如果重点是某个 Ni-Fe 邻位是否向 N₂ 回馈电子，就看差分电荷密度和 N₂ 片段电荷。两类图各回答一个问题，写作时不要互相替代。
式(2) ΔΦ = Φ*N₂/NiFe – ΦNiFe
这个差值表示 N₂ 吸附前后功函数变化。若 ΔΦ 为正，常说明吸附态增加了表面电子逸出难度；若 ΔΦ 为负，常说明吸附诱导的偶极使电子更容易离开表面。这里的符号必须和取面方向一起说明，尤其是非对称 slab。
功函数只能帮你解释趋势
N₂ 活化不能只用“得电子多”来判断。更稳的组合是 N-N 键长、N₂ 片段 Bader 电荷、N₂ 吸附能、*NNH 形成自由能和功函数变化。若某个 Ni-Fe 桥位表现为 Φ 降低、N₂ 得电子、N-N 键拉长，同时 *NNH 形成自由能也下降，才可以说表面电子逸出难度和 N₂ 活化趋势相互支持。若只有功函数降低而 *NNH 仍然很难形成，就不能写成活性提升。
图4 展示了 NRR 与 HER 竞争的质子-电子转移关系。NiFe 表面如果太容易供电子，也可能先稳定 *H，而不是优先稳定 *N₂ 或 *NNH。因此功函数分析要和 ΔG*H 放在一起看：一个表面电子更活跃，不等于对 NRR 更选择性。
图4 DOI: 10.1021/acs.chemrev.4c00155
在写自由能时，*NNH 通常比 N₂ 吸附能更接近起步难度。下面这个式子是简化参照，真正用于电化学条件时还要加入电位项、pH 项和振动校正。
式(3) ΔG*NNH = G*NNH – G*N₂ – 1/2GH₂
这个表达式用来衡量 N₂ 第一步加氢难度。实际电化学条件下还要加入电位项、pH 项和振动校正；但无论采用哪种参照，*NNH 形成自由能都比功函数更接近 NRR 起步难度。
先定同一位点族，再扩展到组成差异
建议对每个 NiFe 位点做三组静态计算。第一组是干净 slab，用来读取 Φ、平面平均静电势和表层原子 Bader 电荷。第二组是 *N₂/NiFe，用来读取 ΔΦ、N₂ 电荷、N-N 键长和差分电荷密度。第三组是 *NNH 及后续中间体，用来判断自由能路径。功函数数据应和具体位点绑定，不能把一个 slab 的 Φ 套到所有吸附构型上。
如果模型数量有限，可以先做同一合金表面的位点对照，而不是急着换很多组分。比如同一个 NiFe(111) slab 上同时比较 Ni-top、Fe-top、Ni-Fe bridge 和缺陷邻位，观察 Φ、ΔΦ、N₂ 电荷和 ΔG*NNH 是否同向变化。这样得到的是位点效应，不会被误读成所有 NiFe 合金组成的普遍规律。
图5 说明 N₂ 与 *H 竞争可以用描述符区分吸附区域。NiFe 合金中，如果 Fe-rich 位点稳定 N₂，但邻近 Ni 位点稳定 *H，就需要把两种位点放进同一张对照表，而不是只挑 N₂ 最稳的构型。表格中至少保留 Φ、ΔΦ、Eads(N₂)、N-N 键长、ΔG*NNH 和 ΔG*H。
图5 DOI: 10.1021/acs.chemrev.4c00155
NiFe 功函数和静电势的结论建议避免绝对化。可以写“较低功函数提示该位点具备更强供电子倾向，并与 N₂ 吸附后电荷转移和 N-N 键拉长相一致”；不要写“功函数降低证明 NRR 活性更高”。若自由能路径不支持，功函数只能说明电子结构趋势，不能替代热力学筛选。
图6 把 NRR 性能、产氨速率和法拉第效率放在同一组结果里。它适合放在结论写法这一节，是因为计算稿最后也要有类似逻辑：电子结构证据解释“为什么可能有利”，自由能路径判断“是否走得通”，而性能或选择性语境负责约束“这种解释有没有现实意义”。
图6 DOI: 10.1021/acs.chemrev.4c00155
还要注意模型边界。不同 NiFe 表面厚度、真空层、磁性初值和偶极修正都会改变 Φ。比较 Ni-rich 和 Fe-rich 表面时，晶胞、层数、固定层、K 点和电势读取方式必须一致。若采用非对称 slab，建议同时报告 LDIPOL、IDIPOL 和真空平台取值位置。
如果要把多个 NiFe 组成放在同一张图里比较，建议先把结构分组。Ni-rich 表面、Fe-rich 表面、Ni-Fe 邻位、台阶位和缺陷位不应直接混成一条散点线。更稳的做法是每组内部比较 Φ 与 ΔG*NNH 的关系，再看不同组之间是否存在共同趋势。若只有某一组呈现相关性，就把结论限定在该结构族中。
静电势图也要避免过度解读颜色。局域电势低并不自动代表 N₂ 会被活化，它只说明电子在空间中的势能环境发生变化。写图注时可以说明取的是平面平均、宏观平均还是三维局域电势；写正文时再把电势变化与 N₂ 端位、侧位、桥位构型对应起来。电势图要服务于构型解释，而不是独立制造机理。
对于已经做过自由能路径的体系，功函数最好放在补充判断中：先用台阶图确定限制步，再回头看限制步前后 Φ 是否发生一致变化。若限制步是 *NNH 形成，关注 *N₂ 到 *NNH 的 ΔΦ；若限制步是 NH₃ 脱附，功函数变化就未必是主要矛盾。这样组织会让电子结构证据和反应步骤一一对应。
NiFe 合金功函数和静电势适合解释表面供电子倾向、界面偶极和 N₂ 吸附诱导的电荷重排，但它们不是 NRR 活性的直接判据。稳妥做法是把 Φ、ΔΦ、LOCPOT 平面平均电势、N₂ Bader 电荷、N-N 键长、ΔG*NNH 和 ΔG*H 放进同一张表。只有电子结构趋势和自由能路径同时支持时，才能说某个 NiFe 位点更有利于 N₂ 活化。
