ionexchangeC 发表于 2024-2-9 22:22
很多时候必须要int=superfine才能消小的虚频，为此统一之前所有的计算内容就很麻烦也没必要。我的处理方法 ...

这是正常处理
为了获得热力学量而做振动分析，我不要求所有体系振动分析的积分格点严格一致，不同体系的opt freq过程我接受ultrafine和superfine和混用

Stars 发表于 2024-2-9 21:45
从2.5里面可以看到，小虚频存在时，其他频率更像是在真实频率的基础上向负方向平移了一定波数，这个在论坛 ...

这等同于QRRHO(Truhlar)并把Shermo的ravib设为50的情况

"其他频率更像是在真实频率的基础上向负方向平移了一定波数" 当前例子中体现了这一点，但也不算是非常普遍的现象，比如一些情况只有个甲基旋转的小虚频，由于与其它自由度耦合很弱，消不消虚频往往对其它低频模式频率没什么影响。正因为有无小虚频结构下低实频的差异有不同的可能性，所以带着小虚频算熵没有唯一理想做法。虽然小虚频当实频同时用QRRHO(Truhlar)可以极大抑制这种差异，但可惜QRRHO(Truhlar)本身却又不够理想。

很多时候必须要int=superfine才能消小的虚频，为此统一之前所有的计算内容就很麻烦也没必要。我的处理方法是先用int=superfine格点优化+算频率得到校正量，然后再拿默认设置（比如int=ultrafine）算个单点再加和。

从2.5里面可以看到，小虚频存在时，其他频率更像是在真实频率的基础上向负方向平移了一定波数，这个在论坛此前也有讨论，不清楚有无将所有频率正移若干波数的处理方式；
Gaussian这种有解析Hessian的还好，此前论坛里有讨论在VASP中遇到虚频的处理方式http://bbs.keinsci.com/thread-13344-1-1.html
那种只有数值方法求频率的情况，出现虚频的原因就可能比较复杂了，到底是数值噪声，还是结构优化不到位，有时还需要反复尝试。常用的vaspkit处理方式不仅不考虑虚频，而且统一将小于某阈值（一般是50）cm-1的实频波数直接按照50对待。其实感觉哪怕不存在虚频，帖子的方法对实频计算的误差也有意义