电催化阳离子效应全解析：双电层分布、作用机制与第一性原理计算物理量

说明：本文华算科技主要介绍电催化界面中的阳离子效应、阳离子在双电层内的分布与作用方式，以及第一性原理计算区分静电场、短程配位和界面水重排时采用的物理量。

什么是阳离子效应？

更换 Li⁺、Na⁺、K⁺ 或 Cs⁺ 后，电极材料没有变化，产物电流、选择性和起始电位却可能明显移动。阳离子既不进入总反应式，也不必留在产物里；它改变的是带电表面附近的溶剂、局域电势与反应物状态。计算若只保留一块真空中的带电表面，便缺少产生这类差异的关键环境。

同一个阳离子在不同电极、阴离子、溶剂和反应路径上可以给出不同排序。Cs⁺ > K⁺ > Na⁺ > Li⁺ 只适用于部分体系，不能当作普适活性序列。阳离子效应的计算任务，是把离子身份转化为可核对的界面数密度、配位构型、电场、表面电荷和自由能变化。

阳离子效应的计算对象是什么？

在电极组成、晶面、阴离子、体相浓度、温度和电位可比时，仅改变电解质阳离子的种类或界面丰度，反应速率、产物分布或中间体自由能随之发生变化，这组关联通常称为阳离子效应。定义里真正受比较的是电极—水—离子共同组成的界面，不是孤立的 Li⁺ 或 K⁺。

负偏压使电极带负电，反离子向界面富集。离子中心仍被水合壳包围，其最近可达位置常用外 Helmholtz 面描述。裸离子半径与水合后的有效尺寸方向相反：Li⁺ 裸离子较小、束缚水较强，Cs⁺ 水合较弱，靠近带电界面时付出的脱水代价通常较低。

体相浓度记录总体配比，反应位点旁还存在随时间涨落的局域浓度。表面电荷、离子水合、阴离子竞争吸附和水分子取向共同决定界面富集。恒浓度实验中换一个电位，阳离子数密度峰值也会移动；恒电位比较中换一个阳离子，表面为维持该电位所需的电子数又可能改变。

固定在电极附近的季铵基团也能形成局域正电荷、改变 Stern 层电势，并抑制 H⁺ 穿过界面，功能上可替代部分可移动金属阳离子。“阳离子参与”不要求发生法拉第反应，它可以在多个电子转移循环中保持自身组成。

朝向溶液侧暴露 O 原子的 ^*CO₂、^*COOH、^*CHO 和 ^*OCCO 容易形成阳离子—氧近邻；吸附较强、负电荷较少或 O 原子埋在表面的物种，响应可能很弱。能否形成稳定近邻取决于瞬时几何与电荷分布，距离轨迹比单个优化构型更能反映溶液中的接触频率。

阳离子进入界面后会改变哪些计算量？

显式溶剂第一性原理分子动力学可记录阳离子沿表面法向的数密度 ρ(z)、阳离子—O 径向分布函数 g(r) 和配位数。ρ(z) 给出离子在哪个距离区间富集，g(r) 的第一峰和积分配位数则区分水合壳与中间体配位。二者来自时间平均，采样长度会影响峰高与误差。

恒电位计算还要记录表面电荷密度 σ。外加电位相同并不保证不同电解质下 σ 相同，因为双电层电容和零电荷电势会随离子分布变化。实验中 CO 生成电流与 C₂H₄ 生成电流可在多种浓度和阳离子下近似随 |σ| 线性变化，表面电荷是连接体相电解质与局域反应速率的中间变量，浓度本身不承担这一角色。

局域电场可由平面平均电势或电荷密度求得，也可通过吸附物振动 Stark 位移间接校准。对于带偶极矩 μ 和极化率 α 的吸附态，场贡献常写成 ΔG_field≈−μ·E−αE²/2。不同中间体对同一电场的响应幅度并不相同，偶极方向相反时，稳定化方向还会翻转。

阳离子靠近 ^*CO₂ 或 ^*OCCO 后，可改变 C—O 键角、吸附高度、Bader 电荷和电荷密度差；含氧物种的未占据轨道下移时，电极电子注入所需的能量降低。自由能降低、稳定接触距离和电荷转移应在同一构型集合中对应，孤立的一张差分电荷图只能描述选定快照。

