3. VASP输入文件INCAR、KPOINTS、POTCAR
VASP的一次计算能否顺利收敛、结果是否可信，很大程度上取决于三个输入文件的配合：INCAR决定“怎么算”，KPOINTS决定“在倒空间哪些点采样”，POTCAR决定“用什么赝势描述离子实与价电子”。POSCAR提供结构，而这三个文件共同定义数值精度与物理近似。理解它们之间的耦合关系，比孤立记忆参数更重要。
INCAR：控制算法与收敛判据
INCAR由一系列 TAG = VALUE 组成，不区分大小写，# 之后为注释。核心参数可分为四类：
- 电子步控制：
ENCUT（平面波截断能，单位eV）、PREC（精度，常用Normal/Accurate）、EDIFF（电子步能量收敛判据，典型 1E-5～1E-6 eV）、NELM（最大电子步数，默认60）、ALGO（算法，Normal/Fast/All/VeryFast）、ISMEAR与SIGMA（展宽方式与宽度）。 - 离子步控制：
IBRION（离子弛豫算法，-1为不动、0为分子动力学、2为共轭梯度）、NSW（最大离子步数）、EDIFFG（离子步收敛判据，力小于该值即停）、ISIF（是否优化晶胞，2为优化离子、3为同时优化离子与晶胞）。 - 自旋与磁性：
ISPIN = 2开启自旋极化，配合MAGMOM给出初始磁矩。 - 输出与任务类型：
LWAVE、LCHARG控制是否写WAVECAR与CHGCAR；ICHARG决定电荷密度来源（0为从初始原子叠加、1为读CHGCAR、11为读CHGCAR且不更新）。
其中 ENCUT 是最需要谨慎的参数。它决定平面波基组的完备程度：截断能过低，总能量与受力都不可靠；过高则计算量急剧上升。工程上通常取 POTCAR 中 ENMAX 最大值的 1.3 倍左右，或直接做一组截断能扫描，观察总能量差小于 1 meV/atom 时即可认为收敛。
KPOINTS：倒空间采样密度
KPOINTS最常用的是自动网格模式（第一行写Auto或0，第二行给网格数，第三行写Gamma或Monkhorst-Pack）。网格密度用 N1×N2×N3 表示，其物理含义是沿三个倒格矢方向的采样点数。对金属体系，K点需足够密才能准确描述费米面；对绝缘体或大带隙体系，较稀的网格即可收敛。
一个实用的经验判据是：K点间距（即倒空间相邻采样点距离）控制在 0.02～0.03 Å⁻¹ 之间。若晶格常数为 a，则 N ≈ a × 0.03 的倒数。例如 a = 4 Å 时，N 约取 8～10。对于二维材料或纳米带，真空方向只需 1 个K点。此外，ISMEAR 与K点密度需匹配：金属用 ISMEAR = 1 或 2，半导体/绝缘体用 ISMEAR = 0 配合较小 SIGMA（如 0.05）。
POTCAR：赝势选择与拼接顺序
POTCAR包含每种元素的赝势数据，其元素顺序必须与POSCAR中元素顺序严格一致，否则会得到完全错误的能量。选择赝势时优先考虑：
- PAW_PBE：通用性最好，适合大多数体系，推荐作为默认选择。
- PAW_LDA：仅在与旧文献对比或特定表面问题时使用。
- PAW_PBE 半芯态版本（如
Ti_pv、Fe_pv）：当体系涉及d电子强关联或需要更高精度时选用，但会显著增加计算量。
拼接POTCAR时可用 cat 命令按POSCAR元素顺序合并，例如 cat POTCAR_A POTCAR_B > POTCAR。务必检查 grep ENMAX POTCAR 确认最大截断能，并据此设置 ENCUT。
自洽与非自洽计算流程
标准流程分为两步：自洽计算（SCF）与非自洽计算（NSCF）。SCF的目标是获得收敛的电荷密度与总能量，此时 ICHARG = 2（默认，从原子叠加开始）或 ICHARG = 0，K点网格较密，EDIFF 取 1E-6。SCF完成后会写出 CHGCAR 与 WAVECAR。
NSCF则固定电荷密度，读取 CHGCAR（ICHARG = 11），用于计算态密度、能带或光学性质。此时K点网格通常更密（尤其做DOS时），且 ISMEAR = -5（四面体方法）配合 LORBIT = 11 输出投影态密度。NSCF不更新电荷密度，因此单步即可完成，但需要保证 NELM 足够大以完成对角化。
常见参数组合与排错要点
下表给出几种典型场景的推荐组合，可作为起步模板，但需根据具体体系微调。
| 场景 | ENCUT | KPOINTS | ISMEAR / SIGMA | IBRION / ISIF |
|---|---|---|---|---|
| 金属体相弛豫 | 1.3×ENMAX | 12×12×12 | 1 / 0.2 | 2 / 3 |
| 半导体能带 | 1.3×ENMAX | 8×8×8 (SCF) → 沿线 (NSCF) | 0 / 0.05 | -1 / 2 |
| 表面吸附 | 1.3×ENMAX | 6×6×1 | 0 / 0.05 | 2 / 2 |
| 分子/团簇 | 1.3×ENMAX | 1×1×1 (Gamma) | 0 / 0.01 | 2 / 2 |
ENCUT低于赝势ENMAX时VASP会警告并可能给出错误结果；金属体系使用ISMEAR=0且SIGMA过小，会导致收敛困难或能量振荡。
ENCUT（如300、400、500、600 eV），观察总能量变化；固定ENCUT、改变K点网格（如6×6×6、8×8×8、10×10×10），观察能量与力。当相邻两档能量差小于1 meV/atom时，可认为参数已收敛。
最后强调一点：INCAR、KPOINTS、POTCAR三者不是独立参数，而是相互制约的。截断能决定基组大小，K点密度决定采样精度，赝势决定价电子描述层次。任何一项设置不当，都会在总能量、力、态密度或能带中留下系统性偏差。建议每次新体系都先做小规模收敛测试，再进入正式计算，这样比事后反复调试更省机时。