1. 吸附能基础概念
做表面计算的人，第一个绕不开的概念就是吸附能。我刚开始接触这个领域时，导师丢给我一句话：「你去算算这个分子在催化剂表面吸附得牢不牢。」当时我一脸懵——什么叫「牢不牢」？后来才明白，这个「牢不牢」就是吸附能。
什么是吸附能
吸附能，说白了就是衡量一个分子（或原子）跟表面结合得有多紧的物理量。你想想看，一个气体分子飘到固体表面，它可能被「粘」住，也可能弹开。这个「粘住」的过程会释放能量，释放得越多，说明结合越强。
我习惯用一个比喻来理解：就像磁铁吸铁钉。磁力越强，你要把铁钉拔下来花的力气就越大。吸附能就是这个「拔下来的力气」，只不过在原子尺度上，我们用量子力学算出来的能量差来表示。
核心定义：吸附能 = 吸附后体系的总能量 − (吸附前表面的能量 + 吸附前自由分子的能量)
通常为负值，越负表示吸附越强。
物理吸附 vs 化学吸附
这里有个关键区分，我在项目中吃过亏，所以特别想讲清楚。
| 特征 | 物理吸附 | 化学吸附 |
|---|---|---|
| 作用力 | 范德华力（分子间作用力） | 化学键（共价键/离子键） |
| 吸附能大小 | 通常 < 0.3 eV | 通常 > 0.5 eV，可达几个eV |
| 吸附距离 | 较远（~3-4 Å） | 较近（~1-2 Å） |
| 电子结构变化 | 几乎不变 | 明显变化，有电荷转移 |
| 可逆性 | 容易脱附 | 较难脱附 |
为什么会这样？物理吸附靠的是范德华力，说白了就是分子间那种「若即若离」的吸引力，很弱。化学吸附则是分子跟表面真正「牵手」了——电子云重叠，形成化学键。
我记得有一次算CO在Pt表面的吸附，一开始用的泛函没加色散修正，算出来吸附能只有0.2 eV。我当时还纳闷：文献上明明说是化学吸附啊？后来加了DFT-D3修正，吸附能直接跳到1.8 eV。嗯，这就是物理吸附和化学吸附的区别在计算中的体现——选不对方法，结果天差地别。
实战经验：判断是物理吸附还是化学吸附，我一般看两个指标：一是吸附能绝对值是否大于0.5 eV，二是吸附后分子键长有没有明显变化。如果分子键长被拉长了0.1 Å以上，基本可以断定是化学吸附。
吸附能的计算公式
公式其实很简单，但细节决定成败。
E_ads = E_total − (E_surface + E_molecule)其中：
- E_total：吸附体系的总能量（表面+吸附质一起算）
- E_surface：干净表面的能量
- E_molecule：自由分子的能量（放在一个大盒子中算）
注意！这里有个坑——基组重叠误差（BSSE）。我曾经算一个弱吸附体系，没做BSSE修正，结果吸附能是正的（表示不吸附），但实验上明明能吸附。后来用counterpoise方法修正后，才得到合理的负值。
避坑指南：我曾经在算H₂O在TiO₂表面的吸附时，发现吸附能总是偏大。排查了三天，最后发现是自由分子的计算盒子大小跟表面计算不一致。记住：自由分子的计算盒子至少要15 Å以上，且要用相同的截断能和k点设置。
吸附能的物理意义
吸附能不仅仅是一个数字。它告诉你：
- 反应活性：吸附能越负，分子越容易被活化
- 覆盖度效应：不同覆盖度下吸附能会变化，这决定了催化剂的活性位点
- 选择性：两种分子在同一表面的吸附能差异，决定了谁优先吸附
我举个例子。在计算CO₂还原时，*CO是一个关键中间体。如果*CO在催化剂表面吸附太强（吸附能太负），它就不容易脱附，催化剂会被毒化；如果吸附太弱，又无法有效活化。所以做催化剂设计时，我们追求的是一个「恰到好处」的吸附能——Sabatier原理说的就是这个道理。
知识体系框架
下面这张图是我自己总结的吸附能知识体系，帮你理清思路：
这张图把吸附能的核心内容串起来了。你从定义出发，搞清楚分类，掌握计算方法，理解物理意义，最后应用到实际问题中。每一步都有坑，但每一步也都有规律可循。
我的建议：刚开始做吸附能计算时，先拿一个简单的体系练手，比如CO在Pt(111)表面。这个体系文献数据多，容易验证。算通了一个，后面就顺了。
好了，吸附能的基础概念就讲到这里。记住三个关键词：能量差、作用力类型、物理意义。后面我们会一步步深入到具体的计算实操中去。
公众号:蓝海资料掘金营，微信deep3321