发现用phonopy算声子谱方法（不用重复提交POSCAR-0*算scf）,而是算SPOSCAR，INCAR里的参数IBRION要设为8（8为计算Hessian矩阵，在vasprun.xml里）,运行phonopy --fc vasprun.xml得到力常数文件FORCE_CONSTRAINTS，然后运行phonopy --dim=“m n l”-c POSCAR-unitcell --readfc -p band.conf可以在线画出声子谱图，最后执行bandplot --gnuplot band.yaml > band.dat得到band.dat,导入origin里画图。

详细的参考https://wenku.baidu.com/view/b8a3ae7b27284b73f2425094.html