第二十章VASP计算实战：静态自洽计算步骤讲解！|2026新版VASP基础教程

一、引言

在上一章《第十九章：VASP二维材料结构优化步骤详解！| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP示例任务。VASP计算一般从结构优化开始，然后对优化后的结构进行性质计算。本章将介绍VASP结构优化计算实例，让大家更好的了解计算流程，具体包括结构优化、半导体性质案例、表面与催化案例。

二、静态自洽计算的物理意义与定位

在 VASP 计算中，静态自洽计算通常指的是 SCF（Self-Consistent Field） 步骤，它是基于已知原子位置（通常是结构优化后的结构）进行的电子结构求解。

位置固定

与结构优化（IBRION ≠ -1）不同，静态自洽计算中原子位置是固定不动的（IBRION = -1）。

电子迭代

仅进行电子迭代（NELM 次），直至电子密度收敛（EDIFF 标准），此时体系达到基态能量极小值。

结果产出

主要输出基态电子密度（CHGCAR）、波函数（WAVECAR）以及能量信息（OUTCAR、OSZICAR）。这些文件是后续性质计算（如能带、态密度）的基础。

三、前期准备：文件与目录结构

在开始静态自洽计算之前，需要准备好 VASP 所需的输入文件。假设你已经完成了结构优化（Relaxation），则静态自洽计算的准备工作如下：

复制优化后的结构文件

结构优化完成后，VASP 会生成 CONTCAR 文件（包含优化后的原子坐标）。静态自洽计算需要使用该文件作为新的 POSCAR。

cp CONTCAR POSCAR 

准备势能文件（POTCAR）

确保你已经为体系中所有原子准备好了正确的势能文件（POTCAR）。这通常需要使用 pseudopotential 脚本或手动拼接。

cat POTCAR_AL POTCAR_O > POTCAR   # 例如对于 Al₂O₃ 

K 点网格（KPOINTS）

静态自洽计算对 K 点网格要求较高，需要足够密集以保证能量收敛。通常使用 Monkhorst-Pack 网格。

Automatic mesh 

0

Monkhorst-Pack

6 6 6 

0 0 0 

四、关键输入文件（INCAR）参数设置

INCAR 是 VASP 的核心控制文件，其中包含了计算类型、收敛标准以及输出选项。以下是静态自洽计算的典型参数设置：

SYSTEM = Al₂O₃ static SCF calculation 

PREC = Accurate      # 计算精度 

ENCUT = 520          # 能量截断，建议比POTCAR推荐值高 30% 左右 

EDIFF = 1E-6         # 电子收敛标准，通常 1E-5 – 1E-6 

NELM = 100           # 最大电子迭代步数 

ISMEAR = 0           # 高斯展宽（针对金属可设 -5） 

SIGMA = 0.05         # 展宽宽度 

IBRION = -1          # 禁止离子运动（固定原子位置） 

NSW = 0              # 步数设为0（不进行离子迭代） 

LWAVE = .TRUE.       # 输出 WAVECAR（波函数文件） 

LCHARG = .TRUE.      # 输出 CHGCAR（电荷密度文件） 

LREAL = Auto         # 实空间投影方式 

参数详解与注意事项

IBRION = -1 与NSW = 0

这两个参数共同作用，确保计算过程仅包含电子迭代，原子位置保持不变。如果缺少这两个设置，VASP 可能会误认为是结构优化或分子动力学。

LWAVE 与LCHARG

通常需要打开这两个选项，以便输出 WAVECAR 和 CHGCAR。后续的非自洽计算（如能带计算）需要读取 WAVECAR（波函数）和 CHGCAR（电荷密度）。

ISMEAR 与SIGMA

对于绝缘体和半导体，推荐使用 ISMEAR = 0（高斯展宽）或 ISMEAR = -5（tetrahedron 方法）。对于金属体系，常用 ISMEAR = 1 或 -1，并设置较小的 SIGMA（如 0.05）。

PREC 与ENCUT

PREC = Accurate 是推荐设置。ENCUT 的选择应基于所使用的 POTCAR 中的最大 ENMAX，通常取该值的 1.3 倍。

五、计算执行与监控

运行 VASP

在准备好 INCAR、POSCAR、POTCAR、KPOINTS 四个文件后，即可在计算节点上提交 VASP 任务。

mpirun -np 32 vasp_std > vasp.out 

监控收敛情况

静态自洽计算的核心是电子收敛。你需要关注以下文件：

OSZICAR：显示每一步的能量变化（E0）和收敛标准（EDIFF）。当 E0 在连续几步内变化非常小（如 1e-6 eV）时，即表示收敛。

OUTCAR：详细记录了每一步的收敛信息。搜索关键词 reached required accuracy 或 converged。

六、计算结束后的结果文件解析

静态自洽计算结束后，VASP 会生成一系列输出文件，其中最关键的三个文件是：

CHGCAR（电荷密度文件）

记录了体系在基态下的电子密度分布。

用途：后续计算（如 Bader 电荷分析、静电势计算）需要读取此文件。

读取：使用 chgsum.pl 或 bader 工具进行后处理。

WAVECAR（波函数文件）

包含了体系的 Kohn-Sham 轨道信息。

用途：非自洽计算（NSCF）和能带计算（Band Structure）必须读取该文件，以便在新的 K 点路径上重新计算能量。

注意：WAVECAR 文件非常大（可能数 GB），复制时需注意磁盘空间。

OUTCAR（输出日志）

记录了详细的计算过程、能量、力、应力张量等信息。

用途：分析体系的总能量（free energy）、费米能级（E-fermi）以及自洽误差。

读取：可以使用脚本（如 grep）提取关键参数。

七、本章要点总结

静态自洽计算知识点梳理

原理与意义：介绍了位置固定与电子迭代

文件与目录结构：介绍了结构文件、K点网络文件、赝势文件

输入输出文件设置：介绍了INCAR、OUTCAR等文件意义

下一步学习建议

下一章将正式引入本次教程的核心—能带结构计算与绘图。我们将从材料结构与输入文件准备、分步进行计算、数据提取与绘图方面详细介绍VASP半导体性质计算案例，敬请期待！