QE计算分子晶体有虚频

破晓时的曙光

我用QE计算一个分子晶体的gamma点频率，优化结构收敛了，但是gamma点频率好几个低频模式都是虚频，请问该如何解决呢

输入文件
&inputph
  tr2_ph=1.0d-12,
  prefix='phonon',
  epsil=.true.,
  amass(1)=12.0107,
  amass(2)=1.00794,
  outdir='./out/',
  fildyn='phonon.dynG',
/
0.0 0.0 0.0

kyuu

先crystal sum下po张频率再说