VASP计算功函数：板模型设计与偶极校正的完整方案

VASP计算功函数的公式只有一步：Φ = V_vacuum – E_Fermi。V_vacuum是真空中静电势的平台值，E_Fermi是费米能级。项目组在计算Pt(111)功函数时，不加偶极校正得到约5.3 eV，加上校正后约5.6 eV——实验值约5.7 eV。差了约0.4 eV，根源是吸附或不对称slab引入的偶极层在真空中产生了额外电场，扭曲了静电势参考点。

功函数的物理本质

功函数是把一个电子从费米能级移到真空静止状态所需的能量。它由两部分组成：体相的费米能级（材料的内禀性质）和表面偶极层（与表面终止方式有关）。同一材料的不同晶面功函数不同——Pt(111)为5.7 eV，Pt(100)为5.6 eV，Pt(110)为5.4 eV——差异来自表面偶极层的强弱。

这个表面依赖性意味着功函数计算必须正确建模表面。slab模型的两个关键参数——板厚度和真空层厚度——直接影响结果的可靠性。

板厚度收敛

板太薄会导致上下表面电子态耦合，功函数不收敛。项目组测试了Pt(111)从3层到10层slab的功函数变化：

不加偶极校正时，5层到10层的变化小于0.05 eV——看似收敛了，但收敛值比实验值低约0.4 eV。加偶极校正后，6层即收敛到约5.6 eV，与实验值5.7 eV偏差约0.1 eV——在DFT精度范围内。

偶极校正如此重要的原因是：即使是”对称”slab（上下表面都是(111)），弛豫后上下表面的原子位移不对称（因为底部固定、顶部弛豫），产生残余偶极矩。这个偶极矩在真空中形成一个电场，使得slab两侧的真空势不相等——V_vacuum_top ≠ V_vacuum_bottom。

偶极校正的具体操作

VASP中偶极校正通过两个参数实现：

LDIPOL = .TRUE.   # 开启偶极校正
IDIPOL = 3         # 校正方向（3=z方向，垂直于表面）

启用LDIPOL后，VASP在每次电子自洽迭代中：

1. 计算slab的偶极矩
2. 在真空区域施加一个补偿电荷平板
3. 重新计算静电势

这样真空中的电场被抵消，V_vacuum恢复到无偶极时的参考值。

项目组踩过一个坑：同时设置了LDIPOL=.TRUE.和DIPOL参数（DIPOL指定偶极矩的中心位置）。DIPOL设置不当会导致校正过冲，功函数反而更差。正确做法是：对称slab不设DIPOL（VASP自动计算），非对称slab（如吸附体系）设置DIPOL在slab中心。

真空层厚度的选择

真空层太薄会导致周期镜像间的静电相互作用，功函数不收敛。项目组的收敛测试：

15 Å已基本收敛，但这是对称slab的情况。对于吸附体系（表面有吸附分子，偶极矩更大），需要20-25 Å真空层。项目组的标准配置是：对称slab 15 Å，吸附slab 20 Å。

功函数计算的后处理

VASP计算完成后，需要从LOCPOT文件提取静电势沿z方向的平均值。项目组用以下Python脚本处理：

from pymatgen.io.vasp import Locpot
import numpy as np

locpot = Locpot.from_file("LOCPOT")
# 获取沿z方向的平均静电势
potential_z = locpot.get_average_along_axis(2)
# 真空势 = 真空区域中的势能平台值
V_vacuum = np.max(potential_z)
# 费米能级从OUTCAR读取
E_fermi = -0.37  # 示例值
work_function = V_vacuum - E_fermi

关键判断：V_vacuum应该是静电势在真空区域中的平台值。如果真空层不够厚，看不到明显的平台，说明需要增加真空层。项目组画势能vs z位置的曲线来直观检查——如果曲线在真空区域是水平的，说明真空层足够；如果还在变化，说明太薄。

吸附体系中的功函数变化

功函数计算的一个重要应用是评估吸附对表面电子性质的影响。CO吸附在Pt(111)上后，功函数降低约0.7-0.8 eV——这是因为CO的5σ轨道向Pt转移电子，在表面形成指向真空的偶极层，降低了电子逸出的能量壁垒。

Fermi能级变化和真空势变化对功函数的贡献需要分开分析：

• ΔE_Fermi 约+0.1 eV（费米能级上移）
• ΔV_vacuum 约-0.6 eV（真空势下降）

功函数变化主要由真空势下降驱动——即表面偶极层的变化。这个分解对于理解吸附引起的电子结构重排很有价值。

反思：功函数的计算精度上限

DFT-PBE对金属功函数的计算精度通常在0.1-0.3 eV。要进一步提升精度需要考虑：

• 杂化泛函（HSE06）：可改善0.05-0.1 eV，但计算量增加10-50倍
• G0W0准粒子修正：可改善0.05-0.15 eV，但只适用于小体系
• 全电子计算（WIEN2k或ELK）：避免赝势引入的误差，改善约0.05 eV

对于大多数应用（如催化剂筛选、异质结能带对齐），PBE+偶极校正的0.1-0.3 eV精度已经足够。但在精确的Schottky势垒计算中，这个偏差可能决定界面是欧姆接触还是肖特基接触——需要更高层级方法。

板模型设计是功函数计算的基础，偶极校正是精度的保障。更多VASP表面计算的实战经验，可以参考VASP/第一性原理栏目，或返回科研学术网首页。