计算化学公社

标题: r2scan-3c优化过渡金属化合物出现的虚频值极大 [打印本页]

作者
Author: neon 时间: 2023-8-4 12:56
标题: r2scan-3c优化过渡金属化合物出现的虚频值极大
大家好，我想用计算过渡金属化合物的spin state splitting energy, 目前用r2scan-3c 优化一45个原子的含Fe的化合物，初始结构是用gfn2-xtb优化的，但是从优化后的结果来看，频率计算含有一个很大的虚频，把.out文件用OfakeG转换后，导入gaussView观察振动模式，发现除了整个分子振动外，箭头所指的O-H键剧烈左右摇摆，相同结构利用b97-3c优化，未发现任何异常，想问一下该如何解释这个虚频，因为测试过很多类似的结构，发现r2scan-3c优化的含铁化合物都含有一个巨大的虚频值，另外，guassview页面底部显示这是一个neutral， singlet，但实际输入文件是+2，5 的high-spin state complex.

作者
Author: yxdd98 时间: 2023-8-4 13:03
本帖最后由 yxdd98 于 2023-8-4 13:07 编辑

没有叫guassview的程序，窗口显示的只是软件对当前体系的判断，如果与实际不符自己在输入文件里修改就成
尝试过消除虚频的各类方法了吗？另外对于r2SCAN-3c可以看ORCA 5.0.4发布 - 量子化学 (Quantum Chemistry) - 计算化学公社 - 第2页 (keinsci.com)的讨论

作者
Author: sobereva 时间: 2023-8-4 15:36
你的截图里明明是B97-3c

作者
Author: neon 时间: 2023-8-4 21:55
谢谢大家的回复，之前有试过sobereva老师帖子的建议，比如提高收敛限，使用精确的hessian计算，但都得到类似的结果，也试过楼上给的建议，使用defgrid3,并做图片中scf设置，但结果还是还有巨大的虚频值。抱歉，放错截图了，但虚频值的确是出现在r2sacn-3c计算中，二者的参数设置完全一样，只是所用方法分别为b97-3c 和r2scan-3c。没有找到图片修改的地方，抱歉无法修改图片错误。

作者
Author: sobereva 时间: 2023-8-4 22:08

neon 发表于 2023-8-4 21:55
谢谢大家的回复，之前有试过sobereva老师帖子的建议，比如提高收敛限，使用精确的hessian计算，但都得到类 ...

如果B97-3c没虚频用B97-3c就完了，对于几何优化目的，r2SCAN-3c不会比B97-3c有明显优势。而如果B97-3c表现很不好的情况（比如对自相互作用误差敏感的体系），用r2SCAN-3c也肯定好不了。我不认为r2SCAN-3c用于opt freq目的比B97-3c有明显额外价值。

作者
Author: chrinide 时间: 2023-8-5 06:46
scan类的所有泛函都有数值问题，对DFT格点要求非常高，如果非要用，建议用更加稳健保守的大格点。而B97-3c是非常数值稳定的泛函，默认格点就可以了。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)