计算化学公社

标题: 如何计算某个反应过程的自由能变？ [打印本页]

作者
Author: 15whyhq 时间: 2019-4-29 16:24
标题: 如何计算某个反应过程的自由能变？
本帖最后由 15whyhq 于 2019-4-29 20:17 编辑

反应物的分子以及产物的分子结构分别在气相下进行opt+freq，计算得到反应的自由能变(产物的减去反应物的） ΔG1 (-52.1 kJ/mol) ;

在气相的结构下进行高精度单点能计算并考虑隐式溶剂模型，计算反应的自由能变(产物的HF值加上Thermal correction to Gibbs Free Energy(气相下的)减去反应物(反应物的HF加上Thermal correction to Gibbs Free Energy(气相下的)) ΔG2 (20.3 kJ/mol)，

ΔG1 与 ΔG2 结果相差很大，是不是算反应中各物质的自由能不能用高精度的单点能HF加上Thermal correction to Gibbs Free Energy进而求反应的自由能变呢？

作者
Author: granvia 时间: 2019-4-29 18:01
本帖最后由 granvia 于 2019-4-29 22:00 编辑

原则上可以的，但主要忧虑是振动配分函数，需假定低精度计算得到的振动频率跟高精度结果差别较小。不过，你算的Delta G2大了好几个数量级肯定是有问题的

作者
Author: DoubeeTwT 时间: 2019-4-29 19:04
http://bbs.keinsci.com/thread-9348-1-1.html
提问没有加求助哦~

理论上应该是可以的，你这G1和G2不但数量级不一样连正负号都变了，反应方向都反过来了。。。可以检查一下第二次的计算哪里有问题

作者
Author: sobereva 时间: 2019-4-30 04:43
100%是你自己计算过程错误“反应中各物质的自由能用高精度的单点能加上Thermal correction to Gibbs Free Energy”可以说几乎是必须这样做，何来不能？

作者
Author: 15whyhq 时间: 2019-4-30 21:29

sobereva 发表于 2019-4-30 04:43
100%是你自己计算过程错误“反应中各物质的自由能用高精度的单点能加上Thermal correction to ...

数据什么的都检查了，没有错啊。我算单点的时候，直接在GView里面打开气相优化好的结构，然后保存坐标算的单点，这个会有问题吗？

作者
Author: pcxlm1 时间: 2019-5-5 09:22

15whyhq 发表于 2019-4-30 21:29
数据什么的都检查了，没有错啊。我算单点的时候，直接在GView里面打开气相优化好的结构，然后保存坐标算 ...

高精度算单点是后HF计算还是DFT计算？

作者
Author: 15whyhq 时间: 2019-5-5 09:23

pcxlm1 发表于 2019-5-5 09:22
高精度算单点是后HF计算还是DFT计算？

DFT计算

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