在上一章《第七章：KPOINTS文件结构解析！| 2026新版VASP基础教程》中，华算科技朱老师详细介绍了VASP的POTCAR文件。VASP输入文件中的KPOINTS提供了计算所需要的数值采样网格。本章将正式介绍KPOINTS，从文件内容介绍、采样方法、采用密度等角度来介绍VASP计算所需要的K点网格。
VASP使用平面波基组展开电子波函数。由于晶体的周期性，电子态在倒空间是连续的，而我们只能离散采样这些点。K点网格密度决定了倒空间积分的精度：
费米面附近电子态变化剧烈，需要极其稠密的K点网格才能准确描述电子分布。
带隙将费米面隔开，电子态分布相对平缓，所需K点密度相对较低。
当原胞被复制成大超胞时，倒空间相应收缩，K点网格可以相应稀疏
最常用的判断标准是总能量变化（Total Energy per Atom）：
随着K点数增加的变化幅度
，即认为K点已收敛。
除了总能量，还可以关注带隙（Band Gap）、态密度（DOS）形状、费米能级（Fermi Energy）以及力（Force）的变化。对于声子谱或磁性计算，磁矩（Magnetic Moment）也是重要指标。
在进行K点测试时，必须保持所有其他计算参数不变，尤其是平面波截断能（ENCUT）。常用的做法是先进行截断能收敛测试，确定一个较高的ENCUT值，然后再进行K点测试。
从一个非常稀疏的网格开始，逐步增加密度。对于立方晶格，常见的测试序列如下：
|
测试序列 |
K点网格 (Monkhorst-Pack) |
|
粗糙测试 |
2x2x2, 3x3x3 |
|
中等密度 |
4x4x4, 5x5x5 |
|
精细收敛 |
6x6x6, 8x8x8, 10x10x10 |
对于低维材料（如石墨烯、薄膜），应保持垂直方向（z方向）K点为1（即1x1x1），只在平面内采样（如6x6x1）。
每完成一次计算，记录以下信息：
KPOINTS：网格大小
E_TOTAL：总能量（OUTCAR文件中）
E_PER_ATOM：每原子的能量
ΔE：相邻两次K点的能量差
Fermi Energy：费米能级（对于金属特别重要）
将 K点数（或总的K点数）作为横轴，总能量（每原子）作为纵轴绘制曲线。您会发现能量随着K点密度增加迅速下降，然后趋于平缓。当曲线在某一点之后基本平直，且能量变化小于1meV/atom时，即可确定收敛的K点网格。
VASP在生成K点网格时会自动考虑晶体对称性，可能会将一个网格简化为对称等价的K点数。例如，3x3x3的网格可能会简化为一个K点（Gamma点）。因此，建议在KPOINTS文件中使用 Gamma 选项或 Monkhorst-Pack 选项明确指定网格，而不是手动计算K点数。
金属计算中，由于电子占据的阶跃性，必须使用平滑技术（Smearing）来辅助收敛。常用的参数为：
ISMEAR = 1 或 2（Methfessel-Paxton）
ISMEAR = -1（Fermi smearing）
ISMEAR = -5（tetrahedron method with Blöchl corrections，适用于DOS计算）
SIGMA：平滑宽度，金属通常取 0.1-0.2 eV。
对于包含数百个原子的超胞模型，通常不需要很大的K点网格。经验法则是：
例如，如果原胞（2个原子）需要 6x6x6 的网格（共 216 个K点），那么 108 个原子的体系可能只需要 2x2x2 的网格即可。
密度泛函理论的态密度（DOS）计算对K点极其敏感。即使总能量已经收敛，为了获得光滑且无噪声的DOS曲线，往往需要显著更密的K点网格（如 20x20x20 或更高）。
对于单层或薄膜材料（如石墨烯、MoS2），应采用二维K点网格：
平面内：使用如 12x12x1 的网格。
垂直方向：设置为 1（避免不必要的计算）。
为什么需要进行K点收敛测试：介绍了金属与半导体体系测试的意义
如何判断K点是否收敛：介绍了收敛标准与检测量
收敛测试的具体步骤：具体步骤为固定参数，增加k点，记录能量曲线
下一章将正式引入本次教程的核心—K点网格设置经验。我们将从基础概念与文件结构、K点密度的初始选择、体系类型对K点设置的影响方面详细介绍VASP KPOINTS文件，以及他们在VASP计算中的应用，敬请期待！