连续介质模型能描述平均介电屏蔽，却没有离散水合壳和特定离子—中间体接触。静态显式水模型补上局域配位，但其结果依赖初始排布。涉及离子交换、脱水和水网络涨落时需要轨迹采样，自由能宜由约束分子动力学、伞形采样或其他增强采样方法沿反应坐标获得。

阳离子怎样参与电子转移和中间体稳定？

CO₂ 接受电子时由近线形分子弯曲为 CO₂^δ−，负电荷主要分布在 O 端。若界面没有足够近的正电中心，溶液中的 CO₂ 可先接受电子，再靠近 Au 表面，属于外球电子转移。阳离子与 O 端形成短程接触后，吸附和电子注入可以协同发生，路径转为内球过程。

外球路径可用约束密度泛函理论结合 Marcus 自由能描述，反应物态和产物态能隙分布给出重组能与电子转移势垒。内球路径则需让 CO₂—Au 距离成为反应坐标，在恒电位分子动力学中计算平均约束力。两条路径对应不同的电子态与溶剂坐标，不能把某个吸附能直接当作完整活化势垒。

阳离子稳定负电中间体时，电子转移与离子靠近彼此耦合，可写作阳离子耦合电子转移。它与质子耦合电子转移的差别在于，M⁺ 调节局域能级与电势，但不提供形成 X—H 键所需的质子。M⁺ 的配位数可随电子注入同步增加，反应后阳离子回到界面液相。

在 Au(111)—水体系中，纯水条件下 CO₂ 靠近表面仍有较高自由能代价；加入 K⁺ 或 Li⁺ 后，内球吸附路径出现可达的自由能曲线。K⁺ 与 Li⁺ 的差别包含配位几何和水合壳重排，曲线差值不宜全部命名为静电场贡献。

真实界面既有宏观电场，也有短程配位。计算分辨二者时，可设置相同表面电荷并改变阳离子位置，或比较固定正电基团与可移动金属离子；还可用受限坐标阻止直接接触，再观察自由能差。对照模型必须保持电位、表面电荷或离子数中的比较基准明确，否则两项效应会一并变化。

为什么 Li+、Na+、K+、Cs+ 的排序会变化？

较弱水合的 K⁺、Cs⁺ 常能靠近外 Helmholtz 面，增强局域电场并稳定带负电的含氧中间体；较强水合的 Li⁺ 更易保留完整第一水合壳。这一趋势成立的前提是反应受带电中间体稳定控制，且界面没有强吸附阴离子占位。

若速率取决于质子到达表面的速度，阳离子改变水网络后会产生另一种响应。酸性 H₂O₂ 电合成中，从 Li⁺ 到 Cs⁺，界面水的 O—H 伸缩带和氢键组分发生变化，质子扩散电流下降，而不带质子传递的 Ru(NH₃)₆³⁺ 探针电流接近不变。阳离子此时主要调节氢键连通与质子迁移，并非只稳定某个吸附态。

路径积分分子动力学把质子的核量子效应纳入水网络采样。Li⁺ 周围紧密且较连续的水合网络，为质子跳跃保留较多供体—受体连接；Cs⁺ 附近水分子排列更松散，长程氢键链减少。水网络差异可用氢键数、供受体距离分布和质子转移自由能描述，平均 O—O 距离无法记录每条质子迁移路径出现的概率。

Au 上 CO₂ 活化重视阳离子—O 配位，Cu 上 C—C 偶联还受 ^*CO 覆盖度影响，Pt 族表面的析氢则对水解离和 H⁺ 输运敏感。同一阳离子序列跨材料失效时，变化的常是控速步骤与界面覆盖度，阳离子的水合化学并未凭空反转。

有限尺寸模型中的一个阳离子可能对应数 mol·L⁻¹ 的界面浓度，周期镜像还会强化离子有序排列。模拟时间仅有数皮秒时，水合壳交换和离子迁移未必充分采样。计算结论应注明电位控制方式、阳离子数、表面面积和采样时间，并把 ρ(z)、配位数、σ、自由能和光谱信号分别与可测量量核对。实验趋势与微观量均获复现后，才可把主导作用归入电场、配位或水网络。