經(jīng)驗分享

公眾號前期推送介紹有關聲子群速計算與數(shù)據(jù)導出內容。生成文件包括phonon_velocity.dat 和 pband_gv.dat 。本篇文章中主要依據(jù)該數(shù)據(jù)文件進行數(shù)據(jù)處理和繪圖。

依據(jù)單層BN結構，計算得到無虛頻聲子譜圖如下，擴包倍數(shù)為4×4×1。

使用phonopy計算得到聲速并借助vaspkit的Sort Phonon Band Structure for Phononpy功能導出數(shù)據(jù)得到數(shù)據(jù)文件。

打開?pband_gv.dat?文件后，將omega列設置為x，

將gv列數(shù)據(jù)依據(jù)導出數(shù)據(jù)空白位置隔斷，并復制為單獨一列，這里只復制了前三段數(shù)據(jù)（聲學支數(shù)據(jù)）

繪圖得到如下，彩色部分為三條聲學支聲子群速。橫坐標單位與聲子譜縱坐標一致。

打開phonon_velocity.dat文件，并直接選中繪圖，結果如下。

曲線為每條聲子譜線聲速，橫坐標與聲子譜波什路徑橫坐標一致。

愿有所成

IEchoQ

引喻失義? ?妄自菲薄