计算化学公社

标题: 求助：VASP大体系SCF不收敛 [打印本页]

作者
Author: diffcwmol 时间: 2022-7-20 10:09
标题: 求助：VASP大体系SCF不收敛
各位同仁，
请教一下VASP里SCF收敛问题。我计算的体系是一个有250个原子，共2000多个电子，x和y方向有真空层，z方向周期性。材料是半导体。ISMEAR =  0 ，SIGMA  =  0.05。

在计算的时候，无论怎么调参数（比如AMIX BMIX），总是出现下边类似的收敛过程。比如下面到21步开始能量就飘了。
尝试过刘博士这个帖子里建议的所有方法（“[VASP] 解决VASP计算SCF/几何优化/过渡态不收敛的方法总结 (待更新)”），怎么着都不收敛。
调试过用不同数量的core，改变NCORE等，都不收敛。

请问各位同仁对此有没有好的解决方案，或者建设性意见？
感谢！

################
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1    0.189539117317E+05 0.18954E+05 -0.94076E+05 19152 0.120E+03
DAV: 2    0.103346226760E+04 -0.17920E+05 -0.17374E+05 21760 0.273E+02
DAV: 3 -0.176724054343E+04 -0.28007E+04 -0.27823E+04 17984 0.995E+01
DAV: 4 -0.187769207482E+04 -0.11045E+03 -0.10879E+03 17752 0.214E+01
DAV: 5 -0.188236618041E+04 -0.46741E+01 -0.46249E+01 20656 0.388E+00 0.182E+02
DAV: 6 -0.175929881483E+04 0.12307E+03 -0.37505E+02 17992 0.156E+01 0.113E+02
DAV: 7 -0.174601209518E+04 0.13287E+02 -0.70641E+01 18968 0.660E+00 0.523E+01
DAV: 8 -0.174360618640E+04 0.24059E+01 -0.25016E+01 19152 0.391E+00 0.188E+01
DAV: 9 -0.174381543849E+04 -0.20925E+00 -0.67185E+00 19272 0.235E+00 0.120E+01
DAV:  10 -0.174338960176E+04 0.42584E+00 -0.13648E+00 18872 0.115E+00 0.644E+00
DAV:  11 -0.174305959989E+04 0.33000E+00 -0.74431E-01 18168 0.775E-01 0.519E+00
DAV:  12 -0.174298027340E+04 0.79326E-01 -0.60390E-01 18768 0.715E-01 0.258E+00
DAV:  13 -0.174294694024E+04 0.33333E-01 -0.16614E-01 18872 0.383E-01 0.920E-01
DAV:  14 -0.174294261764E+04 0.43226E-02 -0.23781E-02 19616 0.159E-01 0.597E-01
DAV:  15 -0.174293994792E+04 0.26697E-02 -0.64633E-03 19304 0.869E-02 0.469E-01
DAV:  16 -0.174293788183E+04 0.20661E-02 -0.39185E-03 19120 0.696E-02 0.350E-01
DAV:  17 -0.174293608018E+04 0.18016E-02 -0.19277E-03 18672 0.481E-02 0.220E-01
DAV:  18 -0.174293570252E+04 0.37767E-03 -0.73810E-04 19056 0.253E-02 0.167E-01
DAV:  19 -0.174293551315E+04 0.18936E-03 -0.43998E-04 18488 0.203E-02 0.110E-01
DAV:  20 -0.174293558518E+04 -0.72025E-04 -0.78437E-04 17664 0.202E-02 0.615E-02
DAV:  21 -0.174293620898E+04 -0.62380E-03 -0.41184E-04 17312 0.133E-02 0.593E-02
DAV:  22 -0.246107449121E+04 -0.71814E+03 -0.22064E+02 19296 0.118E-01 0.215E+01
DAV:  23 -0.508594324239E+06 -0.50613E+06 -0.34303E+06 17800 0.784E+00 0.117E+02
DAV:  24 -0.351035025664E+07 -0.30018E+07 -0.26125E+07 17296 0.353E+01 0.297E+02
DAV:  25 -0.113846776899E+07 0.23719E+07 -0.30747E+04 16800 0.134E+02 0.292E+02

###########################

作者
Author: 七尺贱 时间: 2022-7-21 10:05
虽然我不是很懂vasp，但我觉得你至少应该给个结构的图吧，这样大家都不知道你什么体系，结构合不合理什么的

作者
Author: 后到咸阳斯蒂庞 时间: 2024-10-26 16:03
我倒是用过一种方法，就是先用低k点先进行一次静态计算，结束后在INCAR中设置读取WAVECAR后用高k点再进行一次计算。这种方法有可能会使其收敛。

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