如何计算某个反应过程的自由能变？

15whyhq · 发表于 Post on 2019-4-29 16:24:08

本帖最后由 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进而求反应的自由能变呢？

granvia · 发表于 Post on 2019-4-29 18:01:47

本帖最后由 granvia 于 2019-4-29 22:00 编辑

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

DoubeeTwT · 发表于 Post on 2019-4-29 19:04:01

http://bbs.keinsci.com/thread-9348-1-1.html
提问没有加求助哦~

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

sobereva · 发表于 Post on 2019-4-30 04:43:48

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

15whyhq · 发表于 Post on 2019-4-30 21:29:31

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

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

pcxlm1 · 发表于 Post on 2019-5-5 09:22:12

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

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

15whyhq · 发表于 Post on 2019-5-5 09:23:59

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

DFT计算

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

[Gaussian/gview] 如何计算某个反应过程的自由能变？