生成焓的计算问题

CayJ · 发表于 Post on 2020-5-27 21:10:45

sobereva 发表于 2018-6-10 03:21
如果看了这张幻灯片还不知道怎么算的话，那我也没有能帮你的了

sob老师，苯甲醛的标况生成焓没有重复出来。
计算过程：对苯甲醛、H2、O2和C的输入文件中，关键词输入CBS-QB3，计算完成后读取输出文件最后一行CBS-QB3 Enthalpy数据，然后按照图中公式相减，单位换算成kJ/mol后再加上石墨蒸发焓，结果计算出的数据与实验值不在一个数量级，不知道哪里出错了，请老师指点一下！

zjxitcc · 发表于 Post on 2020-5-27 21:25:22

CayJ 发表于 2020-5-27 21:10
sob老师，苯甲醛的标况生成焓没有重复出来。
计算过程：对苯甲醛、H2、O2和C的输入文件中，关键词输入CB ...

可能的原因很多，而你提供的信息太少，让回答者去猜不合适。例如，我就在论坛上回答过，有人把C(g)误当成单重态了。

CayJ · 发表于 Post on 2020-5-28 11:28:47

zjxitcc 发表于 2020-5-27 21:25
可能的原因很多，而你提供的信息太少，让回答者去猜不合适。例如，我就在论坛上回答过，有人把C(g)误当成 ...

我重复10楼PPT中苯甲醛过程如下：
以苯甲醛为例，输入文件如下：

%chk=H:\benjiaquan.chk
%nprocshared=2
#p CBS-QB3

Title Card Required

0 1
C                -3.22007719 -1.53717998 0.00000000
C                -1.82491719 -1.53717998 0.00000000
C                -1.12737919 -0.32942898 0.00000000
C                -1.82503319 0.87908002 -0.00119900
C                -3.21985819 0.87900202 -0.00167800
C                -3.91745919 -0.32920398 -0.00068200
H                -1.27540919 -2.48969298 0.00131500
H                -0.02769919 -0.32934898 0.00063400
H                -1.27483319 1.83122302 -0.00125800
H                -3.76998019 1.83128302 -0.00263100
H                -5.01706319 -0.32902098 -0.00086200
C                -3.99001256 -2.87089627 0.00063022
H                -3.41462096 -3.61513389 -0.50920787
O                -4.22057712 -3.28754892 1.34901465

焓计算结果读取输出文件中的最后一行：CBS-QB3 Enthalpy=

计算结果	电荷数，自旋多重度	H
C7H6O(g) Hartree	0,1	-344.967512
H2(g) Hartree	0,1	-1.162778
O2(g) Hartree	0,3	-150.161299
C(g) Hartree	0,3	-37.783017
H(C7H6O(g))-3H(H2(g))-0.5H(O2(g))-7H(C(g)) Hartree		-1.9174095
单位换算 kJ/mol		-5034.158642
加上石墨蒸发焓后的最终结果 kJ/mol		-29.15864225

还是和PPT中计算结果-37±4.2 kJ/mol相差好多。这是为什么？

zjxitcc · 发表于 Post on 2020-5-28 12:49:58

本帖最后由 zjxitcc 于 2020-5-28 12:58 编辑

CayJ 发表于 2020-5-28 11:28
我重复10楼PPT中苯甲醛过程如下：
以苯甲醛为例，输入文件如下：

你没算错。ppt中的-37±4.2 kJ/mol是实验值，不是计算值。也不存在你所说的“与实验值不在一个数量级”

CayJ · 发表于 Post on 2020-5-31 01:17:14

zjxitcc 发表于 2020-5-27 00:11
你没写清楚HO-N(SO3Na)2的状态是(g), (l)还是(s)

老师您好！我想计算液相的HO-N(SO3Na)2的标况下的生成焓，我看了论坛的相关一些资料，还是有些拿不准，有以下问题还请老师指点一下：
1. 计算方法采用G4是否合适？
2. 计算液相物质的生成焓，在输入文件中输入什么关键词，输出结果是液相生成焓？
3. 设计的模型为0.5N2(g)+2SO2(g)+1.5O2(g)+2Na(s)+0.5H2(g)=HO-N(SO3Na)2(l)，是否合适？

