24小时热门版块排行榜    

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

zxgxbh2008

铜虫 (小有名气)

[求助] 静态计算

我在做表面问题时,为什么对硅超原胞结构优化时收敛了,但静态计算时怎么不收敛呢,以下是我的INCAR文件,希望各位帮我看看是什么问题,谢谢了
# SCF input for VASP
# Note that VASP uses the FIRST occurence of a keyword
# Extra input added by hand
#RWIGS=1.11
LREAL=Auto
NPAR=1
NELM=200
# End of extra input

SYSTEM = (Si)64  (P1) ~ 265 minimized: (Pearson #1817142)_1 1 0 0  surface (#1)_2x4x1_1 (VASP)

    PREC = Normal
   ENCUT = 400
  IBRION = -1
     NSW = 0
    ALGO = Normal (blocked Davidson)
    NELM = 60
  NELMIN = 2
   EDIFF = 1e-4
  EDIFFG = -0.02
VOSKOWN = 1
  NBLOCK = 1
   ISPIN = 1
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
   LWAVE = .TRUE.
  LCHARG = .TRUE.
ADDGRID = .FALSE.
  ISMEAR = -5
   SIGMA = 0.2
   LREAL = .FALSE.
RWIGS = 1.11
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao7910

金虫 (小有名气)

静态计算里面的是新生成的波函数文件,不知道为什么ismear=-5后,波函数迭代不能使能量收敛。而原来优化时最后一步得到的波函数基本上到了基态了。把那个考过来,程序就能很快迭代收敛。我这么理解的。
9楼2011-07-20 16:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

maomao7910

金虫 (小有名气)

哎,我也遇到,优化时,第一个离子步跑27部电子收敛,后面的离子步7,8步就收敛。做静态ismear= -5就60步也不收敛,dE在 E-01和E-02之间震荡。后面的dfpt还自己算呢。换了ismear有什么好的方法收敛呢?
2楼2011-07-19 19:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao7910

金虫 (小有名气)

【答案】应助回帖

zxgxbh2008(金币+8): 2011-07-20 11:03:00
我找到一个方法,把前面优化时最后收敛那步的wavecar文件考到静态计算文件夹中,istart用默认。这样能很快收敛。
3楼2011-07-20 09:46:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxgxbh2008

铜虫 (小有名气)

引用回帖:
Originally posted by maomao7910 at 2011-07-20 09:46:27:
我找到一个方法,把前面优化时最后收敛那步的wavecar文件考到静态计算文件夹中,istart用默认。这样能很快收敛。

wavecar文件是那个呀,我用的界面是Medea的,我不知道wavecar文件呢
4楼2011-07-20 11:04:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见