计算化学公社

标题: 高斯计算遇到过渡态的难以消除的甲基单键旋转的虚频 [打印本页]

作者
Author: syx123456 时间: 2023-4-23 22:48
标题: 高斯计算遇到过渡态的难以消除的甲基单键旋转的虚频
各位老师好，我在高斯计算过渡态时，遇到了一个难以消除的虚频。计算使用了ωB97XD/6-31G(d)基组，无其他关键词和IOp输入，使用Gaussian 16 A.03版本，结构优化为为opt(TS)+freq，未使用溶剂模型。

从附件中输出文件中可以看到，这是一个氢转移过程，除了过渡态对应的虚频外，还有不想要的甲基单键旋转的虚频。输出文件在附件的压缩包中。

我尝试了将甲基旋转一定角度后，先进行固定优化（此时冻结过渡态中感兴趣的键），再进行opt(TS,calcfc,noeigen)+freq。虚频大小和振动模式不变。
使用GaussView的Manual Displacement功能将不感兴趣的虚频拉动至最大偏移后，重新进行关键词不变的结构优化，均无效。
使用gentor对该甲基进行单键旋转，然后在相同计算级别下对结构进行单点能计算，结果表明该甲基在此振动模式下处于能量极小点。这一步的目的是进行刚性扫描。
由于高斯不支持笛卡尔坐标下进行刚性扫描，因此数据使用python脚本读取SCF能量后，用matplotlib进行画图。
请问我应该如何消除该虚频，或应该对这种表面上看起来的假虚频进行什么解释？

作者
Author: sobereva 时间: 2023-4-24 03:22
加上int=superfine。如果还不行，opt加上tight，如果还不行，加上recalc=3。maxstep=1这么小根本没必要

固定优化毫无意义

作者
Author: syx123456 时间: 2023-4-24 12:22
好的，谢谢老师！

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