计算化学公社

标题: 关于热力学组合方法求生成焓与实验值的对比的提问 [打印本页]

作者
Author: XZDZ 时间: 2017-11-16 11:06
标题: 关于热力学组合方法求生成焓与实验值的对比的提问
老师，请问一下，我用G3，G4以及CBS-APNO三种方法算了苯甲醛的生成焓，还把三种方法求出的焓值做了平均处理，并且与ATCT上查得的苯甲醛生成焓的标准值进行了对比，为什么CBS-APNO方法算出来的误差这么大呢？都已经超过了2kcal/mol？谢谢老师解答。

作者
Author: jjjspring 时间: 2017-11-18 00:46
想知道为什么CBS-APNO误差大需要把CBS-APNO中每一项对焓的贡献分开对比，再对比G3，G4方法中各项对焓的贡献，说不定能看出来是那一项导致的偏差。还有苯甲醛在那么高的温度早就分解了，2000K以下的值更现实一点。

作者
Author: XZDZ 时间: 2017-11-18 08:26

jjjspring 发表于 2017-11-18 00:46
想知道为什么CBS-APNO误差大需要把CBS-APNO中每一项对焓的贡献分开对比，再对比G3，G4方法中各项对焓的贡献 ...

谢谢，也就是说需要分别对比G3，G4和CBS-APNO三种方法当中的平动，转动，振动等的配分函数对总配分函数的贡献吗？

作者
Author: jjjspring 时间: 2017-11-18 12:32

XZDZ 发表于 2017-11-18 08:26
谢谢，也就是说需要分别对比G3，G4和CBS-APNO三种方法当中的平动，转动，振动等的配分函数对总配分函数的 ...

平动没差别，转动是根据几何结构来的差别应该比较小，振动都用校正因子估计差的也不大，主要还是电子结构能量

作者
Author: XZDZ 时间: 2017-11-18 18:29

jjjspring 发表于 2017-11-18 12:32
平动没差别，转动是根据几何结构来的差别应该比较小，振动都用校正因子估计差的也不大，主要还是电子结构 ...

好的，非常感谢您的建议。
这个就得需要去查看提出这三种方法的文献了吧？里面应该都会有介绍？

作者
Author: xylz6188 时间: 2018-1-9 14:21
您好，能否具体指导一下生成焓的计算，我按照gaussian教程来计算（http://gaussian.com/thermo/），更不靠谱，计算结果都让我对计算化学产生怀疑了

作者
Author: sobereva 时间: 2018-1-9 14:28

xylz6188 发表于 2018-1-9 14:21
您好，能否具体指导一下生成焓的计算，我按照gaussian教程来计算（http://gaussian.com/thermo/），更不靠 ...

(, 下载次数 Times of downloads: 34)

输入输出文件
(, 下载次数 Times of downloads: 194)

最终结果
deltaH=(-76.393464+1.164715+150.275349/2)*2625.5=-239.1 kcal/mol

作者
Author: xylz6188 时间: 2018-1-9 15:05

sobereva 发表于 2018-1-9 14:28
输入输出文件

谢谢sob老师，我研究一下。

作者
Author: xylz6188 时间: 2018-1-9 15:19
Sob老师您好！按照你给的文件我计算发现和你的有区别，我想应该是在校正方面的问题，能否点播一下。deltaH=(-76.400910+1.165061+150.325202/2)*2625.5=-192.312624 kJ/mol

作者
Author: sobereva 时间: 2018-1-9 15:26

xylz6188 发表于 2018-1-9 15:19
Sob老师您好！按照你给的文件我计算发现和你的有区别，我想应该是在校正方面的问题，能否点播一下。deltaH= ...

仔细看我的输入文件，焓就在输出文件末尾读G4 Enthalpy=

作者
Author: xylz6188 时间: 2018-1-9 15:53

sobereva 发表于 2018-1-9 15:26
仔细看我的输入文件，焓就在输出文件末尾读G4 Enthalpy=

已经明白，谢谢Sob的耐心解答！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)