| 查看: 3684 | 回復(fù): 19 | ||||
sesy木蟲 (小有名氣)
|
[求助]
請(qǐng)教高斯09算溶劑效應(yīng)及能量校正問題,多送金幣+小紅花
|
|||
我們這級(jí)新接觸Gaussion09,看小木蟲中有很多帖子說09和03做能量校正得到的輸出結(jié)果不一樣, 所以請(qǐng)教高手看看我用09算的方法對(duì)不對(duì),尤其是輸入文件寫的有問題嗎?我先優(yōu)化完結(jié)構(gòu)然后做單點(diǎn)能和溶劑效應(yīng),輸入文件如下,請(qǐng)之前用過09的蟲友們多多指點(diǎn)。 多送金幣+小紅花 希望大家多多關(guān)注,急求解答!!一、使用G09程序包,先用小基組b3lyp/6-31g(d)對(duì)構(gòu)型進(jìn)行優(yōu)化和頻率分析, 輸入文件關(guān)鍵詞 :# opt freq b3lyp/6-31g(d) geom=connectivity 輸出結(jié)果選用的 Thermal correction to Gibbs Free Energy= (用A表示) 二、再用大基組計(jì)算單點(diǎn)能:# b3lyp/6-311++g(d,p) scf=tight 輸出結(jié)果選用的 SCF Done: E(RB3LYP) = -1246.14645414(用B表示) 三、最后考慮到溶劑效應(yīng): 輸入文件# b3lyp/6-311++g(d,p) scf=tight scrf=(PCM,solvent=benzene,read,smd,Dovacuum,ExternalIteration) 。 。 。 Radii=uaks (最后) 輸出結(jié)果選用的 DeltaG (solv) (kcal/mol) = -38.40 (用C表示) 最后自由能表示:自由能=自由能校正(A)+單點(diǎn)能(B)+溶劑化能(C) 請(qǐng)大家看看我這樣做對(duì)不對(duì),還有不對(duì)的話應(yīng)該怎樣寫輸入文件?多謝大家。。。 |
過渡態(tài)計(jì)算 | 高斯 |

金蟲 (小有名氣)
|
一、使用G09程序包,先用小基組b3lyp/6-31g(d)對(duì)構(gòu)型進(jìn)行優(yōu)化和頻率分析, 輸入文件關(guān)鍵詞 :# opt freq b3lyp/6-31g(d) geom=connectivity 輸出結(jié)果選用的 Thermal correction to Gibbs Free Energy= (用A表示) 二、輸入文件# b3lyp/6-311++g(d,p) scf=tight scrf=(solvent=benzene,read,smd,Dovacuum) 。 。 。 Radii=uaks (最后) 輸出結(jié)果選用的 DeltaG (solv) (kcal/mol) = -38.40 (用C表示) 最后自由能表示:自由能=自由能校正(A)+(C) |
金蟲 (小有名氣)
鐵桿木蟲 (正式寫手)
木蟲 (小有名氣)

木蟲 (小有名氣)

木蟲 (小有名氣)

木蟲 (小有名氣)

木蟲 (正式寫手)
|
第二步 二、再用大基組計(jì)算單點(diǎn)能:# b3lyp/6-311++g(d,p) scf=tight 輸出結(jié)果選用的 SCF Done: E(RB3LYP) = -1246.14645414(用B表示) 是沒有必要算的,只算一三步就行。 最后自由能表示:自由能=自由能校正(A)+溶劑化下的單點(diǎn)能 溶劑化下的單點(diǎn)能是g09溶劑化輸出結(jié)果中的SCF Done: E(UB3LYP) = 這一項(xiàng)的能量 |
金蟲 (小有名氣)
木蟲 (小有名氣)

銅蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +7 | mikasa的圍巾 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-28 | 4/200 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 求化學(xué)調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|