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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 用自己生成的赝势算能带效果很好，算声子却一团糟
121/2返回列表12下一页
查看: 2428  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gemucai

木虫 (正式写手)

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

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

» 猜你喜欢

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

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的回帖

goldenfisher

金虫 (著名写手)
★
liliangfang: 金币+1, 谢谢交流 2012-04-08 16:56:18
重金属的semicore很难做，有很多技巧。
声子的计算，如果是高压，这很可能是内壳层电子的影响。你是NCPP的？那要算声子的话，ecutwfc必须足够大才对。
赞一下 回复此楼
5楼2012-04-07 23:43:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
5楼: Originally posted by goldenfisher at 2012-04-07 23:43:48:
重金属的semicore很难做，有很多技巧。
声子的计算，如果是高压，这很可能是内壳层电子的影响。你是NCPP的？那要算声子的话，ecutwfc必须足够大才对。

这个确实好难啊，迄今为止没发现一个能够让程序跑下去的截断半径。
赞一下 回复此楼
6楼2012-04-08 09:51:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
4楼: Originally posted by souledge at 2012-04-07 16:08:37:
Mo的赝势很不好生成，尤其考虑到半芯态时。
这里有两个从CPMD转来的LDA模守恒赝势，06的是没有半芯态，14的是有半芯态的。我近来在做其他的计算，还没时间测试这两个赝势的效果和差别，如果有兴趣，可以测试一下 ...

那两个赝势好像有问题。我测试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.
        (This warning will only be printed once per wavefunction)
      WARNING! Expected number of nodes:   0=     1-  0-  1, number of nodes found:  3.
       Setting wfc to zero for this iteration.
      (This warning will only be printed once per wavefunction)
    Message from routine run_pseudo:
      Errors in PS-KS equations
       WARNING! Expected number of nodes:   0=     2-  1-  1, number of nodes found:  3.
      Setting wfc to zero for this iteration.
      (This warning will only be printed once per wavefunction)

相反，对于Mo的阳离子组态却误差非常小：
n l     nl             e AE (Ry)        e PS (Ry)    De AE-PS (Ry)
3 2     4D   1( 5.00)       -0.78364       -0.78177       -0.00187
1 0     5S   1( 0.00)       -0.74605       -0.74501       -0.00105
2 1     5P   1( 0.00)       -0.44623       -0.44648        0.00024
-----------------------------------------------------------------------------------------
n l     nl             e AE (Ry)        e PS (Ry)    De AE-PS (Ry)
3 2     4D   1( 4.00)       -1.54772       -1.55039        0.00267
1 0     5S   1( 0.00)       -1.32339       -1.31887       -0.00452
2 1     5P   1( 0.00)       -0.94261       -0.94129       -0.00132

至于放在晶体中的情形还没有测试。
赞一下 回复此楼
7楼2012-04-17 11:10:34
已阅   回复此楼   关注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的回帖

tangosnow

木虫 (小有名气)
引用回帖:
8楼: Originally posted by souledge at 2012-04-17 11:59:18
这两个赝势是解析型转过来的，计算时你会看到wfc数为0，这个跟别的赝势不一样。所以最好还是放到晶体里测试一下吧～...

什么是解析型转过来的？
赞一下 回复此楼
9楼2012-07-07 04:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dyc_2008

捐助贵宾 (著名写手)
张辽兄 功力很深厚啊
自己生成偃势
赞一下 回复此楼
一布福田，二修菩提
10楼2012-07-07 21:18:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gemucai 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +3 瞬息宇宙 2026-08-31 6/300 2026-09-02 08:28 by 349506619
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
关闭
请填处理意见
关闭