VASP怎么计算：从输入文件到结果分析的完整入门指南

VASP怎么计算——这个问题几乎是每一个计算材料学初学者面对的第一个门槛。VASP（Vienna Ab-initio Simulation Package）是密度泛函理论框架下最广泛使用的平面波赝势代码之一，但其输入参数众多、报错信息晦涩，新手往往在第一轮计算中就陷入困境。这个项目以硅体的能带和态密度计算为载体，完整记录了从零开始到获得可用数据的全过程，每一步的关键决策和常见陷阱都值得梳理。

VASP的运行依赖四个输入文件：POSCAR、POTCAR、KPOINTS和INCAR。这四个文件的逻辑关系构成了VASP计算的基本框架。

POSCAR定义晶体结构。项目以硅体为例，采用标准菱面胞描述：

“`

Si

5.43

0.0 0.5 0.5

0.5 0.0 0.5

0.5 0.5 0.0

2

Direct

0.00 0.00 0.00

0.25 0.25 0.25

“`

第一行是注释，第二行是晶格常数（Å），随后三行是晶格矢量，接着是原子种类数和坐标。这个格式看起来简单，但项目组初期曾因坐标行前误加空格导致结构读取错误——VASP对格式敏感度极高。

POTCAR是赝势文件，从VASP官方库中拼接获得。硅体选用PBE交换关联泛函下的Si赝势，POTCAR通过命令行拼接生成。赝势版本需与POSCAR中的元素顺序严格对应，顺序错位是最常见的低级错误之一。

KPOINTS文件定义布里渊区采样网格。对于硅体能带计算，项目需要两类KPOINTS：自洽计算用均匀网格，能带计算用高对称路径。

均匀网格采用Gamma中心化方案，从4×4×4开始测试：

– 4×4×4：总能 -10.521 eV

– 8×8×8：总能 -10.546 eV

– 12×12×12：总能 -10.548 eV

– 16×16×16：总能 -10.548 eV

12×12×12以上总能变化进入2 meV平台区。差距不会说谎——稀疏K点在硅体这种小带隙半导体中影响尤为显著，低估采样密度会直接歪曲带隙大小。

ENCUT（截断能）的收敛测试同样不可省略。项目从300 eV扫描至550 eV，400 eV时总能收敛至1 meV以内。对于硅体，ENCUT = 400 eV是PBE赝势ENMAX值（约245 eV）的1.3倍——经验法则建议ENCUT取ENMAX的1.3倍以上，但这个倍率对不同体系并不通用，必须逐体系验证。

INCAR是VASP计算的”控制面板”，参数众多但核心项并不复杂。项目自洽计算的INCAR关键设置如下：

“`

ENCUT = 400

ISMEAR = 0

SIGMA = 0.05

EDIFF = 1E-6

ICHARG = 2

LCHARG = .TRUE.

LWAVE = .TRUE.

“`

ISMEAR的选择有明确的物理依据。半导体和绝缘体用ISMEAR = 0（Gaussian展宽），SIGMA = 0.05 eV；金属体系可用ISMEAR = 1或2（Methfessel-Paxton），SIGMA = 0.2 eV。项目初期误将ISMEAR设为-5（四面体法）做结构弛豫，VASP直接报错——四面体法不支持力的计算，只适用于最终自洽计算。

EDIFF控制电子步收敛。1E-6 eV是能带和态密度计算的推荐值；结构弛豫可用1E-5以提高速度，但最终自洽计算应收紧至1E-6。项目在态密度计算中曾用1E-4，结果费米面附近的能级出现非物理抖动——电子步未充分收敛时，占据数的微小波动会被DOS计算放大。

完整流程分三步走。第一步结构弛豫：ISIF = 3（优化晶胞和原子位置），EDIFFG = -0.01 eV/Å，NSW = 100。弛豫完成后检查CONTCAR中的应力和力是否满足收敛判据。

第二步自洽计算：用弛豫后的CONTCAR替换POSCAR，ICHARG = 2从头做电子自洽，输出CHGCAR。这一步K点网格应与弛豫一致或更密。

第三步能带计算：沿高对称路径（硅体的Γ-X-W-K-Γ-L-U-W-L）设置KPOINTS，ICHARG = 11读取上一步的CHGCAR，LWAVE = .FALSE.关闭波函数输出。项目采用32个K点每段，总路径约200个K点。

能带数据从EIGENVAL文件提取。硅体的计算带隙约0.6 eV，而实验值1.12 eV——这是PBE泛函的系统性低估，并非计算错误。根据Nature Materials的综述讨论，DFT-PBE对半导体带隙的低估源于交换关联势的自相互作用误差，需用GW近似或杂化泛函修正。项目用HSE06重新计算，带隙提升至1.15 eV，与实验高度吻合。

态密度从DOSCAR提取。项目在费米能级附近加密了NEDOS = 2000，确保价带顶和导带底的DOS峰形清晰。SIGMA值在DOS计算中应减小至0.02 eV，过大的展宽会抹平精细结构。

回过头看，VASP怎么计算这个问题的答案并不在于记住一套固定参数模板，而在于深入理解每个参数背后的物理含义，并针对具体体系逐一验证收敛性。项目在硅体上建立的这套计算流程后来成功迁移到了GaAs和ZnO等不同体系，但每次迁移都重新做了ENCUT和K点收敛测试——这是VASP计算结果可信度的根本保障。模板可以复用，但收敛验证永远不可省略。