1. 首先得到力FORCE_SETS，操作为： phonopy -f vasprun.xml-00*

2. 得到态密度和带结构

a. 得到态密度所需文件及操作

 所需文件： mesh.conf FORCE_SETS POSCAR（扩胞前的Unitcell）

 其中mesh.conf内容如下：

ATOM_NAME = Sc

DIM = 3 3 2

MP = 8 8 8

  操作为： phonopy -p --factor=521.471 mesh.conf

b. 得到带结构所需文件及操作

 所需文件：band.conf FORCE_SETS POSCAR（扩胞前的Unitcell）  

 其中band.conf内容如下：

ATOM_NAME = Sc

DIM = 3 3 2

BAND=  0.000  0.000  0.000    0.000  0.000  0.500     -0.333  0.667  0.500    -0.333  0.667  0.000   0.000  0.000  0.000  0.000  0.500  0.000    0.000  0.500  0.500   -0.333  0.667  0.500

操作：1） phonopy -p --factor=521.471 band.conf  得到band.yaml

     2)  bandplot --gnuplot >PhononBAND.dat  将band.yaml转化为可用orgin绘图的数据文件