[AMBER] Reducing positional restraint gradually

From: Rahul Banerjee <rahul.banerjee.chemistry.msu.edu>
Date: Mon, 22 Apr 2013 15:11:44 -0400

Dear Amber users,,

I am running equilibration simulation in presence of a finite positional restraint on all heavy atoms of my system. How can I reduce the restraint force gradually in 1000 steps? My reference structure is the starting structure of the simulation. I am using following script below. At the final step the restraint energy should be zero as the 'k' goes to zero. But it is not so. Am I wrong conceptually?

With best regards,
Rahul
 
Equil
 &cntrl
  imin=0,irest=0,ntx=1,
  nstlim=1000,dt=0.002,
  ntc=2,ntf=2,
  cut=10.0, ntb=2, ntp=1, taup=2.0,
  ntwe=5, ntpr=5, ntwx=5, ntwr=5,
  ntt=3, gamma_ln=2.0, ig=-1
  temp0=300.0, nmropt=1,
  ntr=1, restraintmask=':1-120&!.H=', restraint_wt=8.0,
 /
 &wt
        TYPE='REST', istep1 =0, istep2 = 1000,
        value1 = 8.0, value2 = 0,
 /
 &wt type = 'END'
 /
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Apr 22 2013 - 12:30:02 PDT
Custom Search