计算化学公社

标题: 请问在使用vasp进行MD过程中出现的错误问题 [打印本页]

作者
Author: 白泽666 时间: 2023-1-9 12:04
标题: 请问在使用vasp进行MD过程中出现的错误问题
各位老师同学们我在使用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；

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

作者
Author: 白泽666 时间: 2023-1-9 15:22
好的，谢谢您的提醒；我的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；

作者
Author: 白泽666 时间: 2023-1-9 17:23
本帖最后由白泽666 于 2023-1-9 22:04 编辑

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

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

作者
Author: 白泽666 时间: 2023-2-10 16:43
另外我的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效率区别应该不是很大，是否还需要更改吗

作者
Author: 876449830 时间: 2023-2-10 16:51

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

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

作者
Author: 白泽666 时间: 2023-2-10 17:55

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

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

作者
Author: 876449830 时间: 2023-2-11 11:38

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

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

作者
Author: 白泽666 时间: 2023-2-11 15:09

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

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

作者
Author: honghong 时间: 2023-11-27 20:29
你好，请问解决了么，我也遇到了相同的问题，但是我的晶胞已经优化到了10e-8了，但是还是出现相同的问题。

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