| 查看: 4559 | 回复: 9 | |||
sddtc888至尊木虫 (著名写手)
Hψ=Eψ S=klnΩ dS≥δQ/T
|
[求助]
求助关于过渡态输入文件问题 已有3人参与
|
|
我要找一个过渡态，输入文件如下，但是我一运行Guassian就报错： Old curvilinear step not converged, using linear step: SCX= 7.85D+00 DXMaxT= 1.77-317 SCLim= 0.00D+00 Fact= 0.00D+00 RedCar/ORedCr failed for GTrans. Error termination via Lnk1e in /opt/g09/l101.exe at Tue Oct 21 21:39:44 2014. Job cpu time: 0 days 0 hours 1 minutes 15.0 seconds. File lengths (MBytes): RWF= 5 Int= 0 D2E= 0 Chk= 1 Scr= 1 求助怎么办？是不是我初始结构没设置好？还有关于原子顺序必需一致，我已经做到一致了（我用gview分别画了初始结构和最后结构，原子序号和元素符号我都对好了的） 输入文件如下： %nprocshared=8 %mem=8GB %chk=mg4-h2-qst2.chk #p opt(calcfc,qst2,maxcycle=2000) freq=noraman b3pw91/6-311+g(d,p) empiricaldispersion=gd3 scf(maxcycle=2000,conver=6) nosymm mg4-h2-qst2 0 1 Mg -2.91504837 1.26284766 1.92151075 Mg -6.05287528 0.15350234 0.00000000 Mg -2.91508887 -2.06530526 0.00000000 Mg -2.91504837 1.26284766 -1.92151075 H -2.09425764 -4.95707567 0.32908968 H -1.73598875 -4.83041338 -0.32908968 mg4-h2-qst2 0 1 Mg 1.12634400 -1.03183000 -0.98056000 Mg 1.12757200 1.03115000 0.98044700 Mg -1.12733800 -0.98027800 1.03142900 Mg -1.12657400 0.98095100 -1.03131800 H 2.44087300 -0.00116900 -0.00039100 H -2.44092400 0.00124900 0.00041300 @Civolf 谢谢！！！！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : mg4-h2-qst2.com
2014-10-21 22:17:13, 941 bytes
» 收录本帖的淘帖专辑推荐
 计算化学 |
