应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » IRC 结果分析
162/2返回列表上一页12
查看: 1993  |  回复: 30
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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扫描完以后的结果
回复此楼

» 猜你喜欢
业精于勤荒于嬉，成于思而毁于随
1楼 2012-05-18 14:59:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)
引用回帖:
21楼: Originally posted by dxjdmy at 2012-05-24 14:50:21:
是的，将IRC完了的最后的点，继续优化，看看有什么结果，如果走到反应物与产物了，过渡态就是正确的，如果走不到，过渡态就不正确。

前几天IRC完了，，却出现下面的提示
This type of calculation cannot be archived.


MEMORIES ARE LIKE AN ENGLISH GRAMMER LESSON -
PRESENT TENSE, AND PAST PERFECT.
Job cpu time:  0 days 11 hours 22 minutes 23.3 seconds.
File lengths (MBytes):  RWF=    222 Int=      0 D2E=      0 Chk=     21 Scr=      1
Normal termination of Gaussian 03 at Tue May 22 11:50:59 2012.

我采用的方法是下面的这个
# b3lyp/6-31g scfcyc=300 irc=(calcfc,maxpoints=50,reverse)

necturnal给的建议，现在还没算完，不过总是在第13步死掉，提示l9999错误，加了iop(5/13=1)仍然有。
赞一下 回复此楼
业精于勤荒于嬉，成于思而毁于随
23楼2012-05-25 00:34:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)
引用回帖:
25楼: Originally posted by necturnal at 2012-05-25 08:04:10:
现在看起来，这个情况比较复杂。按理说你的这个过渡态的虚频这么大，IRC应该很容易就找到正确的force方向。
如果你不介意的话，可以把你优化好的过渡态的结构贴出来或者站短我，我看你这个基组不是很大，我可以给

昨天的结果算完了
# b3lyp/6-31g scfcyc=500 iop(5/13=1) irc=(calcfc,maxpoints=100,reverse  ,stepsize=2)

最后也是出现同样的提示
This type of calculation cannot be archived.


YOU KNOW YOU'RE A TEACHER WHEN YOU SAY 2, WRITE 3, AND MEAN 4.
-- RONALD ANSTROM, HIGH SCHOOL TEACHER, UNDERWOOD, N.D. 1974
Job cpu time:  1 days  2 hours 26 minutes 21.8 seconds.
File lengths (MBytes):  RWF=    222 Int=      0 D2E=      0 Chk=     21 Scr=      1
Normal termination of Gaussian 03 at Thu May 24 16:41:53 2012.
等下我把TS的结构贴在下面
赞一下 回复此楼
业精于勤荒于嬉，成于思而毁于随
26楼2012-05-25 09:15:27
已阅   回复此楼   关注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的回帖

ellemiu

铁杆木虫 (正式写手)
引用回帖:
28楼: Originally posted by necturnal at 2012-05-25 15:07:56
我提取了你提供的结构，然后优化了一下，做出了IRC，正常结束。可以很清楚地从我上传的图片上看到H自由基的转移过程。这个IRC还是比较标准的。
你的问题是，你给我的过渡态结构，没有完全converged，距离稳定的过渡 ...

十分感谢necturnal的热心回复。

我给你的过渡态是在我这里计算后收敛的结果，，不知为啥在你那边看没有converged，感觉很奇怪。刚仔细检查了下之前TS的结果，收敛的时候没有出现异常。

我把导出的结果再算一次看看。

TS的结果
赞一下 回复此楼
业精于勤荒于嬉，成于思而毁于随
29楼2012-05-26 10:15:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ellemiu

铁杆木虫 (正式写手)
引用回帖:
28楼: Originally posted by necturnal at 2012-05-25 15:07:56
我提取了你提供的结构，然后优化了一下，做出了IRC，正常结束。可以很清楚地从我上传的图片上看到H自由基的转移过程。这个IRC还是比较标准的。
你的问题是，你给我的过渡态结构，没有完全converged，距离稳定的过渡 ...

你说的这个问题，我已经找到原因了，是因为TS和IRC用的基组不一样造成的，我在ts的是时候用的是6-31g(d,p)，irc的时候在，只是6-31g。所以会造成ts的结果不同。

后来我用6-31g对6-31g(d,p)方法的TS结果再ts了一次，发现大概走了8步才收敛。将得到的结果另存为输入文件后，再走irc，走了两天却还没有收敛，很奇怪
赞一下 回复此楼
业精于勤荒于嬉，成于思而毁于随
30楼2012-06-01 09:17:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ellemiu 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