



尖晶石铁氧体常写成 AB₂O₄,但真正进入 ORR 计算后,A 位和 B 位并不是两个普通标签。四面体位、八面体位、表面终止和阳离子反位会同时改变 可吸附位点和反应中间体 的对应关系。先把这个关系厘清,再比较 O₂、*OOH、*O 和 *OH 的自由能,结论才不会被表面命名带偏。




AB₂O₄ 中 A 位和 B 位是什么?
尖晶石结构由氧骨架和两类阳离子位点组成。A 位通常指四面体配位环境,B 位通常指八面体配位环境;在铁氧体中,Fe、Co、Mn、Ni 等阳离子可能分布在不同位点,也可能出现部分反位。对 ORR 来说,阳离子位置会改变表面金属的配位数、价态趋势、局域磁矩和与含氧中间体的键合方式。
图1 给出 XCo₂O₄ 尖晶石结构和阳离子分布的基本对象。A 位和 B 位 不是晶胞中的静态名称,而是会进入吸附构型、磁态初值和电子结构解释的计算变量。若把尖晶石表面只看成一个普通氧化物 slab,就会忽略 A/B 位带来的模型差异。
为什么 ORR 不只放一个 O₂?
ORR 的四电子路径至少涉及 O₂*、*OOH、*O 和 *OH。O₂ 初始吸附常用于判断分子活化,*OOH 常决定第一步质子-电子转移的热力学难度,*O 与 *OH 则会反映后续断键和脱附倾向。四面体位和八面体位对这些物种的稳定性不一定同向变化,因此一个 O₂ 吸附能不能代表整条 ORR 路径。
以同一晶面上的两个金属位为例,若 O₂ 在 B 位吸附更稳定,但 *OH 在同一位过度稳定,自由能台阶可能仍然不理想。相反,某个 A 位对 O₂ 作用较弱,却能给出更平衡的 *OOH 和 *OH 稳定性。ORR 位点判断要比较 一组中间体,而不是只比较 O₂ 单点吸附能。

图1 XCo₂O₄ 尖晶石氧化物的晶体结构、阳离子分布和 ORR 反应对象。图中可看到四面体位、八面体位和氧骨架的关系。DOI: 10.1038/s41467-019-09503-4。




先选晶面还是先选位点?
尖晶石铁氧体常见低指数晶面包括 (001)、(011)、(111) 等,不同切面暴露的 A 位、B 位和桥氧数量不同。VASP 模型应先确定目标晶面和终止面,再讨论 O₂ 或 *OOH 放在哪里。若先随机放吸附物,再回头命名位点,容易把表面终止差异误读为四面体位和八面体位差异。
图2 展示同一尖晶石系列中成分和 ORR 活性的关联。迁移到计算端,同一晶面上比较 A 位与 B 位 是第一类对照,同一阳离子组成下比较不同终止面是第二类对照。只有固定晶面、终止面和阳离子排布后,位点分配问题才有清晰边界。
磁性和 DFT+U 怎样处理?
Fe、Co、Mn 相关尖晶石通常具有明显局域磁矩。结构优化前,应为不同金属位设置合理 MAGMOM 初值,并比较铁磁、反铁磁或亚铁磁初态对能量和吸附结果的影响。若使用 DFT+U,U 值需要来自文献口径或按研究对象说明测试边界,不能在完整表面和吸附态之间混用不同设置。
一个实操情境是,O₂ 吸附后表面金属磁矩发生变化,同时 O₂ 保留部分自旋。此时只检查总能收敛还不够,还要查看 OUTCAR 中局域磁矩、O-O 键长、吸附物电荷和 DOS。若磁态在不同构型之间跳变,吸附能排序就需要结合磁矩连续性重新判断。

图2 不同 XCo₂O₄ 尖晶石样品的结构表征与电化学 ORR 活性比较。图中体现了成分变化和活性差异需要放在同一模型口径下解释。DOI: 10.1038/s41467-019-09503-4。




