计算化学公社

标题: VASP结构优化收敛问题 [打印本页]

作者
Author: 瓜大小瓜子 时间: 2019-5-3 20:20
标题: VASP结构优化收敛问题
各位老师，请教关于vasp结构优化的能量以及力的收敛问题！
能量如下：
free  energy TOTEN  =    -141.89045799 eV
  free  energy TOTEN  =    -141.89042047 eV
  free  energy TOTEN  =    -141.89050665 eV
  free  energy TOTEN  =    -141.89053585 eV
  free  energy TOTEN  =    -141.89055379 eV
  free  energy TOTEN  =    -141.89054958 eV
  free  energy TOTEN  =    -141.89058764 eV
  free  energy TOTEN  =    -141.89060737 eV
  free  energy TOTEN  =    -141.89062235 eV
  free  energy TOTEN  =    -141.89062464 eV
  free  energy TOTEN  =    -141.89064300 eV
  free  energy TOTEN  =    -141.89065862 eV

我计算的是块体，24个原子，INCAR如下
   ISTART = 0
   ICHARG = 2
   PREC = Accurate
   EDIFF = 1e-4
   EDIFFG  = -0.01
   NPAR =  8
   LPLANE = .TRUE.
   ENCUT = 450
   NELMIN  = 3
   LREAL = Auto
   ISYM = 0
   NWRITE  = 1
   NGX =46
   NGY =46
   NGZ =66
   IVDW=1
   VDW_RADIUS=30.0
   VDW_S6=0.75
      VDW_SR=1.00
   VDW_D=20.0
   VDW_C6=10.800 0.7000  0.1400  10.800 1.7500
   VDW_R0=1.5620 1.3420  1.0010  1.5620 1.4520

  SPIN POLARIZATION SPECIFICATION
   ISPIN = 2
      LORBIT  = 11

IONIC RELAXATION
   NSW    = 500
   NBLOCK  = 1
   IBRION  = 1
   ISIF = 3
   POTIM = 0.5
   LWAVE = .F.
   LCHARG  = .F.
      LDAU= .TRUE.
      LDAUTYPE=2
      LDAUU=3.8 0 0 4.3 0
      LDAUJ=0 0 0 0 0
      LDAUL=2 -1 -1 2 -1
      #LMAXMIX=4

  ELECTRON OCCUPATION
   ISMEAR = 0
   ALGO = fast

  DOS Stuff
   SIGMA = 0.05    check T*S to make sure

优化还在进行，有继续往下降的趋势，这是正常的吗?还是说能量已经收敛，力还没收敛？

作者
Author: 瓜大小瓜子 时间: 2019-5-3 20:25
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.139127979883E+03 0.13711E-03 -0.17145E-03 35920 0.176E-01 0.457E-02
RMM: 2 -0.139128147874E+03 -0.16799E-03 -0.18851E-04 19651 0.674E-02 0.235E-01
RMM: 3 -0.139127967051E+03 0.18082E-03 -0.10005E-04 18342 0.513E-02 0.204E-02
RMM: 4 -0.139127971353E+03 -0.43025E-05 -0.14961E-05 16418 0.191E-02
288 F= -.14189069E+03 E0= -.14188990E+03  d E =-.762699E-05  mag=    4.0000
BRION: g(F)=  0.409E-04 g(S)=  0.907E-06 retain N=  1 mean eig= 0.15
eig: 0.146
bond charge predicted
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.139128099011E+03 -0.13196E-03 -0.55679E-03 33792 0.277E-01 0.259E-02
RMM: 2 -0.139128121947E+03 -0.22936E-04 -0.14827E-04 18393 0.562E-02 0.745E-02
RMM: 3 -0.139128107822E+03 0.14125E-04 -0.32733E-05 17059 0.353E-02
289 F= -.14189072E+03 E0= -.14188993E+03  d E =-.258207E-04  mag=    4.0000
BRION: g(F)=  0.622E-04 g(S)=  0.101E-04 retain N=  1 mean eig= 0.23
eig: 0.234
bond charge predicted
优化已经进行到289步

作者
Author: 啦啦黑还黑 时间: 2019-5-3 23:19

瓜大小瓜子发表于 2019-5-3 20:25
N E dE d eps ncg rms rms(c)
DAV: 1 ...

看这个帖子：
http://bbs.keinsci.com/thread-13065-1-1.html

EDIFF = 1e-4，IVDW=1，都会使得几何优化收敛变慢。

作者
Author: 瓜大小瓜子 时间: 2019-5-4 15:42

啦啦黑还黑发表于 2019-5-3 23:19
看这个帖子：
http://bbs.keinsci.com/thread-13065-1-1.html

好的，谢谢刘老师

作者
Author: 卡开发发 时间: 2019-5-4 22:25
额外补充一点，VASP的CG算法带线搜索确定步长，而RMM-DIIS算法没有，因此对于结构不太好的情况CG相对而言要更加稳定一些。

作者
Author: 瓜大小瓜子 时间: 2019-5-5 19:28

卡开发发发表于 2019-5-4 22:25
额外补充一点，VASP的CG算法带线搜索确定步长，而RMM-DIIS算法没有，因此对于结构不太好的情况CG相对而言要 ...

好的感谢前辈指导

作者
Author: 飞翔的猪 时间: 2019-5-6 08:36
你好，我想请问一下，bulk优化好了之后，得到的结构文件 CONTCAR怎么变成cif？因为后面要切111晶面

作者
Author: 瓜大小瓜子 时间: 2019-5-8 10:54

飞翔的猪发表于 2019-5-6 08:36
你好，我想请问一下，bulk优化好了之后，得到的结构文件 CONTCAR怎么变成cif？因为后面要切111晶面

拖到vesta导出cif

欢迎光临计算化学公社 (http://ccc.keinsci.com/)