关于反应过渡态计算

rwshi

最近，我采用不同的算法和基组[一个是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”这样的警告语句。请大家指教，谢谢！

Shannon

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

brothers

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

rwshi

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

好的  我再试试   谢谢 sobereva  

sobereva

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

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