vasp计算二维磁性材料（在kz方向有厚度的），通过什么vasp计算方法有可能计算出沿kz方向的色散？

vasp计算二维磁性材料，在kz方向有厚度的，是体相切的薄层，通过什么vasp计算方法有可能计算出沿kz方向的色散？可以显示体相在kz方向的性质？因为实验上可以判断，计算上不知道可以有方法判断吗？

图片

问题链接：https://www.zhihu.com/question/1903940826687582989

 

这个问题其实非常典型，很多刚接触二维材料和 VASP 的同学都会混淆“二维体系”“薄层体系”“体相 kz 色散”“ARPES 中的 kz 分辨”这些概念。

如果你的 VASP 模型是真正的二维 slab（即 c 方向加入了大真空层），那么严格来说是不存在真正 kz 色散的。

因为：

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

很多初学者会误以为：

“结构在 z 方向还有厚度，所以应该还能有 kz 色散”

但这里需要区分：

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

例如：

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

虽然它们在 z 方向有实际原子厚度，

但是：

只要模型中加入了真空层（vacuum），
体系就已经不再是一个真正的三维周期晶体。也不存在严格意义上的 (E(kz)) 色散。

你在 VASP 里即使设置：

KPOINTS
0
Gamma
1 1 20

看到的也只是：

  • 数值上的重复采样，
  • 而不是真实物理 kz 色散。

因为真空区域没有电子态传播。

k-path:

Γ-Z
A-L
M-K

这时：
E=E(k_x,k_y,k_z)

中的 (k_z) 才是真实存在的。

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

对于 slab 体系，常见做法有：

方法 1：对比不同厚度 slab

例如：

  • 1 layer
  • 3 layer
  • 5 layer
  • 7 layer

观察：

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

这是最标准的方法。

方法 2：同时计算 bulk 和 slab

这是最推荐的方法。

步骤：

bulk 计算得到：E(k_x,k_y,k_z)
特别关注：

  • Γ-Z
  • A-L
  • K-H

等 kz 路径。

 

 

（2）slab 计算得到：

  • 表面态
  • 有限厚度子带
  • 磁性重构

然后比较：

  • slab 子带是否来源于 bulk kz band folding；
  • 某些态是否是 surface-localized；
  • 是否存在 surface Dirac cone；
  • 是否有层间 AFM 导致的 kz 折叠。

这也是拓扑材料、二维磁体中最常见的分析思路。

 

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

这是很多人最困惑的地方。

ARPES 中的 kz 并不是直接测量得到的。

实验中的：

k_z通常来自通过光电子末态模型推算得到。

因此：

  • bulk 样品能扫描 photon energy 得到 kz dispersion；
  • 真正二维材料通常几乎没有 photon-energy dependence；
  • 如果仍存在明显 kz dependence，
    往往说明体系还具有三维性。

这恰恰是实验判断二维/三维电子结构的重要依据。