phonopy声子群速数据绘图

经验分享

公众号前期推送介绍有关声子群速计算与数据导出内容。生成文件包括phonon_velocity.dat 和 pband_gv.dat 。本篇文章中主要依据该数据文件进行数据处理和绘图。

依据单层BN结构，计算得到无虚频声子谱图如下，扩包倍数为4×4×1。

使用phonopy计算得到声速并借助vaspkit的Sort Phonon Band Structure for Phononpy功能导出数据得到数据文件。

打开 pband_gv.dat 文件后，将omega列设置为x，

将gv列数据依据导出数据空白位置隔断，并复制为单独一列，这里只复制了前三段数据（声学支数据）

绘图得到如下，彩色部分为三条声学支声子群速。横坐标单位与声子谱纵坐标一致。

打开phonon_velocity.dat文件，并直接选中绘图，结果如下。

曲线为每条声子谱线声速，横坐标与声子谱波什路径横坐标一致。

愿有所成

IEchoQ

引喻失义   妄自菲薄