VASP计算吉布斯自由能：金属表面吸附自由能中ZPE、振动熵与平动熵的校正链

一个Cu(211)表面CO₂吸附自由能的计算项目，初始阶段只取了VASP输出的电子能量E_elec做吸附能对比，得到的CO₂吸附自由能是-0.22 eV——暗示CO₂在Cu(211)台阶位上只是弱物理吸附。这个结果和程序升温脱附（TPD）实验中约420 K的脱附峰温度对不上：实验推算的吸附自由能应该在-0.5 eV量级。差距从何而来？

G_ads = E_elec + ZPE + U_vib – TS_vib是标准路径。三项校正中，ZPE和-TS_vib的贡献方向相反但量级相当——叠加之后的”漏掉”并不会互相抵消，反而会累积。

零点能（ZPE = Σhν_i/2）的贡献在CO₂这样的三原子分子上约0.15 eV。VASP通过IBRION=5或DFPT频率分析给出振动本征频率——该项目中CO₂在Cu(211)上的不对称伸缩模（ν_as≈2340 cm⁻¹）贡献了ZPE的主体。吸附后由于Cu-C键的形成，该模红移约80 cm⁻¹，ZPE变化约0.01 eV——量不大但不可忽略，因为在自由能差运算中ZPE从两侧各引入约0.15 eV量级的绝对值，差值误差被控制在0.01 eV以内。

振动的软模问题是这个项目中最难处理的校正项。吸附质在金属表面的平移自由度被冻结为振动模，产物就是一批10-50 cm⁻¹范围的低频软模。谐振近似对这些模式的熵贡献有系统性高估——在298K下，一个10 cm⁻¹的模式贡献约0.09 eV的-TS项，但物理上该模式更接近受阻平移而非谐振振动。这个项目将最软的三个振动模（<30 cm⁻¹）从谐振近似手动替换为二维平移振子模型：用表面扩散势垒计算的截止频率替代谐振频率代入平移配分函数。修正后，-TS吸附态的贡献降低了约0.15 eV，吸附自由能从-0.22 eV变为-0.40 eV——和TPD实验值的差距从0.28 eV缩小到0.10 eV。

平动熵的校正处理涉及标准态问题。气相CO₂在1 atm下的平动熵约0.66 eV（298K），而吸附态的平动自由度被完全冻结。在自由能差ΔG = G_ads – G_gas中，气相平动熵的全部0.66 eV都转化为脱附自由能的驱动力。若采用理想气体近似的标准化学势表达式μ(T,P) = E_elec + ZPE + U_vib – TS_vib + RT ln(P/P₀)，标准态从1 atm换算到实际反应分压（比如10⁻⁶ bar）后，压力修正项RT ln(P/P₀)贡献约-0.34 eV——CO₂吸附自由能进一步移至-0.74 eV。与实验值的残差只剩约0.05 eV，满足热力学精度要求。

准谐振近似通过phonopy接口做声子谱计算是处理固体体系吉布斯自由能的完整方案。phonopy读取VASP的力常数矩阵，通过简正坐标分析得到全布里渊区声子色散，再用配分函数累加得到F(T)曲线。这个项目后来用phonopy做了Cu(211)纯表面的自由能计算作为吸附自由能的基底扣除——体系的完备性由此闭合。