优化过程省略不写了，就是选择出合适的赝势，k-点网格，Ecut等。
第一步建立*.scf.in文件， 运行pw.x <*.scf.in>  *.scf.out生成后面计算所需的波函数；
      此步计算后会生成一个*.scf.out文件。
第二步建立*.ph.in文件， 运行ph.x <*.ph.in>  *.ph.out生成选定q网格点下的动力学矩阵，运算结束会生成 一个*.scf.out文件和 一串*.dyn.0-n文件，*.dyn.0内容是所选q网格会有多少q点，给出了具体坐标，后面的*.dyn.1-n文件分别是n个q点的具体信息文件。
第三步 建立*.q2r.in文件， 运行q2r.x <*.q2r.in> *.q2r.out 傅里叶变换至力常数矩阵，运算结束会生成*.q2r.out和*.fc两个文件。
第四步 建立*.matdyn.in文件，运行matdyn.x <*.matdyn.in > *.matdyn.out生成任意q点下的动力学矩阵，由此可以通过求解久期方程获得格波频率及色散关系。
运行结束会生成*.matdyn.out文件和*.freq文件和*.freq.pq 和*.matdyn.modes四个文件。
第五步 建立*.phdos.in文件，运行matdyn.x <*.phdos.in>*.phdos.out 计算声子态密度，计算结束后会生成*.phdos，*.phdos.out，*.matdyn.freq文件和*.matdyn.freq.pq 四个文件。
******