使用cp2k计算导体声子时，声子谱总是有虚频，qe计算声子谱稳定

honor

最近在学习使用cp2k进行声子相关计算，尝试了几次，计算得出的声子谱总是有虚频。我先优化了MgB2的原胞，优化完后创建了计算能量和力的输入文件8-1.inp,然后使用phonopy --cp2k -c 8-1.inp -d --dim="3 3 2"对输入文件进行扩胞和位移，生成了两个超胞，然后计算超胞的力，后边我使用phonopy --cp2k -f {001..002}/*.xyz生成力场数文件，然后我就使用phonopy --cp2k -c 8-1.inp -p --dim="3 3 2" --pa="auto" --band="1/2 1/2 1/2 0 0 0 1/2 0 1/2" -v，尝试了几次总是有虚频，想请教一下各位老师我的问题出现在哪里，下面是我的输入文件，期待老师的解答