| 查看: 1942 | 回复: 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朵)
» 猜你喜欢
085404(计算机技术)293分 希望调剂到材料相关专业
已经有2人回复
085404(计算机技术)293分 希望调剂到材料相关专业
已经有1人回复
物理化学论文润色/翻译怎么收费?
已经有84人回复
085404(计算机技术)293分 希望调剂到材料相关专业
已经有0人回复
汕头大学化学招调剂,相关内容如下。
已经有4人回复
已知面上挂了,小学校实验方向无平台,感觉太难了
已经有26人回复
求助各位大佬,球球了,在这科萨头
已经有0人回复
» 本主题相关价值贴推荐,对您同样有帮助:
求助高斯文件输入
已经有5人回复
知道赝势基组Hay-Wadt VDZ (n+1) ECP和原子坐标,怎样写高斯的几何优化输入文件
已经有13人回复
怎样在linux系统上提交高斯的输入文件,运算已经查看输出文件
已经有3人回复
请问如何在高斯的结果文件里找到HOMO和LUMO的能量?
已经有8人回复
高斯rwf文件的大小,是随计算增大,还是事先就分配好的?
已经有5人回复
求助高斯输出文件中坐标系怎么定义的?
已经有11人回复
新手求助,高斯计算出现问题~
已经有22人回复
新手求助高斯输入文件写法,求大神帮忙
已经有8人回复
关于高斯PCM模型输入问题
已经有16人回复
同一个输入文件,用windows版的高斯和linux版的高斯分别计算,怎么计算结果不一致?
已经有12人回复
高斯l9999错误 不收敛怎么办?(有图有真相,大家帮看看)
已经有31人回复
优化好的分子结构,用gaussview打开后保存为高斯输入文件,再用该文件做TD,对吗
已经有5人回复
高斯输出文件转化为mol2文件时出现问题 -c resp 与 -c bcc
已经有5人回复
急急急!高斯生成不了wfn文件怎么办
已经有7人回复
晶体数据库找到的分子晶体结构怎么转化成高斯输入文件呢
已经有8人回复
为什么一计算就会出现系统找不到指定文件呢?求高斯大神解答啊....
已经有5人回复
有个求解微分方程用到高斯赛德尔迭代的,算不太明白啊,各位帮忙看看
已经有5人回复
请教各位大侠,高斯技术部的回复。控制电子排布轨道的具体操作如何执行?
已经有4人回复
【资源】如何阅读高斯的计算输出文件
已经有4人回复
【求助】运行高斯-找不到输入文件????
已经有14人回复
【求助】高斯b3lyp/VDZ的输入文件写法
已经有14人回复
【求助】各位大哥,请问怎么看高斯的输出文件?
已经有12人回复

xieshen0321
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 4406.6
- 散金: 137
- 红花: 3
- 帖子: 784
- 在线: 341.4小时
- 虫号: 2052932
- 注册: 2012-10-10
- 性别: GG
- 专业: 理论和计算化学
2楼2014-10-22 21:38:53
meteoric30
木虫 (正式写手)
- 应助: 43 (小学生)
- 金币: 1719.5
- 散金: 12
- 红花: 2
- 帖子: 436
- 在线: 320.7小时
- 虫号: 1046846
- 注册: 2010-06-24
- 性别: GG
- 专业: 催化化学
3楼2014-10-23 03:43:31
阿Q~~
至尊木虫 (文坛精英)
- 应助: 77 (初中生)
- 金币: 33343.6
- 红花: 107
- 帖子: 36186
- 在线: 1028.1小时
- 虫号: 3237212
- 注册: 2014-05-27
- 性别: MM
- 专业: 宇宙学

4楼2014-10-23 06:49:21
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
孙1986: 金币+10, 您的解说受益匪浅 2014-10-23 11:11:18
小木虫: 金币+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
孙1986
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1046.9
- 散金: 1087
- 红花: 19
- 帖子: 455
- 在线: 243.1小时
- 虫号: 1405621
- 注册: 2011-09-17
- 专业: 原子和分子物理

6楼2014-10-23 11:13:28
fatpig8832
铁杆木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 7420.2
- 散金: 31
- 红花: 28
- 帖子: 1438
- 在线: 610.9小时
- 虫号: 770403
- 注册: 2009-05-14
7楼2014-10-23 16:50:19
yongma2008
木虫 (著名写手)
论文?苦!主。
- 应助: 34 (小学生)
- 金币: 2089.6
- 散金: 1670
- 红花: 15
- 帖子: 1610
- 在线: 398.5小时
- 虫号: 1091013
- 注册: 2010-09-06
- 性别: GG
- 专业: 理论和计算化学
8楼2014-10-23 16:52:17
syhu_007012
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 5458.6
- 红花: 3
- 帖子: 348
- 在线: 572.6小时
- 虫号: 68086
- 注册: 2005-05-16
- 专业: 药物化学
9楼2014-10-23 18:52:36













回复此楼
syhu_007012
