| 查看: 1771 | 回复: 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人回复
-
求光催化产过氧化氢的PPT ！！！毕业答辩需要
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有244人回复
-
钠离子硬碳负极扣式半电池组装都没有电流！！！
已经有2人回复
-
Ni元素XPS分析
已经有1人回复
-
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生
已经有0人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有0人回复
-
求一份origin2019以上版本的origin软件压缩包
已经有1人回复
-
湖南大学招收锂离子电池、离子液体、高分子与凝胶材料、电催化方向博士生多名
已经有5人回复
gongxd325
至尊木虫 (文坛精英)
- 应助: 177 (高中生)
- 金币: 37560.1
- 散金: 13187
- 红花: 109
- 帖子: 34003
- 在线: 1843.4小时
- 虫号: 233962
- 注册: 2006-03-30
- 性别: MM
- 专业: 理论和计算化学
12楼2012-12-19 09:24:04
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