计算一步收敛，受力确为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 ： 返回小木虫查看更多