Re: [AMBER] Steered MD

From: Pinky Mazumder <pmazumder67.gmail.com>
Date: Thu, 25 Feb 2021 18:13:19 -0600

Hi Carlos,

I changed it to nmropt=1, and it shows,


 RESTRAINTS:
 Requested file redirections:
  DISANG = ref.RST
  DUMPAVE = ref_vs_t
  LISTIN = POUT
  LISTOUT = POUT
 Restraints will be read from file: ref.RST
Here are comments from the DISANG input file:


* Error: Too many atom ranges need to be stored for center-of-mass
distance restraints. MAXGRP = 1. This needs to be increased.*

How to solve this problem? Please any help would be appreciated.

Sincerely,
Pinky

On Thu, Feb 25, 2021 at 2:27 PM Carlos Simmerling <
carlos.simmerling.gmail.com> wrote:

> are you certain that you need jar=1 here? perhaps nmropt=1 would do what
> you need.
>
> On Thu, Feb 25, 2021 at 3:14 PM Pinky Mazumder <pmazumder67.gmail.com>
> wrote:
>
> > Hi,
> >
> > I am trying to pull one layer (group of atoms) towards another group of
> > atoms using the steered molecular dynamics. I have taken the center of
> mass
> > of both groups. The (.RST) and (pull.in) file is given below:
> >
> > &rst
> > iat=-1,-1,
> > r2=72.67,
> > r2a=0.0,
> > rk2=7.2,
> >
> >
> >
> igr1=2709,2730,2751,2772,2793,2814,2835,2856,2877,2898,2919,2940,2961,2982,3003,3024,3048,3069,3090,3111,3132,3153,3174,3195,3216,3237,3258,3279,3300,3321,3342,3363,3386,3407,3428,3449,3470,3491,3512,3533,3554,3575,3596,3617,3638,3659,3680,3701,3725,3746,3767,3788,3809,3830,3851,3872,3893,3914,3935,3956,3977,3998,4019,4040,
> >
> >
> >
> igr2=1365,1386,1407,1428,1449,1470,1491,1512,1533,1554,1575,1596,1617,1638,1659,1680,1704,1725,1746,1767,1788,1809,1830,1851,1872,1893,1914,1935,1956,1977,1998,2019,2042,2063,2084,2105,2126,2147,2168,2189,2210,2231,2252,2273,2294,2315,2336,2357,2381,2402,2423,2444,2465,2486,2507,2528,2549,2570,2591,2612,2633,2654,2675,2696,
> >
> > &end
> >
> > ASMD Simulation Example MD Input
> > &cntrl
> > imin = 0, nstlim = 20000, dt = 0.002,
> > ntx = 1, temp0 = 300.0,
> > ntt = 3, gamma_ln=5.0
> > ntc = 2, ntf = 2, ntb =1,
> > ntwx = 1000, ntwr = 1000, ntpr = 1000,
> > cut = 10.0, ig=-1, ioutfm=1,
> > irest = 0, jar=1,
> > /
> > &wt type='DUMPFREQ', istep1=1, /
> > &wt type='END', /
> > DISANG=ref.RST
> > DUMPAVE=ref_vs_t
> > LISTIN=POUT
> > LISTOUT=POUT
> >
> > When I am running simulation, it gives following error:
> > Restraints will be read from file: ref.RST
> > Here are comments from the DISANG input file:
> >
> >
> > *jar option running In Jarzynsky runs, there must only one restrain,
> > stopping program*
> >
> > Any help would be appreciated.
> >
> > Thank you.
> >
> > Sincerely,
> > Pinky
> > --
> > Pinky Mazumder
> > PhD Student
> > University of Alabama
> > USA
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>


-- 
Pinky Mazumder
PhD Student
University of Alabama
USA
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Feb 25 2021 - 16:30:03 PST
Custom Search