24小时热门版块排行榜    

查看: 1972  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

newton3915

铁杆木虫 (正式写手)

[求助] pw计算能带图

用pw计算bands的时候经常会得到如图1所示的连线图,会有部分的本征值连接很乱,有时多有时少,但是用点图的话(图2)看起来会好很多,其中scf和nscf的计算都没问题,请教一下各位大侠,这种情况该怎么办呢?是不是哪里出问题了?
pw计算能带图
图1


pw计算能带图-1
图2
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newton3915

铁杆木虫 (正式写手)

引用回帖:
3楼: Originally posted by KalaShayminS at 2013-05-21 00:57:18
QE的bands计算后面用了基于波函数相似性和对称性来判断能级之间的连接顺序,很遗憾这个方法并不是总是可靠的,很多相去甚远的能带也可能是相似的。
最直接的方法就是读取pw.x输出的能量值来画图,而不需要bands.x进 ...

你指的是按照能量值从小到大排列,然后连线成图,不管该能量值在不在所谓一个能带上?
6楼2013-05-22 10:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

xhzha

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流, 2013-05-20 20:58:29
franch: 金币+2, 谢谢回帖交流, 2013-05-20 21:09:01
newton3915: 金币+2, ★★★很有帮助 2013-05-22 11:03:10
这个我算声子谱时极少情况出现过,应该是5.0版本的吧?我猜会不会是个bug,里面是根据简正坐标正交来区分一条条声子带,能带应该类似吧,不知道为什么会出现这种情况,看源程序也没看出问题。我后来想了个办法就是根据一阶导数连续重新排了下数据,很多时候能得到很好的结果。我也期待有高人指点怎么处理这个问题。如果你只想得到好看的曲线,4.3之前的版本应该都没有这个问题,因为它是按能量大小来排布每条带的。
2楼2013-05-20 20:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

KalaShayminS

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
newton3915: 金币+4, ★★★很有帮助 2013-05-22 11:02:51
QE的bands计算后面用了基于波函数相似性和对称性来判断能级之间的连接顺序,很遗憾这个方法并不是总是可靠的,很多相去甚远的能带也可能是相似的。
最直接的方法就是读取pw.x输出的能量值来画图,而不需要bands.x进行处理。
3楼2013-05-21 00:57:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mailoliver

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
newton3915: 金币+4, ★★★很有帮助 2013-05-22 11:03:01
一劳永逸的方法就是自己写个小程序,每个k点的本征值读出来,重新排序一下,写到一个新的文件中,然后用origin画图。
4楼2013-05-21 17:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 10:52 by xSrQNSSjcgoX
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 07:23 by xSrQNSSjcgoX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:33 by xSrQNSSjcgoX
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:38 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 6/300 2026-11-15 02:59 by xSrQNSSjcgoX
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:24 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:13 by Equinoxhua
[考博] 27博士申请(市政/环境) +3 ladf 2026-11-12 3/150 2026-11-13 08:19 by Equinoxhua
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:27 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +7 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 10:03 by HajghKOMynmz
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
信息提示
请填处理意见