24小时热门版块排行榜    

查看: 3393  |  回复: 22

ct_kaoyan

金虫 (小有名气)

[交流] QST2方法找到过渡态后IRC计算出现l202错误

RT,请教高手,如何解决问题,万分感谢!
过渡态是用QST2方法找到的,然后用TS方法优化了QST2找到的过渡态,可以得到唯一虚频,然后再用TS计算得到的过渡态(用gview查看与QST2找到的过渡态结构一样,能量一样)计算IRC找产物和反应物。
QST2找过渡态的基组与方法为:
#p b3lyp/6-31g* fopt(qst2,,maxcyc=150) iop(5/13=1) iop(2/16=1) iop(1/11=1)
TS优化方法与基组:
#p b3lyp/6-31g* fopt(ts,z-matrix,maxcyc=50) iop(5/13=1) iop(2/16=1) freq iop(1/11=1)


IRC-fts计算输入文件如下:
#p b3lyp/6-31g* irc=(maxpoints=50,rcfc,forward,internal) iop(2/16=1) iop(5/13=1) geom=check guess=check

Untitled I Draw Styles-f1ts2

0 1
~
~
~


IRC输出文件的最后几行如下:
   31  H    0.000000
    32  H    0.000000   0.000000
    33  H    0.000000   0.000000   0.000000
Small interatomic distances encountered:
     2    1

     3    1

     3    2

     4    1

     4    2

     4    3

     5    1

     5    2

     5    3

     5    4

     6    1

     6    2

     6    3

     6    4

     6    5

     7    1

     7    2

     7    3

     7    4

     7    5
Problem with the distance matrix.
Error termination via Lnk1e in /home/rwshi/soft/g03/l202.exe at Mon Jun 29 05:54:31 2015.
Job cpu time:  0 days  5 hours 40 minutes 11.9 seconds.
File lengths (MBytes):  RWF=    153 Int=      0 D2E=      0 Chk=     28 Scr=      2
                                                                                                 35419,2       Bot

[ Last edited by ct_kaoyan on 2015-7-1 at 17:19 ]
回复此楼
每天进步一点点!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

ct_kaoyan

金虫 (小有名气)

ct_kaoyan: 回帖置顶 2015-06-30 11:07:12
输入文件如下:
#p b3lyp/6-31g* irc=(maxpoints=50,rcfc,forward,internal) iop(2/16=1) iop(5/13=1) geom=check guess=check

Untitled I Draw Styles-f1ts2

0 1
每天进步一点点!
2楼2015-06-30 11:06:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
ct_kaoyan: 金币+4, 谢谢! 2015-07-03 14:56:42
Problem with the distance matrix.

可能是某些原子距离太近了。把最后一步的结构画出来,也许能发现出错的键长
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
3楼2015-06-30 19:57:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
ct_kaoyan(金币+1): 谢谢参与
4楼2015-07-01 16:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
ct_kaoyan(金币+1): 谢谢参与
5楼2015-07-01 16:35:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
ct_kaoyan(金币+1): 谢谢参与
6楼2015-07-01 16:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ct_kaoyan

金虫 (小有名气)

引用回帖:
3楼: Originally posted by beefly at 2015-06-30 19:57:16
Problem with the distance matrix.

可能是某些原子距离太近了。把最后一步的结构画出来,也许能发现出错的键长

谢谢你的建议,下面是最后一步的坐标,可以看出3号碳原子的Y坐标出现异常,是用星号表示的,gview没法显示结构,请问怎么解决这个问题了?
Atom  N1     Length/X     N2    Alpha/Y     N3     Beta/Z      J
------------------------------------------------------------------------
   1   1  C     0   5.238823          0.438130         -0.116319
   2   2  C     0   5.035136         -0.951357          0.620244
   3   3  C     0   3.780851       ***********          1.596051
   4   4  C     0   2.814175         -0.502438          1.767319
   5   5  C     0   2.966782          0.944417          1.060862
   6   6  C     0   4.220649          1.386053          0.147375
   7   7  N     0   1.372179         -0.846187          2.284808
   8   8  C     0   0.842129          0.415444          2.037471
   9   9  C     0   1.522117          1.800651          1.678391
  10  10  C     0   0.707451         -1.982155          1.823121  
  11  11  C     0  -0.742270         -2.008599          1.595519
  12  12  C     0  -1.408730         -0.794506          1.734006
  13  13  C     0  -0.520675          0.346235          1.816865
  14  14  O     0   0.763312          2.580430         -0.048940
  15  15  C     0  -2.737925         -0.154444          1.722602
  16  16  C     0  -2.479990          1.308359          1.710443
  17  17  N     0  -1.107791          1.559119          1.804029
  18  18  C     0  -4.111902         -0.618518          1.815910
  19  19  C     0  -5.168533          0.345767          1.903303
  20  20  C     0  -4.878766          1.767102          1.869175
  21  21  C     0  -3.544253          2.262634          1.759184
  22  22  H     0   6.224454          0.805132         -0.915509
  23  23  H     0   5.856747         -1.644553          0.408640
  24  24  H     0   3.590293         -2.595137          2.161902
  25  25  H     0   4.398083          2.481576         -0.400283
  26  26  H     0   1.273728         -2.907431          1.636708
  27  27  H     0  -1.222647         -2.989228          1.189731
  28  28  H     0  -0.328120          2.380053          1.413437
  29  29  H     0  -4.352281         -1.711350          1.870470
  30  30  H     0  -6.239101          0.012913          2.038610
  31  31  H     0  -5.726367          2.499518          1.965169
  32  32  H     0  -3.313979          3.355448          1.772589
  33  33  H     0   3.065183          3.615389          3.461502
每天进步一点点!
7楼2015-07-01 16:38:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
ct_kaoyan(金币+1): 谢谢参与
8楼2015-07-01 16:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
ct_kaoyan(金币+1): 谢谢参与
9楼2015-07-01 16:54:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
ct_kaoyan(金币+1): 谢谢参与
10楼2015-07-01 16:54:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ct_kaoyan 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +11 xqy2009 2026-10-13 14/700 2026-10-16 10:28 by goodluck7111
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+4 纳斯达克ML 2026-10-11 5/250 2026-10-16 10:03 by shjy2599
[基金申请] 申请科研项目纯靠实力能成功吗? 200+8 David1982 2026-10-13 12/600 2026-10-16 09:50 by shjy2599
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见