[VASP]求助碲化铋能带计算中的KPOINTS高对称点坐标

Views: 253|回复 Reply: 9

[VASP]

求助碲化铋能带计算中的KPOINTS高对称点坐标

[复制链接 Copy URL]

z57226

z57226

当前离线

UID59174

在线时间0 小时

注册时间2024-2-19

最后登录2025-8-6

好友0

22

帖子0

威望195

eV

Level 3 能力者, 积分 credits 217, 距离下一级还需 283 积分 credits

积分217

Level 3 能力者

发消息 Send PM

电梯直达 Go to

楼主

发表于 Post on 2025-6-19 11:13:50

|

只看该作者 Only view this author

|只看大图 View big picture

|倒序浏览 Reverse view

|阅读模式 Reading model

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

1.jpg

(656.02 KB, 下载次数 Times of downloads: 24)

下载附件 Download

2025-6-19 11:03 上传 Uploaded

收藏 Add to favorites0

回复 Reply

举报 Report

北大-陶豫

北大-陶豫

当前离线

UID19420

在线时间543 小时

注册时间2020-4-24

最后登录2025-8-15

好友11

602

帖子1

威望4788

eV

积分5410

Level 6 (一方通行)

发消息 Send PM

2#

发表于 Post on 2025-6-19 22:25:04

|

只看该作者 Only view this author

感觉没法回答，你对计算过程的描述太粗糙了，建议详细描述一下。

如果你在 vaspkit 的官方群（498683602），去群文件里找一下能带图模版和能带投影图模版，对照着里面的 readme 操作。

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯，我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。

我们怀揣热忱的灵魂天然被赋予对第一性的追求，不屑于单一坐标的约束，钟情于势能面彼端的芬芳。但

回复 Reply

举报 Report

z57226

z57226

当前离线

UID59174

在线时间0 小时

注册时间2024-2-19

最后登录2025-8-6

好友0

22

帖子0

威望195

eV

Level 3 能力者, 积分 credits 217, 距离下一级还需 283 积分 credits

积分217

Level 3 能力者

发消息 Send PM

3#

楼主 Author|

发表于 Post on 2025-6-20 11:36:05

|

只看该作者 Only view this author

北大-陶豫 发表于 2025-6-19 22:25

感觉没法回答，你对计算过程的描述太粗糙了，建议详细描述一下。

如果你在 vaspkit 的官方群（498683602） ...

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

回复 Reply

举报 Report

北大-陶豫

北大-陶豫

当前离线

UID19420

在线时间543 小时

注册时间2020-4-24

最后登录2025-8-15

好友11

602

帖子1

威望4788

eV

积分5410

Level 6 (一方通行)

发消息 Send PM

4#

发表于 Post on 2025-6-20 13:14:10

|

只看该作者 Only view this author

z57226 发表于 2025-6-20 11:36

十分感谢陶老师的回复！我是自己摸索着学的，因此有些地方不专业。我看了您提到的群里的文件里的readme， ...

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

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯，我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。

我们怀揣热忱的灵魂天然被赋予对第一性的追求，不屑于单一坐标的约束，钟情于势能面彼端的芬芳。但

回复 Reply

举报 Report

z57226

z57226

当前离线

UID59174

在线时间0 小时

注册时间2024-2-19

最后登录2025-8-6

好友0

22

帖子0

威望195

eV

Level 3 能力者, 积分 credits 217, 距离下一级还需 283 积分 credits

积分217

Level 3 能力者

发消息 Send PM

5#

楼主 Author|

发表于 Post on 2025-6-20 15:30:21

|

只看该作者 Only view this author

北大-陶豫 发表于 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轴的位置坐标和文献对不上，存在歧义误导阅读了，抱歉。

回复 Reply

举报 Report

北大-陶豫

北大-陶豫

当前离线

UID19420

在线时间543 小时

注册时间2020-4-24

最后登录2025-8-15

好友11

602

帖子1

威望4788

eV

积分5410

Level 6 (一方通行)

发消息 Send PM

6#

发表于 Post on 2025-6-21 10:01:00

|

只看该作者 Only view this author

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

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯，我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。

我们怀揣热忱的灵魂天然被赋予对第一性的追求，不屑于单一坐标的约束，钟情于势能面彼端的芬芳。但

回复 Reply

举报 Report

z57226

z57226

当前离线

UID59174

在线时间0 小时

注册时间2024-2-19

最后登录2025-8-6

好友0

22

帖子0

威望195

eV

Level 3 能力者, 积分 credits 217, 距离下一级还需 283 积分 credits

积分217

Level 3 能力者

发消息 Send PM

7#

楼主 Author|

发表于 Post on 2025-6-21 11:00:13

|

只看该作者 Only view this author

北大-陶豫 发表于 2025-6-21 10:01

只要判断相应的 kpath 路径是否跟文献一样长就行了。不用比横坐标位置是否一样，除非你画能带时的路径也是 ...

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

回复 Reply

举报 Report

z57226

z57226

当前离线

UID59174

在线时间0 小时

注册时间2024-2-19

最后登录2025-8-6

好友0

22

帖子0

威望195

eV

Level 3 能力者, 积分 credits 217, 距离下一级还需 283 积分 credits

积分217

Level 3 能力者

发消息 Send PM

8#

楼主 Author|

发表于 Post on 2025-6-26 14:04:08

|

只看该作者 Only view this author

北大-陶豫 发表于 2025-6-21 10:01

只要判断相应的 kpath 路径是否跟文献一样长就行了。不用比横坐标位置是否一样，除非你画能带时的路径也是 ...

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

2.jpg

(381.7 KB, 下载次数 Times of downloads: 22)

下载附件 Download

2025-6-26 13:52 上传 Uploaded

原胞，convetional cell以及结构优化了的convetional cell的计算结果

回复 Reply

举报 Report

北大-陶豫

北大-陶豫

当前离线

UID19420

在线时间543 小时

注册时间2020-4-24

最后登录2025-8-15

好友11

602

帖子1

威望4788

eV

积分5410

Level 6 (一方通行)

发消息 Send PM

9#

发表于 Post on 2025-6-26 15:26:07

|

只看该作者 Only view this author

你如果用vaspkit生成的kpath，仔细阅读它的输出提示就知道，能带必然是用原胞算的，并且请直接用vaspkit自动生成的PRIMCELL文件。

现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯，我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。

我们怀揣热忱的灵魂天然被赋予对第一性的追求，不屑于单一坐标的约束，钟情于势能面彼端的芬芳。但

回复 Reply

举报 Report

z57226

z57226

当前离线

UID59174

在线时间0 小时

注册时间2024-2-19

最后登录2025-8-6

好友0

22

帖子0

威望195

eV

Level 3 能力者, 积分 credits 217, 距离下一级还需 283 积分 credits

积分217

Level 3 能力者

发消息 Send PM

10#

楼主 Author|

发表于 Post on 2025-6-26 17:02:09

|

只看该作者 Only view this author

北大-陶豫 发表于 2025-6-26 15:26

你如果用vaspkit生成的kpath，仔细阅读它的输出提示就知道，能带必然是用原胞算的，并且请直接用vaspkit自 ...

了解了 谢谢陶老师

回复 Reply

举报 Report