用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文件。

冲冲冲

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

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

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

冲冲冲

这是两个过渡态的震动方式

冲冲冲

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

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

sobereva

不要把6-311G(d,p)写成6311g(d,p)，显得非常缺乏专业素养

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

阿甘

我看了下计算结果，振动方式没变啊……

冲冲冲

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

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

wzkchem5

说明这个过渡态可能能垒比较低，这种情况下一个级别能找到过渡态、另一个级别找不到都是家常便饭，区区一个虚频变了都不算事。