计算化学公社

标题: OHCO体系过渡态搜索问题 [打印本页]

作者
Author: Mikasa 时间: 2016-6-17 22:22
标题: OHCO体系过渡态搜索问题

对于图中OH-CO体系的两个过渡态，用MP2/aug-cc-pVDZ进行过渡态搜索时遇到以下问题：

1. 即使选取与cis-TS1接近的结构进行过渡态搜索，也总是会找到trans-TS1；

2. 冻结HOCO二面角（文献中为零）进行过渡态搜索，得到的结果有两个虚频；

3. 用MP2/aug-cc-pVTZ也是同样的结果

4. 用CCSD(T)/aug-cc-pVTZ计算但还没得到结果

请问这应该怎么解决呢？

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

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

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

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

作者
Author: 北纬18° 时间: 2016-6-17 23:11
哪篇文献？

作者
Author: Mikasa 时间: 2016-6-18 09:04

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

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

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

作者
Author: Mikasa 时间: 2016-6-18 09:05

北纬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组的很多文章都是该体系的。

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

或者买个Molpro。

作者
Author: Mikasa 时间: 2016-6-20 08:42

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

或者 ...

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

作者
Author: sobereva 时间: 2016-6-20 16:09

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

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

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

作者
Author: Mikasa 时间: 2016-6-21 22:05

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算的能量值有问题么？

作者
Author: Mikasa 时间: 2016-6-21 23:10

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

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

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

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