请问在使用vasp进行MD过程中出现的错误问题

白泽666

各位老师同学们我在使用vasp进行分子动力学过程中出现如下报错，请问是什么原因导致的?如何解决这个问题呢？
报错信息：

4829 T=   313. E= -.62209792E+03 F= -.63011612E+03 E0= -.63008082E+03  EK= 0.80182E+01 SP= 0.00E+00 SK= 0.00E+00
Information: wavefunction orthogonal band  816  0.8769
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.150147915858E+08   -0.15014E+08   -0.14925E+08  2136   0.349E+04    0.604E+03
DAV:   2    -0.834123329750E+04    0.15006E+08   -0.13893E+05  2940   0.183E+02    0.260E+03
DAV:   3    -0.343341727973E+06   -0.33500E+06   -0.60295E+06  3036   0.214E+03    0.385E+03
DAV:   4    -0.266033773871E+06    0.77308E+05   -0.22183E+05  2688   0.936E+02    0.286E+03
DAV:   5    -0.733157129703E+05    0.19272E+06   -0.16564E+06  2952   0.613E+02    0.234E+03
DAV:   6    -0.477091142942E+05    0.25607E+05   -0.30425E+05  3264   0.484E+02    0.541E+03
DAV:   7    -0.563786909405E+05   -0.86696E+04   -0.11333E+05  3180   0.365E+02    0.171E+03
DAV:   8    -0.643093937638E+05   -0.79307E+04   -0.85994E+05  2724   0.838E+02    0.209E+03
DAV:   9    -0.352254047937E+05    0.29084E+05   -0.60809E+04  3036   0.221E+02    0.178E+03
DAV:  10    -0.662269911091E+05   -0.31002E+05   -0.13361E+06  3060   0.391E+02    0.242E+03
DAV:  11    -0.595058484349E+05    0.67211E+04   -0.82048E+04  2712   0.304E+02    0.188E+03
DAV:  12    -0.256396464987E+05    0.33866E+05   -0.25073E+05  3156   0.281E+02    0.221E+03
DAV:  13    -0.430661279491E+05   -0.17426E+05   -0.72204E+04  2460   0.265E+02    0.158E+03
DAV:  14    -0.324907882284E+05    0.10575E+05   -0.12555E+06  3132   0.306E+02    0.237E+03
DAV:  15    -0.368068358014E+05   -0.43160E+04   -0.11958E+05  3204   0.269E+02    0.136E+03
DAV:  16    -0.109104438461E+05    0.25896E+05   -0.20241E+05  2652   0.315E+02    0.171E+03
DAV:  17    -0.147340224759E+05   -0.38236E+04   -0.45528E+04  2868   0.155E+02    0.137E+03
DAV:  18    -0.841082704585E+04    0.63232E+04   -0.49998E+04  3024   0.135E+02    0.136E+03
DAV:  19    -0.570180668591E+04    0.27090E+04   -0.29310E+04  2772   0.121E+02    0.122E+03
DAV:  20    -0.463600840482E+04    0.10658E+04   -0.35743E+04  2580   0.124E+02    0.120E+03
DAV:  21    -0.444559384597E+04    0.19041E+03   -0.86845E+04  3012   0.170E+02    0.129E+03
DAV:  22    -0.515107833495E+04   -0.70548E+03   -0.17295E+04  2244   0.903E+01    0.119E+03
DAV:  23    -0.636171570580E+04   -0.12106E+04   -0.15319E+05  3168   0.149E+02    0.120E+03
DAV:  24    -0.424094103800E+04    0.21208E+04   -0.19996E+04  2412   0.104E+02    0.114E+03
DAV:  25    -0.149278553826E+04    0.27482E+04   -0.13008E+04  2388   0.653E+01    0.108E+03
DAV:  26    -0.463237637271E+03    0.10295E+04   -0.20024E+04  2352   0.659E+01    0.105E+03
DAV:  27    -0.881371598117E+03   -0.41813E+03   -0.64917E+03  2088   0.497E+01    0.107E+03
DAV:  28    -0.984056603205E+03   -0.10269E+03   -0.87498E+02  1932   0.245E+01    0.106E+03
DAV:  29    -0.903642099146E+03    0.80415E+02   -0.26610E+02  1932   0.884E+00    0.106E+03
DAV:  30    -0.110092200127E+04   -0.19728E+03   -0.44564E+01  1968   0.448E+00    0.105E+03
DAV:  31    -0.107473377594E+04    0.26188E+02   -0.12802E+01  1968   0.149E+00    0.105E+03
DAV:  32    -0.110727743184E+04   -0.32544E+02   -0.23251E+01  2124   0.931E-01    0.105E+03
DAV:  33    -0.109207107968E+04    0.15206E+02   -0.79370E+00  1944   0.633E-01    0.105E+03
DAV:  34    -0.109494194229E+04   -0.28709E+01   -0.49654E+00  2112   0.956E-01    0.105E+03
DAV:  35    -0.108741676942E+04    0.75252E+01   -0.18568E+00  1824   0.413E-01    0.105E+03
DAV:  36    -0.107582053513E+04    0.11596E+02   -0.46585E+00  2088   0.644E-01    0.105E+03
DAV:  37    -0.106602714685E+04    0.97934E+01   -0.13172E+01  1908   0.532E-01    0.105E+03
DAV:  38    -0.103921081998E+04    0.26816E+02   -0.22555E+01  2028   0.149E+00    0.104E+03
DAV:  39    -0.101316123258E+04    0.26050E+02   -0.35220E+01  2004   0.176E+00    0.103E+03
DAV:  40    -0.965371034414E+03    0.47790E+02   -0.51033E+01  2064   0.389E+00    0.102E+03
DAV:  41    -0.836843042769E+03    0.12853E+03   -0.93374E+01  1992   0.438E+00    0.101E+03
DAV:  42    -0.950924890712E+03   -0.11408E+03   -0.22371E+02  2028   0.113E+01    0.923E+02
DAV:  43    -0.807389641751E+03    0.14354E+03   -0.18700E+02  2040   0.762E+00    0.900E+02
DAV:  44    -0.304534265845E+04   -0.22380E+04   -0.16029E+04  1944   0.269E+01BRMIX: very serious problems
the old and the new charge density differ
old charge density:  1353.01670 new 1370.46654
    0.267E+03
