计算化学公社

标题: 求助碲化铋能带计算中的KPOINTS高对称点坐标 [打印本页]

作者
Author: z57226 时间: 2025-6-19 11:13
标题: 求助碲化铋能带计算中的KPOINTS高对称点坐标
在计算碲化铋能带时，发现vaspkit生成的高对称点路径和文献中常用的路径不一致。因此参考文献中的布里渊区和路径（图2上）对照seekpath生成的布里渊区和高对称点结果（图1），选取了我需要的高对称点（图1下中我框选的）。随后修改KPOINTS，进行计算。结束后读取 KLABELS中高对称点坐标作图。结果发现我得到的高对称点坐标（图2中）和文献（图2上下）中的完全不一致。想请教一下各位大佬是我在找对应的高对称点这一步中找的点错了吗还是我计算的其他设置有问题呢？这个问题困扰很久了，望各位大佬不吝赐教，十分感谢！

作者
Author: 北大-陶豫 时间: 2025-6-19 22:25
感觉没法回答，你对计算过程的描述太粗糙了，建议详细描述一下。
如果你在 vaspkit 的官方群（498683602），去群文件里找一下能带图模版和能带投影图模版，对照着里面的 readme 操作。

作者
Author: z57226 时间: 2025-6-20 11:36

北大-陶豫发表于 2025-6-19 22:25
感觉没法回答，你对计算过程的描述太粗糙了，建议详细描述一下。
如果你在 vaspkit 的官方群（498683602） ...

十分感谢陶老师的回复！我是自己摸索着学的，因此有些地方不专业。我看了您提到的群里的文件里的readme，里面提到了能带作图，从KLABELS文件里可以获得选用的高对称点的坐标（图2中图x轴的数字即是我从KLABELS中得到的）。我现困惑的是我获得的高对称点坐标（图2中）和文献中（图2上下）的不一致（图2中的竖线为参照线）。假设L点坐标是3.171的话，别人的Z点坐标为0.4左右，我的是0.1。我能想到的原因就是高对称点及路径搞错了。我想请教一下这会是其他参数设置错误引起的吗？望陶老师不吝赐教！十分感谢！

作者
Author: 北大-陶豫 时间: 2025-6-20 13:14

z57226 发表于 2025-6-20 11:36
十分感谢陶老师的回复！我是自己摸索着学的，因此有些地方不专业。我看了您提到的群里的文件里的readme， ...

什么叫“假设L点坐标是3.171的话，别人的Z点坐标为0.4左右，我的是0.1”？

作者
Author: z57226 时间: 2025-6-20 15:30

北大-陶豫发表于 2025-6-20 13:14
什么叫“假设L点坐标是3.171的话，别人的Z点坐标为0.4左右，我的是0.1”？

图2中间那张图中0.104,.0.923，3.171这三个值分别是我算出来的Z、F、L高对称点在能带图X轴的位置坐标（从KLABELS中读取的）。我作图以后发现和文献对不上（图2中的竖线为文献中的高对称点在能带图X轴的位置）。我用‘假设......’主要是为了表达我获得的高对称点在能带图X轴的位置坐标和文献对不上，存在歧义误导阅读了，抱歉。

作者
Author: 北大-陶豫 时间: 2025-6-21 10:01
只要判断相应的 kpath 路径是否跟文献一样长就行了。不用比横坐标位置是否一样，除非你画能带时的路径也是 Γ-Z-F-Γ-L。如果你算的路径除了文献的路径外又加了别的路径，可以用群文件里的 BandSimplify 这个脚本。

作者
Author: z57226 时间: 2025-6-21 11:00

北大-陶豫发表于 2025-6-21 10:01
只要判断相应的 kpath 路径是否跟文献一样长就行了。不用比横坐标位置是否一样，除非你画能带时的路径也是 ...

好的了解了，谢谢陶老师的回复！

作者
Author: z57226 时间: 2025-6-26 14:04

北大-陶豫发表于 2025-6-21 10:01
只要判断相应的 kpath 路径是否跟文献一样长就行了。不用比横坐标位置是否一样，除非你画能带时的路径也是 ...

陶老师您好！我后面又研究了一下，发现上面那个问题是因为我用的convetional cell计算引起的。后续我用原胞做了一次计算，高对称点路径在X轴的坐标和文献都对的上了。下图是我对比原胞、convetional cell和结构优化过的convetional cell的晶体结构和KLABELS结果（计算的INCAR等文件都是vaspkit自动生成的，除了KPOINTS是我提供的）。那这样的话，是不是算能带计算最好是在原胞上进行呢？原胞和convetional cell的KLABELS存在明显差别，那后续我算掺杂体系要扩胞的话（针对这个体系），是否只能基于原胞进行扩胞？望陶老师不吝赐教！十分感谢！

作者
Author: 北大-陶豫 时间: 2025-6-26 15:26
你如果用vaspkit生成的kpath，仔细阅读它的输出提示就知道，能带必然是用原胞算的，并且请直接用vaspkit自动生成的PRIMCELL文件。

作者
Author: z57226 时间: 2025-6-26 17:02

北大-陶豫发表于 2025-6-26 15:26
你如果用vaspkit生成的kpath，仔细阅读它的输出提示就知道，能带必然是用原胞算的，并且请直接用vaspkit自 ...

了解了谢谢陶老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)