| 5 | 1/1 | 返回列表 |
| 查看: 1623 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
gn02530640銀蟲(chóng) (小有名氣)
|
[求助]
有沒(méi)有辦法使用 Packmol 取代Material Studio來(lái)縮盒子?(NAMD前處理用)
|
||
|
有沒(méi)有辦法使用 Packmol 取代Material Studio來(lái)縮盒子?(NAMD前處理用,解決初始結(jié)構(gòu)不佳的問(wèn)題) 如題,小弟我目前正在使用NAMD作研究,在進(jìn)行動(dòng)力學(xué)模擬前首先當(dāng)然得先將要模擬的分子和溶劑放入盒子內(nèi)。 在此之前我都是使用Material Studio的 Amorphous Cell Tool 模塊來(lái)達(dá)成,并且也用了縮盒子的方法。關(guān)于縮盒子,按照師兄傳授給我的觀念是:一開(kāi)始將盒子的體積設(shè)得比較大(密度較低),分子會(huì)有比較充裕的空間來(lái)移動(dòng),如此一來(lái)初始結(jié)構(gòu)會(huì)比較正常,之后再慢慢縮小盒子的體積(也就是提高密度或濃度),直到和實(shí)驗(yàn)值吻合為止。 當(dāng)然縮盒子并不是必要的,也是可以在一開(kāi)始直接調(diào)整成和實(shí)驗(yàn)值相同的密度,就把分子放進(jìn)盒子內(nèi)。只是依照我個(gè)人到目前為止的經(jīng)驗(yàn)是,不縮盒子得出的初始結(jié)構(gòu)通常不是很好,在進(jìn)行NAMD模擬前的能量最小化步驟時(shí)就很容易出現(xiàn) ERROR: Constraint failure in RATTLE algorithm for atom 1624! ERROR: Constraint failure; simulation has become unstable. 這樣的錯(cuò)誤( VDW和TOTAL能量值等等都會(huì)變成-9999999999.9999), 即使加大能量最小化的步數(shù)依舊無(wú)法解決。因此在運(yùn)行NAMD動(dòng)力學(xué)模擬前我還是習(xí)慣先縮盒子。 我的問(wèn)題在于:我的體系大約2萬(wàn)多個(gè)原子,受限于個(gè)人計(jì)算器配備的問(wèn)題,每跑一次Amorphous Cell Tool 都要花上數(shù)個(gè)小時(shí),要縮到我想要的密度所耗費(fèi)的時(shí)間實(shí)在太長(zhǎng)。 因此我目前想試著用 Packmol 取代Material Studio來(lái)縮盒子,不過(guò)并不是很順利(有可能是我的方法不對(duì)),想請(qǐng)問(wèn)各位蟲(chóng)友我這個(gè)想法究竟可不可行?非常感謝 |
NAMD |
木蟲(chóng) (正式寫(xiě)手)
|
只要把控壓打開(kāi),盒子就會(huì)自己縮小到合適的密度啊 packmol倒是沒(méi)有什么不妥。。只是我覺(jué)得比較麻煩(不過(guò)我想了想,好像也不是很麻煩) 縮盒子的不好之處在于如果溶劑分子比較復(fù)雜的話,短時(shí)間內(nèi)可能溶劑體系不會(huì)找到最優(yōu)結(jié)構(gòu) 一般的構(gòu)建溶劑盒子的方法是用一個(gè)長(zhǎng)時(shí)間平衡的小盒子通過(guò)周期性復(fù)制得到一個(gè)更大的盒子,這個(gè)可以通過(guò)vmd的solvate模塊實(shí)現(xiàn) |
木蟲(chóng) (正式寫(xiě)手)
銀蟲(chóng) (小有名氣)
|
首先非常感謝您的回復(fù)。 想向你請(qǐng)教一下幾個(gè)問(wèn)題,縮盒子為什么會(huì)是錯(cuò)誤的呢? 又,我是看了這篇 http://www.gaoyang168.com/html/201107/3351593.html 才想試著用PackMol搭建盒子的,請(qǐng)問(wèn)有什么不妥的嗎? 最后,你說(shuō)的「用namd平衡模擬縮盒子」是指逐步縮小cellBasisVector1, cellBasisVector2,cellBasisVector3 的意思嗎? |
銀蟲(chóng) (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +6 | chuanzhu川燭 2026-03-18 | 6/300 |
|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-21 | 4/200 |
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 085600材料與化工 +8 | 安全上岸! 2026-03-16 | 8/400 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|