| 8 | 1/1 | 返回列表 |
| 查看: 1873 | 回復(fù): 7 | ||
樂(lè)逍遙741金蟲(chóng) (小有名氣)
|
[求助]
BSSE校正
|
|
我做離解能時(shí)用的都是Ezpe值,審稿人建議我替換為free energies evaluating thermal and entropic corrections, 在另外一條意見(jiàn)中提到過(guò)最好用BSSE校正。請(qǐng)問(wèn),兩條意見(jiàn)是否矛盾?我查到的 E(校正相互作用能)=E(相互作用能)+BSSE,E(相互作用能)=E(AB)-E(AO)-E(BO), 請(qǐng)問(wèn)E(AB)、E(AO)、E(BO)應(yīng)取什么值?是否可以取審稿人建議的“free energies evaluating thermal and entropic corrections”? 謝謝! |
木蟲(chóng) (正式寫手)
|
不矛盾。BSSE是基組疊加引起的誤差,而free energy則是對(duì)熱效應(yīng)和熵效應(yīng)的校正,完全不同來(lái)源的校正。 E(AB)、E(AO)、E(BO)就是分別優(yōu)化AB, A和B得到的總能。當(dāng)然最好取free energies with thermal and entropic corrections。 |
金蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫手)
|
1. 對(duì)頭。 2. 合理不合理,要看熱力學(xué)和熵校正與BSSE的關(guān)系。 首先HF總能,或MP2總能,相當(dāng)于分子的內(nèi)能U;而熱力學(xué)校正,即焓的計(jì)算是H=U+PV, 自由能G=H-TS,其中熵S是把分子看作剛體計(jì)算轉(zhuǎn)動(dòng)和平動(dòng)熵,并在簡(jiǎn)諧近似下計(jì)算左右的振動(dòng)熵,另外考慮電子熵。其中振動(dòng)熵的貢獻(xiàn)占90%以上。進(jìn)行BSSE,相當(dāng)于是對(duì)焓的內(nèi)能部分進(jìn)行修正;而內(nèi)能修正不影響其他計(jì)算結(jié)果。 當(dāng)然,振動(dòng)熵顯然是基于BSSE校正前的能量進(jìn)行二階微分得到的,這一部分包含BSSE的誤差無(wú)法消除,但是這個(gè)誤差應(yīng)該小到可以忽略。 因此,我認(rèn)為合理。 |
金蟲(chóng) (小有名氣)
|
再問(wèn)您一個(gè)問(wèn)題哈,如何查找free energies with thermal and entropic corrections? 我在高斯輸出文件中尋找free energies with thermal and entropic corrections,發(fā)現(xiàn)沒(méi)有對(duì)應(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 ,請(qǐng)問(wèn)是這兩個(gè)值分別是thermal and entropic 校正值嗎?如果是的話, free energies with thermal and entropic corrections=Etot-28.989-92.327,這樣計(jì)算對(duì)嗎?謝謝您! |
金蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫手)
| 直接取總自由能即可:Sum of electronic and thermal Free Energies= -414.321919 |
金蟲(chóng) (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 5/250 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語(yǔ) 2026-03-20 | 3/150 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|