搜索
VASP计算频率设置IBRION=5计算过渡态，最好是把vtstool的工具编译到vasp中去，好处是仍然可以只用vasp中的NEB,也可以用Cl-NEB，更重要的原因是可以用vtstool中的很多脚本工具来跟踪过渡态计算过程计算过程。计算过程计算过渡态先要摆正心态，不急于下手。步骤如下：（1）做模型，初态IS和终态FS，分别结构优化到基态；（2）线形插入images:nebmake.plPOSCAR.ISPOSCAR.FSNN为image个数。（3）nebmovie.pl,生成movie.xyz。用Xcrysden--xyzmovie.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也是默认的。如果不要climbimage,可以设置LCLIMB=False.（7）收敛判据EDIFFG0。过渡态计算要以力为收敛判据，而不是能量。一般EDIFFG=-0.05就可以接受，-0.02或者-0.01更好。但是作为开始的过渡态计算，可以设置很宽的收敛条件，如EDIFFG=-1.（8）初步过渡态收敛后,修改INCAR中的优化器（IOPT），并修改相应参数（参考vtstool官方论坛），EDIFFG改小（如-0.05），然后运行vfin.pl，这个脚本自动帮你准备在原来的基础上继续运行新的过渡态计算（完成cpCONTCARPOSCAR,保留电荷密度和波函数的操作）。