24小时热门版块排行榜    

查看: 1240  |  回复: 3
当前主题已经存档。

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作

[交流] [求助]计算H3CO --> H2COH过渡态、进行IRC计算的问题

在高斯自带例子e3_03.gjf计算H3CO --> H2COH Reactants的基础上
自己在后面加上了IRC的命令

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


最终IRC优化结束,提示

Optimization completed.
    -- Optimized point #   1 Found.


然后输出优化参数,坐标、键长、键角之类的
但是在输出这些参数以后,出现错误,提示

Framework group  C1[X(CH3O)]
Deg. of freedom     9
Omega: Change in point group or standard orientation.

Old framework group not present or wrong length.
Error termination via Lnk1e in c:\g03w\l202.exe at Fri Oct 24 20:23:23 2008.


使用GView打开输出文件时,点上“read intermediat geometries”也只能看到一步的结构图

请问出现这种错误的原因是什么
应该怎样解决呢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)

好象你的命令行写得有问题,IRC计算分为向前和向后计算两种,用关键词forward或reverse,除此之外还必须指定计算的步长,计算的点数等。如:
#n MP2/6-31G* irc(rcfc,forward,stepsize=5,maxpoints=40) geom=check guess=check iop(2/16=1)
2楼2008-10-24 21:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

houwenjie

铁虫 (初入文坛)



point group出问题的话,在#T UHF/6-31G(d) IRC=RCFC Guess=Read Geom=AllCheck Test一行加上iop(2/16=1):
这是在点群改变情况下,继续irc的
3楼2008-10-24 21:04:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作

好的
多谢大家
补充一下
当不进行
--Link1--
%Chk=H3CO
%NoSave
#T UHF/6-31G(d) IRC=RCFC Guess=Read Geom=AllCheck Test

计算
只进行优化和频率计算时
可以找到一个虚频,数值为-2603.21
但是为什么加上IRC计算就出错呢?
4楼2008-10-24 22:49:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Jldbaobei 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见