| 8 | 1/1 | 返回列表 |
| 查看: 859 | 回復(fù): 7 | ||||
| 當前主題已經(jīng)存檔。 | ||||
tsihu捐助貴賓 (小有名氣)
|
[交流]
[求助] 關(guān)于分子能量的計算與處理
|
|||
|
本人研究的分子有兩種可能構(gòu)象,現(xiàn)在想計算分別的分子能量進行比較以確定哪一種更為穩(wěn)定.感到困擾的是這個分子能量是指構(gòu)型優(yōu)化后的out文件里的SCF Done后的能量還是單點能計算后得到的能量. 還有就是怎么用零點能(zero point energy)進行修正. 望知曉這方面的先行者不吝指教.謝謝. |
仿真建模與計算 |

銀蟲 (小有名氣)
|
既然有兩種可能構(gòu)型,一般都是以這兩種構(gòu)型為初始構(gòu)型,分別進行構(gòu)型優(yōu)化,優(yōu)化后,看最終的構(gòu)型相同還是不相同,如果不相同,那么再比較優(yōu)化后OUT文件中的能量。 如果要進行零點能修正,則在構(gòu)型優(yōu)化后進行頻率計算freq,那么OUT文件中就直接有零點能以及修正后能量的結(jié)果。 建議在計算時,仔細看看Gaussian03手冊中關(guān)于OPT以及Freq頻率計算的相關(guān)內(nèi)容。里面涉及你需要的東西。 |

至尊木蟲 (知名作家)
小木蟲研究生院西南地區(qū)區(qū)長

|
優(yōu)化和單點計算都能給出能量。優(yōu)化對基組和方法的依賴性不是很大,一般用較小基組就可以了,但是能量則不同,往往需要加大基組,比如增加極化函數(shù),彌散函數(shù)。所以LZ的情況,可以直接比較優(yōu)化后的能量(如果可以分清楚的話)也可以在優(yōu)化的基礎(chǔ)上進行更大基組的單點計算已獲得更為準確的能量值 零點能矯正同樓上 頻率計算是為了檢測找到的結(jié)構(gòu)是否為局域極小點,結(jié)構(gòu)是否穩(wěn)定 HOMO-LUMO gap 一般只表示分子的化學(xué)穩(wěn)定性,熱力學(xué)穩(wěn)定性可以計算其鍵能 |
捐助貴賓 (小有名氣)

銀蟲 (小有名氣)
|
各個能量的單位都是Hartree。 構(gòu)型優(yōu)化能否給出相對可靠的分子能量還要看你采用的方法以及基組的大小等因素。進行頻率分析后會得到一些和熱力學(xué)相關(guān)的東西,下面是個具體的例子,看完之后你就能很清楚了。 在分子的熱力學(xué)計算中,包括了零點能的輸出,零點能是對分子的電子能量的矯正,計算了在0K溫度下的分子振動能量.但我們計算的熱力學(xué)數(shù)據(jù)并不總是在OK時,為了計算在較高溫度下的能量,內(nèi)能也要考慮到總能量中,它包括平動能,轉(zhuǎn)動能,和振動能.注意在計算內(nèi)能的時候,已經(jīng)考慮了零點能.下面是一個計算實例 Temperature 298.15 Kelvin. Pressure 1.0000 Atm. Zero-point correction= 0.029201 Thermal correction to Energy= 0.032054 Thermal correction to Enthalpy= 0.032999 Thermal correction to Gibbs Free Energy= 0.008244 Sum of electronic and zero-point Energies= -113.837130 Sum of electronic and thermal Energies= -113.834277 Sum of electronic and thermal Enthalpies= -113.833333 Sum of electronic and thermal Free Energies= -113.858087 后面四行的四個能量分別為E0, E, H, G. 計算為 E0=E(elec) + ZPE E=E0 + E(vib) + E(rol) + E(transl) H=E + RT G=H - TS 上例是在298.15K和1.0Atm下計算的熱力學(xué)狀態(tài)函數(shù),EO即我們的分子體系總能量,而在求298.15K下的內(nèi)能E時,要加上平動能,轉(zhuǎn)動能,和振動能,由統(tǒng)計熱力學(xué)知,這三項處于激發(fā)能級,注意這里的ZPE實際上就是對應(yīng)E(vib)在OK時的能量,但現(xiàn)在并非0K,有關(guān)能量隨溫度變化公式可參考統(tǒng)計熱力著作,最后得到的E就是平常所指的內(nèi)能U,這正是統(tǒng)計熱力學(xué)中對內(nèi)能的定義. |
鐵桿木蟲 (正式寫手)
VASPer

木蟲 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +11 | chaojifeixia 2026-03-19 | 12/600 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | shxchem 2026-03-20 | 5/250 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 一志愿福大288有機化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|