| 查看: 1063 | 回复: 2 | ||
[求助]
江湖救急：在amber里跑md轨迹时，遇到很严重的问题，求大神来指教。 已有2人参与
|
|
我的目的是，跑一个md轨迹，我的输入文件是： dhfr evb &cntrl ievb = 1, irest = 0, ntx=1, ntb = 1, ntt = 3, gamma_ln=100.0, cut = 10.0, nstlim = 40000, dt = 0.0025, ntpr = 40, ntave=1, ntwr = 400, ntwx = 400, ntwe = 400, nscm = 100, temp0 = 300., ntf = 1, ntc = 1, vlimit=20.0, iwrap=1, / &ewald nfft1 = 120, nfft2 = 120, nfft3 = 120, / ~ nstlim*dt 是轨迹模拟总时间，老板的要求是nstlim不变，总模拟时间要跑个100ps，200ps，500ps，1ns。 跑出来的结果用vmd看时 : 导入prmtop文件，选择amber 7 parm，导入mdcrd，选择 AMBER Coordinates with periodic box, 在这个输入文件里，其他变量都不变，我发现 dt =0.001时，在vmd里正常，当 dt=0.002时，以及更大时，在vmd里看 发现会出现乱七八糟的成键信息，水盒子也乱了，模拟的蛋白质也乱成键。请问这个dt设置受温度限制吗？有遇到同类的问题的吗？该怎么解决。 |
回复此楼» 猜你喜欢
河北省自然科学基金
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
论文撤稿了
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
青B发送上会通知了吗
已经有7人回复
-
博士申请
已经有3人回复
【答案】应助回帖
★
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流！ 2014-06-17 15:26:35
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流！ 2014-06-17 15:26:35
|
是不是equilibration没有做好。 按照http://ambermd.org/tutorials/advanced/tutorial8/loop1.htm的步骤重复一次看看。 |
2楼2014-06-07 18:20:45
★
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流！ 2014-06-17 15:26:39
感谢参与，应助指数 +1
月只蓝: 金币+1, 鼓励交流！ 2014-06-17 15:26:39
|
本帖内容被屏蔽 |
3楼2014-06-09 16:11:30
10