| 查看: 2517 | 回复: 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 |
回复此楼» 猜你喜欢
厦门大学侯旭团队招聘青年教师、博士后及博硕士研究生
已经有2人回复
-
哈工大深圳 表界面物理化学方向（胶体与界面，光谱探测研究） 博士招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有283人回复
-
元素硫的化合价和氧化态
已经有1人回复
-
江苏科技大学2026年招聘启事
已经有2人回复
-
International Journal of Nanoscience属于中科院4区吗？
已经有0人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
澳大利亚南昆士兰大学量子点课题组招收博士生
已经有1人回复
-
香港中文大学（深圳）2026年招聘公告
已经有0人回复
-
【接测试】TPD-TPR、TPO、脉冲吸附、原位红外，电化学原位红外，站内联系
已经有1人回复
-
北京工业大学2026年公开招聘公告
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- abinit计算晶体生成焓的问题讨论
已经有7人回复
alexwpch
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2566.2
- 散金: 134
- 红花: 3
- 帖子: 451
- 在线: 705.4小时
- 虫号: 672273
- 注册: 2008-12-11
- 性别: GG
- 专业: 应用有机化学
3楼2012-12-19 00:13:50
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
- 专业: 应用有机化学
|
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
alexwpch
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2566.2
- 散金: 134
- 红花: 3
- 帖子: 451
- 在线: 705.4小时
- 虫号: 672273
- 注册: 2008-12-11
- 性别: GG
- 专业: 应用有机化学
5楼2012-12-19 04:11:54
