计算化学公社

标题: 关于反应过渡态计算 [打印本页]

作者
Author: rwshi 时间: 2015-1-28 13:53
标题: 关于反应过渡态计算
最近，我采用不同的算法和基组[一个是b3lyp/6-311+g(d,p),一个是hf/6-31+g(d)]去计算了同一个自由基反应的过渡态。
结果是采用hf/6-31+g(d)可以得到我所期望的过渡态，但是能垒是很低的；而采用b3lyp/6-311+g(d,p)确得不到我所想要的过渡态。因为实验上人家有做过类似的反应，并且能观察到反应产物。因此，这个过渡态是否存在我还在考虑当中，有可能是不存在的。请教大家有什么好的建议？

另外，在计算的时候电荷和自旋多重度我是都设置的一样的，为什么用gaussview查看b3lyp/6-311+g(d,p)计算的结果会出现“unresolvable inconsistency between charge and multiplicity and orbital occupancies. Orbital data will be ignored”这样的警告语句。请大家指教，谢谢！

作者
Author: sobereva 时间: 2015-1-28 14:23
HF本身会高估势垒，此时势垒又很低，换成HF成份较低的泛函后没有势垒是可以理解的。建议再尝试M062X/6-311+G**，此方法算过渡态和势垒比较理想。

如果确信电荷和自旋多重度都没设错，输出信息中也看不出什么问题，那就是gview判断有问题了，不用管它。

作者
Author: rwshi 时间: 2015-1-28 16:11

sobereva 发表于 2015-1-28 14:23
HF本身会高估势垒，此时势垒又很低，换成HF成份较低的泛函后没有势垒是可以理解的。建议再尝试M062X/6-311+ ...

好的我再试试谢谢 sobereva

作者
Author: brothers 时间: 2015-1-29 15:22
前些天计算过渡态，发现文献报道的一种过渡态用我的方法来算是个极小值......看来还得多试几种方法。
exlploring chemistry with electronic structure method里确实也提到了一个例子，过渡态用不同方法计算会有完全不同的结果

作者
Author: Shannon 时间: 2015-1-29 20:23
我觉得HF的可信度更高点。 GGA泛函找过渡态被坑了好多次，反应势垒低估的太多了，很多时候甚至找不到过渡态。用HF另可错找一千，不可漏过一个。
另外，试过很多半经验方法算过渡态，很多时候 PM3/AM1 能找到过渡态，换成PM6就不行了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)