极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析

说明:本文华算科技主要介绍极限电位与理论过电位的热力学定义、氧化和还原方向的取值方法,以及电位标尺、反应路径和自由能修正对数值的影响。

电催化计算常将多电子总反应写为若干吸附中间体之间的基元步骤,再逐步计算吉布斯自由能。极限电位 UL 描述给定路径在热力学上刚好可以连续进行的电位阈值,理论过电位 ηth 则记录该阈值偏离平衡电位 Ueq 的幅度。两个量共享同一组自由能数据,却对应不同的参照关系。

实验极化曲线中的“极限电流”描述高过电位区受扩散控制的电流平台,符号常写作 jlim;本文的极限电位 UL 是自由能路径给出的热力学量。二者名称相近,物理量、单位和求取方法均不同。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析

极限电位的定义来自哪两个电位?

一个含 n 个电子的总反应具有由反应物和产物化学势决定的平衡电位 Ueq。在相同温度、压力、活度和电位标尺下,Ueq 由总反应自由能确定,路径中插入多少个表面中间体不会改写这个数值。水的四电子氧化与氧还原在标准条件下共享 1.23 V(vs. RHE)的平衡电位,反应方向决定外加电位应向哪一侧移动。

标准状态下的水氧化自由能为 4.92 eV,每转移 4 个电子,对应 4.92/(4e)=1.23 V。中间体吸附会重新分配四个步骤的 ΔGi,但各步之和仍为 4.92 eV。总自由能给出 Ueq,逐步自由能给出 UL,两种取值由此分开。

表面反应路径给出另一组量:*、*OOH、*O、*OH 等状态之间各有一个 ΔGi。同一个总反应可能经过缔合、解离或不同电子数的分支,选定路径及中间体集合决定后面参与取极值的 ΔGi。若遗漏一条能量较高的电化学步骤,算出的 UL 会失去原路径含义。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图1. 氧还原反应的二电子、四电子缔合路径与四电子解离路径。DOI:10.3390/molecules28217264

计算氢电极(CHE)用氢气参照质子—电子对:在 0 V(vs. RHE)时,μ(H++e−)=½G(H2);电位升高 U 后,该对粒子的电化学势降低 eU。这一步把难以在常规周期 DFT 中单独处理的溶液质子和电极电子,换成气相 H2 与电位项。

对于消耗 ni 个质子—电子对的还原步,ΔGi(U)=ΔGi(0)+nieU;对于释放这些粒子的氧化步,符号改为负号。电位改变会按电子数移动每个电化学台阶,未发生电子转移的化学步不会得到同样的线性位移。

在给定路径中逐步移动电位,直至所有电化学步骤都满足 ΔGi(U)≤0,其中至少一个步骤恰好等于 0。这个电位就是极限电位 UL。图中的 0 V、UL 和 Ueq 三组曲线记录同一批中间体在不同电子化学势下的自由能位置。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图2. Mn、Fe、Co 和 Ni–N2O2Cx 位点在 0 V、极限电位和 1.23 V 下的四电子氧还原自由能。DOI:10.3390/molecules28217264
极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析

氧化与还原怎样分别取极限值?

对于向电极释放电子的氧化步骤,升高电位会按该步电子数降低反应自由能。第 i 步达到 ΔGi(U)=0 所需的电位为 ΔGi(0)/(nie),整条路径要满足全部步骤,必须越过数值最高的那个阈值。氧化方向写成UL,ox=max[ΔGi(0)/(nie)],每一步的电子数都要保留。

NiFe-BTC//G 的 OER 路径包含 *OH、*O 和 *OOH。其 Ni 位点上较大的单步自由能为 1.40 eV,对应 UL,ox=1.40 V。同一复合材料中的 Fe 位点与体相 MOF 给出另一组自由能台阶,控制步骤随吸附位点改变,UL 因而属于“材料—位点—路径”这组三元对象。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图3. NiFe-BTC 与石墨烯限域 NiFe-BTC 的电子态、吸附自由能关系、活性图和 OER 自由能路径。DOI:10.1038/s41467-022-33847-z

对于从电极消耗电子的还原步骤,电位降低会按该步电子数使 ΔGi 下移。每一步要求 U≤−ΔGi(0)/(nie),整条路径允许的最高工作电位由这些上限中的最小值控制:UL,red=min[−ΔGi(0)/(nie)]。全部步骤均为单电子转移时,可简写成 −max[ΔGi(0)]/e。

Co–N2O2Cx 的四电子 ORR 例子给出 UL=0.79 V;高于该值时,*OH 的质子化与水生成步骤出现正的 ΔG。触及 ΔG=0 的步骤称为电位决定步骤,这个名称指向热力学台阶。反应速率还受过渡态自由能、覆盖度和静息态分布支配。

分子 HER 路径还会混合电子转移、质子转移和 PCET。若某一步传递两个电子,直接把其 ΔG 除以 e 会把所需电位放大一倍;若某一步不含电子,外加电位在简单 CHE 线性式中不移动该台阶。逐步标注 ni之后,取极值才对应所写反应网络。

