| 查看: 1583 | 回复: 7 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
NAMD跨膜环肽管求助，Periodic cell has become too small for original patch grid!
|
||
|
各位前辈，小妹新手一枚，学习过程中问题多多。。。最近遇到一个棘手的问题，还请大家多多帮忙。我研究的体系是跨膜环肽纳米管，由于要探索纳米管在脂质膜中的倾斜，所以在模拟的过程中对脂质膜和环肽管没有施加任何限制，然后进行超长时间的MD模拟。 由于计算条件限制，我是以20纳秒为间隔，跑完一次再接着续跑20纳秒。但是实际操作后发现，首次二十纳秒进行顺利，在进行续跑，也就是21到40纳秒的模拟时总是出现错误，错误提示是Periodic cell has become too small for original patch grid!。 我之前对照其他帖子给的方法进行修改，增加"margin 3",改大”langevinPistonPeriod 1000. langevinPistonDecay 500“但是最后还是不行。我在这里先贴出前二十纳秒的CONF和续跑的CONF文件，希望各位前辈帮忙看看问题出在哪里，多谢。 |
回复此楼» 猜你喜欢
288求调剂
已经有14人回复
-
一志愿厦大生物学332求调剂
已经有7人回复
-
085600，321分求调剂
已经有17人回复
-
一志愿鲁东大学071000生物学学硕初试分数276求调剂
已经有16人回复
-
材料化工总分334求调剂
已经有11人回复
-
288求调剂，一志愿华南理工大学071005
已经有13人回复
-
复试调剂，一志愿郑州大学材料与化工289分
已经有20人回复
-
一志愿中科大070300化学，314分求调剂
已经有4人回复
-
285求调剂
已经有9人回复
-
材料调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
4楼2013-11-09 18:24:02
|
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
|
接下来是续跑的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
木有大神看到么。。。
5楼2013-11-11 09:50:01
5