| 8 | 1/1 | 返回列表 |
| 查看: 3961 | 回復(fù): 7 | |||
| 本帖產(chǎn)生 1 個 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | |||
zhaoyunm木蟲 (著名寫手)
|
[求助]
金屬配合物穩(wěn)定常數(shù)分子模擬計(jì)算
|
||
| 欲合成一新螯合劑分子,想先模擬計(jì)算該分子與某金屬離子的螯合配位能力(金屬配合物穩(wěn)定常數(shù)),由于沒有學(xué)習(xí)過相關(guān)方面的知識,對這方面一竅不通,請問用什么軟件計(jì)算,該如何計(jì)算?麻煩這方面的高手、前輩們幫幫忙,越詳細(xì)越好,感激不盡! |
精華 |

木蟲 (正式寫手)
|
金屬配合物穩(wěn)定常數(shù)是個熱力學(xué)平衡常數(shù),表示配合反應(yīng)平衡后,配合物濃度與未配合的化合物濃度之間的比例。那么這個平衡比例由什么決定的呢?當(dāng)然是配合前后的自由能差決定的。 分子模擬模擬一個反應(yīng)的自由能差通常有2種方法:cluster方法和MD方法。 1. Cluster方法是用gauss計(jì)算配合能。 你先計(jì)算配合物AB的能量,再減去鰲合物A和配體B的能量就好了。當(dāng)然,三個能量都要做結(jié)構(gòu)優(yōu)化,AB復(fù)合物還要做BSSE校正。然而,這樣得到的是結(jié)合能而不是結(jié)合自由能,怎么辦呢?不要緊,Ab, A, B優(yōu)化完了以后做頻率計(jì)算,會得到每個種類的自由能,用自由能之差即可。那要是這個體系是在水里面而不是氣相Cluster呢?也好辦,結(jié)構(gòu)優(yōu)化和頻率計(jì)算時加用隱式溶劑模型PCM。 2. MD方法 Cluster方法雖好,但隱式溶劑模型有多精確呢?不好說,有的體系很好,有的體系不太好。為什么呢?這要看配體,鰲合物,復(fù)合物與溶劑水之間的作用有多強(qiáng)烈,有多復(fù)雜了。特別是如果成氫鍵,或水分子直接參加配位鍵,那俺老人家就擔(dān)心cluster模型的準(zhǔn)確性了。不要緊,我們可以用顯式溶劑模型,即把鰲合物A和配體B放到一盒真實(shí)的水分子中去跑動力學(xué),讓他們自由結(jié)合到一起,然后通過Kirkwood公式計(jì)算出來整個配合過程的自由能勢能面。然后配合前后的自由能差就能從這個勢能面中直接讀出來。Perfect! 但是,慢著,這個模擬要跑多久?這取決于配合自由能差和配合反應(yīng)的能壘高低,能壘越高越慢。6kcal/mol的能壘在300K溫度下去沖擊能壘,一般要一般2個ns,還不見得一定能沖上去,還要看你初始構(gòu)型。此外,沖上去前后,都需要足夠的采樣來計(jì)算平均勢能,才好由Boltzmann算概率啊,有了概率才能算自由能啊?磥碇苯幽M結(jié)合自由能不行,怎么辦呢?得用加速采樣方法,這就是Umbrella sampling+PMF或用metadynamics。就是這么些了,慢慢玩吧。 |
至尊木蟲 (著名寫手)
木蟲 (正式寫手)
木蟲 (著名寫手)

金蟲 (小有名氣)
銀蟲 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|