VASP如何计算二维材料kz方向色散

这是一个非常常见的问题，许多刚接触二维材料和VASP的用户，容易将“二维体系”、“薄层体系”、“体相kz色散”以及“ARPES中的kz分辨”等概念混淆。

首先给出核心结论：如果VASP模型是真正的二维slab（即在c方向添加了足够大的真空层），那么严格意义上，并不存在真正的kz色散。

背后的原理很简单：

• 在VASP中加入真空层后，周期性在z方向被人为切断。
• 布洛赫波矢k_z不再对应真实晶体中的传播自由度。
• 因此，能带不会具有连续的kz色散。

很多初学者误以为：“结构在z方向还有厚度，所以应该还能有kz色散。”这里需要区分一个关键概念：

一、“有厚度” ≠ “具有kz周期性”

举例说明：

• 单层Fe3GeTe2
• 几层CrI3
• 从体相切出来的NbSe2 slab

这些体系在z方向确实存在原子厚度，但只要模型中加入了真空层，就不再是真正的三维周期晶体。这种情况下，也不存在严格意义上的电子能量随k_z的色散关系。

即使在VASP的KPOINTS文件中将k点设置为：
KPOINTS
0
Gamma
1 1 20

看到的也只是数值上的重复采样，而非真实的物理kz色散——因为真空区域没有电子态传播。

真正需要k_z路径时，应对体相（bulk）结构进行计算，此时k路径可以为：
Γ-Z | A-L | M-K

这种情况下，E = E(kx, ky, kz)中的k_z才是真实存在的。

VASP中如何“间接体现kz信息”

对于slab体系，虽然无法直接计算kz色散，但可以通过以下“曲线救国”方法：

方法1：对比不同厚度的slab

分别计算1层、3层、5层、7层的slab，然后观察：

• 能带数量的增加；
• 子带是否逐渐变密；
• 最终是否逼近bulk的kz色散。

这是最标准且直观的做法。

方法2：同时计算bulk和slab

这是最推荐的方法，操作清晰：

1. bulk计算获得完整的E(kx, ky, kz)，特别关注Γ-Z、A-L、K-H等kz路径。
2. slab计算专门处理表面态、有限厚度子带以及可能的磁性重构。
3. 将两者对比：slab的子带是否来源于bulk kz的能带折叠？某些态是否局域在表面？是否存在表面狄拉克锥？是否有层间反铁磁导致的kz折叠？

这套分析思路在拓扑材料和二维磁体研究中属于常规操作。

实验上的kz分辨，为什么还能看到

这也是很多人困惑的地方。ARPES中的k_z并非直接测量得到。

实际实验中，k_z通常通过光电子末态模型推算出来。因此：

• 对于体相样品，可以通过扫描光子能量得到kz色散；
• 对于真正的二维材料，通常几乎看不到光子能量依赖性；
• 如果实验上仍观察到明显的kz依赖性，反而说明体系具有三维性。

这正是实验上判断二维或三维电子结构的关键依据。