计算化学公社

 找回密码 Forget password
 注册 Register
Views: 422|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

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

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

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

2#
发表于 Post on 2025-9-29 13:43:27 | 只看该作者 Only view this author
建议先把共价键连接的单分子产物结构优化出来,再拉长关键的碳氧单键来找过渡态。B3LYP应当加上D3(BJ)色散校正做这种优化,并且现实的溶剂环境、卡宾前体产生卡宾的方式也得考虑一下。另外有可能该计算级别下此基元反应就是没有电子能垒,改进基组或者更换泛函或许能得到不同结论。
√546=23.36664289109

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-9-30 17:05:36 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-9-29 13:43
建议先把共价键连接的单分子产物结构优化出来,再拉长关键的碳氧单键来找过渡态。B3LYP应当加上D3(BJ)色散 ...

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

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

4#
发表于 Post on 2025-9-30 18:41:34 | 只看该作者 Only view this author
shobuu 发表于 2025-9-30 17:05
感谢您的回复!
我按你的建议试一下。
卡宾前体是卡宾C上有一个=N基团的分子,光照产生卡宾的。

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

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2025-9-30 22:41:05 | 只看该作者 Only view this author
本帖最后由 shobuu 于 2025-9-30 23:28 编辑
Uus/pMeC6H4-/キ 发表于 2025-9-30 18:41
有没有过渡金属催化剂?光照这种重要条件应当在一开始就讲清楚,因为必须考虑电子激发及激发态反应的问题 ...

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

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-9-30 23:35:46 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-9-30 18:41
有没有过渡金属催化剂?光照这种重要条件应当在一开始就讲清楚,因为必须考虑电子激发及激发态反应的问题 ...

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

4

帖子

0

威望

362

eV
积分
366

Level 3 能力者

7#
发表于 Post on 2025-10-7 09:34:00 | 只看该作者 Only view this author
有可能是你的stepsize设的太小了,高斯跑了两步误以为已经到极小点就停下了,可以调成30或更大(默认是10)

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-10-17 22:38:13 | 只看该作者 Only view this author
天水丨伯约 发表于 2025-10-7 09:34
有可能是你的stepsize设的太小了,高斯跑了两步误以为已经到极小点就停下了,可以调成30或更大(默认是10)

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

15

帖子

0

威望

135

eV
积分
150

Level 3 能力者

9#
发表于 Post on 2025-12-22 22:00:17 | 只看该作者 Only view this author
shobuu 发表于 2025-10-17 22:38
之前确实没注意到这个,谢谢解答

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

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2025-12-22 23:40:36 | 只看该作者 Only view this author
leduo 发表于 2025-12-22 22:00
请问你解决了吗?我也是遇到了类似的问题。我对找到的过渡态进行IRC计算,即便将stepsize=30,也只能跑11 ...

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

15

帖子

0

威望

135

eV
积分
150

Level 3 能力者

11#
发表于 Post on 2025-12-23 09:53:17 | 只看该作者 Only view this author
shobuu 发表于 2025-12-22 23:40
我这个是两个分子结合的,能垒很小,后来过渡态找到了,虚频方向也不太对,暂时放弃了。但是算其他反应的 ...

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

134

帖子

1

威望

4105

eV
积分
4259

Level 6 (一方通行)

12#
发表于 Post on 2025-12-23 12:21:49 | 只看该作者 Only view this author
Gaussian 的 IRC 确实走不到底,此时可以考虑用优化代替,也就是用 IRC 两端的结构去优化,最后手动算出完整的 IRC
十八介姑娘一蕾花呀,白白介牙齿、红红介嘴唇,得人惜

68

帖子

0

威望

332

eV
积分
400

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2025-12-27 23:10:08 | 只看该作者 Only view this author
leduo 发表于 2025-12-23 09:53
谢谢。我找的这个过渡态我感觉应该是正确的,但是就是IRC跑不出来,虚频只有-30+,但是振动方向也是对的 ...

楼下回答了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-23 19:49 , Processed in 0.214058 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list