[AMBER] Problem related to the umbrella sampling

From: Sindrila Dutta banik <sindrila.duttabanik.yahoo.com>
Date: Thu, 11 Jul 2013 12:44:35 +0800 (SGT)

Dear all,

I am trying to do umbrella sampling considering two reaction coordinates: increase one N-H bond length (from 1.1A to 1.8A) and decrease one O-H bond (from 1.8 to 1.1A) (which are part of the substrate PPD). To do the same I use the following input for equilibration step. 

10 ps equilibration run
&cntrl
  imin   = 0,
  irest  = 0,
  ntx    = 1,
  ntb    = 1,
  ntc    = 2,
  noshakemask = ":PPD"
  ntf    = 2,
  cut    = 10,
  tempi  = 300.0,
  temp0  = 300.0,
  ntt    = 3,
  gamma_ln = 1.0,
  nstlim = 10000, dt = 0.001
  ntpr = 100, ntwx =100, ntwr =200
  nmropt = 1, ioutfm=1
 /
&wt
  type='END',
 &end
DISANG=comb7

The comb7 file: 

Two restraints N-H and O-H
 &rst
  iat=12230,12245
  r1=0.0, r2=1.7, r3=1.7, r4=10.0,
  rk2=1000, rk3=1000
 &end

 &rst
  iat=12219,12245
  r1=0.0, r2=1.2, r3=1.2, r4=10.0,
  rk2=1000, rk3=1000,
 &end

When I increase N-H bond length from 1.1 upto 1.6, there was no problem. When I set it 1.7 as given in the above input file it gives the following error.   

Coordinate resetting (SHAKE) cannot be accomplished,
     deviation is too large
     NITER, NIT, LL, I and J are :      0      3   4945  10018  10019

     Note: This is usually a symptom of some deeper
     problem with the energetics of the system.

How to solve this problem. Please let me know. 

With best regards
Sindrila
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jul 10 2013 - 22:00:03 PDT
Custom Search