| 查看: 2573 | 回复: 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朵）
syhu_007012» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有184人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有2人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有5人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助高斯文件输入
已经有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 (小学生)
- 金币: 2091.6
- 散金: 1670
- 红花: 15
- 帖子: 1610
- 在线: 398.5小时
- 虫号: 1091013
- 注册: 2010-09-06
- 性别: GG
- 专业: 理论和计算化学
8楼2014-10-23 16:52:17
syhu_007012
木虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 6178.1
- 红花: 3
- 帖子: 353
- 在线: 580.4小时
- 虫号: 68086
- 注册: 2005-05-16
- 专业: 药物化学
9楼2014-10-23 18:52:36
