24小时热门版块排行榜    

查看: 2416  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gemucai

木虫 (正式写手)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +10 wangmingqi 2026-08-28 12/600 2026-08-28 22:06 by water_still
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 17:09 by winsaint
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 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
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +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
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见