求助大佬CL-NEB計(jì)算問(wèn)題
計(jì)算一步收斂,受力確為0。但這也只是針對(duì)某些體系,有的體系能夠正常計(jì)算,并得出結(jié)果。檢查過(guò)中間的結(jié)構(gòu),沒(méi)得問(wèn)題。
主要參數(shù)設(shè)置:
NSW = 200 # number of steps for ionic upd
IBRION = 3 # ionic relaxation 0=MD 1=quasi-New 2=CG
POTIM = 0
ICHAIN = 0
IMAGES = 5
SPRING = -5
LCLIMB = T
LTANGENT = T
LTANGENTOLD = F
LDNEB = F
NEBCELL = F
IOPT = 1
MAXMOVE = 0.2
ILBFGSMEM = 20
LGLOBAL = .TRUE.
INVCURV = 0.01
LLINEOPT = .FALSE
結(jié)果1:
1 F= -.13251778E+03 E0= -.13248957E+03 d E =-.132518E+03
2 F= -.13251783E+03 E0= -.13248962E+03 d E =-.433186E-04
3 F= -.13251784E+03 E0= -.13248963E+03 d E =-.123478E-04
4 F= -.13251785E+03 E0= -.13248964E+03 d E =-.836182E-05
5 F= -.13251785E+03 E0= -.13248964E+03 d E =-.579447E-05
結(jié)果2:
0 0.000000 -128.224300 0.000000
1 0.000000 -132.489600 -4.265300
2 0.000000 -132.340900 -4.116600
3 0.000000 -132.208900 -3.984600
4 0.000000 -132.154000 -3.929700
5 0.000000 -132.165600 -3.941300
6 0.000000 -127.784600 0.439700 :
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
全部為0肯定是錯(cuò)的,除非你的結(jié)構(gòu)原子全Fixed了,
簡(jiǎn)單,全部重新新建重新提交算。。。