计算化学公社

标题: 用b3lyp/631g(d,p)优化的过渡态，再用大基组b3lyp/6311g(d,p)优化后，震动方式变了 [打印本页]

作者
Author: 冲冲冲 时间: 2022-11-1 11:02
标题: 用b3lyp/631g(d,p)优化的过渡态，再用大基组b3lyp/6311g(d,p)优化后，震动方式变了
b3lyp/631g(d,p)优化好的过渡态，再用大基组b3lyp/6311g(d,p)优化后，震动方式变了，虚频从-322降到了-130，这是因为什么呢？
已经尝试了重新摆过渡态结构直接用b3lyp/6311g(d,p)优化，但结果仍未变化，尝试直接用b3lyp/6311g(d,p)柔性扫描，也未找到。
实在不知道怎么办了，用OPT=loose放宽收敛限能行吗？各位老师有没有好的建议呢

附件中放了b3lyp/631g(d,p)优化的log文件，IRC的log文件，以及b3lyp/6311g(d,p)优化的log文件。

作者
Author: wzkchem5 时间: 2022-11-1 13:44
说明这个过渡态可能能垒比较低，这种情况下一个级别能找到过渡态、另一个级别找不到都是家常便饭，区区一个虚频变了都不算事。

作者
Author: 冲冲冲 时间: 2022-11-1 14:55

wzkchem5 发表于 2022-11-1 13:44
说明这个过渡态可能能垒比较低，这种情况下一个级别能找到过渡态、另一个级别找不到都是家常便饭，区区一个 ...

这个过渡态是很关键的一步，请问您知道可以继续尝试什么办法去找到它吗？

作者
Author: 阿甘 时间: 2022-11-1 15:33
我看了下计算结果，振动方式没变啊……

作者
Author: sobereva 时间: 2022-11-1 15:37
不要把6-311G(d,p)写成6311g(d,p)，显得非常缺乏专业素养

更大的基组下死活找不到小基组下能找到的，说明之前基组太小，对势能面描述不正确，本来就不该有那个过渡态（亦或者是基组与理论方法误差抵消的巧合）。如果你觉得那个过渡态就该存在，尝试其它泛函诸如wB97XD，并且恰当考虑溶剂效应等因素。如果在极为合理的计算方式下就是找不到，说明过渡态就是不存在。

作者
Author: 冲冲冲 时间: 2022-11-1 15:43

阿甘发表于 2022-11-1 15:33
我看了下计算结果，振动方式没变啊……

震动方式确实是变了，6311G的结果中芳香环震动的很复杂，IRC结果跑出来也很不一样

作者
Author: 冲冲冲 时间: 2022-11-1 15:44
这是两个过渡态的震动方式

作者
Author: 冲冲冲 时间: 2022-11-1 16:42

sobereva 发表于 2022-11-1 15:37
不要把6-311G(d,p)写成6311g(d,p)，显得非常缺乏专业素养

更大的基组下死活找不到小基组下能找到的，说 ...

谢谢老师，我再试试其他泛函

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