VASP理论计算的核心参数设置与收敛判据

从一个Cu(111)表面计算说起

本项目最近接到一组铜基催化剂的VASP理论计算任务，需要在Cu(111)表面吸附CO分子并分析电荷转移。这个看似常规的任务，实际上对INCAR参数的选取极为敏感——ENCUT偏低会导致力的计算不准，SIGMA过大则让能量展宽偏离物理真实值。本项目中走了不少弯路才摸清各参数的合理边界，这里把经验整理出来，供同行参考。

VASP理论计算的核心流程可以概括为：结构优化→静态自洽计算→性质计算（能带/态密度/电荷密度）。每一步的参数继承关系不是简单复制，而是需要根据前一步的收敛情况做针对性调整。比如结构优化用的EDIFFG=-0.02（力收敛判据0.02 eV/Å），到了静态计算阶段应该收紧到EDIFF=1E-6（能量收敛判据），否则后续的态密度可能出现非物理振荡。

关于[科研计算服务](https://www.keyanxueshu.com/)的整体方案设计，本项目在另一篇中已有详细讨论。本篇聚焦参数层面。

ENCUT截断能：不是越高越好

ENCUT决定了平面波基组的截断能量，直接影响力、应力和能量的精度。VASP手册建议取POTCAR中ENMAX的最大值，但在实际计算中这个值往往不够。

本项目的踩坑经历：计算Cu(111)-CO吸附时，ENCUT=400时吸附能为-1.85 eV，升到500后稳定在-1.42 eV，差异达0.43 eV——这对于判断吸附强弱是致命的。后来看VASP官方论坛的讨论，确认金属体系的ENCUT应至少为ENMAX的1.3倍。这个经验在[DFT计算](https://www.keyanxueshu.com/category/dft/)项目中反复验证过。

KPOINTS网格：密度与计算量的平衡

K点网格的选取遵循”倒空间晶矢越长，K点越少”的原则。对于晶体结构预测类任务，本项目的经验公式是：

• 体相材料：a×b×c (Å) × K_x × K_y × K_z ≈ 30-50（金属）或 15-25（半导体/绝缘体）
• 表面模型（含真空层）：面内方向取体相的值，真空方向取1
• 分子/原子体系：Gamma点即可

具体到Cu(111)的p(2×2)表面模型（晶胞约5.1Å × 5.1Å × 25Å），面内用4×4×1的Monkhorst-Pack网格就足够收敛。但如果要做能带计算，需要用线模式的K点路径，这时Gamma-centered网格更合适。

一个容易忽略的问题：K点密度和ENCUT要同时收敛。只收敛一个而另一个偏低，结果可能比两个都偏低还差——因为误差方向不一致，无法对消。

SIGMA和ISMEAR：展宽方法的选择

金属体系用ISMEAR=1（Methfessel-Paxton）是标准做法，但SIGMA的选择容易被忽视。VASP手册给出一个判据：SIGMA × N_electronic_states_per_K-point 应远小于1 eV，同时检查outcar中的entropy T*S项，该值应小于总能的1 meV/atom。

本项目曾在Ni(111)体系的能带计算中踩过坑：ISMEAR=-5（四面体法）在能带计算中给出非常平滑的曲线，但如果K点不够密，费米面附近的态密度会出现非物理的尖峰。后来改用ISMEAR=0 + SIGMA=0.05配合更密的K点，问题才解决。

收敛判据：EDIFF和EDIFFG的配合

结构优化阶段有两个收敛判据：

• EDIFF：电子步自洽收敛（能量变化），默认1E-4，建议1E-5或1E-6
• EDIFFG：离子步收敛判据，正值表示能量收敛，负值表示力收敛

EDIFFG用负值（力收敛）比正值（能量收敛）更物理，因为力是能量的导数，力收敛意味着体系处于力的平衡点，而能量收敛只说明相邻两步能量差足够小，不代表真的到了极小值。

POTCAR选择：PAW vs Ultrasoft

VASP 5.4+ 默认提供PAW赝势，但同一元素往往有多个版本（如Cu有Cu、Cu_pv、Cu_sv）。选择原则：

• 默认版本（如Cu）：价电子最外层，计算量最小，适合常规体系
• _pv（如Cu_pv）：包含半芯p态，对磁性/光谱计算更准
• _sv（如Cu_sv）：包含半芯s态，最完整但计算量最大

本项目在Cu(111)的磁矩计算中发现，默认Cu赝势给出的磁矩比实验值低15%，改用Cu_pv后误差降到3%以内。但对于纯结构优化，默认版本完全够用，没必要增加计算负担。

常见报错与排查

参数检查清单

每次提交VASP任务前，本项目会过一遍以下清单：

1. ENCUT ≥ POTCAR中最大ENMAX × 1.3
2. KPOINTS在面内方向是否满足收敛（尤其表面模型）
3. ISMEAR与体系类型匹配（金属=1，半导体/绝缘体=0）
4. EDIFF和EDIFFG是否分阶段设置
5. POTCAR中元素顺序与POSCAR一致
6. 真空层≥15Å（表面模型）
7. NCORE或NPAR设置合理（通常NCORE=4-8）

VASP理论计算的参数调优没有万能模板，但理解每个参数的物理含义后，遇到新体系也能快速找到合理的初始值。以上经验均来自本项目的实际计算案例，不同体系可能需要微调，但基本框架是通用的。