VASP+NEB

第一步：
    初末态结构优化建立两个文件夹 ini, fin，在每个文件夹中放vasp计算必备的四个文件（INCAR，POSCAR，KPOINTS，POTCAR），其中的两个POSCAR对应未优化的初末态。确保两个文件夹里面除POSCAR外，其他文件完全一样。对于两个POSCAR，注意每一行的原子一一对应。若有固定位置的原子（比方做表面计算，要固定底部1-2层原子）。
第二步：
  用dist.pl检查两个优化后结构(两个CONTCAR)的相似程度（每个对应原子的初末态距离的平方和，再开根号）。 dist.pl ini/CONTCAR fin/CONTCAR
第三步：
     插点的数目取决于前面dist.pl的返回值，一般插点数目可取（dist.pl返回值/0.8）。 具体插点的命令： nebmake.pl ini/CONTCAR fin/CONTCAR N 最后的N表示插点数目。插入点的算法为线性插值，详情请进前面给的vtst脚本链接。我这里为简单插一个点，执行完命令后出现文件夹00--07。 可用命令nebavoid.pl 1,确保中间插入的点每一个原子间距都大于1A。该命令的参数1表示最小允许间距，可取小数。

System = NEB

PREC = Accurate

#ISYM =0

ENCUT = 650

EDIFF  = 1e-4

EDIFFG = -0.001

KSPACING= 0.5

ISIF   = 2

NSW    = 500

ISMEAR = 0;

SIGMA = 0.05

LCHARG = FALSE

LWAVE  = FALSE

ISTART = 0

IVDW=10

#NPAR = 3

#NEB

IBRION = 3

POTIM = 0

IOPT =   1

ICHAIN = 0

LCLIMB = .TRUE.

SPRING =-5

IMAGES =6

第四步：

    后处理：nebef.pl， nebbarrier.pl 生成neb.dat。

    计算完成使用命令nebresult.pl.完成了nebbarrier.pl, nebspline.pl, nebef.pl, nebmovie.pl, nebjmovie.pl, nebconverge.pl还有对各文件夹中的OUTCAR打包压缩。生成了很多文件。其中mep.eps是以dist.pl距离为横坐标，能量为纵坐标画出的能势垒图.计算过渡态，最好是把vtstool的工具编译到vasp中去，好处是仍然可以只用vasp中的NEB,也可以用Cl-NEB，更重要的原因是可以用vtstool中的很多脚本工具来跟踪过渡态计算过程。
计算过渡态先要摆正心态，不急于下手。步骤如下：
（1）做模型，初态IS和终态FS，分别结构优化到基态；
（2）线形插入images:  nebmake.pl  POSCAR.IS POSCAR.FS  N
N为image个数。
（3）nebmovie.pl,生成movie.xyz。用Xcrysden --xyz movie.xyz 反复观看动画，仔细检查过程的合理性。这里要提醒，POSCAR.IS 和POSCAR.FS中原子坐标列表的顺序必须对应。
（4）写INCAR,选IOPT。注意，最好忘记vasp自带的NEB，而全部改用包含vtstool的vasp. IBRION=3,POTIM=0关闭vasp自带的NEB功能。
（5）过渡态计算第一个离子步最耗时，也最容易出问题，也是模型设计合理性检验的首要环节。所以可以选小一些的ENCUT，可以不用考虑自旋(ISPIN=1)，也不用考虑DFT+U。而且用最快最粗糙的算法（IOPT=3,其他默认）。
（6）带vtstool的vasp-ClNEB(NEB)过渡态计算ICHAIN=0作为入口，这个也是默认的。LCLIMB=TRUE也是默认的。如果不要climb image,可以设置LCLIMB = False.
（7）收敛判据EDIFFG<0。过渡态计算要以力为收敛判据，而不是能量。一般EDIFFG=-0.05就可以接受，-0.02或者-0.01更好。但是作为开始的过渡态计算，可以设置很宽的收敛条件，如EDIFFG=-1.
（8）初步过渡态收敛后,修改INCAR中的优化器（IOPT），并修改相应参数（参考vtstool官方论坛），EDIFFG改小（如-0.05），然后运行vfin.pl，这个脚本自动帮你准备在原来的基础上继续运行新的过渡态计算（完成cp CONTCAR POSCAR, 保留电荷密度和波函数的操作）。<copyright@ muchong>