ms计算自由能

我第一次用 ms计算自由能是把 CO 吸附在 Pt(111) 上，照搬一篇文献的模型跑 DMol³，吸附自由能算出来比人家低 0.28 eV，还以为自己体系更优。后来才发现人家加了隐式溶剂、做了全频率校正，我两样都没做。那次让我明白，ms计算自由能不是把 workflow 点完就交差，热力学链条上每一环都会原样写进最终数字。

一、为什么 ms计算自由能先选对模块

MS 里算自由能主要走两条路：DMol³ 做分子/团簇级、CASTEP 做周期表面。吸附体系我默认用 CASTEP 的周期 slab，单分子或小团簇反应才用 DMol³。曾经我硬用一个 2×2 的表面在 DMol³ 里用大超胞算，k 点没法设，误差直接上 0.2 eV。做 ms计算自由能时，先判断”周期与否”比挑泛函更该优先——表面和体相一定走 CASTEP，孤立分子走 DMol³，别混。

二、核心原理：热力学循环怎么搭

ms计算自由能的主线是热力学循环：ΔG_ads = E(slab+ads) − E(slab) − E(gas) + ΔZPE − TΔS。气相分子的自由能靠频率计算拿 ZPE 和熵，表面态用 DFT 总能量加 ZPE 修正。我用的是标准做法：气相分子单点+频率，表面吸附态频率只算吸附相关模式。关键是所有中间体要锚定同一参考态，否则不同步骤的自由能会整体平移，看不出谁限速。CASTEP 里我习惯用有限位移法算频率，DMol³ 用解析二阶导，二者 ZPE 差一般小于 0.02 eV。

三、关键技术要点：溶剂与色散不能省

液相或湿实验条件做 ms计算自由能，隐式溶剂（COSMO 或 DMol³ 的 solvation）一定要开，否则吸附自由能系统性偏高 0.2~0.4 eV。色散修正同理，Pt 表面 CO 吸附如果不补 DFT-D，弱相互作用算不准。我做过一组对比：关溶剂的 ΔG_ads 是 −0.45 eV，开溶剂变 −0.21 eV，而实验报道约 −0.18 eV。差出来的 0.24 eV 恰好就是溶剂贡献——ms计算自由能里溶剂不是可选项，是默认项。

四、实操流程：从建模到自由能

我的标准动作：确定体系与参考态 → CASTEP 建 slab 并充分弛豫（底层固定两层）→ 气相分子 DMol³ 单点+频率 → 逐个吸附构型优化（至少试 2 个位点）→ 吸附态频率取 ZPE 和熵 → 套热力学循环算 ΔG → 和实验对标。最容易被砍的是吸附态频率，但熵修正的可信度全靠它，我默认保留。

补一个 MOF 吸附甲烷的实算：用 CASTEP 算 ΔG_ads 在 298 K 下是 −0.12 eV（含 ZPE 和熵），只算总能不修正的话是 −0.31 eV，差出 0.19 eV。如果拿没修正的数字去解释吸附等温线，你会误判这个 MOF”吸附很强”，实际在室温下根本吸附不住。这之后我做任何 ms计算自由能的吸附体系，频率和溶剂都是默认打开，不再当作可选项。

五、常见踩坑：参考态不一致

ms计算自由能最常见的错是气相分子和表面态用不同方法、不同基组，导致 ΔG 里混进方法误差。我见过有人气相用 DMol³、表面用 CASTEP，中间差了 0.1 eV 的基组效应。另一个坑是温度，熵随 T 变化，常温算的 ΔG 不能直接套到高温反应。还有人把振动熵用爱因斯坦模型近似，结果在低频模式多的体系误差被放大。

六、复盘：自由能的价值在”可比”

回过头看，ms计算自由能最有用的是让不同体系能在同一把尺子下比较，而不是单给一个绝对值。我接这类项目默认先要参考态和实验条件，再决定算哪些校正项。把模块选型、溶剂、频率三件事做对，ΔG 才敢拿去和吸附实验对话。被证明有用的，是愿意为一个溶剂项多跑半小时的扎实。下次有人拿一个”总能”当自由能报给我，我会先问——你的 ZPE 和熵加了吗。