24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3439  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 338求调剂 +3 我想上岸ii 2026-04-05 3/150 2026-04-05 19:59 by nepu_uu
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 化学0703-一志愿211-338分求调剂 +7 vants 2026-04-05 7/350 2026-04-05 18:17 by cql1109
[考研] 0817化学工程与技术求调剂,一志愿中海洋319 +13 lv945 2026-04-04 13/650 2026-04-05 18:14 by 猪会飞
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 081700,311,求调剂 +12 冬十三 2026-04-04 13/650 2026-04-05 08:43 by qlm5820
[考研] 一志愿北京化工大学,初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 309求调剂 +6 刘刘刘1231 2026-04-02 7/350 2026-04-04 13:41 by liucky
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 334求调剂 +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by 琢珥丶
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 071000生物学调剂 +8 知昭蔓 2026-04-02 8/400 2026-04-03 10:36 by macy2011
[考研] 303求调剂 +3 一色清羽 2026-04-02 4/200 2026-04-03 10:22 by 蓝云思雨
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 学硕化学工程与技术,一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 初试301,代码085701环境工程,本硕一致,四六级已过,有二区一作,共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666!
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
信息提示
请填处理意见