ORR中2e⁻/4e⁻选择性怎么用DFT解释？

说明：本文华算科技主要介绍氧还原两电子与四电子路径的反应分支、自由能和活化能垒，以及电位、溶剂和过氧化氢后续反应对选择性的影响。

两电子与四电子ORR的基本反应是

氧还原反应（ORR）的电子数按每个 O₂ 分子计。在酸性介质中，两电子总反应为 O₂＋2H⁺＋2e⁻→H₂O₂，四电子总反应为 O₂＋4H⁺＋4e⁻→2H₂O。电子数记录总反应计量，不要求电子在同一个基元步骤中一次转移完毕；表面吸附、质子转移和中间体转化可以分步发生。

生成的 H₂O₂ 还可经历 H₂O₂＋2H⁺＋2e⁻→2H₂O。先生成过氧化物、再继续还原的“2＋2”过程，同样总计消耗四个电子。图1分别画出过氧化物离开碳表面和继续反应的情形，并以钠离子调节局域质子分布的模型作对照；最终检测到多少过氧化物，受生成与消耗共同影响。

在常用的缔合机理中，O₂ 通过首个质子—电子转移形成 ^*OOH，其中星号表示表面位点。保留 O—O 键并释放过氧化物可形成两电子产物；若 ^*OOH 发生 O—O 断裂，形成 ^*O 等吸附态，随后经 ^*OH 继续生成水。对先解离氧分子再加氢的表面，初始反应网络另有结构，不能预设所有位点都共用同一个 ^*OOH 分支。

在碱性介质中，两电子总反应常写为 O₂＋H₂O＋2e⁻→HO₂⁻＋OH⁻。水承担质子供体作用，溶液产物可为 HO₂⁻；其与 H₂O₂ 的比例随酸碱平衡变化。DFT 比较酸、碱两套模型时，质子供体、溶液物种和自由能参照必须相互一致。

吸附自由能和火山图解释哪部分选择性

DFT 总能经过零点能、熵及溶剂修正后，才能构造吸附态与溶液物种的自由能。计算氢电极模型在 RHE 标度下用 μ(H⁺＋e⁻)＝½G(H₂)−eU 处理质子—电子对，消耗 n 个质子—电子对的还原步含 ＋neU 项；这里 e 为正的元电荷。若使用 SHE 标度，还要处理质子活度；RHE 已吸收相应电位换算，重复加入同一 pH 校正会改变能量基准。

两电子路径既要形成 ^*OOH，又要使其转化为溶液过氧化物。过弱的 ^*OOH 吸附增加形成中间体的自由能代价，过强的吸附则可能增加释放产物的代价。含氧碳位点的火山图用以水和氢为参照的 ^*OOH 自由能作横轴，用两步均可热力学下坡的最大电位 U_L 作纵轴；峰顶表示这两步之间的平衡。

该火山模型只包含所选两电子反应步骤，O—O 断裂和 H₂O₂ 再还原没有自动进入纵坐标。位于峰顶的位点可能有较小的两电子热力学过电位，却仍有快速的竞争支路。酸性标准条件下，氧还原生成水的平衡电位约为 1.23 V，生成过氧化氢约为 0.70 V；水端更低的总自由能也没有规定反应实际经过哪条路径。

B、P、N、S 掺杂碳的计算在相同电位下比较 ^*OOH、^*O、^*OH 和产物状态，四电子分支的产物端自由能更低；这些材料却表现出过氧化氢选择性。自由能节点只给出中间态相对高低，节点之间尚未包含断键和溶剂重排所需的活化自由能。图3的曲线按每个氧分子的反应计量比较，图中的连线没有给出真实过渡态。

比较两电子和四电子路径时，两组初态须对应相同电位、氧与氢来源及活性位状态。若一组使用中性真空表面，另一组使用带电显式水界面，差值还混入电荷与溶剂模型的差别；把各条曲线的终点任意归零，也会丢失竞争反应原有的驱动力差异。

DFT怎样用活化能垒解释反应分支

