计算化学公社

标题: 活化Gibbs自由能小于零应如何处理？ [打印本页]

作者
Author: Synthcat 时间: 2022-3-30 23:33
标题: 活化Gibbs自由能小于零应如何处理？
本帖最后由 Synthcat 于 2022-3-30 23:35 编辑

各位老师，小弟最近再计算一个质子转移反应的过渡态时，结果显示过渡态的能量比反应物-产物复合物能量要低，请问这种情况有什么解释？或者应该如何处理？计算所用的级别：M062X-D3/def2TZVPP//M062X-D3/def2svp,溶剂模型SMD(甲醇)。
(, 下载次数 Times of downloads: 13)

作者
Author: QuantumicGuy 时间: 2022-3-30 23:48
群里看到的，希望能帮到你

作者
Author: Synthcat 时间: 2022-3-30 23:57

QuantumicGuy 发表于 2022-3-30 23:48
群里看到的，希望能帮到你

非常感谢！我大概明白了，如果遇到这种情况就是说明电子能垒比较小了

作者
Author: sobereva 时间: 2022-3-31 14:37
对于这种情况，如果想分析得比较细致，建议分别考察：
电子能垒
焓的热校正量对自由能能垒的影响
-T*delta_S对自由能垒的影响
此外，还可以用Shermo进一步分解不同振动模式对自由能热校正量的不同成分的影响。

都分解开，自然就能明白导致自由能垒为负的因素。最常见的情况是：计算不当/弄错了数据，或者电子能垒本来就很小。

作者
Author: Synthcat 时间: 2022-3-31 18:44

sobereva 发表于 2022-3-31 14:37
对于这种情况，如果想分析得比较细致，建议分别考察：
电子能垒
焓的热校正量对自由能能垒的影响

感谢社长，我去尝试一下

作者
Author: Synthcat 时间: 2022-3-31 22:16

sobereva 发表于 2022-3-31 14:37
对于这种情况，如果想分析得比较细致，建议分别考察：
电子能垒
焓的热校正量对自由能能垒的影响

社长，抱歉还有一个问题，同样是这个反应，为什么从过渡态到产物的能量变化是正值呢？我检查了一下，应该没有计算错误。

作者
Author: sobereva 时间: 2022-4-1 00:38

Synthcat 发表于 2022-3-31 22:16
社长，抱歉还有一个问题，同样是这个反应，为什么从过渡态到产物的能量变化是正值呢？我检查了一下，应该 ...

可能当前过渡态不是直接连接这个产物的过渡态，做个IRC进行检验。亦可能优化用的级别和单点级别太悬殊，导致能量顺序定性不符

作者
Author: 北大-陶豫 时间: 2022-4-1 09:53
你从D1到TS为啥凭空冒出一个氢来？

作者
Author: 喵星大佬 时间: 2022-4-1 10:09

北大-陶豫发表于 2022-4-1 09:53
你从D1到TS为啥凭空冒出一个氢来？

应该是画错了

作者
Author: Synthcat 时间: 2022-4-1 16:31

北大-陶豫发表于 2022-4-1 09:53
你从D1到TS为啥凭空冒出一个氢来？

这个是chemdraw自动变了，我建模没有错

作者
Author: Synthcat 时间: 2022-4-1 16:32

sobereva 发表于 2022-4-1 00:38
可能当前过渡态不是直接连接这个产物的过渡态，做个IRC进行检验。亦可能优化用的级别和单点级别太悬殊， ...

好的，感谢社长

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