“内建电场”：基本概念及催化应用！

说明：本文华算科技主要介绍内建电场的电荷与电势来源、吸附态和过渡态对电场的响应，以及内建电场在二氧化碳还原、氨合成和单原子电催化中的作用。

内建电场的基本概念对应哪些计算量

内建电场源于体系内部正负电荷中心的空间分离。半导体 p-n 结中的空间电荷区、两种材料接触后的费米能级平衡、铁电晶格的自发极化、Janus 结构两侧的化学环境差异，都会使静电势沿特定方向连续变化。热平衡结构在撤去外部电源后仍保留的场属于内建分量；电催化工作电位、双电层离子和溶剂取向产生的场属于运行分量，实际活性位感受到的是各分量的矢量和。

DFT 自洽结果给出电子密度 n(r) 和静电势。界面法向取为 z 时，三维量经平面平均得到 ρ(z)、φ(z) 与场强 E_z(z)。软件若输出电势 φ，关系式写作 E_z=−dφ/dz；若输出电子势能 U_e=−eφ，则 E_z=(1/e)dU_e/dz。两种曲线的斜率符号相反，报告场方向时应写明坐标正向和输出量定义。

图1对应电催化条件下的总界面场，其中负偏压电极和阳离子分布随工作电位变化。内建分量与电位诱导分量应分别命名：前者可由无外场弛豫结构的电荷密度和电势获得，后者要在给定电荷或给定电位模型中求取。只报告一条电势曲线却省略电位、离子分布和参考态，会使“内建”一词失去明确物理对象。

界面势降 Δφ 除以有效场区宽度 d 可给出平均场强 E_avg≈−Δφ/d，这一数值适用于势能近似线性的区域。原子核附近的势起伏很大，活性位周围的局域场还具有横向变化，平均场强和原子尺度局域场回答不同的空间问题。前者适合比较界面整体极化，后者直接关联吸附键和过渡态电荷分布。

差分电荷密度 Δρ(r) 标出接触或吸附前后的电子积累与耗尽，平面积分可得到面电荷和偶极密度；Bader 电荷给出分区电子数；静电势提供势降和斜率。电荷重排是场的源项，电势梯度是场的直接读数，两类结果在同一坐标和参考态下互相校验。

电场为何会改变吸附态和过渡态能量

吸附能对场强的响应来自偶极矩与极化率

吸附体系沿表面法向具有偶极矩 μ_z 时，外场或内建场会改变其能量。弱场范围可写成 ΔE(F)≈−μ_zF_z−½α_zzF_z²，其中 α_zz 是法向极化率。一次项控制正负场方向造成的非对称响应，二次项描述电子云随场强变形。分子偶极、吸附键取向和表面电荷共同决定 μ_z 的正负。

Ru(0001) 与台阶 Ru(1013) 上的结果显示，N₂^*、H^* 等低偶极吸附态对电场较迟钝，NH_x^* 和 N₂H^* 的响应更强。台阶位会改变吸附构型和电荷中心，因而同一种中间体在两个晶面上具有不同 μ_z。场敏感度由“中间体—位点”整体偶极决定，气相分子的永久偶极只占其中一部分。

过渡态响应取决于初态与过渡态的偶极差

基元步骤的势垒变化由过渡态和初态的响应差决定，可写为 ΔΔG^‡(F)≈−(μ_TS−μ_IS)·F−½(α_TS−α_IS)F²。场向量与偶极差的相对方向决定势垒升降，场强增加并没有固定的活性增益。初态、过渡态和终态若具有相近偶极，整条能量曲线可整体移动，而势垒变化很小。

CO₂ 的线性构型在气相中没有永久偶极，靠近带电表面后会发生弯曲、极化和部分电子占据，吸附态随之获得显著偶极。图3中金属/石墨烯主体 d 态对场的变化较小，CO₂ 相关 s、p 态和两种 O 端取向出现差别。吸附物电子态的场响应可强于催化剂主体能带变化，催化效应由界面复合态承担。

覆盖度升高会引入偶极—偶极相互作用，溶剂分子和阳离子会改变局域介电屏蔽，吸附构型翻转则会改变 μ_z 符号。场强、覆盖度、吸附取向和界面组成应保持同一计算状态，否则两组吸附能之差会混入多种结构变量。

