| 3 | 1/1 | 返回列表 |
| 查看: 589 | 回復(fù): 2 | |||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | |||
hljiang鐵桿木蟲 (著名寫手)
|
[交流]
【求助】能用CPMD判斷一個(gè)小分子能否進(jìn)入孔狀化合物孔道嗎? 已有1人參與
|
||
|
我是做實(shí)驗(yàn)的,我想請教一下: 一個(gè)孔狀化合物,孔道本身和小分子尺寸大小相當(dāng),所以單純從分子尺寸上來說很難判別小分子能否進(jìn)入孔道內(nèi)部,我想通過CPMD計(jì)算來判斷小分子能否進(jìn)入孔道,可以嗎? 哪位能幫個(gè)忙啊? 論文已經(jīng)整理好了,就差這一點(diǎn).我只要判斷小分子能否進(jìn)入孔道,不需要其他結(jié)果. 以合作貢獻(xiàn)論文的形式也可以,送金幣也可以. 多謝先! |

木蟲 (正式寫手)
|
(1) 如果只是考察能否進(jìn)入孔道,問題可以描述為:把分子放到孔道內(nèi)所得到的能量,與分子在孔道外的情況下的能量,是升高,還是降低?如果進(jìn)入孔道后能量升高,且高的多,那么說明孔道禁止分子進(jìn)入,反之,允許進(jìn)入。 這樣的話,只要做一系列能量計(jì)算,得到分子從孔道外進(jìn)入孔道內(nèi)的能量變化曲線。顯示一下能量變化問題即可。 困難點(diǎn)是如何確定模型體系的孔道幾何和分子尺寸。由于分子和孔道原子的作用會(huì)改變它們各自的結(jié)構(gòu),因此,每一步都要優(yōu)化整個(gè)體系的幾何結(jié)構(gòu)。計(jì)算會(huì)比較慢些。 (2) 考慮另外一種情況,假如孔道是柔性的,有不止一種可能的結(jié)構(gòu), 并且有的結(jié)構(gòu)允許分子進(jìn)入,有的不允許。那么上面的計(jì)算就不能全面回答這個(gè)問題。孔道尺寸在特定的溫度下會(huì)會(huì)由于原子熱運(yùn)動(dòng)不斷弛豫。這種情況下,最好要跑MD,計(jì)算一下分子進(jìn)入孔道的自由能刨面(Profile of free energy, 或者Potential of mean force, PMF)。就是把分子用一個(gè)constrain force固定到孔道從外到內(nèi)的一系列點(diǎn)上,分別跑MD,完了后WHAM公式計(jì)算整個(gè)過程的自由能壘,和平均進(jìn)入時(shí)間。如果自由能壘很高,意味著分子不能進(jìn)入,能壘低則能夠進(jìn)入。 和上面的能量計(jì)算相比,自由能計(jì)算的是一種平均的,統(tǒng)計(jì)意義的結(jié)果。這樣的計(jì)算可以充分考慮孔道,或者底物分子的幾何柔性,和熱運(yùn)動(dòng)的影響。 (3) 特別需要注意的是,第一種計(jì)算中,我們假定底物分子是氣態(tài)的。假如你的體系中,底物分子是凝聚態(tài)而不是氣態(tài),那么底物分子之間的作用對進(jìn)入孔道的影響更大,那就必須使用MD的方法,因?yàn)榘吹谝环N方法計(jì)算得到的勢能面可能是錯(cuò)誤的,因?yàn)樗荒苷_反映分子在孔道外的自由能高低。 |
鐵桿木蟲 (著名寫手)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么! 2026-03-15 | 19/950 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|