24小时热门版块排行榜    

查看: 772  |  回复: 7
当前主题已经存档。

baiyang2000

银虫 (正式写手)

[交流] 【求助】求教gaussian过渡态优化问题

最近刚开始研究gaussian,在做gaussian说明书第八章的例题的时候(用程序自带的example),碰到很棘手的问题,当例题中用opt=ts寻找过渡态的时候,每次都出现错误,提示如下:
Opt=TS requires force constants or ModRedundant.
Error termination via Lnk1e in F:\G03W\l1.exe at Sat Jul 25 02:06:12 2009.

但是我加了CalcFC.或者ModRedundant,无任何作用
在厦大bbs有人问类似的问题,g98可以做,但是g03就有类似错误,但是没有人给出解答,例题应该没有错误吧
能有高手分析分析原因么?

注:我目前没有用工作站,只是在用自己的笔记本试做,笔记本的配置是内存2g,硬盘120g,amd turion64,gaussian自带的例题都是非常小的分子,是否电脑的配置问题有影响?

例题 e8_03.gif:
%Chk=formts
#T RHF/6-31G(d) Opt(TS,calcFc,Noeigen) Freq Test

H2 + CO <-> H2CO Opt Freq

0,1
O
C,1,1.13
H,2,1.1,1,164.0
H,3,1.3,2,90.0,1,0.

--Link1--
%Chk=formts
%NoSave
#T RHF/6-31G(d) IRC=RCFC Guess=Read Geom=AllCheck Test

例题,e_8.10:是研究sn2亲核取代反应的,第一步是寻找过渡态,第二步和第三步都是寻找反应中的临时产物离子中性分子混合物,第二,三步单独做,没有问题,有问题的就是第一步
%Chk=sn2
#T RHF/6-31+G(d) Opt=(TS,CalcFC,Noeigentest) Freq Test

Cl...CH3...F TS Opt

  -1   1
C
Cl,1,RCCl
H,1,RCH,2,A3
H,1,RCH,2,A3,3,120.
H,1,RCH,2,A3,3,-120.
X,1,1.,2,90.,3,180.
F,1,R6,6,90.,2,180.
      Variables:
RCCl=2.133
RCH=1.062
R6=2.126
A3=97.3

--Link1--
%Chk=sn2
%NoSave
#T RHF/6-31+G(D) IRC(RCFC,StepSize=20) Geom=AllCheck Guess=Read Test

--Link1--
#T RHF/6-31+G(D) Opt=CalcFC Freq Test

CH3F...Cl Minimum

  -1   1
C
F,1,R2
H,1,R3,2,A3
H,1,R3,2,A3,3,120.
H,1,R3,2,A3,3,-120.
X,1,1.,2,90.,3,180.
Cl,1,R4,6,90.,2,180.
      Variables:
R2=1.398
R3=1.077
R4=3.422
A3=108.4

--Link1--
#T RHF/6-31+G(D) Opt=CalcFC Freq Test

CH3Cl...F Minimum

  -1   1
C
Cl,1,RCCl
H,1,RCH,2,A3
H,1,RCH,2,A3,3,120.
H,1,RCH,2,A3,3,-120.
X,1,1.,2,90.,3,180.
F,1,R6,6,90.,2,180.
      Variables:
RCCl=1.863
RCH=1.07
R6=2.585
A3=107.0

[ Last edited by baiyang2000 on 2009-7-25 at 23:10 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★ ★
baiyang2000(金币+1,VIP+0): 7-25 18:06
lei0736(金币+2,VIP+0):谢谢 7-27 09:54
引用回帖:
Originally posted by baiyang2000 at 2009-7-25 17:53:
最近刚开始研究gaussian,在做gaussian说明书第八章的例题的时候(用程序自带的example),碰到很棘手的问题,当例题中用opt=ts寻找过渡态的时候,每次都出现错误,提示如下:
Opt=TS requires force constants or ModRe ...

试试 opt=(ts,calcfc,noeigentest)
Simpler, stronger, more tolerant
2楼2009-07-25 18:00:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baiyang2000

银虫 (正式写手)

照楼上的方法做了,还是一样的错误提示,难道gaussian的例题有错误?
3楼2009-07-25 18:06:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)


lei0736(金币+1,VIP+0):谢谢 7-27 09:55
引用回帖:
Originally posted by baiyang2000 at 2009-7-25 18:06:
照楼上的方法做了,还是一样的错误提示,难道gaussian的例题有错误?

将完整的输入文件贴上来,否则没法帮你。
Simpler, stronger, more tolerant
4楼2009-07-25 21:46:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

313588908

木虫 (正式写手)


lei0736(金币+1,VIP+0):谢谢 7-27 09:55
输入有误吧。仔细从自带的example检查一下。
5楼2009-07-25 22:43:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baiyang2000

银虫 (正式写手)

已经将出错的两个原文件贴上,请高手指教
6楼2009-07-25 23:11:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★ ★ ★ ★ ★ ★
baiyang2000(金币+4,VIP+0): 7-26 05:58
lei0736(金币+3,VIP+0):谢谢 7-27 09:55
我用我的笔记本(x86-Win32-G03RevB.03)算了你的这个例子:
引用回帖:
例题,e_8.10:是研究sn2亲核取代反应的,第一步是寻找过渡态,第二步和第三步都是寻找反应中的临时产物离子中性分子混合物,第二,三步单独做,没有问题,有问题的就是第一步
%Chk=sn2
#T RHF/6-31+G(d) Opt=(TS,CalcFC,Noeigentest) Freq Test

Cl...CH3...F TS Opt
…………
…………

完全没有问题,正常结束。

“ Job cpu time:  0 days  0 hours  0 minutes 15.0 seconds.
File lengths (MBytes):  RWF=     16 Int=      0 D2E=      0 Chk=      7 Scr=      1
Normal termination of Gaussian 03 at Sat Jul 25 18:14:12 2009.”


是不是你的程序有问题?
Simpler, stronger, more tolerant
7楼2009-07-26 00:17:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baiyang2000

银虫 (正式写手)

确实是程序的问题,我又重新装了程序,现在一切正常了,谢谢楼上的兄弟
8楼2009-07-26 05:58:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 baiyang2000 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见