24小时热门版块排行榜

>论坛更新日志 (1031)
>导师招生 (728)
>虫友互识 (279)
>休闲灌水 (263)
>招聘信息布告栏 (123)
>博后之家 (67)
>论文投稿 (35)
>教师之家 (29)
>硕博家园 (29)
>考研 (22)
>公派出国 (15)
>学术会议 (13)
>考博 (12)
>生物科学 (10)
>基金申请 (8)
>分子生物 (5)

返回列表

陈闹钟

金虫 (初入文坛)

应助: 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

回复此楼