24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2348  |  回复: 6

njdxllwl

至尊木虫 (正式写手)

[求助] QE软件计算电声子耦合,最后一步T_c出现NaN,什么原因呢?已有1人参与

请教大家个问题。
我用QE软件计算电声子耦合,计算得到的elph.0.0.0文件里面从第二行开始的那2个数是负值,进而导致最后处理计算Tc时,omega_log和T_c全是NaN,而如果把elph.0.0.0文件里面从第二行开始的那2个负数改为0,则能算出正的omega_log和T_c,请问这是什么原因导致的呢?如何解决这个负值的问题?
回复此楼

» 猜你喜欢

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

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

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
njdxllwl: 金币+5, 有帮助 2014-01-16 09:02:06
手动改为0即可。实际就应该是0.
真理是一点点接近的!
2楼2014-01-12 03:42:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njdxllwl

至尊木虫 (正式写手)

引用回帖:
2楼: Originally posted by gavinliu7390 at 2014-01-12 03:42:47
手动改为0即可。实际就应该是0.

谢谢回复!但是总感觉直接手动改为0有篡改数据的嫌疑。有没有办法避免为复数呢?
3楼2014-01-12 08:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

这不是篡改数据,是pwscf的bug。
正常就应该是0,只不过由于数值计算的原因,才会是负值。 只要这个负值小于最高频率的5%即可接受。在拟合声子谱的时候,pwscf会把Gamma点三个声学支强制设成0.
真理是一点点接近的!
4楼2014-01-16 06:22:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njdxllwl

至尊木虫 (正式写手)

引用回帖:
4楼: Originally posted by gavinliu7390 at 2014-01-16 06:22:04
这不是篡改数据,是pwscf的bug。
正常就应该是0,只不过由于数值计算的原因,才会是负值。 只要这个负值小于最高频率的5%即可接受。在拟合声子谱的时候,pwscf会把Gamma点三个声学支强制设成0.

谢谢!
5楼2014-01-16 09:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbrother

木虫 (小有名气)

我的是 elph.inp_lambda.1 文件中的第二行的前两个负数  改为0后结果出来了
6楼2014-06-05 07:12:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dx0620

木虫 (正式写手)

引用回帖:
4楼: Originally posted by gavinliu7390 at 2014-01-16 06:22:04
这不是篡改数据,是pwscf的bug。
正常就应该是0,只不过由于数值计算的原因,才会是负值。 只要这个负值小于最高频率的5%即可接受。在拟合声子谱的时候,pwscf会把Gamma点三个声学支强制设成0.

请问一下,对于*.dyn文件,如果出现几个很大的虚频怎么办?如何处理?
7楼2014-09-26 17:26:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 njdxllwl 的主题更新
信息提示
请填处理意见