24小时热门版块排行榜    

CyRhmU.jpeg
查看: 957  |  回复: 6
当前主题已经存档。

loovfnd

至尊木虫 (著名写手)

[交流] 【讨论】关于高斯的能量

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
那总能量是Sum of electronic and thermal Energies= -113.834277这一行吗?

但是吴国是的讲义里又说,
SCF Done: E(RHF) = -74.9606932329 A.U. after 4 cycles (分子总能量)
(1)到底哪个是分子的总能量?
(2)那电子的总能量怎么算呢?
(3)在高斯计算结尾处的能量是什么能量呢?这个能量和自由能什么关系呢?(如图)


(数据是从好几个地方找来的,定性的)

[ Last edited by loovfnd on 2009-12-17 at 17:11 ]
回复此楼
莫强求!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggdh

专家顾问 (正式写手)

★ ★
loovfnd(金币+2,VIP+0):感谢讨论 12-18 08:41
1,E=E0 + E(vib) + E(rol) + E(transl)这个是分子总能量
2E(RHF) = -74.9606932329 A.U. after 4 cycles 是电子总能量
3结尾处就是电子总能量
2楼2009-12-17 17:17:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyan5723

银虫 (正式写手)

★ ★
loovfnd(金币+2,VIP+0):感谢参与讨论! 12-18 08:42
loovfnd(金币+1):感谢 2010-03-12 10:02
你可以看看Thermochemistry in Gaussian这篇文章,Google里可以搜出pdf
最勇敢的事情是认清了生活的真相之后依旧热爱生活。不要害怕欺骗,但要知道世界上存在欺骗。
3楼2009-12-17 21:26:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangyan5723

银虫 (正式写手)

★ ★
loovfnd(金币+2,VIP+0):感谢提供链接 12-18 08:42
loovfnd(金币+2):感谢 2010-03-12 10:01
还可以参考http://muchong.com/html/200907/1438878.html这个帖子
最勇敢的事情是认清了生活的真相之后依旧热爱生活。不要害怕欺骗,但要知道世界上存在欺骗。
4楼2009-12-17 21:59:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)


yjcmwgk(金币+1,VIP+0):感谢转载 12-18 13:35
引用回帖:
Originally posted by ggdh at 2009-12-17 17:17:
1,E=E0 + E(vib) + E(rol) + E(transl)这个是分子总能量
2E(RHF) = -74.9606932329 A.U. after 4 cycles 是电子总能量
3结尾处就是电子总能量

分子总能量如下指出好像是:E0=E(elec) + ZPE

Gaussian计算中分子总能量各项的意义
在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,这正是统计热力学中对内能的定义.
莫强求!
5楼2009-12-18 08:50:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggdh

专家顾问 (正式写手)

★ ★
yjcmwgk(金币+2,VIP+0): 12-18 13:35
loovfnd(金币+1):感谢 2010-03-12 10:02
楼上的意思就是
E=E0 + E(vib) + E(rol) + E(transl)是室温下分子总能量
E0是0K下分子总能量
6楼2009-12-18 12:24:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)

引用回帖:
Originally posted by ggdh at 2009-12-18 12:24:
楼上的意思就是
E=E0 + E(vib) + E(rol) + E(transl)是室温下分子总能量
E0是0K下分子总能量

可以这样理解
莫强求!
7楼2009-12-18 15:26:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 loovfnd 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见