XL-BOMD 全原子模拟显示氮掺杂石墨烯上氧还原随偏压在内球与外球机制间切换

该工作首次把 XL-BOMD 用于多相电催化，在统一框架内同时显式模拟电极与电解质并施加电位偏压，在 622 原子、1840 电子的 NG/水/O2 体系中稳定运行到 20 ps，时间步长 0.1 fs，每个 MD 步约 0.15 秒。 此前静态 DFT 与计算氢电极方法假设质子耦合电子转移、无法描述非吸附型电荷转移，而常规量子分子动力学受 Kohn-Sham DFT 立方标度和自洽场收敛问题限制；XL-BOMD 把电子自由度作为动力学变量，用阴影势避免每步迭代自洽场。 论文报告了 NVE 下 10 ps 无能量漂移、GPMDK 随体系规模的计时标度、以及预条件使初始自洽场收敛呈二阶幂律；DFTB 与 DFT 的吸附能、水分子构型、O2 键长（1.24 Å）与实验值做了对照。

施加到 NG 片的偏压决定 ORR 走内球还是外球路径：较高偏压（0 至 3.0 eV）下 O2 无需化学吸附，电子经溶剂介导转移，4 个电子分两步从 NG 片转移并最终生成 4 个 OH−；较低偏压（-1.25 至 -0.25 eV）下 O2 先吸附到 NG 表面再发生直接电子转移。 论文提出一种通过移动 NG 片原子在轨能量来施加电化学偏压的方法，与恒电位方法不同，体系保持 NVT 正则系综，从而可以系统扫描反应驱动力。 结论基于 58 个最长 20 ps 的产生模拟，内球偏压 10 个、外球偏压 20 个，每个偏压做 3 至 5 次重复；用 Mulliken 电荷与 O-O 间距判定物种与断键，并设了三个对照体系（NG 真空、纯水、纯石墨烯加水），对照中均未观察到 ORR。

模拟得到的反应动力学无法被 Butler-Volmer 和 Marcus-Hush-Chidsey 理论完整描述：内球与外球数据都能拟出两条 Tafel 斜率，但模拟中并未出现对应的速率决定步骤或机制变化，Tafel 斜率呈连续弯曲；用 MHC 理论可以拟合全部数据点并得出看似合理的直线，但模拟明确显示存在两种不同机制。 这提示仅凭宏观动力学拟合推断机制可能产生误导，需要原子尺度的逐步反应图像来补充。 证据来自对 58 个模拟的反应时间统计（反应时间定义为生成全部 OH− 的时间），以及内球 49.2 与 189 mV dec−1、外球 49.3 与 316 mV dec−1 等拟合斜率；论文也指出这些数据范围未覆盖 Tafel 分析通常偏好的一个数量级电流区间。

氮掺杂使石墨烯更活泼的电子结构原因被具体化：与纯石墨烯相比，氮掺杂既移动电子化学势以提高还原电位，又增大费米能级附近的态密度，从而提供电子来源；ORR 后 NG 被氧化产生四个电子空穴，主要定域在与氮直接成键的碳原子上，并可通过共振结构在整个 NG 片上离域。 这把宏观上 NG 优于纯石墨烯的催化活性与可计算的态密度和空穴分布联系起来。 依据是 O2、NG 中氮邻近碳与纯石墨烯的态密度对比图，以及氧化后 0.5 ps 内各近邻碳原子 Mulliken 电荷平均值（第一近邻碳约 0.202，第二近邻约 -0.019）。

这项工作面向电催化理论、计算电化学与量子分子动力学方法研究者，适用于在显式溶剂中研究固-液界面电位依赖反应机制的设定。它使研究者能够在同一模拟中同时处理电极与电解质、施加可控偏压并逐步观察电子转移事件，作者据此提出后续可扩展到稳态电流条件、显式 pH 效应以及缺陷或层状催化剂结构。结果针对的是理想极化电极上的孤立单分子 ORR 事件，作者有意避开覆盖度效应、吸附质相互作用与多位点竞争吸附等复杂因素。

读者仍需留意若干开放问题：偏压到过电位的换算被作者说明需要进一步测试与验证；模拟条件被描述为很可能对应高过电位，且使用了初始 O2 靠近表面与引导分子动力学等加速策略；内球路径的最后一步（吸附氧与水产成两个额外 OH−）在现有模拟中未观察到，作者推测与部分还原或模拟时长有关；瞬态物种的判定依赖 Mulliken 电荷与参考表比对，作者指出 Mulliken 电荷并不正式对应氧化态，且吸附的 O2− 与 O·− 在孤立模拟中电荷相近而难以区分；外球机制中瞬态物种的形成可能与无稳态电流、无电子库、无恒定偏压的模拟条件有关。此外，本文为快速解析文本，图表与支持信息中的具体数值细节需回到原文核对。