DAV:  45    -0.172075091428E+04    0.13246E+04   -0.43506E+03  2172   0.198E+01BRMIX: very serious problems
the old and the new charge density differ
old charge density:  1370.46654 new 1372.13829
    0.820E+02
DAV:  46    -0.145578860020E+06   -0.14386E+06   -0.98337E+03  2304   0.307E+01BRMIX: very serious problems
the old and the new charge density differ
old charge density:  1372.13829 new 1387.79160；

含光君

贴上输入文件和输出文件才便于大家帮你排查
《在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚》（http://bbs.keinsci.com/thread-25787-1-1.html）

白泽666

好的，谢谢您的提醒；我的INCAR 如下：
PREC    =Med            # default: normal (for vasp 5.x), (medium, high low for compatibility)
ISTART  =0                 #job   : 0-new  1-cont  2-samecut
ISPIN   =1
ICHARG  =1                 
#Electronic Relaxation  
NELM   =  200      #maximum of 60 electronic steps
EDIFF  = 1E-04     #accuracy for electronic minimization
LREAL  = A         #default: F, projection done in recipocal space; Auto or A, projection done in real space
ALGO   = N
#Molecular Dynamics
ISYM   = 0
EDIFFG = -0.02   #stopping-criterion for IOM
NSW    = 5000     #number of steps for IOM
IBRION = 0       #ionic relax: 0-MD 1-quasi-New 2-CG
SMASS  = -3       #stress and relaxation
POTIM  = 1     #time-step for ionic-motion  Note:md
NBLOCK = 1
#NCORE   = 7      #if node=1, NPAR=4; node=2, NCORE=7
NPAR   = 12
ISMEAR = 0
SIGMA  = 0.1
#Write flags
NWRITE = 1
LWAVE  = F          #write WAVECAR
LCHARG = F          #write CHGCAR
#LORBIT = 11
因为是动力学所以K取的为1 1 1；

白泽666

这是我作业的报错日志；希望各位老师同学们能够给予帮助，谢谢

白泽666

一个月之后，经过和相关大佬请教，发现问题是我计算时的精度并不够，EDIFF精度太低，建议重新优化结构，大佬推荐的是1e-6，我结构优化时用的是1e-5；MD用的是1e-4，想请问各位老师同学，这种情况下，我是重新优化结构去跑，还是只在MD中提升精度即可

白泽666

另外我的POSCAR如下：
1.0000000000000000
     5.9674090000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000   20.2228150000000007    0.0000000000000000
     0.0000000000000000    0.0000000000000000   80.8912610000000001
只用单K点的话足够吗？
说推荐用PREC=N，但N和M效率区别应该不是很大，是否还需要更改吗

876449830

白泽666 发表于 2023-2-10 16:38
一个月之后，经过和相关大佬请教，发现问题是我计算时的精度并不够，EDIFF精度太低，建议重新优化结构，大 ...

最好在扩胞前优化好，而且EDIFF只是能量值，最好是EDIFFG的值设为-0.02或者-0.01。得到优化结构后，再去跑MD。还有就是你的c轴为什么这么大？

白泽666

876449830 发表于 2023-2-10 16:51
最好在扩胞前优化好，而且EDIFF只是能量值，最好是EDIFFG的值设为-0.02或者-0.01。得到优化结构后，再去 ...

感谢您的指点，我时间较为紧张，所以才请教，我只设置1-E5的话对结构优化处理是不是已经足够了；
我的C轴这么大是因为我构建的晶界体系，两个GB之间需要足够长的距离来得到一个收敛的grain boundary energy

876449830

白泽666 发表于 2023-2-10 17:55
感谢您的指点，我时间较为紧张，所以才请教，我只设置1-E5的话对结构优化处理是不是已经足够了；
我的C ...

我们一般设置EDIFF=1E-6，EDIFFG=-0.02

白泽666

876449830 发表于 2023-2-11 11:38
我们一般设置EDIFF=1E-6，EDIFFG=-0.02

好的，谢谢您。为了保险起见我还是重新优化一下结构看看

honghong

你好，请问解决了么，我也遇到了相同的问题，但是我的晶胞已经优化到了10e-8了，但是还是出现相同的问题。