关于反应能量过程绘制中吉布斯自由能计算的疑问

liustone

各位老师好，本人目前在做一个反应体系的反应过程能量绘制和反应速率计算，在获得物质的自由能时有一些疑问，希望得到大家的指导，
为了减小计算代价，在进行初步优化和振动分析时，我是用了b3lyp/6-311+g**的理论水平计算，读取出一系列能量如下图：

目的是获取吉布斯自用能，这里面对我有用的就是零点能ZPE和最后一行自由能G0，此外还有单点能HF0，另外在http://bbs.keinsci.com/thread-3805-1-1.html获取本理论水平下ZPE校正因子k为0.9881。为了获取更准确的能量，在M062X/def2TZVP下计算得到更高精度的单点能HF1，请问我参考sob老师在http://bbs.keinsci.com/thread-2121-1-1.html中解答的，使用这个公式 G=G0-HF0+HF1-ZPE+k*ZPE 计算得到所需的更准确的自由能是否正确？

请各位老师指导及指正，谢谢

wzkchem5

不完全正确，因为不仅ZPE受校正因子影响，U(T)-U(0)和熵也受校正因子影响

sobereva

仔细看下文，该说的都说了
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

opt freq目的通常毫无必要带弥散，3-zeta基组足矣