用Gaussian搜索二氧化碳环加成反应过渡态，总是搜索不出来

lzu-wyc

各位老师好，我做的二氧化碳环加成反应的反应机理计算，用Gaussian的TS方法搜索二氧化碳插入到开环后的环氧化物这一过程的过渡态时一直搜索不出来正确的过渡态，要不就是虚频和反应方向路径不一样。

Uus/pMeC6H4-/キ

lzu-wyc 发表于 2025-4-11 11:25
可以尝试先优化产物复合物，即先找到C33-O22已经成键、质子已经转移的极小点，再拉长C33-O22搜索过渡态、 ...

那就展示这些计算对应的所有输入输出文件嘛（而且本来应该在1楼就给了），单这样说一句也没法确认当时的计算合不合理，无助于解决问题。

当然也可能设想的反应本就不存在过渡态，论坛里搜“无垒反应”可以看到讨论，这就需要结合更好的计算级别（主要是低档6-31G(d)基组还能按http://sobereva.com/336改进）、足够低的柔性扫描步长上限等等确定。

lzu-wyc

Uus/pMeC6H4-/キ 发表于 2025-4-11 11:05
明白了。可以尝试先优化产物复合物，即先找到C33-O22已经成键、质子已经转移的极小点，再拉长C33-O22搜索 ...

可以尝试先优化产物复合物，即先找到C33-O22已经成键、质子已经转移的极小点，再拉长C33-O22搜索过渡态、跑IRC确认对应于CO2固定的反应。如果不好找，也可结合柔性扫描（http://sobereva.com/474）看看。老师我拿你上面的这些思路跑过，但它最后就优化成起始物或者虚频与反应路径不一样。

Uus/pMeC6H4-/キ

lzu-wyc 发表于 2025-4-11 10:37
有相关的文献，上传的这篇文献第7章就有相关计算

明白了。可以尝试先优化产物复合物，即先找到C33-O22已经成键、质子已经转移的极小点，再拉长C33-O22搜索过渡态、跑IRC确认对应于CO2固定的反应。如果不好找，也可结合柔性扫描（http://sobereva.com/474）看看。还有可能当前体系的质子转移与碳氧成键并不在同一过渡态发生，需要找两步基元过程。

lzu-wyc

Uus/pMeC6H4-/キ 发表于 2025-4-11 10:22
从有机化学常识来说脂肪醇的酸性一般比羧酸要弱，醇与二氧化碳形成的单碳酸酯负离子也倾向于分解，总体来 ...

有相关的文献，上传的这篇文献第7章就有相关计算

Uus/pMeC6H4-/キ

lzu-wyc 发表于 2025-4-11 10:02
好的老师，以后我就不用scf=maxcycle=200这些语言了，就是二氧化碳的33号碳原子插入到22号的氧原子上，然 ...

从有机化学常识来说脂肪醇的酸性一般比羧酸要弱，醇与二氧化碳形成的单碳酸酯负离子也倾向于分解，总体来看设想的反应至少在气相下（因为没使用溶剂模型，不容易稳定反应产生的离子对）并不容易发生；所以首先我想知道有无文献参考或者明确证据支持这样的反应，等确定了再计算过渡态也不迟。

lzu-wyc

Uus/pMeC6H4-/キ 发表于 2025-4-11 09:51
别用opt(maxcyc=200)和scf=maxcycle=200这种没法解决当前问题又徒增耗时的关键词，解决几何优化不收敛的正 ...

好的老师，以后我就不用scf=maxcycle=200这些语言了，就是二氧化碳的33号碳原子插入到22号的氧原子上，然后22号氧原子和21号氢原子的氢键断裂这么一个过程。

Uus/pMeC6H4-/キ

别用opt(maxcyc=200)和scf=maxcycle=200这种没法解决当前问题又徒增耗时的关键词，解决几何优化不收敛的正确方法看http://sobereva.com/164，解决SCF不收敛看http://sobereva.com/61。

仅观察反应物结构不好理解你希望发生的“环加成”“插入到开环”反应是什么，请用键线式等化学语言清晰表达和/或附上可供参考的文献。