内建电场怎样调节催化活性与选择性

二氧化碳还原中的极性中间体

CO₂ 还原的早期步骤会形成弯曲的 ^*CO₂、^*COOH 或 ^*OCHO，这些吸附态具有明显电荷分离。金属/石墨烯模型在文献定义的 −1 V Å⁻¹ 场下改变 CO₂ 吸附和后续自由能，CO 路线与 HCOOH 路线出现不同幅度的能量移动。选择性变化来自各中间体偶极与极化率的差异，也受到金属种类和表面电荷的约束。

图4中的五种金属表面并未得到同一种自由能重排。Ir/石墨烯在该模型场方向下更偏向 CO 生成，部分体系的 HCOOH 中间体或 H^* 响应较弱。催化选择性取决于竞争路线的相对能量差；降低某个吸附态能量后，另一竞争步骤仍可能控制速率，目标产物速率未必升高。

氨合成中的反应序列切换

Ru(1013) 台阶面上的计算把 N₂ 直接解离和 N₂H^* 形成两条反应序列放在相同场强下比较。正场超过约 0.63 V Å⁻¹ 后，N₂H^* 先形成的缔合序列在热力学上占优；NH₂^* 氢化的势垒下降，473 K 下正向速率常数增加。场效应能够改写基元步骤先后顺序，整条自由能曲线和速率常数记录了序列切换，单个 N₂ 吸附能缺少过渡态与后续氢化信息。

纳米曲率产生的局域场调节

球形碳载体的曲率会改变表面电荷分布和局域场强，M–N₄ 单原子位点因而处于可调电场环境。Ni–N₄ 的计算中，^*CO₂→^*COOH 自由能从无场时 0.96 eV 降至 −0.4 V Å⁻¹ 时 0.88 eV，场强增至 −0.8 V Å⁻¹ 后回到 0.90 eV。活性随场强呈现最适区间，过强场会让另一中间体或竞争步骤成为新的能量限制。

Ni、Fe、Co 单原子位点对同一曲率并不呈现相同响应，CO₂ 还原、氧还原、碱性析氢和析氧的最适球径也各不相同。酸性析氢中的 H^* 偶极矩和极化率很小，曲率变化对该反应的影响不明显。中间体偶极谱决定每类反应的场强窗口，材料几何只负责建立可用的局域场。

哪些计算证据能确认催化场效应

场的来源、方向和作用区域应有独立读数

一组可复核结果通常包含四类量：平面平均电荷或差分电荷密度给出源项，静电势斜率给出场方向与平均强度，吸附态偶极和极化率给出能量响应系数，自由能或势垒随场强的斜率给出催化后果。四类量应采用相同的法向、符号约定和结构状态，场从哪侧指向哪侧、作用于哪段空间便有明确答案。

计算得到的电场还可与实验局域探针比较。开尔文探针力显微镜反映表面电势分布，吸附 SCN⁻ 的振动 Stark 位移反映活性位附近场强随电位的变化。Ni-SAC-70 的 Stark 调谐率约为 21 cm⁻¹ V⁻¹，Ni-SAC-350 约为 5 cm⁻¹ V⁻¹，曲率增大对应更强的活性位局域场响应。

模型条件限定场强数值的可比范围

带真空层的非对称 slab 若缺少偶极修正，周期镜像会在真空区引入伪电场；带净电荷的模型还涉及补偿背景和电位零点。真空厚度收敛、偶极修正和电荷补偿属于场强读数的前提，它们会改变 φ(z) 的斜率与平台差。比较两种材料时，晶面、slab 厚度、覆盖度和吸附物取向应分别采用相同取值。

无外场静态 DFT 适合提取极性晶体或异质界面的平衡内建分量；施加锯齿势的定场计算适合测量 μ 和 α；恒电位或带显式溶剂的界面模型用于工作电位下的总局域场。三类模型对应的热力学约束不同，场强数值应随模型类型、参考电位、离子组成和溶剂构型一同报告。

电场引起的吸附能变化若要进入电催化自由能，还应使用同一电位下的零点能、熵、溶剂和覆盖度修正。过渡态偶极差、场强空间分布和恒电位自由能共同限定速率变化；活性位附近的溶剂重排与阳离子配位则决定该场在真实双电层中的幅值。