【摘要】 详解 VASP 计算非对称表面能的公式调整、模型建立与偶极矩校正技巧，明确非对称与对称表面计算的核心差异，科学指南针分享进阶计算要点。

VASP 计算非对称表面能需在对称表面计算基础上进行公式与参数双重调整，公式中需取消弛豫项的除以 2 操作，且因非对称模型会产生偶极矩，需在 VASP 中设置 LDIPOL 与 IDIPOL 参数进行校正，科学指南针详解非对称表面能计算的核心差异与进阶实操方法。

 

非对称表面能与对称表面能的计算有什么核心差异？

非对称表面能与对称表面能计算的核心差异有两点，一是弛豫方式不同，非对称表面仅释放上层原子弛豫、固定底层原子，对称表面则释放上下表面原子；二是公式与参数不同，非对称表面公式弛豫项无需除以 2，且需添加偶极矩校正参数。

对称表面模型的上下表面结构完全一致，弛豫时需同时优化两个表面，因此公式中需除以 2 对两个表面的能量进行平均；而非对称表面模型的上下表面结构存在差异，弛豫时仅优化单个上层表面，底层原子保持固定，因此弛豫项无需除以 2，且非对称的原子分布会导致体系产生偶极矩，若不校正会严重影响能量计算结果。

 

VASP 计算非对称表面能的核心公式如何调整？

VASP 计算非对称表面能的核心公式为 γ = (E_relax - E_unrelax)/A + (E_unrelax - N×E_bulk)/2A，仅将弛豫项的分母 2A 改为 A，非弛豫项保持不变，因弛豫项仅对应单个表面的能量变化，非弛豫项仍由双表面切割产生。

公式调整的核心依据是弛豫的表面数量：非对称表面仅对一个上层表面进行原子弛豫，弛豫项的能量变化仅属于该表面，因此无需除以 2；而非弛豫项是切割晶胞形成 slab 模型时产生的基础能量，无论表面是否对称，切割后都会形成两个表面，因此非弛豫项仍需除以 2，这是公式调整的关键原则，不可混淆。

 

非对称表面模型为什么需要进行偶极矩校正？

非对称表面模型的上下表面原子分布、电子云密度存在显著差异，会导致体系在垂直表面的方向上产生固有偶极矩，偶极矩会引发额外的电场能量，使计算的体系能量出现偏差，因此必须进行偶极矩校正。

非对称表面的原子排列无对称性，上层弛豫后的原子与固定的底层原子在空间分布上形成电荷分离，电子云会向一侧偏移，进而形成偶极矩；该偶极矩产生的电场会与 VASP 计算中的平面波基组相互作用，导致能量计算值偏离实际值，且非对称性越强，偶极矩越大，能量偏差越明显，因此偶极矩校正是非对称表面能计算的必要步骤。

 

VASP 中如何设置参数进行偶极矩校正？

VASP 中对非对称表面模型进行偶极矩校正需在 INCAR 文件中添加两个核心参数，分别是 LDIPOL=.TRUE. 与 IDIPOL=3，LDIPOL 用于开启偶极矩校正功能，IDIPOL=3 代表在 z 方向进行校正，适配 slab 模型的真空层方向。

具体参数设置要点：一是 LDIPOL 参数设为.TRUE.，激活 VASP 的偶极矩校正模块，若设为.FALSE. 则关闭校正功能；二是 IDIPOL 参数设为 3，因非对称 slab 模型的真空层通常设置在 z 方向，偶极矩也沿 z 方向产生，该参数指定校正的空间方向；三是参数设置后需重新进行表面弛豫计算，确保校正后的能量参数为精准值，不可直接使用未校正的能量数据。

 

非对称表面能计算的完整实操流程是什么？

非对称表面能计算的完整流程为晶胞优化→非对称 slab 模型建立→偶极矩参数设置→单表面原子弛豫→能量参数提取→表面能公式计算，核心是模型非对称构建、参数校正与公式弛豫项调整。

第一步，优化晶体原胞得到 E_bulk 与最优晶格参数；第二步，构建非对称的目标晶面 slab 模型，添加真空层，保证上下表面结构差异；第三步，在 INCAR 文件中设置 LDIPOL=.TRUE.、IDIPOL=3，开启偶极矩校正；第四步，固定模型底层原子，释放上层原子进行表面弛豫，得到 E_relax 与 E_unrelax；第五步，统计原子总数 N，计算表面面积 A；第六步，将参数代入调整后的非对称表面能公式，求解得到最终表面能数值。科学指南针拥有专业的材料模拟计算团队，可提供对称与非对称表面能的全流程 VASP 计算服务，解决进阶计算中的各类技术问题。