| 24小時(shí)熱門版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 2785 | 回復(fù): 3 | ||||
gexiaojing木蟲 (著名寫手)
|
[交流]
【求助】熱力學(xué)數(shù)據(jù)計(jì)算 已有1人參與
|
|
高斯計(jì)算后,熱力學(xué)輸出部分如下,請(qǐng)問(wèn)如何得到 H、G、U、S 的數(shù)值?請(qǐng)?jiān)敿?xì)點(diǎn)回答,謝謝! Zero-point correction= 0.099061 (Hartree/Particle) Thermal correction to Energy= 0.104208 Thermal correction to Enthalpy= 0.105152 Thermal correction to Gibbs Free Energy= 0.070422 Sum of electronic and zero-point Energies= -306.393552 Sum of electronic and thermal Energies= -306.388405 Sum of electronic and thermal Enthalpies= -306.387461 Sum of electronic and thermal Free Energies= -306.422191 E (Thermal) CV S KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin Total 65.391 18.370 73.094 Electronic 0.000 0.000 0.000 Translational 0.889 2.981 39.270 Rotational 0.889 2.981 25.997 Vibrational 63.614 12.408 7.828 Vibration 1 0.618 1.903 2.677 Vibration 2 0.640 1.832 2.085 Vibration 3 0.842 1.281 0.694 Vibration 4 0.887 1.178 0.582 |
百寶箱 |

禁蟲 (文壇精英)
密度泛函·小卒
|
金幣退給你吧,這么簡(jiǎn)單的問(wèn)題居然還要問(wèn),好好回去看書去。 http://www.gaoyang168.com/bbs/viewthread.php?tid=1405684 上面這個(gè)帖子中講了初學(xué)者應(yīng)該看什么書。 -------------------------- 以下全部轉(zhuǎn)自郭勇《高斯使用指南》第五章第一節(jié)第二小節(jié)第三部分。你看看他對(duì)這個(gè)問(wèn)題講的多么詳盡。 ------------------------ 據(jù)說(shuō)熱化學(xué)性質(zhì)是根據(jù)正則系綜理想氣體展開得到的。在大話西游論壇上看到過(guò)這樣一種說(shuō)法:氣相下計(jì)算的能量(SCF 能量)是在真空狀態(tài)下計(jì)算的體系的能量,而熱化學(xué)性質(zhì)是在298.15K 一個(gè)大氣壓下條件下計(jì)算得到。 所有頻率計(jì)算都包含有體系的熱化學(xué)分析。默認(rèn)是在298.15K 一個(gè)大氣壓下條件下計(jì)算得到。對(duì)每一種元素都使用含量最高的同位素。下面是熱化學(xué)的輸出部分 - Thermochemistry - Temperature 298.150 Kelvin. Pressure 1.00000 Atm. 計(jì)算的具體條件:溫度,氣壓 Atom 1 has atomic number 35 and mass 78.91834 計(jì)算所用元素的確切質(zhì)量 Atom 2 has atomic number 16 and mass 31.97207 Atom 3 has atomic number 8 and mass 15.99491 可以通過(guò)參數(shù)ReadIsotopes(見Freq 關(guān)鍵詞,矯正因子可以加也可以不加)根據(jù)自己 的具體需要設(shè)定計(jì)算條件 Zero-point correction= 0.318636 (Hartree/Particle) 零點(diǎn)能 Thermal correction to Energy= 0.342790 Thermal correction to Enthalpy= 0.343734 Thermal correction to Gibbs Free Energy= 0.260712 Sum of electronic and zero-point Energies= -4014.081556 Sum of electronic and thermal Energies= -4014.057402 Sum of electronic and thermal Enthalpies= -4014.056458 Sum of electronic and thermal Free Energies= -4014.139480 零點(diǎn)能是對(duì)分子的電子能量的矯正,表明了在0K 溫度下的分子的振動(dòng)能量對(duì)電子能量的影響。當(dāng)比較OK 時(shí)的熱動(dòng)力學(xué)性質(zhì)時(shí),零點(diǎn)能需要被加到總能量中。同頻率計(jì)算一樣,零點(diǎn)能也需要根據(jù)所用的算法采用不同的因子進(jìn)行矯正?捎肦eadIsotopes 命令實(shí)現(xiàn)。據(jù)說(shuō)優(yōu)化所得的能量是在真空下進(jìn)行的,所以有必要進(jìn)行零點(diǎn)能矯正。 在較高的溫度下預(yù)測(cè)體系能量時(shí),要給總能加上內(nèi)能項(xiàng)。內(nèi)能包括特定溫度下的平動(dòng)能,轉(zhuǎn)動(dòng)能和振動(dòng)能,這些因素都受到具體條件的影響。注意計(jì)算內(nèi)能時(shí)已把零點(diǎn)能考慮進(jìn)去了。不要重復(fù)計(jì)算。在考慮內(nèi)能對(duì)總能的影響時(shí),要用Readlsotopes 命令指定矯正因子,否則就要自己手動(dòng)處理,很麻煩。因?yàn)閮?nèi)能中只有零點(diǎn)能需要矯正,如果自己動(dòng)手,就需要將內(nèi)能各項(xiàng)分開,算好后再加合到一起。而用命令,程序自動(dòng)就把一切做好了。 E (Thermal) CV S KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin Total 215.104 90.314 174.735 Electronic 0.000 0.000 0.000 Translational 0.889 2.981 44.059 Rotational 0.889 2.981 37.244 Vibrational 213.326 84.353 93.432 內(nèi)能恒容熱容熵 現(xiàn)在來(lái)看幾個(gè)公式: E0=Eelec+ZPE 表示零K 時(shí)體系的能量等于Ee lec(就是優(yōu)化計(jì)算中所得的SCF 能量)和零點(diǎn)能之和。 E= E0+ Evib+ Erot+ Etransl= E0+ E (Thermal) 表示在某一特定溫度和壓力下體系的能量等于零K 時(shí)體系的能量加上內(nèi)能值。 H=E+RT 表示在某一特定溫度和壓力下體系的焓等于在這一條件下的能量和RT 之和。 G=H-TS 表示在某一特定溫度和壓力下體系的自由能等于在這一條件下的焓和TS 的差。 以上所有的輸出值無(wú)特殊說(shuō)明外都以HARTREE 為單位。 [ Last edited by yjcmwgk on 2009-8-13 at 16:10 ] |
至尊木蟲 (文壇精英)

金蟲 (小有名氣)
|
您好,我想學(xué)習(xí)thermal cal,是個(gè)初學(xué)者,http://www.gaoyang168.com/bbs/viewthread.php?tid=1405684這個(gè)網(wǎng)址打不開了,你能推薦一下初學(xué)者看那些書嗎?萬(wàn)分感激 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +12 | hhhhxn.. 2026-03-23 | 18/900 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 22408 359分調(diào)劑 +3 | Qshers 2026-03-27 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|