计算化学公社

标题: 怎么计算隐式溶剂环境下的反应自由能垒 [打印本页]

作者
Author: ljc2877437407 时间: 2025-8-25 20:02
标题: 怎么计算隐式溶剂环境下的反应自由能垒
各位大佬们好。想请教一下怎么计算在溶剂环境反应的能垒呢，比如A和B产生C，A+B----C。是直接在隐式溶剂模型下计算A、B、C各部分的能量，然后▲E=E(ts)-E(B)-(A)，还是需要分别计算各个部分的溶质在溶剂环境下标况（298.15K、1M）时的自由能，再通过这个公式▲E=E(ts)-E(B)-(A)计算得到嘛

作者
Author: Stardust0831 时间: 2025-8-25 20:40
反应的生成物和反应物，取IRC的两端的结构，继续做opt和freq得到。

作者
Author: ljc2877437407 时间: 2025-8-25 20:45

Stardust0831 发表于 2025-8-25 20:40
反应的生成物和反应物，取IRC的两端的结构，继续做opt和freq得到。

直接在隐式溶剂模型下计算就行了嘛，不用计算各个部分的溶质自由能，公式为（E=真空的能量+溶解自由能+1.89）嘛

作者
Author: Stardust0831 时间: 2025-8-25 21:06

ljc2877437407 发表于 2025-8-25 20:45
直接在隐式溶剂模型下计算就行了嘛，不用计算各个部分的溶质自由能，公式为（E=真空的能量+溶解自由能+1. ...

意思是您不需要把A和B分开处理。直接对A和B的复合物算自由能就可以（结构由对IRC的端点继续优化得到）

作者
Author: ljc2877437407 时间: 2025-8-25 21:07

Stardust0831 发表于 2025-8-25 21:06
意思是您不需要把A和B分开处理。直接对A和B的复合物算自由能就可以（结构由对IRC的端点继续优化得到）

好滴，感谢老师

作者
Author: sobereva 时间: 2025-8-26 05:12
计算G(ts)-G(B)-G(A)得到自由能垒（除非反应复合物很稳定，此时当单分子反应，以AB复合物作为反应物）
每个物质在溶剂下的自由能G都按下文计算就完了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

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