3. 第一性原理方法核心工具箱

表面吸附计算的核心矛盾在于：我们关心的是几个原子在表面上的成键行为，但周期性边界条件迫使我们同时处理数十甚至上百个原子，还要兼顾金属表面的电子离域、吸附质的局域成键以及真空层的静电收敛。第一性原理工具箱中的每一项选择——泛函、赝势、基组、截断能、k 点——都会沿着这条链条放大或衰减误差。本节按“框架—泛函—赝势—基组”的顺序，把每一项的作用位置、可调参数和对吸附能的影响讲清楚。

一、DFT 基本框架：把多体问题变成可解的自洽场

Kohn–Sham 框架把相互作用的电子体系映射为一组无相互作用电子在有效势中运动，有效势包含外势、Hartree 势和交换关联势三部分。实际计算中，我们求解的是 Kohn–Sham 方程：

[-½∇² + V_ext(r) + V_H(r) + V_xc(r)] ψ_i(r) = ε_i ψ_i(r)
ρ(r) = Σ_i f_i |ψ_i(r)|²
V_H(r) = ∫ ρ(r')/|r-r'| dr'
V_xc(r) = δE_xc[ρ]/δρ(r)

求解流程是：给定初始 ρ → 构造 V_H 和 V_xc → 解本征方程得到 ψ_i → 重新占据得到新 ρ → 混合后迭代，直到能量和力收敛。对表面吸附而言，这个自洽循环里有两个关键量：一是总能差 ΔE_ads = E_slab+ads − E_slab − E_ads，二是吸附质与表面原子上的力，用于判断是否达到平衡构型。由于 ΔE_ads 是两个大数相减，总能收敛必须比吸附能精度更严，通常要求总能收敛到 1 meV/atom 量级，才能保证吸附能误差在 10 meV 以内。

二、交换关联泛函：精度与成本的第一次分岔

交换关联泛函是 DFT 中唯一需要近似的部分，也是表面吸附计算中误差的主要来源。常用泛函可分为三个层级：

• LDA：只依赖局域电子密度，对金属体相尚可，但会高估吸附能、低估键长，对分子吸附尤其容易给出过强的结合。
• GGA（PBE、PW91、RPBE）：引入密度梯度，是表面吸附的默认选择。PBE 对金属表面和化学吸附较为均衡；RPBE 通过修正交换增强项，通常给出更弱的吸附能，更接近实验值，但会略微高估键长。
• meta-GGA 与杂化泛函：如 SCAN、HSE06，引入动能密度或精确交换，对 CO、O₂ 等分子的吸附能改善明显，但计算成本显著增加，且对金属体系的收敛更困难。

选择泛函时，应先明确目标：若只比较同一表面上不同吸附位的相对稳定性，PBE 通常足够；若要给出与实验可比的绝对吸附能，需考虑 RPBE 或做 van der Waals 修正。对于物理吸附体系（如苯在石墨烯上），必须加入色散修正（DFT-D3、optB86b-vdW 等），否则结合能可能被低估数倍。

三、赝势与投影缀加波方法：把芯电子“冻结”起来

全电子计算需要处理 1s 等深能级，代价极高。赝势方法将芯电子与核一起冻结为一个有效离子实，只显式处理价电子。投影缀加波（PAW）方法在赝势框架内引入一个线性变换，把赝波函数映射回全电子波函数，从而在保持价电子精度的同时恢复芯区信息。对表面吸附而言，PAW 的优势在于：过渡金属的 d 电子和吸附质的 p 电子都能被准确描述，且磁性和自旋极化处理自然。

使用 PAW 时需注意三个参数：一是价电子配置，例如 Fe 的 3d⁶4s² 与 3d⁷4s¹ 会给出不同的磁矩和吸附能；二是截断能，PAW 数据集通常附带推荐值（如 400–500 eV），但表面体系建议做收敛测试；三是是否包含半芯态（如 3s、3p），对早期过渡金属和碱金属影响较大。一个实用的检查方法是：把吸附能随截断能的变化画出来，当从 400 eV 到 600 eV 的变化小于 10 meV 时，可认为截断能收敛。

四、基组选择：平面波、数值原子轨道与高斯基组

基组决定了波函数的展开方式，也直接决定了计算成本和可处理的体系规模。表面吸附中最常用的是平面波基组，因为它天然满足周期性边界条件，且收敛行为由单一参数——截断能——控制。平面波基组的缺点是：真空层需要足够厚（通常 ≥15 Å），否则波函数在真空中的衰减会被人为截断；同时，过渡金属的局域 d 电子需要较高的截断能。

另一类选择是数值原子轨道基组（如 SIESTA、OpenMX），它把基组局域在原子附近，适合大体系，但存在基组叠加误差（BSSE），且收敛行为不如平面波平滑。高斯基组（如 CRYSTAL）在分子吸附中精度高，但对金属表面的周期性处理较复杂。对大多数表面吸附问题，平面波 + PAW 是稳妥的默认组合；若体系超过 200 个原子，可考虑数值原子轨道基组，但必须做基组收敛和 BSSE 校正。

五、参数选择如何影响吸附能：一个可操作的判断路径

把上述选择串起来，可以得到一条实用的判断路径。第一步，确定泛函：若体系含物理吸附，先加色散修正；若含过渡金属，优先 PBE 或 RPBE。第二步，选择赝势：从赝势库中选取与泛函匹配的 PAW 数据集，检查价电子配置是否包含半芯态。第三步，设定基组和截断能：平面波基组下，先做截断能收敛测试，再固定截断能做 k 点收敛。第四步，检查真空层和偶极修正：若吸附质在表面形成垂直偶极，需加入偶极修正，否则功函数变化和吸附能会有系统偏差。

一个常见的误区是只做单参数收敛。实际上，截断能和 k 点之间存在耦合：k 点不足时，截断能收敛曲线可能提前“假收敛”。建议先取一个较密的 k 点（如 7×7×1）做截断能测试，再固定截断能做 k 点测试。对表面吸附，k 点沿表面方向通常取 3×3×1 到 5×5×1，真空方向取 1。若体系较大，可用 Γ 点近似，但必须验证吸附能变化是否在可接受范围。

最后需要强调：工具箱中的每一项都不是孤立的。泛函决定误差的量级，赝势决定价电子的描述精度，基组决定收敛行为，而截断能和 k 点决定数值误差。只有把这条链上的参数都收敛到吸附能精度以下，后续的吸附位比较、覆盖度效应和反应路径计算才有意义。