计算化学公社

标题: Gaussian做ts时搜索的过渡态是产物结构不是过渡态结构并且出现不收敛问题 [打印本页]

作者
Author: Zyyyy 时间: 2023-8-19 22:39
标题: Gaussian做ts时搜索的过渡态是产物结构不是过渡态结构并且出现不收敛问题
本帖最后由 Zyyyy 于 2023-8-19 23:06 编辑

在做一条反应路径，有两个过渡态，第一步计算的很好 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比较好

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

作者
Author: Zyyyy 时间: 2023-8-28 16:59
本帖最后由 Zyyyy 于 2023-8-28 17:03 编辑

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

作者
Author: ddddnight 时间: 2023-9-22 10:44

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

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

作者
Author: wzkchem5 时间: 2023-9-22 16:00

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

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

作者
Author: ddddnight 时间: 2023-9-22 17:54

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

好的老师，谢谢您

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