卡宾+THF IRC跑两步就停了，求助结构的调整

shobuu · 发表于 Post on 2025-9-29 10:54:36

大家好，计算如图的成键过程的能垒，能搜索到过渡态，也只有一个虚频，但是IRC跑两步就结束了。看卢老师教程，初步判断为过渡态结构不对。可以请教一下大家这个结构可以怎么调整吗？谢谢~

TS1-222 (1).zip (129.51 KB, 下载次数 Times of downloads: 6)

Uus/pMeC6H4-/キ · 发表于 Post on 2025-9-29 13:43:27

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

shobuu · 发表于 Post on 2025-9-30 17:05:36

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

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

Uus/pMeC6H4-/キ · 发表于 Post on 2025-9-30 18:41:34

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

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

shobuu · 发表于 Post on 2025-9-30 22:41:05

本帖最后由 shobuu 于 2025-9-30 23:28 编辑

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

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

shobuu · 发表于 Post on 2025-9-30 23:35:46

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

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

天水丨伯约 · 发表于 Post on 2025-10-7 09:34:00

有可能是你的stepsize设的太小了，高斯跑了两步误以为已经到极小点就停下了，可以调成30或更大（默认是10）

shobuu · 发表于 Post on 2025-10-17 22:38:13

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

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

leduo · 发表于 Post on 2025-12-22 22:00:17

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

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

shobuu · 发表于 Post on 2025-12-22 23:40:36

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

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

leduo · 发表于 Post on 2025-12-23 09:53:17

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

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

王二葛 · 发表于 Post on 2025-12-23 12:21:49

Gaussian 的 IRC 确实走不到底，此时可以考虑用优化代替，也就是用 IRC 两端的结构去优化，最后手动算出完整的 IRC

shobuu · 发表于 Post on 2025-12-27 23:10:08

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

楼下回答了

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 卡宾+THF IRC跑两步就停了，求助结构的调整