回答的不是很完全
给出来的应该默认都是298.15K的热力学量
thermal correction to U,H,G等等都是和前面zero-point零点能相关的校正
一般通过Gibbs自由能可以判断体系的稳定性,以及在反应能垒等等
但U,H,G的计算都包括平动,转动,电子和振动项
前面三项的配分函数很容易处理,但振动配分函数与分子所有的简振模式有关
Gaussian输出里面还有一大堆的所谓vibrational temperature(用#p控制输出)
是和每个振动模相关的
再计算能量的时候开始采用势能面的井底当作参考0点来计算U,H,G等,与实际的振动0点
相差theta/2 (theta是所谓的振动温度),把这一项加到各个热力学函数里面从而进行零
点能校正。
U,H,G的校正值不同,是因为它们的热力学表达式不同
【 在 zhjx (\(^ω^)/) 的大作中提到: 】
: 先优化然后做频率计算
: 在输出文件里会发现这样一段
: Zero-point correction= 0.275565 (Hartree/Particle)
: Thermal correction to Energy= 0.289780
: Thermal correction to Enthalpy= 0.290725
: Thermal correction to Gibbs Free Energy= 0.233918
: Sum of electronic and zero-point Energies= -770.312188
: Sum of electronic and thermal Energies= -770.297973
: Sum of electronic and thermal Enthalpies= -770.297028
: Sum of electronic and thermal Free Energies= -770.353835
: 如果是比较0K的稳定性,就比较zero-point Energies+Zero-point correction
: ...................
--
FROM 166.111.228.*