隐式溶剂模型在催化计算中到底改了什么？

引言

电催化 DFT 计算常从真空 slab 开始，但真实反应发生在电解液界面。把隐式溶剂打开后，很多中间体自由能会变化，尤其是 *OH、*OOH、*COOH、*OCHO 这类极性较强或局部电荷重排明显的物种。隐式溶剂首先改变的是连续介电环境对电荷和偶极的响应，不是自动生成一层真实水分子。

因此，隐式溶剂模型的作用应拆开看：介电屏蔽会稳定极性态，表面电势会改变界面电场的解释，空穴定义会影响溶剂接触区域，中间体自由能则取决于吸附态和清洁表面对溶剂响应的差值。若把这些统称为“溶剂让结果更真实”，就会漏掉关键计算边界。

图1 电催化界面模型的层级示意，从全显式水环境到连续介质近似。图中展示了不同模型对溶剂和电极响应的处理层级。DOI: 10.1002/smtd.202501542。

介电环境究竟改了什么？

它怎样稳定极性中间体？

在隐式模型中，溶剂被看作可极化的连续介质。吸附态产生偶极或局部电荷重排时，周围介质会通过极化响应降低体系自由能。*COOH、*OCHO、*OOH 通常比 *CO 更容易受到这种响应影响，因为前者含有更明显的 C-O、O-H 或 O-O 极性片段。隐式溶剂对不同中间体的修正不是等量平移，而是取决于电荷分布和偶极变化。

ΔG_solv = G_implicit – G_vacuum。式(1)

式(1)定义了一个最直接的溶剂化修正。若 *COOH 的 ΔG_solv 比 *CO 更负，只能说明该吸附态在连续介质中获得更强稳定；它不能单独证明形成了某个特定氢键，也不能说明水分子以某个角度排列。隐式模型给出平均环境响应，不写出局部水结构。

这也是解读图1时应抓住的层级差异：连续介质模型提高了界面环境描述，但它牺牲了水分子的构型自由度。若反应步骤依赖具体质子转移、氢键桥或水解离，隐式溶剂只能作为背景修正，不能替代显式水构型。凡是要讨论具体水参与方式，都需要显式水或采样证据补充。

VASP 参数应怎样理解？

VASPsol 一类隐式溶剂设置通常通过 LSOL 打开，并用溶剂介电常数、空穴参数和离子强度相关参数描述连续环境。不同版本和编译方式支持的关键词会有差异，实际使用前应以当前程序文档和编译说明为准。对水环境，介电常数常按室温水的数量级设置；但空穴边界、表面电荷和离子项仍需按研究对象测试。

参数设置时应先做清洁 slab 对照，再做吸附态对照。若隐式溶剂让清洁表面总能大幅波动，或导致结构优化异常，就不宜直接把吸附态自由能变化当成溶剂效应。隐式溶剂修正必须用同一套模型口径做差，不能只看吸附态单点能。

更清楚的计算步骤是：先在真空中优化清洁表面和吸附态；再用相同几何做隐式静态对照；必要时在隐式环境下重新优化关键吸附态；最后比较真空与隐式环境下的相对自由能。这样可以区分电子环境修正和结构重排贡献。

表面电势为什么会跟着变？

平面平均电势不再只是一个真空平台问题

真空 slab 中常通过 LOCPOT 平面平均静电势寻找真空平台，再计算功函数。加入隐式溶剂后，界面区域不再是简单真空，电势分布会受到介电边界、表面偶极和可能的离子响应影响。此时讨论表面电势，应明确比较的是模型内部电势变化，还是要映射到电化学电位标尺。

图2 展示了显式/隐式混合模型对电极-电解液界面的描述方式。它提示我们：连续介质可以帮助描述长程极化，但局部水层、离子位置和吸附物取向仍会影响界面电势。隐式溶剂能改变电势分布，但不等于已经完成恒电位模拟。

图2 显式溶剂与混合显式/隐式模型示意。图中展示了水分子、离子与连续介质可共同影响局部电场和吸附环境。DOI: 10.1002/smtd.202501542。

如果进一步改变 NELECT 或加入背景电荷，表面电势会同时受到电子数、补偿电荷、介电响应和界面电容影响。没有电位校准时，不能把一个带电计算直接命名为某个 RHE 电位；更合适的表述是该电荷状态下的电子化学势或界面电势发生变化。

它怎样影响吸附自由能？

表面电势改变后，带偶极或电荷转移明显的中间体会受到不同影响。若 *OCHO 的 O 端靠近表面且局部负电性较强，介电环境可能比对 *CO 更明显地稳定它。若 *H 基本呈中性吸附，溶剂修正可能较小。吸附自由能变化来自清洁表面与吸附态对介电环境响应的差值。

ΔG_ads^implicit = G_slab+X^implicit – G_slab^implicit – G_X + ΔG_corr。式(2)

式(2)中的 ΔG_corr 包含零点振动能、熵、电位项和 pH 项等修正，具体按反应体系决定。需要注意的是，气相分子参照通常不在同一个溶剂空穴中计算；若把气相 CO₂、H₂O 或 H₂ 的处理口径改掉，应在自由能表中单独说明。

