自旋态变化会怎样影响反应能垒？

引言

很多反应能垒算不顺，不是 NEB 插点不够，也不是收敛标准太松，而是路径两端的自旋态本来就不在同一条电子态曲线上。尤其遇到 O₂、*OOH、过渡金属单原子或氧化物表面时，如果默认非自旋计算，结果看起来很整齐，物理图像却可能已经错位。

自旋问题麻烦在于它不总是以报错形式出现。计算会正常结束，能垒也会给出一个数，但磁矩、占据和中间体电子结构可能在路径上突然跳变。自旋态没有检查清楚，反应能垒就只能算半个结果。

为什么O₂一上来就不能随便算？

气相 O₂ 是三重态，不是普通闭壳层分子

O₂ 的基态带有两个未成对电子。如果把它当成非自旋体系，参照能就会偏掉；后面所有涉及 O₂、*OOH、*O 或 *OH 的反应自由能都会被这个参照拖着走。OER、ORR 里很多看似“过电位很漂亮”的结果，第一步就可能输在 O₂ 参照上。

在 VASP 里，气相 O₂ 至少要开 ISPIN=2，并给合理 MAGMOM。弛豫后还要看 OUTCAR 或磁矩输出，确认总磁矩和 O 原子局域磁矩没有异常塌掉。O₂ 参照态错了，后面的自由能曲线再平滑也不可靠。

*OOH 是自旋问题和吸附问题的交界处

*OOH 一端连着表面，另一端保留过氧特征，电子结构不像 *OH 那么简单。它可能从表面金属位点取电子，也可能通过 O-O 片段保留一定自由基特征。若金属位点本身有磁矩，*OOH 吸附后常会诱导局域磁矩变化。

因此，不要只比较 *OOH 的吸附能。更好的习惯是同时记录总磁矩、关键原子的局域磁矩、O-O 键长和金属–氧键长。如果 *OOH 的低能构型伴随磁矩重排，机理里就要把自旋态变化写进去，而不是把它当成普通构型差异。

图1 OER中单重态OH⁻与三重态O₂的自旋电子转移示意。 DOI：10.1021/prechem.3c00059

过渡金属位点为什么特别容易出问题？

d 电子占据会跟吸附过程一起变

Fe、Co、Ni、Mn 等位点常有多个接近的自旋态。吸附前是高自旋，吸附后变成中自旋或低自旋，并不罕见。这个变化会改变 M-O 键强度，也会改变过渡态附近的电子转移方式。

如果只用一个初始 MAGMOM 跑到底，VASP 可能收敛到某个局部自旋解。这个解不一定是最低能，也不一定对应真实反应路径。过渡金属位点的 MAGMOM 不是输入一次就结束，而是需要作为构型筛选的一部分。

ΔE^‡ = E_TS − E_IS

当 IS、TS、FS 的自旋态不同，ΔE^‡ 代表的就不只是几何势垒，还混入了电子态切换代价。

看磁矩，不只看总能

同一个结构可以从不同初始磁矩出发，比较最终能量和磁矩分布。若多个磁态能量差小于几十 meV，就不宜轻易说某一个路径绝对占优。对于单原子位点，还要看金属中心和吸附氧上的局域磁矩是否合理。

PDOS 也很有用。若自旋上、下通道在费米能级附近变化明显，说明吸附过程改变了电子占据；若过渡态附近出现突然的自旋翻转，需要考虑路径是否跨了电子态。能垒异常偏低或偏高时，先查磁矩，往往比盲目加插点更快。

图2 自旋催化与磁性催化的概念示意，可用于区分自旋态和外场效应。 DOI：10.1021/prechem.3c00059

VASP里应该怎样做？

先把每个关键态单独跑稳

1.气相 O₂、吸附前表面、*OOH、*O、*OH、过渡态初猜都使用 ISPIN=2。

2.对过渡金属位点设置两到三组 MAGMOM，例如高自旋、中自旋、低自旋初猜，比较最终能量。

3.记录每个关键态的总磁矩、金属局域磁矩、吸附 O 局域磁矩和关键键长。

4.NEB 路径不要只看能量曲线，也要检查每个 image 的磁矩是否连续。

如果某个 image 的磁矩突然跳到另一条曲线，能垒点可能不是几何过渡态，而是电子态切换点。此时可以尝试改变初始磁矩、分段优化路径，或至少在论文里说明该路径涉及自旋重排。

什么时候非自旋结果容易成为无效计算？

遇到闭壳层分子、非磁性金属表面和简单吸附物，非自旋结果有时足够。但只要体系包含 O₂、NO、过渡金属未满 d 壳层、氧空位诱导磁矩、自由基中间体，就不能默认关自旋。

判断标准很朴素：如果开启自旋后总能明显降低，或者出现稳定局域磁矩，非自旋结果就只能作为初筛。非自旋计算不是更干净，只是在某些体系里把磁性自由度冻结掉了。

图3 过渡金属氧化物中低自旋、中自旋和高自旋状态示例。 DOI：10.1021/prechem.3c00059

能垒结果怎样写？

把自旋态作为路径属性写出来

写反应路径时，不要只列 IS、TS、FS 的能量。建议同时给出磁矩变化：例如金属中心从 3.1 μ_B 降到 1.8 μ_B，或 *OOH 吸附后 O 端出现局域自旋。这样别人能看出能垒变化是否与电子态重排有关。

如果你比较了多个初始 MAGMOM，也可以简短说明最低能路径来自哪一组磁态初猜。这个信息不占篇幅，却能显著提高可复查性。

不要把所有能垒差都归因于几何结构

两个催化位点的能垒差，可能来自键长、配位数，也可能来自自旋态不同。若位点 A 在过渡态保持高自旋，位点 B 发生低自旋重排，直接说“A 的几何结构更优”就太简单了。

更稳的写法是把结构因素和电子态因素并列：位点几何决定了吸附姿态，自旋态变化调节了电子转移与过渡态稳定。对含过渡金属和含氧中间体的能垒，磁矩是必须随手记录的证据。

图4 FeN₄/FeN₄O₂位点中 Fe-O 键长与磁矩变化示例。 DOI：10.1021/prechem.3c00059

本章要点总结

要点一：O₂ 是三重态，涉及 OER/ORR 的气相参照不能默认非自旋。

要点二：*OOH 同时牵涉吸附构型和自旋重排，不能只看吸附能。

要点三：过渡金属位点要用多组 MAGMOM 初猜筛查，NEB 路径也要检查 image 磁矩连续性。

要点四：能垒解释要把几何变化和电子态变化分开，避免把自旋贡献误写成结构贡献。