应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gaussian优化l103错误求助
北京石油化工学院2026年研究生招生接收调剂公告
41/1返回列表
查看: 3440  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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错误好像也是输入分子坐标有问题，我现在很困惑我应该怎么修改输入分子？或者是不是应该加其他的一些关键字？
小女子感激不尽！！！！
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-26 21:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luckone2013

银虫 (小有名气)

有情怀的伟大科学家

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2013-09-26 23:20:36
l999错误是优化圈数不够，把out文件保存成gjf，修改后接着优化。这样比较省事。
或者在原来输入文件opt中加入maxcyc=500,把优化圈数加大到500
赞一下(1人) 回复此楼
总得有点奔头！@！@！
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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 288求调剂，一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 一志愿上海大学生物学346 +3 上海大学346调剂 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[考研] 11408，284分，二战真诚求调剂 +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 0856材料与化工调剂，339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 学硕化学工程与技术，一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 348求调剂 +6 吴彦祖24k 2026-04-02 6/300 2026-04-02 14:07 by 给你你注意休息
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by 无际的草原
[考研] 求调剂 +4 图鉴212 2026-03-30 5/250 2026-04-01 15:32 by 图鉴212
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 一志愿大连理工大学材料求调剂 +6 Gymno 2026-03-30 6/300 2026-03-31 07:26 by 无际的草原
信息提示
关闭
请填处理意见
关闭