| 查看: 2515 | 回复: 14 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
alexwpch木虫 (正式写手)
|
[求助]
G2MP2计算生成焓，误差很大，不知道是不是推算方法有问题，请高手帮忙！
|
|
|
生成焓的计算是根据Thermochemistry in Gaussian算的。 原文公式为thermo-1和therom-2。 简化为G2MP2的公式thermo-3（不知道这个简化是不是正确） 以CH4为例子： H的G2MP2 Temperature= 298.150000 Pressure= 1.000000 E(ZPE)= 0.000000 E(Thermal)= 0.001416 E(QCISD(T))= -0.499810 E(Empiric)= -0.000190 DE(MP2)= 0.000000 G2MP2(0 K)= -0.500000 G2MP2 Energy= -0.498584 G2MP2 Enthalpy= -0.497639 G2MP2 Free Energy= -0.510654 C自旋5的G2MP2 Temperature= 298.150000 Pressure= 1.000000 E(ZPE)= 0.000000 E(Thermal)= 0.001416 E(QCISD(T))= -37.623027 E(Empiric)= -0.000760 DE(MP2)= -0.005225 G2MP2(0 K)= -37.629012 G2MP2 Energy= -37.627596 G2MP2 Enthalpy= -37.626652 G2MP2 Free Energy= -37.644040 C自旋1的G2MP2 Temperature= 298.150000 Pressure= 1.000000 E(ZPE)= 0.000000 E(Thermal)= 0.001416 E(QCISD(T))= -37.708117 E(Empiric)= -0.010000 DE(MP2)= -0.016821 G2MP2(0 K)= -37.734938 G2MP2 Energy= -37.733522 G2MP2 Enthalpy= -37.732578 G2MP2 Free Energy= -37.748446 CH4的G2MP2 Temperature= 298.150000 Pressure= 1.000000 E(ZPE)= 0.042658 E(Thermal)= 0.045529 E(QCISD(T))= -40.405888 E(Empiric)= -0.020000 DE(MP2)= -0.026434 G2MP2(0 K)= -40.409664 G2MP2 Energy= -40.406793 G2MP2 Enthalpy= -40.405849 G2MP2 Free Energy= -40.426964 元素的焓值deltaH(0K)采用J. Chem. Phys. 106, 1063 (1997)，见图thermo-4 这样算出来的deltaH（CH4,298K)=-191kJ/mol，跟J. Phys. Chem. 99, 17145 (1995)给出的-74kJ/mol的实验和G2MP2值差很多额。 同样的方法deltaH（NH2NH2,298K)=84.7kJ/mol和维基上查到的50.1kJ/mol也差值很大。 肯定是哪边出问题了，请教下各位高手如何解决额？  thermo-4.JPG  thermo-1.JPG  thermo-2.JPG  thermo-3.JPG |
回复此楼» 猜你喜欢
压裂返排液回用配液：不是处理达标再回用，是按配液要求控制指标
已经有0人回复
-
高级氧化技术怎么选：羟基自由基的"虚假繁荣"与工程现实
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有289人回复
-
苏州国家实验室依托苏州大学，面向2027年秋季入学招收博士生
已经有0人回复
-
生物脱氮的碳源账：甲醇、乙酸钠、葡萄糖怎么选，以及内碳源到底值不值得挖
已经有0人回复
-
含氯化工废水的电催化：活性氯是免费氧化剂还是麻烦制造者？
已经有0人回复
-
求助《生物化学》第三版的章节习题和答案完整版
已经有0人回复
-
求助《生物化学》第三版的章节习题和答案完整版
已经有0人回复
-
石墨电极在化工浓水零排前端的角色：不是主战，但少它不行
已经有0人回复
-
【文献导读】Chem. Rev. 绿氢特刊：1美元/kg目标与SOEC等电解路线要点
已经有0人回复
-
析氧电位Tafel 斜率ECSA这三个看起来很专业的参数,可能一直在骗你
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- abinit计算晶体生成焓的问题讨论
已经有7人回复
alexwpch
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2566.2
- 散金: 134
- 红花: 3
- 帖子: 451
- 在线: 705.4小时
- 虫号: 672273
- 注册: 2008-12-11
- 性别: GG
- 专业: 应用有机化学
9楼2012-12-19 08:55:56
wangth0921
铁虫 (职业作家)
- 应助: 117 (高中生)
- 金币: 1359
- 散金: 200
- 红花: 22
- 帖子: 3740
- 在线: 787.1小时
- 虫号: 1426386
- 注册: 2011-10-03
- 专业: 环境污染化学
2楼2012-12-18 08:24:38
alexwpch
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2566.2
- 散金: 134
- 红花: 3
- 帖子: 451
- 在线: 705.4小时
- 虫号: 672273
- 注册: 2008-12-11
- 性别: GG
- 专业: 应用有机化学
3楼2012-12-19 00:13:50
alexwpch
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2566.2
- 散金: 134
- 红花: 3
- 帖子: 451
- 在线: 705.4小时
- 虫号: 672273
- 注册: 2008-12-11
- 性别: GG
- 专业: 应用有机化学
|
N的0，2算出来 Temperature= 298.150000 Pressure= 1.000000 E(ZPE)= 0.000000 E(Thermal)= 0.001416 E(QCISD(T))= -54.385395 E(Empiric)= -0.010190 DE(MP2)= -0.025526 G2MP2(0 K)= -54.421111 G2MP2 Energy= -54.419695 G2MP2 Enthalpy= -54.418750 G2MP2 Free Energy= -54.435492 N的0，4算出来 Temperature= 298.150000 Pressure= 1.000000 E(ZPE)= 0.000000 E(Thermal)= 0.001416 E(QCISD(T))= -54.491420 E(Empiric)= -0.005570 DE(MP2)= -0.019315 G2MP2(0 K)= -54.516306 G2MP2 Energy= -54.514889 G2MP2 Enthalpy= -54.513945 G2MP2 Free Energy= -54.531341 算出deltaH（NH2NH2,298K)为-403.9和84.7kJ/mol两个，均与50.1kJ/mol相差很大。 |
4楼2012-12-19 00:23:25
