24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2286  |  回复: 14
【奖励】 本帖被评价9次,作者yanhualu2000增加金币 5.25
当前主题已经存档。

yanhualu2000

铁杆木虫 (著名写手)


[资源] 【原创】分享我的一点经验--生成热的计算

最近在学习用Gaussian进行能量计算。经过两个星期的折磨,有一点心得(老虫就不要笑话新人了)。看到很多新人也对能量计算比较迷惑,我就在这里抛砖引玉了。对生成热的计算总结如下:
1. 输出文件中关于能量的信息:
对H2进行计算,# RHF/6-31G Freq Test,在out文件中找到下面的一行:
-------------------
- Thermochemistry -
-------------------
这一行下面就是关于能量的信息。需要注意的是,能量单位为Hartree,而不是kcal/mole或kJ/mol。

下面这一行是焓(我自己这样理解,不知道是否正确)的数值,单位为Hartree。
Sum of electronic and thermal Enthalpies=              -1.111484
前几天比较郁闷的是,我不知道能量单位是Hartree,所以走了不少弯路。

2.下面举例H2O生成热的计算过程。
首先写出反应方程式:
H2 + 1/2O2 → H2O
然后Gaussian计算H2、2O2 、 H2O的焓,计算结果如下:
H2
Sum of electronic and thermal Enthalpies=              -1.111484

O2
Sum of electronic and thermal Enthalpies=            -149.517925

H2O
Sum of electronic and thermal Enthalpies=             -75.956427

H2O的生成热计算值为:
Hf(298.15k) = -75.956427 - ((-1.111484)+ 1/2(-149.517925)) = -0.0859805 Hartree =  -0.0859805 *2625.5kJ/mol = 225.741 kJ/mol

H2O在气态时的生成焓测量值 = -242 kJ/mol。

其他化合物的生成热可以类似的方式计算。

欢迎大家讨论。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★ 三星级,支持鼓励

支持鼓励,呵呵呵,学习
2楼2009-05-18 21:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)


★★★★★ 五星级,优秀推荐


lei0736(金币+1,VIP+0):谢谢 8-11 09:47
是不是还应该考虑零点能(ZPE)校正吧?
3楼2009-05-18 21:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yhl213

金虫 (小有名气)


★★★ 三星级,支持鼓励

楼主对别的化合物计算过吗,这种计算方法是否对所有的都适合啊?
4楼2009-05-19 08:17:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanhualu2000

铁杆木虫 (著名写手)


★ ★
lei0736(金币+2,VIP+0):谢谢 6-5 19:59
引用回帖:
Originally posted by yhl213 at 2009-5-19 08:17:
楼主对别的化合物计算过吗,这种计算方法是否对所有的都适合啊?

我计算过甲醇、乙醛等,HF 6-31G误差在20%以内。
5楼2009-05-19 12:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzfish

金虫 (正式写手)


学习了,多谢楼主
6楼2009-06-05 17:55:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paoge

铁虫 (初入文坛)


楼主你qq号码多少啊
7楼2009-06-30 22:37:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiechong

金虫 (著名写手)


能不能直接计算H2,O2,H2O的单点能  然后再求差值?这样算出来的是不是也是反应热? 像零点能,转动振动能的校正一定要做吗?影响大吗?

[ Last edited by tiechong on 2009-7-1 at 07:45 ]
8楼2009-07-01 07:21:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)


★★★ 三星级,支持鼓励

★ ★
lei0736(金币+2,VIP+0):谢谢 7-17 09:35
首先必须优化!

其实freq里面才有热力学性质
9楼2009-07-01 07:36:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flyskysky

银虫 (正式写手)


★★★ 三星级,支持鼓励

非常好,谢谢!!!
10楼2009-07-16 22:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★ 三星级,支持鼓励

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 7-17 09:36
支持一下新人,鼓励鼓励!做的很好!
其实不建议使用HF。给楼主的建议:DFT计算量不大(略大于HF),但是结论远优于HF,建议使用DFT
tiechong问的也很好。给tiechong的答案:零点能,转动振动能的校正必须做
11楼2009-07-16 22:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsq_phi

木虫 (著名写手)


★★★ 三星级,支持鼓励


lei0736(金币+1,VIP+0):谢谢 8-11 09:48
8楼说的对,必须先优化!按楼主的方法输入应该是:# opt freq RHF/6-31G geom=connectivity
优化后的结果可用。
12楼2009-07-22 21:30:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

penghcp

木虫 (小有名气)


★ 一星级,一般


lei0736(金币+1,VIP+0):谢谢 8-11 09:48
首先必须优化
你的基组 方法都比较低
13楼2009-07-23 20:30:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuyihua

铁杆木虫 (小有名气)


★★★★★ 五星级,优秀推荐

学习了,多谢楼主
14楼2009-08-11 08:34:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)


引用回帖:
Originally posted by yjcmwgk at 2009-7-16 22:30:
支持一下新人,鼓励鼓励!做的很好!
其实不建议使用HF。给楼主的建议:DFT计算量不大(略大于HF),但是结论远优于HF,建议使用DFT
tiechong问的也很好。给tiechong的答案:零点能,转动振动能的校正必须做

零点能,转动振动能的校正如何做呢?
15楼2009-08-11 16:50:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yanhualu2000 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见