| 8 | 1/1 | 返回列表 |
| 查看: 1882 | 回復: 7 | |||
[求助]
BSSE校正
|
|
我做離解能時用的都是Ezpe值,審稿人建議我替換為free energies evaluating thermal and entropic corrections, 在另外一條意見中提到過最好用BSSE校正。請問,兩條意見是否矛盾?我查到的 E(校正相互作用能)=E(相互作用能)+BSSE,E(相互作用能)=E(AB)-E(AO)-E(BO), 請問E(AB)、E(AO)、E(BO)應(yīng)取什么值?是否可以取審稿人建議的“free energies evaluating thermal and entropic corrections”? 謝謝! |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
1. 對頭。 2. 合理不合理,要看熱力學和熵校正與BSSE的關(guān)系。 首先HF總能,或MP2總能,相當于分子的內(nèi)能U;而熱力學校正,即焓的計算是H=U+PV, 自由能G=H-TS,其中熵S是把分子看作剛體計算轉(zhuǎn)動和平動熵,并在簡諧近似下計算左右的振動熵,另外考慮電子熵。其中振動熵的貢獻占90%以上。進行BSSE,相當于是對焓的內(nèi)能部分進行修正;而內(nèi)能修正不影響其他計算結(jié)果。 當然,振動熵顯然是基于BSSE校正前的能量進行二階微分得到的,這一部分包含BSSE的誤差無法消除,但是這個誤差應(yīng)該小到可以忽略。 因此,我認為合理。 |
|
再問您一個問題哈,如何查找free energies with thermal and entropic corrections? 我在高斯輸出文件中尋找free energies with thermal and entropic corrections,發(fā)現(xiàn)沒有對應(yīng)的,如下: Zero-point correction= 0.037424 (Hartree/Particle) Thermal correction to Energy= 0.046197 Thermal correction to Enthalpy= 0.047141 Thermal correction to Gibbs Free Energy= 0.003273 Sum of electronic and zero-point Energies= -414.287768 Sum of electronic and thermal Energies= -414.278995 Sum of electronic and thermal Enthalpies= -414.278051 Sum of electronic and thermal Free Energies= -414.321919 E (Thermal) CV S KCal/Mol Cal/Mol-Kelvin Cal/Mol-Kelvin Total 28.989 28.058 92.327 Electronic 0.000 0.000 0.000 Translational 0.889 2.981 40.186 Rotational 0.889 2.981 26.662 Vibrational 27.212 22.097 25.479 我在下一部分看到有E (Thermal)和S的total 值 28.989和 92.327 ,請問是這兩個值分別是thermal and entropic 校正值嗎?如果是的話, free energies with thermal and entropic corrections=Etot-28.989-92.327,這樣計算對嗎?謝謝您! |
木蟲 (正式寫手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 生物學一志愿985,分數(shù)349求調(diào)劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 南昌大學材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 一志愿南昌大學,327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學大學 ,080500材料科學與工程學碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|