吸附与催化是当前热门研究领域，常见的计算结果包括催化剂结构，磁性，电子性质等，这些性质对催化剂有着决定性作用。本章华算科技朱老师将正式介绍催化领域的表面与低维结构模型构建，磁性体系处理，关键电子性质计算。
功函数的物理定义与计算公式
功函数的定义公式为：
Φ = E_vacuum – E_Fermi
其中：
E_vacuum：真空区域的静电势平均值（从LOCPOT文件提取）
E_Fermi：费米能级（从OUTCAR文件提取）
计算的关键在于准确获取真空能级和费米能级两个参数。
计算前准备：Slab模型构建
表面结构构建
使用Materials Studio、Vesta或其他建模工具构建表面模型。常见步骤包括：
选择晶面：如Cu(111)、Au(100)、Pt(111)等
构建slab模型：通常选择4-7层原子
设置真空层：真空层厚度≥15 Å，以避免上下表面之间的相互作用
结构优化
固定表面底部原子（通常固定下半部分原子层），模拟体相到表面的过渡
优化晶胞参数和原子位置
使用ISIF=2（仅优化原子位置）或ISIF=3（优化晶胞和原子位置）
INCAR参数设置详解
核心参数配置
在INCAR文件中设置以下关键参数：
参数
推荐值
说明
LVTOT
.TRUE.
输出总静电势到LOCPOT文件
LVHAR
.TRUE.
输出Hartree势和离子势
LDIPOL
.TRUE.
开启偶极修正（消除表面不对称）
IDIPOL
3
偶极校正方向（z轴）
PREC
Accurate
提高计算精度
ENCUT
500-520 eV
平面波截断能
ISMEAR
0
金属体系使用高斯展宽
SIGMA
0.05
展宽宽度
ISYM
0
关闭对称性检测（避免影响电势分布）
LCHARG
.TRUE.
保存电荷密度文件
EDIFF
1E-6
电子步收敛精度
IBRION
-1
静态计算模式
NSW
0
不执行离子弛豫
两种参数配置方案
方案一：加极化修正（推荐）
LVTOT = .TRUE.
LDIPOL = .TRUE.
IDIPOL = 3
DIPOL = 0.0 0.0 0.5 # 设置偶极中心
方案二：不加极化修正
LVTOT = .TRUE.
加极化修正适用于表面有明显极化或吸附体系；不加修正适用于对称性良好的无偶极表面。
计算执行步骤
步骤1：结构优化
运行结构优化计算，获得优化的原子位置：
vasp_std
生成CONTCAR文件，后续用于静态计算。
步骤2：静态自洽计算（SCF）
基于优化后的结构进行静态计算：
1.复制CONTCAR到POSCAR
2.修改INCAR文件为静态计算参数
3.设置KPOINTS网格（与优化时保持一致或更密集）
4.运行VASP计算
vasp_std > vasp.out
生成OUTCAR、LOCPOT、CHGCAR等输出文件。
步骤3：检查收敛性
查看OUTCAR中的最终能量收敛情况
检查电子步（EDIFF）是否满足收敛标准
确保SCF计算完全收敛
后处理与数据提取
从OUTCAR提取费米能级
使用以下命令或文本编辑器从OUTCAR中提取费米能级：
grep “Fermi” OUTCAR
或直接在OUTCAR中查找类似：
Fermi energy : 7.12345678 eV
处理LOCPOT文件获取真空能级
LOCPOT文件包含三维静电势分布信息，需进行平面平均处理：
方法一：使用Fortran程序vtotav.f
原vtotav.f程序需向VASP官方申请，也可使用Fortran 95重写的版本处理LOCPOT文件，生成potential.dat文件。
方法二：使用VASPView或Multiwfn
导入LOCPOT文件，软件会自动处理并输出真空能级。
计算功函数
功函数Φ = E_vacuum – E_Fermi
从LOCPOT中提取真空区域的电势平台值作为E_vacuum，从OUTCAR提取E_Fermi，两者相减即得功函数。
本章要点总结
VASP功函数计算知识点梳理
功函数的物理定义：Slab模型构建、INCAR参数设置详解
计算执行步骤：结构优化、静态自洽计算（SCF）、检查收敛性
后处理与数据提取：从OUTCAR提取费米能级、处理LOCPOT文件、计算功函数
下一步学习建议
下一章将正式引入本次教程的核心—静电势分析。我们将从基本物理概念、后处理方法、具体应用场景方面详细介绍静电势分析。
