图片未上传，需要带图版可加wx：lrm______或邮箱liruimin0187@link.tyut.edu.cn 

本教程为使用VASP 5.4.4软件计算扩散能垒教程，请提前安装VESTA（用于建模及转换格式）、materials studio（较大，若已安装VESTA可不安装，用于建模）。

以计算LiF中Ca离子扩散能垒为例，大致步骤为，下载LiF晶体结构 – 扩胞 – 优化 – 构建Ca离子在LiF中初末态构型 – 优化 – 生成中间态构型 – 优化 – 处理数据。

一、获得初末态离子迁移模型

l    同计算吸附能部分从materials project下载LiF (mp-1138) cif文件，扩胞2*2*2，生成VASP计算文件，结构优化后CONTCAR去MS添加Ca在LiF间隙位点，获得初末态离子迁移模型，左图为MS显示图，右图为VESTA图，VESTA图更便于观察配位情况。

二、生成中间态迁移路径模型

l    同计算吸附能部分生成VASP计算文件，更改INCAR中ISIF=2再优化结构，将计算初末态得到的CONTCAR更改为新的POSCAR，分别命名文件夹为start和end，用neb_path.py代码生成9个path的迁移路径

l  此时文件夹下出现11个迁移路径的VASP计算文件夹

l    生成的11个文件夹中，0和10为start和end的POSCAR，仅需计算1-9即可，使用start和end相同的INCAR、POTCAR和KPOINT，逐个计算。

三、生成迁移路径能垒图

l  将服务器计算得到的OUTCAR和CONTCAR保存到各自文件夹下，使用neb_barrier.py代码自动生成Ca迁移路径图，同使用neb_path.py的命令，先cd到当前文件夹路径下，再执行python neb_barrier.py命令，即可自动生成迁移路径能垒pdf文件