目前来说计算声子谱的计算主要还是基于supercell的方法，主要是发展的时间比较长，发展的相对成熟一些。另一种方法是基于线性响应的方法，该方法主要的好处（我个人经验）：计算的结果相对可靠，收敛参数容易控制。只要能量截断和K点收敛了，计算的结果就比较可靠。而且它可以计算任意q点的声子。
但缺点是有时低对称性的结构的低对称性q点不容易收敛。而且对于大体系，计算速度也不是很快。代表软件：abinit， quantum
espresso （pwscf）。
超胞法目前的主流软件是 phonon， phon， fropho， phonopy。 phonon计算的比较准确可靠， 但是要收费。 所以一般人都是选择后几种软件，因为都是免费的。 phon当时是 Dario Alfe 为了计算自由能而写的计算声子谱的软件。但是phon用起来还是比较麻烦的，需要自己写一些脚本。
但天下没有免费的午餐，用免费的软件就是需要一定的经验积累。phon，fropho和phonopy的使用的最大问题就是对于很多体系，会出现离奇的虚频，尤其是gamma点附近，很难算准。 为了克服这个问题，我的建议是：在对原子进行distortion的时候，一定要对原子加上负方向的移动， 这样虽然计算量增加了一倍，但计算会更准确和可靠。
对于phon：
For central differences (LCENTRAL = .T.) the format of the file is the same, but there are twice as many displacements (for each displacement u there is also -u).
在INPHON里要加上 LCENTRAL = .T. ， 这样就可以产生负方向的移动了。
对于fropho：
More di
超胞法目前的主流软件是 phonon， phon， fropho， phonopy。 phonon计算的比较准确可靠， 但是要收费。 所以一般人都是选择后几种软件，因为都是免费的。 phon当时是 Dario Alfe 为了计算自由能而写的计算声子谱的软件。但是phon用起来还是比较麻烦的，需要自己写一些脚本。
但天下没有免费的午餐，用免费的软件就是需要一定的经验积累。phon，fropho和phonopy的使用的最大问题就是对于很多体系，会出现离奇的虚频，尤其是gamma点附近，很难算准。 为了克服这个问题，我的建议是：在对原子进行distortion的时候，一定要对原子加上负方向的移动， 这样虽然计算量增加了一倍，但计算会更准确和可靠。
对于phon：
For central differences (LCENTRAL = .T.) the format of the file is the same, but there are twice as many displacements (for each displacement u there is also -u).
在INPHON里要加上 LCENTRAL = .T. ， 这样就可以产生负方向的移动了。
对于fropho：
More di