当前位置: 首页 > 第一原理 >对Fe110面优化不收敛

对Fe110面优化不收敛

作者 guiwei2018
来源: 小木虫 500 10 举报帖子
+关注

大家好,本人近日在对一包含60个Fe原子的110面基底进行relax时,发现难以收敛,通过降低精度,增加步数也未能解决问题。xcysden打开in文件也没有发现结构错误,希望得到各种大佬的指点,以下给出了in文件。

对Fe110面优化不收敛
Fe.relax.png 返回小木虫查看更多

今日热帖
  • 精华评论
  • guiwei2018

    in文件没上传成功?

  • guiwei2018

    引用回帖:
    2楼: Originally posted by guiwei2018 at 2018-08-31 21:12:49
    in文件没上传成功?

    这是俯视图
    对Fe110面优化不收敛-1
    Fe.relax2.png

  • 张玲海

    增大nbnd 数目,提高degauss,减少conv_thr

  • guiwei2018

    引用回帖:
    4楼: Originally posted by 张玲海 at 2018-09-01 11:32:21
    增大nbnd 数目,提高degauss,减少conv_thr

    大神,我看到手册里说的是对于金属体系,for a metal, 20% more ,意思是超过电子数目一半的20%吗?比如480个电子,nbnd>240*1.2=288,我取了300不知道够不够。另外我把degauss从0.2增加到0.6,conv_thr减小到1.0*10-4

  • 张玲海

    引用回帖:
    5楼: Originally posted by guiwei2018 at 2018-09-01 21:04:49
    大神,我看到手册里说的是对于金属体系,for a metal, 20% more ,意思是超过电子数目一半的20%吗?比如480个电子,nbnd>240*1.2=288,我取了300不知道够不够。另外我把degauss从0.2增加到0.6,conv_thr减小 ...

    大概nbnd 为 380 ,conv_thr减小到1.0*10-5

  • guiwei2018

    引用回帖:
    6楼: Originally posted by 张玲海 at 2018-09-01 23:54:08
    大概nbnd 为 380 ,conv_thr减小到1.0*10-5...

    好的,我再试试!谢谢

  • guiwei2018

    引用回帖:
    6楼: Originally posted by 张玲海 at 2018-09-01 23:54:08
    大概nbnd 为 380 ,conv_thr减小到1.0*10-5...

    确实可以收敛。但是现在有两个问题,一是out文件里出现了SCF correction compared to forces is large: reduce conv_thr to get better values,其实我把conv_thr继续减小到1.0*10-6也可以收敛,但也会出现这样的错误。二是我的课题内容是做吸附,计算结合能,这就要求在优化“基底+吸附物”和‘基底’时,用同一套参数,现在是把优化“基底”的参数改了(增加nbnd=380,更换degauss=0.02为0.06),这样得到的结合能为-0.18与文献里-0.14不符合,请问大神,又该怎么办呢

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