24小时热门版块排行榜    

查看: 1074  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息,有中的吗 +10 archvillain 2026-08-02 12/600 2026-08-02 18:46 by hitlend
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:44 by ca3DVAx4OB19
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 4/200 2026-08-02 18:41 by ca3DVAx4OB19
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 10/500 2026-08-02 18:18 by ca3DVAx4OB19
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 18:14 by ca3DVAx4OB19
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 AS6fX1tr1on4 2026-08-02 12/600 2026-08-02 18:14 by ca3DVAx4OB19
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 18:11 by ca3DVAx4OB19
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 18:11 by ca3DVAx4OB19
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[基金申请] 准备明年的基金了 +7 Tide man 2026-07-27 7/350 2026-08-02 16:43 by gltch
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +18 低垂的野花 2026-07-31 24/1200 2026-08-02 16:03 by swallowwh
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见