怎么计算隐式溶剂环境下的反应自由能垒

ljc2877437407 · 发表于 Post on 2025-8-25 20:02:27

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

Stardust0831 · 发表于 Post on 2025-8-25 20:40:07

反应的生成物和反应物，取IRC的两端的结构，继续做opt和freq得到。

ljc2877437407 · 发表于 Post on 2025-8-25 20:45:21

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

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

Stardust0831 · 发表于 Post on 2025-8-25 21:06:21

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

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

ljc2877437407 · 发表于 Post on 2025-8-25 21:07:11

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

好滴，感谢老师

sobereva · 发表于 Post on 2025-8-26 05:12:24

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 怎么计算隐式溶剂环境下的反应自由能垒

浏览过的版块