sobereva 发表于 2016-6-20 16:09
F12=显式相关计算，看DOI: 10.1002/wcms.68

步长上限先设小了再说calcall的事

老师，我知道我的问题了，我把不同方法算得能量混用了～

sobereva 发表于 2016-6-20 16:09
F12=显式相关计算，看DOI: 10.1002/wcms.68

步长上限先设小了再说calcall的事

老师，我用B3LYP试过了，依然没找到。用M06-2X/aug-cc-pVDZ也是有两个虚频，但是用M06-2X/aug-cc-pVTZ找到了，频率值也与相关文献一致，但是算出的能量值与文献中的差很多，整个势能面结构都不对了。
老师，这说明M06-2X/aug-cc-pVTZ算的能量值有问题么？

Mikasa 发表于 2016-6-20 08:42
嗯呢，原作者就是在Molpro下进行的～不过也没太明白作者的F12是啥意思～

F12=显式相关计算，看DOI: 10.1002/wcms.68

步长上限先设小了再说calcall的事
初猜不好势垒又潜的情况极容易一下就跑其它极小点去了

jjjspring 发表于 2016-6-20 04:01
F12下的TZ相当与正常的QZ或5Z，高斯没有CCSD(T)解析梯度完全做不了这个。先用B3LYP或M062x试试吧。

或者 ...

嗯呢，原作者就是在Molpro下进行的～不过也没太明白作者的F12是啥意思～

F12下的TZ相当与正常的QZ或5Z，高斯没有CCSD(T)解析梯度完全做不了这个。先用B3LYP或M062x试试吧。

或者买个Molpro。

北纬18° 发表于 2016-6-17 23:11
哪篇文献？

Quasi-Classical Trajectory Study of the HO + CO → H + CO2 Reaction on a New ab Initio Based Potential Energy Surface
Guo Hua组的很多文章都是该体系的。

sobereva 发表于 2016-6-17 22:39
若想重现文献数据，先确保所用级别和文献中完全一致，别乱用其它的

把步长上限调小，如maxstep=10，免得 ...

谢谢老师。
我用calcall了，但还是没找到。maxstep没调整过。
文献中是用CCSD(T)-F12/aug-cc-pVTZ计算得到体系的势能面。
MP2/aug-cc-pVDZ是参考其它文献选的，不过好像在这不是那么划算。

哪篇文献？

若想重现文献数据，先确保所用级别和文献中完全一致，别乱用其它的

把步长上限调小，如maxstep=10，免得一下走得太多。还不行就用calcall

CCSD(T)找过渡态就别指望了

PS：这种体系用当前级别极为不划算。B3LYP/TZVP找过渡态足矣
谈谈量子化学中基组的选择
http://sobereva.com/336