计算化学公社

标题: 关于反应能量过程绘制中吉布斯自由能计算的疑问 [打印本页]

作者
Author: liustone 时间: 2023-5-28 10:35
标题: 关于反应能量过程绘制中吉布斯自由能计算的疑问
各位老师好，本人目前在做一个反应体系的反应过程能量绘制和反应速率计算，在获得物质的自由能时有一些疑问，希望得到大家的指导，
为了减小计算代价，在进行初步优化和振动分析时，我是用了b3lyp/6-311+g**的理论水平计算，读取出一系列能量如下图：
(, 下载次数 Times of downloads: 11)
目的是获取吉布斯自用能，这里面对我有用的就是零点能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 计算得到所需的更准确的自由能是否正确？

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

作者
Author: wzkchem5 时间: 2023-5-28 16:39
不完全正确，因为不仅ZPE受校正因子影响，U(T)-U(0)和熵也受校正因子影响

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

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

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