计算化学公社

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

[Gaussian/gview] 氢原子转移反应的过渡态IRC跑三个点就结束了求助

[复制链接 Copy URL]

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

如图所示,计算的是羟自由基进攻葡萄糖1号碳上的氢发生氢原子转移反应的过渡态。首先采用ub3lyp/6-31+g(d,p) scrf=(pcm,solvent=water) opt=modredundant进行羟基的氧原子与葡萄糖1号碳上的氢原子进行键角柔性扫描,得到的能量最高点进行ub3lyp/6-31+g(d,p) scrf=(pcm,solvent=water) opt=(ts,calcfc,noeigen) freq过渡到态搜寻,虚频方向对的。最后采用ub3lyp/6-31+g(d,p) scrf=(pcm,solvent=water) irc(calcfc) guess=read 进行IRC搜索后每次都只跑3步数,采用了tight int=ultrafine跑过渡态,步长改15跑IRC也是三步数停止,步长改20时收敛失败。请问各位老师应该怎么操作。

202409250824478354..png (22.92 KB, 下载次数 Times of downloads: 5)

IRC结果

IRC结果

202409250821564403..png (17.97 KB, 下载次数 Times of downloads: 4)

柔性扫描结果

柔性扫描结果

202409250811474786..png (53.04 KB, 下载次数 Times of downloads: 4)

计算咪表

计算咪表

glucose oh res opt highest p4irc.gjf

2.01 KB, 下载次数 Times of downloads: 9

IRC.gjf

glucose oh res opt highest p4irc.out

211.74 KB, 下载次数 Times of downloads: 3

IRC.out

1063

帖子

0

威望

3277

eV
积分
4340

Level 6 (一方通行)

2#
发表于 Post on 2024-9-25 09:33:23 | 只看该作者 Only view this author
羟自由基进攻葡萄糖1号碳上的氢发生氢原子转移反应的过渡态


你的目标是让葡萄糖1号碳上的氢和 OH 生成 H2O 分子吗?
如果是,那么从你贴的第3张图来看,过渡态可能还没找到,1号碳上 C-H 之间的距离可能还太近,OH 上的 O 和 1号碳上的 H 之间的距离还太远。

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

3#
 楼主 Author| 发表于 Post on 2024-9-25 10:35:49 | 只看该作者 Only view this author
目标是生成H2O分子。老师,这是他的震动的样子,从这张图上来看应该可以生成H2O分子吧

202409251034335144..png (37.98 KB, 下载次数 Times of downloads: 2)

202409251034335144..png

5万

帖子

99

威望

5万

eV
积分
112551

管理员

公社社长

4#
发表于 Post on 2024-9-25 22:08:41 | 只看该作者 Only view this author
该说的都说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400http://bbs.keinsci.com/thread-7736-1-1.html
明确说了,opt和IRC的积分格点设置必须严格一致
另外,牵扯到分子内弱相互作用,最好加上D3
谈谈“计算时是否需要加DFT-D3色散校正?”
http://sobereva.com/413

对于加大步长,可以进一步测试30。SCF不收敛该怎么解决怎么解决
解决SCF不收敛问题的方法
http://sobereva.com/61
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

5

帖子

0

威望

19

eV
积分
24

Level 1 能力者

5#
 楼主 Author| 发表于 Post on 2024-9-27 03:12:59 | 只看该作者 Only view this author
sobereva 发表于 2024-9-25 22:08
该说的都说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400(http://bbs.keinsci.com/t ...

谢谢sob老师!

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 21:06 , Processed in 0.191148 second(s), 24 queries , Gzip On.

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