虎扑体育网

前面我们介绍了使用Phonopy计算声子谱、非解析项校正以及热学性质，今天我们简单介绍下在声子谱计算的基础上，如何通过Phonopy计算声子群速度并利用vaspkit导出相关数据。

使用Phonopy计算声子群速度时，q点处的群速度主要由标签GROUP_VELOCITY控制，群速度被写入yaml文件中，对应于笛卡尔坐标的运行模式。物理单位取决于输入文件的物理单位和频率转换因子，但是如果仅仅使用VASP和默认设置（例如，声子频率为THz），那么物理单位将是Angstrom THz。

操作步骤

(1) 首先对结构高精度优化后扩胞，计算得到没有虚频的声子谱；

(2) 对band.conf文件进行修改，并添加群速度计算相关参数GROUP_VELOCITY = .TRUE.和GV_DELTA_Q = 0.01，程序需要读取本征值和本征矢，所以必须设置 EIGENVECTORS =.TRUE.，此外为了防止vaspkit在导出数据时报错，需要删除BAND_LABELS高对称路径名称，修改后band.conf文件内容如下；

(3) 对声子谱重新进行计算，打开phonopy.yaml文件可以看到，phonopy在计算声子谱的过程中考虑到了群速度；

(4) 打开band.yaml文件中也保存了有关群速度的相关信息；

(5) 随后进行数据提取，以vaspkit1.12版本为例，打开vaspkit，输入73，打开VASP2other Interface 子功能，选择739) Sort Phonon Band Structure for Phononpy，运行完成即可获得group velocity相关数据以及各原子的声子谱数据；