24小时热门版块排行榜    

查看: 1977  |  回复: 30
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ellemiu

铁杆木虫 (正式写手)

[求助] 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扫描完以后的结果
回复此楼
业精于勤荒于嬉,成于思而毁于随
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)

%chk=h3tsirc2.chk
%nprocshared=4
%mem=4GB
# b3lyp/6-31g scfcyc=300 irc=(calcfc,maxpoints=20) iop(5/13=1)

TS-PhOOO-H-OC-Ph

0 2
C              
C                  1            B1
C                  2            B2    1            A1
C                  3            B3    2            A2    1            D1    0
C                  4            B4    3            A3    2            D2    0
C                  1            B5    2            A4    3            D3    0
H                  1            B6    6            A5    5            D4    0
H                  3            B7    2            A6    1            D5    0
H                  4            B8    3            A7    2            D6    0
H                  5            B9    4            A8    3            D7    0
H                  6           B10    1            A9    2            D8    0
C                  2           B11    1           A10    6            D9    0
O                 12           B12    2           A11    1           D10    0
H                 12           B13    2           A12    1           D11    0
C                 12           B14    2           A13    1           D12    0
C                 15           B15   12           A14    2           D13    0
C                 16           B16   15           A15   12           D14    0
C                 17           B17   16           A16   15           D15    0
C                 18           B18   17           A17   16           D16    0
C                 15           B19   12           A18    2           D17    0
H                 15           B20   12           A19    2           D18    0
H                 16           B21   15           A20   12           D19    0
H                 18           B22   17           A21   16           D20    0
H                 19           B23   18           A22   17           D21    0
H                 20           B24   15           A23   12           D22    0
C                 17           B25   16           A24   15           D23    0
O                 26           B26   17           A25   16           D24    0
O                 26           B27   17           A26   16           D25    0
O                 28           B28   26           A27   17           D26    0

   B1             1.40248139
   B2             1.40713742
   B3             1.38886569
   B4             1.40068288
   B5             1.39408101
   B6             1.08550242
   B7             1.08508778
   B8             1.08573132
   B9             1.08625634
   B10            1.08532167
   B11            1.47413846
   B12            1.20211493
   B13            1.27182204
   B14            6.93119797
   B15            1.39328791
   B16            1.40317644
   B17            1.40214419
   B18            1.39168817
   B19            1.39664637
   B20            1.08569148
   B21            1.08351812
   B22            1.08448852
   B23            1.08572741
   B24            1.08614597
   B25            1.48831514
   B26            1.20907486
   B27            1.37983752
   B28            1.40140820
   A1           120.84943922
   A2           119.39420830
   A3           119.88994249
   A4           119.17402935
   A5           121.24181798
   A6           119.21372297
   A7           120.03664471
   A8           119.60455464
   A9           119.78359955
   A10          119.56041105
   A11          127.62224470
   A12          113.30495416
   A13          124.16221422
   A14           19.37602150
   A15          119.73205596
   A16          119.90733064
   A17          120.06475911
   A18          115.54246673
   A19          120.76261940
   A20          120.37715344
   A21          118.60469132
   A22          119.90305320
   A23          119.91500428
   A24          122.84266156
   A25          125.52797146
   A26          110.48184705
   A27          112.78672513
   D1             0.32148514
   D2            -0.29782655
   D3             0.00359667
   D4          -178.60865205
   D5          -179.91132619
   D6           179.76067412
   D7          -179.90730427
   D8           179.55053910
   D9           179.84334744
   D10          173.98818884
   D11           -8.40186142
   D12           41.65658969
   D13          -89.94041521
   D14           81.25330919
   D15            0.07592745
   D16           -0.06612855
   D17           18.83956755
   D18          177.40618261
   D19          -98.76043187
   D20          179.90387643
   D21          179.99714499
   D22          158.69913609
   D23          179.96091803
   D24         -178.34508059
   D25            2.52988933
   D26         -174.30365111


这个是算好的TS的分子信息,我也是被弄的有点焦头烂额,找不到问题所在。necturnal有兴趣可以算下。一起学习。
业精于勤荒于嬉,成于思而毁于随
27楼2012-05-25 09:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 31 个回答

meteoric30

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ellemiu: 金币+2, 有帮助, irc直接拿的TS的结果来继续的,不太明白irc的产物是指什么,你指的是最后一步的状态进行优化吗》?步长设大一点怎样设呢? 2012-05-18 16:19:23
步长设大一点,irc的产物优化了吗?优化的结果看看是不是和ts一样
2楼2012-05-18 15:35:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by meteoric30 at 2012-05-18 15:35:18:
步长设大一点,irc的产物优化了吗?优化的结果看看是不是和ts一样

irc直接拿的TS的结果来继续的,不太明白irc的产物是指什么,你指的是最后一步的状态进行优化吗?

还有步长设大一点怎样设呢?
业精于勤荒于嬉,成于思而毁于随
3楼2012-05-18 16:20:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldzh0531

木虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
ellemiu: 金币+1, 有帮助, 我算出的虚频为-1415,其实我只是想扫势能面看看TS是不是正确的 2012-05-21 10:14:58
IRC只是根据频率计算单点,并不优化。想要得到准确构型要提取出终端点优化。另外,通常氢转移的虚频要大于1000
4楼2012-05-18 17:58:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见