Gaussian做ts时搜索的过渡态是产物结构 不是过渡态结构 并且出现不收敛问题

Zyyyy

在做一条反应路径，有两个过渡态，第一步计算的很好 IRC跑出来很光滑 然后我使用IRC结果的产物 作为下一步TS计算初始构型，并做了一点构型上的调整 更趋向于发生反应 然后进行TS计算
#p opt=(ts,calcfc,noeigen) freq b3lyp/6-311g(d,p) em=gd3bj
然后 跑了200步 报错了 不收敛
Error termination request processed by link 9999.
然后我打开看 发现现在的构型是产物的构型
这步是一个成环反应 是不是用QST2比较好

wzkchem5

建议对所有氧加弥散函数。首先，羧酸根的两个氧带明显负电，必须加弥散；参与亲核进攻的氧会在进攻过程中带有一些负电，即便没有生成游离的醇负离子，也可以预期负电会比较多，所以最好也加弥散；酯的羰基氧在亲核进攻过程中会生成四面体中间体的结构（实际未必是一个中间体，也可能在IRC的坡上，而非一个驻点），也会带明显负电，所以也要加弥散；最后剩下那个氧自始至终都不带明显负电，但产物的碳酸酯结构里的两个五元环氧几乎等价，已经有一个五元环氧加了弥散的话，另一个不加会导致高估两个氧的不等价性，所以最好也加弥散。
解决这个问题以后，如果仍然找不到过渡态，应当考虑第二步亲核进攻和第一步是协同异步反应的可能性，方法是IRC跑出来的中间体先不要急着用来找过渡态，而是做opt优化到收敛，如果自动变成了最终产物，就说明第二步没有单独的过渡态

Zyyyy

非常感谢老师的帮助 加了弥散之后 经过几次计算 已经算出了过渡态
然后进行IRC 发现路径结果不是很好 有一个方向好像没算完 可能是算到步数了 我在计算时候关键词写的 maxpoints=200, stepsize=5 但是算出来的结果是和我预期的结果是一致的 想问问老师 这样的结果可以使用吗

ddddnight

wzkchem5 发表于 2023-8-19 23:56
建议对所有氧加弥散函数。首先，羧酸根的两个氧带明显负电，必须加弥散；参与亲核进攻的氧会在进攻过程中带 ...

老师，看到您回答的，结合卢老师关于是否加弥散函数的博文，对于这种含有带负电比较多的O的体系，算过渡态都应该对O加弥散。请问，是不是对于这种体系，简单结构优化的时候要不要对O加弥散函数呢？

wzkchem5

ddddnight 发表于 2023-9-22 03:44
老师，看到您回答的，结合卢老师关于是否加弥散函数的博文，对于这种含有带负电比较多的O的体系，算过渡 ...

也要加。这里有两个额外因素：（1）结构优化所需的精度不如单点能所需精度高，导致加弥散的必要性降低；（2）结构优化用的基组比较小，同等条件下小基组对弥散函数的需求比大基组大，导致加弥散的必要性增加。总的来说，结构优化加弥散的必要性和单点能不会有太大差别

ddddnight

wzkchem5 发表于 2023-9-22 16:00
也要加。这里有两个额外因素：（1）结构优化所需的精度不如单点能所需精度高，导致加弥散的必要性降低； ...

好的老师，谢谢您