求助大佬CL-NEB计算问题
计算一步收敛,受力确为0。但这也只是针对某些体系,有的体系能够正常计算,并得出结果。检查过中间的结构,没得问题。
主要参数设置:
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
结果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
结果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 : 返回小木虫查看更多
全部为0肯定是错的,除非你的结构原子全Fixed了,
简单,全部重新新建重新提交算。。。