计算化学公社

标题: 求助：Gaussian优化复合物总是出现很小的虚频 [打印本页]

作者
Author: suosuosky 时间: 2019-4-18 10:50
标题: 求助：Gaussian优化复合物总是出现很小的虚频
请问各位老师，我在使用g16优化复合物的构型是，总是出现-5~-4 大小的虚频，反复微调构型多次后，最后还是消不掉，请问应如何解决呢?
我目前是采用B3LYP-D3优化构型，然后使用B2PLYP/jun-cc-pVTZ计算结合能。大部分构型优化没有问题，只有其中1-2个出现上述消不掉小虚频的情况。

输入文件如下：
#p opt freq b3lyp/6-311+g(d,p) empiricaldispersion=gd3bj output=wfn

作者
Author: suosuosky 时间: 2019-4-18 13:38
本帖最后由 suosuosky 于 2019-4-18 17:09 编辑

把GD3BJ改成GD3，用opt=tight int=ultrafine 都试过还是一直有 -5 左右的小虚频，用m062x/6-311+G**优化相同的构型，则没有虚频，分子间的整体构型想死，但是构型有部分差异，包括分子间的距离、取代基的角度等等。

（1）请问这种情况下，我是否可以使用m062x/6-311+G**优化好的构型来计算BSSE校正的结合能，这个数值能否与之前其他构型的B3LYP-D3/6-311+G**基础上得到的BSSE校正的结合能进行大小比较？因为其他复合物都已经用B3LYP-D3/6-311+G**优化并计算了能量，如果全部再用m062x/6-311+G**，需要重新计算的体系就很多了。。。

（2）此外，我还需要比较复合物和单体的某些键长的变化，m062x/6-311+G**得到的复合物构型是不是必须与m062x/6-311+G**的单体进行比较，B3LYP-D3/6-311+G**得到的复合物构型是不是必须与B3LYP-D3/6-311+G**的单体进行比较。还是说小虚频的构型也可以直接使用，与单体比较键长？

谢谢！

作者
Author: sobereva 时间: 2019-4-18 17:29
用m062x/6-311+G**或许这次碰巧没虚频，但你用多了这个级别，以后会碰到更多虚频，到时候就知道痛苦了

去掉弥散函数再优化，弥散函数对于弱相互作用体系优化意义不大，还容易导致零碎的数值问题

作者
Author: suosuosky 时间: 2019-4-18 18:08
本帖最后由 suosuosky 于 2019-4-18 18:13 编辑

sobereva 发表于 2019-4-18 17:29
用m062x/6-311+G**或许这次碰巧没虚频，但你用多了这个级别，以后会碰到更多虚频，到时候就知道痛苦了

...

谢谢sob老师，我试了一下，还是有虚频。应该如何解决呢？
计算的是卤键作用(计算体系小分子为BrCF3)，log文件已经上传。

作者
Author: sobereva 时间: 2019-4-19 03:33

suosuosky 发表于 2019-4-18 18:08
谢谢sob老师，我试了一下，还是有虚频。应该如何解决呢？
计算的是卤键作用(计算体系小分子为BrCF3)，lo ...

直接上calcall结合opt=tight

作者
Author: suosuosky 时间: 2019-4-19 14:12

sobereva 发表于 2019-4-19 03:33
直接上calcall结合opt=tight

谢谢老师，我试了，计算结果还是与小的虚频。附件是log文件。

我能否以现在的构型比较构型，并且用双杂化泛函计算结合能。

作者
Author: sobereva 时间: 2019-4-19 22:12

suosuosky 发表于 2019-4-19 14:12
谢谢老师，我试了，计算结果还是与小的虚频。附件是log文件。

我能否以现在的构型比较构型，并且用双 ...

用int=superfine再试，如果还不行按照虚频轻微调节结构再试
不管你最后用什么算能量，虚频还是得消（虽然对计算结合能计算其实没有影响。如果你不消的话别让审稿人知道）

作者
Author: 李总把头 时间: 2021-4-6 21:47
楼主问下是怎么优化复合物的方便说下嘛谢谢

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