经典高斯总结.doc

一般地,优化所得驻点的性质(极小点还是过渡态)要靠频率来确定;而对过渡态,要确定反应路径(即到底是哪个反应的过渡态)必需要做 IRC 了,不然靠不住的(往往用 QST 找到的过渡态并不一定就是连接输入反应物和产物的过渡态)。在我们用 QST2 或 QST3 来优化过渡态时,需输入反应物和产物,实际上反应物和产物的输入顺序是没有关系的。就是说,先输反应物后输产物和先输产物后输反应物得到的是同样的过渡态。这也好理解, QST2 里对过渡态的初始猜测实际上是程序自动将输入的反应物和产物的各变量取个平均,所以输入顺序是没有关系的。对 QST3 和 TS ,需人为指定过渡态的初始猜测。上面说的反应物和产物的输入顺序没有关系,有个前提条件,就是反应物和产物的自旋多重度一致。对 QST3 , 因为是人为指定过渡态的初始猜测,所以没有影响;而对 QST2 , 过渡态的自旋多重度默认和后面输入的一致。如果反应物和产物的自旋多重度一致,那随便先输哪个都没关系;而如果反应物和产物的自旋多重度不一致,这时该怎么办??? 当反应物基态多重度为 1,产物基态多重度为 3时,到底将过渡态的多重度定为 1还是 3?还是两个都试,哪个能量低取哪个?假如取 1,在下来做 IRC 验证时, Reverse 还好办, Forward 却是按多重度为 1做的, 这样怎么能和多重度为3的产物连接起来? 要将其关联起来,需有一个为激发态,这样才能在自旋多重度上保持一致。可这种处理对吗?还有更好的思想吗?另外,我觉得 IRC 难用得很,很难控制, 也不知自己钻进了哪条死胡同!! 一直对此很迷茫,希望各位大侠援助啊! 对于这点确实很迷茫,我觉得好像得用 cas 解决,有看过类似的文献,上面有用 cas 得出过反应物基态经激发态再回到基态得出产物的反映路径。我没有试过, 对于 cas 我很是头疼,理论和实践都一无所知。作 IRC 能说明哪些问题? irc做完,就说明我们完成了这个反应路径地计算, 在 SUMMARY OF REACTION PATH FOLLOWING: , 我们可以看到分子过渡态的键长键角与能量随着反应坐标的变化而变化,如果你将反应坐标与能量作图, 就可以得到一条过渡态曲线,由于 irc 计算是结构沿着反应路径的方向,在每个点进行优化的,所以如果你找的过渡态是正确且优化是成功的话, TS 确实是连接两个 minimum 的。一般 IRC 的两个终点并不恰好就是反应物和生成物,除非不断加大 Maxpoints 。然而此时 IRC 又往往会出错, 或者并不能达到指定的 Maxpoints ,例如,我想正反两个方向各算 20 各点, IRC 计算正常完成,但结果中并没有 41 个构型。如果这时作能量-反应坐标曲线,如何确定反应物和生成物的反应坐标? 如果你的过渡态确实做得漂亮,那么就一定可以成功地由过渡态找到与之对应的反应物、生成物。可这只是理论上的东东,现实中,由于反应的势能面实在太复杂, 你千辛万苦所得到的过渡态是很难天衣无缝地连接反应物、生成物的。不断加大 Maxpoints 简直是 Kill time 。一般地,我们可以将 IRC 最后得到的两个结构的所有数据分别与你做 QST 时所用的两个结构进行对比,只要结构数据的差别、能量差别比较小,就可以基本认定(因为反应的势能面太复杂)你 IRC 是成功的了。说得好!取点越多越好,但我发现往往作不到这一点,所以我也是将 IRC 最后得到的两个结构的数据分别生产物态和产物态的两个结构进行对比。而且我发现,作 IRC 时频率计算给出的讯息很重要,呵呵! 过渡态寻找小结刚刚做了一段时间的过渡态,期间碰到了许多的困难。寻找过渡态不是一件容易的事(对于我和大多数刚涉及量化的人来说) ,因此我希望通过写这个经验小结能对大家有些帮助。 ,用哪种方法来寻找过渡态? GAUSSIAN 提供的方法是 QSTN 和 TSN 方法。两种方法各有优点和缺点。 QSTN 方法特别 QST3 方法要求输入反应物,过渡态的猜测结构,产物这三者的结构。特别麻烦。但很管用,一般不会出现不收敛的情况。对于 TSN (对应关键词为 OPT=TS )方法,只要求输入过渡态的初始结构,但这个初始结构非常的关键,如果结构不好,则很容易出现不收敛的情况。所以我建议,如果是刚开始做过度态的话,用 QSTN 方法是好的选择,等有了“感觉”之后,再用 TS N 方法。 ? 在找过度态的时候, 经常碰到的一些问题就是不收敛(1).,有一个错误的本征值( 错误信息为: there isa wrong sign eigenvalue in hessian matrix.....) (2), 和 LINK9999 错误导致退出。(