| 查看: 3061 | 回复: 15 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【讨论】大家好，加水后跑完分子动力学，在log中发现了问题 已有4人参与
|
||||
|
加水后跑完分子动力学，在log中发现了一些问题，有很多warning和一些error，希望大家能帮我看看是不是我参数设计的问题还是我加的水时不合适造成的体系问题，谢谢各位了，期待大家的回复： 这是conf文件： ############################################################# ## JOB DESCRIPTION ## ############################################################# # Minimization and Equilibration of # 1kimz in a Water Box ############################################################# ## ADJUSTABLE PARAMETERS ## ############################################################# structure ../common/1kimz_wb.psf coordinates ../common/1kimz_wb.pdb set temperature 310 set outputname 1kimz_wb_eq firsttimestep 0 ############################################################# ## SIMULATION PARAMETERS ## ############################################################# # Input paraTypeCharmm on parameters ../common/par_all27_prot_lipid.inp temperature $temperature # Force-Field Parameters exclude scaled1-4 1-4scaling 1.0 cutoff 10. switching on switchdist 8. pairlistdist 12 # Integrator Parameters timestep 2.0 ;# 1fs/step rigidBonds all ;# needed for 2fs steps nonbondedFreq 1 fullElectFrequency 2 stepspercycle 10 # Constant Temperature Control langevin on ;# do langevin dynamics langevinDamping 5 ;# damping coefficient (gamma) of 5/ps langevinTemp $temperature langevinHydrogen off ;# don't couple langevin bath to hydrogens # Periodic Boundary Conditions cellBasisVector1 80. 0. 0. cellBasisVector2 0. 70. 0. cellBasisVector3 0. 0 75. cellOrigin 29. 85. 42. wrapAll on # PME (for full-system periodic electrostatics) PME yes PMEGridSizeX 82 PMEGridSizeY 71 PMEGridSizeZ 78 # Constant Pressure Control (variable volume) useGroupPressure yes ;# needed for rigidBonds useFlexibleCell no useConstantArea no langevinPiston on langevinPistonTarget 1.01325 ;# in bar -> 1 atm langevinPistonPeriod 100. langevinPistonDecay 50. langevinPistonTemp $temperature # Output outputName $outputname restartfreq 500 ;# 500steps = every 1ps dcdfreq 250 xstFreq 250 outputEnergies 100 outputPressure 100 ############################################################# ## EXTRA PARAMETERS ## ############################################################# ############################################################# ## EXECUTION SCRIPT ## ############################################################# # Minimization minimize 1000 reinitvels $temperature run 5000 ;# 100ps |
回复此楼» 收录本帖的淘帖专辑推荐
 精华帖转载 |
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教大家氧化硅果冻状凝胶加水后还能稀释成溶胶吗
已经有6人回复
- 水性聚氨酯合成加水乳化后出现絮状物，没有得带乳液
已经有7人回复
- 【求助】加水后 碳纳米管不见了（已解决）
已经有4人回复
- 【求助】将聚异丁烯丁二酸酐溶于二甲苯中，加水洗涤立马有白色状沉淀是为什么？
已经有5人回复
- 【求助/交流】蒽酮试剂加水后有沉淀……
已经有6人回复
- 【求助】聚氨酯乳化先加水还是先加乙二胺
已经有9人回复
7楼2010-12-12 19:32:50
|
弄到凌晨了，问题也没解决。按照坛子里出现类似问题的帖子进行尝试，发现： （1）把步长变小，行不通，还是出现同样的问题： warning：Bad global exclusion count, possible error warning: Increasing cutoff during minimization may avoid this. ERROR: Constraint failure in RATTLE algorithm for atom 787!.... Constraint failure :simulation has become unstable （2）按照出现的错误： warning：Bad global exclusion count, possible error warning: Increasing cutoff during minimization may avoid this. 我把cutoff的值变大，NAMD的光标不删了，log也不更新，也就是NAMD直接罢工了（选了两个值试都不行） 我把cutoff的值变小，也不行同样出现： warning：Bad global exclusion count, possible error warning: Increasing cutoff during minimization may avoid this. （3）其他的原因：psf文件有问题？？我是菜鸟也看不出来，在产生psf文件的log里没发现有任何问题。 （4）原子重叠，哎，我都不知道怎样查看，笨。 我把我的psf与pdb放到网盘里了（咱们坛子传不上来）：http://u.115.com/?ac=my#ct=frame ... id=1&cid=998954 用户名：tonme 密码：002725 （5）期待大家的指点，谢谢了。 |
2楼2010-12-12 17:43:56
|
最新进展： 缩小步长不行，还是出现同样的问题。 现在我把cutoff改为2或1，发现每当算到第528步的时候就算不下去了，log文件不更新了。在linux系统终端的那个光标还在闪，不知道这是为什么？？？ 另外： conf文件里的cellBasisVector1/2/3,我一直是按照这个公式进行的： cellBasisVector1 = |min(x)-max(x)|,0,0 cellBasisVector2 = 0,|min(y)-max(y)|,0 cellBasisVector3 = 0,0,|min(z)-max(z)| cellOrigin是写的是center的值， PME我是按照baoding说的那样，cellBasisVector/PMEGridSize<=1,且是2 or 3 or 5 的倍数 我也延长了minimize的时间，可没用，只算到第528步的时候log文件就停了，知识终端的光标还在闪。 大家帮我分析一下，这是什么问题啊？？还有哪里我需要修改的？？先谢谢各位了。 psf pdb 以及conf log文件均在：网盘里了（咱们坛子传不上来）： http://u.115.com/?ac=my#ct=frame ... id=1&cid=998954 用户名：tonme 密码：002725 |
3楼2010-12-12 17:44:09
★ ★
zh1987hs(金币+2):thank you 2010-12-12 18:44:56
zh1987hs(金币+2):thank you 2010-12-12 18:44:56
|
现在知道为什么总是算到第528步就算不下去的原因了，是因为cellOrigin是我写的是center的值，当我帮它改成0,0,0时就可以完全算完了。下一步我想改cutoff试试，能否消去 warning：Bad global exclusion count, possible error warning: Increasing cutoff during minimization may avoid this. ERROR: Constraint failure in RATTLE algorithm for atom 787!.... Constraint failure :simulation has become unstable 这个问题 |
4楼2010-12-12 17:59:28
20
投票:
