求助 在计算频率的时候出现了比较小的虚频 求教怎么消除虚频

nosugar

sobereva 发表于 2023-6-27 23:28
得弄清楚造成的原因
诸如优化收敛已经足够精确，但由于CUTOFF太低等数值因素造成Hessian太糙而产生虚频 ...

好的 非常谢谢社长！

红烧肉盖饭

老师您好，我目前遇到了消虚频问题。通过jmol软件查看，我的孤立体系计算也存在虚频。具体数据在我的最新论坛博文里。
我一开始想通过jmol可视化，将结构按照虚频旋转的方向保存进行重新优化。但是，我尝试保存失败了。
后续，我又尝试vaspkit的507程序，但是显示说“no imaginary-freq”。告诉我不存在虚频，与vaspkit论坛里，写的也不相同。
或许老师，有没有什么好的消虚频的方式？

我的体系中设置EDIFFG = -0.03，EDIFF = 1E-7，ENCUT = 500

hellozhaoming

sobereva 发表于 2019-5-27 04:02
分子计算，50、100的虚频根本不能随随便便忽略，除非对自己的数据根本不负责任。

1 一方面这已经可能意 ...

我的个人看法：
波数越大的振动越接近局域的谐振子，波数越小的振动谐振子的力常数越接近零，越接近自由平动。波数非常接近零的振动自由度，已经接近于理想自由平动模式了，用简谐近似不合理，但直接忽略也不合理。
如果振动模式是由H、N、O、F元素贡献的，且如果用隐式溶剂模型，比如水溶液，计算振动时并没有考虑振动模式与溶剂氢键网络的耦合。这些接近平动的模式实际上不应该存在，而是与氢键网络的振动耦合，（除非该元素与氢键网络没有相互作用），所以把小波数的振动手动改成水的振动波数可能更合理一些。当然这只是推测，没有严格证明。

yika

granvia 发表于 2019-5-28 01:22
我认真看了一下那个VASP教程，其实没有问题，因为他计算得到乙醇的频率有3*9=27个值，其中三个是出现虚数 ...

老师你好，那这么说是不是不对应振动的虚频无论多大都可以忽略？

yika

granvia 发表于 2019-5-28 01:22
我认真看了一下那个VASP教程，其实没有问题，因为他计算得到乙醇的频率有3*9=27个值，其中三个是出现虚数 ...

老师，我明白了，因为最后几个频率对应的是平动和转动自由度，所以对应的能量应该远小于振动能级，所以只要虚频在平动和转动自由度上，且其能量相较于振动能量占比很少即可忽略，上述理解不知道是否有问题，谢谢老师