Improved Dimer Method计算处理求助

user0203

各位老师同学好！
我现在在算过渡态，之前用CI-NEB 方法跑出来不收敛，于是继续跑了IDM 计算，但是网络上很少有基于该算法的结果分析和后处理，能否请各位友友给点意见和帮助？
目前关于我的体系下，IDM的INCAR 参数这样设置，算完后结果文件里没有出现收敛标志，但是查看OUTCAR 的力应该是收敛了，这种情况下的CONTCAR结构可以作为过渡态吗？

卡开发发

不好说，没停止的原因可能是你的curvature非负，你可以看这里：

The dimer is rotated such that its axis is parallel with the direction of the maximal negative curvature.

user0203

卡开发发 发表于 2023-12-30 03:23
不好说，没停止的原因可能是你的curvature非负，你可以看这里：

你好，IDM 运算结束后怎么查看结果如Torque呢？网上关于这部分的描述还挺少的。用什么后处理方法呢？

卡开发发

user0203 发表于 2024-1-5 21:36
你好，IDM 运算结束后怎么查看结果如Torque呢？网上关于这部分的描述还挺少的。用什么后处理方法呢？

没有Torque的信息，IDM的信息可以在OUTCAR中查看，类似于这样的字段：

  
DIMER METHOD: (5)
  -----------------------------------------------------------
    max. gradient (eV/A):   0.07770164
    dimer rotated by (deg.):   21.5967
    curvature along the dimer direction:   -5.8685

user0203

卡开发发 发表于 2024-1-5 23:19
没有Torque的信息，IDM的信息可以在OUTCAR中查看，类似于这样的字段：

老师好，我去查了下OUTCAR 文件，最后四次移动的值均为负，具体变化情况可以看下我画这张图。另外贴出了OUTCAR的最后一页，还请指点

卡开发发

user0203 发表于 2024-1-6 21:52
老师好，我去查了下OUTCAR 文件，最后四次移动的值均为负，具体变化情况可以看下我画这张图。另外贴出了O ...

那就不是curvature的问题，你看你最终的力是多少，你上面给出的力是个3*N的矢量F=(Fx,Fy,Fz)，但对VASP而言实际力收敛的标准是max(|F|)=max(sqrt(Fx^2+Fy^2+Fz^2))。。如上面字段中DIMER METHOD:这里往下一点：

max.  gradient (eV/A):

或者如果你装了VASP那个VTST插件的话，可以用grep FORCES: OUTCAR来看力变化的情况，这是比较简单的做法。当然其他也可以自己写一些脚本来处理。

卡开发发

user0203 发表于 2024-1-6 21:52
老师好，我去查了下OUTCAR 文件，最后四次移动的值均为负，具体变化情况可以看下我画这张图。另外贴出了O ...

你看你顶楼OUTCAR中力的输出倒数第十个数据，分量都大于>0.05了（那个-0.051805），另外你的EDIFFG设置的是5e-3。