反应路径中参照态怎么选？

引言

自由能路径里最隐蔽的错误，往往不是中间体结构，而是参照态。H₂、H₂O、N₂、NO₃⁻、CO₂ 各自有不同自旋、溶剂、电荷和热力学校正问题。参照态一旦不统一，后面的每一个台阶都会被同一个错误拖着走。

气相参照最容易错在哪里？

H2 和 N2 为什么不能随手算？

H₂ 常用于 CHE，把 H⁺ + e⁻ 和 1/2H₂ 联系起来。N₂ 三键很强，盒子大小、自旋、能量收敛和热校正都会影响 NRR 台阶。气相分子要单独建大盒子，不能从吸附态里截出来当参照。

H₂ 的问题常被低估。它看起来只是一个双原子分子，但在 CHE 里相当于把 H⁺ + e^– 的化学势折算到 1/2H₂ 上，所以 H₂ 的能量、零点能和熵修正会直接影响每一个质子电子转移步骤。

N₂ 则是另一类麻烦。它三键很强，能量收敛通常不难，但一旦拿它去做 NRR 或含氮物种参照，后面每一个 *NNH、*NH、*NH₂ 相关台阶都会继承这个零点。气相 N₂ 不要从吸附结构里截取，也不要和不同设置下的气相目录混用。

气相参照要独立收敛，不能借用吸附结构里的片段能量。

实操时建议给气相分子单独建目录，固定一个足够大的盒子、相同泛函和截断能，并把 OUTCAR、OSZICAR、振动频率和热修正来源一起留存。以后自由能表出了异常，才能回头查到参照态到底来自哪里。

式(1)  H⁺ + e⁻ ↔ 1/2H₂

O2、CO2 这类分子有什么额外问题？

O₂ 是三重态，CO₂ 是线性分子，二者对自旋、熵和振动修正都敏感。CO₂RR 路径还常把 CO₂、CO、HCOOH、CH₃OH 等多个产物参照放在一起，电子数和 H 数必须严格配平。

O₂ 最容易被忘的是自旋。若把三重态 O₂ 当成普通闭壳层分子处理，OER、ORR 中含氧中间体的自由能零点会整体偏掉。很多研究会避开直接使用 O₂ 总能，改用 H₂O/H₂ 组合或经验修正，这些选择都要在方法里讲清楚。

CO₂ 参照看似简单，但 CO₂RR 往往会同时牵涉 CO、HCOOH、CH₃OH、CH₄ 等产物。产物越多，越要先把 C、O、H 和电子数配平，否则同一张路径图里可能混进不同零点。

图1展示 Sabatier 最优区间与反应自由能关系，适合提醒参照态要和每一步反应式逐项对应。

如果图里的某一步突然比相邻步骤低很多，不要先急着解释成材料优势。先查这一步用到的气相分子、H 数、电子数和热修正是否与前后步骤一致，往往比盯着吸附构型更快找到问题。

图1  Sabatier 最优区间与自由能路径关系。DOI: 10.1002/anie.3566312。

水和带电物种为什么更麻烦？

H₂O 参照和 pH 怎样连起来？

在碱性 HER、OER 或 ORR 里，H₂O、OH⁻、H⁺ 和 e⁻ 的关系要通过 CHE 或显式反应式处理。不能一会儿用酸性写法，一会儿用碱性写法。

pH 的作用不是在自由能图末尾随便加一个修正项。酸性写法里常用 H⁺ + e^–，碱性写法里则经常把 H₂O 和 OH^– 放进反应式。两种写法可以互相换算，但不能在同一条路径里随手混用。

尤其是碱性 HER 和 OER，水既可能是反应物，也可能参与稳定界面氢键网络。若模型只算了干净 slab 上的单个吸附物，方法里就要说明这是热力学近似，而不是把水层效应也一并包含了。

图2可以用来说明反应路径中不同中间体的自由能变化；只要牵涉质子电子转移，就要把电位和 pH 条件写清楚。

这类图的价值不只是看哪一步最高，还要看每一步使用了多少电子和质子。只要电位修正或 pH 修正漏了一步，自由能线条就会看起来很顺，但化学计量已经不对。

式(2)  ΔG(U) = ΔG(0) – neU

图2  反应自由能与吸附强度关系。DOI: 10.1002/anie.3566312。

