24小时热门版块排行榜    

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

play900325

铁虫 (小有名气)

[求助] gaussian优化l103错误求助

各位高人,我现在正在学习用qst2,qst3还有ts找过渡态的方法,以下是我用于qst2文件的一个输入分子,想先把它优化一下,输入文件是:
%nprocshared=2
%mem=200MW
%chk=product.chk
#p opt=(gdiis,cartesian) b3lyp/6-31g geom=connectivity

qst2

1 1
C                 -1.24478188   -0.99808527    0.66678566
H                 -0.65028357   -1.67056742    0.08434294
H                 -0.65039672   -0.15738417    0.95803803
H                 -2.07701494   -0.66185324    0.08434581
Br                -0.22816152    0.43983571   -1.82346861
Cl                -1.83142182   -1.82783639    2.10378484

1 2 1.0 3 1.0 4 1.0 6 1.0
2
3
4
5
6
我的分子中有一个一氯甲烷,旁边是一个溴离子,这个是产物分子;还有一个反应物的分子是一溴甲烷和氯离子。上述输入文件计算结果的末尾部分:
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GDIIS optimizer.
Eigenvalue    13 is   1.05D-04 should be less than     0.000000 Eigenvector:
                          A11       A5        A6        A10       R7
   1                    0.64185  -0.41187   0.40853  -0.30067  -0.22819
                          R8        D3        D2        A2        D6
   1                    0.20959  -0.12634  -0.12190  -0.09768   0.08561
Eigenvalue    14 is   3.01D-07 should be less than     0.000000 Eigenvector:
                          A10       A3        A11       R6        A6
   1                    0.68522  -0.42654   0.34444  -0.23518   0.20748
                          A5        D4        R8        A4        R7
   1                    0.20372   0.14807   0.11267  -0.10801   0.10775
NTrRot=    -1 NTRed=    26 NAtoms=     6 NSkip=    14 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /apps/g09/l103.exe at Thu Sep 26 21:07:38 2013.
Job cpu time:  0 days  0 hours  0 minutes 19.5 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1
我查了一下,l103错误一般是输入分子的坐标有问题,将Cl-C-Br之间夹角稍作改变(原来是180度)会出现l9999错误:
Mulliken charges with hydrogens summed into heavy atoms:
              1
     1  C    0.069956
     5  Br   0.632050
     6  Cl   0.297994
Sum of Mulliken charges with hydrogens summed into heavy atoms =   1.00000
Electronic spatial extent (au):  <R**2>=            583.9483
Charge=              1.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=             -0.5497    Y=              0.8657    Z=              1.3687  Tot=              1.7103
Quadrupole moment (field-independent basis, Debye-Ang):
   XX=            -29.4137   YY=            -29.7419   ZZ=            -33.3349
   XY=              2.9438   XZ=             -0.8125   YZ=             -1.9354
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
   XX=              1.4165   YY=              1.0883   ZZ=             -2.5047
   XY=              2.9438   XZ=             -0.8125   YZ=             -1.9354
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=            -25.2953  YYY=             -0.6801  ZZZ=              1.0026  XYY=             -8.5894
  XXY=             -5.0319  XXZ=             -0.1328  XZZ=            -10.4852  YZZ=              0.8477
  YYZ=              0.9751  XYZ=              2.5991
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=           -485.9089 YYYY=            -73.9002 ZZZZ=            -43.1648 XXXY=             -5.8307
XXXZ=             -0.5355 YYYX=              4.7087 YYYZ=             -3.6091 ZZZX=              0.5800
ZZZY=             -1.3081 XXYY=            -92.6112 XXZZ=            -94.7711 YYZZ=            -16.1956
XXYZ=             -1.9851 YYXZ=             -1.3529 ZZXY=              3.2027
N-N= 2.219478992611D+02 E-N=-7.753497419430D+03  KE= 3.059649242168D+03
No NMR shielding tensors so no spin-rotation constants.
Leave Link  601 at Thu Sep 26 21:23:24 2013, MaxMem=  209715200 cpu:       0.3
(Enter /apps/g09/l9999.exe)


HE THAT FOLLOWS NATURE IS NEVER LOST.

                             -- BACON
Error termination request processed by link 9999.
Error termination via Lnk1e in /apps/g09/l9999.exe at Thu Sep 26 21:23:25 2013.
Job cpu time:  0 days  0 hours  4 minutes 47.5 seconds.
File lengths (MBytes):  RWF=      5 Int=      0 D2E=      0 Chk=      1 Scr=      1

而l9999错误好像也是输入分子坐标有问题,我现在很困惑我应该怎么修改输入分子?或者是不是应该加其他的一些关键字?
小女子感激不尽!!!!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

play900325

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by luckone2013 at 2013-09-26 22:56:40
l999错误是优化圈数不够,把out文件保存成gjf,修改后接着优化。这样比较省事。
或者在原来输入文件opt中加入maxcyc=500,把优化圈数加大到500

谢谢这位师兄~计算是算完了,但是现在有另一个问题,优化出来的结果是一个溴化氢和氯甲基离子,这跟我要输入qst2的分子差别有点大啊,这样的分子输入进去以后会找不到正确的过渡态吧~~这种情况您遇到过吗?怎么解决呢?
3楼2013-09-27 10:21:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

luckone2013

银虫 (小有名气)

有情怀的伟大科学家

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
gmy1990: 金币+3 2013-09-26 23:20:36
l999错误是优化圈数不够,把out文件保存成gjf,修改后接着优化。这样比较省事。
或者在原来输入文件opt中加入maxcyc=500,把优化圈数加大到500
总得有点奔头!@!@!
2楼2013-09-26 22:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yypcnbe

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2013-09-28 09:36:54
优化过渡态的时候得加opt=noeigentest,不然会出现l9999错误,具体看手册
4楼2013-09-27 17:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见