计算化学公社

标题: 分子优化过程中出现虚频 [打印本页]

作者
Author: ZMINER 时间: 2019-4-15 19:59
标题: 分子优化过程中出现虚频
各位大牛，本人小白，发帖为了求助我在优化分子的过程中出现以下问题（关于虚频处理问题上）软件G09 B.01：
1.我在优化分子时，在B3LYP/6-311+G 水平下做freq计算没有出现虚频（见压缩包中文件1.log），但在B3LYP/6-311+G(d,p) 的水平下做freq计算却出现虚频（见压缩包中文件2.log，如下图），请问这是怎么回事？

(, 下载次数 Times of downloads: 37)

2.在检查问题1中的2.log文件，发现最后一个步却是四个标准YES：

Maximum Force          0.000057    0.000450    YES
RMS    Force          0.000020    0.000300    YES
Maximum Displacement    0.000852    0.001800    YES
RMS    Displacement    0.000249    0.001200    YES

和opt文件是一样的，那为什么会出现虚频，这四个收敛标准值和虚频出现有啥联系？

3.在处理B3LYP/6-311+G(d,p) 的水平下出现的虚频，采用sob老师发的帖子：http://sobereva.com/278 中  做法一，虚频没有解决，接着采用做法二：opt=calcall 成功解决虚频（见压缩包中文件3.log），但是算出来的能量和直接使用opt（文件2.log）算出来的能量有差异，为啥能量是有差异的？

压缩包文件：
(, 下载次数 Times of downloads: 21)

以上疑问希望您能给予解答，谢谢您！

作者
Author: happyknighthawk 时间: 2019-4-15 20:05
见Sob老师的帖子：http://sobereva.com/278

作者
Author: sobereva 时间: 2019-4-15 23:52
B3LYP/6-311+G计算根本毫无意义，Cs对称性下这个级别没出现虚频正说明这个级别太烂，连真正的极小点结构都体现不出来

2.log本来就有虚频，显然没虚频的3.log和这个能量不可能相同

根本就用不着昂贵的calcall，博文里说了，虚频破坏对称性的情况直接按照虚频调结构重新优化就完了。

反复看此文10遍
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈弥散函数和“月份”基组
http://sobereva.com/119
别没事乱加弥散。记住不加极化函数计算结果毫无实际意义，数据根本无法发表

作者
Author: ZMINER 时间: 2019-4-17 10:18

happyknighthawk 发表于 2019-4-15 20:05
见Sob老师的帖子：http://sobereva.com/278

谢谢您的指点

作者
Author: ZMINER 时间: 2019-4-17 10:18

sobereva 发表于 2019-4-15 23:52
B3LYP/6-311+G计算根本毫无意义，Cs对称性下这个级别没出现虚频正说明这个级别太烂，连真正的极小点结构都 ...

谢谢老师的解答。

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