24小时热门版块排行榜    

查看: 3600  |  回复: 6

ludeng8710

新虫 (小有名气)

[求助] gaussian09过渡态优化出现l103错误,求解决

gaussian09过渡态优化出现l103错误,这是在找到过渡态,后面算完频率的l1110和l1002之后出现的错误。求解决,或者说我可以提取这个输出文件的吉布斯自由能作为过渡态的能量从而求得反应能垒么?谢谢

%mem=2000mb
%rwf=1,2000mb,2,2000mb,3,2000mb,4,2000mb,5,2000mb,6,2000mb,7,-1
%nosave
%nprocshared=4
#p opt=(Cartesian,calcfc,ts,noeigen,maxcyc=400,loose) freq=noraman scf=(maxcyc=500) ub3lyp/6-31g iop(1/8=1)

Title Card Required

0 1
C                  0.70633000    2.73491100   -1.52012200
C                  0.42185800    3.68426600    0.35426500
C                  2.62543500    3.24215900    1.33093200
C                  3.53396200    2.53624600    2.08422400
C                  3.13211000    1.36304800    2.78728100
C                  1.75952500    0.98489700    2.66692600
....

输出文件:

   104        9           0.000008645    0.000034546   -0.000152872
    105        9           0.000227423    0.000071027    0.000080418
    106        9          -0.000023744    0.000053105    0.000042804
-------------------------------------------------------------------
Cartesian Forces:  Max     0.002460157 RMS     0.000380694
Leave Link  716 at Sun Jul 15 21:01:15 2012, MaxMem=  262144000 cpu:       5.8
(Enter /home/gan/g09/l103.exe)

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Eigenvalue   312 is   0.00D+00 should be greater than     0.000001 Eigenvector:
                          R8        R2        R29       R1        A72
   1                    0.37495   0.33425  -0.30739  -0.27016   0.26566
                          A73      A159       R7        A6        A64
   1                    0.26021  -0.22688  -0.21478   0.19047  -0.17755
NTrRot=    -1 NTRed=   559 NAtoms=   106 NSkip=   247 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /home/gan/g09/l103.exe at Sun Jul 15 21:01:16 2012.
Job cpu time:  1 days 17 hours 17 minutes 23.8 seconds.
File lengths (MBytes):  RWF=  11795 Int=      0 D2E=      0 Chk=     41 Scr=      1

求解决
谢谢~

[ Last edited by ludeng8710 on 2012-7-15 at 10:11 ]
回复此楼

» 猜你喜欢

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

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

左手疗伤

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
103可能是程序找不到minima,你限定的步长是不是太小了,去掉iop试试
非淡泊无以明志,非宁静无以致远
2楼2012-07-15 10:54:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
听说这个是安装盗版g09的一个bug,我也没遇过.
Error in internal coordinate system,按道理加上cartesian关键字就可以避免那个bug了,不懂.
3楼2012-07-15 11:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

astringent

铜虫 (著名写手)

帮顶并同问,我在优化时也经常遇到这种情况,不知道怎么解决。请高手解答
4楼2012-07-15 17:37:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZZU2011

木虫 (正式写手)

【答案】应助回帖

换成z-matrix试一下,不行的话,就得调整结构了,建议楼主先用低基组优化一下,然后再升高基组
5楼2012-09-01 10:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaogaiyun

铜虫 (小有名气)

【答案】应助回帖

我觉得应该用QST3来算,就是再加一个你的猜测结构
6楼2013-04-19 10:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marson

银虫 (正式写手)

【答案】应助回帖

我曾经遇到过这样的问题。两种解决方法:
通过G03版本优化去解决,可能是G09的一个bug;
局部稍微改动构型,这种内坐标错误一般是因为构型中的一些二面角之类的在优化过程中超出G09识别了,我的是179°然后再增加的话就到181°,这样一般G09就报错了~建议还是改动局部构型吧~
思想有多远,我们就能走多远~
7楼2013-04-19 15:16:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ludeng8710 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-31 02:11 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-31 01:16 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 01:07 by hZiFeudZyoGR
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-31 00:56 by hZiFeudZyoGR
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-31 00:42 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 15/750 2026-08-31 00:10 by hZiFeudZyoGR
[基金申请] 国社科又开始会评了,不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 4FFAWE8HcgUD 2026-08-29 5/250 2026-08-30 20:41 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 4FFAWE8HcgUD 2026-08-29 3/150 2026-08-30 20:10 by hZiFeudZyoGR
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 18:28 by hZiFeudZyoGR
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 18:17 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-30 17:56 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
信息提示
请填处理意见