求助：VASP大体系SCF不收敛

diffcwmol

各位同仁，
请教一下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

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

七尺贱

虽然我不是很懂vasp，但我觉得你至少应该给个结构的图吧，这样大家都不知道你什么体系，结构合不合理什么的

后到咸阳斯蒂庞

我倒是用过一种方法，就是先用低k点先进行一次静态计算，结束后在INCAR中设置读取WAVECAR后用高k点再进行一次计算 。这种方法有可能会使其收敛。