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

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

瓜大小瓜子

      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步

啦啦黑还黑

瓜大小瓜子 发表于 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，都会使得几何优化收敛变慢。

瓜大小瓜子

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

好的，谢谢刘老师

卡开发发

额外补充一点，VASP的CG算法带线搜索确定步长，而RMM-DIIS算法没有，因此对于结构不太好的情况CG相对而言要更加稳定一些。

瓜大小瓜子

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

好的 感谢前辈指导

飞翔的猪

你好，我想请问一下，bulk优化好了之后，得到的结构文件 CONTCAR怎么变成cif？因为后面要切111晶面

瓜大小瓜子

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

拖到vesta导出cif