24小时热门版块排行榜    

查看: 3590  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[基金申请] 固定端突然变了,今天 +3 archvillain 2026-08-06 6/300 2026-08-06 20:11 by olivermiaoer
[基金申请] 听说今天filecode变了 +12 布布和一二 2026-08-06 25/1250 2026-08-06 19:37 by wdc28
[基金申请] 基金中了 +4 laoda193707 2026-08-06 4/200 2026-08-06 18:56 by olivermiaoer
[基金申请] filecode +7 布布和一二 2026-08-06 10/500 2026-08-06 16:22 by janehuiminli
[基金申请] filecode +8 等待解的谜 2026-08-06 12/600 2026-08-06 16:11 by janehuiminli
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见