1、先确定初末结构，初末结构不仅要做结构优化，还要做频率计算，确定是否为稳定结构
2、将初末结构复制到过渡态计算的文件夹下，分别命名为fi in文件夹，然后用nebmake.pl脚本创建插入结构，插入结构数量为dist.pl脚本返回值除以0.8，具体情况根据自己服务器的计算能力而定。
3、过渡态计算参数，我使用的是vtst脚本，所以计算参数额外有一些；
IBRION=3;POTIM=0;IOPT=1;ICHAIN=0;LCLIMB=.T.;IMGES=3;SPRING=-5
此外LWAVE=.T.建议输出波函数，因为后面可能会用到。ISIF=2等于3会莫名停掉。
4、计算过程中或者完成后用nebefs.pl查看受力情况和能垒，然后确定过渡态，计算频率，如果有大于一个虚频就要进行消虚频；做法如下
频率计算完成后；使用grep f/i OUTCAR 搜索文件，然后确定要消几号虚频；使用grep -As “x f/i” OUTCAR 输出虚频的振动模式，如下图
前面x y z 表示原子位置，dx dy dz后面表示振动方向。将前面的X=x+0.1dx 换算成新的坐标 然后复制到过渡态的POSCAR中，其他的插入点用CONTCAR替换POSCAR。然后重新计算过渡态，计算频率，消虚频，直到有且只有一个虚频。
5、参考看这里