计算化学公社

标题: 卡宾+THF IRC跑两步就停了，求助结构的调整 [打印本页]

作者
Author: shobuu 时间: 2025-9-29 10:54
标题: 卡宾+THF IRC跑两步就停了，求助结构的调整
大家好，计算如图的成键过程的能垒，能搜索到过渡态，也只有一个虚频，但是IRC跑两步就结束了。看卢老师教程，初步判断为过渡态结构不对。可以请教一下大家这个结构可以怎么调整吗？谢谢~
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 6)

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-9-29 13:43
建议先把共价键连接的单分子产物结构优化出来，再拉长关键的碳氧单键来找过渡态。B3LYP应当加上D3(BJ)色散校正做这种优化，并且现实的溶剂环境、卡宾前体产生卡宾的方式也得考虑一下。另外有可能该计算级别下此基元反应就是没有电子能垒，改进基组或者更换泛函或许能得到不同结论。

作者
Author: shobuu 时间: 2025-9-30 17:05

Uus/pMeC6H4-/キ发表于 2025-9-29 13:43
建议先把共价键连接的单分子产物结构优化出来，再拉长关键的碳氧单键来找过渡态。B3LYP应当加上D3(BJ)色散 ...

感谢您的回复！
我按你的建议试一下。
卡宾前体是卡宾C上有一个=N基团的分子，光照产生卡宾的。

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-9-30 18:41

shobuu 发表于 2025-9-30 17:05
感谢您的回复！
我按你的建议试一下。
卡宾前体是卡宾C上有一个=N基团的分子，光照产生卡宾的。

有没有过渡金属催化剂？光照这种重要条件应当在一开始就讲清楚，因为必须考虑电子激发及激发态反应的问题，只像一楼的计算那样还按基态/闭壳层单重态处理是不妥的。这方面我不熟，或许可以请教如wzkchem5老师等人

作者
Author: shobuu 时间: 2025-9-30 22:41
本帖最后由 shobuu 于 2025-9-30 23:28 编辑

Uus/pMeC6H4-/キ发表于 2025-9-30 18:41
有没有过渡金属催化剂？光照这种重要条件应当在一开始就讲清楚，因为必须考虑电子激发及激发态反应的问题 ...

谢谢您的指正！体系里是没有过渡金属催化剂的。我们认为光照可能只对卡宾的产生起作用。卡宾和THF结合的过程猜测是基态反应。后续有更多结果我再补充进来。

作者
Author: shobuu 时间: 2025-9-30 23:35

Uus/pMeC6H4-/キ发表于 2025-9-30 18:41
有没有过渡金属催化剂？光照这种重要条件应当在一开始就讲清楚，因为必须考虑电子激发及激发态反应的问题 ...

抱歉，前面发了点错误信息00~

作者
Author: 天水丨伯约 时间: 2025-10-7 09:34
有可能是你的stepsize设的太小了，高斯跑了两步误以为已经到极小点就停下了，可以调成30或更大（默认是10）

作者
Author: shobuu 时间: 2025-10-17 22:38

天水丨伯约发表于 2025-10-7 09:34
有可能是你的stepsize设的太小了，高斯跑了两步误以为已经到极小点就停下了，可以调成30或更大（默认是10）

之前确实没注意到这个，谢谢解答

作者
Author: leduo 时间: 2025-12-22 22:00

shobuu 发表于 2025-10-17 22:38
之前确实没注意到这个，谢谢解答

请问你解决了吗？我也是遇到了类似的问题。我对找到的过渡态进行IRC计算，即便将stepsize=30，也只能跑11步，优化也得不到对应中间体。我也参考了sob老师的IRC常见问题的帖子，我认为可能是过渡态没有找准确，可能找到了过渡态旁边红点这种位置（sob老师IRC常见问题帖子里那张图，我不知道怎么上传不了图片）。但是有个疑问是：正确的过渡态结构是与原本找到的过渡态结构有较小的区别，还是说可能有较大的区别呢？

作者
Author: shobuu 时间: 2025-12-22 23:40

leduo 发表于 2025-12-22 22:00
请问你解决了吗？我也是遇到了类似的问题。我对找到的过渡态进行IRC计算，即便将stepsize=30，也只能跑11 ...

我这个是两个分子结合的，能垒很小，后来过渡态找到了，虚频方向也不太对，暂时放弃了。但是算其他反应的时候，我会按照文献的过渡态结构，摆成和文献差不多的，然后往其中一边偏一点，做柔性势能面扫描，找到最高点，再从两边找两个点，做QST3的计算。就能做成功了。没有文献的话，可以从反应物，产物出发，各自做一个键长拉大的柔性势能面扫描？

作者
Author: leduo 时间: 2025-12-23 09:53

shobuu 发表于 2025-12-22 23:40
我这个是两个分子结合的，能垒很小，后来过渡态找到了，虚频方向也不太对，暂时放弃了。但是算其他反应的 ...

谢谢。我找的这个过渡态我感觉应该是正确的，但是就是IRC跑不出来，虚频只有-30+，但是振动方向也是对的，不知道为什么。等过了今天我发个帖子问问

作者
Author: 王二葛 时间: 2025-12-23 12:21
Gaussian 的 IRC 确实走不到底，此时可以考虑用优化代替，也就是用 IRC 两端的结构去优化，最后手动算出完整的 IRC

作者
Author: shobuu 时间: 2025-12-27 23:10

leduo 发表于 2025-12-23 09:53
谢谢。我找的这个过渡态我感觉应该是正确的，但是就是IRC跑不出来，虚频只有-30+，但是振动方向也是对的 ...

楼下回答了

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