24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

陈闹钟

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 894.5
帖子: 19
在线: 131.1小时
虫号: 2325240
注册: 2013-03-07
性别: MM
专业: 理论和计算化学

[求助] AMBER跑MD想同时加多个restraint,碰到的问题？已有1人参与

我在用AMBER跑MD，想加多个restraint,请问有人知道怎么加吗？我做的输入文件如下：

我按照我这个输入文件加restraint，我看输出的out文件，系统只认最后一行的restraint，没有前几行的restraint。谢谢！

input  file:
&cntrl
  imin = 0, irest = 1, ntx = 7,
  ntb = 2, pres0 = 1.0, ntp = 1,
  taup = 2.0,
  cut = 12.0, ntr = 1,
  ntc = 2, ntf = 2,
  tempi = 300.0, temp0 = 300.0,
  ntt = 3, gamma_ln = 1.0,
  nstlim = 250000, dt = 0.002,
  ntpr = 1000, ntwx = 1000, ntwr = 1000,
restraintmask = '(:62 >@10.0) & !@H='  restraint_wt = 12.0
restraintmask = '((:62 >@9.0) & (:62 <@10.0)) & !@H='  restraint_wt = 10.0
restraintmask = '((:62 >@8.0) & (:62 <@9.0)) & !@H='  restraint_wt = 8.0
restraintmask = '((:62 >@7.0) & (:62 <@8.0)) & !@H='  restraint_wt = 6.0
restraintmask = '((:62 >@6.0) & (:62 <@7.0)) & !@H='  restraint_wt = 4.0
restraintmask = '((:62 >@5.0) & (:62 <@6.0)) & !@H='  restraint_wt = 2.0,
/
Hold the 8HQ fixed
0.0
RES 1 225
END
END

回复此楼