常见问题

过渡态不收敛

过渡态收敛要比几何优化困难的多，经常会碰到各种问题导致过渡态优化不收敛，原因非常多，再另外的博文中具体分析，所有用来解决几何优化不收敛的方法同样使用与过渡态。

（1）用nebmake.pl产生的初始点是线性插值的，很可能不合理。用idpp.py脚本非线型插点，可以避免不合理的初始结构。

（2）CI-NEB插点的数目不合理，（NEB是插点越多越精确，但是CI-NEB不是越多越精确！！！因为能量最高的点会自动climb到鞍点，如果CI-NEB计算正确结束，用1个点和用10个点找出来的过渡态是一样的。）

用dist.pl ini/CONTCAR fin/CONTCAR的返回值除以0.8大概就是要插的点的数目。

（3）先用NEB或者CINEB进行粗收敛：

EDIFF = 1E-5；EDIFFG = -0.5；IBRION = 3；POTIM = 0 ；IOPT = 7；ICHAIN = 0；LCLIMB = .TRUE.

再用DIMER或者更高精度NEB精收敛

（4）再用DIMER精收敛的方法：

在NEB文件夹依次运行nebresults.pl和neb2dim.pl产生新的dim文件夹下进行Dimer运算，修改INCAR：

EDIFF = 1E-7；EDIFFG = -0.03；IBRION = 3；POTIM = 0；ICHAIN = 2；IOPT = 2

注意：这里新生成的POSCAR使用的原NEB文件夹里的两个能量最高的结构的中点，如果粗收敛用的NEB是没问题的，但是粗收敛用的CINEB则需要把原CI-NEB文件夹里能量最高的点所对应的结构复制过来， 比如:cp ../03/CONTCAR ./POSCAR

MODECAR可以用neb2dim.pl自动产生的， 也可以先给新POSCAR算个频率，看一下振动方向，把自己认为最可能是过渡态的虚频的振动模复制出来做成MODECAR的格式。（该振动一般是频率最高的一个虚频） 也可以直接把所有原子的振动模都写成0，只把最关键的一个或几个原子(发生转移的)的模手动设置。（该方法只适用于发生转移的原子明显，且其他原子几乎不动的情况）

（5）再用CI-NEB精收敛的方法：

运行vfin.pl保存粗收敛的计算结果，再把所有CONTCAR复制成POSCAR重新计算，修改INCAR：

EDIFF = 1E-7；EDIFFG = -0.03；IBRION = 3；POTIM = 0 ；IOPT = 1或2；ICHAIN = 0；LCLIMB = .TRUE.

或者用VASP自带的DIIS收敛方法（此法往往比vtst自带的收敛算法更易收敛）：

EDIFF = 1E-7；EDIFFG = -0.03；IBRION = 1；POTIM = 0.1 ；IOPT = 0；ICHAIN = 0；LCLIMB = .TRUE.