» 猜你喜欢
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有1人回复
-
中科院西安光机所——促进学术交流，推动产学研融合，分享领域成果，加强校企合作。
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有78人回复
-
求助：三乙酰丙酮铁（Fe(acac)₃）标准红外(IR)与XRD谱图/数据
已经有1人回复
-
Comsol 6.0 桶内肥料产热产湿过程模拟求助
已经有1人回复
-
分子动力学 Plumed教程求助
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
申请2027年材料/化学类博士
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何写两个NO分解成一个N2和O2的过渡态，和反应历程，能量变化？
已经有3人回复
- 过渡态优化问题求助
已经有13人回复
- 新手求助高斯输入文件写法，求大神帮忙
已经有8人回复
- 为什么会出现l9999错误及过渡态初始结构怎么设定？
已经有14人回复
- gaussian优化l103错误求助
已经有3人回复
- 计算过渡态报l999错误
已经有12人回复
- GaussianView内坐标格式的输入文件解释和QST2问题?
已经有7人回复
- shapefile文件属性表中添加属性出错
已经有7人回复
- 怎么寻找过渡态和中间体
已经有11人回复
- vasp计算过渡态问题
已经有13人回复
- 求助SiH2 + H2 --> SiH4过渡态的计算
已经有12人回复
- Gaussian 03计算过渡态问题
已经有4人回复
- 模拟键的断裂是不是不需要做过渡态的计算？
已经有10人回复
- 【求助】opt=QST2输入方法
已经有10人回复
- 【求助】在gaussian过渡态优化计算出现的问题~请教高手~
已经有3人回复
- 【求助】用vasp搜寻过渡态的请进
已经有16人回复
- 【求助】怎样把ChemBio 3D 文件转化为 Gaussian 输入文件
已经有6人回复
- 【求助】请教一个关于gaussian的问题
已经有17人回复
- 【求助】求助高手一个做过渡态输入文件的问题
已经有6人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
sddtc888: 金币+200, ★★★★★最佳答案, 谢谢！ 2014-10-25 13:52:22
感谢参与，应助指数 +1
sddtc888: 金币+200, ★★★★★最佳答案, 谢谢！ 2014-10-25 13:52:22
|
该问题是由于反应物和产物的结构相差过大，无法产生合适的内坐标而引起的。 解决方法是添加合适的虚拟键，使得反应物和产物具有相同的connectivity，这样有助于程序产生合适的内坐标 另外，maxcycle=2000不推荐，通常不能解决问题，还浪费计算时间；DFT缺省freq=noraman，所以只需写freq #p opt=(calcfc,qst2) freq b3pw91/6-311+g(d,p) nosymm geom=connectivity empiricaldispersion=gd3 mg4-h2-qst2 0 1 Mg -2.91504837 1.26284766 1.92151075 Mg -6.05287528 0.15350234 0.00000000 Mg -2.91508887 -2.06530526 0.00000000 Mg -2.91504837 1.26284766 -1.92151075 H -2.09425764 -4.95707567 0.32908968 H -1.73598875 -4.83041338 -0.32908968 1 2 0.5 5 0.5 3 0.5 4 0.5 2 5 0.5 4 0.5 3 4 0.5 6 0.5 4 6 0.5 5 6 0.5 6 mg4-h2-qst2 0 1 Mg 1.12634400 -1.03183000 -0.98056000 Mg 1.12757200 1.03115000 0.98044700 Mg -1.12733800 -0.98027800 1.03142900 Mg -1.12657400 0.98095100 -1.03131800 H 2.44087300 -0.00116900 -0.00039100 H -2.44092400 0.00124900 0.00041300 1 2 1.0 5 0.5 3 0.5 2 5 0.5 4 0.5 3 0.5 3 4 1.0 6 0.5 4 6 0.5 5 6 0.5 6 |
2楼2014-10-22 06:29:06
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
3楼2014-10-22 07:40:41
4楼2014-10-22 07:56:54
感谢参与，应助指数 +1
|
本帖内容被屏蔽 |
5楼2014-10-22 08:27:15
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
6楼2014-10-22 09:06:00
|
那就对了，你说的就是 “@Civolf”，楼主应该是希望Civolf这个ID给他回答吧 “@”在高斯中的意思就是在该处插入额外文件 http://www.gaussian.com/g_tech/g_ur/m_input.htm The contents of an external file may be included within a Gaussian 09 input file using the following syntax: @filename. This causes the entire file to be placed at the current location in the input stream. Appending /N to such commands will prevent the included file’s contents from being echoed at the start of the output file. |
7楼2014-10-22 10:14:18
8楼2014-10-22 10:16:35
sddtc888
至尊木虫 (著名写手)
Hψ=Eψ S=klnΩ dS≥δQ/T
- 应助: 181 (高中生)
- 贵宾: 2.508
- 金币: 10753.4
- 散金: 1948
- 红花: 29
- 沙发: 1
- 帖子: 1609
- 在线: 1165.1小时
- 虫号: 1410152
- 注册: 2011-09-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
我按照您的方法，结果算出来了，你看一下对不对！ 还有，我想问一下，有了log文件，怎么看过渡态的结构，还有，怎么从log文件中看过渡态的频率（主要想看看是不是只有一个虚频），还有，过渡态每一步的能量怎么从log文件中提取 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : mg4h2-qst2.log
2014-10-23 14:01:39, 1.58 M
9楼2014-10-23 14:01:13
10楼2014-10-25 05:36:42
15