O₂ 有哪些初始构型?
O₂ 可以采用端位、侧向和桥位构型:端位通常让一个 O 原子靠近金属中心,侧向构型让 O-O 键平行贴近表面,桥位构型则连接两个邻近金属或金属-氧环境。对尖晶石表面,A 位、B 位和桥氧附近都应进入初始构型矩阵,优化后再按最终键长和吸附位点重新归类。
O₂ 吸附能可以用于初筛构型,其中 Eslab-O₂、Eslab 和 EO₂ 必须来自同一泛函、同一自旋设置和同一赝势口径。可连接到 *OOH 的端点 是 O₂ 构型筛选的核心目标,不是只寻找最低吸附能。由于气相 O₂ 的三重态和 GGA 误差会影响绝对值,ORR 中更常把自由能台阶与 CHE 参照合并讨论。
Eads(O₂*) = Eslab-O₂ – Eslab – EO₂ 式(1)
*OOH 为什么是关键中间体?
*OOH 同时包含 O-O 键和 O-H 键,几何上比 O₂ 更能暴露位点选择问题。若 *OOH 通过端位 O 连接金属中心,另一个 O 与 H 指向溶液侧,常对应较典型的缔合 ORR 路径;若 *OOH 在优化中断裂成 *O 和 *OH,则该构型已经偏离目标中间体,需要作为另一条路径单独处理。
图3 和图4 中的 ORR 活性与描述符关联说明,*OOH 的稳定性会改变理论过电位和速控步骤。对比 A 位与 B 位时,应记录 O-O 键长、M-O 键长、H 位置、ΔG*OOH 和后续 *O/*OH 的自由能。这里的判断重点是,哪个位点能让 *OOH 形成而不把后续 *OH 稳定得过强。

图3 XCo₂O₄ 尖晶石氧化物的 ORR 电化学活性和动力学结果。性能差异提示阳离子占位和中间体吸附需要联动分析。DOI: 10.1038/s41467-019-09503-4。




CHE 公式怎样落到 ORR?
在酸性 ORR 的 CHE 近似中,H⁺ + e⁻ 可用 1/2H₂ 代替;碱性口径也可以通过 H₂O/OH⁻ 参照转换。对同一尖晶石表面,四电子 ORR 可写成 O₂ → *OOH → *O → *OH → H₂O 的自由能序列。每一步 ΔG 都需要包含 DFT 能量差、零点能、熵修正、电位项和 pH 项。
自由能修正表达的是某一基元步从 DFT 总能进入电化学尺度的转换。实际比较 A 位与 B 位时,几何端点、自由能台阶和电子态旁证 要同时保留,不能只输出最终过电位。每个中间体所在的几何位点和磁矩状态也应随结果一起记录。
ΔG = ΔEDFT + ΔZPE – TΔS – neU + kBT ln10 × pH 式(2)
哪些对比会改写位点判断?
*OOH 最稳定的位点并不一定给出最低过电位。若同一位点上的 *OH 明显过稳,最后一步形成 H₂O 的自由能会升高;此时位点优势就从“能形成 *OOH”转向“能否完成 *OH 后续处理”。这个情境能提醒我们,A 位或 B 位的优劣不能停在单个中间体的最低能构型上。
还有一类情况来自 O₂ 入口本身。某个位点上的 O₂ 几乎不吸附,但 *OOH 可通过氢键或邻近桥氧稳定,这时 O₂* 不是唯一入口,模型需要加入显式水或邻近质子环境再比较。图4 将描述符与活性关联起来,说明位点优劣不能只从几何图中判断。建议在代表构型上输出 DOSCAR、CHGCAR、AECCAR 和必要的 Bader 电荷,检查金属 d 态、吸附物 O 2p 态和差分电荷密度是否支持自由能趋势。

图4 尖晶石氧化物 ORR 活性与反应中间体吸附强度之间的关联。火山关系提示 *OOH、*O 和 *OH 不能只看单一物种。DOI: 10.1038/s41467-019-09503-4。




最少需要哪些算例?
一套可执行任务可以从清洁 slab 开始:完整优化目标晶面和终止面,固定底层并保持 15 Å 量级真空层作为初始测试;随后对 A 位、B 位和桥位分别放置 O₂、*OOH、*O、*OH。每个吸附态优化完成后,再做静态计算和电子结构分析。截断能、k 点、自旋初值、U 值和色散修正必须在同一批比较中保持一致或明确测试边界。
若目标是“谁更适合吸附 O₂ 与 *OOH”,最终表格至少应包含位点名称、最终吸附构型、关键键长、ΔG*OOH、ΔG*OH、局域磁矩和电子结构旁证。尖晶石铁氧体 ORR 的 VASP 结果应以 位点矩阵 为主线,而不是以单个最低能构型收尾。
哪些结果不能混在一起?
不同晶面、不同终止面、不同 U 值和不同磁态初值下得到的吸附能不能直接混合排序。若 B 位模型来自 (001) 面,而 A 位模型来自 (111) 面,结论只能说明两个表面模型的差异,不能说明 A/B 位的本征差异。若吸附态优化后表面发生明显重构,还要判断该结构是否仍代表原始位点。
下一步可以把代表位点放入显式水层或隐式溶剂模型中,检查 *OOH 与 *OH 是否受界面氢键改变。对于铁氧体这类磁性氧化物,结构、磁态、自由能和电子结构四类结果需要 在同一个模型族内闭合,这样 A 位和 B 位的比较才有明确边界。




要点:尖晶石铁氧体的 A 位和 B 位会改变表面配位、磁态和含氧中间体键合方式。
要点:O₂、*OOH、*O、*OH 应按同一晶面和终止面构建位点矩阵,不能只比较一个 O₂ 吸附能。
要点:ORR 自由能台阶需要和最终构型、局域磁矩、PDOS 与差分电荷密度一起判断。
要点:不同晶面、U 值或磁态口径下的结果不能直接用于 A/B 位本征排序。
