电化学界面计算里，K⁺、Na⁺ 和一层水不是简单的背景环境。它们会改变局域电场、氢键方向、吸附物偶极和质子-电子转移端点。显式阳离子模型要回答的是中间体周围真实出现了哪些近邻相互作用，尤其是 *COOH 和 *OOH 这类带有明显偶极与 O-H 键的中间体。
显式阳离子放进 slab 模型时，不能把它当成新的催化位点随意靠近表面。K⁺ 或 Na⁺ 更多代表电解质双电层里的局域正电中心，通常需要与水分子、吸附物氧端和表面电荷分布共同定义。若阳离子直接与金属表面成键，模型含义会从界面离子效应变成共吸附原子效应。
对 CO₂RR，可把 *COOH 的 O 端附近作为候选位点，同时保留 K⁺/Na⁺ 与水分子的配位距离。图1 的标度关系偏移显示，共吸附阳离子能改变 C₁ 中间体的吸附能关系。阳离子位置必须和吸附物构型一起收敛，否则无法判断能量变化来自局域静电稳定化，还是来自不合理的接触结构。
图1 C₁ 中间体在不同共吸附阳离子下的吸附标度关系偏移。该图用于说明 K、Na、Mg 与吸附物相互作用可改变 *COOH 等中间体稳定性。 DOI: 10.1021/acscatal.4c00727。
界面水层可从单个水、二到四个水分子的小团簇，或短时 AIMD 代表帧开始。对 K⁺，水合壳层较松散，常需要较大的 K-O 距离；对 Na⁺，局域水配位更紧，靠近 *COOH 或 *OOH 时可能改变 O-H 方向。所有初构型都要在同一 slab、同一覆盖度和同一电荷设置下比较。
若只放阳离子而没有水，局域电场会被过度理想化；若只放水而没有阳离子，又会丢失电解质中特定离子的稳定化差异。水层负责给出氢键与介电屏蔽，阳离子负责提供局域正电和配位约束。二者分开建一组对照，才能看出各自贡献。
图2 中 Na、Mg 与 K 的稳定化关系并不完全相同，说明不同阳离子不能只按电荷符号等同处理。VASP 优化后应记录 M-Owater、M-Oads、O-H 方向和吸附物-表面键长，避免把几何差异压缩成单个吸附能数值。
图2 Na 与 Mg 提供的中间体稳定化能和 K 稳定化能的对应关系。该图用于比较不同显式阳离子对 CO₂RR 中间体的非等价影响。 DOI: 10.1021/acscatal.4c00727。
*COOH 的阳离子效应可通过稳定化能拆出来。把同一表面、同一水层、同一阳离子坐标作为参照，分别计算含 *COOH 和不含 *COOH 的端点，再与无阳离子端点比较。这个差值比直接比较 Eads 更容易隔离 K⁺/Na⁺ 对中间体的局域作用。
实际建模时，空白界面端点不能随意重优化成另一种水层排布。若去掉 *COOH 后水分子重新占据吸附位，差分能量会混入水重排贡献；更稳妥的做法是同时保留优化结构和约束对照，用几何变化解释能量变化。
ΔEstab(M⁺) = [E*COOH+M⁺+H₂O/slab – EM⁺+H₂O/slab] – [E*COOH+H₂O/slab – EH₂O/slab] 式(1)
式(1)中 M⁺ 代表 K⁺ 或 Na⁺，H₂O 可代表固定数量的界面水分子。公式前后两组差值必须来自相同 slab 尺寸、k 点、ENCUT、偶极修正和优化标准。若加入阳离子后总电荷设置改变，还要同时说明 NELECT 或背景电荷处理方式。
当 *COOH 的 O 端靠近 K⁺ 时，C-O 键、O-H 键和金属-C/O 接触距离都可能变化。ΔEstab 只在几何与参照一致时才具有清楚含义。若优化后 K⁺ 脱离水层或直接落到表面空位，应该把该构型单独归为强接触态，而不是和外层阳离子态混合比较。
图3 给出 CO₂ 到 CO 路径中 *COOH 描述符随真空、水和阳离子条件变化的趋势。它说明同一个金属表面，在加入界面物种后，*COOH 相关自由能可能移动到火山图的不同区域。VASP 后处理时，ΔG*COOH 应统一包含 ZPE、熵、电位项和界面模型差异。
如果目标是比较 K⁺ 与 Na⁺，应优先固定表面和水层数量，再改变阳离子种类；如果目标是比较水层数量，应固定阳离子并改变 nH₂O。两类对照不能混在同一个横坐标下；对 *COOH 的趋势判断，应同时保留中间体自由能和阳离子近邻结构，再讨论活性变化。
图4 的 Au-water-2K 界面展示了 K⁺ 与水层对 CO₂ 活化路径的影响。对于短时 AIMD 或 cDFT-MD 结果，静态 VASP 可以先取若干代表帧做单点能与局域结构分析，再决定是否继续做更长时间采样。
图3 Au、Ag、Cu 和 Pd 表面经 *COOH 通往 CO 的火山图，图中比较真空、水和不同阳离子条件。该图对应 *COOH 受界面环境调控的后处理方式。 DOI: 10.1021/acscatal.4c00727。
图4 Au-water 与 Au-water-2K 界面中外球电子转移的能量分布和关键结构。该图用于说明 K⁺ 与水层位置会改变界面电场和 CO₂ 活化。 DOI: 10.1038/s41467-023-43300-4。
*OOH 同时包含 O-O 键和 O-H 键，靠近阳离子或水层时会受到局域电场、氢键和表面键合共同影响。对 ORR 或 OER，*OOH 的自由能常和限制电位、H₂O₂ 选择性或四电子路径分叉相连，因此界面坐标改变会直接进入反应趋势。
图5 中 *OOH 在不同碳缺陷与官能团附近的构型响应展示了局域电场效应。若外加电场或阳离子让 *OOH 的末端 O 更靠近水层，O-O 键长和 O-H 取向可能同步改变。此时应把 *OOH 与表面的键长、O-O 距离、Bader 电荷和偶极矩变化一起核对。
*OOH 不能只按单一吸附能判断界面效应。若 *OOH 能量降低但 O-O 键显著拉长，可能意味着路径选择性也在变化；若能量变化很小但 O-H 方向改变，后续质子转移端点仍可能需要重新定义。
图5 不同碳缺陷与官能团上 *OOH 的优化构型及局域电场影响。该图用于说明 *OOH 对界面电场和局域极化的响应。 DOI: 10.1021/acsenergylett.3c02743。
K⁺ 和 Na⁺ 的差异不只是半径。Na⁺ 与水氧的配位更紧，K⁺ 的配位壳层更宽松；在有限水层模型中，这会改变阳离子相对 *COOH 或 *OOH 的可达位置。若两种离子优化后处在不同高度，能量差应解释为离子-水-吸附物三者共同作用，而不是单独的离子大小效应。
图6 的 ORR 界面水和 OOHad/OHad 自由能路径说明，含氧中间体与水层、电解质环境之间存在耦合。同一组阳离子比较至少要保留无阳离子、K⁺、Na⁺ 三个端点，并记录水层构型是否保持同一配位模式。
图6 碱性 ORR 中界面水、OOHad/OHad 和不同电解质环境下的自由能路径。该图对应阳离子-水层改变含氧中间体稳定性的判断。 DOI: 10.1002/anie.8900381。
slab 通常需要足够真空层、固定底层和偶极修正。非对称界面可设置 LDIPOL = .TRUE.、IDIPOL = 3，具体方向要和表面法向一致。金属表面可用适合金属的展宽进行优化，再做更严格的静态能量对照；半导体或绝缘表面需检查费米能级和带边位置。
参考 INCAR 示例可包含 ENCUT、EDIFF、EDIFFG、ISMEAR、SIGMA、LREAL、LDIPOL、IDIPOL 和 LASPH。参数值只作为起始设置，最终以体系收敛测试为准。对带净电荷模型，还要说明背景电荷、真空层和能量修正边界。
若使用隐式溶剂和显式阳离子共同描述界面，介电常数、空穴边界和显式水团簇要固定在同一组对照中。隐式溶剂改变长程屏蔽，显式阳离子改变局域配位；二者同时出现时，不能再把能量变化只归因于 K⁺ 或 Na⁺。若只想看局域离子作用，可先在真空加水层模型中完成 K⁺/Na⁺ 对照，再把隐式溶剂作为第二层边界测试。
结构目录可分为 clean-water、clean-water-K、clean-water-Na、COOH-water、COOH-water-K、COOH-water-Na、OOH-water、OOH-water-K 和 OOH-water-Na。每组先优化，再做静态输出 CHGCAR、LOCPOT、DOSCAR 和 Bader 所需电荷材料。
同一张记录表中可列出 Eads、ΔEstab、ΔG*COOH、ΔG*OOH、M-Owater、M-Oads、关键 O-H/O-O 键长、Bader 电荷、偶极矩和真空平台。几何量与能量量放在同一组记录中，才能判断阳离子效应来自静电稳定化、氢键改变还是直接接触。
如果 K⁺ 与 Na⁺ 给出相反趋势，应优先检查优化后的水合环境是否仍可比较。若 Na⁺ 形成更紧的水合壳层，而 K⁺ 更靠近吸附物氧端，二者改变 *COOH 或 *OOH 的方式并不相同。
对含 *COOH 的 CO₂RR 算例，还应把 *H 竞争吸附放在同一水层和阳离子环境下计算；对含 *OOH 的 ORR/OER 算例，则要同时保留 *OH 和 *O 端点。这样得到的不是孤立中间体的能量变化，而是阳离子对路径分叉和限制步骤的影响；阳离子效应只有进入同一反应路径时才具有催化含义，不能停在单个端点。
最终判断中，显式阳离子模型的可靠性来自成套对照，而不是单个含离子构型，并应同时说明表面模型、水层数量、阳离子位置、吸附物构型、能量读取口径和自由能修正方式。
要点：显式阳离子模型要同时定义阳离子、水层、吸附物和表面位置。
要点：*COOH 可用稳定化能拆分 K⁺/Na⁺ 对中间体的局域作用。
要点：*OOH 对 O-H 方向、O-O 键长和局域电场更敏感，需结合几何和电荷结果判断。
要点：无阳离子、K⁺ 和 Na⁺ 端点应在同一 slab、同一水层数量和同一能量口径下比较。
