VASP计算单点能 — 静态能量计算的精度控制与工程规范

单点能计算——固定原子坐标不动，做一步自洽迭代得到体系的总能量——是VASP中使用频率最高、看起来最简单的操作。但”简单”带来的问题是：很多人直接拿别人INCAR文件里的参数跑，从不做收敛性测试，然后拿着精度在0.1 eV量级的数据去论证一个只有0.05 eV差异的物理结论。VASP计算单点能的精度控制不是玄学，而是一套有明确测试流程的工程规范。这篇复盘把每一步的精度来源拆清楚。

一、截断能收敛测试：最基础也最容易被跳过

截断能（ENCUT）是控制平面波基组完备性的核心参数。ENCUT越高，基组越完备，总能量越低（变分原理）——但计算量也越大。关键不是选一个绝对值，而是确认你选的ENCUT能让目标物理量的计算误差小于你关心的物理效应大小。

收敛测试的标准做法：取一个典型的输入结构，固定k点和所有其他参数，从300 eV开始以25 eV步长递增ENCUT直到600-700 eV，记录每个ENCUT下的总能量。画出E_tot vs ENCUT曲线，找到能量变化小于1 meV/atom的ENCUT值。

不同元素的PAW势函数推荐的默认ENCUT不同（如O是400 eV，W是223 eV的ENMAX×1.3≈290 eV），但默认值只是”可用”级别——如果你想做结合能计算（精度要求<0.01 eV/atom），至少取ENMAX×1.5。对于含第一周期过渡金属（3d元素）的体系，ENCUT通常需要550-600 eV才能收敛到1 meV/atom。

二、k点收敛测试：总能量对k点的敏感度比你想象的低

VASP计算单点能时，k点密度影响总能量的方式与截断能不同。总能量对k点密度的收敛行为不是单调的——加密k点不保证总能量始终下降，而是趋于一个饱和值。

一个常见误解是认为k点越密精度越高，所以无脑用很高的k点（如15×15×15）。实际上，对于大多数周期性固体，总能量在k点间距<0.03 Å⁻¹时就收敛到1 meV/atom以内。以晶格常数4 Å的立方晶胞为例，6×6×6 (k点间距0.025 Å⁻¹)通常已经够用，更高密度的收益递减。

对于表面slab模型的计算，k点只在与表面平行的方向（k_x、k_y）做采样，垂直方向（k_z）用1即可（如8×8×1）。这是因为slab模型在垂直方向有真空层，布里渊区在该方向的色散几乎为零。

收敛测试做法：固定ENCUT和其他参数，从4×4×4到16×16×16测试k点对总能量的影响。如果8×8×8和10×10×10的E_tot差异<1 meV/atom，选8×8×8就够了。

三、展宽方法与SIGMA值的选择

展宽方法（ISMEAR）和展宽宽度（SIGMA）的选择在单点能计算中非常重要——因为在展宽不为零时，VASP输出的”Total free energy”不是严格的基态能量，它包含一个熵项：

F = E₀ − TS

其中F是VASP输出的”free energy TOTEN”，E₀是外推到零展宽的”energy without entropy”（即”energy(sigma→0)”）。对于需要精确总能量的应用（结合能、吸附能、形成能），应该使用”energy(sigma→0)”而非”free energy TOTEN”。

ISMEAR选择指南：

• 半导体/绝缘体：ISMEAR=-5（四面体方法），是最精确的，不引入任何展宽熵。但需要足够的k点（>4×4×4）。
• 金属：ISMEAR=1（Methfessel-Paxton，N=1阶），SIGMA=0.1-0.2。SIGMA太大会让”free energy”与”energy(sigma→0)”差距过大（>0.01 eV/atom），SIGMA太小可能导致电子步收敛困难。
• ISMEAR=0（Gaussian展宽）是通用的安全选项，但对金属不如Methfessel-Paxton精确，对绝缘体不如四面体方法。

四、自旋极化与磁矩初始化

许多体系在非磁性（ISPIN=1）计算中能正常收敛，但关闭自旋极化可能让你漏掉能量上更稳定的磁性基态。对含Fe、Co、Ni、Mn及大多数3d过渡金属的体系，默认必须开ISPIN=2。

初始磁矩的设定（MAGMOM标签）对收敛路径有显著影响。如果初始磁矩设反了方向（如本该铁磁耦合的原子设成了反铁磁），VASP可能收敛到错误的局域极小（亚稳态），给出的总能量比真实基态高0.1-0.5 eV/atom。

好的做法是：对过渡金属原子，初始磁矩设为其原子磁矩（如Fe=5, Co=3, Ni=2），对非磁原子设0。对反铁磁体系，相邻同种原子的初始磁矩设相反符号。在输出中检查磁矩是否稳定且合理，如果不合理，换初始猜测重跑。

五、单点能与结构弛豫的先后关系

一个经典错误：用低精度参数做结构优化（省算力），然后用高精度参数在同一几何结构上做单点能计算。这会导致”精度不匹配”——你优化的结构是在低精度势能面上的极小点，但用高精度势能面评估这个点时，它可能不再是极小点。

正确的做法是：用与最终单点能计算相同或相近的精度参数做结构优化。如果最终要用HSE06做单点能，结构优化也应该用HSE06（或至少在PBE优化的结构上确认原子受力在HSE06下仍然接近零）。

对于大多数应用，PBE优化+HSE06单点的组合是可以接受的，只要确认HSE06单点计算中原子受力<0.05 eV/Å。如果受力超标，需要在HSE06下继续弛豫。

六、VASP单点能计算的完整检查清单

一套标准的VASP计算单点能工程checklist：

1. ☐ ENCUT至少为所有元素ENMAX值×1.3（推荐1.5），已做收敛测试
2. ☐ k点间距<0.03 Å⁻¹，已做收敛测试
3. ☐ ISMEAR选择与体系匹配（半导体=-5，金属=1）
4. ☐ 如果需要精确能量，使用energy(sigma→0)而非free energy
5. ☐ ISPIN=2用于含磁性元素的体系，MAGMOM设置了合理的初始磁矩
6. ☐ 结构弛豫用的参数精度与单点能计算一致（或至少确认高精度下的受力仍然足够小）
7. ☐ PREC=Accurate（对于高精度需求，这会启用更密集的实空间网格）
8. ☐ 自洽计算的电子步收敛到EDIFF=1E-6 eV（而非默认的1E-4）用于高精度能量
9. ☐ LREAL=.FALSE.（避免实空间投影带来的额外数值误差，对高精度单点能）

其中第8条EDIFF=1E-6经常被忽略。将EDIFF从默认1E-4收紧到1E-6，总能量的变化可以到0.01 eV/atom——对于结合能在几个eV量级的计算，这个误差占比远小于1%，但对于需要区分meV级别能量差异的计算（如相变能量差、磁各向异性能），1E-4的默认精度是不够的。

配图建议：

1. 截断能收敛测试曲线（ALT：”VASP计算单点能截断能收敛测试总能量vs ENCUT收敛行为1meV/atom”）：展示E_tot随ENCUT变化的收敛曲线。
2. k点收敛测试图（ALT：”VASP单点能k点收敛测试总能量vs k点网格密度”）：展示E_tot随k点密度变化的收敛行为。
3. 精度参数对比柱状图（ALT：”VASP单点能不同精度参数设置PBE PBE+U HSE06能量对比”）：对比不同泛函/参数下的计算能量差异。