MS计算吸附能：建模到分析的完整经验指南

吸附能是表面科学和催化研究中最核心的计算量。无论是评估催化剂活性、筛选材料还是理解反应机理，吸附能计算都是绕不开的第一步。用Materials Studio的CASTEP模块算吸附能，流程上不难，但要把结果做准、做可靠，每一步都有经验可循。本文从大量实际项目中提炼出最实用的操作经验。

一、吸附能基础

1.1 公式

ΔE_ads = E(surface+adsorbate) – E(surface) – E(adsorbate)

• 负值表示吸附放热（有利）
• 正值表示吸附吸热（不利）
• 通常催化中有意义的范围：-0.5 ~ -3.0 eV

1.2 常见吸附质参考能量

吸附能计算需要单独优化吸附质分子的能量。以下是常见小分子的优化经验值（PBE/450eV）：

重要经验：O2和N2必须开自旋极化计算！分子态O2是三重态，如果用单态计算，能量偏高约1.5 eV，直接导致吸附能偏差巨大。

二、建模经验

2.1 表面模型选择

超胞选择经验：

• 覆盖率1/4 ML → (2×2)，最常用
• 覆盖率1/9 ML → (3×3)，低覆盖
• 覆盖率1/2 ML → (√2×√2)R45°，高覆盖
• 分子吸附（如CO、H2O）→ 至少(2×2)，避免周期镜像相互作用

2.2 吸附位点

以FCC金属(111)面为例，常见吸附位点：

经验：小原子（H、O、N）倾向空位吸附，分子（CO、H2O）倾向顶位或桥位。但不要假设——所有位点都算一遍才知道最稳定的。

2.3 初始吸附距离

吸附质到表面第一层原子的初始距离：

踩坑经验：初始距离太近（<1.2Å）会导致优化发散，太远（>4Å）则可能优化不到化学吸附态。我一般设1.8-2.0Å，基本万无一失。

三、CASTEP计算参数

3.1 优化参数

Task: Geometry Optimization
Quality: Fine
Energy cutoff: 450 eV
k-point spacing: 0.03 1/Å
SCF tolerance: 1.0E-6 eV/atom
Max force: 0.02 eV/Å
Max iterations: 300
Spin polarized: 看体系

3.2 色散校正

在MS中开启DFT-D：

• CASTEP → Calculation → Electronic → 勾选 “TS method” (Grimme)
• 对弱吸附体系（苯/有机分子在金属表面）影响0.3-0.8 eV
• 对强化学吸附（O/H在金属表面）影响<0.1 eV

经验判断：

• 有π电子体系 → 必须开D3
• 纯金属+小原子 → 可不开
• MOF/COF体系 → 必须开D3

四、完整计算流程

Step 1: 优化干净表面

1. 建slab模型，固定底层
2. CASTEP优化，记录E(surface)

Step 2: 优化吸附质分子

1. 在真空中建分子模型（15×15×15Å盒子）
2. 只需Γ点
3. 记录E(adsorbate)

Step 3: 构建吸附构型

1. 在MS中将分子放到吸附位点上
2. 调整初始距离和取向
3. 固定底层原子
4. CASTEP优化，记录E(surface+adsorbate)

Step 4: 计算吸附能

ΔE_ads = E(surface+adsorbate) – E(surface) – E(adsorbate)

Step 5: 验证与校正

1. 检查优化后的构型是否合理（键长、角度）
2. 频率计算验证无虚频（过渡态需有且仅有1个虚频）
3. ZPE校正：ΔE_ZPE ≈ 0.05-0.15 eV（通常使吸附能变正）
4. 色散校正（如未在计算中开启）

五、经验值参考

5.1 金属表面吸附能（PBE，eV）

5.2 常见问题判断

六、进阶分析

6.1 差分电荷密度

Δρ = ρ(surface+adsorbate) – ρ(surface) – ρ(adsorbate)

在MS中：

1. 对优化后的吸附体系做单点能计算
2. 对相同位置的干净表面和分子分别做单点能
3. 用MS的Volume Visualization分析电荷密度差

经验：差分电荷密度图能直观显示电荷从催化剂转移到吸附质（或反之），是论文中的重要配图。

6.2 Bader电荷分析

MS不直接提供Bader分析，但可以：

1. 导出CHGCAR格式
2. 用Henkelman组的Bader程序分析
3. 得到每个原子的精确电荷转移量

6.3 DOS分析

对比吸附前后的PDOS：

• 吸附质轨道与表面d带杂化 → 化学吸附
• 无明显杂化 → 物理吸附
• d带中心偏移 → 电子结构调控

七、项目报价参考

结语

吸附能计算的关键在于：一致的参数、合理的模型、全面的位点扫描。建议每次计算前先查文献值做参考，计算后与文献对标验证。如有需求，欢迎联系我们提供从建模到分析的完整技术服务。