 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1550 | 回復(fù): 7 | ||
[求助]
NAMD跨膜環(huán)肽管求助，Periodic cell has become too small for original patch grid!
|
|
各位前輩，小妹新手一枚，學(xué)習(xí)過程中問題多多。。。最近遇到一個(gè)棘手的問題，還請(qǐng)大家多多幫忙。我研究的體系是跨膜環(huán)肽納米管，由于要探索納米管在脂質(zhì)膜中的傾斜，所以在模擬的過程中對(duì)脂質(zhì)膜和環(huán)肽管沒有施加任何限制，然后進(jìn)行超長(zhǎng)時(shí)間的MD模擬。 由于計(jì)算條件限制，我是以20納秒為間隔，跑完一次再接著續(xù)跑20納秒。但是實(shí)際操作后發(fā)現(xiàn)，首次二十納秒進(jìn)行順利，在進(jìn)行續(xù)跑，也就是21到40納秒的模擬時(shí)總是出現(xiàn)錯(cuò)誤，錯(cuò)誤提示是Periodic cell has become too small for original patch grid!。 我之前對(duì)照其他帖子給的方法進(jìn)行修改，增加"margin 3",改大”langevinPistonPeriod 1000. langevinPistonDecay 500“但是最后還是不行。我在這里先貼出前二十納秒的CONF和續(xù)跑的CONF文件，希望各位前輩幫忙看看問題出在哪里，多謝。 |
回復(fù)此樓» 猜你喜歡
293求調(diào)劑
已經(jīng)有12人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有4人回復(fù)
-
268求調(diào)劑
已經(jīng)有9人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有5人回復(fù)
-
301求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有9人回復(fù)
-
304求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有5人回復(fù)
-
344求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】namd 出錯(cuò)Periodic cell has become too small for original patch grid!
已經(jīng)有10人回復(fù)
|
structure 8pnt_mem_water.psf coordinates eq.pdb paraTypeCharmm on parameters par_all27_prot_lipid.inp parameters par_extraterms.inp set temperature 310 temperature $temperature if {0} { set inputname eq_r binCoordinates $inputname.restrart.coor binVelocities $inputname.vel ;# remove the "temperature" entry if you use this! #extendedSystem $inputname.xsc } firsttimestep 0 if {1} { cellBasisVector1 58. 0. 0. cellBasisVector2 0. 58. 0. cellBasisVector3 0. 0. 94. cellOrigin 0. 0. 0. } wrapAll on PME yes PMEGridSizeX 64 PMEGridSizeY 64 PMEGridSizeZ 100 exclude scaled1-4 cutoff 10. switching on switchdist 8. pairlistdist 12. timestep 1. nonbondedFreq 2 fullElectFrequency 4 stepspercycle 20 langevin on langevinDamping 5 langevinTemp $temperature langevinPiston on langevinPistonTarget 1. langevinPistonPeriod 1000. langevinPistonDecay 500. langevinPistonTemp $temperature useFlexibleCell yes useConstantRatio yes if {0} { constraints on consref constr.pdb conskfile constr.pdb conskcol B constraintscaling 10.0 } binaryoutput no outputName eq_1 dcdfreq 1000 xstFreq 1000 restartname eq_r_1 restartfreq 1000 binaryrestart yes outputEnergies 1000 outputtiming 1000 if {1} { minimize 10000 reinitvels $temperature } run 10000000 |
2樓2013-11-09 16:34:50
|
接下來是續(xù)跑的CONF structure 8pnt_mem.psf coordinates eq.pdb paraTypeCharmm on parameters par_all27_prot_lipid.inp parameters par_extraterms.inp set temperature 310 if {1} { set inputname eq_r_1 binCoordinates $inputname.coor binVelocities $inputname.vel ;# remove the "temperature" entry if you use this! #extendedSystem $inputname.xsc } firsttimestep 25723000 if {0} { cellBasisVector1 58. 0. 0. cellBasisVector2 0. 58. 0. cellBasisVector3 0. 0. 94. cellOrigin 0. 0. 0. PME yes PMEGridSizeX 64 PMEGridSizeY 64 PMEGridSizeZ 100 } wrapAll on margin 5 exclude scaled1-4 cutoff 10. switching on switchdist 8. pairlistdist 12. timestep 1. nonbondedFreq 2 fullElectFrequency 4 stepspercycle 20 langevin on langevinDamping 5 langevinTemp $temperature langevinPiston on langevinPistonTarget 1. langevinPistonPeriod 1000. langevinPistonDecay 500. langevinPistonTemp $temperature useFlexibleCell yes useConstantRatio yes if {0} { constraints on consref constr.pdb conskfile constr.pdb conskcol B constraintscaling 10.0 } binaryrestart yes binaryoutput no outputName eq_1_xu dcdfreq 1000 xstFreq 1000 restartname eq_r_1_xu restartfreq 1000 outputEnergies 1000 outputtiming 1000 if {0} { minimize 2000 reinitvels $temperature } run 4307000 |
3樓2013-11-09 16:35:49
4樓2013-11-09 18:24:02
木有大神看到么。。。
5樓2013-11-11 09:50:01
|
本帖內(nèi)容被屏蔽 |
6樓2013-11-12 08:45:47
7樓2013-11-12 10:34:16
ruthxu
金蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1109.2
- 散金: 62
- 紅花: 2
- 帖子: 319
- 在線: 138小時(shí)
- 蟲號(hào): 1159022
- 注冊(cè): 2010-11-30
- 性別: MM
- 專業(yè): 高分子科學(xué)
8樓2014-08-05 10:54:18
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yylcx 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277調(diào)劑 +4 | 自由煎餅果子 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 341求調(diào)劑 +4 | 番茄頭--- 2026-03-10 | 4/200 |
|
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|