24小时热门版块排行榜    

查看: 2570  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 青基如何破局 +3 lyz123lyz 2026-08-27 4/200 2026-08-27 13:08 by kulium
[基金申请] 基金不中,共勉 +8 eulota 2026-08-26 8/400 2026-08-27 12:34 by lice_1987
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+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
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见