24小时热门版块排行榜    

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

gemucai

木虫 (正式写手)

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

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

dyc_2008

捐助贵宾 (著名写手)

张辽兄 功力很深厚啊
自己生成偃势
一布福田,二修菩提
10楼2012-07-07 21:18:10
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见