在上一章《第二十四VASP计算实战：吸附能计算流程讲解！ | 2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP表面与催化计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。
在第一性原理计算中，吉布斯自由能G 通常通过以下公式计算：
：VASP计算得到的基态电子能量（通常为SCF收敛后OUTCAR中的
）。
ZPE (Zero Point Energy)：零点振动能，是分子在0K时仍存在的能量。
：热焓修正，包含了从0K到设定温度（如298.15K）过程中的振动、转动、平动贡献（Pv项对于固体通常可以忽略）。
TS：熵项，包含振动、转动、平动及电子熵。
注意：对于固体（如吸附在表面的分子、晶体材料），通常只考虑振动贡献；而对于气体分子，则需要考虑平动、转动和振动三方面的熵与焓。
整个流程主要分为结构优化、力常数矩阵计算（频率计算）和热力学后处理三大步骤。
这是所有计算的基础，确保计算的是系统的基态结构。
输入文件设置：
INCAR：PREC = Accurate，ENCUT设置足够高（如500 eV），EDIFF = 1E-5（能量收敛），EDIFFG = -0.02（力收敛），ISIF = 3（同时优化离子和晶格）。
KPOINTS：对于体相材料，使用密集的k点网格（如
）；对于表面或分子，使用Gamma点或少量k点。
POTCAR：选择合适的赝势（如PAW-PBE）。
POSCAR：输入初始原子坐标。
运行：使用标准的vasp_std或vasp_gam执行。
目标：得到优化后的能量
和对应的基态结构文件（CONTCAR）。
为了得到更精确的电子能量，通常在优化结构的基础上进行一次更高精度的静态计算。
设置：
NSW = 0（关闭离子步进，固定结构）。
ISMEAR = -5（Gaussian Smearing，适用于半导体/绝缘体），或 ISMEAR = 0（适用于金属）。
LCHARG = .TRUE.，LWAVE = .TRUE.（保存电荷密度和波函数以备后用）。
目的：获得最终的基态能量
，用于后续的热力学计算。
这是吉布斯自由能计算中最关键的一步，决定了ZPE、振动熵和热焓的精确度。
计算方法：
VASP内部的Finite Displacement (FD) 方法：最常用。通过在每个原子上施加微小位移（通常为0.01 Å），计算受扰动后的力。
关键设置：
INCAR：IBRION = 5（Finite Displacement），ISIF = 2（仅优化离子），POTIM = 0.015（位移步长），NFREE = 2（双边位移），EDIFF = 1E-8（高精度）。
KPOINTS：通常使用Gamma点即可。
这是计算吉布斯自由能的核心环节，需要将振动频率转化为热力学量。
计算公式（基于统计力学）：
ZPE：
振动自由能：
振动熵：
对于气体分子（如
等），需要考虑理想气体的平动和转动贡献。
经验公式：
平动：
转动：
焓修正：
（理想气体）
频率计算是最耗时的环节，因为需要进行3N 次（N 为原子数）的单点计算。
加速策略：并行计算：使用 vasp_gibbs 脚本时，可通过参数 -n 指定CPU核数，或通过 -m 指定MPI执行文件。
对称性分析：自动识别分子对称性（如
），剔除冗余位移，大幅减少计算量。
理论基础：介绍了吉布斯自由能计算公式中各项物理意义
计算流程：介绍结构优化，静态计算，频率计算等
关键参数：介绍了气象物种处理与计算加速
下一章将正式引入本次教程的核心—迁移路径构建。我们将从理论基础、迁移路径构建、输入文件准备方面详细介绍VASP 电池材料案例，敬请期待！
