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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 用自己生成的赝势算能带效果很好，算声子却一团糟
121/2返回列表12下一页
查看: 2425  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 中青基了要发朋友圈吗？ +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
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
关闭
请填处理意见
关闭