| 查看: 5378 | 回复: 32 | |||
| 【奖励】 本帖被评价21次，作者604gq增加金币 19.5 个 | |||
| 当前主题已经存档。 | |||
[资源]
Gaussian寻找过渡态经验小结
|
|||
|
寻找过渡态不是一件容易的事（对于我和大多数刚涉及量化的人来说），因此我希望通过写这个经验小结能对大家有些帮助。 1.首先遇到的问题是，用哪种方法来寻找过渡态？ GAUSSIAN提供的方法是QSTN和TSN方法。两种方法各有优点和缺点。QSTN方法特别QST3方法要求输入反应物，过渡态的猜测结构，产物这三者的结构。特别麻烦。但很管用，一般不会出现不收敛的情况。对于TSN（对应关键词为OPT=TS）方法，只要求输入过渡态的初始结构，但这个初始结构非常的关键，如果结构不好，则很容易出现不收敛的情况。所以我建议，如果是刚开始做过度态的话，用QSTN方法是好的选择，等有了“感觉”之后，再用TSN方法。 2. 怎么解决经常出现的错误？ 在找过度态的时候，经常碰到的一些问题就是（1）不收敛，（2）有一个错误的本征值（错误信息为：there is a wrong sign eigenvalue in hessian matrix.....)，（3）和LINK9999错误导致退出。对于不收敛的情况，可以分为两类，比如提示信息里的CONVERGENCE FAILER 提醒收敛到了10（-5），而此时你设定的SCF循环次数也仅仅是64步，那么完全有希望通过加大SCF循环次数来达到收敛的目的。倘若只收敛到10（-3）或10（-2），此时加大循环次数可能就没用了。结果还是CONVERGE FAILER。此时可采用SCF=QC，来达到强制收敛的目的 。因为SCF=QC（LINK508）的计算量比默认的L502要大，所以不到万不得以就不用它了。出现第二个错误可以直接用关键词OPT=NOEIGEN 来实现。LINK9999出错是因为已经走完了默认的步数，但还未完成。系统会自动跳出。出现这种情况大多数就是因为优化步数和SCF步数超过了默认值。可用OPT（MAXCYCLE=100）和SCF（MAXCYCLE=300）来改错。 3.怎么样控制过渡态的优化，使得过渡态不至于收敛到其他的分子结构中去？ 我用GAUSS VIEW 可以解决这个问题 ，当刚开始运行GAUSSIAN时，你用GVIEW去打开输出文件时，你可以看到你的过渡态的初始输入结构，当一个循环过后（从上一个LINK502到下一个LINK502），你再打开输出文件，你就可以清晰地看到优化一步后分子的构型，这样就可以随时监控过度态分子的结构，倘若已经有收敛到其他分子构型的趋势时，你就可以把它给KILL了，而不至于需要等全部工作结束后，打开输出文件才知道已经不是想要的过渡态了。 如果收敛到其它的构型上去，可以考虑缩小OPT的步长.iop(1/8=2或3）即可。 4. 还需要加其它的关键词吗？ 建议在OPT中加入CALCFC。这样可以加大找到过渡态的几率。本人深有体会！ 5. 如何控制寻找过渡态的步长？ 这个可以用IOP来实现。具体相应的IOP是iop(1/8)=2，此语句说明在找过度态的时候，以2A（A为基本单位长度）为单位来寻找过渡态。 |
回复此楼» 收录本帖的淘帖专辑推荐
 过渡态 |
» 猜你喜欢
南京大学软物质计算模拟方向硕博连读/直博生招生
已经有0人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有169人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
美国杰克逊州立大学周剑歌教授的化学计算与AI实验室诚招博士后研究员
已经有3人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Tc数据拟合
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2008-05-23 11:59:41
3楼2008-05-24 09:49:07
4楼2008-05-26 09:22:35
5楼2008-05-29 11:50:34
6楼2008-06-02 11:23:45
7楼2008-06-03 10:17:51
8楼2008-06-03 15:18:51
9楼2008-06-03 15:40:51
10楼2009-01-01 20:59:24
11楼2009-01-03 12:18:21
12楼2009-01-04 10:50:28
13楼2009-01-18 21:09:52
14楼2009-01-20 10:15:20
15楼2009-01-24 20:39:31
16楼2009-02-11 21:42:49
17楼2009-02-23 15:40:54
18楼2009-02-23 15:42:57
19楼2009-02-23 20:33:04
21楼2009-03-01 11:38:22
23楼2009-05-18 17:51:13
24楼2009-05-18 21:35:16
简单回复
caipingl20楼
2009-02-24 03:59
回复
谢谢分享
zhaooming22楼
2009-05-18 17:46
回复
100