关于能带计算的一些疑问

yoghurt117

最近一直在初学cp2k的能带计算，因为网上用cp2k计算的例子不多，结合查到的例子和QE等软件的关于石墨烯例子，感觉有些疑问，想请教一下群里的高人，

1. 因为查到的例子输入的结构参数都不一样，但石墨烯通常有60度的结构和120度的结构，左边结构和右边结构选取的基矢量不同，我按照两种结构算了一下倒格子，得到了k点分别列在图下面，计算能带的话，是不是 应该是按照不同的输入结构，来输入不同的k点路径吧？




2. 第二个疑问是，我要做掺杂或者缺陷能带计算的话，就得构建超胞模型，但是超胞构建好之后，做完结构优化，k点路径的选择还是按照上面两种方式来选择k点路径么？因为涉及到能带折叠问题，所以对于超胞的计算不是很明白；

3. 第三个 问题，我做了一个 磷掺杂的超胞优化，结果 如下图 ，

我把结构提交给了 seekpath，结果系统给出了一堆k点 ，因为结构优化结果导致对称性和完美石墨烯不同，所以系统给出的k点多了，这个可以理解，但是看有的文献依然用的是上面4个k点来算的 ，所以我不太明白是必须要按照seekpath给出的来算还是说可以按照完美石墨烯的上面gamma-->M-->K-->gamma来算 。

谢谢！

阮文升-河内国大

1. 是的。你可以利用 VASPKIT，自动调整到标准取向（甚至能帮你自动识别出超胞中的原胞）然后自动生成 kpath，见：
cp2k利用纯泛函进行能带计算的简单方法
http://bbs.keinsci.com/forum.php ... 4&fromuid=89463
(出处: 计算化学公社)

2. 我不清楚。

3. 磷的偏心位移会降低对称性，降低到何种程度取决于磷的偏心位移方向。例如，如果只在z方向上有偏心位移，就不会影响xy方向上的对称性。当然，kpath 的选取，也取决于你是否需要计算 z 方向的能带（一般不需要）。

yoghurt117

阮文升-河内国大 发表于 2026-7-26 21:17
1. 是的。你可以利用 VASPKIT，自动调整到标准取向（甚至能帮你自动识别出超胞中的原胞）然后自动生成 kpat ...

好的，太感谢了，
我回头查一下multiwfn手册