24小时热门版块排行榜    

查看: 2129  |  回复: 5

边城雪

新虫 (初入文坛)

[求助] QE输出的矩阵元有很多的NaN已有1人参与

大家好,我用qe-6.3里的bands.x执行文件计算matrix elements of the momentum operator p between conduction and valence bands(需设置lp=.true.),
计算步骤是
pw.x <scf.in>scf.out
pw.x <band.in>band.out
bands.x <bands.in>bands.out
在第二步band.in的输入文件中,如果我使用的K点比较密的话,比如:
K_POINTS {crystal_c}
3
0.00    0.00    0.00      0
1.00    0.00    0.00      45
0.00    1.00    0.00      45
输出的matrix elements就会有很多的NaN(not a number);
但如果我用较疏的k-mesh,比如
K_POINTS {crystal_c}
3
0.00    0.00    0.00      0
1.00    0.00    0.00      15
0.00    1.00    0.00      15
输出矩阵元素都可以正常显示,我试了不同的服务器,结果都是如此,我也查询了一些资料,大多认为是数学库对角化的问题。但我不明白为什么会和设置的k点疏密有关?大家有计算这个或者遇到此类问题的吗?希望大家给一些建议,任何指点建议都可以,谢谢~
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天飘雪

版主 (知名作家)

认真做事,踏实做人

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
边城雪: 金币+30, 有帮助 2019-01-27 00:17:20

siestafonda

至尊木虫 (著名写手)

边城雪

新虫 (初入文坛)

漫天飘雪: 回帖置顶 2019-01-28 15:47:55

边城雪

新虫 (初入文坛)

边城雪

新虫 (初入文坛)

相关版块跳转 我要订阅楼主 边城雪 的主题更新
信息提示
请填处理意见