对同一反应路径，隐式溶剂可能改变限制步。例如 CO₂RR 中 *COOH 被稳定后，后续 *CO 脱附或进一步氢化可能成为主要台阶；OER 中 *OOH 若被更强稳定，过电位可能下降。判断时不能只报单个中间体降低多少，还要回到完整自由能路径。

中间体稳定性为什么会重新排序？

极性、暴露面积和取向都会参与

隐式溶剂修正的大小与中间体暴露给溶剂的表面积、偶极方向、局部电荷分布和表面屏蔽能力有关。同一种 *COOH，若 C 端吸附和 O 端吸附的暴露方向不同，溶剂响应也可能不同。隐式环境下的构型排序可能不同于真空下的构型排序。

图3 和图4 来自 Ag/H₂O 界面 CO₂ 吸附研究，显示界面水密度、取向和自由能面会随反应环境变化。虽然这是显式界面研究，但它能帮助理解隐式模型的边界：连续介质能描述平均极化，无法直接给出水密度峰和具体取向。

图3 Ag/H₂O 界面水密度和取向分布。图中展示了电位条件下界面水结构的变化。DOI: 10.1038/s41467-025-65630-1。

图4 CO₂ 吸附过程中的自由能面与水化结构变化。图中展示了显式界面对反应坐标和中间体稳定性的影响。DOI: 10.1038/s41467-025-65630-1。

如果隐式结果和显式水结果方向一致，可以说明极化稳定和局部水结构都支持该趋势；如果二者不一致，则说明具体水构型可能主导该步骤。此时应检查氢键数量、水分子取向、吸附物暴露方向和自由能差，而不是简单选择一个更有利的结果。

如何设置一组清楚的对照？

建议按四层对照推进：真空 slab、隐式溶剂静态、隐式溶剂优化、显式水或显式水加隐式背景。每增加一层，只改变一个主要因素。这样可以把介电修正、结构重排和局部水作用拆开，而不是把所有变化混成一个溶剂修正值。

同一组对照必须统一 slab 层数、真空层、偶极修正、覆盖度、吸附位点、磁性初值和收敛标准。对含过渡金属位点的体系，还要同步固定 U 值和 LASPH 设置。溶剂模型改变时，电子结构参数不能同时漂移。

若显式水层参与比较，水分子数也要保持一致，或用明确反应式处理水化学势。否则，ΔG 差值会混入额外水吸附和水层重排贡献。

结果应怎样落到催化判断？

先看完整路径，再看单个中间体

隐式溶剂最常见的误读，是看到某个中间体变低，就直接判断催化活性提高。正确顺序应是：先计算完整反应路径，再看限制步是否改变，最后判断过电位、限制电位或竞争反应是否同步变化。单个中间体稳定不等于整条路径改善。

ΔG(U,pH) = ΔG(0,0) – neU – k_BT ln(10) × pH × m。式(3)

式(3)表示电位和 pH 项进入自由能的一般形式，其中 n 是电子转移数，m 与质子参与数及参照标尺有关。若使用 RHE 标尺，pH 项的处理方式会与 SHE 标尺不同。隐式溶剂只改变环境能量项，不能自动替代电位和 pH 修正。

图5 展示介电环境对 OER 中间体自由能和过电位的影响。读这类图时应关注相邻台阶是否被不等量改变，而不是只看所有柱子有没有整体下降。介电修正真正改变催化判断的方式，是改写相邻中间体之间的自由能差。

图5 介电环境对 OER 中间体自由能和过电位的影响。图中展示了不同介电处理下自由能台阶与过电位变化。DOI: 10.1038/s42004-025-01508-z。

哪些结论需要加边界？

如果只做隐式溶剂，应把结论限定为连续介质环境下的趋势；若要讨论氢键、质子穿梭或水解离路径，就需要显式水结构或动力学采样补充。若体系带电或需要恒电位，应说明电荷补偿、电位校准和双电层处理方式。

对于材料筛选，隐式溶剂适合比较同一口径下大量候选的趋势；对于精细机理，尤其是 CO₂RR、碱性 HER、OER 中 *OOH 形成等步骤，最好加入代表水构型或多构型检查。隐式溶剂提高了环境描述层级，但没有消除界面模型的不确定性。

最终交付给自由能图的，不应只是“加溶剂后某态降低”。更有用的是列出真空与隐式环境下每个中间体的 ΔG、限制步、过电位或限制电位，并说明是否改变材料排序。

本章要点总结

要点一：隐式溶剂主要描述连续介电屏蔽和极化响应，不直接给出具体水分子构型。

要点二：吸附自由能变化来自清洁表面和吸附态对溶剂环境响应的差值，不能只看单个吸附态能量。

要点三：隐式溶剂可能改变极性中间体排序，但涉及氢键、质子转移和水解离时仍需显式水或采样证据。

要点四：电位、pH、带电模型和介电环境是不同层级的修正，不能互相替代。