| 查看: 2357 | 回复: 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 |
回复此楼» 猜你喜欢
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
关于合成MIL-125的求助
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有221人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
投稿咨询
已经有7人回复
-
自科面上 散金币 祈福
已经有61人回复
-
求问：LATP煅烧应该使用什么容器？
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- abinit计算晶体生成焓的问题讨论
已经有7人回复
gongxd325
至尊木虫 (文坛精英)
- 应助: 177 (高中生)
- 金币: 38306.1
- 散金: 13187
- 红花: 109
- 帖子: 34003
- 在线: 1847.4小时
- 虫号: 233962
- 注册: 2006-03-30
- 性别: MM
- 专业: 理论和计算化学
10楼2012-12-19 09:09:36
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
