24小时热门版块排行榜    

查看: 2665  |  回复: 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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见