24小时热门版块排行榜    

查看: 2569  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

孙1986

金虫 (正式写手)

[交流] 高斯大侠留步,看看这个高斯输入文件有何问题?

#p fopt=(ts z-matrix,calcfc) iop(5/13=1,1/11=1) freq optcyc=500 scf(maxcyc=300) ub3lyp/6-31g(d,p) geom=connectivity

NO+CH2 reaction

0 2
N                  0.00000000    0.00000000   -0.61797712
O                  0.00000000    0.00000000    0.54072998
C                  1.31900152   -1.04941686   -1.22943436
H                  1.31900152   -0.05250786   -1.64471336
H                  1.31900152   -2.04632586   -1.64471336

1 2 2.0
2
3 4 1.0 5 1.0
4
5
自己为高斯菜鸟,现在初学阶段,模拟别人的结果,输入这样一个文件,有错误,请哪位大侠指导一二
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

临渊羡鱼,不如退而结网
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高斯技术支持

新虫 (初入文坛)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
孙1986: 金币+10, 您的解说受益匪浅 2014-10-23 11:11:18
该例子有一定的典型性,有很多高斯初学者会遇到的问题,在此详细讨论一下:
首先,这个输入文件可以运行,但是容易导致错误意想不到的结果
(1)当给出opt=z-matrix时,高斯会去读取输入文件中的Z-matrix以及变量。可是楼主给出的是直角坐标(cartesian),所以程序会按照opt=cartesian来进行优化。参考http://www.gaussian.com/g_tech/g_ur/k_opt.htm
(When a Z-matrix without any variables is used for the molecule specification, and Opt=Z-matrix is specified, then the optimization will actually be performed in Cartesian coordinates.)
使用该方法(opt=cartesian)收敛速度远远低于默认的冗余内坐标(redundant internal coordinate)
(2) fopt是配合opt=z-matrix使用的,要求对所有变量进行优化。这个选项现在已经很少人用了,通常直接使用opt即可。
(3)iop(5/13=1,1/11=1)。iop(5/13)=1要求哪怕SCF不收敛都继续计算,该选项是用来调试用的,在实际应用中不应该出现。iop(1/11=1)即noeigentest,建议在结构已经很接近过渡态时使用,否则容易得到其他意想不到的虚频,比如甲基旋转。一般的建议是,能使用关键词尽量避免使用IOP。首先IOP在多步计算中(比如opt freq)只能适用于第一步,而使用关键词对所有步骤都适用。此外,在不同版本的高斯程序中,IOP参数可能会发生改变。
(4)optcyc=N, 即opt(maxcyc=N),还有scf(maxcyc=300)。如果结构或者SCF不收敛,增大默认的周期数通常不会有改善,反而耗费更多的计算时间。
(5) 如果optcyc写成opt(maxcyc=N),同时又使用opt=(ts,z-matrix),同一个opt出现两次,这种情况应该避免,因为程序在生成非标准输入(即输出文件开始部分生成的一大堆IOP选项)时会对相同选项进行加和,可能导致意想不到的问题。
(6)总的来说,初学高斯,建议从最少的选项开始,完成一个基本计算,然后添加适当的选项以实现额外的功能。直接从其他人那里拷贝过来输入文件,可能适用于特定体系,但不一定在任何情况都适用。尤其是比如IOP(5/13=1),在正式计算中绝对应该避免出现
5楼2014-10-23 08:57:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

xieshen0321

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你不把错误信息贴出来估计神仙也难解啊!!!
2楼2014-10-22 21:38:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meteoric30

木虫 (正式写手)

#p fopt?  多个f吧。。。
3楼2014-10-23 03:43:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿Q~~

至尊木虫 (文坛精英)


小木虫: 金币+0.5, 给个红包,谢谢回帖
帮顶一下,不懂 ~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
4楼2014-10-23 06:49:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 我不理解! +8 Edward_pc 2026-08-26 13/650 2026-08-26 15:37 by Tide man
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见