生成焓的计算问题

CayJ

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

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

zjxitcc

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

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

CayJ

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

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

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

CayJ

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

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

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

CayJ

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

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

zjxitcc

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

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

yaoyuan0711

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

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

sobereva

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

S8

是

yaoyuan0711

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

是

谢谢sob老师！

wzkchem5

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

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

yaoyuan0711

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

好的，谢谢！

yaoyuan0711

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

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

wzkchem5

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)生成焓来换算