??圖片未上傳，需要帶圖版可加wx：lrm______或郵箱liruimin0187@link.tyut.edu.cn?

本教程為使用VASP 5.4.4軟件計算擴散能壘教程，請?zhí)崆鞍惭bVESTA（用于建模及轉(zhuǎn)換格式）、materials studio（較大，若已安裝VESTA可不安裝，用于建模）。

以計算LiF中Ca離子擴散能壘為例，大致步驟為，下載LiF晶體結(jié)構(gòu) – 擴胞 – 優(yōu)化 – 構(gòu)建Ca離子在LiF中初末態(tài)構(gòu)型 – 優(yōu)化 – 生成中間態(tài)構(gòu)型 – 優(yōu)化 – 處理數(shù)據(jù)。

一、獲得初末態(tài)離子遷移模型

l? ??同計算吸附能部分從materials project下載LiF (mp-1138) cif文件，擴胞2*2*2，生成VASP計算文件，結(jié)構(gòu)優(yōu)化后CONTCAR去MS添加Ca在LiF間隙位點，獲得初末態(tài)離子遷移模型，左圖為MS顯示圖，右圖為VESTA圖，VESTA圖更便于觀察配位情況。

二、生成中間態(tài)遷移路徑模型

l? ??同計算吸附能部分生成VASP計算文件，更改INCAR中ISIF=2再優(yōu)化結(jié)構(gòu)，將計算初末態(tài)得到的CONTCAR更改為新的POSCAR，分別命名文件夾為start和end，用neb_path.py代碼生成9個path的遷移路徑

l? 此時文件夾下出現(xiàn)11個遷移路徑的VASP計算文件夾

l? ??生成的11個文件夾中，0和10為start和end的POSCAR，僅需計算1-9即可，使用start和end相同的INCAR、POTCAR和KPOINT，逐個計算。

三、生成遷移路徑能壘圖

l? 將服務(wù)器計算得到的OUTCAR和CONTCAR保存到各自文件夾下，使用neb_barrier.py代碼自動生成Ca遷移路徑圖，同使用neb_path.py的命令，先cd到當(dāng)前文件夾路徑下，再執(zhí)行python neb_barrier.py命令，即可自動生成遷移路徑能壘pdf文件

?

?

?