24小时热门版块排行榜    

查看: 3613  |  回复: 9

秋雨子

金虫 (小有名气)

[交流] 生成焓的计算值与实验值得差别 已有5人参与

请教大家是怎么计算一个化合物在298K,1atm下的生成焓的?
  我通过两种方法计算发现与实验值相差很大,不知道是什么原因,很奇怪
       以HCL 为例,一种方法是用gaussian opt+freq计算HCL,H2,Cl2,取thermochemistry 中Sum of electronic and thermal Enthalpies的值,用HCL-1/2(H2+Cl2)计算求得-12,用HF或者MP2,计算区别不大,仍然和实验值相差悬殊。
       另一种方法是利用gaussian官方网站上 thermochemstry in gaussian 计算求的-32,发现与实验值-90都相差很远,另外发现官网上http://www.gaussian.com/g_whitepap/thermo/thermo.pdf取得例子C2H6的生成焓计算是5.7,可是查出来实验值是-86左右,怎么也相差很远,让我有点匪夷所思,如我我计算有错误,可是怎么解释gaussian自己的例子和实验值还差这么多?还请各位高手指教一二
       (单位都已经换算成KJ/MOL)
回复此楼

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

量化学习资料 Gaussian

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

To the time to life, rather than to life in time
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jintongyin

木虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖
童鞋,标准生成焓不是这样子算滴,你这算出来的只是反应焓,或者连反应焓都不是,是四不像
仔细查两篇文献,好好看看吧,需要用一种等化学键法
2楼2011-12-29 09:57:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lstone03

木虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖
能给个具体计算的例子不?
感觉还不错!
3楼2011-12-29 12:43:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主


小木虫(金币+0.5):给个红包,谢谢回帖
能量计算很依赖方法,而传统的DFT都有这方面的缺陷
我猜你的方法选用错误
Acc Chem Res 2008, 41, 157
Org Biomol Chem 2007, 5, 741
J Chem Phys, 2010, 132, 154104
4楼2011-12-29 14:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋雨子

金虫 (小有名气)

谢谢各位高手的意见,我发现自己有一个错误是kcal/mol和kj/mol的错误,但是gaussian网站上http://muchong.com/bbs/url.php?s=http%3A%2F%2Fwww.gaussian.com%2Fg_whitepap%2Fthermo%2Fthermo.pdf用的原子的方法计算生成焓是不对的,或者他采用的机组不对。
我又用exploring chemistry with electronic structure method P181的等构反应计算了一下,生成焓和实验值相差32kj/mol,是不是有点太大了,我用的am1机组,机组是不太好,可是相差也有点太大,不知道实验和理论相差多少kj/mol的时候可以接受呢?
To the time to life, rather than to life in time
5楼2011-12-30 09:24:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

n370470948

银虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我用dmol模块算出的生成焓居然是正数。。。。不知道你现在知道怎么算了吗?指教下小弟。
6楼2012-05-09 22:53:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

n370470948

银虫 (正式写手)

那标准焓如何算?比如计算乙烷气体,画出其3D 球棍模型,然后呢?用Dmol模块可以算吗?但在勾选properties标签栏中勾选frequency后,计算出来的是正值啊。。。郁闷死了
7楼2012-05-09 22:54:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhudgyg

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by n370470948 at 2012-05-09 22:54:41
那标准焓如何算?比如计算乙烷气体,画出其3D 球棍模型,然后呢?用Dmol模块可以算吗?但在勾选properties标签栏中勾选frequency后,计算出来的是正值啊。。。郁闷死了

偶用castep计算出来的也是正值,郁闷呢,请问您怎么解决的?
8楼2013-03-09 23:54:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

n370470948

银虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
8楼: Originally posted by zhuzhudgyg at 2013-03-09 23:54:19
偶用castep计算出来的也是正值,郁闷呢,请问您怎么解决的?...

最后不了了之了....惭愧啊
9楼2013-03-11 08:49:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzhudgyg

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
9楼: Originally posted by n370470948 at 2013-03-11 08:49:54
最后不了了之了....惭愧啊...

没什么惭愧的,做计算的都很纠结,互相理解吧
10楼2013-03-11 18:59:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 秋雨子 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:24 by vZfe6xYu34yj
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 QUjhNVAcOSff 2026-09-13 7/350 2026-09-14 21:21 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
请填处理意见