计算化学公社
标题: 等键反应计算生成焓的一些问题 [打印本页]
作者Author: jianghuaiyu 时间: 2019-5-10 12:18
标题: 等键反应计算生成焓的一些问题
等键反应计算生成焓用的是下面的公式：
ΔH298=ΔE+ΔZPE+ΔHT
高斯优化结果如下：
Zero-point correction= 0.009332 (Hartree/Particle)
Thermal correction to Energy= 0.011692
Thermal correction to Enthalpy= 0.012636
Thermal correction to Gibbs Free Energy= -0.007075
Sum of electronic and zero-point Energies= -100.473051
Sum of electronic and thermal Energies= -100.470691
Sum of electronic and thermal Enthalpies= -100.469747
Sum of electronic and thermal Free Energies= -100.489458
我想问下如下对应关系是否正确：
Zero-point correction对应ZPE
Thermal correction to Enthalpy对应ZPE+HT
Sum of electronic and zero-point Energies对应E+ZPE
Sum of electronic and thermal Enthalpies对应E+ZPE+HT
如果对应关系正确，那么计算时是否只需要用sum of的第三项，即Sum of electronic and thermal Enthalpies代入公式中一次性直接求出ΔE+ΔZPE+ΔHT的值？
作者Author: liyuanhe211 时间: 2019-5-10 12:37
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327
作者Author: jianghuaiyu 时间: 2019-5-10 12:42
这个好像和我问的内容关系不大。。。
作者Author: jianghuaiyu 时间: 2019-5-10 12:47
看了“自己写的一个简易的计算物质气相生成焓的ppt”这个帖子，感觉正确方法是Sum of electronic and zero-point Energies和Sum of electronic and thermal Enthalpies（粗鄙之人算HT。。。）这两项加起来是E+ZPE+HT，这样对吗？
作者Author: liyuanhe211 时间: 2019-5-10 12:50
包含。
作者Author: jianghuaiyu 时间: 2019-5-10 13:10
我看的贴子里计算生成焓时，都没有提到和吉布斯自由能相关的信息，您推荐的网站里写的“(1)巨懒无敌懒人：# B3LYP/6-31G* opt freq，读Sum of electronic and thermal Free Energies=后面的值。”，不懂为什么要加上吉布斯自由能的校正
作者Author: granvia 时间: 2019-5-10 14:18
http://gaussian.com/thermo/ 官网解释得很清楚
作者Author: jianghuaiyu 时间: 2019-5-10 14:53
这个我也看过，按照这上面说的Sum of electronic and thermal enthalpies=Ɛ0+Hcorr=Ɛ0+Etot+KbT，即第三个sum of包含了总的电子能，零点能和温度矫正，我就是想问下能否用所有生成物的这个值的和减去所有反应物的这个值的和来直接计算ΔE+ΔZPE+ΔHT
作者Author: sobereva 时间: 2019-5-10 16:23
别当粗鄙之人
(, 下载次数 Times of downloads: 74)
起码用一般人的做法计算焓，然后再拿产物的减反应物的
作者Author: sobereva 时间: 2019-5-10 16:25
相加毫无意义，大错特错。搞懂焓的热校正量和ZPE是什么关系再去计算
作者Author: sobereva 时间: 2019-5-10 16:25
我博文里说的是自由能的计算
算焓举一反三
作者Author: sobereva 时间: 2019-5-10 16:27
问问题时候把各项都交代清楚，E、HT都是什么东西
作者Author: jianghuaiyu 时间: 2019-5-10 18:50
谢谢老师，现在搞明白了，下次提问会注意
作者Author: davi 时间: 2023-3-2 15:43
老师，看了您的这页PPT，想请教一个问题，就是反应物如果是单质的化，比如C、H2，也需要先opt+freq然后高精度计算单点能吗
作者Author: chands 时间: 2023-3-2 16:07
单原子直接freq就行了，H2当然要opt freq。
作者Author: davi 时间: 2023-3-2 16:21
那氢气还需要在高精度下重新算能量吗，老师
作者Author: chands 时间: 2023-3-2 16:36
当然需要在高级别下计算单点能。
作者Author: davi 时间: 2023-3-3 11:30
等键反应计算生成焓，就是在opt+freq输出文件中读取第三项Thermal correction to Enthalpy，然后在高精度下计算单点能，两者相加，再产物减去生成物吗？
作者Author: chands 时间: 2023-3-3 13:53
是。注意生成焓的零点是一系列单质，你要设计反应将单质和你要算的物质连接起来。比如算CO2的生成焓，要设计反应 C(s, 石墨）+ O2(g) = CO2 (g)
作者Author: davi 时间: 2023-3-3 14:01
那老师，他这样和原子化焓计算不就是一样的吗？感觉要被原子化含和等键反应来求生成焓搞昏了
作者Author: chands 时间: 2023-3-3 14:08
不一样呀，原子化焓是拆成原子的焓变。C(s, 石墨）+ O2(g) = CO2 (g)， O2(g)又不是原子。
作者Author: davi 时间: 2023-3-3 14:29
我看原子化焓来算生成热，一些公式里面又有0K，一些又只有298K，是为啥呢
作者Author: chands 时间: 2023-3-3 14:32
生成焓一般是298 K和 1 bar下的，如果不是，要换算到298 K和1 bar。你找本普通化学瞄一下。
作者Author: davi 时间: 2023-3-3 15:14
老师，我用更高精度来计算单点能，这是输入文件，为什么报错呢
作者Author: chands 时间: 2023-3-3 21:10
关键词写错，写B2PLYPD3/def2tzvp。以后问这种问题要贴出输出文件。
http://sobereva.com/210
http://gaussian.com/basissets/
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |