在上一章《第二章VASP计算实战：MgO结构优化步骤详解！| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP结构优化计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。
在对二维材料进行结构优化时，最容易忽视的细节是真空层（Vacuum Layer）和C 轴（垂直方向）的处理。
由于 VASP 使用三维周期性边界条件（PBC），二维材料必须在垂直方向（通常是 C 轴）添加真空层，以避免不同层之间的相互作用。
在优化过程中，如果 C 轴不固定，真空层的厚度可能会被误认为是应力大的区域，导致盒子在 C 方向“坍塌”，甚至把真空层压缩到 0 Å，导致层与层之间相互作用，结果毫无物理意义。
最常见的做法是通过 Materials Project 或 The Materials Project 等数据库下载材料的 .cif 文件。
使用 VESTA、Materials Studio 或 ASE 等工具打开 .cif 文件，将 C 轴拉伸到 15 Å – 20 Å，确保层与层之间足够分离。
确认原子顺序（例如：Mo、S、S），然后保存为 VASP 可读的 POSCAR 文件。
在 VASP 计算中，四个文件（POSCAR、POTCAR、KPOINTS、INCAR）决定了计算的全部参数。
确保使用与元素对应的赝势，并放在同一目录下。
对于二维材料，KPOINTS 的设置通常是：
X-Y 平面：密度较高，如 12 12 1 或 15 15 1（根据体系大小调整）。
Z 方向：固定为 1，因为有真空层，Z 方向的周期性并不重要。
这是最关键的控制文件，以下参数必须特别关注：
# 计算类型
IBRION = 2 # 使用共轭梯度法进行离子弛豫
ISIF = 2 # 只优化离子位置，固定晶格参数
NSW = 100 # 最大离子步数
EDIFF = 1E-5 # 能量收敛标准
# 电子迭代
ENCUT = 500 # 截断能，确保足够大以适应真空层
EDIFFG = -0.02 # 力收敛标准，单位为 eV/Å
# 体系特性
ISYM = 0 # 关闭对称性（对二维材料常见）
LCHARG = .FALSE. # 不生成 CHGCAR
LWAVE = .FALSE. # 不生成 WAVECAR
特别提示： ISIF=2 是关键，它指示 VASP 固定格点（包括 C 轴长度），只优化原子位置。
在开始正式优化之前，建议先进行一次快速预优化：
使用较低的 ENCUT（如 350 eV）。
使用较粗的 KPOINTS（如 6 6 1）。
目的是快速消除大的力（Force）或应力，节省时间。
完成预优化后，使用更高精度的参数（如 ENCUT=500，KPOINTS=12 12 1）进行正式优化。
能量变化（ECHANGE）：通常要求低于 1E-5 eV。
最大力（FORCE）：要求低于 0.02 eV/Å。
检查 CONTCAR：计算结束后，VASP 会生成 CONTCAR 文件，包含优化后的原子坐标。
构建POSCAR：介绍了获取结构、添加真空层、保存文件
输入文件设置：介绍了INCAR、KPOINTS、POTCAR文件
计算流程：介绍了预优化、正式优化、监控收敛
下一章将正式引入本次教程的核心—静态自洽计算。我们将从原理与意义、文件与目录结构、输入输出文件设置方面详细介绍VASP半导体性质计算案例，敬请期待！
