计算化学公社

标题: 如下情况是否需要加入色散校正？ [打印本页]

作者
Author: 小强 时间: 2017-7-20 11:47
标题: 如下情况是否需要加入色散校正？
本帖最后由小强于 2017-7-20 11:47 编辑

有一体系我在b3lyp/6-31g(d)计算水平下顺利得到其过渡态(TS1, 经IRC验证)及中间体(IM1)；由于体系存在芳环，故想初略考察下色散作用对相应过渡态及中间体构型及能量的影响，因此我以上述所得的TS1和IM1为初猜结构，在b3lyp-D3/6-31g(d)计算水平下重新优化并得到过渡态(TS2)和中间体(IM2)。结果发现: TS1与TS2的构型差异很小，能量相差约0.1Bohr; 而IM1与IM2是两种完全不同的构型，能量相差也约为0.1Bohr。(ps: 我计算目的主要是考察势垒大小。) 请问各位老师和同学:
看了sob老师写的关于计算色散作用的博文，有个疑惑就是: 关键词empirical dispersion=gd3/gd3(bj) 与弥散函数是不是总是同时出现？我采用上述未加弥散函数的基组所得结果是否有参考价值？
考虑到加入弥散函数后会使计算量猛增且易发生震荡，对于上述过渡态及中间体的优化考虑色散作用是否是必须的？谢谢大家了！

作者
Author: shx 时间: 2017-7-20 17:55
结果没有问题的，加上dispersion得到的结构一般来讲更为合理一些。

作者
Author: 小强 时间: 2017-7-20 18:49

shx 发表于 2017-7-20 17:55
结果没有问题的，加上dispersion得到的结构一般来讲更为合理一些。

首先谢谢你的回复！另外再请教下你，如果加上dispersion的话，那么弥散函数在优化时需不需要加呢？我知道计算色散作用弥散函数是很重要的，但是不知道弥散函数应该是在计算单点时加上还是凡是关于色散作用的计算都应该加上？

作者
Author: 1130240115 时间: 2017-7-20 19:09
一般中性体系优化时可以不加弥散函数，只在计算单点能时加弥散函数即可。见：谈谈弥散函数和月份基组http://sobereva.com/119

作者
Author: 冰释之川 时间: 2017-7-20 19:36
首先楼主得先确定研究体系是否存在弱相互作用，并且是否要研究弱相互作用。如果不存在，那加不加色散对能垒影响不大（当然加色散肯定会有一点好处的）。研究静电主导的弱相互作用体系，在优化时可以不用带弥散。但是如果要计算弱相互作用，基组建议6-31+G(d)起步，也就是要带上弥散函数。6-31+G(d)是计算弱相互作用的底限，有些人拿6-31G(d)计算弱相互作用，纯属搞笑。如果降格到6-31G(d)，前提条件是必须用counterpoise或者gCP方式考虑BSSE问题。 具体参考《乱谈DFT-D》http://sobereva.com/83

作者
Author: 小强 时间: 2017-7-20 20:31

1130240115 发表于 2017-7-20 19:09
一般中性体系优化时可以不加弥散函数，只在计算单点能时加弥散函数即可。见：谈谈弥散函数和月份基组http:/ ...

谢谢指教，漏看了这篇博文

作者
Author: 小强 时间: 2017-7-20 21:08

冰释之川发表于 2017-7-20 19:36
首先楼主得先确定研究体系是否存在弱相互作用，并且是否要研究弱相互作用。如果不存在，那加不加色散对能垒 ...

大概明白了，我再消化消化，谢谢您的指点！

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