4. 吸附位点识别：高对称位点与判据

吸附位点的识别，说白了就是搞清楚分子到底“坐”在表面的哪个位置。这步要是搞错了，后面的吸附能计算全是白费功夫。我刚开始做表面吸附时，就吃过这个亏——选错了位点，算出来的能量差了好几个eV，还以为是自己的计算参数设错了。

今天咱们就聊聊怎么识别这些位点，以及怎么判断它们是不是真的稳定。

4.1 高对称位点：顶位、桥位、空位

在晶体表面，原子排列是有规律的。分子吸附时，通常会选择几个“特殊”的位置。这些位置因为对称性高，能量上往往更有利。常见的就三种：顶位、桥位、空位。

顶位（Top Site）：吸附原子正对着表面某个原子的正上方。说白了，就是“骑”在一个原子上。

桥位（Bridge Site）：吸附原子在两个表面原子的中间。你可以想象成“跨”在两个原子上。

空位（Hollow Site）：吸附原子在三个或四个表面原子围成的空穴中心。这个位置通常是最稳定的，因为配位数最高。

我举个例子，拿最常见的fcc(111)面来说：

• 顶位：正对着一个表面原子
• 桥位：在两个相邻表面原子之间
• fcc空位：在三个原子围成的三角形中心，但下方第二层没有原子
• hcp空位：同样在三个原子围成的三角形中心，但下方第二层有原子

你想想看，fcc和hcp空位就差那么一点点，但吸附能可能差不少。我有个项目就是研究CO在Pt(111)上的吸附，fcc空位比hcp空位稳定了大概0.2 eV。别小看这0.2 eV，在催化反应里，这足以改变整个反应路径。

4.2 几何判据：距离与配位

怎么判断一个位点是不是真的“吸附”了？几何上就看两点：距离和配位。

距离判据：吸附原子与最近表面原子的距离，应该接近它们之间的共价键长或离子键长。如果距离太远，那根本就没成键，只是物理吸附。

我个人习惯用这个经验值：吸附原子与表面原子的距离，一般比它们之间的共价键半径之和大个0.1-0.3 Å。比如O原子在Pt表面，Pt-O键长大约2.0 Å左右。如果算出来是3.0 Å，那肯定不对。

配位判据：吸附原子周围有几个表面原子？配位数越高，吸附越强。顶位配位数是1，桥位是2，空位是3或4。一般来说，空位吸附最强，桥位次之，顶位最弱。

但这里有个坑——不是所有体系都遵循这个规律。我记得有一次算H在Pd(111)上的吸附，结果发现顶位反而比桥位稳定。后来一查文献，原来是因为H原子太小，在桥位上反而有应力。嗯，这里要注意，小原子（H、Li等）经常不按套路出牌。

4.3 电子判据：电荷密度与态密度

几何判据只能告诉我们“有没有吸附”，但电子判据能告诉我们“为什么吸附”。

电荷密度差：这是我最常用的方法。公式很简单：

Δρ = ρ(吸附体系) - ρ(表面) - ρ(孤立分子)

如果Δρ在吸附原子和表面原子之间出现明显的电荷积累（红色区域），说明形成了化学键。如果只是微弱的变形，那可能是物理吸附。

态密度（DOS）分析：吸附前后，表面原子的d轨道态密度会发生变化。如果吸附后，d轨道在费米能级附近出现新的峰，或者原有峰发生劈裂，说明有成键。

我举个例子，CO在过渡金属表面的吸附。CO的5σ轨道会向金属的d轨道捐赠电子，而金属的d轨道会向CO的2π*轨道反馈电子。在DOS图上，你会看到CO的5σ峰往低能级移动，而金属的d轨道在费米能级附近出现新的反键态。

为什么会这样？说白了，这就是“捐赠-反馈”机制。你想想看，CO和金属之间不是单方向的电子转移，而是双向的。这个机制解释了为什么CO在Pt、Pd等贵金属上吸附特别强。

4.4 实战：如何系统识别吸附位点

好了，理论说完了，咱们来点实际的。我一般按这个流程走：

1. 建表面模型：切出合适的表面，比如fcc(111)、bcc(110)等。记得加真空层，至少15 Å。
2. 标记高对称位点：用VESTA或ASE标出顶位、桥位、空位的位置坐标。
3. 放置吸附原子：在每个位点上方1.5-2.0 Å处放一个原子。
4. 结构优化：固定下面几层原子，只优化吸附原子和表面顶层。
5. 分析结果：看最终位置、吸附能、电荷密度差、DOS。

这里有个小技巧：我习惯先用DFT+U算一遍，再用杂化泛函验证关键位点。因为杂化泛函太贵，全算一遍不现实。但关键位点的能量差，杂化泛函会更准。

下面是我画的一张流程图，帮你理清思路：

这个流程我用了好几年，基本没出过问题。当然，具体体系还得具体分析。比如二维材料（石墨烯、MoS₂）的吸附位点就和金属表面不太一样，但思路是相通的。

---

公众号:蓝海资料掘金营，微信deep3321