吸附自由能计算：从吸附能到表面反应自由能的完整路径

吸附自由能计算是电催化（OER/ORR/HER/CO₂RR）理论研究的基石。用DFT优化出一个吸附构型→算出吸附能→代入计算氢电极（CHE）模型拿到反应自由能台阶图→预测过电位——这条标准路径看起来皮实，但实际上每一步都有隐式假设，而吸附熵、覆盖度效应和溶剂化修正三条腿的任何一条不稳，整条自由能链就偏了。

吸附能≠吸附自由能：熵贡献的三种来源

DFT算出来的吸附能是0K静态电子能之差：E_ads = E(slab+ads) – E(slab) – E(ads_gas)。但室温下的吸附自由能需要三个额外的熵修正项。

分子平动熵损失是最大头。气相分子在气相中有完整的平动和转动自由度，对应标准摩尔熵可达150-200 J/mol·K。吸附在表面上后，这些平移和转动被限制为二维扩散/受阻转动和低频振动——熵损失约100-150 J/mol·K，对应-TΔS贡献约+30到+45 kJ/mol。这意味着就算吸附是放热的（ΔH<0），ΔG可能只勉强为负甚至为正——在低覆盖度下，分子的吸附是一个自发性被熵减”抵消”放热的过程。

表面振动熵增益通常被忽略，因为从slab→slab+ads新增的低频振动模式（<200 cm⁻¹，对应吸附质-表面的平动/振动耦合）会产生约+5到+15 J/mol·K的熵增益。虽然量级只有平动熵损失的1/10，但它部分缓冲了净熵减。

构型熵在大分子（C3+有机物）吸附中不可忽略。一个丁烷分子在Pt(111)表面可能有多个几乎能量简并的吸附构型（gauche/反式、不同吸附位点组合）——每个构型贡献k_B·ln(N)的熵。如果存在3个~等能构型，构型熵贡献约9 J/mol·K——不大但积少成多。

覆盖度效应：从稀释极限到高压工况

CHE模型的基础假设是表面吸附质之间无相互作用（Langmuir型），因此ΔG_ads不随覆盖度变化。这个假设在低覆盖度下（θ<0.1 ML）成立，但实际催化反应常在高覆盖度下运行——ORR在Pt(111)上的工作覆盖度接近0.5-1.0 ML OH*/O*。

覆盖度对吸附自由能的影响来自两个渠道：吸附质之间的直接静电/空间排斥（使吸附能降低）和间接的电子效应（吸附诱导的表面电子态d带中心位移）。以Pt(111)上O的吸附为例：0.11 ML覆盖度下ΔG_ads(O) = 0.95 eV；0.25 ML时降到0.68 eV；0.50 ML时仅剩0.42 eV。0.27 eV的覆盖度效应不能忽略——它足以把ORR的过电位预测从0.45 V推高到0.72 V。

正确处理覆盖度效应需要在每个目标覆盖度下重新做覆盖度依赖的自由能计算——即构建不同大小的超胞、放不同数量的吸附质→逐一优化→提取覆盖度相关的ΔG_ads(θ)。

溶剂化修正：水分子不只是介质

电催化环境中的吸附自由能计算还有一个核心挑战：界面水分子网络。水分子以氢键网络的形式覆盖在电极表面，它们不是被动的”溶剂背景”——水分子直接参与中间体的稳定化（如通过氢键稳定OH和OOH），有时直接参与质子转移。

处理界面水的策略有三层：第一层，在隐式溶剂模型（VASPsol）下计算吸附自由能，可以捕获连续介质的静电溶剂化效应。VASPsol对OH和OOH的稳定化能量在0.1-0.3 eV级别，这对自由能台阶的能级排布有关键影响。

第二层，加入1-2层有序的显式水分子（冰状水双层），用DFT优化它们的位置和取向。这个方案对HER的酸性机理修正不大（H*的溶剂化效应弱），但对OER——尤其是O-O偶联步骤——显式水的氢键稳定可以影响键形成/断裂的能垒0.2-0.3 eV。

第三层，在显式水层的基础上跑AIMD（从头算分子动力学）来采样水分子网络的动态重排——这是精度最高的方案，但计算成本巨大（每自由能点需~1000核时）。

自由能台阶图：误差从哪来一目了然

将吸附自由能代入CHE模型，画出反应自由能台阶图，每个台阶的ΔG值与理想催化剂（所有台阶等距，ΔG=1.23 eV for ORR）的偏差直接决定了过电位的理论极限。但自由能台阶图中每个台阶的误差会在台阶之间累积——如果*OOH的吸附自由能高了0.2 eV，OER的过电位预测就低了0.2 V——这个灵敏度要求每一步吸附自由能计算的不确定度都必须压在<0.1 eV以内。

吸附自由能计算的最终检验标准，是与实验Tafel斜率和起始电位的对比——不是看一个点对不对，而是看整条自由能链的行为在覆盖度-电势坐标系下能否复现实验的动力学特征。

更多内容请访问 https://www.keyanxueshu.com/