24小时热门版块排行榜    

查看: 1786  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gongxd325

至尊木虫 (文坛精英)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
alexwpch: 金币+25, ★★★很有帮助 2012-12-20 09:08:07

alexwpch

木虫 (正式写手)

普通回帖

wangth0921

铁虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
alexwpch: 金币+20, ★★★很有帮助 2012-12-20 09:08:41

alexwpch

木虫 (正式写手)

alexwpch

木虫 (正式写手)

alexwpch

木虫 (正式写手)

alexwpch

木虫 (正式写手)

gongxd325

至尊木虫 (文坛精英)

【答案】应助回帖

alexwpch

木虫 (正式写手)

gongxd325

至尊木虫 (文坛精英)

【答案】应助回帖

相关版块跳转 我要订阅楼主 alexwpch 的主题更新
信息提示
请填处理意见