24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2511  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见