最近一直在折腾qe的事情，总结几点经验：
1. phonopy的超胞法和位移法在胞不是足够大的时候，可能有一些电荷密度波的贡献算不出来。所以最好有时间用qe折腾一遍，尤其是一些比较零碎的组分和低对称性结构。
2. qe有时候会出现gamma点处的声学波虚频，有的还比较大（>1THz）同时会连带着周围的q点也有一点虚频，消除方法是增加ph.in里面的tr2_ph=1e-14或-15
3. 对于qe结构弛豫不收敛的情况，增加 ecutrho可以解决（大多数情况增加ecutrho并不明显增加计算时间）。但比较奇怪的是大多数情况增加ecutrho对力和应力的改变非常有限。