计算化学公社

标题: 过渡态搜索中的频率问题 [打印本页]

作者
Author: jing-quan63 时间: 2019-5-24 16:07
标题: 过渡态搜索中的频率问题
大家好，最近仿照文献计算一个体系（甲基提取正戊铜端部氢原子）的过渡态，采用M062X/6-311G(2df,2p)方法进行过渡态搜素，出现两个多余虚频，分别对应甲基和正戊铜的整体震动，尝试过以下方式：
   1、在b3lyp/6-31g(d,p)优化成功，改用文献方法继续优化
   2、在b3lyp/6-311g(2df,2p)优化成功,改用文献方法继续优化
   3 、m062x/6-311g(d,p)优化成功，再采用文献方法
以上均未能消除虚频（int=ultrafine格点精度）
   最后，采用M062X/6-311(2df,2p)  结合int=superfine ,才成功消除多余虚频，问题如下：
1：对于m062x泛函，积分格点精度是否必须在整个几何优化和走IRC过程中一致，要是我用int=superfine消除虚频，难道之前int=ultrafine的优化结果都要作废吗？
2：如何根据体系确定积分格点，比如我觉得不肯上来就用int=superfine，毕竟耗时。
3：特别地，针对我上面解决的过程，有无建议意见，感谢提出
谢谢大家！

作者
Author: liyuanhe211 时间: 2019-5-24 21:49
可以用你的路径中ultrafine和superfine都能收敛到正确的过渡态的结构测试一下使用superfine的几何结构计算ultrafine下的单点、然后将superfine的频率计算得出的热力学校正量加上去的做法会造成多大的误差。

作者
Author: sobereva 时间: 2019-5-25 02:52
6-311G(2df,2p)是很垃圾的基组，大小相近的def2-TZVP比其好得多得多，仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

对于优化过渡态目的，根本没必要用6-311G(2df,2p)这个档次的，用def-TZVP或6-311(2d,p)就足够了。

M06-2X一般都用ultrafine，有诸如你这种情况搞不定的时候再试superfine
总的来说对于稍大的体系而且还涉及一批体系的研究不太建议在opt freq IRC这些阶段用M06-2X，否则这类问题比较烦人。用B3LYP-D3(BJ)一般就够了，速度还快很多，对积分格点只需要fine程度就完全够用了。

作者
Author: jing-quan63 时间: 2019-5-25 09:18

sobereva 发表于 2019-5-25 02:52
6-311G(2df,2p)是很垃圾的基组，大小相近的def2-TZVP比其好得多得多，仔细看
谈谈量子化学中基组的选择
h ...

谢谢老师，主要是我后期想采用kisthelp来计算反应速率常数，感觉优化的方法要是精度不够的话，可能热力学校正也不准，不知道我理解有没有错，或者老师有无好的建议？

作者
Author: jing-quan63 时间: 2019-5-25 09:25

liyuanhe211 发表于 2019-5-24 21:49
可以用你的路径中ultrafine和superfine都能收敛到正确的过渡态的结构测试一下使用superfine的几何结构计算u ...

好的我尝试一下，谢谢李老师

作者
Author: sobereva 时间: 2019-5-25 09:41

jing-quan63 发表于 2019-5-25 09:18
谢谢老师，主要是我后期想采用kisthelp来计算反应速率常数，感觉优化的方法要是精度不够的话，可能热力学 ...

仔细看
浅谈为什么优化和振动分析不需要用大基组
http://sobereva.com/387（http://bbs.keinsci.com/thread-6600-1-1.html）

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