24小时热门版块排行榜    

查看: 1223  |  回复: 6

ct_kaoyan

金虫 (小有名气)

[求助] 找一个互变异构反应过渡态结构 已有1人参与

我现要要计算下面反应的过渡态,我分别用TS,QST2和QST3方法找过渡态。找到两种过渡态ts1和ts2,不不知道正确与否,ts1用gview查看负虚频震动方向是H原子在20号碳原子和21号氧原子之间震动,如下图ts1;ts2是用gview查看负虚频震动方向氢原子在20号碳原子和8号碳原子之间震动,如下图ts2。请问哪种正确?或者两种都不正确,正确的过渡态的氢原子是不是应该在碳8和氧21之间震动?请有经验的量化高手解答,非常感谢!!

找一个互变异构反应过渡态结构
反应.png


找一个互变异构反应过渡态结构-1
ts1.png


找一个互变异构反应过渡态结构-2
ts2.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

每天进步一点点!
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ct_kaoyan

金虫 (小有名气)

怎么没人呢?自己顶一下!
每天进步一点点!
2楼2015-07-03 21:18:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

oinkmasta

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
ct_kaoyan: 金币+50, ★★★很有帮助, 谢谢你的答复! 2015-07-08 20:10:15
两个看起来都似乎是反应路径中的一环,所以可以试试以下几步:
1)如果ts1和ts2各自虚频都正常(有且仅有一负),可以各做irc,看看是不是ts1的反应物向和ts2的生成物向有结构一致。如果有,那就是中间态了。
2)如果irc找不到一致性,或者说很可能没有中间态,则ts1对/ts2对/ts1,2都不对。
3)假设ts1和ts2有一个正确,只要比较能量取高即可。
4)都不对也有可能?那就需要柔性扫描了,扫H和8号C的距离,或者扫H和O的距离,应该都可以。有条件也可以尝试二维柔性扫描。
16届本科在北美
3楼2015-07-05 00:17:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ct_kaoyan

金虫 (小有名气)

引用回帖:
3楼: Originally posted by oinkmasta at 2015-07-05 00:17:57
两个看起来都似乎是反应路径中的一环,所以可以试试以下几步:
1)如果ts1和ts2各自虚频都正常(有且仅有一负),可以各做irc,看看是不是ts1的反应物向和ts2的生成物向有结构一致。如果有,那就是中间态了。
2) ...

首先谢谢你的回复。按照你的方法对两过渡态各做了irc,可是得不到ts1的反应物向和ts2的生成物,总是出错。我想可能是两个过渡态都不对,想尝试重新寻找。
每天进步一点点!
4楼2015-07-08 20:09:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

oinkmasta

金虫 (小有名气)

【答案】应助回帖

引用回帖:
4楼: Originally posted by ct_kaoyan at 2015-07-08 07:09:34
首先谢谢你的回复。按照你的方法对两过渡态各做了irc,可是得不到ts1的反应物向和ts2的生成物,总是出错。我想可能是两个过渡态都不对,想尝试重新寻找。...

出错是指高斯不能完成irc?还是得到的构型不对?你算出的虚频还挺合理的,我猜不应该是后者呀。
如果是高斯报错,尝试把rcfc换成calcfc,重新计算force constants。我最近从ts虚频读到irc的时候也碰到过这个情况。

下面是我一般默认选取的irc命令,供参考用。步长和步数当然都可以调节。
# irc=(calcfc,stepsize=5,maxpoints=30)
16届本科在北美
5楼2015-07-08 21:15:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

oinkmasta

金虫 (小有名气)

更正一下,上面我写的第三步不正确,它对应的是有中间态的情况。
没有中间态的时候比较ts1和ts2,应当看哪一个的irc能连接反应物和生成物。
16届本科在北美
6楼2015-07-08 22:53:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ct_kaoyan

金虫 (小有名气)

引用回帖:
5楼: Originally posted by oinkmasta at 2015-07-08 21:15:26
出错是指高斯不能完成irc?还是得到的构型不对?你算出的虚频还挺合理的,我猜不应该是后者呀。
如果是高斯报错,尝试把rcfc换成calcfc,重新计算force constants。我最近从ts虚频读到irc的时候也碰到过这个情况。 ...

不好意识,这几天出差没有及时回复。出错是指不能完成IRC,会出现“Problem with the distance matrix.”,'l202.exe'错误或者出现“ The SCF is confused.”,l502.exe错误。我用的irc命令是“#p b3lyp/6-31g* irc=(maxpoints=50,rcfc,forward,internal) iop(2/16=1) iop(5/13=1)”,先正在尝试把rcfc换成calcfc,重新计算force constants。
每天进步一点点!
7楼2015-07-13 20:35:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ct_kaoyan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 303求调剂 +7 元夕元 2026-03-20 8/400 2026-03-26 20:38 by 不吃魚的貓
[考研] 求调剂 +8 Auroracx 2026-03-22 8/400 2026-03-26 19:55 by 不吃魚的貓
[考研] 266分求材料化工冶金矿业等专业的调剂 +3 哇呼哼呼哼 2026-03-26 3/150 2026-03-26 19:16 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工(085600)296求调剂 +9 稻妻小编 2026-03-26 9/450 2026-03-26 16:16 by 不吃魚的貓
[考研] 资源与环境 调剂申请(333分) +9 holy J 2026-03-21 9/450 2026-03-26 15:47 by 161765490
[考研] 材料与化工考研调剂 +10 孅華 2026-03-22 10/500 2026-03-26 15:40 by zzll406
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 材料学硕,求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 289材料与化工(085600)B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
信息提示
请填处理意见