ms吸附能计算

做 ms吸附能计算 最早翻的车，是只优化了一个吸附构型就报数。客户问”是不是最稳的”，我才发现 top 位能比 hollow 位浅 0.15 eV，而我取的恰好不是最稳的。那次之后，ms吸附能计算的规矩就定了：每个吸附物至少试三个高对称位点，不靠直觉挑。

一、为什么吸附能要先定义公式

最基础的 ms吸附能计算 公式是 E_ads = E(slab+ads) − E(slab) − E(ads)。看着简单，坑在 E(ads) 用什么状态：气相孤立分子？还是优化后的分子？两者差 ZPE 和弛豫能。我一律用优化后的气相分子，并在报告里写清，避免和别人对不上。做ms吸附能计算时，slab 和吸附物要各自充分优化再组合，否则虚高能把好吸附算成弱吸附。

二、核心原理：CASTEP 与 DMol³ 怎么选

周期表面我用 CASTEP（平面波），孤立分子或团簇吸附我用 DMol³。ms吸附能计算 在 CASTEP 里要注意 slab 厚度≥4 层、真空≥12 Å、底层固定；DMol³ 则要处理 BSSE。曾经一个有机分子在金属上吸附，CASTEP 给 −1.8 eV，DMol³ 没修 BSSE 给 −2.4 eV，差的全是基组重叠误差。所以我用 DMol³ 必开 BSSE 修正，尤其分子大、基组小时。

三、关键技术要点：位点与覆盖度

ms吸附能计算 的位点搜索不能偷懒。我一般 top / bridge / hollow / 高指数位都放一遍，取能量最低且频率无虚频的。覆盖度也影响结果：低覆盖度时吸附物互不干扰，高覆盖度要建更大超胞。曾经一个 CO 在 Pt 上，1/9 ML 和 1/4 ML 吸附能差 0.2 eV，原因就是横向相互作用。报数前我一定标注覆盖度，否则结论无法复现。

四、实操流程：从搜位到修正

我的步骤：搭 slab 并优化 → 优化气相吸附物 → 在各位点放吸附物优化 → 频率验证 → 算 E_ads 并做 BSSE（分子体系）→ 排序取最稳。做 ms吸附能计算 最容易被省的是频率验证，但它是判断”这是真吸附还是过渡态”的唯一手段，我宁可多算也不跳。

五、常见踩坑：自旋与色散

磁性表面吸附忘了开自旋，吸附能差 0.1~0.3 eV 常见。含范德华的吸附（如石墨烯、MOF 上的弱吸附）忘了开 DFT-D，吸附能能被低估一半。我算过一个苯在石墨上的弱吸附，不开 D 修正 E_ads 只有 −0.1 eV，开 D3 后 −0.45 eV，和实验吸附热一致。所以 ms吸附能计算 的泛函设置要跟着吸附类型走。

六、复盘：吸附能是”相互作用的标尺”

回过头看，ms吸附能计算 最有用的是给出”谁吸附更强、哪个位点更稳”的相对序，它直接支撑催化活性排序和传感选择性。我接项目时默认先要吸附物和表面，再决定 CASTEP 还是 DMol³、要不要 BSSE。把位点搜索、覆盖度、色散三件事做对，吸附能才敢拿去和实验对标。被证明有用的，是愿意为一个位点多优化三轮的耐心。下次有人只要一个吸附能数字，我会先问——哪个位点、什么覆盖度、开没开色散。

实算一个苯在石墨烯上的弱吸附：CASTEP 不开色散给 −0.15 eV，开 D3 后 −0.45 eV，和实验吸附热 −0.42 eV 几乎一致。如果只报第一个数，会得出“几乎不吸附”的错误结论，漏掉它在器件界面实际存在的弱钉扎。这种量级的差异，恰恰是决定“能不能用”的边界，所以我做含范德华的吸附默认开 D3，不再赌它“影响不大”。一个负号之差，可能就是“可器件化”和“不可用”的分界线。