在HFSCF方法计算中,体系采取固定核近似后,体系的哈密顿算符分为电子哈密顿和核哈密顿两部分,电子哈密顿分为三部分,电子动能,电子间斥能,及核电排斥,通过单电子近似,得到单电子的HF方程,自洽场的结果得到了电子的能量.核哈密顿分为两部分,核动能及排斥,由于核固定动能为零,只剩下一项,因此体系的能量就是电子能量加上核斥能,,这就是我们计算得到的分子总能量,因为分子在OK时仍在平衡位置振动,要加上零点校正.最后结果即例子中的:
Sum of electronic and zero-point Energies= -113.837130
在分子的热力学计算中,包括了零点能的输出,零点能是对分子的电子能量的矫正,计算了在0K温度下的分子振动能量.但我们计算的热力学数据并不总是在OK时,为了计算在较高温度下的能量,内能也要考虑到总能量中,它包括平动能,转动能,和振动能.注意在计算内能的时候,已经考虑了零点能.下面是一个计算实例
Temperature 298.15 Kelvin. Pressure 1.0000 Atm.
Zero-point correction= 0.029201
Thermal corection to Energy= 0.032054
Thermal correction to Enthalpy= 0.032999
Thermal correction to Gibbs Free Energy= 0.008244
Sum of electronic and zero-point Energies= -113.837130
Sum of electronic and thermal Energies= -113.834277
Sum of electronic and thermal Enthalpies= -113.833333
Sum of electronic and thermal Free Energies= -113.858087
后面四行的四个能量分别为E0, E, H, G. 计算为
E0=E(elec) + ZPE
E=E0 + E(vib) + E(rol) + E(transl)
H=E + RT
G=H - TS
上例是在298.15K和1.0Atm下计算的热力学状态函数,EO即我们的分子体系总能量,而在求298.15K下的内能E时,要加上平动能,转动能,和振动能,由统计热力学知,这三项处于激发能级,注意这里的ZPE实际上就是对应E(vib)在OK时的能量,但现在并非0K,有关能量随温度变化公式可参考统计热力著作,最后得到的E就是平常所指的内能U,这正是统计热力学中对内能的定义.对H,G状态函数,不再叙述.
由于本版无法上传公式,还请见谅,有关内容请参看相关著作