用PES得到过渡态初猜结构使用TS寻找过渡态优化不收敛和振动虚频不在反应方向上

吸血的蚊子

各位老师好，我使用gaussian计算喹啉N夺取乙硒醇的质子的自由能垒，首先需要找到过渡态。寻找过渡态的过程如下：首先分别优化质子转移前后的化合物结构，然后以N,H两个原子的距离为扫描变量，使用柔性扫描找可能的过渡态结构。1.柔性扫描,关键词为：opt=modredundant b3lyp/6-31+g(d,p) scrf em=gd3bj ；步数为10，步长为-0.11埃。
2.查找过渡态，扫描结果出现两个峰，以最高点作为过渡态的初猜结构，使用opt=TS查找过渡态，关键词为opt=(ts,calcfc,noeigen,tight) freq b3lyp/6-31+g(d,p) scrf em=gd3bj int=ultrafine，第一次优化不收敛，查看结果优化并未出现震荡，应该是初猜过渡态结构不合理。
3.关键词同上，再次使用柔性扫描，缩短步长至-0.05埃，仍旧是两个峰，以最高峰作为过渡态初猜结构再次优化，结果收敛但对应的虚频不是在反应坐标方向上.
请问是我势能面扫描不够准确的原因吗？应该如何解决呢？谢谢。

wzkchem5

试试优化较低的那个峰对应的结构。你选的那个峰虽然较高，但是可能属于扫描坐标选取不合理导致的假最高点（因为势能曲线看起来在这个点不连续），参见DOI: 10.1002/wcms.34的Figure 3

dnlx

另外，质子转移反应有可能不存在过渡态。
http://bbs.keinsci.com/thread-26076-1-1.html

吸血的蚊子

wzkchem5 发表于 2023-4-28 19:20
试试优化较低的那个峰对应的结构。你选的那个峰虽然较高，但是可能属于扫描坐标选取不合理导致的假最高点（ ...

谢谢老师，根据您的建议进行了尝试，找到的过渡态优化收敛，且虚频方向在反应坐标上。

吸血的蚊子

dnlx 发表于 2023-4-29 11:50
另外，质子转移反应有可能不存在过渡态。
http://bbs.keinsci.com/thread-26076-1-1.html

好的，多谢提醒