版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

[求助] G2MP2計算生成焓，誤差很大，不知道是不是推算方法有問題，請高手幫忙！

生成焓的計算是根據(jù)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

回復此樓

» 猜你喜歡

材料/化學相關專業(yè)2026級學術(shù)/專業(yè)型碩士研究生已經(jīng)有19人回復
求催化作用基礎第三版（甄開吉) 課件已經(jīng)有5人回復
物理化學論文潤色/翻譯怎么收費? 已經(jīng)有213人回復
26年材料/化學相關專業(yè)博士研究生招生溫州大學化學與材料工程學院(26年秋季入學) 已經(jīng)有0人回復
求助電催化劑制備，怎么在碳布表面生長PdCo合金已經(jīng)有0人回復
求助電催化劑制備，怎么在碳布表面生長PdCo合金已經(jīng)有1人回復
304求調(diào)劑已經(jīng)有3人回復
上海華東理工大學機械與動力工程學院先進材料與制造研究組招聘公告已經(jīng)有0人回復
304求調(diào)劑已經(jīng)有6人回復
304求調(diào)劑已經(jīng)有1人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

abinit計算晶體生成焓的問題討論已經(jīng)有7人回復

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

gongxd325

至尊木蟲 (文壇精英)

應助: 177 (高中生)
金幣: 37753.1
散金: 13187
紅花: 109
帖子: 34003
在線: 1844.8小時
蟲號: 233962
注冊: 2006-03-30
性別: MM
專業(yè): 理論和計算化學

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
alexwpch: 金幣+25, ★★★很有幫助 2012-12-20 09:08:07

(1) 建議不要用原子化能求生成焓
(2) 自旋多重度根據(jù)原子或分子的電子排布情況決定，無不成對電子為1，有1個不成對電子為2，2個不成對電子為3

贊一下(1人)

回復此樓

一個人做點好事并不難難的是一輩子做好事!!!

6樓2012-12-19 07:53:02

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

引用回帖:

12樓: Originally posted by gongxd325 at 2012-12-19 09:24:04
我的意思你沒完全理解，你用G2計算等鍵反應得到生成焓結(jié)果更好，原子化反應只是作為對照。C和O都有2個不成對電子，多重度都是3，N有3個不成對電子，多重度是4

N的多重度是4的情況下算出來的NH2NH2的生成焓是84.7kJ/mol與查到的50kJ/mol相差比較大。

贊一下(1人)

回復此樓

13樓2012-12-19 09:35:01

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

wangth0921

鐵蟲 (職業(yè)作家)

應助: 117 (高中生)
金幣: 1359
散金: 200
紅花: 22
帖子: 3740
在線: 787.1小時
蟲號: 1426386
注冊: 2011-10-03
專業(yè): 環(huán)境污染化學

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
alexwpch: 金幣+20, ★★★很有幫助 2012-12-20 09:08:41

C的計算錯了，S=1，計算時設：
# G2MP2

C-atom

0,3
C

贊一下

回復此樓

2樓2012-12-18 08:24:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

引用回帖:

2樓: Originally posted by wangth0921 at 2012-12-18 08:24:38
C的計算錯了，S=1，計算時設：
# G2MP2

C-atom

0,3
C

那我的NH2NH2應該也是N的自旋多重度錯了？N的自旋多重度應該是多少？輸入寫0,2還是0,4？

贊一下

回復此樓

3樓2012-12-19 00:13:50

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

引用回帖:

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

贊一下

回復此樓

4樓2012-12-19 00:23:25

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

引用回帖:

2樓: Originally posted by wangth0921 at 2012-12-18 08:24:38
C的計算錯了，S=1，計算時設：
# G2MP2

C-atom

0,3
C

還有O的應該是多少，自旋多重度為3？

贊一下

回復此樓

5樓2012-12-19 04:11:54

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

引用回帖:

6樓: Originally posted by gongxd325 at 2012-12-19 07:53:02
(1) 建議不要用原子化能求生成焓
(2) 自旋多重度根據(jù)原子或分子的電子排布情況決定，無不成對電子為1，有1個不成對電子為2，2個不成對電子為3

那用G2MP2算的結(jié)果，如何來求生成焓？
我用b3pw91/6-311++g(3df,3pd)構(gòu)建等鍵反應算的生成焓，審稿人說我的體系比較小，直接用G2級別的算。

贊一下

回復此樓

7樓2012-12-19 08:09:33

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

gongxd325

至尊木蟲 (文壇精英)

應助: 177 (高中生)
金幣: 37753.1
散金: 13187
紅花: 109
帖子: 34003
在線: 1844.8小時
蟲號: 233962
注冊: 2006-03-30
性別: MM
專業(yè): 理論和計算化學

【答案】應助回帖

既然b3pw91/6-311++g(3df,3pd)計算等鍵反應會用，為何G2不會用？不過就是改用G2計算焓變

贊一下

回復此樓

一個人做點好事并不難難的是一輩子做好事!!!

8樓2012-12-19 08:52:03

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

alexwpch

木蟲 (正式寫手)

應助: 15 (小學生)
金幣: 2566.2
散金: 134
紅花: 3
帖子: 451
在線: 705.4小時
蟲號: 672273
注冊: 2008-12-11
性別: GG
專業(yè): 應用有機化學

引用回帖:

8樓: Originally posted by gongxd325 at 2012-12-19 08:52:03
既然b3pw91/6-311++g(3df,3pd)計算等鍵反應會用，為何G2不會用？不過就是改用G2計算焓變

審稿人原話如下：
“Why use an isodesmic procedure to calculate heats of formation? These molecules are small enough that it could be done by high level ab initio methods such as G2 or CCSD(T).”

贊一下

回復此樓

9樓2012-12-19 08:55:56

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

gongxd325

至尊木蟲 (文壇精英)

應助: 177 (高中生)
金幣: 37753.1
散金: 13187
紅花: 109
帖子: 34003
在線: 1844.8小時
蟲號: 233962
注冊: 2006-03-30
性別: MM
專業(yè): 理論和計算化學

【答案】應助回帖

等鍵反應計算生成焓精度高，是很好的辦法，你可以不用DFT而用G2計算，作為比較，你也可以用G2和原子化能計算（這就是審稿的意思），也就是借助于原子的實驗生成焓計算你的分子的生成焓

贊一下

回復此樓

一個人做點好事并不難難的是一輩子做好事!!!

10樓2012-12-19 09:09:36

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉(zhuǎn) 我要訂閱樓主 alexwpch 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 機械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 一志愿中國石油大學（華東）本科齊魯工業(yè)大學 +3	石能偉 2026-03-17	3/150	2026-03-21 02:22 by JourneyLucky
[考研] 279分求調(diào)劑一志愿211 +11	chaojifeixia 2026-03-19	12/600	2026-03-21 01:49 by 星空星月
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 北科281學碩材料求調(diào)劑 +5	tcxiaoxx 2026-03-20	5/250	2026-03-20 21:35 by laoshidan
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 281求調(diào)劑（0805） +14	煙汐憶海 2026-03-16	25/1250	2026-03-20 15:47 by yuncha
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 0817調(diào)劑 +3	沒有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 302求調(diào)劑 +4	小賈同學123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學123
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879

24小時熱門版塊排行榜

[求助] G2MP2計算生成焓，誤差很大，不知道是不是推算方法有問題，請高手幫忙！

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

[求助] G2MP2計算生成焓，誤差很大，不知道是不是推算方法有問題，請高手幫忙！

» 本主題相關價值貼推薦，對您同樣有幫助: