| 4 | 1/1 | 返回列表 |
| 查看: 2787 | 回復(fù): 3 | ||||
gexiaojing木蟲 (著名寫手)
|
[交流]
【求助】熱力學(xué)數(shù)據(jù)計(jì)算 已有1人參與
|
|
高斯計(jì)算后,熱力學(xué)輸出部分如下,請(qǐng)問如何得到 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)單的問題居然還要問,好好回去看書去。 http://www.gaoyang168.com/bbs/viewthread.php?tid=1405684 上面這個(gè)帖子中講了初學(xué)者應(yīng)該看什么書。 -------------------------- 以下全部轉(zhuǎn)自郭勇《高斯使用指南》第五章第一節(jié)第二小節(jié)第三部分。你看看他對(duì)這個(gè)問題講的多么詳盡。 ------------------------ 據(jù)說熱化學(xué)性質(zhì)是根據(jù)正則系綜理想氣體展開得到的。在大話西游論壇上看到過這樣一種說法:氣相下計(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 可以通過參數(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ù)說優(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)在來看幾個(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 的差。 以上所有的輸出值無特殊說明外都以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é)者看那些書嗎?萬分感激 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 0703化學(xué) +6 | 妮妮ninicgb 2026-03-27 | 6/300 |
|
|
[考研] 292求調(diào)劑 +13 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 14/700 |
|
|
[考研] 269專碩求調(diào)劑 +10 | 金恩貝 2026-03-21 | 10/500 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 5/250 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|