本文主要是介绍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
引喻失义 妄自菲薄
这篇关于phonopy声子群速数据绘图的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!