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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 用自己生成的赝势算能带效果很好，算声子却一团糟
汕头大学海洋科学接受调剂
51/1返回列表
查看: 2100  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gemucai

木虫 (正式写手)

[求助] 用自己生成的赝势算能带效果很好，算声子却一团糟

我用的是PWSCF,能带计算和声子计算都比较熟悉了。但最近我自己生成赝势之后，算出的能带与别人文献比较几乎无差别，但gamma点的声子频率却都是负的几千1/cm。我做过了优化，力和压力都几乎为0了，但还是有如此大的虚频。为何电子和声子计算会差别如此之大？有可能是我的赝势哪里没做好才导致了DFPT过程的异常？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-04-07 08:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
引用回帖:
7楼: Originally posted by gemucai at 2012-04-17 11:10:34:
那两个赝势好像有问题。我测试4p5 5s1组态的时候，ld1.x报错
WARNING! Expected number of nodes:   0=     3-  2-  1, number of nodes found:  3.
       Setting wfc to zero for this iteration.
       ...

这两个赝势是解析型转过来的，计算时你会看到wfc数为0，这个跟别的赝势不一样。所以最好还是放到晶体里测试一下吧～
赞一下 回复此楼
思想重于技巧，内涵重于表象
8楼2012-04-17 11:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+2, 专家考核, 谢谢指教 2012-04-07 10:11:38
gemucai: 金币+3, ★有帮助, 谢谢提醒！我又重启了去年没做完的课题，又要用Mo的赝势了。 2012-04-07 14:06:31
啥的赝势？就是你在邮件列表里说的那个C的全电子赝势么？一般能带计算不过是total energy的计算，对于体系的弛豫与否不是那么敏感~所以测试赝势时一般都是至少测试晶格常数、弹性常数的~至于声子，那是相当的敏感~
赞一下 回复此楼
思想重于技巧，内涵重于表象
2楼2012-04-07 10:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
2楼: Originally posted by souledge at 2012-04-07 10:05:19:
啥的赝势？就是你在邮件列表里说的那个C的全电子赝势么？一般能带计算不过是total energy的计算，对于体系的弛豫与否不是那么敏感~所以测试赝势时一般都是至少测试晶格常数、弹性常数的~至于声子，那是相当的敏感~

之所以用碳是想使问题简化。其实真正要生成的是钼的赝势，我没有做semicore，影响会这么大么？我测电子组态的时候误差已经只有几个mRy了。

当然没做semicore的原因是当我考虑4S和4P的时候，不管我用什么样的Rcut，ld1.x都会报错，想先用一个简单点的看看行不。
赞一下 回复此楼
3楼2012-04-07 14:02:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+2, 谢谢分享 2012-04-07 16:13:37
gemucai: 金币+7, ★★★很有帮助, 我又做了一个稍微好点的，会跟这两个比较。先谢谢了！ 2012-04-07 16:16:17
引用回帖:
3楼: Originally posted by gemucai at 2012-04-07 14:02:42:
之所以用碳是想使问题简化。其实真正要生成的是钼的赝势，我没有做semicore，影响会这么大么？我测电子组态的时候误差已经只有几个mRy了。

当然没做semicore的原因是当我考虑4S和4P的时候，不管我用什么样的R ...

Mo的赝势很不好生成，尤其考虑到半芯态时。
这里有两个从CPMD转来的LDA模守恒赝势，06的是没有半芯态，14的是有半芯态的。我近来在做其他的计算，还没时间测试这两个赝势的效果和差别，如果有兴趣，可以测试一下～
顺便，如果测试了，还希望早点发上测试结果图来，也有利于其他使用QE的人～免得赝势不好让太多人浪费了时间，或者免得大家都得自己跑测试～
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : Mo.LDA06.UPF
    • 2012-04-07 16:08:30, 295.5 K
  • 附件 2 : Mo.LDA14.UPF
    • 2012-04-07 16:08:34, 295.5 K
思想重于技巧，内涵重于表象
4楼2012-04-07 16:08:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研求调剂 +5 ban班小七 2026-04-11 5/250 2026-04-14 10:58 by dick_runner
[考研] 本科211，报考085601-310分 +13 ararak 2026-04-13 13/650 2026-04-14 09:28 by 852137818
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 37/1850 2026-04-14 07:39 by Abskk
[考研] 求调剂，985材料与化工348分 +9 涵竹刘 2026-04-11 14/700 2026-04-13 22:26 by 涵竹刘
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[材料工程] 材料调剂推荐 +9 蛋糕x2 2026-04-07 9/450 2026-04-13 09:07 by lhj2009
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 346，工科0854求调剂，专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 一志愿211，化学310分，本科重点双非，求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 314求调剂 +14 weltZeng 2026-04-09 14/700 2026-04-09 23:14 by wolf97
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 085404，334分，求调剂 +5 sunjie8888 2026-04-08 8/400 2026-04-09 07:26 by sunjie8888
[考研] 软件工程求调剂22软工296分求调剂，接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 生物学学硕，初试351分，求调剂 +4 …～、王…～ 2026-04-08 5/250 2026-04-08 21:49 by limeifeng
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
信息提示
关闭
请填处理意见
关闭