文章介绍了如何使用phonopy进行声子群速计算，并通过vaspkit的SortPhononBandStructureforPhononpy功能处理数据。重点展示了从pband_gv.dat和phonon_velocity.dat文件中提取信息，绘制声子谱和声速曲线的过程，特别关注了声学支的声子群速。
经验分享
公众号前期推送介绍有关声子群速计算与数据导出内容。生成文件包括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文件，并直接选中绘图，结果如下。
曲线为每条声子谱线声速，横坐标与声子谱波什路径横坐标一致。
愿有所成
1930
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
