| 5 | 1/1 | 返回列表 |
| 查看: 2788 | 回復(fù): 3 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
gexiaojing木蟲 (著名寫手)
|
[交流]
【求助】熱力學(xué)數(shù)據(jù)計算 已有1人參與
|
|||
|
高斯計算后,熱力學(xué)輸出部分如下,請問如何得到 H、G、U、S 的數(shù)值。空堅敿(xì)點回答,謝謝! 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 |
百寶箱 |

至尊木蟲 (文壇精英)

禁蟲 (文壇精英)
密度泛函·小卒
|
金幣退給你吧,這么簡單的問題居然還要問,好好回去看書去。 http://www.gaoyang168.com/bbs/viewthread.php?tid=1405684 上面這個帖子中講了初學(xué)者應(yīng)該看什么書。 -------------------------- 以下全部轉(zhuǎn)自郭勇《高斯使用指南》第五章第一節(jié)第二小節(jié)第三部分。你看看他對這個問題講的多么詳盡。 ------------------------ 據(jù)說熱化學(xué)性質(zhì)是根據(jù)正則系綜理想氣體展開得到的。在大話西游論壇上看到過這樣一種說法:氣相下計算的能量(SCF 能量)是在真空狀態(tài)下計算的體系的能量,而熱化學(xué)性質(zhì)是在298.15K 一個大氣壓下條件下計算得到。 所有頻率計算都包含有體系的熱化學(xué)分析。默認(rèn)是在298.15K 一個大氣壓下條件下計算得到。對每一種元素都使用含量最高的同位素。下面是熱化學(xué)的輸出部分 - Thermochemistry - Temperature 298.150 Kelvin. Pressure 1.00000 Atm. 計算的具體條件:溫度,氣壓 Atom 1 has atomic number 35 and mass 78.91834 計算所用元素的確切質(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è)定計算條件 Zero-point correction= 0.318636 (Hartree/Particle) 零點能 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 零點能是對分子的電子能量的矯正,表明了在0K 溫度下的分子的振動能量對電子能量的影響。當(dāng)比較OK 時的熱動力學(xué)性質(zhì)時,零點能需要被加到總能量中。同頻率計算一樣,零點能也需要根據(jù)所用的算法采用不同的因子進行矯正?捎肦eadIsotopes 命令實現(xiàn)。據(jù)說優(yōu)化所得的能量是在真空下進行的,所以有必要進行零點能矯正。 在較高的溫度下預(yù)測體系能量時,要給總能加上內(nèi)能項。內(nèi)能包括特定溫度下的平動能,轉(zhuǎn)動能和振動能,這些因素都受到具體條件的影響。注意計算內(nèi)能時已把零點能考慮進去了。不要重復(fù)計算。在考慮內(nèi)能對總能的影響時,要用Readlsotopes 命令指定矯正因子,否則就要自己手動處理,很麻煩。因為內(nèi)能中只有零點能需要矯正,如果自己動手,就需要將內(nèi)能各項分開,算好后再加合到一起。而用命令,程序自動就把一切做好了。 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)在來看幾個公式: E0=Eelec+ZPE 表示零K 時體系的能量等于Ee lec(就是優(yōu)化計算中所得的SCF 能量)和零點能之和。 E= E0+ Evib+ Erot+ Etransl= E0+ E (Thermal) 表示在某一特定溫度和壓力下體系的能量等于零K 時體系的能量加上內(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,是個初學(xué)者,http://www.gaoyang168.com/bbs/viewthread.php?tid=1405684這個網(wǎng)址打不開了,你能推薦一下初學(xué)者看那些書嗎?萬分感激 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求化學(xué)調(diào)劑 +5 | wulanna 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +10 | jhybd 2026-03-23 | 15/750 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +9 | 景隨風(fēng) 2026-03-29 | 14/700 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +10 | 芋泥小鈴鐺 2026-03-28 | 10/500 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 求調(diào)劑 +4 | 青春裁為三截 2026-03-29 | 4/200 |
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運來上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|