计算化学公社

标题: 高对称点的kpath不连续时如何作出断开的能带图 [打印本页]

作者
Author: 珊珊来迟 时间: 2019-5-29 21:31
标题: 高对称点的kpath不连续时如何作出断开的能带图
(, 下载次数 Times of downloads: 78)
请问画这种kpath不连续的能带图时，X和M处的能带是断开的，有什么好的方法吗？我现在只能想到作图时找到X和M这两个高对称点对应的坐标，然后保持纵坐标不变，横坐标改为一样的，画出来的能带图就是断开的。但是对于路径中有多处断开的情况，一条条能带找着去改这个横坐标，工作量太大了。。想问问大家有什么好的处理的方法吗？
ps：在选择高对称点路径时我参考的是这篇文献Computational Materials Science 49 (2010) 299–312 (, 下载次数 Times of downloads: 56) ，这篇文献中给出的kpath很多都是像上图这种不连续的情况，但文献中很少见到这么复杂的能带图，请问有必要用这么复杂的kpath吗？用简单的kpath又担心不能完整的表现出晶体能带结构的基本信息，想请教大家如何看待这一问题？

作者
Author: 卡开发发 时间: 2019-5-29 21:59
本帖最后由卡开发发于 2019-5-29 22:08 编辑

我不太理解你的意思，你是计算过程中断开的能带不知道怎么做还是仅仅只是画图？如果你在计算的时候不知道断开的能带怎么做，按照你给的图，KPOINTS可以这么写：

KPOINTS
number of intersections
L
rec
x.xxx x.xxx x.xxx !Gamma
x.xxx x.xxx x.xxx !X
#你可以看到X连起来的部分是连续出现两次的
x.xxx x.xxx x.xxx !X
x.xxx x.xxx x.xxx !M
x.xxx x.xxx x.xxx !M
x.xxx x.xxx x.xxx !Gamma
x.xxx x.xxx x.xxx !Gamma
x.xxx x.xxx x.xxx !R
x.xxx x.xxx x.xxx !R
x.xxx x.xxx x.xxx !X
#X断开的部分X只出现一次
x.xxx x.xxx x.xxx !M
x.xxx x.xxx x.xxx !R

复制代码

这样的图做出来直接就是断开的。

作者
Author: granvia 时间: 2019-5-29 23:44
我以前用java写过对此进行自动处理的程序，但只适用于CRYSTAL程序。其实，这只是在表达上似乎好看一些罢了，对结果没有实际影响。最客观省事的做法是，不要用连线图画能带，直接用散点图即可

作者
Author: 珊珊来迟 时间: 2019-5-30 19:07
本帖最后由珊珊来迟于 2019-5-30 19:10 编辑

卡开发发发表于 2019-5-29 21:59
我不太理解你的意思，你是计算过程中断开的能带不知道怎么做还是仅仅只是画图？如果你在计算的时候不知道断 ...

老师您好，我是这样设置的kpoints，但我用origin作图时，X和M是连起来的，因为这两个点横坐标不同，比如某一条能带x（1,2），M（2,4），x和M两个点会有线连起来。除非我把X点坐标改为（2,2）这样才能画出我帖子中给出的那个图的效果（图中x和M是同一个横坐标）

作者
Author: 珊珊来迟 时间: 2019-5-30 19:09

granvia 发表于 2019-5-29 23:44
我以前用java写过对此进行自动处理的程序，但只适用于CRYSTAL程序。其实，这只是在表达上似乎好看一些罢了 ...

谢谢您的解答，确实对结果没有实际影响，但想要画出好看的图，看来不会写程序是硬伤

作者
Author: s541040349 时间: 2021-4-9 09:57

卡开发发发表于 2019-5-29 21:59
我不太理解你的意思，你是计算过程中断开的能带不知道怎么做还是仅仅只是画图？如果你在计算的时候不知道断 ...

卡卡老师，如果想把断开的图做成连续的这样反向处理可以吗

作者
Author: SSAW 时间: 2022-10-28 13:09

珊珊来迟发表于 2019-5-30 19:09
谢谢您的解答，确实对结果没有实际影响，但想要画出好看的图，看来不会写程序是硬伤

你好，请问你的问题解决了吗？我也遇到相同问题

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