24小时热门版块排行榜    

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

gemucai

木虫 (正式写手)

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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 个回答

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

goldenfisher

金虫 (著名写手)


liliangfang: 金币+1, 谢谢交流 2012-04-08 16:56:18
重金属的semicore很难做,有很多技巧。
声子的计算,如果是高压,这很可能是内壳层电子的影响。你是NCPP的?那要算声子的话,ecutwfc必须足够大才对。
5楼2012-04-07 23:43:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 本子写完了,给DS兄弟看了,得了92分 +3 Doma 2026-03-01 7/350 2026-03-02 00:00 by jnzsy
[考研] 一志愿郑大材料学硕298分,求调剂 +5 wsl111 2026-03-01 5/250 2026-03-01 23:45 by 暮雨星晴
[考研] 江苏省农科院招调剂1名 +3 Qwertyuop 2026-03-01 3/150 2026-03-01 23:18 by aaadim
[考研] 0854复试调剂 276 +3 wmm9 2026-03-01 3/150 2026-03-01 23:13 by 热情沙漠
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 8/400 2026-03-01 22:50 by jian_
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 0805总分292,求调剂 +7 幻想之殇 2026-03-01 7/350 2026-03-01 21:22 by 公瑾逍遥
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 12/600 2026-03-01 19:48 by 好好好1233
[考研] 272求调剂 +6 材紫有化 2026-02-28 6/300 2026-03-01 18:58 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 290求调剂 +9 材料专硕调剂; 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
信息提示
请填处理意见