| 查看: 1084 | 回复: 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设置受温度限制吗？有遇到同类的问题的吗？该怎么解决。 |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
【答案】应助回帖
★
感谢参与，应助指数 +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
