VASP入門到精通：[20]計算聲子譜2?

之前介紹到怎麼去計算聲子譜，計算完成後，可以去分析聲子譜

工具/原料

VASP phonopy

方法/步驟

之前小編給大家介紹了怎麼樣去計算聲子譜，在14節有介紹，沒有看到的小夥伴可以去看一下

注意INCAR裡不要多加入一些沒用的參數設置，這樣容易出錯，只要必須的參數設置就可以了，如小編上一節介紹的

計算完成後結果如圖所示：下圖是OSZICAR數據：

下圖是OUTCAR最後的數據

注意要有 ELASTIC MODULI IONIC CONTR (kBar)

這一行的數據

vasprun.xml裡面要有力學Hessian矩陣才能繪製聲子譜，步驟如下：

運行命令：

phonopy --fc vasprun.xml

然後會看到生成力學文件：FORCE_CONSTRAINS

在band.conf文件中加入一行FORCE_CONSTRAINS=READ

注意修改裡面的原包基矢和高對稱點

然後運行phonopy --dim="m n l" -c POSCAR-unitcell band.conf

接下來繪製，bandplot --gnuplot> PHONO.dat

這個PHONO.dat文件便是聲子譜，把它放到origin裡面去處理吧，大家就能得到想要的聲子譜了！！！

有問題的歡迎來問哦，求評價，如果覺得好，給個贊吧