24小时热门版块排行榜    

查看: 3642  |  回复: 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的回帖

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的回帖

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的回帖

yypcnbe

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
gkf高: 金币+2 2013-09-28 09:36:54
优化过渡态的时候得加opt=noeigentest,不然会出现l9999错误,具体看手册
4楼2013-09-27 17:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 play900325 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:03 by LH5NkK5Ud4bh
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:27 by DgNGHc3h5tPl
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:11 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 03:59 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-16 02:38 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +6 两三岁ss 2026-09-14 6/300 2026-09-15 20:17 by 只吃肉葡萄
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
信息提示
请填处理意见