zjxitcc · 发表于 Post on 2020-6-1 11:33:23

本帖最后由 zjxitcc 于 2020-6-1 11:34 编辑

CayJ 发表于 2020-5-31 01:17
老师您好！我想计算液相的HO-N(SO3Na)2的标况下的生成焓，我看了论坛的相关一些资料，还是有些拿不准，有 ...

1.合适
2.“液相生成焓”没办法直接通过简单的关键词获得，同理固相也是。(l)->(g)的焓变（一般称作汽化焓或蒸发焓，或者反过来找液化放热也行）一般通过查询物化手册/数据库得到（就像ppt里的石墨蒸发焓）。计算上可以通过MD模拟获得，耗时较长，精度没法保证，不如找数据来得好。显然，也并不是所有物质的数据都查得到、都有人测过。
3.你确定HO-N(SO3Na)2通常是以液体形式存在？这不是盐么？

CayJ · 发表于 Post on 2020-6-1 13:15:44

zjxitcc 发表于 2020-6-1 11:33
1.合适
2.“液相生成焓”没办法直接通过简单的关键词获得，同理固相也是。(l)->(g)的焓变（一般称作汽化 ...

对，HO-N(SO3Na)2是盐，是要计算一个反应在水溶液下的反应热

zjxitcc · 发表于 Post on 2020-6-1 13:24:29

CayJ 发表于 2020-6-1 13:15
对，HO-N(SO3Na)2是盐，是要计算一个反应在水溶液下的反应热

如果你是盐在水中，应该用(aq)状态。(s)或(l)->(aq)也是有能量变化的。

yaoyuan0711 · 发表于 Post on 2021-6-16 09:54:23

sobereva 发表于 2018-6-10 03:21
如果看了这张幻灯片还不知道怎么算的话，那我也没有能帮你的了

请问sob老师，计算生成焓时，基态S是计算S单原子还是计算S8团簇的能量/8？S单原子是3重态吧？

sobereva · 发表于 Post on 2021-6-16 10:10:34

yaoyuan0711 发表于 2021-6-16 09:54
请问sob老师，计算生成焓时，基态S是计算S单原子还是计算S8团簇的能量/8？S单原子是3重态吧？

S8

是

yaoyuan0711 · 发表于 Post on 2021-6-16 10:22:38

sobereva 发表于 2021-6-16 10:10
S8

是

谢谢sob老师！

wzkchem5 · 发表于 Post on 2021-6-16 14:40:17

yaoyuan0711 发表于 2021-6-16 02:54
请问sob老师，计算生成焓时，基态S是计算S单原子还是计算S8团簇的能量/8？S单原子是3重态吧？

补充一下，算出来的S8是气相结果，需要用实验的升华焓换算成固相S8的焓

yaoyuan0711 · 发表于 Post on 2021-6-16 14:49:54

wzkchem5 发表于 2021-6-16 14:40
补充一下，算出来的S8是气相结果，需要用实验的升华焓换算成固相S8的焓

好的，谢谢！

yaoyuan0711 · 发表于 Post on 2021-7-1 15:51:02

wzkchem5 发表于 2021-6-16 14:40
补充一下，算出来的S8是气相结果，需要用实验的升华焓换算成固相S8的焓

再请教一下您：我在CCCBDB网站上查到的S8(g)生成焓是100.42kJ/mol, S原子(g)是277.17kJ/mol。应该用S8的生成焓除以8？

wzkchem5 · 发表于 Post on 2021-7-1 16:16:24

yaoyuan0711 发表于 2021-7-1 08:51
再请教一下您：我在CCCBDB网站上查到的S8(g)生成焓是100.42kJ/mol, S原子(g)是277.17kJ/mol。应该用S8的 ...

对，要么用S8(g)生成焓不经除以8直接用来换算S8的焓，要么你把计算的S8的焓除以8以后再用除过8的S8(g)生成焓来换算

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 生成焓的计算问题