从同一个 ^*OOH 前体出发，过氧化物释放与继续生成 ^*O 具有不同过渡构型。两条支路的活化自由能须从共同前体量起，并保持电位、温度和质子条件一致。显式水模型可保留水分子重排与质子转移；真空静态 NEB 只描述相应势能路径，直接把它的能垒当作水界面的活化自由能，会遗漏环境贡献。

在 B 掺杂碳的恒电位显式水模拟中，^*OOH 脱离为水相过氧化物前体的自由能垒约为 0.53 eV，形成 ^*O 的分支约为 0.70 eV。较低的产物释放能垒与较高的 O—O 断裂能垒使两电子反应获得动力学优势，尽管四电子产物端更低。图4保留了水界面与不同阶段的原子构型；这两个能垒仅适用于论文采用的模型条件。

对同一前体的两个竞争基元步骤，过渡态理论给出 k₂/k₄≈(A₂/A₄)exp[−(ΔG₂^‡−ΔG₄^‡)/(k_BT)]。能垒差通过指数影响分支速率，不同反应前因子仍可改变两者的比值；若 ^*OOH 长时间占据大部分位点，新的氧分子吸附会受抑制。把速率常数比直接写成实验过氧化氢百分数，会省去位点占据、后续反应和产物损失。

双吡咯氮位点的碳模型中，^*OOH 在水界面的朝向和周围氢键影响其脱离与断裂。沿 C—O 与 O—O 键伸长分别计算自由能，可比较表面—吸附物分离和分子内部断键。图5同时给出平均约束力与积分得到的自由能；波动的约束力曲线和红色自由能曲线具有不同物理量与纵轴。

该碱性显式水模型得到的两电子与四电子分支能垒分别约为 0.54 和 0.92 eV，局域构型对两种过渡态的稳定程度不同，才产生速率差异。若只给出 Bader 电荷、态密度峰或吸附后 O—O 键长，尚缺这两类反应坐标上的活化自由能；单个电子结构变化也没有给出产物选择率。

电位/界面环境/产物去向怎样改变选择率

电位改变界面电荷、吸附物偶极与水分子排列，过渡态对电位的响应可与吸附初态不同。在反应路径上维持恒电位，需要随构型调整电子数；仅给反应终态加上计算氢电极校正，尚未计算过渡态的带电与溶剂响应。将电位依赖的速率常数写入微观动力学方程后，可求解空位、^*OOH 和 ^*OH 各自所占的位点比例与稳态产物流量。

酸性碳电极上的钠离子研究把界面质子分布、^*OOH 分支与 H₂O₂ 进一步还原分别计算或测试。局域质子供给同时影响初始选择与产物存留：即使已经生成 H₂O₂，后续还原仍会减少收集到的产物。图6中的水界面模拟、活化能和过氧化氢还原电流分别对应这些环节；体相溶液 pH 相近时，局域质子环境仍可能不同。

旋转环盘电极中，环电流记录从盘面离开并到达环的过氧化物，盘电流还包含生成水的电荷消耗。图7将盘环电流、过氧化氢选择率与平均电子数 n 分别列出。按只有两种终产物、收集效率已校准且没有额外损失的模型，若每个氧分子进入两电子终产物的比例为 S₂，则 n＝4−2S₂。

在上述理想两产物模型中，过氧化氢法拉第效率为 2S₂/(4−2S₂)，分母按两条路径各自消耗的电荷计数，与分子比例 S₂ 使用不同权重。实验存在析氢、过氧化物分解或收集损失时，还需计入相应通量；平均电子数接近 4，也可能含先生成过氧化物再还原的贡献。

将 DFT 结果映射到产物流量时，过氧化物再吸附、继续还原和离开催化层的速率与初始分支共同参与稳态质量守恒。位点间距离、催化层厚度和传质改变停留时间，同一个局域位点模型可对应不同的收集选择率；产物脱离后的溶液浓度及其向相邻位点扩散，也进入过氧化物继续消耗的速率。