应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 电声耦合计算中，alpha2f.dat里面全是NaN
111/2返回列表12下一页
查看: 4251  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Local_Spin

铜虫 (小有名气)

[求助] 电声耦合计算中，alpha2f.dat里面全是NaN

各位，你们有没有遇到过这样的情况，elph.*.*.*文件里都是正常的，也没有虚频，到最后lambda.x那一步，得到的alpha2f.dat里面都是NaN，lambda.out里面的logw和TC也都是NaN。以前都没有遇到过这样的情况，很费解。

不知道各位有没有这样的经验，是怎么解决的？
回复此楼

» 猜你喜欢
1楼 2013-12-19 10:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

星期五的天气

金虫 (小有名气)
楼主你怎么解决的？我最近算的和您出现一样的情况呢。
赞一下(1人) 回复此楼
8楼2016-05-30 20:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gemucai

木虫 (正式写手)
★
sunyang1988: 金币+1, 谢谢提示 2014-01-03 10:07:25
确定没虚频？
回复此楼
2楼2013-12-19 11:01:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Local_Spin

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by gemucai at 2013-12-19 11:01:16
确定没虚频？

确定没虚频。lambda.out里面的第一列，也就是对各个q点lambda直接求和那一列，数字都是正常的，后面的列因为要涉及到a2f，于是就变成NaN了。很神奇
赞一下 回复此楼
3楼2013-12-19 11:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastfreedom

木虫 (小有名气)
★
sunyang1988: 金币+1, 谢谢交流 2014-01-03 10:07:33
引用回帖:
3楼: Originally posted by Local_Spin at 2013-12-19 11:42:00
确定没虚频。lambda.out里面的第一列，也就是对各个q点lambda直接求和那一列，数字都是正常的，后面的列因为要涉及到a2f，于是就变成NaN了。很神奇...

楼主可以看看自己的体系是不是绝缘体，如果费米能级附近的DOS为0
的话，在后面的lambda以及a2f 求和当中会出现除0的情况，这样可能就会出现NAN
赞一下 回复此楼
4楼2013-12-19 16:07:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Local_Spin

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by lastfreedom at 2013-12-19 16:07:16
楼主可以看看自己的体系是不是绝缘体，如果费米能级附近的DOS为0
的话，在后面的lambda以及a2f 求和当中会出现除0的情况，这样可能就会出现NAN...

谢谢你的回答。不过的确是金属。暂时没法编辑主楼，在此贴出lambda.out的内容。可以看到，第一列的lambda是有数值的，后面的N（EF）也不为0。然而，logw却是NaN。

=====lambda.out========


lambda = 2.158147 (       NaN)  <log w>=      NaNK  N(Ef)=19.138134 at degauss=0.005
     lambda = 1.850261 (       NaN)  <log w>=      NaNK  N(Ef)=19.177627 at degauss=0.010
     lambda = 1.622056 (       NaN)  <log w>=      NaNK  N(Ef)=18.059126 at degauss=0.015
     lambda = 1.440294 (       NaN)  <log w>=      NaNK  N(Ef)=16.482164 at degauss=0.020
     lambda = 1.286550 (       NaN)  <log w>=      NaNK  N(Ef)=14.947632 at degauss=0.025
     lambda = 1.155481 (       NaN)  <log w>=      NaNK  N(Ef)=13.614817 at degauss=0.030
     lambda = 1.045781 (       NaN)  <log w>=      NaNK  N(Ef)=12.522326 at degauss=0.035
     lambda = 0.956206 (       NaN)  <log w>=      NaNK  N(Ef)=11.654469 at degauss=0.040
     lambda = 0.884078 (       NaN)  <log w>=      NaNK  N(Ef)=10.974005 at degauss=0.045
     lambda = 0.826542 (       NaN)  <log w>=      NaNK  N(Ef)=10.439546 at degauss=0.050
lambda        omega_log          T_c
   2.15815           NaN                NaN
   1.85026           NaN                NaN
   1.62206           NaN                NaN
   1.44029           NaN                NaN
   1.28655           NaN                NaN
   1.15548           NaN                NaN
   1.04578           NaN                NaN
   0.95621           NaN                NaN
   0.88408           NaN                NaN
   0.82654           NaN                NaN
赞一下 回复此楼
5楼2013-12-19 20:26:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njdxllwl

至尊木虫 (正式写手)
引用回帖:
5楼: Originally posted by Local_Spin at 2013-12-19 20:26:22
谢谢你的回答。不过的确是金属。暂时没法编辑主楼，在此贴出lambda.out的内容。可以看到，第一列的lambda是有数值的，后面的N（EF）也不为0。然而，logw却是NaN。

=====lambda.out========


lambda = 2.158 ...

请问楼主知道什么原因了吗？偶也想知道。
赞一下 回复此楼
6楼2014-01-02 23:28:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastfreedom

木虫 (小有名气)
漫天飘雪: 回帖置顶 2018-04-10 08:32:36
引用回帖:
5楼: Originally posted by Local_Spin at 2013-12-19 20:26:22
谢谢你的回答。不过的确是金属。暂时没法编辑主楼，在此贴出lambda.out的内容。可以看到，第一列的lambda是有数值的，后面的N（EF）也不为0。然而，logw却是NaN。

=====lambda.out========


lambda = 2.158 ...

楼主在lambda.x 的输入文件第一行中：
emax  degaussq  ngaussq

ngaussq 取的是否为1？

如果是你去0  (gauss) 或者 -99 (fermi-dirac) 试试
在求a2f的积分中有个delta函数，用大于0阶mp方法可能会出现负值，
赞一下 回复此楼
7楼2014-01-04 18:41:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peoplemotu

新虫 (初入文坛)
请问楼主当时是怎么解决的？能贴一下方法吗？我遇到了相同的问题
赞一下 回复此楼
9楼2018-10-08 17:02:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哒哒虫

新虫 (初入文坛)
请问楼主当时是怎么解决的？我遇到了相同的问题
赞一下 回复此楼
10楼2018-11-15 09:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Local_Spin 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 怎么查啊 +4 huang1991js 2026-08-26 4/200 2026-08-26 11:26 by panjy.cn
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 能否退出参与的面上项目解除限项 +22 koalala 2026-08-24 25/1250 2026-08-26 11:14 by dxcharlary
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 8/400 2026-08-26 11:03 by zjlp
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +24 木水思豆 2026-08-25 27/1350 2026-08-26 08:26 by yanchen918
[基金申请] 范进中举一文的中心思想 +8 炎黄贵胄 2026-08-22 9/450 2026-08-26 07:47 by WASM
[基金申请] 2026国自然函评费到账 +21 羊腰板 2026-08-21 24/1200 2026-08-26 07:20 by finnigan
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
关闭
请填处理意见
关闭