| 5 | 1/1 | 返回列表 |
| 查看: 1621 | 回復(fù): 4 | |||
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ě)手)
銀蟲(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) (正式寫(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) (小有名氣)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過(guò) +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +3 | vv迷 2026-03-22 | 3/150 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|