应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】热力学数据计算
41/1返回列表
查看: 2626  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gexiaojing

木虫 (著名写手)

[交流] 【求助】热力学数据计算已有1人参与

高斯计算后，热力学输出部分如下，请问如何得到 H、G、U、S 的数值啊？请详细点回答，谢谢！


Zero-point correction=                           0.099061 (Hartree/Particle)
Thermal correction to Energy=                    0.104208
Thermal correction to Enthalpy=                  0.105152
Thermal correction to Gibbs Free Energy=         0.070422
Sum of electronic and zero-point Energies=           -306.393552
Sum of electronic and thermal Energies=              -306.388405
Sum of electronic and thermal Enthalpies=            -306.387461
Sum of electronic and thermal Free Energies=         -306.422191

                     E (Thermal)             CV                S
                      KCal/Mol        Cal/Mol-Kelvin    Cal/Mol-Kelvin
Total                   65.391             18.370             73.094
Electronic               0.000              0.000              0.000
Translational            0.889              2.981             39.270
Rotational               0.889              2.981             25.997
Vibrational             63.614             12.408              7.828
Vibration  1             0.618              1.903              2.677
Vibration  2             0.640              1.832              2.085
Vibration  3             0.842              1.281              0.694
Vibration  4             0.887              1.178              0.582
回复此楼

» 收录本帖的淘帖专辑推荐

百宝箱

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
我的心是个集体宿舍，里面住着宅女·妖精·泼妇·大妈·萝莉·熟女，以及天使和恶魔，她们轮流当宿舍管理员。
1楼2009-08-13 15:44:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★ ★ ★ ★ ★ ★
lei0736(金币+6,VIP+0):谢谢 8-14 08:34
金币退给你吧，这么简单的问题居然还要问，好好回去看书去。
http://muchong.com/bbs/viewthread.php?tid=1405684
上面这个帖子中讲了初学者应该看什么书。
--------------------------
以下全部转自郭勇《高斯使用指南》第五章第一节第二小节第三部分。你看看他对这个问题讲的多么详尽。
------------------------
据说热化学性质是根据正则系综理想气体展开得到的。在大话西游论坛上看到过这样一种说法：气相下计算的能量（SCF 能量）是在真空状态下计算的体系的能量，而热化学性质是在298.15K 一个大气压下条件下计算得到。
所有频率计算都包含有体系的热化学分析。默认是在298.15K 一个大气压下条件下计算得到。对每一种元素都使用含量最高的同位素。下面是热化学的输出部分
- Thermochemistry -
Temperature 298.150 Kelvin. Pressure 1.00000 Atm. 计算的具体条件：温度，气压
Atom 1 has atomic number 35 and mass 78.91834 计算所用元素的确切质量
Atom 2 has atomic number 16 and mass 31.97207
Atom 3 has atomic number 8 and mass 15.99491
可以通过参数ReadIsotopes（见Freq 关键词，矫正因子可以加也可以不加）根据自己
的具体需要设定计算条件
Zero-point correction= 0.318636 (Hartree/Particle) 零点能
Thermal correction to Energy= 0.342790
Thermal correction to Enthalpy= 0.343734
Thermal correction to Gibbs Free Energy= 0.260712
Sum of electronic and zero-point Energies= -4014.081556
Sum of electronic and thermal Energies= -4014.057402
Sum of electronic and thermal Enthalpies= -4014.056458
Sum of electronic and thermal Free Energies= -4014.139480
零点能是对分子的电子能量的矫正，表明了在0K 温度下的分子的振动能量对电子能量的影响。当比较OK 时的热动力学性质时，零点能需要被加到总能量中。同频率计算一样，零点能也需要根据所用的算法采用不同的因子进行矫正。可用ReadIsotopes 命令实现。据说优化所得的能量是在真空下进行的，所以有必要进行零点能矫正。
在较高的温度下预测体系能量时，要给总能加上内能项。内能包括特定温度下的平动能，转动能和振动能，这些因素都受到具体条件的影响。注意计算内能时已把零点能考虑进去了。不要重复计算。在考虑内能对总能的影响时，要用Readlsotopes 命令指定矫正因子，否则就要自己手动处理，很麻烦。因为内能中只有零点能需要矫正，如果自己动手，就需要将内能各项分开，算好后再加合到一起。而用命令，程序自动就把一切做好了。
E (Thermal) CV S
KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin
Total 215.104 90.314 174.735
Electronic 0.000 0.000 0.000
Translational 0.889 2.981 44.059
Rotational 0.889 2.981 37.244
Vibrational 213.326 84.353 93.432
内能恒容热容熵
现在来看几个公式：
E0=Eelec+ZPE
表示零K 时体系的能量等于Ee lec（就是优化计算中所得的SCF 能量）和零点能之和。
E= E0+ Evib+ Erot+ Etransl= E0+ E (Thermal)
表示在某一特定温度和压力下体系的能量等于零K 时体系的能量加上内能值。
H=E+RT
表示在某一特定温度和压力下体系的焓等于在这一条件下的能量和RT 之和。
G=H-TS
表示在某一特定温度和压力下体系的自由能等于在这一条件下的焓和TS 的差。
以上所有的输出值无特殊说明外都以HARTREE 为单位。

[ Last edited by yjcmwgk on 2009-8-13 at 16:10 ]
赞一下(12人) 回复此楼
2楼2009-08-13 16:09:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gongxd325

至尊木虫 (文坛精英)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+3,VIP+0):谢谢 8-14 08:35
U=Sum of electronic and thermal Energies
H=Sum of electronic and thermal Enthalpies
G=Sum of electronic and thermal Free Energies
S=S  (73.094)
赞一下(13人) 回复此楼
一个人做点好事并不难难的是一辈子做好事!!!
3楼2009-08-14 06:38:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

417233518

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by yjcmwgk at 2009-08-13 16:09:06:
金币退给你吧，这么简单的问题居然还要问，好好回去看书去。
http://muchong.com/bbs/viewthread.php?tid=1405684
上面这个帖子中讲了初学者应该看什么书。
--------------------- ...

您好，我想学习thermal cal，是个初学者，http://muchong.com/bbs/viewthread.php?tid=1405684这个网址打不开了，你能推荐一下初学者看那些书吗？万分感激
赞一下(1人) 回复此楼
4楼2011-03-21 15:35:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gexiaojing 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