VASP+NEB

2024-03-06 约 646 字预计阅读 2 分钟

前言

简述VASP + VTST + NEB计算缺陷迁移动力学、反应路径、过渡态等。

前期准备

¹²³⁴

计算缺陷迁移（反应路径、过渡态等类似）的始态和终态，relax计算，ini和fin两个文件夹。
运行dist.pl ini/CONTCAR fin/CONTCAR，结果/0.8为需要插入的结构数，下文结构数以5为例。
运行nebmake.pl ini/CONTCAR fin/CONTCAR 5生成00~06这7个文件夹。接着，运行nebavoid.pl 1，确保没有原子间距小于1 Å的结构。
运行cp ini/OUTCAR 00/ && cp fin/OUTCAR 06/将ini和fin的OUTCAR复制到对应的00和06文件夹内。
修改nebmovie.pl文件，

1
2
3


#if($xyzflag==0){
#    unlink "movie.xyz";
#}

运行nebmovie.pl，可得到movie和movie.xyz文件。

movie.xyz可以通过ase gui movie.xyz查看，也可以通过Jmol查看，我倾向后者。

在relax的INCAR基础上修改，

1
2
3
4
5
6
7
8
9


EDIFF = 1e-5
EDIFFG = -0.01
IBRION = 3
POTIM = 0
IOPT = 3
ICHAIN = 0
LCLIMB = .T.
IMAGES = 5
SPRING = -5

注意

job.sh中的vasp_std使用VTST编译的版本。

运行nebef.pl。

运行nebbarrier.pl，生成neb.dat用于作图。