VASP+Phonopy计算声子谱

在第一性原理计算过程中，研究体系的稳定性是经常遇到的一个问题，我们可以从很多方面来解释，其中可以使用声子谱研究体系的动力学稳定性。

在固体理论中，声子是晶格振动的简正模能量量子，声子用来描述晶格的简谐振动。在量子力学中，固体内存在原子核之间的相互作用、电子间的相互作用还有原子核与电子间的相互作用。其中，电子的运动规律用密度泛函理论得到，而原子核的运动规律则用声子来描述。

目前，声子谱是研究材料热力学性质的一个很好的切入点，对于三维块体材料，声子谱分光学波（高）和声学波频率（低），当声子谱全部在0点以上，说明材料没有出现虚频，也就是说材料是相对稳定存在的。

那么，如何计算声子谱是我们关注的重点。目前计算声子谱的方法有两种，分别是直接法即有限位移法和密度泛函微扰理论。

1.有限位移法，或称Finite displacement方法

通过在优化后的平衡结构中引入原子位移，计算作用在原子上的Hellmann-Feynman力，进而由动力学矩阵算出声子色散曲线。

2.密度泛函微扰理论

密度泛函微扰理论或称DFPT，通过计算系统能量对外场微扰的响应来求出晶格动力学性质，直接计算出原子的移动而导致的势场变化，再进一步构造出动力学矩阵，进而算出声子谱。

这两种方法的计算方式略有不同&#x