在上一章《第七章：如何选择VASP执行文件？| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了Linux系统基础—脚本编写。VASP需要4个输入文件INCAR、KPOINTS、POSCAR、POTCAR来控制计算任务。本章将正式介绍VASP执行流程，从输入与输出文件角度来介绍VASP运行方法。
VASP计算的核心是四个输入文件，它们共同定义了计算的“硬件配置”（赝势）、“软件设置”（参数）和“待测物体”（结构）。
描述原子体系的几何结构，是计算的“原料”。
包含每种原子的赝势信息，决定了计算的精度和方法。
控制计算的各种开关，类似于操作系统的“系统设置”。
定义电子波矢的采样网格，类似于“取样频率”。
这四个文件必须放在同一工作目录下，且文件名必须全部大写，无扩展名，文件内容必须是纯文本。
POSCAR 是 VASP计算中最直观的文件，它描述了你要计算的物质是什么。
编写顺序与格式
1. 标题行（任意字符串，用于描述体系，如 Si or SiO2）
2. 缩放因子（通常设置为 1.0，除非你想整体放大/缩小体系）
3. 晶格向量（3 行，每行 3 个数，单位为 Å，定义了晶胞的形状和大小）
4. 原子种类（如 Si O）
5. 每种原子的个数（对应上面的顺序，如 2 4）
6. 坐标系类型（选择 Direct（分数坐标）或 Cartesian（笛卡尔坐标））
7. 原子坐标（每个原子对应一行坐标，若为 Direct 则为相对晶格向量的比例，若为 Cartesian 则为 Å 坐标）
编写技巧
坐标系选择：Direct 坐标是最常用的，因为它与晶格参数直接关联，便于晶胞形变时的坐标转换。
结构来源：POSCAR 可以手动输入（适合小体系），也可以通过 Materials Studio 导出为 CIF/XSD 格式，再使用 VESTA、VESTA2 或 VASPKIT 等工具转换生成。
高对称性：如果体系具有高对称性（如金刚石、石墨），请务必使用晶格的标准格（primitive cell），以减少计算量。
POTCAR 包含了原子核电荷与价电子之间的有效相互作用势，它是VASP计算最重要的“核心”文件。
获取方式
POTCAR 不能随意编写，必须使用 VASP 官方提供的 potcar 命令或下载的赝势库（POTCAR files）进行合并。
合并方法
假设你的体系中有 Si 和 O 两种元素，需要先下载对应元素的赝势文件（如 Si 和 O 文件夹），然后在工作目录下运行：
cat /path_to_potentials/Si/POTCAR /path_to_potentials/O/POTCAR > POTCAR
注意事项
赝势版本：VASP 提供了不同的赝势版本（如 PBE、LDA、PAW_GGA 等），请确保你选择的赝势与计算的交换相关函数（XC functional）一致。
赝势检查：合并完成后，建议使用 grep 命令查看 POTCAR 中的关键参数（如 ENMAX），确保所有元素的能量截断要求（ENCUT）相近。
INCAR 是 VASP计算中最灵活、也是最复杂的文件，它决定了计算是要“优化结构”还是“计算能带”，以及计算的精度。
核心参数解释
以下参数是最基础且必不可少的设置：
SYSTEM = My_Si_Diamond # 计算体系的名称，纯粹用于区分计算
ISTART = 0 # 计算开始方式：0表示新计算，1表示接着上一次计算
ICHARG = 2 # 电荷密度的初始化方式：2表示自洽初始化
ENCUT = 520 # 平面波基组截断能（eV），通常取最大 ENMAX 的 1.3~1.5 倍
PREC = Accurate # 计算精度：Normal/Accurate/Low
EDIFF = 1E-5 # 电子自洽收敛标准（总能量变化阈值）
ISMEAR = 0 # 积分方式：0用于绝缘体，-5用于金属，1用于金属
SIGMA = 0.05 # 高斯展宽（eV），通常与 ISMEAR 搭配使用
IBRION = 2 # 离子弛豫算法：2表示共轭梯度法（CG）
NSW = 100 # 离子最大步数
ISIF = 3 # 是否弛豫晶格：3表示弛豫原子和晶格
KPOINTS 决定了电子波矢的采样密度，直接影响计算结果的收敛性。
两种主要模式
1、自动模式（Automatic Generation）
这是最常用的方式，适用于结构优化、能量计算等。
Automatic mesh
0
Gamma
11 11 11
0 0 0
第一行是描述行。
第二行设置为 0（表示自动生成）。
第三行选择网格类型：Gamma（Gamma centered）或 Monkhorst-Pack。
第四行是网格密度（如 11x11x11），密度越高计算越慢但越精确。
第五行是位移（通常设为 0）。
2、显式模式（Explicit List）
适用于能带计算或自定义 k 点路径。
K-points along high symmetry lines
8
Line-mode
Reciprocal
0.0 0.0 0.0 ! Gamma
0.5 0.0 0.0 ! X
0.5 0.5 0.0 ! M
0.0 0.0 0.0 ! Gamma
第三行设为 Line-mode。
第四行选择坐标系：Reciprocal（倒空间）。
之后的每两行定义了一条 k 轨迹。
收敛性测试
在正式计算前，务必进行 KPOINTS 的收敛性测试。一般来说，体心立方（BCC）或面心立方（FCC）体系的 k 网格密度可以设得较低（如 6x6x6），而低对称性或大超胞体系则需要更稀疏的网格（如 2x2x2）。
INCAR：计算内容，收敛标准，计算方法
POSCAR：计算模型
KPOINTS：布里渊区数值积分采用网格
POTCAR：元素赝势
下一章将正式引入本次教程的核心—执行文件vasp_std/vasp_gam。我们将从vasp_std功能与适用范围、vasp_std功能与适用范围，vasp_std功能与适用范围方面详细介绍执行文件vasp_std/vasp_gam，以及他们在VASP计算中的应用，敬请期待！
