24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

alexwpch

木虫 (正式写手)

应助: 15 (小学生)
金币: 2566.2
散金: 134
红花: 3
帖子: 451
在线: 705.4小时
虫号: 672273
注册: 2008-12-11
性别: GG
专业: 应用有机化学

[求助] 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

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

abinit计算晶体生成焓的问题讨论已经有7人回复

1楼2012-12-18 03:23:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

alexwpch

木虫 (正式写手)

应助: 15 (小学生)
金币: 2566.2
散金: 134
红花: 3
帖子: 451
在线: 705.4小时
虫号: 672273
注册: 2008-12-11
性别: GG
专业: 应用有机化学

引用回帖:

2楼: Originally posted by wangth0921 at 2012-12-18 08:24:38
C的计算错了，S=1，计算时设：
# G2MP2

C-atom

0,3
C

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相差很大。