| 查看: 1036 | 回复: 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设置受温度限制吗？有遇到同类的问题的吗？该怎么解决。 |
回复此楼» 猜你喜欢
265求调剂
已经有9人回复
-
085600材料与化工调剂
已经有20人回复
-
专硕 351 086100 也是考的材科基 本科也是材料
已经有6人回复
-
085600专硕材料与化工348分求调剂
已经有10人回复
-
085600 295分求调剂
已经有21人回复
-
285求调剂
已经有5人回复
-
一志愿安徽大学0817化学工程与技术，求调剂
已经有9人回复
-
一志愿0817化学工程与技术，求调剂
已经有8人回复
-
271分求调剂学校
已经有3人回复
-
生物学308分求调剂（一志愿华东师大）
已经有7人回复
【答案】应助回帖
★
感谢参与，应助指数 +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
