24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ellemiu

铁杆木虫 (正式写手)

应助: 3 (幼儿园)
金币: 5123
散金: 1596
红花: 2
帖子: 701
在线: 309.6小时
虫号: 923559
注册: 2009-12-09
性别: GG
专业: 化学反应工程

[求助] IRC 结果分析

TS寻找过渡态，完了以后用IRC扫描势能面，顺利收敛，，在分析结果的时候却意外的发现，irc的势能面扫描过程中，开始点和结束点的几何构型几乎一样，反应产物，感觉只是键在来回震荡。

好吧，问题出来了
1、这样的结果说明什么？如果这是一个正常的结果，是不是意味着我的反应不能发生

2、如果不是正常的结果，，改如何避免这样的情况呢？
拜请各位大神不吝赐教。

本人计算的是一个是H原子的迁移过程，涉及自由基。
IRC输入文件如下：
%chk=h3tsirc2.chk
%nprocshared=4
%mem=4GB
# b3lyp/6-31g scfcyc=300 irc=(calcfc,maxpoints=20)

TS-Ph-H-OC-Ph

0 2

第一张，irc扫描完以后的结果

irc扫描完以后的结果

回复此楼

» 猜你喜欢

中国科大电池方向任晓迪课题组招收2026级博士生-电解液/电池安全性/人工智能方向已经有20人回复
26年秋季博士申请已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有241人回复
推荐给英语教学者的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有58人回复
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有25人回复
如何掌握一门外语？已经有0人回复
企业为什么一定要做“专精特新”认定？如何成为专精特新企业？已经有5人回复
求助ASIQwin新版本的软件已经有0人回复
核磁分析软件MestReNova打开文件时报错已经有0人回复

业精于勤荒于嬉，成于思而毁于随

1楼 2012-05-18 14:59:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

necturnal

铜虫 (初入文坛)

应助: 4 (幼儿园)
金币: 266.5
红花: 1
帖子: 38
在线: 39.6小时
虫号: 534084
注册: 2008-03-27
专业: 化学动力学

我提取了你提供的结构，然后优化了一下，做出了IRC，正常结束。可以很清楚地从我上传的图片上看到H自由基的转移过程。这个IRC还是比较标准的。
你的问题是，你给我的过渡态结构，没有完全converged，距离稳定的过渡态结构，还是有微小的结构差别的。我提取了你给我的结构，在你的所用的基组和方法下，走了7步之后，才完全converged了。这可能是就是你为什么做IRC，总走到一个同一个结构的原因，因为你提供的结构没有收敛，所用IRC在找force的方向的时候，总会偏向于你的过渡态所在的势能面鞍点的方向。
下边是我的input的route section还有优化好的结构,for you refence
#B3LYP/6-31G  irc(maxpoints=14,calcfc)

CHECK

0 2
C                -2.63111700 -0.70748700 0.56895500
C                -2.97700400 0.42177700 -0.19740600
C                -4.18867800 0.45474300 -0.92437400
C                -5.04835900 -0.63912900 -0.87442100
C                -4.70445400 -1.76609100 -0.10952400
C                -3.50058500 -1.80090500 0.60443000
H                -1.68296300 -0.73042600 1.09608400
H                -4.43441600 1.33479400 -1.50904000
H                -5.98215900 -0.62085900 -1.42632500
H                -5.37748500 -2.61700800 -0.07481100
H                -3.23705600 -2.67754400 1.18566100
C                -2.07578400 1.57485900 -0.23763000
O                -2.15712500 2.58760900 -0.91821100
H                -1.06959100 1.50288200 0.65255700
C                4.94303000 0.75826500 -0.43868100
C                3.63962900 0.92450800 0.03522000
C                2.78501900 -0.18898700 0.12927900
C                3.24515100 -1.46246100 -0.25093100
C                4.54860200 -1.62076400 -0.72284200
C                5.39878500 -0.51083000 -0.81800400
H                5.60270300 1.61654500 -0.51072800
H                3.28158800 1.90214900 0.33244300
H                2.57117500 -2.30744100 -0.16795600
H                4.90192000 -2.60395800 -1.01529700
H                6.41266300 -0.63478500 -1.18519600
C                1.39365900 -0.08144600 0.62397600
O                0.61239400 -1.02219600 0.78588700
O                1.05326900 1.25529300 0.86960100
O                -0.25275500 1.39171200 1.54951400