计算化学公社

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

[新手求助] 求助C-S键断键的过渡态IRC验证均为反应物

[复制链接 Copy URL]

87

帖子

0

威望

1101

eV
积分
1188

Level 4 (黑子)

拟考察C-S键断键的过渡态,跑IRC时总是遇到过渡态虚频,震荡趋势都正确,但是跑irc时能量都不在最高点,且两端都是反应物的情况。

压缩包文件及其关键词:
TS24: #p opt=(calcall,tight,ts,noeigen) freq ub3lyp/6-311+g(d,p) scrf=(solvent=water) nosymm em=gd3bj integral=ultrafine 成功收敛且只有一个虚频,振动方向也是C-S键断键方向。
IRC24:#p IRC(LQA,calcfc,recalc=10,stepsize=30,maxpoints=50) ub3lyp/6-311+g(d,p) scrf=(solvent=water) nosymm guess=read geom=allcheck em=gd3bj integral=ultrafine oldchk文件是TS的chk文件,先前也尝试过 IRC(calcfc,ralcall) IRC(calcfc,ralcall,LQA) IRC(calcfc,ralcall,GS2) IRC(LQA,calcfc,ralcall,stepsize=40) 但都是类似的结果 当加入forward或reverse关键词时,forward可以产生曲线且为成键方向,但reverse第二个点能量略升高呈断键趋势,而后开始下降呈成键方向。
请教老师,这种情况要怎么才能解决呢?我已经困在这问题里很久了,恳请老师指点。

202510041259353228..png (47.67 KB, 下载次数 Times of downloads: 0)

forward方向

forward方向

202510041258284628..png (54.36 KB, 下载次数 Times of downloads: 0)

reverse方向

reverse方向

TS_IRC.rar

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

1209

帖子

6

威望

2442

eV
积分
3771

Level 5 (御坂)

傻傻的木瓜

2#
发表于 Post on 2025-10-4 13:28:51 | 只看该作者 Only view this author
NTf2阴离子在水溶液中发生“C-S键断键”的反应是否的确为单分子单步基元反应,而不分多步发生、没有溶剂水分子或氢氧根等参与?如果是的话,是均裂为F3CSO2NSO2·自由基和·CF3自由基,还是异裂为阴阳离子?确认不是论坛里讨论极多的无垒反应么?
√546=23.36664289109

87

帖子

0

威望

1101

eV
积分
1188

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-10-8 03:42:47 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-10-4 13:28
NTf2阴离子在水溶液中发生“C-S键断键”的反应是否的确为单分子单步基元反应,而不分多步发生、没有溶剂水 ...

谢谢老师的回复,学生参考博文474《详谈使用Gaussian做势能面扫描》,先断键至3.5处 guess=mix stable=opt ub3lyp/6-311+g(d,p) scrf=(solvent=water) nosymm em=gd3bj integral=ultrafine
获得稳定波函数 然后关键词为opt=modredundant pop=always guess=(mix,always) ub3lyp/6-311+g(d,p) scrf=(solvent=water) nosymm em=gd3bj integral=ultrafine 进行了柔性扫描 所得曲线如图 请问老师能否通过柔性扫描的结果判断为此条件下断键没有能垒?

202510080330312560..png (50.8 KB, 下载次数 Times of downloads: 0)

柔性扫描

柔性扫描

1209

帖子

6

威望

2442

eV
积分
3771

Level 5 (御坂)

傻傻的木瓜

4#
发表于 Post on 2025-10-8 11:28:51 | 只看该作者 Only view this author
xxtnenu 发表于 2025-10-8 03:42
谢谢老师的回复,学生参考博文474《详谈使用Gaussian做势能面扫描》,先断键至3.5处 guess=mix stable=op ...

这个趋势是有点像,不过扫描可以再精细些(一方面缩短点步长,一方面也调调上下限,毕竟短于平衡键长的区域没必要扫),此外最好把各点提取出来另外做振动分析获得热力学校正量、算高精度单点能(记得检验各点对称破缺态的波函数稳定性),最后拿自由能曲线说事。这种处理方法也在别的帖子里讨论过:

http://bbs.keinsci.com/thread-25056-1-1.html
http://bbs.keinsci.com/thread-26668-1-1.html

√546=23.36664289109

87

帖子

0

威望

1101

eV
积分
1188

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-10-8 21:42:56 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-10-8 11:28
这个趋势是有点像,不过扫描可以再精细些(一方面缩短点步长,一方面也调调上下限,毕竟短于平衡键长的区 ...

谢谢老师的解答和建议,学生再进一步尝试

本版积分规则 Credits rule

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

GMT+8, 2026-1-25 11:18 , Processed in 0.492470 second(s), 23 queries , Gzip On.

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