NO₃⁻ 为什么不能直接当普通气相分子？

NO₃⁻ 是溶液离子，直接放到真空盒子里算会牵涉带电体系、背景电荷和溶剂稳定问题。NO₃RR 常用实验化学势、溶液参照或约定自由能修正，必须在方法里说明。

NO₃^– 不能像 CO₂、N₂ 那样简单放进真空盒子里直接取总能。带电体系会引入均匀背景电荷，真空层、盒子大小和偶极修正都会影响能量；更关键的是，真实 NO₃^– 在溶液中被水合稳定。

因此 NO₃RR 更常见的做法是使用约定化学势、实验自由能或文献修正，把溶液物种和表面吸附态接到同一套热力学标尺上。这个标尺不交代清楚，后面对 *NO₃、*NO₂、*NO、*NH_x 的讨论都会缺一个可信零点。

图3显示含氮反应中 Sabatier 最优区间的概念示意，带电溶液物种的参照态必须交代来源，不能和中性气相分子同等处理。

看到含氧氮物种路径时，最好顺手检查每一步脱氧、加氢和质子电子转移是否写完整。NO₃RR 的路径分支很多，参照态一乱，最容易把不同加氢方向算成同一条线。

图3  含氮反应中 Sabatier 最优区间示意。DOI: 10.1002/anie.3566312。

不同反应体系的参照态到底差在哪里？

HER、OER和ORR为什么不能只换一个分子名？

HER 最常见的参照是 CHE，很多时候只需要把 H⁺ + e^– 与 1/2H₂ 对上；但到了 OER 和 ORR，参照态就会牵涉 H₂O、OH^–、O₂ 以及电位和 pH。看起来只是多了氧物种，实际是每一步的 H、O、电子来源都变复杂了。

所以同一张自由能图里，不能把酸性 OER 的写法、碱性 ORR 的写法和 HER 的 CHE 写法随手拼在一起。哪怕每个中间体总能都算得很收敛，只要参照规则混了，台阶高度就会被人为改写。

CO₂RR 和 NO₃RR 为什么更容易乱？

CO₂RR 的麻烦在于碳产物太多，*CO、HCOOH、CH₃OH、CH₄ 对应的氢数、电子数和氧去向都不同。NO₃RR 的麻烦则在于 NO₃^– 本身是溶液离子，不能像中性气相分子那样直接丢进真空盒子后拿总能当标准参照。

写这两类路径时，建议先按产物分支整理“碳/氮来源、氧来源、氢来源、电子数”四列。这个底稿不一定放进正文，但能帮你提前发现某一步少了一个 H₂O，或者把 NO₃^– 当成了普通中性分子。

自由能图画出来后先复查哪几类异常？

整条线整体偏移，通常不是材料突然变强

如果某条路径上所有含 H 的步骤都一起变低，先查 H₂、H₂O 或 CHE 写法；如果所有含 O 的物种都异常，优先查 O₂、H₂O/OH^– 的参照；如果含 N 物种整体奇怪，就回到 N₂、NO、NO₃^– 的处理方式。参照态错误常常不是让一个点离谱，而是让一组点一起漂移。

图4可以作为简单参照：HER 里 H* 自由能看似只牵涉一个吸附态，实际也离不开 H₂ 和 CHE 零点。

图4  H 吸附自由能与覆盖度关系示例。DOI: 10.1021/acsomega.6c03086。

单个中间体离谱时，先看结构还是参照？

如果只有某一个吸附态突然异常，通常要同时看两件事：一是 CONTCAR 里中间体有没有断键、迁移或变成另一种构型；二是这个中间体对应的反应式有没有少配一个 H⁺ + e^–、H₂O 或 OH^–。结构跑偏和参照配错都会制造“漂亮但不真实”的台阶。

最后再决定是否重算。若是结构问题，就回到初始构型和约束；若是参照问题，通常需要统一整条路径的自由能表，而不是只修一个数字。自由能图的可靠性，更多来自这类复查，而不是图画得多顺滑。

本章要点总结

要点1：参照态错误会系统性影响整条自由能路径。

要点2：气相分子要独立收敛，带电溶液物种要说明化学势来源。

要点3：酸性和碱性写法不能混用。

要点4：每一步都应检查原子数、电荷数和电子数守恒。