若路径含 ΔG>0 且 ni=0 的纯化学步骤,线性 CHE 中任何电位都不会消除该自由能代价。此时 UL 仍可由电化学步骤计算,整条路径却保留一个电位无法调平的热力学上坡步。这类反应网络需另查化学步骤的驱动力与活化自由能。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图4. 四种分子催化剂在 U=−1.26 V、pH=7 时含或不含 PCET 的 HER 自由能路径。DOI:10.1021/acsomega.5c10226
极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析

理论过电位怎样从电位差得到?

平衡电位描述总反应达到热力学平衡的位置,极限电位描述选定表面路径中全部电化学台阶刚好顺行的位置。两者之间的电位差定义理论过电位。氧化方向需要把电位升到 UL,写作ηth,ox=UL,ox−Ueq;还原方向需要把电位降到 UL,写作ηth,red=Ueq−UL,red。

用上面的 OER 数值,1.40 V 与 1.23 V 的差得到 ηth=0.17 V;Co–N2O2Cx 的 ORR 数值则为 1.23−0.79=0.44 V。同一个 1.23 V 基准在两种反应方向中采用相反的减法顺序,否则会得到负的“过电位”。

理想四电子 OER 把总计 4.92 eV 的自由能平均分配到四个单电子步骤,每步均为 1.23 eV,此时 UL=Ueq,ηth=0。真实表面对 *OH、*O 和 *OOH 的稳定幅度不同,最大的逐步阈值偏离 1.23 V,理论过电位量化的是这种台阶不均衡,数值越小表示所选热力学路径越接近均匀分配。

碱性 ORR 若使用 SHE 标尺,平衡电位会随 pH 移动。图中的 pH=14、Ueq=0.40 V(vs. SHE)与酸性条件下常见的 1.23 V(vs. RHE)属于两套参照写法。电位标尺和 pH 必须成对报告,把 0.40 V 与 1.23 V 直接相减没有同一热力学基准。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图5. Fe@NC、Cu@NC 和 Fe–Cu@NC 在 pH=14 时的 0 V、平衡电位与极限电位 ORR 自由能。DOI:10.1038/s41598-024-61223-y

理论过电位读取的是基元步骤自由能,不含跨越过渡态所需的活化自由能。电流密度与 Tafel 斜率还受电荷转移系数、覆盖度、双电层、传质和有效位点数影响;实验过电位还混入欧姆降、气泡遮挡与电极结构。ηth 适合比较同一计算口径下的热力学电位需求,实验速率需由动力学计算和电化学数据另行约束。

部分文献把 UL,ox 写作理论起始电位 Uonset,th,或把 max[ΔGi]/e 直接称为 OER 的限制电位。术语核对要回到公式和反应方向:数值是与 Ueq 相减前的电位,还是已经相减后的 ηth,两列数据的单位虽同为 V,含义不同。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析

哪些计算条件会移动极限电位?

RHE 标尺把理想质子—电子对的 pH 位移吸收到电位定义中,SHE 标尺则保留约 −0.059 V×pH 的能量修正(298 K)。路径若含单独质子转移、单独电子转移或溶液物种酸碱平衡,斜率会改变,控制步骤还可能在某个 pH 切换。电位—pH 图上的斜线、水平段和突变分别对应 PCET、纯电子转移和优势路径变化。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图6. 四种分子催化剂使 HER 路径达到零能量跨度的电位–pH 关系。DOI:10.1021/acsomega.5c10226

溶剂会按中间体的偶极矩、氢键受体能力和带电程度产生不同稳定作用。Fe@NC、Cu@NC 与 Fe–Cu@NC 的气相和溶剂自由能曲线差别集中在 *OOH、*O 和 *OH,并非整条曲线刚性平移。某一步的 ΔGi 改变后,取极值的步骤及 UL都可能换成另一台阶。

极限电位与理论过电位:热力学定义、取值方法及计算条件影响解析
图7. Fe@NC、Cu@NC 和 Fe–Cu@NC 的气相与溶剂模型 ORR 自由能路径。DOI:10.1038/s41598-024-61223-y

零点能、振动熵、气体压力和溶液活度进入各物种的 G;交换相关泛函、自旋态、色散修正和构型搜索改变 DFT 能量;覆盖度、共吸附、界面电场与表面电荷改变工作态吸附强度。这些修正应逐项进入同一套 ΔGi,混用不同模型得到的单步能量会破坏取极值所需的共同基准。

CHE 把每个 PCET 步骤近似为随 U 线性变化,真实界面中的吸附偶极、差分电容和表面电荷会引入非线性响应。恒电位计算可检验 ΔGi(U) 的斜率,特别适合吸附态偶极差大、界面场强或电荷数随电位显著变化的表面。

极限电位由某一个逐步阈值控制,该步的自由能误差会直接传到 UL。若两个候选控制步骤只差 0.05 eV,而溶剂修正或泛函差异约为 0.2 eV,控制步骤排序应写成存在竞争,不宜给出唯一机理归属。相同口径下的相对趋势通常比跨论文比较两位小数更有意义。

一组可复核的结果应列出反应方向、电位标尺、pH、温度、所选路径、每步电子数和自由能修正,并给出全部 ΔGi、控制步骤、UL、Ueq 与 ηth。实验对照需再注明电流密度、iR 校正、传质条件和催化剂负载量,热力学电位差与测得电极响应才具有可比口径。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
赞 (0)
上一篇 3天前
下一篇 2小时前

相关推荐