什么是吸附能公式？
在表面化学、催化科学及材料科学领域，吸附能（Adsorption Energy）是一个至关重要的物理量。它量化了吸附质（Adsortate）与吸附剂（Adsorbent/Surface）之间相互作用的强弱。理解吸附能公式不仅是进行第一性原理计算（如DFT）的基础，也是解释实验吸附等温线、预测催化剂活性的关键。
简单来说，吸附能描述了当气体分子或液体分子附着在固体表面时，系统能量的变化过程。这一过程通常伴随着能量的释放（放热）或吸收（吸热）。对于绝大多数催化反应而言，反应物在催化剂表面的吸附是反应发生的第一步，因此吸附能的大小直接决定了反应的难易程度。
⚙️ 核心定义
吸附能是指单位面积或单位摩尔吸附质在表面形成吸附层时所释放或吸收的能量。它是衡量吸附作用强度的直接指标。
⚡ 重要意义
在火山型曲线（Volcano Plot）中，吸附能是描述催化剂活性趋势的核心描述符（Descriptor）。Sabatier原理指出，最佳催化剂的吸附能应处于适中范围。
? 应用领域
广泛应用于氢气存储材料开发、CO2捕获、电催化析氢（HER）、析氧反应（OER）以及环境污染物去除等领域。
吸附能公式的深度解析
在密度泛函理论（DFT）计算中，吸附能公式的表达形式最为常见且基础。它是连接微观电子结构与宏观物理性质的桥梁。以下是标准吸附能公式及其各项的物理意义。
Eads：吸附能。若结果为负值，表示吸附过程放热，吸附稳定；若为正值，表示吸热，吸附不稳定。
Etotal(Surface + Adsorbate)：吸附体系（表面+吸附质）的总能量，单位通常为 eV。
Etotal(Surface)：纯净表面模型的总能量，单位通常为 eV。
Etotal(Adsorbate)：孤立吸附质分子（或原子）的总能量，单位通常为 eV。
实验吸附热的关联
需要注意的是，DFT计算得到的吸附能通常是在0 K下的电子基态能量差。而在实验条件下（如298 K），我们需要考虑热力学修正。实验测得的吸附热（Qst）与计算吸附能的关系如下：
其中 Δ(H - TS) 代表了焓变和熵变的贡献。对于气体吸附，熵的损失通常较大，因此实验测得的吸附热往往比单纯计算的电子能量差要小（绝对值）。
多分子吸附的修正公式
当表面覆盖度（Coverage, θ）较高时，吸附质之间存在相互作用，此时简单的吸附能公式需要进行修正。对于 N 个吸附质分子：
这种修正确保了吸附能与覆盖度 θ 的对应关系，是研究高压吸附条件下的关键步骤。
物理吸附 vs 化学吸附
根据吸附能的大小及作用机制的不同，吸附现象主要分为物理吸附（Physisorption）和化学吸附（Chemisorption）。理解两者的区别对于选择合适的吸附能公式模型至关重要。
| 特征 | 物理吸附 | 化学吸附 |
|---|---|---|
| 作用力 | 范德华力（Van der Waals forces） | 化学键力（共价键、离子键、配位键等） |
| 吸附能范围 | 4 - 40 kJ/mol (≈ 0.04 - 0.4 eV) | 40 - 800 kJ/mol (≈ 0.4 - 8.0 eV) |
| 选择性 | 无选择性，任何气体均可被任何固体吸附 | 具有高度选择性，仅特定气体-固体组合发生 |
| 层数 | 可形成多层吸附 | 通常仅为单层吸附 |
| 温度影响 | 低温有利于吸附，高温脱附 | 可能需要活化能，高温下更易发生 |
| 可逆性 | 高度可逆 | 往往不可逆，或难以完全脱附 |
网友还关心：如何区分两者？
在实际计算中，可以通过观察吸附质与表面原子间的电荷转移（Bader Charge Analysis）和键长变化来辅助判断。如果吸附能较小且无显著电荷转移，多为物理吸附；若吸附能较大且伴随明显的电子重排和键长缩短，则为化学吸附。近年来，针对物理吸附中范德华力（vdW）的修正（如DFT-D2, DFT-D3方法）成为研究热点，因为标准泛函往往严重低估物理吸附能。
如何正确计算吸附能？
计算吸附能看似简单，但在实际操作中极易出错。以下是基于DFT计算的标准化流程指南。
模型构建要点
- 表面模型：使用超胞（Supercell）方法构建表面，确保真空层厚度足够（通常 > 15 Å），以避免周期性镜像相互作用。
- 吸附位点：需测试多种可能的吸附位点（如Top, Bridge, Hollow等），选择能量最低（吸附能最负）的构型作为稳定构型。
- 吸附质：确保孤立分子的自旋极化设置正确，特别是对于O2、NO等顺磁性分子。
参数设置关键
- 截断能（Cut-off Energy）：三个计算（表面、分子、复合系统）必须使用完全相同的截断能，以确保能量相减的准确性。
- K点网格：表面计算需使用密集的K点网格，而孤立分子计算通常只需Γ点或极稀疏的网格。注意在VASP等软件中，需对分子计算进行修正或使用大盒子。
- 收敛标准：能量收敛标准建议设为 10-5 eV，力收敛标准设为 0.01 eV/Å，以确保几何结构充分弛豫。
误差校正方法
- 零点能校正（ZPE）：振动频率计算后获得的零点振动能需加到电子能量上，公式：Ecorrected = EDFT + 0.5·Σhνi。
- 基组叠加误差（BSSE）：对于小基团或弱相互作用体系，需使用Counterpoise方法校正BSSE，否则吸附能会被高估。
- 熵校正：若需与实验298K数据对比，需通过振动频率计算获得焓（H）和熵（S）的贡献，计算 Gibbs 自由能变化 ΔG。
影响吸附能的关键因素
理解哪些因素调控吸附能的大小，是设计高效吸附剂和催化剂的核心。主要因素包括电子结构、几何结构及外部环境。
① d带中心理论
对于过渡金属，吸附能与金属d带中心（d-band center, εd）密切相关。d带中心上移，反键轨道填充减少，吸附能增强（结合更牢）；d带中心下移，吸附能减弱。
② 晶面效应
不同晶面具有不同的原子配位数和表面能。通常，低指数密排面（如Pt(111)）吸附能较弱，而高指数台阶面（如Pt(211)）吸附能较强，活性位点更多。
③ 缺陷工程
空位（Vacancy）、掺杂（Doping）和应变（Strain）能显著改变局部电子结构，从而调控吸附能。例如，石墨烯中的氮掺杂可增强对O2的吸附能。
④ 溶剂化效应
在电催化中，溶剂分子（如水）会与吸附质竞争表面位点。隐式溶剂模型（如VASPsol）可模拟溶剂对吸附能的屏蔽和稳定作用，使计算更贴近实验。
吸附能研究前沿动态
近年来，随着计算能力的提升和机器学习的发展，吸附能的研究进入了一个新的阶段。
2010s 初期：d带中心理论的成熟
Nørskov等人提出的d带中心模型成为预测过渡金属催化剂吸附能和活性的标准工具，极大地加速了催化剂筛选过程。
2015s 中期：范德华力的精确描述
DFT-D2/D3等色散校正方法的广泛应用，解决了传统泛函无法准确描述物理吸附能的难题，使得多孔材料（如MOFs, COFs）的气体吸附模拟成为可能。
2020s 至今：机器学习与高通量筛选
利用机器学习势函数（Machine Learning Potentials）和描述符，研究者可以在毫秒级时间内预测数百万种材料的吸附能，实现了从“试错法”到“理性设计”的转变。
未来趋势：动态吸附与原位模拟
从静态吸附能计算向动态过程模拟发展，结合原位光谱数据，揭示真实反应条件下吸附能的动态演化规律。
网友最关心的吸附能公式问题
我们在日常科研交流中收集了以下高频问题，并结合吸附能公式原理进行详细解答。
正值表示吸附过程在能量上是不利的（吸热）。可能的原因包括：1. 吸附位点不是全局或局部能量最低点，尝试其他吸附构型；2. 表面模型存在缺陷或张力，导致表面能过高；3. 吸附质分子本身不稳定；4. 未考虑熵效应，在高温下熵项可能主导，导致吉布斯自由能为正，但电子能量差仍可能为负。建议检查几何优化是否充分收敛。
吸附能（Eads）通常指0 K下电子基态的能量差，而吸附热（Qst）是实验测量的热力学量，包含内能变化和PV功等。在恒温恒压下，两者关系涉及焓变。对于固体表面吸附，通常近似认为 Qst ≈ -Eads，但精确计算需加入热校正项（ZPE, H, TS）。
随着覆盖度增加，吸附质之间的相互作用（排斥或吸引）会增强。通常，高覆盖度下吸附质间存在排斥作用，导致吸附能的绝对值减小（结合变弱）。因此，在报告吸附能时，必须注明对应的覆盖度 θ。
自相互作用误差在局域泛函（如LDA, GGA）中较明显，会导致电子过度离域，影响吸附能准确性。对于强关联体系（如含d/f电子的材料），建议使用杂化泛函（如HSE06）或DFT+U方法，以修正电子自相互作用，获得更可靠的吸附能。