VASP 怎么计算对称表面能？分步实操教程

VASP 计算对称表面能需遵循 “晶胞优化 - 模型建立 - 原子弛豫 - 能量计算 - 结果求解” 的标准化流程，核心采用专用拆分公式，需注意弛豫时固定中间层原子、释放上下表面原子，且公式中需除以 2 消除双表面影响，科学指南针结合 Cu (001) 面实例，详解每一步实操要点。

VASP 计算对称表面能的核心公式是什么？

VASP 计算对称表面能的总公式为 γ = (E_relax - E_unrelax)/2A + (E_unrelax - N×E_bulk)/2A，公式分为弛豫项与非弛豫项，除以 2 是因为对称表面模型切割后会产生两个表面，需对能量进行平均拆分。

其中，E_relax 为表面弛豫后的体系总能量，E_unrelax 为表面弛豫前的 slab 模型能量，此为公式的弛豫项，代表表面弛豫所减小的能量；E_unrelax 为弛豫前 slab 能量，E_bulk 为晶体最小晶格中单个原子的能量，N 为 slab 模型中的原子总数，此为公式的非弛豫项，代表切割表面产生的基础能量；A 为模型的表面面积，是能量归一化的关键参数。

如何建立对称表面的 slab 模型？

建立对称表面 slab 模型需先优化晶体晶胞得到最优参数，再根据目标晶面进行切面操作，模型层数建议选择奇数层，且需添加足够厚度的真空层，避免上下表面间的相互作用。

以 Cu (001) 面为例，首先通过 VASP 优化 Cu 的原胞，得到晶格常数、原胞能量等最优参数；再基于最优参数对 Cu 晶胞进行 (001) 面切面，建立七层的 Cu (001) slab 模型，七层为奇数层，便于后续固定中间层原子；最后在模型 z 方向添加 15-20 Å 的真空层，完全消除上下表面的相互作用，保证计算结果的准确性。

图注：展示 Cu (001) 面七层 slab 模型的三维结构，标注表层原子、中间层原子与真空层的位置及厚度，明确模型建立的关键参数。

VASP 计算对称表面能的分步实操流程是什么？

VASP 计算对称表面能的实操流程分为五步，依次是晶体晶胞优化、对称 slab 模型建立、表面原子弛豫、能量参数提取、表面能最终计算，每一步都需严格设置 VASP 参数并验证结果合理性。

第一步，优化晶胞：设置 VASP 弛豫参数，对 Cu 原胞进行结构优化，得到最优 E_bulk 与晶格参数；第二步，建立模型：基于最优参数构建 Cu (001) 面七层 slab 模型，添加真空层；第三步，原子弛豫：固定模型中间三层原子，释放上下两层表面原子，设置弛豫参数进行表面弛豫，得到 E_relax 与 E_unrelax；第四步，参数提取：统计 slab 模型中的原子总数 N，计算模型的表面面积 A；第五步，计算结果：将所有参数代入对称表面能核心公式，求解得到 Cu (001) 面的表面能数值。

对称表面能计算的关键注意事项有哪些？

对称表面能计算的核心注意事项有三点，一是弛豫时必须固定中间层原子，避免整体晶格畸变；二是公式必须除以 2，匹配双表面的模型特征；三是真空层厚度需足够，消除表面间相互作用。

若弛豫时未固定中间层原子，晶胞会发生整体变形，导致计算的弛豫能量与实际表面弛豫能量偏差过大；若公式未除以 2，会将双表面的能量计入单表面，使表面能结果偏大一倍；若真空层厚度不足，上下表面的原子会产生相互作用，破坏表面的对称性，导致所有能量参数失真，科学指南针提供专业的 VASP 计算服务，可规避实操误差，保障计算结果精准。