能垒计算时是否需要加弥散函数

wwp

如老师所说的：一种粗略但比较省事的计算溶剂下焓或自由能的方法如下：
(1)用opt freq优化结构同时获得焓或自由能的热校正量（获得焓的热校正量时建议用ZPE校正因子）。溶剂模型是否需要用看情况
(2)在(1)的结构下结合溶剂模型用高级别计算单点
然后将1和2的结果相加，再加上1.89 kcal/mol即可。

体系是一个阴离子，使用的方法和基组是b3lyp-d3（bj）/6-311g**。我想请问一下在（1）的几何优化步骤中得到的自由能的热校正量是否需要添加弥散函数？使用的方法和基组是b3lyp-d3（bj）/6-311g**，体系是一个阴离子。因为在添加弥散函数之后优化出来的结果和没有添加的结果相差巨大。

liyuanhe211

对“负电荷比较多的”原子加弥散。
如果没有计算错误，以加了弥散的为准
但是你说“加/不加弥散，优化的结果差别巨大”仍然听起来不正常。怎么个巨大法

wwp

liyuanhe211 发表于 2020-3-7 19:52
对“负电荷比较多的”原子加弥散。
如果没有计算错误，以加了弥散的为准
但是你说“加/不加弥散，优化的 ...

李老师在没加弥散下 -OOH中的氢与氧相连接；加了弥散后-OOH中的氢却跑到了氮上面去，氮是负一价氮

wwp

wwp 发表于 2020-3-7 20:39
李老师在没加弥散下 -OOH中的氢与氧相连接；加了弥散后-OOH中的氢却跑到了氮上面去，氮是负一价氮

李老师我在一篇帖子看到：
• 加弥散对改进结果有益，计算条件允许可以加弥散：
计算反应势垒
优化阴离子体系或弱相互作用体系的结构、对其进行振动分析

里面的计算反应势垒是指哪些方面呐？包括那些自由能热校正的结果吗？

liyuanhe211

wwp 发表于 2020-3-7 20:39
李老师在没加弥散下 -OOH中的氢与氧相连接；加了弥散后-OOH中的氢却跑到了氮上面去，氮是负一价氮

键连方式都不一样了，显然不是常规意义上的“弥散有助于改善势垒”的范畴，那是描述同一个过程时的比较。直观上这似乎更多的是优化算法的事情。虽然质子化在氮上而不是过氧醇上确实从直观上更稳定，但不至于N为负离子、过氧醇中性的极小点就不存在。可以尝试优化的时候把步长减小，并且调整初猜的二面角不要让O-H和C=N在同一侧

liyuanhe211

wwp 发表于 2020-3-7 20:46
李老师我在一篇帖子看到：
• 加弥散对改进结果有益，计算条件允许可以加弥散：
计算反应势垒

这里一般指电子能。当然加了热力学量完了之后电子能的改变也还在里头。

wwp

liyuanhe211 发表于 2020-3-8 02:34
这里一般指电子能。当然加了热力学量完了之后电子能的改变也还在里头。

李老师 我还想请问一下：就按着上面的那个阴离子体系我能否利用在b3lyp-d3（bj）/6-311g**优化后得到的自由能热校正量+在高精度下得到的单点能，来作为体系的自由能呢？谢谢老师了。

liyuanhe211

wwp 发表于 2020-3-8 07:54
李老师 我还想请问一下：就按着上面的那个阴离子体系我能否利用在b3lyp-d3（bj）/6-311g**优化后得到的自 ...

高精度单点要给负电荷集中的原子带弥散。优化不带弥散通常没有大问题，但是如果你已经知道加、不加弥散对优化结果有显著影响，就应该搞清楚这个影响是怎么回事儿。

wwp

liyuanhe211 发表于 2020-3-8 12:19
高精度单点要给负电荷集中的原子带弥散。优化不带弥散通常没有大问题，但是如果你已经知道加、不加弥散对 ...

好的李老师，我看看哪里出现了问题。那对于自由能热校正量来说不带弥散也是可以的吗？