| 查看: 3173 | 回復(fù): 13 | ||
初衷。新蟲 (初入文壇)
|
[求助]
解離能怎么做 已有2人參與
|
|
|
我在做一個(gè)HgX2(x是鹵素)的解離能 不太清楚具體怎么用高斯來計(jì)算 已經(jīng)做了分子的opt+freq 接下來是不是分成片段Hg+X2來opt+freq然后再算前后的能差呢 求各位前輩大神教教我 都困擾了好幾天 網(wǎng)上也沒看到別人具體怎么做的 發(fā)自小木蟲Android客戶端 |
版主 (著名寫手)
|
conventional解離能可以直接opt算出complex的能量 和每一個(gè)monomer的能量 然后用樓上樓上給的公式BDE=E(complex)-E(monomer1)-E(monomer2)。 也可以考慮BSSE并且用counterpoise方法進(jìn)行校正。 關(guān)鍵詞就是counterpoise input寫法大概就是 # B3LYP/SDD Counterpoise=2 NoSymm Opt HgX + X, optimization with counterpoise correction using ECP basis 0 1 0 2 0 2 ! charge and multiplicity for complex, monomer1, and monomer 2, respectively Hg(Fragment=1) 0.00000000 0.00000000 2.000 X(Fragment=1) 0.00000000 0.00000000 1.000 X(Fragment=2) 0.00000000 0.00000000 3.000 正常算完的話找到輸出文件里面類似這樣的話 Counterpoise corrected energy = xxxxx BSSE energy = xxxxx sum of fragments = xxxxxx complexation energy = xxxxx kcal/mole (raw) complexation energy = xxxxx kcal/mole (corrected) 基本上就能得到解離能了 |
專家顧問 (著名寫手)
新蟲 (初入文壇)
新蟲 (初入文壇)
|
太感謝了 這個(gè)真的特詳細(xì) 我還在sob的帖子里看到另一種感覺是不同的方法 先算整個(gè)分子HgX2的能量 然后減去Hg原子個(gè)X2分子的單點(diǎn)能 完了后兩者想加減分子能量再加上熱化學(xué)修正的量 這樣算是不是也可以 那兩種方法哪種的更可靠呢 發(fā)自小木蟲Android客戶端 |
版主 (著名寫手)
|
零點(diǎn)校正 是應(yīng)該加進(jìn)去的 既然分解前后振動(dòng)差距很大。做法就是 分別opt freq complex和兩個(gè)monomer 得到零點(diǎn)校正以后的能量(zero-point corrected energy) 然后再用公式BDE=E(complex_corr)-E(mono1_corr)-E(mono2_corr)就可以了.這個(gè)和BSSE校正是不同類型的校正 BSSE校正主要是修正計(jì)算中人為加入的誤差 發(fā)自小木蟲IOS客戶端 |
新蟲 (初入文壇)
木蟲 (著名寫手)
木蟲 (著名寫手)
木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|