优化结构之后对于Gaussian和vasp计算频率的疑问

薛定谔的老猫

诸位大佬好，有一个疑问还请大佬指教，在采用Gaussian优化结构之后，进行频率计算是必须的，以确定分子是不是真的处在势能面的极小值点，而且在采用Gaussian计算的文章中一般都会强调这一点；但是vasp在几何结构优化之后，很少看到有进行频率计算的，文章中也很少提及这一点。当然vasp一般计算晶体，通过声子谱来看结构的稳定性，可是很多文章中除非专门研究声子谱的，也很少会在结构优化之后计算声子谱。于是我的疑问就是，vasp优化结构之后难道不要进行频率或者声子谱计算看一下有无虚频？如果需要，为啥很少有文章提及？感谢！

卡开发发

1、Gaussian或ORCA真的做上百原子的体系恐怕算频率的人也就很少了，最主要的原因还是计算量太大。如果真的按照原则，周期体系的计算是应该检验声子谱。对于周期体系，即便是小的晶格可能要牵扯到声子的BZ采样，要么就是需要造超晶格，这样计算速度也就很慢。就算是同样做，稍大点的体系LCAO的速度总是比平面波快多了。
3、对于固体，牵扯到变晶格优化，数值方面不是很稳定；对表面体系，周期性边界会引起一些旋转不变的破缺，吸附的分子可能会相对晶格发生一些内转动。这些类似的因素都会导致计算容易出现虚频，而这些虚频消除一般也比较麻烦。
4、做过渡态的情况还是有人算，不过很多人也不是算全频率，这种做法可靠性值得商榷。另外，分析工具确实也不及计算小分子方面的工具来的全面。

薛定谔的老猫

卡开发发 发表于 2020-3-23 14:47
1、Gaussian或ORCA真的做上百原子的体系恐怕算频率的人也就很少了，最主要的原因还是计算量太大。如果真的 ...

哈哈，原来是这样，感谢！